Numerical study of loss of hyperbolicity using a cold plasma model

This paper proposes a new implicit numerical method in Euler variables to solve one-dimensional cold plasma equations with density-dependent collision coefficients, effectively overcoming computational challenges associated with the loss of hyperbolicity while confirming theoretical predictions regarding solution smoothness.

The Gist

The Big Picture: A Crowd of Runners

Imagine a stadium filled with runners (electrons) who are all supposed to stay in their lanes. In a "cold plasma," these runners are packed so tightly that they move together like a single fluid.

Usually, these runners oscillate (run back and forth) in a smooth, rhythmic wave. However, if they run too fast or start too close together, the wave can "break." In physics terms, this is called a singularity or breaking effect. It's like a traffic jam where the cars suddenly pile up so high that the density becomes infinite. At this point, the mathematical rules that describe their motion stop working (the system "loses hyperbolicity"), and standard computer simulations crash or give nonsense results.

The Problem: Friction That Changes the Rules

Scientists have known for a long time that if you add "friction" (collisions between electrons and ions) to this system, it can smooth things out.

• Constant Friction: Imagine if every runner had the same amount of friction, no matter how crowded the track was. This helps, but it doesn't always stop the traffic jam from forming if the runners start too aggressively.
• Variable Friction (The New Idea): The paper looks at a more realistic scenario where friction depends on how crowded the track is. If the runners bunch up (high density), the friction gets stronger. It's like a crowd that gets harder and harder to push through the more people are in it.

The Catch: While this "crowd-dependent friction" is physically realistic, it breaks the math. It changes the type of equations from a stable "hyperbolic" system (like a predictable wave) to a tricky "non-hyperbolic" system (like a Jordan block). Standard computer tools designed for waves fail here because the math becomes unstable and prone to exploding with errors.

The Solution: A New Way to Calculate

The authors, Chizhonkov and Rozanova, built a new computer algorithm (a set of instructions for a computer) to handle this tricky math.

• The Old Way: Think of the old method as taking a snapshot of the runners, guessing where they will be next, and then correcting the guess. This works great for smooth waves but fails when the friction changes based on density.
• The New Way: They created an implicit method. Imagine instead of just guessing the future, the computer solves a puzzle where it figures out the future state and the current state simultaneously. It's like solving a maze by looking at the exit and the entrance at the same time. This approach is much more stable and prevents the computer from crashing, even when the math gets weird.

What They Found: The Results

They tested this new method on two scenarios: slow runners (non-relativistic) and super-fast runners (relativistic).

1. Smoothing the Waves: When they used the "crowd-dependent friction" (where friction increases with density), the waves didn't break as easily. The friction acted like a shock absorber that got stronger exactly when the runners started to bunch up.
2. Stopping the Break: In many cases, this variable friction completely stopped the "traffic jam" (the singularity) from forming, even when the runners started with enough energy to cause a crash in a frictionless world.
3. The Threshold: They found a "tipping point." If the friction is strong enough (specifically, if it grows faster than linearly with density), the waves stay smooth forever. If the friction is just a constant number, the waves might still break.
4. Relativity: Even when the runners were moving near the speed of light, the new method worked perfectly. It showed that while collisions delay the crash, they don't always stop it unless the friction is strong enough.

The Takeaway

The paper doesn't just say "collisions are good." It says: "If you model collisions correctly (where friction grows with density), you can prevent the mathematical breakdown of the system."

However, the authors also warn that this "fix" isn't magic. In some extreme cases, the waves can still break, but the new computer method allows scientists to see exactly when and how that happens without the simulation crashing. They successfully proved that their new "implicit" calculator is the right tool for the job, matching all known theoretical predictions.

In short: They built a better calculator for a specific type of physics problem that usually breaks computers, and they used it to show that "crowd-dependent friction" is a powerful way to keep plasma waves from crashing.

Technical Summary

Technical Summary: Numerical Study of Loss of Hyperbolicity Using a Cold Plasma Model

Problem Statement
The paper addresses the numerical simulation of one-dimensional cold plasma oscillations, specifically focusing on the system of hydrodynamic and Maxwell's equations that include electron-ion collisions. A central mathematical challenge arises when the collision frequency $\nu$ is modeled as a linear function of the electron density $N$ (i.e., $\nu = \nu_0 N + \nu_1$). While a constant collision coefficient preserves the hyperbolic nature of the system (albeit non-strictly), introducing a density-dependent coefficient causes the system to lose hyperbolicity. The Jacobian matrix of the system becomes a Jordan block with coinciding eigenvalues but only one eigenvector. This structural change renders traditional explicit difference schemes, designed for hyperbolic conservation laws, unstable over long time intervals due to the growth of computational errors. Furthermore, while constant collisions are known to damp oscillations, the specific impact of a linear density-dependent collision term on the "breaking" effect (the formation of singularities in electron density) requires rigorous numerical investigation, particularly in the threshold case where smoothness might still be lost.

Methodology
To overcome the instability associated with the loss of hyperbolicity, the authors propose a new implicit second-order solution method in Euler variables. This method is a generalization of the well-known explicit MacCormack scheme, adapted to handle the specific structure of the cold plasma equations with density-dependent collisions.

• Scheme Construction: The algorithm employs a predictor-corrector approach.
  • Predictor Step: Calculates an intermediate state using forward difference operators and the matrix $A(V, \gamma)$, which represents the system's derivatives.
  • Corrector Step: Refines the solution using backward difference operators.
  • Implicit Treatment: Crucially, the scheme includes an implicit term controlled by a parameter $\lambda$. When $\lambda = 0$, the scheme reduces to the explicit MacCormack scheme. For the non-hyperbolic case (where $\nu_0 \neq 0$), the authors utilize an implicit formulation (e.g., $\lambda = 1$) to ensure stability over long time intervals.
• Stability Analysis: The paper provides a spectral stability analysis demonstrating that for systems with a second-order Jordan block, explicit schemes (or those with $\lambda=0$) become unstable as the number of time steps increases, regardless of grid parameters. The implicit scheme is shown to be necessary for stable long-term integration.
• Model Scope: The method is applied to both the full relativistic system (using the Lorentz factor $\gamma$) and the non-relativistic approximation.
• Initial Conditions: Simulations utilize localized initial electric field perturbations (Gaussian-like profiles) to excite plasma oscillations, with zero initial velocity and momentum.

Key Results
The computational experiments, conducted for both non-relativistic and relativistic cases, yield the following findings:

1. Validation of Theory: The numerical results for constant collision coefficients ($\nu_0 = 0, \nu_1 > 0$) are in full agreement with existing theoretical results. Specifically, constant collisions lead to the damping of oscillation amplitudes and, in the relativistic case, can delay or prevent the formation of density singularities depending on the magnitude of $\nu_1$.
2. Impact of Linear Dependence ($\nu = \nu_0 N$):
   • Non-Relativistic Case: For initial data that would otherwise lead to global smooth solutions, the linear dependence enhances damping. For data leading to breaking in the collisionless case, the linear dependence prevents the breaking effect entirely for the tested parameters, suggesting an expansion of the set of initial data leading to global smooth solutions.
   • Relativistic Case: The linear dependence significantly slows down the breaking process of multi-period oscillations. In cases where breaking occurs with constant collisions, the variable coefficient delays the singularity formation. In some scenarios, it completely eliminates the breaking effect, preventing the electron density from becoming infinite.
3. Nature of Singularity: In collisionless relativistic cases, the breaking effect is confirmed to be a "gradient catastrophe," where the momentum and electric field remain bounded while their derivatives (and consequently the density) become infinite.
4. Threshold Behavior: The study identifies that while a linear dependence ($\nu \propto N$) acts as a regularizer, it is a "threshold" case. Complete suppression of breaking in the general case (beyond affine solutions) theoretically requires a dependence stronger than linear ($\nu \propto N^\gamma, \gamma > 1$). However, the numerical experiments show that the linear case is highly effective at delaying or preventing breaking for the specific localized initial conditions tested.

Significance and Claims
The paper claims significance primarily in the development of a robust numerical tool for a class of equations that are mathematically challenging due to the loss of hyperbolicity.

• Numerical Method: The primary contribution is the construction and validation of an implicit MacCormack-type scheme capable of handling non-strictly hyperbolic systems (specifically those with Jordan block structures) over long time intervals where explicit methods fail.
• Physical Insight: The study provides numerical evidence that modeling electron-ion collisions as a linear function of electron density (a physically motivated approach) qualitatively changes the solution behavior by enhancing damping and potentially preventing the breaking of plasma oscillations.
• Applicability: The method is presented as a viable tool for numerical modeling of laser-plasma interactions, particularly in regimes where collision effects are significant and the system approaches or loses hyperbolicity. The authors note that their results support the use of such regularization in cold plasma hydrodynamics, though they acknowledge that complete regularization for general initial data may require stronger dependencies than linear.

The authors maintain a modest tone regarding the theoretical justification of the linear dependence preventing breaking in the general case, noting that while their results are consistent with affine solution theories, a complete theoretical proof for general initial data remains an open area. The work is framed as a successful numerical study that bridges the gap between theoretical predictions of singularity formation and the practical need for stable computational algorithms in non-hyperbolic regimes.