An experimental investigation of moderate reynolds number flow in a T-Channel

An experimental investigation of water flow in a T-shaped channel with rectangular cross section (20 × 20 mm inlet ID and 20 × 40 mm outlet ID) has been conducted for a Reynolds number Re range of 56–422, based on inlet diameter. Dynamical conditions and the T-channel geometry of the current study are applicable to the microscale. 2-D planar particle imaging velocimetry (PIV) and laser-induced fluorescence (LIF) were used in multiple locations of the T-channel to investigate local dynamical behaviors. Steady symmetric and asymmetric flow regimes predicted in the literature, which is largely numerical, are experimentally verified. Unsteady flow regimes, which are numerically predicted to occur at higher Re but have not yet been experimentally characterized, are also examined, and real-time LIF results illuminate the evolution of unsteady structure. Experimental data of the present resolution and scope are not presently available for unsteady flow regimes. Time scales are presented for unsteady flow regimes, which are found to exhibit periodic behavior and to occur for Re  ≥ 195. An unsteady symmetrical regime is identified for Re ≥ 350 that is detrimental to mixing. Momentum fields and dynamical behaviors of all flow regimes are characterized in detail, such that published mixing trends may be better understood. Results of all experimental trials were used to construct a regime map. A symmetric topology is found to be dominant for Re from 56 to 116, when flow is steady, and 350 to 422, when flow is characterized by unsteady stagnation-point oscillation in the T-channel junction. Asymmetric flow, which is positively indicated for mixing, is dominant for Re between 142 and 298, and the fluid interface exhibits both steady (two standing vortices) and unsteady (shear-layer type roll-up) behaviors. This result is based on multiple experiments and suggests a practical operating range of 142  ≤ Re ≤ 298 where asymmetric flow is highly likely to experimentally occur. The identification of an upper limit on Re,  beyond which mixing appears negatively impacted by a more symmetrical momentum field, is practically important as pressure drops on the microscale are significant.     

The effect of chamber length and Reynolds number on jet precession

The jet axial velocity field exiting from a nozzle/chamber configuration with an expansion ratio of 5 is investigated using Stereo-PIV for a range of chamber lengths and Reynolds (Re) numbers. The jet can exit the chamber in axial jet (AJ) mode with the maximum velocity near the chamber axis or precessing jet (PJ) mode with the maximum velocity near the chamber wall and rotating or precessing about the chamber axis. Algorithms were developed to determine the jet mode from exit conditions and allow conditional averaging of the velocity field in PJ mode. The probability of the jet in PJ mode was found to be a strong function of chamber length, L/D and only a mild function of Re for Re > 10,000. High precession probability was found for chambers of length in the range 2 < L/D < 2.75 for all cases for Re > 10,000. An abrupt reduction in precession probability occurred for chamber lengths L/D~3. For increasing chamber lengths, an increase in precession probability was observed. The ratio of entrainment-into-the-chamber of surrounding fluid to jet exit fluid was found not to be a function of Re or jet mode (AJ or PJ) but only a function of L/D. A maximum ratio entrainment-into-the-chamber was observed to occur in the range 2 < L/D < 2.5. Conditionally averaged velocity profiles also showed the exiting jet to be a strong function of L/D and with only a mild effect of Re for all cases of Re > 10,000.     

Flow behaviour of negatively buoyant jets in immiscible ambient fluid

In this paper we investigate experimentally the injection of a negatively buoyant jet into a homogenous immiscible ambient fluid. Experiments are carried out by injecting a jet of dyed fresh water through a nozzle in the base of a cylindrical tank containing rapeseed oil. The fountain inlet flow rate and nozzle diameter were varied to cover a wide range of Richardson Ri (8 × 10⁻⁴ < Ri < 1.98), Reynolds Re (467 < Re < 5,928) and Weber We (2.40 < We < 308.56) numbers. Based on the Re, Ri and We values for the experiments, we have determined a regime map to define how these values may control the occurrence of the observed flow types. Whereas Ri plays a stronger role when determining the maximum penetration height, the effect of the Reynolds number is stronger predicting the flow behaviour for a specific nozzle diameter and injection velocity.     

Experimental investigation of a scaled-up passive micromixer with uneven interdigital inlet and teardrop obstruction elements

Micromixers are vital components in micro total analysis systems. It is desirable to develop micromixers which are capable of rapidly mixing two or more fluids in a small footprint area, while minimizing mechanical losses. A novel planar scaled-up passive micromixer is experimentally investigated in this study. The design incorporates a 7-substream uneven interdigital inlet which supplies two liquid species in a parallel arrangement and promotes diffusion along the side walls. Forty-eight staggered teardrop-shaped obstruction elements located along the channel length combined with 32 side walls protrusions increase the two-fluid interfacial area while converging the flow due to periodic reductions in cross-sectional area. The scaled-up micromixer has a mixing channel length of 110 mm with a mixing channel height and width of 2 and 5 mm, respectively. Experimental investigations are carried out at four locations along the channel length and at Reynolds numbers equal to 1, 5, 10, 25, 50, and 100, where the Reynolds number is calculated based on total two-fluid flow and the mixing channel hydraulic diameter. Flow visualization is employed to study flow patterns, while induced fluorescence (IF), using de-ionized water and low concentration Rhodamine 6G solutions, provides mixing efficiency data. Results show a change in dominant mixing mechanism from mass diffusion to mass advection, with a critical Reynolds number of 25. At high Reynolds numbers, the formation of additional lamellae is observed, as is the formation of Dean vortices in the vicinity of the teardrop obstructions. Of the tested cases, the highest outlet mixing efficiency, 68.5%, is achieved at a Reynolds number of 1, where mass diffusion dominates. At low Reynolds numbers, superior mixing efficiency is due primarily to the implementation of the uneven interdigital inlet. A comparable mixing length is proposed to allow for reasonable comparison with published studies.     

Investigation of slip boundary conditions in the T-shaped microchannel

Flow regimes and mixing performance in a T-type micromixer at high Reynolds numbers were studied by numerical solution of the Navier–Stokes equations. The Reynolds number was varied from 1 to 1000. The cross section of the mixing channel was 100 μm × 200 μm, and its length was 1400 μm. The transverse inlet channels were symmetric about the mixing channel, and their cross-section was 100 μm × 100 μm, and the total length was 800 μm. Five different flow regimes were identified: (i) stationary vortex-free flow (Re < 5); (ii) stationary symmetric vortex flow with two horseshoe vortices (5 < Re < 150); (iii) stationary asymmetric vortex flow (150 < Re < 240); (iv) non-stationary periodic flow (240 < Re < 400); and (v) stochastic flow (Re > 400). Maximum mixing efficiency is obtained for nonstationary asymmetric vortex flow. In this case, an S-shaped vortex structure is formed in the flow field. The effect of the slip conditions on the flow pattern and mixing efficiency is studied. The slip length varied from 1 to 70 μm in the calculations. It was shown that the mixing can be controlled by hydrophobic coating.     

Temperature fields induced by direct contact condensation of steam in a cross-flow in a channel

The temperature fields in the center plane of a channel with a square cross-section have been measured. Steam injected at relatively low mass fluxes through a small hole in one of the walls of the channel condensed intermittently in a small area close to the inlet. The upstream temperature of the liquid cross-flow, T _L , the momentum ratio, J, and the Prandtl number proved to be important for the single-phase temperature field induced in the jet further away from the steam inlet. Jet centerlines of velocity and temperature are measured and positions are compared. Different locations for J < 100 and low T _L are explained from dependencies on Reynolds and Prandtl numbers. Next to the jet centerline a second high-temperature zone was found to occur, close to the wall and downstream of the steam inlet. The importance of capillary forces is investigated with the aid of 3D CFD computations.     

Numerical investigation of turbulent free jet flows issuing from rectangular nozzles: the influence of small aspect ratio

In this research, the fluid and thermal characteristics of a rectangular turbulent jet flow is studied numerically. The results of three-dimensional jet issued from a rectangular nozzle are presented. A numerical method employing control volume approach with collocated grid arrangement was employed. Velocity and pressure fields are coupled with SIMPLEC algorithm. The turbulent stresses are approximated using k–e{\varepsilon} model with two different inlet conditions. The velocity and temperature fields are presented and the rates of their decay at the jet centerline are noted. The velocity vectors of the main flow and the secondary flow are illustrated. Also, effect of aspect ratio on mixing in rectangular cross-section jets is considered. The aspect ratios that were considered for this work were 1:1 to 1:4. The results showed that the jet entrains more with smaller AR. Special attention has been drawn to the influence of the Reynolds number (based on hydraulic diameter) as well as the inflow conditions on the evolution of the rectangular jet. An influence on the jet evolution is found for smaller Re, but the jet is close to a converged state for higher Reynolds numbers. The inflow conditions have considerable influence on the jet characteristics.     

High-resolution imaging of turbulence structures in jet flames and non-reacting jets with laser Rayleigh scattering

A comparative study of the length scales and morphology of dissipation fields in turbulent jet flames and non-reacting jets provides a quantitative analysis of the effects of heat release on the fine-scale structure of turbulent mixing. Planar laser Rayleigh scattering is used for highly resolved measurements of the thermal and scalar dissipation in the near fields of CH₄/H₂/N₂ jet flames (Re _d  = 15,200 and 22,800) and non-reacting propane jets (Re _d  = 7,200–21,700), respectively. Heat release increases the dissipation cutoff length scales in the reaction zone of the flames such that they are significantly larger than the cutoff scales of non-reacting jets with comparable jet exit Reynolds numbers. Fine-scale anisotropy is enhanced in the reaction zone. At x/d = 10, the peaks of the dissipation angle PDFs in the Re _d  = 15,200 and 22,800 jet flames exceed those of non-reacting jets with corresponding jet exit Reynolds numbers by factors of 2.3 and 1.8, respectively. Heat release significantly reduces the dissipation layer curvature in the reaction zone and in the low-temperature periphery of the jet flames. These results suggest that the reaction zone shields the outer regions of the jet flame from the highly turbulent flow closer to the jet axis.     

Experimental investigations of a swirling jet in both stationary and rotating surroundings

The ‘plug’ flow emerging from a long rotating tube into a large stationary reservoir was used in the experimental investigation of swirling jets with Reynolds numbers, Re = 600, 1,000 and 2,000, and swirl numbers, S = ΩR/U, in the range 0–1.1, to cover flow regimes from the non-rotating jet to vortex breakdown. Here Ω is the nozzle rotation rate, R is the radius of the nozzle exit, and U is the mean mass axial velocity. The jet was more turbulent and eddies shed faster at larger Re. However the flow criticality and shear layer morphology remained unchanged with Re. After the introduction of sufficient rotation, co-rotating and counter-winding helical waves replaced vortex rings to become the dominant vortex structure. The winding direction of the vortex lines suggests that Kelvin–Helmholtz and generalized centrifugal instability dominated the shear layer. A quantitative visualization study has been carried out for cases where the reservoir was rotating independently with S _a  = Ω _a R/U = ±0.35, ±0.51 and ±0.70 at Re = 1,000 and 2000, where Ω _a is the rotation rate of the reservoir. The criterion for breakdown was found to be mainly dependent on the absolute swirl number of the jet, S. This critical swirl number was slightly different in stationary and counter-swirl surroundings but obviously smaller when the reservoir co-rotated, i.e. S _c  = 0.88, 0.85 and 0.70, respectively. These results suggest that the flow criticality depends mainly on the velocity distributions of the vortex core, while instabilities resulting from the swirl difference between the jet and its ambient seem to have only a secondary effect.     

Modified flapping jet for increased jet spreading using synthetic jets

The present paper is an experimental investigation, using a PIV system, on modified rectangular jet flow co-flowing with a pair of synthetic jets placed symmetrically with respect to the geometric centerline of the main flow. The objective was to determine the optimal forcing conditions that would result in jet spreading beyond what would be obtained in a simple flapped jet. The main jet had an exit Re_h = 36,000, based on the slot height, h. The synthetic jets were operated in a periodic manner with a periodic momentum coefficient of about 3.3% and at a frequency of the main jet preferred mode. A short, wide angle diffuser of half angle of about 45° was attached to the main jet. Generally for the vectored jet, much of the flow features found here resembled those reported in the literature except that the deflection angle in this study increased with downstream distances inside the diffuser and then remained roughly unchanged thereafter. Larger jet spreading was achieved when the main jet was subjected to simultaneous actuation of the synthetic jets but the flow did not achieve the initial jet spreading that was observed in the vectored jet. Further jet spreading was achieved when the synthetic jets were alternately actuated in which each synthetic jet was actuated for a number of cycles before switching. This technique allowed the jet to flap across the flow between transverse positions larger than what would be obtained in a simple flip-flop jet. Under the present flow geometry and Reynolds number, it was found that when the ratio f_s/f_al, where f_s is the synthetic actuation frequency and f_al is the alternating frequency, was larger than 10, the mean streamwise velocity of the main jet had two peaks symmetrically placed with respect to the jet axis and the jet had the appearance of flowing into two streams each moving nearly parallel to the diffuser wall. For a value of f_s/f_al of about 10, the optimal value in this study, the desired flow properties were achieved in that, the mean velocity was nearly uniform with an increase in the jet width compared to the simultaneous actuations, and the jet flapping was more effective in redistributing and homogenizing the turbulent kinetic energy across the main jet.     

Near-field instability of variable property jet in normal gravity and microgravity fields

The near-field instability of variable property jets of air, CO₂, and He, issued into the ambient air, has been investigated experimentally within normal gravity and microgravity fields. The density ratio to the ambient air is unity for air jets, more than unity (1.53) for CO₂ jets, and less than unity (0.14) for He jets, respectively. The ratio of kinematic viscosity to the ambient air is unity for air jets, less than unity for CO₂ jets (0.53), and more than unity for He jets (7.75), respectively. The jet velocity is varied from 0.4 to 1.8 m/s and then the jet Reynolds number varies from 60 for Helium jet to 2,000 for CO₂ jet, while the Richardson number varies from negative to positive values. The motion of the jet is visualized using a laser tomographic method and recorded by a high-speed digital video camera with 250 frames/s. The result shows that the instability of the jet is intensified when Re > 800 while it is weakened at Re < 800 at the microgravity field, indicating that the viscosity plays an important role in weakening the instability. Under a normal gravity field, the buoyancy also becomes important. In order to quantify the instability criteria, the quantity of the instability is introduced, which consists of the Kelvin–Helmholtz instability, buoyancy effect and viscous effect. When the ratio of the sum of Kelvin–Helmholtz and buoyancy forces to viscous force exceeds a certain value, around 12 in the present study, the jet becomes unstable even when Re < 800. These results reveal that the instability of variable property jets is influenced by the Kelvin–Helmholtz instability, the viscous effect and the buoyancy effect.     

The effect of sidewalls on rectangular jets

An experimental study is presented regarding the influence of sidewalls on the turbulent free jet flow issuing from a smoothly contracting rectangular nozzle of aspect ratio 15. “Sidewalls” are two parallel plates, flush with each of the slots’ short sides, practically establishing bounding walls extending the nozzle sidewalls in the downstream direction. Measurements of the streamwise and lateral velocity mean and turbulent characteristics have been accomplished, with an x-sensor hot wire anemometer, up to an axial distance of 35 nozzle widths, for jets with identical inlet conditions with and without sidewalls. Centreline measurements for both configurations have been collected for three Reynolds numbers, Re_D = 10,000, 20,000 and 30,000. For Re_D = 20,000 measurements in the transverse direction were collected at 13 different downstream locations in the range, x = 0–35 nozzle widths, and in the spanwise direction at three different downstream locations, x = 2, 6 and 25 nozzle widths.Results indicate that, the two jet configurations (with and without sidewalls) produce statistically different flow fields. Sidewalls do not lead to the production of a 2D flow field as undulations in the spanwise mean velocity distribution indicate. They do increase the two-dimensionality of the jet increasing the longevity of 2D spanwise rollers structures formed in the initial stages of entrainment, which are responsible for the convection of longitudinal momentum towards the outer field, establishing larger streamwise mean velocities at the jet edges. In the near field, up to 25 nozzle widths, lower outward lateral velocities in the presence of the sidewalls are held responsible for the decrease of turbulent terms including rms of velocity fluctuations and Reynolds stresses. Skewness factors increase monotonically across the shear layers from negative values to positive forming sharp peaks at the outer edges of the jet, illustrative of the presence of well defined 2D roller structures in the jet with sidewalls.     

Solitary waves on inclined films: their characteristics and the effects on wall shear stress

The properties of solitary waves, developing from inlet disturbances of controlled frequency along an inclined film flow, are systematically studied experimentally and computationally. Time-variations of film height and wall shear stress are measured, using respectively a capacitance probe and an electrodiffusion sensor. Computational data are provided from simulations performed by a Galerkin finite element scheme. The height and spacing of solitary humps, their phase velocity and the wavelength of the preceding capillary ripples are reported as functions of the Reynolds number (10<Re<100) and the inlet frequency (0.5 Hz< f<2.5 Hz). The wall shear stress modulation imposed by the passage of solitary waves is studied experimentally and computationally as a function of Re. Distinct nonlinear characteristics are noted, including a steep maximum and a negative minimum, with the effects intensifying at intermediate Re. All computer predictions are found to be in good quantitative agreement with the experimental data.     

The entrainment behavior of a turbulent axisymmetric jet in a viscous host fluid

Although turbulent jets have been studied extensively, one configuration that has not received much attention is the viscosity-stratified jet, wherein a turbulent jet of lower viscosity issues into a density-matched host liquid of higher viscosity. We present experimental data for scalar dispersion and two-dimensional velocity measurements in the axial plane of a turbulent axisymmetric jet with a Reynolds number (Re) of 2,000 issuing into a viscous host liquid at viscosity ratios (m) ranging from 1 to 55. The presence of a strong viscosity discontinuity across the jet edge results in a significant decrease in the scalar spread rate. We attribute this to the rapid reduction in turbulence intensity and the suppression of large engulfing eddies at the jet edge. The velocity profile, on the other hand, indicates that the velocity width and mass flux reduce with increasing m up to about 20, but then increase for higher values of m. This non-monotonic variation is explained by the growing influence of viscous stress for m>20. The scalar spread rate, the velocity spread rate, the centerline velocity decay rate, and the jet mass flux are all minimized for m20 for Re=2,000.

---

The effects of acceleration in jets: kinematics of the near field vortices

Direct and large-eddy simulations (DNS/LES) of accelerating round jets are used to analyze the effects of acceleration on the kinematics of vortex rings in the near field of the jet (x/D < 12). The acceleration is obtained by increasing the nozzle jet velocity with time, in a previously established (steady) jet, and ends once the inlet jet velocity is equal to twice its initial value. Several acceleration rates (α = 0.02–0.6) and Reynolds numbers (Re _D = 500–20000) were simulated. Acceleration maps were used to make a detailed study of the kinematics of vortex rings in accelerating jets. One of the effects of the acceleration is to cause a number of new primary and secondary vortex merging events that are absent from steady jets. As the acceleration rate α increases, both the number of primary merging events between rings and the axial position where these take place decreases. The statistics for the speed of the starting ring that forms at the start of the acceleration phase for each simulation, agree well with the statistics for the “front” speed observed by Zhang and Johari (Phys Fluids 8:2185–2195, 1996). Acceleration maps and flow visualizations show that during the acceleration phase the near field coherent vortices become smaller and are formed at an higher frequency than in the steady jet, and their (mean) shedding frequency increases linearly with the acceleration rate. Finally, it was observed that the acceleration decreases the spreading rate of the jet, in agreement with previous experimental works.      

Secondary flow patterns and mixing in laminar pulsating flow through a curved pipe

Mixing by secondary flow is studied by particle image velocimetry (PIV) in a developing laminar pulsating flow through a circular curved pipe. The pipe curvature ratio is η = r ₀/r _c  = 0.09, and the curvature angle is 90°. Different secondary flow patterns are formed during an oscillation period due to competition among the centrifugal, inertial, and viscous forces. These different secondary-flow structures lead to different transverse-mixing schemes in the flow. Here, transverse mixing enhancement is investigated by imposing different pulsating conditions (Dean number, velocity ratio, and frequency parameter); favorable pulsating conditions for mixing are introduced. To obviate light-refraction effects during PIV measurements, a T-shaped structure is installed downstream of the curved pipe. Experiments are carried out for the Reynolds numbers range 420 ≤ Re_st ≤ 1,000 (Dean numbers 126.6 ≤ Dn ≤ 301.5) corresponding to non-oscillating flow, velocity component ratios 1 ≤ (β = U _max,osc/U _m,st) ≤ 4 (the ratio of velocity amplitude of oscillations to the mean velocity without oscillations), and frequency parameters 8.37 < (α = r ₀(ω/ν)^0.5) < 24.5, where α² is the ratio of viscous diffusion time over the pipe radius to the characteristic oscillation time. The variations in cross-sectional average values of absolute axial vorticity (|ζ|) and transverse strain rate (|ε|) are analyzed in order to quantify mixing. The effects of each parameter (Re_st, β, and α) on transverse mixing are discussed by comparing the dimensionless vorticities (|ζ _P |/|ζ _S |) and dimensionless transverse strain rates (|ε _P |/|ε _S |) during a complete oscillation period.     

Moderate Reynolds number axisymmetric jet development downstream an extended conical diffuser: Influence of extension length

Low-velocity (bulk velocity of 4.4 m/s) and moderate Reynolds (7350) axisymmetrical jet development is studied by hot-film single sensor anemometry. The jet issues from a conical convergent-divergent diffuser with uniform extension (diameter 25 mm). Decreasing the length-to-diameter ratio of the extension tube from 20 down to 0.4 is shown to alter severely the mean velocity profile at the tube outlet from Blasius to top-hat whereas turbulence intensities increases from 6 up to 50%. Next, the influence of the initial velocity profile at the tube outlet on axisymmetrical jet development is assessed. The velocity development exhibits a self-similar far field characteristic for axisymmetrical jet development. Although, the jet centerline decay constant increases and the jet spreading rate decreases as length-to-diameter ratios are increased from 0.4 up to 7.2 for which the initial centerline velocity decreases. Therefore, scaling of the centerline decay constant and inverse scaling of the spreading rate with initial centerline velocity U₀ or initial velocity Reynolds number Re₀ reported for moderate Reynolds numbers and low initial turbulence levels [18,22] does not hold as the turbulence level exceeds a threshold value in the range 12–27%. In addition, the influence of initial conditions on near and far field turbulence properties is shown. A transition in near field behaviour is observed for length-to-diameters around 3.6. Flow and geometrical configurations under study are relevant to e.g. upper airway flow.     

PIV–LES analysis of channel flow rotating about the streamwise axis

The flow field of a channel rotating about the streamwise axis is analyzed experimentally and numerically. The current investigations were carried out at a bulk velocity based Reynolds number of Re_m = 2850 and a friction velocity based Reynolds number of Re_τ = 180, respectively. Particle-image velocimetry (PIV) measurements are compared with large-eddy simulation data to show earlier direct numerical simulation findings to generate too large a reverse flow region in the center region of the spanwise flow. The development of the mean spanwise velocity distribution and the influence of the rotation on the turbulent properties, i.e., the Reynolds stresses and the two-point correlations of the flow, are confirmed in both investigations. The rotation primarily influences those components of the Reynolds shear stresses, which contain the spanwise velocity component. The size of the correlation areas and thus the length scales of the flow generally grow in all three coordinate directions leading to longer structures. Furthermore, experimental results of the same channel flow at a significantly lower bulk Reynolds number of Re_{m, l} = 665, i.e., a laminar flow in a non-rotating channel, are introduced. The experiments show the low Reynolds number flow to become turbulent under rotation and to develop the same characteristics as the high Reynolds number flow.