Finite element method for the Navier-Stokes equation with local conservation properties

If you’re solving the incompressible Navier-Stokes equations, I’d recommend looking into divergence conforming Discontinuous Galerkin (DG) methods. These methods are non-conforming and are similar in some ways to finite volume methods. Upwinding can be incorporated naturally through the specification of the numerical flux. Provided the finite element spaces are chosen such that the divergence operator maps the velocity space into the pressure space, mass is conserved exactly. Exact mass conservation is advantageous for a number of reasons, one of which being that the error estimates for the velocity field are independent of the pressure field, and thus the pressure approximation cannot pollute the velocity approximation.

I am planning on adding a FEniCSx demo of this type of method soon. In the meantime, you may find these references helpful:

8 Likes