Quantification of spurious dissipation and mixing – Discrete variance decay in a Finite-Volume framework

Abstract It is well known that in numerical models the advective transport relative to fixed or moving grids needs to be discretised with sufficient accuracy to minimise the spurious decay of tracer variance (spurious mixing). In this paper a general analysis of discrete variance decay (DVD) caused by advective and diffusive fluxes is established. Lacking a general closed derivation for the local DVD rate, two non-invasive methods to estimate local DVD during model runtime are discussed. Whereas the first was presented recently by Burchard and Rennau (2008), the second is a newly proposed alternative. This alternative analysis method is argued to have a more consistent foundation. In particular, it recovers a physically sound definition of discrete variance in a Finite-Volume cell. The diagnosed DVD can be separated into physical and numerical (spurious) contributions, with the latter originating from discretisation errors. Based on the DVD analysis, a 3D dissipation analysis is developed to quantify the physically and numerically induced loss of kinetic energy. This dissipation analysis provides a missing piece of information to assess the discrete energy conservation of an ocean model. Analyses are performed and evaluated for three test cases, with complexities ranging from idealised 1D advection to a realistic ocean modelling application to the Western Baltic Sea. In all test cases the proposed alternative DVD analysis method is demonstrated to provide a reliable diagnostic tool for the local quantification of physically and numerically induced dissipation and mixing.

[1]  Mehmet Ilicak,et al.  Spurious dianeutral mixing and the role of momentum closure , 2012 .

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

[3]  Hans Burchard,et al.  Second-order turbulence closure models for geophysical boundary layers. A review of recent work , 2005 .

[4]  K. Bolding,et al.  Three-dimensional modelling of estuarine turbidity maxima in a tidal estuary , 2004 .

[5]  Julia Getzlaff,et al.  Diagnostics of diapycnal diffusion in z-level ocean models , 2008 .

[6]  R. Feistel,et al.  Vertical mixing in the Baltic Sea and consequences for eutrophication - A review , 2009 .

[7]  O. Reynolds.  On the dynamical theory of incompressible viscous fluids and the determination of the criterion , 1995, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences.

[8]  Gareth P. Williams.  Friction Term Formulation and Convective Instability in a Shallow Atmosphere , 1972 .

[9]  H. Burchard,et al.  Quantitative analysis of numerically induced mixing in a coastal model application , 2009 .

[10]  James C. McWilliams,et al.  Quasi-Monotone Advection Schemes Based on Explicit Locally Adaptive Dissipation , 1998 .

[11]  K. Nagel,et al.  Hydrographisch-chemische Zustandseinschätzung der Ostsee 2004 , 2003 .

[12]  K. Bolding,et al.  Dynamics of medium-intensity dense water plumes in the Arkona Basin, Western Baltic Sea , 2005 .

[13]  U. Gräwe,et al.  The future of the western Baltic Sea: two possible scenarios , 2013, Ocean Dynamics.

[14]  S. Schimmels,et al.  On the effect of structure-induced resistance and mixing on inflows into the Baltic Sea: A numerical model study , 2012 .

[15]  K. Taylor.  Summarizing multiple aspects of model performance in a single diagram , 2001 .

[16]  Jean-Marie Beckers,et al.  Non-uniform adaptive vertical grids for 3D numerical ocean models , 2010 .

[17]  James C. McWilliams,et al.  Anisotropic horizontal viscosity for ocean models , 2003 .

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

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

[20]  M. Prather.  Numerical advection by conservation of second-order moments. [for trace element spatial distribution and chemical interaction in atmosphere] , 1986 .

[21]  U. Gräwe,et al.  Storm surges in the Western Baltic Sea: the present and a possible future , 2012, Climate Dynamics.

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

[23]  U. Gräwe,et al.  Deep‐water dynamics and boundary mixing in a nontidal stratified basin: A modeling study of the Baltic Sea , 2014 .

[24]  Stephen M. Griffies,et al.  Spurious Diapycnal Mixing Associated with Advection in a z-Coordinate Ocean Model , 2000 .

[25]  Arthur E. P. Veldman.  "Missing" Boundary Conditions? Discretize First, Substitute Next, and Combine Later , 1990, SIAM J. Sci. Comput..

[26]  Knut Klingbeil,et al.  Implementation of a direct nonhydrostatic pressure gradient discretisation into a layered ocean model , 2013 .

[27]  J. Pietrzak,et al.  The Use of TVD Limiters for Forward-in-Time Upstream-Biased Advection Schemes in Ocean Modeling , 1998 .

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

[29]  M. A. Morales Maqueda,et al.  Second-order moment advection scheme applied to Arctic Ocean simulation , 2006 .

[30]  Andreas Oschlies,et al.  Diagnostics of diapycnal diffusivity in z-level ocean models part I: 1-Dimensional case studies , 2010 .

[31]  Robert L. Street,et al.  Reducing numerical diffusion in interfacial gravity wave simulations , 2005 .

[32]  R. C. Wajsowicz.  A consistent formulation of the anisotropic stress tensor for use in models of the large-scale ocean circulation , 1993 .

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

[34]  Quantification of Numerical and Physical Mixing in Coastal Ocean Model Applications , 2013 .

[35]  K. Bolding,et al.  Model simulations of dense bottom currents in the Western Baltic Sea , 2009 .

[36]  P. Linden,et al.  The stability of vortices in a rotating, stratified fluid , 1981, Journal of Fluid Mechanics.

[37]  R. Feistel,et al.  State and Evolution of the Baltic Sea, 1952–2005 , 2008 .

[38]  Burchard Hans,et al.  GETM, a General Estuarine Transport Model , 2002 .

[39]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[40]  L. Umlauf,et al.  Internal wave mixing in the Baltic Sea: Near-inertial waves in the absence of tides , 2011 .

[41]  H. Burchard,et al.  Comparative quantification of physically and numerically induced mixing in ocean models , 2008 .

[42]  Hans Burchard.  Quantification of numerically induced mixing and dissipation in discretisations of shallow water equations , 2012 .

[43]  K. Ruddick,et al.  A coastal ocean model intercomparison study for a three-dimensional idealised test case , 1998 .

[44]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[45]  A. Kirwan.  Formulation of constitutive equations for large‐scale turbulent mixing , 1969 .

[46]  M. S. Dubovikov,et al.  Ocean Turbulence I: One-Point Closure Model Momentum and Heat Vertical Diffusivities , 2001 .

[47]  Ivane Pairaud,et al.  Energy conservation issues in sigma-coordinate free-surface ocean models , 2008 .

[48]  J. Beckers,et al.  Realistic modelling of the exceptional inflows into the central Baltic Sea in 2003 using terrain-following coordinates , 2011 .