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.     

Development and characterization of a variable turbulence generation system

Experimental turbulent combustion studies require systems that can simulate the turbulence intensities [u′/U ₀ ~ 20–30% (Koutmos and McGuirk in Exp Fluids 7(5):344–354, 1989)] and operating conditions of real systems. Furthermore, it is important to have systems where turbulence intensity can be varied independently of mean flow velocity, as quantities such as turbulent flame speed and turbulent flame brush thickness exhibit complex and not yet fully understood dependencies upon both U ₀ and u′. Finally, high pressure operation in a highly pre-heated environment requires systems that can be sealed, withstand high gas temperatures, and have remotely variable turbulence intensity that does not require system shut down and disassembly. This paper describes the development and characterization of a variable turbulence generation system for turbulent combustion studies. The system is capable of a wide range of turbulence intensities (10–30%) and turbulent Reynolds numbers (140–2,200) over a range of flow velocities. An important aspect of this system is the ability to vary the turbulence intensity remotely, without changing the mean flow velocity. This system is similar to the turbulence generators described by Videto and Santavicca (Combust Sci Technol 76(1):159–164, 1991) and Coppola and Gomez (Exp Therm Fluid Sci 33(7):1037–1048, 2009), where variable blockage ratio slots are located upstream of a contoured nozzle. Vortical structures from the slots impinge on the walls of the contoured nozzle to produce fine-scale turbulence. The flow field was characterized for two nozzle diameters using three-component Laser Doppler velocimetry (LDV) and hotwire anemometry for mean flow velocities from 4 to 50 m/s. This paper describes the key design features of the system, as well as the variation of mean and RMS velocity, integral length scales, and spectra with nozzle diameter, flow velocity, and turbulence generator blockage ratio.     

Interaction of turbulent natural convection and surface thermal radiation in inclined square enclosures

Two-dimensional numerical studies of flow and temperature fields for turbulent natural convection and surface radiation in inclined differentially heated enclosures are performed. Investigations are carried out over a wide range of Rayleigh numbers from 10⁸ to 10¹², with the angle of inclination varying between 0° and 90°. Turbulence is modeled with a novel variant of the k–ε closure model. The predicted results are validated against experimental and numerical results reported in literature. The effect of the inclination of the enclosure on pure turbulent natural convection and the latter’s interaction with surface radiation are brought out. Profiles of turbulent kinetic energy and effective viscosity are studied to observe the net effect on the intensity of turbulence caused by the interaction of natural convection and surface radiation. The variations of local Nusselt number and average Nusselt number are presented for various inclination angles. Marked change in the convective Nusselt number is found with the orientation of enclosure. Also analyzed is the influence of change in emissivity on the flow and heat transfer. A correlation relevant to practical applications in the form of average Nusselt number, as a function of Rayleigh number, Ra, radiation convection parameter, N _RC and inclination angle of the enclosure, φ is proposed.     

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).     

'T-RANS' Simulation of Deterministic Eddy Structure in Flows Driven by Thermal Buoyancy and Lorentz Force

The paper reports on the application of the Time-dependent Reynolds-Averaged Navier–Stokes (T-RANS) approach to analysing the effects of magnetic force and bottom-wall configuration on the reorganisation of a large coherent structure and its role in the transport processes in Rayleigh–Bénard convection. The large-scale deterministic motion is fully resolved in time and space, whereas the unresolved stochastic motion is modelled by a `subscale' model for which the conventional algebraic stress/flux expressions were used, closed with the low-Re number (k)-(ε)-(θ²) three-equation model. The applied method reproduces long-term averaged mean flow properties, turbulence second moments, and all major features of the coherent roll/cell structure in classic Rayleigh–Bénard convection in excellent agreement with the available DNS and experimental results. Application of the T-RANS approach to Rayleigh–Bénard convection with wavy bottom walls and a superimposed magnetic field yielded the expected effects on there organisation of the eddy structure and consequent modifications of the mean and turbulence parameters and wall heat transfer. This revised version was published online in August 2006 with corrections to the Cover Date.     

Stability for the 3D Navier–Stokes Equations with Nonzero far Field Velocity on Exterior Domains

Concerning to the non-stationary Navier–Stokes flow with a nonzero constant velocity at infinity, just a few results have been obtained, while most of the results are for the flow with the zero velocity at infinity. The temporal stability of stationary solutions for the Navier–Stokes flow with a nonzero constant velocity at infinity has been studied by Enomoto and Shibata (J Math Fluid Mech 7:339–367, 2005), in L ^p spaces for p ≥ 3. In this article, we first extend their result to the case \frac32 < p{\frac{3}{2} < p} by modifying the method in Bae and Jin (J Math Fluid Mech 10:423–433, 2008) that was used to obtain weighted estimates for the Navier–Stokes flow with the zero velocity at infinity. Then, by using our generalized temporal estimates we obtain the weighted stability of stationary solutions for the Navier–Stokes flow with a nonzero velocity at infinity.     

Numerical investigation of turbulent flow,heat transfer and entropy generation in a helical coiled tube with larger curvature ratio

Three-dimensional turbulent forced convective heat transfer and flow characteristics, and the non-dimensional entropy generation number in a helical coiled tube subjected to uniform wall temperature are simulated using the k–ε standard turbulence model. A finite volume method is employed to solve the governing equations. The effects of Reynolds number, curvature ratio, and coil pitch on the average friction factor and Nusselt number are discussed. The results presented in this paper cover a Reynolds number range of 2 × 10⁴ to 6 × 10⁴, a pitch range of 0.1–0.2 and a curvature ratio range of 0.1–0.3. The results show that the coil pitch, curvature ratio and Reynolds number have different effects on the average friction factor and Nusselt number at different cross-sections. In addition, the flow and heat transfer characteristics in a helical coiled tube with a larger curvature ratio for turbulent flow are different from that of smaller curvature ratio for laminar and turbulent flow in certain ways. Some new features that are not obtained in previous researches are revealed. Moreover, the effects of Reynolds number, curvature ratio, and coil pitch on the non-dimensional entropy generation number of turbulent forced convection in a helical coiled tube are also discussed.     

Characterizing developing adverse pressure gradient flows subject to surface roughness

An experimental study was conducted to examine the effects of surface roughness and adverse pressure gradient (APG) on the development of a turbulent boundary layer. Hot-wire anemometry measurements were carried out using single and X-wire probes in all regions of a developing APG flow in an open return wind tunnel test section. The same experimental conditions (i.e., T _∞, U _ref, and C _p) were maintained for smooth, k ⁺ = 0, and rough, k ⁺ = 41–60, surfaces with Reynolds number based on momentum thickness, 3,000 < Re _θ < 40,000. The experiment was carefully designed such that the x-dependence in the flow field was known. Despite this fact, only a very small region of the boundary layer showed a balance of the various terms in the integrated boundary layer equation. The skin friction computed from this technique showed up to a 58% increase due to the surface roughness. Various equilibrium parameters were studied and the effect of roughness was investigated. The generated flow was not in equilibrium according to the Clauser (J Aero Sci 21:91–108, 1954) definition due to its developing nature. After a development region, the flow reached the equilibrium condition as defined by Castillo and George (2001), where Λ = const, is the pressure gradient parameter. Moreover, it was found that this equilibrium condition can be used to classify developing APG flows. Furthermore, the Zagarola and Smits (J Fluid Mech 373:33–79, 1998a) scaling of the mean velocity deficit, U _∞δ*/δ, can also be used as a criteria to classify developing APG flows which supports the equilibrium condition of Castillo and George (2001). With this information a ‘full APG region’ was defined.     

Dynamical characteristics of coherent structure in the outer region of turbulence boundary layer

In the present paper the coherent structures in the outer region of turbulent boundary layer were investigated experimentally and analytically. From the observation of the flow field over smooth wall, rough wall and sand wave wall, it was found that the direct effect of wall on the flow structure can reachy ⁺¹≈100, and both lateral and vertical vortices exist in the outer region, but the coherent structures in the outer region are mainly the formation, development and decay of the large-scale lateral vortices. By experimental and dynamical analysis, some influence factors and their relations associated with the dynamical process of lateral vortices were deduced. The project supported by the National Natural Science Foundation of China     

Boundary layer structure of oscillatory open-channel shallow flows over smooth and rough beds

The boundary layer structure of oscillatory shallow open channel flows has been studied in a wide flume. Fluorescence solution was released at a porous rough bed through a diffuser covered by gravel of 0.5 cm grain size. A planar laser-induced fluorescence (PLIF) system was used to visualise the dye plumes in both vertical and horizontal planes for a qualitative understanding of the roles of large-scale flow structures in mass transport. A variety of tests were conducted for a range of oscillatory periods (30–240 s), water depths (3–16 cm) and velocity amplitudes (0.027–0.325 m/s), which cover a wide range of oscillatory flows with Reynolds numbers Re _a varied from 0.3 × 10⁴ (laminar) to 2.1 × 10⁶ (fully turbulent). For quantitative investigation, a novel technique, namely combined laser-induced fluorescence (LIF) and 2D laser Doppler velocimetry (LDV) (LIF/LDV), was developed and used to measure the velocity and solute concentration simultaneously in a vertical plane over 50 cycles. From the dye plumes revealed by the PLIF in transitional flows, there are different patterns of flow structure and solute transport with three representative stages of acceleration, deceleration and flow reversal. In the acceleration stage, turbulence was suppressed with dye layers adhering to the surface with little vertical mass transport. In the deceleration stage, flame-like turbulent structures occurred when turbulence generation was prominent. This was investigated quantitatively by recording the percentage occurrence of the adhered smooth layers per cycle. For those smooth bed cases with Re _a < 1.8 × 10⁵, the adhered smooth dye layers type of boundary layer occupied 100% of the oscillation period. Over a sufficiently high Re _a , a rough bed can generate fully turbulent oscillatory flows without the appearance of adhering dye layers. Between these two extremes, a transitional flow regime occurs in a wide range of flow conditions: Re _a > 2.7 × 10⁴ over the rough bed and Re _a > 8.3 × 10⁶ over a smooth bed.     

Large-scale flow structures and heat transport of turbulent forced and mixed convection in a closed rectangular cavity

Results of an experimental study of flow structure formation and heat transport in turbulent forced and mixed convection are presented. The experiments were conducted in a rectangular cavity with a square cross section, which has an aspect ratio between length and height of Γ_xz = 5. Air at atmospheric pressure was used as working fluid. The air inflow was supplied through a slot below the ceiling, while exhausting was provided by another slot, which is located directly above the floor. Both vents extend over the whole length of the cell. In order to induce thermal convection the bottom of the cell is heated while the ceiling is maintained at a constant temperature. This configuration allows to generate and study mixed convection under well defined conditions. Results of forced convection at Re = 1.07 × 10⁴ as well as mixed convection at 1.01 × 10⁴ ? Re ? 3.4 × 10⁴ and Ra = 2.4 × 10⁸ (3.3 ? Ar ? 0.3), which were obtained by means of Particle Image Velocimetry and local temperature measurements, are presented. For purely forced convection a 2D mean wind, which can be approximated by a solid body rotation, is found. With increasing Archimedes number this structure becomes unstable, leading to a transition of the solid body rotation into additional smaller convection rolls. Proper orthogonal decomposition of the instantaneous velocity fields has been performed for further analysis of these coherent large-scale structures. Their fingerprint is found in the spatial temperature distribution of the out flowing air at the end of the outlet channel, which reveals a temporally stable profile with two maxima over the length of the outlet. Moreover a maximum in the global enthalpy transport by the fluid is found at Ar ≈ 0.6.     

Statistically Steady Incompressible DNS to Validate a New Correlation for Turbulent Burning Velocity in Turbulent Premixed Combustion

Incompressible 3-D DNS is performed in non-decaying turbulence with single step chemistry to validate a new analytical expression for turbulent burning velocity. The proposed expression is given as a sum of laminar and turbulent contributions, the latter of which is given as a product of turbulent diffusivity in unburned gas and inverse scale of wrinkling at the leading edge. The bending behavior of U _T at higher u′ was successfully reproduced by the proposed expression. It is due to decrease in the inverse scale of wrinkling at the leading edge, which is related with an asymmetric profile of FSD with increasing u′. Good agreement is achieved between the analytical expression and the turbulent burning velocities from DNS throughout the wrinkled, corrugated and thin reaction zone regimes. Results show consistent behavior with most experimental correlations in literature including those by Bradley et al. (Philos Trans R Soc Lond A 338:359–387, 1992), Peters (J Fluid Mech 384:107–132, 1999) and Lipatnikov et al. (Progr Energ Combust Sci 28:1–74, 2002).     

Turbulence beneath finite amplitude water waves

Babanin and Haus (J Phys Oceanogr 39:2675–2679, 2009) recently presented evidence of near-surface turbulence generated below steep non-breaking deep-water waves. They proposed a threshold wave parameter a ²ω/ν = 3,000 for the spontaneous occurrence of turbulence beneath surface waves. This is in contrast to conventional understanding that irrotational wave theories provide a good approximation of non-wind-forced wave behaviour as validated by classical experiments. Many laboratory wave experiments were carried out in the early 1960s (e.g. Wiegel 1964). In those experiments, no evidence of turbulence was reported, and steep waves behaved as predicted by the high-order irrotational wave theories within the accuracy of the theories and experimental techniques at the time. This contribution describes flow visualisation experiments for steep non-breaking waves using conventional dye techniques in the wave boundary layer extending above the wave trough level. The measurements showed no evidence of turbulent mixing up to a value of a ²ω/ν = 7,000 at which breaking commenced in these experiments. These present findings are in accord with the conventional understandings of wave behaviour.     

Turbulent Flow Produced by Piston Motion in a Spark-ignition Engine

Turbulence produced by the piston motion in spark-ignition engines is studied by 2D axisymmetric numerical simulations in the cylindrical geometry as in the theoretical and experimental work by Breuer et al. (Flow Turbul Combust 74:145, 2005). The simulations are based on the Navier–Stokes gas-dynamic equations including viscosity, thermal conduction and non-slip at the walls. Piston motion is taken into account as a boundary condition. The turbulent flow is investigated for a wide range of the engine speed, 1,000–4,000 rpm, assuming both zero and non-zero initial turbulence. The turbulent rms-velocity and the integral length scale are investigated in axial and radial directions. The rms-turbulent velocity is typically an order-of-magnitude smaller than the piston speed. In the case of zero initial turbulence, the flow at the top-dead-center may be described as a combination of two large-scale vortex rings of a size determined by the engine geometry. When initial turbulence is strong, then the integral turbulent length demonstrates self-similar properties in a large range of crank angles. The results obtained agree with the experimental observations of Breuer et al. (Flow Turbul Combust 74:145, 2005).     

Large-Eddy Simulation of Dispersion from Line Sources in a Turbulent Channel Flow

Large-eddy simulations of the dispersion from scalar line sources at various locations within a fully developed turbulent channel flow at Re = uh/ν = 10,400 are presented. Both mean and fluctuating scalar quantities are compared with those from the single available set of experimental data (Lavertu and Mydlarski, J Fluid Mech 528:135–172, 2005) and differences are highlighted and discussed. The results are also discussed in the context of scalar dispersion in other kinds of turbulent flows, e.g. homogeneous shear-flow. Initial computations at a much lower Reynolds number are also reported and compared with the two available direct numerical simulation data sets.     

Testing of Model Equations for the Mean Dissipation using Kolmogorov Flows

Direct Numerical Simulations (DNS) of Kolmogorov flows are performed at three different Reynolds numbers Re _λ between 110 and 190 by imposing a mean velocity profile in y-direction of the form U(y) = F sin(y) in a periodic box of volume (2π)³. After a few integral times the turbulent flow turns out to be statistically steady. Profiles of mean quantities are then obtained by averaging over planes at constant y. Based on these profiles two different model equations for the mean dissipation ε in the context of two-equation RANS (Reynolds Averaged Navier–Stokes) modelling of turbulence are compared to each other. The high Reynolds number version of the k-ε-model (Jones and Launder, Int J Heat Mass Transfer 15:301–314, 1972), to be called the standard model and a new model by Menter et al. (2006), to be called the Menter–Egorov model, are tested against the DNS results. Both models are solved numerically and it is found that the standard model does not provide a steady solution for the present case, while the Menter–Egorov model does. In addition a fairly good quantitative agreement of the model solution and the DNS data is found for the averaged profiles of the kinetic energy k and the dissipation ε. Furthermore, an analysis based on flow-inherent geometries, called dissipation elements (Wang and Peters, J Fluid Mech 608:113–138, 2008), is used to examine the Menter–Egorov ε model equation. An expression for the evolution of ε is derived by taking appropriate moments of the equation for the evolution of the probability density function (pdf) of the length of dissipation elements. A term-by-term comparison with the model equation allows a prediction of the constants, which with increasing Reynolds number approach the empirical values.     

The nature of near-wall convection velocity in turbulent channel flow

A novel notion of turbulent structure the local cascade structure-is introduced to study the convection phenomenon in a turbulent channel flow. A space-time cross-correlation method is used to calculate the convection velocity. It is found that there are two characteristic convection speeds near the wall, one associated with small-scale streaks of a lower speed and another with streamwise vortices and hairpin vortices of a higher speed. The new concept of turbulent structure is powerful to illustrate the dominant role of coherent structures in the near-wall convection, and to reveal also the nature of the convection-the propagation of patterns of velocity fluctuations-which is scale-dependent.     

Coherent Vortex Simulation (CVS), A Semi-Deterministic Turbulence Model Using Wavelets

In the spirit of Ha Minh's semi-deterministic model, we propose a new method for computing fully-developed turbulent flows, called Coherent Vortex Simulation (CVS). It is based on the observation that turbulent flows contain both an organized part, the coherent vortices, and a random part, the incoherent background flow. The separation into coherent and incoherent contributions is done using the wavelet coefficients of the vorticity field and the Biot–Savart kernel to reconstruct the coherent and incoherent velocity fields. The evolution of the coherent part is computed using a wavelet basis, adapted at each time step to resolve the regions of strong gradients, while the incoherent part is discarded during the flow evolution, which models turbulent dissipation. The CVS method is similar to LES, but it uses nonlinear multiscale band-pass filters, which depend on the instantaneous flow realization, while LES uses linear low-pass filters, which do not adapt to the flow evolution. As example, we apply the CVS method to compute a time developing two-dimensional mixing layer and a wavelet forced two-dimensional homogeneous isotropic flow. We also demonstrate how walls or obstacles can be taken into account using penalization and compute a two-dimensional flow past an array of cylinders. Finally, we perform the same segmentation into coherent and incoherent components in a three-dimensional homogeneous isotropic turbulent flow. We show that the coherent components correspond to vortex tubes, which exhibit non-Gaussian statistics and long-range correlation, with the same k ^−5/3power-law energy spectrum as the total flow. In contrast, the incoherent components correspond to an homogeneous random background flow which does not contain organized structures and presents an energy equipartition together with a Gaussian PDF of velocity. This justifies their elimination during the CVS computation to model turbulent dissipation. This revised version was published online in August 2006 with corrections to the Cover Date.     

Very high Reynolds number boundary layers over 3D sparse roughness and obstacles: the mean flow

Hot-wire and oil-film interferometry measurements are taken for 3D rough wall boundary layers at very high Reynolds numbers (61,000 < Re θ < 120,000) with low blockage ratios, 10 < δ/H < 135, and high roughness, 100 < H ⁺ < 4,900. The results cover flows over both rough walls and over obstacles and are compared with and provide extension to recent lower Reynolds number results. The validity of the Townsend ‘wall similarity hypothesis’ in relation to consistently increasing 3D roughness is interrogated. In agreement with recent work, Schultz and Flack (J Fluid Mech 580:381–405, 2007) and Castro (J Fluid Mech 585:469–485, 2007) found that, for relatively low roughness, Townsend’s hypothesis holds for the mean velocity field. With increasing roughness, the equilibrium layer diminishes and gradually vanishes. The viscous component of the wall shear stress decreases, while the turbulent component increases as the roughness effects extend across the boundary layer.