Staggering domino-like blast front motion in a one-dimensional cold gas

TL;DR

This paper investigates a one-dimensional particle system modeling blast wave dynamics, revealing a special set of mass ratios where the system exhibits a domino-like, ballistic shock front motion with explicit solutions.

Contribution

It provides analytical and numerical analysis of the model for various mass ratios, discovering a family of ratios causing a drastic change to domino-like, ballistic shock front behavior with exact formulas.

Findings

For m/μ=2, results align with hydrodynamic predictions.

Identifies a family of mass ratios where only triplets move, leading to ballistic shock fronts.

Derives explicit formulas for particle velocities and positions at these ratios.

Abstract

One-dimensional alternating particle systems are widely used to study interconnections between the hydrodynamics of blast waves in a gas-like medium and the Newtonian dynamics of its corpuscular constituents. We study the model in which point particles with masses $m,\mu, m,\mu,\dots, (m\geq\mu)$ are distributed on the positive half-line $\mathbb{R}_{+}$. Their dynamics are initiated by giving a positive velocity to the leftmost particle; in its course, the particles undergo elastic collisions. For this model with $m/\mu=2$, it has previously been established that the dynamics that start from random initial positions are consistent with predictions based on Euler's hydrodynamic equation. In particular, they have the following properties: (i) the position of the rightmost particle (shock front) evolves as $t^\delta$ with $\delta<1$; (ii) recoiled particles behind the front enter the negative half-axis; (iii) particles with locations $x\leq0$ move ballistically and eventually take over the total energy of the system. In this paper, we present numerical and analytical results for the dynamics of this model with nonrandom (typically equidistant) initial positions and various values of $m/\mu$. For $m/\mu=2$ and equidistant initial positions, our results qualitatively agree with those just mentioned. At the same time, we found an infinite family of numbers $\{\mathcal{M}_k\}$ such that, for $m/\mu=\mathcal{M}_k$, the hydrodynamic behavior mentioned changes drastically to the following. At each moment, only a single triplet $m,\mu, m$ is in motion, whereas all other particles are at rest. As a result, the shock front moves ballistically with an average velocity equal to the initial one. Such a `staggering domino-like' picture is obtained as an exact solution, which yields, in particular, explicit formulas for $\mathcal{M}_k$ and the particle velocities and positions.