Numerical Bow Shock Instabilities in Inert Polyatomic Gases

This paper identifies and characterizes a previously undocumented class of numerically induced carbuncle-type instabilities in hypersonic bow shock simulations of inert polyatomic gases with low specific heat ratios, highlighting the critical need to distinguish these computational artifacts from genuine physical shock oscillations observed in real thermochemically relaxing flows.

The Gist

The Big Picture: When Computers "Hallucinate" Instability

Imagine you are a weather forecaster trying to predict a hurricane. You have a supercomputer running a simulation. Suddenly, the computer shows a massive, swirling storm front that is wobbling and shaking violently.

You might think, "Wow, nature is chaotic today!" But what if the storm isn't real? What if the shaking is just a glitch in your computer code caused by the way you drew your map?

This is exactly what G. S. Sidharth and A. Dwivedi discovered in their study. They found that when simulating high-speed gas flowing around a sphere (like a spacecraft entering an atmosphere), the computer can invent fake, wavy instabilities that look like real physics but are actually just numerical errors.

The Setting: The "Heavy" Gas and the "Bumpy" Map

To understand the problem, we need two ingredients:

1. The "Heavy" Gas (Low Gamma):
   Think of the gas molecules as dancers. In normal air, they are light and energetic. In the gases studied here (like propane or CO2 at high speeds), the molecules are "heavy" and sluggish. When a shockwave hits them, they get squished together very tightly. This creates a huge jump in density, like a traffic jam where cars go from moving fast to bumper-to-bumper instantly.

   • The Analogy: Imagine a highway where cars suddenly turn into giant, heavy boulders. The pile-up is massive.

2. The "Bumpy" Map (Coarse Grid):
   Computers don't see a smooth curve; they see a grid of tiny squares (pixels). To save time, scientists sometimes use a "coarse" grid where the squares are big.

   • The Analogy: Imagine trying to draw a perfect circle using only large Lego blocks. The circle looks jagged and bumpy.

The Problem: The "Ghost Wave"

When the researchers combined the "heavy gas" (huge density jump) with the "bumpy map" (coarse grid), something strange happened.

In real physics, if you shoot a gas at a sphere, the shockwave (the invisible wall of compressed gas) should sit still and look smooth. But in their computer simulation, the shockwave started wiggling and traveling around the sphere like a snake.

The Twist:
Usually, scientists know about a computer error called the "Carbuncle." This is a static, jagged glitch that sits right in the center of the shock. It's like a permanent bump on the nose of the sphere.

• The Discovery: This paper found something new. Instead of a static bump, the computer created traveling waves. The shockwave started rippling and moving side-to-side, mimicking a real physical instability that scientists have seen in experiments with real gases.

Why This is Dangerous (The "Fake News" of Physics)

This is tricky because real gases do sometimes become unstable. In real life, when gas gets super hot, the molecules break apart (dissociate), which absorbs energy and causes the shockwave to wobble. This is a real physical phenomenon.

However, the researchers showed that you can get the exact same wobbling in a computer simulation just by using a "heavy" gas model and a "bumpy" grid, even if the gas is chemically inert (nothing is breaking apart).

• The Metaphor: It's like a musician playing a note.
  • Real Instability: The musician is actually playing a wobbly note because their instrument is broken.
  • Numerical Instability: The musician is playing a perfect note, but the microphone is picking up a weird hum from the electrical wiring.
  • The Danger: If you don't check your equipment, you might think the musician is playing a wobbly note and blame the instrument, when it was just the wiring.

The Solution: Sharpening the Pencil

The researchers tested how to stop these fake waves:

1. Refine the Grid (Zoom In): When they made the "Lego blocks" smaller (finer grid), the fake waves got smaller and eventually disappeared.
   • Lesson: If the wave vanishes when you zoom in, it was probably a computer glitch, not real physics.
2. The "Entropy Fix" (The Stabilizer): They used a mathematical trick (an eigenvalue limiter) to force the computer to be a little more "dissipative" (smoother) right at the shock. This acted like a shock absorber, stopping the fake waves from forming.

The Takeaway: How to Tell Real from Fake

The paper concludes with a warning for scientists: Just because you see a wiggling shockwave in a simulation, it doesn't mean it's real.

To know if an instability is real, you must check for Grid Convergence:

• Real Physics: If you zoom in (make the grid finer), the wave stays there, but it gets sharper. The frequency of the wobble stays the same.
• Computer Glitch: If you zoom in, the wave gets smaller and eventually disappears. The frequency changes wildly.

Summary in One Sentence

This paper warns us that when simulating high-speed gas flows, a computer can invent fake, traveling shockwaves that look exactly like real physics, but these "ghost waves" are actually just artifacts of using a low-resolution grid, and they vanish when you sharpen your digital tools.

Technical Summary

1. Problem Statement

The paper addresses a critical ambiguity in computational fluid dynamics (CFD) simulations of hypersonic flow over blunt bodies: distinguishing between genuine physical instabilities and numerical artifacts.

• Context: In experiments involving polyatomic gases (e.g., $CO_2$, $CF_4$, propane) at high Mach numbers, bow shocks often exhibit large-amplitude oscillations and unsteady spatial structures. These are physically attributed to thermochemical relaxation (dissociation) and high density ratios ($\rho_2/\rho_1$) across the shock, which drive baroclinic vorticity generation and Kelvin-Helmholtz instabilities.
• The Issue: When simulating these flows, researchers often model the gas as "inert" (calorically perfect) but with an effectively low specific heat ratio ($\gamma \approx 1.1–1.2$) to mimic the high density ratios observed in real gases.
• The Gap: It was unclear whether the unsteadiness observed in such simulations was a physical phenomenon or a numerical error. The authors investigate whether inviscid, inert gas simulations with high density ratios and coarse meshes can spontaneously generate traveling-wave instabilities purely due to numerical factors, specifically the carbuncle phenomenon and grid stretching.

2. Methodology

The authors conducted a systematic numerical investigation using the following setup:

• Governing Equations: Axisymmetric Euler equations (inviscid) and Navier-Stokes equations (viscous) were solved.
• Flow Conditions: Two cases were simulated to represent high density ratios:
  1. Mars Science Laboratory entry conditions: $M_\infty = 28$ in $CO_2$ ($\gamma \approx 1.15$, density ratio $\approx 14$).
  2. Propane conditions: $M_\infty = 11$ with $\gamma \approx 1.08$ (density ratio $> 20$).
• Numerical Scheme: A second-order finite-volume method with a Roe-type Riemann solver and MUSCL reconstruction. No artificial dissipation or entropy fixes were initially applied to allow instabilities to manifest naturally.
• Grid Design: Structured, body-fitted grids with stretched spacing in the wall-normal direction (to resolve the boundary layer).
  • Grid resolutions tested: Coarse ($\Delta s \approx 0.1$), Medium ($\Delta s \approx 0.05$), and Fine ($\Delta s \approx 0.025$), where $\Delta s$ is the cell size at the shock normalized by standoff distance.
  • Comparison between uniform grids and stretched grids.
• Control Mechanism: The study utilized an eigenvalue limiter (entropy fix) applied to the dissipative term of the flux. This limiter prevents the entropy-wave eigenvalue from vanishing at the stagnation point, a known trigger for carbuncle instabilities. The parameter $\varepsilon_0$ (ranging from 0.2 to 0.4) controlled the strength of this limiter.

3. Key Contributions

This work provides the first documentation of well-organized traveling-wave patterns arising purely from numerical artifacts in the absence of chemistry or real-gas effects. Key contributions include:

1. Identification of a New Instability Mode: Demonstrating that inert, low-$\gamma$ flows on stretched grids can produce traveling wave instabilities that mimic physical shock oscillations.
2. Grid-Dependence Characterization: Proving that these instabilities are grid-driven. Unlike physical instabilities (anchored by chemical induction lengths), these numerical waves vanish or change scale significantly with mesh refinement.
3. Role of Grid Stretching: Revealing that grid stretching (common in aerothermodynamics) qualitatively changes the instability from a steady "kink" (carbuncle) to a propagating traveling wave due to spatially varying truncation errors.
4. Diagnostic Framework: Establishing that grid convergence of oscillation frequency is the primary metric to distinguish physical shock unsteadiness from numerical noise.

4. Key Results

• Instability Manifestation:

  • On coarse grids ($\Delta s \approx 0.1$) with low $\gamma$, the bow shock developed visible corrugations and unsteady traveling waves propagating circumferentially.
  • These waves were absent in high-$\gamma$ ($\gamma=1.4$) cases on the same grid, confirming the dependence on high density ratios.

• Effect of Eigenvalue Limiter:

  • Low Limiter ($\varepsilon_0 = 0.2$): Resulted in a steady carbuncle deformation (localized kink) at the stagnation point.
  • Intermediate Limiter ($\varepsilon_0 = 0.3$): The instability transitioned to a traveling-wave regime, with waves propagating along the shock layer.
  • High Limiter ($\varepsilon_0 = 0.4$): The shock layer was largely stabilized, with only residual perturbations.
  • Conclusion: The instability is controllable via standard numerical stabilization techniques used for carbuncles.

• Grid Resolution Dependence:

  • Coarse Grid: Instability amplitude was high and persisted across all limiter values.
  • Fine Grid: Instability amplitude dropped by several orders of magnitude. The traveling waves became smaller in amplitude and wavelength.
  • Spectral Analysis: The frequency of the traveling waves decreased by a factor of ~2 when the grid was refined, confirming the waves are seeded by the grid rather than physical flow physics.

• Viscous vs. Inviscid:

  • Navier-Stokes simulations showed no significant change in instability behavior compared to inviscid Euler simulations.
  • This confirms the instability is fundamentally inviscid, driven by the interaction of large density gradients, shock curvature, and coarse mesh resolution, rather than boundary layer effects.

• Uniform vs. Stretched Grids:

  • Uniform Grids: Produced a steady-state carbuncle (kink) without propagation.
  • Stretched Grids: Produced traveling waves. The spatial variation in cell size creates an inhomogeneous truncation error that acts as a forcing mechanism, driving the wave propagation.

5. Significance and Implications

• Simulation Caution: The presence of shock-layer corrugations in a simulation does not automatically confirm a physical instability. Researchers must be cautious when modeling polyatomic gases as inert fluids with low $\gamma$, as they may inadvertently trigger numerical artifacts that mimic physical phenomena.
• Validation Protocol: The paper emphasizes that grid convergence studies are essential. Physical instabilities (like those in shock-induced combustion) maintain a characteristic frequency and scale anchored by chemistry, changing by only a few percent with mesh refinement. Numerical instabilities, conversely, change drastically or vanish upon refinement.
• Future Directions: The authors suggest that active molecular rate processes (dissociation/combustion) observed in real experiments likely play a role beyond just raising the density ratio, potentially providing the necessary physical "anchor" that prevents these numerical artifacts from dominating. Future work will focus on adjoint-based sensitivity analysis to identify the "wavemaker" of these instabilities.

In summary, this paper serves as a critical warning to the hypersonics community: numerical artifacts can masquerade as complex physical instabilities in high-density-ratio flows, and rigorous grid-convergence testing is required to validate simulation results.