Stability and consistency of nonhydrostatic free‐surface models using the semi‐implicit θ‐method

SUMMARY The θ-method is a popular semi-implicit finite-difference method for simulating free-surface flows. Problem stiffness, arising because of the presence of both fast and slow timescale processes, is easily handled by the θ-method. In most ocean, coastal, and estuary modeling applications, stiffness is caused by fast surface gravity wave timescales imposed on slower timescales of baroclinic variability. The method is well known to be unconditionally stable for shallow water (hydrostatic) models when 12≤ θ ≤1, where θ is the implicitness parameter. In this paper, we demonstrate that the method is also unconditionally stable for nonhydrostatic models, when 12≤ θ ≤1 for both pressure projection and pressure correction methods. However, the methods result in artificial damping of the barotropic mode. In addition to investigating stability, we also estimate the form of artificial damping induced by both the free surface and nonhydrostatic pressure solution methods. Finally, this analysis may be used to estimate the damping or growth associated with a particular wavenumber and the overall order of accuracy of the discretization. Copyright © 2012 John Wiley & Sons, Ltd.

[1]  O. Fringer,et al.  Physical vs. numerical dispersion in nonhydrostatic ocean modeling , 2011 .

[2]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[3]  Alexander F. Shchepetkin,et al.  The regional oceanic modeling system (ROMS): a split-explicit, free-surface, topography-following-coordinate oceanic model , 2005 .

[4]  G. Stelling,et al.  On a momentum conservative z-layer unstructured C-grid ocean model with flooding , 2012 .

[5]  G. Dahlquist.  A special stability problem for linear multistep methods , 1963 .

[6]  P. Moin.  Fundamentals of Engineering Numerical Analysis , 2001 .

[7]  Chin H. Wu,et al.  A two‐dimensional vertical non‐hydrostatic σ model with an implicit method for free‐surface flows , 2004 .

[8]  António M. Baptista,et al.  SELFE: A semi-implicit Eulerian–Lagrangian finite-element model for cross-scale ocean circulation , 2008 .

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

[10]  Yefei Bai,et al.  Depth‐integrated free‐surface flow with a two‐layer non‐hydrostatic formulation , 2012 .

[11]  Time Accuracy of Pressure Methods for Nonhydrostatic Free-Surface Flows , 2006 .

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

[13]  L. Perelman,et al.  A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers , 1997 .

[14]  B. Hodges.  Numerical Techniques in CWR-ELCOM (code release v.1) , 2000 .

[15]  S. Mitran,et al.  An approximated method for the solution of elliptic problems in thin domains: Application to nonlinear internal waves , 2008 .

[16]  L. Mysak,et al.  Waves in the ocean , 1978 .

[17]  Vincenzo Casulli,et al.  A SEMI-IMPLICIT FINITE DIFFERENCE METHOD FOR NON-HYDROSTATIC, FREE-SURFACE FLOWS , 1999 .

[18]  Alexander F. Shchepetkin,et al.  Algorithm for non-hydrostatic dynamics in the Regional Oceanic Modeling System , 2007 .

[19]  M. Gerritsen,et al.  An unstructured-grid, finite-volume, nonhydrostatic, parallel coastal ocean simulator , 2006 .

[20]  R. Falconer,et al.  Three‐dimensional numerical modelling of free surface flows with non‐hydrostatic pressure , 2002 .

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

[22]  Geoffrey W. Cowles,et al.  A nonhydrostatic version of FVCOM: 1. Validation experiments , 2010 .

[23]  Y. Kanarska,et al.  A non-hydrostatic numerical model for calculating free-surface stratified flows , 2003 .

[24]  V. Casulli,et al.  Stability, accuracy and efficiency of a semi-implicit method for three-dimensional shallow water flow☆ , 1994 .

[25]  G. Stelling,et al.  Numerical simulation of the vertical structure of discontinuous flows , 2001 .

[26]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[27]  J. Crank,et al.  A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type , 1947 .

[28]  G. Stelling,et al.  Development and validation of a three-dimensional morphological model , 2004 .

[29]  H. Hasumi,et al.  Developments in ocean climate modelling , 2000 .

[30]  Luca Bonaventura,et al.  Consistency with continuity in conservative advection schemes for free‐surface models , 2002 .

[31]  A three‐dimensional non‐hydrostatic vertical boundary fitted model for free‐surface flows , 2008 .

[32]  J. Koseff,et al.  A non-staggered grid, fractional step method for time-dependent incompressible Navier-Stokes equations in curvilinear coordinates , 1994 .

[33]  Roger Alexander Falconer,et al.  An implicit numerical algorithm for solving non-hydrostatic free-surface flow problems , 2001 .

[34]  A. Chorin.  Numerical Solution of the Navier-Stokes Equations* , 1989 .

[35]  Vincenzo Casulli,et al.  Semi-implicit numerical modeling of nonhydrostatic free-surface flows for environmental problems , 2002 .

[36]  Xinjian Chen.  A fully hydrodynamic model for three‐dimensional, free‐surface flows , 2003 .

[37]  Robert L. Street,et al.  An analysis and comparison of the time accuracy of fractional‐step methods for the Navier–Stokes equations on staggered grids , 2002 .

[38]  Joseph Oliger,et al.  A nonhydrostatic mesoscale ocean model. Part II: Numerical implementation , 1996 .

[39]  R. T. Cheng,et al.  SEMI-IMPLICIT FINITE DIFFERENCE METHODS FOR THREE-DIMENSIONAL SHALLOW WATER FLOW , 1992 .

[40]  H. H. Rachford,et al.  The Numerical Solution of Parabolic and Elliptic Differential Equations , 1955 .

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

[42]  J. Kan.  A second-order accurate pressure correction scheme for viscous incompressible flow , 1986 .

[43]  A. Blumberg,et al.  A Description of a Three‐Dimensional Coastal Ocean Circulation Model , 2013 .

[44]  Marcel Zijlema,et al.  Further experiences with computing non‐hydrostatic free‐surface flows involving water waves , 2005 .