The boundary layer equations for plane, incompressible, and steady flow are 

$$\matrix{ {u{{\partial u} \over {\partial x}} + v{{\partial u} \over {\partial y}} = - {1 \over \varrho }{{\partial p} \over {\partial x}} + v{{{\partial ^2}u} \over {\partial {y^2}}},} \cr {0 = {{\partial p} \over {\partial y}},} \cr {{{\partial u} \over {\partial x}} + {{\partial v} \over {\partial y}} = 0.} \cr }$$

Boundary Layer Theory

Coherent structures in wall turbulence transport momentum and provide a means of producing turbulent kinetic energy. Above the viscous wall layer, the hairpin vortex paradigm of Theodorsen coupled with the quasistreamwise vortex paradigm have gained considerable support from multidimensional visualization using particle image velocimetry and direct numerical simulation experiments. Hairpins can autogenerate to form packets that populate a significant fraction of the boundary layer, even at very high Reynolds numbers. The dynamics of packet formation and the ramifications of organization of coherent structures (hairpins or packets) into larger-scale structures are discussed. Evidence for a large-scale mechanism in the outer layer suggests that further organization of packets may occur on scales equal to and larger than the boundary layer thickness.

Hairpin vortex organization in wall turbulencea)

A survey is given of the major developments in three-dimensional velocity field measurements using the tomographic particle image velocimetry (PIV) technique. The appearance of tomo-PIV dates back seven years from the present review (Elsinga et al 2005a 6th Int. Symp. PIV (Pasadena, CA)) and this approach has rapidly spread as a versatile, robust and accurate technique to investigate three-dimensional flows (Arroyo and Hinsch 2008 Topics in Applied Physics vol 112 ed A Schroder and C E Willert (Berlin: Springer) pp 127–54) and turbulence physics in particular. A considerable number of applications have been achieved over a wide range of flow problems, which requires the current status and capabilities of tomographic PIV to be reviewed. The fundamental aspects of the technique are discussed beginning from hardware considerations for volume illumination, imaging systems, their configurations and system calibration. The data processing aspects are of uppermost importance: image pre-processing, 3D object reconstruction and particle motion analysis are presented with their fundamental aspects along with the most advanced approaches. Reconstruction and cross-correlation algorithms, attaining higher measurement precision, spatial resolution or higher computational efficiency, are also discussed. The exploitation of 3D and time-resolved (4D) tomographic PIV data includes the evaluation of flow field pressure on the basis of the flow governing equation. The discussion also covers a-posteriori error analysis techniques. The most relevant applications of tomo-PIV in fluid mechanics are surveyed, covering experiments in air and water flows. In measurements in flow regimes from low-speed to supersonic, most emphasis is given to the complex 3D organization of turbulent coherent structures.

Tomographic PIV: principles and practice

A model-based description of the scaling and radial location of turbulent fluctuations in turbulent pipe flow is presented and used to illuminate the scaling behaviour of the very large scale motions. The model is derived by treating the nonlinearity in the perturbation equation (involving the Reynolds stress) as an unknown forcing, yielding a linear relationship between the velocity field response and this nonlinearity. We do not assume small perturbations. We examine propagating helical velocity response modes that are harmonic in the wall-parallel directions and in time, permitting comparison of our results to experimental data. The steady component of the velocity field that varies only in the wall-normal direction is identified as the turbulent mean profile. A singular value decomposition of the resolvent identifies the forcing shape that will lead to the largest velocity response at a given wavenumber–frequency combination. The hypothesis that these forcing shapes lead to response modes that will be dominant in turbulent pipe flow is tested by using physical arguments to constrain the range of wavenumbers and frequencies to those actually observed in experiments. An investigation of the most amplified velocity response at a given wavenumber–frequency combination reveals critical-layer-like behaviour reminiscent of the neutrally stable solutions of the Orr–Sommerfeld equation in linearly unstable flow. Two distinct regions in the flow where the influence of viscosity becomes important can be identified, namely wall layers that scale with R+1/2 and critical layers where the propagation velocity is equal to the local mean velocity, one of which scales with R+2/3 in pipe flow. This framework appears to be consistent with several scaling results in wall turbulence and reveals a mechanism by which the effects of viscosity can extend well beyond the immediate vicinity of the wall. The model reproduces inner scaling of the small scales near the wall and an approach to outer scaling in the flow interior. We use our analysis to make a first prediction that the appropriate scaling velocity for the very large scale motions is the centreline velocity, and show that this is in agreement with experimental results. Lastly, we interpret the wall modes as the motion required to meet the wall boundary condition, identifying the interaction between the critical and wall modes as a potential origin for an interaction between the large and small scales that has been observed in recent literature as an amplitude modulation of the near-wall turbulence by the very large scales.

/pdf/a-critical-layer-framework-for-turbulent-pipe-flow-31bl89j0fy.pdf

A critical-layer framework for turbulent pipe flow

Pipe flow is a prominent example among the shear flows that undergo transition to turbulence without mediation by a linear instability of the laminar profile. Experiments on pipe flow, as well as plane Couette and plane Poiseuille flow, show that triggering turbulence depends sensitively on initial conditions, that between the laminar and the turbulent states there exists no intermediate state with simple spatial or temporal characteristics, and that turbulence is not persistent, i.e., it can decay again, if the observation time is long enough. All these features can consistently be explained on the assumption that the turbulent state corresponds to a chaotic saddle in state space. The goal of this review is to explain this concept, summarize the numerical and experimental evidence for pipe flow, and outline the consequences for related flows.

Turbulence transition in pipe flow

Finite amplitude coherent structures with a reflection symmetry in the spanwise direction of a parallel boundary layer flow are reported together with a preliminary analysis of their stability. The search for the solutions is based on the self-sustaining process originally described by Waleffe (Phys. Fluids 9, 883 (1997)). This requires adding a body force to the Navier-Stokes equations; to locate a relevant nonlinear solution it is necessary to perform a continuation in the nonlinear regime and parameter space in order to render the body force of vanishing amplitude. Some states computed display a spanwise spacing between streaks of the same length scale as turbulence flow structures observed in experiments (S.K. Robinson, Ann. Rev. Fluid Mech. 23, 601 (1991)), and are found to be situated within the buffer layer. The exact coherent structures are unstable to small amplitude perturbations and thus may be part of a set of unstable nonlinear states of possible use to describe the turbulent transition. The nonlinear solutions survive down to a displacement thickness Reynolds number Re
 * = 496 , displaying a 4-vortex structure and an amplitude of the streamwise root-mean-square velocity of 6% scaled with the free-stream velocity. At this Re* the exact coherent structure bifurcates supercritically and this is the point where the laminar Blasius flow starts to cohabit the phase space with alternative simple exact solutions of the Navier-Stokes equations.

Unstable flow structures in the Blasius boundary layer.

The flow in a square duct is considered. Finite amplitude approximate traveling wave solutions, obtained using the self-sustaining-process approach introduced by Waleffe [Phys. Fluids 9, 883 (1997)], are obtained at low to moderate Reynolds numbers and used as initial conditions in direct numerical simulations. The ensuing dynamics is analyzed in a suitably defined phase space. Only one among the traveling wave solutions found is capable of surviving for a long time, with the flow trajectory forming quasiregular loops in phase space. Eventually, also this trajectory escapes along the manifold of a chaotic saddle and relaminarization ensues.

http://www.dicat.unige.it/wedin/PhysFluid2007.pdf

Coherent flow states in a square duct

The nonlinear stability of the asymptotic suction boundary layer is studied numerically, searching for finite-amplitude solutions that bifurcate from the laminar flow state. By changing the boundary conditions for disturbances at the plate from the classical no-slip condition to more physically sound ones, the stability characteristics of the flow may change radically, both for the linearized as well as the nonlinear problem. The wall boundary condition takes into account the permeability K of the plate; for very low permeability, it is acceptable to impose the classical boundary condition (K=0). This leads to a Reynolds number of approximately Re(c)=54400 for the onset of linearly unstable waves, and close to Re(g)=3200 for the emergence of nonlinear solutions [F. A. Milinazzo and P. G. Saffman, J. Fluid Mech. 160, 281 (1985); J. H. M. Fransson, Ph.D. thesis, Royal Institute of Technology, KTH, Sweden, 2003]. However, for larger values of the plate's permeability, the lower limit for the existence of linear and nonlinear solutions shifts to significantly lower Reynolds numbers. For the largest permeability studied here, the limit values of the Reynolds numbers reduce down to Re(c)=796 and Re(g)=294. For all cases studied, the solutions bifurcate subcritically toward lower Re, and this leads to the conjecture that they may be involved in the very first stages of a transition scenario similar to the classical route of the Blasius boundary layer initiated by Tollmien-Schlichting (TS) waves. The stability of these nonlinear solutions is also investigated, showing a low-frequency main unstable mode whose growth rate decreases with increasing permeability and with the Reynolds number, following a power law Re(-ρ), where the value of ρ depends on the permeability coefficient K. The nonlinear dynamics of the flow in the vicinity of the computed finite-amplitude solutions is finally investigated by direct numerical simulations, providing a viable scenario for subcritical transition due to TS waves.

/pdf/effect-of-plate-permeability-on-nonlinear-stability-of-the-252ewe9e9w.pdf

Effect of plate permeability on nonlinear stability of the asymptotic suction boundary layer.

Saddle-node bifurcations of traveling waves in the asymptotic suction boundary layer flow

http://iopscience.iop.org/article/10.1088/0169-5983/48/6/061411/pdf

Håkan Wedin

Papers

Unstable flow structures in the Blasius boundary layer.

Coherent flow states in a square duct

Effect of plate permeability on nonlinear stability of the asymptotic suction boundary layer.

Saddle-node bifurcations of traveling waves in the asymptotic suction boundary layer flow

Permeability models affecting nonlinear stability in the asymptotic suction boundary layer: the Forchheimer versus the Darcy model