A Helmholtz-Leray projection method with variational multiscale stabilization for the Navier-Stokes equations

Abstract

The Galerkin finite element formulation of the incompressible Navier-Stokes equations presents two principal challenges: maintaining stable velocity-pressure coupling and controlling instability in advection-dominated regimes. Moreover, the monolithic formulation produces a coupled nonlinear saddle-point problem. In this work, we present a residual-based variational multiscale (VMS) stabilization of an incremental Helmholtz-Leray projection method that replaces this coupled saddle-point problem with a nonlinear velocity predictor, a pressure Poisson equation, and a velocity projection. The multiscale decomposition is applied only to the predicted velocity; neither the pressure nor the corrected, weakly divergence-free velocity is decomposed into coarse and fine scales. The modeled velocity fine scale contributes consistently to all three subproblems, introducing SUPG-like stabilization in the momentum predictor and a PSPG-like residual contribution in the pressure Poisson equation. We provide a formal error decomposition that separates the BDF2 time-discretization, projection-splitting, and spatial-VMS errors and, under stated stability and spatial-approximation assumptions, yields a combined velocity error estimate with second-order temporal accuracy. Numerical results for manufactured solutions, lid-driven cavity flow, flow past a cylinder, and the Taylor-Green vortex agree closely with established reference data. Comparisons with monolithic VMS indicate that omitting the pressure fine scale reduces drag overprediction and excess modeled dissipation, at the cost of increased divergence error. In the Taylor-Green tests, the projection formulation also reduces the average solution time per step by factors ranging from approximately $1.3\times$ to $2.7\times$ under identical solver settings.