Mass, momentum and energy conserving (MaMEC) discretizations on general grids for the compressible Euler and shallow water equations

The paper explains a method by which discretizations of the continuity and momentum equations can be designed, such that they can be combined with an equation of state into a discrete energy equation. The resulting 'MaMEC' discretizations conserve mass, momentum as well as energy, although no explicit conservation law for the total energy is present. Essential ingredients are (i) discrete convection that leaves the discrete energy invariant, and (ii) discrete consistency between the thermodynamic terms. Of particular relevance is the way in which finite volume fluxes are related to nodal values. The method is an extension of existing methods based on skew-symmetry of discrete operators, because it allows arbitrary equations of state and a larger class of grids than earlier methods. The method is first illustrated with a one-dimensional example on a highly stretched staggered grid, in which the MaMEC method calculates qualitatively correct results and a non-skew-symmetric finite volume method becomes unstable. A further example is a two-dimensional shallow water calculation on a rectilinear grid as well as on an unstructured grid. The conservation of mass, momentum and energy is checked, and losses are found negligible up to machine accuracy.

[1]  Arthur Veldman,et al.  Spectro-Consistent Discretization of Navier-Stokes: a Challenge to RANS and LES , 1998 .

[2]  Jason Frank,et al.  Conservation Properties of Smoothed Particle Hydrodynamics Applied to the Shallow Water Equation , 2001 .

[3]  Arthur Veldman,et al.  A Volume-of-Fluid based simulation method for wave impact problems , 2005 .

[4]  Luca Bonaventura,et al.  A semi‐implicit method conserving mass and potential vorticity for the shallow water equations on the sphere , 2005 .

[5]  J. Blair Perot,et al.  Discrete calculus methods for diffusion , 2007, J. Comput. Phys..

[6]  Roel Verstappen,et al.  A new symmetry‐preserving Cartesian‐grid method for computing flow past arbitrarily shaped objects , 2005 .

[7]  A. Gorobets,et al.  Parameter-free symmetry-preserving regularization modeling of a turbulent differentially heated cavity , 2010 .

[8]  Bülent Karasözen,et al.  Cosymmetry preserving finite-difference methods for convection equations in a porous medium , 2005 .

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

[10]  Vincent Guinot,et al.  High-Order Fluxes for Conservative Skew-Symmetric-like Schemes in Structured Meshes , 2000 .

[11]  Jason Frank,et al.  Statistical relevance of vorticity conservation in the Hamiltonian particle-mesh method , 2009, J. Comput. Phys..

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

[13]  Jian-Guo Liu,et al.  Energy and helicity preserving schemes for hydro- and magnetohydro-dynamics flows with symmetry , 2004 .

[14]  Maarten van Reeuwijk,et al.  A mimetic mass, momentum and energy conserving discretization for the shallow water equations , 2011 .

[15]  Antony Jameson,et al.  The Construction of Discretely Conservative Finite Volume Schemes that Also Globally Conserve Energy or Entropy , 2008, J. Sci. Comput..

[16]  Arthur Veldman,et al.  Symmetry-preserving upwind discretization of convection on non-uniform grids , 2008 .

[17]  J. Blair Perot,et al.  Higher-order mimetic methods for unstructured meshes , 2006, J. Comput. Phys..

[18]  Roel Verstappen.  On restraining the production of small scales of motion in a turbulent channel flow , 2008 .

[19]  Oleg V. Vasilyev.  High Order Finite Difference Schemes on Non-uniform Meshes with Good Conservation Properties , 2000 .

[20]  Yohei Morinishi,et al.  Skew-symmetric form of convective terms and fully conservative finite difference schemes for variable density low-Mach number flows , 2010, J. Comput. Phys..

[21]  Stanly Steinberg.  A Discreate Calculus with Applications of High-Order Discretizations to Boundary-Value Problems , 2004 .

[22]  Thomas S. Lund,et al.  Kinetic energy conservation issues associated with the collocated mesh scheme for incompressible flow , 2006, J. Comput. Phys..

[23]  B. Strand.  Summation by parts for finite difference approximations for d/dx , 1994 .

[24]  S. Steinberg.  A Discrete Calculus with Applications of High-Order Discretizations to Boundary-Value Problems , 2004 .

[25]  Pramod K. Subbareddy,et al.  A fully discrete, kinetic energy consistent finite-volume scheme for compressible flows , 2009, J. Comput. Phys..

[26]  M. Shashkov,et al.  The Construction of Compatible Hydrodynamics Algorithms Utilizing Conservation of Total Energy , 1998 .

[27]  A. Veldman,et al.  Symmetry-preserving discretization of turbulent flow , 2003 .

[28]  J. C. Kok,et al.  A high-order low-dispersion symmetry-preserving finite-volume method for compressible flow on curvilinear grids , 2009, J. Comput. Phys..

[29]  Gareth P. Williams,et al.  Conservation properties of convection difference schemes , 1970 .

[30]  P. Moin,et al.  Fully Conservative Higher Order Finite Difference Schemes for Incompressible Flow , 1998 .

[31]  A. Arakawa,et al.  Numerical methods used in atmospheric models , 1976 .

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

[33]  Rik Wemmenhove,et al.  Numerical simulation of two-phase flow in offshore environments , 2008 .

[34]  Roel Luppes,et al.  The numerical simulation of liquid sloshing on board spacecraft , 2007, J. Comput. Phys..

[35]  Joel H. Ferziger,et al.  Numerical simulation of a compressible, homogeneous, turbulent shear flow , 1981 .

[36]  J. Blair Perot,et al.  Discrete Conservation Properties of Unstructured Mesh Schemes , 2011 .

[37]  Arthur Veldman,et al.  Direct Numerical Simulation of Turbulence at Lower Costs , 1997 .

[38]  Guus S. Stelling,et al.  A staggered conservative scheme for every Froude number in rapidly varied shallow water flows , 2003 .