A high-order, pressure-robust, and decoupled finite difference method for the Stokes problem

Qiwei Feng, Bin Han and Michael Neilan

Mathematics and Computers in Simulation (MATCOM), 2026, vol. 241, issue PB, 634-649

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 A1uh(1)=b1 and A2uh(2)=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 first approximate the pressure gradient locally, and then calculate the pressure itself locally with sixth-order accuracy, both 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, pressure gradient, and pressure 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 ν.

Keywords: Stokes problems; Pressure- and viscosity-robust FDMs; Sixth-order convergence rates; The decouple property; A novel biharmonic equation; Explicit formulas of FDMs (search for similar items in EconPapers)
Date: 2026
References: View complete reference list from CitEc
Citations:

Downloads: (external link)
http://www.sciencedirect.com/science/article/pii/S0378475425004628
Full text for ScienceDirect subscribers only

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:eee:matcom:v:241:y:2026:i:pb:p:634-649

DOI: 10.1016/j.matcom.2025.10.033

Access Statistics for this article

Mathematics and Computers in Simulation (MATCOM) is currently edited by Robert Beauwens

More articles in Mathematics and Computers in Simulation (MATCOM) from Elsevier
Bibliographic data for series maintained by Catherine Liu ().