Entropy-stable discontinuous Galerkin approximation with summation-by-parts property for the incompressible Navier-Stokes equations with variable density and artificial compressibility

Abstract We present a provably stable discontinuous Galerkin spectral element method for the incompressible Navier–Stokes equations with artificial compressibility and variable density. Stability proofs, which include boundary conditions, that follow a continuous entropy analysis are provided. We define a mathematical entropy function that combines the traditional kinetic energy and an additional energy term for the artificial compressibility, and derive its associated entropy conservation law. The latter allows us to construct a provably stable split–form nodal Discontinuous Galerkin (DG) approximation that satisfies the summation–by–parts simultaneous–approximation–term (SBP–SAT) property. The scheme and the stability proof are presented for general curvilinear three–dimensional hexahedral meshes. We use the exact Riemann solver and the Bassi–Rebay 1 (BR1) scheme at the inter–element boundaries for inviscid and viscous fluxes respectively, and an explicit low storage Runge–Kutta RK3 scheme to integrate in time. We assess the accuracy and robustness of the method by solving a manufactured solution, the Kovasznay flow, a lid driven cavity, the inviscid Taylor–Green vortex, and the Rayleigh–Taylor instability.

[1]  Gregor Gassner,et al.  Split form nodal discontinuous Galerkin schemes with summation-by-parts property for the compressible Euler equations , 2016, J. Comput. Phys..

[2]  David A. Kopriva,et al.  Metric Identities and the Discontinuous Spectral Element Method on Curvilinear Meshes , 2006, J. Sci. Comput..

[3]  D. Jacqmin.  Regular Article: Calculation of Two-Phase Navier–Stokes Flows Using Phase-Field Modeling , 1999 .

[4]  Esteban Ferrer,et al.  A high order Discontinuous Galerkin - Fourier incompressible 3D Navier-Stokes solver with rotating sliding meshes , 2012, J. Comput. Phys..

[5]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[6]  Travis C. Fisher,et al.  High-order entropy stable finite difference schemes for nonlinear conservation laws: Finite domains , 2013, J. Comput. Phys..

[7]  Chunlei Liang,et al.  A high-order solver for unsteady incompressible Navier-Stokes equations using the flux reconstruction method on unstructured grids with implicit dual time stepping , 2015, J. Comput. Phys..

[8]  Esteban Ferrer,et al.  An interior penalty stabilised incompressible discontinuous Galerkin-Fourier solver for implicit large eddy simulations , 2017, J. Comput. Phys..

[9]  Christopher A. Kennedy,et al.  Reduced aliasing formulations of the convective terms within the Navier-Stokes equations for a compressible fluid , 2008, J. Comput. Phys..

[11]  Gonzalo Rubio,et al.  Entropy-stable discontinuous Galerkin approximation with summation-by-parts property for the incompressible Navier-Stokes/Cahn-Hilliard system , 2019, J. Comput. Phys..

[12]  J. Peraire,et al.  Subgrid-scale modeling and implicit numerical dissipation in DG-based Large-Eddy Simulation , 2017 .

[13]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[14]  Esteban Ferrer,et al.  A free-energy stable nodal discontinuous Galerkin approximation with summation-by-parts property for the Cahn-Hilliard equation , 2019, J. Comput. Phys..

[15]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[16]  J. Williamson.  Low-storage Runge-Kutta schemes , 1980 .

[17]  Esteban Ferrer,et al.  The Bassi Rebay 1 scheme is a special case of the Symmetric Interior Penalty formulation for discontinuous Galerkin discretisations with Gauss-Lobatto points , 2018, J. Comput. Phys..

[18]  Steven H. Frankel,et al.  Entropy Stable Spectral Collocation Schemes for the Navier-Stokes Equations: Discontinuous Interfaces , 2014, SIAM J. Sci. Comput..

[19]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[20]  Eitan Tadmor,et al.  Entropy stability theory for difference approximations of nonlinear conservation laws and related time-dependent problems , 2003, Acta Numerica.

[21]  L. Kovasznay.  Laminar flow behind a two-dimensional grid , 1948 .

[22]  G. Karniadakis,et al.  Spectral/hp Element Methods for Computational Fluid Dynamics , 2005 .

[23]  Gregor Gassner,et al.  On the use of kinetic energy preserving DG-schemes for large eddy simulation , 2017, J. Comput. Phys..

[24]  S. Rebay,et al.  A High-Order Accurate Discontinuous Finite Element Method for the Numerical Solution of the Compressible Navier-Stokes Equations , 1997 .

[25]  Spencer J. Sherwin,et al.  Stability of Projection Methods for Incompressible Flows Using High Order Pressure-Velocity Pairs of Same Degree: Continuous and Discontinuous Galerkin Formulations , 2014 .

[26]  Jie Shen,et al.  PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS , 2007 .

[27]  Nikolaus A. Adams,et al.  A physically consistent weakly compressible high-resolution approach to underresolved simulations of incompressible flows , 2013 .

[28]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .

[29]  J. Sethian,et al.  A Fast Level Set Method for Propagating Interfaces , 1995 .

[30]  L. Quartapelle,et al.  A projection FEM for variable density incompressible flows , 2000 .

[31]  Gonzalo Rubio,et al.  On the role of numerical dissipation in stabilising under-resolved turbulent simulations using discontinuous Galerkin methods , 2018 .

[32]  Florian J. Hindenlang,et al.  Stability of Wall Boundary Condition Procedures for Discontinuous Galerkin Spectral Element Approximations of the Compressible Euler Equations , 2019, Lecture Notes in Computational Science and Engineering.

[33]  Héctor D. Ceniceros,et al.  Computation of multiphase systems with phase field models , 2002 .

[34]  Alessandro Colombo,et al.  Artificial compressibility Godunov fluxes for variable density incompressible flows , 2017, Computers & Fluids.

[35]  C. W. Hirt,et al.  SOLA-VOF: a solution algorithm for transient fluid flow with multiple free boundaries , 1980 .

[36]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[37]  Gregor Gassner,et al.  The BR1 Scheme is Stable for the Compressible Navier–Stokes Equations , 2017, J. Sci. Comput..