Simulations of Turbulent Flow in a Plane Asymmetric Diffuser

Large-eddy simulations (LES) of a planar, asymmetric diffuser flow have been performed. The diverging angle of the inclined wall of the diffuser is chosen as 8.5°, a case for which recent experimental data are available. Reasonable agreement between the LES and the experiments is obtained. The numerical method is further validated for diffuser flow with the diffuser wall inclined at a diverging angle of 10°, which has served as a test case for a number of experimental as well as numerical studies in the literature (LES, RANS). For the present results, the subgrid-scale stresses have been closed using the dynamic Smagorinsky model. A resolution study has been performed, highlighting the disparity of the relevant temporal and spatial scales and thus the sensitivity of the simulation results to the specific numerical grids used. The effect of different Reynolds numbers of the inflowing, fully turbulent channel flow has been studied, in particular, Re _b  = 4,500, Re _b  = 9,000 and Re _b  = 20,000 with Re _b being the Reynolds number based on the bulk velocity and channel half width. The results consistently show that by increasing the Reynolds number a clear trend towards a larger separated region is evident; at least for the studied, comparably low Reynolds-number regime. It is further shown that the small separated region occurring at the diffuser throat shows the opposite behaviour as the main separation region, i.e. the flow is separating less with higher Re _b . Moreover, the influence of the Reynolds number on the internal layer occurring at the non-inclined wall described in a recent study has also been assessed. It can be concluded that this region close to the upper, straight wall, is more distinct for larger Re _b . Additionally, the influence of temporal correlations arising from the commonly used periodic turbulent channel flow as inflow condition (similar to a precursor simulation) for the diffuser is assessed.     

Scanning PIV investigation of the laminar separation bubble on a SD7003 airfoil

A laminar separation bubble occurs on the suction side of the SD7003 airfoil at an angle of attack α =  4–8° and a low Reynolds number less than 100,000, which brings about a significant adverse aerodynamic effect. The spatial and temporal structure of the laminar separation bubble was studied using the scanning PIV method at α =  4° and Re = 60,000 and 20,000. Of particular interest are the dynamic vortex behavior in transition process and the subsequent vortex evolution in the turbulent boundary layer. The flow was continuously sampled in a stack of parallel illuminated planes from two orthogonal views with a frequency of hundreds Hz, and PIV cross-correlation was performed to obtain the 2D velocity field in each plane. Results of both the single-sliced and the volumetric presentations of the laminar separation bubble reveal vortex shedding in transition near the reattachment region at Re = 60,000. In a relatively long distance vortices characterized by paired wall-normal vorticity packets retain their identities in the reattached turbulent boundary layer, though vortices interact through tearing, stretching and tilting. Compared with the restricted LSB at Re = 60,000, the flow at Re = 20,000 presents an earlier separation and a significantly increased reversed flow region followed by “huge” vortical structures.     

Stabilization and evolution of parameters of a symmetric laminar flow in a plane channel with sudden expansion

A problem is formulated for computing the fields of parameters of a stationary laminar symmetric flow. A two-dimensional flow in a channel with a sudden change in the cross-sectional area is computed. The evolution of a three-dimensional perturbation inserted into the channel at the initial stage of computations is analyzed. It is demonstrated that the parameters of a two-dimensional flow in the channel at a Reynolds number Re = 50 become stabilized at a dimensionless time t > 20, whereas the steady state is reached under the same conditions at t ≈ 100. At a distance of approximately 10h (h is the channel width at the entrance), the flow becomes one-dimensional, but the streamwise component of the velocity vector remains a function of the streamwise coordinate owing to flow compressibility. __________ Translated from PrikladnayaMekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 35–42, January–February, 2007.     

Experimental investigation of mixed convection in a rectangular duct with a heated plate in the middle of cross section

The flow and heat transfer in an inclined and horizontal rectangular duct with a heated plate longitudinally mounted in the middle of cross section was experimentally investigated. The heated plate and rectangular duct were both made of highly conductive materials, and the heated plate was subjected to a uniform heat flux. The heat transfer processes through the test section were under various operating conditions: Pr ≈ 0.7, inclination angle ϕ = −60° to +60°, Reynolds number Re = 334–1,911, Grashof number Gr = 5.26 × 10²–5.78 × 10⁶. The experimental results showed that the average Nusselt number in the entrance region was 1.6–2 times as large as that in the fully developed region. The average Nusselt numbers and pressure drops increased with the Reynolds number. The average Nusselt numbers and pressure drops decreased with an increase in the inclination angle from −60° to +60° when the Reynolds number was less than 1,500. But when the Reynolds number increased to over about 1,800, the heat transfer coefficients and pressure drops were independent of inclination angles.     

Correlation of the separation region length in shock wave/channel boundary layer interaction

 Results of an experimental investigation of the characteristics of a separation region induced by the interaction of an externally generated oblique shock with the turbulent boundary layer formed in a rectangular half channel are discussed. The experiments were carried out in the supersonic wind tunnel of the Institute of Theoretical and Applied Mechanics SB RAS at a free-stream Mach number M _∞=3.01 over a range of Reynolds numbers Re ₁=(9.7–47.5)×10⁶ m^-1 and at zero incidence and zero yaw of the model. Particular attention is paid to the size of the zone of the upstream propagation of disturbances (upstream influence region) under different experimental conditions: with varied values of the shock wave strength, half channel width, and Reynolds number. It is shown, in particular, that the normalized upstream influence region length as a function of inclination angle of the shock generator in a rectangular half channel is readily approximated by a simple exponential function. In support of the known reference data obtained for supersonic numbers M _∞ and moderate Re in other configurations, it is also shown that the upstream influence region length decreases with increasing Reynolds number. Generalization of experimental data on the length of the upstream influence region formed in similar geometric configurations is possible using an additional reference linear scale which is the distance from the leading edge of the shock generator to the exposed surface. A substantial dependence of the reference dimensions of separation region on the half channel width is also established. Received: 20 January 1997/Accepted: 30 June 1997     

Time resolved analysis of steady and oscillating flow in the upper human airways

In this experimental study a thorough analysis of the steady and unsteady flow field in a realistic transparent silicone lung model of the first bifurcation of the upper human airways will be presented. To determine the temporal evolution of the flow time-resolved particle-image velocimetry recordings were performed for a Womersley number range 3.3 ≤ α ≤ 5.8 and Reynolds numbers of Re _D = 1,050, 1,400, and 2,100. The results evidence a highly three-dimensional and asymmetric character of the velocity field in the upper human airways, in which the influence of the asymmetric geometry of the realistic lung model plays a significant role for the development of the flow field in the respiratory system. At steady inspiration, the flow shows independent of the Reynolds number a large zone with embedded counter-rotating vortices in the left bronchia ensuring a continuous streamwise transport into the lung. At unsteady flow the critical Reynolds number, which describes the onset of vortices in the first bifurcation, is increased at higher Womersley number and decreased at higher Reynolds number. At expiration the unsteady and steady flows are almost alike.     

Investigation of the vortex induced unsteadiness of a separation bubble via time-resolved and scanning PIV measurements

A transitional separation bubble on the suction side of an SD7003 airfoil is considered. The transition process that forces the separated shear layer to reattach seems to be governed by Kelvin–Helmholtz instabilities. Large scale vortices are formed due to this mechanism at the downstream end of the bubble. These vortices possess a three-dimensional structure and detach from the recirculation region, while other vortices are formed within the bubble. This separation of the vortex is a highly unsteady process, which leads to a bubble flapping. The structure of these vortices and the flapping of the separation bubble due to these vortices are temporally and spatially analyzed at angles of attack from 4° to 8° and chord-length based Reynolds numbers Re _c = 20,000–60,000 using time-resolved PIV measurements in a 2D and a 3D set-up, i.e., stereo-scanning PIV measurements are done in the latter case. These measurements complete former studies at a Reynolds number of Re _c = 20,000. The results of the time-resolved PIV measurements in a single light-sheet show the influence of the angle of attack and the Reynolds number. The characteristic parameters of the separation bubble are analyzed focusing on the unsteadiness of the separation bubble, e.g., the varying size of the main recirculation region, which characterizes the bubble flapping, and the corresponding Strouhal number are investigated. Furthermore, the impact of the freestream turbulence is investigated by juxtaposing the current and former results. The stereo-scanning PIV measurements at Reynolds numbers up to 60,000 elucidate the three-dimensional character of the vortical structures, which evolve at the downstream end of the separation bubble. It is shown that the same typical structures are formed, e.g., the c-shape vortex and the screwdriver vortex at each Reynolds number and angle of attack investigated and the occurrence of these patterns in relation to Λ-structures is discussed. To evidence the impact of the freestream turbulence, these results are compared with findings of former measurements.     

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

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

Effect of magnetic Reynolds number on the two‐dimensional hydromagnetic flow around a cylinder

Numerical experiments have been conducted to study the effect of magnetic Reynolds number on the steady, two‐dimensional, viscous, incompressible and electrically conducting flow around a circular cylinder. Besides usual Reynolds number Re, the flow is governed by the magnetic Reynolds number R_m and Alfvén number β. The flow and magnetic field are uniform and parallel at large distances from the cylinder. The pressure Poisson equation is solved to find the pressure fields in the entire flow region. The effects of the magnetic field and electrical conductivity on the recirculation bubble, drag coefficient, standing vortex and pressure are presented and discussed. For low interaction parameter (N<1), the suppression of the flow‐separation is nearly independent of the conductivity of the fluid, whereas for large interaction parameters, the conductivity of the fluid strongly influences the control of flow‐separation. Copyright © 2008 John Wiley & Sons, Ltd.     

Flow Bifurcations in a Thin Gap Between Two Rotating Spheres

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

Forced convection flow over a microsphere

Forced convection heat transfer characteristics around a microsphere subjected to uniform heat flux boundary condition is numerically investigated in this study. Moderate to high values of Reynolds number and a wide range of Prandtl number are considered. The analysis assumes that the continuity assumption is valid and hence the Navier–Stokes equations are solved for the range of Knudsen number of 0.001 ≤ Kn ≤ 0.1. The appropriate boundary conditions at the surface of the microsphere; the velocity slip and temperature jump are applied. The effect of the flow parameters: Re, Pr and Kn on the velocity and temperature distribution is presented and hence a better control on the boundary layer thickness can be achieved in the microscale level. Furthermore, the effect of the controlling parameters on the delay of flow separation, reduced shear stress, drag coefficient and on the Nusselt number profiles is also presented in the results.     

Flow past a square cylinder at low Reynolds numbers

Results are presented for the flow past a stationary square cylinder at zero incidence for Reynolds number, Re ? 150. A stabilized finite‐element formulation is employed to discretize the equations of incompressible fluid flow in two‐dimensions. For the first time, values of the laminar separation Reynolds number, Re_s, and separation angle, θ_s, at Re_s are predicted. Also, the variation of θ_s with Re is presented. It is found that the steady separation initiates at Re = 1.15. Contrary to the popular belief that separation originates at the rear sharp corners, it is found to originate from the base point, i.e. θ_s=180^° at Re = Re_s. For Re > 5, θ_s approaches the limit of 135 ^°. The length of the separation bubble increases approximately linearly with increasing Re. The drag coefficient varies as Re^?0.66. Flow characteristics at Re ? 40 are also presented for elliptical cylinders of aspect ratios 0.2, 0.5, 0.8 and 1 (circle) having the same characteristic dimension as the square and major axis oriented normal to the free‐stream. Compared with a circular cylinder, the flow separates at a much lower Re from a square cylinder leading to the formation of a bigger wake (larger bubble length and width). Consequently, at a given Re, the drag on a square cylinder is more than the drag of a circular cylinder. This suggests that a cylinder with square section is more bluff than the one with circular section. Among all the cylinder shapes studied, the square cylinder with sharp corners generates the largest amount of drag. Copyright © 2010 John Wiley & Sons, Ltd.     

An experimental investigation of moderate reynolds number flow in a T-Channel

An experimental investigation of water flow in a T-shaped channel with rectangular cross section (20 × 20 mm inlet ID and 20 × 40 mm outlet ID) has been conducted for a Reynolds number Re range of 56–422, based on inlet diameter. Dynamical conditions and the T-channel geometry of the current study are applicable to the microscale. 2-D planar particle imaging velocimetry (PIV) and laser-induced fluorescence (LIF) were used in multiple locations of the T-channel to investigate local dynamical behaviors. Steady symmetric and asymmetric flow regimes predicted in the literature, which is largely numerical, are experimentally verified. Unsteady flow regimes, which are numerically predicted to occur at higher Re but have not yet been experimentally characterized, are also examined, and real-time LIF results illuminate the evolution of unsteady structure. Experimental data of the present resolution and scope are not presently available for unsteady flow regimes. Time scales are presented for unsteady flow regimes, which are found to exhibit periodic behavior and to occur for Re  ≥ 195. An unsteady symmetrical regime is identified for Re ≥ 350 that is detrimental to mixing. Momentum fields and dynamical behaviors of all flow regimes are characterized in detail, such that published mixing trends may be better understood. Results of all experimental trials were used to construct a regime map. A symmetric topology is found to be dominant for Re from 56 to 116, when flow is steady, and 350 to 422, when flow is characterized by unsteady stagnation-point oscillation in the T-channel junction. Asymmetric flow, which is positively indicated for mixing, is dominant for Re between 142 and 298, and the fluid interface exhibits both steady (two standing vortices) and unsteady (shear-layer type roll-up) behaviors. This result is based on multiple experiments and suggests a practical operating range of 142  ≤ Re ≤ 298 where asymmetric flow is highly likely to experimentally occur. The identification of an upper limit on Re,  beyond which mixing appears negatively impacted by a more symmetrical momentum field, is practically important as pressure drops on the microscale are significant.     

Computation of Vortex Shedding and Radiated Sound for a Circular Cylinder: Subcritical to Transcritical Reynolds Numbers

The Lighthill acoustic analogy combined with Reynolds-averaged Navier–Stokes flow computations are used to investigate the ability of existing technology to predict the tonal noise generated by vortex shedding from a circular cylinder for a range of Reynolds numbers (100 < Re < 5 million). Computed mean drag, mean coefficient of pressure, Strouhal number, and fluctuating lift are compared with experiment. Two-dimensional calculations produce a Reynolds number trend similar to experiment but incorrectly predict many of the flow quantities. Different turbulence models give inconsistent results in the critical Reynolds number range (Re≈ 100000). The computed flow field is used as input for noise prediction. Two-dimensional inputs overpredict both noise amplitude and frequency; however, if an appropriate correlation length is used, predicted noise amplitudes agree with experiment. Noise levels and frequency content agree much better with experiment when three-dimensional flow computations are used as input data. Received 5 May 1998 and accepted 28 September 1998     

The wake of falling disks at low Reynolds numbers

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

Parametric study of unsteady laminar flows

Results of a parametric study of unsteady laminar flows are analyzed. Three-dimensional unsteady equations of hydromechanics for a compressible medium are solved. The range of the characteristic Reynolds number Re = 400–900 is considered. It is demonstrated that the laminar flow in a plane channel ceases to be steady at Re = 415. As the Reynolds number increases, the unsteady processes become more intense, disturbances penetrate inward the channel, and separation zones lose their stability. In the vicinity of the channel exit, however, the flow tends to stabilize, though it remains unsteady. No transition to a turbulent flow occurs in the examined range of Reynolds numbers.     

Effect of bulk viscosity on Kelvin-Helmholtz instability

An energy functional leading to a resolvable variational problem for determining the critical Reynolds number of laminar-turbulent transition Re _cr is constructed within the framework of the nonlinear energy stability theory of compressible flows. Asymptotic estimates containing the characteristic dependence in the main order are obtained for the stability of various modes of Couette compressible gas flow. The asymptotics considered are long-wave approximations. This suggests that the obtained dependence describes the effect of bulk viscosity on the large-scale vortex structures characteristic of Kelvin-Helmholtz instability. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 73–84, May–June, 2008.     

High-fidelity simulations of moving and flexible airfoils at low Reynolds numbers

The present paper highlights results derived from the application of a high-fidelity simulation technique to the analysis of low-Reynolds-number transitional flows over moving and flexible canonical configurations motivated by small natural and man-made flyers. This effort addresses three separate fluid dynamic phenomena relevant to small fliers, including: laminar separation and transition over a stationary airfoil, transition effects on the dynamic stall vortex generated by a plunging airfoil, and the effect of flexibility on the flow structure above a membrane airfoil. The specific cases were also selected to permit comparison with available experimental measurements. First, the process of transition on a stationary SD7003 airfoil section over a range of Reynolds numbers and angles of attack is considered. Prior to stall, the flow exhibits a separated shear layer which rolls up into spanwise vortices. These vortices subsequently undergo spanwise instabilities, and ultimately breakdown into fine-scale turbulent structures as the boundary layer reattaches to the airfoil surface. In a time-averaged sense, the flow displays a closed laminar separation bubble which moves upstream and contracts in size with increasing angle of attack for a fixed Reynolds number. For a fixed angle of attack, as the Reynolds number decreases, the laminar separation bubble grows in vertical extent producing a significant increase in drag. For the lowest Reynolds number considered (Re _c  = 10⁴), transition does not occur over the airfoil at moderate angles of attack prior to stall. Next, the impact of a prescribed high-frequency small-amplitude plunging motion on the transitional flow over the SD7003 airfoil is investigated. The motion-induced high angle of attack results in unsteady separation in the leading edge and in the formation of dynamic-stall-like vortices which convect downstream close to the airfoil. At the lowest value of Reynolds number (Re _c  = 10⁴), transition effects are observed to be minor and the dynamic stall vortex system remains fairly coherent. For Re _c  = 4 × 10⁴, the dynamic-stall vortex system is laminar at is inception, however shortly afterwards, it experiences an abrupt breakdown associated with the onset of spanwise instability effects. The computed phased-averaged structures for both values of Reynolds number are found to be in good agreement with the experimental data. Finally, the effect of structural compliance on the unsteady flow past a membrane airfoil is investigated. The membrane deformation results in mean camber and large fluctuations which improve aerodynamic performance. Larger values of lift and a delay in stall are achieved relative to a rigid airfoil configuration. For Re _c = 4.85 × 10⁴, it is shown that correct prediction of the transitional process is critical to capturing the proper membrane structural response.     

Some Reynolds number effects on two- and three-dimensional turbulent boundary layers

Experimental data for a two-dimensional (2-D) turbulent boundary layer (TBL) flow and a three-dimensional (3-D) pressure-driven TBL flow outside of a wing/body junction were obtained for an approach Reynolds number based on momentum thickness of Re _θ =23,200. The wing shape had a 3:2 elliptical nose, NACA 0020 profiled tail, and was mounted on a flat wall. Some Reynolds number effects are examined using fine spatial resolution (Δy ⁺=1.8) three-velocity-component laser-Doppler velocimeter measurements of mean velocities and Reynolds stresses at nine stations for Re _θ =23,200 and previously reported data for a much thinner boundary layer at Re _θ =5,940 for the same wing shape. In the 3-D boundary layers, while the stress profiles vary considerably along the flow due to deceleration, acceleration, and skewing, profiles of the parameter correlate well and over available Reynolds numbers. The measured static pressure variations on the flat wall are similar for the two Reynolds numbers, so the vorticity flux and the measured mean velocities scaled on wall variables agree closely near the wall. The stresses vary similarly for both cases, but with higher values in the outer region of the higher Re _θ case. The outer layer turbulence in the thicker high Reynolds number case behaves similarly to a rapid distortion of the flow, since stream-wise vortical effects from the wall have not diffused completely through the boundary layer at all measurement stations. Received: 9 June 2000/Accepted: 26 January 2001     

LES and RANS for Turbulent Flow over Arrays of Wall-Mounted Obstacles

Large-eddy simulation (LES) has been applied to calculate the turbulent flow over staggered wall-mounted cubes and staggered random arrays of obstacles with area density 25%, at Reynolds numbers between 5 × 10³ and 5 10⁶, based on the free stream velocity and the obstacle height. Re = 5 × 10³ data were intensively validated against direct numerical simulation (DNS) results at the same Re and experimental data obtained in a boundary layer developing over an identical roughness and at a rather higher Re. The results collectively confirm that Reynolds number dependency is very weak, principally because the surface drag is predominantly form drag and the turbulence production process is at scales comparable to the roughness element sizes. LES is thus able to simulate turbulent flow over the urban-like obstacles at high Re with grids that would be far too coarse for adequate computation of corresponding smooth-wall flows. Comparison between LES and steady Reynolds-averaged Navier-Stokes (RANS) results are included, emphasising that the latter are inadequate, especially within the canopy region.