Sub-Kolmogorov Intermittency and Multifractal Dissipation in Multiphase Turbulence

Through direct numerical simulations, this study reveals that in multiphase turbulence, interface breakup and coalescence drive a distinct multifractal organization of dissipation, causing intense energy dissipation events to extend deep into the sub-Kolmogorov range and significantly broaden the local dissipative cutoff compared to single-phase turbulence.

The Gist

Imagine a pot of water boiling on a stove. If you just have water (single-phase), the bubbles and swirls are chaotic, but they follow a somewhat predictable pattern of size and energy. Now, imagine adding oil to that water and stirring it vigorously. You get a messy mix of droplets, streams, and bubbles constantly forming, merging, and splitting apart. This is multiphase turbulence.

This paper investigates what happens at the tiniest, most invisible levels of that chaotic mix. The researchers wanted to understand why the "smallest swirls" in a mix of liquids behave so differently—and more violently—than in a single liquid.

Here is the story of their discovery, broken down into simple concepts:

1. The "Safety Net" That Doesn't Exist

In normal fluid physics, there is a theoretical "safety net" called the Kolmogorov scale. Think of this as the smallest size a whirlpool can get before the fluid's natural stickiness (viscosity) smooths it out and kills the energy. In a single liquid, the energy stops there.

However, the researchers found that in a mix of liquids (like oil and water), this safety net is broken.

• The Analogy: Imagine a trapeze artist (the energy) swinging. In a single liquid, they stop swinging at a specific height. In a liquid mix, the trapeze artist keeps swinging way lower, deep into a zone where physics says they should have stopped.
• The Finding: The "dissipative cutoff" (the point where energy dies) doesn't just stop at the usual limit; it stretches deep into a "sub-Kolmogorov" range. The energy fluctuations become much more intense and extreme than anyone expected.

2. The Culprits: Breaking and Merging

Why does this happen? The paper identifies the specific "crime scenes" where this extreme energy is generated.

• The Analogy: Think of a crowd of people moving randomly. If two people bump into each other and merge, or if a group splits apart, it causes a sudden, chaotic jolt.
• The Finding: The most intense, tiny-scale energy bursts happen specifically at the interfaces where the liquids meet. Specifically, they occur when droplets break apart (breakup) or smash together and merge (coalescence).
• These events create sharp curves and sudden changes in speed that the fluid can't smooth out easily, forcing the energy to go deeper and deeper into the microscopic realm.

3. The "Fractal" Geometry of Chaos

The researchers used a mathematical tool called multifractal analysis.

• The Analogy: Imagine looking at a coastline. From far away, it looks like a line. Up close, it's jagged. Up closer, it's full of bays and rocks. A "fractal" is a shape that looks complex at every level of zoom.
• The Finding: In a single liquid, the "roughness" of the energy distribution is fairly consistent. But in a liquid mix, the geometry of the chaos changes completely at the smallest scales.
  • The "roughness" becomes much more extreme.
  • The most violent energy events aren't spread out evenly; they are concentrated on very thin, thread-like structures (like the neck of a droplet right before it snaps).
  • The paper describes these intense events as being supported on "sparse structures," meaning they are rare, isolated, and incredibly sharp, rather than a general fog of turbulence.

4. The Prediction Machine

The researchers didn't just observe this; they proved they could predict it.

• They used the mathematical "shape" of the chaos (the singularity spectrum) to predict exactly how often these extreme, tiny events would happen.
• The Result: When they looked at the "near" and "sub-Kolmogorov" zones (the deep, tiny scales), their predictions matched the computer simulations perfectly. This confirms that the strange, extreme behavior is a direct result of the liquid interfaces breaking and merging.

The Bottom Line

The paper concludes that when you mix two liquids, the turbulence doesn't just get "a little messier." The act of droplets breaking and merging fundamentally rewrites the rules of the smallest scales. It creates a new, distinct type of chaotic geometry where the most violent energy events are locked into thin, thread-like regions at the interface.

In short: Breaking and merging droplets don't just disturb the flow; they imprint a unique, extreme, and highly organized pattern of chaos onto the very smallest scales of the fluid.

Technical Summary

Technical Summary: Sub-Kolmogorov Intermittency and Multifractal Dissipation in Multiphase Turbulence

Problem Statement
Multiphase turbulent flows exhibit stronger intermittency than their single-phase counterparts, yet the physical origin and geometrical organization of their most intense small-scale fluctuations remain poorly understood. While it is established that surface tension mediates a scale-dependent exchange between turbulent kinetic energy and interfacial energy—modifying the classical energy cascade—the specific impact of topology-changing interfacial events (breakup and coalescence) on the statistics of the dissipative field and the local dissipative cutoff is unclear. A critical open question is how these events influence the smallest dynamically active scales, where a significant fraction of mixing and dissipation occurs, and whether they alter the multifractal organization of the flow.

Methodology
The authors employ Direct Numerical Simulations (DNS) of the incompressible Navier–Stokes equations coupled with surface tension. The interface dynamics are captured using a Volume-of-Fluid (VOF) method implemented in the FLUTAS code.

• Simulation Setup: Simulations are conducted in a triply periodic box ($L=2\pi$) with a grid resolution of $512^3$ points. Turbulence is maintained by large-scale forcing.
• Configurations: Two cases are compared: a single-phase reference flow ($Re_\lambda \approx 50$, $\eta_K \approx 6.8 \Delta x$) and a multiphase flow with identical density and viscosity for carrier and dispersed phases ($\psi = 0.1$, $We(d_v) = 690$). The matching of density and viscosity isolates the effects of interfacial dynamics and surface tension.
• Analytical Framework: The study combines an analysis of fluctuating dissipative scales with a multifractal formalism.
  • Local Dissipative Scale: Defined via a scale-dependent local Reynolds number, $Re(x, \eta) \equiv \eta |\delta u_\eta(x)| / \nu \sim 1$, to characterize the probability density function (PDF), $Q(\eta)$, of the local dissipative cutoff.
  • Multifractal Analysis: The dissipation field is treated as a singular measure. The study computes generalized R'enyi dimensions ($D_q$) and singularity spectra ($f(\alpha)$) over two distinct scale intervals: Above Kolmogorov (AK, $2.4 < r/\eta_K < 20$) and Below Kolmogorov (BK, $0.3 < r/\eta_K < 2.4$).
  • Validation: The multifractal spectrum is used to theoretically reconstruct the distribution of cutoff lengths, $Q(\eta)$, to verify the predictive capability of the approach.

Key Results

1. Broadening of the Dissipative Cutoff: In multiphase turbulence, the PDF of the local dissipative scale, $Q(\eta)$, is markedly broader than in single-phase turbulence. Specifically, the left tail is substantially enhanced, indicating a high probability of events where the local dissipative scale $\eta$ is significantly smaller than the global Kolmogorov scale ($\eta \ll \eta_K$). This confirms that interfaces amplify the intermittency of dissipative dynamics, allowing intense fluctuations to penetrate deep into the sub-Kolmogorov range.
2. Spatial Concentration of Events: The intense sub-Kolmogorov events are not randomly distributed but are spatially concentrated around topology-changing interfacial regions, specifically breakup and coalescence. These regions are characterized by large interface curvature and extreme small-scale velocity fluctuations.
3. Multifractal Spectrum Divergence:
   • Above Kolmogorov (AK): The singularity spectrum for the multiphase case remains close to the single-phase case for moderate and large $\alpha$, suggesting that weakly dissipative regions are only mildly affected by the interface.
   • Below Kolmogorov (BK): The multiphase case exhibits a qualitatively different behavior. The spectrum develops a markedly broader singularity spectrum with strong left tails, supporting very small values of $\alpha$ on sets with finite fractal dimension. The value $f(\alpha \approx 0) \approx 1$ suggests that the most singular dissipative structures are supported on approximately one-dimensional sets (e.g., thinning ligaments and reconnecting threads).
4. Predictive Capability: Using the BK-derived singularity spectrum, the theoretical reconstruction of $Q(\eta)$ matches the DNS data with excellent agreement. Conversely, using the AK spectrum fails to capture the sub-Kolmogorov events, demonstrating that the multifractal approach is predictive only when the appropriate scale intervals (specifically those capturing interfacial effects) are selected.

Significance and Claims
The paper establishes that breakup and coalescence in multiphase turbulence do not merely perturb the flow locally; they imprint a distinct multifractal geometry on the dissipation field. The findings demonstrate that interfacial dynamics sustain a strongly multifractal dissipative regime in the vicinity of and below the Kolmogorov scale, which is fundamentally different from the inertial-range intermittency of single-phase turbulence. By linking the local Reynolds number fluctuations to the singularity spectrum, the study provides a direct connection between interfacial topology changes and the extreme small-scale activity of multiphase turbulence. This work extends the multifractal formalism to multiphase flows, offering a framework to describe how interface-mediated energy transfer reshapes the statistical geometry of dissipation.