A gPAV-Based Unconditionally Energy-Stable Scheme for Incompressible Flows with Outflow/Open Boundaries
arXiv:1908.01852 · doi:10.1016/j.cma.2020.112969
Abstract
We present an unconditionally energy-stable scheme for approximating the incompressible Navier-Stokes equations on domains with outflow/open boundaries. The scheme combines the generalized Positive Auxiliary Variable (gPAV) approach and a rotational velocity-correction type strategy, and the adoption of the auxiliary variable simplifies the numerical treatment for the open boundary conditions. The discrete energy stability of the proposed scheme has been proven, irrespective of the time step sizes. Within each time step the scheme entails the computation of two velocity fields and two pressure fields, by solving an individual de-coupled Helmholtz (including Poisson) type equation with a constant pre-computable coefficient matrix for each of these field variables. The auxiliary variable, being a scalar number, is given by a well-defined explicit formula within a time step, which ensures the positivity of its computed values. Extensive numerical experiments with several flows involving outflow/open boundaries in regimes where the backflow instability becomes severe have been presented to test the performance of the proposed method and to demonstrate its stability at large time step sizes.
34 pages, 18 figures, 6 tables
References in corpus (8)
- 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
- Wall-bounded multiphase flows of N immiscible incompressible fluids: consistency and contact-angle boundary condition
- 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
- Energy-Stable Boundary Conditions Based on a Quadratic Form: Applications to Outflow/Open-Boundary Problems in Incompressible Flows
- A Three-Dimensional Hybrid Spectral Element-Fourier Spectral Method for Wall-Bounded Two-Phase Flows