土木学会論文集 B2( 海岸工学 ) Vol. 70, No. 2, 2014, I_016-I_020 非線形長波モデルと流体粒子法による津波シミュレータの開発 Development of a Tsunami Simulator Integrating the Smoothed-Particle Hydrodynamics Method and the Nonlinear Shallow Water Wave Model Tamon SUWA, Fumihiko IMAMURA, and Daisuke SUGAWARA We develop a tsunami simulator integrating a 3-D fluid simulation technology that runs on large-scale parallel computers using smoothed-particle hydrodynamics (SPH) method, together with a 2-D tsunami propagation simulation technique using a nonlinear shallow water wave model. We use the 2-D simulation to calculate tsunami propagation of scale of about 1000km from epicenter to near shore. The 3-D SPH method can be used to calculate the force that a tsunami can exert on a building, and to simulate flooding patterns in urban area of at most km scale. By applying the processing power of computers to the technologies resulting from this research, we seek to contribute to improved disaster preparedness and disaster mitigation through a better understanding of tsunami's mechanisms. 1. 2011 3 11 TUNAMI-N2 Smoothed Particle Hydrodynamics SPH 2006 Boussinesq Moving Particle Semi-implicit MPS Boussinesq 2. 1 2 3 h M N x y D n Manning TUNAMI-N2 TUNAMI-N2 Goto 1997 SPH Suwa 2013.
非線形長波モデルと流体粒子法による津波シミュレータの開発 I_17 4 6 5 m ρ v p h g a b a 2h b r ab a b 4 5 5 Fang 2009 2011 5 W r ab h 5 Wendland 1995 q= r ab /h a d W r ab h 1 3 21/ 16π h 3 p r 7 r 0 1000kgm -3 p 0 2.25 10 7 Pa p 0 SPH SPH 2 SPH SPH 2 図 -1a - d 1350m 450m 150m 50m e SPH d 図 -2 SPHx y SPH
I_18 土木学会論文集 B2( 海岸工学 ),Vol. 70,No. 2,2014 dx 2 SPH 2 SPH SPH 0 r 0 SPH V SPH 2 v S n dt V 8 v 2 v=0 SPH V 0 SPH V 0 -V 0 V 0 -V 0 +V 0 SPH 2 v r 0 3. 1 1350m 450m 150m 50m 図 -1 図 -1 d 0.5km 1.7km SPH 図 -1 e 図 -3 SPH 図 -1 e SPH 図 -4 SPH
非線形長波モデルと流体粒子法による津波シミュレータの開発 I_19 SPH 3,360 1m SPH 図 -2 図 -4 図 -2 SPH 0sec 20 図 -3 SPH 図 -1 e SPH SPH 図 -4 SPH 図 -1 e SPH 230 40 336 Fujitsu PRIMERGY RX200S3 8 64CPU 20 図 -5 SPH 2 SPH 図 -5 1/10 15m 500m 図 -5 2 1 2 SPH 9 x 0 =500m λ=75m A 3.0m 15.0m 2 図 -6 x=100m 0m A=15.0m A=3.0m A=15m SPH A=3m SPH 1/10 SPH A=3m x=150m SPH x=150m SPH 図 -6 SPHA=15m A=3.0m m
I_20 土木学会論文集 B2( 海岸工学 ),Vol. 70,No. 2,2014 図 -6 SPH x=100m SPH 4. 3 SPH 図 -4 8 160 0.5km 1.7km 1m 3 SPH 160 4 2 0.1m 0.5m 10km 100 参考文献 2006 Boussinesq 53 pp. 11-15. 2011 B3 Vol. 67, pp. 268-273. Fang, J., A. Parriaux, M. Rentschler, and C. Ancey (2009) : Improved SPH methods for simulating free surface flows of viscous fluids, Applied numerical mathematics, Vol. 59, pp. 251-271. Goto, C., Y. Ogawa, N. Shuto, and F. Imamura (1997) : Numerical method of tsunami simulation with the leap-frog scheme (IUGG/IOC Time Project), IOC Manual, UNESCO, No. 35, 126 p. Suwa, T., T. Nakagawa, and K. Murakami (2013) : A Study of the Wave Transformation Passing over an Artificial Reef using SPH Method, Journal of computational science and technology, Vol. 7, No. 2, pp. 126-133. Wendland, H. (1995): Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree, Advances in computational mathematics, Vol. 4, pp. 389-396.