Pith. sign in

REVIEW 3 major objections 4 minor 1 cited by

Asymptotic-Preserving Dynamical Low-Rank Method for the Stiff Nonlinear Boltzmann Equation

T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read A dynamical low-rank integrator solves the stiff Boltzmann equation with only $r^2$ collision evaluations per time step and time steps independent of the Knudsen number.

desk verdict Two genuinely new DLR integrators with a clean AP proof for a BGK model; the abstract's Boltzmann AP claim is an extrapolation from that model and Maxwellian-only tests, but the paper is honest and deserves refereeing. read the letter →

arxiv 2502.08951 v1 pith:JEMJKRLL submitted 2025-02-13 math.NA cs.NA

Lukas Einkemmer , Jingwei Hu , Shiheng Zhang This is my paper · ORCID
classification math.NAcs.NA MSC 35Q2065F55
keywords Boltzmannequationdynamicallow-rankapproximationasymptotic-preservingschemestiffcollisionoperatorIMEXpenalizationBGKintegratorKnudsennumber
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Numerically solving the Boltzmann equation is prohibitively expensive because the collision operator is a high-dimensional integral evaluated at every spatial grid point and every time step, and it becomes stiff in the fluid regime. This paper claims that advancing a dynamical low-rank representation of the distribution function with two new integrators, called XL and sXL, removes that bottleneck: the Boltzmann collision operator is evaluated only $r^2$ times per step, where the rank $r$ is far below the number of spatial grid points. The same construction is what makes the schemes asymptotic-preserving, so time steps can be chosen independently of the Knudsen number $\varepsilon$ and the method reduces to a consistent discretization of the compressible Euler equations as $\varepsilon \to 0$. A reader should care because this offers a concrete route to computing Boltzmann dynamics at fluid-regime cost while retaining kinetic detail when it matters.

What carries the argument

The central object is the XL integrator, a factored-step scheme for a matrix differential equation that updates the spatial basis through one forward differential equation, updates the velocity coefficients through a second, and then truncates the rank back to $r$ by singular value decomposition; it is exact on exact rank-$r$ solutions. Its specialization, the sXL integrator, exploits the separation property $F(KV^T) = \mathcal{K}(K)\mathcal{V}(V)^T$ satisfied by the Boltzmann collision operator to replace the first differential equation with a single algebraic augmentation of the basis by products of existing basis functions, so only one differential equation is solved per step. The stiff part of the split system is written as $(Q(f,f)+\lambda f-\lambda f)/\varepsilon$ with $\lambda$ chosen near the collision loss rate, making the implicit term a cheap scalar multiple of $f$ inside the low-rank factors.

What would settle it

A direct test is to run the proposed XL scheme on the full Boltzmann equation starting from non-equilibrium initial data at small Knudsen number; if $\|f^n - M^n\|$ does not drop to $O(\Delta t)$ within a few steps while a full-tensor AP scheme does, then the central AP claim beyond the BGK/equilibrium setting fails.

Watch

Extended reading notes

Core claim

The central claim is that the stiff nonlinear Boltzmann equation can be advanced on a low-rank manifold $f = \sum_{i,j} X_i(x,t) S_{ij}(t) V_j(v,t)$ using two new integrators, the XL integrator and its separable-problem specialization sXL, that are cheaper than existing projector-splitting or basis-update-and-Galerkin alternatives. The key algorithmic ingredients are an operator splitting that treats advection explicitly and treats the stiff collision part with a BGK-type penalty $\lambda f$ handled implicitly, avoiding both implicit evaluation of the Boltzmann collision operator and costly projection of the Maxwellian. As a result the collision operator is evaluated $r^2$ times per step instead of once per spatial grid point. The paper proves the asymptotic-preserving property for a special one-dimensional BGK equation whose Maxwellian is a rank-3 function, assuming equilibrium initial data and rank at least six, and demonstrates with sine-wave and shock-tube experiments in 1D+2V that the schemes match a full-tensor asymptotic-preserving solver for $\varepsilon$ from 1 down to $10^{-6}$ with time steps independent of $\varepsilon$.

Load-bearing premise

The asymptotic-preserving property is proved only for a simplified rank-3 BGK collision model with equilibrium initial data, and the paper assumes the same behavior transfers to the full nonlinear Boltzmann collision operator and to general initial data; the paper states that proving AP for the full Boltzmann equation is challenging.

Editorial extensions

If this is right

  • The collision operator is evaluated $O(r^2)$ times per time step instead of once per spatial grid point, so the dominant cost no longer multiplies by the number of grid points.
  • Time steps can be taken independently of $\varepsilon$ in the stiff fluid regime, subject only to the advective CFL condition, as the sine and shock-tube tests at $\varepsilon = 10^{-6}$ show.
  • In the limit $\varepsilon \to 0$ with equilibrium initial data, the scheme becomes a first-order consistent discretization of the compressible Euler equations, so one solver covers kinetic and fluid regimes without switching.
  • The sXL variant reduces a low-rank step to one differential equation for separable kinetic operators, a class that includes advection, Fokker-Planck terms, and special BGK-type terms with low-rank Maxwellians.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper: if the contraction argument around $\lambda$ were extended to the full Boltzmann operator with non-equilibrium data, the asymptotic-preserving claim would apply generally; a direct test is to start from a strongly non-Maxwellian $f_0$ at small $\varepsilon$ and monitor $\|f^n - M^n\|$.
  • Beyond the paper: the sXL mechanism appears transferable to other separated kinetic operators, such as radiation transport, Fokker-Planck models, and stochastic-Galerkin uncertainty quantification systems, where one differential equation per step would replace the full update.
  • Beyond the paper: the rank needed for accuracy in the fluid regime is likely tied to the polynomial complexity of the local Maxwellian, which suggests principled choices of $r$ or adaptive tolerances instead of manual ranks.
  • Beyond the paper: since the scheme is first order and the authors note the underlying full-tensor AP scheme limits the accuracy, a second-order IMEX time integrator with the same penalization is a natural extension the paper does not test.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper proposes two dynamical low-rank integrators, XL and sXL, for matrix differential equations, and applies them to the stiff nonlinear Boltzmann equation with a penalty-based IMEX splitting. The XL integrator augments the spatial basis by a K-step and then performs one L-step, while the sXL integrator exploits a separability property of the right-hand side to replace the K-step by direct basis augmentation. In the kinetic application, advection is handled by an explicit KSL-type splitting and collision by the penalized IMEX update (4.19)–(4.20), which requires only r^2 evaluations of the Boltzmann collision operator per time step. The authors prove an exactness property for the generic XL integrator, prove an asymptotic-preserving (AP) property for a special 1D1V BGK model with equilibrium initial data, and present 1+2 dimensional numerical tests for sine-type and shock-tube initial data showing accuracy comparable to the full-tensor scheme for epsilon from 1 down to 1e-6.

Significance. If the advertised claims hold, the paper would make a useful contribution to kinetic simulation: the O(r^2) collision-evaluation count is a concrete and well argued complexity improvement, and epsilon-independent time stepping in the stiff regime is practically important. The generic XL integrator and the separable sXL construction are elegant, and the exactness proof in Theorem 2.2 is transparent. The AP analysis for the special BGK model is nontrivial and goes beyond a formal consistency check. However, the central AP claim for the nonlinear Boltzmann equation is not actually proved: Theorem 5.1 covers only a rank-3 BGK operator with Maxwellian initial data, and the contraction condition (4.9) needed for the Boltzmann operator is left open. The numerical experiments all start from equilibrium initial data, so the non-equilibrium contraction mechanism is neither proved nor tested. This gap is substantial and should be addressed before the abstract and conclusion claim an AP scheme for the stiff nonlinear Boltzmann equation.

major comments (3)
  1. [§4 (Eqs. (4.19)–(4.20)) and §5 (Theorem 5.1)] The AP property is proven only for the special rank-3 BGK operator (5.1)–(5.2) and only under the assumption f0 = M0. The paper itself states at the start of Section 5 that proving AP for the nonlinear Boltzmann operator is challenging. For the Boltzmann collision step, the epsilon-to-zero limit of (4.19)–(4.20) is not an exact projection onto the Maxwellian; it gives an explicit correction f^{n+1} = f^* + (Delta t / lambda) Q(f^*, f^*) in the augmented basis, and AP depends on the contraction inequality (4.9), whose validity for the Boltzmann operator is explicitly acknowledged just before (4.13) to be not easy to show. The abstract's unconditional statement that the proposed low-rank schemes are asymptotic-preserving is therefore an extrapolation from the BGK model. The authors should either extend the proof to the Boltzmann setting under stated conditions, or weaken the abstract and conclusion to say that AP is established for the special BGK model and only formally expected, with numerical evidence, for the full Boltzmann operator.
  2. [§5 (Theorem 5.1) and §6 (Figs. 1–6)] Even for the special BGK model, Theorem 5.1 assumes equilibrium initial data f0 = M0. The proof uses this assumption to obtain f^n = M^n inductively and therefore does not exercise the contraction transient (4.11) that is the stated mechanism for AP from non-equilibrium data. All numerical experiments in Section 6 also start from Maxwellian initial data, namely (6.1) for the sine problem and (6.3) for the shock tube. Consequently, the non-equilibrium regime is neither proved nor numerically tested. At minimum, please add a non-Maxwellian initial-data experiment and state clearly that the AP theorem is conditional on equilibrium initial data.
  3. [§4 (Approaches 1 and 2) and §5 (Theorem 5.1)] The AP proof for the sXL integrator uses the full augmentation (5.9), which in the Boltzmann setting corresponds to augmenting by all r^2 products {X_k^* X_l^*}. The two heuristic rank-selection variants DLR-sXL-1 and DLR-sXL-2, which select subsets of these products by coefficient magnitude or by a tolerance, are used in the numerical experiments but are not covered by the theorem. No error or AP analysis is given for these selected subsets, so the numerical results for DLR-sXL-1 and DLR-sXL-2 are not tied to the proven AP statement. Please either analyze these variants under the stated assumptions or present them clearly as heuristic and not part of the AP guarantee.
minor comments (4)
  1. [§2 and §4] There are typographical artifacts 'T runcation' and 'T runcation-step' in Section 2 and in Step 6 of Section 4; these should be corrected to 'Truncation'.
  2. [§6] In the setup of Section 6, 'V elocity domain' should read 'Velocity domain'.
  3. [References] References [25] and [26] are the same paper (Hu and Wang, J. Sci. Comput. 92 (2022), p. 75) and should be consolidated to avoid duplication.
  4. [§4] The step numbering is confusing: 'Step 4' is used for both the XL X-step in (4.19) and the sXL augmentation step in (4.22). Please label the variants distinctly, for example Step 4a and Step 4b, so that the algorithm is unambiguous.

Circularity Check

0 steps flagged · score 1.0 of 10

No material circularity: the AP proof is self-contained, though the full-Boltzmann AP claim is extrapolated from a BGK analog with Maxwellian initial data.

full rationale

The paper's central derivation is not circular. The AP property is proven in Section 5 for a special 1D1V BGK equation by direct manipulation of the discrete low-rank equations: equation (5.15) computes the macroscopic flux change exactly, equations (5.16)-(5.26) bound the intermediate updates using the discretization itself (e.g., S^(1)_ij - S^(2)_ij = O(Delta t) from (4.17)), and equations (5.27)-(5.33) show that in the epsilon -> 0 limit the low-rank factors collapse to the Maxwellian. No fitted constant, no assumed Euler limit, and no imported uniqueness theorem is used. The O(r^2) collision-evaluation count is a structural property of the scheme, not a prediction. Self-citations appear (e.g., [24], [25]) but they supply building blocks (full-tensor AP splitting, nonstiff DLR method) that the paper re-derives or verifies in its own proof; they are not load-bearing in a circular sense. The main weakness is scope, not circularity: Theorem 5.1 proves AP only for the simplified BGK operator with rank-3 Maxwellian and equilibrium initial data, while the abstract claims AP for the stiff nonlinear Boltzmann equation. That is an extrapolation or a limitation, not a reduction of the result to its inputs, so it does not raise the circularity score above the minor self-citation range.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The central claim rests on the low-rank ansatz, the penalty contraction condition, and the transfer of the AP proof from a simplified BGK model to the full Boltzmann operator. No experimental constants are fitted; the free parameters are algorithmic choices such as rank, tolerance, and the penalty parameter.

free parameters (3)
  • Penalty parameter lambda = 1.1
    Introduced in the BGK penalized splitting in Section 4. For the special BGK equation the AP proof needs lambda greater than 1/2, while for the full Boltzmann operator the existence of lambda satisfying contraction condition (4.9) is not established. The value 1.1 is a heuristic choice taken from prior AP work.
  • Low-rank r and augmented tilde_r = 6, 10, 14, 20 depending on test; tilde_r = r
    The rank is chosen by hand for each numerical test, and the AP theorem assumes r at least 6. The central claim that r is much smaller than the number of spatial grid points is an assumption of the method rather than a proven property.
  • Augmentation tolerance tol = 1, 0.05, 0.01, etc. depending on test
    Tolerance controls how many product basis functions are added in the DLR-sXL-2 variant. It is selected empirically in Section 6 and plays no role in the AP proof.
assumptions (4)
  • domain assumption The separability condition (2.7) applies to the Boltzmann collision operator, so F(KV^T) = K(K)V(V)^T and the augmented x-basis can be built from products X_k X_l.
    Used in Sections 2 and 4 to derive the sXL integrator. It is true for the bilinear, velocity-local collision operator, but the rank of the augmented basis is assumed not to explode in the fluid regime.
  • ad hoc to paper There exists a penalty parameter lambda such that sup|D*|/lambda < 1 for the full Boltzmann collision operator, as required by (4.9).
    The contraction condition underlies the AP argument in equations (4.8) through (4.11). The paper states that showing this for the Boltzmann operator is not easy, and only numerical evidence with lambda = 1.1 is provided.
  • domain assumption The special BGK Maxwellian (5.2) is a rank-3 function and the Hermite recurrence (5.13) holds for the orthonormalized velocity basis.
    Used in Theorem 5.1 to track moments through the splitting steps. It is true for the weakly compressible model but does not transfer to the full Maxwellian, which is not low-rank in general.
  • domain assumption SVD truncation and orthonormalization preserve accuracy and no rank explosion occurs as epsilon tends to 0.
    Assumed throughout Sections 5 and 6. The AP proof shows the rank stays 3 for the special BGK, but for the full Boltzmann operator rank growth is only controlled numerically by the chosen tolerances.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Asymptotic-Preserving Dynamical Low-Rank Method for the Stiff Nonlinear Boltzmann Equation." pith.science (2026). https://pith.science/paper/JEMJKRLL

@misc{pith2026250208951,
  author       = {Pith},
  title        = {Pith review of: Asymptotic-Preserving Dynamical Low-Rank Method for the Stiff Nonlinear Boltzmann Equation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JEMJKRLL}},
  note         = {Machine review of arXiv:2502.08951}
}
abstract

In kinetic theory, numerically solving the full Boltzmann equation is extremely expensive. This is because the Boltzmann collision operator involves a high-dimensional, nonlinear integral that must be evaluated at each spatial grid point and every time step. The challenge becomes even more pronounced in the fluid (strong collisionality) regime, where the collision operator exhibits strong stiffness, causing explicit time integrators to impose severe stability restrictions. In this paper, we propose addressing this problem through a dynamical low-rank (DLR) approximation. The resulting algorithm requires evaluating the Boltzmann collision operator only $r^2$ times, where $r$, the rank of the approximation, is much smaller than the number of spatial grid points. We propose a novel DLR integrator, called the XL integrator, which reduces the number of steps compared to the available alternatives (such as the projector splitting or basis update & Galerkin (BUG) integrator). For a class of problems including the Boltzmann collision operator which enjoys a separation property between physical and velocity space, we further propose a specialized version of the XL integrator, called the sXL integrator. This version requires solving only one differential equation to update the low-rank factors. Furthermore, the proposed low-rank schemes are asymptotic-preserving, meaning they can capture the asymptotic fluid limit in the case of strong collisionality. Our numerical experiments demonstrate the efficiency and accuracy of the proposed methods across a wide range of regimes, from non-stiff (kinetic) to stiff (fluid).

Figures

Figures reproduced from arXiv: 2502.08951 by the authors.

Figure 1
Figure 1. Comparison of density, bulk velocity (first component), and temperature for the sine problem with [PITH_FULL_IMAGE:figures/full_fig_p020_1.png] view at source ↗
Figure 2
Figure 2. Comparison of density, bulk velocity (first component), and temperature for the sine problem with [PITH_FULL_IMAGE:figures/full_fig_p021_2.png] view at source ↗
Figure 3
Figure 3. Comparison of DLR-sXL-2 with rank r = 10 under three different tolerances, tol ∈ {0.2, 0.1, 0.05}, for ε = 10−6 . Solid line: full tensor method, red circle: DLR-sXL-2 with tol = 0.2, green star: DLR￾sXL-2 with tol = 0.1, blue diamond: DLR-sXL-2 with tol = 0.05. In Figure 3c, one can observe a slight discrepancy between DLR-sXL-2-0.2 and others. The additional ranks refers to the number of extra basis functions intr… view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Comparison of density, bulk velocity (first component), and temperature for the shock tube problem [PITH_FULL_IMAGE:figures/full_fig_p023_4.png]
Figure 5
Figure 5. Figure 5: Comparison of density, bulk velocity (first component), and temperature for the shock tube problem [PITH_FULL_IMAGE:figures/full_fig_p024_5.png]
Figure 6
Figure 6. Figure 6: Comparison of DLR-sXL-2 with rank r = 20 under three different tolerances, tol ∈ {0.1, 0.05, 0.01}, for ε = 10−6 . Solid line: full tensor method, red circle: DLR-sXL-2 with tol = 0.1, green star: DLR-sXL-2 with tol = 0.05, blue diamond: DLR-sXL-2 with tol = 0.01. The …

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. A Galerkin Alternating Projection Method for Kinetic Equations in the Diffusive Limit

    math.NA 2025-05 conditional novelty 6.0 of 10

    The GAP scheme is a new dynamical low-rank integrator for the radiative transfer equation that provably preserves the diffusive limit and avoids CFL restrictions.

Reference graph

Works this paper leans on

36 extracted references · 25 canonical work pages · cited by 1 Pith paper

  1. [1]

    Bardos, F

    C. Bardos, F. Golse, and D. Levermore, Fluid dynamic limits of kinetic equations. I. Formal derivations, J. Stat. Phys., 63 (1991), pp. 323–344

    work page 1991

  2. [2]

    Energy stable and conservative dynamical low-rank approximation for the Su-Olson problem

    L. Baumann, L. Einkemmer, C. Klingenberg, and J. Kusch, Energy stable and conservative dynamical low-rank approximation for the Su-Olson problem , (2023), https://doi.org/10.48550/ ARXIV.2307.07538

    work page Pith review arXiv doi:10.48550/arxiv.2307.07538 2023

  3. [3]

    Cercignani, The Boltzmann Equation and Its Applications , Springer-Verlag, New York, 1988

    C. Cercignani, The Boltzmann Equation and Its Applications , Springer-Verlag, New York, 1988

    work page 1988

  4. [4]

    Ceruti, J

    G. Ceruti, J. Kusch, and C. Lubich, A rank-adaptive robust integrator for dynamical low-rank approximation, BIT Numerical Mathematics, 62 (2022), pp. 1149–1174

    work page 2022

  5. [5]

    Ceruti, J

    G. Ceruti, J. Kusch, and C. Lubich, A Parallel Rank-Adaptive Integrator for Dynamical Low- Rank Approximation, SIAM J. Sci. Comput., 46 (2024), pp. B205–B228, https://doi.org/10.1137/ 23m1565103

    work page 2024

  6. [6]

    Ceruti and C

    G. Ceruti and C. Lubich, An unconventional robust integrator for dynamical low-rank approximation, BIT Numerical Mathematics, 62 (2022), pp. 23–44

    work page 2022

  7. [7]

    Chen and G

    S. Chen and G. D. Doolen, Lattice Boltzmann method for fluid flows , Annual review of fluid mechanics, 30 (1998), pp. 329–364

    work page 1998

  8. [8]

    Coughlin and J

    J. Coughlin and J. Hu, Efficient dynamical low-rank approximation for the Vlasov-Amp` ere-Fokker- Planck system , J. Comput. Phys., 470 (2022), p. 111590, https://doi.org/10.1016/j.jcp.2022. 111590

    doi:10.1016/j.jcp.2022 2022

Show all 36 references
  1. [9]

    Coughlin, J

    J. Coughlin, J. Hu, and U. Shumlak, Robust and conservative dynamical low-rank methods for the Vlasov equation via a novel macro-micro decomposition , J. Comput. Phys., 509 (2024), p. 113055

    2024

  2. [10]

    Dektor and L

    A. Dektor and L. Einkemmer, Interpolatory dynamical low-rank approximation for the 3+3d boltzmann-bgk equation, arXiv:2411.15990, (2024), https://doi.org/10.48550/ARXIV.2411.15990

    2024 doi arXiv

  3. [11]

    Z. Ding, L. Einkemmer, and Q. Li, Dynamical Low-Rank Integrator for the Linear Boltzmann Equation: Error Analysis in the Diffusion Limit , SIAM J. Numer. Anal., 59 (2021), https://doi.org/ 10.1137/20M1380788

    2021 doi

  4. [12]

    Einkemmer, A low-rank algorithm for weakly compressible flow , SIAM Journal on Scientific Computing, 41 (2019), pp

    L. Einkemmer, A low-rank algorithm for weakly compressible flow , SIAM Journal on Scientific Computing, 41 (2019), pp. A2795–A2814. 26

    2019

  5. [13]

    Einkemmer, Accelerating the simulation of kinetic shear Alfv´ en waves with a dynamical low-rank approximation, J

    L. Einkemmer, Accelerating the simulation of kinetic shear Alfv´ en waves with a dynamical low-rank approximation, J. Comput. Phys., 501 (2024), p. 112757, https://doi.org/10.1016/j.jcp.2024. 112757

    2024 doi

  6. [14]

    Einkemmer, J

    L. Einkemmer, J. Hu, and J. Kusch, Asymptotic-Preserving and Energy Stable Dynamical Low- Rank Approximation, SIAM Journal on Numerical Analysis, 62 (2024), pp. 73–92, https://doi.org/ 10.1137/23m1547603

    2024 doi

  7. [15]

    Einkemmer, J

    L. Einkemmer, J. Hu, and Y. W ang, An asymptotic-preserving dynamical low-rank method for the multi-scale multi-dimensional linear transport equation , J. Comput. Phys., 439 (2021)

    2021

  8. [16]

    Einkemmer, J

    L. Einkemmer, J. Hu, and L. Ying, An Efficient Dynamical Low-Rank Algorithm for the Boltzmann- BGK Equation Close to the Compressible Viscous Flow Regime , SIAM J. Sci. Comput., 43 (2021), pp. B1057–B1080, https://doi.org/10.1137/21m1392772

    2021 doi

  9. [17]

    Einkemmer, K

    L. Einkemmer, K. Kormann, J. Kusch, R. McClarren, and J.-M. Qiu, A review of low-rank methods for time-dependent kinetic simulations , arXiv:2412.05912, (2024), https://arxiv.org/abs/ 2412.05912

    2024 arXiv

  10. [18]

    Einkemmer and C

    L. Einkemmer and C. Lubich, A low-rank projector-splitting integrator for the Vlasov–Poisson equation, SIAM J. Sci. Comput., 40 (2018), pp. B1330–B1360,https://doi.org/10.1137/18M116383X

    2018 doi

  11. [19]

    Einkemmer, A

    L. Einkemmer, A. Ostermann, and C. Piazzola, A low-rank projector-splitting integrator for the Vlasov–Maxwell equations with divergence correction , J. Comput. Phys., 403 (2020), p. 109063, https://doi.org/10.1016/j.jcp.2019.109063

    2020

  12. [20]

    Filbet and S

    F. Filbet and S. Jin, A class of asymptotic-preserving schemes for kinetic equations and related problems with stiff sources , J. Comput. Phys., 229 (2010), pp. 7625–7648

    2010

  13. [21]

    Gamba, J

    I. Gamba, J. Haack, C. Hauck, and J. Hu, A fast spectral method for the Boltzmann collision operator with general collision kernels , SIAM J. Sci. Comput., 39 (2017), pp. B658–B674

    2017

  14. [22]

    W. Guo, J. F. Ema, and J.-M. Qiu, A local Macroscopic Conservative (LoMaC) Low Rank Tensor Method with the Discontinuous Galerkin Method for the Vlasov Dynamics , Commun. Appl. Math. Comput., 6 (2024), pp. 550–575

    2024

  15. [23]

    J. Hu, S. Jin, and Q. Li, Asymptotic-preserving schemes for multiscale hyperbolic and kinetic equations, in Handbook of Numerical Methods for Hyperbolic Problems: Applied and Modern Issues, R. Abgrall and C.-W. Shu, eds., North-Holland, 2017, ch. 5, pp. 103–129

    2017

  16. [24]

    J. Hu, S. Jin, and R. Shu, On stochastic Galerkin approximation of the nonlinear Boltzmann equation with uncertainty in the fluid regime , Journal of Computational Physics, 397 (2019), p. 108838

    2019

  17. [25]

    Hu and Y

    J. Hu and Y. W ang, An adaptive dynamical low rank method for the nonlinear Boltzmann equation , Journal of Scientific Computing, 92 (2022), p. 75

    2022

  18. [26]

    Hu and Y

    J. Hu and Y. W ang, An adaptive dynamical low rank method for the nonlinear Boltzmann equation , J. Sci. Comput., 92 (2022), p. 75. 27

    2022

  19. [27]

    Jin, Asymptotic preserving (AP) schemes for multiscale kinetic and hyperbolic equations: a review , Riv

    S. Jin, Asymptotic preserving (AP) schemes for multiscale kinetic and hyperbolic equations: a review , Riv. Mat. Univ. Parma, 3 (2012), pp. 177–216

    2012

  20. [28]

    Kusch and P

    J. Kusch and P. Stammer, A robust collision source method for rank adaptive dynamical low- rank approximation in radiation therapy , ESAIM: Mathematical Modelling and Numerical Analysis, 57 (2023), pp. 865–891

    2023

  21. [29]

    Lubich, From quantum to classical molecular dynamics: reduced models and numerical analysis , European Mathematical Society, Z¨ urich, 2008

    C. Lubich, From quantum to classical molecular dynamics: reduced models and numerical analysis , European Mathematical Society, Z¨ urich, 2008

    2008

  22. [30]

    Lubich and I

    C. Lubich and I. V. Oseledets , A projector-splitting integrator for dynamical low-rank approximation, BIT Numerical Mathematics, 54 (2014), pp. 171–188

    2014

  23. [31]

    Meyer, F

    H.-D. Meyer, F. Gatti, and G. A. Worth, Multidimensional quantum dynamics , John Wiley & Sons, 2009

    2009

  24. [32]

    Mouhot and L

    C. Mouhot and L. Pareschi, Fast algorithms for computing the Boltzmann collision operator , Math. Comp., 75 (2006), pp. 1833–1852

    2006

  25. [33]

    Nakao, J.-M

    J. Nakao, J.-M. Qiu, and L. Einkemmer , Reduced Augmentation Implicit Low-rank (RAIL) integrators for advection-diffusion and Fokker-Planck models , arXiv:2311.15143, (2023)

    2023 arXiv

  26. [34]

    Pareschi and G

    L. Pareschi and G. Russo, Numerical solution of the Boltzmann equation I: spectrally accurate approximation of the collision operator , SIAM J. Numer. Anal., 37 (2000), pp. 1217–1245

    2000

  27. [35]

    Z. Peng, R. McClarren, and M. Frank, A low-rank method for two-dimensional time-dependent radiation transport calculations , J. Comput. Phys., 421 (2020), https://doi.org/10.1016/j.jcp. 2020.109735

    2020

  28. [36]

    V an Leer, Towards the ultimate conservative difference scheme

    B. V an Leer, Towards the ultimate conservative difference scheme. V. A second-order sequel to Godunov’s method, Journal of computational Physics, 32 (1979), pp. 101–136. 28

    1979

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.