Pith. sign in

REVIEW 3 minor 2 cited by

The polytopal discontinuous Galerkin method for coupled non-Newtonian Stokes-Darcy systems is well-posed and stable with established error bounds.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.3

2026-06-27 09:11 UTC pith:ZOUTXSUS

load-bearing objection The paper sets up a polytopal DG scheme for non-Newtonian Stokes-Darcy coupling and gives a full a-priori analysis via generalized inf-sup, with numerics to match.

arxiv 2606.11935 v1 pith:ZOUTXSUS submitted 2026-06-10 math.NA cs.NA

Polytopal Discontinuous Galerkin Discretizations of Coupled Non-Newtonian Stokes-Darcy Systems

Paola F. Antonietti , Michele Botti , Nicola Parolini , Valentina Pederzoli , Marco Verani This is my paper
classification math.NA cs.NA
keywords polytopal discontinuous Galerkinnon-Newtonian Stokes-Darcygeneralized inf-sup theoryshear-dependent viscosityvelocity-dependent viscositya priori error estimatescoupled flow systemsporous media flow
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper introduces a polytopal discontinuous Galerkin discretization for modeling the interaction between a non-Newtonian fluid in free flow and in a porous medium. It provides a complete a-priori analysis for shear-dependent viscosity in the free-flow region and velocity-dependent viscosity in the porous region. Well-posedness, stability, and error estimates are proven using generalized inf-sup theory, which supports the method's use on meshes with complex geometries. This matters because it enables accurate numerical simulations of such coupled systems without limitations from mesh type or order.

Core claim

The proposed polytopal discontinuous Galerkin method for the coupled non-Newtonian Stokes-Darcy system achieves well-posedness, stability, and optimal error bounds in the framework of generalized inf-sup theory for both shear-dependent and velocity-dependent non-Newtonian viscosity models.

What carries the argument

Polytopal discontinuous Galerkin discretization analyzed through generalized inf-sup theory for the coupled system.

Load-bearing premise

The generalized inf-sup theory framework applies directly to the coupled non-Newtonian Stokes-Darcy system with the chosen viscosity models and polytopal mesh assumptions.

What would settle it

A numerical experiment on a benchmark coupled system where the computed solution fails to converge at the predicted rate or violates stability would falsify the error bounds and well-posedness claims.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • The discretization is suitable for configurations with complex geometries due to its geometric flexibility.
  • Arbitrary-order accuracy is supported by the method.
  • Error estimates are derived for the chosen non-Newtonian models.
  • Numerical results confirm the theoretical error bounds.

Where Pith is reading between the lines

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

  • The framework could be tested on real-world applications like filtration processes or groundwater flow with non-Newtonian fluids.
  • Extensions to time-dependent problems or other interface conditions might follow from the analysis.
  • Implementation on general polytopal meshes could improve efficiency in high-performance computing settings.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

0 major / 3 minor

Summary. The paper proposes a polytopal discontinuous Galerkin discretization for coupled non-Newtonian Stokes-Darcy systems modeling free-flow and porous-medium interactions. It claims a complete a-priori analysis establishing well-posedness, stability, and error bounds via generalized inf-sup theory, covering shear-dependent viscosity in the free-flow region and velocity-dependent viscosity in the porous region, with numerical confirmation of the error estimates.

Significance. If the analysis holds, the work provides a geometrically flexible, arbitrary-order method for a challenging class of multiphysics non-Newtonian flows. The use of polytopal meshes and the extension of generalized inf-sup theory to this coupled setting with distinct viscosity models would be a useful contribution to the numerical analysis of Stokes-Darcy problems.

minor comments (3)
  1. The abstract states that error estimates are confirmed by numerical results, but the manuscript should include a dedicated section or table explicitly comparing observed convergence rates against the predicted orders for both viscosity models.
  2. Notation for the viscosity functions (shear-dependent vs. velocity-dependent) should be introduced with explicit definitions and assumptions on their growth and monotonicity properties early in the analysis section to aid readability.
  3. The polytopal mesh assumptions and the precise form of the DG numerical fluxes at the Stokes-Darcy interface should be stated with a reference to the relevant equation or definition for clarity.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our work on polytopal discontinuous Galerkin methods for non-Newtonian Stokes-Darcy coupling and for recommending minor revision. No major comments were provided in the report.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper presents a standard a-priori analysis establishing well-posedness, stability and error bounds for the polytopal DG scheme on the coupled non-Newtonian Stokes-Darcy system via generalized inf-sup theory. This is an independent mathematical derivation from the discrete formulation and viscosity model assumptions, with no self-definitional reductions, fitted inputs renamed as predictions, or load-bearing self-citation chains that collapse the claims to the paper's own inputs. The analysis relies on established inf-sup frameworks applied to the given setting without circular renaming or ansatz smuggling.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review provides no details on free parameters, axioms, or invented entities; cannot populate ledger entries.

pith-pipeline@v0.9.1-grok · 5658 in / 914 out tokens · 17804 ms · 2026-06-27T09:11:48.333257+00:00 · methodology

0 comments
read the original abstract

We propose and analyze a polytopal discontinuous Galerkin method for the numerical approximation of a coupled non-Newtonian Stokes-Darcy system modeling the interaction between a non-Newtonian free-flow fluid and a non-Newtonian flow through a porous medium. Due to its geometric flexibility and arbitrary-order accuracy, the proposed discretization scheme is well-suited to configurations with complex geometries. We provide a complete a-priori analysis that considers shear-dependent and velocity-dependent non-Newtonian viscosity models for the free-flow and porous media regions, respectively. The well-posedness, stability, and error bounds of the method are established in the framework of generalized inf-sup theory. Error estimates are confirmed by numerical results.

Figures

Figures reproduced from arXiv: 2606.11935 by Marco Verani, Michele Botti, Nicola Parolini, Paola F. Antonietti, Valentina Pederzoli.

Figure 1
Figure 1. Figure 1: Here, ΩS is the free-flow region modeled by Stokes equations, and ΩD is the porous region modeled by Darcy’s law. We split the boundary of the domain in ΓS = ∂ΩS ∩∂Ω and ΓD = ∂ΩD ∩∂Ω. Moreover, we define nS, nD the unit normal vectors to ΓS and ΓD, respectively. We define n S Γ the normal unit vector to the interface Γ directed towards ΩD, and n D Γ = −n S Γ the unit normal vector to 3 [PITH_FULL_IMAGE:fi… view at source ↗
Figure 1
Figure 1. Figure 1: The computational domain. The non-Newtonian Stokes–Darcy coupled system reads: given fS,fD, find (uS, pS) and (uD, pD) such that:    −∇ · (gS(|D(u S )|)D(u S ) − p S I) = fS in ΩS, ∇ · u S = 0 in ΩS, u S = 0 on ΓS, K−1 gD(|u D|)u D + ∇p D = fD in ΩD, ∇ · u D = 0 in ΩD, u D · nD = 0 on ΓD, u S · n S Γ + u D · n D Γ = 0 on Γ, p S − (gS(|D(u S )|)D(u S )n S Γ )… view at source ↗

discussion (0)

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

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

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

    math.NA 2026-07 accept novelty 6.0

    A C1–C0 lowest-order virtual element method for the Stokes–Darcy interface problem in stream-function–pressure form is well-posed, optimally convergent, and mesh-flexible.

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

    math.NA 2026-07 unverdicted novelty 5.0

    A stream function-pressure virtual element method is introduced for the coupled Stokes-Darcy system on polygonal meshes with interface conditions enforced.

Reference graph

Works this paper leans on

48 extracted references · cited by 1 Pith paper

  1. [1]

    P. F. Antonietti, F. Bonaldi, and I. Mazzieri. A high-order discontinuous Galerkin approach to the elasto-acoustic problem.Computer Methods in Applied Mechanics and Engineering, 358:112634, 2020

    2020

  2. [2]

    P. F. Antonietti, S. Bonetti, M. Botti, M. Corti, I. Fumagalli, and I. Mazzieri. lymph: Dis- continuous Polytopal Methods for Multi-Physics Differential Problems.ACM Transaction on Mathematical Software, 51(1), 2025

    2025

  3. [3]

    P. F. Antonietti, A. Cangiani, J. Collis, Z. Dong, E. H. Georgoulis, S. Giani, and P. Houston. Review of discontinuous Galerkin finite element methods for partial differential equations on complicated domains. InBuilding Bridges: Connections and Challenges in Modern Approaches to Numerical Partial Differential Equations, volume 114 ofLecture Notes in Comp...

    2016

  4. [4]

    P. F. Antonietti, C. Facciol` a, A. Russo, and M. Verani. Discontinuous Galerkin approximation of flows in fractured porous media on polytopic grids.SIAM Journal on Scientific Computing, 41(1):A109–A138, 2019. 20

    2019

  5. [5]

    P. F. Antonietti, C. Facciol` a, and M. Verani. Polytopic discontinuous Galerkin methods for the numerical modelling of flow in porous media with networks of intersecting fractures.Journal of Computational and Applied Mathematics, 393:113506, 2021

    2021

  6. [6]

    P. F. Antonietti, S. Giani, and P. Houston.hp-Version Composite Discontinuous Galerkin Meth- ods for Elliptic Problems on Complicated Domains.SIAM Journal on Scientific Computing, 35(3):A1417–A1439, 2013

    2013

  7. [7]

    P. F. Antonietti, L. Mascotto, M. Verani, and S. Zonca. Stability analysis of polytopic dis- continuous Galerkin approximations of the Stokes problem with applications to fluid–structure interaction problems.Journal of Scientific Computing, 90(23), 2022

    2022

  8. [8]

    Arnold, F

    D. Arnold, F. Brezzi, B. Cockburn, and L. Marini. Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems.SIAM Journal on Numerical Analysis, 39(5):1749–1779, 2002

    2002

  9. [9]

    Badia, A

    S. Badia, A. Quaini, and A. Quarteroni. Coupling Biot and Navier-Stokes problems for fluid- poroelastic structure interaction.Technical report, Universitat Politcnica de Catalunya, 2008

    2008

  10. [10]

    Bassi, L

    F. Bassi, L. Botti, A. Colombo, D. Di Pietro, and P. Tesini. On the flexibility of agglomeration based physical space discontinuous Galerkin discretizations.Journal of Computational Physics, 231(1):45–65, 2012

    2012

  11. [11]

    Beavers and D

    G. Beavers and D. Joseph. Boundary conditions at a naturally permeable wall.Journal of Fluid Mechanics, 30(1):197–207, 1967

    1967

  12. [12]

    Becker and M

    R. Becker and M. Braack. The concept of mapped coercivity for nonlinear operators in Banach spaces.Journal of Functional Analysis, 289(3):110893, 2025

    2025

  13. [13]

    Bonaldi, K

    F. Bonaldi, K. Brenner, J. Droniou, and R. Masson. Gradient discretization of two-phase flows coupled with mechanical deformation in fractured porous media.Computers and MAthematics with Applications, 98:40–68, 2021

    2021

  14. [14]

    Botti and L

    M. Botti and L. Mascotto. Sobolev–Poincar´ e inequalities for piecewiseW1,p functions over general polytopic meshes.SIAM Journal on Numerical Analysis, 2026. Accepted for publication

    2026

  15. [15]

    J. Boyd, J. Buick, and S. Green. Analysis of the Casson and Carreau-Yasuda non-Newtonian blood models in steady and oscillatory flows using the lattice Boltzmann method.Physics of Fluids, 19(9):093103, 2007

    2007

  16. [16]

    Burman and P

    E. Burman and P. Hansbo. A unified stabilized method for Stokes’ and Darcy’s equations.Journal of Computational and Applied Mathematics, 198(1):35–51, 2007

    2007

  17. [17]

    Cangiani, Z

    A. Cangiani, Z. Dong, E. Georgoulis, and P. Houston.hp-version discontinuous Galerkin methods on polygonal and polyhedral meshes.Mathematical Models and Methods in Applied Sciences, 24(10):2009–2041, 2014

    2009

  18. [18]

    Cangiani, Z

    A. Cangiani, Z. Dong, E. H. Georgoulis, and P. Houston.hp-Version Discontinuous Galerkin Methods on Polygonal and Polyhedral Meshes. SpringerBriefs in Mathematics. Springer, Cham, 2017

    2017

  19. [19]

    Y. Cao, M. Gunzburger, F. Hua, and X. Wang. Coupled Stokes-Darcy model with Beavers-Joseph interface boundary condition.Communications in Mathematical Sciences, 8(1):1 – 25, 2010

    2010

  20. [20]

    Cesmelioglu, S

    A. Cesmelioglu, S. Rhebergen, and G. Wells. An embedded–hybridized discontinuous Galerkin method for the coupled Stokes–Darcy system.Journal of Computational and Applied Mathemat- ics, 367:112476, 2020

    2020

  21. [21]

    Chhabra and S

    R. Chhabra and S. A. Patel.Non-Newtonian Flow and Applied Rheology, Third Edition. Elsevier, 2025. 21

    2025

  22. [22]

    Chow and G

    S.-S. Chow and G. F. Carey. Numerical approximation of generalized Newtonian fluids using Powell–Sabin–Heindl elements: I. theoretical estimates.International Journal for Numerical Methods in Fluids, 41(10):1085–1118, 2003

    2003

  23. [23]

    Di Pietro and J

    D. Di Pietro and J. Droniou. A hybrid high-order method for Leray-Lions elliptic equations on general meshes.Mathematics of Computations, 86(307):2159–2191, 2017

    2017

  24. [24]

    Di Pietro and J

    D. Di Pietro and J. Droniou.The Hybrid High-Order Method for Polytopal Meshes. Springer, 2020

    2020

  25. [25]

    Di Pietro and A

    D. Di Pietro and A. Ern.Mathematical Aspects of Discontinuous Galerkin Methods. Springer, 2012

    2012

  26. [26]

    Discacciati, E

    M. Discacciati, E. Miglio, and A. Quarteroni. Mathematical and numerical models for coupling surface and groundwater flows.Applied Numerical Mathematics, 43(1):57–74, 2002

    2002

  27. [27]

    Discacciati, A

    M. Discacciati, A. Quarteroni, et al. Navier-Stokes/Darcy coupling: modeling, analysis, and numerical approximation.Revista Matem´ atica Complutense, 22(2):315–426, 2009

    2009

  28. [28]

    Distaccati

    M. Distaccati. Domain decomposition methods for the coupling of surface and groundwater flows. Ph.D. thesis, Ecole Polytechnique Federale de Sausanne, Sausanne, Switzerland, 2024

    2024

  29. [29]

    Dunn and K

    J. Dunn and K. Rajagopal. Fluids of differential type: Critical review and thermodynamic analysis.International Journal of Engineering Science, 33(5):689–729, 1995

    1995

  30. [30]

    Ervin, E

    V. Ervin, E. Jenkins, and S. Sun. Coupling nonlinear Stokes and Darcy flow using mortar finite elements.Applied Numerical Mathematics, 61(11):1198–1222, 2011

    2011

  31. [31]

    V. J. Ervin, E. W. Jenkins, and S. Sun. Coupled Generalized Nonlinear Stokes Flow with Flow through a Porous Medium.SIAM Journal on Numerical Analysis, 47(2):929–952, 2009

    2009

  32. [32]

    Girault and B

    V. Girault and B. 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

    2052

  33. [33]

    Girault and M

    V. Girault and M. F. Wheeler. Numerical discretization of a Darcy–Forchheimer model.Nu- merische Mathematik, 110(2):161–198, 2008

    2008

  34. [35]

    Hanspal, A

    N. Hanspal, A. Waghode, V. Nassehi, and R. Wakeman. Numerical Analysis of Coupled Stokes/Darcy Flows in Industrial Filtrations.Transport in Porous Media, 64:73–101, 2006

    2006

  35. [36]

    Hanspal, A

    N. Hanspal, A. Waghode, V. Nassehi, and R. Wakeman. Development of a predictive mathemati- cal model for coupled Stokes/Darcy flows in cross-flow membrane filtration.Chemical Engineering Journal, 149(1):132–142, 2009

    2009

  36. [37]

    Jingyan and Z

    H. Jingyan and Z. Guanyu. The Well-Posedness of Discontinuous Galerkin Approximation for the Non-Newtonian Stokes-Darcy-Forchheimer Coupling System.Journal of Scientific Computing, 97, 2023

    2023

  37. [38]

    J¨ ager and A

    W. J¨ ager and A. Mikeli´ c. On the Interface Boundary Condition of Beavers, Joseph, and Saffman. SIAM Journal on Applied Mathematics, 60(4):1111–1127, 2000

    2000

  38. [39]

    Kanschat and B

    G. Kanschat and B. Rivi` ere. A strongly conservative finite element method for the coupling of stokes and darcy flow.Journal of Computational Physics, 229(17):5933–5943, 2010. 22

    2010

  39. [40]

    Kanschat and B

    G. Kanschat and B. Rivi` ere. A strongly conservative finite element method for the coupling of stokes and darcy flow.International Journal of Numerical Analysis and Modeling, 11(3):496–524, 2014

    2014

  40. [41]

    R. Li, Y. Gao, C.-S. Zhang, and Z. Chen. A Stokes–Darcy–Darcy model and its discontinuous Galerkin method on polytopic grids.Journal of Computational Physics, 501:112780, 2024

    2024

  41. [42]

    Lipnikov, D

    K. Lipnikov, D. Vassilev, and I. Yotov. Discontinuous Galerkin and mimetic finite difference methods for coupled Stokes–darcy flows on polygonal and polyhedral grids.Numerical Methods for Partial Differential Equations, 30(2):515–543, 2014

    2014

  42. [43]

    Quaini.Algorithms for Fluid-Structure Interaction Problems Arising in Hemodynamics

    A. Quaini.Algorithms for Fluid-Structure Interaction Problems Arising in Hemodynamics. PhD thesis, Ecole Polytechnique F´ ed´ erale de Lausanne, Switzerland, 2008

    2008

  43. [44]

    Rivi` ere

    B. Rivi` ere. Analysis of a Discontinuous Finite Element Method for the Coupled Stokes and Darcy Problems.J Sci Comput, 22:479–500, 2005

    2005

  44. [45]

    Rivi` ere.Discontinuous Galerkin Methods for Solving Elliptic and Parabolic Equations

    B. Rivi` ere.Discontinuous Galerkin Methods for Solving Elliptic and Parabolic Equations. Society for Industrial and Applied Mathematics, 2008

    2008

  45. [46]

    Rivi` ere and I

    B. Rivi` ere and I. Yotov. Locally Conservative Coupling of Stokes and Darcy Flows.SIAM Journal on Numerical Analysis, 42(5):1959–1977, 2005

    1959

  46. [47]

    P. G. Saffman. On the Boundary Condition at the Surface of a Porous Medium.Studies in Applied Mathematics, 50(2):93–101, 1971

    1971

  47. [48]

    Santesarti, M

    G. Santesarti, M. Marino, F. Viola, R. Verzicco, and G. Vairo. A quasi-analytical solution for Carreau–Yasuda-type shear-thinning flows in slightly tapered pipes using a truncated power-law model.Journal of Non-Newtonian Fluid Mechanics, 349:105569, 2026

    2026

  48. [49]

    Ye and S

    X. Ye and S. Zhang. A conforming discontinuous Galerkin finite element method for the Stokes problem on polytopal meshes.Journal of Scientific Computing, 87:1–26, 2021. A Appendix In this section, we provide a complete proof of the continuity ofb h(·,·), i.e., |bh(vh, ph)|≲| | |vh| | |div| | |ph| | |Mh ∀vh ∈X h,∀p h ∈M h. 23 Proof.From the definitions of ...

    2021