pith. sign in

arxiv: 2604.15509 · v1 · submitted 2026-04-16 · 🧮 math.NA · cs.NA

A Correction Function-based KFBI Method for Brinkman Interface Problems

Han Zhou , Wenjun Ying This is my paper

Pith reviewed 2026-05-10 09:46 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords interface problemsboundary integral methodcorrection functioncollocation methodStokes equationsBrinkman equationsmoving interfacesnumerical discretization
0
0 comments X

The pith

A correction-function-based kernel-free boundary integral method solves Stokes and Brinkman interface problems by recasting them as boundary integrals and using collocation on local jump problems.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper develops a method for fluid interface problems in which coefficients jump across a boundary by converting the original discontinuous-coefficient equations into equivalent boundary integral equations. These integrals involve potential functions that satisfy simpler problems without discontinuities inside the domain. The integrals are discretized with a corrected Marker-and-Cell scheme, while jumps in the solution are captured by a local correction function that satisfies a Cauchy problem near the interface. The Cauchy problem is solved by a collocation method for which the authors give criteria on the minimal number of points and prove solvability. Numerical tests on both fixed and moving interfaces confirm that the resulting scheme is accurate and efficient.

Core claim

We propose a correction-function-based kernel-free boundary integral (CF-KFBI) method for solving Stokes- and Brinkman-type interface problems. We begin by recasting the original interface problem with discontinuous coefficients as boundary integral equations, in which the integral operators can be interpreted as boundary data for potential functions that satisfy simpler interface problems without coefficient discontinuities. Each such interface problem is discretized using a corrected Marker-and-Cell (MAC) scheme. Within a narrow band around the interface, we introduce a local correction function that represents the solution jump, leading to a local Cauchy problem. This problem is solved by

What carries the argument

The local correction function that encodes the solution jump across the interface and is obtained by solving a local Cauchy problem with a collocation method whose solvability is proven for a minimal set of points.

If this is right

  • The method produces accurate numerical solutions for Stokes and Brinkman interface problems while avoiding direct treatment of discontinuous coefficients inside the domain.
  • The approach extends without change to both fixed and moving interfaces.
  • Explicit criteria for the minimal number of collocation points guarantee solvability of the local correction problem under the stated smoothness assumptions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same local-correction structure could be applied to other elliptic interface problems whose jumps can be isolated in a narrow band.
  • Replacing the current MAC discretization with higher-order schemes might raise the global convergence rate while keeping the collocation step unchanged.
  • Automatic selection of collocation points based on local geometry could make the method more robust for highly curved or deforming interfaces.

Load-bearing premise

The local Cauchy problem for the correction function is well-posed and is solved to sufficient accuracy by the chosen collocation points, which requires the interface and solution to be sufficiently smooth.

What would settle it

A computation on an interface problem with deliberately reduced smoothness at the interface that produces instability or loss of accuracy in the collocation solve for the correction function.

Figures

Figures reproduced from arXiv: 2604.15509 by Han Zhou, Wenjun Ying.

Figure 1
Figure 1. Figure 1: Schematic of the computational domain B, the subdomains Ω±, and the interface Γ. Let B ⊂ R 2 be a rectangular domain and let Ω+ ⊂ B be a complex-shaped open subdomain embedded in B with a smooth boundary Γ = ∂Ω+, see [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Time evolution of the elastic interface and velocity field in the Brinkman flow. Starting from a highly [PITH_FULL_IMAGE:figures/full_fig_p021_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Time evolution of the enclosed area, area error, and GMRES iteration count for [PITH_FULL_IMAGE:figures/full_fig_p021_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Time evolution of the elastic interface and velocity field in unsteady Stokes flow. Starting from a highly [PITH_FULL_IMAGE:figures/full_fig_p022_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Time evolution of the enclosed area, area error, and GMRES iteration count for [PITH_FULL_IMAGE:figures/full_fig_p022_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Time evolution of the elastic interface in Brinkman shear flow for [PITH_FULL_IMAGE:figures/full_fig_p023_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Time evolution of the elastic interface in Brinkman shear flow for [PITH_FULL_IMAGE:figures/full_fig_p023_7.png] view at source ↗
read the original abstract

In this work, we propose a correction-function-based kernel-free boundary integral (CF-KFBI) method for solving Stokes- and Brinkman-type interface problems. We begin by recasting the original interface problem with discontinuous coefficients as boundary integral equations, in which the integral operators can be interpreted as boundary data for potential functions that satisfy simpler interface problems without coefficient discontinuities. Each such interface problem is discretized using a corrected Marker-and-Cell (MAC) scheme. Within a narrow band around the interface, we introduce a local correction function that represents the solution jump, leading to a local Cauchy problem. This problem is solved with a collocation method, for which we provide criteria for a minimal choice of collocation points and prove solvability. Several numerical experiments, including both fixed- and moving-interface problems, are presented to demonstrate the accuracy and efficiency of the proposed method.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The paper proposes a correction-function-based kernel-free boundary integral (CF-KFBI) method for Stokes- and Brinkman-type interface problems. It recasts the original discontinuous-coefficient problems as boundary integral equations whose kernels act as boundary data for simpler continuous-coefficient interface problems; these are discretized by a corrected Marker-and-Cell (MAC) scheme. Jumps are handled by introducing a local correction function that satisfies a Cauchy problem near the interface, which is solved by collocation; criteria for a minimal number of collocation points are given and solvability of the resulting linear system is proved under smoothness assumptions. Numerical experiments on both fixed- and moving-interface problems are presented to demonstrate accuracy and efficiency.

Significance. If the numerical accuracy holds and the local solvability result extends to a global convergence statement, the method would supply a practical, high-order approach for interface problems with discontinuous viscosity or permeability. The explicit criteria and proof for the collocation system constitute a concrete theoretical contribution, and the combination of kernel-free boundary integrals with a corrected MAC scheme and local correction functions is a natural extension of prior KFBI work. The absence of a global error analysis, however, currently confines the significance assessment to computational evidence.

major comments (1)
  1. [Theoretical analysis of the composite scheme] The manuscript proves solvability of the local collocation system for the correction function (see the section on the local Cauchy problem and collocation) but supplies no global a priori error estimate for the composite CF-KFBI scheme. Consequently the headline claims of accuracy rest entirely on the reported numerical experiments; the interaction of truncation errors arising from the boundary-integral representation, the corrected MAC discretization, and the collocation correction is not analyzed. This is load-bearing for the central claim of an accurate and efficient method.
minor comments (2)
  1. [Abstract] The abstract states that 'criteria for a minimal choice of collocation points' are provided; the precise statement of these criteria (e.g., dependence on interface curvature or solution smoothness) should be highlighted already in the abstract for clarity.
  2. [Section introducing the correction function] Notation for the correction function and the local Cauchy problem is introduced in the main text; a short dedicated subsection or table summarizing the minimal collocation-point conditions would improve readability.

Simulated Author's Rebuttal

1 responses · 1 unresolved

We are grateful to the referee for the thorough review and valuable feedback on our manuscript. We have carefully considered the major comment regarding the theoretical analysis and provide our response below.

read point-by-point responses

  1. Referee: The manuscript proves solvability of the local collocation system for the correction function (see the section on the local Cauchy problem and collocation) but supplies no global a priori error estimate for the composite CF-KFBI scheme. Consequently the headline claims of accuracy rest entirely on the reported numerical experiments; the interaction of truncation errors arising from the boundary-integral representation, the corrected MAC discretization, and the collocation correction is not analyzed. This is load-bearing for the central claim of an accurate and efficient method.

    Authors: We acknowledge that a global a priori error estimate for the composite scheme is not derived in the manuscript. The focus of the theoretical contribution is on establishing the solvability of the local collocation system for the correction function under the given smoothness assumptions. The accuracy of the overall method is supported by extensive numerical experiments on both fixed and moving interface problems, which show the expected convergence rates. Analyzing the global error would require a detailed study of how the truncation errors from the kernel-free boundary integral representation, the corrected MAC scheme, and the local collocation interact, particularly in the presence of the interface. Such an analysis is technically involved and lies beyond the scope of the current work, which prioritizes the development and validation of the CF-KFBI method. In the revised version, we have added a paragraph in the introduction and conclusion to explicitly state the scope of the theoretical results and to discuss the error sources heuristically based on the numerical observations. revision: partial

standing simulated objections not resolved
  • Deriving a global a priori error estimate for the full CF-KFBI scheme.

Circularity Check

0 steps flagged

Minor self-citation of prior KFBI/MAC work; central local collocation solvability proof is independent

full rationale

The derivation recasts the interface problem as boundary integrals, applies a corrected MAC scheme, and solves a local Cauchy problem for the correction function via collocation, with explicit criteria and a solvability proof under smoothness assumptions. This local proof and the composite scheme do not reduce by construction to fitted parameters or self-referential definitions. Prior KFBI and MAC references appear but are not load-bearing for the new correction-function component or the collocation analysis, which stands on its own equations and assumptions. No global error estimate is claimed as a derivation, so no circular reduction occurs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The approach rests on standard assumptions from numerical PDE theory and introduces a computational correction function; no free parameters or new physical entities are described.

axioms (1)
  • domain assumption The interface is sufficiently smooth and the solution sufficiently regular for the local Cauchy problem to be well-posed.
    Invoked when setting up the collocation method and proving solvability.
invented entities (1)
  • local correction function no independent evidence
    purpose: Represents the jump in the solution across the interface inside a narrow band.
    Introduced as a new computational device to convert the jump into a local solvable problem.

pith-pipeline@v0.9.0 · 5438 in / 1270 out tokens · 78424 ms · 2026-05-10T09:46:07.958060+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

54 extracted references · 54 canonical work pages

  1. [1]

    Durlofsky, J

    L. Durlofsky, J. Brady, Analysis of the brinkman equation as a model for flow in porous media, Physics of Fluids 30 (11) (1987) 3329

    work page 1987

  2. [2]

    Auriault, On the domain of validity of brinkman’s equation, Transport in porous media 79 (2) (2009) 215–223

    J.-L. Auriault, On the domain of validity of brinkman’s equation, Transport in porous media 79 (2) (2009) 215–223

    work page 2009

  3. [3]

    H. C. Brinkman, A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles, Applied Scientific Research 1 (1) (1949) 27–34

    work page 1949

  4. [4]

    D. A. Nield, A. Bejan, Convection in Porous Media, 4th Edition, Springer, Cham, 2017

    work page 2017

  5. [5]

    Pozrikidis, Boundary Integral and Singularity Methods for Linearized Viscous Flow, Cam- bridge University Press, Cambridge, 1992

    C. Pozrikidis, Boundary Integral and Singularity Methods for Linearized Viscous Flow, Cam- bridge University Press, Cambridge, 1992

    work page 1992

  6. [6]

    Barthès-Biesel, Motion and deformation of elastic capsules and vesicles in flow, Annual Review of Fluid Mechanics 48 (2016) 25–52

    D. Barthès-Biesel, Motion and deformation of elastic capsules and vesicles in flow, Annual Review of Fluid Mechanics 48 (2016) 25–52

    work page 2016

  7. [7]

    Lesinigo, C

    M. Lesinigo, C. D’Angelo, A. Quarteroni, A multiscale darcy–brinkman model for fluid flow in fractured porous media, Numerical Methods for Partial Differential Equations 27 (5) (2011) 112–136

    work page 2011

  8. [8]

    Kahshan, D

    M. Kahshan, D. Lu, N. H. Abu-Hamdeh, A. Golmohammadzadeh, M. Rahimi-Gorji, Darcy– brinkman flow of a viscous fluid through a porous duct: Application in blood filtration process, Journal of the Taiwan Institute of Chemical Engineers 117 (2020) 223–230

    work page 2020

  9. [9]

    Rohan, J

    E. Rohan, J. C. Turjanicová, V. Lukeš, Multiscale modelling and simulations of tissue perfusion using the biot–darcy–brinkman model, Computers & Structures 251 (2021) 106404. 24

    work page 2021

  10. [10]

    A.-R. A. Khaled, K. Vafai, The role of porous media in modeling flow and heat transfer in biological tissues, International Journal of Heat and Mass Transfer 46 (26) (2003) 4989–5003

    work page 2003

  11. [11]

    Babuška, The finite element method for elliptic equations with discontinuous coefficients, Computing 5 (3) (1970) 207–213

    I. Babuška, The finite element method for elliptic equations with discontinuous coefficients, Computing 5 (3) (1970) 207–213

    work page 1970

  12. [12]

    Guyomarc’h, C.-O

    G. Guyomarc’h, C.-O. Lee, K. Jeon, A discontinuous galerkin method for elliptic interface problems with application to electroporation, Communications in Numerical Methods in Engi- neering 25 (10) (2009) 991–1008

    work page 2009

  13. [13]

    J. H. Bramble, J. T. King, A finite element method for interface problems in domains with smooth boundaries and interfaces, Advances in Computational Mathematics 6 (1) (1996) 109– 138

    work page 1996

  14. [14]

    W. C. Wang, A jump condition capturing finite difference scheme for elliptic interface problems, SIAM Journal on Scientific Computing 25 (5) (2004) 1479–1496

    work page 2004

  15. [15]

    X. Xie, J. Xu, G. Xue, Uniformly-stable finite element methods for darcy-stokes-brinkman models, Journal of Computational Mathematics (2008) 437–455

    work page 2008

  16. [16]

    J. Li, J. M. Melenk, B. Wohlmuth, J. Zou, Optimal a priori estimates for higher order finite elements for elliptic interface problems, Applied Numerical Mathematics 60 (1–2) (2010) 19–37

    work page 2010

  17. [17]

    S. Groß, A. Reusken, Numerical Methods for Two-Phase Incompressible Flows, Vol. 40 of Springer Series in Computational Mathematics, Springer, 2011

    work page 2011

  18. [18]

    C. S. Peskin, The immersed boundary method, Acta Numerica 11 (2002) 479–517

    work page 2002

  19. [19]

    R. J. LeVeque, Z. Li, The immersed interface method for elliptic equations with discontinuous coefficients and singular sources, SIAM Journal on Numerical Analysis 31 (4) (1994) 1019–1044

    work page 1994

  20. [20]

    R. J. LeVeque, Z. Li, Immersed interface methods for Stokes flow with elastic boundaries or surface tension, SIAM Journal on Scientific Computing 18 (3) (1997) 709–735

    work page 1997

  21. [21]

    Li, An overview of the immersed interface method and its applications, Taiwanese Journal of Mathematics 7 (1) (2003) 1–49

    Z. Li, An overview of the immersed interface method and its applications, Taiwanese Journal of Mathematics 7 (1) (2003) 1–49

    work page 2003

  22. [22]

    Z. Li, K. Ito, The Immersed Interface Method: Numerical Solutions of PDEs Involving Inter- faces and Irregular Domains, Vol. 33 of Frontiers in Applied Mathematics, SIAM, 2006

    work page 2006

  23. [23]

    Y. Zhou, S. Zhao, M. Feig, G.-W. Wei, High order matched interface and boundary (MIB) methods for elliptic equations with discontinuous coefficients and singular sources, Journal of Computational Physics 213 (1) (2006) 1–30

    work page 2006

  24. [24]

    S. Yu, Y. Zhou, G.-W. Wei, Matched interface and boundary (MIB) method for elliptic prob- lems with sharp-edged interfaces, Journal of Computational Physics 224 (2) (2007) 729–756

    work page 2007

  25. [25]

    Xia, G.-W

    K. Xia, G.-W. Wei, MIB galerkin method for elliptic interface problems, Journal of Computa- tional and Applied Mathematics 272 (2014) 195–220

    work page 2014

  26. [26]

    A. N. Marques, J.-C. Nave, R. R. Rosales, A correction function method for Poisson problems with interface jump conditions, Journal of Computational Physics 230 (20) (2011) 7567–7597. 25

    work page 2011

  27. [27]

    R. P. Fedkiw, T. Aslam, B. Merriman, S. Osher, A non-oscillatory Eulerian approach to inter- faces in multimaterial flows (the ghost fluid method), Journal of Computational Physics 152 (2) (1999) 457–492

    work page 1999

  28. [28]

    Johansen, P

    H. Johansen, P. Colella, A Cartesian grid embedded boundary method for Poisson’s equation on irregular domains, Journal of Computational Physics 147 (1) (1998) 60–85

    work page 1998

  29. [29]

    Albright, Y

    J. Albright, Y. Epshteyn, M. Medvinsky, Q. Xia, High-order numerical schemes based on dif- ference potentials for 2D elliptic problems with material interfaces, Applied Numerical Math- ematics 116 (2017) 141–163

    work page 2017

  30. [30]

    N. Moes, J. Dolbow, T. Belytschko, A finite element method for crack growth without remesh- ing, International Journal for Numerical Methods in Engineering 46 (1) (1999) 131–150

    work page 1999

  31. [31]

    Belytschko, R

    T. Belytschko, R. Gracie, G. Ventura, A review of extended/generalized finite element methods for material modeling, Modelling and Simulation in Materials Science and Engineering 17 (4) (2009) 043001

    work page 2009

  32. [32]

    Z. Li, T. Lin, X. Wu, New Cartesian grid methods for interface problems using the finite element formulation, Numerische Mathematik 96 (1) (2003) 61–98

    work page 2003

  33. [33]

    Z. Li, T. Lin, R. C. Rogers, An immersed finite element space and its approximation capability, Numerical Methods for Partial Differential Equations 20 (3) (2004) 338–367

    work page 2004

  34. [34]

    Z. Li, X. Yang, An immersed finite element method for elasticity equations with interfaces, in: Recent Advances in Adaptive Computation, Vol. 383 of Contemporary Mathematics, American Mathematical Society, Providence, RI, 2005, pp. 285–298

    work page 2005

  35. [35]

    X. He, T. Lin, Y. Lin, X. Zhang, Immersed finite element methods for parabolic equations with moving interface, Numerical Methods for Partial Differential Equations 29 (2) (2013) 619–646

    work page 2013

  36. [36]

    Power, G

    H. Power, G. Miranda, Second kind integral equation formulation of Stokes flows past a particle of arbitrary shape, SIAM Journal on Applied Mathematics 47 (4) (1987) 689–698

    work page 1987

  37. [37]

    Greengard, M

    L. Greengard, M. C. Kropinski, An integral equation approach to the two-dimensional Stokes flow in doubly periodic domains, Journal of Computational Physics 197 (2) (2004) 344–368

    work page 2004

  38. [38]

    Pozrikidis, Boundary Integral and Singularity Methods for Linearized Viscous Flow, Vol

    C. Pozrikidis, Boundary Integral and Singularity Methods for Linearized Viscous Flow, Vol. 8 of Cambridge Texts in Applied Mathematics, Cambridge University Press, 1992

    work page 1992

  39. [39]

    Ahmadi, R

    E. Ahmadi, R. Cortez, H. Fujioka, Boundary integral formulation for flows containing an in- terface between two porous media, Journal of Fluid Mechanics 816 (2017) 71–93

    work page 2017

  40. [40]

    Tlupova, R

    S. Tlupova, R. Cortez, Boundary integral solutions of coupled stokes and darcy flows, Journal of Computational Physics 228 (2009) 158–179

    work page 2009

  41. [41]

    Boubendir, S

    Y. Boubendir, S. Tlupova, Domain decomposition methods for solving stokes–darcy problems with boundary integrals, SIAM Journal on Scientific Computing 35 (2013) B82–B106

    work page 2013

  42. [42]

    Tlupova, A domain decomposition solution of the stokes-darcy system in 3d based on bound- ary integrals, Journal of Computational Physics 450 (2022) 110824

    S. Tlupova, A domain decomposition solution of the stokes-darcy system in 3d based on bound- ary integrals, Journal of Computational Physics 450 (2022) 110824

    work page 2022

  43. [43]

    T. Y. Hou, J. S. Lowengrub, M. J. Shelley, Boundary integral methods for multicomponent fluids and multiphase materials, Journal of Computational Physics 169 (2) (2001) 302–362. 26

    work page 2001

  44. [44]

    L. Ying, G. Biros, D. Zorin, A high-order 3d boundary integral equation solver for elliptic PDEs in smooth domains, Journal of Computational Physics 219 (1) (2006) 247–275

    work page 2006

  45. [45]

    D. M. Ambrose, M. Siegel, S. Tlupova, A small-scale decomposition for 3d boundary integral computations with surface tension, Journal of Computational Physics 247 (2013) 168–191

    work page 2013

  46. [46]

    W. Ying, C. S. Henriquez, A kernel-free boundary integral method for elliptic boundary value problems, Journal of computational physics 227 (2) (2007) 1046–1074

    work page 2007

  47. [47]

    Y. Xie, W. Ying, A fourth-order kernel-free boundary integral method for implicitly defined surfaces in three space dimensions, Journal of Computational Physics 415 (2020) 109526

    work page 2020

  48. [48]

    H. Zhou, J. Yang, W. Ying, A kernel-free boundary integral method for the nonlinear poisson- boltzmann equation, Journal of Computational Physics 493 (2023) 112423

    work page 2023

  49. [49]

    H.Zhou, W.Ying, Acorrectionfunction–basedkernel-freeboundaryintegralmethodforelliptic pdes with implicitly defined interfaces, Journal of Computational Physics 496 (2024) 112545

    work page 2024

  50. [50]

    G. C. Hsiao, W. L. Wendland, Boundary Integral Equations, Vol. 164 of Applied Mathematical Sciences, Springer, New York, 2008

    work page 2008

  51. [51]

    McLean, Strongly Elliptic Systems and Boundary Integral Equations, Cambridge University Press, Cambridge, 2000

    W. McLean, Strongly Elliptic Systems and Boundary Integral Equations, Cambridge University Press, Cambridge, 2000

    work page 2000

  52. [52]

    Beale, A

    T. Beale, A. Layton, On the accuracy of finite difference methods for elliptic problems with interfaces, Communications in Applied Mathematics and Computational Science 1 (1) (2006) 91–119

    work page 2006

  53. [53]

    H. Dong, Z. Zhao, S. Li, W. Ying, J. Zhang, Second Order Convergence of a Modified MAC Scheme for Stokes Interface Problems, Journal of Scientific Computing 96 (1) (2023) 27

    work page 2023

  54. [54]

    Trottenberg, C

    U. Trottenberg, C. W. Oosterlee, A. Schüller, Multigrid, Academic Press, San Diego, 2001. 27

    work page 2001