The disappearance of viscous and laminar wakes in complex flows

The singular effects of steady large-scale external strain on the viscous wake generated by a rigid body and the overall flow field are analysed. In an accelerating flow strained at a positive rate, the vorticity field is annihilated owing to positive and negative vorticity either side of the wake centreline diffusing into one another and the volume flux in the wake decreases with downwind distance. Since the wake disappears, the far-field flow changes from monopolar to dipolar. In this case, the force on the body is no longer proportional to the strength of the monopole, but is proportional to the strength of the far field dipole. These results are extended to the case of strained turbulent wakes and this is verified against experimental wind tunnel measurements of Keffer (1965) and Elliott & Townsend (1981) for positive and negative strains. The analysis demonstrates why the total force acting on a body may be estimated by adding the viscous drag and inviscid force due to the irrotational straining field. Applying the analysis to the wake region of a rigid body or a bubble shows that the wake volume flux decreases even in uniform flows owing to the local straining flow in the near-wake region. While the wake volume flux decreases by a small amount for the flow over streamline and bluff bodies, for the case of a clean bubble the decrease is so large as to render Betz's (1925) drag formula invalid. To show how these results may be applied to complex flows, the effects of a sequence of positive and negative strains on the wake are considered. The average wake width is much larger than in the absence of a strain field and this leads to diffusion of vorticity between wakes and the cancellation of vorticity. The latter mechanism leads to a net reduction in the volume flux deficit downstream which explains why in calculations of the flow through groups of moving or stationary bodies the wakes of upstream bodies may be ignored even though their drag and lift forces have a significant effect on the overall flow field

The formation of shear and density layers in stably stratified trubulent flows: linear processes

The initial evolution of the momentum and buoyancy fluxes in a freely decaying, stably stratified homogeneous turbulent flow with r.m.s. velocity u′0 and integral lengthscale l0 is calculated using a weakly inhomogeneous and unsteady form of the rapid distortion theory (RDT) in order to study the growth of small temporal and spatial perturbations in the large-scale mean stratification N(z, t) and mean velocity profile u(z, t) (here N is the local Brunt–Väisälä frequency and u is the local velocity of the horizontal mean flow) when the ratio of buoyancy forces to inertial forces is large, i.e. Nl0/u′0[dbl greater-than sign]1. The lengthscale L of the perturbations in the mean profiles of stratification and shear is assumed to be large compared to l0 and the presence of a uniform background mean shear can be taken into account in the model provided that the inertial shear forces are still weaker than the buoyancy forces, i.e. when the Richardson number Ri = (N/[partial partial differential]zu)2[dbl greater-than sign]1 at each height. When a mean shear perturbation is introduced initially with no uniform background mean shear and uniform stratification, the analysis shows that the perturbations in the mean flow profile grow on a timescale of order N-1. When the mean density profile is perturbed initially in the absence of a background mean shear, layers with significant density gradient fluctuations grow on a timescale of order N−10 (where N0 is the order of magnitude of the initial Brunt–Väisälä frequency) without any associated mean velocity gradients in the layers. These results are in good agreement with the direct numerical simulations performed by Galmiche et al. (2002) and are consistent with the earlier physically based conjectures made by Phillips (1972) and Posmentier (1977). The model also shows that when there is a background mean shear in combination with perturbations in the mean stratification, negative shear stresses develop which cause the mean velocity gradient to grow in the density layers. The linear analysis for short times indicates that the scale on which the mean perturbations grow fastest is of order u′0/N0, which is consistent with the experiments of Park et al. (1994). We conclude that linear mechanisms are widely involved in the formation of shear and density layers in stratified flows as is observed in some laboratory experiments and geophysical flows, but note that the layers are also significantly influenced by nonlinear and dissipative processes at large times

Structure of unsteady stably stratified turbulence with mean shear.

The statistics of unsteady turbulence with uniform stratification N (Brunt–Väisälä frequency) and shear α(=dU1/dx3) are analysed over the entire time range (00 and \it Ri>0.25 respectively, oscillatory momentum and positive and negative density fluxes develop. Above a critical value of \it Ri\scriptsize\it crit(∼0.3), their average values are persistently countergradient. This structural change in the turbulence is the primary mechanism whereby stable stratification reduces the fluxes and the production of variances. It is quite universal and differs from the energy and stability mechanisms of Richardson (1926) and Taylor (1931). The long-time asymptotics of the energy ratio ER(=\it PE/VKE) of the potential energy to the vertical kinetic energy generally decreases with \it Ri(≥0.25), reaching the smallest value of 3/2 when there is no shear (\it Ri→∞). For strong mean shear (\it Ri<0.25), RDT significantly overestimates ER since (as in unstratified shear flow) it underestimates the vertical kinetic energy VKE. The RDT results show that the asymptotic values of the energy ratio ER and the normalized vertical density flux are independent of the initial value of ER, in agreement with DNS. This independence of the initial condition occurs because the ratios of the contributions from the initial values PE0 and KE0 are the same for PE and VKE and can be explained by the linear processes. Stable stratification generates buoyancy oscillations in the direction of the energy propagation of the internal gravity wave and suppresses the generation of turbulence by mean shear. Because the shear distorts the wavenumber fluctuations, the low-wavenumber spectrum of the vertical kinetic energy has the general form E33(k)∝(αtk)−1, where (LXαt)−1≪k≪L−1X (LX: integral scale). The viscous decay is controlled by the shear, so that the components of larger streamwise wavenumber k1 decay faster. Then, combined with the spectrum distortion by the shear, the energy and the flux are increasingly dominated by the small-k1 components as time elapses. They oscillate at the buoyancy period π/N because even in a shear flow the components as k1→0 are weakly affected by the shear. The effects of stratification N and shear α at small scales are to reduce both VKE and PE. Even for the same \it Ri, larger N and α reduce the high-wavenumber components of VKE and PE. This supports the applicability of the linear assumption for large N and α. At large scales, the stratification and shear effects oppose each other, i.e. both VKE and PE decrease due to the stratification but they increase due to the shear. We conclude that certain of these unsteady results can be applied directly to estimate the properties of sheared turbulence in a statistically steady state, but others can only be applied qualitatively

The disappearance of laminar and turbulent wakes in complex flows

The singular effects of steady large-scale external strain on the viscous wake generated by a rigid body and the overall flow field are analysed. In an accelerating flow strained at a positive rate, the vorticity field is annihilated owing to positive and negative vorticity either side of the wake centreline diffusing into one another and the volume flux in the wake decreases with downwind distance. Since the wake disappears, the far-field flow changes from monopolar to dipolar. In this case, the force on the body is no longer proportional to the strength of the monopole, but is proportional to the strength of the far field dipole. These results are extended to the case of strained turbulent wakes and this is verified against experimental wind tunnel measurements of Keffer (1965) and Elliott & Townsend (1981) for positive and negative strains. The analysis demonstrates why the total force acting on a body may be estimated by adding the viscous drag and inviscid force due to the irrotational straining field.Applying the analysis to the wake region of a rigid body or a bubble shows that the wake volume flux decreases even in uniform flows owing to the local straining flow in the near-wake region. While the wake volume flux decreases by a small amount for the flow over streamline and bluff bodies, for the case of a clean bubble the decrease is so large as to render Betz's (1925) drag formula invalid.To show how these results may be applied to complex flows, the effects of a sequence of positive and negative strains on the wake are considered. The average wake width is much larger than in the absence of a strain field and this leads to diffusion of vorticity between wakes and the cancellation of vorticity. The latter mechanism leads to a net reduction in the volume flux deficit downstream which explains why in calculations of the flow through groups of moving or stationary bodies the wakes of upstream bodies may be ignored even though their drag and lift forces have a significant effect on the overall flow field

Adjustment of a turbulent boundary layer to a 'canopy' of roughness elements

A model is developed for the adjustment of the spatially averaged time-mean flow of a deep turbulent boundary layer over small roughness elements to a canopy of larger three-dimensional roughness elements. Scaling arguments identify three stages of the adjustment. First, the drag and the finite volumes of the canopy elements decelerate air parcels; the associated pressure gradient decelerates the flow within an impact region upwind of the canopy. Secondly, within an adjustment region of length of order Lc downwind of the leading edge of the canopy, the flow within the canopy decelerates substantially until it comes into a local balance between downward transport of momentum by turbulent stresses and removal of momentum by the drag of the canopy elements. The adjustment length, Lc, is proportional to (i) the reciprocal of the roughness density (defined to be the frontal area of canopy elements per unit floor area) and (ii) the drag coefficient of individual canopy elements. Further downstream, within a roughness-change region, the canopy is shown to affect the flow above as if it were a change in roughness length, leading to the development of an internal boundary layer. A quantitative model for the adjustment of the flow is developed by calculating analytically small perturbations to a logarithmic turbulent velocity profile induced by the drag due to a sparse canopy with L/Lc≪1, where L is the length of the canopy. These linearized solutions are then evaluated numerically with a nonlinear correction to account for the drag varying with the velocity. A further correction is derived to account for the finite volume of the canopy elements. The calculations are shown to agree with experimental measurements in a fine-scale vegetation canopy, when the drag is more important than the finite volume effects, and a canopy of coarse-scale cuboids, when the finite volume effects are of comparable importance to the drag in the impact region. An expression is derived showing how the effective roughness length of the canopy, \z0eff, is related to the drag in the canopy. The value of \z0eff varies smoothly with fetch through the adjustment region from the roughness length of the upstream surface to the equilibrium roughness length of the canopy. Hence, the analysis shows how to resolve the unphysical flow singularities obtained with previous models of flow over sudden changes in surface roughness

Stratified separated flow around a mountain with an inversion layer below the mountain top

This paper presents analytical and numerical results for separated stratified inviscid flow over and around an isolated mountain in the limit of small Fronde number. The vertical density profile consists of a lower strongly stratified layer whose depth is just less than that of the mountain. It is separated from a semi-infinite upper stably stratified layer by a thin, highly stable, inversion layer. The paper aims to provide, for this particular profile, a thorough analysis of the three-dimensional separated flow over a mountain top with strong stratification. The Froude numbers F and F, of the lower layer and the interface are small with F-1 << F << 1, but the upper-layer Froude number is arbitrary. The flow at each height in the lower layer is governed by the two-dimensional Euler equations and moves horizontally around the mountain. It is given by a modification of a previous model using Kirchhoff free-streamline theory for the separated flow region downstream of the mountain. The pressure variations associated with the lower-layer flow are of the same order as the dynamic head and induce significant displacements of the inversion layer. When the inversion is near the top of the mountain these deflections are of the same order as the height of the projecting part of the mountain top and combine with the flow over the mountain top to excite vertically propagating internal waves in the upper layer. The resultant pressure field, vertical stream surface displacements, and surface streamlines in the upper layer are described consistently in the hydrostatic limit. Many of the features of the upper flow, including the perturbations of the critical dividing streamlines, are similar to those in flows with uniform stable stratification at low Froude number. Comparisons are made with experiments and approximate models for these summit flows based on the assumption that the dividing streamlines have small vertical displacement

Micro structure and Lagrangian statistcs of the scalar field with a mean gradient in isotropic turbulence

This paper presents an analysis and numerical study of the relations between the small-scale velocity and scalar fields in fully developed isotropic turbulence with random forcing of the large scales and with an imposed constant mean scalar gradient. Simulations have been performed for a range of Reynolds numbers from Reλ = 22 to 130 and Schmidt numbers from Sc = 1/25 to 144. The simulations show that for all values of Sc [gt-or-equal, slanted] 0.1 steep scalar gradients are concentrated in intermittently distributed sheet-like structures with a thickness approximately equal to the Batchelor length scale η/Sc[fraction one-half] with η the Kolmogorov length scale. We observe that these sheets or cliffs are preferentially aligned perpendicular to the direction of the mean scalar gradient. Due to this preferential orientation of the cliffs the small-scale scalar field is anisotropic and this is an example of direct coupling between the large- and small-scale fluctuations in a turbulent field. The numerical simulations also show that the steep cliffs are formed by straining motions that compress the scalar field along the imposed mean scalar gradient in a very short time period, proportional to the Kolmogorov time scale. This is valid for the whole range of Sc. The generation of these concentration gradients is amplified by rotation of the scalar gradient in the direction of compressive strain. The combination of high strain rate and the alignment results in a large increase of the scalar gradient and therefore in a large scalar dissipation rate. These results of our numerical study are discussed in the context of experimental results (Warhaft 2000) and kinematic simulations (Holzer & Siggia 1994). The theoretical arguments developed here follow from earlier work of Batchelor & Townsend (1956), Betchov (1956) and Dresselhaus & Tabor (1991)

Effects of rotation and sloping terrain on fronts of density current fronts

The initial stage of the adjustment of a gravity current to the effects of rotation with angular velocity f/2 is analysed using a short time analysis where Coriolis forces are initiated in an inviscid von Kármán–Benjamin gravity current front at tF=0. It is shown how, on a time-scale of order 1/f, as a result of ageostrophic dynamics, the slope and front speed UF are much reduced from their initial values, while the transverse anticyclonic velocity parallel to the front increases from zero to O(NH0), where N=g′/H0−−−−−√ is the buoyancy frequency, and g′=gΔρ/ρ0 is the reduced acceleration due to gravity. Here ρ0 is the density and Δρ and H0 are the density difference and initial height of the current. Extending the steady-state theory to account for the effect of the slope σ on the bottom boundary shows that, without rotation, UF has a maximum value for σ=\upi/6, while with rotation, UF tends to zero on any slope. For the asymptotic stage when ftF≫1, the theory of unsteady waves on the current is reviewed using nonlinear shallow-water equations and the van der Pol averaging method. Their motions naturally split into a ‘balanced’ component satisfying the Margules geostrophic relation and an equally large ‘unbalanced’ component, in which there is horizontal divergence and ageostrophic vorticity. The latter is responsible for nonlinear oscillations in the current on a time scale f−1, which have been observed in the atmosphere and field experiments. Their magnitude is mainly determined by the initial potential energy in relation to that of the current and is proportional to the ratio \it Bu−−−−−√=LR/R0, where LR=NH0/f is the Rossby deformation radius and R0 is the initial radius. The effect of slope friction also prevents the formation of a steady front. From the analysis it is concluded that a weak mean radial flow must be driven by the ageostrophic oscillations, preventing the mean front speed UF from halting sharply at ftF∼1. Depending on the initial value of LR/R0, physical arguments show that UF decreases slowly in proportion to (ftF)−1/2, i.e. UF/UF0=F(ftF,\it Bu). Thus the front only tends to the geostrophic asymptotic state of zero radial velocity very slowly (i.e. as ftF→∞) for finite values of LR/R0. However, as LR/R0→0, it reaches this state when ftF∼1. This analysis of the overall nonlinear behaviour of the gravity current is consistent with two two-dimensional non-hydrostatic (Navier–Stokes) and axisymmetric hydrostatic (shallow-water) Eulerian numerical simulations of the varying form of the rotating gravity current. When the effect of surface friction is considered, it is found that the mean movement of the front is significantly slowed. Furthermore, the oscillations with angular frequency f and the slow growth of the radius, when ftF≥1, are consistent with recent experiments

The influence of the thermal diffusivity of the lower boundary on eddy motion in convection

The paper presents new concepts and results for the eddy structure of turbulent convection in a horizontal fluid layer of depth h which lies above a solid base with thickness hb. The fluid parameters are the kinematic viscosity ν, the thermal diffusivity κ, which is taken to be comparable with ν, the density ρ, the specific heat cp and the expansion parameter β. The thermal diffusivity of the solid is κb. The results are an extension of the more commonly studied cases, where a constant heat flux or constant temperature is applied at the interface between the fluid and the base. The buoyancy forces induce eddy motions with a typical velocity w∗∼(gβFθh)1/3 where ρcpFθ is the average heat flux and Fθ the covariance of the fluctuations of the temperature and of the vertical velocity. At moderate Reynolds numbers (Re=w∗h/ν), say less than about 103, an order-of-magnitude analysis shows that for the case of high diffusivity of the base (i.e. κb≫κ) elongated ‘plumes’ form at the surface and extend to the top of the fluid layer. When the base diffusivity is low (i.e. κb≤κ) the surface cools below the developing ‘plume’ and either the plume breaks up into elongated puffs or, if κb≪κ, horizontal pressure gradients form so that only small-scale puffs can form near the surface. At very high Reynolds numbers, approximately greater than 104, the surface boundary layer below each puff/plume is highly turbulent with a local logarithmic velocity and temperature profile. An approximate analysis indicates for this case that there is insufficient buoyancy flux from the base, irrespective of its diffusivity, to maintain plumes, because of the high turbulent heat transfer. So puffs dominate high-Reynolds-number thermal convection as numerical simulations and field experiments demonstrate. However, when the surface heat flux is uniform, for example as a result of radiant heat transfer or by forcing with a constant heat flux below a very thin conducting base, plumes are the dominant form of eddy motion, as is commonly observed. In the numerical solutions presented here, where Re∼3×102 and the slab thickness hb=h, it is shown that the spatial scales of eddy structures in the fluid layer close to the surface become significantly smaller as κb/κ is reduced from 100 to 0.1. At the same time in the core of the convective layer the change in the autocorrelation and spatial correlation function indicates that there is a transition from long-duration plumes into shorter-duration and smaller-length-scale elongated puffs. The simulations show that the largest temperature fluctuations near the surface occur when a constant heat flux is applied at the bottom of the fluid layer. The smallest temperature fluctuations are associated with the constant-temperature boundary condition. The finite base diffusivity cases lie in between these limits, with the largest fluctuations occurring when the thermal diffusivity of the base is small. The hypothesis introduced above has been tested qualitatively in a laboratory set-up when the effective diffusivity of the base was varied. The flow structure was observed as it changed from being characterized by nearly steady plumes, into unsteady plumes and finally into puffs when the thickness of the conducting base was first increased and then the diffusivity was decreased

Displacement of inviscid fluid by a sphere moving away from a wall

We develop a theoretical analysis of the displacement of inviscid fluid particles and material surfaces caused by the unsteady flow around a solid body that is moving away from a wall. The body starts at position hs from the wall, and the material surface is initially parallel to the wall and at distance hL from it. A volume of fluid Df+ is displaced away from the wall and a volume Df- towards the wall. Df+ and Df- are found to be sensitive to the ratio hL/hs. The results of our specific calculations for a sphere can be extended in general to other shapes of bodies. When the sphere moves perpendicular to the wall the fluid displacement and drift volume Df+ are calculated numerically by computing the flow around the sphere. These numerical results are compared with analytical expressions calculated by approximating the flow around the sphere as a dipole moving away from the wall. The two methods agree well because displacement is an integrated effect of the fluid flow and the largest contribution to displacement is produced when the sphere is more than two radii away from the wall, i.e. when the dipole approximation adequately describes the flow. Analytic expressions for fluid displacement are used to calculate Df+ when the sphere moves at an acute angle α away from the wall. In general the presence of the wall reduces the volume displaced forward and this effect is still significant when the sphere starts 100 radii from the wall. A sphere travelling perpendicular to the wall, α = 0, displaces forward a volume Df+(0) = 4πa3hL/33/2hS when the marked surface starts downstream, or behind the sphere, and displaces a volume Df+(0) [similar] 2πa3/3 forward when it is marked upstream or in front of the body. A sphere travelling at an acute angle away from the wall displaces a volume Df+(α) [similar] Df+(0) cos α forward when the surface starts downstream of the sphere. When the marked surface is initially upstream of the sphere, there are two separate regions displaced forward and a simple cosine dependence on α is not found. These results can all be generalized to calculate material surfaces when the sphere moves at variable speed, displacements no longer being expressed in terms of time, but in relation to the distance travelled by the sphere