: LES Energy RANS Kawai & Lele JCP (2010) Wavenumber Accuracy LES, DNS RANS 1970 s Efficiency
|
|
|
- えいじろう つちかね
- 7 years ago
- Views:
Transcription
1 MHD (ISAS), (JAXA) [email protected] hp://flab.eng.isas.jaxa.jp/member/kawai JAXA
2 : LES Energy RANS Kawai & Lele JCP (2010) Wavenumber Accuracy LES, DNS RANS 1970 s Efficiency
3 FDS, FVS, AUSM: 1970 s Q n+1 j = Q n+1 j + x E j+1/2 E j 1/2 E j+1/2 = 1 2 [E j+1 + E j A ave (Q j+1 Q j )] Turbulence: K-H insabiliy w/u Exac FDS High-order cenral x/r c
4 : non-dissipaive : FDS, FVS, AUSM 5 4 Densiy Reference FDS x WENO (Jiang & Shu JCP 1996) WCNS (Deng & Maekawa JCP 1997) Universiy of Chicago LAD (Kawai & Lele JCP 2008,2010) Soshi Kawai, ISAS/JAXA
5 ESA NASA B =0 NASA
6 MHD B =0 Eigh-wave formulaion Powell NASA repor (1994) - Projecion Brackbill & Barnes JCP (1980) - B =0 Consrained Transpor - 2 Evans & Hawley APJ (1988) B =0
7 LAD Kawai e al. JCP (2008, 2010) u + ( uu + p ar ( u) ) =0 u u
8 LAD Kawai e al. JCP (2008, 2010) u + ( uu + p ar ( u) ) =0 u u Kawai & Lele AIAA J. (2010) B B B B
9 B B B B B x 0 A z A y B y = RHS phyb + A z 0 A x B z A y A x 0 A = A x A y = ar B = ar J A z
10 B B B B B x B y = RHS phyb + B z : : ( B) ( B) =0 0 A z A y A z 0 A x A y A x 0 2 (A z A z ) x y A = ar B = ar J + 2 (A y A y ) x z = J = RHS phyj + 2 (A) ( A) + 2 (A x A x ) y z J = RHS phyj + 2 ( ar J) Soshi Kawai, ISAS/JAXA
11 ar = C ar 1 a 4 J x 4 r J 2 x r x r+3 B x B y = RHS phyb + B z J A z A y A z 0 A x A y A x 0 A = ar B = ar J x J ar PDE ( B) =0 (D B) n+1 =(D B)n =0
12 B (n[i,j,k]) ( B) n = 1 x DL Bx,i [B x]+ 1 y DL By,j [B y]+ 1 z DL Bz,k [B z] B x (n[i,j,k]) B n+1 x, = Bn x, y DL RHS,j [RHS]+D L Az,j [A z] + z DL RHS,k [RHS]+D L Ay,k [A y] 1x DL Bx,i Inducion eqs.: B x B y = RHS phyb + B z 0 A z A y A z 0 A x A y A x 0 1x DL Bx,i B n+1 x, = 1 x DL Bx,i B n x, + x y DL Bx,i x z DL Bx,i D L RHS,j [RHS]+D L Az,j [A z] D L RHS,k [RHS]+D L Ay,k [A y]
13 B (n[i,j,k]) ( B) n+1 =( B)n + x y {F xy} + x z {F xz} + y z {F yz} if F xy = F xz = F yz =0 ( B) n+1 =( B)n =( B)0 F xy = D L Bx,i +D L By,j D L RHS,j [RHS]+D L Az,j [A z] D L RHS,i [RHS]+D L Az,i [A z] if D L Bx,i DL Az,j = DL By,j DL Az,i F xy = F xz = F yz =0
14 D L Bx,i DL Az,j = DL By,j DL Az,i D L Bx,i f = N a, p=m a, a p, f (i+p+1,j,k) f (i+p,j,k) D L Az,j f = D L By,j f = D L Az,i f = N b, q=m b, b q, f (i,j+q+1,k) f (i,j+q,k) N c, q=m c, c q, f (i,j+q+1,k) f (i,j+q,k) N d, p=m d, d p, f (i+p+1,j,k) f (i+p,j,k) D L Bx,i DL Az,j f = N a, p=m a, N b,(i+p+1,j,k) q=m b,(i+p+1,j,k) a p, b q,(i+p+1,j,k) f (i+p+1,j+q+1,k) f (i+p+1,j+q,k) D L By,j DL Az,i f = N c, q=m c, N d,(i,j+q+1,k) p=m d,(i,j+q+1,k) c q, d p,(i,j+q+1,k) f (i+p+1,j+q+1,k) f (i+p,j+q+1,k) a m, = d m,(i,j+n,k) k i f ( OK!) 6 j i Soshi Kawai, ISAS/JAXA
15 PDE B x B y = RHS phyb + B z 0 A z A y A z 0 A x A y A x 0 4 J A = ar B = ar J ar x 4 PDE OK 6 k i f j Soshi Kawai, ISAS/JAXA i
16 : C 1) ar = C 1 a r J 2 x r x r+3 B x B y = RHS phyb + B z 0 A z A y A z 0 A x A y A x 0 A = ar B = ar J 2) : Alfven 3) : 4) : Orszag-Tang
17 C 1D slow swich-off shock problem Falle e al. MNRAS (1998) ar = C 1 a r J 2 x r x r+3 J = RHS phyj + 2 ( ar J) By x C =0 C =1 C = 10 C = 100 Wiggles ampliude(%) Thickness Wiggles C Shock hickness (grid poins)
18 Alfven Exac N=4 N=8 N=16 = 45 2N B N x N=8 B = N=16 Soshi Kawai, ISAS/JAXA Toh JCP (2000)
19 h-order slope Error Y L error + L 2 error X
20 Dai-Woodward 1.6 = 45 Dai & Woodward JCP (1994) Lef condiions, u,u,u z,p,b,b,b z L =(1.08, 1.2, 0.01, 0.5, 0.95, Righ condiions 2/ 4, 3.6/ 4, 2/ 4 ), u,u,u z,p,b,b,b z R =(1, 0, 0, 0, 1, 2/ 4, 4/ 4, 2/ 4 ) B D (N=4000) 1D HLLD (N=800) Miyoshi & Kusano JCP (2005) 2D presen (N=800) x x B 0 Soshi Kawai, ISAS/JAXA
21 Orszag-Tang Orszag & Tang JFM (1979) Tempraure Soshi Kawai, ISAS/JAXA
22 Orszag-Tang Orszag & Tang JFM (1979) Temperaure Arificial diffusion A = ar J
23 Orszag-Tang Orszag & Tang JFM (1979) Temperaure Divergence B field B = 10 12
24 y =0.5 1 N=1600 N=200 N=400 N=800 Temperaure T = p/ ar Arificial diffusion x x Soshi Kawai, ISAS/JAXA
25 MHD - - PDE - =0 pressure u-velociy X1 X2 X3 X1 X2 X3 MHD =20 =40 =60 =80 =100 Soshi Kawai, ISAS/JAXA
untitled
20 7 1 22 7 1 1 2 3 7 8 9 10 11 13 14 15 17 18 19 21 22 - 1 - - 2 - - 3 - - 4 - 50 200 50 200-5 - 50 200 50 200 50 200 - 6 - - 7 - () - 8 - (XY) - 9 - 112-10 - - 11 - - 12 - - 13 - - 14 - - 15 - - 16 -
untitled
19 1 19 19 3 8 1 19 1 61 2 479 1965 64 1237 148 1272 58 183 X 1 X 2 12 2 15 A B 5 18 B 29 X 1 12 10 31 A 1 58 Y B 14 1 25 3 31 1 5 5 15 Y B 1 232 Y B 1 4235 14 11 8 5350 2409 X 1 15 10 10 B Y Y 2 X 1 X
untitled
CFD JAXA CFD CFD Navier-Stokes -1- -2-1 CFD CFD Navier-Stokes 1. 2. 2,500,000 CFD Copyright Boeing CFD Boeing B777 HP -3- -4-2 CFD European Transonic Wind tunnel (-163 ) JAXA 2mx2m CFD (Computational Fluid
宿泊産業活性化のための実証実験
121 32 10 12 12 19 2 15 59 40 33 34 35 36 37 38 3637 20 39 12 19 OFF 2008/12/19 2008/12/25 3 1 1 72,000 2008/12/19 2008/12/26 2 1 1 36,000 2008/12/28 2009/1/5 2 1 1 24,000 2009/1/6 2009/1/16 3 1 1 25,200
901 902 2 40 5 786 30 2 2 100 10100200 903 904 2 3 2 12 905 6765 30 3 61016 1 10162532 253240 2 2 1 2 100 24 45 545 1 2 2 510 1515 1010 50 300 0 10 2942 560 2 1 1 2 24 15 2565 2 10 2942 560 3 3 56 03 18
1037 1038 2 40 5 876 30 2 2 100 10100200 1039 1040 2 3 2 12 1041 6765 30 1 1 2 2 1 2 100 24 45 545 1 2 2 510 1515 1010 50 300 0 10 2942 560 2 1 3 2 10 2942 560 3 61016 1 10162532 253240 1 2 24 15 2565
MHD) MHD MHD 1977 MHD ApJ MHD simulation Abstract/Keywords ADS simulation ApJ MHD simulation 1982
MHD) MHD MHD 1977 MHD ApJ MHD simulation Abstract/Keywords 1982 86 5 23 1987 91 57 1992 1996 127 1997 2001 232 ADS simulation ApJ MHD simulation 1982 2001 MHD simulation ApJ 39 ( 71 ) MHD simulation M2
.....Z...^.[.......\..
15 10 16 42 55 55 56 60 62 199310 1995 134 10 8 15 1 13 1311 a s d f 141412 2 g h j 376104 3 104102 232 4 5 51 30 53 27 36 6 Y 7 8 9 10 8686 86 11 1310 15 12 Z 13 14 15 16 102193 23 1712 60 27 17 18 Z
応力とひずみ.ppt
in [email protected] 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
0 (1 ) 0 (1 ) 01 Excel Excel ( ) = Excel Excel =5+ 5 + 7 =5-5 3 =5* 5 10 =5/ 5 5 =5^ 5 5 ( ), 0, Excel, Excel 13E+05 13 10 5 13000 13E-05 13 10 5 0000
1 ( S/E) 006 7 30 0 (1 ) 01 Excel 0 7 3 1 (-4 ) 5 11 5 1 6 13 7 (5-7 ) 9 1 1 9 11 3 Simplex 1 4 (shadow price) 14 5 (reduced cost) 14 3 (8-10 ) 17 31 17 3 18 33 19 34 35 36 Excel 3 4 (11-13 ) 5 41 5 4
(5) 75 (a) (b) ( 1 ) v ( 1 ) E E 1 v (a) ( 1 ) x E E (b) (a) (b)
(5) 74 Re, bondar laer (Prandtl) Re z ω z = x (5) 75 (a) (b) ( 1 ) v ( 1 ) E E 1 v (a) ( 1 ) x E E (b) (a) (b) (5) 76 l V x ) 1/ 1 ( 1 1 1 δ δ = x Re x p V x t V l l (1-1) 1/ 1 δ δ δ δ = x Re p V x t V
TM
NALTR-1390 TR-1390 ISSN 0452-2982 UDC 533.6.013.1 533.6.013.4 533.6.69.048 NAL TECHNICAL REPORT OF NATIONAL AEROSPACE LABORATORY TR-1390 e N 1999 11 NATIONAL AEROSPACE LABORATORY ... 1 e N... 2 Orr-Sommerfeld...
1 NASA 2 PLANET-C 2 ISAS No.320 2007.11
ISSN 0285-2861 2007.11 No. 320 Rocket Girls20079276 ISAS No.320 2007.11 1 1 NASA 2 PLANET-C 2 ISAS No.320 2007.11 ISAS No.320 2007.11 3 3 PLANET-C 5 ISAS BOGUSLAWSKY 72.943.2 4 ISAS No.320 2007.11 ISAS
ac b 0 r = r a 0 b 0 y 0 cy 0 ac b 0 f(, y) = a + by + cy ac b = 0 1 ac b = 0 z = f(, y) f(, y) 1 a, b, c 0 a 0 f(, y) = a ( ( + b ) ) a y ac b + a y
01 4 17 1.. y f(, y) = a + by + cy + p + qy + r a, b, c 0 y b b 1 z = f(, y) z = a + by + cy z = p + qy + r (, y) z = p + qy + r 1 y = + + 1 y = y = + 1 6 + + 1 ( = + 1 ) + 7 4 16 y y y + = O O O y = y
( ) x y f(x, y) = ax
013 4 16 5 54 (03-5465-7040) [email protected] hp://lecure.ecc.u-okyo.ac.jp/~nkiyono/inde.hml 1.. y f(, y) = a + by + cy + p + qy + r a, b, c 0 y b b 1 z = f(, y) z = a + by + cy z = p + qy
: u i = (2) x i Smagorinsky τ ij τ [3] ij u i u j u i u j = 2ν SGS S ij, (3) ν SGS = (C s ) 2 S (4) x i a u i ρ p P T u ν τ ij S c ν SGS S csgs
![: u i = (2) x i Smagorinsky τ ij τ [3] ij u i u j u i u j = 2ν SGS S ij, (3) ν SGS = (C s ) 2 S (4) x i a u i ρ p P T u ν τ ij S c ν SGS S csgs](/thumbs/91/106610562.jpg ": u i = (2) x i Smagorinsky τ ij τ [3] ij u i u j u i u j = 2ν SGS S ij, (3) ν SGS = (C s ) 2 S (4) x i a u i ρ p P T u ν τ ij S c ν SGS S csgs")
15 C11-4 Numerical analysis of flame propagation in a combustor of an aircraft gas turbine, 4-6-1 E-mail: [email protected], 2-11-16 E-mail: [email protected], 4-6-1 E-mail: [email protected],
2 ISAS2006.3
ISAS2006.3 1 ISSN 0285-2861 2006.3 2 ISAS2006.3 ISAS2006.3 3 4 ISAS2006.3 ISAS2006.3 5 6 ISAS2006.3 ISAS2006.3 7 2 8 ISAS2006.3 ISAS2006.3 9 ISAS 2006.3 ISSN 0285-2861 229-8510 3-1-1 TEL042-759-8008 [email protected]
t Z
t Z 1 11 1 1 1 1 1 1 1 1 1 1 1 1 UU 1 as a s 1 1 1 1 1 1 G q w e r q q w e r a s 1 1 1 1 1 1 1 a s a gg g q w e r 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
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
1 a b cc b * 1 Helioseismology * * r/r r/r a 1.3 FTD 9 11 Ω B ϕ α B p FTD 2 b Ω * 1 r, θ, ϕ ϕ * 2 *
448 8542 1 e-mail: [email protected] 1. 400 400 1.1 10 1 1 5 1 11 2 3 4 656 2015 10 1 a b cc b 22 5 1.2 * 1 Helioseismology * 2 6 8 * 3 1 0.7 r/r 1.0 2 r/r 0.7 3 4 2a 1.3 FTD 9 11 Ω B ϕ α B
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
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
870727_ガイドブック2016_vol1.indd
VOL.1 VOL.2 VOL.3 2016 I N D E X C h e c k! 3 C h e c k! 4 5 1 2 3 6 4 C h e c k! 5 7 8 9 10 11 12 Q Q Q A A A Q A C h e c k! 13 14 1,531,513 72,364 3,737,465 1,110,000 3,337,105 1,119,421 C h e c k! 774,000
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
北アルプス_燕岳～穂高_－2.doc
1 ( )22 2001 2226 30 1300m 27kg 2 JR 1000 22 10 3 500-400 au 300 N 2763m 4 N N 10 N 500m PM9:00 5 ( 2 ) 7 23 2001 2226 10 N N 40 10 30 6 5m N N N 12 3 49 30 20 2922m 7 K ( ) U S N ( ) 50 15 N 35 (2758m)
1 n =3, 2 n 3 x n + y n = z n x, y, z 3 a, b b = aq q a b a b b a b a a b a, b a 0 b 0 a, b 2
n =3, 200 2 10 1 1 n =3, 2 n 3 x n + y n = z n x, y, z 3 a, b b = aq q a b a b b a b a a b a, b a 0 b 0 a, b 2 a, b (a, b) =1a b 1 x 2 + y 2 = z 2, (x, y) =1, x 0 (mod 2) (1.1) x =2ab, y = a 2 b 2, z =
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
B line of mgnetic induction AB MN ds df (7.1) (7.3) (8.1) df = µ 0 ds, df = ds B = B ds 2π A B P P O s s Q PQ R QP AB θ 0 <θ<π
8 Biot-Svt Ampèe Biot-Svt 8.1 Biot-Svt 8.1.1 Ampèe B B B = µ 0 2π. (8.1) B N df B ds A M 8.1: Ampèe 107 108 8 0 B line of mgnetic induction 8.1 8.1 AB MN ds df (7.1) (7.3) (8.1) df = µ 0 ds, df = ds B
.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,
50 2 I SI MKSA r q r q F F = 1 qq 4πε 0 r r 2 r r r r (2.2 ε 0 = 1 c 2 µ 0 c = m/s q 2.1 r q' F r = 0 µ 0 = 4π 10 7 N/A 2 k = 1/(4πε 0 qq
49 2 I II 2.1 3 e e = 1.602 10 19 A s (2.1 50 2 I SI MKSA 2.1.1 r q r q F F = 1 qq 4πε 0 r r 2 r r r r (2.2 ε 0 = 1 c 2 µ 0 c = 3 10 8 m/s q 2.1 r q' F r = 0 µ 0 = 4π 10 7 N/A 2 k = 1/(4πε 0 qq F = k r
宇宙天気予報と磁気流体力学
2008 年 7 月 17 日 ( 木 ), 宇宙プラズマ物理学, 東京大学 宇宙天気予報と磁気流体力学 片岡龍峰 理化学研究所 目次 宇宙嵐 フレア放射線 オーロラ嵐 バンアレン帯 磁気流体力学 Frozen-in, 保存形式, RH 関係式 有限体積法による数値解法 数値宇宙天気予報 太陽風 バンアレン帯のモデリング研究 確率予報 数値予報 宇宙天気図 1. 宇宙嵐 代表的な宇宙嵐について解説する
4 B 67 654 (00-) Eec o Pulsaing Gas Injecion on Granular Flow in Gas-Solid Fluidized Beds Asushi MIYOSHI, Toshihiro KAWAGUCHI, Toshisugu TANAKA and Yuaka TSUJI Osaka Universiy, De. o Machanohysics Engineering
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)...