A robust unstructured grid discretization for 3-dimensional hydrostatic flows in spherical geometry: A new numerical structure for ocean general circulation modeling

Current models of the general circulation of the global oceans employ a spatial discretization of the relevant hydrodynamic fields on Cartesian rectilinear grids. For many applications, significant benefit would be expected to accrue from the versatility offered by unstructured grids. However, until very recently, available numerical methods for performing integrations on unstructured grids could not conserve discrete dynamical invariants, a numerical model characteristic that is important for robust ocean simulations over large space and time scales (as needed, for instance, in climate modeling applications). Our purpose in this paper is to describe such a conservative discretization scheme for rotating hydrostatic Boussinesq fluid dynamics on general triangular tessellations of the sphere, and to demonstrate its properties in a number of simulations that incorporate realistic ocean basin geometry. Several different implicit time discretizations are possible, each of which exhibits a form of exact numerical energy conservation in the absence of dissipation. The properties of this new numerical methodology are validated through analysis of a sequence of unforced nonlinear dynamical problems, which clearly demonstrate the capacity of the model to resolve the geostrophic adjustment process and the onset of baroclinic instability in collapsing density fronts. As a final test, a number of ocean modeling experiments with realistic climatological and wind-stress forcing are performed in order to investigate the manner in which different mesh structures and resolutions influence the simulated phenomenology. As the theoretical properties of the numerical methodology suggest, it is thereby shown to be both robust and stable. The further work that will be required to implement this structure in a state-of-the-art oceanic general circulation model, as well as other potential applications of the techniques, are discussed in the concluding section of the paper.

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

[2]  Stephen M. Griffies,et al.  Fundamentals of Ocean Climate Models , 2004 .

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

[4]  J. Côté,et al.  A Lagrange multiplier approach for the metric terms of semi‐Lagrangian models on the sphere , 1988 .

[5]  Vincenzo Casulli,et al.  A Semi-Implicit Finite Difference Model for Three-Dimensional Tidal Circulation , 1992 .

[6]  P. Frederickson,et al.  Icosahedral Discretization of the Two-Sphere , 1985 .

[7]  Edward P. Myers,et al.  A cross-scale model for 3D baroclinic circulation in estuary-plume-shelf systems: I , 2004 .

[8]  K. Bryan,et al.  The Circulation of the World Ocean: A Numerical Study. Part I, A Homogeneous Model , 1972 .

[9]  James C. McWilliams,et al.  The Formation of Oceanic Eddies in Symmetric and Asymmetric Jets. Part I: Early Time Evolution and Bulk Eddy Transports , 1995 .

[10]  Deborah Greaves,et al.  Hierarchical tree-based finite element mesh generation , 1999 .

[11]  R. Dmowska,et al.  International Geophysics Series , 1992 .

[12]  D. Schmidt,et al.  Accuracy and conservation properties of a three-dimensional unstructured staggered mesh scheme for fluid dynamics , 2002 .

[13]  Robert M. Chervin,et al.  Ocean general circulation from a global eddy‐resolving model , 1992 .

[14]  J. Dukowicz,et al.  Implicit free‐surface method for the Bryan‐Cox‐Semtner ocean model , 1994 .

[15]  S. Gorshkov,et al.  World ocean atlas , 1976 .

[16]  Luc Giraud,et al.  A Set of GMRES Routines for Real and Complex Arithmetics , 1997 .

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

[18]  R. Pacanowski,et al.  Parameterization of Vertical Mixing in Numerical Models of Tropical Oceans , 1981 .

[19]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

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

[21]  R. Chervin,et al.  A simulation of the global ocean circulation with resolved eddies , 1988 .

[22]  James C. McWilliams,et al.  The formation of oceanic eddies in symmetric and asymmetric jets. Part II: Late time evolution and coherent vortex formation , 1996 .

[23]  G. R. Stuhne,et al.  Vortex Erosion and Amalgamation in a New Model of Large Scale Flow on the Sphere , 1996 .

[24]  W. R. Peltier,et al.  The climate of the Earth at Last Glacial Maximum: statistical equilibrium state and a mode of internal variability , 2004 .

[25]  Ross J. Murray,et al.  Explicit Generation of Orthogonal Grids for Ocean Models , 1996 .

[26]  R. Michael Clancy,et al.  A Note on Finite Differencing of the Advection-Diffusion Equation , 1981 .

[27]  W. R. Peltier,et al.  The Structure and Nonlinear Evolution of Synoptic Scale Cyclones: Life Cycle Simulations with a Cloud-Scale Model , 1990 .

[28]  Krishnan Mahesh,et al.  A numerical method for large-eddy simulation in complex geometries , 2004 .

[29]  Rainald Löhner,et al.  An unstructured grid-based, parallel free surface solver , 1997 .

[30]  A. Arakawa.  Computational design for long-term numerical integration of the equations of fluid motion: two-dimen , 1997 .

[31]  G. R. Stuhne,et al.  New icosahedral grid-point discretization of the shallow water equation on the sphere , 1999 .

[32]  Herbert E. Huppert,et al.  Axisymmetric gravity currents in a rotating system: experimental and numerical investigations , 2001, Journal of Fluid Mechanics.

[33]  R. A. Nicolaides.  Incompressible Computational Fluid Dynamics: The Covolume Approach to Computing Incompressible Flows , 1993 .

[34]  N. Wedi,et al.  Extending Gal-Chen and Somerville terrain-following coordinate transformation on time-dependent curvilinear boundaries , 2004 .

[35]  K. Bryan,et al.  A Nonlinear Model of an Ocean Driven by Wind and Differential Heating: Part I. Description of the Three-Dimensional Velocity and Density Fields , 1968 .

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

[37]  Sol Hellerman,et al.  Normal Monthly Wind Stress Over the World Ocean with Error Estimates , 1983 .

[38]  S. H. Lo,et al.  Finite element mesh generation over analytical curved surfaces , 1996 .

[39]  C. Reason,et al.  A Curvilinear Ocean Model Using a Grid Regionally Compressed in the South Indian Ocean , 2001 .

[40]  S. Manabe,et al.  A global ocean-atmosphere climate model with seasonal variation for future studies of climate sensitivity , 1979 .

[41]  M. Cox,et al.  A Nonlinear Model of an Ocean Driven by Wind and Differential Heating: Part II. An Analysis of the Heat, Vorticity and Energy Balance , 1968 .

[42]  Pascale Bouruet-Aubertot,et al.  Nonlinear adjustment of density fronts. Part 1. The Rossby scenario and the experimental reality , 2004, Journal of Fluid Mechanics.

[43]  Douglas K. Lilly,et al.  ON THE COMPUTATIONAL STABILITY OF NUMERICAL SOLUTIONS OF TIME-DEPENDENT NON-LINEAR GEOPHYSICAL FLUID DYNAMICS PROBLEMS , 1965 .

[44]  B. Perot.  Conservation Properties of Unstructured Staggered Mesh Schemes , 2000 .

[45]  Malcolm L. Spaulding,et al.  Estuarine and coastal modeling : proceedings of the 2nd international conference , 1992 .

[46]  Y. H. Yamazaki,et al.  The Existence of Subsynoptic-Scale Baroclinic Instability and the Nonlinear Evolution of Shallow Disturbances , 2001 .

[47]  P. Swarztrauber,et al.  A standard test set for numerical approximations to the shallow water equations in spherical geometry , 1992 .

[48]  S. Manabe,et al.  A Global Ocean-Atmosphere Climate Model. Part II. The Oceanic Circulation , 1975 .