Turbulence closure in the light of phase transition

✨

This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are standing by a river, watching a smooth stream of water pour out of a pipe into a calm lake. At first, the water coming out is a tight, organized beam. But as it travels further, it starts to swirl, mix, and churn, eventually blending completely with the lake water. This chaotic mixing is called turbulence.

For over a century, scientists have struggled to write a perfect "rulebook" (mathematical equations) to predict exactly how this swirling water behaves. The old rulebooks worked okay, but they were like using a blurry map to navigate a stormy sea.

This paper introduces a brand new rulebook based on a surprising idea: Turbulence isn't just chaos; it's actually a "phase transition," just like water turning into ice or a magnet suddenly snapping into alignment.

Here is the breakdown of the paper's big ideas, translated into everyday language:

1. The Big Idea: Turbulence is Like a "Light Switch"

In physics, a "phase transition" is when a material suddenly changes its state.

Analogy: Think of a magnet. When it's hot, the tiny magnetic atoms inside are jiggling randomly (disordered). As it cools down to a specific "critical temperature," they suddenly snap into perfect alignment (ordered). That moment of snapping is a phase transition.

The author, Mohammed Azim, suggests that turbulence is the same thing.

The "Off" State: Smooth, calm water (laminar flow) is like the hot, disordered magnet.
The "On" State: Chaotic, swirling water (turbulent flow) is like the cold, ordered magnet.
The Switch: When water speeds up past a certain point, it doesn't just get "a little messy." It undergoes a sudden, fundamental shift in its nature, similar to water freezing into ice.

2. The Old Way vs. The New Way

The Old Way (The Eddy Viscosity):
Scientists used to treat turbulence like a thick, sticky syrup. They assumed the swirling water acted like a single, giant "eddy" (a big swirl) that slowed things down.

Metaphor: Imagine trying to predict traffic by assuming every car is just a giant, slow-moving blob. It's a rough guess, but it misses the details of individual cars swerving.

The New Way (The Phase Transition):
This paper treats turbulence as a complex system with its own "free energy" (a measure of how much the system wants to change).

Metaphor: Instead of seeing the traffic as a blob, the author sees the cars as individual drivers reacting to a "critical point." When the road gets crowded enough, the drivers suddenly switch from driving calmly to driving frantically. The math now captures that switch rather than just the average speed.

3. The "Symmetry" Secret

One of the coolest findings in the paper is about symmetry.

In the new math, the "stress" (the force the swirling water exerts) behaves in a very specific, predictable pattern.
The Analogy: Imagine a seesaw.
- Even Symmetry: If you push down on the left, the right goes up in a perfect mirror image. This happens with the "normal" forces in the water.
- Odd Symmetry: If you push the left side, the right side reacts in a way that flips the sign (like a mirror image that is also inverted). This happens with the "shear" forces (the twisting forces).
The paper found that the math of turbulence naturally creates these perfect mirror patterns, just like the physics of magnets or ice. This proves that turbulence follows the same deep laws of nature as other phase transitions.

4. Putting It to the Test: The Jet Engine

To see if this new theory works, the author built a computer simulation of a plane jet (like a jet of air shooting out of a nozzle).

The Setup: They shot a stream of air and watched how it spread out and mixed with the surrounding air.
The Result: The computer used the new "Phase Transition" rules to calculate the flow.
The Outcome: The results matched real-world experiments and high-level supercomputer simulations almost perfectly.
- It correctly predicted how fast the jet slows down.
- It correctly predicted how wide the jet spreads.
- It correctly predicted how the swirling forces behave.

5. Why Does This Matter?

Think of this like upgrading from a paper map to a GPS.

Current Models: Good enough for driving to the grocery store, but they might get you lost in a complex city or a storm.
This New Model: Because it understands the fundamental nature of the switch from calm to chaotic, it is more accurate and robust.

In a nutshell:
This paper argues that we should stop trying to force turbulence into old, clumsy boxes. Instead, we should recognize that turbulence is a natural, sudden shift in the state of matter, governed by the same elegant laws that turn water to ice or make magnets work. By using these laws, we can predict how fluids move with much greater precision, which helps engineers design better airplanes, cars, and weather models.

1. Problem Statement

Turbulent flows are governed by the Navier-Stokes equations, but direct numerical simulation (DNS) of all scales is computationally prohibitive for high Reynolds numbers in complex geometries. Consequently, Reynolds-Averaged Navier-Stokes (RANS) equations are used, which introduce unknown Reynolds stresses (correlations of fluctuating velocities).

The Core Issue: The "closure problem" requires modeling these Reynolds stresses to close the system of equations.
Limitations of Current Models: Traditional models rely on the Boussinesq hypothesis, which assumes a linear relationship between Reynolds stresses and mean velocity gradients via a scalar eddy viscosity ( $\nu_T$ ). This approach often fails to capture complex anisotropic behaviors and the fundamental physics of the transition from laminar to turbulent flow.
The Gap: There is a need for a closure model that treats turbulence not just as a statistical phenomenon but as a physical process analogous to continuous (second-order) phase transitions, utilizing concepts like order parameters and free energy symmetry.

2. Methodology

The study proposes a novel theoretical framework and numerical implementation:

A. Theoretical Derivation (Phase Transition Analogy)

Analogy: The author draws a parallel between the onset of turbulence and a continuous phase transition (CPT). In CPT, a system undergoes qualitative changes in macroscopic properties as it crosses a critical point, governed by an order parameter ( $\Phi$ ).
Order Parameter: The mean velocity ( $u$ ) is identified as the control parameter (order parameter) for the fluid system.
Free Energy Formulation: The Reynolds stresses ( $\overline{u_i' u_j'}$ ) are treated as the free energy of the turbulence. The study derives closure equations by applying the condition that at the phase transition point, the derivative of free energy with respect to the order parameter is zero ( $\partial F / \partial \Phi = 0$ ).
Derivation Steps:
1. Start with the incompressible Navier-Stokes equations.
2. Apply the divergence theorem to relate stress correlations to free energy equilibrium.
3. Formulate implicit functions for Reynolds stresses based on the symmetry of free energy in CPT.
4. Derive specific algebraic expressions for normal stresses ( $\overline{u'u'}$ , $\overline{v'v'}$ ) and shear stress ( $\overline{u'v'}$ ) as functions of normalized velocity.
5. Key Innovation: Unlike the scalar eddy viscosity in Boussinesq's hypothesis, the new model treats turbulent viscosity as a tensor ( $\nu_{T,ij}$ ), allowing for anisotropic stress modeling.

B. Numerical Implementation

Test Case: A plane turbulent jet (free shear flow) is selected as the benchmark.
Governing Equations: The derived closure equations (Eq. 18) are coupled with the RANS equations (Continuity and Momentum).
Numerical Scheme:
- Method: Finite Volume Method (FVM) on a structured rectangular grid.
- Algorithm: SIMPLE (Semi-Implicit Method for Pressure Linked Equations) for pressure-velocity coupling.
- Discretization: Second-order upwind for convective terms and second-order central difference for diffusive terms.
- Solver: Iterative line-by-line Tridiagonal Matrix Algorithm (TDMA).
Parameters: Simulations were conducted for an air jet with Reynolds number $Re = 3.2 \times 10^4$ , orifice height $h=0.04$ m, and exit velocity $u_o = 12$ m/s.

3. Key Contributions

New Closure Equations: Derivation of closed-form RANS equations where Reynolds stresses are expressed as functions of the mean velocity and its derivatives, rooted in the thermodynamics of phase transitions.
Tensorial Viscosity: Moving beyond the scalar eddy viscosity assumption, the model utilizes a turbulent viscosity tensor, better capturing the directional nature of turbulent stresses.
Symmetry Manifestation: The study demonstrates that turbulent stresses exhibit odd and even symmetries relative to the normalized mean velocity, directly reflecting the free energy symmetry of continuous phase transitions.
- Shear stress is modeled as an odd function (polynomial of degree 3 or 4).
- Normal stresses are modeled as even functions (derived via integration of the shear stress function).
Analytical Constants: The paper provides specific values for scaling constants ( $C_{ij}$ , $C_s$ , $C_n$ , $\alpha$ ) determined by matching the theoretical model with established experimental and DNS data.

4. Results

The numerical solution of the new closure equations for a plane jet yielded the following findings:

Velocity Profiles:
- The mean axial velocity profiles ( $u/u_o$ ) and transverse velocity profiles ( $v/u_o$ ) showed excellent agreement with existing experimental data (Ramaprian & Chandrasekhara) and DNS data (Klein et al.).
- The jet exhibited self-similarity in the far-field region ( $x/h \geq 15$ ), collapsing well onto a Gaussian curve.
Jet Growth and Decay:
- The growth of the jet half-width ( $y_{1/2}$ ) and the decay of the centerline velocity ( $u_c$ ) followed linear trends consistent with theoretical predictions for free jets.
- The model successfully captured the effects of a small co-flow on jet growth rates.
Reynolds Stresses:
- The normalized normal stresses ( $\overline{u'u'}$ , $\overline{v'v'}$ ) and shear stress ( $\overline{u'v'}$ ) achieved self-similarity in the far field ( $x/h \geq 20$ ).
- The predicted stress profiles matched closely with DNS and experimental measurements, validating the new closure relations.
Stability: The numerical scheme remained stable using specific constants ( $b_1=3.1, b_2=1.7$ ) and under-relaxation factors, confirming the robustness of the derived equations.

5. Significance and Conclusion

Phenomenological Validation: The study provides strong evidence that the transition from laminar to turbulent flow can be effectively modeled as a continuous phase transition. The success of the closure equations supports the hypothesis that turbulence dynamics obey the same symmetry and scaling laws as thermodynamic phase transitions.
Improved Modeling: By treating turbulence as a phase transition and using a viscosity tensor, the model offers a more physically grounded alternative to empirical eddy viscosity models. It captures the intrinsic symmetries of the flow without relying solely on curve-fitting.
Future Implications: This approach opens a new avenue for turbulence modeling, potentially leading to more accurate predictions for complex flows where traditional RANS models struggle, by leveraging the mathematical framework of Landau's theory of phase transitions.

In summary, Azim successfully derived and validated a new class of turbulence closure equations that bridge fluid dynamics and statistical thermodynamics, demonstrating that treating turbulence as a continuous phase transition yields accurate and physically consistent results for plane jet flows.