Towards a full-scale version of Yakhot's model of strong turbulence

This paper presents a parameter-free, closed-form extension of Yakhot's model that successfully describes second and third-order structure functions across the full range of turbulent scales, from the smallest dissipative regimes to the system scale, by integrating an empirically observed relation with Kolmogorov's four-fifth law.

The Gist

Imagine you are standing by a rushing river. You see big, rolling waves (the large scales) and tiny, frantic ripples (the small scales). For over 80 years, physicists have been trying to write a single "rulebook" that explains how the water moves from the giant waves down to the tiniest ripples. This is the problem of turbulence.

This paper by Christoph Renner is a significant step toward writing that rulebook. Here is the story of what he did, explained without the heavy math.

The Problem: A Broken Bridge

For a long time, scientists had two different rulebooks:

1. The "Big Wave" Book: This explained how the large, rolling waves behave.
2. The "Tiny Ripple" Book: This explained how the smallest, dissipating ripples behave.

The problem was that these two books didn't talk to each other. There was a "bridge" in the middle (called the inertial range) where the rules changed, and no one had a single equation that could smoothly walk you from the biggest waves to the tiniest ripples without tripping over a gap.

A scientist named Yakhot had a great model for the middle part, but it was missing two things:

• It didn't know how to stop at the very top (the big waves).
• It didn't know how to handle the very bottom (the tiny ripples where friction kills the motion).

The Solution: Building the Full Bridge

Renner's paper is about taking Yakhot's model and extending it to cover the entire river, from the source to the sea.

1. The "Magic Glue" (The New Discovery)

To connect the tiny ripples to the rest of the river, Renner looked at experimental data (measurements from a high-speed helium gas jet). He was looking for a pattern in how the "energy" of the water changes as you get smaller.

He found a surprising "secret handshake" between the data.

• Imagine you have a set of measuring cups of different sizes (representing different "orders" of turbulence).
• Renner discovered that if you look at how fast the water level changes in one cup, it is directly related to the size of the next bigger cup, multiplied by a specific "magic number" (a constant).
• This relationship held true from the tiniest ripples all the way up to the middle of the river. It was like finding a single rule that governed how the water behaved at almost every size.

2. The "Transition Zone" (The New Ruler)

In the old models, there was a sharp line between "friction" (where water stops moving) and "turbulence" (where it swirls). Renner introduced a new concept: a transition zone.

Think of it like a ramp instead of a cliff.

• Old view: You are driving on a highway (turbulence), and suddenly you hit a wall of mud (friction) and stop instantly.
• Renner's view: The highway gently slopes down into a muddy field. There is a specific point on that slope where the car starts to feel the mud.

Renner calculated exactly where this "slope" begins. He found that this point depends on the Reynolds number (a fancy way of saying "how fast and how big the river is"). The faster and bigger the river, the smaller this transition zone becomes, but it always exists.

3. The Final Result: One Equation to Rule Them All

By combining:

• Yakhot's original model for the middle,
• A recent fix for the big waves,
• And his new "magic glue" for the tiny ripples...

Renner created a single, closed-form equation.

• No free parameters: You don't need to guess numbers to make it fit. You just plug in the size of the river and how fast it's flowing.
• Perfect fit: When he tested this equation against real data, it matched perfectly from the tiniest ripples to the biggest waves.

Why This Matters

Think of turbulence like a complex puzzle. For decades, we had the corner pieces (big waves) and the center pieces (tiny ripples), but the middle was a mess.

Renner didn't just find a few more pieces; he found the instruction manual that tells you exactly how all the pieces fit together. He showed us that the transition from "swirling chaos" to "frictional stop" isn't a mystery; it's a predictable slope that can be calculated.

In a nutshell:
This paper takes a famous theory of turbulence, fixes the holes at the very top and very bottom, and connects them with a new, simple rule discovered in the data. The result is a complete, self-contained map of how fluids move, from the largest storms to the smallest eddies, without needing to guess any numbers along the way.

Technical Summary

1. Problem Statement

Turbulent fluid motion, particularly fully developed homogeneous and isotropic turbulence, remains a fundamental challenge in physics. While Yakhot's model (based on Castaing's thermodynamic analogy) successfully describes the transition from inertial to large scales and captures the behavior of odd-order structure functions, it has two critical limitations:

1. Lack of Dissipative Range Description: The original and previously extended Yakhot models do not account for viscous dissipation effects, which dominate at small scales (the dissipation range).
2. Incomplete Small-Scale Dynamics: The models fail to correctly describe the convergence of even-order structure functions to their large-scale limits and do not capture the specific scaling laws observed in the dissipation range (where structure functions scale as $r^n$).

The goal of this paper is to extend Yakhot's framework to cover the full scale of turbulence, from the smallest dissipative scales up to the system scale, by deriving closed-form expressions for second and third-order structure functions ($S_2$ and $S_3$) without free parameters.

2. Methodology

The research combines theoretical modeling with empirical analysis of experimental data from a cryogenic axisymmetric helium gas jet ($Re \approx 2.3 \times 10^5$).

• Theoretical Framework: The study builds upon Yakhot's probability density function (PDF) equation for velocity increments. It incorporates a recent large-scale extension (introducing parameters $C$ and $D$) and seeks to add a term to account for small-scale dissipation.
• Dimensionless Formulation: All equations are expressed in dimensionless variables: $r = l/L$ (separation normalized by system scale), $u = v/\sigma$ (velocity increment normalized by RMS), and $S_n(r)$.
• Hypothesis Testing:
  • Hypothesis 1 (H1): The function $d(r)$ (associated with the large-scale symmetry term) is independent of the order $n$.
  • Hypothesis 2 (H2): The function $d(r)$ retains its large-scale value of 1 for most scales but approaches zero rapidly enough to have no influence on small-scale dynamics.
• Residual Analysis: The authors analyzed the residuals ($R_n$) between the theoretical differential equations and experimental data. They sought a universal relationship between these residuals and higher-order structure functions.
• Closure via Kolmogorov's 4/5 Law: Since the system is under-determined for general $n$, the authors utilized Kolmogorov's four-fifths law (relating $S_3$ to the energy dissipation rate $\epsilon$ and the derivative of $S_2$) to close the system specifically for orders $n=2$ and $n=3$.

3. Key Contributions

A. Discovery of a Universal Scaling Relation (Eq. 22)

The most significant empirical finding is a relation between the spatial derivative of even-order structure functions and the next higher odd-order structure function. For even orders $n$, the residual $R_n(r)$ follows a power law:
$$-n \frac{R_n(r)}{S_{n+1}(r)} = \frac{\tau}{r^2}$$
where $\tau \approx 0.026$ is a constant derived from data. This relation holds from the dissipation range well into the inertial range.

B. Introduction of a Characteristic Length Scale ($\rho$)

The extension introduces a new length scale parameter, $\rho$, which marks the crossover point from the dissipative regime to the inertial regime.

• Definition: $\rho^2 = 3\tau / Re$.
• Physical Meaning: It represents the scale where the scaling behavior transitions from the viscous $r^2$ law to the inertial $r^{\zeta_2}$ law.
• Reynolds Number Dependence: The paper derives an expression linking $\rho$ to the Reynolds number ($Re$) and energy dissipation ($\epsilon$), showing $\rho \propto Re^{-3/4}$, similar to the Kolmogorov microscale but representing the transition point rather than the point of total dissipation dominance.

C. Derivation of Closed-Form Full-Scale Models

Using the new scaling relation and the 4/5 law, the author derives explicit, closed-form solutions for $S_2(r)$ and $S_3(r)$ that are valid across all scales.

• $S_2(r)$: The model combines the large-scale Yakhot solution with a small-scale term involving $\rho$. It correctly reproduces the $r^2$ behavior at $r \to 0$ and the $r^{\zeta_2}$ behavior in the inertial range.
• $S_3(r)$: Derived by combining the 4/5 law with the new $S_2$ model, ensuring consistency with the $r^3$ small-scale limit.

4. Results

• Agreement with Experimental Data: The derived full-scale models for $S_2$ and $S_3$ show excellent agreement with experimental data across the entire range of length scales (from the smallest dissipative scales to the system scale $L$).
• No Free Parameters: The models are fully determined by the large-scale characteristics of the flow (system scale $L$, Reynolds number $Re$, and dimensionless energy dissipation $\epsilon$). The parameter $\tau$ (or $\rho$) is not a free fitting parameter but is derived from the Reynolds number.
• Validation of Hypotheses:
  • The assumption that $d(r) \to 0$ at small scales (H2) is validated by the model's ability to reproduce the correct small-scale asymptotics ($S_2 \propto r^2$).
  • The scaling law (Eq. 22) holds for high Reynolds numbers ($Re \approx 2.3 \times 10^5$).
• Limitations at Low Reynolds Numbers: Preliminary analysis of lower Reynolds number datasets ($Re \approx 7 \times 10^3$) shows deviations from the simple $r^{-2}$ power law in the residual scaling, suggesting the model is an approximation that becomes more accurate at higher $Re$.
• Odd vs. Even Asymmetry: The study notes an asymmetry where even-order residuals collapse into a universal power law, while odd-order residuals do not follow the same simple pattern, suggesting the need for a more general framework involving longitudinal and transverse increments for a complete theory.

5. Significance

This work represents a substantial step toward a complete analytical theory of fully developed turbulence.

1. Bridging the Gap: It successfully unifies the description of turbulence across the dissipative, inertial, and integral scales within a single mathematical framework.
2. Predictive Power: By eliminating free parameters, the model becomes a predictive tool based solely on macroscopic flow properties ($Re$, $\epsilon$).
3. Refinement of Yakhot's Model: It resolves the long-standing issue of Yakhot's model failing at small scales, providing a rigorous extension that respects both Kolmogorov's 4/5 law and the viscous scaling limits.
4. New Scaling Insight: The identification of the $\rho$ scale and the specific $r^{-2}$ scaling of the residuals offers new physical insights into the transition mechanisms between turbulence regimes.

In conclusion, Renner presents a robust, parameter-free model that accurately describes the structure functions of second and third order from the smallest dissipative scales to the system scale, effectively extending Yakhot's model to a "full-scale" version.