Turbulent/non-turbulent Interface Research Articles

The effect of subgrid-scale models on the entrainment of a passive scalar in a turbulent planar jet

Classical large-eddy simulation (LES) modelling assumes that the passive subgrid-scale (SGS) models do not influence large-scale quantities, even though there is now ample evidence of this in many flows. In this work, direct numerical simulation (DNS) and large-eddy simulations of turbulent planar jets at Reynolds number ReH = 6000 including a passive scalar with Schmidt number Sc = 0.7 are used to study the effect of several SGS models on the flow integral quantities e.g. velocity and scalar jet spreading rates. The models analysed are theSmagorinsky, dynamic Smagorinsky, shear-improved Smagorinsky and the Vreman. Detailed analysis of the thin layer bounding the turbulent and non-turbulent regions – the so-called turbulent/non-turbulent interface (TNTI) – shows that this region raises new challenges for classical SGS models. The small scales are far from equilibrium and contain a high fraction of the total kinetic energy and scalar variance, but the situation is worse for the scalar than for the velocity field. Both a-priori and a-posteriori (LES) tests show that the dynamic Smagorinsky and shear-improved models give the best results because they are able to accurately capture the correct statistics of the velocity and passive scalar fluctuations near the TNTI. The results also suggest the existence of a critical resolution Δx, of the order of the Taylor scale λ, which is needed for the scalar field. Coarser passive scalar LES i.e. Δx ≥ λ results in dramatic changes in the integral quantities. This fact is explained by the dynamics of the small scales near the jet interface.

Characteristics of the turbulent/nonturbulent interface in boundary layers, jets and shear-free turbulence

The characteristics of turbulent/nonturbulent interfaces (TNTI) from boundary layers, jets and shear-free turbulence are compared using direct numerical simulations. The TNTI location is detected by assessing the volume of turbulent flow as function of the vorticity magnitude and is shown to be equivalent to other procedures using a scalar field. Vorticity maps show that the boundary layer contains a larger range of scales at the interface than in jets and shear-free turbulence where the change in vorticity characteristics across the TNTI is much more dramatic. The intermittency parameter shows that the extent of the intermittency region for jets and boundary layers is similar and is much bigger than in shear-free turbulence, and can be used to compute the vorticity threshold defining the TNTI location. The statistics of the vorticity jump across the TNTI exhibit the imprint of a large range of scales, from the Kolmogorov micro-scale to scales much bigger than the Taylor scale. Finally, it is shown that contrary to the classical view, the low-vorticity spots inside the jet are statistically similar to isotropic turbulence, suggesting that engulfing pockets simply do not exist in jets.

Multiscale analysis of fluxes at the turbulent/non-turbulent interface in high Reynolds number boundary layers

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.

The turbulence boundary of a temporal jet

AbstractWe examine the structure of the turbulence boundary of a temporal plane jet at$\mathit{Re}= 5000$using statistics conditioned on the enstrophy. The data is obtained by direct numerical simulation and threshold values span 24 orders of magnitude, ranging from essentially irrotational fluid outside the jet to fully turbulent fluid in the jet core. We use two independent estimators for the local entrainment velocity${v}_{n} $based on the enstrophy budget. The data show clear evidence for the existence of a viscous superlayer (VSL) that envelopes the turbulence. The VSL is a nearly one-dimensional layer with low surface curvature. We find that both its area and viscous transport velocity adjust to the imposed rate of entrainment so that the integral entrainment flux is independent of threshold, although low-Reynolds-number effects play a role for the case under consideration. This threshold independence is consistent with the inviscid nature of the integral rate of entrainment. A theoretical model of the VSL is developed that is in reasonably good agreement with the data and predicts that the contribution of viscous transport and dissipation to interface propagation have magnitude$2{v}_{n} $and$- {v}_{n} $, respectively. We further identify a turbulent core region (TC) and a buffer region (BR) connecting the VSL and the TC. The BR grows in time and inviscid enstrophy production is important in this region. The BR shows many similarities with the turbulent–non-turbulent interface (TNTI), although the TNTI seems to extend into the TC. The average distance between the TC and the VSL, i.e. the BR thickness is about 10 Kolmogorov length scales or half a Taylor length scale, indicating that intense turbulent flow regions and viscosity-dominated regions are in close proximity.

Lagrangian statistics across the turbulent-nonturbulent interface in a turbulent plane jet

Lagrangian statistics from millions of particles are used to study the turbulent entrainment mechanism in a direct numerical simulation of a turbulent plane jet at Re(λ) ≈ 110. The particles (tracers) are initially seeded at the irrotational region of the jet near the turbulent shear layer and are followed as they are drawn into the turbulent region across the turbulent-nonturbulent interface (TNTI), allowing the study of the enstrophy buildup and thereby characterizing the turbulent entrainment mechanism in the jet. The use of Lagrangian statistics following fluid particles gives a more correct description of the entrainment mechanism than in previous works since the statistics in relation to the TNTI position involve data from the trajectories of the entraining fluid particles. The Lagrangian statistics for the particles show the existence of a velocity jump and a characteristic vorticity jump (with a thickness which is one order of magnitude greater than the Kolmogorov microscale), in agreement with previous results using Eulerian statistics. The particles initially acquire enstrophy by viscous diffusion and later by enstrophy production, which becomes "active" only deep inside the turbulent region. Both enstrophy diffusion and production near the TNTI differ substantially from inside the turbulent region. Only about 1% of all particles find their way into pockets of irrotational flow engulfed into the turbulent shear layer region, indicating that "engulfment" is not significant for the present flow, indirectly suggesting that the entrainment is largely due to "nibbling" small-scale mechanisms acting along the entire TNTI surface. Probability density functions of particle positions suggests that the particles spend more time crossing the region near the TNTI than traveling inside the turbulent region, consistent with the particles moving tangent to the interface around the time they cross it.

Effects of mean shear on the local turbulent entrainment process

AbstractWe report on effects of mean shear on the turbulent entrainment process, focusing in particular on their relation to small-scale processes in the proximity of the turbulent/non-turbulent interface (TNTI). Three-dimensional particle tracking velocimetry (3D-PTV) measurements of an axisymmetric jet are compared to data from a direct numerical simulation (DNS) of a zero-mean-shear (ZMS) flow. First, conditional statistics relative to the interface position are investigated in a pseudo-Eulerian view (i.e. in a fixed frame relative to the interface position) and in a Lagrangian view. We find that in a pseudo-Eulerian frame of reference, both vorticity fluctuations and mean shear contribute to the vorticity jump at the boundary between irrotational and turbulent regions. In contrast, the Lagrangian evolution of enstrophy along trajectories crossing the entrainment interface is almost exclusively dominated by vorticity fluctuations, at least during the first Kolmogorov time scales after passing the interface. A mapping between distance to the instantaneous interface versus conditional time along the trajectory shows that entraining particles remain initially close to the TNTI and therefore attain lower average enstrophy values. The ratio between the rate of change of enstrophy in the two frames of references defines the local entrainment velocity ${v}_{n} = - (\mathrm{D} {\omega }^{2} / \mathrm{D} t)/ (\partial {\omega }^{2} / \partial {\hat {x} }_{n} )$, where ${\omega }^{2} $ is enstrophy and ${\hat {x} }_{n} $ is the coordinate normal to the TNTI. The quantity ${v}_{n} $ is decomposed into mean and fluctuating components and it is found that mean shear enhances the local entrainment velocity via inviscid and viscous effects. Further, the analysis substantiates that for all investigated flow configurations the local entrainment velocity depends considerably on the geometrical shape of the interface. Depending on the surface shape, different small-scale mechanisms are dominant for the local entrainment process, i.e. viscous effects for convex shapes and vortex stretching for concave shapes, looking from the turbulent region towards the convoluted boundary. Moreover, turbulent fluctuations display a stronger dependence on the shape of the interface than mean shear effects.

Multiscale Geometry and Scaling of the Turbulent-Nonturbulent Interface in High Reynolds Number Boundary Layers

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.

Entrainment and mixing of water droplets across a shearless turbulent interface with and without gravitational effects

We describe experiments of the entrainment and mixing of water (sub-Kolmogorov scale) droplets across a turbulent–non-turbulent interface (TNI) as well a turbulent–turbulent interface (TTI) in shearless grid turbulence, over a time scale in which evaporation is insignificant. The flow is produced by means of a splitter plate with an active grid and water sprays on one side and screens or an active grid on the other side. The Taylor microscale Reλ on the turbulent side is 275 and the average dissipation scale Stokes number, Stη ≈ 0.2, and based on the integral scale, Stl ≈ 0.003. By changing the orientation of the grid system, gravitational effects may be excluded or included. We show that in the absence of gravity, for the Stokes number range studied (0.06 ≤ Stη ≤ 1.33), the droplet distribution does not change across the interface. With gravity, the larger drops are selectively mixed and this is more pronounced for the TNI than for the TTI. The particle concentration distribution is an error function for the TTI but departs significantly for the TNI due to the intermittency in the flow. In terms of particle concentration, the entrainment is most efficient for the TTI with gravity. The results are related to droplet entrainment in clouds.

Acceleration, pressure and related quantities in the proximity of the turbulent/non-turbulent interface

This paper presents an analysis of flow properties in the proximity of the turbulent/non-turbulent interface (TNTI), with particular focus on the acceleration of fluid particles, pressure and related small scale quantities such as enstrophy,ω2=ωiωi, and strain,s2=sijsij. The emphasis is on the qualitative differences between turbulent, intermediate and non-turbulent flow regions, emanating from the solenoidal nature of the turbulent region, the irrotational character of the non-turbulent region and the mixed nature of the intermediate region in between. The results are obtained from a particle tracking experiment and direct numerical simulations (DNS) of a temporally developing flow without mean shear. The analysis reveals that turbulence influences its neighbouring ambient flow in three different ways depending on the distance to the TNTI: (i) pressure has the longest range of influence into the ambient region and in the far region non-local effects dominate. This is felt on the level of velocity as irrotational fluctuations, on the level of acceleration as local change of velocity due to pressure gradients, Du/Dt≃ ∂u/∂t≃ −∇p/ρ, and, finally, on the level of strain due to pressure-Hessian/strain interaction, (D/Dt)(s2/2) ≃ (∂/∂t)(s2/2) ≃ −sijp,ij> 0; (ii) at intermediate distances convective terms (both for acceleration and strain) as well as strain production −sijsjkski> 0 start to set in. Comparison of the results at Taylor-based Reynolds numbersReλ= 50 andReλ= 110 suggests that the distances to the far or intermediate regions scale with the Taylor microscale λ or the Kolmogorov length scale η of the flow, rather than with an integral length scale; (iii) in the close proximity of the TNTI the velocity field loses its purely irrotational character as viscous effects lead to a sharp increase of enstrophy and enstrophy-related terms. Convective terms show a positive peak reflecting previous findings that in the laboratory frame of reference the interface moves locally with a velocity comparable to the fluid velocity fluctuations.

Experiments on the spreading of shear-free turbulence under the influence of confinement and rotation

From Lie-group (symmetry) analysis of the multi-point correlation equation Oberlack and Gunther (Fluid Dyn Res 33:453-476, 2003) found three different solutions for the behavior of shear-free turbulence: (i) a diffusion like solution, in which turbulence diffuses freely into the adjacent calm fluid, (ii) a deceleration wave like solution when there is an upper bound for the integral length scale and (iii) a finite domain solution for the case when rotation is applied to the system. This paper deals with the experimental validation of the theory. We use an oscillating grid to generate turbulence in a water tank and Particle Image Velocimetry (PIV) to determine the two-dimensional velocity and out-of-plane vorticity components. The whole setup is placed on a rotating table. After the forcing is initiated, a turbulent layer develops which is separated from the initially irrotational fluid by a sharp interface, the so-called turbulent/non-turbulent interface (TNTI). The turbulent region grows in time through entrainment of surrounding fluid. We measure the propagation of the TNTI and find quantitative agreement with the predicted spreading laws for case one and two. For case three (system rotation), we observe that there is a sharp transition between a 3D turbulent flow close to the source of energy and a more 2D-like wavy flow further away. We measure that the separation depth becomes constant and in this sense, we confirm the theoretical finite domain solution.

On turbulent entrainment and dissipation in dilute polymer solutions

We present a comparative experimental study of a turbulent flow developing in clear water and dilute polymer solutions (25 and 50 wppm polyethylene oxide). The flow is forced by a planar grid that oscillates vertically with stroke S and frequency f in a square container of initially still fluid. Two-component velocity fields are measured in a vertical plane passing through the center of the tank by using time resolved particle image velocimetry. After the forcing is initiated, a turbulent layer develops that is separated from the initially irrotational fluid by a sharp interface, the so-called turbulent/nonturbulent interface (TNTI). The turbulent region grows in time through entrainment of surrounding fluid until the fluid in the whole container is in turbulent motion. From the comparison of the experiments in clear water and polymer solutions we conclude: (i) Polymer additives modify the large scale shape of the TNTI. (ii) Both, in water and in the polymer solution the mean depth of the turbulent layer, H(t), follows the theoretical prediction for Newtonian fluids H(t)∝Kt, where K∝S2f is the “grid action.” (iii) We find a larger grid action for dilute polymer solutions than for water. As a consequence, the turbulent kinetic energy of the flow increases and the rate of energy input becomes higher. (iv) The entrainment rate β=ve/vrms (where ve=dH/dt is the interface propagation velocity and vrms is the root mean square of the vertical velocity) is lower for polymers (βp≈0.7) than for water (βw≈0.8). The measured values for β are in good agreement with similarity arguments, from which we estimate that in our experiment about 28% of the input energy is dissipated by polymers.

Effect of large-scale intermittency and mean shear on scaling-range exponents in a turbulent jet

The present study investigates the combined impact of the intermittency associated with the turbulent-nonturbulent interface and the mean shear rate in an axisymmetric jet on the structure of turbulence in the scaling range, where the spectrum exhibits a power-law behavior. Second-order structure functions, autocorrelations of the dissipation rate, and spectra of both the longitudinal velocity fluctuation and the passive temperature fluctuation are measured at a distance of 40 diameter downstream from the nozzle exit. All the scaling range exponents are influenced by the large-scale intermittency and the mean shear. The scalar fluctuation is much more sensitive to the variation in large-scale intermittency than the velocity fluctuation.

Conditional Sampling in Turbulence Measurement

(a) Shear flows which exhibit a turbulent-nonturbulent interface; (b) Shear layers perturbed by interaction with another field of turbu­ lence; (c) Quasi-periodic or periodic flows (these include flows behind pass­ ing blades in rotating machinery, boundary layers over a solid or liquid surface with stationary or progressive periodic waves, and more generally flows subjected to either internal or external peri­ odic perturbations); (d) Coherent structures in different shear flows (these structures are currently receiving close attention by the turbulence community in view of their importance to the flow dynamics).

Transition to turbulence in starting pipe flows.

An experimental investigation was performed to study the transition to turbulence in pipe flow started impulsively with a supercritical Reynolds number. The results show that a turbulentnon-turbulent interface propagates downstream by consuming the non-turbulent region where laminar boundary-layer flow is developing with time. For the startup of a completely quieted fluid, the propagation velocity is equal to the maximum velocity of the steady-state turbulent flow. For restarting the fluid flow after a brief shutoff, however, the propagation velocity is larger than the maximum velocity, and the difference increases with shorter quieting times. In the process of interface propagation, a wave of very low frequency appears ahead of the interface centring around the radial position where the smallest change in velocity is observed at the instant of propagation. A study of the mechanism of interface propagation shows that the final break of the abovementioned wave continuously triggers downstream propagation of the turbulent-flow region. In addition, random jumping of the interface in turbulence-decaying flow is concluded to be the main cause for the increase in propagation velocity for short quieting times.

Turbulent mixing in a two-dimensional jet

An ensemble of instantaneous temperature profiles across a two-dimensional heated jet was obtained by shooting a fine platinum resistance thermometer across the jet at speeds much higher than the local velocity of the jet. Each recording of the ensemble was shifted to have a common center, and ensemble averages were obtained counting only the profiles for which the spatial location in question was inside the jet. This conditional averaging produced mean and mean-square fluctuations which are quite flat across the jet, indicating that the fluid inside the jet is well-mixed. Visual inspection of instantaneous records showed no cold spots in the jet. Ensemble averages were also measured by aligning the edge or beginning of the instantaneous profiles. This led to a sharper rise in the mean temperature profile and an overshoot in the mean-square temperature fluctuations. This indicates that mixing across the turbulent-nonturbulent interface is stronger than in the main body of the jet.

Turbulent-nonturbulent interface detector

An analog circuit is presented for the detection in real time of the turbulent-nonturbulent interface at the free stream edge of turbulent shear flows. The output of the circuit is a square wave whose value is unity when the detector probe is in turbulent flow, and is zero when in nonturbulent flow. The effectiveness of several detection schemes is evaluated. A method for calibrating the detector circuit is given by using a ``synthetic,'' intermittently turbulent simulated signal of known intermittency properties. A circuit for generating the simulated signal is also presented.

Detection of the turbulent-nonturbulent interface in slightly heated turbulent shear flows

The detection of the interface between turbulent and nonturbulent zones in slightly heated turbulent shear flows on the basis of a temperature signal is discussed. Two detection parameters, a threshold level or gate δθg and a hold time τh, are used. A complete description of the detection process and of the procedure used to select nominal values for the detection param eters is presented. It is shown that the detection process is credible in the sense that the statistics of the intermittency function are nearly independent of the values of the detection parameters. Comparison of the intermittency function statistics with previous results is made.

Simple Phenomenological Theory of Turbulent Shear Flows

A rate equation is proposed to govern the variation of the effective turbulent viscosity. The effects of generation, convection, diffusion, and decay are each represented by appropriate terms leaving only two empirical constants to be determined by experiment. This rate equation together with the equations of motion form a closed system applicable to quasiparallel turbulent shear flows. For an incompressible turbulent boundary layer with zero pressure gradient, solutions were obtained by assuming local similarity and a linear growth of the boundary-layer thickness. Another problem, the turbulent-nonturbulent interface at the outer edge of the boundary layer was treated by using the further assumption that the large scale motion of the interface has no significant contribution to the Reynolds stress. It can be shown that for a nearly homogeneous domain, Prandtl's mixing length theory is a limiting case of the present theory.