A high‐order triangular discontinuous Galerkin oceanic shallow water model

A high-order triangular discontinuous Galerkin (DG) method is applied to the two-dimensional oceanic shallow water equations. The DG method can be characterized as the fusion of finite elements with finite volumes. This DG formulation uses high-order Lagrange polynomials on the triangle using nodal sets up to 15th order. Both the area and boundary integrals are evaluated using order 2N Gauss cubature rules. The use of exact integration for the area integrals leads naturally to a full mass matrix; however, by using straight-edged triangles we eliminate the mass matrix completely from the discrete equations. Besides obviating the need for a mass matrix, triangular elements offer other obvious advantages in the construction of oceanic shallow water models, specifically the ability to use unstructured grids in order to better represent the continental coastlines for use in tsunami modeling. In this paper, we focus primarily on testing the discrete spatial operators by using six test cases—three of which have analytic solutions. The three tests having analytic solutions show that the high-order triangular DG method exhibits exponential convergence. Furthermore, comparisons with a spectral element model show that the DG model is superior for all polynomial orders and test cases considered. Copyright © 2007 John Wiley & Sons, Ltd.

[1]  Francis X. Giraldo.  A SPECTRAL ELEMENT SEMI-LAGRANGIAN METHOD FOR THE SHALLOW WATER EQUATIONS ON UNSTRUCTURED GRIDS , 1998 .

[2]  J. Remacle,et al.  Shock detection and limiting with discontinuous Galerkin methods for hyperbolic conservation laws , 2004 .

[3]  Jean-François Remacle,et al.  An adaptive discretization of shallow‐water equations based on discontinuous Galerkin methods , 2006 .

[4]  T. Koornwinder.  Two-Variable Analogues of the Classical Orthogonal Polynomials , 1975 .

[5]  E. Toro.  Shock-Capturing Methods for Free-Surface Shallow Flows , 2001 .

[6]  Jan S. Hesthaven,et al.  From Electrostatics to Almost Optimal Nodal Sets for Polynomial Interpolation in a Simplex , 1998 .

[7]  M. Piggott,et al.  A Nonhydrostatic Finite-Element Model for Three-Dimensional Stratified Oceanic Flows. Part II: Model Validation , 2004 .

[8]  Ethan J. Kubatko,et al.  hp Discontinuous Galerkin methods for advection dominated problems in shallow water flow , 2006 .

[9]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[10]  K. Anastasiou,et al.  SOLUTION OF THE 2D SHALLOW WATER EQUATIONS USING THE FINITE VOLUME METHOD ON UNSTRUCTURED TRIANGULAR MESHES , 1997 .

[11]  Spencer J. Sherwin,et al.  A triangular spectral/hp discontinuous Galerkin method for modelling 2D shallow water equations , 2004 .

[12]  Stephan T. Grilli,et al.  Source Constraints and Model Simulation of the December 26 , 2007 .

[13]  Zygmunt Kowalik,et al.  NUMERICAL MODELING OF THE GLOBAL TSUNAMI : Indonesian Tsunami of 26 December 2004 , 2006 .

[14]  Hong Li,et al.  The discontinuous Galerkin finite element method for the 2D shallow water equations , 2001 .

[15]  Jörn Behrens,et al.  Comparison of approaches with unstructured meshes and nested grids , 2022 .

[16]  Moshe Dubiner.  Spectral methods on triangles and other domains , 1991 .

[17]  H. Stommel,et al.  The westward intensification of wind‐driven ocean currents , 1948 .

[18]  David A. Kopriva,et al.  Metric Identities and the Discontinuous Spectral Element Method on Curvilinear Meshes , 2006, J. Sci. Comput..

[19]  Dirk Schwanenberg,et al.  A Discontinuous Galerkin Method for the Shallow Water Equations with Source Terms , 2000 .

[20]  Francis X. Giraldo,et al.  A diagonal-mass-matrix triangular-spectral-element method based on cubature points , 2007 .

[21]  Of references. , 1966, JAMA.

[22]  Mark A. Taylor,et al.  An Algorithm for Computing Fekete Points in the Triangle , 2000, SIAM J. Numer. Anal..

[23]  B. Dick,et al.  The International Federation of Red Cross and Red Crescent Societies , 1992, Tropical doctor.

[24]  A. K. Lal,et al.  The Great Tsunami of 26 December 2004: A description based on tide-gauge data from the Indian subcontinent and surrounding areas , 2006 .

[25]  R. Cools.  Monomial cubature rules since “Stroud”: a compilation—part 2 , 1999 .

[26]  Natalja Rakowsky,et al.  amatos: Parallel adaptive mesh generator for atmospheric and oceanic simulation , 2005 .

[27]  Dale B. Haidvogel,et al.  A three-dimensional spectral element model for the solution of the hydrostatic primitive equations , 2003 .

[28]  Vasily Titov,et al.  Implementation and testing of the Method of Splitting Tsunami (MOST) model , 1997 .

[29]  Extracting the 2004 Indian Ocean tsunami signals from sea surface height data observed by satellite altimetry , 2008 .

[30]  Ian Parsons,et al.  Surface deformation due to shear and tensile faults in a half-space , 1986 .

[31]  Jonathan Richard Shewchuk,et al.  Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator , 1996, WACG.

[32]  J. Boyd,et al.  A staggered spectral element model with application to the oceanic shallow , 1995 .

[33]  Traian Iliescu,et al.  Submitted to: Journal of Physical Oceanography , 2022 .

[34]  M. Gerritsen,et al.  An unstructured-grid, finite-volume, nonhydrostatic, parallel coastal ocean simulator , 2006 .

[35]  V. Titov,et al.  Sumatra tsunami: lessons from modeling , 2006 .

[36]  Francis X. Giraldo,et al.  High-order triangle-based discontinuous Galerkin methods for hyperbolic equations on a rotating sphere , 2006, J. Comput. Phys..

[37]  John P. Boyd,et al.  Equatorial Solitary Waves. Part I: Rossby Solitons , 1980 .

[38]  F. Dupont,et al.  The Adaptive Spectral Element Method and Comparisons with More Traditional Formulations for Ocean Modeling , 2004 .

[39]  H. Kitamura.  BookThe emerging Japanese superstate: Challenge and response: by Herman Kahn 274 pages. $7.95. Prentice-Hall, Inc, Englewood Cliffs, NJ , 1971 .

[40]  Francis X. Giraldo,et al.  Semi‐implicit time‐integrators for a scalable spectral element atmospheric model , 2005 .

[41]  Hong Ma,et al.  A spectral element basin model for the shallow water equations , 1993 .

[42]  S. Lentz,et al.  Classical tidal harmonic analysis including error estimates in MATLAB using T_TIDE , 2002 .

[43]  A. Stroud.  Approximate calculation of multiple integrals , 1973 .

[44]  Francis X. Giraldo,et al.  The Lagrange-Galerkin Method for the Two-dimensional Shallow Water Equations on Adaptive Grids , 2000 .

[45]  Daniel R. Lynch,et al.  Comprehensive coastal circulation model with application to the Gulf of Maine , 1996 .

[46]  L. Kantha,et al.  Numerical models of oceans and oceanic processes , 2000 .

[47]  John P. Boyd,et al.  Equatorial Solitary Waves. Part 3: Westward-Traveling Modons , 1985 .

[48]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[49]  Jan S. Hesthaven,et al.  Nodal high-order discontinuous Galerkin methods for the spherical shallow water equations , 2002 .

[50]  Ronald Cools,et al.  A survey of numerical cubature over triangles , 1993 .

[51]  Francis X. Giraldo,et al.  A nodal triangle-based spectral element method for the shallow water equations on the sphere , 2005 .

[52]  Richard E. Thomson,et al.  The 26 December 2004 Sumatra Tsunami: Analysis of Tide Gauge Data from the World Ocean Part 1. Indian Ocean and South Africa , 2007 .

[53]  Alain Vande Wouwer,et al.  Upwinding in the method of lines , 2001 .

[54]  R. Cools,et al.  Monomial cubature rules since “Stroud”: a compilation , 1993 .

[55]  Vasily Titov,et al.  The Global Reach of the 26 December 2004 Sumatra Tsunami , 2005, Science.

[56]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[57]  Clinton N Dawson,et al.  A discontinuous Galerkin method for two-dimensional flow and transport in shallow water , 2002 .

[58]  Roy A. Walters,et al.  A semi‐implicit finite element model for non‐hydrostatic (dispersive) surface waves , 2005 .

[59]  Fred F. Pollitz,et al.  Implications of the 26 December 2004 Sumatra-Andaman earthquake on tsunami forecast and assessment models for great subduction-zone earthquakes , 2007 .