Skeleton-stabilized divergence-conforming B-spline discretizations for incompressible flow problems of high Reynolds number

Abstract

We consider a stabilization method for divergence-conforming B-spline discretizations of the incompressible Navier–Stokes problem wherein jumps in high-order normal derivatives of the velocity field are penalized across interior mesh facets. We prove that this method is pressure robust, consistent, and energy stable, and we show how to select the stabilization parameter appearing in the method so that excessive numerical dissipation is avoided in both the cross-wind direction and in the diffusion-dominated regime. We examine the efficacy of the method using a suite of numerical experiments, and we find the method yields optimal L2 and H1 convergence rates for the velocity field, eliminates spurious small-scale structures that pollute Galerkin approximations, and is effective as an Implicit Large Eddy Simulation (ILES) methodology.