Strongly enstrophy-stable integrators for the incompressible Navier-Stokes equations
Abstract
We propose a mixed finite element discretisation for the incompressible Navier-Stokes equations that preserves the evolution laws of both energy and enstrophy, in a stronger sense than previous discretisations. In two dimensions, the evolution law for enstrophy only permits dissipation for thermodynamically isolated systems, leading to a Reynolds-number-independent bound on the velocity gradient that naturally stabilises the scheme, even on severely under-resolved meshes. In three dimensions, the scheme preserves both dissipation and the generation of enstrophy through vortex stretching. We enforce these evolution laws by systematically introducing auxiliary variables into the discretisation. While conforming implementations of these schemes require discrete Stokes complexes with enhanced regularity, we introduce both (i) equivalent reparametrisations and (ii) penalty formulations that require only the typical curl- and div-conforming spaces from the standard discrete de Rham complex. The scheme handles different types of boundary conditions and curved domains. The robust stabilisation properties of the proposed scheme are demonstrated through numerical simulations of a shear flow, a spherical vortex, and flow past an obstacle. We observe numerically that preserving the discrete evolution of enstrophy in this way has a strong stabilising effect on the numerical solution, especially in two dimensions.
Related
- High-order Energy-stable and Charge-conservative Lagrangian FEM for 3D Incompressible Inductionless MHD equations with Variable Density
- Noncommutative Lp-integrators and stochastic differential equations in Lp
- Global Regularity for the 3D Incompressible Navier–Stokes Equations via the Frohmanian Symplectic Tether: Lean 4 Formalization and Conditional Cubic Comparison Framework (v2.1)
- Nonuniqueness of Solutions to 2D Incompressible Euler Equations with Unbounded Vorticity
- A dual-field structure-preserving mixed finite element discretization for incompressible Hall MHD equations
BibTeX
@misc{andrews2026strongly,
title = {{Strongly enstrophy-stable integrators for the incompressible Navier-Stokes equations}},
author = {Boris D. Andrews and Matin Shams and Patrick E. Farrell},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.15520},
eprint = {2609.15520},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.15520},
}
RIS
TY - UNPB
TI - Strongly enstrophy-stable integrators for the incompressible Navier-Stokes equations
AU - Boris D. Andrews
AU - Matin Shams
AU - Patrick E. Farrell
PY - 2026
DA - 2026/09/14
PB - arXiv
DO - 10.48550/arxiv.2609.15520
UR - https://arxiv.org/abs/2609.15520
AB - We propose a mixed finite element discretisation for the incompressible Navier-Stokes equations that preserves the evolution laws of both energy and enstrophy, in a stronger sense than previous discretisations. In two dimensions, the evolution law for enstrophy only permits dissipation for thermodynamically isolated systems, leading to a Reynolds-number-independent bound on the velocity gradient that naturally stabilises the scheme, even on severely under-resolved meshes. In three dimensions, the scheme preserves both dissipation and the generation of enstrophy through vortex stretching. We enforce these evolution laws by systematically introducing auxiliary variables into the discretisation. While conforming implementations of these schemes require discrete Stokes complexes with enhanced regularity, we introduce both (i) equivalent reparametrisations and (ii) penalty formulations that require only the typical curl- and div-conforming spaces from the standard discrete de Rham complex. The scheme handles different types of boundary conditions and curved domains. The robust stabilisation properties of the proposed scheme are demonstrated through numerical simulations of a shear flow, a spherical vortex, and flow past an obstacle. We observe numerically that preserving the discrete evolution of enstrophy in this way has a strong stabilising effect on the numerical solution, especially in two dimensions.
ER -