Mechanical Engineering - Research Publications
Permanent URI for this collection
Search Results
1 - 9 of 9
-
Item
-
ItemMultiscale analysis of fluxes at the turbulent/non-turbulent interface in high Reynolds number boundary layers(AIP Publishing, 2014-01)Analysis of fluxes across the turbulent/non-turbulent interface (TNTI) of turbulent boundary layers is performed using data from two-dimensional particle image velocimetry (PIV) obtained at high Reynolds numbers. The interface is identified with an iso-surface of kinetic energy, and the rate of change of total kinetic energy (K) inside a control volume with the TNTI as a bounding surface is investigated. Features of the growth of the turbulent region into the non-turbulent region by molecular diffusion of K, viscous nibbling, are examined in detail, focussing on correlations between interface orientation, viscous stress tensor elements, and local fluid velocity. At the level of the ensemble (Reynolds) averaged Navier-Stokes equations (RANS), the total kinetic energy K is shown to evolve predominantly due to the turbulent advective fluxes occurring through an average surface which differs considerably from the local, corrugated, sharp interface. The analysis is generalized to a hierarchy of length-scales by spatial filtering of the data as used commonly in Large-Eddy-Simulation (LES) analysis. For the same overall entrainment rate of total kinetic energy, the theoretical analysis shows that the sum of resolved viscous and subgrid-scale advective flux must be independent of scale. Within the experimental limitations of the PIV data, the results agree with these trends, namely that as the filter scale increases, the viscous resolved fluxes decrease while the subgrid-scale advective fluxes increase and tend towards the RANS values at large filter sizes. However, a definitive conclusion can only be made with fully resolved three-dimensional data, over and beyond the large dynamic spatial range presented here. The qualitative trends from the measurement results provide evidence that large-scale transport due to the energy-containing eddies determines the overall rate of entrainment, while viscous effects at the smallest scales provide the physical mechanism ultimately responsible for entrainment. Data spanning over a decade in Reynolds number suggest that the fluxes (or the entrainment velocity) scale with the friction velocity (or equivalently the local turbulent fluctuating velocity), whereas Taylor microscale and boundary-layer thickness are the appropriate length scales at small and large filter sizes, respectively.
-
ItemSpatial averaging of velocity measurements in wall-bounded turbulence: single hot-wires(IOP Publishing Ltd, 2013-11)
-
ItemSpatial averaging of streamwise and spanwise velocity measurements in wall-bounded turbulence using ν- and x-probes(IOP PUBLISHING LTD, 2013-11)
-
ItemMultiscale Geometry and Scaling of the Turbulent-Nonturbulent Interface in High Reynolds Number Boundary Layers(AMER PHYSICAL SOC, 2013-07-24)The scaling and surface area properties of the wrinkled surface separating turbulent from nonturbulent regions in open shear flows are important to our understanding of entrainment mechanisms at the boundaries of turbulent flows. Particle image velocimetry data from high Reynolds number turbulent boundary layers covering three decades in scale are used to resolve the turbulent-nonturbulent interface experimentally and, for the first time, determine unambiguously whether such surfaces exhibit fractal scaling. Box counting of the interface intersection with the measurement plane exhibits power-law scaling, with an exponent between -1.3 and -1.4. A complementary analysis based on spatial filtering of the velocity fields also shows power-law behavior of the coarse-grained interface length as a function of filter width, with an exponent between -0.3 and -0.4. These results establish that the interface is fractal-like with a multiscale geometry and fractal dimension of Df≈2.3-2.4.
-
ItemScaling of the turbulent/non-turbulent interface in boundary layers(CAMBRIDGE UNIV PRESS, 2014-07)Abstract Scaling of the interface that demarcates a turbulent boundary layer from the non-turbulent free stream is sought using theoretical reasoning and experimental evidence in a zero-pressure-gradient boundary layer. The data-analysis, utilising particle image velocimetry (PIV) measurements at four different Reynolds numbers ($\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\delta u_{\tau }/\nu =1200\mbox{--}14\, 500$), indicates the presence of a viscosity dominated interface at all Reynolds numbers. It is found that the mean normal velocity across the interface and the tangential velocity jump scale with the skin-friction velocity$u_{\tau }$and are approximately$u_{\tau }/10$and$u_{\tau }$, respectively. The width of the superlayer is characterised by the local vorticity thickness$\delta _{\omega }$and scales with the viscous length scale$\nu /u_{\tau }$. An order of magnitude analysis of the tangential momentum balance within the superlayer suggests that the turbulent motions also scale with inner velocity and length scales$u_{\tau }$and$\nu /u_{\tau }$, respectively. The influence of the wall on the dynamics in the superlayer is considered via Townsend’s similarity hypothesis, which can be extended to account for the viscous influence at the turbulent/non-turbulent interface. Similar to a turbulent far-wake the turbulent motions in the superlayer are of the same order as the mean velocity deficit, which lends to a physical explanation for the emergence of the wake profile in the outer part of the boundary layer.
-
ItemThe turbulent/non-turbulent interface and entrainment in a boundary layer(CAMBRIDGE UNIV PRESS, 2014-03)Abstract The turbulent/non-turbulent interface in a zero-pressure-gradient turbulent boundary layer at high Reynolds number ($\mathit{Re}_\tau =14\, 500$) is examined using particle image velocimetry. An experimental set-up is utilized that employs multiple high-resolution cameras to capture a large field of view that extends $2\delta \times 1.1\delta $ in the streamwise/wall-normal plane with an unprecedented dynamic range. The interface is detected using a criteria of local turbulent kinetic energy and proves to be an effective method for boundary layers. The presence of a turbulent/non-turbulent superlayer is corroborated by the presence of a jump for the conditionally averaged streamwise velocity across the interface. The steep change in velocity is accompanied by a discontinuity in vorticity and a sharp rise in the Reynolds shear stress. The conditional statistics at the interface are in quantitative agreement with the superlayer equations outlined by Reynolds (J. Fluid Mech., vol. 54, 1972, pp. 481–488). Further analysis introduces the mass flux as a physically relevant parameter that provides a direct quantitative insight into the entrainment. Consistency of this approach is first established via the equality of mean entrainment calculations obtained using three different methods, namely, conditional, instantaneous and mean equations of motion. By means of ‘mass-flux spectra’ it is shown that the boundary-layer entrainment is characterized by two distinctive length scales which appear to be associated with a two-stage entrainment process and have a substantial scale separation.
-
ItemReynolds-number-dependent turbulent inertia and onset of log region in pipe flows(CAMBRIDGE UNIV PRESS, 2014-10)Abstract A detailed analysis of the ‘turbulent inertia’ (TI) term (the wall-normal gradient of the Reynolds shear stress,$\mathrm{d} \langle -uv\rangle /\mathrm{d} y $), in the axial mean momentum equation is presented for turbulent pipe flows at friction Reynolds numbers$\delta ^{+} \approx 500$, 1000 and 2000 using direct numerical simulation. Two different decompositions for TI are employed to further understand the mean structure of wall turbulence. In the first, the TI term is decomposed into the sum of two velocity–vorticity correlations ($\langle v \omega _z \rangle + \langle - w \omega _y \rangle $) and their co-spectra, which we interpret as an advective transport (vorticity dispersion) contribution and a change-of-scale effect (associated with the mechanism of vorticity stretching and reorientation). In the second decomposition, TI is equivalently represented as the wall-normal gradient of the Reynolds shear stress co-spectra, which serves to clarify the accelerative or decelerative effects associated with turbulent motions at different scales. The results show that the inner-normalised position,$y_m^{+}$, where the TI profile crosses zero, as well as the beginning of the logarithmic region of the wall turbulent flows (where the viscous force is leading order) move outwards in unison with increasing Reynolds number as$y_m^{+} \sim \sqrt{\delta ^{+}}$because the eddies located close to$y_m^{+}$are influenced by large-scale accelerating motions of the type$\langle - w \omega _y \rangle $related to the change-of-scale effect (due to vorticity stretching). These large-scale motions of$O(\delta ^{+})$gain a spectrum of larger length scales with increasing$\delta ^{+}$and are related to the emergence of a secondary peak in the$-uv$co-spectra. With increasing Reynolds number, the influence of the$O(\delta ^{+})$motions promotes viscosity to act over increasingly longer times, thereby increasing the$y^{+}$extent over which the mean viscous force retains leading order. Furthermore, the TI decompositions show that the$\langle v \omega _z \rangle $motions (advective transport and/or dispersion of vorticity) are the dominant mechanism in and above the log region, whereas$\langle - w \omega _y \rangle $motions (vorticity stretching and/or reorientation) are most significant below the log region. The motions associated with$\langle - w \omega _y \rangle $predominantly underlie accelerations, whereas$\langle v \omega _z \rangle $primarily contribute to decelerations. Finally, a description of the structure of wall turbulence deduced from the present analysis and our physical interpretation is presented, and is shown to be consistent with previous flow visualisation studies.
-
ItemSelf-similarity in the inertial region of wall turbulence(AMER PHYSICAL SOC, 2014-12-24)The inverse of the von Kármán constant κ is the leading coefficient in the equation describing the logarithmic mean velocity profile in wall bounded turbulent flows. Klewicki [J. Fluid Mech. 718, 596 (2013)] connects the asymptotic value of κ with an emerging condition of dynamic self-similarity on an interior inertial domain that contains a geometrically self-similar hierarchy of scaling layers. A number of properties associated with the asymptotic value of κ are revealed. This is accomplished using a framework that retains connection to invariance properties admitted by the mean statement of dynamics. The development leads toward, but terminates short of, analytically determining a value for κ. It is shown that if adjacent layers on the hierarchy (or their adjacent positions) adhere to the same self-similarity that is analytically shown to exist between any given layer and its position, then κ≡Φ(-2)=0.381966..., where Φ=(1+√5)/2 is the golden ratio. A number of measures, derived specifically from an analysis of the mean momentum equation, are subsequently used to empirically explore the veracity and implications of κ=Φ(-2). Consistent with the differential transformations underlying an invariant form admitted by the governing mean equation, it is demonstrated that the value of κ arises from two geometric features associated with the inertial turbulent motions responsible for momentum transport. One nominally pertains to the shape of the relevant motions as quantified by their area coverage in any given wall-parallel plane, and the other pertains to the changing size of these motions in the wall-normal direction. In accord with self-similar mean dynamics, these two features remain invariant across the inertial domain. Data from direct numerical simulations and higher Reynolds number experiments are presented and discussed relative to the self-similar geometric structure indicated by the analysis, and in particular the special form of self-similarity shown to correspond to κ=Φ(-2).