Energy conservation issues in sigma-coordinate free-surface ocean models

This paper focuses on the energy conservation properties of a hydrostatic, Boussinesq, coastal ocean model using a classic finite difference method. It is shown that the leapfrog time-stepping scheme, combined with the sigma-coordinate formalism and the motions of the free surface, prevents the momentum advection from exactly conserving energy. Because of the leapfrog scheme, the discrete form of the kinetic energy depends on the product of velocities at odd and even time steps and thus appears to be possibly negative when high-frequency modes develop. Besides, the study of the energy balance clarifies the numerical choices made for the computation of mixing processes. The time-splitting technique used to reduce the computation costs associated to the resolution of surface waves leads to the well-known external and internal mode equations. We show that these equations do not conserve energy if the coupling of these two modes is forward in time. Even if non-linear terms are negligible, this shortcoming can be significant regarding the pressure gradient term ‘frozen' over a baroclinic time step. An alternative energy-conserving time-splitting technique is proposed in this paper. Discussion and conclusions are conducted in the light of a set of numerical experiments dedicated to surface and internal gravity waves.

[1]  P. Marsaleix,et al.  Secondary flows induced by wind forcing in the Rhône region of freshwater influence , 2004 .

[2]  Robert L. Higdon,et al.  A two-level time-stepping method for layered ocean circulation models: further development and testing , 2005 .

[3]  G. J. Haltiner.  Numerical Prediction and Dynamic Meteorology , 1980 .

[4]  Dale B. Haidvogel,et al.  Numerical Simulation of Flow around a Tall Isolated Seamount. Part I: Problem Formulation and Model Accuracy , 1993 .

[5]  P. Marsaleix,et al.  Application of an Inverse Method to Coastal Modeling , 2000 .

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

[7]  Rüdiger Gerdes,et al.  A primitive equation ocean circulation model using a general vertical coordinate transformation: 1. Description and testing of the model , 1993 .

[8]  I. D. James,et al.  Advection schemes for shelf sea models , 1996 .

[9]  Richard Asselin,et al.  Frequency Filter for Time Integrations , 1972 .

[10]  James C. McWilliams,et al.  A method for computing horizontal pressure‐gradient force in an oceanic model with a nonaligned vertical coordinate , 2003 .

[11]  G. Mellor,et al.  Development of a turbulence closure model for geophysical fluid problems , 1982 .

[12]  Y. Song,et al.  A General Pressure Gradient Formulation for Ocean Models. Part II: Energy, Momentum, and Bottom Torque Consistency , 1998 .

[13]  A. Arakawa,et al.  Vertical Differencing of the Primitive Equations in Sigma Coordinates , 1983 .

[14]  Patrick Marsaleix,et al.  Considerations on Open Boundary Conditions for Regional and Coastal Ocean Models , 2006 .

[15]  S. Griffies,et al.  Tracer Conservation with an Explicit Free Surface Method for z-Coordinate Ocean Models , 2001 .

[16]  N. Heaps,et al.  Three-dimensional coastal ocean models , 1987 .

[17]  U. C. Mohanty,et al.  Simulation of storm surges using a three‐dimensional numerical model: An application to the 1977 Andhra cyclone , 1983 .

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

[19]  Philippe Gaspar,et al.  A simple eddy kinetic energy model for simulations of the oceanic vertical mixing: Tests at Station Papa and long-term upper ocean study site , 1990 .

[20]  Kraig B. Winters.  Available potential energy and mixing in density-stratified fluids , 1996 .

[21]  I. Shulman,et al.  Optimized Boundary Conditions and Data Assimilation with Application to the M2 Tide in the Yellow Sea , 1998 .

[22]  Y. Song,et al.  A General Pressure Gradient Formulation for Ocean Models. Part I: Scheme Design and Diagnostic Analysis , 1998 .

[23]  P. Delecluse,et al.  A Three-Dimensional Numerical Study of Deep-Water Formation in the Northwestern Mediterranean Sea , 1991 .

[24]  G. Mellor.  USERS GUIDE for A THREE-DIMENSIONAL, PRIMITIVE EQUATION, NUMERICAL OCEAN MODEL , 1998 .

[25]  Robert L. Haney,et al.  On the Pressure Gradient Force over Steep Topography in Sigma Coordinate Ocean Models , 1991 .

[26]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[27]  Eric Deleersnijder,et al.  Numerical Mass Conservation in a Free-surface Sigma-coordinate Marine Model With Mode Splitting , 1993 .

[28]  A. E. Gill.  Atmosphere-Ocean Dynamics , 1982 .

[29]  Hans Burchard.  Energy-conserving discretisation of turbulent shear and buoyancy production , 2002 .

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

[31]  Steven A. Orszag,et al.  Large Eddy Simulation of Complex Engineering and Geophysical Flows , 2010 .