A divergence-free velocity reconstruction for incompressible flows

Abstract

In incompressible flows with vanishing normal velocities at the boundary, irrotational forces in the momentum equations should be balanced completely by the pressure gradient. Unfortunately, nearly all available discretizations for incompressible flows violate this property. The origin of the problem is that discrete velocities are usually not divergence-free. Hence, the use of divergence-free velocity reconstructions is proposed wherever an L2 scalar product appears in the discrete variational formulation. The approach is illustrated and applied to a nonconforming MAC-like discretization for unstructured Delaunay grids. It is numerically demonstrated that a divergence-free velocity reconstruction based on the lowest-order Raviart–Thomas element increases the robustness and accuracy of an existing convergent discretization, when irrotational forces appear in the momentum equations.