Stream function -- pressure virtual element methods for the Stokes--Darcy interface problem

Franco Dassi; Rekha Khot; David Mora; Andres E. Rubiano; Ricardo Ruiz-Baier

Stream function -- pressure virtual element methods for the Stokes--Darcy interface problem

Numerical Analysis 2026-07-02 v1

Authors: Franco Dassi , Rekha Khot , David Mora , Andres E. Rubiano , Ricardo Ruiz-Baier

Abstract

This paper introduces a novel Virtual Element Method (VEM) for the coupled Stokes--Darcy system in primal-primal form. In the free-flow Stokes domain, we implement a stream function formulation that inherently satisfies the incompressibility constraint and reduces computational cost. Across the interface, mass conservation, normal stress balance, and the Beavers--Joseph--Saffman slip condition are enforced to couple the biharmonic stream function equation with the Darcy's pressure equation. Leveraging VEM's ability to handle general polygonal meshes, the proposed method naturally accommodates irregular interface geometries without requiring remeshing or adaptive refinement. The accuracy of the method is validated through several numerical simulations that include applications to dead-end filtration, and network flow in bioartificial organs.

Cite

@article{arxiv.2607.01622,
  title  = {Stream function -- pressure virtual element methods for the Stokes--Darcy interface problem},
  author = {Franco Dassi and Rekha Khot and David Mora and Andres E. Rubiano and Ricardo Ruiz-Baier},
  journal= {arXiv preprint arXiv:2607.01622},
  year   = {2026}
}