Turbulent transport in an inclined jet in crossflow

The present study experimentally investigates a turbulent jet in crossflow relevant to film cooling applications. The jet is inclined at 30°, and its mean velocity is the same as the crossflow. Magnetic resonance imaging is used to obtain the full three-dimensional velocity and concentration fields, whereas Reynolds stresses are obtained along selected planes by Particle Image Velocimetry. The critical role of the counter-rotating vortex pair in the mixing process is apparent from both velocity and concentration fields. The jet entrainment is not significantly higher than in an axisymmetric jet without crossflow, because the proximity of the wall inhibits the turbulent transport. Reynolds shear stresses correlate with velocity and concentration gradients, consistent with the fundamental assumptions of simple turbulence models. However the eddy viscosity is strongly anisotropic and non-homogeneous, being especially low along the leeward side of the jet close to injection. Turbulent diffusion acts to decouple mean velocity and concentration fields, as demonstrated by the drop in concentration flux within the streamtube issued from the hole. Volume-averaged turbulent diffusivity is calculated using a mass–flux balance across the streamtube emanating from the jet hole, and it is found to vary slowly in the streamwise direction. The data are compared with Reynolds-Averaged Navier–Stokes simulations with standard k − ε closure and an optimal turbulent Schmidt number. The computations underestimate the strength of the counter-rotating vortex pair, due to an overestimated eddy viscosity. On the other hand the entrainment is increasingly underpredicted downstream of injection. To capture the correct macroscopic trends, eddy viscosity and eddy diffusivity should vary spatially in different ways. Therefore a constant turbulent Schmidt number formulation is inadequate for this flow.     

The Flow Structure Produced by Pulsed-jet Vortex Generators in a Turbulent Boundary Layer in an Adverse Pressure Gradient

The effect of pulsed jet vortex generators on the structure of an adverse pressure gradient turbulent boundary layer flow was investigated. Two geometrically optimised vortex generator configurations were used, co-rotating and counter-rotating. The duty cycle and pulse frequency were both varied and measurements of the skin friction (using hot films) and flow structure (using stereo PIV) were performed downstream of the actuators. The augmentation of the mean wall shear stress was found to be dependent on the net mass flow injected by the actuators. A quasi steady flow structure was found to develop far downstream of the injection location for the highest pulse frequency tested. The actuator near field flow structure was observed to respond very quickly to variations in the jet exit velocity.     

Large Eddy Simulation of Compressible Subsonic Turbulent Jet Starting From a Smooth Contraction Nozzle

Compressible subsonic turbulent starting jet with a relatively large Reynolds number of significant practical importance is investigated using large eddy simulation (LES), starting from a smooth contraction nozzle. The computational domain of truncated conical shape is determined through the comparison of the time-averaged numerical solution with the particle imaging velocimetry measurements for the steady jet. It is shown that the starting jet consists of a leading vortex ring followed by a quasi-steady jet, and the instantaneous velocity field exhibits contraction and expansion zones, corresponding to the high pressure (HP) and low pressure (LP) regions formed by the convecting vortex rings, and are related to the Kelvin-Helmholtz instability. The thin boundary layer inside the smooth contraction nozzle evolves into a shear layer at the nozzle exit and develops with the downstream penetration of the jet. Using λ ₂ criterion, the formation and evolution of the vortical structures are temporally visualized, illustrating distortion of vortex rings into lobed shapes prior to break-down. Rib-shape streamwise vortex filaments exist in the braid region between a pair of consecutive vortex rings due to secondary instabilities. Finally, formation and dynamics of hairpin vortices in the shear layer is identified.     

Large-Scale Streamwise Vortices in Turbulent Channel Flow Induced by Active Wall Actuations

Direct numerical simulations of turbulent flow in a plane channel using spanwise alternatively distributed strips (SADS) are performed to investigate the characteristics of large-scale streamwise vortices (LSSVs) induced by small-scale active wall actuations, and their role in suppressing flow separation. SADS control is obtained by alternatively applying out-of-phase control (OPC) and in-phase control (IPC) to the wall-normal velocity component of the lower channel wall, in the spanwise direction. Besides the non-controlled channel flow simulated as a reference, four controlled cases with 1, 2, 3 and 4 pairs of OPC/IPC strips are studied at M =?0.2 and R e =?6,000, based on the bulk velocity and the channel half height. The case with 2 pairs of strips, whose width is Δz ⁺ =?264 based on the friction velocity of the non-controlled case, is the most effective in terms of generating large-scale motions. It is also found that the OPC (resp. IPC) strips suppress (resp. enhance) the coherent structures and that leads to the creation of a vertical shear layer, which is responsible for the LSSVs presence. They are in a statistically steady state and their cores are located between two neighbouring OPC and IPC strips. These motions contribute significantly to the momentum transport in the wall-normal and spanwise directions showing potential for flow separation suppression.     

Effect of Multilateral Jet Mixing on Stability and Structure of Turbulent Partially-Premixed Flames

An experimental study was conducted to investigate the effects of multilateral jet mixing, using both three and four side-jets, on the structure and stability of turbulent partially-premixed flames. Particle Image Velocimetry and OH*-chemiluminescence were used to study the effects of geometry and operating conditions on the resulting flow-field and reaction zone structures, respectively. These effects were compared under varying ratios of side-jet to primary flow momentum, whilst keeping the bulk flow constant. It was found that the mixing regimes upstream of the nozzle exit affect the flame characteristics, i.e. an impinging regime is likely to generate a lifted flame whilst a backflow regime is likely to generate an attached flame. Unlike the 4 side-jets cases, the OH* images and v _{r m s} profiles for the 3 side-jets cases show distinct asymmetry, with intense OH* and low velocity fluctuations on the opposite sides of the fuel injection. It was also found that the flow and scalar fields become independent of the upstream conditions, for both 3 and 4 side-jets, after one diameter downstream of the nozzle exit.     

Predicting Turbulent Spectra in Drag-reduced Flows

In the present work we describe how turbulent skin-friction drag reduction obtained through near-wall turbulence manipulation modifies the spectral content of turbulent fluctuations and Reynolds shear stress with focus on the largest scales. Direct Numerical Simulations (DNS) of turbulent channels up to Re _τ = 1000 are performed in which drag reduction is achieved either via artificially removing wall-normal turbulent fluctuations in the vicinity of the wall or via streamwise-travelling waves of spanwise wall velocity. This near-wall turbulence manipulation is shown to modify turbulent spectra in a broad range of scales throughout the whole channel. Above the buffer layer, the observed changes can be predicted, exploiting the vertical shift of the logarithmic portion of the mean streamwise velocity profile, which is a classic performance measure for wall roughness or drag-reducing riblets. A simple model is developed for predicting the large-scale contribution to turbulent fluctuation and Reynolds shear stress spectra in drag-reduced turbulent channels in which a flow control acts at the wall. Any drag-reducing control that successfully interacts with large scales should deviate from the predictions of the present model, making it a useful benchmark for assessing the capability of a control to affect large scales directly.     

Turbulent shear stress profiles in a bubbly channel flow assessed by particle tracking velocimetry

Particle tracking velocimetry (PTV) is applied to a bubbly two-phase turbulent flow in a horizontal channel at Re = 2 × 10⁴ to investigate the turbulent shear stress profile which had been altered by the presence of bubbles. Streamwise and vertical velocity components of liquid phase are obtained using a shallow focus imaging method under backlight photography. The size of bubbles injected through a porous plate in the channel ranged from 0.3 to 1.5 mm diameter, and the bubbles show a significant backward slip velocity relative to liquid flow. After bubbles and tracer particles are identified by binarizing the image, velocity of each phase and void fraction are profiled in a downstream region. The turbulent shear stress, which consists of three components in the bubbly two-phase flow, is computed by analysis of PTV data. The result shows that the fluctuation correlation between local void fraction and vertical liquid velocity provides a negative shear stress component which promotes frictional drag reduction in the bubbly two-phase layer. The paper also deals with the source of the negative shear stress considering bubble’s relative motion to liquid.     

Measurement of local forcing on a turbulent boundary layer using PIV

An experimental study was carried out to investigate the effect of periodic blowing and suction on a turbulent boundary layer. Particle image velocimetry (PIV) was used to probe the characteristics of the flow. Local forcing was introduced to the boundary layer via a sinusoidally-oscillating jet issuing from a thin spanwise slot. Three forcing frequencies (f⁺=0.44, 0.66 and 0.88) with a fixed forcing amplitude (A⁺=0.6) were employed at Re _θ =690. The effect of three different forcing angles (α=60°, 90° and l20°) was investigated under a fixed forcing frequency (f⁺=0.088). The PIV results showed that the wall-region velocity decreases on imposition of the local forcing. Inspection of the phase-averaged velocity profiles revealed that spanwise large-scale vortices are generated downstream of the slot and persist farther downstream. The highest reduction in skin friction was achieved at the highest forcing frequency (f⁺=0.088) and a forcing angle of α=120°. The spatial fraction of the vortices was examined to analyze the skin friction reduction.     

Vortex-enhanced mixing through active and passive flow control methods

This study aims to understand the underlying physics of vortex-enhanced mixing through active and passive flow control methods. To find a best flow control method that enhances turbulent mixing through the generation of streamwise vortices, an experimental investigation was carried out to compare active and passive flow control methods of an incompressible axisymmetric jet. For active flow control, the lip of the circular jet was equipped with a single small flap deflected away from the jet stream at an angle of 30° to the jet axis. The flap incorporated a flow control slot through which steady and oscillatory suction were implemented. The active flow control methods require power input to the suction devices. For passive flow control, the lip of the circular jet was equipped with a single small delta tab deflected into the jet stream at an angle of 30° to the jet axis. The chord lengths of the flap and delta tab were one-sixth of the jet diameter. The momentum of jet increased in the case of active flow control by entraining the ambient fluid, whereas momentum decreased in the case of passive flow control. The effect of steady suction saturated for volumetric suction coefficient values greater than 0.82 %. The strength of streamwise vortices generated by the flap were greater than those generated by the delta tab. Steady suction produced positive pressures just downstream of the flow control slot in the central portion of the flap and negative pressures at the flap edges. Oscillatory suction was highly dependent on dimensionless frequency (F ⁺) based on the distance from the flow control slot to the flap trailing edge; the pressures on the central portion of the flap increased for F ⁺ ≤ 0.11 and then decreased for greater F ⁺; finally attained negative pressures at F ⁺ = 0.44. The increase in jet momentum and turbulence intensity, combined with the induced streamwise vorticity, makes steady suction a potential concept for increasing propulsion efficiency through vortex-enhanced mixing. The flow control methods modify the jet flow, which in turn would alter the jet noise spectra.     

Mean Flow and Turbulence Measurements in a Turbulent Free Cruciform Jet

The results of measurements of all three components of the mean velocity vector, the Reynolds normal and primary shear stresses and the mean static pressure in a turbulent free jet, issuing from a sharp-edged cruciform orifice, are presented in this paper. The measurements were made with an x-array hot-wire probe and a pitot-static tube in the near flow field of the jet. The Reynolds number, based upon the equivalent diameter of the orifice, was 1.70 × 10⁵. In addition to the quantities measured directly, the mean streamwise centreline velocity decay, the jet half-velocity widths, the jet spreading rate, the mean streamwise vorticity, the mass entrainment rate, the integral momentum flux and the one-dimensional energy spectra have been derived from the measured data. The results show that the mean streamwise centreline velocity decay rate of the cruciform jet is higher than that of a round jet issuing from an orifice with the same exit area as that of the cruciform orifice. The mean streamwise velocity field changed shape continuously from a cruciform close to the orifice exit plane to circular at 12 and half equivalent diameters downstream. The mean streamwise vorticity field, up to about three equivalent diameters downstream of the orifice exit plane, consists of four pairs of counter-rotating cells, which are aligned with the four edges in the centre of the cruciform orifice.     

Turbulent Scalar Mixing in a Skewed Jet in Crossflow: Experiments and Modeling

Turbulent mixing of an inclined, skewed jet injected into a crossflow is investigated using MRI-based experiments and a high-fidelity LES of the same configuration. The MRI technique provides three-dimensional fields of mean velocity and mean jet concentration. The 30° skew of the jet relative to the crossflow produces a single dominant vortex which introduces spanwise asymmetries to the velocity and concentration fields. The turbulent scalar transport of the skewed jet is investigated in further detail using the LES, which is validated against the experimental measurements. Mixing is found to be highly anisotropic throughout the jet region. Isotropic turbulent diffusivity and viscosity are used to calculate an optimal value of the turbulent Schmidt number, which varies widely over the jet region and lies mostly outside of the typically accepted range 0.7 ≤ S c _t ≤ 0.9. Finally, three common scalar flux models of increasing complexity are evaluated based on their ability to capture the anisotropy and predict the scalar concentration field of the present configuration. The higher order models are shown to better represent the turbulent scalar flux vector, leading to more accurate calculations of the concentration field. While more complex models are better able to capture the turbulent mixing, optimization of model constants is shown to significantly affect the results.     

Numerical Study of Turbulent Jet Ignition in a Lean Premixed Configuration

Direct numerical simulations (DNS) of a hot combustion product jet interacting with a lean premixed hydrogen-air coflow are conducted to fundamentally investigate turbulent jet ignition (TJI) in a three-dimensional configuration. TJI is an efficient method for initiating and controlling combustion in ultra-lean combustion systems. Fully compressible gas dynamics and species equations are solved with high order finite difference methods. The hydrogen-air reaction is simulated with a reliable detailed chemical kinetics mechanism. The physical processes involved in the TJI-assisted combustion are investigated by considering the flame heat release, temperature, species concentrations, vorticity, and Baroclinc torque. The complex turbulent flame and flow structures are delineated in three main: i) hot product jet, ii) burned-mixed, and iii) flame zones. In the TJI-assisted combustion, the flow structures and the flame features such as flame speed, temperature, and species distribution are found to be quite different than those in “standard” turbulent premixed combustion due to the existence of a high energy turbulent hot product jet. The flow structures and statistics are also found to be different than those normally seen in non-isothermal non-reacting jets.     

Helical Structures in the Near Field of a Turbulent Pipe Jet

We perform a finely resolved Large-eddy simulation to study coherent vortical structures populating the initial (near-nozzle) zone of a pipe jet at the Reynolds number of 5300. In contrast to ‘top-hat’ jets featured by Kelvin-Helmholtz rings with the non-dimensional frequency S t≈0.3?0.6, no high-frequency dominant mode is observed in the near field of a jet issuing from a fully-developed pipe flow. Instead, in shear layers we observe a relatively wide peak in the power spectrum within the low-frequency range (S t≈0.14) corresponding to the propagating helical waves entering with the pipe flow. This is confirmed by the Fourier transform with respect to the azimuthal angle and the Proper Orthogonal Decomposition complemented with the linear stability analysis revealing that this low-frequency motion is not connected to the Kelvin-Helmholtz instability. We demonstrate that the azimuthal wavenumbers m=1?5 contain the most of the turbulent kinetic energy and that a common form of an eigenmode is a helical vortex rotating around the axis of symmetry. Small and large timescales are identified corresponding to “fast” and “slow” rotating modes. While the “fast” modes correspond to background turbulence and stochastically switch from co- to counter-rotation, the “slow” modes are due to coherent helical structures which are long-lived and have low angular velocities, in agreement with the previously described spectral peak at low S t.     

Direct Numerical Simulation of Subsonic Round Turbulent Jet

Direct numerical simulation(DNS) of spatially developing round turbulent jet flow with Reynolds number 4,700 was carried out. Over 20 million grid points were used in this simulation. Fully compressible three-dimensional Navier–Stokes equations were solved. High order explicit spatial difference schemes and Runge–Kutta time integration scheme were used to calculate derivatives and time marching, respectively. Non-reflecting boundary conditions and exit zone techniques were adopted. Some refined computational grids were used in order to capture the smallest turbulent structures near the centerline of the jet. Low level disturbance were imposed on the jet inflow velocity to trigger the developing of turbulence. Turbulent statistics such as mean velocity, Reynolds stresses, third order velocity moments were obtained and compared with experimental data. One-dimensional velocity autospectra was also calculated. The inertial region where the spectra decays according to the k^− 5/3 was observed. The quantitative profiles of mean velocity and all of the third order velocity moments which were difficult to measure via experimental techniques were presented here in detail. The jet flow was proven to be close to fully self-similar around 19 jet diameters downstream of jet exit. The statistic data and revealed flow feature obtained in this paper can provide valuable reference for round turbulent jet research.     

Pressure-Gradient Turbulent Boundary Layers Developing Around a Wing Section

A direct numerical simulation database of the flow around a NACA4412 wing section at R e _c = 400,000 and 5^° angle of attack (Hosseini et al. Int. J. Heat Fluid Flow 61, 117–128, 2016), obtained with the spectral-element code Nek5000, is analyzed. The Clauser pressure-gradient parameter β ranges from ? 0 and 85 on the suction side, and from 0 to ? 0.25 on the pressure side of the wing. The maximum R e _?? and R e _τ values are around 2,800 and 373 on the suction side, respectively, whereas on the pressure side these values are 818 and 346. Comparisons between the suction side with zero-pressure-gradient turbulent boundary layer data show larger values of the shape factor and a lower skin friction, both connected with the fact that the adverse pressure gradient present on the suction side of the wing increases the wall-normal convection. The adverse-pressure-gradient boundary layer also exhibits a more prominent wake region, the development of an outer peak in the Reynolds-stress tensor components, and increased production and dissipation across the boundary layer. All these effects are connected with the fact that the large-scale motions of the flow become relatively more intense due to the adverse pressure gradient, as apparent from spanwise premultiplied power-spectral density maps. The emergence of an outer spectral peak is observed at β values of around 4 for λ _z ? 0.65δ ₉₉, closer to the wall than the spectral outer peak observed in zero-pressure-gradient turbulent boundary layers at higher R e _?? . The effect of the slight favorable pressure gradient present on the pressure side of the wing is opposite the one of the adverse pressure gradient, leading to less energetic outer-layer structures.     

Numerical Simulation of Two Phase Turbulent Flow of Nanofluids in Confined Slot Impinging Jet

In this article, a numerical investigation is performed on flow and heat transfer of confined impinging slot jet, with a mixture of water and Al₂O₃ nanoparticles as the working fluid. Two-dimensional turbulent flow is considered and a constant temperature is applied on the impingement surface. The k ? ω turbulence model is used for the turbulence computations. Two-phase mixture model is implemented to study such a flow field. The governing equations are solved using the finite volume method. In order to consider the effect of obstacle angle on temperature fields in the channel, the numerical simulations were performed for different obstacle angles of 0° ? 60°. Also different geometrical parameters, volume fractions and Reynolds numbers have been considered to study the behavior of the system in terms of stagnation point, average and local Nusselt number and stream function contours. The results showed that the intensity and size of the vortex structures depend on jet- impingement surface distance ratio (H/W) and volume fraction. The maximum Nusselt number occurs at the stagnation point with the highest values at about H/W = 1. Increasing obstacle angle, from 15° to 60°, enhances the heat transfer rate. It was also revealed that the minimum value of average Nusselt number occurs in higher H/W ratios with decreasing the channel length.     

Turbulent Energy Scale-Budget Equations for nearly Homogeneous Sheared Turbulence

For moderate Reynolds numbers, the isotropic relation between second-order and third-order moments for velocity increments (Kolmogorov's equation) is not respected, reflecting a non-negligible correlation between the scales responsible for the injection, transfer and dissipation of the turbulent energy. For (shearless) grid turbulence, there is only one dominant large-scale phenomenon, which is the non-stationarity of statistical moments resulting from the decay of energy downstream of the grid. In this case, the extension of Kolmogorov's analysis, as carried out by Danaila, Anselmet, Zhou and Antonia, J. Fluid Mech. 391, 1999 359-369) is quite straightforward. For shear flows, several large-scale phenomena generally coexist with similar amplitudes. This is particularly the case for wall-bounded flows, where turbulent diffusion and shear effects can present comparable amplitudes. The objective of this work is to quantify, in a fully developed turbulent channel flow and far from the wall, the influence of these two effects on the scale-by-scale energy budget equation. A generalized Kolmogorov equation is derived. Relatively good agreement between the new equation and hot-wire measurements is obtained in the outer region (40 < x ⁺ ₃ < 150) of the channel flow, for which the turbulent Reynolds number is R _λ≈ 36.     

An Experimental Study of the Rotational Effects on Separated Turbulent Flow During Stall Delay

Three-dimensional velocity fields were measured using tomographic particle image velocimetry (Tomo-PIV) on a model of the blade of a small-scale horizontal axis wind turbine (HAWT) to study the effects of rotation on separated turbulent flows during stall delay at a global tip speed ratio (TSR) of 3 and a Reynolds number of 4800. The flow fields on a static airfoil were also measured at a similar angle-of-attack (AOA) and Reynolds number for comparison. It was observed that the blade’s rotation in the streamwise direction significantly affected both the mean flow and the turbulence statistics over the suction surface. The mean velocity fields revealed that, different from the airfoil flow at large AOA, the recirculation region with reversed flow did not exist on the suction surface of the blade and the flow was rather attached. Mean spanwise flow from blade’s root to its tip was also generated by the rotation. The mean vorticity vector of the blade flow was found to be tilted in the rotational direction of the blade, as well as in the wall-normal direction. Of particular effects of the rotation on Reynolds stresses were the enhancement of 〈w ²〉 and the creation of strong 〈v w〉. The production of Reynolds stresses was also affected by blade’s rotation directly through the rotational production terms and indirectly by dramatically changing the fluctuating velocity fields. The distribution of enstrophy was observed to be modified by rotation, too.