Projection-based stabilization for high-order incompressible flow solvers

This work presents a novel stabilization strategy for the Galerkin formulation of the incompressible Navier-Stokes equations, developed to achieve high accuracy while ensuring convergence and compatibility with high-order elements on unstructured meshes. The numerical algorithm employs a fractional step method with carefully defined boundary conditions to obtain a consistent pressure field, enabling high-order temporal accuracy. The proposed stabilization is seamlessly integrated into the algorithm and shares the same underlying principle as the natural stabilization inherent in the fractional step method, both rely on the difference between the gradient operator and its projection. The numerical dissipation associated to the stabilization term is found to diminish with increasing polynomial order of the elements. Numerical test cases confirm the effectiveness of the method, demonstrating convergence under mesh refinement and increasing polynomial order.