A discrete Boltzmann model with state-dependent power-law relaxation time for nonequilibrium transport in compressible flows
Pith reviewed 2026-05-20 00:30 UTC · model grok-4.3
The pith
A discrete Boltzmann model with power-law relaxation time depending on density and temperature recovers nonequilibrium transport across shocks and discontinuities.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The DTRT-DBM incorporates a relaxation time of the form τ = τ₀ (ρ/ρ₀)^a (T/T₀)^b. This formulation extends the discrete Boltzmann framework to flows with spatially varying nonequilibrium intensity. Validation against the Sod shock tube and analytical solutions for viscous stress and heat flux shows that the model recovers macroscopic wave structures and nonequilibrium quantities across shock waves, rarefaction waves, and contact discontinuities. Phase diagrams of viscous stress and heat flux reveal exponential dependence of extrema on a and b, with asymmetric sensitivity: more to a for density-dominated gradients and more to b for temperature-dominated ones.
What carries the argument
The state-dependent relaxation time τ=τ₀(ρ/ρ₀)^a(T/T₀)^b inside the discrete Boltzmann collision operator, which adjusts the rate at which the distribution relaxes toward equilibrium according to local density and temperature.
Load-bearing premise
A simple power-law dependence on density and temperature alone is enough to represent how nonequilibrium intensity changes from place to place in the flows examined.
What would settle it
Measurements or high-resolution simulations of viscous stress and heat flux in a compressible flow featuring strong chemical reactions or ionization, where the model's predictions would deviate noticeably from the data if the power-law form fails to capture the true state dependence.
Figures
read the original abstract
Thermodynamic nonequilibrium effects play a central role in momentum and energy transport in compressible flows. In conventional BGK kinetic models, the relaxation time $\tau$ is taken as a constant, which neglects the dependence of the relaxation process on local macroscopic states. To overcome this limitation, we develop a discrete Boltzmann model with a density- and temperature-dependent power-law relaxation time, termed DTRT-DBM, in which $\tau=\tau_0(\rho/\rho_0)^a(T/T_0)^b$. This formulation extends the discrete Boltzmann framework to flows with spatially varying nonequilibrium intensity. The model is validated by the Sod shock tube and by analytical solutions for viscous stress and heat flux, demonstrating accurate recovery of both macroscopic wave structures and nonequilibrium quantities across shock waves, rarefaction waves, and contact discontinuities. On this basis, phase diagrams of viscous stress and heat flux are constructed to examine how these quantities depend on the power-law exponents $a$ and $b$. The extrema of these quantities depend exponentially on the model parameters and exhibit regime-dependent behaviour. The roles of $a$ and $b$ are not symmetric: the nonequilibrium response is more sensitive to $a$ when density gradients dominate, but more sensitive to $b$ when temperature gradients dominate. Within the parameter range and flow configurations examined here, higher-order viscous stress increases the growth rate of the total viscous-stress extremum, whereas higher-order heat flux reduces the growth rate of the total heat-flux extremum. These results show that the proposed model can capture different higher-order nonequilibrium responses in compressible flows and provides a framework for the modelling and analysis of multiscale nonequilibrium processes.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a discrete Boltzmann model (DTRT-DBM) that replaces the constant relaxation time of standard BGK-type models with a power-law form τ = τ₀(ρ/ρ₀)^a (T/T₀)^b. The model is tested on the Sod shock-tube problem and against analytical expressions for viscous stress and heat flux; the authors report accurate recovery of both macroscopic wave structures and nonequilibrium moments across shocks, rarefactions, and contacts. Phase diagrams are constructed to map the dependence of extremal viscous stress and heat flux on the exponents a and b, revealing asymmetric sensitivity (a more influential under density gradients, b under temperature gradients) and differing effects of higher-order moments.
Significance. If the variable-τ formulation is shown to be consistent with the discrete velocity discretization and the cited analytical benchmarks, the work would supply a compact, tunable extension of discrete Boltzmann methods for flows with spatially varying nonequilibrium intensity. The phase-diagram analysis supplies concrete, falsifiable predictions about how nonequilibrium extrema scale with a and b.
major comments (2)
- [Model formulation and Chapman-Enskog recovery] § on model formulation and Chapman-Enskog recovery: the power-law τ(ρ,T) necessarily produces additional terms proportional to ∇τ (and therefore to ∇ρ and ∇T) in the first-order stress and heat-flux expressions. These terms are absent from the constant-τ analytical solutions used for validation. The manuscript must either derive the modified transport coefficients explicitly or demonstrate that the discrete velocity set and collision operator cancel the extra contributions for arbitrary a,b; otherwise the reported agreement with analytical viscous stress/heat flux is limited to the hydrodynamic limit or to special (a,b) pairs.
- [Validation against analytical viscous stress and heat flux] Validation against analytical viscous stress and heat flux (Sod-shock-tube section): the strongest claim—that the model recovers nonequilibrium quantities across discontinuities—rests on direct comparison with analytical expressions that assume constant τ. Without an accompanying Chapman-Enskog verification that includes the ∇τ contributions, the quantitative match cannot be taken as general evidence for the power-law form.
minor comments (3)
- Specify the numerical values and ranges of a and b examined in the phase diagrams, together with any physical constraints (e.g., positivity of τ, realizability of moments).
- Add quantitative error measures (L2 norms, maximum deviations) for the Sod-shock-tube profiles rather than qualitative statements of agreement.
- Clarify whether the discrete velocity set is the same as in prior DBM papers or modified to accommodate the variable-τ operator.
Simulated Author's Rebuttal
We thank the referee for the thorough review and valuable comments on our manuscript. The points raised regarding the Chapman-Enskog recovery and validation are well-taken, and we have prepared detailed responses and revisions to address them.
read point-by-point responses
-
Referee: [Model formulation and Chapman-Enskog recovery] the power-law τ(ρ,T) necessarily produces additional terms proportional to ∇τ (and therefore to ∇ρ and ∇T) in the first-order stress and heat-flux expressions. These terms are absent from the constant-τ analytical solutions used for validation. The manuscript must either derive the modified transport coefficients explicitly or demonstrate that the discrete velocity set and collision operator cancel the extra contributions for arbitrary a,b; otherwise the reported agreement with analytical viscous stress/heat flux is limited to the hydrodynamic limit or to special (a,b) pairs.
Authors: We agree that the variable relaxation time introduces additional terms involving ∇τ in the Chapman-Enskog expansion. In the revised manuscript, we will derive the modified transport coefficients explicitly, accounting for these contributions. We will also examine whether the specific discrete velocity discretization leads to any cancellations for the chosen power-law form. This will provide a more rigorous foundation for the validation against analytical expressions. revision: yes
-
Referee: [Validation against analytical viscous stress and heat flux] the strongest claim—that the model recovers nonequilibrium quantities across discontinuities—rests on direct comparison with analytical expressions that assume constant τ. Without an accompanying Chapman-Enskog verification that includes the ∇τ contributions, the quantitative match cannot be taken as general evidence for the power-law form.
Authors: We acknowledge the limitation in the original validation. We will include a detailed Chapman-Enskog analysis that incorporates the ∇τ terms in the revised manuscript. The comparisons in the Sod shock-tube section will be updated to use the corrected analytical expressions for viscous stress and heat flux, thereby providing stronger evidence for the power-law relaxation time model. revision: yes
Circularity Check
No significant circularity; modeling choice validated against external benchmarks
full rationale
The paper introduces the state-dependent relaxation time τ=τ₀(ρ/ρ₀)^a(T/T₀)^b explicitly as a modeling extension to conventional constant-τ BGK models, rather than deriving it from prior equations or fitting it to the target nonequilibrium quantities. Validation proceeds via direct comparison to the independent Sod shock-tube problem and to analytical expressions for viscous stress and heat flux; these benchmarks lie outside the model construction and are not obtained by re-expressing the same inputs. No load-bearing step reduces a claimed prediction to a self-citation chain, a fitted parameter renamed as output, or a definitional equivalence. The derivation chain is therefore self-contained against external references.
Axiom & Free-Parameter Ledger
free parameters (2)
- a
- b
axioms (1)
- domain assumption BGK-type single-relaxation-time collision operator remains valid when relaxation time is made state-dependent.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
τ = τ₀(ρ/ρ₀)^a (T/T₀)^b ... Chapman–Enskog multiscale analysis ... Δ∗(1)₂ = −μ[∇u + (∇u)ᵀ − 2/(n+2)I∇·u]
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
D2V25 discrete velocity set ... moment constraints M0–M5,3
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
- [1]
-
[2]
C. D. Lin and K. H. Luo , journal=. 2019 , publisher=
work page 2019
-
[3]
M. Perucho and M. Hanasz and J. M. Mart. Phys. Rev. E , volume=. 2007 , publisher=
work page 2007
-
[4]
H. Yang and Z. Xu and E. K. Lim and S. Kim and K. S. Cho and Y. H. Kim and J. Chae and K. Cho and K. Ji , journal=. 2018 , publisher=
work page 2018
-
[5]
V. N. Gamezo and A. M. Khokhlov and E. S. Oran and A. Y. Chtchelkanova and R. O. Rosenberg , journal=. 2003 , publisher=
work page 2003
- [6]
-
[7]
S. Boccelli and W. Kaufmann and T. E. Magin and J. G. McDonald , journal=. 2024 , publisher=
work page 2024
-
[8]
F. Chen and A. G. Xu and J. H. Song and Y. B. Gan and Y. D. Zhang and N. Guan , journal=. 2024 , publisher=
work page 2024
-
[9]
H. B. Ding and Y. Zhang and C. Q. Sun and Y. Yang and C. Wen , journal=. 2022 , publisher=
work page 2022
-
[10]
M. Bernardini and D. Modesti and F. Salvadore and S. Sathyanarayana and G. Della Posta and S. Pirozzoli , journal=. 2023 , publisher=
work page 2023
- [11]
-
[12]
M. Kl. J. Chem. Theory Comput. , volume=. 2018 , publisher=
work page 2018
-
[13]
Y. X. Li and M. Hekmatifar and Y. L. Sun and others , journal=
- [14]
- [15]
- [16]
-
[17]
Y. J. Zhu and C. W. Zhong and K. Xu , journal=. 2017 , publisher=
work page 2017
-
[18]
Z. L. Guo and K. Xu , journal=. 2021 , publisher=
work page 2021
-
[19]
L. M. Yang and C. Shu and W. M. Yang and J. Wu , journal=. 2019 , publisher=
work page 2019
-
[20]
L. M. Yang and C. Shu and J. Wu and Y. Y. Liu and X. Shen , journal=. 2022 , publisher=
work page 2022
-
[21]
S. B. Q. Tran and F. Y. Leong and Q. T. Le and D. V. Le , journal=. 2022 , publisher=
work page 2022
-
[22]
J. Latt and C. Coreixas and J. B. Phil. Trans. R. Soc. A , volume=. 2020 , publisher=
work page 2020
-
[23]
C. J. Zhan and X. Liu and Z. H. Chai and B. C. Shi , journal=. 2025 , publisher=
work page 2025
-
[24]
Z. D. Wang and Y. K. Wei and Y. H. Qian , journal=. 2022 , publisher=
work page 2022
-
[25]
L. L. Fei and D. Derome and J. Carmeliet , journal=. 2024 , publisher=
work page 2024
-
[26]
Y. Yang and J. Tu and M. L. Shan and Z. J. Zhang and C. Chen and H. X. Li , journal=. 2025 , publisher=
work page 2025
-
[27]
G. X. Hou and L. L. Fei and C. D. Lin and W. W. Yan , journal=. 2025 , publisher=
work page 2025
-
[28]
D. J. Zhang and A. G. Xu and Y. D. Zhang and Y. B. Gan and Y. J. Li , journal=. 2022 , publisher=
work page 2022
-
[29]
G. L. Sun and Y. B. Gan and A. G. Xu and Q. F. Shi , journal=. 2024 , publisher=
work page 2024
-
[30]
J. H. Song and A. G. Xu and L. Miao and F. Chen and Z. P. Liu and L. F. Wang and N. F. Wang and X. Hou , journal=. 2024 , publisher=
work page 2024
-
[31]
J. H. Song and L. Miao and F. Chen and Y. B. Gan and A. G. Xu and L. X. Li , journal=. 2025 , publisher=
work page 2025
-
[32]
Y. F. Li and C. D. Lin , journal=. 2024 , publisher=
work page 2024
- [33]
-
[34]
Y. K. Wei and H. Yang and H. S. Dou and Z. Lin and Z. D. Wang and Y. H. Qian , journal=. 2018 , publisher=
work page 2018
-
[35]
B. E. Morgan and B. J. Olson and J. E. White and J. A. McFarland , journal=. 2017 , publisher=
work page 2017
-
[36]
H. S. Zhang and R. Betti and R. Yan and H. Aluie , journal=. 2020 , publisher=
work page 2020
-
[37]
Z. D. Wang and Y. K. Wei and Y. H. Qian , journal=. 2020 , publisher=
work page 2020
-
[38]
N. Sawant and B. Dorschner and I. V. Karlin , journal=. 2022 , publisher=
work page 2022
-
[39]
H. Y. Huang and K. Jin and K. Li and H. Li and X. J. Zheng , journal=. 2024 , publisher=
work page 2024
-
[40]
L. Vienne and A. Giauque and E. L. Phys. Fluids , volume=. 2024 , publisher=
work page 2024
-
[41]
S. A. Hosseini and P. Boivin and D. Th. Prog. Energy Combust. Sci. , volume=. 2024 , publisher=
work page 2024
-
[42]
F. Chen and J. H. Song and Y. B. Gan and A. G. Xu and W. K. He , journal=. 2025 , publisher=
work page 2025
-
[43]
W. H. Huang and C. D. Lin and X. L. Su and J. Li , journal=. 2025 , publisher=
work page 2025
- [44]
- [45]
-
[46]
H. W. Li and A. G. Xu and G. Zhang and Y. M. Shan , journal=. 2022 , publisher=
work page 2022
- [47]
-
[48]
C. D. Lin and K. H. Luo and A. G. Xu and Y. B. Gan and H. L. Lai , journal=. 2021 , publisher=
work page 2021
-
[49]
C. D. Lin , booktitle=. 2023 , organization=
work page 2023
-
[50]
Q. H. Guo and Y. B. Gan and B. Yang and Y. H. Wu and H. L. Lai and A. G. Xu , journal=. 2025 , publisher=
work page 2025
-
[51]
Z. Y. He and Y. B. Gan and B. Yang and D. M. Li and H. L. Lai and A. G. Xu , journal=. 2025 , publisher=
work page 2025
-
[52]
Y. H. Wu and Y. B. Gan and A. G. Xu and B. Yang , journal=. 2025 , publisher=
work page 2025
-
[53]
Q. B. Wu and C. D. Lin and H. L. Lai , journal=. 2025 , publisher=
work page 2025
-
[54]
A. G. Xu and G. C. Zhang and Y. B. Gan and F. Chen and X. J. Yu , journal=. 2012 , publisher=
work page 2012
-
[55]
A. G. Xu and G. C. Zhang and Y. D. Zhang , journal=. 2018 , publisher=
work page 2018
-
[56]
A. G. Xu and Y. D. Zhang , year=
-
[57]
A. G. Xu and D. J. Zhang and Y. B. Gan , journal=. 2024 , publisher=
work page 2024
-
[58]
Y. B. Gan and A. G. Xu and G. C. Zhang and Y. D. Zhang and S. Succi , journal=. 2018 , publisher=
work page 2018
-
[59]
Y. B. Gan and A. G. Xu and G. C. Zhang and C. D. Lin and H. L. Lai and Z. P. Liu , journal=. 2019 , publisher=
work page 2019
- [60]
-
[61]
L. F. Wang and W. H. Ye and Y. J. Li , journal=
-
[62]
L. F. Wang and W. H. Ye and W.-S. Don and Z. M. Sheng and Y. J. Li and X. T. He , journal=. 2010 , publisher=
work page 2010
-
[63]
L. F. Wang and W. H. Ye and X. T. He , journal=. 2012 , publisher=
work page 2012
-
[64]
J. D. Sadler and S. Green and S. T. Li and Y. Zhou and K. A. Flippo and H. Li , journal=. 2022 , publisher=
work page 2022
-
[65]
C. Almarcha and P. Clavin and L. Duchemin and J. Sanz , journal=. 2007 , publisher=
work page 2007
-
[66]
L. F. Wang and W. H. Ye and W. Y. Zhang and X. T. He , journal=. 2013 , publisher=
work page 2013
-
[67]
R. C. Millikan and D. R. White , journal=. 1963 , publisher=
work page 1963
-
[68]
J. W. Streicher and A. Krish and R. K. Hanson , journal=. 2020 , publisher=
work page 2020
-
[69]
W. Q. Hu and Z. H. Li and A. P. Peng and X. Y. Jiang , journal=
- [70]
-
[71]
Z.H. Li and W.Q. Hu and A.P. Peng and J.L. Wu and C.H. Lee , journal=. 2020 , publisher=
work page 2020
-
[72]
J. L. Wu and Z. H. Li and Z. B. Zhang and A. P. Peng , journal=. 2021 , publisher=
work page 2021
-
[73]
L. H. Holway , journal=. 1966 , publisher=
work page 1966
-
[74]
E. M. Shakhov , journal=. 1968 , publisher=
work page 1968
-
[75]
G. J. Liu , journal=. 1990 , publisher=
work page 1990
-
[76]
V. A. Rykov and V. A. Titarev and E. M. Shakhov , journal=. 2008 , publisher=
work page 2008
-
[77]
J. H. Shi and B. Zhao and X. X. Niu and Q. Xin and H. Z. Xu and X. Q. Lu , journal=. 2024 , publisher=
work page 2024
-
[78]
Y. B. Gan and A. G. Xu and G. C. Zhang and H. L. Lai , journal=. 2018 , publisher=
work page 2018
-
[79]
G. S. Jiang and C. W. Shu , journal=. 1996 , publisher=
work page 1996
-
[80]
U. M. Ascher and S. J. Ruuth and R. J. Spiteri , journal=. 1997 , publisher=
work page 1997
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.