A High-Order, Pressure-Robust, and Decoupled Finite Difference Method for the Stokes Problem
Abstract
In this paper, we consider the Stokes problem with Dirichlet boundary conditions and the constant kinematic viscosity in an axis-aligned domain . We decouple the velocity u and pressure p by deriving a novel biharmonic equation in and third-order boundary conditions on ∂. In contrast to the fourth-order streamfunction approach, our formulation does not require to be simply connected. For smooth velocity fields u in two dimensions, we explicitly construct a finite difference method (FDM) with sixth-order consistency to approximate u at all relevant grid points: interior points, boundary side points, and boundary corner points. The resulting scheme yields two linear systems A1u(1)h=b1 and A2u(2)h=b2, where A1,A2 are constant matrices, and b1,b2 are independent of the pressure p and the kinematic viscosity . Thus, the proposed method is pressure- and viscosity-robust. To accommodate velocity fields with less regularity, we modify the FDM by removing singular terms in the right-hand side vectors. Once the discrete velocity is computed, we apply a sixth-order finite difference operator to approximate the pressure gradient locally, without solving any additional linear systems. In our numerical experiments, we test both smooth and non-smooth solutions ( u,p) in a square domain, a triply connected domain, and an L-shaped domain in two dimensions. The results confirm sixth-order convergence of the velocity and pressure gradient in the ∞-norm for smooth solutions. For non-smooth velocity fields, our method achieves the expected lower-order convergence. Moreover, the observed velocity error \| uh- u\|∞ is independent of the pressure p and viscosity .
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.