Quantifying wall turbulence via a symmetry approach. Part I. A Lie group theory

TL;DR

This paper introduces a symmetry-based Lie group approach to analytically derive mean velocity profiles in wall-bounded turbulence, providing new invariant solutions validated by DNS data.

Contribution

It develops a novel symmetry-based framework using Lie group analysis to derive analytic expressions for mean velocity profiles in wall turbulence, including a new core layer.

Findings

Derived invariant solutions for different flow layers

Validated solutions with DNS data

Unified multi-layer flow formula

Abstract

First principle based prediction of mean flow quantities of wall-bounded turbulent flows (channel, pipe, and turbulent boundary layer - TBL) is of great importance from both physics and engineering standpoints. Here (Part I), we present a symmetry-based approach which derives analytic expressions governing the mean velocity profile (MVP) from an innovative Lie-group analysis. The new approach begins by identifying a set of order functions (e.g. stress length, shear-induced eddy length), in analogy with the order parameter in Landau's mean-field theory, which aims at capturing symmetry aspects of the fluctuations (e.g. Reynolds stress, dissipation). The order functions are assumed to satisfy a dilation group invariance - representing the effects of the wall on fluctuations - which allows us to postulate three new kinds of invariant solutions of the mean momentum equation (MME), focusing on group invariants of the order functions (rather than those of the mean velocity as done in previous studies). The first - a power law solution - gives functional forms for the viscous sublayer, the buffer layer, the log-layer, and a newly identified central `core' (for channel and pipe, but non-existent for TBL). The second - a defect power law of form $1-r^{m}$ ($r$ being the distance from the center line) - describes the `bulk zone' (the region of balance between production and dissipation). The third - a relation between the group invariants of the stress length function and its first derivative - describes scaling transition between adjacent layers. A combination of these three expressions yields a multi-layer formula covering the entire flow domain, identifying three important parameters: scaling exponent, layer thickness, and transition sharpness. All three kinds of invariant solutions are validated, individually and in combination, by data from direct numerical simulations (DNS).