2D–3D hybrid stabilized finite element method for tsunami runup simulations

This paper presents a two-dimensional (2D)–three-dimensional (3D) hybrid stabilized finite element method that enables us to predict a propagation process of tsunami generated in a hypocentral region, which ranges from offshore propagation to runup to urban areas, with high accuracy and relatively low computational costs. To be more specific, the 2D shallow water equation is employed to simulate the propagation of offshore waves, while the 3D Navier–Stokes equation is employed for the runup in urban areas. The stabilized finite element method is utilized for numerical simulations for both of the 2D and 3D domains that are independently discretized with unstructured meshes. The multi-point constraint and transmission methods are applied to satisfy the continuity of flow velocities and pressures at the interface between the resulting 2D and 3D meshes, since neither their spatial dimensions nor node arrangements are consistent. Numerical examples are presented to demonstrate the performance of the proposed hybrid method to simulate tsunami behavior, including offshore propagation and runup to urban areas, with substantially lower computation costs in comparison with full 3D computations.

[1]  Koji Fujima,et al.  Development of the 2D/3D Hybrid Model for Tsunami Numerical Simulation , 2002 .

[2]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[3]  Nicolas G. Wright,et al.  A robust 2D shallow water model for solving flow over complex topography using homogenous flux method , 2013 .

[4]  Koji Fujima,et al.  NUMERICAL ANALYSIS OF TSUNAMI BY USING 2D/3D HYBRID MODEL , 2001 .

[5]  Tayfun E. Tezduyar,et al.  Stabilized-® nite-element / interface-capturing technique for parallel computation of unsteady ̄ ows with interfaces q , 2000 .

[6]  Graham F. Carey,et al.  A symmetric formulation and SUPG scheme for the shallow-water equations , 1996 .

[7]  Tayfun E. Tezduyar,et al.  Parallel finite element methods for large-scale computation of storm surges and tidal flows , 1997 .

[8]  Tayfun E. Tezduyar,et al.  Space–time SUPG finite element computation of shallow-water flows with moving shorelines , 2011 .

[9]  T. Umetsu,et al.  Stabilized Bubble Function Method for Shallow Water Long Wave Equation , 2003 .

[10]  Tayfan E. Tezduyar,et al.  Stabilized Finite Element Formulations for Incompressible Flow Computations , 1991 .

[11]  Ethan J. Kubatko,et al.  A wetting and drying treatment for the Runge-Kutta discontinuous Galerkin solution to the shallow water equations , 2009 .

[12]  Kwok Fai Cheung,et al.  Depth‐integrated, non‐hydrostatic model with grid nesting for tsunami generation, propagation, and run‐up , 2010 .

[13]  Gui-Rong Liu,et al.  The Finite Element Method , 2007 .

[14]  Mourad Heniche,et al.  A two-dimensional finite element drying-wetting shallow water model for rivers and estuaries , 2000 .

[15]  Jalal Abedi,et al.  Simulation technique for wave generation , 2003 .

[16]  The Application of a Hybrid 2D/3D Numerical Tsunami Inundation-Propagation Flow Model to the 2011 off the Pacific Coast of Tohoku Earthquake Tsunami , 2013 .