pith. sign in

arxiv: 2604.21752 · v1 · submitted 2026-04-23 · 🧮 math.NA · cs.NA

Stable and asymptotic preserving space-time discretizations of a linear kinetic transport equation in diffusive scaling

Anita Gjesteland , Sigrun Ortleb , Salim Elghawi , David C. Del Rey Fern\'andez This is my paper

Pith reviewed 2026-05-09 21:02 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords linear kinetic transportmicro-macro decompositionsummation-by-partsenergy stabilityasymptotic preservationdiffusive scalingsimultaneous approximation terms
0
0 comments X

The pith

Space-time discretization of linear kinetic transport is unconditionally energy stable and asymptotically preserving

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

The paper develops a tensor-product space-time discretization for a linear kinetic transport equation in diffusive scaling that is unconditionally energy stable and preserves the asymptotic limit as the scaling parameter vanishes. It uses micro-macro decomposition to split the particle distribution into equilibrium and non-equilibrium parts, then proves stability and preservation properties hold for any spatial and temporal schemes that satisfy the summation-by-parts property. A new boundary treatment based on simultaneous approximation terms enforces Dirichlet conditions while retaining energy stability. Readers would care because transport models like radiative transfer become stiff in the diffusive regime, and this approach guarantees reliable behavior on coarse grids without artificial instabilities or loss of the physical limit.

Core claim

For the micro-macro decomposed linear kinetic transport equation, general spatial and temporal discretizations satisfying the summation-by-parts property, together with a simultaneous approximation term treatment for Dirichlet boundaries, produce a fully discrete scheme that is unconditionally energy stable and asymptotically preserving in the diffusive scaling limit.

What carries the argument

The micro-macro decomposition of the distribution function into equilibrium and non-equilibrium components, paired with summation-by-parts operators in space and time plus simultaneous approximation terms for boundary conditions.

Load-bearing premise

The spatial and temporal discretizations must satisfy the summation-by-parts property and the micro-macro decomposition must remain valid for the linear kinetic transport equation under diffusive scaling.

What would settle it

Numerical experiments on a non-summation-by-parts discretization that exhibit energy growth or fail to approach the diffusion solution as the scaling parameter goes to zero would disprove the stability and asymptotic preservation claims.

Figures

Figures reproduced from arXiv: 2604.21752 by Anita Gjesteland, David C. Del Rey Fern\'andez, Salim Elghawi, Sigrun Ortleb.

Figure 1
Figure 1. Figure 1: Example grid with one element in the spatial direction and one slab in the temporal direction. [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: figure 2. For the sake of clarity in the presentation of the scheme and the readability of the proofs [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Example grid with multiple elements in space and two slabs in time. [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Example with variable scattering frequency and homogeneous Dirichlet boundary conditions. [PITH_FULL_IMAGE:figures/full_fig_p017_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Numerical solutions at times t “ 0.1, 0.4, 1.0, 1.6, 4.0 for inhomogeneous Dirichlet boundary conditions in the kinetic regime ε “ 1. Multi-element spatial domain with K “ 10 cells and N “ 2 (left) vs N “ 3 (right) nodes per cell. 17 [PITH_FULL_IMAGE:figures/full_fig_p017_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Numerical solutions at times t “ 0.1, 0.4, 1.0, 1.6, 4.0 for inhomogeneous Dirichlet boundary conditions in the diffusive regime ε “ 10´8 . Multi-element spatial domain with K “ 10 cells and N “ 2 (left) vs N “ 3 (right) nodes per cell. 5 Conclusions and future work We developed an unconditionally energy-stable tensor-product space-time discretization framework based on the combination of spatial and tempo… view at source ↗
read the original abstract

We develop an unconditionally energy-stable tensor-product space-time discretization framework for the solution of a linear kinetic transport equation in one space dimension. The kinetic equation is a simplified model of radiative transfer formulated as a hyperbolic balance law in diffusive scaling for a particle distribution function of the independent variables space, time and velocity. Our numerical discretization is based on the well-known technique of micro-macro decomposition which results in a system of balance laws for equilibrium and non-equilibrium quantities and facilitates preservation of the asymptotic limit for vanishing scaling parameters at the discrete level. We prove fully discrete stability and asymptotic preservation for general spatial and temporal discretizations having the summation-by-parts property. A new provably energy-stable Dirichlet boundary treatment for the micro-macro decomposed system is developed based on the introduction of simultaneous approximation terms. Numerical results show convergence for smooth problems and demonstrate energy stability of the proposed boundary treatment.

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

2 major / 2 minor

Summary. The paper develops an unconditionally energy-stable tensor-product space-time discretization for a one-dimensional linear kinetic transport equation in diffusive scaling. It employs micro-macro decomposition to obtain a system of balance laws, proves fully discrete energy stability and asymptotic preservation for general summation-by-parts (SBP) spatial and temporal discretizations, introduces a new simultaneous approximation term (SAT) Dirichlet boundary treatment that is provably energy-stable, and reports numerical results on convergence for smooth problems and energy stability of the boundary treatment.

Significance. If the proofs of stability and asymptotic preservation hold, including for the new SAT boundary treatment, the work supplies a general framework for high-order, stable, and limit-preserving discretizations of kinetic equations that applies to arbitrary SBP operators. This is valuable for multiscale radiative transfer simulations, where preserving the diffusion limit at the discrete level on bounded domains is essential. The tensor-product space-time approach combined with SBP and SAT for boundaries represents a technical advance over existing methods.

major comments (2)
  1. [Boundary treatment and AP sections (e.g., the SAT formulation and the fully discrete AP proof)] The proof of asymptotic preservation (AP) must explicitly include the SAT boundary terms to confirm consistency with the diffusion limit. It is not clear from the abstract or the stated claims whether the SAT penalty parameters remain consistent (i.e., enforce the correct boundary conditions on the equilibrium variable) as the scaling parameter vanishes; an O(1) inconsistency in the boundary treatment would render the scheme stable but not AP on bounded domains.
  2. [Energy stability proof for the boundary treatment] In the energy stability analysis for the micro-macro system with SAT boundaries, the estimates should be re-derived in the diffusive limit to verify that the boundary contributions do not prevent the scheme from recovering the correct diffusion equation and boundary conditions at the discrete level.
minor comments (2)
  1. [Abstract and §1] Clarify in the abstract and introduction whether 'unconditionally energy-stable' holds for arbitrary time steps or requires a CFL-type restriction independent of the scaling parameter.
  2. [Numerical results] The numerical results section should include a test case that specifically checks asymptotic preservation (e.g., convergence to the diffusion solution as ε→0) with the new SAT boundaries active.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments regarding the boundary treatment and asymptotic preservation analysis. We address each major comment point by point below. We agree that the proofs require explicit extension to the SAT terms and will revise the manuscript accordingly to strengthen the presentation.

read point-by-point responses

  1. Referee: [Boundary treatment and AP sections (e.g., the SAT formulation and the fully discrete AP proof)] The proof of asymptotic preservation (AP) must explicitly include the SAT boundary terms to confirm consistency with the diffusion limit. It is not clear from the abstract or the stated claims whether the SAT penalty parameters remain consistent (i.e., enforce the correct boundary conditions on the equilibrium variable) as the scaling parameter vanishes; an O(1) inconsistency in the boundary treatment would render the scheme stable but not AP on bounded domains.

    Authors: We agree that the AP proof must explicitly incorporate the SAT boundary terms for a complete demonstration on bounded domains. The SAT penalties in our formulation are constructed to weakly impose the Dirichlet conditions on the equilibrium variable in a manner consistent with the micro-macro decomposition. However, the fully discrete AP analysis in the manuscript primarily addresses the interior discretization and volume terms. We will revise Section 4 (and related boundary sections) to include the boundary SAT contributions in the limit process as ε → 0. This extension will verify that the penalty parameters scale appropriately to enforce the correct boundary conditions of the limiting diffusion equation, ruling out any O(1) inconsistency. We believe this will fully confirm the AP property for the complete scheme including boundaries. revision: yes

  2. Referee: [Energy stability proof for the boundary treatment] In the energy stability analysis for the micro-macro system with SAT boundaries, the estimates should be re-derived in the diffusive limit to verify that the boundary contributions do not prevent the scheme from recovering the correct diffusion equation and boundary conditions at the discrete level.

    Authors: The energy stability proof for the SAT boundary treatment is derived in a general form that holds for arbitrary positive scaling parameters. To address this point, we will add an explicit limiting analysis in the revised manuscript: we take the diffusive limit (ε → 0) directly on the energy estimates, including all boundary SAT terms. This will show that the boundary contributions reduce consistently to those of the discrete diffusion equation with the appropriate boundary conditions, without preventing recovery of the limit. Since unconditional stability holds for all ε > 0, the limit preserves the stability property, but we will provide the detailed verification as requested. revision: yes

Circularity Check

0 steps flagged

No circularity: proofs build on standard SBP properties and micro-macro decomposition from prior literature

full rationale

The paper establishes fully discrete stability and asymptotic preservation via mathematical proofs that invoke the summation-by-parts (SBP) property of the underlying discretizations and the standard micro-macro decomposition of the linear kinetic transport equation. These are external, well-established techniques cited from prior independent literature rather than derived within the paper or fitted to its own outputs. The new simultaneous approximation term (SAT) Dirichlet boundary treatment is analyzed directly for energy stability, with no reduction of the asymptotic-preservation claim to a self-definition, renamed fit, or self-citation chain. The derivation chain remains self-contained and non-tautological against external benchmarks such as the diffusive limit and energy estimates.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the summation-by-parts property of the chosen discretizations and the applicability of micro-macro decomposition to the given linear equation; no free parameters or new entities are introduced.

axioms (2)
  • domain assumption Spatial and temporal discretizations possess the summation-by-parts property
    Invoked to prove fully discrete energy stability and asymptotic preservation.
  • domain assumption Micro-macro decomposition is valid for the linear kinetic transport equation in diffusive scaling
    Forms the basis for splitting into equilibrium and non-equilibrium components.

pith-pipeline@v0.9.0 · 5462 in / 1296 out tokens · 27047 ms · 2026-05-09T21:02:38.002955+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

32 extracted references · 32 canonical work pages

  1. [1]

    Bennoune, M

    M. Bennoune, M. Lemou, and L. Mieussens. Uniformly stable numerical schemes for the Boltzmann equation preserving the compressible Navier–Stokes asymptotics.Journal of Com- putational Physics, 227(8):3781–3803, 2008

    work page 2008

  2. [2]

    Boscarino, L

    S. Boscarino, L. Pareschi, and G. Russo. Implicit-explicit Runge–Kutta schemes for hyperbolic systems and kinetic equations in the diffusion limit.SIAM Journal on Scientific Computing, 35(1):A22–A51, 2013

    work page 2013

  3. [3]

    Bourgat, P

    J. Bourgat, P. Tallec, B. Perthame, and Y. Qiu. Coupling Boltzmann and Euler equations without overlapping.Contemp. Math., 157, 1994

    work page 1994

  4. [4]

    M. H. Carpenter, D. Gottlieb, and S. Abarbanel. Time-stable boundary conditions for finite- difference schemes solving hyperbolic systems: Methodology and application to high-order compact schemes.Journal of Computational Physics, 111(2):220–236, 1994

    work page 1994

  5. [5]

    Coron and B

    F. Coron and B. Perthame. Numerical passage from kinetic to fluid equations.SIAM Journal on Numerical Analysis, 28(1):26–42, 1991

    work page 1991

  6. [6]

    Crean, J

    J. Crean, J. E. Hicken, D. C. D. R. Fern´ andez, D. W. Zingg, and M. H. Carpenter. Entropy- stable summation-by-parts discretizations of the euler equations on general curved elements. Journal of Computational Physics, 356:410–438, 2018

    work page 2018

  7. [7]

    Degond and S

    P. Degond and S. Jin. A smooth transition model between kinetic and diffusion equations. SIAM Journal on Numerical Analysis, 42(6):2671–2687, 2005. 22

    work page 2005

  8. [8]

    D. C. Del Rey Fern´ andez, P. D. Boom, M. H. Carpenter, and D. W. Zingg. Extension of tensor-product generalized and dense-norm summation-by-parts operators to curvilinear coordinates.Journal of Scientific Computing, 80(4):1957–1996, 2019

    work page 1957

  9. [9]

    D. C. Del Rey Fern´ andez, J. E. Hicken, and D. W. Zingg. Review of summation-by-parts oper- ators with simultaneous approximation terms for the numerical solution of partial differential equations.Computer & Fluids, 95:171–196, 2014

    work page 2014

  10. [10]

    Dimarco and L

    G. Dimarco and L. Pareschi. Asymptotic preserving implicit-explicit runge–kutta methods for nonlinear kinetic equations.SIAM Journal on Numerical Analysis, 51(2):1064–1087, 2013

    work page 2013

  11. [11]

    I. M. Gamba, S. Jin, and L. Liu. Micro-macro decomposition based asymptotic-preserving numerical schemes and numerical moments conservation for collisional nonlinear kinetic equa- tions.Journal of Computational Physics, 382:264–290, 2019

    work page 2019

  12. [12]

    J. E. Hicken, D. C. Del Rey Fern´ andez, and D. W. Zingg. Multidimensional summation- by-parts operators: General theory and application to simplex elements.SIAM Journal on Scientific Computing, 38(4):A1935–A1958, 2016

    work page 2016

  13. [13]

    J. E. Hicken, D. C. Del Rey Fern´ andez, and D. W. Zingg. Multidimensional summation- by-parts operators: general theory and application to simplex elements.SIAM Journal on Scientific Computing, 38(4):A1935–A1958, 2016

    work page 2016

  14. [14]

    J. Jang, F. Li, J.-M. Qiu, and T. Xiong. Analysis of asymptotic preserving DG-IMEX schemes for linear kinetic transport equations in a diffusive scaling.SIAM Journal on Numerical Analysis, 52(4):2048–2072, 2014

    work page 2048

  15. [15]

    J. Jang, F. Li, J.-M. Qiu, and T. Xiong. High order asymptotic preserving DG-IMEX schemes for discrete-velocity kinetic equations in a diffusive scaling.Journal of Computational Physics, 281:199–224, 2015

    work page 2015

  16. [16]

    S. Jin, L. Pareschi, and G. Toscani. Diffusive relaxation schemes for multiscale discrete- velocity kinetic equations.SIAM Journal on Numerical Analysis, 35(6):2405–2439, 1998

    work page 1998

  17. [17]

    Kreiss and G

    H.-O. Kreiss and G. Scherer. Finite element and finite difference methods for hyperbolic partial differential equations. In C. d. Boor, editor,Mathematical Aspects of Finite Elements in Partial Differential Equations, pages 195–212. Academic Press, 1974

    work page 1974

  18. [18]

    Lemou and L

    M. Lemou and L. Mieussens. A new asymptotic preserving scheme based on micro-macro formulation for linear kinetic equations in the diffusion limit.SIAM Journal on Scientific Computing, 31(1):334–368, 2008

    work page 2008

  19. [19]

    Lundquist and J

    T. Lundquist and J. Nordstr¨ om. The SBP-SAT technique for initial value problems.Journal of Computational Physics, 270:86–194, 2014

    work page 2014

  20. [20]

    Mattsson

    K. Mattsson. Diagonal-norm upwind sbp operators.Journal of Computational Physics, 335:283–310, 2017

    work page 2017

  21. [21]

    Naldi and L

    G. Naldi and L. Pareschi. Numerical schemes for kinetic equations in diffusive regimes.Applied Mathematics Letters, 11(2):29–35, 1998

    work page 1998

  22. [22]

    Naldi and L

    G. Naldi and L. Pareschi. Numerical schemes for hyperbolic systems of conservation laws with stiff diffusive relaxation.SIAM Journal on Numerical Analysis, 37(4):1246–1270, 2000

    work page 2000

  23. [23]

    Nordstr¨ om, J

    J. Nordstr¨ om, J. Gong, E. van der Weide, and M. Sv¨ ard. A stable and conservative high order multi-block method for the compressible navier-stokes equations.Journal of Computational Physics, 228(24):9020–9035, 2009

    work page 2009

  24. [24]

    Nordstr¨ om and T

    J. Nordstr¨ om and T. Lundquist. Summation-by-parts in time.Journal of Computational Physics, 251:487–499, 2013

    work page 2013

  25. [25]

    P. Olsson. Summation by parts, projections, and stability. i.Mathematics of Computation, 64(211):1035–1065, 1995

    work page 1995

  26. [26]

    P. Olsson. Summation by parts, projections, and stability. ii.Mathematics of Computation, 64(212):1473–1493, 1995. 23

    work page 1995

  27. [27]

    S. Ortleb. On unconditional stability of upwind generalized summation-by-parts schemes for kinetic equations in the diffusion limit.PAMM, 24(4):e202400119, 2024

    work page 2024

  28. [28]

    Z. Peng, Y. Cheng, J.-M. Qiu, and F. Li. Stability-enhanced AP IMEX1-LDG method: Energy-based stability and rigorous AP property.SIAM Journal on Numerical Analysis, 59(2):925–954, 2021

    work page 2021

  29. [29]

    Peng and F

    Z. Peng and F. Li. Asymptotic preserving imex-dg-s schemes for linear kinetic transport equations based on schur complement.SIAM Journal on Scientific Computing, 43(2):A1194– A1220, 2021

    work page 2021

  30. [30]

    B. Strand. Summation by parts for finite difference approximations for d/dx.Journal of Computational Physics, 110(1):47–67, 1994

    work page 1994

  31. [31]

    Sv¨ ard and J

    M. Sv¨ ard and J. Nordstr¨ om. Review of summation-by-parts schemes for initial-boundary- value problems.Journal of Computational Physics, 268:17–38, 2014

    work page 2014

  32. [32]

    Tiwari and A

    S. Tiwari and A. Klar. An adaptive domain decomposition procedure for Boltzmann and Euler equations.Journal of Computational and Applied Mathematics, 90(2):223–237, 1998. 24

    work page 1998