Direct Numerical Simulation of Transitional Noncircular Buoyant Reactive Jets

The near field dynamics of transitional buoyant reactive jets established on noncircular geometries, including a rectangular nozzle with an aspect ratio of 2:1 and a square nozzle with the same cross-sectional area, are investigated by three-dimensional spatial direct numerical simulations. Without applying external perturbations at the inflow boundary, large vortical structures develop naturally in the flow field due to buoyancy effects. Simulation results and analysis describe the details and clarify mechanisms of vortex dynamics of the noncircular buoyant reactive jets. The interaction between density gradients and gravity initiates the flow vorticity. Among the major vorticity transport terms, the gravitational term mainly promotes flow vorticity in the cross-streamwise direction. For the baroclinic torque, it can either create or destroy flow vorticity depending on the local flow structure. The vortex stretching term has different effects on the streamwise and cross-streamwise vorticity. Streamwise vorticity is mainly created by vortex stretching, while this term can either create or destroy cross-streamwise vorticity. Under the coupling effects of buoyancy and noncircular nozzle geometry, three-dimensional vortex interactions lead to the transitional behavior of the reactive jets. Simulations also show that the rectangular jet is more vortical than the square jet. The rectangular jet has a stronger tendency of transition to turbulence at the downstream due to the aspect ratio effect. Mean flow property calculations show that the rectangular buoyant reactive jet has a higher entrainment rate than its square counterpart. Received 13 December 2000 and accepted 24 July 2001     

Some features of a turbulent wing–body junction vortical flow

Laser-Doppler velocimeter measurements of a wing/body junction flow field made within a plane to the side of the wing/wall junction and perpendicular both to a 3:2 elliptical nose—NACA 0020 tail wing, and a flat wall are presented. Reynolds number of the approach boundary layer was, Re_θ = 5940, and free-stream air velocity was, U_ref = 27.5 m/s. A large vortical structure residing in the outer region redirects the low-turbulence free-stream flow to the vicinity of the wing/wall junction, resulting in thin boundary layers with velocity magnitudes higher than free-stream flow. Lateral pressure gradients result in a three-dimensional separation on the uplifting side of the vortex. Additionally, a high vorticity vortical structure with opposite sense to the outer-layer vortex forms beneath the outer-layer vortex. Normal and shear stresses increase to attain values an order of magnitude larger compared to values measured in a three-dimensional boundary layer just outside the junction vortex. Bimodal histograms of the w fluctuating velocity occur under the outer-layer vortex near the wall due to the time-dependent nature of the horseshoe vortex. In such a flow the shear-stress angle (SSA) highly lags the flow-gradient angle (FGA), and the turbulence diffusion is highly altered due to presence of vortical structures.     

Three-dimensional numerical study of jet-in-crossflow characteristics at low Reynolds number

A numerical study of a square jet in a cross flow is carried out at a Reynolds number of 100. The flow field and heat transfer characteristic downstream of the jet have been explored by solving three-dimensional unsteady Navier–Stokes equations and energy equation using higher order spatial and temporal discretization. The projection of vortical structure on a plane is seen to give the component of vortex normal to the plane. Four combinations of velocity profile namely (1) uniform crossflow and uniform jet, (2) laminar boundary layer crossflow and uniform jet, (3) uniform crossflow and parabolic jet profile, and (4) laminar boundary layer crossflow and parabolic jet are compared at same phase to see their effect on the flow field and heat transfer characteristic. All the four cases are seen to exhibit unsteadiness but the jet with parabolic profile is seen to show stronger unsteadiness. The instantaneous vortical structures of all the cases at the same phase show that the structures are more complex for the jet with parabolic velocity profile. The temperature field is seen to be correlated with the vortical structures. Comparison of the time averaged flow field reveals that the jet penetration is the highest for the jet having parabolic profile and boundary layer crossflow. The adiabatic effectiveness is observed to be more for the jet with uniform velocity profile and uniform crossflow and was least for the jet with parabolic velocity profile and boundary layer crossflow.     

Impingement of a low Reynolds number turbulent circular jet onto a flat plate at normal incidence

The velocity field of a circular water jet impinging onto a flat plate has been measured using particle image velocimetry, or PIV. The velocity field has been recorded at several instants in time, producing thousands of simultaneous two-dimensional velocity measurements for each realization. The instantaneous velocity, vorticity and rate-of-strain fields reveal the interaction of vortices near the impinging wall within the radial wall jet downstream from the stagnation point. An ensemble average of the instantaneous fields produces a mean velocity field of the jet flow, which reveals many of the processes leading to boundary layer separation and vortex breakaway within the wall jet. The PIV system extracts the velocity measurements using a two-dimensional autocorrelation method, and can obtain thousands of highly accurate velocity measurements within a few minutes. The structure found in these experiments may be similar to the ground level structure of atmospheric microburst phenomena.A version of this paper was presented at the 11th Symposium on Turbulence, University of Missouri-Rolla, 17–19 October 1988     

Experimental investigation of the interaction between a plane wall jet and a parallel offset jet

The dual-jet flow generated by a plane wall jet and a parallel offset jet at an offset ratio of d/w = 1.0 has been investigated using Particle Image Velocimetry (PIV). The particle images are captured, processed, and subsequently used to characterize the flow in terms of the 2D velocity and vorticity distributions. Statistical characteristics of the flow are obtained through ensemble averaging of 360 instantaneous velocity fields. Also presented is a time series of instantaneous flow fields to illustrate the dynamic interaction between the two jets. Results reveal that the near field of the flow is characterized by a periodic large-scale Karman-like vortex shedding similar to what would be expected in the wake of a bluff body. The existence of the Karman-like vortices results in periodic interactions between the two jets; in addition, these vortices produce noticeable impact on the jet outer layers, i.e., the free shear layer of the offset jet and the wall boundary layer of the wall jet. A schematic of vortex/shear layer interaction is proposed to illustrate the flow pattern.     

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-Eddy Simulations of Forced Three-Dimensional Impinging Jets

Large-eddy simulations of the flow field around twin three-dimensional impinging jets were carried out to simulate the near-ground hover configuration of a vertical takeoff and landing (VTOL) aircraft. Both the impinging jet and the upwash caused by the collision of the wall jets are modeled in this study. The evolution of the vortical structures in the impinging jet flow field, due to the introduction of axisymmetric and azimuthal perturbations at the jet exit, has been investigated. The vortical structures formed in the jet shear layer due to azimuthal forcing, show significant three-dimensional vortex stretching effects when compared to the structures formed during axisymmetric forcing. Breakdown of the large-scale structures into smaller vortices also occurs much earlier during azimuthal forcing. When compared to the upwash formed during axisymmetric forcing, the azimuthally perturbed jet forms an upwash that is less coherent and results in a weaker upload or lift-off force on the aircraft undersurface. Comparison with available experimental data indicates good agreement for the centerline velocity decay, the wall pressure variation and the phase speed of the vortical structures.     

Turbulent structure during transition to self-similarity in a round jet

 The developing turbulent region of a round jet was investigated using an improved implementation of digital particle image velocimetry (DPIV). The two-dimensional flow field in planes normal and parallel to the axial velocity was measured at locations between 15 and 30 diameters downstream, for two Reynolds numbers of 5500 and 16,000. The study consisted of instantaneous snapshots of the velocity and vorticity fields as well as measurements of velocity correlations up to third order. In this regime, the Reynolds number had a significant effect on both the instantaneous flow structure and the profiles of mean velocity across the jet. Coherent streamwise structures were present in the jet core for the lower Reynolds number. Additional structures whose evolution was governed by time scales two orders of magnitude larger than the convective scale inside the jet were observed in the entrainment field. The velocity correlations provided further support for the validity of DPIV turbulence measurements. The data was consistent with the equations of motion and momentum was conserved. DPIV measurements of turbulent kinetic energy components agreed with the hot-wire measurements of previous studies. Received: 27 November 1996/Accepted: 14 July 1997     

Instantaneous and time-averaged flow fields of multiple vortices in the tip region of a ducted propulsor

The instantaneous and time-averaged flow fields in the tip region of a ducted marine propulsor are examined. In this flow, a primary tip-leakage vortex interacts with a secondary, co-rotating trailing edge vortex and other co- and counter-rotating vorticity found in the blade wake. Planar particle imaging velocimetry (PIV) is used to examine the flow in a plane approximately perpendicular to the mean axis of the primary vortex. An identification procedure is used to characterize multiple regions of compact vorticity in the flow fields as series of Gaussian vortices. Significant differences are found between the vortex properties from the time-averaged flow fields and the average vortex properties identified in the instantaneous flow fields. Variability in the vortical flow field results from spatial wandering of the vortices, correlated fluctuations of the vortex strength and core size, and both correlated and uncorrelated fluctuations in the relative positions of the vortices. This variability leads to pseudo-turbulent velocity fluctuations. Corrections for some of this variability are performed on the instantaneous flow fields. The resulting processed flow fields reveal a significant increase in flow variability in a region relatively far downstream of the blade trailing edge, a phenomenon that is masked through the process of simple averaging. This increased flow variability is also accompanied by the inception of discrete vortex cavitation bubbles, which is an unexpected result, since the mean flow pressures in the region of inception are much higher than the vapor pressure of the liquid. This suggests that unresolved fine-scale vortex interactions and stretching may be occurring in the region of increased flow variability.     

PIV investigation of the shear layer vortices in the near wake of a circular cylinder

Particle image velocimetry measurements are performed in the near wake of a circular cylinder at a Reynolds number of 12,500. Attention is focused on the shear layer that develops just downstream of the separation point from the cylinder surface to investigate the possible existence of a preferred spatio-temporal organization in this flow region and the possible occurrence of the vortex pairing phenomenon. Eddy structures are identified in instantaneous velocity maps in order to investigate their spatial relationships. For that purpose a vortex extraction procedure is designed, based on the wavelet transform of instantaneous maps of the swirling strength. This algorithm allows not only the detection of the vortical structures from the instantaneous velocity fields, giving access to their instantaneous location, but also the estimation of their main characteristics such as their radius, intensity and convection velocity. The vortex population detected in the shear layer is found to be of small diameter compared to that of the von Kármán vortex and of rather high intensity, in agreement with the existence of a thin shear layer. The strong flapping motion of the shear layer and its complex spatial development is also confirmed. By employing conditional analysis of the computed data and their proper scaling, the surrounding of the detected vortex cores is investigated. A preferred spatial vortex separation is detected and is shown to vary with the longitudinal distance from the origin of the shear layer, in agreement with the qualitative behavior of a turbulent plane mixing layer. Evidence of the occurrence of the vortex pairing or amalgamation mechanisms in the shear layer is also demonstrated.     

Reynolds stress and vorticity in turbulent wall flows

Experimental measurements in a boundary layer and a large-eddy simulation of plane channel flow have been used to study the dynamics of vorticity and mass transport in the nearwall region. It was found that Reynolds stress generation occurs in the vicinity of quasi-streamwise vortices, and that smoke particles tend to be ejected from the wall near these vortical structures.     

Perturbed vortical layers and shear sheltering

New theoretical results and physical interpretations are presented concerning the interactions between different types of velocity fields that are separated by thin interfacial layers, where there are dynamically significant variations of vorticity across the layers and, in some cases within them. It is shown how, in different types of complex engineering and environmental flow, the strengths of these interactions vary from the weakest kind of superposition to those where they determine the flow structure, for example by mutual exclusion of velocity fields from the other region across the interface, or by local resonance near the interface. We focus here on the excluding kinds of interactions between, on the one hand, elongated and compact regions containing vortical flows and large variations in velocity, and on the other hand various kinds of weak perturbation in the surrounding external flow region: rotational, irrotational; time-varying, steady; large, small; coplanar, non-coplanar; non-diffusive, diffusive. It is shown how all these kinds of external disturbances can be wholly, or partially, ‘blocked’ at the interface with the vortical region, so that beyond a certain sheltering distance into the interior of this region the fluctuations can be very small. For the special case of quasi-parallel co-planar external straining motions outside non-directional shear flows, weak sheltering occurs if the mean velocity of the shear flow increases – otherwise the perturbations are amplified. For non-parallel flows, the sheltering effect can be greater when the vorticity is distributed in thin vortex sheets. The mechanism whereby the vortical flow induces ‘blocking’ and ‘shear-sheltering’ effects can be quantitatively explained in terms of the small adjustments of the vorticity in the vortical layers, and in some cases by the change in impulse of these layers. If the vorticity in the outer part of the vortical region is weak, it can be ‘stripped away’ by the external disturbances until the remaining vorticity is strong enough to ‘block’ the disturbances and shelter the inner flow of the vortical region. The mechanisms presented here appear to explain on the one hand some aspects of the observed robustness of vortical structures and jet or plume like shear flows in turbulent and geophysical flows, and on the other hand the levels of external perturbation needed to erode or breakdown turbulent shear flows.     

An experimental study of a radial wall jet formed by the normal impingement of a round synthetic jet

The flow field of a radial wall jet created by the impingement of a round synthetic jet normal to a flat surface was characterized using hot-wire anemometry. In the synthetic wall jets the width of the outer layer was observed to increase linearly with the radial distance along the wall, while the local maximum velocity varied inversely. The synthetic wall jet exhibits self-similar behavior as distinguished by the collapse of the mean and rms velocity profiles when normalized by the outer layer scaling variables. Increasing the actuator driving amplitude at a fixed frequency (i) increased the growth rate of the outer layer, and (ii) decreased the decay rate of the local velocity maximum. The flow field of the synthetic wall jet was dominated by vortical structures associated with the actuator driving frequency, and harmonics connected with the interaction of the produced vortex structures. For the actuator conditions investigated, neither the classical laminar nor fully turbulent analytical solutions for continuous wall jets were amenable to modeling the synthetic wall jet profile due to the transitional and unsteady nature of the synthetic wall jet.     

A bypass wake induced laminar/turbulent transition

The process of laminar to turbulent transition induced by a von Karman vortex street wake, was studied for the case of a flat plate boundary layer. The boundary layer developed under zero pressure gradient conditions. The vortex street was generated by a cylinder positioned in the free stream. An X-type hot-wire probe located in the boundary layer, measured the streamwise and normal to the wall velocity components. The measurements covered two areas; the region of transition onset and development and the region where the wake and the boundary layer merged producing a turbulent flow. The evolution of Reynolds stresses and rms-values of velocity fluctuations along the transition region are presented and discussed. From the profiles of the Reynolds stress and the mean velocity profile, a ‘negative' energy production region along the transition region, was identified. A quadrant splitting analysis was applied to the instantaneous Reynolds stress signals. The contributions of the elementary coherent structures to the total Reynolds stress were evaluated, for several x-positions of the near wall region. Distinct regions in the streamwise and normal to the wall directions were identified during the transition.     

Active Control of Swirling Coaxial Jet Mixing with Manipulation of Large-Scale Vortical Structures

Large-scale vortical structures and associated mixing in methane/air swirling coaxial jets are actively controlled by manipulating the outer shear layer of the outer swirling coaxial jet with miniature flap actuators. In order to investigate the control mechanisms, stereoscopic particle image verocimetry (stereo-PIV) and plannar laser-induced fluorescence (PLIF) techniques are employed. It is found that intense vortex rings are produced in the outer shear layer in phase with the periodic flap motion regardless of the swirl number examined. The vortical structures in the inner shear layer, however, are strongly dependent on the swirl rate. This is because the central methane jet is accelerated by the negative axial pressure gradient, of which strength is determined by the swirl. As a result, the inner vortex formation is significantly suppressed at a higher swirl rate. On the other hand, at a relatively low swirl rate, the inner vortices are shed continuously and the methane jet is pinched off. This particular mode promotes the mixing of methane and air, so that the flammable mixture can be formed at an earlier stage of the jet flow development. In addition, the evolution of secondary streamwise vortices is prompted by the combination of the periodic vortex ring shedding and the swirl. They also contribute to the mixing enhancement in the downstream region.     

PIV analysis of near-wake behind a sphere at a subcritical Reynolds number

The vortical structure of near-wake behind a sphere is investigated using a PIV technique in a circulating water channel at Re = 11,000. The measured velocity fields show a detailed vortical structure in the recirculation region such as recirculation vortices, reversed velocity zone, and out-of-plane vorticity distribution. The vorticity distribution of the sphere wake shows waviness in cross-sectional planes. The time-averaged turbulent structures are consistent with the visualized flow showing the onset of shear layer instability. The spatial distributions of turbulent intensities provide turbulent statistics for validating numerical predictions.     

Planar velocity measurements of a two-dimensional compressible wake

The present study describes the application of particle image velocimetry (PIV) to investigate the compressible flow in the wake of a two-dimensional blunt base at a freestream Mach number M_X=2. The first part of the study addresses specific issues related to the application of PIV to supersonic wind tunnel flows, such as the seeding particle flow-tracing fidelity and the measurement spatial resolution. The seeding particle response is assessed through a planar oblique shock wave experiment. The measurement spatial resolution is enhanced by means of an advanced image-interrogation algorithm. In the second part, the experimental results are presented. The PIV measurements yield the spatial distribution of mean velocity and turbulence. The mean velocity distribution clearly reveals the main flow features such as expansion fans, separated shear layers, flow recirculation, reattachment, recompression and wake development. The turbulence distribution shows the growth of turbulent fluctuations in the separated shear layers up to the reattachment location. Increased velocity fluctuations are also present downstream of reattachment outside of the wake due to unsteady flow reattachment and recompression. The instantaneous velocity field is analyzed seeking coherent flow structures in the redeveloping wake. The instantaneous planar velocity and vorticity measurements return evidence of large-scale turbulent structures detected as spatially coherent vorticity fluctuations. The velocity pattern consistently shows large masses of fluid in vortical motion. The overall instantaneous wake flow is organized as a double row of counter-rotating structures. The single structures show vorticity contours of roughly elliptical shape in agreement with previous studies based on spatial correlation of planar light scattering. Peak vorticity is found to be five times higher than the mean vorticity value, suggesting that wake turbulence is dominated by the activity of large-scale structures. The unsteady behavior of the reattachment phenomenon is studied. Based on the instantaneous flow topology, the reattachment is observed to fluctuate mostly in the streamwise direction suggesting that the unsteady separation is dominated by a pumping-like motion.     

The Vortical Structures in the Rear Separation and Wake Produced by a Supersonic Micro-Ramp

The vortical structures in the rear separation and wake region produced by a micro-ramp that immersed in a supersonic turbulent boundary layer are investigated. The small scale separation close to the trailing edge was revealed and this confirms the previous experimental observation. Between the reverse region and surrounding fast moving flow, a three-dimensional shear layer was formed, and vortices are generated. By using vortex line method, the spiral points were illustrated as the cross-sections of the Ω-shaped vortices that follow the shape of the separation. The vortical structure was analogous to that in the wake region, where similar Ω-shaped vortex which follows the deficit region caused by the micro-ramp can be observed. Finally, the revealed flow topology was conceived new and beneficial to the studying of wall bounded turbulence which involves similar vortical structures but in a smaller scale, compared with the vortical pattern in the current micro-ramp wake.     

Time resolved flow-field measurements of a turbulent mixing layer over a rectangular cavity

High Reynolds number, low Mach number, turbulent shear flow past a rectangular, shallow cavity has been experimentally investigated with the use of dual-camera cinematographic particle image velocimetry (CPIV). The CPIV had a 3 kHz sampling rate, which was sufficient to monitor the time evolution of large-scale vortices as they formed, evolved downstream and impinged on the downstream cavity wall. The time-averaged flow properties (velocity and vorticity fields, streamwise velocity profiles and momentum and vorticity thickness) were in agreement with previous cavity flow studies under similar operating conditions. The time-resolved results show that the separated shear layer quickly rolled-up and formed eddies immediately downstream of the separation point. The vortices convect downstream at approximately half the free-stream speed. Vorticity strength intermittency as the structures approach the downstream edge suggests an increase in the three-dimensionality of the flow. Time-resolved correlations reveal that the in-plane coherence of the vortices decays within 2–3 structure diameters, and quasi-periodic flow features are present with a vortex passage frequency of ~1 kHz. The power spectra of the vertical velocity fluctuations within the shear layer revealed a peak at a non-dimensional frequency corresponding to that predicted using linear, inviscid instability theory.     

A vorticity stretching diagnostic for turbulent and transitional flows

Vorticity stretching in wall-bounded turbulent and transitional flows has been investigated by means of a new diagnostic measure, denoted by Γ, designed to pick up regions with large amounts of vorticity stretching. It is based on the maximum vorticity stretching component in every spatial point, thus yielding a three-dimensional scalar field. The measure was applied in four different flows with increasing complexity: (a) the near-wall cycle in an asymptotic suction boundary layer (ASBL), (b) K-type transition in a plane channel flow, (c) fully turbulent channel flow at Re _τ = 180 and (d) a complex turbulent three-dimensional separated flow. Instantaneous data show that the coherent structures associated with intense vorticity stretching in all four cases have the shape of flat ‘pancake’ structures in the vicinity of high-speed streaks, here denoted ‘h-type’ events. The other event found is of ‘l-type’, present on top of an unstable low-speed streak. These events (l-type) are further thought to be associated with the exponential growth of streamwise vorticity in the turbulent near-wall cycle. It was found that the largest occurrence of vorticity stretching in the fully turbulent wall-bounded flows is present at a wall-normal distance of y ⁺?=?6.5, i.e. in the transition between the viscous sublayer and buffer layer. The associated structures have a streamwise length of ~200–300 wall units. In K-type transition, the Γ-measure accurately locates the regions of interest, in particular the formation of high-speed streaks near the wall (h-type) and the appearance of the hairpin vortex (l-type). In the turbulent separated flow, the structures containing large amounts of vorticity stretching increase in size and magnitude in the shear layer upstream of the separation bubble but vanish in the backflow region itself. Overall, the measure proved to be useful in showing growing instabilities before they develop into structures, highlighting the mechanisms creating high shear region on a wall and showing turbulence creation associated with instantaneous separations.