How far does the influence of the free surface extend in turbulent open channel flow?

TL;DR

This study uses direct numerical simulations to analyze the multi-layer structure of turbulence near the free surface in open channel flow, revealing how surface influence extends through the entire channel height at high Reynolds numbers.

Contribution

It provides detailed scaling laws and identifies a four-layer structure that describes the influence of the free surface in turbulent open channel flow.

Findings

The free surface influence extends throughout the entire channel height at Re_τ ≥ 400.

A four-layer structure describes the turbulence near the free surface.

The Kolmogorov and viscous sublayers scale with specific length scales.

Abstract

Turbulent open channel flow is known to feature a multi-layer structure near the free surface. In the present work we employ direct numerical simulations considering Reynolds numbers up to $\mathrm{Re}_\tau=900$ and domain sizes large enough ($L_x=12 \pi h$, $L_z=4 \pi h$) to faithfully capture the effect of very-large-scale motions in order to test the proposed scaling laws and ultimately answer the question: How far does the influence of the free surface extend? In the region near the free surface, where fluctuation intensities of velocity and vorticity become highly anisotropic, we observe the previously documented triple-layer structure, consisting of a wall-normal velocity damping layer that scales with the channel height $h$, and two sublayers that scale with the near-surface viscous length scale $\ell_V=\mathrm{Re}_b^{-1/2}h$ and with the Kolmogorov length scale $\ell_K=\mathrm{Re}_b^{-3/4}h$, respectively. The Kolmogorov sublayer measures $\delta_K \approx 20 \ell_K$ and the layer, where the wall-normal turbulence intensity decreases linearly to zero near the free surface, scales with $\ell_V$ and the corresponding near-surface viscous sublayer measures $\delta_V \approx \ell_V$. Importantly, the streamwise turbulence intensity profile for $\mathrm{Re}_\tau \ge 400$ suggests that the influence of the free-slip boundary penetrates essentially all the way down to the solid wall through the appearance of enhanced very-large-scale motions ($\delta_{SIL}\approx h$). In contrast, the layer where the surface-normal turbulence intensity is damped to zero is restricted to the free surface ($\delta_{NVD}\approx 0.3h$). As a consequence, the partitioning of the surface-influenced region has to be expanded to a four-layer structure that spans the entire channel height $h$.