pith. sign in

arxiv: 2601.03472 · v2 · submitted 2026-01-07 · ❄️ cond-mat.soft · cond-mat.stat-mech

Kinetic theory of dilute weakly charged granular gases with hard-core and inverse power-law interactions under uniform shear flow

Yuria Kobayashi , Makoto R. Kikuchi , Shunsuke Iizuka , Satoshi Takada This is my paper

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

classification ❄️ cond-mat.soft cond-mat.stat-mech
keywords kinetic theorygranular gasuniform shear flowinverse power lawGrad moment methodrheologyrestitution coefficientDSMC
0
0 comments X

The pith

A kinetic theory using Grad's moment expansion accurately describes the rheology of dilute granular gases with inverse power-law interactions under uniform shear flow.

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

The paper develops a framework to study the steady rheology of dilute gases with repulsive hard-core and inverse power-law potentials under uniform shear flow. It starts from the Boltzmann equation incorporating a restitution coefficient that varies with impact velocity and derives evolution equations for the stress tensor using Grad's moment expansion. Collisional rates and transport coefficients are then approximated by simple analytical functions of temperature over wide shear rate ranges. These expressions show excellent quantitative agreement with direct simulation Monte Carlo results for key quantities like shear stress, temperature anisotropy, and shear viscosity. The analysis further indicates that the velocity distribution function stays close to Maxwellian even at strong shear rates.

Core claim

Starting from the Boltzmann equation for particles with a velocity-dependent restitution coefficient determined by the hard-core plus inverse power-law potential, Grad's moment method yields closed evolution equations for the stress tensor components. The collisional contributions are evaluated and fitted to analytical forms that capture the temperature dependence, producing explicit predictions for the shear stress and viscosity as functions of the shear rate. These predictions match the steady-state values from direct simulation Monte Carlo simulations to high accuracy across the tested range of parameters.

What carries the argument

Grad's moment expansion applied to the Boltzmann equation with a velocity-dependent restitution coefficient for the combined hard-core and inverse power-law interaction potential.

If this is right

  • The shear viscosity and normal stress differences can be computed from simple analytical expressions fitted to the collisional rates.
  • Temperature anisotropy in the velocity distribution is quantitatively predicted by the moment equations.
  • The system maintains a nearly Maxwellian velocity distribution under strong uniform shear.
  • Transport coefficients show consistent temperature dependence that holds over a broad range of shear rates.

Where Pith is reading between the lines

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

  • The approach may extend naturally to other interaction potentials by adjusting the restitution coefficient form.
  • Confirmation of near-Maxwellian behavior suggests that simpler hydrodynamic descriptions could work well for these dilute systems under shear.
  • Similar kinetic-theory derivations could be applied to study transient flows or more complex geometries.
  • The analytical fits might facilitate incorporation into larger-scale simulations of granular flows.

Load-bearing premise

Grad's moment expansion provides a sufficiently accurate representation of the velocity distribution function under strong shear, and the model for the velocity-dependent restitution coefficient correctly represents the physics of the inverse power-law plus hard-core potential.

What would settle it

Observation of significant deviations between the predicted shear viscosity and direct simulation results at high shear rates, particularly in the temperature anisotropy or stress components, would indicate that the theory fails to capture the distribution function accurately.

Figures

Figures reproduced from arXiv: 2601.03472 by Makoto R. Kikuchi, Satoshi Takada, Shunsuke Iizuka, Yuria Kobayashi.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic of our system. Monodisperse particles are [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The interparticle potential used in this paper. The [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Dependence of the collision angle [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Temperature dependence of Ω [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Shear rate dependence of (left) the temperature, (mi [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Shear rate dependence of the the viscosity for [PITH_FULL_IMAGE:figures/full_fig_p006_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Comparison between Ω [PITH_FULL_IMAGE:figures/full_fig_p006_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Relative errors of the viscosity obtained from the [PITH_FULL_IMAGE:figures/full_fig_p007_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Deviation of the dimensionless velocity distribut [PITH_FULL_IMAGE:figures/full_fig_p007_10.png] view at source ↗
read the original abstract

We develop a kinetic-theory framework to investigate the steady rheology of a dilute gas interacting via a repulsive potential under uniform shear flow. Starting from the Boltzmann equation with a restitution coefficient that depends on the impact velocity and potential strength, we derive evolution equations for the stress tensor based on Grad's moment expansion. The resulting expressions for the collisional rates and transport coefficients are fitted with simple analytical functions that capture their temperature dependence over a wide range of shear rates. Comparison with direct simulation Monte Carlo (DSMC) results shows excellent quantitative agreement for the shear stress, temperature anisotropy, and shear viscosity. We also analyze the velocity distribution functions, revealing that the system remains nearly Maxwellian even under strong shear.

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 / 1 minor

Summary. The manuscript develops a kinetic-theory framework for the steady rheology of dilute granular gases with hard-core plus inverse power-law interactions under uniform shear flow. Starting from the Boltzmann equation with a velocity-dependent restitution coefficient, the authors close the moment hierarchy at the stress-tensor level via Grad's expansion, derive evolution equations for the stress tensor, and supply analytical fitting functions for the collisional rates and transport coefficients that capture their temperature dependence. Direct comparison with DSMC simulations is reported to show excellent quantitative agreement for shear stress, temperature anisotropy, and shear viscosity, while the velocity distribution function is stated to remain nearly Maxwellian even at strong shear.

Significance. If the central results hold, the work supplies a practical analytical route to the rheology of weakly charged granular gases with composite potentials, together with DSMC-validated fitting expressions that can be used in larger-scale modeling. The explicit validation against independent particle simulations for multiple rheological observables is a clear strength.

major comments (2)
  1. [Velocity distribution function analysis and collisional-rate fitting] The Grad closure assumes that higher-order Hermite coefficients remain negligible when evaluating the velocity-dependent restitution integrals. The manuscript states that the VDF remains nearly Maxwellian and reports excellent DSMC agreement for second-order moments, yet provides no explicit comparison of fourth-order moments or the high-velocity tails (which dominate inelastic collisions) between the truncated theory and the DSMC data. A direct check of these quantities, for example in the strong-shear regime, is needed to confirm that the fitted collisional rates do not systematically misrepresent the rheology.
  2. [Analytical fitting functions for collisional rates] The analytical fitting functions for the collisional rates contain free parameters calibrated to DSMC data. While this is acceptable for practical use, the manuscript should state the explicit functional forms, the fitting procedure, and the range of validity (shear rate and restitution parameters) in a dedicated section or table so that the expressions can be reproduced without re-fitting.
minor comments (1)
  1. [Introduction and model definition] Notation for the inverse-power-law exponent and the hard-core diameter should be introduced once and used consistently; a short table of symbols would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and the positive recommendation for minor revision. We address the two major comments point by point below.

read point-by-point responses

  1. Referee: [Velocity distribution function analysis and collisional-rate fitting] The Grad closure assumes that higher-order Hermite coefficients remain negligible when evaluating the velocity-dependent restitution integrals. The manuscript states that the VDF remains nearly Maxwellian and reports excellent DSMC agreement for second-order moments, yet provides no explicit comparison of fourth-order moments or the high-velocity tails (which dominate inelastic collisions) between the truncated theory and the DSMC data. A direct check of these quantities, for example in the strong-shear regime, is needed to confirm that the fitted collisional rates do not systematically misrepresent the rheology.

    Authors: We agree that an explicit check of fourth-order moments would provide stronger support for the validity of the Grad closure and the near-Maxwellian assumption. In the revised manuscript we will add a new figure comparing selected fourth-order moments (including those related to the high-velocity tails) obtained from DSMC simulations against the predictions of the truncated theory, with particular attention to the strong-shear regime. This addition will confirm that the deviations remain small enough that the fitted collisional rates are not systematically biased. revision: yes

  2. Referee: [Analytical fitting functions for collisional rates] The analytical fitting functions for the collisional rates contain free parameters calibrated to DSMC data. While this is acceptable for practical use, the manuscript should state the explicit functional forms, the fitting procedure, and the range of validity (shear rate and restitution parameters) in a dedicated section or table so that the expressions can be reproduced without re-fitting.

    Authors: We thank the referee for this useful suggestion. In the revised version we will insert a dedicated subsection (new Section 3.3) that lists the explicit analytical expressions for all fitting functions, describes the least-squares procedure used to calibrate the free parameters against the DSMC data, and provides a table summarizing the fitting coefficients together with the ranges of shear rate and restitution coefficient for which the fits are valid. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation from Boltzmann equation via Grad's expansion is independent and externally validated

full rationale

The paper starts from the Boltzmann equation with a velocity-dependent restitution coefficient, applies Grad's moment expansion to obtain closed evolution equations for the stress tensor, derives the collisional rates and transport coefficients (as integrals over the distribution), and fits those derived expressions with simple analytical functions solely to capture their temperature dependence for practical use. The resulting model is then compared to independent DSMC simulations for shear stress, anisotropy, and viscosity, showing quantitative agreement. No step reduces by construction to its own inputs: the fitting is post-derivation convenience, not a redefinition of the target observables; DSMC provides external falsification; no self-citations, uniqueness theorems, or ansatzes are invoked as load-bearing. The chain remains self-contained against the external benchmark.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The framework rests on standard kinetic-theory assumptions plus a small number of fitting parameters introduced for practical use of the collisional integrals.

free parameters (1)
  • parameters in analytical fits for collisional rates and transport coefficients
    Simple analytical functions are fitted to capture the temperature dependence of rates over wide shear-rate ranges.
axioms (2)
  • domain assumption Boltzmann equation applies to dilute gas with the given interaction potential
    Standard starting point for kinetic theory of dilute systems.
  • domain assumption Grad's moment expansion truncates the velocity distribution sufficiently for stress-tensor evolution
    Approximation invoked to close the moment equations.

pith-pipeline@v0.9.0 · 5431 in / 1176 out tokens · 40267 ms · 2026-05-16T17:23:53.366530+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

38 extracted references · 38 canonical work pages

  1. [1]

    Case for α = 4 We first consider the case α = 4. Let us consider the following equation: 1 − x2 − (x ˜b )4 = 0, (A1) which admits two solutions, x2 1 ≡ ˜b2 ( − ˜b2 + √ ˜b4 + 4 ) 2 , x 2 2 ≡ ˜b2 ( ˜b2 + √ ˜b4 + 4 ) 2 , (A2) with x1 ≤ x2. Using these, Eq. (A1) can be factorized as 1 − x2 − (x ˜b )4 = 1 ˜b4 (x2 1 − x2)(x2 2 + x2). (A3) Accordingly, the collis...

    work page

  2. [2]

    In this case, the turning point is determined by the cubic equation in x2: 1 − x2 − (x ˜b )6 = 0

    Case for α = 6 We next consider the case α = 6. In this case, the turning point is determined by the cubic equation in x2: 1 − x2 − (x ˜b )6 = 0. (A12) 9 This equation can in principle be solved by the cubic formula or by applying Vieta’s formula. For convenience, we define F (X) ≡ 1 − X − X 3 ˜b6 , (A13) which is a strictly decreasing function. Since F (0...

    work page

  3. [3]

    Complex plas- mas and colloidal dispersions: particle-resolved studies of classical liquids and solids

    A. Ivlev , G. Morfill, and C. P. Royall, “Complex plas- mas and colloidal dispersions: particle-resolved studies of classical liquids and solids” (World Scientific, 2012)

    work page 2012

  4. [4]

    Dusty plasmas

    V. E. Fortov, A. G. Khrapak, S. A. Khrapak, V. I. Molotkov, and O. F. Petrov, “Dusty plasmas” Physics- Uspekhi 47, 447–492 (2004)

    work page 2004

  5. [5]

    Contact electrifica- tion of insulating materials

    D. J. Lacks and R. M. Sankaran, “Contact electrifica- tion of insulating materials ” J. Phys. D: Appl. Phys. 44, 453001 (2011)

    work page 2011

  6. [6]

    Experimental gen- eration of volcanic lightning

    C. Cimarelli, M. A. Alatorre-Ibarg¨ uengoitia, U. Kuep- pers, B. Scheu, and D. B. Dingwell, “Experimental gen- eration of volcanic lightning” Geol. 42, 79–82, (2014)

    work page 2014

  7. [7]

    Kinetic theory for plane flows of a dense gas of identical, rough, inelastic, circular disks

    J. T. Jenkins and M. W. Richman “Kinetic theory for plane flows of a dense gas of identical, rough, inelastic, circular disks ” Phys. Fluids. 28, 3485–3494, (1985)

    work page 1985

  8. [8]

    Hy- drodynamics for granular flow at low density

    J. J. Brey, J. W. Dufty, C. S. Kim, and A. Santos, “Hy- drodynamics for granular flow at low density” Phys. Rev. E 58, 4638 (1998)

    work page 1998

  9. [9]

    Dense fluid transport for in- elastic hard spheres

    V. Garz´ o and J. W. Dufty, “Dense fluid transport for in- elastic hard spheres” Phys. Rev. E 59, 5895–5911 (1999)

    work page 1999

  10. [10]

    Rapid granular flows

    I. Goldhirsch, “Rapid granular flows” Annu. Rev. Fluid Mech. 267–293,35 , (2003)

    work page 2003

  11. [11]

    Non-Newtonian stress, collisional dissipation and heat flux in the shear flow of inelastic disks: a reduction via Grad’s moment method

    S. Saha and M. Alam, “Non-Newtonian stress, collisional dissipation and heat flux in the shear flow of inelastic disks: a reduction via Grad’s moment method” J. Fluid Mech. 757, 251–296 (2014)

    work page 2014

  12. [12]

    Velocity correlations in dense granular shear flows: Effects on energy dissipation and normal stress

    N. Mitarai and H. Nakanishi, “Velocity correlations in dense granular shear flows: Effects on energy dissipation and normal stress” Phys. Rev. E 75, 031305 (2007)

    work page 2007

  13. [13]

    A modified kinetic theory for frictional granular flows in dense and dilute regimes

    S. Chialvo and S. Sundaresan, “A modified kinetic theory for frictional granular flows in dense and dilute regimes” Phys. Fluids 25, 070603 (2013)

    work page 2013

  14. [14]

    Collision rates in charged granular gases

    T. Scheffler and D. E. Wolf, “Collision rates in charged granular gases” Granul. Matter 4, 103 (2002)

    work page 2002

  15. [15]

    Long- time behavior of granular gases with impact-velocity de- pendent coefficient of restitution

    T. P¨ oschel, N. V. Brilliantov, and T. Schwager, “Long- time behavior of granular gases with impact-velocity de- pendent coefficient of restitution” Physica A 325, 274 (2003)

    work page 2003

  16. [16]

    Homogeneous cooling state of dilute granular gases of charged parti- cles

    S. Takada, D. Serero, and T. P¨ oschel, “Homogeneous cooling state of dilute granular gases of charged parti- cles” Phys. Fluids 29, 083303 (2017)

    work page 2017

  17. [17]

    Early-stage aggregation in three-dimensional charged granular gas

    C. Singh and M. G. Mazza, “Early-stage aggregation in three-dimensional charged granular gas” Phys. Rev. E 97, 022904 (2018)

    work page 2018

  18. [18]

    Electrification in granular gases leads to constrained fractal growth

    C. Singh and M. G. Mazza, “Electrification in granular gases leads to constrained fractal growth” Sci. Rep. 9, 9049 (2019)

    work page 2019

  19. [19]

    Transport coeffi- cients for granular gases of electrically charged particle s

    S. Takada, D. Serero, and T. P¨ oschel, “Transport coeffi- cients for granular gases of electrically charged particle s” J. Fluid Mech. 935, A38 (2022)

    work page 2022

  20. [20]

    Rheology of dilute granular gas mixtures where the grains interact via a square shoulder and well potential

    K. Yoshii, S. Takada, K. Kurosawa, and T. P¨ oschel, “Rheology of dilute granular gas mixtures where the grains interact via a square shoulder and well potential” Phys. Fluids 35, 013327 (2023)

    work page 2023

  21. [21]

    Kinetic theory of dilute granular gases hav- ing an inverse power law potential

    S. Takada, “Kinetic theory of dilute granular gases hav- ing an inverse power law potential” arXiv:2406.10864

    work page arXiv

  22. [22]

    Rheology of dilute granular gases with hard-core and inverse power- law potentials

    Y. Kobayashi, S. Iizuka, and S. Takada, “Rheology of dilute granular gases with hard-core and inverse power- law potentials” EPJ Web. Conf. 340, 03003 (2025)

    work page 2025

  23. [23]

    Classical Mechan- ics Third Edition

    H. Goldstein, C. Poole, and J. Safko, “Classical Mechan- ics Third Edition” (Addison Wesley, Bostion, 2001)

    work page 2001

  24. [24]

    The Mathematical Theory of Non Uniform Gases

    S. Chapman and T. G. Cowling, “The Mathematical Theory of Non Uniform Gases” (Cambridge Univ. Press, Cambridge, 1970)

    work page 1970

  25. [25]

    Kinetic Theory of Granular Gases

    N. V. Brilliantov and T. P¨ oschel, “Kinetic Theory of Granular Gases” (Oxford Univ. Press, Oxford, 2004)

    work page 2004

  26. [26]

    Granular Gaseous Flows —A Kinetic The- ory Approach to Granular Gaseous Flows—

    V. Garz´ o, “Granular Gaseous Flows —A Kinetic The- ory Approach to Granular Gaseous Flows—” (Springer, 2019)

    work page 2019

  27. [27]

    Inherent rheology of a granular fluid in uniform shear flow

    A. Santos, V. Garz´ o, J. W. Dufty, “Inherent rheology of a granular fluid in uniform shear flow” Phys. Rev. E 69, 061303 (2004)

    work page 2004

  28. [28]

    Kinetic theory of discon- tinuous rheological phase transition for a dilute iner- tial suspension

    H. Hayakawa and S. Takada, “Kinetic theory of discon- tinuous rheological phase transition for a dilute iner- tial suspension” Prog. Theor. Exp. Phys. 2019, 083J01 (2019)

    work page 2019

  29. [29]

    Kinetic the- ory of shear thickening for a moderately dense gas-solid suspension: From discontinuous thickening to continuous thickening

    H. Hayakawa, S. Takada, and V. Garz´ o, “Kinetic the- ory of shear thickening for a moderately dense gas-solid suspension: From discontinuous thickening to continuous thickening” Phys. Rev. E 96, 042903 (2017), [Erratum] 10 101, 069904 (2020)

    work page 2017

  30. [30]

    En- skog kinetic theory of rheology for a moderately dense inertial suspension

    S. Takada, H. Hayakawa, A. Santos, and V. Garz´ o, “En- skog kinetic theory of rheology for a moderately dense inertial suspension” Phys. Rev. E 102, 022907 (2020)

    work page 2020

  31. [31]

    On the kinetic theory of rarefied gases

    H. Grad, “On the kinetic theory of rarefied gases” Com- mun. Pure Appl. Math. 2, 331–407 (1949)

    work page 1949

  32. [32]

    Rheological properties in a low-density granular mixture

    M. Montanero and V. Garz´ o, “Rheological properties in a low-density granular mixture” Physica A 310, 17–38 (2002)

    work page 2002

  33. [33]

    Energy nonequipartition in a sheared granular mixture

    M. Montanero and V. Garz´ o, “Energy nonequipartition in a sheared granular mixture” Mol. Simul. 29, 357–362 (2003)

    work page 2003

  34. [34]

    Molecular Gas Dynamics and the Direct Simulation of Gas Flows

    G. A. Bird, “Molecular Gas Dynamics and the Direct Simulation of Gas Flows” (Oxford Univ. Press, New York, 1994)

    work page 1994

  35. [35]

    Computational Granular Dynamics

    T. P¨ oschel and T. Schwager, “Computational Granular Dynamics” (Springer, Berlin, 2005)

    work page 2005

  36. [36]

    DSMC eval- uation of the Navier-Stokes shear viscosity of a granular fluid

    J. M. Montanero, A. Santos, and V. Garz´ o, “DSMC eval- uation of the Navier-Stokes shear viscosity of a granular fluid” AIP Conf. Proc. 762, 797–802 (2005)

    work page 2005

  37. [37]

    Rheology of a dilute binary mixture of inertial suspension under simple shear flow

    S. Takada, H. Hayakawa, V. Garz´ o, “Rheology of a dilute binary mixture of inertial suspension under simple shear flow” Prog. Theor. Exp. Phys. 2023, 113J01 (2023)

    work page 2023

  38. [38]

    Handbook of Mathe- matical Functions: With Formulas, Graphs, and Mathe- matical Tables

    M. Abramowitz and I. A. Stegun, “Handbook of Mathe- matical Functions: With Formulas, Graphs, and Mathe- matical Tables” (Dover, New York, 1964)

    work page 1964