pith. sign in

arxiv: 1907.03416 · v1 · pith:CV6WZDF7new · submitted 2019-07-08 · ❄️ cond-mat.stat-mech

Shock propagation in the hard sphere gas in two dimensions: comparison between simulations and hydrodynamics

Jilmy P. Joy , R. Rajesh This is my paper

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

classification ❄️ cond-mat.stat-mech
keywords hard sphere gasshock propagationhydrodynamicsmolecular dynamicsblast wavetwo dimensionsnon-equilibrium
0
0 comments X

The pith

Large-scale simulations in two dimensions show that hydrodynamics does not describe shock propagation in a hard sphere gas well.

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

The paper uses event-driven molecular dynamics simulations to follow a two-dimensional hard sphere gas initially at rest after localized kinetic energy injection. Radial profiles of density, pressure, temperature and flow velocity become scale invariant at long times. These profiles are compared against the hydrodynamic equations obtained from the continuity relations for mass, momentum and energy, the same equations that describe a classical blast wave. The simulations show a clear mismatch with the hydrodynamic predictions. The paper verifies that local equilibrium holds and an equation of state exists yet the discrepancy remains.

Core claim

Contrary to earlier smaller-scale reports of agreement in two dimensions, large-scale simulations demonstrate that the radial distributions from the hard sphere gas do not match the solutions of the hydrodynamic continuity equations for mass, momentum and energy, just as observed in three dimensions. The mismatch persists after the authors check the assumptions of local equilibrium, the existence of an equation of state, and the neglect of heat conduction and viscosity.

What carries the argument

Hydrodynamic description from the continuity equations for mass, momentum and energy, tested against event-driven molecular dynamics simulations of hard spheres.

If this is right

  • Hydrodynamics fails to match simulation data in both two and three dimensions.
  • Verification of local equilibrium and the equation of state does not remove the mismatch.
  • Neglect of heat conduction and viscosity in the hydrodynamic treatment contributes to the observed discrepancy.
  • Scale-invariant growth appears in the simulations but is not reproduced by the hydrodynamic solutions.

Where Pith is reading between the lines

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

  • Earlier reports of agreement in two dimensions likely resulted from insufficient system sizes that masked the mismatch.
  • The result raises the question whether hydrodynamic descriptions require explicit dissipative terms for any dimension when modeling strong shocks in dilute gases.
  • Similar large-scale simulations with different interaction potentials could test whether the mismatch is specific to hard spheres.

Load-bearing premise

The hydrodynamic continuity equations are assumed to capture the dynamics once local equilibrium is reached and an equation of state is available.

What would settle it

Quantitative comparison of the simulated radial density or velocity profiles at late times against the corresponding hydrodynamic solutions would show whether they agree within statistical fluctuations.

Figures

Figures reproduced from arXiv: 1907.03416 by Jilmy P. Joy, R. Rajesh.

Figure 1
Figure 1. Figure 1: Moving (red) and stationary (blue) particles at times (a) t = 1 × 105 , (b) t = 1.5 × 105 , (c) t = 2.0 × 105 and (d) t = 2.5 × 105 , after the initial injection of four energetic particles at the center. At time t = 2.5 × 105 , there are 919407 moving particles. The data are for the ambient number density ρ0 = 0.15. of the particles i and j undergoing collisions, δt is the time duration of measurement, an… view at source ↗
Figure 2
Figure 2. Figure 2: Simulation results for the temporal variation of (a) the radius of the shock R(t) and (b) the number of moving particles N(t). The solid lines are power laws (a) √ t and (b) t. The data are for the ambient number density ρ0 = 0.15. energy can be replaced by that for entropy. Since the flow is isotropic, the different thermodynamic quantities cannot depend on the angle. Thus, in radial coordinates, the cont… view at source ↗
Figure 3
Figure 3. Figure 3: (Color online) The scaling functions (a) R(ξ), (b) V (ξ), (c) E(ξ), and (d) P(ξ) corresponding to density, velocity, temperature and pressure respectively versus ξ obtained from hydrodynamic equations for ambient number density (a)–(d) ρ0 = 0.15 and (e)–(h) ρ0 = 0.382. n refers to the number of terms that is retained in the virial expansion (n = 1 is ideal gas). The insets show the plots on a log-log scale… view at source ↗
Figure 4
Figure 4. Figure 4: (Color online) The variation of the scaling functions (a) R(ξ), (b) V (ξ), (c) E(ξ) and (d) P(ξ) corresponding to non-dimensionalised density, velocity, temperature and pressure [see Eq. (11)] with scaled distance ξ. The data are shown for 2 different initial densities ρ0 = 0.15 and 0.382. For ρ0 = 0.15, the different times are t = 100000, 150000, 200000, 250000, and for ρ = 0.382, t = 50000, 100000, 13000… view at source ↗
Figure 5
Figure 5. Figure 5: (Color online) The data in [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: (Color online) The variation of χ(n) [see Eq. (19)] with ξ for n = 2, 4, 6, 8, 10. The data are for times 150000 and 250000 and for ambient number density ρ0 = 0.15. For large n, χ(n) converges to one. 5. Verifying the assumptions of the TvNS theory We now numerically check the different assumptions of the TvNS theory. 5.1. Equation of state One of the key assumptions of the TvNS theory is an EOS relates t… view at source ↗
Figure 7
Figure 7. Figure 7: (Color online) The variation of ζ, the ratio of thermal energies in the radial and transverse directions [see Eq. (20)] with the scaled distance ξ. The data is for four different times with keys as in [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: (Color online) The variation with scaled distance ξ of (a) the kurtosis κr for the radial velocity fluctuations. (b) the kurtosis κ⊥ for the velocity fluctuations in the θ direction and (c) skewness S for the radial velocity fluctuations. The data are for ρ0 = 0.15 and ρ0 = 0.382 and for four different times with keys as in [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: (Color online) The variation of Eflow [see Eq. (23)] with time. The data are for ambient number density ρ0 = 0.15. compared with the fit to a gaussian. Clearly, the distribution deviates from a gaussian, is asymmetric, and is skewed towards the larger positive fluctuations. 5.4. Energy of mean flow The total energy of the system can be divided into two parts: one is from the mean flow velocity and the othe… view at source ↗
read the original abstract

We study the radial distribution of pressure, density, temperature and flow velocity fields at different times in a two dimensional hard sphere gas that is initially at rest and disturbed by injecting kinetic energy in a localized region through large scale event driven molecular dynamics simulations. For large times, the growth of these distributions are scale invariant. The hydrodynamic description of the problem, obtained from the continuity equations for the three conserved quantities -- mass, momentum, and energy -- is identical to those used to describe the hydrodynamic regime of a blast wave propagating through a medium at rest, following an intense explosion, a classic problem in gas dynamics. Earlier work showed that the results from simulations matched well with the predictions from hydrodynamics in two dimensions, but did not match well in three dimensions. To resolve this contradiction, we perform large scale simulations in two dimensions, and show that like in three dimensions, hydrodynamics does not describe the simulation data well. To account for this discrepancy, we check in our simulations the different assumptions of the hydrodynamic approach like local equilibrium, existence of an equation of state, neglect of heat conduction and viscosity.

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

1 major / 0 minor

Summary. The paper performs large-scale event-driven molecular dynamics simulations of a 2D hard-disk gas with localized kinetic energy injection. It reports that the resulting scale-invariant radial profiles of density, flow velocity, pressure and temperature fail to match the predictions of the standard hydrodynamic blast-wave solution obtained from the three continuity equations, in contrast to earlier 2D claims but consistent with 3D results. The authors explicitly check the hydrodynamic assumptions of local equilibrium, existence of an equation of state, and negligibility of heat conduction and viscosity.

Significance. If the reported mismatch is quantitatively robust, the work demonstrates that conventional hydrodynamics can fail to describe blast-wave propagation in two-dimensional gases even when the usual assumptions appear satisfied, with potential implications for the validity of hydrodynamic descriptions in low-dimensional non-equilibrium systems.

major comments (1)
  1. [Abstract] Abstract: the central claim that 'hydrodynamics does not describe the simulation data well' is not supported by any quantitative measure of mismatch (e.g., integrated squared deviation, point-wise relative errors, or statistical significance) or by explicit data-exclusion criteria; without these the strength of the discrepancy cannot be assessed.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting the need for quantitative support of our central claim. We address the single major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim that 'hydrodynamics does not describe the simulation data well' is not supported by any quantitative measure of mismatch (e.g., integrated squared deviation, point-wise relative errors, or statistical significance) or by explicit data-exclusion criteria; without these the strength of the discrepancy cannot be assessed.

    Authors: We agree that the strength of the reported discrepancy would be clearer with explicit quantitative measures. In the revised manuscript we will add (i) the integrated squared deviation between the scaled simulation profiles and the hydrodynamic solution over the self-similar region, (ii) point-wise relative errors at selected radii and times, and (iii) an explicit statement of the radial interval used for each comparison (thereby removing any ambiguity about data-exclusion criteria). These additions will be placed both in the abstract and in a new subsection of the results. revision: yes

Circularity Check

0 steps flagged

No significant circularity; independent MD vs. standard hydrodynamics

full rationale

The central comparison uses event-driven molecular dynamics trajectories that are generated independently of the hydrodynamic model. The hydrodynamic blast-wave solution is obtained directly from the three continuity equations for mass, momentum and energy with no parameters fitted from the present simulations. The paper explicitly verifies local equilibrium and the equation of state in the data but does not insert any simulation-derived quantities back into the hydrodynamic prediction; the reported mismatch is therefore an external test rather than a self-referential construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard continuity equations for mass, momentum, and energy together with the validity of event-driven molecular dynamics; no additional free parameters or invented entities appear in the abstract.

axioms (1)
  • standard math Continuity equations for mass, momentum, and energy conservation govern the hydrodynamic description
    Explicitly invoked in the abstract to obtain the blast-wave equations.

pith-pipeline@v0.9.0 · 5727 in / 1212 out tokens · 27340 ms · 2026-05-25T01:18:50.461348+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

44 extracted references · 44 canonical work pages · 1 internal anchor

  1. [1]

    Whitham G 1974 Linear and Nonlinear Waves (New York: Wiley)

    work page 1974

  2. [2]

    Barenblatt G 1987 Scaling, Self-similarity, and Intermediate Asymptotics: Dimensional Analysis and Intermediate Asymptotics (Cambridge: Cambridge University Press)

    work page 1987

  3. [3]

    Taylor G 1950 Proc. R. Soc. Lond. A 201 159

    work page 1950

  4. [4]

    Taylor G 1950 Proc. R. Soc. Lond. A 201 175–186 Shock propagation in hard sphere gas 17

    work page 1950

  5. [5]

    von Neumann J 1963 Collected Works (Oxford: Pergamon Press) p 219

    work page 1963

  6. [6]

    Sedov L 1993 Similarity and Dimensional Methods in Mechanics 10th ed (Florida: CRC Press)

    work page 1993

  7. [7]

    Sedov L 1946 J. Appl. Math. Mech. 10 241

    work page 1946

  8. [8]

    Woltjer L 1972 Ann. Rev. Astron. Astrophys. 10 129–158

    work page 1972

  9. [9]

    Gull S 1973 Mon. Not. R. Astr. Soc. 161 47–69

    work page 1973

  10. [10]

    Cioffi D F, Mckee C F and Bertschinger E 1988 The Astrophysical Journal 334 252–265

    work page 1988

  11. [11]

    Ostriker J P and McKee C F 1988 Rev. Mod. Phys. 60(1) 1–68

    work page 1988

  12. [12]

    Zel’dovich Y B and Raizer Y P 2002 Physics of Shock Waves and High Temperature Hydrodynamic Phenomena (New York: Dover Publications, Inc.)

    work page 2002

  13. [13]

    Edwards M J, MacKinnon A J, Zweiback J, Shigemori K, Ryutov D, Rubenchik A M, Keilty K A, Liang E, Remington B A and Ditmire T 2001 Phys. Rev. Lett. 87(8) 085004

    work page 2001

  14. [14]

    Plasmas 11(11) 4968–4972

    Edens A, Ditmire T, Hansen J, Edwards M, Adams R, Rambo P, Ruggles L, Smith I and Porter J 2004 Phys. Plasmas 11(11) 4968–4972

    work page 2004

  15. [15]

    Plasmas 12(05) 052707–1–052707–7

    Moore A S, Symes D R and Smith R A 2005 Phys. Plasmas 12(05) 052707–1–052707–7

    work page 2005

  16. [16]

    Dokuchaev V I 2002 Astronomy and Astrophysics 395(3) 1023–1029

    work page 2002

  17. [17]

    and Astrophys

    Falle S 1975 Aston. and Astrophys. 43 323–336

    work page 1975

  18. [18]

    Fluid Mech 117 473–491

    Ghoniem A, Kamel M, Berger S and Oppenheim A 1982 J. Fluid Mech 117 473–491

    work page 1982

  19. [19]

    Abdel-Raouf A and Gretler W 1991 Fluid Dyn. Res. 8 273–285

    work page 1991

  20. [20]

    Fluids 6 2154

    Steiner H and Gretler W 1994 Phys. Fluids 6 2154

    work page 1994

  21. [21]

    VonNeumann J and Richtmyer R 1950 Journal of Applied Physics 21 232

    work page 1950

  22. [22]

    Latter R 1955 Journal of Applied Physics 26 954

    work page 1955

  23. [23]

    Brode H L 1955 Journal of Applied Physics 26 766

    work page 1955

  24. [24]

    Fluids 13 2665

    Plooster M N 1970 The Phys. Fluids 13 2665

    work page 1970

  25. [25]

    Walsh A M, Holloway K E, Habdas P and de Bruyn J R 2003 Phys. Rev. Lett. 91(10) 104301

    work page 2003

  26. [26]

    Proc 1145 767

    Metzger P T, Latta R C, Schuler J M and Immer C D 2009 AIP Conf. Proc 1145 767

    work page 2009

  27. [27]

    Grasselli Y and Herrmann H J 2001 Gran Matt 3 201–204

    work page 2001

  28. [28]

    Boudet J F, Cassagne J and Kellay H 2009 Phys. Rev. Lett. 103(22) 224501

    work page 2009

  29. [29]

    Jabeen Z, Rajesh R and Ray P 2010 Eur. Phys. Lett. 89 34001

    work page 2010

  30. [30]

    Pathak S N, Jabeen Z, Ray P and Rajesh R 2012 Phys. Rev. E 85(6) 061301

    work page 2012

  31. [31]

    Cheng X, Xu L, Patterson A, Jaeger H M and Nagel S R 2008 Nat Phys 4 234

    work page 2008

  32. [32]

    Sandnes B, Knudsen H A, M˚ aløy K J and Flekkøy E G 2007 Phys. Rev. Lett. 99(3) 038001 URL http://link.aps.org/doi/10.1103/PhysRevLett.99.038001

    work page doi:10.1103/physrevlett.99.038001 2007

  33. [33]

    Pinto S F, Couto M S, Atman A P F, Alves S G, Bernardes A T, de Resende H F V and Souza E C 2007 Phys. Rev. Lett. 99(6) 068001

    work page 2007

  34. [34]

    Johnsen O, Toussaint R, M˚ aløy K J and Flekkøy E G 2006Phys. Rev. E 74(1) 011301

    work page

  35. [35]

    Huang H, Zhang F and Callahan P 2012 Phys. Rev. Lett. 108(25) 258001

    work page 2012

  36. [36]

    Joy J P, Pathak S N, Dibyendu D and Rajesh R 2017 Phys. Rev. E 96(3) 032908

    work page 2017

  37. [37]

    Barbier M, Villamaina D and Trizac E 2015 Phys. Rev. Lett. 115(21) 214301

    work page 2015

  38. [38]

    Fluids 28 083302

    Barbier M, Villamaina D and Trizac E 2016 Phys. Fluids 28 083302

    work page 2016

  39. [39]

    Antal T, Krapivsky P L and Redner S 2008 Phys. Rev. E 78(3) 030301

    work page 2008

  40. [40]

    Joy J P, Pathak S N and Rajesh R 2018 arXiv:1812.03638 [cond-mat.stat-mech]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  41. [41]

    Rapaport D C 2004 The art of molecular dynamics simulations (Cambridge: Cambridge University Press)

    work page 2004

  42. [42]

    Isobe M 2016 Molecular Simulation 42(16) 1317–1329

    work page 2016

  43. [43]

    Landau L and Lifshitz E 1987 Course of Theoretical Physics- Fluid Mechanics (Oxford: Butterwrth-Heinemann)

    work page 1987

  44. [44]

    McCoy B M 2009 Advanced Statistical Mechanics (Oxford: Oxford Science Publications)

    work page 2009