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
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.
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
- 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
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.
Referee Report
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)
- [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.
- [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)
- [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.
- [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
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
-
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
-
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
Robustness of condensed preconditioner relies on authors' 2025 self-citation framework
specific steps
-
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
axioms (1)
- domain assumption The hybridizable discontinuous Galerkin discretization of the Stokes-Darcy system is uniformly well-posed with respect to physical parameters.
Reference graph
Works this paper leans on
-
[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
work page 1967
-
[2]
Daniele Boffi, Franco Brezzi, and Michel Fortin.Mixed finite element methods and applications, vol- ume 44. Springer, 2013
work page 2013
-
[3]
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
work page 2045
-
[4]
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]
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
work page 2022
-
[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
work page 2004
-
[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
work page 2005
-
[8]
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
work page 2015
-
[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
work page 2024
-
[10]
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
work page 2023
-
[11]
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
work page 2020
-
[12]
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
work page 2017
-
[13]
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
work page 2009
-
[14]
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
work page 2002
-
[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
work page 2021
-
[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
work page 2023
-
[17]
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
work page 2018
-
[18]
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
work page 1911
-
[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
work page 2017
-
[20]
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
work page 2052
-
[21]
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
work page 2025
-
[22]
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]
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
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[24]
Karl Erik Holter, Miroslav Kuchta, and Kent-Andre Mardal. Robust preconditioning of monolithically coupled multiphysics problems.arXiv preprint arXiv:2001.05527, 2020
-
[25]
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
work page 2021
-
[26]
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]
Iury Igreja. A scalar/vector stabilized hybrid mixed method for Darcy–Stokes flows.Journal of Com- putational and Applied Mathematics, page 116879, 2025
work page 2025
-
[28]
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
work page 2018
-
[29]
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
work page 2010
-
[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
work page 2002
-
[31]
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
work page 2018
-
[32]
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
work page 2011
-
[33]
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
work page 2015
-
[34]
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
work page 1982
-
[35]
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
work page 1936
-
[36]
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
work page 1959
-
[37]
Philip Geoffrey Saffman. On the boundary condition at the surface of a porous medium.Studies in applied mathematics, 50(2):93–101, 1971
work page 1971
-
[38]
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
work page 1997
-
[39]
Joachim Sch¨ oberl. C++ 11 implementation of finite elements in NGSolve.Institute for analysis and scientific computing, Vienna University of Technology, 30, 2014
work page 2014
-
[40]
H. Van Emden and U. M. Yang. BoomerAMG: A parallel algebraic multigrid solver and preconditioner. Applied Numerical Mathematics, 41(1):155–177, 2002
work page 2002
-
[41]
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...
work page 2011
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.