Revisiting the mixing length scaling in pressure-gradient turbulent boundary layers via a symmetry approach

This paper presents a unified, symmetry-based analytical model that accurately predicts full-profile mixing length, mean velocity, and Reynolds shear stress in equilibrium adverse pressure-gradient turbulent boundary layers, thereby clarifying the transition from logarithmic to half-power-law scaling and determining key layer thicknesses without ad hoc fitting.

The Gist

Imagine a river flowing smoothly over a flat rock. This is easy to predict: the water moves in neat, parallel layers. But what happens when that river hits a hill? The water has to fight against gravity (or in this case, a "pressure gradient") to keep moving. It gets turbulent, chaotic, and starts swirling in complex ways.

For a hundred years, scientists have tried to write a single rulebook to predict exactly how this chaotic water moves, especially when it's being pushed hard against a wall. This paper, written by Wei-Tao Bi from Peking University, offers a new, unified rulebook for a specific type of chaotic flow called an Adverse Pressure Gradient (APG) Turbulent Boundary Layer.

Here is the breakdown of what the paper does, using simple analogies:

1. The Problem: The "Mixing Length" Mystery

To understand turbulence, scientists use a concept called "Mixing Length." Think of this as the average distance a swirling eddy (a tiny whirlpool of water) travels before it bumps into another one and loses its energy.

• The Old Rule: A century ago, a scientist named Prandtl said, "The mixing length is just a straight line that gets longer the further you get from the wall." This worked great for calm rivers (Zero Pressure Gradient).
• The Problem: When the river hits a hill (Adverse Pressure Gradient), that straight-line rule breaks. The water behaves differently, and scientists have been arguing for decades about how to fix the rulebook. Some say the mixing length stays constant; others say it changes shape.

2. The Solution: A "Symmetry" Approach

Instead of just guessing numbers to fit the data, this author uses a symmetry approach.

• The Analogy: Imagine a piece of clay. If you squeeze it from the sides (pressure), it doesn't just get shorter; it bulges out in a specific, predictable way based on the laws of physics. The author argues that turbulence has a hidden "symmetry" or a specific shape it must take when squeezed by pressure.
• By finding this hidden shape, the author builds a mathematical model that describes the entire flow profile, from the wall all the way to the edge of the flow, without needing to patch it together with different rules for different parts.

3. The Key Discoveries

A. The "Critical Tipping Point" (The Beta Number)
The paper identifies a specific "tipping point" in the strength of the pressure pushing back on the flow.

• Below the Tipping Point: The flow still has a "logarithmic" zone (a region where the speed increases in a predictable, steady way).
• Above the Tipping Point: The pressure is so strong that the "logarithmic" zone gets squashed and disappears. The flow transitions to a new rule called the "Half-Power Law."
• The Finding: The author calculates this tipping point to be a specific number (around 6.2). If the pressure is stronger than this, the old rules stop working, and the new "Half-Power" rules take over.

B. The "Universal Constant" (The Karman Constant)
Scientists have long argued about a specific number (called the Kármán constant) that appears in the math of these flows. Some say it changes depending on the flow; others say it's always the same.

• The Paper's Claim: The author argues that this number is always the same (0.45) if you look at the entire flow profile correctly. The reason it seems to change in experiments is that scientists were only looking at a small slice of the flow. When you look at the whole picture, the number is invariant (unchanging).

C. The "Self-Adjusting" Layers
The model automatically figures out how thick the different layers of the flow are (like the sticky layer right next to the wall vs. the chaotic layer further out).

• The Analogy: Think of the flow as a multi-layer cake. As the pressure increases, the bottom layers (the sticky ones) get squished thinner, and the top layers (the wake) get bigger. The author's math calculates exactly how much they squish without needing to manually measure them for every single experiment.

4. How They Tested It

The author didn't just write equations; they tested them against a massive library of data.

• They compared their model against Direct Numerical Simulations (super-computer simulations of water molecules), Large Eddy Simulations, and real-world wind tunnel experiments.
• The data covered a huge range: from gentle flows to flows so strong they are about to stop moving entirely (separation).
• The Result: The model matched the data incredibly well across the board, predicting the speed of the wind/water and the swirling forces (Reynolds stress) with high accuracy.

5. Why This Matters (According to the Paper)

• It Unifies the Rules: It connects the calm flow rules with the chaotic, high-pressure flow rules into one single, smooth mathematical formula.
• It Solves the "Log Law" Debate: It explains why and when the famous "Log Law" breaks down under strong pressure, replacing it with the "Half-Power Law."
• It Removes Guesswork: Unlike previous models that required scientists to tweak numbers to fit specific experiments, this model only needs one small correction factor (based on the maximum stress) and then predicts everything else automatically.

Summary

In short, this paper says: "We found the hidden symmetry in how turbulent water behaves when pushed hard against a wall. We found the exact point where the old rules stop working and new rules take over. And we proved that a fundamental constant of nature remains the same, provided you look at the whole picture."

It's a new, unified map for navigating the most chaotic parts of fluid flow, validated by decades of data and super-computer simulations.

Technical Summary

1. Problem Statement

The mixing length ($\ell_m$), introduced by Prandtl a century ago, remains a cornerstone of turbulence modeling. While Prandtl's linear scaling hypothesis ($\ell_m = \kappa y$) works well for zero-pressure-gradient (ZPG) flows, its applicability to Adverse Pressure Gradient (APG) turbulent boundary layers (TBLs) is debated.

• The Conflict: In APG flows, the total shear stress (TSS) is no longer linear, and the logarithmic law of the wall begins to break down as the flow approaches separation.
• The Gap: Existing models either rely on empirical fitting parameters that vary with flow conditions or fail to provide a unified, physically consistent description of the full mixing-length profile (from the viscous sublayer to the wake region). Specifically, the transition from the logarithmic law to the half-power (Stratford) law under strong APG lacks a rigorous analytical derivation.
• Objective: To develop a symmetry-based, analytical model that predicts the full-profile mixing length in equilibrium APG TBLs without ad hoc fitting, unifying the inner layer, log region, transition zone, and wake region.

2. Methodology

The authors extend the Structural Ensemble Dynamics (SED) theory, which uses Lie-group symmetry analysis to describe wall turbulence, and couple it with a recently developed symmetry-based TSS model.

• Theoretical Framework:

  • Symmetry Approach: The model assumes that the spatial evolution of physical quantities (mixing length and TSS) obeys dilation symmetry, broken by the wall and the pressure gradient.
  • Crossover Scaling Ansatz: The authors use a specific functional form, $[1 + (y/y_0)^\zeta]^{\Delta/\zeta}$, to smoothly transition between different scaling laws (e.g., from linear to power-law) across layer boundaries.
  • Two-Layer TSS Model: They utilize a TSS model ($\tau^+$) that accounts for the shear stress overshoot induced by APG, characterized by a critical thickness $y_p^*$.

• Key Hypotheses & Derivations:

  1. Logarithmic Region: The mixing length scales as $\ell_m = \kappa y \sqrt{\tau^+}$ to preserve the log law for mean velocity, where $\tau^+$ is the dimensionless TSS. This introduces a $\sqrt{\tau^+}$ correction to Prandtl's hypothesis.
  2. Wake Region: The mixing length becomes constant ($\ell_m = \lambda \delta$). The parameter $\lambda$ is derived analytically as a function of the Clauser parameter ($\beta$), decreasing from the ZPG value ($\kappa/4$) to a lower asymptotic limit.
  3. Half-Power Law Transition: A critical Clauser parameter ($\beta_c \approx 6.2$) is identified. Above this value, the logarithmic region shrinks, and the flow transitions to a half-power law scaling ($\ell_m \propto \sqrt{y}$).
  4. Inner Layer: The thicknesses of the viscous sublayer ($y_s^+$) and buffer layer ($y_b^+$) are determined self-consistently by matching the model to Nickels' universal inner-layer scaling, dependent on the pressure gradient parameter $p^+$.

• Model Formulation:
  The final full-profile model (Equation 26) combines these elements into a single compact expression containing only one empirical parameter ($\sigma_p$ or $\tau_{max}$) to account for finite-Reynolds-number effects.

3. Key Contributions

• Unified Full-Profile Model: The paper provides the first analytical formulation that unifies the viscous sublayer, buffer layer, log-law region, half-power transition, and wake region for equilibrium APG TBLs.
• Identification of Critical $\beta_c$: The study identifies a critical Clauser parameter $\beta_c \approx 6.2$. Below this, the log law holds with a shrinking upper bound; above this, the log law breaks down progressively, transitioning to the half-power law.
• Invariant Kármán Constant: The model posits a global, intrinsic Kármán constant $\kappa = 0.45$ (consistent with SED theory), distinguishing it from the locally fitted $\kappa$ values that appear to vary with APG strength in traditional analyses.
• Analytical Wake Parameter: An explicit formula for the wake-region mixing-length parameter $\lambda(\beta)$ is derived, replacing empirical constants with a physics-based function of the pressure gradient.
• Self-Consistent Inner Layers: The model determines inner-layer thicknesses ($y_s^+, y_b^+$) without ad hoc fitting, linking them directly to the pressure gradient and Reynolds number.

4. Results and Validation

The model was validated against a comprehensive database of Direct Numerical Simulations (DNS), Large Eddy Simulations (LES), and experimental data, covering:

• Reynolds Numbers: $Re_\theta = 1250 \sim 50,980$.
• Pressure Gradients: $\beta = 0 \sim 39$ (from ZPG to near-separation).

Key Findings:

• Mixing Length Profiles: The model accurately predicts the full mixing-length profile across all regimes, capturing the multi-layer scaling behavior. It outperforms recent semi-empirical models (e.g., Ma et al.) which fail in high-Reynolds-number, strong-APG cases.
• Mean Velocity & Reynolds Shear Stress: By integrating the mixing length and TSS models, the authors accurately predict mean velocity and Reynolds shear stress profiles.
  • The model correctly captures the "overshoot" in shear stress and the shift of the peak location away from the wall under strong APG.
  • It successfully reproduces the transition from the log law to the half-power law for cases with $\beta > 6.2$.
• Robustness: The model maintains high fidelity across the entire parameter space without case-by-case tuning, unlike competing models that require flow-dependent parameters ($b, n$) to compensate for TSS formulation errors.
• Diagnostic Functions: Analysis of diagnostic functions ($y^+ du^+/dy^+$, $\sqrt{y^+} du^+/dy^+$, etc.) confirms that $\beta_c = 6.2$ is the physically correct threshold for the log-law breakdown, whereas lower values (e.g., $\beta_c = 2$) lead to premature transitions inconsistent with DNS data.

5. Significance

• Resolution of Long-standing Debates: The work resolves the century-old debate on mixing-length scaling in APG flows by providing a physically consistent mechanism for the log-law breakdown and the emergence of the half-power law.
• Universality of $\kappa$: It supports the hypothesis that the intrinsic Kármán constant is invariant ($\kappa = 0.45$) and that apparent variations in fitted $\kappa$ are artifacts of finite Reynolds number effects and the breakdown of the log layer, not a change in the fundamental constant.
• Engineering Application: The model offers a robust, parameter-free (except for one measurable peak stress) closure for Reynolds-Averaged Navier-Stokes (RANS) turbulence models, significantly improving predictions for aerodynamic flows involving strong pressure gradients and near-separation conditions.
• Future Direction: The symmetry-based framework is shown to be extensible to non-equilibrium flows, suggesting a path toward a unified theory for complex wall-bounded turbulence.

In summary, this paper presents a rigorous, symmetry-based analytical framework that successfully unifies the description of mixing length in pressure-gradient turbulent boundary layers, offering a significant advancement over empirical and semi-empirical approaches.