pith. sign in

arxiv: 2605.23290 · v1 · pith:MGGODUVDnew · submitted 2026-05-22 · 🧮 math.NA · cs.NA

High-order, long-time stable and parallel decoupled GBDFk SAV ensemble schemes for the Navier--Stokes--Darcy flow with random hydraulic conductivity tensors

Wei-Wei Han , Fukeng Huang , Changxin Qiu This is my paper

Pith reviewed 2026-05-25 03:58 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Navier-Stokes-Darcy flowensemble schemesGSAVGBDFkhigh-order accuracylong-time stabilitydecoupled schemesrandom parameters
0
0 comments X

The pith

GSAV-GBDFk ensemble schemes achieve high-order temporal accuracy and long-time stability for Navier-Stokes-Darcy flows with random parameters.

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

The paper develops ensemble numerical schemes for the unsteady Navier-Stokes-Darcy system with uncertainties in initial conditions, forces, conductivity tensors, and interface conditions. By combining partitioned decoupling with the generalized scalar auxiliary variable approach and generalized BDFk time stepping, the methods achieve high-order accuracy while maintaining stability over long times without restricting the time step size. They also allow explicit handling of nonlinear terms and reuse the same coefficient matrix for multiple realizations, improving efficiency.

Core claim

The GSAV-GBDF$k$-Ensemble schemes integrate partitioned decoupling, the generalized scalar auxiliary variable approach, and generalized BDF$k$ discretizations for the unsteady Navier--Stokes--Darcy system with uncertain initial conditions, forcing terms, hydraulic conductivity tensors, and interface conditions. This yields high-order temporal accuracy, long-time stability, explicit nonlinear treatment, uniform-in-time bounds without step restrictions, optimal error estimates, and shared coefficient matrices across realizations.

What carries the argument

The GSAV-GBDF$k$-Ensemble scheme combining partitioned decoupling with generalized scalar auxiliary variable and generalized BDF$k$ discretizations.

If this is right

  • High-order temporal accuracy is attained.
  • Numerical solutions satisfy uniform-in-time bounds independent of time-step size.
  • Optimal-order error estimates hold.
  • The nonlinear term is treated explicitly while preserving stability.
  • A single coefficient matrix is shared across realizations.

Where Pith is reading between the lines

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

  • The shared-matrix structure may reduce computation time when realizations run in parallel.
  • The uniform bounds could guide adaptive time-step choices in long simulations of uncertain porous-media flow.

Load-bearing premise

The random hydraulic conductivity tensors and Lions-Beavers-Joseph-Saffman interface conditions allow the partitioned decoupling strategy combined with the GSAV approach to preserve uniform-in-time bounds and keep the shared coefficient matrix well-conditioned across realizations.

What would settle it

A computation showing that the observed temporal convergence rate falls below the expected order k for a chosen random conductivity realization, or that solution norms grow unbounded for large fixed time steps, would falsify the claims.

Figures

Figures reproduced from arXiv: 2605.23290 by Changxin Qiu, Fukeng Huang, Wei-Wei Han.

Figure 1
Figure 1. Figure 1: Domain Ωp for the porous media and Ωf for the free flow. A = (0, 1), B = (0, 3/4), C = (1/2, 1/4), D = (1/2, 0), E = (3/4, 0), F = (3/4, 1/4), G = (1, 1/4), H = (1, 1/2), I = (3/4, 1/2) and J = (1/4, 1). 21 [PITH_FULL_IMAGE:figures/full_fig_p021_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Mean velocity from GSAV-GBDF3-Ensemble algorithm for [PITH_FULL_IMAGE:figures/full_fig_p022_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Mean velocity from tradition algorithm for [PITH_FULL_IMAGE:figures/full_fig_p022_3.png] view at source ↗
read the original abstract

We develop and analyze high-order ensemble schemes for the unsteady Navier--Stokes--Darcy system with uncertain initial conditions, forcing terms, hydraulic conductivity tensors, and Lions-Beavers-Joseph-Saffman interface conditions. The proposed schemes which are called GSAV-GBDF$k$-Ensemble schemes integrate a partitioned decoupling strategy, the generalized scalar auxiliary variable (GSAV) approach, and generalized BDF$k$ discretizations. This framework achieves high-order temporal accuracy and long-time stability, permits explicit treatment of the nonlinear term, and facilitates an efficient ensemble implementation for multiple parameter realizations by sharing a single, unified coefficient matrix at each time step. Moreover, the numerical solutions are shown to satisfy uniform-in-time bounds without time-step restrictions. Owing to the ensemble formulation, the resulting linear systems share common coefficient matrices, which significantly improves computational efficiency. We further establish optimal-order error estimates for the proposed high-order schemes. Numerical results are included to confirm the theoretical analysis and to illustrate the accuracy, stability, and efficiency of the proposed methods.

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 develops GSAV-GBDF$k$-Ensemble schemes for the unsteady Navier-Stokes-Darcy system with random initial conditions, forcing, hydraulic conductivity tensors, and Lions-Beavers-Joseph-Saffman interface conditions. The schemes combine partitioned decoupling, the generalized scalar auxiliary variable (GSAV) approach, and generalized BDF$k$ time discretizations to achieve high-order temporal accuracy, long-time stability without time-step restrictions, explicit treatment of the nonlinear term, uniform-in-time bounds, shared coefficient matrices across realizations for efficiency, and optimal-order error estimates, with numerical results confirming the analysis.

Significance. If the stability and error claims hold, the work would deliver a computationally efficient framework for ensemble simulations of coupled free/porous-media flows with uncertainty in conductivity tensors. The combination of GSAV stabilization with ensemble sharing of a single matrix at each step, together with unconditional uniform bounds, addresses a practical need in uncertainty quantification for Navier-Stokes-Darcy problems; the provision of optimal error estimates and high-order accuracy strengthens the contribution relative to existing first-order or conditionally stable methods.

major comments (2)
  1. [§4.2, Theorem 4.1] §4.2, Theorem 4.1: the uniform-in-time bound is stated to hold without time-step restrictions for arbitrary random conductivity tensors satisfying the given positivity and boundedness assumptions, but the proof sketch relies on the GSAV auxiliary variable absorbing the interface terms; it is unclear whether the constant in the bound remains independent of the realization when the conductivity tensor varies strongly across the ensemble (the abstract claims this independence, yet the estimate appears to retain a dependence on the essential supremum of the tensor eigenvalues).
  2. [§5.1, Eq. (5.3)] §5.1, Eq. (5.3): the error estimate for the velocity and pressure is claimed to be optimal order in time and space, but the analysis of the decoupled interface terms under Lions-Beavers-Joseph-Saffman conditions introduces an additional consistency error whose order is not explicitly tracked; if this term is only first-order, it would cap the overall temporal accuracy below the advertised BDF$k$ order for $k>1$.
minor comments (2)
  1. The notation for the ensemble index set and the shared matrix is introduced without a dedicated table or list of symbols; adding one would improve readability.
  2. Figure 6.3 caption does not specify the mesh size or the number of realizations used; this information should be added for reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading, positive assessment of the contribution, and constructive comments. We address each major point below with clarifications based on the analysis in the manuscript.

read point-by-point responses

  1. Referee: [§4.2, Theorem 4.1] §4.2, Theorem 4.1: the uniform-in-time bound is stated to hold without time-step restrictions for arbitrary random conductivity tensors satisfying the given positivity and boundedness assumptions, but the proof sketch relies on the GSAV auxiliary variable absorbing the interface terms; it is unclear whether the constant in the bound remains independent of the realization when the conductivity tensor varies strongly across the ensemble (the abstract claims this independence, yet the estimate appears to retain a dependence on the essential supremum of the tensor eigenvalues).

    Authors: The uniform-in-time bound of Theorem 4.1 is independent of individual realizations. Under the problem assumptions, all conductivity tensors satisfy the same uniform positivity and boundedness constants (i.e., there exist λ_min > 0 and λ_max < ∞ independent of the realization such that λ_min |ξ|^2 ≤ K(x,ω)ξ·ξ ≤ λ_max |ξ|^2 a.e.). The GSAV auxiliary variable absorbs the interface contributions, and the resulting discrete energy estimate depends only on these uniform constants (together with data norms that are uniform across the ensemble). The essential supremum is therefore taken over the entire ensemble, not per realization. We will add an explicit remark after the theorem statement clarifying this uniformity. revision: partial

  2. Referee: [§5.1, Eq. (5.3)] §5.1, Eq. (5.3): the error estimate for the velocity and pressure is claimed to be optimal order in time and space, but the analysis of the decoupled interface terms under Lions-Beavers-Joseph-Saffman conditions introduces an additional consistency error whose order is not explicitly tracked; if this term is only first-order, it would cap the overall temporal accuracy below the advertised BDF$k$ order for $k>1$.

    Authors: The consistency error arising from the decoupled treatment of the Lions-Beavers-Joseph-Saffman interface conditions is estimated in the proof of the error bound (5.3) and is shown to be O(Δt^k) by using the same BDFk extrapolation and GSAV stabilization that preserve the temporal order. The interface terms are handled via the partitioned scheme and absorbed into the higher-order truncation error without order reduction. We agree that the tracking of this term can be made more transparent and will insert a short auxiliary lemma (or expanded paragraph) that isolates the interface consistency error and verifies its order explicitly. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper develops and analyzes GSAV-GBDFk ensemble schemes for the Navier-Stokes-Darcy system. Claims of high-order accuracy, uniform-in-time bounds, and optimal error estimates rest on direct analysis of the partitioned decoupling, GSAV stabilization, and generalized BDFk discretizations applied to the random-tensor problem with Lions-Beavers-Joseph-Saffman conditions. No steps reduce by construction to fitted parameters, self-definitional loops, or load-bearing self-citations whose validity depends on the present work. The framework is self-contained against external benchmarks of numerical PDE analysis.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Only the abstract is available; no explicit free parameters, ad-hoc axioms, or invented entities are described beyond standard domain assumptions on the flow system.

axioms (1)
  • domain assumption The Lions-Beavers-Joseph-Saffman interface conditions hold for the coupled system.
    Stated in the abstract as part of the problem setup.

pith-pipeline@v0.9.0 · 5727 in / 1263 out tokens · 24994 ms · 2026-05-25T03:58:24.288194+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

62 extracted references · 62 canonical work pages

  1. [1]

    SIAM Journal on Numerical Analysis , volume=

    DG approximation of coupled Navier--Stokes and Darcy equations by Beaver--Joseph--Saffman interface condition , author=. SIAM Journal on Numerical Analysis , volume=. 2009 , publisher=

    work page 2009

  2. [2]

    Computer Methods in Applied Mechanics and Engineering , volume=

    On the solution of the coupled Navier--Stokes and Darcy equations , author=. Computer Methods in Applied Mechanics and Engineering , volume=. 2009 , publisher=

    work page 2009

  3. [3]

    Studies in applied mathematics , volume=

    On the boundary condition at the surface of a porous medium , author=. Studies in applied mathematics , volume=. 1971 , publisher=

    work page 1971

  4. [4]

    Journal of fluid mechanics , volume=

    Boundary conditions at a naturally permeable wall , author=. Journal of fluid mechanics , volume=. 1967 , publisher=

    work page 1967

  5. [5]

    Long time stability and convergence for fully discrete nonlinear

    Shen, Jie , journal=. Long time stability and convergence for fully discrete nonlinear. 1990 , publisher=

    work page 1990

  6. [6]

    Stability and error analysis of a new class of higher-order consistent splitting schemes for the

    Huang, Fukeng and Shen, Jie , journal=. Stability and error analysis of a new class of higher-order consistent splitting schemes for the

    work page

  7. [7]

    On a new class of

    Huang, Fukeng and Shen, Jie , journal=. On a new class of. 2024 , publisher=

    work page 2024

  8. [8]

    Stability and error analysis of

    Li, Xiaoli and Wang, Weilong and Shen, Jie , journal=. Stability and error analysis of. 2022 , publisher=

    work page 2022

  9. [9]

    Stability and error analysis of a second-order consistent splitting scheme for the

    Huang, Fukeng and Shen, Jie , journal=. Stability and error analysis of a second-order consistent splitting scheme for the. 2023 , publisher=

    work page 2023

  10. [10]

    Analysis of divergence-free

    Han, Yongbin and Hou, Yanren and Zhang, Min , journal=. Analysis of divergence-free

    work page

  11. [11]

    Error analysis of the

    Li, Xiaoli and Shen, Jie , journal=. Error analysis of the. 2020 , publisher=

    work page 2020

  12. [12]

    Numerical approximation of incompressible

    Lin, Lianlei and Yang, Zhiguo and Dong, Suchuan , journal=. Numerical approximation of incompressible. 2019 , publisher=

    work page 2019

  13. [13]

    A Class of New Linear, Efficient and High-Order Implicit-Explicit Methods for the Unsteady

    Wang, Xinhui and Guo, Xu and Li, Xiaoli , journal=. A Class of New Linear, Efficient and High-Order Implicit-Explicit Methods for the Unsteady. 2024 , publisher=

    work page 2024

  14. [14]

    Decoupled Stabilized

    Jia, Xiaofeng and Feng, Hui , journal=. Decoupled Stabilized. 2022 , publisher=

    work page 2022

  15. [15]

    A decoupling method with different subdomain time steps for the non-stationary

    Jia, Huiyong and Shi, Peilin and Li, Kaitai and Jia, Hongen , journal=. A decoupling method with different subdomain time steps for the non-stationary. 2017 , publisher=

    work page 2017

  16. [16]

    A second order in time incremental pressure correction finite element method for the

    Wang, Yunxia and Li, Shishun and Si, Zhiyong , journal=. A second order in time incremental pressure correction finite element method for the

    work page

  17. [17]

    An efficient and long-time accurate third-order algorithm for the

    Chen, Wenbin and Gunzburger, Max and Sun, Dong and Wang, Xiaoming , journal=. An efficient and long-time accurate third-order algorithm for the. 2016 , publisher=

    work page 2016

  18. [18]

    A decoupling method with different subdomain time steps for the nonstationary

    Shan, Li and Zheng, Haibiao and Layton, William J , journal=. A decoupling method with different subdomain time steps for the nonstationary. 2013 , publisher=

    work page 2013

  19. [19]

    Partitioned time stepping method for fully evolutionary

    Shan, Li and Zheng, Haibiao , journal=. Partitioned time stepping method for fully evolutionary. 2013 , publisher=

    work page 2013

  20. [20]

    SIAM Journal on Numerical Analysis , volume=

    Analysis of long time stability and errors of two partitioned methods for uncoupling evolutionary groundwater--surface water flows , author=. SIAM Journal on Numerical Analysis , volume=. 2013 , publisher=

    work page 2013

  21. [21]

    Efficient and long-time accurate second-order methods for the

    Chen, Wenbin and Gunzburger, Max and Sun, Dong and Wang, Xiaoming , journal=. Efficient and long-time accurate second-order methods for the. 2013 , publisher=

    work page 2013

  22. [22]

    Decoupled schemes for a non-stationary mixed

    Mu, Mo and Zhu, Xiaohong , journal=. Decoupled schemes for a non-stationary mixed

    work page

  23. [23]

    A new class of implicit--explicit

    Huang, Fukeng and Shen, Jie , journal=. A new class of implicit--explicit. 2022 , publisher=

    work page 2022

  24. [24]

    Stability and Error Analysis of a Class of High-Order

    Huang, Fukeng and Shen, Jie , journal=. Stability and Error Analysis of a Class of High-Order. 2021 , publisher=

    work page 2021

  25. [25]

    SIAM Journal on Scientific Computing , volume=

    A highly efficient and accurate new scalar auxiliary variable approach for gradient flows , author=. SIAM Journal on Scientific Computing , volume=. 2020 , publisher=

    work page 2020

  26. [26]

    SIAM Review , volume=

    A new class of efficient and robust energy stable schemes for gradient flows , author=. SIAM Review , volume=. 2019 , publisher=

    work page 2019

  27. [27]

    The scalar auxiliary variable (

    Shen, Jie and Xu, Jie and Yang, Jiang , journal=. The scalar auxiliary variable (. 2018 , publisher=

    work page 2018

  28. [28]

    Numerical Methods for Partial Differential Equations , volume=

    Time relaxation algorithm for flow ensembles , author=. Numerical Methods for Partial Differential Equations , volume=. 2016 , publisher=

    work page 2016

  29. [29]

    Numerical Methods for Partial Differential Equations , volume=

    Numerical analysis of two ensemble eddy viscosity numerical regularizations of fluid motion , author=. Numerical Methods for Partial Differential Equations , volume=. 2015 , publisher=

    work page 2015

  30. [30]

    A high order ensemble algorithm for dual-porosity-

    Qiu, Changxin and Hou, Jiangyong and Xia, Yinhua and Shan, Li , journal=. A high order ensemble algorithm for dual-porosity-. 2025 , publisher=

    work page 2025

  31. [31]

    Highly efficient ensemble algorithms for computing the

    Jiang, Nan and Yang, Huanhuan , journal=. Highly efficient ensemble algorithms for computing the. 2024 , publisher=

    work page 2024

  32. [32]

    An artificial compressibility

    Jiang, Nan and Li, Ying and Yang, Huanhuan , journal=. An artificial compressibility. 2021 , publisher=

    work page 2021

  33. [33]

    2021 , publisher=

    Jiang, Nan and Yang, Huanhuan , journal=. 2021 , publisher=

    work page 2021

  34. [34]

    An artificial compressibility ensemble algorithm for a stochastic

    He, Xiaoming and Jiang, Nan and Qiu, Changxin , journal=. An artificial compressibility ensemble algorithm for a stochastic. 2020 , publisher=

    work page 2020

  35. [35]

    An efficient ensemble algorithm for numerical approximation of stochastic

    Jiang, Nan and Qiu, Changxin , journal=. An efficient ensemble algorithm for numerical approximation of stochastic. 2019 , publisher=

    work page 2019

  36. [36]

    An efficient algorithm for computation of

    Mohebujjaman, Muhammad and Rebholz, Leo G , journal=. An efficient algorithm for computation of. 2017 , publisher=

    work page 2017

  37. [37]

    SIAM Journal on Numerical Analysis , volume=

    A second order ensemble timestepping algorithm for natural convection , author=. SIAM Journal on Numerical Analysis , volume=. 2018 , publisher=

    work page 2018

  38. [38]

    SIAM Journal on Scientific Computing , volume=

    Stabilized scalar auxiliary variable ensemble algorithms for parameterized flow problems , author=. SIAM Journal on Scientific Computing , volume=. 2021 , publisher=

    work page 2021

  39. [39]

    IMA Journal of Numerical Analysis , volume=

    An efficient algorithm for simulating ensembles of parameterized flow problems , author=. IMA Journal of Numerical Analysis , volume=. 2019 , publisher=

    work page 2019

  40. [40]

    Journal of Scientific Computing , volume=

    A higher order ensemble simulation algorithm for fluid flows , author=. Journal of Scientific Computing , volume=. 2015 , publisher=

    work page 2015

  41. [41]

    International Journal for Uncertainty Quantification , volume=

    An algorithm for fast calculation of flow ensembles , author=. International Journal for Uncertainty Quantification , volume=. 2014 , publisher=

    work page 2014

  42. [42]

    SIAM/ASA Journal on Uncertainty Quantification , volume=

    Comparison of surrogate-based uncertainty quantification methods for computationally expensive simulators , author=. SIAM/ASA Journal on Uncertainty Quantification , volume=. 2017 , publisher=

    work page 2017

  43. [43]

    Combustion and flame , volume=

    Uncertainty quantification in reacting-flow simulations through non-intrusive spectral projection , author=. Combustion and flame , volume=. 2003 , publisher=

    work page 2003

  44. [44]

    Computer methods in applied mechanics and engineering , volume=

    Stochastic collocation and mixed finite elements for flow in porous media , author=. Computer methods in applied mechanics and engineering , volume=. 2008 , publisher=

    work page 2008

  45. [45]

    SIAM Journal on Numerical Analysis , volume=

    A sparse grid stochastic collocation method for partial differential equations with random input data , author=. SIAM Journal on Numerical Analysis , volume=. 2008 , publisher=

    work page 2008

  46. [46]

    SIAM Journal on Numerical Analysis , volume=

    A stochastic collocation method for elliptic partial differential equations with random input data , author=. SIAM Journal on Numerical Analysis , volume=. 2007 , publisher=

    work page 2007

  47. [47]

    SIAM Journal on Scientific Computing , volume=

    High-order collocation methods for differential equations with random inputs , author=. SIAM Journal on Scientific Computing , volume=. 2005 , publisher=

    work page 2005

  48. [48]

    Kuo, Frances Y and Schwab, Christoph and Sloan, Ian H , journal=. Quasi-. 2012 , publisher=

    work page 2012

  49. [49]

    Multilevel

    Barth, Andrea and Lang, Annika , journal=. Multilevel. 2012 , publisher=

    work page 2012

  50. [50]

    doi:10.1137/140951758 , journal =

    An Adaptive Shifted Power Method for Computing Generalized Tensor Eigenpairs , author =. doi:10.1137/140951758 , journal =

    work page doi:10.1137/140951758

  51. [51]

    Nick Higham , title =

    work page

  52. [52]

    Kolda and Ali Pinar , eprint =

    Chengbin Peng and Tamara G. Kolda and Ali Pinar , eprint =. Accelerating Community Detection by Using

    work page

  53. [53]

    and Zhang, Shanrong and Merritt, Matthew E

    Woessner, Donald E. and Zhang, Shanrong and Merritt, Matthew E. and Sherry, A. Dean , title =. Magnetic Resonance in Medicine , doi =

    work page

  54. [54]

    2003 , eid =

    Properties of Highly Clustered Networks , author =. 2003 , eid =. doi:10.1103/PhysRevE.68.026121 , journal =

    work page doi:10.1103/physreve.68.026121 2003

  55. [55]

    Clawpack Software , author =

    work page

  56. [56]

    : A Document Preparation System

    Leslie Lamport. : A Document Preparation System. 1986

    work page 1986

  57. [57]

    Frank Mittlebach and Michel Goossens , title =

    work page

  58. [58]

    and Van Loan, Charles F

    Golub, Gene H. and Van Loan, Charles F. , title =

    work page

  59. [59]

    Paul Dawkins , title =

    work page

  60. [60]

    User's Guide for the

    work page

  61. [61]

    Michael Downes , title =

    work page

  62. [62]

    Christian Feuers\"anger , title =

    work page