An experimental investigation of flow past a dual step cylinder

Flow development in the wake of a dual step cylinder has been investigated experimentally using Laser Doppler Velocimetry and flow visualization. The dual step cylinder model is comprised of a large diameter cylinder (D) mounted at the mid-span of a small diameter cylinder (d). The experiments have been performed for a Reynolds number (Re _D ) of 1,050, a diameter ratio (D/d) of 2, and a range of large cylinder aspect ratios (L/D). The results show that the flow development is highly dependent on L/D. The following four distinct flow regimes can be identified based on vortex dynamics in the wake of the large cylinder: (1) for L/D ≥ 15, three vortex shedding cells form in the wake of the large cylinder, one central cell bounded by two cells of lower frequency, (2) for 8 < L/D ≤ 14, a single vortex shedding cell forms in the wake of the large cylinder, (3) for 2 < L/D ≤ 6, vortex shedding from the large cylinder is highly three-dimensional. When spanwise vortices are shed, they deform substantially and attain a hairpin shape in the near wake, (4) for 0.2 ≤ L/D ≤ 1, the large cylinder induces vortex dislocations between small cylinder vortices. The results show that for Regimes I to III, on the average, the frequency of vortex shedding in the large cylinder wake decreases with L/D, which is accompanied by a decrease in coherence of the shed vortices. In Regime IV, small cylinder vortices connect across the large cylinder wake, but these connections are interrupted by vortex dislocations. With decreasing L/D, the frequency of dislocations decreases and the dominant frequency in the large cylinder wake increases toward the small cylinder shedding frequency.     

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.     

Application of digital particle image velocimetry to insect aerodynamics: measurement of the leading-edge vortex and near wake of a Hawkmoth

Some insects use leading-edge vortices to generate high lift forces, as has been inferred from qualitative smoke visualisations of the flow around their wings. Here we present the first Digital Particle Image Velocimetry (DPIV) data and quantitative analysis of an insect’s leading-edge vortex and near wake at two flight speeds. This allows us to describe objectively 2D slices through the flow field of a tethered Tobacco Hawkmoth (Manduca sexta). The near-field vortex wake appears to braodly resemble elliptical vortex loops. The presence of a leading-edge vortex towards the end of the downstroke is found to coincide with peak upward force production measured by a six-component force–moment balance. The topology of Manduca’s leading-edge vortex differs from that previously described because late in the downstroke, the structure extends continuously from wingtip across the thorax to the other wingtip.     

Effect of Schmidt number on the structure and propagation of density currents

The results of a numerical study of two- and three-dimensional Boussinesq density currents are described. They are aimed at exploring the role of the Schmidt number on the structure and dynamics of density driven currents. Two complementary approaches are used, namely a spectral method and a finite-volume interface capturing method. They allow for the first time to describe density currents in the whole range of Schmidt number 1 ≤ Sc ≤ ∞ and Reynolds number 10² ≤ Re ≤ 10⁴. The present results confirm that the Schmidt number only weakly influences the structure and dynamics of density currents provided the Reynolds number of the flow is large, say of O(10⁴) or more. On the contrary low- to moderate-Re density currents are dependant on Sc as the structure of the mixing region and the front velocities are modified by diffusion effects. The scaling of the characteristic density thickness of the interface has been confirmed to behave as (ScRe)^−1/2. Three-dimensional simulations suggest that the patterns of lobes and clefts are independent of Sc. In contrast the Schmidt number is found to affect dramatically (1) the shape of the current head as a depression is observed at high-Sc, (2) the formation of vortex structures generated by Kelvin–Helmholtz instabilities. A criterion is proposed for the stability of the interface along the body of the current based on the estimate of a bulk Richardson number. This criterion, derived for currents of arbitrary density ratio, is in agreement with present computed results as well as available experimental and numerical data.      

Dual leading-edge vortex structure for flow over a simplified butterfly model

The dye visualization experiments show that a dual leading-edge vortex (LEV) structure exists on the suction side of a simplified butterfly model of Papilio ulysses at α = 8°−12°. Furthermore, the results of particle image velocimetry (PIV) measurement indicate that the axial velocity of the primary (outer) vortex core reaches the lower extreme value while a transition from a “wake-like” to a “jet-like” axial velocity profile occurs. The work reveals for the first time the existence of dual LEV structure on the butterfly-like forward-sweep wing configuration.     

The formation of vortex rings in a strongly forced round jet

The periodic formation of vortex rings in the developing region of a round jet subjected to high-amplitude acoustic forcing is investigated with High-Speed Particle Image Velocimetry. Harmonic velocity oscillations ranging from 20 to 120% of the mean exit velocity of the jet was achieved at several forcing frequencies determined by the acoustic response of the system. The time-resolved history of the formation process and circulation of the vortex rings are evaluated as a function of the forcing conditions. Overall, high-amplitude forcing causes the shear layers of the jet to breakup into a train of large-scale vortex rings, which share many of the features of starting jets. Features of the jet breakup such as the roll-up location and vortex size were found to be both amplitude and frequency dependent. A limiting time-scale of t/T ≈ 0.33 based on the normalized forcing period was found to restrict the growth of a vortex ring in terms of its circulation for any given arrangement of jet forcing conditions. In sinusoidally forced jets, this time-scale corresponds to a kinematic constraint where the translational velocity of the vortex ring exceeds the shear layer velocity that imposes pinch-off. This kinematic constraint results from the change in sign in the jet acceleration between t = 0 and t = 0.33T. However, some vortex rings were observed to pinch-off before t = 0.33T suggesting that they had acquired their maximum circulation. By invoking the slug model approximations and defining the slug parameters based on the experimentally obtained time- and length-scales, an analytical model based on the slug and ring energies revealed that the formation number for a sinusoidally forced jet is L/D ≈ 4 in agreement with the results of Gharib et al. (J Fluid Mech 360:121–140, 1998).     

Flow past two in-tandem airfoils undergoing sinusoidal oscillations

The interaction of the wake, generated behind an upstream oscillating NACA 0012 airfoil, with the downstream NACA 0012 airfoil, oscillated at the same conditions but with ϕ = 0° and 180° different phases (relative to the upstream airfoil), was investigated by particle image velocimetry and surface pressure measurements. The results show that the axial spacing and phase difference determined the strength of the undesirable interference effects and, subsequently, the behavior of the dynamic-load loops of the downstream airfoil. The boundary-layer events on the downstream airfoil were persistently different from those observed on the baseline oscillating airfoil. The downwash induced by the upstream airfoil disrupted leading-edge vortex (LEV) formation on the downstream airfoil. The absence of LEV-induced transient effects also led to a significantly lowered aerodynamic loading and C _l -hysteresis compared to the baseline airfoil. The aerodynamic performance of the ϕ = 180° case, however, outperformed that of the ϕ = 0° case.     

Numerical study of self-sustained oscillatory flows and heat transfer in channels with a tandem of transverse vortex generators

A numerical investigation was conducted into channel flows with a tandem of transverse vortex generators in the form of rectangular cylinders. The oscillatory behavior of the flow is studied. Data for heat transfer and flow losses are presented for 100≤Re≤400 and cylinder separation distances 1≤S/H≤4. The results are obtained by numerical solution of the full Navier-Stokes equations and the energy equation. Self-sustained flow oscillations are found for Re>100. Alternate and dynamic shedding of large vortex structures from the cylinders is observed by visualization of the numerically determined flow field. A heat transfer enhancement up to a factor 1.78 compared to plane channel flow is observed. Received on 16 July 1997     

Non-Darcy buoyancy driven flows in a fluid saturated porous medium: the use of asymptotic computational fluid dynamics (ACFD) approach

Natural convection in a fluid saturated porous medium has been numerically investigated using a generalized non-Darcy approach. The governing equations are solved by using Finite Volume approach. First order upwind scheme is employed for convective formulation and SIMPLE algorithm for pressure velocity coupling. Numerical results are presented to study the influence of parameters such as Rayleigh number (10⁶ ≤Ra ≤10⁸), Darcy number (10⁻⁵ ≤ Da ≤ 10⁻²), porosity (0.4 ≤ ɛ ≤ 0.9) and Prandtl number (0.01 ≤ Pr ≤ 10) on the flow behavior and heat transfer. By combining the method of matched asymptotic expansions with computational fluid dynamics (CFD), so called asymptotic computational fluid dynamics (ACFD) technique has been employed to generate correlation for average Nusselt number. The technique is found to be an attractive option for generating correlation and also in the analysis of natural convection in porous medium over a fairly wide range of parameters with fewer simulations for numerical solutions.     

Non-Darcian and Anisotropic Effects on Free Convection in a Porous Enclosure

The free convective flow and heat transfer, within the framework of Boussinesq approximation, in an anisotropic fluid filled porous rectangular enclosure subjected to end-to-end temperature difference have been investigated using Brinkman extended non-Darcy flow model. The studies involve simultaneous consideration of hydrodynamic and thermal anisotropy. The flow and temperature fields in general are governed by, Ra, the Rayleigh number, AR, the aspect ratio of the slab, K*, the permeability ratio and k*, the thermal conductivity ratio, and Da, Darcy number. Numerical solutions employing the successive accelerated replacement (SAR) scheme have been obtained for 100 ≤ Ra ≤ 1000, 0.5 ≤ AR ≤ 5, 0.5 ≤ K* ≤ 5, 0.5 ≤ k* ≤ 5, and 0 ≤ Da ≤ 0.1. It has been found that [`(Nu)]{\overline {Nu}}, average Nusselt number increases with increase in K* and decreases as k* increases. However, the magnitude of the change in [`(Nu)]{\overline {Nu}} depends on the parameter Da, characterizing the Brinkman extended non-Darcy flow.     

Numerical study of transient laminar natural convection heat transfer over a sphere subjected to a constant heat flux

Transient laminar natural convection over a sphere which is subjected to a constant heat flux has been studied numerically for high Grashof numbers (10⁵ ≤ Gr ≤ 10⁹) and a wide range of Prandtl numbers (Pr = 0.02, 0.7, 7, and 100). A plume with a mushroom-shaped cap forms above the sphere and drifts upward continuously with time. The size and the level of temperature of the transient cap and plume stem decrease with increasing Gr and Pr. Flow separation and an associated vortex may appear in the wake of the sphere depending on the magnitude of Gr and Pr. A recirculation vortex which appears and grows until “steady state” is attained was found only for the very high Grashof numbers (10⁵ ≤ Gr ≤ 10⁹) and the lowest Prandtl number considered (Pr = 0.02). The appearance and subsequent disappearance of a vortex was observed for Gr = 10⁹ and Pr = 0.7. Over the lower hemisphere, the thickness of both the hydrodynamic (δ_H) and the thermal (δ_T) boundary layers remain nearly constant and the sphere surface is nearly isothermal. The surface temperature presents a local maximum in the wake of the sphere whenever a vortex is established in the wake of the sphere. The surface pressure recovery in the wake of the sphere increases with decreasing Pr and with increasing Gr. For very small Pr, unlike forced convection, the ratio δ_T/δ_H remains close to unity. The results are in good agreement with experimental data and in excellent agreement with numerical results available in the literature. A correlation has also been presented for the overall Nusselt number as a function of Gr and Pr.     

Numerical investigation of natural convection of Al<Subscript>2</Subscript>O<Subscript>3</Subscript> nanofluid in vertical annuli

This paper presents the numerical study of internal free convection of Al₂O₃ water nanofluid in vertical annuli. Vertical walls are maintained at constant temperatures and horizontal walls are adiabatic. Results are validated by experimental data. Effect of nanofluids on natural convection is investigated as a function of geometrical and physical parameters and particle fractions for aspect ratio of 1 ≤ H/L ≤ 5, Grashof number of 10³ ≤ Gr ≤ 10⁵ and concentration of 0 ≤ ϕ ≤ 0.06. More than 330 different numerical cases are investigated to develop a new correlation for the Nusselt number. This correlation is presented as a function of Nusselt number of base fluid and particle fraction which is a linear decreasing function of particle fraction. The developed correlation for annuli is also valid for the natural convection of Al₂O₃ water nanofluid in a square cavity. Furthermore, the effect of the viscosity and conductivity models on the Nusselt number of nanofluids in cylindrical cavities are discussed.     

Investigation of the different base fluid effects on the nanofluids heat transfer and pressure drop

A numerical study of laminar forced convective flows of three different nanofluids through a horizontal circular tube with a constant heat flux condition has been performed. The effect of Al₂O₃ volume concentration 0 ≤ φ ≤ 0.09 in the pure water, water-ethylene glycol mixture and pure ethylene glycol as base fluids, and Reynolds number of 100 ≤ Re ≤ 2,000 for different power inputs in the range of 10 ≤ Q(W) ≤ 400 have been investigated. In this study, all of the nanofluid properties are temperature and nanoparticle volume concentration dependent. The governing equations have been solved using finite volume approach with the SIMPLER algorithm. The results indicate an increase in the averaged heat transfer coefficient with increasing the mass of ethylene glycol in the water base fluid, solid concentration and Reynolds number. From the investigations it can be inferred that, the pressure drop and pumping power in the nanofluids at low solid volumetric concentration (φ < 3%) is approximately the same as in the pure base fluid in the various Reynolds numbers, but the higher solid nanoparticle volume concentration causes a penalty drop in the pressure. Moreover, this study shows it is possible to achieve a higher heat transfer rate with lower wall shear stress with the use of proper nanofluids.     

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     

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.     

Characterization of fluid flow patterns and heat transfer in horizontal channel mixed convection

Two mechanisms of roll initiation are highlighted in a horizontal channel flow, uniformly heated from below, at constant heat flux (Γ = 10, Pr = 7, 50 ≤ Re ≤ 100, 0 ≤ Ra ≤ 10⁶). The first mechanism is the classical one, it occurs for low Rayleigh numbers and is initiated by the lateral wall effect. The second occurs for higher Rayleigh numbers and combines the previous effect with a supercritical vertical temperature gradient in the lower boundary layer, which simultaneously triggers pairs of rolls in the whole zone in between the two lateral rolls. We have found that in the present configuration, the transition between the two roll initiation mechanisms occurs for Ra/Re ² ≈ 18. Consequently, the heat transfer is significantly enhanced compared to the pure forced convection case owing to the flow pattern responsible of the continuous flooding the heated wall with cold fluid.     

Buoyancy and cross flow effects on heat transfer of multiple impinging slot air jets cooling a flat plate at different orientations

The present article reports on heat transfer characteristics associated with multiple laminar impinging air jet cooling a hot flat plat at different orientations. The work aims to study the interactions of the effects of cross flow, buoyancy induced flow, orientation of the hot surface with respect to gravity, Reynolds numbers and Rayleigh numbers on heat transfer characteristics. Experiments have been carried out for different values of jet Reynolds number, Rayleigh number and cross flow strength and at different orientations of the air jet with respect to the target hot plate. In general, the effective cooling of the plate has been observed to be increased with increasing Reynolds number and Rayleigh number. The results concluded that the hot surface orientation is important for optimum performance in practical applications. It was found that for Re ≥ 400 and Ra ≥ 10,000 (these ranges give 0.0142 ≤ Ri ≤ 1.59 the Nusselt number is independent on the hot surface orientation. However, for Re ≤ 300 and Ra ≥ 100,000 (these ranges give 1.59 ≤ Ri ≤ 42.85): (i) the Nusselt number for horizontal orientation with hot surface facing down is less that that of vertical orientation and that of horizontal orientation with hot surface facing up, and (ii) the Nusselt number of vertical orientation is approximately the same as that of horizontal orientation with hot surface facing up. For all surfaces orientations and for the entire ranges of Re and Ra, it was found that increasing the cross flow strength decreases the effective cooling of the surface.     

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

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

Convection in a Lid-Driven Heat-Generating Porous Cavity with Alternative Thermal Boundary Conditions

Mixed convection flow in a two-sided lid-driven cavity filled with heat-generating porous medium is numerically investigated. The top and bottom walls are moving in opposite directions at different temperatures, while the side vertical walls are considered adiabatic. The governing equations are solved using the finite-volume method with the SIMPLE algorithm. The numerical procedure adopted in this study yields a consistent performance over a wide range of parameters that were 10⁻⁴ ≤ Da ≤ 10⁻¹ and 0 ≤ Ra _I ≤ 10⁴. The effects of the parameters involved on the heat transfer characteristics are studied in detail. It is found that the variation of the average Nusselt number is non-linear for increasing values of the Darcy number with uniform or non-uniform heating condition.     

Optimal design of angled rib turbulators in a cooling channel

In the present study, an optimal design for enhancement of heat transfer and thermal performance for a stationary channel with angled rib turbulators was investigated to find the most suitable rib geometry. Among various design parameters, two design variables, rib angle of attack (α) and pitch-to-rib height (p/e), were chosen. The ranges of two design variables were set as 30° ≤ α ≤ 80° and 3.0 ≤ p/e ≤ 15.0. Approximations for design of the best rib turbulators were obtained using the advanced response surface method with functional variables. The second-order response surfaces (or correlations) within the ranges of two design variables were completed by this method. As for the optimized results, maximum averaged heat transfer value was obtained at α = 53.31° and p/e = 6.50, while the highest thermal performance value was presented at α = 54.67° and p/e = 6.80.