Size: px
Start display at page:

Download ""

Transcription

1

2 (MHD ) GFV (Galium Field Visualizer) GFV OpenGL GFV GFV GFV

3 GFV GFV GFV

4 1 MHD (Magnetohydrodynamics) MHD MHD MHD Ultrasonic Velocity Profiler, UVP UVP CG (Visualization) CG MHD MHD MHD 1

5 2 2.1 CG CG API OpenGL Open Graphics Library OpenGL D GPU OpenGL CG OpenGL GLU GLUT GLU OpenGL Utility Library OpenGL CG GLUT OpenGL Utility Toolkit GLUT OpenGL gltranslate gluperspective glutmousefunc [1] 2.2 [2] φ X(x,y,z) Fig. 1 X x y z [i][j][k] (x i, y j, z k ) φ ijk 2

6 φ w ijk w ijk = f x i f y j f z k f x i = x i+1 x = x i+1 x x i+1 x i x f x i+1 = 1 f x i f y j = y j+1 y = y j+1 y y j+1 y j y f y j+1 = 1 f y j fk z = z k+1 z = z k+1 z z k+1 z k z f z k+1 = 1 f z k Fig. 1 X φ φ(x, y, z) = w 000 φ w 001 φ w 010 φ w 011 φ w 100 φ w 101 φ w 110 φ w 111 φ [3] dx dt = v(x) 3

7 Fig. 1: Location of φ. x x n t x x n+1 x 0 = x n, k 1 = t v(x 0 ), k 2 = t v(x 0 + k 1 2 ), k 3 = t v(x 0 + k 2 2 ), k 4 = t v(x 0 + k 3 ), x n+1 = x (k 1 + 2k 2 + 2k 3 + k 4 ). dx dt = v x(x, y, z) dy dt = v y(x, y, z) dz dt = v z(x, y, z) 4

8 x 0 = x n, y 0 = y n, z 0 = z n k1 x = t v x (x 0, y 0, z 0 ) k y 1 = t v y (x 0, y 0, z 0 ) k1 z = t v z (x 0, y 0, z 0 ) k x 2 = t v x (x 0 + kx 1 2, y 0 + ky 1 2, z 0 + kz 1 2 ) k y 2 = t v y (x 0 + kx 1 2, y 0 + ky 1 2, z 0 + kz 1 2 ) k z 2 = t v z (x 0 + kx 1 2, y 0 + ky 1 2, z 0 + kz 1 2 ) k x 3 = t v x (x 0 + kx 2 2, y 0 + ky 2 2, z 0 + kz 2 2 ) k y 3 = t v y (x 0 + kx 2 2, y 0 + ky 2 2, z 0 + kz 2 2 ) k z 3 = t v z (x 0 + kx 2 2, y 0 + ky 2 2, z 0 + kz 2 2 ) k x 4 = t v x (x 0 + k x 3, y 0 + k y 3, z 0 + k z 3) k y 4 = t v y (x 0 + k x 3, y 0 + k y 3, z 0 + k z 3) k z 4 = t v z (x 0 + k x 3, y 0 + k y 3, z 0 + k z 3) x n+1 = x (kx 1 + 2k x 2 + 2k x 3 + k x 4) y n+1 = y (ky 1 + 2k y 2 + 2k y 3 + k y 4) z n+1 = z (kz 1 + 2k z 2 + 2k z 3 + k z 4) 5

9 3 MHD C OpenGL GFV (Galium Field Visualizer) 3.1 GFV x y z GFV Fig. 2 x y z GFV NX,NY,NZ X,Y,Z XMIN,XMAX X YMIN,YMAX Y ZMIN,ZMAX Z PARTICLEMAX 6

10 Fig. 2: Bounding box. x y z 7

11 MESH 3.2 Fig. 3 [Slice] [Arrow] [ParticleTracer] [SelectSliceData] [SelectArrowData] Reset View End Fig. 3: Example of the menu of GFV. 8

12 Fortran90 C C Fortran x y z (NX,NY,NZ) C [0][0][0] [0][0][1] [0][0][2] Fortran (1,1,1) (2,1,1) (3,1,1) Fortran C C [NZ][NY][NX] C C Fortran NX*NY*NZ Fig. 4: Multidimensional arrays in Fortran and C

13 0x123456AB 16 [ AB] [AB ] CPU (union) union { } data ; f l o a t f ; unsigned char adr [ s i z e o f ( f l o a t ) ] ; adr data.adr[n] 1 float GFV x y z [Slice] [SelectSliceData] Velocity( ) Pressure( ) Enstrophy( ) Helicity( ) xy Shift + yz Ctrl + zx Shift + 10

14 OpenGL glbegin(gl TRIANGLE STRIP) glend() xy Fig. 5 y x OpenGL xy yz zx (Fig. 6) Fig. 5: A slice plane by OpenGL s triangle strip. xy 11

15 Fig. 6: A snapshot of slice planes that visualize flow velocity amplitude. 12

16 5 GFV 5.1 (Fig. 7) xy yz zx x y y z z x [Arrow] Slicexy xy Sliceyz yz Slicezx zx xy x y 50%, 50% 5.2 x y z (Fig.8 Fig. 9) [Arrow] print3darrow 13

17 x y z 50% 50% atan2 x y xy z 70% 5.3 [Particle Tracer] set xy z zx y yz x [Particle Tracer] start particle.txt x y z [ParticleTracer] input set [ParticleTracer] output particle.txt 14

18 Fig. 7: A snapshot of velocity vector visualization by 2D arrows. Fig. 8: A snapshot of velocity vector visualization by 3D arrows. Fig. 7 Fig. 8 15

19 Fig. 9: Another snapshot of velocity vector visualization by 3D arrows. 2.3 x y z 2.2 start GLUT glutidlefunc glutidlefunc 16

20 Fig. 10: A snapshot of velocity vector visualization by animated particle tracers. 17

21 6 GFV 6.1 (UVP) Fig mm 200mm 200mm Fig. 12 x 6.2 [5] Fig. 13 5:5:1 x y z xz yz (xy ) Ra Ek Ra = (g/t 0)βd 4 (κ/ρ 0 C p )(ν/ρ 0 ) Ek = ν ΩH 2 [4] NX=251 NY=251 NZ=51 Ra = Ek =

22 Fig. 11: Liquid galium contained in a convection vessel. 6.3 GFV GFV NX=251 NY=251 NZ=51 x XMIN XMAX 1.25 YMIN YMAX 1.25 ZMIN ZMAX Fig. 14 x yz x Fig. 15 xy x x Fig. 16 yz yz x=0 Fig

23 Fig. 12: Helmholtz coils for horizontal magnetic field applied to the convection vessel. Fig. 13: Simulation model. 20

24 Fig. 18 Fig. 18 x Fig. 17 Fig. 19 Fig. 19 Fig. 20 GFV Fig. 21 x Fig. 22 x x Fig. 22 x Fig. 23 ω = v (ω ω) [6] Fig. 23 Fig. 24 v ω [7] 0 Fig. 24 y y Fig.25 GFV Fig

25 Fig. 14: A snapshot of slice planes for velocity field amplitude of liquid galium convection. x 22

26 Fig. 15: Velocity vectors by 2D arrows. xy yz 23

27 Fig. 16: Velocity vectors visualized by 2D and 3D arrows. yz x 24

28 Fig. 17: Distribution of magnetic field amplitude in the liquid galium convection. Fig. 18: Magnetic field vectors visualized by 2D and 3D arrows. 25

29 Fig. 19: Combined visualization of 3D arrows and slice planes. 26

30 Fig. 20: Magnetic field concentration that is visualized by the red stripes in the horizontal slice planes. 27

31 Fig. 21: Helical flows structure visualized by particle tracers. Fig. 22: Torus structure of the flow in convection rolls. 28

32 Fig. 23: The enstrophy distribution on slice planes with helical flow lines. Fig. 24: Two Different kinds of helical flow lines; right-winding and left winding. The color shows the helicity density distribution; positive in red and negative in blue. 29

33 Fig. 25: A pair of opposite windings of helical flow, or a pair of opposite signs of helicity density, in each convection roll. Fig. 26: Pressure on slice planes. 30

34 7 MHD GFV x GFV x Fig. 16 x x x x GFV GFV 31

35 [1], OpenGL CG, (2003) [2],, (2005) [3] W.H.Press, B.P.Flannery, S.A.Teukolsky and W.T.Vetterling, Numerical Recipes in C : The Art of Scientific Computing, Second Edition,Cambridge Press, pp (1992) [4] D.O.Gough,D.R.Moore,E.A.Spiegel,N.O.Weiss, CONVECTIVE INSTABIL- ITY IN A COMPRESSIBLE ATOMOSPHERE. II, THE ASTROPHYSI- CAL JOURNAL, 206, pp (1976) [5],, (2010) [6],,,pp (1982) [7],,,p.162(1998) MHD 32

応力とひずみ.ppt

応力とひずみ.ppt in yukawa@numse.nagoya-u.ac.jp 2 3 4 5 x 2 6 Continuum) 7 8 9 F F 10 F L L F L 1 L F L F L F 11 F L F F L F L L L 1 L 2 12 F L F! A A! S! = F S 13 F L L F F n = F " cos# F t = F " sin# S $ = S cos# S S

More information

A

A A05-132 2010 2 11 1 1 3 1.1.......................................... 3 1.2..................................... 3 1.3..................................... 3 2 4 2.1............................... 4 2.2

More information

変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy,

変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy, 変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy, z + dz) Q! (x + d x + u + du, y + dy + v + dv, z +

More information

13,825,228 3,707,995 26.8 4.9 25 3 8 9 1 50,000 0.29 1.59 70,000 0.29 1.74 12,500 0.39 1.69 12,500 0.55 10,000 20,000 0.13 1.58 30,000 0.00 1.26 5,000 0.13 1.58 25,000 40,000 0.13 1.58 50,000 0.00 1.26

More information

II A A441 : October 02, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka )

II A A441 : October 02, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka ) II 214-1 : October 2, 214 Version : 1.1 Kawahira, Tomoki TA (Kondo, Hirotaka ) http://www.math.nagoya-u.ac.jp/~kawahira/courses/14w-biseki.html pdf 1 2 1 9 1 16 1 23 1 3 11 6 11 13 11 2 11 27 12 4 12 11

More information

.1 z = e x +xy y z y 1 1 x 0 1 z x y α β γ z = αx + βy + γ (.1) ax + by + cz = d (.1') a, b, c, d x-y-z (a, b, c). x-y-z 3 (0,

.1 z = e x +xy y z y 1 1 x 0 1 z x y α β γ z = αx + βy + γ (.1) ax + by + cz = d (.1') a, b, c, d x-y-z (a, b, c). x-y-z 3 (0, .1.1 Y K L Y = K 1 3 L 3 L K K (K + ) 1 1 3 L 3 K 3 L 3 K 0 (K + K) 1 3 L 3 K 1 3 L 3 lim K 0 K = L (K + K) 1 3 K 1 3 3 lim K 0 K = 1 3 K 3 L 3 z = f(x, y) x y z x-y-z.1 z = e x +xy y 3 x-y ( ) z 0 f(x,

More information

III No (i) (ii) (iii) (iv) (v) (vi) x 2 3xy + 2 lim. (x,y) (1,0) x 2 + y 2 lim (x,y) (0,0) lim (x,y) (0,0) lim (x,y) (0,0) 5x 2 y x 2 + y 2. xy x2 + y

III No (i) (ii) (iii) (iv) (v) (vi) x 2 3xy + 2 lim. (x,y) (1,0) x 2 + y 2 lim (x,y) (0,0) lim (x,y) (0,0) lim (x,y) (0,0) 5x 2 y x 2 + y 2. xy x2 + y III No (i) (ii) (iii) (iv) (v) (vi) x 2 3xy + 2. (x,y) (1,0) x 2 + y 2 5x 2 y x 2 + y 2. xy x2 + y 2. 2x + y 3 x 2 + y 2 + 5. sin(x 2 + y 2 ). x 2 + y 2 sin(x 2 y + xy 2 ). xy (i) (ii) (iii) 2xy x 2 +

More information

食品工学.indb

食品工学.indb , Vol. 11, No. 1, pp. 51-58, Mar. 2010 IH The Visualization and Quantification of the Flow in the Pan During Induction Heating and Gas Heating Haruna KAWAKAMI, Zensyu TOU, Mika FUKUOKA, and Noboru SAKAI

More information

II 1 3 2 5 3 7 4 8 5 11 6 13 7 16 8 18 2 1 1. x 2 + xy x y (1 lim (x,y (1,1 x 1 x 3 + y 3 (2 lim (x,y (, x 2 + y 2 x 2 (3 lim (x,y (, x 2 + y 2 xy (4 lim (x,y (, x 2 + y 2 x y (5 lim (x,y (, x + y x 3y

More information

液晶の物理1:連続体理論(弾性,粘性)

液晶の物理1:連続体理論(弾性,粘性) The Physics of Liquid Crystals P. G. de Gennes and J. Prost (Oxford University Press, 1993) Liquid crystals are beautiful and mysterious; I am fond of them for both reasons. My hope is that some readers

More information

18 2 F 12 r 2 r 1 (3) Coulomb km Coulomb M = kg F G = ( ) ( ) ( ) 2 = [N]. Coulomb

18 2 F 12 r 2 r 1 (3) Coulomb km Coulomb M = kg F G = ( ) ( ) ( ) 2 = [N]. Coulomb r 1 r 2 r 1 r 2 2 Coulomb Gauss Coulomb 2.1 Coulomb 1 2 r 1 r 2 1 2 F 12 2 1 F 21 F 12 = F 21 = 1 4πε 0 1 2 r 1 r 2 2 r 1 r 2 r 1 r 2 (2.1) Coulomb ε 0 = 107 4πc 2 =8.854 187 817 10 12 C 2 N 1 m 2 (2.2)

More information

平成 22 年度 ( 第 32 回 ) 数学入門公開講座テキスト ( 京都大学数理解析研究所, 平成 ~8 22 月年 58 日開催月 2 日 ) V := {(x,y) x n + y n 1 = 0}, W := {(x,y,z) x 3 yz = x 2 y z 2

平成 22 年度 ( 第 32 回 ) 数学入門公開講座テキスト ( 京都大学数理解析研究所, 平成 ~8 22 月年 58 日開催月 2 日 ) V := {(x,y) x n + y n 1 = 0}, W := {(x,y,z) x 3 yz = x 2 y z 2 3 90 2006 1. V := {(x,y) x n + y n 1 = 0}, W := {(x,y,z) x 3 yz = x 2 y z 2 = xz y 2 = 0} V (x,y) n = 1 n = 2 (x,y) V n = 1 n = 2 (3/5,4/5),(5/13,12/13)... n 3 V (0,±1),(±1,0) ( ) n 3 x n + y n = z n,

More information

x () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x

x () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x [ ] IC. f(x) = e x () f(x) f (x) () lim f(x) lim f(x) x + x (3) lim f(x) lim f(x) x + x (4) y = f(x) ( ) ( s46). < a < () a () lim a log xdx a log xdx ( ) n (3) lim log k log n n n k=.3 z = log(x + y ),

More information

1.2 y + P (x)y + Q(x)y = 0 (1) y 1 (x), y 2 (x) y 1 (x), y 2 (x) (1) y(x) c 1, c 2 y(x) = c 1 y 1 (x) + c 2 y 2 (x) 3 y 1 (x) y 1 (x) e R P (x)dx y 2

1.2 y + P (x)y + Q(x)y = 0 (1) y 1 (x), y 2 (x) y 1 (x), y 2 (x) (1) y(x) c 1, c 2 y(x) = c 1 y 1 (x) + c 2 y 2 (x) 3 y 1 (x) y 1 (x) e R P (x)dx y 2 1 1.1 R(x) = 0 y + P (x)y + Q(x)y = R(x)...(1) y + P (x)y + Q(x)y = 0...(2) 1 2 u(x) v(x) c 1 u(x)+ c 2 v(x) = 0 c 1 = c 2 = 0 c 1 = c 2 = 0 2 0 2 u(x) v(x) u(x) u (x) W (u, v)(x) = v(x) v (x) 0 1 1.2

More information

ma22-9 u ( v w) = u v w sin θê = v w sin θ u cos φ = = 2.3 ( a b) ( c d) = ( a c)( b d) ( a d)( b c) ( a b) ( c d) = (a 2 b 3 a 3 b 2 )(c 2 d 3 c 3 d

ma22-9 u ( v w) = u v w sin θê = v w sin θ u cos φ = = 2.3 ( a b) ( c d) = ( a c)( b d) ( a d)( b c) ( a b) ( c d) = (a 2 b 3 a 3 b 2 )(c 2 d 3 c 3 d A 2. x F (t) =f sin ωt x(0) = ẋ(0) = 0 ω θ sin θ θ 3! θ3 v = f mω cos ωt x = f mω (t sin ωt) ω t 0 = f ( cos ωt) mω x ma2-2 t ω x f (t mω ω (ωt ) 6 (ωt)3 = f 6m ωt3 2.2 u ( v w) = v ( w u) = w ( u v) ma22-9

More information

. ev=,604k m 3 Debye ɛ 0 kt e λ D = n e n e Ze 4 ln Λ ν ei = 5.6π / ɛ 0 m/ e kt e /3 ν ei v e H + +e H ev Saha x x = 3/ πme kt g i g e n

. ev=,604k m 3 Debye ɛ 0 kt e λ D = n e n e Ze 4 ln Λ ν ei = 5.6π / ɛ 0 m/ e kt e /3 ν ei v e H + +e H ev Saha x x = 3/ πme kt g i g e n 003...............................3 Debye................. 3.4................ 3 3 3 3. Larmor Cyclotron... 3 3................ 4 3.3.......... 4 3.3............ 4 3.3...... 4 3.3.3............ 5 3.4.........

More information

II (10 4 ) 1. p (x, y) (a, b) ε(x, y; a, b) 0 f (x, y) f (a, b) A, B (6.5) y = b f (x, b) f (a, b) x a = A + ε(x, b; a, b) x a 2 x a 0 A = f x (

II (10 4 ) 1. p (x, y) (a, b) ε(x, y; a, b) 0 f (x, y) f (a, b) A, B (6.5) y = b f (x, b) f (a, b) x a = A + ε(x, b; a, b) x a 2 x a 0 A = f x ( II (1 4 ) 1. p.13 1 (x, y) (a, b) ε(x, y; a, b) f (x, y) f (a, b) A, B (6.5) y = b f (x, b) f (a, b) x a = A + ε(x, b; a, b) x a x a A = f x (a, b) y x 3 3y 3 (x, y) (, ) f (x, y) = x + y (x, y) = (, )

More information

k m m d2 x i dt 2 = f i = kx i (i = 1, 2, 3 or x, y, z) f i σ ij x i e ij = 2.1 Hooke s law and elastic constants (a) x i (2.1) k m σ A σ σ σ σ f i x

k m m d2 x i dt 2 = f i = kx i (i = 1, 2, 3 or x, y, z) f i σ ij x i e ij = 2.1 Hooke s law and elastic constants (a) x i (2.1) k m σ A σ σ σ σ f i x k m m d2 x i dt 2 = f i = kx i (i = 1, 2, 3 or x, y, z) f i ij x i e ij = 2.1 Hooke s law and elastic constants (a) x i (2.1) k m A f i x i B e e e e 0 e* e e (2.1) e (b) A e = 0 B = 0 (c) (2.1) (d) e

More information

JKR Point loading of an elastic half-space 2 3 Pressure applied to a circular region Boussinesq, n =

JKR Point loading of an elastic half-space 2 3 Pressure applied to a circular region Boussinesq, n = JKR 17 9 15 1 Point loading of an elastic half-space Pressure applied to a circular region 4.1 Boussinesq, n = 1.............................. 4. Hertz, n = 1.................................. 6 4 Hertz

More information

1 No.1 5 C 1 I III F 1 F 2 F 1 F 2 2 Φ 2 (t) = Φ 1 (t) Φ 1 (t t). = Φ 1(t) t = ( 1.5e 0.5t 2.4e 4t 2e 10t ) τ < 0 t > τ Φ 2 (t) < 0 lim t Φ 2 (t) = 0

1 No.1 5 C 1 I III F 1 F 2 F 1 F 2 2 Φ 2 (t) = Φ 1 (t) Φ 1 (t t). = Φ 1(t) t = ( 1.5e 0.5t 2.4e 4t 2e 10t ) τ < 0 t > τ Φ 2 (t) < 0 lim t Φ 2 (t) = 0 1 No.1 5 C 1 I III F 1 F 2 F 1 F 2 2 Φ 2 (t) = Φ 1 (t) Φ 1 (t t). = Φ 1(t) t = ( 1.5e 0.5t 2.4e 4t 2e 10t ) τ < 0 t > τ Φ 2 (t) < 0 lim t Φ 2 (t) = 0 0 < t < τ I II 0 No.2 2 C x y x y > 0 x 0 x > b a dx

More information

A

A A 2563 15 4 21 1 3 1.1................................................ 3 1.2............................................. 3 2 3 2.1......................................... 3 2.2............................................

More information

II 2014 2 (1) log(1 + r/100) n = log 2 n log(1 + r/100) = log 2 n = log 2 log(1 + r/100) (2) y = f(x) = log(1 + x) x = 0 1 f (x) = 1/(1 + x) f (0) = 1

II 2014 2 (1) log(1 + r/100) n = log 2 n log(1 + r/100) = log 2 n = log 2 log(1 + r/100) (2) y = f(x) = log(1 + x) x = 0 1 f (x) = 1/(1 + x) f (0) = 1 II 2014 1 1 I 1.1 72 r 2 72 8 72/8 = 9 9 2 a 0 1 a 1 a 1 = a 0 (1+r/100) 2 a 2 a 2 = a 1 (1 + r/100) = a 0 (1 + r/100) 2 n a n = a 0 (1 + r/100) n a n a 0 2 n a 0 (1 + r/100) n = 2a 0 (1 + r/100) n = 2

More information

pdf

pdf http://www.ns.kogakuin.ac.jp/~ft13389/lecture/physics1a2b/ pdf I 1 1 1.1 ( ) 1. 30 m µm 2. 20 cm km 3. 10 m 2 cm 2 4. 5 cm 3 km 3 5. 1 6. 1 7. 1 1.2 ( ) 1. 1 m + 10 cm 2. 1 hr + 6400 sec 3. 3.0 10 5 kg

More information

, = = 7 6 = 42, =

, = = 7 6 = 42, = http://www.ss.u-tokai.ac.jp/~mahoro/2016autumn/alg_intro/ 1 1 2016.9.26, http://www.ss.u-tokai.ac.jp/~mahoro/2016autumn/alg_intro/ 1.1 1 214 132 = 28258 2 + 1 + 4 1 + 3 + 2 = 7 6 = 42, 4 + 2 = 6 2 + 8

More information

i

i i 3 4 4 7 5 6 3 ( ).. () 3 () (3) (4) /. 3. 4/3 7. /e 8. a > a, a = /, > a >. () a >, a =, > a > () a > b, a = b, a < b. c c n a n + b n + c n 3c n..... () /3 () + (3) / (4) /4 (5) m > n, a b >, m > n,

More information

6. Euler x

6. Euler x ...............................................................................3......................................... 4.4................................... 5.5......................................

More information

i

i 009 I 1 8 5 i 0 1 0.1..................................... 1 0.................................................. 1 0.3................................. 0.4........................................... 3

More information

n Y 1 (x),..., Y n (x) 1 W (Y 1 (x),..., Y n (x)) 0 W (Y 1 (x),..., Y n (x)) = Y 1 (x)... Y n (x) Y 1(x)... Y n(x) (x)... Y n (n 1) (x) Y (n 1)

n Y 1 (x),..., Y n (x) 1 W (Y 1 (x),..., Y n (x)) 0 W (Y 1 (x),..., Y n (x)) = Y 1 (x)... Y n (x) Y 1(x)... Y n(x) (x)... Y n (n 1) (x) Y (n 1) D d dx 1 1.1 n d n y a 0 dx n + a d n 1 y 1 dx n 1 +... + a dy n 1 dx + a ny = f(x)...(1) dk y dx k = y (k) a 0 y (n) + a 1 y (n 1) +... + a n 1 y + a n y = f(x)...(2) (2) (2) f(x) 0 a 0 y (n) + a 1 y

More information

δ ij δ ij ˆx ˆx ŷ ŷ ẑ ẑ 0, ˆx ŷ ŷ ˆx ẑ, ŷ ẑ ẑ ŷ ẑ, ẑ ˆx ˆx ẑ ŷ, a b a x ˆx + a y ŷ + a z ẑ b x ˆx + b

δ ij δ ij ˆx ˆx ŷ ŷ ẑ ẑ 0, ˆx ŷ ŷ ˆx ẑ, ŷ ẑ ẑ ŷ ẑ, ẑ ˆx ˆx ẑ ŷ, a b a x ˆx + a y ŷ + a z ẑ b x ˆx + b 23 2 2.1 n n r x, y, z ˆx ŷ ẑ 1 a a x ˆx + a y ŷ + a z ẑ 2.1.1 3 a iˆx i. 2.1.2 i1 i j k e x e y e z 3 a b a i b i i 1, 2, 3 x y z ˆx i ˆx j δ ij, 2.1.3 n a b a i b i a i b i a x b x + a y b y + a z b

More information

DVIOUT

DVIOUT A. A. A-- [ ] f(x) x = f 00 (x) f 0 () =0 f 00 () > 0= f(x) x = f 00 () < 0= f(x) x = A--2 [ ] f(x) D f 00 (x) > 0= y = f(x) f 00 (x) < 0= y = f(x) P (, f()) f 00 () =0 A--3 [ ] y = f(x) [, b] x = f (y)

More information

p = mv p x > h/4π λ = h p m v Ψ 2 Ψ

p = mv p x > h/4π λ = h p m v Ψ 2 Ψ II p = mv p x > h/4π λ = h p m v Ψ 2 Ψ Ψ Ψ 2 0 x P'(x) m d 2 x = mω 2 x = kx = F(x) dt 2 x = cos(ωt + φ) mω 2 = k ω = m k v = dx = -ωsin(ωt + φ) dt = d 2 x dt 2 0 y v θ P(x,y) θ = ωt + φ ν = ω [Hz] 2π

More information

1

1 PalmGauss SC PGSC-5G Instruction Manual PalmGauss SC PGSC-5G Version 1.01 PalmGauss SC PGSC5G 1.... 3 2.... 3 3.... 3 3.1... 3 3.2... 3 3.3 PalmGauss... 4 3.4... 4 3.4.1 (Fig. 4)... 4 3.4.2 (Fig. 5)...

More information

y π π O π x 9 s94.5 y dy dx. y = x + 3 y = x logx + 9 s9.6 z z x, z y. z = xy + y 3 z = sinx y 9 s x dx π x cos xdx 9 s93.8 a, fx = e x ax,. a =

y π π O π x 9 s94.5 y dy dx. y = x + 3 y = x logx + 9 s9.6 z z x, z y. z = xy + y 3 z = sinx y 9 s x dx π x cos xdx 9 s93.8 a, fx = e x ax,. a = [ ] 9 IC. dx = 3x 4y dt dy dt = x y u xt = expλt u yt λ u u t = u u u + u = xt yt 6 3. u = x, y, z = x + y + z u u 9 s9 grad u ux, y, z = c c : grad u = u x i + u y j + u k i, j, k z x, y, z grad u v =

More information

第62巻 第1号 平成24年4月/石こうを用いた木材ペレット

第62巻 第1号 平成24年4月/石こうを用いた木材ペレット Bulletin of Japan Association for Fire Science and Engineering Vol. 62. No. 1 (2012) Development of Two-Dimensional Simple Simulation Model and Evaluation of Discharge Ability for Water Discharge of Firefighting

More information

I, II 1, A = A 4 : 6 = max{ A, } A A 10 10%

I, II 1, A = A 4 : 6 = max{ A, } A A 10 10% 1 2006.4.17. A 3-312 tel: 092-726-4774, e-mail: hara@math.kyushu-u.ac.jp, http://www.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html Office hours: B A I ɛ-δ ɛ-δ 1. 2. A 1. 1. 2. 3. 4. 5. 2. ɛ-δ 1. ɛ-n

More information

第3章 OpenGL の基礎

第3章 OpenGL の基礎 3 OpenGL April 20, 2012 1 / 23 31 ( ) OpenGL OpenGL 2 / 23 32 OpenGL OpenGL OpenGL (Open Graphics Library) Silicon Graphics, Inc 2 3 API (Application Program Interface) [4] UNIX OS Windows Macintosh CAD

More information

(3) (2),,. ( 20) ( s200103) 0.7 x C,, x 2 + y 2 + ax = 0 a.. D,. D, y C, C (x, y) (y 0) C m. (2) D y = y(x) (x ± y 0), (x, y) D, m, m = 1., D. (x 2 y

(3) (2),,. ( 20) ( s200103) 0.7 x C,, x 2 + y 2 + ax = 0 a.. D,. D, y C, C (x, y) (y 0) C m. (2) D y = y(x) (x ± y 0), (x, y) D, m, m = 1., D. (x 2 y [ ] 7 0.1 2 2 + y = t sin t IC ( 9) ( s090101) 0.2 y = d2 y 2, y = x 3 y + y 2 = 0 (2) y + 2y 3y = e 2x 0.3 1 ( y ) = f x C u = y x ( 15) ( s150102) [ ] y/x du x = Cexp f(u) u (2) x y = xey/x ( 16) ( s160101)

More information

第3章 OpenGL の基礎

第3章 OpenGL の基礎 3 OpenGL April 11, 2017 1 / 28 3.1 ( ) OpenGL OpenGL 2 / 28 3.2 OpenGL OpenGL OpenGL (Open Graphics Library) Silicon Graphics, Inc. 2 3 API (Application Program Interface) [4] UNIX OS Windows Macintosh

More information

, 3, STUDY ON IMPORTANCE OF OPTIMIZED GRID STRUCTURE IN GENERAL COORDINATE SYSTEM 1 2 Hiroyasu YASUDA and Tsuyoshi HOSHINO

, 3, STUDY ON IMPORTANCE OF OPTIMIZED GRID STRUCTURE IN GENERAL COORDINATE SYSTEM 1 2 Hiroyasu YASUDA and Tsuyoshi HOSHINO , 3, 2012 9 STUDY ON IMPORTANCE OF OPTIMIZED GRID STRUCTURE IN GENERAL COORDINATE SYSTEM 1 2 Hiroyasu YASUDA and Tsuyoshi HOSHINO 1 950-2181 2 8050 2 950-2181 2 8050 Numerical computation of river flows

More information

120 9 I I 1 I 2 I 1 I 2 ( a) ( b) ( c ) I I 2 I 1 I ( d) ( e) ( f ) 9.1: Ampère (c) (d) (e) S I 1 I 2 B ds = µ 0 ( I 1 I 2 ) I 1 I 2 B ds =0. I 1 I 2

120 9 I I 1 I 2 I 1 I 2 ( a) ( b) ( c ) I I 2 I 1 I ( d) ( e) ( f ) 9.1: Ampère (c) (d) (e) S I 1 I 2 B ds = µ 0 ( I 1 I 2 ) I 1 I 2 B ds =0. I 1 I 2 9 E B 9.1 9.1.1 Ampère Ampère Ampère s law B S µ 0 B ds = µ 0 j ds (9.1) S rot B = µ 0 j (9.2) S Ampère Biot-Savart oulomb Gauss Ampère rot B 0 Ampère µ 0 9.1 (a) (b) I B ds = µ 0 I. I 1 I 2 B ds = µ 0

More information

Gmech08.dvi

Gmech08.dvi 145 13 13.1 13.1.1 0 m mg S 13.1 F 13.1 F /m S F F 13.1 F mg S F F mg 13.1: m d2 r 2 = F + F = 0 (13.1) 146 13 F = F (13.2) S S S S S P r S P r r = r 0 + r (13.3) r 0 S S m d2 r 2 = F (13.4) (13.3) d 2

More information

1.1 1 A

1.1 1 A . A..2 2 2. () (xyz) ( xyz) ( xy z) = (x x)yz ( xy z) = yz ( xy z) = y(z ( x z)) = y((z x)(z z)) = y( x z) (2) (3) M aj (x, y, M aj ( x, ȳ, z)) = xy ȳm aj ( x, ȳ, z) M aj ( x, ȳ, z)x M aj (x, y, z) x =

More information

7 27 7.1........................................ 27 7.2.......................................... 28 1 ( a 3 = 3 = 3 a a > 0(a a a a < 0(a a a -1 1 6

7 27 7.1........................................ 27 7.2.......................................... 28 1 ( a 3 = 3 = 3 a a > 0(a a a a < 0(a a a -1 1 6 26 11 5 1 ( 2 2 2 3 5 3.1...................................... 5 3.2....................................... 5 3.3....................................... 6 3.4....................................... 7

More information

x, y x 3 y xy 3 x 2 y + xy 2 x 3 + y 3 = x 3 y xy 3 x 2 y + xy 2 x 3 + y 3 = 15 xy (x y) (x + y) xy (x y) (x y) ( x 2 + xy + y 2) = 15 (x y)

x, y x 3 y xy 3 x 2 y + xy 2 x 3 + y 3 = x 3 y xy 3 x 2 y + xy 2 x 3 + y 3 = 15 xy (x y) (x + y) xy (x y) (x y) ( x 2 + xy + y 2) = 15 (x y) x, y x 3 y xy 3 x 2 y + xy 2 x 3 + y 3 = 15 1 1977 x 3 y xy 3 x 2 y + xy 2 x 3 + y 3 = 15 xy (x y) (x + y) xy (x y) (x y) ( x 2 + xy + y 2) = 15 (x y) ( x 2 y + xy 2 x 2 2xy y 2) = 15 (x y) (x + y) (xy

More information

untitled

untitled 2004/12/21 2/2 (11/16) DT-MRI (11/30) /OpenGL 12/7 12/14 (12/21) 1/11 (1/18) OpenGL ~ ~ OpenGL Silicon Graphics, OpenGL ~ ~ OpenGL OpenGL Utility Library (GLU) OpenGL. OpenGL. OpenGL Utility Toolkit (GLUT)

More information

7. y fx, z gy z gfx dz dx dz dy dy dx. g f a g bf a b fa 7., chain ule Ω, D R n, R m a Ω, f : Ω R m, g : D R l, fω D, b fa, f a g b g f a g f a g bf a

7. y fx, z gy z gfx dz dx dz dy dy dx. g f a g bf a b fa 7., chain ule Ω, D R n, R m a Ω, f : Ω R m, g : D R l, fω D, b fa, f a g b g f a g f a g bf a 9 203 6 7 WWW http://www.math.meiji.ac.jp/~mk/lectue/tahensuu-203/ 2 8 8 7. 7 7. y fx, z gy z gfx dz dx dz dy dy dx. g f a g bf a b fa 7., chain ule Ω, D R n, R m a Ω, f : Ω R m, g : D R l, fω D, b fa,

More information

,. Black-Scholes u t t, x c u 0 t, x x u t t, x c u t, x x u t t, x + σ x u t, x + rx ut, x rux, t 0 x x,,.,. Step 3, 7,,, Step 6., Step 4,. Step 5,,.

,. Black-Scholes u t t, x c u 0 t, x x u t t, x c u t, x x u t t, x + σ x u t, x + rx ut, x rux, t 0 x x,,.,. Step 3, 7,,, Step 6., Step 4,. Step 5,,. 9 α ν β Ξ ξ Γ γ o δ Π π ε ρ ζ Σ σ η τ Θ θ Υ υ ι Φ φ κ χ Λ λ Ψ ψ µ Ω ω Def, Prop, Th, Lem, Note, Remark, Ex,, Proof, R, N, Q, C [a, b {x R : a x b} : a, b {x R : a < x < b} : [a, b {x R : a x < b} : a,

More information

『共形場理論』

『共形場理論』 T (z) SL(2, C) T (z) SU(2) S 1 /Z 2 SU(2) (ŜU(2) k ŜU(2) 1)/ŜU(2) k+1 ŜU(2)/Û(1) G H N =1 N =1 N =1 N =1 N =2 N =2 N =2 N =2 ĉ>1 N =2 N =2 N =4 N =4 1 2 2 z=x 1 +ix 2 z f(z) f(z) 1 1 4 4 N =4 1 = = 1.3

More information

() n C + n C + n C + + n C n n (3) n C + n C + n C 4 + n C + n C 3 + n C 5 + (5) (6 ) n C + nc + 3 nc n nc n (7 ) n C + nc + 3 nc n nc n (

() n C + n C + n C + + n C n n (3) n C + n C + n C 4 + n C + n C 3 + n C 5 + (5) (6 ) n C + nc + 3 nc n nc n (7 ) n C + nc + 3 nc n nc n ( 3 n nc k+ k + 3 () n C r n C n r nc r C r + C r ( r n ) () n C + n C + n C + + n C n n (3) n C + n C + n C 4 + n C + n C 3 + n C 5 + (4) n C n n C + n C + n C + + n C n (5) k k n C k n C k (6) n C + nc

More information

75 unit: mm Fig. Structure of model three-phase stacked transformer cores (a) Alternate-lap joint (b) Step-lap joint 3 4)

75 unit: mm Fig. Structure of model three-phase stacked transformer cores (a) Alternate-lap joint (b) Step-lap joint 3 4) 3 * 35 (3), 7 Analysis of Local Magnetic Properties and Acoustic Noise in Three-Phase Stacked Transformer Core Model Masayoshi Ishida Kenichi Sadahiro Seiji Okabe 3.7 T 5 Hz..4 3 Synopsis: Methods of local

More information

21 2 26 i 1 1 1.1............................ 1 1.2............................ 3 2 9 2.1................... 9 2.2.......... 9 2.3................... 11 2.4....................... 12 3 15 3.1..........

More information

日立金属技報 Vol.34

日立金属技報 Vol.34 Influence of Misorientation Angle between Adjacent Grains on Magnetization Reversal in Nd-Fe-B Sintered Magnet Tomohito Maki Rintaro Ishii Mitsutoshi Natsumeda Takeshi Nishiuchi Ryo Uchikoshi Masaaki Takezawa

More information

MUFFIN3

MUFFIN3 MUFFIN - MUltiFarious FIeld simulator for Non-equilibrium system - ( ) MUFFIN WG3 - - JCII, - ( ) - ( ) - ( ) - (JSR) - - MUFFIN sec -3 msec -6 sec GOURMET SUSHI MUFFIN -9 nsec PASTA -1 psec -15 fsec COGNAC

More information

2012専門分科会_new_4.pptx

2012専門分科会_new_4.pptx d dt L L = 0 q i q i d dt L L = 0 r i i r i r r + Δr Δr δl = 0 dl dt = d dt i L L q i q i + q i i q i = q d L L i + q i i dt q i i q i = i L L q i L = 0, H = q q i L = E i q i i d dt L q q i i L = L(q

More information

2.4 ( ) ( B ) A B F (1) W = B A F dr. A F q dr f(x,y,z) A B Γ( ) Minoru TANAKA (Osaka Univ.) I(2011), Sec p. 1/30

2.4 ( ) ( B ) A B F (1) W = B A F dr. A F q dr f(x,y,z) A B Γ( ) Minoru TANAKA (Osaka Univ.) I(2011), Sec p. 1/30 2.4 ( ) 2.4.1 ( B ) A B F (1) W = B A F dr. A F q dr f(x,y,z) A B Γ( ) I(2011), Sec. 2. 4 p. 1/30 (2) Γ f dr lim f i r i. r i 0 i f i i f r i i i+1 (1) n i r i (3) F dr = lim F i n i r i. Γ r i 0 i n i

More information

D v D F v/d F v D F η v D (3.2) (a) F=0 (b) v=const. D F v Newtonian fluid σ ė σ = ηė (2.2) ė kl σ ij = D ijkl ė kl D ijkl (2.14) ė ij (3.3) µ η visco

D v D F v/d F v D F η v D (3.2) (a) F=0 (b) v=const. D F v Newtonian fluid σ ė σ = ηė (2.2) ė kl σ ij = D ijkl ė kl D ijkl (2.14) ė ij (3.3) µ η visco post glacial rebound 3.1 Viscosity and Newtonian fluid f i = kx i σ ij e kl ideal fluid (1.9) irreversible process e ij u k strain rate tensor (3.1) v i u i / t e ij v F 23 D v D F v/d F v D F η v D (3.2)

More information

2005

2005 20 30 8 3 190 60 A,B 67,2000 98 20 23,600 100 60 10 20 1 3 2 1 2 1 12 1 1 ( ) 340 20 20 30 50 50 ( ) 6 80 5 65 17 21 5 5 12 35 1 5 20 3 3,456,871 2,539,950 916,921 18 10 29 5 3 JC-V 2 ( ) 1 17 3 1 6

More information

II 2 II

II 2 II II 2 II 2005 yugami@cc.utsunomiya-u.ac.jp 2005 4 1 1 2 5 2.1.................................... 5 2.2................................. 6 2.3............................. 6 2.4.................................

More information

7-12.dvi

7-12.dvi 26 12 1 23. xyz ϕ f(x, y, z) Φ F (x, y, z) = F (x, y, z) G(x, y, z) rot(grad ϕ) rot(grad f) H(x, y, z) div(rot Φ) div(rot F ) (x, y, z) rot(grad f) = rot f x f y f z = (f z ) y (f y ) z (f x ) z (f z )

More information

II No.01 [n/2] [1]H n (x) H n (x) = ( 1) r n! r!(n 2r)! (2x)n 2r. r=0 [2]H n (x) n,, H n ( x) = ( 1) n H n (x). [3] H n (x) = ( 1) n dn x2 e dx n e x2

II No.01 [n/2] [1]H n (x) H n (x) = ( 1) r n! r!(n 2r)! (2x)n 2r. r=0 [2]H n (x) n,, H n ( x) = ( 1) n H n (x). [3] H n (x) = ( 1) n dn x2 e dx n e x2 II No.1 [n/] [1]H n x) H n x) = 1) r n! r!n r)! x)n r r= []H n x) n,, H n x) = 1) n H n x) [3] H n x) = 1) n dn x e dx n e x [4] H n+1 x) = xh n x) nh n 1 x) ) d dx x H n x) = H n+1 x) d dx H nx) = nh

More information

微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.

微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 初版 1 刷発行時のものです. 微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. ttp://www.morikita.co.jp/books/mid/00571 このサンプルページの内容は, 初版 1 刷発行時のものです. i ii 014 10 iii [note] 1 3 iv 4 5 3 6 4 x 0 sin x x 1 5 6 z = f(x, y) 1 y = f(x)

More information

EGunGPU

EGunGPU Super Computing in Accelerator simulations - Electron Gun simulation using GPGPU - K. Ohmi, KEK-Accel Accelerator Physics seminar 2009.11.19 Super computers in KEK HITACHI SR11000 POWER5 16 24GB 16 134GFlops,

More information

2 1,384,000 2,000,000 1,296,211 1,793,925 38,000 54,500 27,804 43,187 41,000 60,000 31,776 49,017 8,781 18,663 25,000 35,300 3 4 5 6 1,296,211 1,793,925 27,804 43,187 1,275,648 1,753,306 29,387 43,025

More information

x ( ) x dx = ax

x ( ) x dx = ax x ( ) x dx = ax 1 dx = a x log x = at + c x(t) = e at C (C = e c ) a > 0 t a < 0 t 0 (at + b ) h dx = lim x(t + h) x(t) h 0 h x(t + h) x(t) h x(t) t x(t + h) x(t) ax(t) h x(t + h) x(t) + ahx(t) 0, h, 2h,

More information

No δs δs = r + δr r = δr (3) δs δs = r r = δr + u(r + δr, t) u(r, t) (4) δr = (δx, δy, δz) u i (r + δr, t) u i (r, t) = u i x j δx j (5) δs 2

No δs δs = r + δr r = δr (3) δs δs = r r = δr + u(r + δr, t) u(r, t) (4) δr = (δx, δy, δz) u i (r + δr, t) u i (r, t) = u i x j δx j (5) δs 2 No.2 1 2 2 δs δs = r + δr r = δr (3) δs δs = r r = δr + u(r + δr, t) u(r, t) (4) δr = (δx, δy, δz) u i (r + δr, t) u i (r, t) = u i δx j (5) δs 2 = δx i δx i + 2 u i δx i δx j = δs 2 + 2s ij δx i δx j

More information

Optical Lenses CCD Camera Laser Sheet Wind Turbine with med Diffuser Pitot Tube PC Fig.1 Experimental facility. Transparent Diffuser Double Pulsed Nd:

Optical Lenses CCD Camera Laser Sheet Wind Turbine with med Diffuser Pitot Tube PC Fig.1 Experimental facility. Transparent Diffuser Double Pulsed Nd: *1 *2 *3 PIV Measurement of Field of the Wind Turbine with a med Diffuser Kazuhiko TOSHIMITSU *4, Koutarou NISHIKAWA and Yuji OHYA *4 Department of Mechanical Engineering, Matsue National Collage of Technology,

More information

PDF

PDF 1 1 1 1-1 1 1-9 1-3 1-1 13-17 -3 6-4 6 3 3-1 35 3-37 3-3 38 4 4-1 39 4- Fe C TEM 41 4-3 C TEM 44 4-4 Fe TEM 46 4-5 5 4-6 5 5 51 6 5 1 1-1 1991 1,1 multiwall nanotube 1993 singlewall nanotube ( 1,) sp 7.4eV

More information

ii

ii ii iii 1 1 1.1..................................... 1 1.2................................... 3 1.3........................... 4 2 9 2.1.................................. 9 2.2...............................

More information

Appendix A BASIC BASIC Beginner s All-purpose Symbolic Instruction Code FORTRAN COBOL C JAVA PASCAL (NEC N88-BASIC Windows BASIC (1) (2) ( ) BASIC BAS

Appendix A BASIC BASIC Beginner s All-purpose Symbolic Instruction Code FORTRAN COBOL C JAVA PASCAL (NEC N88-BASIC Windows BASIC (1) (2) ( ) BASIC BAS Appendix A BASIC BASIC Beginner s All-purpose Symbolic Instruction Code FORTRAN COBOL C JAVA PASCAL (NEC N88-BASIC Windows BASIC (1 (2 ( BASIC BASIC download TUTORIAL.PDF http://hp.vector.co.jp/authors/va008683/

More information

( ) ( )

( ) ( ) 20 21 2 8 1 2 2 3 21 3 22 3 23 4 24 5 25 5 26 6 27 8 28 ( ) 9 3 10 31 10 32 ( ) 12 4 13 41 0 13 42 14 43 0 15 44 17 5 18 6 18 1 1 2 2 1 2 1 0 2 0 3 0 4 0 2 2 21 t (x(t) y(t)) 2 x(t) y(t) γ(t) (x(t) y(t))

More information

1 1 sin cos P (primary) S (secondly) 2 P S A sin(ω2πt + α) A ω 1 ω α V T m T m 1 100Hz m 2 36km 500Hz. 36km 1

1 1 sin cos P (primary) S (secondly) 2 P S A sin(ω2πt + α) A ω 1 ω α V T m T m 1 100Hz m 2 36km 500Hz. 36km 1 sin cos P (primary) S (secondly) 2 P S A sin(ω2πt + α) A ω ω α 3 3 2 2V 3 33+.6T m T 5 34m Hz. 34 3.4m 2 36km 5Hz. 36km m 34 m 5 34 + m 5 33 5 =.66m 34m 34 x =.66 55Hz, 35 5 =.7 485.7Hz 2 V 5Hz.5V.5V V

More information

福岡大学人文論叢47-3

福岡大学人文論叢47-3 679 pp. 1 680 2 681 pp. 3 682 4 683 5 684 pp. 6 685 7 686 8 687 9 688 pp. b 10 689 11 690 12 691 13 692 pp. 14 693 15 694 a b 16 695 a b 17 696 a 18 697 B 19 698 A B B B A B B A A 20 699 pp. 21 700 pp.

More information

, 1.,,,.,., (Lin, 1955).,.,.,.,. f, 2,. main.tex 2011/08/13( )

, 1.,,,.,., (Lin, 1955).,.,.,.,. f, 2,. main.tex 2011/08/13( ) 81 4 2 4.1, 1.,,,.,., (Lin, 1955).,.,.,.,. f, 2,. 82 4.2. ζ t + V (ζ + βy) = 0 (4.2.1), V = 0 (4.2.2). (4.2.1), (3.3.66) R 1 Φ / Z, Γ., F 1 ( 3.2 ). 7,., ( )., (4.2.1) 500 hpa., 500 hpa (4.2.1) 1949,.,

More information

空力騒音シミュレータの開発

空力騒音シミュレータの開発 41 COSMOS-V, an Aerodynamic Noise Simulator Nariaki Horinouchi COSMOS-V COSMOS-V COSMOS-V 3 The present and future computational problems of the aerodynamic noise analysis using COSMOS-V, our in-house

More information

6kg 1.1m 1.m.1m.1 l λ ϵ λ l + λ l l l dl dl + dλ ϵ dλ dl dl + dλ dl dl 3 1. JIS 1 6kg 1% 66kg 1 13 σ a1 σ m σ a1 σ m σ m σ a1 f f σ a1 σ a1 σ m f 4

6kg 1.1m 1.m.1m.1 l λ ϵ λ l + λ l l l dl dl + dλ ϵ dλ dl dl + dλ dl dl 3 1. JIS 1 6kg 1% 66kg 1 13 σ a1 σ m σ a1 σ m σ m σ a1 f f σ a1 σ a1 σ m f 4 35-8585 7 8 1 I I 1 1.1 6kg 1m P σ σ P 1 l l λ λ l 1.m 1 6kg 1.1m 1.m.1m.1 l λ ϵ λ l + λ l l l dl dl + dλ ϵ dλ dl dl + dλ dl dl 3 1. JIS 1 6kg 1% 66kg 1 13 σ a1 σ m σ a1 σ m σ m σ a1 f f σ a1 σ a1 σ m

More information

( ) Loewner SLE 13 February

( ) Loewner SLE 13 February ( ) Loewner SLE 3 February 00 G. F. Lawler, Conformally Invariant Processes in the Plane, (American Mathematical Society, 005)., Summer School 009 (009 8 7-9 ) . d- (BES d ) d B t = (Bt, B t,, Bd t ) (d

More information

211 kotaro@math.titech.ac.jp 1 R *1 n n R n *2 R n = {(x 1,..., x n ) x 1,..., x n R}. R R 2 R 3 R n R n R n D D R n *3 ) (x 1,..., x n ) f(x 1,..., x n ) f D *4 n 2 n = 1 ( ) 1 f D R n f : D R 1.1. (x,

More information

? FPGA FPGA FPGA : : : ? ( ) (FFT) ( ) (Localization) ? : 0. 1 2 3 0. 4 5 6 7 3 8 6 1 5 4 9 2 0. 0 5 6 0 8 8 ( ) ? : LU Ax = b LU : Ax = 211 410 221 x 1 x 2 x 3 = 1 0 0 21 1 2 1 0 0 1 2 x = LUx = b 1 31

More information

.5 z = a + b + c n.6 = a sin t y = b cos t dy d a e e b e + e c e e e + e 3 s36 3 a + y = a, b > b 3 s363.7 y = + 3 y = + 3 s364.8 cos a 3 s365.9 y =,

.5 z = a + b + c n.6 = a sin t y = b cos t dy d a e e b e + e c e e e + e 3 s36 3 a + y = a, b > b 3 s363.7 y = + 3 y = + 3 s364.8 cos a 3 s365.9 y =, [ ] IC. r, θ r, θ π, y y = 3 3 = r cos θ r sin θ D D = {, y ; y }, y D r, θ ep y yddy D D 9 s96. d y dt + 3dy + y = cos t dt t = y = e π + e π +. t = π y =.9 s6.3 d y d + dy d + y = y =, dy d = 3 a, b

More information

i 18 2H 2 + O 2 2H 2 + ( ) 3K

i 18 2H 2 + O 2 2H 2 + ( ) 3K i 18 2H 2 + O 2 2H 2 + ( ) 3K ii 1 1 1.1.................................. 1 1.2........................................ 3 1.3......................................... 3 1.4....................................

More information

( 12 ( ( ( ( Levi-Civita grad div rot ( ( = 4 : 6 3 1 1.1 f(x n f (n (x, d n f(x (1.1 dxn f (2 (x f (x 1.1 f(x = e x f (n (x = e x d dx (fg = f g + fg (1.2 d dx d 2 dx (fg = f g + 2f g + fg 2... d n n

More information

(u(x)v(x)) = u (x)v(x) + u(x)v (x) ( ) u(x) = u (x)v(x) u(x)v (x) v(x) v(x) 2 y = g(t), t = f(x) y = g(f(x)) dy dx dy dx = dy dt dt dx., y, f, g y = f (g(x))g (x). ( (f(g(x)). ). [ ] y = e ax+b (a, b )

More information

Fermat s Last Theorem Hajime Mashima November 19, 2018 Abstract About 380 years ago, Pierre de Fermat wrote the following idea to Diophantus s Arithme

Fermat s Last Theorem Hajime Mashima November 19, 2018 Abstract About 380 years ago, Pierre de Fermat wrote the following idea to Diophantus s Arithme Fermat s Last Theorem Hajime Mashima November 19, 2018 Abstract About 380 years ago, Pierre de Fermat wrote the following idea to Diophantus s Arithmetica. Cubum autem in duos cubos, aut quadratoquadratum

More information

入試の軌跡

入試の軌跡 4 y O x 4 Typed by L A TEX ε ) ) ) 6 4 ) 4 75 ) http://kumamoto.s.xrea.com/plan/.. PDF) Ctrl +L) Ctrl +) Ctrl + Ctrl + ) ) Alt + ) Alt + ) ESC. http://kumamoto.s.xrea.com/nyusi/kumadai kiseki ri i.pdf

More information

96 7 1m =2 10 7 N 1A 7.1 7.2 a C (1) I (2) A C I A A a A a A A a C C C 7.2: C A C A = = µ 0 2π (1) A C 7.2 AC C A 3 3 µ0 I 2 = 2πa. (2) A C C 7.2 A A

96 7 1m =2 10 7 N 1A 7.1 7.2 a C (1) I (2) A C I A A a A a A A a C C C 7.2: C A C A = = µ 0 2π (1) A C 7.2 AC C A 3 3 µ0 I 2 = 2πa. (2) A C C 7.2 A A 7 Lorentz 7.1 Ampère I 1 I 2 I 2 I 1 L I 1 I 2 21 12 L r 21 = 12 = µ 0 2π I 1 I 2 r L. (7.1) 7.1 µ 0 =4π 10 7 N A 2 (7.2) magnetic permiability I 1 I 2 I 1 I 2 12 21 12 21 7.1: 1m 95 96 7 1m =2 10 7 N

More information

dx dt = f x,t ( ) t

dx dt = f x,t ( ) t MATLAB 1. 2. Runge-Kutta 3. 1. 2. 3. dx dt = f x,t ( ) t dx( t) dt = lim Δt x( t + Δt) x( t) Δt t = nδt n = 0,1,2,3,4,5, x = x( nδt) n Δt dx ( ) dt = f x,t x n +1 x n Δt = f ( x,nδt) n 1 x n = x 0 n =

More information

73

73 73 74 ( u w + bw) d = Ɣ t tw dɣ u = N u + N u + N 3 u 3 + N 4 u 4 + [K ] {u = {F 75 u δu L σ (L) σ dx σ + dσ x δu b δu + d(δu) ALW W = L b δu dv + Aσ (L)δu(L) δu = (= ) W = A L b δu dx + Aσ (L)δu(L) Aσ

More information

(1.2) T D = 0 T = D = 30 kn 1.2 (1.4) 2F W = 0 F = W/2 = 300 kn/2 = 150 kn 1.3 (1.9) R = W 1 + W 2 = = 1100 N. (1.9) W 2 b W 1 a = 0

(1.2) T D = 0 T = D = 30 kn 1.2 (1.4) 2F W = 0 F = W/2 = 300 kn/2 = 150 kn 1.3 (1.9) R = W 1 + W 2 = = 1100 N. (1.9) W 2 b W 1 a = 0 1 1 1.1 1.) T D = T = D = kn 1. 1.4) F W = F = W/ = kn/ = 15 kn 1. 1.9) R = W 1 + W = 6 + 5 = 11 N. 1.9) W b W 1 a = a = W /W 1 )b = 5/6) = 5 cm 1.4 AB AC P 1, P x, y x, y y x 1.4.) P sin 6 + P 1 sin 45

More information

WebGL OpenGL GLSL Kageyama (Kobe Univ.) Visualization / 57

WebGL OpenGL GLSL Kageyama (Kobe Univ.) Visualization / 57 WebGL 2014.04.15 X021 2014 3 1F Kageyama (Kobe Univ.) Visualization 2014.04.15 1 / 57 WebGL OpenGL GLSL Kageyama (Kobe Univ.) Visualization 2014.04.15 2 / 57 WebGL Kageyama (Kobe Univ.) Visualization 2014.04.15

More information

course pptx

course pptx ParaView () 20105141CAE OPENFOAM (R) is a registered trade mark of OpenCFD Limited, the producer of the OpenFOAM software and owner of the OPENFOAM (R) and OpenCFD (R) trade marks. This offering is not

More information

() (, y) E(, y) () E(, y) (3) q ( ) () E(, y) = k q q (, y) () E(, y) = k r r (3).3 [.7 ] f y = f y () f(, y) = y () f(, y) = tan y y ( ) () f y = f y

() (, y) E(, y) () E(, y) (3) q ( ) () E(, y) = k q q (, y) () E(, y) = k r r (3).3 [.7 ] f y = f y () f(, y) = y () f(, y) = tan y y ( ) () f y = f y 5. [. ] z = f(, y) () z = 3 4 y + y + 3y () z = y (3) z = sin( y) (4) z = cos y (5) z = 4y (6) z = tan y (7) z = log( + y ) (8) z = tan y + + y ( ) () z = 3 8y + y z y = 4 + + 6y () z = y z y = (3) z =

More information

I ( ) 1 de Broglie 1 (de Broglie) p λ k h Planck ( Js) p = h λ = k (1) h 2π : Dirac k B Boltzmann ( J/K) T U = 3 2 k BT

I ( ) 1 de Broglie 1 (de Broglie) p λ k h Planck ( Js) p = h λ = k (1) h 2π : Dirac k B Boltzmann ( J/K) T U = 3 2 k BT I (008 4 0 de Broglie (de Broglie p λ k h Planck ( 6.63 0 34 Js p = h λ = k ( h π : Dirac k B Boltzmann (.38 0 3 J/K T U = 3 k BT ( = λ m k B T h m = 0.067m 0 m 0 = 9. 0 3 kg GaAs( a T = 300 K 3 fg 07345

More information

210 資料 TI 89 (1) TI 89 2nd ON HOME ( ) ( ) HOME =! ENTER ( ) = (10) ENTER ( ) [ ] { } ( )! 2 =! ( ) ( ) 2 3x ( 2y + yz) ( ) 3x ( ( ) 2y + y z)

210 資料 TI 89 (1) TI 89 2nd ON HOME ( ) ( ) HOME =! ENTER ( ) = (10) ENTER ( ) [ ] { } ( )! 2 =! ( ) ( ) 2 3x ( 2y + yz) ( ) 3x ( ( ) 2y + y z) 210 資料 TI 89 (1) TI 89 2nd ON HOME () () HOME =! ENTER () = 3 10 3 (10) ENTER ( ) [ ] { } ( )! 2 =! ( ) () 2 3x ( 2y + yz) ( ) 3x ( ( ) 2y + y z) ENTER () 2nd 9 2nd 9) ENTER ( ) 2nd 7) ENTER 7 7 ) ENTER

More information

M3 x y f(x, y) (= x) (= y) x + y f(x, y) = x + y + *. f(x, y) π y f(x, y) x f(x + x, y) f(x, y) lim x x () f(x,y) x 3 -

M3 x y f(x, y) (= x) (= y) x + y f(x, y) = x + y + *. f(x, y) π y f(x, y) x f(x + x, y) f(x, y) lim x x () f(x,y) x 3 - M3............................................................................................ 3.3................................................... 3 6........................................... 6..........................................

More information