pith. sign in

arxiv: 2604.02890 · v1 · submitted 2026-04-03 · 🧮 math.NA · cs.NA

A posteriori error estimates for mixed finite element discretization of the multigroup Neutron Simplified Transport equations with Robin boundary condition

Patrick Ciarlet (POEMS) , Minh-Hieu Do (SERMA) , Mario Gervais (SERMA) , Fran\c{c}ois Madiot (SERMA) This is my paper

Pith reviewed 2026-05-13 18:48 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords a posteriori error estimatesmixed finite elementsSPN equationsRobin boundary conditiondomain decompositionadaptive mesh refinementneutron transport
0
0 comments X

The pith

Guaranteed and locally efficient a posteriori error estimators are derived for mixed finite element discretizations of the multigroup SPN equations with Robin boundary conditions.

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

The paper derives a posteriori error estimates for the mixed finite element discretization of the multigroup neutron simplified transport equations when Robin boundary conditions are imposed. These estimates are shown to be guaranteed upper bounds on the discretization error while also being locally efficient, which supports reliable adaptive mesh refinement. A dedicated estimator component is constructed to control the contribution from the Robin boundary term, and the same framework is extended to a domain decomposition method that incorporates L2 jump terms across subdomain interfaces. The results apply to both simplicial and Cartesian meshes and cover cases with mixed boundary conditions of Dirichlet, Neumann, or Fourier type. If correct, the estimates give practitioners a practical tool to drive mesh adaptation in neutron transport calculations without requiring over-refinement of the entire domain.

Core claim

For the mixed finite element approximation of the multigroup SPN equations with Robin boundary conditions, residual-based estimators are constructed that deliver a guaranteed upper bound on the error in the natural energy norm and locally efficient lower bounds on the same quantity. The estimators include volume residuals, inter-element flux jumps, and an explicit boundary residual term adapted to the Robin condition; the analysis is carried out on both simplicial and Cartesian meshes. The same estimator construction is then applied inside a non-overlapping domain decomposition scheme with L2 jump penalization, yielding an adaptive refinement strategy that remains reliable across subdomain 0

What carries the argument

Residual-based a posteriori error estimators for the mixed formulation that include a dedicated Robin-boundary residual term.

If this is right

  • The estimators directly drive adaptive mesh refinement that reduces the global error at optimal rates.
  • The same estimator construction remains reliable when the domain is split into subdomains with Cartesian meshes and L2 interface jumps.
  • Separate residual terms allow the theory to cover any combination of Dirichlet, Neumann, or Robin boundary conditions on different parts of the boundary.
  • Local efficiency guarantees that the estimator does not over-estimate the error on any single element or subdomain.

Where Pith is reading between the lines

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

  • The estimator construction could be reused for other first-order transport models that admit a similar mixed variational formulation.
  • Combining the DD+L2 jumps framework with the estimators suggests a route to scalable parallel adaptive solvers for large reactor geometries.
  • The reliance on standard approximation properties implies the method should extend without change to higher-order mixed elements once those properties are verified.
  • Numerical experiments in the paper already illustrate the approach; repeating them with a known exact solution would provide an independent check of the constants appearing in the bounds.

Load-bearing premise

The continuous SPN problem with Robin boundary conditions is assumed to be well-posed and the mixed finite element spaces are assumed to satisfy standard approximation properties on the given meshes.

What would settle it

A concrete numerical counter-example on a single-element or low-refinement mesh where the computed estimator fails to upper-bound the true error measured in the energy norm would falsify the guarantee.

Figures

Figures reproduced from arXiv: 2604.02890 by Fran\c{c}ois Madiot (SERMA), Mario Gervais (SERMA), Minh-Hieu Do (SERMA), Patrick Ciarlet (POEMS).

Figure 1
Figure 1. Figure 1: Description of the AMR process. 25 [PITH_FULL_IMAGE:figures/full_fig_p025_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The benchmark geometry. Source Reflector Core Control rod S 1 f 0. 9.09319 × 10−3 0. S 2 f 0. 2.90183 × 10−1 0 [PITH_FULL_IMAGE:figures/full_fig_p027_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Subdomains and initial mesh for the multi-domain approach. [PITH_FULL_IMAGE:figures/full_fig_p028_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Relative error in the ∥ · ∥S norm (left) and maximum of the total estimator (right) as a function of the total number of mesh elements. Figures 5, 6 and 7 respectively show the final meshes of the mono-domain, DDM-1 and DDM-2 multi-domain configurations. We observe that refinement mostly takes place near the material interfaces. The DDM-based refinement is able to focus on this interface-focused refinement… view at source ↗
Figure 5
Figure 5. Figure 5: Final mesh for the mono-domain configurations. [PITH_FULL_IMAGE:figures/full_fig_p030_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Final mesh for the DDM-1 multi-domain configuration. [PITH_FULL_IMAGE:figures/full_fig_p031_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Final mesh for the DDM-2 multi-domain configuration. [PITH_FULL_IMAGE:figures/full_fig_p031_7.png] view at source ↗
read the original abstract

We analyse a posteriori error estimates for the discretization with mixed finite elements on simplicial or Cartesian meshes of the multigroup neutron simplified transport (SPN ) equations, in the case where a Robin (or Fourier type) boundary condition is imposed on the boundary. This boundary condition is of particular importance in neutronics, since it corresponds to the well-known vacuum boundary condition. We provide guaranteed and locally efficient estimators. In particular, a specific estimator is designed to handle the Robin boundary condition. We also develop the theory in the case of mixed imposed boundary conditions, of Dirichlet, Neumann or Fourier type. The approach is further extended to a Domain Decomposition Method, the so-called DD+L 2 jumps method. In this framework, the adaptive mesh refinement strategy is implemented for a discretization using Cartesian meshes on each subdomain. Numerical experiments illustrate the theory.

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

1 major / 2 minor

Summary. The manuscript derives guaranteed and locally efficient a posteriori error estimators for mixed finite element discretizations of the multigroup simplified neutron transport (SPN) equations subject to Robin boundary conditions on simplicial or Cartesian meshes. It also treats mixed Dirichlet/Neumann/Fourier boundary conditions and extends the estimators to a domain decomposition method (DD + L2 jumps) with adaptive refinement on Cartesian subdomains, illustrated by numerical experiments.

Significance. If the reliability proof is complete, the work supplies a practical, computable error control tool for neutronics applications where vacuum (Robin) boundary conditions are standard. The DD extension adds value for large-scale simulations. The contribution is incremental rather than foundational, as it builds on standard residual-based a posteriori techniques for mixed formulations.

major comments (1)
  1. [derivation of the a posteriori upper bound] The reliability (guaranteed upper bound) of the estimators is obtained from the stability of the continuous mixed variational formulation. The manuscript invokes well-posedness of the multigroup SPN problem with Robin boundary conditions without providing a proof, a specific reference, or an explicit stability constant for the coupled system (scattering and fission cross-sections link the groups). Standard single-group arguments do not automatically carry over; if the continuous problem is not stable for some physically relevant data, the reliability constant becomes uncontrolled and the 'guaranteed' claim fails.
minor comments (2)
  1. [Abstract] The abstract introduces the DD+L2 jumps method without a brief parenthetical definition or forward reference; expand on first use for readability.
  2. [Problem formulation] Notation for the multigroup flux and cross-section matrices should be introduced once in a dedicated notation subsection rather than piecemeal.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the detailed review and for highlighting the need to explicitly justify the well-posedness assumption underlying the reliability proof. We address this point below and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [derivation of the a posteriori upper bound] The reliability (guaranteed upper bound) of the estimators is obtained from the stability of the continuous mixed variational formulation. The manuscript invokes well-posedness of the multigroup SPN problem with Robin boundary conditions without providing a proof, a specific reference, or an explicit stability constant for the coupled system (scattering and fission cross-sections link the groups). Standard single-group arguments do not automatically carry over; if the continuous problem is not stable for some physically relevant data, the reliability constant becomes uncontrolled and the 'guaranteed' claim fails.

    Authors: We agree that the continuous multigroup SPN problem with Robin boundary conditions requires explicit justification for the reliability result to be fully rigorous. The well-posedness follows from standard coercivity and continuity arguments for the mixed formulation under the usual physical assumptions (non-negative, bounded cross-sections with the fission operator satisfying a spectral radius condition less than one). In the revised version we will (i) cite a specific reference establishing well-posedness for the multigroup SPN system with Robin conditions (e.g., the analysis in the neutronics literature for the diffusion approximation of the transport equation), and (ii) add a short appendix sketching the derivation of the stability constant for the coupled system by viewing the scattering/fission terms as a compact perturbation and applying a Neumann-series argument. This will make the dependence of the reliability constant on the data explicit and controlled for all physically admissible cross-section data. No change to the estimator construction or numerical results is required. revision: yes

Circularity Check

0 steps flagged

No significant circularity; estimators constructed from residuals

full rationale

The derivation constructs guaranteed a posteriori estimators directly from the residuals of the discrete mixed variational formulation plus explicit Robin boundary terms. Reliability follows from the assumed well-posedness of the continuous multigroup SPN problem and standard mixed-FEM approximation properties on simplicial/Cartesian meshes; these are external inputs rather than quantities defined or fitted inside the estimator itself. No self-definitional loops, fitted parameters renamed as predictions, or load-bearing self-citations appear in the central chain. The DD+L2-jumps extension likewise uses the same residual structure. The only minor deduction is the unverified invocation of continuous well-posedness for the coupled multigroup case, which does not create a circular reduction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard properties of mixed finite element spaces and the well-posedness of the SPN problem with Robin conditions; no free parameters, ad-hoc axioms, or invented entities are introduced in the abstract.

axioms (2)
  • domain assumption Well-posedness of the continuous multigroup SPN problem with Robin boundary conditions
    Invoked to guarantee existence of the exact solution against which the discrete error is measured.
  • standard math Approximation properties of mixed finite element spaces on simplicial and Cartesian meshes
    Used to obtain local efficiency and reliability constants independent of mesh size.

pith-pipeline@v0.9.0 · 5471 in / 1352 out tokens · 49471 ms · 2026-05-13T18:48:35.147471+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

35 extracted references · 35 canonical work pages

  1. [1]

    Arbogast and Z

    T. Arbogast and Z. Chen. On the implementation of mixed methods as nonconforming methods for second-order elliptic problems.Math. Comp., 64(211):943–972, 1995

    work page 1995

  2. [2]

    Baudron and J.-J

    A.-M. Baudron and J.-J. Lautard.SPN core calculations in the APOLLO3 System. In Mathematics and Computational Methods Applied to Nuclear Science and Engineering (M&C 2011). Latin American Section (LAS) / American Nuclear Society (ANS), 2011

    work page 2011

  3. [3]

    Boffi, F

    D. Boffi, F. Brezzi, and M. Fortin.Mixed and hybrid finite element methods and applications. Springer-Verlag, 2013

    work page 2013

  4. [4]

    Carstensen

    C. Carstensen. A posteriori error estimate for the mixed finite element method. Math. Comp., 66(218):465–476, 1997

    work page 1997

  5. [5]

    Ciarlet.The finite element method for elliptic problems, volume 40 of Classics in Applied Mathmetics

    P.G. Ciarlet.The finite element method for elliptic problems, volume 40 of Classics in Applied Mathmetics. SIAM, 2002

    work page 2002

  6. [6]

    Ciarlet Jr

    P. Ciarlet Jr. T-coercivity: a practical tool for the study of variational formulations in Hilbert spaces. Technical Report hal-05421231v1, HAL, 2025

    work page 2025

  7. [7]

    Ciarlet, Jr., M.-H

    P. Ciarlet, Jr., M.-H. Do, M. Gervais, and F. Madiot. A posteriori error estimates for the DD+L2 jumps method on the neutron diffusion equations. Computers & Mathematics with Applications, 195:349–365, 2025

    work page 2025

  8. [8]

    Ciarlet, Jr., M

    P. Ciarlet, Jr., M. H. Do, and F. Madiot. A posteriori error estimates for mixed finite element discretizations of the Neutron Diffusion equa- tions. ESAIM: Mathematical Modelling and Numerical Analysis, 57(1):1– 27, 2023

    work page 2023

  9. [9]

    Ciarlet, Jr., M.-H

    P. Ciarlet, Jr., M.-H. Do, and F. Madiot. Adaptive mesh refinement on Cartesian meshes applied to the mixed finite element discretization of the 32 multigroup neutron diffusion equations. In Error Control, Adaptive Dis- cretizations, and Applications, Part 3, volume 60 ofAdvances in Applied Mechanics, pages 201–234. Elsevier, 2025

    work page 2025

  10. [10]

    Ciarlet, Jr., L

    P. Ciarlet, Jr., L. Giret, E. Jamelot, and F. D. Kpadonou. Numerical analysis of the mixed finite element method for the neutron diffusion eigen- problem with heterogeneous coefficients. ESAIM: Math. Modell. Numer. Anal., 52:2003–2035, 2018

    work page 2003

  11. [11]

    Ciarlet, Jr., E

    P. Ciarlet, Jr., E. Jamelot, and F. D. Kpadonou. Domain decomposition methods for the diffusion equation with low-regularity solution.Comput. Math. Applic., 74:2369–2384, 2017

    work page 2017

  12. [12]

    J. J. Duderstadt and L. J. Hamilton.Nuclear reactor analysis. John Wiley & Sons, Inc., 1976

    work page 1976

  13. [13]

    Ern and J.-L

    A. Ern and J.-L. Guermond.Finite Elements I, volume 72 ofTexts Appl. Math. Springer, 2021

    work page 2021

  14. [14]

    Ern and J.-L

    A. Ern and J.-L. Guermond.Finite Elements II, volume 73 ofTexts Appl. Math. Springer, 2021

    work page 2021

  15. [15]

    E. M. Gelbard. Application of spherical harmonics method to reactor prob- lems. Bettis Atomic Power Laboratory, West Mifflin, PA, Technical Report No. WAPD-BT-20, 1960

    work page 1960

  16. [16]

    M. Gervais. Estimateurs a posteriori pour une méthode de décomposition de domaine non-conforme pour l’équation de diffusion neutronique. PhD thesis, Université Paris Saclay, 2026

    work page 2026

  17. [17]

    Giret.Non-conforming domain decomposition for the multigroup neutron SPN equations

    L. Giret.Non-conforming domain decomposition for the multigroup neutron SPN equations. PhD thesis, Université Paris Saclay, 2018

    work page 2018

  18. [18]

    Jamelot, A.-M

    E. Jamelot, A.-M. Baudron, and J.-J. Lautard. Domain decomposition for the SPN solver MINOS. Transport Theory and Statistical Physics, 41(7):495–512, 2012

    work page 2012

  19. [19]

    Jamelot and P

    E. Jamelot and P. Ciarlet, Jr. Fast non-overlapping schwarz domain de- composition methods for solving the neutron diffusion equation.Journal of Computational Physics, 241:445–463, 2013

    work page 2013

  20. [20]

    Jamelot and F

    E. Jamelot and F. Madiot. Numerical analysis of the neutron multigroup SPN equations. Comptes Rendus. Mathématique, 359(5), 2021

    work page 2021

  21. [21]

    Jamelot and P

    E. Jamelot and P. Ciarlet, Jr. Fast non-overlapping Schwarz domain de- compositionmethodsforsolvingtheneutrondiffusionequation. J. Comput. Phys., 241:445–463, 2013

    work page 2013

  22. [22]

    Könnö, D

    J. Könnö, D. Schötzau, and R. Stenberg. Mixed finite element methods for problems with Robin boundary conditions. SIAM Journal on Numerical Analysis, 49(1):285 – 308, 2011. 33

    work page 2011

  23. [23]

    M. G. Larson and A. Målqvist. A posteriori error estimates for mixed finite element approximations of elliptic problems. Numer. Math., 108(3):487– 500, 2008

    work page 2008

  24. [24]

    J. J. Lee. A posteriori error estimates of Darcy flows with Robin-type jump interface conditions. Computers & Mathematics with Applications, 174:417–430, 2024

    work page 2024

  25. [25]

    Lovadina and R

    C. Lovadina and R. Stenberg. Energy norm a posteriori error estimates for mixed finite element methods.Math. Comp., 75(256):1659–1674, 2006

    work page 2006

  26. [26]

    Marchuk and V.I

    G. Marchuk and V.I. Lebedev.Numerical methods in the theory of neutron transport. Harwood Academic Pub, New York, NY, 1986

    work page 1986

  27. [27]

    J.-C. Nédélec. Mixed finite elements inR3. Numer. Math., 35(3):315–341, 1980

    work page 1980

  28. [28]

    P. Oswald. On a BPX-preconditioner for P1 elements. Computing, 51(2):125–133, 1993

    work page 1993

  29. [29]

    G. V. Pencheva, M. Vohralík, M. F. Wheeler, and T. Wildey. Robust a Posteriori Error Control and Adaptivity for Multiscale, Multinumerics, and Mortar Coupling. SIAM Journal on Numerical Analysis, 51(1):526–554, 2013

    work page 2013

  30. [30]

    Raviart and J.-M

    P.-A. Raviart and J.-M. Thomas. A mixed finite element method for second order elliptic problems. In Mathematical aspects of finite element meth- ods, volume 606 ofLecture Notes in Mathematics, pages 292–315. Springer, 1977

    work page 1977

  31. [31]

    Takeda and H

    T. Takeda and H. Ikeda. 3-D neutron transport benchmarks.Journal of Nuclear Science and Technology, 28(7):656–669, 1991

    work page 1991

  32. [32]

    Vohralík

    M. Vohralík. A posteriori error estimates for lowest-order mixed finite element discretizations of convection-diffusion-reaction equations. SIAM J. Numer. Anal., 45(4):1570–1599, 2007

    work page 2007

  33. [33]

    Vohralík

    M. Vohralík. Unified primal formulation-based a priori and a posteriori erroranalysisofmixedfiniteelementmethods. Math. Comp., 79(272):2001– 2032, 2010

    work page 2001

  34. [34]

    M. F. Wheeler and I. Yotov. A posteriori error estimates for the mortar mixed finite element method. SIAM J. Numer. Anal., 43(3):1021–1042, 2005

    work page 2005

  35. [35]

    Wohlmuth and R

    B. Wohlmuth and R. Hoppe. A comparison of a posteriori error estimators for mixed finite element discretizations by Raviart-Thomas elements.Math. Comp., 68(228):1347–1378, 1999. 34 A A model with mixed boundary conditions In this section, we describe how a posteriori estimation theory can be extended to the case where mixed boundary conditions are imposed...

    work page 1999