An Energy-Stable Scheme for Incompressible Navier-Stokes Equations with Periodically Updated Coefficient Matrix
arXiv:1909.08395 · doi:10.1016/j.jcp.2020.109624
Abstract
We present an energy-stable scheme for simulating the incompressible Navier-Stokes equations based on the generalized Positive Auxiliary Variable (gPAV) framework. In the gPAV-reformulated system the original nonlinear term is replaced by a linear term plus a correction term, where the correction term is put under control by an auxiliary variable. The proposed scheme incorporates a pressure-correction type strategy into the gPAV procedure, and it satisfies a discrete energy stability property. The scheme entails the computation of two copies of the velocity and pressure within a time step, by solving an individual de-coupled linear equation for each of these field variables. Upon discretization the pressure linear system involves a constant coefficient matrix that can be pre-computed, while the velocity linear system involves a coefficient matrix that is updated periodically, once every time steps in the current work, where is a user-specified integer. The auxiliary variable, being a scalar-valued number, is computed by a well-defined explicit formula, which guarantees the positivity of its computed values. It is observed that the current method can produce accurate simulation results at large (or fairly large) time step sizes for the incompressible Navier-Stokes equations. The impact of the periodic coefficient-matrix update on the overall cost of the method is observed to be small in typical numerical simulations. Several flow problems have been simulated to demonstrate the accuracy and performance of the method developed herein.
arXiv comment34 pages, 11 figures, 10 tables
References in corpus (6)
- Numerical Approximation of Incompressible Navier-Stokes Equations Based on an Auxiliary Energy Variable
- A Roadmap for Discretely Energy-Stable Schemes for Dissipative Systems Based on a Generalized Auxiliary Variable with Guaranteed Positivity
- An Unconditionally Energy-Stable Scheme Based on an Implicit Auxiliary Energy Variable for Incompressible Two-Phase Flows with Different Densities Involving Only Precomputable Coefficient Matrices
- A Convective-like Energy-Stable Open Boundary Condition for Simulations of Incompressible Flows
- A pressure correction scheme for generalized form of energy-stable open boundary conditions for incompressible flows
- A gPAV-Based Unconditionally Energy-Stable Scheme for Incompressible Flows with Outflow/Open Boundaries