Non-linear dynamics of semi-dilute polydisperse polymer solutions in microfluidics: A study of a benchmark flow problem

The complex flow behaviour of semi-dilute (15 < c/c^* < 22.5) polydisperse polyethylene oxide (PEO) aqueous solutions flowing through a planar microfluidic geometry with an 8:1:8 contraction-expansion is systematically studied. The molecular weight and distribution of the PEO samples are analysed by Gel Permeation Chromatography (GPC). Full rheometric characterizations using various techniques including piezoelectric axial vibrator (PAV) measurements at frequencies as high as 6700 Hz are carried out for one semi-dilute PEO solution. Complex flows over a wide range of elasticity numbers (20 ? El ? 120), Weissenberg numbers (7 ? Wi ? 121) and Reynolds numbers (0.08 ? Re ? 4.5) are characterized using micro-particle image velocimetry (μ-PIV) and pressure drop measurements. The evolution of vortex formation and dynamics has been visualized through a step-flow-rate experiment. The effect of El on vortex stability has been studied. Various flow dynamics regimes have been quantified and are presented in a Wi-Re diagram. The experimental results reveal that the elastic behaviour of polymer solutions is very sensitive to high molecular weight polymer in the polydisperse polymer samples, and the contraction ratio and the aspect ratio of flow geometry are the important design parameters in controlling the non-linear dynamics of semi-dilute polymer solutions in microfluidics.     

Reynolds Number Effects on the Coherent Dynamics of the Turbulent Horseshoe Vortex System

The adverse pressure gradient induced by a surface-mounted obstacle in a turbulent boundary layer causes the approaching flow to separate and form a dynamically rich horseshoe vortex system (HSV) in the junction of the obstacle with the wall. The Reynolds number of the flow (Re) is one of the important parameters that control the rich coherent dynamics of the vortex, which are known to give rise to low-frequency, bimodal fluctuations of the velocity field (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). We carry out detached eddy simulations (DES) of the flow past a circular cylinder mounted on a rectangular channel for Re = 2.0 × 10⁴ and 3.9 × 10⁴ (Dargahi, Exp Fluids 8:1–12, 1989) in order to systematically investigate the effect of the Reynolds number on the HSV dynamics. The computed results are compared with each other and with previous experimental and computational results for a related junction flow at a much higher Reynolds number (Re = 1.15 × 10⁵) (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). The computed results reveal significant variations with Re in terms of the mean-flow quantities, turbulence statistics, and the coherent dynamics of the turbulent HSV. For Re = 2.0 × 10⁴ the HSV system consists of a large number of necklace-type vortices that are shed periodically at higher frequencies than those observed in the Re = 3.9 × 10⁴ case. For this latter case the number of large-scale vortical structures that comprise the instantaneous HSV system is reduced significantly and the flow dynamics becomes quasi-periodic. For both cases, we show that the instantaneous flowfields are dominated by eruptions of wall-generated vorticity associated with the growth of hairpin vortices that wrap around and disorganize the primary HSV system. The intensity and frequency of these eruptions, however, appears to diminish rapidly with decreasing Re. In the high Re case the HSV system consists of a single, highly energetic, large-scale necklace vortex that is aperiodically disorganized by the growth of the hairpin mode. Regardless of the Re, we find pockets in the junction region within which the histograms of velocity fluctuations are bimodal as has also been observed in several previous experimental studies.     

Dynamics of hairpin vortices generated by a mixing tab in a channel flow

To better understand mixing by hairpin vortices, time-series particle image velocimetry (PIV) was applied to the wake of a trapezoidal-shaped passive mixing tab mounted at the bottom of a square turbulent channel (Re _h =2,080 based on the tab height). Instantaneous velocity/vorticity fields were obtained in sequences of 10 Hz in the tab wake in the center plane (x–y) and in a plane (x–z) parallel to the wall. Periodically-shed hairpin vortices were clearly identified and seen to rise as they advected downstream. Experimental evidence shows that the vortex-induced ejection of the near-wall viscous fluid to the immediate upstream is important to the dynamics of hairpin vortices. It can increase the strength of the hairpin vortices in the near tab region and cause generation of secondary hairpin vortices further downstream when the hairpin heads are farther away from the wall. Measurements also reveal the existence of a type of new secondary vortice with the opposite-sign spanwise vorticity. The distribution of vortex loci in the x–y plane shows that the hairpin vortices and the reverse vortices are spatially segregated in distinct layers. Turbulence statistics, including mean velocity profiles, Reynolds stresses, and turbulent kinetic energy dissipation rate distributions, were obtained from the PIV data. These statistical quantities clearly reveal imprints of the identified vortex structures and provide insight into mixing effectiveness. Received: 24 February 2000/Accepted: 24 October 2000     

Effect of acceleration on the wavy Taylor vortex flow

In this paper, we use a laser optical technique to investigate the characteristics of a wavy Taylor vortex flow between two concentric cylinders, with the inner cylinder subjected to a wide range of predetermined acceleration and the outer one at rest. We focus on the inner/outer radius ratio of 0.894, with an acceleration (dRe/dt*) from 0.1123 to 2,247, and Reynolds number from Re/Re _c =1.0 to 36. The results show that, with increasing Reynolds number, there is an initial increase in the wavelength of the wavy vortex flow (λ), and a decrease in the wave speed (c) before they asymptote to a constant value, which is a function of the acceleration. As for the wave amplitude (A), it is found that the effect of acceleration is significant only in a very narrow range of Reynolds numbers. Received: 21 August 2001 / Accepted: 22 November 2001     

Modelling the three-dimensional flow of a semi-dilute polymer solution in microfluidics—on the effect of aspect ratio

The flow of polymer solutions in microfluidic devices is inherently three-dimensional, especially in the non-linear flow regime, and often results in flow phenomena that might not even be encountered in macro-devices. Using a multi-mode Phan-Thien–Tanner model, three-dimensional (3-D) simulations of a semi-dilute polyethylene oxide (PEO) solution through 8:1 planar contraction micro-channels with various depths have been carried out to systematically study the effect of the aspect ratio on the flow fields. Vortex dynamics in the upstream flow section and excess pressure drop are quantified in detail. A transition from a salient-corner vortex mechanism to a lip vortex mechanism is observed as the aspect ratio is varied from 1 to 1/4, which corresponds to the elasticity numbers El = 36.3 to 48.3. The numerical results show that varying the aspect ratio of microfluidic channels has similar effects to varying other parameters, such as fluid properties, which influence the elasticity number. Thus, our results support the view that vortex growth mechanism is determined by the elasticity number, which is fixed for a given fluid and geometry. The principle is of significance to the design of new microfluidic chips for a wide range of applications.     

The ultra-low Reynolds number airfoil wake

Lift force and the near wake of an NACA 0012 airfoil were measured over the angle (α) of attack of 0°–90° and the chord Reynolds number (Re _c ), 5.3 × 10³–5.1 × 10⁴, with a view to understand thoroughly the near wake of the airfoil at low- to ultra-low Re _c . While the lift force is measured using a load cell, the detailed flow structure is captured using laser-Doppler anemometry, particle image velocimetry, and laser-induced fluorescence flow visualization. It has been found that the stall of an airfoil, characterized by a drop in the lift force, occurs at Re _c  ≥ 1.05 × 10⁴ but is absent at Re _c  = 5.3 × 10³. The observation is connected to the presence of the separation bubble at high Re _c but absence of the bubble at ultra-low Re _c , as evidenced in our wake measurements. The near-wake characteristics are examined and discussed in detail, including the vortex formation length, wake width, spanwise vorticity, wake bubble size, wavelength of K–H vortices, Strouhal numbers, and their dependence on α and Re _c .     

Experiments on the effect of acceleration on the drag of tapwater bubbles

Three-dimensional (3D) trajectories of spherical air bubbles passing through a converging part of a rectangular channel have been measured. Bubble diameters, d _b, were less than 1 mm and the Reynolds numbers, Re _b, for stagnant tapwater and for mean liquid velocity vˉ _L=0.25 m/s were in about same range. Received: 15 January 2001 / Accepted: 12 June 2001     

The wake of falling disks at low Reynolds numbers

We visualized the wake structure of circular disks falling vertically in quiescent water.The evolution of the wake was shown to be similar to the flow patterns behind a fixed disk.The Reynolds number,Re = Ud/ν,is in the range of 40 200.With the ascension of Reynolds numbers,a regular bifurcation occurred at the first critical Reynolds number Re c 1,leading to a transition from an axisymmetric wake structure to a plane symmetric one;A Hopf bifurcation took place at the second critical Reynolds number Re c 2,as the wake structure became unsteady.Plane symmetry of the wake structure was first lost as periodic vortex shedding appeared,but recovered at higher Reynolds number.The difference between the two critical Reynolds numbers was found to be shape-dependent,as we compared our results for thin discs with those for other falling bodies,such as spheres and cones.This observation could be understood in terms of the instability mechanism of the vortical structure.     

Experimental investigation into vortex structure and pressure drop across microcavities in 3D integrated electronics

Hydrodynamics in microcavities with cylindrical micropin fin arrays simulating a single layer of a water-cooled electronic chip stack is investigated experimentally. Both inline and staggered pin arrangements are investigated using pressure drop and microparticle image velocimetry (μPIV) measurements. The pressure drop across the cavity shows a flow transition at pin diameter–based Reynolds numbers (Re _d ) ~200. Instantaneous μPIV, performed using a pH-controlled high seeding density of tracer microspheres, helps visualize vortex structure unreported till date in microscale geometries. The post-transition flow field shows vortex shedding and flow impingement onto the pins explaining the pressure drop increase. The flow fluctuations start at the chip outlet and shift upstream with increasing Re _d . No fluctuations are observed for a cavity with pin height-to-diameter ratio h/d = 1 up to Re _d ~330; however, its pressure drop was higher than for a cavity with h/d = 2 due to pronounced influence of cavity walls.     

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.     

A rotating elliptic airfoil in fluid at rest and in a parallel freestream

This paper reports results of DPIV measurements on a two-dimensional elliptic airfoil rotating about its own axis of symmetry in a fluid at rest and in a parallel freestream. In the former case, we examined three rotating speeds (Re _c,Ω = 400, 1,000 and 2,000), and in the later case, four rotating speeds (Ro _c,Ω = 2.4, 1.2, 0.6 and 0.4), together with two freestream velocities (Re _c,u  = 200 and 1,000) and two starting configurations of the airfoil (i.e., chord parallel to (α ₀ = 0°) or normal (α ₀ = 90°) to the freestream). Results show that a rotating airfoil in a stationary fluid produces two distinct types of vortex structures depending on the Reynolds number. The first type occurs at the lowest Reynolds number (Re _c,Ω = 400), where vortices shed from the two edges or tips of the airfoil dissipated quickly, resulting in the airfoil rotating in a layer of diffused vorticity. The second type occurs at higher Reynolds numbers (i.e., Re _c,Ω = 1,000 and 2,000), where the corresponding vortices rotated together with the airfoil. Due to the vortex suction effect, the torque characteristics are likely to be heavily damped for the first type because of the rapidly subsiding vortex shedding, and more oscillatory for the second type due to persistent presence of tip vortices. In a parallel freestream, increasing the tip-speed ratio (V/U) of the airfoil (i.e., decreasing the Rossby number, Ro _c,Ω) transformed the flow topology from periodic vortex shedding at Ro _c,Ω = 2.4 to the generation of a “hovering vortex” at Ro _c,Ω = 0.6 and 0.4. The presence of the hovering vortex, which has not been reported in literature before, is likely to enhance the lift characteristics of the airfoil. Freestream Reynolds number is found to have minimal effect on the vortex formation and shedding process, although it enhances shear layer instability and produces more small-scale flow structures that affect the dynamics of the hovering vortex. Likewise, initial starting configuration of the airfoil, while affecting the flow transient during the initial phase of rotation, has insignificant effect on the overall flow topology. Unfortunately, technical constraint of our apparatus prevented us from carrying out complimentary force measurements; nevertheless, the results presented herein, which are more extensive than those computed by Lugt and Ohring (1977), will provide useful benchmark data, from which more advanced numerical calculations can be carried out to ascertain the corresponding force characteristics, particularly for those conditions with the presence of hovering vortex.     

Heat transfer from a surface fitted with rectangular blocks at different orientation angle

 An experimental investigation was carried out to study the enhancement of the heat transfer from a heated flat plate fitted with rectangular blocks of 1 × 2 × 2 cm³ dimensions in a channel flow as a function of Reynolds number (Re_h), spacing (S _y ) of blocks in the flow direction, and the block orientation angle (α) with respect to the main flow direction. The experiments were performed in a channel of 18 cm width and 10 cm height, with air as the working fluid. For fixed S _x =3.81 cm, which is the space between the blocks in transverse to the flow direction, the experimental ranges of the parameters were S _y =3.33–4.33 cm, α=0–45°, Re_h=7625–31550 based on the hydraulic diameter and the average velocity at the beginning of the test section in the channel. Correlations for Nusselt number were developed, and the ratios of heat transfer with blocks to those with no blocks were given. The results indicated that the heat transfer could be enhanced or reduced depending on the spacing between blocks, and the block orientation angle. The maximum heat transfer rate was obtained at the orientation angle of 45°. Received on 13 December 2000 / Published online: 29 November 2001     

Measurement and visualization of impingement cooling in narrow channels

Experimental measurement techniques such as naphthalene sublimation, liquid crystal thermography and real-time holographic interferometry are standard. Their application in narrow channels causes problems and is therefore limited. The channel width must not change too much because the naphthalene sublimation and the liquid crystal coating necessary for the thermography may cause non-negotiable variations. The interferometry fails in turbulent flow area. The diffraction along the channel edges is an additional difficulty. A comparison of the results obtained from the application of all three techniques, which has not been considered in earlier publications, is made here. The methods were used to measure and visualize the heat transfer characteristics of an array of 1.2 mm diameter impinging jets in an enclosed channel (≥2.2 mm) with single-sided flow-off at Reynolds numbers of about Re _z  ≈ 20,000. Scale-up ratios as low as 2.4 have been used in order to maintain similarity as it has not been previously reported. The naphthalene technique provided a high spatially resolved measurement of the Sherwood number along a downstream line. The liquid crystal thermography technique provided 2D contours of the Nusselt number. The temperature distribution within dead water zones was visualized with holographic interferometry. The cross-flow effects caused a shift in the stagnation point and a monotone decrease in the Nusselt number in the downstream direction. Received: 21 April 2000/Accepted: 6 July 2000     

Flow Bifurcations in a Thin Gap Between Two Rotating Spheres

This paper presents the use of a parameter continuation method and a test function to solve the steady, axisymmetric incompressible Navier–Stokes equations for spherical Couette flow in a thin gap between two concentric, differentially rotating spheres. The study focuses principally on the prediction of multiple steady flow patterns and the construction of bifurcation diagrams. Linear stability analysis is conducted to determine whether or not the computed steady flow solutions are stable. In the case of a rotating inner sphere and a stationary outer sphere, a new unstable solution branch with two asymmetric vortex pairs is identified near the point of a symmetry-breaking pitchfork bifurcation which occurs at a Reynolds number equal to 789. This solution transforms smoothly into an unstable asymmetric 1-vortex solution as the Reynolds number increases. Another new pair of unstable 2-vortex flow modes whose solution branches are unconnected to previously known branches is calculated by the present two-parameter continuation method. In the case of two rotating spheres, the range of existence in the (Re ₁ , Re ₂ ) plane of the one and two vortex states, the vortex sizes as a function of both Reynolds numbers are identified. Bifurcation theory is used to discuss the origin of the calculated flow modes. Parameter continuation indicates that the stable states are accompanied by certain unstable states. Received 26 November 2001 and accepted 10 May 2002 Published online 30 October 2002 Communicated by M.Y. Hussaini     

Flow of dilute to semi-dilute polystyrene solutions through a benchmark 8:1 planar abrupt micro-contraction

We have studied the flow of thermodynamically ideal solutions of a high molecular weight (M_w = 6.9 MDa) atactic polystyrene in the θ solvent dioctyl phthalate (aPS in DOP) through a micro-fabricated 8:1 planar abrupt contraction geometry. The channel is much deeper than most micro-scale geometries, providing an aspect ratio of 16:1 and a good approximation to 2D flow in the narrow channel. The solutions span a range of concentration 0.03 wt.% < c < 0.6 wt.%, encompassing the dilute to semi-dilute regimes and providing a range of fluid viscosities and relaxation times such that we achieve a range of Weissenberg numbers (8.7 < Wi < 1562) and Reynolds numbers (0.07 < Re < 11.2), giving elasticity numbers between 31 < El < 295. We study the flow using a combination of micro-particle image velocimetry (μ-PIV) to characterize the flow field, pressure measurements to evaluate the non-Newtonian viscosity, and birefringence measurements to assess macromolecular strain. Flow field observations reveal three broad flow regimes characterized by Newtonian-like flow, unstable flow and vortex growth in the upstream salient corners. Transitions between the flow regimes scale with Wi, independent of El, indicating the dominance of elastic over inertial effects in all the fluids. Transitions in the flow field are also reflected by transitions in the relative viscosity (determined from the pressure drop) and the macromolecular strain (determined from birefringence measurements). The strain through the 8:1 contraction saturates at a value of ～4–5 at high Wi. The results of these experiments broaden the limited set of literature on flow through micro-fluidic planar contractions and should be of significant value for optimizing lab-on-a-chip design and for comparison with modeling studies with elasticity dominated fluids.     

Effect of natural ventilation on the boundary layer separation and near-wake vortex shedding characteristics of a sphere

Experiments were conducted in water and wind tunnels on spheres in the Reynolds number range 6 × 10³ to 6.5 × 10⁵ to study the effect of natural ventilation on the boundary layer separation and near-wake vortex shedding characteristics. In the subcritical range of Re (<2 × 10⁵), ventilation caused a marginal downstream shift in the location of laminar boundary layer separation; there was only a small change in the vortex shedding frequency. In the supercritical range (Re > 4 × 10⁵), ventilation caused a downstream shift in the mean locations of boundary layer separation and reattachment; these lines showed significant axisymmetry in the presence of venting. No distinct vortex shedding frequency was found. Instead, a dramatic reduction occurred in the wake unsteadiness at all frequencies. The reduction of wake unsteadiness is consistent with the reduction in total drag already reported. Based on the present results and those reported earlier, the effects of natural ventilation on the flow past a sphere can be categorized in two broad regimes, viz., weak and strong interaction regimes. In the weak interaction regime (subcritical Re), the broad features of the basic sphere are largely unaltered despite the large addition of mass in the near wake. Strong interaction is promoted by the closer proximity of the inner and outer shear layers at supercritical Re. This results in a modified and steady near-wake flow, characterized by reduced unsteadiness and small drag. Received: 8 September 1998 / Accepted: 1 January 2000     

PIV measurements of flow past a confined cylinder

This study investigates the flow past a confined circular cylinder built into a narrow rectangular duct with a Reynolds number range of 1,500 ≤ Re _d ≤ 6,150, by employing the particle image velocimetry technique. In order to better explain the 3-D flow behaviour in the juncture regions of the lower and upper plates and the cylinder, respectively, as well as the dynamics of the horseshoe vortex system, both time-averaged and instantaneous flow data are presented for regions upstream and downstream of the cylinder. The size, intensity and interaction of the vortex systems vary substantially with the Reynolds number. Although the narrow rectangular duct with a single built-in cylinder is a geometrically symmetrical arrrangement, instantaneous flow data have revealed that the flow structures in both the lower and upper plate–cylinder junction regions are not symmetrical with respect to the centreline of the flow passage. The vortical flow structures obtained in side-view planes become dominant sometimes in the lower juncture region and sometimes in the upper juncture region in unsteady mode.     

Role of the shear layer instability in the near wake behavior of two side-by-side circular cylinders

Wakes, and their interaction behind two parallel cylinders lying in a plane perpendicular to the flow, have been investigated experimentally in the sub-critical Reynolds number regime. The experiments were performed in a water channel using laser Doppler velocimetry. The gap between the two cylinders was less than the cylinder diameter, a geometry referred to as strong interaction configuration. In this case the blockage is strong and a gap-jet appears between the cylinders. Two flow regimes of the near wake region have been identified: one below a critical Reynolds number Re _c ]1000;1700[, where the gap jet is stably deflected to one side and the double near-wake becomes asymmetric; the other, above Re _c, where the gap-jet deflection is unstable and a random flopping phenomenon takes place. When Re<Re _c, two different Strouhal numbers are identified, related to the Kármán vortex shedding behind each cylinder. When Re>Re _c, a third frequency appears in the near wake, related to the development of Kelvin-Helmholtz vortices in the separated shear layer of the cylinders [Prasad A, Williamson CHK (1997) J Fluid Mech 333:375]. The observed flopping behavior is attributed to the birth of these Kelvin-Helmholtz instabilities and their intermittent nature. Further downstream, beyond about five cylinder diameters, the random flopping flow phenomena disappear while a slightly asymmetric single wake persists. It is characterized by a Strouhal number St=0.13, a value that one would normally measure behind a single cylinder of twice its diameter.     

Effect of local forcing on a turbulent boundary layer

An experimental study is performed to analyze flow structures behind local suction and blowing in a flat-plate turbulent boundary layer. The local forcing is given to the boundary layer flow by means of a sinusoidally oscillating jet issuing from a thin spanwise slot at the wall. The Reynolds number based on the momentum thickness is about Re _θ =1700. The effects of local forcing are scrutinized by altering the forcing frequency (0.011 ≤ f ⁺≤ 0.044). The forcing amplitude is fixed at A ₀=0.4. It is found that a small local forcing reduces the skin friction and the skin friction reduction increases with the forcing frequency. A phase-averaging technique is employed to capture the large-scale vortex evolution. An organized spanwise vortical structure is generated by the local forcing. The cross-sectional area of vortex and the time fraction of vortex are examined by changing the forcing frequency. An investigation of the random fluctuation components reveals that turbulent energy is concentrated near the center of vortical structures. Received: 17 March 2000/Accepted: 3 April 2001     

Effect of an Exit-Wedge Angle on Pinch-off and Mass Entrainment of Vortex Rings in Air

The effects of exit-wedge angle on evolution, formation, pinch-off, propagation and diffusive mass entrainment of vortex rings in air were studied using digital particle image velocimetry. Vortex rings were generated by passing a solenoid-valve-controlled air jet through a cylindrical nozzle. Experiments were performed over a wide range of exit-wedge angles (10° ≤ α ≤ 90°) of the cylindrical nozzle, initial Reynolds numbers (450 ≤ Re ≤ 4,580) and length-to-diameter ratios (0.9 ≤ L/D ≤ 11) of the air jet. For sharp edges (α ≤ 10°), a secondary ring may emerge at high Reynolds numbers, which tended to distort the vortex ring if ingested into it. For blunt edges (α ≥ 45°), by contrast, stable vortex rings were produced. The formation phase of a vortex ring was found to be closely related to its evolution pattern. An exit-wedge angle of 45° was found to be optimal for rapid pinch-off and faster propagation and better stability of a vortex ring. Diffusive mass entrainment was found to be between 35% and 40% in the early stages of a vortex ring propagation and it gradually increased throughout the course of vortex ring propagation. Entrainment fraction was found to be sensitive to the L/D ratio of the initial jet and decreases when the L/D ratio is increased.