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.     

The effect of incidence angle on the mean wake of surface-mounted finite-height square prisms

The mean wake of a surface-mounted finite-height square prism was studied experimentally in a low-speed wind tunnel to explore the combined effects of incidence angle (α) and aspect ratio (AR). Measurements of the mean wake velocity field were made with a seven-hole pressure probe for finite square prisms of AR = 9, 7, 5 and 3, at a Reynolds number of Re = 3.7 × 10⁴, for incidence angles from α = 0° to 45°. The relative thickness of the boundary layer on the ground plane, compared to the prism width, was δ/D = 1.5. As the incidence angle increases from α = 0° to 15°, the mean recirculation zone shortens and the mean wake shifts in the direction opposite to that of the mean lift force. The downwash is also deflected to this side of the wake and the mean streamwise vortex structures in the upper part of the wake become strongly asymmetric. The shortest mean recirculation zone, and the greatest asymmetry in the mean wake, is found at the critical incidence angle of α_critical ≈ 15°. As the incidence angle increases from α = 15° to 45°, the mean recirculation zone lengthens and the mean streamwise vortex structures regain their symmetry. These vortices also elongate in the wall-normal direction and become contiguous with the horseshoe vortex trailing arms. The mean wake of the prism of AR = 3 has some differences, such as an absence of induced streamwise vorticity near the ground plane, which support its classification as lying below the critical aspect ratio for the present flow conditions.     

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.     

Influence of rounding corners on the wake of a finite-length cylinder: An experimental study

This study aims to investigate experimentally the influence of rounding corners (r) as well as aspect ratio (AR) on the flow structures of a surface-mounted finite cylinder. The cylinders with sharp (r* = r/D = 0) and rounded corners (r*=0.167, 0.25 and 0.5) and aspect ratio or height-to-width/diameter ratio (AR = H/D) between 2 and 7 are utilized. The experiments are based on the five-hole probe and hot-wire measurements as well as the oil flow visualization. Wake measurements are made in an open return wind tunnel at the Reynolds number, Re = 1.6 × 10⁴, where Re is defined based on the side width/diameter (D) of the cylinder cross-section and the freestream velocity. It is found that r* and AR have significant effects on the flow structure from the perspective of wake topology, strength of streamwise vortices, and vortex shedding frequency. For all r* considered, the wake is characterized by a quadrupole type (both the tip and base vortices are present) at AR = 7, while a dipole type occurs for AR = 2 and 4 (the base vortices are absent). The strength (circulation) of the streamwise vortex structures is affected by r*. For all AR examined in the present study, the strengths of tip and base vortex structures decrease with increasing r*. The oil flow visualization demonstrates that the features of the horseshoe vortex are sensitive to r* and AR. With increasing r*, the location of the separation line moves downstream and the distance between horseshoe vortex legs decreases. Velocity measurements reveal that the downwash flow enhances with increasing r*. It is also found that the Strouhal number increases progressively by 60% as r* increases from 0 to 0.5, regardless of AR.     

The effect of a splitter plate on the flow around a finite prism

The effect of a wake-mounted splitter plate on the flow around a surface-mounted finite-height square prism was investigated experimentally in a low-speed wind tunnel. Measurements of the mean drag force and vortex shedding frequency were made at Re=7.4×10⁴ for square prisms of aspect ratios AR=9, 7, 5 and 3. Measurements of the mean wake velocity field were made with a seven-hole pressure probe at Re=3.7×10⁴ for square prisms of AR=9 and 5. The relative thickness of the boundary layer on the ground plane was δ/D=1.5–1.6 (where D is the side length of the prism). The splitter plates were mounted vertically from the ground plane on the wake centreline, with a negligible gap between the leading edge of the plate and rear of the prism. The splitter plate heights were always the same as the heights of prisms, while the splitter plate lengths ranged from L/D=1 to 7. Compared to previously published results for an “infinite” square prism, a splitter plate is less effective at drag reduction, but more effective at vortex shedding suppression, when used with a finite-height square prism. Significant reduction in drag was realized only for short prisms (of AR≤5) when long splitter plates (of L/D≥5) were used. In contrast, a splitter plate of length L/D=3 was sufficient to suppress vortex shedding for all aspect ratios tested. Compared to previous results for finite-height circular cylinders, finite-height square prisms typically need longer splitter plates for vortex shedding suppression. The effect of the splitter plate on the mean wake was to narrow the wake width close to the ground plane, stretch and weaken the streamwise vortex structures, and increase the lateral entrainment of ambient fluid towards the wake centreline. The splitter plate has little effect on the mean downwash. Long splitter plates resulted in the formation of additional streamwise vortex structures in the upper part of the wake.     

Separated flow structures around a cylindrical obstacle in a narrow channel

A detailed experimental study is performed on the separated flow structures around a low aspect-ratio circular cylinder (pin-fin) in a practical configuration of liquid cooling channel. Distinctive features of the present arrangement are the confinement of the cylinder at both ends, water flow at low Reynolds numbers (Re = 800, 1800, 2800), very high core flow turbulence and undeveloped boundary layers at the position of the obstacle. The horseshoe vortex system at the junctions between the cylinder and the confining walls and the near wake region behind the obstacle are deeply investigated by means of Particle Image Velocimetry (PIV). Upstream of the cylinder, the horseshoe vortex system turns out to be perturbed by vorticity bursts from the incoming boundary layers, leading to aperiodical vortex oscillations at Re = 800 or to break-away and secondary vorticity eruptions at the higher Reynolds numbers. The flow structures in the near wake show a complex three-dimensional behaviour associated with a peculiar mechanism of spanwise mass transport. High levels of free-stream turbulence trigger an early instabilization of the shear layers and strong Bloor–Gerrard vortices are observed even at Re = 800. Coalescence of these vortices and intense spanwise flow inhibit the alternate primary vortex shedding for time periods whose length and frequency increase as the Reynolds number is reduced. The inhibition of alternate vortex shedding for long time periods is finally related to the very large wake characteristic lengths and to the low velocity fluctuations observed especially at the lowest Reynolds number.     

Aerodynamic forces and three-dimensional flow structures in the mean wake of a surface-mounted finite-height square prism

Different flow models have been proposed for the flow around surface-mounted finite-height square prisms, but there is still a lack of consensus about the origin and connection of the streamwise tip vortices with the other elements of the wake. This numerical study was performed to address this gap, in addition to clarifying the relationship of the near-wake structures with the far wake and the near-wall flow, which is associated with the fluid forces. A large-eddy simulation approach was adopted to solve the flow around a surface-mounted finite-height square prism with an aspect ratio of AR = 3 and a Reynolds number Re = 500. The mean drag and normal forces and the bending moment for the prism were quantitatively compared in terms of skin-friction and pressure contributions, and related to the near-wall flow. Both three-dimensional visualizations and planar projections of the time-averaged flow field were used to identify, qualitatively, the main structures of the wake, including the horseshoe vortex, corner vortices and regions of high streamwise vorticity in the upper part of the wake. These features showed the same qualitative behavior as reported in high Reynolds number studies. It was found that some regions of high streamwise vorticity magnitude, like the tip vortices, are associated with the three-dimensional bending of the flow, and the tip vortices did not continuously extend to the free end of the prism. The three-dimensional flow analysis, which integrated different observations of the flow field around surface-mounted finite-height square prisms, also revealed that the mean near-wake structure is composed of two sections of different origin and location of dominance.     

Tip gap height effects on flow structure and heat/mass transfer over plane tip of a high-turning turbine rotor blade

The effects of tip gap height-to-chord ratio, h/c, on the flow structure and heat/mass transfer over the plane tip surface of a large-scale high-turning turbine rotor blade have been investigated for h/c = 1.0%, 2.0%, 3.0% and 4.0%. For near-wall tip gap flow visualizations, a high-resolution oil film method is employed, and the naphthalene sublimation technique is used for local heat/mass transfer rate measurements. From the tip surface visualizations, a pair of vortices named “tip gap vortices” is identified in the leading edge region within the tip gap. The overall tip gap flow is characterized not only by the tip gap vortices but also by the flow separation/re-attachment process along the pressure-side tip edge. Within the separation bubble, there exist complicated near-wall flows moving toward a mid-chord flow converging area. With increasing h/c, the tip gap vortices, the flow separation/re-attachment, and the converging flows within the separation bubble tend to be intensified. In general, higher thermal load is found along the loci of the tip gap vortices and along the re-attachment line, while lower thermal load is observed behind the tip gap vortex system and near the mid-chord flow converging area. Heat/mass transfer characteristics with the variation of h/c are discussed in detail in conjunction with the tip gap flow features. Based on the flow visualizations and heat/mass transfer data, new realistic tip gap flow models have been proposed for h/c = 1.0 and 4.0%.     

Effect of side ratio on fluid flow and heat transfer from rectangular cylinders using the PANS method

A computational study of heat transfer from rectangular cylinders is carried out. Rectangular cylinders are distinguished based on the ratio of the length of streamwise face to the height of the cross-stream face (side ratio, R). The simulations were performed to understand the heat transfer in a flow field comprising separation, reattachment, vortex shedding and stagnation. The Partially-Averaged Navier–Stokes (PANS) modeling approach is used to solve the turbulent flow physics associated and the wall resolve approach is used for the near wall treatment because of the flow separation involved. The simulations were performed using a finite volume based opensource software, OpenFOAM, at Reynolds number (Re) = 22,000 for rectangular cylinder at constant temperature kept in an air stream. Two critical side ratios were obtained, R = 0.62 and 3.0. At R = 0.62, the maximum value of the drag coefficient (C_d) = 2.681 was observed which gradually reduced by 54% at R = 4.0. The base pressure coefficient and global Nusselt number also attained the maximum value at R = 0.62 and from R = 2.5 to 3.0 a sharp discontinuous increase by 140% in the Strouhal number was observed. At R = 0.62, it was observed that the separated flow reattaches at the trailing edge after rolling over the side face and therefore increases the overall Nusselt number. The phase averaging was also performed to analyze the unsteady behavior of heat transfer.     

Investigation of flow field of clap and fling motion using immersed boundary coupled lattice Boltzmann method

This paper deals with the investigation of flow field due to clap and fling mechanism using immersed boundary coupled with lattice Boltzmann method. The lattice Boltzmann method (LBM), an alternative to Navier–Stokes solver, is used because of its simplicity and computational efficiency in solving complex moving boundary problems. Benchmark problems are simulated to validate the code, which is then used for simulating flow over two elliptic wing of aspect ratio 5 performing clap and fling flapping motion for different flow parameters such as Reynolds number (Re=75, 100, 150), advance ratio (J=10E−3.10E−2,0.2) and frequency (f=0.05 Hz, 0.25 Hz). Numerical simulation is able to capture typical low Reynolds number unsteady phenomena such as, ׳wake vortex wing interaction׳, ׳Kramer effect׳ and ׳delayed stall׳. The results are both qualitatively and quantitatively consistent with experimental observation. The parametric study involving different combinations of Re, f and J depict distinctly different aerodynamic performances providing physical insights into the flow physics. It is observed that a combination of low f, low J and high Re flow results in better aerodynamic performance. Pronounced lift enhancement via leading edge vortices are obtained in unsteady regime (J<1) compared to quasi-steady regime (J>1). The role of leading edge vortices in enhancing lift are investigated by studying the size and strength of these vortices for different flow conditions. For a given Re, the magnitude of maximum lift coefficient decreases with increasing f irrespective of the value of J; while the same is enhanced with the increasing Re.     

URANS vs. experiments of flow induced motions of multiple circular cylinders with passive turbulence control

Two-dimensional Unsteady Reynolds-Average Navier–Stokes equations with the Spalart–Allmaras turbulence model are used to simulate the flow induced motions of multiple circular cylinders with passive turbulence control (PTC) in steady uniform flow. Four configurations with 1, 2, 3, and 4 cylinders in tandem are simulated and studied at a series of Reynolds numbers in the range of 30 000<Re<120 000. Simulation results are verified by experimental data measured in the Marine Renewable Energy Laboratory. Good agreement was observed between the values of vorticity, amplitude ratio, and frequency ratio predicted by numerical simulations and experimental measurements. The amplitude and frequency response show the initial and upper branches in vortex induced vibration (VIV), transition from VIV to galloping, and galloping branch for all PTC-cylinders. The maximum amplitude of 2.9 diameters for the first cylinder is achieved at Re=104 356 in the numerical results. Compared with the first cylinder, the VIV initial branch starts at higher Re for the downstream cylinders due to the presence of the upstream cylinder(s). 2P and 2P+2S vortex patterns are observed at Re=62 049 and Re=90 254 for the single PTC-cylinder. Furthermore, the shed vortices of the downstream cylinders are strongly disrupted and modified by the vortices shed from the upstream one in the cases of multiple PTC-cylinders.     

The flow field around a micropillar confined in a microchannel

The flow field over a low aspect ratio (AR) circular pillar (L/D = 1.5) in a microchannel was studied experimentally. Microparticle image velocimetry (μPIV) was employed to quantify flow parameters such as flow field, spanwise vorticity, and turbulent kinetic energy (TKE) in the microchannel. Flow regimes of cylinder-diameter-based Reynolds number at 100 ⩽ Re_D ⩽ 700 (i.e., steady, transition from quasi-steady to unsteady, and unsteady flow) were elucidated at the microscale. In addition, active flow control (AFC), via a steady control jet (issued from the pillar itself in the downstream direction), was implemented to induce favorable disturbances to the flow in order to alter the flow field, promote turbulence, and increase mixing. Together with passive flow control (i.e., a circular pillar), turbulent kinetic energy was significantly increased in a controllable manner throughout the flow field.     

Simulation of a valveless pump with an elastic tube

A valveless pump consisting of a pumping chamber with an elastic tube was simulated using an immersed boundary (IB) method. The interaction between the motion of the elastic tube and the pumping chamber generated a net flow toward the outlet throughout a full cycle of the pump. The net flow rate of the valveless pump was examined by varying the stretching coefficient (ϕ), bending coefficient (γ), the aspect ratio (l/d) of the elastic tube, and the frequency (f) of the pumping chamber. As the stretching and bending coefficients of the elastic tube increased, the net flow through the valveless pump decreased. Elastic tubes with aspect ratios in the range of 2 ⩽ l/d ⩽ 3 generated a higher flow rate than that generated for tubes with aspect rations of l/d = 1 or 4. As the frequency of the pumping chamber increased, the net flow rate of the pump for l/d = 2 increased. However, the net flow rate for l/d = 3 was nonlinearly related to the pumping frequency due to the complexity of the wave motions. Snapshots of the fluid velocity vectors and the wave motions of the elastic tube were examined over one cycle of the pump to gain a better understanding of the mechanism underlying the valveless pump. The relationship between the average gap in the elastic tube and the average flow rate of the pump was analyzed. A smaller gap in the elastic tube during the expansion mode and a wider gap in the elastic tube during the contraction mode played a dominant role in generating a high average flow rate in the pump, regardless of the stretching coefficient (ϕ), the aspect ratio (l/d) of the elastic tube, or the pumping frequency of the pumping chamber (f).     

Experimental study on the effect of pulsating inflow to an enclosure for improved mixing

Optimal control of inlet jet flows is of broad interest for enhanced mixing in ventilated rooms. The general approach in mechanical ventilation is forced convection by means of a constant flow rate supply. However, this type of ventilation may cause several problems such as draught and appearance of stagnation zones, which reduces the ventilation efficiency. A potential way to improve the ventilation quality is to apply a pulsating inflow, which has been hypothesised to reduce the stagnation zones due to enhanced mixing. The present study aims at testing this hypothesis, experimentally, in a small-scale two-dimensional water model using Particle Image Velocimetry with an in-house vortex detection program. We are able to show that for an increase in pulsation frequency or alternatively in the flow rate the stagnation zones are reduced in size and the distribution of vortices becomes more homogeneous over the considered domain. The number of vortices (all scales) increases by a factor of four and the swirl-strength by about 50% simply by turning on the inflow pulsation. Furthermore, the vortices are well balanced in terms of their rotational direction, which is validated by the symmetric Probability Density Functions of vortex circulation (Γ) around Γ = 0. There are two dominating vortex length scales in the flow, namely 0.6 and 0.8 inlet diameters and the spectrum of vortex diameters become broader by turning on the inflow pulsation. We conclude that the positive effect for enhanced mixing by increasing the flow rate can equally be accomplished by applying a pulsating inflow.     

Direct numerical simulation of a NACA0012 in full stall

This work aims at investigating the mechanisms of separation and the transition to turbulence in the separated shear-layer of aerodynamic profiles, while at the same time to gain insight into coherent structures formed in the separated zone at low-to-moderate Reynolds numbers. To do this, direct numerical simulations of the flow past a NACA0012 airfoil at Reynolds numbers Re = 50,000 (based on the free-stream velocity and the airfoil chord) and angles of attack AOA = 9.25° and AOA = 12° have been carried out. At low-to-moderate Reynolds numbers, NACA0012 exhibits a combination of leading-edge/trailing-edge stall which causes the massive separation of the flow on the suction side of the airfoil. The initially laminar shear layer undergoes transition to turbulence and vortices formed are shed forming a von Kármán like vortex street in the airfoil wake. The main characteristics of this flow together with its main features, including power spectra of a set of selected monitoring probes at different positions on the suction side and in the wake of the airfoil are provided and discussed in detail.     

Impact of aspect ratio on flow boiling of water in rectangular microchannels

In this paper we focus on the impact of varying the aspect ratio of rectangular microchannels, on the overall pressure drop involving water boiling. An integrated system comprising micro-heaters, sensors and microchannels has been realized on (1 1 0) silicon wafers, following CMOS compatible process steps. Rectangular microchannels were fabricated with varying aspect ratios (width [W] to depth [H]) but constant hydraulic diameter of 142 ± 2 μm and length of 20 mm. The invariant nature of the hydraulic diameter is confirmed through two independent means: physical measurements using profilometer and by measuring the pressure drop in single-phase fluid flow. The experimental results show that the pressure drop for two-phase flow in rectangular microchannels experiences minima at an aspect ratio of about 1.6. The minimum is possibly due to opposing trends of frictional and acceleration pressure drops, with respect to aspect ratio. In a certain heat flux and mass flux range, it is observed that the two-phase pressure drop is lower than the corresponding single-phase value. This is the first study to investigate the effect of aspect ratio in two-phase flow in microchannels, to the best of our knowledge. The results are in qualitative agreement with annular flow model predictions. These results improve the possibility of designing effective heat-sinks based on two-phase fluid flow in microchannels.     

LES of turbulent jet in cross-flow: Part 1 – A numerical validation study

The paper presents results of a LES based numerical simulation of the turbulent jet-in-cross-flow (JICF) flowfield, with Reynolds number based on cross-flow velocity and jet diameter Re = 2400 and jet-to-cross-flow velocity ratio of R = 3.3. The JICF flow case has been investigated in great detail, involving conduction of two independent precursor simulations, prior to the main JICF simulation, as the considered case has turbulent inflow conditions on both jet and cross-stream side. The LES results are directly compared to pointwise Laser Doppler Anemometry (LDA) measurements, showing a very good agreement on the level of various statistical quantities in all flow regions but the immediate jet-to-cross-flow exhaustion zone. Several LES computations involving grids of up to 15 million grid points have been conducted, showing no improvement in the agreement between numerical results and measurements, possibly indicating a LDA measurement problem in this particular region.     

Flow pattern and molecular visualization of DNA solutions through a 4:1 planar micro-contraction

The present work explores unusual flow behavior of entangled fluids in an abrupt contraction flow device. Fluorescent imaging was carried out on four different entangled DNA solutions with concentrations ranging from 0.1 to 1.0% (with a wide range of entanglements per chain Z = 7–55). For weakly entangled solutions (Z < 30), vortex flow was dominant at high flow rates. However, for well-entangled DNA solutions (Z ≥ 30), unusual time dependant shear banding was observed at the contraction entrance. Upon reducing the slip length by adding sucrose to the well-entangled DNA solution, vortex flow became dominant again. In vortex flow, most DNA chains remained coiled at the corner in regular recirculation. However, when jerky-shear-banding flow developed, significant stable stretching of DNA chains occurred at the center-line, with quasi-periodic switching between stretching and recoil at the corner.     

Experimental study on turbulent mixing of spray droplets in crossflow

Centrifugal spray injected at various angles in gas crossflow has been studied experimentally using PIV visualization system and image-processing techniques. Experiments were carried out inside a rectangular duct (95 mm × 95 mm in cross-section) at ambient temperature and pressure, with different gas Reynolds numbers (vary from 12,900 to 45,000) and three injection angles (60°, 90° and 120°). The spray angle of the centrifugal nozzle is 80°, with D₃₂ of 80 μm. The instantaneous images of droplets distribution and the values of the line-averaged D₃₂ at different positions on the cross-sections along the flow field for each condition were obtained, and their flow field configurations were achieved. Quantitative assessments of mixing degree between two phases for different injection angles were determined using a spatial unmixedness parameter. It is found that the addition of droplets into the gas crossflow enhanced the turbulence intensity of the gas crossflow and caused different-scale vortices. The flow field structure, to a great extent, is dependent on the injection angle. The entrainment and centrifugal force of large vortex lead to uneven droplet distribution and moreover influence the mixing of droplets and gas crossflow. A better mixing result can be obtained with the injection angle of nozzles of 60°.