pith. sign in

arxiv: 2602.17936 · v2 · submitted 2026-02-20 · 🧮 math.NA · cs.NA

Optimal error estimate of an isoparametric upwind discontinuous Galerkin method for radiation transport equation on curved domains

Changhui Yao , Yunpan Ma , Lingxiao Li This is my paper

Pith reviewed 2026-05-15 21:12 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords discontinuous Galerkinisoparametric mappingradiation transport equationcurved domainsoptimal error estimateupwind schemehyperbolic PDEgeometric approximation
0
0 comments X

The pith

An isoparametric upwind discontinuous Galerkin scheme achieves optimal convergence in the DG norm for the radiation transport equation on curved domains.

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

The paper constructs an isoparametric upwind discontinuous Galerkin method that approximates both the solution and the curved geometry of the domain using high-order mappings. It addresses extra difficulties from non-affine transformations by introducing an auxiliary operator that restores continuity of the bilinear form with respect to the DG norm. Precise bounds are derived for the geometric error on the inflow boundary and for the mismatch between discrete and continuous normal vectors. These steps together prove that the scheme retains its optimal convergence rate in the DG norm. The result matters for accurate high-order modeling of radiation transport where mesh boundaries must follow physical curves without sacrificing accuracy or requiring excessive refinement.

Core claim

With the aid of an isoparametric auxiliary operator, the bilinear form remains continuous in the DG norm when its first argument is the isoparametric projection error. The geometric approximation error on the inflow boundary of the original domain is bounded at the optimal order, and the difference between discrete and continuous normal vectors is shown to be sufficiently small. These estimates together establish that the isoparametric upwind discontinuous Galerkin method for the radiation transport equation converges at the optimal rate in the DG norm on bounded domains whose boundary is piecewise C^{k+1} smooth.

What carries the argument

The isoparametric mapping that approximates the curved domain, paired with an auxiliary operator that restores continuity of the discrete bilinear form for the projection error.

If this is right

  • Optimal rates hold in the DG norm, controlling both solution values and inter-element fluxes simultaneously.
  • The method applies directly to high-order simulations on realistic curved geometries without order reduction from geometric approximation.
  • Two- and three-dimensional tests confirm that the theoretical rates are attained in practice.
  • Geometric errors on the inflow boundary remain controlled at the same order as the interior approximation error.

Where Pith is reading between the lines

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

  • The same auxiliary-operator technique may extend optimal-order proofs to other first-order hyperbolic equations on curved domains.
  • Mesh-generation requirements near curved interfaces can be relaxed when boundary smoothness matches the chosen polynomial degree.
  • Radiation-transport codes in engineering or medical physics could reduce the number of elements needed near curved boundaries while keeping global accuracy.

Load-bearing premise

The physical domain boundary must be piecewise C^{k+1} smooth so that the isoparametric mapping supplies geometric approximation order matching the polynomial degree.

What would settle it

Numerical experiments on a domain whose boundary has smoothness strictly less than C^{k+1} that exhibit convergence rates strictly below the optimal order predicted by the analysis.

Figures

Figures reproduced from arXiv: 2602.17936 by Changhui Yao, Lingxiao Li, Yunpan Ma.

Figure 1
Figure 1. Figure 1: Left: A curved tetrahedron K. Right: A sphere domain divided by curved tetrahedra Th. Moreover, we denote the boundary of Dh by Γh and the inflow boundary by Γ− h , namely Γh = ∂Dh, Γ − h = {y ∈ Γh : Ω · nh(y) < 0}. Here nh(y) is the unit outward normal vector of Γh at point y. We also define the inflow (outflow) boundary for each element K ∈ Th : ∂−(+)K = {x ∈ ∂K : Ω · nh(x) < (>)0}. Combining the assumpt… view at source ↗
Figure 2
Figure 2. Figure 2: Isoparametric mapping and numerical results in 2D. [PITH_FULL_IMAGE:figures/full_fig_p013_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Mesh and domain evolution through isoparametric mapping in 3D. [PITH_FULL_IMAGE:figures/full_fig_p014_3.png] view at source ↗
read the original abstract

In recent years, high-order finite element methods on high-order meshes have attracted considerable attention. This work investigates the isoparametric upwind discontinuous Galerkin method for the radiation transport equation on a bounded domain with a piecewise $C^{k+1}$ smooth curved boundary. We use the isoparametric mapping to approximate the curved domain and construct a curved upwind discontinuous Galerkin scheme. The first-order hyperbolic nature and the complexity introduced by non-affine transformation, lead to additional difficulties for geometric approximation, numerical stability and the optimal error estimate. To address these issues, with the help of an isoparametric auxiliary operator, we first prove that the bilinear form is continuous with respect to the DG norm when its first argument is the isoparametric projection error. Then the geometric approximation error of inflow boundary of original domain is precisely estimated. The error order between discrete normal vectors and the continuous ones are also proven. Finally, the rigorous analysis yields an optimal convergence rate in the DG norm. Two- and three-dimensional numerical tests are conducted to support the theoretical results.

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

0 major / 2 minor

Summary. The manuscript develops an isoparametric upwind discontinuous Galerkin method for the radiation transport equation on bounded domains with piecewise C^{k+1} smooth curved boundaries. It approximates the domain via isoparametric mappings, constructs the corresponding curved upwind DG scheme, proves continuity of the bilinear form on the projection error via an auxiliary operator, bounds geometric approximation errors on the inflow boundary, establishes the order of consistency between discrete and continuous normal vectors, and derives an optimal convergence rate in the DG norm, with supporting 2D and 3D numerical tests.

Significance. If the central claims hold, the work is significant for numerical analysis of hyperbolic problems on complex geometries: it supplies a rigorous optimal-error theory for high-order isoparametric DG schemes that accounts for non-affine transformations and geometric approximation, a setting that arises frequently in radiation transport but has lacked such estimates.

minor comments (2)
  1. Abstract: the phrase 'optimal convergence rate' is used without stating the precise order (e.g., O(h^{k+1/2}) in the DG norm); adding the explicit rate would improve immediate readability.
  2. Section describing the numerical tests: the presentation of mesh families, refinement strategy, and exact error-measurement procedure is terse; expanding these details would make the computational validation easier to reproduce.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful summary of our manuscript and for the positive assessment of its significance for high-order DG methods on curved domains. We appreciate the recommendation of minor revision. No specific major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The derivation proceeds by constructing an isoparametric auxiliary operator, proving continuity of the bilinear form on the projection error in the DG norm, bounding geometric approximation errors on the inflow boundary at the required order, and establishing consistency between discrete and continuous normal vectors. These steps close the argument for the optimal DG-norm convergence rate under the stated piecewise C^{k+1} boundary regularity. No equation reduces the target error bound to a fitted quantity defined by the result itself, no self-citation chain is load-bearing for the central claim, and the proof relies on standard DG stability plus geometric estimates that remain independent of the final rate. The analysis is therefore self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The analysis rests on standard finite-element assumptions for hyperbolic problems plus domain smoothness; no free parameters or new entities are introduced.

axioms (1)
  • domain assumption The domain boundary is piecewise C^{k+1} smooth
    Invoked to guarantee that the isoparametric mapping achieves the geometric approximation order required for the optimal DG error bound.

pith-pipeline@v0.9.0 · 5498 in / 1218 out tokens · 23925 ms · 2026-05-15T21:12:20.846228+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

  1. [1]

    R. A. Adams and J. J. F. Fournier. Sobolev spaces (2nd ed.).Academic Press, (2003)

    work page 2003

  2. [2]

    D. N. Arnold and S. W. Walker. The Hellan–Herrmann–Johnson method with curved elements.SIAM J. Numer. Anal., 58(5) (2020), pp. 2829–2855

    work page 2020

  3. [3]

    Aylwin and C

    R. Aylwin and C. Jerez-Hanckes. Finite-element domain approximation for Maxwell varia- tional problems on curved domains.SIAM J. Numer. Anal., 61(3) (2023), pp. 1139–1171

    work page 2023

  4. [4]

    M. A. Badri, P. Jolivet, B. Rousseau and Y. Favennec. High performance computation of radiative transfer equation using the finite element method.J. Comput. Phys., 360 (2018), pp. 74–92

    work page 2018

  5. [5]

    Bertrand, S

    F. Bertrand, S. Munzenmaier and G. Starke. First-order system least squares on curved boundaries: Higher-order Raviart–Thomas elements.SIAM J. Numer. Anal., 52(6) (2014), pp. 3165–3180

    work page 2014

  6. [6]

    Bertrand and G

    F. Bertrand and G. Starke. Parametric Raviart–Thomas elements for mixed methods on domains with curved surfaces.SIAM J. Numer. Anal., 54(6) (2016), pp. 3648–3667

    work page 2016

  7. [7]

    J. H. Bramble and J. T. King. A robust finite element method for nonhomogeneous Dirich- let problems in domains with curved boundaries.Math. Comp., 63(207) (1994), pp. 1–17. 15

    work page 1994

  8. [8]

    S. C. Brenner and L. R. Scott. The mathematical theory of finite element methods (3nd ed.).New York: Springer, (2008)

    work page 2008

  9. [9]

    Brezzi, L

    F. Brezzi, L. D. Marini and E. S¨ uli. Discontinuous Galerkin methods for first-order hyper- bolic problems.Math. Models Methods Appl. Sci., 14(12) (2004), pp. 1893–1903

    work page 2004

  10. [10]

    Cheng and C

    J. Cheng and C. W. Shu. A third order conservative Lagrangian type scheme on curvilinear meshes for the compressible Euler equations.Commun. Comput. Phys., 4(5) (2008), pp. 1008–1024

    work page 2008

  11. [11]

    P. G. Ciarlet and P. A. Raviart. The combined effect of curved boundaries and numerical integration in isoparametric finite element methods.Academic Press, (1972)

    work page 1972

  12. [12]

    P. G. Ciarlet. The finite element method for elliptic problems. Amsterdam: North-Holland, 1978

    work page 1978

  13. [13]

    V. A. Dobrev, T. V. Kolev and R. N. Rieben. High-order curvilinear finite element methods for Lagrangian hydrodynamics.SIAM J. Sci. Comput., 34 (2012), pp. B606-B641

    work page 2012

  14. [14]

    J. Fu, J. Cheng, W. Li, T. Xiong and Y. Wang. An asymptotic-preserving IMEX PN method for the gray model of the radiative transfer equation.J. Sci. Comput., 103(2) (2025), pp. 71

    work page 2025

  15. [15]

    Garcke, R

    H. Garcke, R. N¨ urnberg, S. Praetorius and G. Zhang. Isoparametric finite element methods for mean curvature flow and surface diffusion.J. Comput. Phys., (2025), pp. 114248

    work page 2025

  16. [16]

    J. L. Guermond and G. Kanschat. Asymptotic analysis of upwind discontinuous Galerkin approximation of the radiative transport equation in the diffusive limit.SIAM J. Numer. Anal., 48(1) (2010), pp. 53–78

    work page 2010

  17. [17]

    J. S. Hesthaven and T. Warburton. Nodal discontinuous Galerkin methods: Algorithms, analysis, and applications.New York: Springer, (2008)

    work page 2008

  18. [18]

    Houston, C

    P. Houston, C. Schwab and E. S¨ uli. Stabilized hp-finite element methods for first-order hyperbolic problems.SIAM J. Numer. Anal., 37(5) (2000), pp. 1618–1643

    work page 2000

  19. [19]

    Houston, C

    P. Houston, C. Schwab and E. S¨ uli. Discontinuous hp-finite element methods for advection- diffusion-reaction problems.SIAM J. Numer. Anal., 39(6) (2002), pp. 2133–2163

    work page 2002

  20. [20]

    Johnson and J

    C. Johnson and J. Pitk¨ aranta. An analysis of the discontinuous Galerkin method for a scalar hyperbolic equation.Math. Comp., 46(173) (1986), pp. 1–26

    work page 1986

  21. [21]

    Kashiwabara

    T. Kashiwabara. Finite element analysis of a generalized Robin boundary value problem in curved domains based on the extension approach.IMA J. Numer. Anal., 45(5) (2025), pp. 2723–2755

    work page 2025

  22. [22]

    K. D. Lathrop. Use of discrete-ordinates methods for solution of photon transport problems. Nucl. Sci. Eng., 24(4) (1966), pp. 381–388

    work page 1966

  23. [23]

    J. M. Lee. Smooth manifolds.Introduction to smooth manifolds, pp. 1–29. New York: Springer, 2013

    work page 2013

  24. [24]

    Lehrenfeld and A

    C. Lehrenfeld and A. Reusken. Analysis of a high-order unfitted finite element method for elliptic interface problems.IMA J. Numer. Anal., 38(3) (2018), pp. 1351–1387

    work page 2018

  25. [25]

    Lehtikangas, T

    O. Lehtikangas, T. Tarvainen, A. D. Kim and S. R. Arridge. Finite element approximation of the radiative transport equation in a medium with piece-wise constant refractive index. J. Comput. Phys., 282 (2015), pp. 345–359

    work page 2015

  26. [26]

    M. Lenoir. Optimal isoparametric finite elements and error estimates for domains involving curved boundaries.SIAM J. Numer. Anal., 23(3) (1986), pp. 562–580

    work page 1986

  27. [27]

    Lesaint and P

    P. Lesaint and P. A. Raviart. On a finite element method for solving the neutron transport equation.Publ. S´ emin. Math. Inform. Rennes, (S4) (1974), pp. 1–40. 16

    work page 1974

  28. [28]

    B. Li, Y. Xia and Z. Yang. Optimal convergence of arbitrary Lagrangian-Eulerian iso- parametric finite element methods for parabolic equations in an evolving domain.IMA J. Numer. Anal., 43(1) (2023), pp. 501–534

    work page 2023

  29. [29]

    B. Li, W. Qiu, Y. Xie and W. Yu. Weak discrete maximum principle of isoparametric finite element methods in curvilinear polyhedra.Math. Comp., 93(345) (2024), pp. 1–34

    work page 2024

  30. [30]

    L. Li, H. Su, H. Zhang and W. Zheng. A divergence-free parametric finite element method for 3D Stokes equations on curved domains.arXiv preprintarXiv:2512.16216, (2025)

    work page arXiv 2025

  31. [31]

    M. F. Modest and S. Mazumder. Radiative heat transfer. Academic Press, 2021

    work page 2021

  32. [32]

    Olivier and T

    S. Olivier and T. S. Haut. High-order finite element second moment methods for linear transport.Nucl. Sci. Eng., 198(6) (2024), pp. 1179–1214

    work page 2024

  33. [33]

    Pazner and T

    W. Pazner and T. Haut. A short note on the accuracy of the discontinuous Galerkin method with reentrant faces.J. Comput. Phys., 443 (2021), pp. 110448

    work page 2021

  34. [34]

    T. E. Peterson. A note on the convergence of the discontinuous Galerkin method for a scalar hyperbolic equation.SIAM J. Numer. Anal., 28(1) (1991), pp. 133–140

    work page 1991

  35. [35]

    M. Tang, L. Wang and X. Zhang. Accurate front capturing asymptotic preserving scheme for nonlinear gray radiative transfer equation.SIAM J. Sci. Comput., 43(3) (2021), pp. B759–B783

    work page 2021

  36. [36]

    D. N. Woods and T. S. Palmer. RZ geometry discrete ordinates radiation transport us- ing higher-order finite element spatial discretizations on meshes with curved surfaces.J. Comput. Theor. Transport, 48(3) (2019), pp. 109–128

    work page 2019

  37. [37]

    L. B. Zhang. A parallel algorithm for adaptive local refinement of tetrahedral meshes using bisection.Numer. Math. Theory Methods Appl., 2(1) (2009), pp. 65–89

    work page 2009

  38. [38]

    Z. W. Zong, M. S. Cheng, Y. C. Yu and Z. M. Dai. A multithreaded parallel upwind sweep algorithm for the SN transport equations discretized with discontinuous finite elements. Nucl. Sci. Tech., 34(12) (2023), pp. 200. 17

    work page 2023