pith. sign in

arxiv: 2604.22641 · v1 · submitted 2026-04-24 · 🧮 math.NA · cs.NA

Preconditioning of a hybridizable discontinuous Galerkin method for the coupled Stokes--Darcy system

Pith reviewed 2026-05-08 10:20 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords hybridizable discontinuous GalerkinStokes-Darcy systemparameter-robust preconditionerstatic condensationoperator preconditioninguniform well-posednessnumerical linear algebra
0
0 comments X

The pith

Parameter-robust preconditioners are proven for the statically condensed linear systems arising from hybridizable discontinuous Galerkin discretizations of the Stokes-Darcy system.

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

The paper shows how to build preconditioners that remain effective no matter how the physical parameters of the Stokes-Darcy problem vary. First the hybridizable discontinuous Galerkin scheme is shown to be uniformly well-posed with respect to those parameters, allowing an operator-preconditioning framework to produce a robust solver for the full discrete system. A second framework is then used to prove that the same preconditioner stays robust after static condensation reduces the system size. Readers care because the Stokes-Darcy model describes groundwater flow, filtration, and many other interface problems in which permeability and viscosity ratios can differ by many orders of magnitude and cause conventional solvers to fail.

Core claim

Proving uniform well-posedness of the hybridizable discontinuous Galerkin discretization of the coupled Stokes-Darcy system with respect to all physical parameters permits construction of a parameter-robust preconditioner for the non-condensed scheme via the operator-preconditioning framework. The resulting preconditioner is shown to remain robust on the statically condensed reduced system by applying the condensation-robustness framework. Numerical examples confirm that iteration counts stay bounded independently of mesh size and parameter values.

What carries the argument

Uniform well-posedness of the hybridizable discontinuous Galerkin scheme for the Stokes-Darcy coupling, which enables operator preconditioning whose robustness survives static condensation.

If this is right

  • The number of iterations required by the preconditioned solver remains bounded independently of mesh size and all physical parameters.
  • Large-scale three-dimensional simulations of coupled free-flow and porous-media problems become feasible without manual parameter tuning.
  • The same design strategy applies directly to other hybridizable discontinuous Galerkin formulations of interface problems once uniform well-posedness is established.
  • Static condensation can be performed without sacrificing the parameter robustness already achieved on the full system.

Where Pith is reading between the lines

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

  • The approach suggests that similar two-step preconditioning arguments could be used for other multiphysics couplings discretized with hybridizable methods.
  • If the uniform well-posedness result extends to time-dependent or nonlinear versions of the Stokes-Darcy system, the same preconditioners would immediately become available for those problems.
  • Implementation in existing finite-element libraries would allow direct testing on heterogeneous media with spatially varying parameters.

Load-bearing premise

The hybridizable discontinuous Galerkin discretization of the Stokes-Darcy system is uniformly well-posed independent of the physical parameters.

What would settle it

Numerical experiments in which the preconditioned iteration count grows without bound as the permeability-to-viscosity ratio tends to zero or infinity would disprove the claimed robustness.

Figures

Figures reproduced from arXiv: 2604.22641 by Esteban Henr\'iquez, Jeonghun J. Lee, Miroslav Kuchta, Sander Rhebergen.

Figure 1
Figure 1. Figure 1: Permeability, computed velocity field, and computed pressure field for the test case described in view at source ↗
read the original abstract

We propose parameter-robust preconditioners for the statically condensed linear system arising from a hybridizable discontinuous Galerkin discretization of the coupled Stokes--Darcy system. The design strategy relies on first applying the operator-preconditioning framework [Numer. Linear Algebra Appl., 18(1):1--40, 2011] to construct a preconditioner for the non-condensed discretization. This is done by proving uniform well-posedness of the scheme. Next, we prove robustness of the resulting condensed preconditioner applied to the reduced linear system using the framework we proposed in [SIAM J. Sci. Comput., 47(6):A3212--A3238, 2025]. Numerical examples demonstrate robustness of the proposed preconditioners.

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

2 major / 2 minor

Summary. The manuscript proposes parameter-robust preconditioners for the statically condensed linear system obtained from a hybridizable discontinuous Galerkin (HDG) discretization of the coupled Stokes-Darcy system. The strategy first establishes uniform well-posedness of the non-condensed HDG formulation to construct a preconditioner via the 2011 operator-preconditioning framework, then transfers robustness to the condensed system using the authors' 2025 SIAM J. Sci. Comput. framework, with numerical examples provided to illustrate the robustness.

Significance. If the uniform well-posedness holds with constants independent of viscosity, permeability, and interface parameters, the work offers a valuable contribution to robust solvers for multiphysics interface problems. The combination of HDG discretization, operator preconditioning, and static condensation extends the authors' prior framework in a concrete setting, and the inclusion of numerical validation strengthens the practical relevance for applications such as porous-media flow.

major comments (2)
  1. [Uniform well-posedness analysis] The uniform inf-sup and continuity estimates for the HDG scheme (particularly those involving the hybrid variables and the Stokes-Darcy transmission conditions) are load-bearing for both the initial preconditioner construction and the subsequent condensation argument. These estimates must remain independent of viscosity, permeability tensor anisotropy, and interface coupling parameters; any deterioration would invalidate the robustness claim. The manuscript should state the explicit parameter independence in the relevant theorem.
  2. [Condensed preconditioner analysis] When applying the 2025 SIAM J. Sci. Comput. framework to obtain robustness of the condensed preconditioner, the manuscript must confirm that static condensation does not introduce new parameter dependence (e.g., in the Schur complement associated with the hybrid variables). The transfer step should be verified explicitly for this HDG discretization.
minor comments (2)
  1. [Notation and parameters] Clarify the precise definition of the physical parameters (viscosity, permeability tensor, interface coefficients) and their ranges in both the analysis and the numerical experiments to avoid ambiguity.
  2. [Numerical examples] Ensure that all figures in the numerical section clearly label the parameter values tested and include reference to the corresponding theorem establishing robustness.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and constructive suggestions. We address the two major comments below and will revise the manuscript accordingly to improve clarity on the parameter independence.

read point-by-point responses
  1. Referee: [Uniform well-posedness analysis] The uniform inf-sup and continuity estimates for the HDG scheme (particularly those involving the hybrid variables and the Stokes-Darcy transmission conditions) are load-bearing for both the initial preconditioner construction and the subsequent condensation argument. These estimates must remain independent of viscosity, permeability tensor anisotropy, and interface coupling parameters; any deterioration would invalidate the robustness claim. The manuscript should state the explicit parameter independence in the relevant theorem.

    Authors: We agree that explicit statement of parameter independence strengthens the presentation. Theorem 3.1 establishes the uniform well-posedness of the HDG formulation, with the inf-sup and continuity constants independent of the viscosity, the permeability tensor (including anisotropy), and the interface coupling parameters. To address the referee's request, we will insert a short remark immediately after Theorem 3.1 that explicitly records this independence and references the relevant estimates in the proof. This change will be made in the revised version. revision: yes

  2. Referee: [Condensed preconditioner analysis] When applying the 2025 SIAM J. Sci. Comput. framework to obtain robustness of the condensed preconditioner, the manuscript must confirm that static condensation does not introduce new parameter dependence (e.g., in the Schur complement associated with the hybrid variables). The transfer step should be verified explicitly for this HDG discretization.

    Authors: We appreciate this point. Section 4 applies the 2025 framework directly to the statically condensed system, and the robustness follows from the uniform well-posedness of the non-condensed formulation together with the algebraic structure of the HDG hybridization. Nevertheless, to make the transfer explicit, we will add a brief paragraph in Section 4.2 that verifies the absence of new parameter dependence in the Schur complement for the hybrid variables, using the block structure of the HDG system and the fact that the condensation operator is parameter-independent. This verification will be included in the revision. revision: yes

Circularity Check

1 steps flagged

Robustness of condensed preconditioner relies on authors' 2025 self-citation framework

specific steps
  1. self citation load bearing [Abstract]
    "Next, we prove robustness of the resulting condensed preconditioner applied to the reduced linear system using the framework we proposed in [SIAM J. Sci. Comput., 47(6):A3212--A3238, 2025]."

    The load-bearing step establishing parameter-robustness after static condensation is justified solely by citation to the authors' prior work rather than a self-contained derivation or external theorem; the 2025 framework is invoked as the mechanism that allows the preconditioner to survive condensation.

full rationale

The paper proves uniform well-posedness of the HDG Stokes-Darcy scheme (new content) to apply the external 2011 operator-preconditioning framework, then transfers robustness to the statically condensed system exclusively via the authors' own 2025 SIAM J. Sci. Comput. framework. This constitutes moderate self-citation load-bearing on the central condensed-preconditioner claim, but the specific application, well-posedness estimates, and numerical validation remain independent. No self-definitional reductions, fitted predictions, or ansatz smuggling are present.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the uniform well-posedness of the HDG scheme (treated as a domain assumption) and on the applicability of the cited operator-preconditioning frameworks; no free parameters or new invented entities are introduced in the abstract.

axioms (1)
  • domain assumption The hybridizable discontinuous Galerkin discretization of the Stokes-Darcy system is uniformly well-posed with respect to physical parameters.
    Invoked to justify application of the operator-preconditioning framework to the non-condensed system.

pith-pipeline@v0.9.0 · 5437 in / 1344 out tokens · 23632 ms · 2026-05-08T10:20:29.539406+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

41 extracted references · 41 canonical work pages · 1 internal anchor

  1. [1]

    Boundary conditions at a naturally permeable wall.Journal of fluid mechanics, 30(1):197–207, 1967

    Gordon S Beavers and Daniel D Joseph. Boundary conditions at a naturally permeable wall.Journal of fluid mechanics, 30(1):197–207, 1967

  2. [2]

    Springer, 2013

    Daniele Boffi, Franco Brezzi, and Michel Fortin.Mixed finite element methods and applications, vol- ume 44. Springer, 2013

  3. [3]

    A parameter-robust iterative method for Stokes–Darcy problems retaining local mass conservation.ESAIM: Mathematical Modelling and Numerical Analysis, 54(6):2045–2067, 2020

    Wietse M Boon. A parameter-robust iterative method for Stokes–Darcy problems retaining local mass conservation.ESAIM: Mathematical Modelling and Numerical Analysis, 54(6):2045–2067, 2020

  4. [4]

    Parameter-robust preconditioners for the Stokes-Darcy coupled problem without fractional operators.arXiv preprint arXiv:2501.06369, 2025

    Wietse M Boon, Xiaozhe Hu, and Xue Wang. Parameter-robust preconditioners for the Stokes-Darcy coupled problem without fractional operators.arXiv preprint arXiv:2501.06369, 2025. 17

  5. [5]

    Robust monolithic solvers for the Stokes–Darcy problem with the Darcy equation in primal form.SIAM Journal on Scientific Computing, 44(4):B1148–B1174, 2022

    Wietse M Boon, Timo Koch, Miroslav Kuchta, and Kent-Andre Mardal. Robust monolithic solvers for the Stokes–Darcy problem with the Darcy equation in primal form.SIAM Journal on Scientific Computing, 44(4):B1148–B1174, 2022

  6. [6]

    Korn’s inequalities for piecewiseH 1 vector fields.Mathematics of Computation, pages 1067–1087, 2004

    Susanne C Brenner. Korn’s inequalities for piecewiseH 1 vector fields.Mathematics of Computation, pages 1067–1087, 2004

  7. [7]

    Stabilized Crouzeix-Raviart element for the Darcy-Stokes problem

    Erik Burman and Peter Hansbo. Stabilized Crouzeix-Raviart element for the Darcy-Stokes problem. Numerical Methods for Partial Differential Equations: An International Journal, 21(5):986–997, 2005

  8. [8]

    New fully-mixed finite element methods for the Stokes–Darcy coupling.Computer Methods in Applied Me- chanics and Engineering, 295:362–395, 2015

    Jessika Camano, Gabriel N Gatica, Ricardo Oyarz´ ua, Ricardo Ruiz-Baier, and Pablo Venegas. New fully-mixed finite element methods for the Stokes–Darcy coupling.Computer Methods in Applied Me- chanics and Engineering, 295:362–395, 2015

  9. [9]

    Aycil Cesmelioglu, Jeonghun J Lee, and Sander Rhebergen. A strongly conservative hybridizable discon- tinuous Galerkin method for the coupled time-dependent Navier–Stokes and Darcy problem.ESAIM: Mathematical Modelling and Numerical Analysis, 58(1):273–302, 2024

  10. [10]

    A hybridizable discontinuous Galerkin method for the coupled Navier–Stokes and Darcy problem.Journal of Computational and Applied Mathematics, 422:114923, 2023

    Aycil Cesmelioglu and Sander Rhebergen. A hybridizable discontinuous Galerkin method for the coupled Navier–Stokes and Darcy problem.Journal of Computational and Applied Mathematics, 422:114923, 2023

  11. [11]

    An embedded–hybridized discontinuous Galerkin method for the coupled Stokes–Darcy system.Journal of Computational and Applied Mathe- matics, 367:112476, 2020

    Aycil Cesmelioglu, Sander Rhebergen, and Garth N Wells. An embedded–hybridized discontinuous Galerkin method for the coupled Stokes–Darcy system.Journal of Computational and Applied Mathe- matics, 367:112476, 2020

  12. [12]

    Convergence of IPDG for coupled time-dependent Navier–Stokes and Darcy equations.Journal of Computational and Applied Mathematics, 324:25–48, 2017

    Nabil Chaabane, Vivette Girault, Charles Puelz, and Beatrice Riviere. Convergence of IPDG for coupled time-dependent Navier–Stokes and Darcy equations.Journal of Computational and Applied Mathematics, 324:25–48, 2017

  13. [13]

    A hybridizable discontinuous Galerkin method for steady-state convection-diffusion-reaction problems.SIAM Journal on Scientific Computing, 31(5):3827–3846, 2009

    Bernardo Cockburn, Bo Dong, Johnny Guzm´ an, Marco Restelli, and Riccardo Sacco. A hybridizable discontinuous Galerkin method for steady-state convection-diffusion-reaction problems.SIAM Journal on Scientific Computing, 31(5):3827–3846, 2009

  14. [14]

    Mathematical and numerical models for coupling surface and groundwater flows.Applied Numerical Mathematics, 43(1-2):57–74, 2002

    Marco Discacciati, Edie Miglio, and Alfio Quarteroni. Mathematical and numerical models for coupling surface and groundwater flows.Applied Numerical Mathematics, 43(1-2):57–74, 2002

  15. [15]

    Guosheng Fu. Uniform auxiliary space preconditioning for HDG methods for elliptic operators with a parameter dependent low order term.SIAM Journal on Scientific Computing, 43(6):A3912–A3937, 2021

  16. [16]

    Guosheng Fu and Wenzheng Kuang. Uniform block-diagonal preconditioners for divergence-conforming HDG methods for the generalized stokes equations and the linear elasticity equations.IMA Journal of Numerical Analysis, 43(3):1718–1741, 2023

  17. [17]

    A strongly conservative hybrid DG/mixed FEM for the coupling of Stokes and Darcy flow.Journal of Scientific Computing, 77(3):1605–1620, 2018

    Guosheng Fu and Christoph Lehrenfeld. A strongly conservative hybrid DG/mixed FEM for the coupling of Stokes and Darcy flow.Journal of Scientific Computing, 77(3):1605–1620, 2018

  18. [18]

    Analysis of fully-mixed finite element methods for the Stokes-Darcy coupled problem.Mathematics of Computation, 80(276):1911–1948, 2011

    Gabriel Gatica, Ricardo Oyarz´ ua, and Francisco-Javier Sayas. Analysis of fully-mixed finite element methods for the Stokes-Darcy coupled problem.Mathematics of Computation, 80(276):1911–1948, 2011

  19. [19]

    Analysis of the HDG method for the Stokes–Darcy coupling

    Gabriel N Gatica and Fil´ ander A Sequeira. Analysis of the HDG method for the Stokes–Darcy coupling. Numerical Methods for Partial Differential Equations, 33(3):885–917, 2017

  20. [20]

    DG approximation of coupled Navier–Stokes and Darcy equations by Beaver–Joseph–Saffman interface condition.SIAM Journal on Numerical Analysis, 47(3):2052–2089, 2009

    Vivette Girault and B´ eatrice Rivi` ere. DG approximation of coupled Navier–Stokes and Darcy equations by Beaver–Joseph–Saffman interface condition.SIAM Journal on Numerical Analysis, 47(3):2052–2089, 2009. 18

  21. [21]

    Parameter-robust preconditioning for hybridizable symmetric discretizations.SIAM Journal on Scientific Computing, 47(6):A3212–A3238, 2025

    Esteban Henr´ ıquez, Jeonghun J Lee, and Sander Rhebergen. Parameter-robust preconditioning for hybridizable symmetric discretizations.SIAM Journal on Scientific Computing, 47(6):A3212–A3238, 2025

  22. [22]

    Preconditioning of a hybridizable discon- tinuous Galerkin method for Biot’s consolidation model.arXiv preprint arXiv:2508.12991, 2025

    Esteban Henr´ ıquez, Jeonghun J Lee, and Sander Rhebergen. Preconditioning of a hybridizable discon- tinuous Galerkin method for Biot’s consolidation model.arXiv preprint arXiv:2508.12991, 2025

  23. [23]

    Robust preconditioning for an HDG discretization of the time-dependent Stokes equations

    Esteban Henr´ ıquez, Jeonghun J Lee, and Sander Rhebergen. Robust preconditioning for an HDG discretization of the time-dependent Stokes equations.arXiv preprint arXiv:2604.07202, 2026

  24. [24]

    Robust preconditioning of monolithically coupled multiphysics problems.arXiv preprint arXiv:2001.05527, 2020

    Karl Erik Holter, Miroslav Kuchta, and Kent-Andre Mardal. Robust preconditioning of monolithically coupled multiphysics problems.arXiv preprint arXiv:2001.05527, 2020

  25. [25]

    Robust preconditioning for coupled Stokes– Darcy problems with the Darcy problem in primal form.Computers & Mathematics with Applications, 91:53–66, 2021

    Karl Erik Holter, Miroslav Kuchta, and Kent-Andre Mardal. Robust preconditioning for coupled Stokes– Darcy problems with the Darcy problem in primal form.Computers & Mathematics with Applications, 91:53–66, 2021

  26. [26]

    Parameter-robust preconditioner for Stokes-Darcy coupled problem with Lagrange multiplier.arXiv preprint arXiv:2503.19261, 2025

    Xiaozhe Hu, Miroslav Kuchta, Kent-Andre Mardal, and Xue Wang. Parameter-robust preconditioner for Stokes-Darcy coupled problem with Lagrange multiplier.arXiv preprint arXiv:2503.19261, 2025

  27. [27]

    A scalar/vector stabilized hybrid mixed method for Darcy–Stokes flows.Journal of Com- putational and Applied Mathematics, page 116879, 2025

    Iury Igreja. A scalar/vector stabilized hybrid mixed method for Darcy–Stokes flows.Journal of Com- putational and Applied Mathematics, page 116879, 2025

  28. [28]

    A stabilized hybrid mixed DGFEM naturally coupling Stokes–Darcy flows.Computer Methods in Applied Mechanics and Engineering, 339:739–768, 2018

    Iury Igreja and Abimael FD Loula. A stabilized hybrid mixed DGFEM naturally coupling Stokes–Darcy flows.Computer Methods in Applied Mechanics and Engineering, 339:739–768, 2018

  29. [29]

    A strongly conservative finite element method for the coupling of Stokes and Darcy flow.Journal of Computational Physics, 229(17):5933–5943, 2010

    Guido Kanschat and B´ eatrice Riviere. A strongly conservative finite element method for the coupling of Stokes and Darcy flow.Journal of Computational Physics, 229(17):5933–5943, 2010

  30. [30]

    Coupling fluid flow with porous media flow

    William J Layton, Friedhelm Schieweck, and Ivan Yotov. Coupling fluid flow with porous media flow. SIAM Journal on Numerical Analysis, 40(6):2195–2218, 2002

  31. [31]

    Quasi-optimality of a pressure-robust nonconforming finite element method for the Stokes-problem.Mathematics of Compu- tation, 87(312):1543–1566, 2018

    Alexander Linke, Christian Merdon, Michael Neilan, and Felix Neumann. Quasi-optimality of a pressure-robust nonconforming finite element method for the Stokes-problem.Mathematics of Compu- tation, 87(312):1543–1566, 2018

  32. [32]

    Preconditioning discretizations of systems of partial differ- ential equations.Numerical Linear Algebra with Applications, 18(1):1–40, 2011

    Kent-Andre Mardal and Ragnar Winther. Preconditioning discretizations of systems of partial differ- ential equations.Numerical Linear Algebra with Applications, 18(1):1–40, 2011

  33. [33]

    Strong coupling of finite element meth- ods for the Stokes–Darcy problem.IMA Journal of Numerical Analysis, 35(2):969–988, 2015

    Antonio M´ arquez, Salim Meddahi, and Francisco-Javier Sayas. Strong coupling of finite element meth- ods for the Stokes–Darcy problem.IMA Journal of Numerical Analysis, 35(2):969–988, 2015

  34. [34]

    Analysis of a hybridized/interface stabilized finite element method for the Stokes equations.SIAM Journal on Numerical Analysis, 55(4):1982–2003, 2017

    Sander Rhebergen and Garth N Wells. Analysis of a hybridized/interface stabilized finite element method for the Stokes equations.SIAM Journal on Numerical Analysis, 55(4):1982–2003, 2017

  35. [35]

    Preconditioning of a hybridized discontinuous Galerkin finite element method for the Stokes equations.Journal of Scientific Computing, 77(3):1936–1952, 2018

    Sander Rhebergen and Garth N Wells. Preconditioning of a hybridized discontinuous Galerkin finite element method for the Stokes equations.Journal of Scientific Computing, 77(3):1936–1952, 2018

  36. [36]

    Locally conservative coupling of Stokes and Darcy flows.SIAM Journal on Numerical Analysis, 42(5):1959–1977, 2005

    B´ eatrice Rivi` ere and Ivan Yotov. Locally conservative coupling of Stokes and Darcy flows.SIAM Journal on Numerical Analysis, 42(5):1959–1977, 2005

  37. [37]

    On the boundary condition at the surface of a porous medium.Studies in applied mathematics, 50(2):93–101, 1971

    Philip Geoffrey Saffman. On the boundary condition at the surface of a porous medium.Studies in applied mathematics, 50(2):93–101, 1971

  38. [38]

    NETGEN an advancing front 2D/3D-mesh generator based on abstract rules.Com- puting and visualization in science, 1(1):41–52, 1997

    Joachim Sch¨ oberl. NETGEN an advancing front 2D/3D-mesh generator based on abstract rules.Com- puting and visualization in science, 1(1):41–52, 1997. 19

  39. [39]

    C++ 11 implementation of finite elements in NGSolve.Institute for analysis and scientific computing, Vienna University of Technology, 30, 2014

    Joachim Sch¨ oberl. C++ 11 implementation of finite elements in NGSolve.Institute for analysis and scientific computing, Vienna University of Technology, 30, 2014

  40. [40]

    Van Emden and U

    H. Van Emden and U. M. Yang. BoomerAMG: A parallel algebraic multigrid solver and preconditioner. Applied Numerical Mathematics, 41(1):155–177, 2002

  41. [41]

    Analysis of an interface stabilized finite element method: the advection-diffusion- reaction equation.SIAM journal on numerical analysis, 49(1):87–109, 2011

    Garth N Wells. Analysis of an interface stabilized finite element method: the advection-diffusion- reaction equation.SIAM journal on numerical analysis, 49(1):87–109, 2011. A Proof of inf-sup condition eq. (15) We will prove that for anyq d h ∈Q d h there existsv h ∈V d h such that|||v h|||v,d ≤c −1 3 |||qd h|||q,d andb d h(qd h, vh) = |||qd h|||2 q,d. Fi...