On the removal of the barotropic condition in helicity studies of the compressible Euler and ideal compressible MHD equations

This paper introduces new definitions of helicity and cross-helicity densities for non-barotropic compressible Euler and ideal MHD equations that remove the restrictive barotropic pressure assumption, revealing that while global conservation is lost, the rate of change of these quantities depends solely on potential vorticity and kinetic energy, thereby enabling the derivation of an inverse resolution length scale bounded by initial potential vorticity.

The Gist

Imagine a fluid, like air or water, as a giant, invisible dance floor. On this floor, invisible threads called vortex lines swirl and twist. Sometimes, these threads get tangled into knots or linked together like chains.

For a long time, scientists studying these fluids had a major rule they couldn't break: they had to assume the fluid was "barotropic." In plain English, this means they had to pretend that the fluid's pressure and its density (how crowded the molecules are) were perfectly locked together, like two dancers holding hands so tightly they can never move apart. If the pressure changed, the density changed in a predictable way. This made the math easy, but it wasn't very realistic for real-world weather or stars, where pressure and density often act independently.

The Problem
The scientists wanted to measure a specific property of these swirling threads called helicity. Think of helicity as a score that tells you how "knotted" or "twisted" the fluid is.

• The Old Way: Under the strict "barotropic" rule, this helicity score was perfectly conserved. It was like a bank account where the total money never changed, no matter how the dancers moved.
• The Issue: When you remove that strict rule (because real fluids don't always follow it), the helicity score starts to fluctuate wildly. The math gets messy because the pressure term acts like a wild card that ruins the conservation law. It's like trying to balance a checkbook when someone keeps randomly adding or subtracting money without telling you.

The New Solution
The authors of this paper, Daniel Boutros and John Gibbon, came up with a clever new way to define the helicity score. Instead of just looking at the velocity of the fluid, they decided to weigh the measurement by the fluid's density.

• The Analogy: Imagine you are counting how many people are dancing.
  • Old Method: You just count the number of people moving.
  • New Method: You count the "mass" of the movement. If a heavy person (high density) moves, it counts for more than a light person.
  • By defining their new helicity density as (Density × Velocity) · (Density × Vorticity), they created a new score that behaves much better.

What They Found
Even though this new score isn't perfectly conserved (it doesn't stay exactly the same forever), the authors found a beautiful pattern in how it changes:

1. The Pressure Problem is Solved: In their new equation, the messy pressure terms are tucked away inside a "flux" (a flow of information). If you look at the whole system (like a closed room), these pressure terms cancel each other out. The pressure stops being a wild card.
2. The Real Culprit: The only thing that actually changes the total helicity score is something called Potential Vorticity. Think of this as a measure of how the fluid's spin interacts with its density changes.
3. The "Speed Limit": Because this potential vorticity is a "material constant" (it travels with the fluid like a passenger on a train and doesn't change its value), the authors proved that the rate at which the helicity score can change is strictly limited. It can't grow infinitely fast.

The "Resolution Length" (The Ruler)
Using this new understanding, the authors invented a new kind of ruler, which they call $\lambda_H$.

• Imagine you are looking at a blurry photo of the swirling fluid.
• This new ruler tells you the smallest detail you can possibly see before the "knots" and "twists" of the fluid start to break down or become too chaotic to track.
• They showed that this smallest detail is directly related to how "knotted" the fluid was at the very beginning. If the fluid starts with very complex knots, this ruler gets smaller, meaning the fluid can get very detailed and chaotic.

In Summary
The paper says: "We found a new way to measure the 'knottedness' of fluids that works even when the fluid isn't perfectly simple. By weighting our measurement by density, we can ignore the confusing effects of pressure and focus on the true driver of change: the interaction between spin and density. This allows us to put a strict limit on how fast the fluid's topological structure can change and gives us a new way to measure the smallest scales of fluid chaos."

They also applied this same logic to Magnetohydrodynamics (MHD), which is the study of electrically conducting fluids (like plasma in stars) interacting with magnetic fields, showing that the same "density-weighted" trick works there too.

Technical Summary

Technical Summary: On the removal of the barotropic condition in helicity studies of the compressible Euler and ideal compressible MHD equations

Problem Statement
The helicity $H = \int_V \mathbf{u} \cdot \boldsymbol{\omega} \, dV$ is a fundamental topological invariant in ideal fluid dynamics, constraining the evolution of vortex lines and linking structures. However, in the context of compressible flows, the conservation of standard helicity relies heavily on the barotropic assumption, where pressure $P$ is a function of density $\rho$ alone ($P=P(\rho)$). Under this assumption, the gradients $\nabla P$ and $\nabla \rho$ are parallel, causing the baroclinic term $\nabla(\rho^{-1}) \times \nabla P$ in the vorticity equation to vanish.

In non-barotropic scenarios (e.g., atmospheric flows with temperature gradients or stellar interiors), this term is non-zero, breaking the conservation of standard helicity. Previous attempts to define conserved quantities in non-barotropic flows (e.g., Mobbs 1981; Salmon 1988) introduced generalized helicities that depend on non-local variables, specifically the Lagrangian time integral of temperature. The authors identify a gap in the literature: the lack of a local definition of helicity for non-barotropic compressible flows that facilitates the analysis of topological dynamics without requiring global conservation or non-local variables.

Methodology
The authors propose a modification to the definition of helicity density to accommodate non-barotropic conditions. Instead of the canonical density $h = \mathbf{u} \cdot \boldsymbol{\omega}$, they introduce a new density:
$$h_\rho = (\rho \mathbf{u}) \cdot \text{curl}(\rho \mathbf{u}) = (\rho \mathbf{u}) \cdot (\rho \boldsymbol{\omega})$$
This definition is applied to three distinct systems:

1. Inhomogeneous Incompressible Euler Equations: Where density varies but velocity is divergence-free ($\nabla \cdot \mathbf{u} = 0$).
2. Fully Compressible Euler Equations: Including specific internal energy $e$ and allowing for compressibility ($\nabla \cdot \mathbf{u} \neq 0$).
3. Ideal Compressible Magnetohydrodynamics (MHD): Extending the approach to cross-helicity.

The core analytical method involves deriving the evolution equation for the new density $h_\rho$. The authors manipulate the momentum and vorticity equations to express the time derivative $\partial_t h_\rho$ in the form of an entropy-type relation (in the sense of hyperbolic conservation laws):
$$\partial_t h_\rho + \nabla \cdot \mathbf{J}_\rho = \sigma_\rho$$
where $\mathbf{J}_\rho$ is a flux vector and $\sigma_\rho$ is a source term. Crucially, the authors demonstrate that all pressure terms are contained within the flux $\mathbf{J}_\rho$. Consequently, when integrated over a periodic domain (or a domain with suitable boundary conditions), the pressure contribution vanishes, leaving the rate of change of total helicity dependent only on local kinematic and thermodynamic variables.

Key Contributions and Results

1. Removal of the Barotropic Constraint: The paper proves that by redefining helicity density as $h_\rho$, the strict requirement for a barotropic equation of state is removed. While the total integral $H = \int_V h_\rho \, dV$ is no longer strictly conserved in the non-barotropic case, its evolution is governed by a closed local equation.

2. Entropy-Type Evolution Laws:

   • Inhomogeneous Incompressible Euler: The evolution is given by $\partial_t h_\rho + \nabla \cdot \mathbf{J}_\rho = \sigma_\rho$, where the source term is $\sigma_\rho = -q \rho |\mathbf{u}|^2$. Here, $q = \boldsymbol{\omega} \cdot \nabla \rho$ is the potential vorticity.
   • Fully Compressible Euler: The source term includes an additional term involving the divergence of velocity: $\tilde{\sigma}_\rho = \sigma_\rho - 2h_\rho \nabla \cdot \mathbf{u}$.
   • Ideal MHD: A similar non-barotropic cross-helicity density $h_c = \rho \mathbf{u} \cdot \mathbf{B}$ is introduced, satisfying a similar evolution law involving the magnetic potential vorticity $q_c = \mathbf{B} \cdot \nabla \rho$.

3. Boundedness and Length Scale Estimation:

   • For the inhomogeneous incompressible case, the potential vorticity $q$ is shown to be a material constant ($Dq/Dt = 0$). This implies that $\|q(\cdot, t)\|_\infty \leq \|q_0\|_\infty$.
   • Using this bound, the authors derive an upper bound for the rate of change of helicity: $|dH/dt| \leq 2 \|q_0\|_\infty E_0$, where $E_0$ is the total kinetic energy.
   • This bound allows for the definition of an inverse resolution length scale $\lambda_H^{-1}$, analogous to the Kolmogorov scale, which characterizes the smallest scale of significant topological variation. The paper establishes the bound:
     $$\lambda_H^{-1} \leq \left( \frac{4}{E_0 \rho_0} \right)^{1/7} \|q_0\|_\infty^{2/7}$$
     This result suggests that the growth of topological complexity is bounded by the initial potential vorticity and energy.

4. Generalization to MHD: The methodology successfully extends to ideal compressible MHD, providing a local cross-helicity density whose evolution is similarly decoupled from explicit pressure terms in the integral sense.

Significance and Claims
The authors claim that this work provides a local framework for analyzing helicity dynamics in non-barotropic compressible flows, overcoming the limitations of previous non-local generalizations. The primary significance lies in:

• Analytical Tractability: The new helicity density obeys an equation where pressure terms are relegated to a flux, making the global rate of change independent of the specific equation of state (barotropic or otherwise).
• Topological Insight: The dynamics of the new helicity are shown to be driven by the potential vorticity (and velocity divergence in the compressible case), offering a clearer physical interpretation of how baroclinicity affects topological structures.
• Quantitative Bounds: For inhomogeneous incompressible flows, the work provides a rigorous upper bound on the inverse length scale of helicity variations, linking topological evolution directly to the initial potential vorticity and energy.

The paper explicitly notes that these results are valid only for time intervals where sufficiently smooth solutions exist, acknowledging that finite-time singularities (shocks or blow-up) could invalidate the manipulations required. The authors do not claim to solve the global existence problem but rather to provide a robust diagnostic tool for the dynamics of smooth solutions in non-barotropic regimes.