Pith. sign in

REVIEW 3 minor 1 cited by

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

T0 review · 0 major / 3 minor · reviewed 2026-06-27 · grok-4.3

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

Pith's one-line read The polytopal discontinuous Galerkin method for coupled non-Newtonian Stokes-Darcy systems is well-posed and stable with established error bounds.

desk verdict 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. read the letter →

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

Paola F. Antonietti , Michele Botti , Nicola Parolini , Valentina Pederzoli , Marco Verani This is my paper · ORCID
classification math.NAcs.NA
keywords polytopaldiscontinuousGalerkinnon-NewtonianStokes-Darcygeneralizedinf-suptheoryshear-dependentviscosityvelocity-dependentapriorierrorestimatescoupledflowsystemsporousmedia
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

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.

What carries the argument

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

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

Extended reading notes

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.

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.

Editorial extensions

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.

Reading between the lines

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.
Share X Bluesky LinkedIn Reddit HN

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 · score 0.0 of 10

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.

Assumptions & free parameters 0 free parameters · 0 assumptions · 0 invented entities

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

how reviews work

0 comments
Cite this review

Pith. "Pith review of Polytopal Discontinuous Galerkin Discretizations of Coupled Non-Newtonian Stokes-Darcy Systems." pith.science (2026). https://pith.science/paper/ZOUTXSUS

@misc{pith2026260611935,
  author       = {Pith},
  title        = {Pith review of: Polytopal Discontinuous Galerkin Discretizations of Coupled Non-Newtonian Stokes-Darcy Systems},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZOUTXSUS}},
  note         = {Machine review of arXiv:2606.11935}
}
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 the authors.

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:figures/full… view at source ↗
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 Γ ) · n S Γ =… view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

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

    math.NA 2026-07 unverdicted novelty 6.0 of 10

    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

Show all 48 references
  1. [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

  2. [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

  3. [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

  4. [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

  5. [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

  6. [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

  7. [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

  8. [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

  9. [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

  10. [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

  11. [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

  12. [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

  13. [21]

    Chhabra and S

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

    2025

  14. [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

  15. [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

  16. [24]

    Di Pietro and J

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

    2020

  17. [25]

    Di Pietro and A

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

    2012

  18. [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

  19. [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

  20. [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

  21. [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

  22. [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

  23. [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

  24. [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

  25. [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

  26. [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

  27. [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

  28. [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

  29. [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

  30. [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

  31. [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

  32. [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

  33. [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

  34. [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

  35. [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

  36. [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

  37. [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

  38. [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

  39. [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

  40. [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)|≲| |...

    2021

Pith tools

Reviewed June 27, 2026 · model on record in the stance chip above.