Local conservation laws of continuous Galerkin method for the incompressible Navier--Stokes equations in EMAC form

Maxim A. Olshanskii,Leo G. Rebholz

DOI: https://doi.org/10.1016/j.cma.2023.116583

2024-08-23

Abstract:We consider local balances of momentum and angular momentum for the incompressible Navier-Stokes equations. First, we formulate new weak forms of the physical balances (conservation laws) of these quantities, and prove they are equivalent to the usual conservation law formulations. We then show that continuous Galerkin discretizations of the Navier-Stokes equations using the EMAC form of the nonlinearity preserve discrete analogues of the weak form conservation laws, both in the Eulerian formulation and the Lagrangian formulation (which are not equivalent after discretizations). Numerical tests illustrate the new theory.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: How to ensure that the continuous Galerkin method can maintain the local conservation properties of momentum and angular momentum when solving the incompressible Navier - Stokes equations (NSE). Specifically, the author has studied whether using the nonlinear terms in the EMAC (energy, momentum and angular momentum conservation) form can accurately maintain these physical conservation laws at the discrete level.

### Background and Problem Description

The incompressible Navier - Stokes equations (NSE) are the fundamental equations in fluid mechanics and are used to describe the motion of viscous incompressible fluids. Their general form is:

\[
\begin{cases}
\frac{\partial \mathbf{u}}{\partial t}+(\mathbf{u} \cdot \nabla) \mathbf{u}-\text{div} \, \sigma = 0 \\
\text{div} \, \mathbf{u} = 0
\end{cases} \quad \text{in } \Omega
\]

Here, $\mathbf{u}$ is the velocity field, $\sigma$ is the Cauchy stress tensor, and for Newtonian fluids, $\sigma = 2\nu D(\mathbf{u})-pI$, where $D(\mathbf{u})=\frac{1}{2}(\nabla \mathbf{u}+(\nabla \mathbf{u})^T)$ is the deformation rate tensor and $p$ is the pressure field.

It is known that smooth solutions satisfy a series of conservation laws, such as the conservation of momentum, energy and vorticity. However, numerical methods usually break these conservation laws during the discretization process. Therefore, developing numerical methods that can maintain as many conservation laws as possible at the discrete level is a long - term challenge in the field of computational fluid dynamics.

### The Core Problem of the Paper

The main problem of this paper is to verify whether the continuous Galerkin method can maintain the local conservation of momentum and angular momentum at the discrete level when using the nonlinear terms in the EMAC form. Specifically, the author proves this by introducing the so - called "diffusive volume form" of the weak conservation law and shows the effectiveness of this form in the Eulerian and Lagrangian frameworks.

### Main Contributions

1. **Theoretical Proof**: The author has proven that the continuous Galerkin method can accurately maintain the local conservation of momentum and angular momentum when using the EMAC form.
2. **Numerical Verification**: The theoretical results are verified through numerical experiments, showing the superiority of the EMAC form in long - time integration.
3. **Comparison with Other Forms**: Other common nonlinear forms (such as SKEW, CONV and ROT) are compared, and it is proven that they cannot maintain local conservation like EMAC.

In summary, this paper aims to provide theoretical and numerical evidence to prove that the continuous Galerkin method in the EMAC form can maintain the local conservation of momentum and angular momentum at the discrete level, thus providing a more accurate and stable solution for numerical simulation.