pith. sign in

arxiv: 2407.05894 · v2 · submitted 2024-07-08 · 🧮 math.NA · cs.NA· math.AP

On Nonlinear Closures for Moment Equations Based on Orthogonal Polynomials

Eda Yilmaz , Georgii Oblapenko , Manuel Torrilhon This is my paper

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

classification 🧮 math.NA cs.NAmath.AP
keywords moment closureorthogonal polynomialsGram matrixgas kinetic theoryhyperbolicityrealizabilityGrad closuremaximum entropy closure
0
0 comments X

The pith

The Gramian closure, based on orthogonal polynomials from Gram matrices, provides accurate results with attractive properties for moment equations in gas kinetic theory.

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

This paper introduces a method for closing moment equations by constructing orthogonal polynomials from the Gram matrices associated with the moment basis. The approach is applied to the moment closure problem in gas kinetic theory, where it is shown to possess several desirable mathematical properties. Numerical experiments on model particle distributions demonstrate that the Gramian closure yields highly accurate approximations across a broad range of cases. Comparisons are made with established techniques such as Grad's closure and the maximum-entropy closure. Readers interested in reducing kinetic models to fluid-like equations while maintaining fidelity would find this relevant because it offers both theoretical advantages and practical accuracy.

Core claim

An approach to the moment closure problem is proposed on the basis of orthogonal polynomials derived from Gram matrices. In the context of gas kinetic theory, this Gramian closure is proven to have multiple attractive mathematical properties. Numerical studies for model gas particle distributions show that it provides very accurate results for a wide range of distribution functions when compared to Grad's closure and the maximum-entropy method.

What carries the argument

The Gramian closure, defined by using orthogonal polynomials constructed from Gram matrices of the moment basis to determine the closure relation.

If this is right

  • It provides accurate approximations for a wide range of distribution functions.
  • It possesses multiple attractive mathematical properties such as preserving hyperbolicity and realizability.
  • It does not introduce new instabilities for distributions arising in kinetic theory.
  • It serves as a competitive alternative to Grad's closure and the maximum-entropy method.

Where Pith is reading between the lines

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

  • The method could be adapted to moment problems in other physical systems beyond gases, such as in radiation transport or semiconductor modeling.
  • Further analysis might reveal how the choice of basis affects the closure's performance in high-dimensional cases.
  • Testing the closure in full simulations of shock waves or boundary layers could validate its practical utility.

Load-bearing premise

The polynomials derived from the Gram matrices will always produce a closure that maintains the necessary hyperbolicity and realizability properties for the relevant distribution functions without causing instabilities.

What would settle it

Observation of a specific distribution function in kinetic theory for which the Gramian closure leads to unphysical negative densities or oscillatory instabilities in the moment equations, unlike the maximum-entropy method.

Figures

Figures reproduced from arXiv: 2407.05894 by Eda Yilmaz, Georgii Oblapenko, Manuel Torrilhon.

Figure 1
Figure 1. Figure 1: Family of distribution functions for test problems in subsection [PITH_FULL_IMAGE:figures/full_fig_p019_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Relative error er(uclosure) of the next higher moment calculated for different number of moments using different closure techniques is shown for the Mott-Smith shock wave distribution [39]. We consider Ma = 4, γ = 5/3 and compute moments with the distribution fMS for each x values from -10 to 10 with step size 0.25. For the maximum entropy closure, we discretize velocity domain [−6, 9] with 1000 points. No… view at source ↗
Figure 3
Figure 3. Figure 3: Relative error between uclosure and uM+1 with fixed parameters v0 = 1.5, β = −0.05 for the electron hole distribution [5] with different electrostatic potentials 0 ≤ ϕ ≤ 2 with step size 0.04. Even (left) and odd (right) case scenarios were examined and are shown separately due to differences in the definition of the Gramian and extended Gramian closures. For the maximum entropy closure, the velocity domai… view at source ↗
Figure 4
Figure 4. Figure 4: Relative error er(uclosure) of the next higher moment for a range of widths w in the test of the realizability boundaries. The left part of the figure shows the relative error for even numbers of moments M = 4, 6 while the right one shows odd cases M = 5, 7. For the maximum entropy closure, computational domain is taken with c ∈ [−4, 5] with 1000 grid points. of distribution function fEH can be seen in the… view at source ↗
Figure 5
Figure 5. Figure 5: Condition numbers of different closure methods as a function of the order of the moment [PITH_FULL_IMAGE:figures/full_fig_p023_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: The relative errors of the next higher moment for the Mott-Smith distribution from subsection [PITH_FULL_IMAGE:figures/full_fig_p024_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: The relative error of different closure methods for the electron-hole test problem in subsection [PITH_FULL_IMAGE:figures/full_fig_p026_7.png] view at source ↗
read the original abstract

In the present work, an approach to the moment closure problem on the basis of orthogonal polynomials derived from Gram matrices is proposed. Its properties are studied in the context of the moment closure problem arising in gas kinetic theory, for which the proposed approach is proven to have multiple attractive mathematical properties. Numerical studies are carried out for model gas particle distributions and the approach is compared to other moment closure methods, such as Grad's closure and the maximum-entropy method. The proposed ``Gramian'' closure is shown to provide very accurate results for a wide range of distribution functions.

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 paper proposes a Gramian closure for the moment closure problem in gas kinetic theory, constructed via orthogonal polynomials derived from the Gram matrix of the moment basis. It proves that the resulting closure preserves realizability, yields a hyperbolic system, and satisfies entropy dissipation (Sections 3–4). Numerical experiments in Section 5 compare the method to Grad’s closure and the maximum-entropy closure on several model distributions, reporting high accuracy across a wide range of test cases.

Significance. If the stated proofs and numerical results hold, the work supplies a new nonlinear closure with multiple mathematically attractive properties (realizability preservation, hyperbolicity, entropy dissipation) that are explicitly constructed and verified in Sections 3–4, together with reproducible numerical evidence in Section 5. This combination of theoretical guarantees and demonstrated accuracy on kinetic-theory distributions constitutes a substantive contribution to moment methods.

minor comments (2)
  1. [§2] §2: the precise definition of the Gram matrix and the associated inner product could be restated with an explicit formula to aid readers who are not already familiar with the construction.
  2. [§5] §5: the tables and figures would benefit from a brief statement of the precise error norm and the range of Knudsen numbers or moments retained in each test.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript, for the positive summary and significance assessment, and for the recommendation to accept. We are pleased that the combination of theoretical guarantees (realizability preservation, hyperbolicity, and entropy dissipation) and the numerical comparisons were viewed as a substantive contribution.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper constructs the Gramian closure explicitly from the Gram matrix of the moment basis using orthogonal polynomials, then derives and proves properties (realizability preservation, hyperbolicity, entropy dissipation) in Sections 3–4 via direct mathematical arguments on the resulting system. Numerical comparisons in Section 5 are performed against independent baselines (Grad, maximum-entropy) on model distributions without any fitted parameters or self-referential definitions. No load-bearing step reduces to a self-citation, fitted input renamed as prediction, or ansatz smuggled via prior work; the central claims rest on an independent construction whose proofs are self-contained within the manuscript.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review provides no explicit free parameters, axioms, or invented entities; all details of the Gram-matrix construction and proofs are absent.

pith-pipeline@v0.9.0 · 5625 in / 1007 out tokens · 19793 ms · 2026-05-23T23:01:06.547456+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

52 extracted references · 52 canonical work pages · 2 internal anchors

  1. [1]

    M. R. Abdelmalik, Z. Cai, and T. Pichard , Moment methods for the radiative trans- fer equation based on φ-divergences, Comput. Methods Appl. Mech. Eng., 417 (2023), p. 116454

    work page 2023

  2. [2]

    Aheizer and M

    N. Aheizer and M. Krein , Some Questions in the Theory of Moments , vol. 2 of Trans- lations of Mathematical Monographs, American Mathematical Society, 1962

    work page 1962

  3. [3]

    G. W. Alldredge, C. D. Hauck, and A. L. Tits , High-order entropy-based closures for linear transport in slab geometry II: A computational study of the optimization problem , SIAM Journal on Scientific Computing, 34 (2012), pp. B361–B391

    work page 2012

  4. [4]

    B¨ohmer and M

    N. B¨ohmer and M. Torrilhon , Entropic quadrature for moment approximations of the Boltzmann-BGK equation, J. Comput. Phys., 401 (2020), p. 108992

    work page 2020

  5. [5]

    Bujarbarua and H

    S. Bujarbarua and H. Schamel , Theory of finite-amplitude electron and ion holes , Journal of Plasma Physics, 25 (1981), pp. 515–529

    work page 1981

  6. [6]

    Z. Cai, Y. Fan, and R. Li, A framework on moment model reduction for kinetic equation , SIAM Journal on Applied Mathematics, 75 (2015), pp. 2001–2023

    work page 2015

  7. [7]

    Cercignani , The Boltzmann Equation and its Applications , Applied Mathematical Sciences, Springer, New York, 1988

    C. Cercignani , The Boltzmann Equation and its Applications , Applied Mathematical Sciences, Springer, New York, 1988

    work page 1988

  8. [8]

    Chalons, R

    C. Chalons, R. Fox, and M. Massot, A multi-gaussian quadrature method of moments for gas-particle flows in a les framework , in Proceedings of the Summer Program, 2010, pp. 347–358

    work page 2010

  9. [9]

    Beyond pressureless gas dynamics: Quadrature-based velocity moment models

    C. Chalons, D. Kah, and M. Massot , Beyond pressureless gas dynamics: quadrature- based velocity moment models, arXiv preprint arXiv:1011.2974, (2010)

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  10. [10]

    Curto, Recursiveness, positivity and truncated moment problems , Houston Journal of Mathematics, 17 (1991), pp

    R. Curto, Recursiveness, positivity and truncated moment problems , Houston Journal of Mathematics, 17 (1991), pp. 603–635

    work page 1991

  11. [11]

    Desjardins, R

    O. Desjardins, R. O. Fox, and P. Villedieu , A quadrature-based moment method for dilute fluid-particle flows , Journal of Computational Physics, 227 (2008), pp. 2514–2539

    work page 2008

  12. [12]

    Dreyer, Maximisation of the entropy in non-equilibrium , Journal of Physics A: Math- ematical and General, 20 (1987), p

    W. Dreyer, Maximisation of the entropy in non-equilibrium , Journal of Physics A: Math- ematical and General, 20 (1987), p. 6505. 28 Yilmaz, Oblapenko, Torrilhon (2024) Nonlinear Closures

    work page 1987

  13. [13]

    R. O. Fox , A quadrature-based third-order moment method for dilute gas-particle flows , Journal of Computational Physics, 227 (2008), pp. 6313–6350

    work page 2008

  14. [14]

    R. O. Fox and F. Laurent , Hyperbolic quadrature method of moments for the one- dimensional kinetic equation , SIAM J. Appl. Math., 82 (2022), pp. 750–771

    work page 2022

  15. [15]

    R. O. Fox, F. Laurent, and A. Passalacqua , The generalized quadrature method of moments, Journal of Aerosol Science, 167 (2023), p. 106096

    work page 2023

  16. [16]

    R. O. Fox, F. Laurent, and A. Vi ´e, Conditional hyperbolic quadrature method of mo- ments for kinetic equations , Journal of Computational Physics, 365 (2018), pp. 269–293

    work page 2018

  17. [17]

    Frank, B

    M. Frank, B. Dubroca, and A. Klar , Partial moment entropy approximation to ra- diative heat transfer, J. Comput. Phys., 218 (2006), pp. 1–18

    work page 2006

  18. [18]

    Gautschi , On generating orthogonal polynomials , SIAM J

    W. Gautschi , On generating orthogonal polynomials , SIAM J. Sci. Stat. Comput., 3 (1982), pp. 289–317

    work page 1982

  19. [19]

    Gautschi, Orthogonal polynomials: applications and computation , Acta numerica, 5 (1996), pp

    W. Gautschi, Orthogonal polynomials: applications and computation , Acta numerica, 5 (1996), pp. 45–119

    work page 1996

  20. [20]

    C. S. Gillespie, Moment-closure approximations for mass-action models , IET Syst. Biol., 3 (2009), pp. 52–58

    work page 2009

  21. [21]

    Grad , On the kinetic theory of rarefied gases , Comm

    H. Grad , On the kinetic theory of rarefied gases , Comm. Pure Appl. Math., 2 (1949), pp. 331–407

    work page 1949

  22. [22]

    H. L. Hamburger, Hermitian transformations of deficiency-index (1, 1), Jacobi matrices and undetermined moment problems, American Journal of Mathematics, 66 (1944), pp. 489– 522

    work page 1944

  23. [23]

    Huang, S

    Q. Huang, S. Li, and W.-A. Yong, Stability analysis of quadrature-based moment meth- ods for kinetic equations , SIAM Journal on Applied Mathematics, 80 (2020), pp. 206–231

    work page 2020

  24. [24]

    Junk , Domain of definition of Levermore’s five-moment system , J

    M. Junk , Domain of definition of Levermore’s five-moment system , J. Stat. Phys., 93 (1998), pp. 1143–1167

    work page 1998

  25. [25]

    D. S. Kershaw, Flux limiting natures own way – a new method for numerical solution of the transport equation , tech. rep., Lawrence Livermore National Lab.(LLNL), Livermore, CA (United States), 1976

    work page 1976

  26. [26]

    Koellermeier, R

    J. Koellermeier, R. P. Schaerer, and M. Torrilhon , A framework for hyperbolic approximation of kinetic equations using quadrature-based projection methods., Kinet. Relat. Models, 7 (2014)

    work page 2014

  27. [27]

    M. N. Kogan and M. N. Kogan , The equations of the kinetic theory of gases , Rarefied Gas Dynamics, (1969), pp. 29–104

    work page 1969

  28. [28]

    Kowalewski, Einf¨ uhrung in die Determinantentheorie, 1925

    G. Kowalewski, Einf¨ uhrung in die Determinantentheorie, 1925

    work page 1925

  29. [29]

    Kuehn, Moment closure – a brief review, in Control of self-organizing nonlinear systems, E

    C. Kuehn, Moment closure – a brief review, in Control of self-organizing nonlinear systems, E. Sch¨ oll, S. Klapp, and P. H¨ ovel, eds., Springer, 2016, pp. 253–271

    work page 2016

  30. [30]

    Laplante and C

    J. Laplante and C. Groth , Comparison of maximum entropy and quadrature-based moment closures for shock transitions prediction in one-dimensional gaskinetic theory , in AIP Conference Proceedings, vol. 1786, AIP Publishing, 2016

    work page 2016

  31. [31]

    C. D. Levermore , Moment closure hierarchies for kinetic theories , J. Stat. Phys., 83 (1996), pp. 1021–1065

    work page 1996

  32. [32]

    , Entropy-based moment closures for kinetic equations , Transport Theory and Statis- tical Physics, 26 (1997), pp. 591–606

    work page 1997

  33. [33]

    McDonald and M

    J. McDonald and M. Torrilhon, Affordable robust moment closures for CFD based on the maximum-entropy hierarchy, J. Comput. Phys., 251 (2013), pp. 500–523

    work page 2013

  34. [34]

    McGraw , Description of aerosol dynamics by the quadrature method of moments , Aerosol science and Technology, 27 (1997), pp

    R. McGraw , Description of aerosol dynamics by the quadrature method of moments , Aerosol science and Technology, 27 (1997), pp. 255–265. 29 Nonlinear Closures Yilmaz, Oblapenko, Torrilhon (2024)

    work page 1997

  35. [35]

    Milbrandt and M

    J. Milbrandt and M. Yau , A multimoment bulk microphysics parameterization. part II: A proposed three-moment closure and scheme description , J. Atmos. Sci., 62 (2005), pp. 3065–3081

    work page 2005

  36. [36]

    Mitrovi´c and M

    D. Mitrovi´c and M. Nedeljkov , Delta shock waves as a limit of shock waves , J. Hy- perbolic Differ. Equ., 4 (2007), pp. 629–653

    work page 2007

  37. [37]

    Higher order minimum entropy approximations in radiative transfer

    P. Monreal and M. Frank, Higher order minimum entropy approximations in radiative transfer, arXiv preprint arXiv:0812.3063, (2008)

    work page internal anchor Pith review Pith/arXiv arXiv 2008

  38. [38]

    , Moment realizability and Kershaw closures in radiative transfer , 2013

    work page 2013

  39. [39]

    H. M. Mott-Smith , The solution of the boltzmann equation for a shock wave , Physical Review, 82 (1951), p. 885

    work page 1951

  40. [40]

    Pichard, A moment closure based on a projection on the boundary of the realizability domain: 1d case , Kinetic and related models, 13 (2020), pp

    T. Pichard, A moment closure based on a projection on the boundary of the realizability domain: 1d case , Kinetic and related models, 13 (2020), pp. 1243–1280

    work page 2020

  41. [41]

    , A moment closure based on a projection on the boundary of the realizability domain: Extension and analysis , Kinetic and Related Models, 15 (2022), p. 793

    work page 2022

  42. [42]

    C. R. J. Roger A. Horn , Matrix Analysis, Cambridge University Press, 2013

    work page 2013

  43. [43]

    Schm¨udgen, The Moment Problem , Graduate Texts in Mathematics, Springer Inter- national Publishing, 2017

    K. Schm¨udgen, The Moment Problem , Graduate Texts in Mathematics, Springer Inter- national Publishing, 2017

    work page 2017

  44. [44]

    Schneider, Kershaw closures for linear transport equations in slab geometry I: Model derivation, Journal of Computational Physics, 322 (2016), pp

    F. Schneider, Kershaw closures for linear transport equations in slab geometry I: Model derivation, Journal of Computational Physics, 322 (2016), pp. 905–919

    work page 2016

  45. [45]

    Shohat and J

    J. Shohat and J. Tamarkin, The problem of moments , AMS, Providence, 1945

    work page 1945

  46. [46]

    Struchtrup, Macroscopic Transport Equations for Rarefied Gas Flows, Interaction of Mechanics and Mathematics, Springer, New York, 2005

    H. Struchtrup, Macroscopic Transport Equations for Rarefied Gas Flows, Interaction of Mechanics and Mathematics, Springer, New York, 2005

    work page 2005

  47. [47]

    Szeg¨o, Orthogonal Polynomials, AMS, Providence, 1939

    G. Szeg¨o, Orthogonal Polynomials, AMS, Providence, 1939

    work page 1939

  48. [48]

    P.-Y. C. Taunay and M. E. Mueller , Quadrature-based moment methods for kinetic plasma simulations, Journal of Computational Physics, 473 (2023), p. 111700

    work page 2023

  49. [49]

    Torrilhon, Modeling nonequilibrium gas flow based on moment equations , Annu

    M. Torrilhon, Modeling nonequilibrium gas flow based on moment equations , Annu. Rev. Fluid Mech., 48 (2016), pp. 429–458

    work page 2016

  50. [50]

    Van Cappellen, M

    M. Van Cappellen, M. R. Vetrano, and D. Laboureur , Higher order hyperbolic quadrature method of moments for solving kinetic equations , Journal of Computational Physics, 436 (2021), p. 110280

    work page 2021

  51. [51]

    C. Yuan, F. Laurent, and R. Fox , An extended quadrature method of moments for population balance equations, J. Aerosol Sci., 51 (2012), pp. 1–23

    work page 2012

  52. [52]

    Zhang, Y

    R. Zhang, Y. Chen, Q. Huang, and W.-A. Yong , Dissipativeness of the hyper- bolic quadrature method of moments for kinetic equations, arXiv preprint arXiv:2406.13931, (2024). 30

    work page arXiv 2024