This program implements an HLLC Riemann solver for the full shallow water equations, coded from scratch based on the methodology described by Hou et al. (2015). The case study in the paper simulates the 1993 tsunami that struck Okushiri Island, off the west coast of Hokkaido, Japan. The focus is Monai Valley on the island’s southwestern coast, where an extreme vertical runup of 30 metres was recorded. The current simulation predicts a maximum runup of approximately 32 metres (to scale), closely matching field observations. For more details, visit the NOAA Centre for Tsunami Research.
animation.mp4
| File | Purpose |
|---|---|
| main.f90 | Main program containing functions and subroutines for flow solver. |
| monai.py* | Pre-processing program to set up initial conditions for Monai Valley case study. |
| read.py | Post-processing program to create flow animation with Plotly. |
| b_raw.txt | Raw bathymetry data for Monai Valley case study used in monai.py. |
| ht_raw.txt | Raw inflow height data for Monai Valley case study used in monai.py. |
| monai.geo | File containing geometric information to generate mesh in Gmsh. |
| monai.su2 | File containing mesh information to input to main program. |
* Run monai.py if running the case study for the first time, or when there is a change in mesh size or time step.
[1] Hou et al. (2015). An Efficient Unstructured MUSCL Scheme for Solving the 2D Shallow Water Equations.
[2] Song et al. (2011). A Robust Well-Balanced Finite Volume Model for Shallow Water Flows with Wetting and Drying over Irregular Terrain.