The properties of the mean momentum balance in turbulent boundary layer, pipe and channel flows are explored both experimentally and theoretically. Available highquality data reveal a dynamically relevant four-layer description that is a departure from the mean profile four-layer description traditionally and nearly universally ascribed to turbulent wall flows. Each of the four layers is characterized by a predominance of two of the three terms in the governing equations, and thus the mean dynamics of these four layers are unambiguously defined. The inner normalized physical extent of three of the layers exhibits significant Reynolds-number dependence. The scaling properties of these layer thicknesses are determined. Particular significance is attached to the viscous/Reynolds-stress-gradient balance layer since its thickness defines a required length scale. Multiscale analysis (necessarily incomplete) substantiates the four-layer structure in developed turbulent channel flow. In particular, the analysis verifies the existence of at least one intermediate layer, with its own characteristic scaling, between the traditional inner and outer layers. Other information is obtained, such as (i) the widths (in order of magnitude) of the four layers, (ii) a flattening of the Reynolds stress profile near its maximum, and (iii) the asymptotic increase rate of the peak value of the Reynolds stress as the Reynolds number approaches infinity. Finally, on the basis of the experimental observation that the velocity increments over two of the four layers are unbounded with increasing Reynolds number and have the same order of magnitude, there is additional theoretical evidence (outside traditional arguments) for the asymptotically logarithmic character of the mean velocity profile in two of the layers; and (in order of magnitude) the mean velocity increments across each of the four layers are determined. All of these results follow from a systematic train of reasoning, using the averaged momentum balance equation together with other minimal assumptions, such as that the mean velocity increases monotonically from the wall.

The growth of the two-dimensional single-mode Rayleigh-Taylor instability (RTI) at low Atwood number (A=0.04) is investigated using Direct Numerical Simulations. The main result of the paper is that, at long times and sufficiently high Reynolds numbers, the bubble acceleration becomes stationary, indicating mean quadratic growth. This is contrary to the general belief that single-mode Rayleigh-Taylor instability reaches a constant bubble velocity at long times. At unity Schmidt number, the development of the instability is strongly influenced by the perturbation Reynolds number, defined as Rep≡λsqrt[Agλ/(1+A)]/ν. Thus, the instability undergoes different growth stages at low and high Rep. A new stage, chaotic development, was found at sufficiently high Rep values, after the reacceleration stage. During the chaotic stage, the instability experiences seemingly random acceleration and deceleration phases, as a result of complex vortical motions, with strong dependence on the initial perturbation shape (i.e., wavelength, amplitude, and diffusion thickness). Nevertheless, our results show that the mean acceleration of the bubble front becomes constant at late times, with little influence from the initial shape of the interface. As Rep is lowered to small values, the later instability stages, chaotic development, reacceleration, potential flow growth, and even the exponential growth described by linear stability theory, are subsequently no longer reached. Therefore, the results suggest a minimum Reynolds number and a minimum development time necessary to achieve all stages of single-mode RTI development, requirements which were not satisfied in the previous studies of single-mode RTI.

Elements of the first-principles-based theory of Wei et al. (J. Fluid Mech., vol. 522, 2005, p. 303), Fife et al. (Multiscale Model. Simul., vol. 4, 2005a, p. 936; J. Fluid Mech., vol. 532, 2005b, p. 165) and Fife, Klewicki & Wei (J. Discrete Continuous Dyn. Syst., vol. 24, 2009, p. 781) are clarified and their veracity tested relative to the properties of the logarithmic mean velocity profile. While the approach employed broadly reveals the mathematical structure admitted by the time averaged Navier–Stokes equations, results are primarily provided for fully developed pressure driven flow in a two-dimensional channel. The theory demonstrates that the appropriately simplified mean differential statement of Newton's second law formally admits a hierarchy of scaling layers, each having a distinct characteristic length. The theory also specifies that these characteristic lengths asymptotically scale with distance from the wall over a well-defined range of wall-normal positions, y. Numerical simulation data are shown to support these analytical findings in every measure explored. The mean velocity profile is shown to exhibit logarithmic dependence (exact or approximate) when the solution to the mean equation of motion exhibits (exact or approximate) self-similarity from layer to layer within the hierarchy. The condition of pure self-similarity corresponds to a constant leading coefficient in the logarithmic mean velocity equation. The theory predicts and clarifies why logarithmic behaviour is better approximated as the Reynolds number gets large. An exact equation for the leading coefficient (von Kármán coefficient κ) is tested against direct numerical simulation (DNS) data. Two methods for precisely estimating the leading coefficient over any selected range of y are presented. These methods reveal that the differences between the theory and simulation are essentially within the uncertainty level of the simulation. The von Kármán coefficient physically exists owing to an approximate self-similarity in the flux of turbulent force across an internal layer hierarchy. Mathematically, this self-similarity relates to the slope and curvature of the Reynolds stress profile, or equivalently the slope and curvature of the mean vorticity profile. The theory addresses how, why and under what conditions logarithmic dependence is approximated relative to the specific mechanisms contained within the mean statement of dynamics.

No abstract

Recent studies by the present authors have empirically and analytically explored the properties and scaling behaviours of the Reynolds averaged momentum equation as applied to wall-bounded flows. The results from these efforts have yielded new perspectives regarding mean flow structure and dynamics, and thus provide a context for describing flow physics. A physical model of the turbulent boundary layer is constructed such that it is consonant with the dynamical structure of the mean momentum balance, while embracing independent experimental results relating, for example, to the statistical properties of the vorticity field and the coherent motions known to exist. For comparison, the prevalent, well-established, physical model of the boundary layer is briefly reviewed. The differences and similarities between the present and the established models are clarified and their implications discussed.

This paper provides a multiscale analytical study of steady incompressible turbulent flow through a channel of either Couette or pressure-driven Poiseuille type. Mathematically, the paper's two most novel features are that (1) the analysis begins with an underdetermined singular perturbation problem, namely the Reynolds averaged mean momentum balance equation, and (2) it leads to the existence of an infinite number of length scales. (These two features are probably linked, but the linkage will not be pursued.) The paper develops a credible assumption of a mathematical nature which, when added to the initial underdetermined problem, results in a knowledge of almost the complete layer (scaling) structure of the mean velocity and Reynolds stress profiles. This structure in turn provides a lot of other important information about those profiles. The possibility of almostlogarithmic sections of the mean velocity profile is given special attention. The sense in which the length scales are asymptotically proportional to the distance from the wall is determined. Most traditional theoretical analyses of these wall-bounded flows are based ultimately on either the classical overlap hypothesis, mixing length concepts, or similarity arguments. The present paper avoids those approaches and their attendant assumptions. Empirical data are also not used, except that the Reynolds stress takes on positive values. Instead, reasonable criteria are proposed for recognizing scaling layers in the flow, and they are then used to determine the scaling structure and much more information.

No abstract

A clear and consistent framework for the analysis of the outer region of adverse-pressure-gradient turbulent boundary layers is established in this paper based on basic principles and theory, and the help of six adverse-pressure-gradient turbulent boundary layer databases and a zero-pressure-gradient one. Outer velocity and length scales for the mean velocity defect and the Reynolds stresses are discussed first. The conditions of validity of four velocity scales are determined in terms of the shape factor, since one scale is restricted to small velocity-defect boundary layers (the friction velocity $u_{\unicode[STIX]{x1D70F}}$), one to large-defect ones (the pressure-gradient velocity $U_{po}$), while the two others are proper scales for all velocity-defect conditions (the Zagarola–Smits velocity $U_{zs}$ and the mixing-layer-type velocity $U_{m}$). The turbulent boundary layer equations are then used to bring out, in a consistent manner and without assuming any self-similar behaviour, a set of non-dimensional parameters characterizing the outer region of turbulent boundary layers with arbitrary pressure gradients. In terms of a generic outer length scale $L_{o}$ and velocity scale $U_{o}$, these non-dimensional parameters are the pressure-gradient parameter $\unicode[STIX]{x1D6FD}_{o}=L_{o}/(\unicode[STIX]{x1D70C}U_{o}^{2})\,\text{d}p_{e}/\text{d}x$, the Reynolds number $Re_{o}=U_{o}L_{o}/\unicode[STIX]{x1D708}(U_{o}/U_{e})$ and the inertial parameter $\unicode[STIX]{x1D6FC}_{o}=U_{e}V_{e}/U_{o}^{2}$, where $U_{e}$ and $V_{e}$ are respectively the streamwise and wall-normal components of mean velocity at the boundary layer edge. These parameters have a clear physical meaning: they are ratios of the order of magnitude of forces, with the Reynolds shear stress gradient (apparent turbulent force) as the reference force – inertial to apparent turbulent forces for $\unicode[STIX]{x1D6FC}_{o}$, pressure to apparent turbulent forces for $\unicode[STIX]{x1D6FD}_{o}$ and apparent turbulent to viscous forces for $Re_{o}$. We discuss at length their significance and determine under what conditions they retain their meaning depending on the outer velocity scale that is considered. The seven boundary layer databases are analysed and compared using the established framework. An astonishing and impressive result is obtained: by choosing $U_{o}=U_{zs}$, the streamwise evolution of the three ratios of forces in the outer region can be accurately followed with $\unicode[STIX]{x1D6FD}_{zs}$, $\unicode[STIX]{x1D6FC}_{zs}$ and $Re_{zs}$ in all these flows – not just the order of magnitude of these ratios. This cannot be achieved with $u_{\unicode[STIX]{x1D70F}}$ and $U_{po}$, and only imperfectly with $U_{m}$. Consequently, $\unicode[STIX]{x1D6FD}_{zs}$, $\unicode[STIX]{x1D6FC}_{zs}$ and $Re_{zs}$ can be used to follow – in a global sense – the streamwise evolution of the streamwise mean momentum balance in the outer region.