Advection schemes for unstructured grid ocean modelling

We study advection schemes for unstructured grid ocean models. Four linear advection schemes are investigated by solving a scalar transport equation. Schemes under consideration include continuous, nonconforming and discontinuous finite elements and finite volumes. A comprehensive derivation of the numerical schemes is presented and conservation and dispersion properties are discussed. An assessment is made by performing the test problem introduced by Hecht et al. [J. Geophys. Res. 100 (1995) 20763] in which a passive scalar field is advected through an analytical Stommel gyre. It is found that continuous finite elements and finite volumes have some difficulties to represent accurately solutions with steep gradients. As a result they are prone to generate unphysical oscillations. On the other hand, discontinuous and non-conforming finite element schemes perform better. This is due to their higher flexibility that makes them better suited to highly sheared flows. (C) 2003 Elsevier Ltd. All rights reserved.

[1]  Thomas J. R. Hughes,et al.  The Continuous Galerkin Method Is Locally Conservative , 2000 .

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

[3]  R. C. Berger,et al.  Discrete Fluxes and Mass Balance in Finite Elements , 2002 .

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

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

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

[7]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

[8]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

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

[10]  海老 豊,et al.  F. Thomasset: Implementation of Finite Element Methods for Navier-Stokes Equations, Springer-Verlag, New York and Heidelberg, 1981, viii+162ページ, 9,270円 (Springer Series in Computational Physics). , 1983 .

[11]  Beth A. Wingate,et al.  A better, more discriminating test problem for ocean tracer transport , 2000 .

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

[13]  Thomas J. R. Hughes,et al.  Consistent finite element calculations of boundary and internal fluxes , 1998 .

[14]  J. Rae,et al.  Implementation of finite element Methods for Navier-Stokes equations , 1982 .

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

[16]  P. Smolarkiewicz.  A Fully Multidimensional Positive Definite Advection Transport Algorithm with Small Implicit Diffusion , 1984 .

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

[18]  Akira Okubo,et al.  Oceanic diffusion diagrams , 1971 .

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

[20]  Philip J. Rasch,et al.  Upwind‐weighted advection schemes for ocean tracer transport: An evaluation in a passive tracer context , 1995 .

[21]  Robert L. Lee,et al.  THE CONSISTENT GALERKIN FEM FOR COMPUTING DERIVED BOUNDARY QUANTITIES IN THERMAL AND/OR FLUIDS PROBLEMS , 1987 .

[22]  Jiayan Yang.  On the importance of resolving the western boundary layer in wind-driven ocean general circulation models , 2003 .

[23]  J. W. Eastwood,et al.  Springer series in computational physics Editors: H. Cabannes, M. Holt, H.B. Keller, J. Killeen and S.A. Orszag , 1984 .

[24]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[25]  Ivo Babuška,et al.  The post-processing approach in the finite element method—part 1: Calculation of displacements, stresses and other higher derivatives of the displacements , 1984 .

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

[27]  D. B. Simons,et al.  Forms of bed roughness in alluvial channels , 1961 .

[28]  T. Hughes,et al.  The Galerkin/least-squares method for advective-diffusive equations , 1988 .