Creation of very-low-Reynolds-number chaotic fluid motions in microchannels using viscoelastic surfactant solution

Solutions of flexible high-molecular-weight polymers or some kinds of surfactant are viscoelastic fluids. The elastic stress is induced in such viscoelastic fluid flows and grows nonlinearly with the flow-rate resulting in many particular flow phenomena, including purely elastic instability. The purely elastic instability can even result in a kind of chaotic fluid motion, the so-called elastic turbulence, which is a recently discovered flow phenomenon and arises at arbitrarily small Reynolds number. By using viscoelastic surfactant solution, we attempted to create the peculiar chaotic fluid motions in several specially designed microchannels in which flows with curvilinear streamlines can be generated. The viscoelastic working fluids were aqueous solutions of surfactant, CTAC/NaSal (cetyltrimethyl ammonium chloride/sodium salicylate). CTAC solutions with weight concentration of 200 ppm (part per million) and 1000 ppm, respectively, at room temperature were tested. For comparison, water flows in the same microchannels were also visualized. The Reynolds numbers for all the microchannel flows were quite small (for solution flows, the Reynolds numbers were the order of or smaller than one) and the flow should be definitely laminar for Newtonian fluid. It was found that the regular laminar flow patterns for low-Reynolds-number Newtonian fluid flow in different microchannels were strongly deformed in solution flows: either asymmetrical flow structures or time-dependent vortical fluid motions appeared. These chaotic flow phenomena were considered to be induced by the viscoelasticity of the CTAC solutions. Discussions about the potential applications using such kind of chaotic fluid motions were also made.     

矩形槽道内表面活性剂减阻流体流场特性

采用粒子图像测速仪对矩形槽道内表面活性减阻流体在流动方向(x方向)与壁面垂直方向(y方向)所在平面的流场进行了测量,分析了速度、涡量、速度脉动相关量在流场内的瞬态分布,以及对500幅相同工况的流场进行了统计平均. 结果显示: 与牛顿流体相比, 表面活性剂减阻流体接近于层流流动,横向速度脉动被大幅减弱,导致湍流输运减弱,雷诺应力远远小于水. 减阻流体流向速度脉动呈条带特征,沿流动方向发展,反映了减阻流体不同于水的湍流输运特征.      

Flow of drag-reducing surfactant solutions in rough pipes

The paper reports the results of experimental study of the flow of hexadecyltrimethylammonium chloride (CTAC) solutions with addition of sodium salicylate (NaSal) in the rough pipes. Measurements were performed in the range of the surfactant concentration from 200 to 400 ppm at a constant molar ratio CTAC/NaSal of 1:2. Five pipes of the relative roughness k/D varying from 1.2 × 10^?2 to 5.6 × 10^?2, obtained by the covering of inner surface of the pipes with glued silicon carbide particles of different size, were studied. The roughness was observed to increase the drag of flow of CTAC/NaSal solutions already at Reynolds numbers higher than 800. With increasing relative roughness k/D, the critical value of Reynolds number, at which the drag reduction disappears, was found to decrease. However, no influence of the roughness on the critical shear stress was noted. The ratio of the critical Reynolds number for rough pipes to that of hydraulically smooth pipes was independent of the surfactant concentration. The degree of drag reduction by the flow of surfactants was greater in rough pipes than in smooth pipes.     

Flow pattern transition in liquid-liquid flows with a transverse cylinder

The effect of a cylindrical bluff body on the interface characteristics of stratified two-phase, oil-water, pipe flows is experimentally investigated with high speed Particle Image Velocimetry (PIV). The motivation was to study the feasibility of flow pattern map actuation by using a transverse cylinder immersed in water in the stratified pattern, and particularly the transition from separated to dispersed flows. The cylinder has a diameter of 5 mm and is located at 6.75 mm from the bottom of the pipe in a 37 mm ID acrylic test section. Velocity profiles were obtained in the middle plane of the pipe. For reference, single phase flows were also investigated for Reynolds numbers from 1550 to 3488. It was found that the flow behind the cylinder was similar to the two dimensional cases, while the presence of the lower pipe wall diverted the vorticity layers towards the top. In two-phase flows, the Froude number (from 1.4 to 1.8) and the depth of the cylinder submergence below the interface affected the generation of waves. For high Froude numbers and low depths of submergence the counter rotating von Karman vortices generated by the cylinder interacted with the interface. In this case, the vorticity clusters from the top of the cylinder were seen to attach at the wave crests. At high depths of submergence, a jet like flow appeared between the top of the cylinder and the interface. High speed imaging revealed that the presence of the cylinder reduced to lower mixture velocities the transition from separated to dual continuous flows where drops of one phase appear into the other.     

Temperature and concentration influence on drag reduction of very low concentrated CTAC/NaSal aqueous solution in turbulent pipe flow

Some additives, like surfactants or polymers, in aqueous solutions lead to dramatic drag reducing effects. Classically the range of concentration of surfactant additive is greater than 500 ppm. It was recently demonstrated that an aqueous solution of CTAC/NaSal (cetyltrimethyl ammonium chloride and sodium salicylate) with a concentration lower than 100 ppm has been found to reduce drag. In this work, the influence of temperature (10–50 °C) and concentration (25–150 ppm) on the drag reduction is studied with an experimental set-up. A cross analysis is made from those results. It reinforces the link between the structure of the micelles (rod-like shape) and the drag reduction rate.     

Analysis of an organized turbulent structure using a pattern recognition technique in a drag-reducing surfactant solution flow

This work presents the investigation for an organized turbulent structure in a drag-reducing flow of dilute surfactant solution by utilizing a particle image velocimetry system to perform the pattern recognition technique on a trajectory in four quadrants of streamwise and wall-normal velocity fluctuations. The pattern recognition is added to a new algorithm in order to directly capture the spatial rotation motion. The Reynolds number based on the channel height and bulk mean velocity was set to 1.5 × 10⁴. Surfactant solution with a weight concentration of 150 ppm was employed and the drag reduction rate was 65%. In the drag-reducing flow, we observe increased frequencies of occurrence of the flow events that correspond to a meandering motion in the wall-normal direction of the high-and low-speed regions. Three findings from investigation of the ensemble-averaged Reynolds shear stress and vortex structure are as follows: (i) the Reynolds shear stress in the large fluctuation range occurs in the narrow region; (ii) Size, strength, arrangement and inclination in the spatial vortex structure in the drag-reducing flow differ from those of the water; and (iii) all trajectory contributions for the wall friction coefficient decrease. Finally, we interpreted that the viscoelasticity characterizing the viscoelastic stress and relaxation time in rheological properties of the flow changes specific elementary vortex for the drag-reducing flow, and the trajectories of each flow pattern change drastically.     

Investigation of turbulence structures in a drag-reduced turbulent channel flow with surfactant additive by stereoscopic particle image velocimetry

In the present study, we employed stereoscopic particle image velocimetry (PIV) to investigate the characteristics of turbulence structures in a drag-reduced turbulent channel flow with addition of surfactant. The tested drag-reducing fluid was a CTAC/NaSal/Water (CTAC: cetyltrimethyl ammonium chloride; NaSal: sodium salicylate) system at 25°C. The weight concentration of CTAC was 30 ppm. Stereoscopic PIV measurement was performed for a water flow (Re=1.1×10⁴) and a CTAC solution flow (Re=1.5×10⁴ with 54% drag reduction) in both the streamwise–spanwise and wall-normal-spanwise planes, respectively. The three-dimensionality of hairpin vortex structures in the near-wall region for wall-bounded turbulent flow was reproduced by conditionally averaging the stereoscopic two-dimensional-three-component velocity fields. A series of wall-normal vortex cores were found to align with the near-wall low-speed streaks with opposite vorticity signals at both sides of the streaks and with the vorticity decreased on average by about one order of magnitude in CTAC solution flow compared with water flow; the spanwise spacing between the near-wall low-speed streaks in the solution flow is increased by about 46%. The streamwise vorticity of the vortex cores appearing in the wall-normal-spanwise plane was also decreased by the use of drag-reducing surfactant additives.     

Instantaneous flow field above the free end of finite-height cylinders and prisms

The flow above the free ends of surface-mounted finite-height circular cylinders and square prisms was studied experimentally using particle image velocimetry (PIV). Cylinders and prisms with aspect ratios of AR = 9, 7, 5, and 3 were tested at a Reynolds number of Re = 4.2 × 10⁴. The bodies were mounted normal to a ground plane and were partially immersed in a turbulent zero-pressure-gradient boundary layer, where the boundary layer thickness relative to the body width was δ/D = 1.6. PIV measurements were made above the free ends of the bodies in a vertical plane aligned with the flow centreline. The present PIV results provide insight into the effects of aspect ratio and body shape on the instantaneous flow field. The recirculation zone under the separated shear layer is larger for the square prism of AR = 3 compared to the more slender prism of AR = 9. Also, for a square prism with low aspect ratio (AR = 3), the influence of the reverse flow over the free end surface becomes more significant compared to that for a higher aspect ratio (AR = 9). For the circular cylinder, a cross-stream vortex forms within the recirculation zone. As the aspect ratio of the cylinder decreases, the reattachment point of the separated flow on the free end surface moves closer to the trailing edge. For both the square prism and circular cylinder cases, the instantaneous velocity vector field and associated in-plane vorticity field revealed small-scale structures mostly generated by the separated shear layer.     

Air–water flow characteristics in high-velocity free-surface flows with 50% void fraction

High-velocity free-surface flows are complex two-phase flows and limited information is available about the interactions between air and water for void fractions of about 50%. Herein a detailed experimental study was conducted in the intermediate flow region (C ∼ 50%) on a stepped spillway and the microscopic air–water flow characteristics were investigated. The results showed differences in water and droplet chord times with comparatively larger number of air chord times (0–2 ms), and larger number of water chord times (2–6 ms). A monotonic decrease of particle chord modes was observed with increasing bubble count rates. Several characteristic time scales were identified based upon inter-particle arrival time analyses of characteristic chord time classes as well as spectral analyses of the instantaneous void fraction signal. Chord times of 3–5 ms appeared to be characteristic time scales of the intermediate flow region having similar time scales compared to the local correlation and integral turbulent time scales and to time scales associated with bubble break-up and turbulent velocity fluctuations. A further characteristic time scale of 100 ms was identified in a frequency analysis of instantaneous void fraction. This time scale was of the same order of magnitude as free-surface auto-correlation time scales suggesting that the air–water flow structure was affected by the free-surface fluctuations.     

Drag reduction for liquid flow through micro-apertures

Pressure drops in the flow through micro-orifices and capillaries were measured for silicone oils, aqueous solutions of polyethylene glycol (PEG), and surfactant aqueous solutions. The diameter of micro-orifices ranged from 5 μm to 400 μm. The corresponding length/diameter ratio was from 4 to 0.05 and capillary diameters were 105 μm and 450 μm. The following results were obtained: silicone oils of 10^?6 m²/s and 10^?5 m²/s in kinematic viscosity generated a reduction of pressure drop (RPD), that is, drag reduction, similar to the RPD of water and a glycerol/water mixture reported in the previous paper by the present authors. When RPD occurred, the pressure drop (PD) of silicone oils of 10^?6 m²/s and 10^?5 m²/s had nearly the same magnitude. Namely, the difference in viscosity did not influence RPD. A 10³ ppm aqueous solution of PEG20000 provided almost the same PD as that of PEG8000 for the 400 μm to 15 μm orifices, but a greater PD than that of PEG8000 for the 10 μm to 5 μm orifices. A non-ionic surfactant and a cationic surfactant were highly effective in RPD compared with anionic surfactants: the non-ionic and cationic surfactant solutions had PD one order of magnitude lower than that of water under some flow conditions in the concentration range from 1 ppm to 10⁴ ppm, but the anionic surfactant solutions did not generate RPD except in the case of the smallest orifice of 5 μm in diameter. The PD of the non-ionic surfactant solution showed a steep rise at a Reynolds number (Re_t) for 400 μm to 15 μm orifices. The Re_t provides the relationship Re_t = K/D, where D is the orifice diameter, and K is a constant of 2 × 10^?2 m for the 100–20 μm orifices irrespective of liquid concentration. Capillary flow experiment revealed that the PEG, non-ionic and cationic surfactant solutions generated RPD also in a laminar flow through the capillary of 105 μm in diameter, but not in the flow through the capillary of 450 μm in diameter. In order to clarify the cause of RPD, an additional experiment was carried out by changing the orifice material from metal to acrylic resin. The result gave a different appearance of RPD, suggesting that RPD is related to an interfacial phenomenon between the liquid and wall. The large RPDs found in the present experiment are very interesting from both academic and practical viewpoints.     

Experimental study on the effect of drag reducing polymer on flow patterns and drag reduction in a horizontal oil–water flow

In this study, a HMW anionic co-polymer of 40:60 wt/wt NaAMPS/acrylamide was used as a drag reducing polymer (DRP) for oil–water flow in a horizontal 25.4 mm ID acrylic pipe. The effect of polymer concentration in the master solution and after injection in the main water stream, oil and water velocities, and pipe length on drag reduction (DR) was investigated. The injected polymer had a noticeable effect on flow patterns and their transitions. Stratified and dual continuous flows extended to higher superficial oil velocities while annular flow changed to dual continuous flow. The results showed that as low as 2 ppm polymer concentration was sufficient to create a significant drag reduction across the pipe. DR was found to increase with polymer concentration increased and reached maximum plateau value at around 10 ppm. The results showed that the drag reduction effect tends to increase as superficial water velocity increased and eventually reached a plateau at U_sw of around 1.3 m/s. At U_sw > 1.0 m/s, the drag reduction decreased as U_so increased while at lower water velocities, drag reduction is fluctuating with respect to U_so. A maximum DR of about 60% was achieved at U_so = 0.14 m/s while only 45% was obtained at U_so = 0.52 m/s. The effectiveness of the DRP was found to be independent of the polymer concentration in the master solution and to some extent pipe length. The friction factor correlation proposed by Al-Sarkhi et al. (2011) for horizontal flow of oil–water using DRPs was found to underpredict the present experimental pressure gradient data.     

Numerical evaluation of turbulence models for dense to dilute gas–solid flows in vertical conveyor

A two-fluid model (TFM) of multiphase flows based on the kinetic theory and small frictional limit boundary condition of granular flow was used to study the behavior of dense to dilute gas–solid flows in vertical pneumatic conveyor. An axisymmetric 2-dimensional, vertical pipe with 5.6 m length and 0.01 m internal diameter was chosen as the computation domain, same to that used for experimentation in the literature. The chosen particles are spherical, of diameter 1.91 mm and density 2500 kg/m³. Turbulence interaction between the gas and particle phases was investigated by Simonin's and Ahmadi's models and their numerical results were validated for dilute to dense conveying of particles. Flow regimes transition and pressure drop were predicted. Voidage and velocity profiles of each phase were calculated in radial direction at different lengths of the conveying pipe. It was found that the voidage has a minimum, and gas and solid velocities have maximum values along the center line of the conveying pipe and pressure drop has a minimum value in transition from dense slugging to dilute stable flow regime. Slug length and pressure fluctuation reduction were predicted with increasing gas velocity, too. It is shown that solid phase turbulence plays a significant role in numerical prediction of hydrodynamics of conveyor and the capability of particles turbulence models depends on tuning parameters of slip-wall boundary condition.     

Influence of wall proximity on characteristics of wake behind a square cylinder: PIV measurements and POD analysis

Influence of wall proximity on characteristics of the wake behind a two-dimensional square cylinder was experimentally studied in the present work. A low-speed recirculation water channel was established for the experiment; the Reynolds number based on the free-stream velocity and cylinder width (D) was kept at Re_D = 2250. Four cases with different gap width, e.g., G/D = 0.1, 0.2, 0.4 and 0.8, were chosen for comparison. Two experimental techniques, e.g., the standard PIV with high image-density CCD camera and TR-PIV with a high-speed camera were employed in measuring the wake field, enabling a comprehensive view of the time-averaged wake pattern at high spatial resolution and the instantaneous flow field at high temporal resolution, respectively. For the four cases, the difference in spatial characteristics of the wake in the vicinity of the plane wall was analyzed in terms of the time-averaged quantities measured by the standard PIV, e.g., the streamline pattern, the vector field, the streamwise velocity fluctuation intensity and the reverse-flow intermittency. The proper orthogonal decomposition (POD) method was extensively used to decompose the TR-PIV measurements, giving a close-up view of the energetic POD modes buried in the wake. The low-order flow model of the wake at G/D = 0.8 and 0.4 was constructed by using the linear combination of the first two POD modes and the time-mean flow field, which reflected well the vortex shedding process in the sense of the phase-dependent patterns. The intermittent appearance of the weakly separated region near the wall was found at G/D = 0.4. On going from G/D = 0.8 to 0.4, the remarkable variation of the instantaneous wake in the longitudinal direction confirmed that the wall constraint stretches the vortices in the plane of the wall and transfers the energy to the longitudinal component at the expense of the lateral one.     

Instabilities in viscoelastic flow through an axisymmetric sudden contraction

Kinematics and dynamics of the viscoelastic flow in an axisymmetric 4 : 1 sudden contraction geometry are studied for a highly elastic polyisobutylene (PIB) based polymer solution (referred to as PIB-Boger fluid). The critical conditions for the onset of the elastic instabilities and the dynamics of the resulting secondary flows are measured for various flow rates. The spatio-temporal characteristics of the flow are determined by instantaneous pressure measurements and streakline photography. The nonlinear dynamics of the global flow field both upstream and downstream of the contraction plane are systematically examined. New dynamic flow behavior and elastic instabilities downstream of the contraction plane are reported. It is shown that the instantaneous pressure measurements along with flow visualization can be used as an effective tool to characterize viscoelastic flows in complex geometries.     

Finite-length effects on dynamical behavior of rod-like particles in wall-bounded turbulent flow

Combined Particle Image Velocimetry (PIV) and Particle Tracking Velocimetry (PTV) measurements have been performed in dilute suspensions of rod-like particles in wall turbulence. PIV results for the turbulence field in the water table flow apparatus compared favorably with data from Direct Numerical Simulations (DNS) of channel flow turbulence and the universality of near-wall turbulence justified comparisons with DNS of fiber-laden channel flow. In order to examine any shape effects on the dynamical behavior of elongated particles in wall-bounded turbulent flow, fibers with three different lengths but the same diameter were used. In the logarithmic part of the wall-layer, the translational fiber velocity was practically unaffected by the fiber length l. In the buffer layer, however, the fiber dynamics turned out to be severely constrained by the distance z to the wall. The short fibers accumulated preferentially in low-speed areas and adhered to the local fluid speed. The longer fibers (l/z > 1) exhibited a bi-modal probability distribution for the fiber velocity, which reflected an almost equal likelihood for a long fiber to reside in an ejection or in a sweep. It was also observed that in the buffer region, high-speed long fibers were almost randomly oriented whereas for all size cases the slowly moving fibers preferentially oriented in the streamwise direction. These phenomena have not been observed in DNS studies of fiber suspension flows and suggested l/z to be an essential parameter in a new generation of wall-collision models to be used in numerical studies.     

Flow evolution of a turbulent submerged two-dimensional rectangular free jet of air. Average Particle Image Velocimetry (PIV) visualizations and measurements

The paper presents average flow visualizations and measurements, obtained with the Particle Image Velocimetry (PIV) technique, of a submerged rectangular free jet of air in the range of Reynolds numbers from Re = 35,300 to Re = 2200, where the Reynolds number is defined according to the hydraulic diameter of a rectangular slot of height H. According to the literature, just after the exit of the jet there is a zone of flow, called zone of flow establishment, containing the region of mixing fluid, at the border with the stagnant fluid, and the potential core, where velocity on the centerline maintains a value almost equal to the exit one. After this zone is present the zone of established flow or fully developed region. The goal of the paper is to show, with average PIV visualizations and measurements, that, before the zone of flow establishment is present a region of flow, never mentioned by the literature and called undisturbed region of flow, with a length, L_U, which decreases with the increase of the Reynolds number. The main characteristics of the undisturbed region is the fact that the velocity profile maintains almost equal to the exit one, and can also be identified by a constant height of the average PIV visualizations, with length, L_CH, or by a constant turbulence on the centerline, with length L_CT. The average PIV velocity and turbulence measurements are compared to those performed with the Hot Film Anemometry (HFA) technique. The average PIV visualizations show that the region of constant height has a length L_CH which increases from L_CH = H at Re = 35,300 to L_CH = 4–5H at Re = 2200. The PIV measurements on the centerline of the jet show that turbulence remains constant at the level of the exit for a length, L_CT, which increases from L_CT = H at Re = 35,300 to L_CT = 4–5H at Re = 2200. The PIV measurements show that velocity remains constant at the exit level for a length, L_U, which increases from L_U = H at Re = 35,300 to L_U = 6H at Re = 2200 and is called undisturbed region of flow. In turbulent flow the length L_U is almost equal to the lengths of the regions of constant height, L_CH, and constant turbulence, L_CT. In laminar flow, Re = 2200, the length of the undisturbed region of flow, L_U, is greater than the lengths of the regions of constant height and turbulence, L_CT = L_CH = 4–5H. The average PIV and HFA velocity measurements confirm that the length of potential core, L_P, increases from L_P = 4–5H at Re = 35,300 to L_P = 7–8H at Re = 2200, and are compared to the previous experimental and theoretical results of the literature in the zone of mixing fluid and in the fully developed region with a good agreement.     

Technique for quantitative mapping of three-dimensional liquid–gas phase boundaries in microchannel flows

A diagnostic technique capable of characterizing interfaces between transparent, immiscible fluids is developed and demonstrated by investigating the morphology of liquid–gas interfaces in an adiabatic two-phase flow through a microchannel of 500 μm × 500 μm square cross section. Water seeded with 0.5 μm-diameter fluorescent polystyrene particles is pumped through the channel, and the desired adiabatic two-phase flow regime is achieved through controlled air injection. The diagnostic technique relies on obtaining particle position data through epifluorescent imaging of the flow at excitation and emission wavelengths of 532 nm and 620 nm, respectively. The particle position data are then used to resolve interface locations to within ±1 μm in the focal plane. By mapping the interface within individual focal planes at various depths within the channel, it is possible to determine the complete liquid–gas interface geometry across the channel cross section in a dynamic flow environment. Utilizing this approach, the liquid–gas phase boundaries of annular flows within a microchannel have been successfully characterized.     

Large Eddy Simulations of turbulent particle-laden channel flow

This paper scrutinises the Large Eddy Simulation (LES) approach to simulate the behaviour of inter-acting particles in a turbulent channel flow. A series of simulations that are fully (four-way), two-way and one-way coupled are performed in order to investigate the importance of the individual physical phenomena occurring in particle-laden flows. Moreover, the soft sphere and hard sphere models, which describe the interaction between colliding particles, are compared with each other and the drawbacks and advantages of each algorithm are discussed. Different models to describe the sub-grid scale stresses with LES are compared. Finally, simulations accounting for the rough walls of the channel are compared to simulations with smooth walls. The results of the simulations are discussed with the aid of the experimental data of Kussin J. and Sommerfeld M., 2002, Experimental studies on particle behaviour and turbulence modification in horizontal channel flow with different wall roughness, Exp. in Fluids, 33, pp. 143–159 of Reynolds number 42,000 based on the full channel height. The simulations are carried out in a three-dimensional domain of 0.175 m × 0.035 m  × 0.035 m where the direction of gravity is perpendicular to the flow. The simulation results demonstrate that rough walls and inter-particle collisions have an important effect in redistributing the particles across the channel, even for very dilute flows. A new roughness model is proposed which takes into account the fact that a collision in the soft sphere model is fully resolved and it is shown that the new model is in very good agreement with the available experimental data.     

Influence of aspect ratio on the flow above the free end of a surface-mounted finite cylinder

The flow above the free end of a surface-mounted finite-height cylinder was studied in a low-speed wind tunnel using particle image velocimetry (PIV). Velocity measurements were made in vertical and horizontal measurement planes above the free end of finite cylinders of aspect ratios AR = 9, 7, 5 and 3, at a Reynolds number of Re = 4.2 × 10⁴. The relative thickness of the boundary layer on the ground plane was δ/D = 1.7. Flow separating from the leading edge formed a prominent recirculation zone on the free-end surface. The legs of the mean arch vortex contained within the recirculation zone terminate on the free-end surface on either side of the centreline. Separated flow from the leading edge attaches onto the upper surface of the cylinder along a prominent attachment line. Local separation downstream of the leading edge is also induced by the reverse flow and arch vortex circulation within the recirculation zone. As the cylinder aspect ratio is lowered from AR = 9 to AR = 3, the thickness of the recirculation zone increases, the arch vortex centre moves downstream and higher above the free-end surface, the attachment position moves downstream, and the termination points of the arch vortex move upstream. A lowering of the aspect ratio therefore results in accentuated curvature of the arch vortex line. Changes in aspect ratio also influence the vorticity generation in the near-wake region and the shape of the attachment line.     

PIV study of large-scale flow organisation in slot jets

The paper reports on particle image velocimetry (PIV) measurements in turbulent slot jets bounded by two solid walls with the separation distance smaller than the jet width (5–40%). In the far-field such jets are known to manifest features of quasi-two dimensional, two component turbulence. Stereoscopic and tomographic PIV systems were used to analyse local flows. Proper orthogonal decomposition (POD) was applied to extract coherent modes of the velocity fluctuations. The measurements were performed both in the initial region close to the nozzle exit and in the far fields of the developed turbulent slot jets for Re ⩾ 10,000. A POD analysis in the initial region indicates a correlation between quasi-2D vortices rolled-up in the shear layer and local flows in cross-stream planes. While the near-field turbulence shows full 3D features, the wall-normal velocity fluctuations day out gradually due to strong wall-damping resulting in an almost two-component turbulence. On the other hand, the longitudinal vortex rolls take over to act as the main agents in wall-normal and spanwise mixing and momentum transfer. The quantitative analysis indicates that the jet meandering amplitude was aperiodically modulated when arrangement of the large-scale quasi-2D vortices changed between asymmetric and symmetric pattern relatively to the jet axis. The paper shows that the dynamics of turbulent slot jets are more complex than those of 2D, plane and rectangular 3D jets. In particular, the detected secondary longitudinal vortex filaments and meandering modulation is expected to be important for turbulent transport and mixing in slot jets. This issue requires further investigations.