71 特集 オープンソースの大きな流れ Nonlinear Sloshing Analysis in a Three-dimensional Rectangular Pool Ken UZAWA, The Center for Computational Sciences and E-systems, Japan Atomic Energy Agency 1 1.1 ( ) (RIST) (ORNL/RSICC) 8 8 (DOE) MCNP( ) RELAP( ) MCNP RELAP TRACE &(erification & alidation ) erification alidation 77 8587 5-1-5 E-mail: uzawa.ken@jaea.go.jp erification 1) 1. ) 198 R.Stallman 199 Linux Apache, MySQL GNU OS 3 7) & &
7 OpenFOAM 6) OpenFOAM (Fujitsu PRIMERGY BX9, TFLOPS) Fluent 8) ( ) 9, 1) 11 13) OpenFOAM - realizable k-ε 1) Launder-Gibson 15) OpenFOAM 1.6 CFD ( ) 16 18) 3 3.1 1 (L, H, W )[m] (x,d )[m] 1[kg/m 3 ],.1[Pa s] 1[kg/m 3 ], 1 [Pa s] g 9.81[m/s ] (N x,n y,n z ) d L [m] W [m] y H [m] [m] x z x [m] 1 3. Martin 16)
73 (x,d ) = (.5715,.5715) [m] (x,d ) = (.8575,.5715) [m] (L,H,W)=(d, 15x,d ) (N x,n y,n z ) = (1, 75, 5) Martin Martin ( ) RKE LG LG RKE t =.35[s] t =.75[s] LG RKE LG LG [m] [m].7 (Martin).6.5 RKE LG..3..1.1..3..5...18.16.1.1.6 (Martin) RKE LG..5.1.15 [Pa] 15 1 5 experiment (Hu) laminar realizable RKE k-ε LG LG -5...6.8 1. 3 A t =.5[s] LG LG 3.3 Hu 19) Hu (x,d,l,h,w) = (.68,.1, 1.18,.,.1) [m] (N x,n y,n z )=(7, 96, 8) 3 A (x, y, z) = (1.18,.1,.6) [m] 8 ( ) RKE LG t =[s] t =.5[s] 3. - K = P + G +, (1)
7 K = P + G +, () K = P + G +. (3) P G K K d 1 = dt ρu i d, () P p = u i d, x i (5) G = u i ρgδ i d, (6) = u i µd ij d. x j (7) K,K d 1 = dt ρū i d, (8) P,P = p ū i d, x i (9) G,G = ū i ρgδ i d, (1) = = ū i ū i (µ t + µ)d ij d, (11) x j ( ρτ ij +µd ij )d. (1) x j ρ u i p g δ ij µ D ij = 1 ( ui ) x j + u j x i τ ij = ũ i ũ j 5 RKE (c)lg 5 RKE LG t=.35[s] Smagorinsky LES LG LES (N x,n y,n z )=(9, 19, 96) 5(d) LG LES LG RKE 5 K G P...6 K P G...6 G K P -...6 5 3 1 K (LES) K (LES) (d) (c).1..3..5 ( )
75.1 6 (L, H, W )= (3.,.,.5)[m] (d,h)= (., 1.5)[m] d =1.5[m] 3 K-Net N-S ( 7) (N x, N y, N z ) = (16, 1, 5) 1). 8 RKE LG LES 8 t =5[s] RKE LG LES LG LES 1-1 [cm] - -3 L [m] W [m] - 1 3 5 [s] 5 6 [cm/s ] z 1 - y x x d [m] [m] h [m] H [m] 8 6 - -6-8 -1 5 1 15 5 [s] 8 t =5[s].3 9 RKE LG LES RKE LG LES 1 7 (3 K-Net N-S )
76 [Pa s]. 9 5 3 1 5 1 3 5 5 6 38 3 3 5 5 t =5[s] OpenFOAM 6 ( ) CAE 11 CAE ( ) K-NET K-NET 1) eye (, 1) 9-55. ) http://www.opensource.org/docs/osd 3) http://www.code-aster.org/ ) http://adventure.sys.t.u-tokyo.ac.jp/jp/ 5) http://www.csc.fi/english/pages/elmer 6) http://www.openfoam.com/ 7) http://www.ciss.iis.u-tokyo.ac.jp/rss1/theme/multi/fluid/ 8) Fire Service Act, Firefighting Hazards, No.1. 9) Muto, K., Kasai, Y., Nakahara, M. & Ishida, Y.: Experimental Tests on Sloshing Response of a Water Pool with Submergeb Blocks, Proc. Pressure essels and Piping Conference (1985) 9-1. 1) :,, 33B (1987) 79-9. 11) Huerta, A., Liu, W.K. & Gvildys, J.: Large amplitude sloshing with submerged blocks, Proc. Pressure essels and Piping Conference (1989) 13-18. 1) Sakai, M., Higashi, S., Sato, K. & Tanaka. N.: Nonlinear sloshing response evaluation method for rectangular tank with overflowing water, CRIEPI Report N631 (7). 13) Eguchi Y.: Development of SMART-slg to simulate sloshing in an oil tank with a floating roof, CRIEPI Report N75 (8). 1) Shih, T. -H., Liou, W. W., Shabbir, A., Yang, Z. & Zhu, J.: A New k ϵ Eddy iscosity Model for High Reynolds Number Turbulent Flows, Computers Fluids, (1995) 7-38. 15) Gibson, G. G. & Launder, B. E.: Ground effects on pressure fluctuations in the atmospheric boundary layer, J. Fluid. Mech., 85 (1978) 91-511. 16) Martin, J.C. & Moyce, W.J.: PartI. An experimental study of the collapse of liquid columns on a rigid horizontal plane, Phil. Trans. R. Soc. Lond. A, (195) 31-3. 17) Koshizuka, S., Tamako, H. & Oka, Y.: A particle method for incompressible viscous flow with fluid fragmentation, Computational Fluid Mechanics Journal, 113 (1995) 13-17. 18) Chanson, H.: Analytical Solutions of Laminar and Turbulent Dam Break Wave, Proc. Intl Conf. Fluvial Hydraulics River Flow 6, A3 (6) 65-7. 19) Hu, C.: A CIP-based method for numerical simulations of violent free surface flows, J. Mar. Sci. Technol, 9 () 13-157.