Turbulent Kinetic Energy Budget

Updated 29 January 2026

Turbulent kinetic energy budget is a framework that quantifies energy transfer, production, and dissipation in various turbulent flows.
It employs multi-scale and spectral decompositions to analyze cascade processes and diagnose interactions in wall-bounded, stratified, and multiphase environments.
The approach informs improved model closures and predictive tools for environmental, industrial, and plasma turbulence applications.

Turbulent kinetic energy (TKE) budgets provide an exact diagnostic of the transfer, production, and dissipation of kinetic energy in turbulent flows. They have central significance for theoretical, computational, and applied turbulence research. The foundations and detailed structures of TKE budgets vary across canonical flows, including wall-bounded turbulence, stably- and unstably-stratified environments, geophysical flows, multiphase systems, and kinetic plasmas. The following exposition covers exact budget forms, underlying mechanisms, multi-scale and spectral decompositions, model closures, and representative physical findings—specializing in technical detail for the research community.

1. Fundamental Formulation: The General TKE Budget

For an incompressible fluid, the instantaneous velocity field is decomposed into mean and fluctuation parts $u_i = U_i + u_i'$ , and the pressure is similarly decomposed. The Reynolds-averaged TKE per unit mass is $k = \frac{1}{2} \overline{u_i' u_i'}$ . The general TKE transport equation, neglecting body forces and density fluctuations, reads: $\frac{\partial k}{\partial t} + U_j \frac{\partial k}{\partial x_j} = -\overline{u_i' u_j'}\,\frac{\partial U_i}{\partial x_j} - \nu \overline{\frac{\partial u_i'}{\partial x_j} \frac{\partial u_i'}{\partial x_j}} - \frac{\partial}{\partial x_j} \left[ \frac{1}{2} \overline{u_i' u_i' u_j'} + \frac{1}{\rho} \overline{p' u_j'} - \nu \frac{\partial k}{\partial x_j} \right]$ where, in order, the right-hand terms are shear production ( $P$ ), viscous dissipation ( $\varepsilon$ ), and combined turbulent transport, pressure diffusion, and viscous diffusion. Extensions incorporate buoyancy effects, compressibility, stratification, multiphase physics, and other physics as needed (Wagner et al., 2020, Rogachevskii et al., 2021, Kleeorin et al., 2021, Kleeorin et al., 2018).

Classical Term Decomposition

Symbol	Physical Meaning	Canonical Expression
$P$	Shear production	$-\overline{u_i' u_j'} \frac{\partial U_i}{\partial x_j}$
$\varepsilon$	Dissipation (viscous)	$\nu \overline{\frac{\partial u_i'}{\partial x_j} \frac{\partial u_i'}{\partial x_j}}$
$T^T$	Turbulent transport (triple correlations)	$-\frac{1}{2} \frac{\partial}{\partial x_j} \overline{u_i' u_i' u_j'}$
$T^p$	Pressure diffusion	$-\frac{1}{\rho} \frac{\partial}{\partial x_j} \overline{p' u_j'}$
$D$	Viscous diffusion	$\nu \frac{\partial^2 k}{\partial x_j^2}$
$P_b$	Buoyancy production/consumption	$\beta \overline{u_z' \theta'}$ (Boussinesq; sign depends on stability)

The specific structure and importance of each term is domain-dependent.

2. Multi-Scale and Spectral Budgets

2.1 Scale-by-Scale Energy Transfer

The two-point or scale-filtered TKE budget generalizes the one-point budget by considering correlations between points separated by a vector $\mathbf{r}$ or by filtering at scale $\ell$ (Kirubel et al., 2024, Zhang et al., 2021, Thiesset et al., 15 Sep 2025, Adhikari et al., 16 Oct 2025).

Kármán–Howarth–Monin Equation for Homogeneous Isotropic Turbulence

For homogeneous isotropic flows, the Kármán–Howarth–Monin (KHM) equation governs two-point statistics: $0 = -\frac{1}{4r^4} \frac{d}{dr}\left[r^4 S_3(r)\right] + 2 \nu \frac{1}{r^4} \frac{d}{dr} \left[r^4 \frac{d F}{dr}\right] + I(r)$ where $S_3(r)$ is the third-order velocity structure function and $I(r)$ large-scale energy injection. In the inertial range, the dominant scale-by-scale nonlinear transfer (cascade) of energy ( $\Pi(r)$ ) equals the mean dissipation ( $\varepsilon$ ), independently of viscosity (Sy et al., 2020).

Scale-Space and Inhomogeneous Flows

For inhomogeneous and anisotropic flows, two-point frameworks include additional terms for interscale transport due to both mean-flow (SST $_m$ ) and fluctuating-field (SST $_f$ ) inhomogeneity, as well as explicit scale-space viscous diffusion and dissipation (Kirubel et al., 2024).

Multiphase and Compressible Extensions

In multiphase turbulence, surface tension mediates transfer across a pivotal scale ( $r_{\mathcal{S}}$ ), separating deformation and restoration mechanisms; pressure transport exhibits phase-dependent sign reversal (Thiesset et al., 15 Sep 2025, Perlekar, 2018). In compressible or kinetic plasma turbulence, scale-filtered pressure–strain interactions (notably the Pi–D shear component) dominate flow-to-internal-energy conversion at kinetic scales (Cao et al., 2020, Adhikari et al., 16 Oct 2025).

2.2 Spectral Decomposition

A spectral TKE budget expresses covariance and transfer as functions of wavenumber components: $\frac{\partial E_K(k_x, y, k_z)}{\partial t} = P + D + \Pi^s + W + T^\parallel - \varepsilon = 0$ Each term corresponds to the above physical processes, but as a function of scale and wall-normal position. Production occurs predominantly at large scales (outer VLSMs, streamwise streaks), while dissipation redistributes energy nearly isotropically at smaller scales. Nonlinear transfers (interscale and wall-normal) mediate cascades and energy redistribution among flow components (Lee et al., 2018).

3. Stratification, Buoyancy, and Nonlocal Forcing

3.1 Stratified and Wave-Driven Flows

In stratified environments (both stably and unstably stratified), buoyancy modifies the energy exchange between kinetic and potential reservoirs (Rogachevskii et al., 2021, Zilitinkevich et al., 2011, Kleeorin et al., 2021, Kirubel et al., 2024, Kleeorin et al., 2018). The EFB theory gives an augmented TKE budget,

$\frac{D E_K}{D t} + \partial_z \Phi_K = -\,\tau_{iz}\,\partial_z U_i + \beta F_z - \varepsilon_K$

and establishes parameterizations connecting Richardson numbers, turbulent Prandtl/Schmidt numbers, and dissipation timescales. Critically, even in very stable stratification, shear production ensures finite $E_K$ at large Ri—contrary to critical-Ri "shutoff" in older closures. Additional coupling to internal gravity waves provides a source for TKE under strongly damped, wave-driven conditions (Kleeorin et al., 2018).

3.2 Convectively Driven Flows

In strongly convective boundary layers and Rayleigh–Bénard convection, the governing TKE balance becomes

$\frac{\partial E_K}{\partial t} + \nabla \cdot \Phi_K = P_b + P_s - \varepsilon_K$

with buoyancy-dominated structures and pronounced anisotropy in the upper boundary. Multi-scale and filtered analyses (using horizontal filters or two-point statistics) have revealed that most TKE from buoyancy and large-scale rolls is transferred down to small scales, with only a small proportion dissipated directly at the largest scales (Green et al., 2019, Kirubel et al., 2024).

4. Specific Flows and Applications

4.1 Wall-Bounded Shear Flows

Oscillatory pipe flow: DNS demonstrates that production primarily occurs near the wall during deceleration; dissipation always peaks at the wall, and diffusive terms are generally secondary. The competition between production and dissipation across the Stokes layer controls turbulence onset, amplification, and decay. Scaling is set by the critical parameter $Re_\delta$ (Wagner et al., 2020).
Turbulent channel and boundary layer: The major TKE sink in the near-wall region is dissipation, which forms a structural "envelope" for both production and viscous diffusion; in outer regions, turbulent transport and pressure–velocity correlations dominate. Classical coherent structures (streaks, vortices, Q events) do not fully capture the mechanisms underpinning each TKE budget term; data-driven, explainable deep learning methods (SHAP analysis) offer more precise mappings (Alcántara-Ávila et al., 27 Jan 2026).

4.2 Environmental and Geophysical Flows

Oceanic mesoscale (West Coast): The surface TKE budget is dominated by wind-stress work and turbulent dissipation in narrow coastal zones (eddy-kine viscosity $\nu_e \sim 10^{-2}$ m $^2$ /s), with energy and vorticity generated nearshore then advected offshore, where local wind-work is negligible. Close balance between input and dissipation yields quasi-stationary mesoscale fields (Tóth et al., 2017).
Stratified atmospheric boundary layer: Double decomposition allows isolation of wake-added TKE behind wind turbines; shear production exceeds buoyant production in most conditions, and the wake TKE peak location is shifted downstream and enhanced in stable stratification due to nonlinear base–deficit interactions—core mechanisms underrepresented in empirical wake models (Klemmer et al., 2024).

4.3 Multiphase and Binary Fluid Turbulence

Multiphase (liquid–gas): The scale-by-scale TKE budget includes a surface-tension transfer term, which induces a two-pronged energy cascade. Above a pivotal scale $r_\mathcal{S}$ (close to the Hinze scale), kinetic energy is stored in interface deformation; below it, energy is re-injected by interface restoration. The total scale-wise transfer is generally larger than in single-phase turbulence and is nonmonotonic in phase fraction and density ratio (Thiesset et al., 15 Sep 2025). For phase-separating binary fluids, Cahn–Hilliard–Navier–Stokes simulations reveal a direct interfacial energy transfer channel that alters the classical inertial-range dynamics, reducing large-scale content and modifying cascade properties at the droplet scale (Perlekar, 2018).

4.4 Kinetic Plasmas

Collisionless kinetic scales: At sub-ion kinetic scales, the pressure–strain (particularly its shear Pi–D component) dominates flow-to-internal-energy conversion; the normal and shear contributions are strongly anti-correlated. The amplitude of these terms scales with the temperature ratio of electrons and ions, directly governing species-specific heating in space plasmas (Adhikari et al., 16 Oct 2025).

5. Model Closure and Parameterization

TKE budget closure is essential for practical modeling (LES, RANS) and is underpinned by physically and empirically grounded relationships:

EFB models employ algebraic or prognostic balances for all major TKE budget terms, including detailed dependence of turbulent Prandtl/Schmidt numbers and vertical anisotropy on stratification parameters (Zilitinkevich et al., 2011, Rogachevskii et al., 2021).
Scale-aware subgrid models for compressible turbulence must incorporate non-negligible pressure-dilatation, dilatational dissipation, and non-local transport, especially at high Mach number (Cao et al., 2020).
Wall-bounded turbulence models should enforce near-wall hierarchy, with dissipation envelopes subsuming production and diffusion at $y^+<10$ (Alcántara-Ávila et al., 27 Jan 2026).

6. Key Physical Insights, Scaling, and Implications

Hierarchy and spatial segregation: Near-wall turbulent energy dynamics are nested—dissipation structures "envelope" production and diffusion structures at small wall distances; away from boundaries and in outer layers, energy redistribution becomes dominated by transport and component exchange.
Multi-cascade and bottleneck effects: In multiphase and stratified turbulence, a two-pronged cascade (classical inertial plus surface tension or interfacial transfer) is manifest; phase interface and stratification impact not just the magnitude but the functional form and distribution of scale-by-scale kinetic energy transfer.
Self-similar scaling: In sufficiently extended domains and high Reynolds number, the spectral energy flows and cascades exhibit self-similar collapse with wall-normal coordinate or characteristic scale (e.g., $k_x y = \mathcal{O}(1)$ ), confirming the consistency with Townsend/Millikan scaling for outer/inner layer interactions (Lee et al., 2018).

7. Outlook and Open Directions

Modern approaches, including scale-resolved (DNS, spectral, and two-point) analyses, data-driven structure discovery, and detailed phase-conditioned budgets, continue to refine both physical understanding and closure strategies. Key challenges remain in quantifying scale-space cascades in inhomogeneous and multiphase environments, improving closure and subgrid representations for compressible, stratified, and kinetic flows, and incorporating these insights into operational predictive models across a range of scientific and engineering applications (Alcántara-Ávila et al., 27 Jan 2026, Kirubel et al., 2024, Thiesset et al., 15 Sep 2025, Adhikari et al., 16 Oct 2025).

Cited arXiv resources:

(Wagner et al., 2020, Alcántara-Ávila et al., 27 Jan 2026, Sy et al., 2020, Cao et al., 2020, Tóth et al., 2017, Kirubel et al., 2024, Rogachevskii et al., 2021, Thiesset et al., 15 Sep 2025, Kleeorin et al., 2018, Klemmer et al., 2024, Green et al., 2019, Adhikari et al., 16 Oct 2025, Kleeorin et al., 2021, Zhang et al., 2021, Perlekar, 2018, Zilitinkevich et al., 2011, Calvino et al., 2023, Lee et al., 2018)