Low Reynolds number flow over a square cylinder with a detached flat plate

A numerical study of the alteration of a square cylinder wake using a detached downstream thin flat plate is presented. The wake is generated by a uniform flow of Reynolds number 150 based on the side length of the cylinder, D. The sensitivity of the near wake structure to the downstream position of the plate is investigated by varying the gap distance (G) along the wake centerline in the range 0 ⩽ G ⩽ 7D for a constant plate length of L = D. A critical gap distance is observed to occur at G_c ≈ 2.3D that indicates the existence of two flow regimes. Regime I is characterised by vortex formation occurring downstream of the gap while for regime II, formation occurs within the gap. By varying the plate length and gap distance, a condition is found where significant unsteady total lift reduction can occur. The root mean square lift reduction is limited by an unsteady stall process on the plate.     

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.     

Comparison of flow structures in the downstream region of a cylinder and sphere

An experimental investigation of flow structures downstream of a circular cylinder and sphere immersed in a free-stream flow is performed for Re = 5000 and 10,000 using qualitative and quantitative flow visualization techniques. The obtained results are presented in terms of time-averaged velocity vectors, patterns of streamlines, vorticity, Reynolds stress correlations and turbulent kinetic energy distributions. Flow data reveal that the size of wake flow region, the location of singular and double points, the peak values of turbulence quantities, such as Reynolds stress correlations, vorticity fluctuations and turbulent kinetic energy vary as a function of models’ geometry and Reynolds Numbers. The concentration of small scale vortices is more dominant in the wake of the sphere than that of the cylinder. The maximum value of turbulent kinetic energy (TKE) occurs close to the saddle point for the cylinder case while two maximum values of TKE occur along shear layers for the sphere one because of the 3-D flow behavior.     

Comparison of flow structures in the downstream region of a cylinder and sphere

An experimental investigation of flow structures downstream of a circular cylinder and sphere immersed in a free-stream flow is performed for Re = 5000 and 10,000 using qualitative and quantitative flow visualization techniques. The obtained results are presented in terms of time-averaged velocity vectors, patterns of streamlines, vorticity, Reynolds stress correlations and turbulent kinetic energy distributions. Flow data reveal that the size of wake flow region, the location of singular and double points, the peak values of turbulence quantities, such as Reynolds stress correlations, vorticity fluctuations and turbulent kinetic energy vary as a function of models’ geometry and Reynolds Numbers. The concentration of small scale vortices is more dominant in the wake of the sphere than that of the cylinder. The maximum value of turbulent kinetic energy (TKE) occurs close to the saddle point for the cylinder case while two maximum values of TKE occur along shear layers for the sphere one because of the 3-D flow behavior.     

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.     

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.     

The influence of pipe length on thermal statistics computed from DNS of turbulent heat transfer

We present results from direct numerical simulation of turbulent heat transfer in pipe flow at a bulk flow Reynolds number of 5000 and Prandtl numbers ranging from 0.025 to 2.0 in order to examine the effect of streamwise pipe length (πδ ≡ πD/2 ⩽ L ⩽ 12πδ) on the convergence of thermal turbulence statistics. Various lower and higher order thermal statistics such as mean temperature, rms of fluctuating temperature, turbulent heat fluxes, two-point auto and cross-correlations, skewness and flatness were computed and it is found that the value of L required for convergence of the statistics depends on the Prandtl number: larger Prandtl numbers requires comparatively shorter pipe length for convergence of most of the thermal statistics.     

Turbulent flow past a bubble and an ellipsoid using shadow-image and PIV techniques

An experimental investigation on flow around an oscillating bubble and solid ellipsoid with a flat bottom was conducted. A single air bubble (equivalent diameter D_e=9.12 mm) was attached to a small disk (∼1 mm) at the end of a needle and suspended across a vertical square channel (100 mm) by wire wherein water flowed downward at a constant flowrate. The solid ellipsoid (D_e∼9.1 mm) was suspended across the square channel in the same manner. The equivalent diameter-based Reynolds and Eotvos number range, 1950<Re<2250 and 11<Eo<11.5, placed the bubble in the ‘wobbly’ regime while the flow in its wake was turbulent. A constant flowrate and one bubble size was used such that flow in the wake was turbulent. Velocity measurements of the flow field around the bubble or solid were made using a one CCD camera Digital Particle Image Velocimetry (DPIV) system enhanced by Laser Induced Fluorescence (LIF). The shape of the bubble or solid was simultaneously recorded along with the velocity using a second CCD camera and an Infrared Shadow Technique (IST). In this way both the flow-field and the boundary of the bubble (solid) were measured. The velocity vector plots of flow around and in the wake of a bubble/solid, supplemented by profiles and contours of the average and root-mean-square velocities, vorticity, Reynolds stress and turbulent kinetic energy, revealed differences in the wake flow structure behind a bubble and solid. One of the significant differences was in the inherent, oscillatory motion of the bubble which not only produced vorticity in the near-wake, but as a result of apparent vorticity stretching distributed the turbulent kinetic energy associated with this flow more uniformly on its wake, in contrast to the solid.     

Dynamic interaction of fixed dual spheres for several configurations and inflow conditions

The changes in force characteristics as well as the shedding patterns for various dual sphere configurations are studied. The Reynolds numbers considered are 300, 600 and two different inflow conditions are used: steady and pulsating. The sphere formations are defined by the separation distance D₀ between the spheres and the angle between the line connecting the centres of the spheres and the main flow direction, γ. The position of one of the spheres is varied in the range 0°–90° using a 15° increment. Two separation distances are studied; 1.5D and 3D. The method used for the simulations is the Volume of Solid (VOS) approach, a method based on Volume of Fluid (VOF). A major conclusion from this work is that the sphere interaction alters the wake dynamics by obstructing the vortex shedding (generating a steady wake or a wake with lower Strouhal number) and by changing the direction of the lift force so that it in most cases is directed in the plane containing the sphere centres. The results also show that changing the inflow condition gives the same relative change in drag and lift as for a single sphere. The drag is substantially reduced by placing the sphere downstream in a tandem arrangement and slightly increased in a side-by-side arrangement. However, the effect is decreased by increasing separation distance and increasing Reynolds number.     

Hydrodynamic forces on a rotating sphere

The wake dynamics of a rotating sphere with prescribed rotation axis angles are quantitatively analysed by carrying out numerical simulations at Reynolds numbers of Re = 100, 250 and 300, non-dimensional rotational rates Ω¹ = 0–1 and rotation axis angles α = 0, π/6, π/3 and π/2 measured from the free stream axis. These parameters are the same as those in an earlier study (Poon et al., 2010, Int. J. Heat Fluid Flow) where the instantaneous flow structures were discussed qualitatively. This study extends the findings of the earlier study by employing phase diagrams (C_Lx, C_Ly) and (C_D, C_L) to provide a quantitative analysis of the time-dependent behaviour of the flow structures. At Re = 300 and Ω¹ = 0.05, the phase diagrams (C_Lx, C_Ly) show ‘saw tooth’ patterns for both α = 0 and π/6. The ‘saw tooth’ pattern indicates that the flow structures comprise a higher frequency oscillation component at a Reynolds number of 300 which is not observed until Re ≈ 800 for a stationary sphere. This ‘saw tooth’ pattern disappears as Ω¹ increases. The employment of the phase diagrams also reveals that different flow structures induce different oscillation amplitudes on both lateral force coefficients. With the exception of the vortices formed from a shear layer instability, all other flow regimes show larger fluctuations in C_L than C_D.     

Flow-induced force fluctuations on a sphere at high Strouhal number

An empirical model is developed to estimate the broadband unsteady force spectrum induced on a rigid sphere in a nominally steady, uniform flow. The Reynolds number is sub-critical, and the frequency range considered is above the low-mode Strouhal shedding frequency of the sphere (0.5⩽fd/U₀⩽100, where f is the frequency, d is the diameter, and U₀ is the mean flow speed). The model uses the separation of variables assumption for the cross-power spectral densities of the surface pressure fluctuations. The assumption is shown to be a proper engineering approximation except in the lower part of the considered frequency range. In addition, the flow-induced unsteady lift and drag forces are measured independently of each other using towed spheres in a basin of water. Both estimations, from the empirical model and the data measured in the tow tank, show that the dimensionless power spectral densities of broadband unsteady lift and drag forces are constant for fd/U₀<1, and (fd/U₀)⁻³ dependent for 1⩽fd/U₀⩽100. The model predicts that the broadband spectral density of the unsteady lift force is about 5 dB higher than that of the unsteady drag force, while the measured data show the level difference between 3 and 7 dB. The empirical model presented here has application in predicting the flow-induced noise of underwater hydrophones that sense acoustic particle velocity or acceleration.     

CFD based mini- vs. micro-system delineation in elongated bubble flow regime

Delineation of mini- and micro-scale channels with respect to two-phase flow has been the subject of many research papers. There is no consensus on when the small channel can be characterized as a mini-channel or micro-channel. The idea proposed by this paper is to use the normalized bubble nose radius, liquid film thickness top over bottom ratio, and bubble shape contour, which are found under normal gravity conditions in slug flow through a horizontal adiabatic channel, as the delineation criteria. The input parameters are bubble nose radius and bubble nose velocity as the characteristic length scale and characteristic velocity scale respectively. 3D numerical simulation with ANSYS FLUENT was used to obtain the necessary data. Following CFD practice, a mesh independence study and a numerical model validation against published experimental data were both conducted. Analysis of the numerical simulation results showed that channels with D ⩽ 100 μm can be characterized as a micro-system, while channels with D ⩾ 400 μm belong to mini-systems. The region 200 μm ⩽ D ⩽ 300 μm represents a transition from the micro-scale to mini-scale.     

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.     

Flow structure around perforated cylinders in shallow water

The experimental investigations were carried out in order to have detailed information on the flow structure around perforated cylinders using high-image density Particle Image Velocimetry technique in shallow water flow. The depth-averaged free-stream velocity was kept constant as U_∞=100 mm/s corresponding to the Reynolds number of Re=10 000 based on the perforated cylinder diameter. In order to analyze the effect of porosity, β on the flow structure, the porosities in the range of 0.1≤β≤0.8 with an increment of 0.1 were used and the results were compared with the bare cylinder case by means of velocity and vorticity contours, turbulent kinetic energy, Reynolds shear stress and streamline topologies. It was concluded that the porosity, β had a substantial effect on the control of large-scale vortical structures downstream of the cylinder in which the shear layers were elongated, fluctuations were significantly attenuated and formation of Karman Vortex Street was successfully prevented by the use of perforated cylinders.     

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.     

Flow structures and heat transfer on dimples in a staggered arrangement

Vortex structures and heat transfer enhancement mechanism of turbulent flow over a staggered array of dimples in a narrow channel have been investigated using Large Eddy Simulation (LES), Laser Doppler Velocimetry (LDV) and pressure measurements for Reynolds numbers Re_H = 6521 and Re_H = 13,042.The flow and temperature fields are calculated by LES using dynamic mixed model applied both for the velocity and temperature. Simulations have been validated with experimental data obtained for smooth and dimpled channels and empiric correlations. The flow structures determined by LES inside the dimple are chaotic and consist of small eddies with a broad range of scales where coherent structures are hardly to detect. Proper Orthogonal Decomposition (POD) method is applied on resolved LES fields of pressure and velocity to identify spatial–temporal structures hidden in the random fluctuations. For both Reynolds numbers it was found that the dimple package with a depth h to diameter D ratio of h/D = 0.26 provides the maximum thermo-hydraulic performance. The heat transfer rate could be enhanced up to 201% compared to a smooth channel.     

The effects of freestream turbulence on the drag coefficient of a sphere

The effects of freestream turbulence intensity and integral length scale as freestream turbulent parameters on the drag coefficient of a sphere were experimentally investigated in a closed circuit wind tunnel. The Reynolds number, Re = Ud/ν, was varied from 2.2 × 10⁴ to 8 × 10⁴ by using spheres with diameter d of 20, 51 and 102 mm in addition to altering the freestream velocity, U. The freestream turbulence intensity Tu and flow integral length scale Λ were manipulated by the utilization of orifice perforated plates. The proper combination of orifice perforated plate hole diameter, sphere size, and sphere location along the center line of the wind tunnel enabled the independent alterations of turbulence intensity and relative integral length scale (Λ/d) from 1.8% to 10.7% and from 0.1 to 2.6, respectively, at each studied Reynolds number. Results show that over the range of conditions studied, the drag always decreases with increasing Tu and, the critical Reynolds number at which the drag coefficient is dramatically reduced is decreased by increasing Tu. Most interestingly, the drag at any particular Re and Tu may be significantly lowered by reducing Λ/d; this is particularly the case at high Re and Tu.     

Influence of aspect ratio on the mean flow field of a surface-mounted finite-height square prism

The flow around surface-mounted, finite-height square prisms at a Reynolds number of Re = 4.2 × 10⁴ was investigated experimentally in a low-speed wind tunnel using particle image velocimetry. The thickness of the boundary layer on the ground plane relative to the width of the prism was δ/D = 1.5. Four prism aspect ratios were tested, AR = 9, 7, 5, and 3, to study how the aspect ratio influences the flow field close to the prism. Upstream of the prism, lowering the aspect ratio from AR = 9 to AR = 3 causes the stagnation point on the upstream face to move closer to the free end, but there is no influence on the location and strength of the horseshoe vortex. Lowering the aspect ratio from AR = 9 to AR = 3 causes the cross-stream vortices in the upper and lower halves of the wake to move downstream and upstream, respectively; the latter vortex is absent for AR = 3, suggesting this prism sits below the critical aspect ratio. Above the free end of the prism, within the region of separated flow, lowering the aspect ratio from AR = 9 to AR = 3 shifts the location of the cross-stream vortex farther downstream. For the prism of AR = 3, reverse flow above the free end is stronger yet more unsteady compared to the more slender prisms, while the streamwise edge vortices are smaller and weaker.     

Theoretical stress analysis of intersecting cylindrical shells subjected to external loads transmitted through branch pipes

A thin shell theoretical solution of two normally intersecting cylindrical shells subjected to thrust-out force and three kinds of moments transmitted through branch pipes is presented in this paper. The solutions of modified Morley equation, which can be applicable up to ρ₀ = d/D ⩽ 0.8 and λ = d/(DT)^1/2 ⩽ 8 and the order of accuracy is raised to O(T/D), for the four loading cases are given. The accurate continuity conditions of generalized forces and displacements at the intersecting curve of two cylindrical shells for the four loading cases and the condition of the uniqueness of displacements are derived in this paper. The presented results are verified by experimental and numerical results successfully. They are in agreement with WRC Bulletin 297 when d/D is small.     

URANS of flow and endwall heat transfer in a pinned passage relevant to gas-turbine blade cooling

This paper presents some results of URANS study of flow and heat transfer in a matrix of wall-bounded 8 × 8 round pins, mimicking internal cooling passage of gas-turbine blades. The focus is on flow unsteadiness, its role in heat transfer and the capabilities of RANS models to reproduce these features in a set-up of industrial relevance. The results for two Reynolds numbers, 10 000 and 30 000, are compared with the available experiments and LES. It is shown that the elliptic-relaxation eddy-viscosity model, ζ-f captures vortex shedding and the consequent gross effects on the flow development. However, a closer look at flow details reveals discrepancies, especially around the first three pin rows, where the unsteadiness reproduced by URANS shows much weaker amplitudes as compared with LES. Only further downstream the succession of forcing from a series of pins produced unsteadiness akin to those captured by LES. The comparison suggests that smaller structures undetected by URANS need to be resolved to capture properly the separation and wake characteristics of each row. At Re = 10 000, the average endwall Nusselt number agrees well with the LES, both being about 20% lower than in the experiment. For Re = 30 000 the URANS Nusselt is within 10% of the experimental value.