Turbulent flow in a circular separationless diffuser at Reynolds numbers smaller than 2000

The results of experimental and numerical investigation of flow in a circular conical diffuser with a small conicity angle ensuring separationless flow are presented. The measurements are carried out in an air flow with the Reynolds number Re₂ in the diffuser exit section ranging from 600 to 3000. A considerable effect of the channel expansion on the flow pattern is found to exist. It is shown that, as distinct from a tube, in which only laminar flow can be realized as steady for Re < 2000, in the exit section of a diffuser with the generator slope of 0.3° and a length equal to 70 entry diameters a developed turbulent flow is formed for Re₂ > 1000. For Re₂ > 1300 this flow is steady, that is, almost independent of the turbulence level at the entry, and is determined by the Reynolds number Re₂ in the exit section. For Re₂ ≈ 1000 the turbulent flow continuously goes over into a laminar flow. The flow parameters measured at the diffuser exit correspond to calculations in accordance with the threeequation turbulence model.     

Wall modeling for implicit large-eddy simulation and immersed-interface methods

We propose and analyze a wall model based on the turbulent boundary layer equations (TBLE) for implicit large-eddy simulation (LES) of high Reynolds number wall-bounded flows in conjunction with a conservative immersed-interface method for mapping complex boundaries onto Cartesian meshes. Both implicit subgrid-scale model and immersed-interface treatment of boundaries offer high computational efficiency for complex flow configurations. The wall model operates directly on the Cartesian computational mesh without the need for a dual boundary-conforming mesh. The combination of wall model and implicit LES is investigated in detail for turbulent channel flow at friction Reynolds numbers from Re _τ  = 395 up to Re _τ =100,000 on very coarse meshes. The TBLE wall model with implicit LES gives results of better quality than current explicit LES based on eddy viscosity subgrid-scale models with similar wall models. A straightforward formulation of the wall model performs well at moderately large Reynolds numbers. A logarithmic-layer mismatch, observed only at very large Reynolds numbers, is removed by introducing a new structure-based damping function. The performance of the overall approach is assessed for two generic configurations with flow separation: the backward-facing step at Re _h = 5,000 and the periodic hill at Re _H = 10,595 and Re _H = 37,000 on very coarse meshes. The results confirm the observations made for the channel flow with respect to the good prediction quality and indicate that the combination of implicit LES, immersed-interface method, and TBLE-based wall modeling is a viable approach for simulating complex aerodynamic flows at high Reynolds numbers. They also reflect the limitations of TBLE-based wall models.     

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.     

Experiments on the effect of acceleration on the drag of tapwater bubbles

Three-dimensional (3D) trajectories of spherical air bubbles passing through a converging part of a rectangular channel have been measured. Bubble diameters, d _b, were less than 1 mm and the Reynolds numbers, Re _b, for stagnant tapwater and for mean liquid velocity vˉ _L=0.25 m/s were in about same range. Received: 15 January 2001 / Accepted: 12 June 2001     

Influence of the surface roughness on the third-order moments of velocity fluctuations

Roughness wall effects in a zero pressure gradient turbulent boundary layers were investigated using hot-wire anemometry. The skewness and diffusion factors of u and v, the longitudinal and normal velocity fluctuations, were measured and represented using wall variables. The results indicate that the wall roughness removes the crossover point between sweep and ejection events to the outer region of the layer for a single Reynolds number Re _θ  > 3,000. This behaviour exhibits that the roughness surface favours the maintaining of sweep events obtained by a quadrant analysis. These results show that communication between the wall region and outer region of a turbulent boundary layer exists and the wall similarity hypothesis for a rough wall is questionable. The effect of the wall roughness on the position of the point crossover from sweep to ejection motions with respect to the wall seems to be the same as that obtained when the Reynolds number is higher. Received: 8 March 2000/Accepted: 15 May 2000     

Scalar transport from a point source in flows over wavy walls

Simultaneous measurements of the velocity and concentration field in fully developed turbulent flows over a wavy wall are described. The concentration field originates from a low-momentum plume of a passive tracer. PLIF and digital particle image velocimetry are used to make spatially resolved measurements of the structure of the scalar distribution and the velocity. The measurements are performed at three different Reynolds numbers of Re _b = 5,600, Re _b = 11,200 and Re _b = 22,400, respectively, based on the bulk velocity u _b and the total channel height 2h. The velocity field and the scalar field are investigated in a water channel with an aspect ratio of 12:1, where the bottom wall of the test section consists of a train of sinusoidal waves. The wavy wall is characterized by the amplitude to wavelength ratio α = 0.05 and the ratio β between the wave amplitude and the half channel height where β = 0.1. The scalar is released from a point source at the wave crest. For the concentration measurements, Rhodamine B is used as tracer dye. At low to moderate Reynolds number, the flow field is characterized through a recirculation zone which develops after the wave crest. The recirculation zone induces high intensities of the fluctuations of the streamwise velocity and wall-normal velocity. Furthermore, large-scale structures are apparent in the flow field. In previous investigations it has been shown that these large-scale structures meander laterally in flows over wavy bottom walls. The investigations show a strong effect of the wavy bottom wall on the scalar mixing. In the vicinity of the source, the scalar is transported by packets of fluid with a high scalar concentration. As they move downstream, these packets disintegrate into filament-like structures which are subject to strong gradients between the filaments and the surrounding fluid. The lateral scale of the turbulent plume is smaller than the lateral scale of the large-scale structures in the flow field and the plume dispersion is dominated by the structures in the flow field. Due to the lateral meandering of the large-scale structures of the flow field, also the scalar plume meanders laterally. Compared to turbulent plumes in plane channel flows, the wavy bottom wall enhances the mixing effect of the turbulent flow and the spreading rate of the scalar plume is increased.     

2D and 3D Turbulent Fluctuations in Open Channel Flow with Re τ = 590 Studied by Large Eddy Simulation

Well-resolved 3D Large Eddy Simulations (LES) are presented for open channel flow at a Reynolds number Re _τ  = 590 based on friction velocity u _τ and water depth h. The results are depth-averaged and thereby information is obtained on the 2D horizontal fluctuations in the channel. The total turbulence is decomposed into 2D and 3D fluctuations and the energy content of these as well as their spectral distribution is studied. It is found that only 15% of the fluctuating energy is contained in the 2D fluctuations and that these are mostly of scales larger than the water depth while the 3D fluctuations are restricted by the limited vertical extent of the water body and have scales smaller than the water depth. Information is obtained on the dispersion terms arising from the depth-averaging procedure, and scalar transport due to a vertical line source of tracer is studied thereby investigating the contribution of the 2D and 3D fluctuations to the transverse mixing of the scalar.     

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     

PIV–LES analysis of channel flow rotating about the streamwise axis

The flow field of a channel rotating about the streamwise axis is analyzed experimentally and numerically. The current investigations were carried out at a bulk velocity based Reynolds number of Re_m = 2850 and a friction velocity based Reynolds number of Re_τ = 180, respectively. Particle-image velocimetry (PIV) measurements are compared with large-eddy simulation data to show earlier direct numerical simulation findings to generate too large a reverse flow region in the center region of the spanwise flow. The development of the mean spanwise velocity distribution and the influence of the rotation on the turbulent properties, i.e., the Reynolds stresses and the two-point correlations of the flow, are confirmed in both investigations. The rotation primarily influences those components of the Reynolds shear stresses, which contain the spanwise velocity component. The size of the correlation areas and thus the length scales of the flow generally grow in all three coordinate directions leading to longer structures. Furthermore, experimental results of the same channel flow at a significantly lower bulk Reynolds number of Re_{m, l} = 665, i.e., a laminar flow in a non-rotating channel, are introduced. The experiments show the low Reynolds number flow to become turbulent under rotation and to develop the same characteristics as the high Reynolds number flow.     

Large eddy simulation of compressible channel flow

Astract The present study is a contribution to the analysis of wall-bounded compressible flows, including a special focus on wall modeling for compressible turbulent boundary layer in a plane channel. large eddy simulation (LES) of fully developed isothermal channel flows at Re = 3,000 and Re = 4,880 with a sufficient mesh refinement at the wall are carried out in the Mach number range 0.3 ≤ M ≤ 3 for two different source term formulations: first the classical extension of the incompressible configuration by Coleman et al. (J. Fluid Mech. 305:159–183, 1995), second a formulation presently derived to model both streamwise pressure drop and streamwise internal energy loss in a spatially developed compressible channel flow. It is shown that the second formulation is consistent with the spatial problem and yields a much stronger cooling effect at the wall than the classical formulation. Based on the present LES data bank, compressibility and low Reynolds number effects are analysed in terms of coherent structure and statistics. A study of the universality of the structure of the turbulence in non-hypersonic compressible boundary layers (M≤5) is performed in reference to Bradshaw (Annu. Rev. Fluid. Mech. 9:33–54, 1977). An improvement of the van Driest transformation is proposed; it accounts for both density and viscosity changes in the wall layer. Consistently, a new integral wall scaling (y ^c+) which accounts for strong temperature gradients at the wall is developed for the present non-adiabatic compressible flow. The modification of the strong Reynolds analogy proposed by Huang et al. (J. Fluid Mech. 305:185–218, 1995) to model the correlation between velocity and temperature for non-adiabatic wall layers is assessed on the basis of a Crocco–Busemann relation specific to channel flow. The key role of the mixing turbulent Prandtl number Pr _m is pointed out. Results show very good agreement for both source formulations although each of them involve a very different amount of energy transfer at the wall. The present work was performed within the framework of the French–German research initiative “large eddy simulation of complex flows’ (UR 507). The computing resources were provided by IDRIS-France. The authors gratefully acknowledge the financial support from the Centre National de la Recherche Scientifique (CNRS), the Centre d’été Mathématique de Recherche Avancée en Calcul Scientifique (CEMRACS) and the Direction Générale de l’Armement (DGA/D4S).     

Oscillating flow in a 2-D diffuser

Separating oscillating flows in an internal, adverse pressure gradient geometry are studied experimentally. Simultaneous velocity and pressure measurements demonstrate that the minor losses associated with oscillating flow in an adverse pressure gradient geometry can be smaller or larger than those for steady flow. Separation is found to begin high in the diffuser and propagate downward. The flow is able to remain attached further into the diffuser with larger Reynolds numbers, small displacement amplitudes, and smaller diffuser angles. The extent of separation grows with L ₀/h. The minor losses grow with increasing displacement amplitude in the measured range 10 < L ₀/h < 40. Losses decrease with increasing Re _δ in the measured range of 380 < Re _δ < 740. It is found that the losses increase with increasing diffuser angle over the measured range of 12° < θ < 30°. The nondimensional acoustic power dissipation increases with Reynolds number in the measured range and decreases with displacement amplitude.     

Effects of solid and slotted baffles on the convection characteristics of backward-facing step flow in a channel

This study investigates the enhancement of the laminar forced convection characteristics of backward-facing step flow in a two-dimensional channel through the installation of solid and slotted baffles onto the channel wall. The effects of the height of baffle H _b, inclination of baffle installation ϕ_b, height of slot in baffle H _t, inclination of slot in baffle ϕ_t, and distance between the backward-facing step and baffle D on the flow structure, temperature distribution and Nusselt number variation for the system at various Re are numerically explored. Results show that a slotted baffle can enhance the average Nusselt number for the heating section of channel plate by the maximum 190% when Pr=0.7, H _s=0.5, L=5, H _b ≤ 0.3, W _b ≤ 0.2, 0.1 ≤ D ≤ 0.5, 0° ≤ ϕ_b ≤ 45°, H _t ≤ 0.1, 0° ≤ ϕ_t ≤ 45° and 50 ≤ Re ≤ 400. As for the solid baffle, the enhancement may be up by 230%. The solid baffle might cause the re-separation of main stream, and consequently result in poor local heat transfer coefficient in the end region of heating section. This disadvantage can be obviously improved as the baffle is slotted. Besides the penalty of increase in pressure drop due to the baffle installation is much higher for the situation with solid baffle.     

Turbulent channel flow with 2D wedges of random height on one wall

Direct numerical simulations (DNS) of a turbulent channel flow with 2D wedges of random height on the bottom wall have been performed. In addition, two other simulations have been carried out to assess the effect of the geometry on the overlying flow. In the first simulation, the four smallest elements were removed while in the other, a uniform distribution of wedges with the same area was used. Two Reynolds numbers were studied, Re_b=2500 and Re_b=5000 which correspond in case of smooth walls to Re_τ=180 and 300, respectively. Roughness on the wall induces separated regions, the reattachment occurring on the walls of the wedges or on the bottom wall. The pressure gradients on the walls increase the ejections and inrushes towards the wall. As a consequence the flow is more isotropic. The mechanism inducing an improved isotropy has been explained in term of the spectra and budgets of Reynolds stress. The comparison of the 3 surfaces has shown that near the wall, the uniformly distributed roughness represents only a poor approximation of the surface with wedges of random height. The Reynolds stresses, pressure distribution and spectra on the modified wall agree well with those on the random surface. Energy spectra show the pitch to height ratio of the largest elements to be the more appropriate geometrical parameter to describe the geometry.     

Large Eddy Simulations: How to evaluate resolution

The present paper gives an analysis of fully developed channel flow at Reynolds number of Re=u_τδ/ν=4000 based on the friction velocity, u_τ, and half the channel height, δ. Since the Reynolds number is high, the LES is coupled to a URANS model near the wall (hybrid LES–RANS) which acts as a wall model. It it found that the energy spectra is not a good measure of LES resolution; neither is the ratio of the resolved turbulent kinetic energy to the total one (i.e. resolved plus modelled turbulent kinetic energy). It is suggested that two-point correlations are the best measures for estimating LES resolution. It is commonly assumed that SGS dissipation takes place at high wavenumbers. Energy spectra of the fluctuating velocity gradients show that this is not true; the major part of the SGS dissipation takes place at low to midrange wavenumbers. Furthermore, the energy spectra of the fluctuating velocity gradients reveals that the accuracy of the predicted velocity gradients at the highest resolved wavenumbers is very poor.     

Turbulent flow in a channel at a low Reynolds number

A two-component laser Doppler velocimeter with high spatial and temporal resolution was used to obtain measurements for fully developed turbulent flow of water through a channel with an aspect ratio of 12 : 1 at Re=5700 (based on the centerline velocity and the half-height of the channel). Statistical quantities that were determined are the mean streamwise velocity, the root-mean-square of the fluctuations of the streamwise and the normal velocities, the Reynolds shear stress and higher order moments. Turbulence production is calculated from these quantities. Turbulence statistics obtained from experiments are compared with results from a direct numerical simulation at the same Reynolds number. The good agreement validates a recent DNS, at Re=5700, which is approximately twice as large as used in most previous studies. Received: 12 May 1997 / Accepted: 8 April 1998     

Development of a turbulent boundary layer after a step from smooth to rough surface

The flow developing downstream of a step change from smooth to rough surface condition is studied in the light of Townsend’s wall similarity hypothesis. Previous studies seem to support the hypothesis for channel and pipe flows, but there are considerable controversies about its application to boundary layers and in particular to surface roughness formed by spanwise bars. It has been suggested that this controversy arises from insufficient separation of scales between the boundary layer thickness and the roughness length scale. An experimental investigation has therefore been undertaken where the flow evolves from a fully developed smooth wall boundary layer at high Reynolds numbers over a step in surface roughness (Re _θ = 13,400 at the step). The flow is mapped through the development of the internal layer until the flow is fully developed over the rough wall. The internal layer is found to grow as δ ∼ X ^0.73, and after about 15 boundary layer thicknesses at the step, the internal layer has reached the outer edge of the incoming layer. At the last rough wall measurement station, the Reynolds number has grown to Re _θ ≈ 32,600 and the ratio of boundary layer to roughness length scales is δ/k ≈ 140. The outer layer differences between the smooth and the rough wall data were found to be sufficiently small to conclude that for this setup the Townsend’s wall similarity hypothesis appears to hold.     

Three-dimensional Turbulent Boundary Layer in a Shrouded Rotating System

A thre-dimensional direct numerical simulation is combined with a laboratory study to describe the turbulent flow in an enclosed annular rotor-stator cavity characterized by a large aspect ratio G = (b − a)/h = 18.32 and a small radius ratio a/b = 0.152, where a and b are the inner and outer radii of the rotating disk and h is the interdisk spacing. The rotation rate Ω considered is equivalent to the rotational Reynolds number Re = Ωb ²/ν= 9 .5 × 10⁴ (ν the kinematic viscosity of water). This corresponds to a value at which experiment has revealed that the stator boundary layer is turbulent, whereas the rotor boundary layer is still laminar. Comparisons of the computed solution with velocity measurements have given good agreement for the mean and turbulent fields. The results enhance evidence of weak turbulence by comparing the turbulence properties with available data in the literature (Lygren and Andersson, J Fluid Mech 426:297–326, 2001). An approximately self-similar boundary layer behavior is observed along the stator. The wall-normal variations of the structural parameter and of characteristic angles confirm that this boundary layer is three-dimensional. A quadrant analysis (Kang et al., Phys Fluids 10:2315–2322, 1998) of conditionally averaged velocities shows that the asymmetries obtained are dominated by Reynolds stress-producing events in the stator boundary layer. Moreover, Case 1 vortices (with a positive wall induced velocity) are found to be the major source of generation of special strong events, in agreement with the conclusions of Lygren and Andersson (J Fluid Mech 426:297–326, 2001).     

An efficient very large eddy simulation model for simulation of turbulent flow

Among the various hybrid methodologies, Speziale's very large eddy simulation (VLES) is one that was proposed very early. It is a unified simulation approach that can change seamlessly from Reynolds Averaged Navier–Stokes (RANS) to direct numerical simulation (DNS) depending on the numerical resolution. The present study proposes a new improved variant of the original VLES model. The advantages are achieved in two ways: (i) RANS simulation can be recovered near the wall which is similar to the detached eddy simulation concept; (ii) a LES subgrid scale model can be reached by the introduction of a third length scale, that is, the integral turbulence length scale. Thus, the new model can provide a proper LES mode between the RANS and DNS limits. This new methodology is implemented in the standard k ? ? model. Applications are conducted for the turbulent channel flow at Reynolds number of Re_τ = 395, periodic hill flow at Re = 10,595, and turbulent flow past a square cylinder at Re = 22,000. In comparison with the available experimental data, DNS or LES, the new VLES model produces better predictions than the original VLES model. Furthermore, it is demonstrated that the new method is quite efficient in resolving the large flow structures and can give satisfactory predictions on a coarse mesh. Copyright © 2012 John Wiley & Sons, Ltd.     

Non-linear dynamics of semi-dilute polydisperse polymer solutions in microfluidics: effects of flow geometry

The non-linear dynamics of a semi-dilute (c/c* = 15) polydisperse polyethylene oxide (PEO) solution in microfluidics are studied experimentally using benchmark contraction–expansion flow geometries with three contraction–expansion ratios (4:1:4, 8:1:8 and 16:1:16) and two narrow channel lengths (L _c/D _h = 53 and 5.3, where L _c is the length of the narrow channel and D _h is its hydraulic diameter). Complex flows over a range of elasticity numbers (El), Weissenberg numbers (Wi) and Reynolds numbers (Re) are characterized using micro-particle image velocimetry ( \upmu\upmu-PIV) and pressure drop measurements. The evolution of vortex formation and dynamics has been visualized through a step-flow-rate experiment. Various flow dynamics regimes have been quantified and are presented in a Wi–Re diagram. The experimental results reveal that the contraction ratio can result in qualitatively different vortex dynamics of semi-dilute polymer solutions in microfluidics, whereas the length of the narrow channel merely affects the dynamics at a quantitative level. A single elasticity number, if defined by the size of the narrow channel, is not sufficient to account for the effects of contraction ratio on the non-linear vortex dynamics.     

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