Longtime behavior for homoenergetic solutions in the collision dominated regime for hard potentials

This paper establishes the long-time convergence of collision-dominated homoenergetic solutions for hard potentials with initially high temperatures to a Maxwellian distribution with diverging temperature, deriving precise asymptotic formulas via a Hilbert-type expansion ansatz for both cutoff and non-cutoff kernels.

The Gist

Imagine a vast, invisible dance floor filled with billions of tiny, energetic dancers. These dancers are gas molecules. In a calm room, they bump into each other randomly, eventually settling into a predictable, average pattern of movement known as a Maxwellian distribution (think of it as the "chill mode" where everyone is moving at a comfortable, average speed).

But what happens if the dance floor itself starts to move? What if the floor stretches, shrinks, or gets twisted by a giant invisible hand? This is the world of homoenergetic solutions studied in this paper.

Here is the story of what the author, Bernhard Kepka, discovered, explained simply:

The Setup: The Stretching Dance Floor

The paper looks at a specific type of gas behavior where the gas isn't just sitting still; it's being acted upon by a "drift." Imagine the gas is inside a box that is being:

1. Sheared: Like pulling the top of a deck of cards to the side while holding the bottom still.
2. Dilated: Like stretching a rubber sheet in all directions.
3. A mix of both.

Usually, when you stretch or shear a gas, you might expect it to cool down (because the molecules spread out) or behave chaotically. However, this paper focuses on "Hard Potentials." Think of these as gas molecules that are like super-bouncy, hard rubber balls. When they hit each other, they don't just gently nudge; they collide with significant force.

The Big Surprise: The Gas Gets Hotter and Hotter

The main discovery is counter-intuitive. If you start with a gas that is already very hot (the dancers are moving fast), and you start shearing or stretching the container:

• The gas does not cool down.
• Instead, the temperature goes to infinity.

The "dance" gets more and more frantic. The faster the gas moves, the more violent the collisions become, which pumps even more energy into the system. It's a runaway effect where the gas heats up faster and faster as time goes on.

The Secret Weapon: The "Hilbert Expansion"

How did the author prove this? He used a mathematical trick called a Hilbert-type expansion.

Imagine trying to describe the chaotic motion of a million dancers. It's impossible to track everyone individually. So, the author uses a "layered" approach:

1. Layer 1 (The Base): He assumes the gas is mostly in that calm, average "Maxwellian" state (the chill mode).
2. Layer 2 (The Wobble): He adds a tiny "wobble" or disturbance to account for the stretching and shearing forces.
3. Layer 3 (The Correction): He adds even tinier corrections to fix the math.

By stacking these layers, he could show that the "wobble" caused by the shear is so strong that it overpowers the cooling effect of the stretching. The math proves that the "wobble" keeps feeding energy into the system, driving the temperature up forever.

The Three Scenarios

The author looked at three specific ways the "dance floor" could move, and found that in all three, the gas gets hotter, but at different speeds:

1. Simple Shear: The floor is pulled sideways at a constant speed. The temperature grows steadily, like a car accelerating at a constant rate.
2. Shear with Decaying Stretch: The floor is pulled sideways, but the stretching effect fades away over time. The temperature still grows, but the math gets a bit more complex (like a car accelerating but with a slightly slipping tire).
3. Combined Orthogonal Shear: The floor is twisted in two different directions at once. This is the most chaotic scenario. Here, the temperature explodes even faster, growing like the cube of time (if you wait twice as long, it gets 8 times hotter, not just 2 times).

Why Does This Matter?

You might ask, "Who cares about gas in a stretching box?"

This isn't just about boxes. This math helps us understand:

• Aerospace: How gases behave in the thin atmosphere of space where planes fly at high speeds.
• Astrophysics: How gas clouds behave when they are being stretched by the gravity of stars or black holes.
• Engineering: How to model fluids in extreme conditions where standard "cooling" assumptions fail.

The Takeaway

The paper is a rigorous proof that for certain types of bouncy, hard gas molecules, if you twist or stretch them hard enough, they won't settle down. Instead, they will enter a state of perpetual heating, getting hotter and hotter forever, driven by the very collisions that usually keep things stable.

The author didn't just guess this; he built a mathematical bridge (using the "Hilbert expansion") to prove that the gas stays close to a predictable pattern even as it heats up to infinity, giving us precise formulas for exactly how fast it gets hot. It's a victory for understanding the chaotic dance of the universe.

Technical Summary

1. Problem Statement

The paper investigates the long-time asymptotic behavior of a specific class of solutions to the inhomogeneous Boltzmann equation, known as homoenergetic solutions. These solutions describe the dynamics of a dilute gas under the combined action of particle collisions and external flow fields (specifically shear, dilation, or a combination).

The governing equation for the distribution function $f(t,x,v)$ is:
$$\partial_t f + v \cdot \nabla_x f = Q(f,f)$$
Homoenergetic solutions take the form $f(t,x,v) = g(t, v - L(t)x)$, where $L(t)$ is a time-dependent matrix satisfying $\dot{L} + L^2 = 0$. This reduces the problem to a kinetic equation for $g(t,w)$ involving a drift term and the collision operator $Q$.

The Core Challenge:
The paper focuses on the collision-dominated regime for hard potentials ($\gamma > 0$). In this regime, the collision operator dominates the drift term as $t \to \infty$. Previous work (specifically [22]) had formulated conjectures regarding the long-time behavior of these solutions using formal Hilbert-type expansions, suggesting that solutions converge to a Maxwellian distribution with a temperature that grows to infinity. However, these previous results lacked rigorous mathematical proofs. The goal of this paper is to provide a rigorous justification of these conjectures, including precise asymptotic formulas for the temperature and convergence rates.

2. Methodology

The author employs a multi-step analytical strategy combining formal asymptotic analysis with rigorous functional analysis in weighted Sobolev spaces.

A. Rescaling and Ansatz

The analysis begins by rescaling the distribution function to normalize mass, momentum, and temperature. Let $f_t(v)$ be the rescaled distribution. The author introduces a Hilbert-type expansion ansatz:
$$f_t(v) = \mu(v) + \bar{\mu}_t(v) + h_t(v)$$
where:

• $\mu(v)$ is the standard Maxwellian equilibrium.
• $\bar{\mu}_t$ is the first-order correction term derived from the linearized collision operator $L$.
• $h_t$ is the remainder (error term).

The inverse temperature $\beta_t = T_t^{-1}$ is shown to satisfy a differential equation driven by the trace of the shear matrix and the interaction between the drift and the distribution.

B. Functional Framework

The proofs rely on estimates in weighted Sobolev spaces $H^1_{p}$ and $L^1_m$.

• Non-cutoff kernels: The paper primarily treats non-cutoff kernels (grazing collisions), where the angular part of the collision kernel has a singularity. This requires the use of anisotropic norms (specifically the norm $\|\cdot\|_{H^{s,*}_p}$ introduced in [21]) to handle the fractional regularity induced by the singularity.
• Linearized Operator: The linearized collision operator $L$ is utilized. For hard potentials ($\gamma > 0$), $L$ possesses a spectral gap in appropriate weighted spaces, allowing for exponential decay estimates.

C. Key Technical Tools

1. Povzner Estimates: Used to establish moment bounds and uniqueness of weak solutions.
2. Coercivity Estimates: Leveraging the spectral gap of $L$ to control the dissipation of the error term $h_t$.
3. Continuation Argument: A bootstrap argument is used. The author assumes a priori bounds on the error term and the temperature growth, proves that these bounds are self-consistent for small initial perturbations, and then extends the result to all time $t \geq 0$.
4. Time Rescaling: For cases where the shear matrix $L(t)$ grows unbounded (e.g., combined orthogonal shear), a change of time variable is introduced to map the problem to a setting with bounded coefficients, allowing the application of the main stability theorem.

3. Key Contributions

• Rigorous Proof of Conjectures: The paper provides the first rigorous proof of the long-time asymptotics for homoenergetic solutions with hard potentials, confirming the formal computations found in [22].
• Precise Asymptotic Formulas: It derives exact asymptotic rates for the temperature $T_t$ (or inverse temperature $\beta_t$) for three distinct types of shear flows:
  1. Simple Shear: $T_t \sim t^{2/\gamma}$.
  2. Simple Shear with Decaying Dilatation: $T_t \sim t^{2/\gamma}$ (with a specific logarithmic correction factor depending on the trace).
  3. Combined Orthogonal Shear: $T_t \sim t^{6/\gamma}$ (due to the linear growth of the shear matrix).
• Stability in Weighted Spaces: The results establish the local stability of the Maxwellian distribution in weighted Sobolev spaces $H^1_{p_0}$, proving that solutions starting close to equilibrium remain close and converge.
• Extension to Cutoff Kernels: The methodology is adapted to cover both non-cutoff (grazing collisions) and cutoff kernels, demonstrating the robustness of the approach.

4. Main Results

The paper establishes the following theorems under the assumption that the initial data is sufficiently close to a Maxwellian and the initial temperature is sufficiently high (ensuring the collision term remains dominant).

Theorem 1.1 (Convergence and Temperature Asymptotics):
For matrices $L(t)$ corresponding to simple shear, simple shear with dilatation, or combined orthogonal shear:

• The rescaled solution $f_t$ converges to the Maxwellian $\mu$ in the $H^1_{p_0}$ norm as $t \to \infty$.
• The inverse temperature $\beta_t$ decays according to specific power laws:
  • Simple Shear: $\lim_{t\to\infty} \beta_t^{-\gamma/2} / t = \text{const}$.
  • Shear with Dilatation: $\lim_{t\to\infty} \beta_t^{-\gamma/2} / t^2 = \text{const}$.
  • Combined Shear: $\lim_{t\to\infty} \beta_t^{-\gamma/2} / t^3 = \text{const}$.

Theorem 1.2 (Quantitative Decay Rates):
The paper provides explicit decay rates for the perturbation $h_t = f_t - \mu - \bar{\mu}_t$:

• For simple shear, $\|h_t\| \sim O((1+t)^{-2})$.
• For combined orthogonal shear, the decay is faster, $\|h_t\| \sim O((1+t)^{-4})$, due to the stronger driving mechanism of the shear.
• The error term is bounded by a combination of the initial perturbation size and the inverse of the temperature growth parameter $\zeta_t$.

5. Significance

• Mathematical Rigor: This work bridges the gap between formal asymptotic expansions and rigorous analysis for non-equilibrium kinetic equations. It validates the use of Hilbert expansions in the context of time-dependent, non-stationary flows.
• Physical Insight: The results clarify the competition between dilatation (which cools the gas) and shear (which heats the gas via viscous dissipation). The paper proves that for hard potentials, the shear effect dominates, leading to an unbounded increase in temperature (energy) over time, preventing the system from reaching a static equilibrium.
• Methodological Advancement: The successful application of anisotropic norms and weighted Sobolev estimates to the non-cutoff Boltzmann equation in a time-dependent setting provides a template for analyzing other non-equilibrium kinetic problems where the collision operator dominates the transport.
• Generality: By handling both cutoff and non-cutoff kernels, the results apply to a wide range of physical interactions, including inverse power-law potentials and hard spheres.

In summary, Kepka's paper resolves a significant open problem in kinetic theory by proving that homoenergetic solutions for hard potentials converge to a time-dependent Maxwellian with a temperature growing polynomially in time, driven by the dominance of the shear flow over the expansion.