The Origin of Da Scaling: Suppressed Cooling in Fast-Cooling Mixing Layers

This paper explains the transition in radiative cooling scaling from $\dot{E}_{\rm cool} \propto {\rm Da}^{1/2}$ to $\dot{E}_{\rm cool} \propto {\rm Da}^{1/4}$ in fast-cooling turbulent mixing layers as a result of ram pressure from inflowing gas suppressing the turbulent folding and fractal structure of the interface.

The Gist

Imagine two giant rivers of gas flowing past each other in space: one is a hot, thin river, and the other is a cold, thick river. Where they meet, they don't just slide past; they churn, mix, and create a turbulent "mixing layer." As these gases mix, they get hot enough to glow and radiate away energy as light. This process is called a Turbulent Radiative Mixing Layer (TRML).

For a long time, scientists thought they understood how fast this energy is lost. They believed that if the gas cooled down very quickly (a "fast-cooling" regime), the amount of light emitted would follow a specific mathematical rule. However, new simulations by Lachlan Lancaster and his team have uncovered a twist: the rule changes, and the reason is surprisingly physical.

Here is the story of their discovery, explained simply.

The Two Regimes: Stirring vs. Folding

To understand the discovery, imagine you are trying to mix a drop of dye into a glass of water.

1. The Slow-Cooling Regime (The "Stirred Reactor"):
   If the dye takes a long time to disappear (cool down), the swirling water has plenty of time to mix it thoroughly. The turbulence acts like a giant spoon, smoothing out the boundary between the hot and cold gas. In this scenario, the faster the turbulence swirls, the more energy is radiated. The relationship is straightforward: more turbulence equals more cooling.

2. The Fast-Cooling Regime (The "Fractal Fold"):
   Now, imagine the dye disappears almost instantly. The water swirls, but before it can smooth things out, the dye vanishes. In this case, the turbulence doesn't smooth the surface; instead, it crumples and folds it, like a piece of paper being crumpled into a ball. This creates a massive amount of surface area (a "fractal" structure) where the hot and cold gases touch. Because there is so much surface area, the gas cools very efficiently.

Scientists expected that even in this "fast-cooling" regime, the cooling rate would keep increasing in a predictable way as the turbulence got stronger. But the simulations showed something different: the cooling rate started to grow much more slowly than expected.

The Discovery: The "Wind" Stops the Folding

The paper asks: Why does the cooling rate slow down when the gas cools very fast?

The authors found the answer lies in the inflow of gas. To keep the mixing layer going, hot gas must constantly flow in to replace the gas that has cooled and fallen out.

• The Analogy: Imagine a strong wind blowing against a pile of dry leaves.
  • When the wind is gentle (Low "Damköhler number"): The wind isn't strong enough to stop the leaves from tumbling and folding over each other. The pile stays messy and has a huge surface area.
  • When the wind is a hurricane (High "Damköhler number"): The wind is so powerful that it smashes the leaves flat against the ground. It suppresses the tumbling and folding. The pile becomes smooth and flat, losing all that extra surface area.

In the paper's language:

• The "wind" is the ram pressure of the inflowing hot gas.
• The "tumbling leaves" are the turbulent folds of the mixing layer.
• When the cooling is extremely fast, the inflow of gas becomes so violent that its pressure crushes the turbulent folds. The interface between the hot and cold gas stops being a crumpled, high-surface-area fractal and becomes a smoother, flatter surface.

Because the surface area shrinks, the gas has less "skin" to radiate energy from, so the total cooling rate drops below what scientists previously predicted.

The "Damköhler Number" (The Speedometer)

The paper uses a specific number called the Damköhler number (Da) to measure this. Think of it as a speedometer comparing two things:

1. How fast the turbulence swirls (the eddy turnover time).
2. How fast the gas cools down (the cooling time).

• Low Da: Cooling is slow; turbulence wins and smooths the surface.
• High Da: Cooling is fast; turbulence tries to fold the surface, but the inflow pressure wins and flattens it.

The authors show that the transition where the cooling rate changes its behavior happens exactly when the pressure of the incoming gas becomes stronger than the pressure of the turbulence itself.

What This Means for the Math

Previously, some theories suggested that the change in cooling rate was due to complex changes in how heat diffuses through the gas. The authors argue that this is incorrect.

Instead, they propose a new, simpler explanation:

1. The cooling rate depends on how much surface area exists between the hot and cold gas.
2. In the fast-cooling regime, the incoming gas acts like a heavy hand, pressing down on the turbulence.
3. This pressure suppresses the "fractal" (crumpled) nature of the surface, reducing the area available for cooling.
4. This suppression perfectly explains why the cooling rate follows a new, slower mathematical rule (scaling with the 1/4 power instead of the 1/2 power).

Summary

In short, the paper reveals that in the universe's most energetic mixing layers, you can't have your cake and eat it too. If the gas cools down too quickly, the force required to keep feeding that cooling process (the inflow) becomes so strong that it smashes the turbulent folds. This flattens the interface, reduces the surface area, and slows down the total energy loss. The "fast-cooling" regime isn't just about speed; it's about the suppression of chaos by the sheer force of the inflowing gas.

Technical Summary

Technical Summary: The Origin of Da Scaling: Suppressed Cooling in Fast-Cooling Mixing Layers

Problem Statement
In numerical simulations of Turbulent Radiative Mixing Layers (TRMLs), a distinct change in the scaling behavior of the total radiative cooling rate ($\dot{E}_{cool}$) is observed as the Damköhler number ($Da \equiv t_{eddy}/t_{cool}$) increases. In the "slow-cooling" regime ($Da < 1$), the cooling rate scales as $\dot{E}_{cool} \propto Da^{1/2}$. However, in the "fast-cooling" regime ($Da > 1$), which is relevant to many astrophysical environments, the scaling transitions to $\dot{E}_{cool} \propto Da^{1/4}$. While this $Da^{1/4}$ scaling is empirically established in both astrophysical and turbulent combustion literature, the physical mechanism driving this transition has been debated. Previous explanations have relied on concepts such as an "effective" cooling time (Tan et al. 2021) or a balance between turbulent wrinkling and a specific fractal dimension and Kolmogorov scaling (Fielding et al. 2020). This paper seeks to identify the precise physical origin of this scaling transition.

Methodology
The authors utilize a suite of high-resolution hydrodynamical simulations of TRMLs performed with the GPU-accelerated AthenaK code. The simulations model the interface between hot, low-density gas and cold, high-density gas in relative motion, subject to a radiative cooling function that dominates at intermediate temperatures.

• Parameters: The study varies the density contrast ($\chi = \rho_{cold}/\rho_{hot}$), the Mach number of the shear flow ($M$), and the cooling strength (parameterized by $\xi$, which is closely related to $Da$).
• Measurements: The authors measure the total cooling rate ($\dot{E}_{cool}$), the bulk inflow velocity ($v_{bulk}$), the integral-scale turbulent velocity ($v_t(L_{int})$), and the interface surface area ($A_{int}$) across different scales.
• Analysis: They analyze the scaling of these quantities with $Da$ and investigate the fractal dimension ($d$) of the interface and the anisotropy of the turbulent velocity field. The study explicitly checks for resolution independence and compares simulation measurements against theoretical models.

Key Results

1. Transition Condition: The transition from $\dot{E}_{cool} \propto Da^{1/2}$ to $\dot{E}_{cool} \propto Da^{1/4}$ occurs precisely when the bulk inflow velocity required to balance cooling ($v_{bulk}$) exceeds the integral-scale turbulent velocity ($v_t(L_{int})$).
2. Suppression of Interface Area: In the fast-cooling regime ($Da \gg 1$), the total cooling rate is determined by the suppression of the interface surface area ($A_{int}$). While $A_{int}$ remains roughly constant in the slow-cooling regime, it becomes progressively suppressed at high $Da$.
3. Physical Mechanism: The suppression of $A_{int}$ is caused by the ram pressure of the inflowing gas ($\rho_{hot}v_{bulk}^2$) disrupting the turbulent folding of the interface. When $v_{bulk} \gtrsim v_t(L_{int})$, the ram pressure dominates the turbulent pressure, preventing the turbulence from wringing the interface into a multi-scale, fractal structure on large scales.
4. Fractal Dimension and Anisotropy: The study finds that the excess fractal dimension ($d$) of the interface, which is approximately $1/2$ on small scales, is suppressed on large scales in the fast-cooling regime. This suppression correlates directly with the anisotropy of the velocity field induced by the inflow. Furthermore, the turbulent velocity structure function does not follow the Kolmogorov $v_t(\ell) \propto \ell^{1/3}$ scaling assumed in some previous models; instead, it follows a shallower power law ($\propto \ell^{1/5}$).
5. Rejection of Previous Models: The authors demonstrate that the Fielding et al. (2020) model, which relies on Kolmogorov scaling and a constant fractal dimension across all scales, does not hold in their simulations. Similarly, the Tan et al. (2021) model, which posits a specific "effective" cooling time, does not account for the observed suppression of $A_{int}$.

New Theoretical Framework
The authors propose a new interpretation based on the suppression of the fractal structure. They derive a relation where the scale at which fractal structure is realized ($L_{frac}$) is suppressed relative to the outer turbulent scale ($L_{int}$) by the ratio of the inflow ram pressure to the turbulent pressure.

• By assuming the fractal dimension $d=1/2$ applies only below the suppressed scale $L_{frac}$, and incorporating the ram-pressure scaling ($L/L_{frac} \propto v_{bulk}^2/v_t^2$), they derive:
  $$\frac{\dot{E}_{cool}}{v_t(L_{int})} \propto Da^{\frac{1+2d}{2}}$$
• Substituting $d=1/2$ recovers the observed $\dot{E}_{cool} \propto Da^{1/4}$ scaling.
• The model predicts that deviations from $Da^{1/4}$ should occur if the fractal dimension $d$ differs from $1/2$ or if the scaling of $v_t$ with scale differs from the assumptions made.

Significance
The paper claims to resolve the origin of the $Da^{1/4}$ scaling in fast-cooling mixing layers by identifying the ram-pressure suppression of turbulent folding as the governing mechanism. This finding challenges previous interpretations that relied on effective reaction times or constant fractal dimensions across all scales. The authors conclude that the $Da^{1/4}$ scaling is a consequence of the inflow dynamically suppressing the large-scale fractal structure of the interface, a result that self-consistently explains the transition observed in numerical experiments. They note that while the $Da^{1/4}$ scaling is robust in the tested regime, its applicability to extreme parameters ($Da \gg 10^3$) requires further verification, as the measured power-law indices in some simulations slightly deviate from $1/4$.