ALBERT

Description: When air blows over water the wind exerts a stress at the interface thereby inducing in the water a sheared turbulent drift current. We present scaling arguments showing that, if a wind suddenly starts blowing, then the sheared drift current grows in depth on a timescale that is larger than the wave period, but smaller than a timescale for wave growth. This argument suggests that the drift current can influence growth of waves of wavelength λ that travel parallel to the wind at speed c. In narrow ‘inner’ regions either side of the interface, turbulence in the air and water flows is close to local equilibrium; whereas above and below, in ‘outer’ regions, the wave alters the turbulence through rapid distortion. The depth scale, la, of the inner region in the air flow increases with c/u*a (u*a is the unperturbed friction velocity in the wind). And so we classify the flow into different regimes according to the ratio la/λ. We show that different turbulence models are appropriate for the different flow regimes. When (u*a + c)/UB(λ) ≪ 1 (UB(z) is the unperturbed wind speed) la is much smaller than λ. In this limit, asymptotic solutions are constructed for the fully coupled turbulent flows in the air and water, thereby extending previous analyses of flow over irrotational water waves. The solutions show that, as in calculations of flow over irrotational waves, the air flow is asymmetrically displaced around the wave by a non-separated sheltering effect, which tends to make the waves grow. But coupling the air flow perturbations to the turbulent flow in the water reduces the growth rate of the waves by a factor of about two. This reduction is caused by two distinct mechanisms. Firstly, wave growth is inhibited because the turbulent water flow is also asymmetrically displaced around the wave by non-separated sheltering. According to our model, this first effect is numerically small, but much larger erroneous values can be obtained if the rapid-distortion mechanism is not accounted for in the outer region of the water flow. (For example, we show that if the mixing-length model is used in the outer region all waves decay!) Secondly, non-separated sheltering in the air flow (and hence the wave growth rate) is reduced by the additional perturbations needed to satisfy the boundary condition that shear stress is continuous across the interface. In a companion paper, we develop a numerical model for the coupled air-water flow with waves of arbitrary speed and in another we examine the detailed energy budget of the wave motions. © 1994, Cambridge University Press. All rights reserved.

Description: We investigate the changes to a fully developed turbulent boundary layer caused by the presence of a two-dimensional moving wave of wavelength L = 2π/k and amplitude a. Attention is focused on small slopes, ak, and small wave speeds, c, so that the linear perturbations are calculated as asymptotic sequences in the limit (U*+ c)/UB(L) →0 (U*is the unperturbed friction velocity and UB(L) is the approach-flow mean velocity at height L). The perturbations can then be described by an extension of the four-layer asymptotic structure developed by Hunt, Leibovich & Richards (1988) to calculate the changes to a boundary layer passing over a low hill. When (u*+ c)/UB(L) is small, the matched height, zm (the height where UBequals c), lies within an inner surface layer, where the perturbation Reynolds shear stress varies only slowly. Solutions across the matched height are then constructed by considering an equation for the shear stress. The importance of the shear-stress perturbation at the matched height implies that the inviscid theory of Miles (1957) is inappropriate in this parameter range. The perturbations above the inner surface layer are not directly influenced by the matched height and the region of reversed flow below zm: they are similar to the perturbations due to a static undulation, but the ‘effective roughness length’ that determines the shape of the unperturbed velocity profile is modified to zm = z0exp (KC/U*). The solutions for the perturbations to the boundary layer are used to calculate the growth rate of waves, which is determined at leading order by the asymmetric pressure perturbation induced by the thickening of the perturbed boundary layer on the leeside of the wave crest. At first order in (u*+ c)/UB(L), however, there are three new effects which, numerically, contribute significantly to the growth rate, namely: the asymmetries in both the normal and shear Reynolds stresses associated with the leeside thickening of the boundary layer, and asymmetric perturbations induced by the varying surface velocity associated with the fluid motion in the wave; further asymmetries induced by the variation in the surface roughness along the wave may also be important. © 1993, Cambridge University Press

Description: We investigate, using theoretical and computational techniques, the processes that lead to the drag force on a rigid surface that has two-dimensional undulations of length L and height H (with H/L 〈 1) caused by the flow of a turbulent boundary layer of thickness h. The recent asymptotic analyses of Sykes (1980) and Hunt, Leibovich & Richards (1988) of the linear changes induced in a turbulent boundary layer that flows over an undulating surface are extended in order to calculate the leading-order contribution to the drag. It is assumed that L is much less than the natural lengthscale h* = hUJu* over which the boundary layer evolves (w* is the unperturbed friction velocity and t/0a mean velocity scale in the approach flow). At leading order, the perturbation to the drag force caused by the undulations arises from a pressure asymmetry at the surface that is produced by the thickening of the perturbed boundary layer in the lee of the undulation. This we term non-separated sheltering to distinguish it from the mechanism proposed by Jeffreys (1925). Order of magnitude estimates are derived for the other mechanisms that contribute to the drag; the next largest is shown to be smaller than the non-separated sheltering effect by 0(w*/t/o). The theoretical value of the drag induced by the non-separated sheltering effect is in good agreement with both the values obtained by numerical integration of the nonlinear equations with a second-order-closure model and experiments. Although the analytical solution is developed using the mixing-length model for the Reynolds stresses, this model is used only in the inner region, where the perturbation shear stress has a significant effect on the mean flow. The analytical perturbation shear stresses are approximately equal to the results from a higher-order closure model, except where there is strong acceleration or deceleration. The asymptotic theory and the results obtained using the numerical model show that the perturbations to the Reynolds stresses in the outer region do not directly contribute a significant part of the drag. This explains why several previous analyses and computations that use the mixing-length model inappropriately throughout the flow lead to values of the drag force that are too large by up to 100 %. © 1993, Cambridge University Press

Description: An asymptotic analysis is developed for turbulent boundary layers in strong adverse pressure gradients. It is found that the boundary layer divides into three distinguishable regions: these are the wall layer, the wake layer and a transition layer. This structure has two key differences from the zero-pressure-gradient boundary layer: the wall layer is not eXponentially thinner than the wake; and the wake has a large velocity deficit, and cannot be linearized. The mean velocity profile has a y1/2 behaviour in the overlap layer between the wall and transition regions. The analysis is done in the conteXt of eddy viscosity closure modelling. It is found that k-ε-type models are suitable to the wall region, and have a power-law solution in the y1/2 layer. The outer-region scaling precludes the usual ε-equation. The Clauser, constant-viscosity model is used in that region. An asymptotic eXpansion of the mean flow and matching between the three regions is carried out in order to determine the relation between skin friction and pressure gradient. Numerical calculations are done for self-similar flow. It is found that the surface shear stress is a double-valued function of the pressure gradient in a small range of pressure gradients. © 1992, Cambridge University Press. All rights reserved.

Description: A body moves at uniform speed in an unbounded inviscid fluid. Initially, the body is infinitely far upstream of an infinite plane of marked fluid; later, the body moves through and distorts the plane and, finally, the body is infinitely far downstream of the marked plane. Darwin (1953) suggested that the volume between the initial and final positions of the surface of marked fluid (the drift volume) is equal to the volume of fluid associated with the ‘added-mass’ of the body.We re-examine Darwin's (1953) concept of drift and, as an illustration, we study flow around a sphere. Two lengthscales are introduced: ρmax, the radius of a circular plane of marked particles; and x0, the initial separation of the sphere and plane. Numerical solutions and asymptotic expansions are derived for the horizontal Lagrangian displacement of fluid elements. These calculations show that depending on its initial position, the Lagrangian displacement of a fluid element can be either positive – a Lagrangian drift – or negative – a Lagrangian reflux. By contrast, previous investigators have found only a positive horizontal Lagrangian displacement, because they only considered the case of infinite x0. For finite x0, the volume between the initial and final positions of the plane of marked fluid is defined to be the ‘partial drift volume’, which is calculated using a combination of the numerical solutions and the asymptotic expansions. Our analysis shows that in the limit corresponding to Darwin's study, namely that both x0 and ρmax become infinite, the partial drift volume is not well-defined: the ordering of the limit processes is important. This explains the difficulties Darwin and others noted in trying to prove his proposition as a mathematical theorem and indicates practical, as well as theoretical, criteria that must be satisfied for Darwin's result to hold.We generalize our results for a sphere by re-considering the general expressions for Lagrangian displacement and partial drift volume. It is shown that there are two contributions to the partial drift volume. The first contribution arises from a reflux of fluid and is related to the momentum of the flow; this part is spread over a large area. It is well-known that evaluating the momentum of an unbounded fluid is problematic since the integrals do not converge; it is this first term which prevented Darwin from proving his proposition as a theorem. The second contribution to the partial drift volume is related to the kinetic energy of the flow caused by the body: this part is Darwin's concept of drift and is localized near the centreline. Expressions for partial drift volume are generalized for flow around arbitrary-shaped two- and three-dimensional bodies. The partial drift volume is shown to depend on the solid angles the body subtends with the initial and final positions of the plane of marked fluid. This result explains why the proof of Darwin's proposition depends on the ratio ρmax/x0.An example of drift due to a sphere travelling at the centre of a square channel is used to illustrate the differences between drift in bounded and unbounded flows.