An efficient Eulerian finite element method for the shallow water equations

The accuracy and efficiency of an Eulerian method is assessed by solving the non-linear shallow water equations and compared with the performances of an existing semi-Lagrangian method. Both methods use a linear non-conforming finite element discretization for velocity and a linear conforming finite element discretization for surface elevation. This finite element pair is known to be computationally efficient and free of pressure modes. The model equations are carefully derived and a comparison is performed by simulating the propagation of slow Rossby waves in the Gulf of Mexico. Simulations show that the Eulerian model performs well and gives results comparable to high order semi-Lagrangian schemes using kriging interpolators. C 2004 Elsevier Ltd. All rights reserved.

[1]  Andrew Staniforth,et al.  An accurate interpolating scheme for semi-Lagrangian advection on an unstructured mesh for ocean modelling , 1997 .

[2]  A. Staniforth,et al.  Semi-Lagrangian integration schemes for atmospheric models - A review , 1991 .

[3]  D. Krige.  A statistical approach to some basic mine valuation problems on the Witwatersrand, by D.G. Krige, published in the Journal, December 1951 : introduction by the author , 1951 .

[4]  Clive Temperton,et al.  An Efficient Two‐Time‐Level Semi‐Lagrangian Semi‐Implicit Integration Scheme , 1987 .

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

[6]  G. Matheron.  The intrinsic random functions and their applications , 1973, Advances in Applied Probability.

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

[8]  Jens Schröter,et al.  A diagnostic stabilized finite-element ocean circulation model , 2003 .

[9]  Janusz A. Pudykiewicz,et al.  Some properties and comparative performance of the semi‐Lagrangian method of Robert in the solution of the advection‐diffusion equation , 1984 .

[10]  P. Raviart,et al.  On a Finite Element Method for Solving the Neutron Transport Equation , 1974 .

[11]  G. Christakos.  On the Problem of Permissible Covariance and Variogram Models , 1984 .

[12]  A. Robert.  A Semi-Lagrangian and Semi-Implicit Numerical Integration Scheme for the Primitive Meteorological Equations , 1982 .

[13]  Andrew Staniforth,et al.  A Semi-implicit Semi-Lagrangian Finite-Element Shallow-Water Ocean Model , 2000 .

[14]  A. D. Kirwan,et al.  Genesis of a Gulf of Mexico ring as determined from kinematic analyses , 1987 .

[15]  André Robert,et al.  A stable numerical integration scheme for the primitive meteorological equations , 1981 .

[16]  J. R. Bates,et al.  Improving the Estimate of the Departure Point Position in a Two-Time Level Semi-Lagrangian and Semi-Implicit Scheme , 1987 .

[17]  Francois Thomasset,et al.  A noise-free finite element scheme for the two-layer shallow water equations , 1984 .

[18]  D. Purnell.  Solution of the Advective Equation by Upstream Interpolation with a Cubic Spline , 1976 .

[19]  A. Robert,et al.  A semi-Lagrangian and semi-implicit numerical integration scheme for multilevel atmospheric models , 1985 .

[20]  Roy A. Walters,et al.  An unstructured grid, three‐dimensional model based on the shallow water equations , 2000 .

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

[22]  F. Trochu.  A contouring program based on dual kriging interpolation , 1993, Engineering with Computers.

[23]  Jens Schröter,et al.  A finite-element ocean model: principles and evaluation , 2004 .

[24]  Eric Deleersnijder,et al.  Delaunay mesh generation for an unstructured-grid ocean general circulation model , 2000 .

[25]  Eric Blayo,et al.  A Comparison of Two Numerical Methods for Integrating a Quasi-geostrophic Multilayer Model of Ocean Circulations , 1994 .

[26]  A Quantitative Analysis of the Dissipation Inherent in Semi-Lagrangian Advection , 1988 .

[27]  E. Süli,et al.  Discontinuous hp-finite element methods for advection-diffusion problems , 2000 .

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

[29]  P. Raviart,et al.  Conforming and nonconforming finite element methods for solving the stationary Stokes equations I , 1973 .

[30]  A. Snider,et al.  Atmospheric and ocean modeling with an adaptive finite element solver for the shallow-water equations , 1986 .

[31]  Changsheng Chen,et al.  An Unstructured Grid, Finite-Volume, Three-Dimensional, Primitive Equations Ocean Model: Application to Coastal Ocean and Estuaries , 2003 .

[32]  Eric Deleersnijder,et al.  A comparison of three finite elements to solve the linear shallow water equations , 2003 .

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

[34]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[35]  Andrew J. Weaver,et al.  A Diagnostic Barotropic Finite-Element Ocean Circulation Model , 1995 .

[36]  Eric Deleersnijder,et al.  Advection schemes for unstructured grid ocean modelling , 2004 .