410 likes | 424 Views
This article discusses the simulation of tsunamis generated by earthquakes in Japan's Nankai Trough and the flow of groundwater through porous media. It showcases the capabilities of the Earth Simulator in simulating these phenomena.
E N D
Tsunami Simulation • Supposed to be generated by Earthquake at Nankai Trough in the South-West part of Japan. • Quasi 3D Shallow-Water Equations • Kawahara-Kashiyama Scheme • 3-Step Explicit Time Marching • Linear Triangle Elements
Nankai Trough Cabinet Office/ Government of Japan
Nankai Trough Nagoya University
Tsunami Simulation8.80105nodes, 1.76106elem's Mesh Size = 1.25 km
Demonstration onthe Earth Simulator • 64 PEs on Earth Simulator • 512 GFLOPS peak • 1,200 km 1,200 km Region • 923,521 nodes, 1,843,200 elements • 42,255 steps • Parallel Programming Models • Flat MPI • OpenMP/MPI Hybrid
Demonstration onthe Earth Simulator64 PEs, 8 SMP Nodes, 42,255 steps ●Flat MPI ●Hybrid Parallel Speed-Up
Demonstration onthe Earth Simulator64 PEs, 8 SMP Nodes, 42,255 steps ●Flat MPI 222 GFLOPS (43% of the Peak) GFLOPS
Background • Deep Geological Disposal of HLW • High-Level Radioactive Waste • Natural/Engineered Barriers • Long-term simulation (106~108 y) for nuclide transportation • Random Heterogeneous Field • Sedimentary Rock : Porous Media • Granite : Fracture
V.Test Hole V.Test Hole Sedimentary Rock ANDRA Tunnels Granite Fault Deep Geological Disposalof HLWDisposal Site JNC, Japan.
Demonstration on the Earth Simulator • Coupled Flow Equations • Saturated Darcy Flow & Energy Equations • Primitive Variable Approach • ICCG/Multigrid for Poisson’s Eqn. • Up to 2M Meshes on ES (1PE) • 2.4 GFLOPS (30% of the Peak)
Problem Setting Uniform Temperature/Pressure at Z=0 surface Uniform Initial Flow Heat Generation
Results on EST=0.50, Temperature Uini=0.00 Uini=0.02
Results on EST=1.00, Temperature Uini=0.00 Uini=0.02
Results on EST=1.50, Temperature Uini=0.00 Uini=0.02
Results on EST=2.00, Temperature Uini=0.00 Uini=0.02
Results on EST=2.50, Temperature Uini=0.00 Uini=0.02
Results on EST=3.00, Temperature Uini=0.00 Uini=0.02
Results on EST=3.50, Temperature Uini=0.00 Uini=0.02
Results on EST=4.00, Temperature Uini=0.00 Uini=0.02
Results on EST=4.50, Temperature Uini=0.00 Uini=0.02
Results on EST=5.00, Temperature Uini=0.00 Uini=0.02
Results on EST=5.50, Temperature Uini=0.00 Uini=0.02
Results on EST=6.00, Temperature Uini=0.00 Uini=0.02
Results on EST=6.50, Temperature Uini=0.00 Uini=0.02
Results on EST=7.00, Temperature Uini=0.00 Uini=0.02
Results on EST=7.50, Temperature Uini=0.00 Uini=0.02
Results on EST=8.00, Temperature Uini=0.00 Uini=0.02
Results on EST=8.50, Temperature Uini=0.00 Uini=0.02
Results on EST=9.00, Temperature Uini=0.00 Uini=0.02
Results on EST=9.50, Temperature Uini=0.00 Uini=0.02
Results on EST=10.0, Temperature Uini=0.00 Uini=0.02
Performance on 1 PEGFLOPS (Over All) ●ICCG ●MG
Thermal Stress Analysis of Nuclear Experimental Facility ( Tube Sheet )
Demonstration on the Earth Simulator • FEA-Visualization Concurrent Coupling • 200M DOF on ES (64 SMP nodes) • 2.8 GFLOPS / PE (35% of the Peak)
Thermal Stress Analysis of Tube Sheet • 全体を均等加熱 • PVRによるミーゼス応力表示