REVIEW 2 major objections 2 minor 90 references
An unfitted finite-element scheme for Darcy flow keeps exact local mass conservation even when cut cells are arbitrarily small.
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 →
T0 review · grok-4.5
2026-07-13 17:41 UTC pith:3ID3TCS6
load-bearing objection Package mismatch: abstract is unfitted H(div) Darcy FEM, full text is an unrelated soft X-ray transient paper; claims cannot be checked. the 2 major comments →
Divergence-free unfitted finite element discretisations for the Darcy problem
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
An unfitted compatible discretisation of the Darcy problem, based on H(div)-conforming flux spaces and discontinuous pressure spaces, can be stabilised by a combination of L2 flux stabilisation and a mixed-term ghost penalty so that pointwise discrete mass conservation is preserved while all stability and a-priori error constants remain independent of the cut configuration.
What carries the argument
The mixed-term ghost-penalty stabilisation (realised either cell-wise or face-based) that restores pressure control without destroying the local conservation structure of the H(div)-discontinuous-pressure pair.
Load-bearing premise
That the two stabilisation terms can be scaled so they simultaneously restore pressure control and leave the local conservation structure intact for every possible cut configuration.
What would settle it
A single numerical experiment on a sequence of meshes that produce arbitrarily small cut cells, with pure pressure boundary data, in which either the discrete mass residual fails to reach solver tolerance or the flux error constant grows unboundedly as the cut size tends to zero.
If this is right
- Unfitted meshes can be used for Darcy-type porous-media problems without sacrificing cell-wise mass conservation.
- Conditioning and error constants remain cut-independent, so adaptive or immersed-boundary geometries become practical.
- Pressure-robust flux estimates hold under pure pressure boundary conditions, improving accuracy for velocity-driven quantities of interest.
- The augmented-Lagrangian variant supplies a well-conditioned algebraic system that admits efficient preconditioning strategies.
Where Pith is reading between the lines
- The same mixed-term stabilisation idea may transfer to other mixed formulations (Stokes, Brinkman, poroelasticity) that require both H(div) conformity and cut robustness.
- Because mass conservation is preserved up to solver tolerance, the scheme is a natural candidate for multiphase or transport problems that couple tightly to the Darcy velocity.
- Face-based versus bulk ghost-penalty realisations give implementers a practical trade-off between sparsity and ease of assembly that can be tuned per application.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The abstract of arXiv:2603.26212 claims an unfitted compatible finite-element scheme for the Darcy problem that uses H(div)-conforming fluxes and discontinuous pressures, preserves pointwise discrete mass conservation on arbitrarily small cut cells, and obtains cut-independent stability and a priori error estimates via an L2 flux stabilisation combined with a mixed-term (bulk or face-based ghost-penalty) stabilisation. Mixed boundary conditions are weakly imposed, pressure-robust flux bounds are asserted for pure pressure BCs, and an augmented-Lagrangian variant is proposed for better constraint control and preconditioning. The numerical claims include optimal rates, cut-independent conditioning and mass conservation up to solver tolerance. The full manuscript text supplied for review, however, is an unrelated astrophysics article on a soft X-ray dirty fireball (arXiv-style 2603.26213); none of the stated proofs, lemmas, cut-configuration analysis or Darcy numerics appear.
Significance. If the abstract claims were substantiated by a correct analysis and supporting experiments, the work would be a useful contribution to unfitted mixed methods for porous-media flow: cut-robust, locally conservative H(div)–L2 schemes with pressure-robust flux estimates remain of practical interest for complex geometries. The combination of L2 flux stabilisation with a mixed-term ghost penalty that is claimed not to destroy local conservation, together with an AL variant amenable to preconditioning, would be a concrete technical advance. Because the supplied full text contains none of that material, the significance of the actual submission cannot be assessed.
major comments (2)
- The review package is misassembled: the abstract and paper_id describe an unfitted H(div)–discontinuous Darcy FEM (arXiv:2603.26212), yet the full manuscript text is an astrophysics paper on an energetic dirty fireball detected in soft X-rays. Consequently every load-bearing claim—pointwise discrete mass conservation under the proposed stabilisations, cut-independent stability and a priori constants, pressure-robust flux bounds for pure pressure BCs, and the AL variant—cannot be verified. A correct manuscript matching the abstract must be supplied before any scientific evaluation is possible.
- Even from the abstract alone, the central technical assertion (that L2 flux stabilisation plus the mixed-term bulk/face ghost-penalty restores pressure control for every cut configuration without destroying local conservation) is parameter-dependent and typically requires careful scaling of the free stabilisation coefficients. Without the proofs, parameter ranges, or cut-configuration experiments that the abstract promises, this claim remains uncheckable and is therefore a blocker for acceptance.
minor comments (2)
- Once a correct manuscript is submitted, the free parameters listed in the abstract (L2 flux coefficient, mixed-term/ghost-penalty parameters for bulk and face realisations, AL penalty) should be stated with their admissible ranges and scaling with mesh size and cut fraction.
- The abstract mentions both bulk and face-based ghost-penalty realisations; the eventual paper should make clear which (if either) is preferred for conditioning and for the pressure-robust flux estimate.
Circularity Check
No circularity detectable: abstract describes a standard FEM analysis paper; supplied full text is an unrelated astrophysics manuscript, so no derivation chain exists to reduce.
full rationale
The abstract of arXiv:2603.26212 claims an unfitted H(div)–discontinuous Darcy scheme that preserves pointwise discrete mass conservation, is robust to small cut cells via L2 flux plus mixed-term (bulk or face) ghost-penalty stabilisation, admits cut-independent stability and a priori estimates, and yields pressure-robust flux bounds under pure pressure BCs, plus an augmented-Lagrangian variant. These are ordinary numerical-analysis claims of the form ‘scheme X is stable/convergent under assumptions Y’; they do not recycle fitted parameters as predictions, invoke self-citation uniqueness theorems, or rename known empirical patterns. The CACHEABLE full-text block, however, is an entirely different paper (soft X-ray dirty fireball, arXiv-style 2603.26213). Consequently no equations, proofs, or self-citations from the Darcy FEM work are available to inspect. With no derivation chain present, no circular reduction can be exhibited. Score 0 is therefore the only honest outcome under the hard rules (no speculation, quote-required evidence only).
Axiom & Free-Parameter Ledger
free parameters (3)
- L2 flux stabilisation coefficient
- mixed-term / ghost-penalty stabilisation parameters (bulk and face-based)
- augmented Lagrangian penalty parameter
axioms (5)
- domain assumption Continuous Darcy problem (flux in H(div), pressure in L², constitutive law relating flux to pressure gradient and permeability) is well-posed under standard mixed boundary conditions.
- domain assumption Unfitted meshes may produce arbitrarily small cut cells; analysis must be independent of cut configuration.
- standard math H(div)-conforming flux spaces paired with discontinuous pressure spaces yield a discrete complex that supports pointwise (cell-wise) mass conservation when the divergence of the flux space lands in the pressure space.
- domain assumption Ghost-penalty / bulk stabilisation can control jumps or gradients across cut faces without destroying consistency at the optimal rate.
- ad hoc to paper Weak imposition of flux and pressure traces on unfitted boundaries is consistent and stable when combined with the interior stabilisations.
read the original abstract
We develop an unfitted compatible finite element discretisation for the Darcy problem based on $H(\mathrm{div})$-conforming flux spaces and discontinuous pressure spaces. The method is designed to preserve pointwise discrete mass conservation while remaining robust in the presence of arbitrarily small cut cells arising from unfitted meshes. Robustness is achieved by combining an $L^2$-stabilisation of the flux with an additional mixed-term stabilisation that enhances pressure control without destroying the local conservation structure. We consider both cell-wise (bulk) and face-based ghost-penalty realisations of the stabilisation. Mixed boundary conditions are handled by weak imposition of both flux and pressure traces on unfitted boundaries. We prove stability and a priori error estimates with constants independent of the cut configuration, and establish pressure-robust flux error bounds in the case of pure pressure boundary conditions. We also introduce an augmented Lagrangian variant that improves control of the conservation constraint and is amenable to efficient preconditioning strategies. Numerical experiments for a range of cut configurations, boundary-condition regimes and parameter choices confirm the theoretical results, demonstrating optimal convergence, cut-independent conditioning and mass conservation up to solver tolerance.
Reference graph
Works this paper leans on
-
[1]
S. E. Woosley, J. S. Bloom,Annu. 剥v. 䅳瑲on. 䅳瑲潰桹献44, 507 (2006)
2006
-
[2]
R. W. Klebesadel, I. B. Strong, R. A. Olson,䅳瑲潰桹献 J. 䱥琮182, L85 (1973)
1973
-
[3]
Zhang,The Physics 潦 Lauua-Ray Bur獴猬 ISBN: 㤷㠭ㄭㄳ㤭㈲㘵㌭〮 䍡浢物摧e 啮椭 veristy 偲ess(2018)
B. Zhang,The Physics 潦 Lauua-Ray Bur獴猬 ISBN: 㤷㠭ㄭㄳ㤭㈲㘵㌭〮 䍡浢物摧e 啮椭 veristy 偲ess(2018)
2018
-
[4]
Lithwick, R
Y . Lithwick, R. Sari,䅳瑲潰桹献 J.555, 540 (2001)
2001
-
[5]
J. L. Racusin,整 al.,䅳瑲潰桹献 J.738, 138 (2011)
2011
-
[6]
Campana,整 al.,Nature442, 1008 (2006)
S. Campana,整 al.,Nature442, 1008 (2006)
2006
-
[7]
Sakamoto,整 al.,䅳瑲潰桹献 J.679, 570 (2008)
T. Sakamoto,整 al.,䅳瑲潰桹献 J.679, 570 (2008)
2008
- [8]
-
[9]
Amati,整 al.,䅳瑲on
L. Amati,整 al.,䅳瑲on. 䅳瑲潰桹献390, 81 (2002)
2002
-
[10]
Paczy ´nski,䅳瑲潰桹献 J
B. Paczy ´nski,䅳瑲潰桹献 J. 䱥琮494, L45 (1998)
1998
-
[11]
C. D. Dermer, J. Chiang, M. B ¨ottcher,䅳瑲潰桹献 J.513, 656 (1999)
1999
-
[12]
Y . F. Huang, Z. G. Dai, T. Lu,䵯渮 乯琮 刮 䅳瑲on. 卯挮332, 735 (2002)
2002
-
[13]
J. E. Rhoads,䅳瑲潰桹献 J.591, 1097 (2003)
2003
-
[14]
Zhang, S
W. Zhang, S. E. Woosley, A. Heger,䅳瑲潰桹献 J.608, 365 (2004)
2004
-
[15]
S. B. Cenko,整 al.,䅳瑲潰桹献 J.769, 130 (2013)
2013
-
[16]
A. Y . Q. Ho,整 al.,䅳瑲潰桹献 J.938, 85 (2022)
2022
- [17]
-
[18]
Liu,整 al.,Nature Astronouy(2025)
Y . Liu,整 al.,Nature Astronouy(2025). 16
2025
-
[19]
Burns, M
E. Burns, M. E. Ravasio, P. G. Jonker, Fermi-GBM Team,LRB 䍯潲dinates Network 38238, 1 (2024)
2024
-
[20]
Materials and methods are available as supplementary materials
-
[21]
Gruber,整 al.,䅳瑲潰桹献 J
D. Gruber,整 al.,䅳瑲潰桹献 J. Supp.211, 12 (2014)
2014
-
[22]
Band,整 al.,䅳瑲ophys
D. Band,整 al.,䅳瑲ophys. J.413, 281 (1993)
1993
-
[23]
Kaneko,整 al.,䅳瑲ophys
Y . Kaneko,整 al.,䅳瑲ophys. J. Supp.166, 298 (2006)
2006
-
[24]
L. Nava, G. Ghirlanda, G. Ghisellini, A. Celotti,䅳瑲on. 䅳瑲潰桹献530, A21 (2011)
2011
-
[25]
Z. Y . Liu,整 al.,LRB 䍯潲dinat敳 Network38211, 1 (2024)
2024
-
[26]
Zhang,整 al.,䅳瑲ophys
B. Zhang,整 al.,䅳瑲ophys. J.642, 354 (2006)
2006
-
[27]
J. A. Nousek,整 al.,Astr潰桹献 J.642, 389 (2006)
2006
-
[28]
Rossi,整 al.,LRB Coordinates Network38233, 1 (2024)
A. Rossi,整 al.,LRB Coordinates Network38233, 1 (2024)
2024
-
[29]
J. Wei,整 al.,arXiv e-printsp. arXiv:1610.06892 (2016)
Pith/arXiv arXiv 2016
-
[30]
J. B. Oke,整 al.,Publ. Astron. 卯挮 Pac.107, 375 (1995)
1995
- [31]
-
[32]
Tagliaferri,整 慬.,Nature436, 985 (2005)
G. Tagliaferri,整 慬.,Nature436, 985 (2005)
2005
-
[33]
Zhang, E.-W
B.-B. Zhang, E.-W. Liang, B. Zhang,䅳瑲潰桹献 J.666, 1002 (2007)
2007
-
[34]
Kumar, A
P. Kumar, A. Panaitescu,䅳瑲潰桹献 J. 䱥琮541, L51 (2000)
2000
-
[35]
Zhang, B
B.-B. Zhang, B. Zhang, E.-W. Liang, X.-Y . Wang,䅳瑲潰桹献 J. 䱥琮690, L10 (2009)
2009
-
[36]
Zhang, P
B. Zhang, P. M ´esz´aros,䅳瑲潰桹献 J.571, 876 (2002)
2002
-
[37]
Rossi, D
E. Rossi, D. Lazzati, M. J. Rees,䵯渮 乯琮 刮 䅳瑲on. 卯挮332, 945 (2002)
2002
-
[38]
Z. G. Dai, T. Lu,䅳瑲on. 䅳瑲潰桹献333, L87 (1998). 17
1998
-
[39]
R. Shen, C. D. Matzner,䅳瑲潰桹献 J.744, 36 (2012)
2012
-
[40]
Dereli-B ´egu´e,整 al.,Nature 䍯浭畮楣慴楯湳13, 5611 (2022)
H. Dereli-B ´egu´e,整 al.,Nature 䍯浭畮楣慴楯湳13, 5611 (2022)
2022
-
[41]
Z. G. Dai, T. Lu,䵯渮 乯琮 刮 䅳瑲on. 卯挮298, 87 (1998)
1998
-
[42]
R. A. Chevalier, Z.-Y . Li,䅳瑲潰桹献 J. 䱥琮520, L29 (1999)
1999
-
[43]
A. I. MacFadyen, S. E. Woosley,䅳瑲潰桹献 J.524, 262 (1999)
1999
-
[44]
Narayan, T
R. Narayan, T. Piran, P. Kumar,䅳瑲潰桹献 J.557, 949 (2001)
2001
-
[45]
R. D. Blandford, R. L. Znajek,䵯渮 乯琮 刮 䅳瑲on. 卯挮179, 433 (1977)
1977
-
[46]
W.-H. Lei, B. Zhang, E.-W. Liang,Astr潰桹献 J.765, 125 (2013)
2013
-
[47]
Zhang, Z
D. Zhang, Z. G. Dai,䅳瑲ophys. J.703, 461 (2009)
2009
-
[48]
Zhang, S
W. Zhang, S. E. Woosley, A. I. MacFadyen,䅳瑲潰桹献 J.586, 356 (2003)
2003
-
[49]
Zhang,整 al.,䅳瑲ophys
C. Zhang,整 al.,䅳瑲ophys. J. 䱥琮941, L2 (2022)
2022
-
[50]
Cheng,整 al.,Experiuental 䅳瑲onouy57, 10 (2024)
H. Cheng,整 al.,Experiuental 䅳瑲onouy57, 10 (2024)
2024
-
[51]
Chen,整 al.,Space Telescopes and Instruuentation 2020: 啬瑲aviolet 瑯 Lauua Ray, J.-W
Y . Chen,整 al.,Space Telescopes and Instruuentation 2020: 啬瑲aviolet 瑯 Lauua Ray, J.-W. A. den Herder, S. Nikzad, K. Nakazawa, eds. (2020), vol. 11444 ofSociety 潦 偨潴漭 Optical Instruuentation 䕮杩湥敲s ⡓偉䔩 䍯湦敲敮捥 卥物敳, p. 114445B
2020
-
[52]
A. Goldstein,整 al.,arXiv 攭灲楮瑳p. arXiv:1903.12597 (2019)
Pith/arXiv arXiv 1903
-
[53]
Connaughton,整 al.,䅳瑲潰桹献 J
V . Connaughton,整 al.,䅳瑲潰桹献 J. Supp.216, 32 (2015)
2015
-
[54]
Poolakkil,整 al.,䅳瑲潰桹献 J.913, 60 (2021)
S. Poolakkil,整 al.,䅳瑲潰桹献 J.913, 60 (2021)
2021
-
[55]
Vianello,䅳瑲潰桹献 J
G. Vianello,䅳瑲潰桹献 J. Supp.236, 17 (2018)
2018
-
[56]
Pedichini,整 al.,䥮獴牵浥湴 Design and P敲景牭慮捥 for Optical/Infr慲敤 䝲潵湤ⵢ慳敤 Telescopes, M
F. Pedichini,整 al.,䥮獴牵浥湴 Design and P敲景牭慮捥 for Optical/Infr慲敤 䝲潵湤ⵢ慳敤 Telescopes, M. Iye, A. F. M. Moorwood, eds. (2003), vol. 4841 ofSociety 潦 Photo-Optical Instruuentation 䕮杩湥敲s ⡓偉䔩 䍯湦敲敮捥 卥物敳, pp. 815–826. 18
2003
-
[57]
Fontana,整 al.,䅳瑲on
A. Fontana,整 al.,䅳瑲on. 䅳瑲潰桹献570, A11 (2014)
2014
-
[58]
I. S. McLean,整 al.,䝲潵湤ⵢ慳敤 and 䅩牢潲湥 Instruuentation for 䅳瑲onouy 䥖, I. S. McLean, S. K. Ramsay, H. Takami, eds. (2012), vol. 8446 ofSociety 潦 Photo-Optical Instruuentation 䕮杩湥敲s ⡓偉䔩 䍯湦敲敮捥 卥物敳, p. 84460J
2012
-
[59]
D. Lang, D. W. Hogg, K. Mierle, M. Blanton, S. Roweis,䅳瑲on. J.139, 1782 (2010)
2010
-
[60]
M. F. Skrutskie,整 慬.,䅳瑲on. J.131, 1163 (2006)
2006
-
[61]
䅳瑲潰桹献674, A1 (2023)
Gaia Collaboration,整 al.,䅳瑲on. 䅳瑲潰桹献674, A1 (2023)
2023
-
[62]
K. C. Chambers,整 慬.,arXiv 攭灲楮瑳p. arXiv:1612.05560 (2016)
Pith/arXiv arXiv 2016
-
[63]
Prochaska,整 al.,T桥 J潵牮慬 潦 佰敮 卯畲捥 Software5, 2308 (2020)
J. Prochaska,整 al.,T桥 J潵牮慬 潦 佰敮 卯畲捥 Software5, 2308 (2020)
2020
-
[64]
V . S. Dhillon,整 al.,䵯渮 乯琮 刮 䅳瑲on. 卯挮507, 350 (2021)
2021
-
[65]
A. C. Carnall, R. J. McLure, J. S. Dunlop, R. Dav ´e,䵯渮 乯琮 刮 䅳瑲on. 卯挮480, 4379 (2018)
2018
-
[66]
A. C. Carnall,整 al.,䅳瑲潰桹献 J.873, 44 (2019)
2019
-
[67]
A. Wang, T. An, C. Dai, X.-Y . Wang, EP radio follow-up Team,LRB 䍯潲dinates Network 38659, 1 (2024)
2024
-
[68]
Yang,整 al.,Nature612, 232 (2022)
J. Yang,整 al.,Nature612, 232 (2022)
2022
-
[69]
Yang,整 al.,䅳瑲潰桹献 J
J. Yang,整 al.,䅳瑲潰桹献 J. 䱥琮947, L11 (2023)
2023
-
[70]
J. M. Dickey, F. J. Lockman,Annu. 剥v. 䅳瑲on. 䅳瑲潰桹献28, 215 (1990)
1990
-
[71]
P. M. W. Kalberla,整 慬.,䅳瑲on. 䅳瑲潰桹献440, 775 (2005)
2005
-
[72]
Willingale, R
R. Willingale, R. L. C. Starling, A. P. Beardmore, N. R. Tanvir, P. T. O’Brien,䵯渮 乯琮 刮 䅳瑲on. 卯挮431, 394 (2013)
2013
-
[73]
䅳瑲潰桹献594, A116 (2016)
HI4PI Collaboration,整 al.,䅳瑲on. 䅳瑲潰桹献594, A116 (2016). 19
2016
-
[74]
M. E. Ravasio,整 al.,䅳瑲on. 䅳瑲潰桹献613, A16 (2018)
2018
-
[75]
M. E. Ravasio, G. Ghirlanda, L. Nava, G. Ghisellini,䅳瑲on. 䅳瑲潰桹献625, A60 (2019)
2019
-
[76]
de la Cruz-Dombriz, P
´A. de la Cruz-Dombriz, P. K. S. Dunsby, O. Luongo, L. Reverberi,J潵牮慬 潦 Cosuolo杹 and 䅳瑲oparticle Physics2016, 042 (2016)
2016
-
[77]
K. P. Burnham, D. R. Anderson,Sociological Methods & 剥獥慲ch33, 261 (2004)
2004
-
[78]
T. P. Li, Y . Q. Ma,䅳瑲潰桹献 J.272, 317 (1983)
1983
-
[79]
Zhang,整 al.,䅳瑲ophys
B. Zhang,整 al.,䅳瑲ophys. J.703, 1696 (2009)
2009
-
[80]
B. B. Zhang,整 al.,Nature 䅳瑲onouy2, 69 (2018)
2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.