pith. sign in

arxiv: 2604.05299 · v1 · submitted 2026-04-07 · 🧮 math.NA · cs.NA

Multi-Scale Kinetic Simulation: Asymptotic Preserving IMEX-BDF-DG Schemes with Three Implicit-Explicit Partitionings

Kimberly Matsuda , Fengyan Li This is my paper

Pith reviewed 2026-05-10 20:01 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords asymptotic preserving methodsIMEX-BDF schemesdiscontinuous Galerkinkinetic transportmicro-macro decompositionradiative transfermulti-scale simulation
0
0 comments X

The pith

Three families of IMEX-BDF-DG schemes preserve the asymptotic diffusive limit of kinetic transport models at the discrete level.

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

The paper develops asymptotic preserving numerical methods for a linear kinetic transport model in diffusive scaling. It starts from the micro-macro decomposition and combines discontinuous Galerkin discretization in space with implicit-explicit BDF time integrators using three different partitionings of the terms. The resulting schemes are shown to recover the correct macroscopic diffusion behavior uniformly even on meshes that do not resolve the small scale. Systematic analysis and tests compare the three partitionings for stability, accuracy, computational cost, and preservation of the limit, and contrast them with related Runge-Kutta versions.

Core claim

The authors establish that the proposed IMEX-BDF-DG schemes with three implicit-explicit partitionings are asymptotic preserving for the simplified radiative transfer equation. This means the discrete solutions approach the solution of the diffusion equation as the scaling parameter tends to zero, without requiring the mesh size to shrink proportionally, while retaining high-order accuracy in space and time.

What carries the argument

The micro-macro decomposition of the kinetic distribution function, paired with three specific ways to split the resulting equations into implicit and explicit parts for BDF time stepping.

Load-bearing premise

The micro-macro decomposition remains valid and the chosen partitionings keep the asymptotic limit and stability without extra conditions on mesh size or scaling parameter.

What would settle it

A computation on a fixed under-resolved mesh where the discrete solution deviates from the diffusion-equation solution as the scaling parameter approaches zero would show the asymptotic preserving property fails.

Figures

Figures reproduced from arXiv: 2604.05299 by Fengyan Li, Kimberly Matsuda.

Figure 1
Figure 1. Figure 1: , with the scheme defined and analyzed in Appendix A. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 1.1
Figure 1.1. Figure 1.1: Telegraph equation: f(x, v = 1, T = 1) on [−π, π] computed by the first order upwind finite difference (also P 0 upwind discontinuous Galerkin) method in space with the backward Euler method in time for (1.1). Here σs = 1, σa = 0, G = 0 and the exact solution f(x, v, t) = 1 r e rt sin(x)+ εvert cos(x), r = −2 1+√ 1−4ϵ 2 , h = ∆t = π/40. The discrepancy in the computed and exact solutions for ε = 10−6 evi… view at source ↗
Figure 3.1
Figure 3.1. Figure 3.1: Stability regions of IMEX-BDFr-DGr with Strategy 1 (and 2) (top row) and Strategy 3 (bottom row) for r = 1, 2, 3; White: Stable, Black: Unstable; α = ε/(σs,mh), β = ∆t/(εh). From the stability plots, we further extract numerically the explicit formulas for the time step conditions, ∆t ≤ ∆t (k) CF Lr. These formulas are in similar ansatzes as in Theorem 3.1, except that with Strategy 3, we set the time st… view at source ↗
Figure 3.2
Figure 3.2. Figure 3.2: Stability curves of IMEX-BDF3-DG3, with the stable region below each curve. (a): [PITH_FULL_IMAGE:figures/full_fig_p014_3_2.png] view at source ↗
Figure 3.3
Figure 3.3. Figure 3.3: Maximum ∆t allowed for stability of IMEX-BDF1-DG1-S1 [PITH_FULL_IMAGE:figures/full_fig_p015_3_3.png] view at source ↗
Figure 3.4
Figure 3.4. Figure 3.4: is for IMEX-BDF1-DG1-S3. Similar observations as in [PITH_FULL_IMAGE:figures/full_fig_p015_3_4.png] view at source ↗
Figure 5.1
Figure 5.1. Figure 5.1: Initialization for IMEX-BDF3-DG3-Sk by methods of the same family of lower order temporal accuracy. Adaptive time-stepping region: by IMEX-BDF1-DG3-Sk (blue), IMEX-BDF2- DG3-Sk (purple); transitional region (red) and uniform time-stepping region (black) by IMEX￾BDF3-DG3-Sk (black). domain (i.e. the uniform time-stepping region). Our choice of the hyper-parameter ν ∈ (1, 2] is for the consideration of acc… view at source ↗
Figure 7.1
Figure 7.1. Figure 7.1: Example 1: Smooth example with constant material properties. Convergence study for [PITH_FULL_IMAGE:figures/full_fig_p022_7_1.png] view at source ↗
Figure 7.2
Figure 7.2. Figure 7.2: Example 1: Smooth example with constant material properties. To compare IMEX [PITH_FULL_IMAGE:figures/full_fig_p024_7_2.png] view at source ↗
Figure 7.3
Figure 7.3. Figure 7.3: Example 3: Two-material problem. The computed [PITH_FULL_IMAGE:figures/full_fig_p027_7_3.png] view at source ↗
Figure 7.4
Figure 7.4. Figure 7.4: Example 3: Two-material problem. Absolute error in [PITH_FULL_IMAGE:figures/full_fig_p028_7_4.png] view at source ↗
Figure 7.5
Figure 7.5. Figure 7.5: Example 4: Riemann problem for telegraph equation. [PITH_FULL_IMAGE:figures/full_fig_p029_7_5.png] view at source ↗
Figure 7.6
Figure 7.6. Figure 7.6: Example 4: Riemann problem for telegraph equation. [PITH_FULL_IMAGE:figures/full_fig_p030_7_6.png] view at source ↗
Figure 7.7
Figure 7.7. Figure 7.7: Example 4: Riemann problem for telegraph equation. [PITH_FULL_IMAGE:figures/full_fig_p031_7_7.png] view at source ↗
read the original abstract

Kinetic transport models are mesoscopic mathematical descriptions of the transport of particles as well as their interactions with the background media or among themselves, and they have wide applications in many areas of mathematical physics such as nuclear and biomedical engineering, rarefied gas dynamics, and plasma physics. They are often multi-scale, with different characteristics (e.g. hyperbolic, diffusive) depending on the material properties. As our continuing effort to design and analyze numerical methods for accurate and robust simulation of the multi-scale kinetic transport models, in this work, we consider a linear kinetic transport model, a simplified radiative transfer equation, in a diffusive scaling, and propose and analyze three families of asymptotic preserving (AP) methods. Numerical methods with the AP property, that is to preserve the asymptotic behavior of the models at the discrete level on under-resolved meshes, can work uniformly well to simulate multi-scale models across a wide range of scales. The proposed methods start from the micro-macro decomposition of the model, and involve discontinuous Galerkin (DG) methods in space, the discrete ordinates method (i.e. $S_N$ method) in velocity, and implicit-explicit (IMEX) BDF methods in time, with three different IMEX partitionings. A systematic study, both analytically and computationally, is presented regarding their difference in stability, accuracy, computational complexity and AP property. These methods, with multi-step time integrators, are also compared in terms of their accuracy and efficiency with the ones that only differ in using certain IMEX Runge-Kutta methods in time. Together with our previous developments, the present work further contributes to high order DG AP methods for multi-scale kinetic simulation, especially by utilizing the structure of the micro-macro decomposition of the models.

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 proposes and analyzes three families of asymptotic preserving (AP) IMEX-BDF-DG schemes for a linear kinetic transport model (simplified radiative transfer equation) in diffusive scaling. Starting from the micro-macro decomposition, the methods combine discontinuous Galerkin spatial discretization, discrete ordinates velocity discretization, and IMEX BDF time stepping with three distinct partitionings. It provides a systematic analytical and computational study of stability, accuracy, computational complexity, and the AP property, while comparing the BDF-based schemes to analogous IMEX Runge-Kutta methods.

Significance. If the AP property is verified to hold uniformly for all three partitionings on under-resolved meshes, the work strengthens the toolkit for high-order multi-scale kinetic simulations in applications such as radiative transfer and plasma physics. The explicit comparison between multi-step BDF and Runge-Kutta integrators, together with the focus on different IMEX partitionings, supplies practical guidance on trade-offs in stability and efficiency. The combination of micro-macro decomposition with DG and established IMEX techniques is a natural extension of prior developments in the field.

major comments (2)
  1. [analysis of AP property (likely §4)] The load-bearing claim is that each of the three IMEX-BDF partitionings, when combined with the micro-macro decomposition and DG discretization, yields an AP scheme whose ε→0 limit (fixed h, Δt) recovers a consistent discretization of the diffusion equation. The manuscript must explicitly pass to the limit inside the fully discrete equations for each partitioning separately, showing that explicit contributions vanish and the micro component relaxes to zero without residual inconsistent fluxes. This verification is required in the analysis section to support the uniform AP assertion.
  2. [numerical results section (likely §5)] Numerical tables demonstrating the AP property should report errors or solution profiles for successively smaller ε with fixed mesh size and time step, for all three partitionings. Without such data, it is difficult to confirm that the discrete limit is indeed the correct diffusion discretization rather than a perturbed or inconsistent scheme.
minor comments (2)
  1. [method formulation] Notation for the three IMEX partitionings should be introduced with explicit formulas for the implicit and explicit terms in each case to facilitate comparison of their stability and limit behavior.
  2. [throughout] The abstract states a 'systematic study' of stability, accuracy, complexity, and AP property; the corresponding sections should include a summary table or clear cross-references so readers can quickly compare the three partitionings.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and valuable suggestions on clarifying the asymptotic preserving (AP) analysis and strengthening the numerical evidence. We address the two major comments point by point below. Both points can be addressed by targeted revisions that make the existing analysis and results more explicit without altering the core contributions.

read point-by-point responses

  1. Referee: [analysis of AP property (likely §4)] The load-bearing claim is that each of the three IMEX-BDF partitionings, when combined with the micro-macro decomposition and DG discretization, yields an AP scheme whose ε→0 limit (fixed h, Δt) recovers a consistent discretization of the diffusion equation. The manuscript must explicitly pass to the limit inside the fully discrete equations for each partitioning separately, showing that explicit contributions vanish and the micro component relaxes to zero without residual inconsistent fluxes. This verification is required in the analysis section to support the uniform AP assertion.

    Authors: We agree that an explicit, side-by-side passage to the ε→0 limit inside the fully discrete equations for each of the three partitionings would improve clarity. The current analysis in §4 derives the AP property from the micro-macro decomposition and shows that the limiting scheme is consistent with a DG discretization of the diffusion equation; however, the limiting process is presented in a unified manner rather than expanded separately for each IMEX-BDF partitioning. In the revised manuscript we will insert a new subsection (or expanded paragraphs) in §4 that performs this limit explicitly for each partitioning, confirming that the explicit flux terms vanish, the micro component relaxes to zero, and no inconsistent residual fluxes remain. revision: yes

  2. Referee: [numerical results section (likely §5)] Numerical tables demonstrating the AP property should report errors or solution profiles for successively smaller ε with fixed mesh size and time step, for all three partitionings. Without such data, it is difficult to confirm that the discrete limit is indeed the correct diffusion discretization rather than a perturbed or inconsistent scheme.

    Authors: We appreciate this concrete recommendation. Section 5 already contains numerical experiments that include small-ε regimes to illustrate the AP behavior and compares the three BDF-based schemes with their Runge-Kutta counterparts. To make the uniform convergence to the diffusion limit fully transparent, we will add a dedicated table (or set of tables) in §5 that reports L² errors (or solution profiles) against a reference diffusion solution for a sequence of decreasing ε values while keeping h and Δt fixed, for each of the three partitionings. This will directly demonstrate that the schemes recover the correct diffusion discretization uniformly. revision: yes

Circularity Check

0 steps flagged

No significant circularity; schemes derived from standard micro-macro and IMEX components with independent analysis

full rationale

The derivation begins from the established micro-macro decomposition of the linear kinetic model, applies standard DG spatial discretization and S_N velocity discretization, then introduces three explicit IMEX partitionings into BDF time stepping. The AP property is asserted to follow from passing to the ε→0 limit inside the fully discrete equations for each partitioning separately, with the explicit terms vanishing and the macro equation closing to a consistent diffusion discretization. This limit passage is a direct algebraic verification on the scheme equations rather than a self-definition or fitted-parameter renaming. Self-citations to prior DG-AP work appear only as context for the continuing effort and do not supply the uniqueness or correctness of the three new partitionings. No step reduces the claimed AP property, stability, or accuracy to an input by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard numerical analysis assumptions for consistency and stability of DG and IMEX schemes plus the validity of the micro-macro decomposition for the linear kinetic model.

axioms (1)
  • standard math Standard assumptions on stability and consistency of discontinuous Galerkin methods and IMEX time integrators for hyperbolic and diffusive regimes.
    Invoked throughout the analysis of the proposed schemes.

pith-pipeline@v0.9.0 · 5627 in / 1199 out tokens · 59819 ms · 2026-05-10T20:01:02.950430+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

42 extracted references · 42 canonical work pages

  1. [1]

    Discontinuous finite element transport solutions in thick diffusive prob- lems

    Marvin L Adams. Discontinuous finite element transport solutions in thick diffusive prob- lems. Nuclear science and engineering, 137(3):298–333, 2001

    work page 2001

  2. [2]

    Fast iterative methods for discrete-ordinates particle transport calculations

    Marvin L Adams and Edward W Larsen. Fast iterative methods for discrete-ordinates particle transport calculations. Progress in nuclear energy, 40(1):3–159, 2002

    work page 2002

  3. [3]

    Implicit-explicit methods for time- dependent partial differential equations

    Uri M Ascher, Steven J Ruuth, and Brian TR Wetton. Implicit-explicit methods for time- dependent partial differential equations. SIAM Journal on Numerical Analysis, 32(3):797– 823, 1995. 32

    work page 1995

  4. [4]

    Implicit-explicit Runge- Kutta schemes for hyperbolic systems and kinetic equations in the diffusion limit

    Sebastiano Boscarino, Lorenzo Pareschi, and Giovanni 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

  5. [5]

    Coupling Boltzmann and Euler equations without overlapping

    JF Bourgat, Patrick Le Tallec, B Perthame, and Y Qiu. Coupling Boltzmann and Euler equations without overlapping. Contemporary Mathematics, 157:377–398, 1994

    work page 1994

  6. [6]

    Asymp- totically complexity diminishing schemes (ACDS) for kinetic equations in the diffusive scaling

    Ana¨ ıs Crestetto, Nicolas Crouseilles, Giacomo Dimarco, and Mohammed Lemou. Asymp- totically complexity diminishing schemes (ACDS) for kinetic equations in the diffusive scaling. Journal of Computational Physics, 394:243–262, 2019

    work page 2019

  7. [7]

    A smooth transition model between kinetic and diffusion equations

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

    work page 2005

  8. [8]

    Implicit-explicit linear multistep methods for stiff kinetic equations

    Giacomo Dimarco and Lorenzo Pareschi. Implicit-explicit linear multistep methods for stiff kinetic equations. SIAM Journal on Numerical Analysis, 55(2):664–690, 2017

    work page 2017

  9. [9]

    Asymptotic-preserving Monte Carlo methods for transport equations in the diffusive limit

    Giacomo Dimarco, Lorenzo Pareschi, and Giovanni Samaey. Asymptotic-preserving Monte Carlo methods for transport equations in the diffusive limit. SIAM Journal on Scientific Computing, 40(1):A504–A528, 2018

    work page 2018

  10. [10]

    An asymptotic-preserving dynamical low- rank method for the multi-scale multi-dimensional linear transport equation

    Lukas Einkemmer, Jingwei Hu, and Yubo Wang. An asymptotic-preserving dynamical low- rank method for the multi-scale multi-dimensional linear transport equation. Journal of Computational Physics, 439:110353, 2021

    work page 2021

  11. [11]

    A hierarchy of hybrid numerical methods for multiscale kinetic equations

    Francis Filbet and Thomas Rey. A hierarchy of hybrid numerical methods for multiscale kinetic equations. SIAM Journal on Scientific Computing, 37(3):A1218–A1247, 2015

    work page 2015

  12. [12]

    The convergence of numerical transfer schemes in diffusive regimes i: Discrete-ordinate method

    Fran¸ cois Golse, Shi Jin, and C David Levermore. The convergence of numerical transfer schemes in diffusive regimes i: Discrete-ordinate method. SIAM Journal on Numerical Analysis, 36(5):1333–1369, 1999

    work page 1999

  13. [13]

    Asymptotic analysis of upwind discontinuous Galerkin approximation of the radiative transport equation in the diffusive limit

    Jean-Luc Guermond and Guido Kanschat. Asymptotic analysis of upwind discontinuous Galerkin approximation of the radiative transport equation in the diffusive limit. SIAM Journal on Numerical Analysis, 48(1):53–78, 2010

    work page 2010

  14. [14]

    A low rank tensor representation of linear transport and non- linear Vlasov solutions and their associated flow maps

    Wei Guo and Jing-Mei Qiu. A low rank tensor representation of linear transport and non- linear Vlasov solutions and their associated flow maps. Journal of Computational Physics, 458:111089, 2022

    work page 2022

  15. [15]

    Solving Ordinary Differential Equations I: Nonstiff Problems

    Ernst Hairer, Syvert Paul Nørsett, and Gerhard Wanner. Solving Ordinary Differential Equations I: Nonstiff Problems. Springer-Verlag, 1987

    work page 1987

  16. [16]

    Nodal discontinuous Galerkin methods: algorithms, analysis, and applications

    Jan S Hesthaven and Tim Warburton. Nodal discontinuous Galerkin methods: algorithms, analysis, and applications. Springer Science & Business Media, 2007

    work page 2007

  17. [17]

    Analysis of asymptotic preserving DG-IMEX schemes for linear kinetic transport equations in a diffusive scaling

    Juhi Jang, Fengyan Li, Jing-Mei Qiu, and Tao 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

  18. [18]

    High order asymptotic preserving DG-IMEX schemes for discrete-velocity kinetic equations in a diffusive scaling

    Juhi Jang, Fengyan Li, Jing-Mei Qiu, and Tao 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

  19. [19]

    Asymptotic preserving (AP) schemes for multiscale kinetic and hyperbolic equa- tions: a review

    Shi Jin. Asymptotic preserving (AP) schemes for multiscale kinetic and hyperbolic equa- tions: a review. Lecture notes for summer school on methods and models of kinetic theory (M&MKT), Porto Ercole (Grosseto, Italy), pages 177–216, 2010

    work page 2010

  20. [20]

    Uniformly accurate diffusive relax- ation schemes for multiscale transport equations

    Shi Jin, Lorenzo Pareschi, and Giuseppe Toscani. Uniformly accurate diffusive relax- ation schemes for multiscale transport equations. SIAM Journal on Numerical Analysis, 38(3):913–936, 2000

    work page 2000

  21. [21]

    On numerical solutions of transport problems in the diffusion limit

    Edward W Larsen. On numerical solutions of transport problems in the diffusion limit. Nuclear Science and Engineering, 83(1):90–99, 1983

    work page 1983

  22. [22]

    Asymptotic solutions of numerical transport problems in optically thick, diffusive regimes

    Edward W Larsen, Jim E Morel, and Warren F Miller Jr. Asymptotic solutions of numerical transport problems in optically thick, diffusive regimes. Journal of Computational Physics, 69(2):283–324, 1987

    work page 1987

  23. [23]

    Micro-macro schemes for kinetic equations includ- ing boundary layers

    Mohammed Lemou and Florian M´ ehats. Micro-macro schemes for kinetic equations includ- ing boundary layers. SIAM Journal on Scientific Computing, 34(6):B734–B760, 2012. 33

    work page 2012

  24. [24]

    A new asymptotic preserving scheme based on micro-macro formulation for linear kinetic equations in the diffusion limit

    Mohammed Lemou and Luc 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

  25. [25]

    Computational Methods of Neutron Transport

    Elmer Eugene Lewis and Warren F Miller. Computational Methods of Neutron Transport. John Wiley and Sons, Inc., New York, NY, 1984

    work page 1984

  26. [26]

    Boltzmann equation: micro-macro decompositions and positivity of shock profiles

    Tai-Ping Liu and Shih-Hsien Yu. Boltzmann equation: micro-macro decompositions and positivity of shock profiles. Communications in mathematical physics, 246(1):133–179, 2004

    work page 2004

  27. [27]

    Multi-scale simulation and model order reduction for the radiative transfer equation

    Kimberly Matsuda. Multi-scale simulation and model order reduction for the radiative transfer equation. PhD thesis, Rensselaer Polytechnic Institute, 2025

    work page 2025

  28. [28]

    Efficient asymptotic preserving deterministic meth- ods for the Boltzmann equation

    Lorenzo Pareschi and Giovanni Russo. Efficient asymptotic preserving deterministic meth- ods for the Boltzmann equation. Models and Computational Methods for Rarefied Flows, AVT-194 RTO AVT/VKI, Rhode St. Genese, Belgium, page 34, 2011

    work page 2011

  29. [29]

    OEDG: Oscillation-eliminating discontinuous Galerkin method for hyperbolic conservation laws

    Manting Peng, Zheng Sun, and Kailiang Wu. OEDG: Oscillation-eliminating discontinuous Galerkin method for hyperbolic conservation laws. Mathematics of Computation, 2024

    work page 2024

  30. [30]

    Structure-preserving discontinuous Galerkin methods for multi-scale kinetic transport equations and nonlinear optics models

    Zhichao Peng. Structure-preserving discontinuous Galerkin methods for multi-scale kinetic transport equations and nonlinear optics models. PhD thesis, Rensselaer Polytechnic In- stitute, 2020

    work page 2020

  31. [31]

    Asymptotic and positiv- ity preserving methods for Kerr-Debye model with Lorentz dispersion in one dimension

    Zhichao Peng, Vrushali A Bokil, Yingda Cheng, and Fengyan Li. Asymptotic and positiv- ity preserving methods for Kerr-Debye model with Lorentz dispersion in one dimension. Journal of Computational Physics, 402:109101, 2020

    work page 2020

  32. [32]

    A micro-macro decomposed re- duced basis method for the time-dependent radiative transfer equation.Multiscale Modeling & Simulation, 22(1):639–666, 2024

    Zhichao Peng, Yanlai Chen, Yingda Cheng, and Fengyan Li. A micro-macro decomposed re- duced basis method for the time-dependent radiative transfer equation.Multiscale Modeling & Simulation, 22(1):639–666, 2024

    work page 2024

  33. [33]

    Stability-enhanced AP IMEX- LDG schemes for linear kinetic transport equations under a diffusive scaling

    Zhichao Peng, Yingda Cheng, Jing-Mei Qiu, and Fengyan Li. Stability-enhanced AP IMEX- LDG schemes for linear kinetic transport equations under a diffusive scaling. Journal of Computational Physics, 415:109485, 2020

    work page 2020

  34. [34]

    Stability-enhanced AP IMEX1-LDG method: energy-based stability and rigorous AP property

    Zhichao Peng, Yingda Cheng, Jing-Mei Qiu, and Fengyan 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

  35. [35]

    Asymptotic preserving IMEX-DG-S schemes for linear kinetic transport equations based on Schur complement

    Zhichao Peng and Fengyan 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

  36. [36]

    Operator bounds and time step conditions for the DG and central DG methods

    Matthew A Reyna and Fengyan Li. Operator bounds and time step conditions for the DG and central DG methods. Journal of Scientific Computing, 62(2):532–554, 2015

    work page 2015

  37. [37]

    Reduced order modeling applied to the discrete ordinates method for radiation heat transfer in participating media

    John Tencer, Kevin Carlberg, Roy Hogan, and Marvin Larsen. Reduced order modeling applied to the discrete ordinates method for radiation heat transfer in participating media. In Heat Transfer Summer Conference, volume 50336, page V002T15A011. American Society of Mechanical Engineers, 2016

    work page 2016

  38. [38]

    Variable step-size implicit-explicit linear multistep methods for time-dependent partial differential equations

    Dong Wang and Steven J Ruuth. Variable step-size implicit-explicit linear multistep methods for time-dependent partial differential equations. Journal of Computational Mathematics, pages 838–855, 2008

    work page 2008

  39. [39]

    High order asymptotic preserv- ing nodal discontinuous Galerkin IMEX schemes for the BGK equation

    Tao Xiong, Juhi Jang, Fengyan Li, and Jing-Mei Qiu. High order asymptotic preserv- ing nodal discontinuous Galerkin IMEX schemes for the BGK equation. Journal of Computational Physics, 284:70–94, 2015

    work page 2015

  40. [40]

    The L2-norm stability analysis of Runge–Kutta discontinuous Galerkin methods for linear hyperbolic equations

    Yuan Xu, Qiang Zhang, Chi-wang Shu, and Haijin Wang. The L2-norm stability analysis of Runge–Kutta discontinuous Galerkin methods for linear hyperbolic equations. SIAM Journal on Numerical Analysis, 57(4):1574–1601, 2019

    work page 2019

  41. [41]

    Adaptive time-stepping schemes for the solution of the Poisson-Nernst-Planck equations

    David Yan, Mary C Pugh, and Francis P Dawson. Adaptive time-stepping schemes for the solution of the Poisson-Nernst-Planck equations. Applied Numerical Mathematics, 163:254–269, 2021

    work page 2021

  42. [42]

    Asymptotic preserving and uniformly un- conditionally stable finite difference schemes for kinetic transport equations

    Guoliang Zhang, Hongqiang Zhu, and Tao Xiong. Asymptotic preserving and uniformly un- conditionally stable finite difference schemes for kinetic transport equations. SIAM Journal on Scientific Computing, 45(5):B697–B730, 2023. 34 Appendix A An example of non-AP implicit methods Withα(x, t) =f(x, v= 1, t),β(x, t) =f(x, v=−1, t), the telegraph equation is g...

    work page 2023