We investigate the temporal accuracy of two generalized‐ schemes for the incompressible Navier‐Stokes equations. In a widely‐adopted approach, the pressure is collocated at the time step
New SAV-pressure correction methods for the Navier-Stokes equations: stability and error analysis
We construct new first- and second-order pressure correctionschemes using the scalar auxiliary variable approach for the Navier-Stokes equations. These schemes are linear, decoupled and only require solving a sequence of Poisson type equations at each time step. Furthermore, they are unconditionally energy stable. We also establish rigorous error estimates in the two dimensional case for the velocity and pressure approximation of the first-order scheme without any condition on the time step.
more »
« less
- Award ID(s):
- 2012585
- NSF-PAR ID:
- 10329244
- Date Published:
- Journal Name:
- Mathematics of Computation
- Volume:
- 91
- Issue:
- 333
- ISSN:
- 0025-5718
- Page Range / eLocation ID:
- 141 to 167
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract t n + 1second‐order accurate in time . We developed a suite of numerical code using inf‐sup stable higher‐order non‐uniform rational B‐spline (NURBS) elements for spatial discretization. In doing so, we are able to achieve high spatial accuracy and to investigate asymptotic temporal convergence behavior. Numerical evidence suggests that onlyfirst‐order accuracy is achieved, at least for the pressure, in this aforesaid temporal discretization approach. On the other hand, evaluating the pressure at the intermediate time step recovers second‐order accuracy, and the numerical implementation is simplified. We recommend this second approach as the generalized‐ scheme of choice when integrating the incompressible Navier‐Stokes equations. -
null (Ed.)Pressure Poisson equation (PPE) reformulations of the incompressible Navier–Stokes equations (NSE) replace the incompressibility constraint by a Poisson equation for the pressure and a suitable choice of boundary conditions. This yields a time-evolution equation for the velocity field only, with the pressure gradient acting as a nonlocal operator. Thus, numerical methods based on PPE reformulations are representatives of a class of methods that have no principal limitations in achieving high order. In this paper, it is studied to what extent high-order methods for the NSE can be obtained from a specific PPE reformulation with electric boundary conditions (EBC). To that end, implicit–explicit (IMEX) time-stepping is used to decouple the pressure solve from the velocity update, while avoiding a parabolic time-step restriction; and mixed finite elements are used in space, to capture the structure imposed by the EBC. Via numerical examples, it is demonstrated that the methodology can yield at least third order accuracy in space and time.more » « less
-
Abstract In this work we introduce semi-implicit or implicit finite difference schemes for the continuity equation with a gradient flow structure. Examples of such equations include the linear Fokker–Planck equation and the Keller–Segel equations. The two proposed schemes are first-order accurate in time, explicitly solvable, and second-order and fourth-order accurate in space, which are obtained via finite difference implementation of the classical continuous finite element method. The fully discrete schemes are proved to be positivity preserving and energy dissipative: the second-order scheme can achieve so unconditionally while the fourth-order scheme only requires a mild time step and mesh size constraint. In particular, the fourth-order scheme is the first high order spatial discretization that can achieve both positivity and energy decay properties, which is suitable for long time simulation and to obtain accurate steady state solutions.more » « less
-
more » « less
Abstract The one‐leg, two‐step time discretization proposed by Dahlquist, Liniger and Nevanlinna is second order and variable step G‐stable. G‐stability for systems of ordinary differential equations (ODEs) corrresponds to unconditional, long time energy stability when applied to the Navier–Stokes equations (NSEs). In this report, we analyze the method of Dahlquist, Liniger and Nevanlinna as a variable step, time discretization of the Navier–Stokes equations. We prove that the kinetic energy is bounded for variable time‐steps, show that the method is second‐order accurate, characterize its numerical dissipation and prove error estimates. The theoretical results are illustrated by several numerical tests.
-
Abstract This paper concerns the existence of global weak solutions à la Leray for compressible Navier–Stokes equations with a pressure law which depends on the density and on time and space variables t and x . The assumptions on the pressure contain only locally Lipschitz assumption with respect to the density variable and some hypothesis with respect to the extra time and space variables. It may be seen as a first step to consider heat-conducting Navier–Stokes equations with physical laws such as the truncated virial assumption. The paper focuses on the construction of approximate solutions through a new regularized and fixed point procedure and on the weak stability process taking advantage of the new method introduced by the two first authors with a careful study of an appropriate regularized quantity linked to the pressure.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1090/mcom/3651
