A mass, energy, vorticity, and potential enstrophy conserving lateral fluid-land boundary scheme for the shallow water equations

A numerical scheme for treating fluid-land boundaries in inviscid shallow water flows is derived that conserves the domain-summed mass, energy, vorticity, and potential enstrophy in domains with arbitrarily shaped boundaries. The boundary scheme is derived from a previous scheme that conserves all four domain-summed quantities only in periodic domains without boundaries. It consists of a method for including land in the model along with evolution equations for the vorticity and extrapolation formulas for the depth at fluid-land boundaries. Proofs of mass, energy, vorticity, and potential enstrophy conservation are given. Numerical simulations are carried out demonstrating the conservation properties and accuracy of the boundary scheme for inviscid flows and comparing its performance with that of four alternative boundary schemes. The first of these alternatives extrapolates or finite-differences the velocity to obtain the vorticity at boundaries; the second enforces the free-slip boundary condition; the third enforces the super-slip condition; and the fourth enforces the no-slip condition. Comparisons of the conservation properties demonstrate that the new scheme is the only one of the five that conserves all four domain-summed quantities, and it is the only one that both prevents a spurious energy cascade to the smallest resolved scales and maintains the correct flow orientation with respect to an external forcing. Comparisons of the accuracy demonstrate that the new scheme generates vorticity fields that have smaller errors than those generated by any of the alternative schemes, and it generates depth and velocity fields that have errors about equal to those in the fields generated by the most accurate alternative scheme.

[1]  H. Hurlburt,et al.  A Numerical Study of Loop Current Intrusions and Eddy Shedding , 1980 .

[2]  Richard J. Greatbatch,et al.  Four-Gyre Circulation in a Barotropic Model with Double-Gyre Wind Forcing , 2000 .

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

[4]  Tomonori Matsuura,et al.  Two Different Aperiodic Phases of Wind-Driven Ocean Circulation in a Double-Gyre, Two-Layer Shallow-Water Model , 2006 .

[5]  Russell L. Elsberry,et al.  Interactions between a Barotropic Vortex and an Idealized Subtropical Ridge. Part I: Vortex Motion , 1991 .

[6]  R. Salmon,et al.  A General Method for Conserving Energy and Potential Enstrophy in Shallow-Water Models , 2007 .

[7]  Jian Wang,et al.  Balanced Models and Dynamics for the Large- and Mesoscale Circulation , 1997 .

[8]  François W. Primeau,et al.  Multiple Equilibria of a Double-Gyre Ocean Model with Super-Slip Boundary Conditions , 1998 .

[9]  R. Sadourny.  The Dynamics of Finite-Difference Models of the Shallow-Water Equations , 1975 .

[10]  G. Parisi,et al.  : Multiple equilibria , 2022 .

[11]  A. Adcroft,et al.  Representation of Topography by Shaved Cells in a Height Coordinate Ocean Model , 1997 .

[12]  Zavisa Janjic,et al.  Nonlinear Advection Schemes and Energy Cascade on Semi-Staggered Grids , 1984 .

[13]  G. Tripoli.  A Nonhydrostatic Mesoscale Model Designed to Simulate Scale Interaction , 1992 .

[14]  Jian Wang,et al.  Emergence of Fofonoff states in inviscid and viscous ocean circulation models , 1994 .

[15]  Todd D. Ringler,et al.  The ZM Grid: An Alternative to the Z Grid , 2002 .

[16]  Akio Arakawa,et al.  Computational Design of the Basic Dynamical Processes of the UCLA General Circulation Model , 1977 .

[17]  Rick Salmon.  Poisson-Bracket Approach to the Construction of Energy- and Potential-Enstrophy- Conserving Algorithms for the Shallow-Water Equations , 2004 .

[18]  D. Randall,et al.  A Potential Enstrophy and Energy Conserving Numerical Scheme for Solution of the Shallow-Water Equations on a Geodesic Grid , 2002 .

[19]  Akio Arakawa,et al.  Integration of the Nondivergent Barotropic Vorticity Equation with AN Icosahedral-Hexagonal Grid for the SPHERE1 , 1968 .

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

[21]  Frank Lunkeit,et al.  Regime Transitions in a Stochastically Forced Double-Gyre Model , 2001 .

[22]  K. W. Morton,et al.  Vorticity-Preserving Lax-Wendroff-Type Schemes for the System Wave Equation , 2001, SIAM J. Sci. Comput..

[23]  R. Salmon,et al.  A general method for conserving quantities related to potential vorticity in numerical models , 2005 .

[24]  Frank Abramopoulos.  A New Fourth-Order Enstrophy and Energy Conserving Scheme , 1991 .

[25]  Henk A. Dijkstra,et al.  On the Physics of the Agulhas Current: Steady Retroflection Regimes , 2001 .

[26]  A. Arakawa,et al.  A Potential Enstrophy and Energy Conserving Scheme for the Shallow Water Equations , 1981 .

[27]  Roger Temam,et al.  Low-Frequency Variability in Shallow-Water Models of the Wind-Driven Ocean Circulation. Part I: Steady-State Solution* , 2003 .

[28]  Michael Ghil,et al.  Multiple Equilibria, Periodic, and Aperiodic Solutions in a Wind-Driven, Double-Gyre, Shallow-Water Model , 1995 .

[29]  Jonas Nycander,et al.  Nonuniform Upwelling in a Shallow-Water Model of the Antarctic Bottom Water in the Brazil Basin* , 2004 .

[30]  T. Ringler,et al.  Analysis of Discrete Shallow-Water Models on Geodesic Delaunay Grids with C-Type Staggering , 2005 .

[31]  Slobodan Nickovic,et al.  The Step-Mountain Coordinate: Model Description and Performance for Cases of Alpine Lee Cyclogenesis and for a Case of an Appalachian Redevelopment , 1988 .

[32]  Francis J. Poulin,et al.  The Influence of Topography on the Stability of Jets , 2005 .

[33]  Miodrag Ic.  Fourth-Order Horizontal Advection Schemes on the Semi-staggered Grid , 1988 .

[34]  Fedor Mesinger,et al.  Horizontal Advection Schemes of a Staggered Grid—An Enstrophy and Energy-Conserving Model , 1981 .

[35]  George F. Carnevale,et al.  Dynamic Boundary Conditions Revisited , 2001 .

[36]  R. E. Hart,et al.  Simulations of Dual-Vortex Interaction within Environmental Shear , 1999 .

[37]  A. Hollingsworth,et al.  An internal symmetric computational instability , 1983 .

[38]  Frank Abramopoulos.  Generalized Energy and Potential Enstrophy Conserving Finite Difference Schemes for the Shallow Water Equations , 1988 .

[39]  Lynne D. Talley,et al.  Generalizations of Arakawa's Jacobian , 1989 .

[40]  Kenji Takano,et al.  A Fourth Order Energy and Potential Enstrophy Conserving Difference Scheme , 1982 .