Visualization of turbulent flow around a sphere at subcritical reynolds numbers

The shear layer evolution and turbulent structure of near-wake behind a sphere atRe= 11,000 and 5,300 were investigated using a smoke-wire visualization method. A laminar flow separation was found to occur near the equator. The smooth laminar shear layers appeared to be axisymmetrically stable to the downstream location of aboutx/d=1.0 atRe=11,000 andx/d= 1.7∼1.8 atRe=5,300, respectively. At Re=11,000, the vortex ring-shaped protrusions were observed with the onset of shear layer instability. Moreover, the transition from laminar to turbulence in the separated flow region occurred earlier at the hiher Reynolds number ofRe=11,000 than atRe=5,300. The PIV measurements in the streamwise and cross-sectional planes atRe=11,000 clearly revealed the turbulent structures of the sphere wake such as recirculating flow, shear layer instability, vortex roll-up, and small-scale turbulent eddies.     

Strengthened opposition control for skin-friction reduction in wall-bounded turbulent flows

An opposition control scheme with strengthened control input is proposed and tested in turbulent channel flows at friction Reynolds number Re_τ = 180 by direct numerical simulations. When the detection plane is located at less than 20 wall units, the drag reduction rate can be greatly enhanced by increasing the control amplitude parameter. The maximum drag reduction rate achieved in the present study is around 33%, which is much higher than the best value of 25% reported in literature. The strengthened control can be more efficient to attain a given drag reduction rate. Based on the total shear stress at the virtual wall established between the real wall and the detection plane by the control, a new friction velocity is proposed and the corresponding coordinate transform is made. Scaled by the proposed friction velocity, the wall-normal velocity fluctuation and the Reynolds shear stress of the controlled flows are collapsed well with those of the uncontrolled flow in the new coordinate. Based on the similarity, a relation between drag reduction rate and the effectiveness of the virtual wall is deduced, which disclosed that the elevation and residual Reynolds shear stress at the virtual wall are the key parameters to determine the drag reduction rate. The conclusion are also validated at Re_τ = 395 and 590. The decrease of the drag reduction rate with the increase of the Reynolds number is attributed to the enhanced residual Reynolds shear stress at the virtual wall.     

Effect of vertical large eddy breakup devices on the drag of a flat plate

Experiments are performed to study the possibility of decreasing the net drag of a flat plate with the use of streamwise-oriented vertical elements mounted normal to the surface in an incompressible equilibrium turbulent boundary layer. The Reynolds number based on the boundary-layer momentum thickness in the section where the vertical elements are placed is 820. It is demonstrated that vertical large eddy breakup elements with the geometry used do not reduce the drag of a flat plate in the major part of the range of Reynolds numbers Re _x examined. It is only at extremely low values of Re _x that a certain gain in the net drag is reached, as compared with the value for a non-modified flow.     

雷诺数对圆形渐缩喷嘴湍流射流的影响

本文研究雷诺数(Re)对圆形渐缩喷嘴湍流射流的影响.实验在射流出口雷诺数为 Re = 4050—20100 的范围内进行,分别测量了射流出口、中心线的平均及湍流流场以及部分径向剖面速度分布.所有测量均采用单热线恒温热线风速仪进行高频采样,所测流场范围在轴向上为 0—30d(这里d为射流出口直径).虽然出口速度分布均为"平顶帽"形,但测量结果依然反映出Re对射流出口以及下游流场有强烈的影响.当Re小于临界值(～10000)时 关键词： 雷诺数 圆形射流 热线风速仪     

Viscous Interaction Between a Vortex Tube and a Rotating Spherical Particle

The transient advection of a cylindrical vortex tube in a viscous incompressible flow field and its interaction with a rotating/spinning spherical particle has been investigated numerically at Reynolds numbers in the range of 20≤ Re≤200 for angular velocities of 0≤Ω≤0.5. The effects of vortex parameters such as size, circulation strength and initial position relative to the particle, on the temporal behavior of the lift and drag forces are studied. Vortex‐sphere interactions bring about major changes in the flow field particularly when coupled with particle rotation. It is observed that the forces acting on the particle are significantly influenced during the time that the vortex core is in the vicinity of the particle. The extent of these local changes are about ±30% in the drag coefficient and about ±200% in the lift coefficient as compared to flow over a rotating solid sphere with no vortex interaction. It is also found that a vortex with core radius between one and two particle diameters creates the strongest temporal variations in the lift and drag coefficients. Furthermore, maximum lift variations occur for the vortex‐particle head on collision, while a vortex with an offset distance of about one diameter from the principal flow axis generates the maximum drag variations.     

Investigating asymptotic suction boundary layers using a one-dimensional stochastic turbulence model

The turbulent asymptotic suction boundary layer is studied using a one-dimensional turbulence (ODT) model. ODT is a fully resolved, unsteady stochastic simulation technique. While flow properties reside on a one-dimensional domain, turbulent advection is represented using mapping events whose occurrences are governed by a random process. Due to its reduced spatial dimensionality, ODT achieves major cost reductions compared to three-dimensional (3D) simulations. A comparison to recent direct numerical simulation (DNS) data at moderate Reynolds number (Re = u_∞ / v₀ = 333, where u_∞ and v₀ are the free stream and suction velocity, respectively) suggests that the ODT model is capable of reproducing several velocity statistics, i.e. mean velocity and turbulent kinetic energy budgets, while peak turbulent stresses are under-estimated by ODT. Variation of the Reynolds number in the range Re ∈ [333,400,500,1000] shows that ODT can reproduce various trends observed as a result of increased suction in turbulent asymptotic suction boundary layers, i.e. the reduction of Reynolds stresses and enhanced skin friction. While up to Re = 500 our results can be directly compared to recent LES data, the simulation at Re = 1000 is currently not feasible through full 3D simulations, hence ODT may assist the design of future DNS or LES simulations at larger Reynolds numbers.     

Harbingers and latecomers – the order of appearance of exact coherent structures in plane Poiseuille flow

The transition to turbulence in plane Poiseuille flow (PPF) is connected with the presence of exact coherent structures. We here discuss a variety of different structures that are relevant for the transition, compare the critical Reynolds numbers and optimal wavelengths for their appearance, and explore the differences between flows operating at constant mass flux or at constant pressure drop. The Reynolds numbers quoted here are based on the mean flow velocity and refer to constant mass flux. Reynolds numbers based on constant pressure drop are always higher. The Tollmien–Schlichting (TS) waves bifurcate subcritically from the laminar profile at Re = 5772 at wavelength 6.16 and reach down to Re = 2610 at a different optimal wave length of 4.65. Their streamwise localised counter part bifurcates at the even lower value Re = 2334. Three-dimensional exact solutions appear at much lower Reynolds numbers. We describe one exact solutions that has a critical Reynolds number of 316. Streamwise localised versions of this state require higher Reynolds numbers, with the lowest bifurcation occurring near Re = 1018. The analysis shows that the various branches of TS-waves cannot be connected with transition observed near Re ≈ 1000 and that the exact coherent structures related to downstream vortices come in at lower Reynolds numbers and prepare for the transition.     

Generation of vortices by a streamwise oscillating cylinder

The wake of a streamwise oscillating cylinder is presently investigated. The Reynolds number investigated is 300, based on the cylinder diameterd. The cylinder oscillates at an amplitude of 0.5d and a frequencyf _e/f_s=1.8, wheref _e is the cylinder oscillating frequency andf _s is the natural vortex shedding frequency of a stationary cylinder. Under these conditions the flow is essentially two dimensional. A two-dimensional direct numerical simulation (DNS) scheme has been developed to calculate the flow. The DNS results display a street of binary vortices, each containing two counter-rotating vortical structures, symmetrical about the centerline, which is in excellent agreement with measurements. The drag and lift on the cylinder have been examined. The time averaged drag and lift are 1.4 and 0, respectively, which are the same as those on a stationary cylinder at the sameRe. However, the fluctuating drag was high, about 2.68. It has been found that, being symmetrically formed about the centerline, the binary vortices induce an essentially zero fluctuating lift, which may have a profound implication in flow control and engineering.     

Vorticity,backscatter and counter-gradient transport predictions using two-level simulation of turbulent flows

The two-level simulation (TLS) method evolves both the large-and the small-scale fields in a two-scale approach and has shown good predictive capabilities in both isotropic and wall-bounded high Reynolds number (Re) turbulent flows in the past. Sensitivity and ability of this modelling approach to predict fundamental features (such as backscatter, counter-gradient turbulent transport, small-scale vorticity, etc.) seen in high Re turbulent flows is assessed here by using two direct numerical simulation (DNS) datasets corresponding to a forced isotropic turbulence at Taylor’s microscale-based Reynolds number Re_λ ≈ 433 and a fully developed turbulent flow in a periodic channel at friction Reynolds number Re_τ ≈ 1000. It is shown that TLS captures the dynamics of local co-/counter-gradient transport and backscatter at the requisite scales of interest. These observations are further confirmed through a posteriori investigation of the flow in a periodic channel at Re_τ = 2000. The results reveal that the TLS method can capture both the large- and the small-scale flow physics in a consistent manner, and at a reduced overall cost when compared to the estimated DNS or wall-resolved LES cost.     

Effect of the number of grooves on flow characteristics around a circular cylinder with triangular grooves

In the flow around a circular cylinder, a sudden decrease in the mean drag coefficient occurs at a high Reynolds number, but the same phenomenon occurs at a lower Reynolds number in the case where there exist grooves or roughness on the cylinder surface. In this paper, in order to make clear the flow characteristics around a cylinder with 20, 26 and 32 triangular grooves, the mean drag coefficient, pressure distribution, velocity distribution and turbulence intensity distribution were measured. Moreover, the flow around the cylinder was analyzed by applying the RNGk − ɛ turbulent model, and the surface flow pattern was investigated using the oil-film technique. From these results, it was found that a sudden decrease in the mean drag coefficient of a cylinder with 32 triangular grooves occurs at a lower Reynolds number compared with 20 and 26 triangular grooves.     

Scaling in large Prandtl number turbulent thermal convection

We study the scaling properties of heat transfer Nu in turbulent thermal convection at large Prandtl number Pr using a quasi-linear theory. We show that two regimes arise, depending on the Reynolds number Re. At low Reynolds number, NuPr ^-1/2 and Re are a function of RaPr ^-3/2. At large Reynolds number NuPr ^1/3 and RePr are function only of RaPr ^2/3 (within logarithmic corrections). In practice, since Nu is always close to Ra ^1/3, this corresponds to a much weaker dependence of the heat transfer in the Prandtl number at low Reynolds number than at large Reynolds number. This difference may solve an existing controversy between measurements in SF₆ (large Re) and in alcohol/water (lower Re). We link these regimes with a possible global bifurcation in the turbulent mean flow. We further show how a scaling theory could be used to describe these two regimes through a single universal function. This function presents a bimodal character for intermediate range of Reynolds number. We explain this bimodality in term of two dissipation regimes, one in which fluctuation dominate, and one in which mean flow dominates. Altogether, our results provide a six parameters fit of the curve Nu(Ra, Pr) which may be used to describe all measurements at Pr≥0.7. Received 27 February 2002 / Received in final form 29 May 2002 Published online 31 July 2002     

Large-eddy simulations of temporally accelerating turbulent channel flow

The unsteady turbulent channel flow subject to the temporal acceleration is considered in this study. Large-eddy simulations were performed to study the response of the turbulent flow to the temporal acceleration. The simulations were started with the fully developed turbulent channel flow at an initial Reynolds number of Re₀ = 3500 (based on the channel half-height and the bulk-mean velocity), and then a constant temporal acceleration was applied. During the acceleration, the Reynolds number of the channel flow increased linearly from the initial Reynolds number to the final Reynolds number of Re₁ = 22,600. The effect of grid resolution, domain size, time step size on the simulation results was assessed in a preliminary study using simulations of the accelerating turbulent flow as well as simulations of the steady turbulent channel flow at various Reynolds numbers. Simulation parameters were carefully chosen from the preliminary study to ascertain the accuracy of the simulation. From the accelerating turbulent flow simulations, the delays in the response of various flow properties to the temporal acceleration were measured. The distinctive features of the delays responsible for turbulence production, energy redistribution, and radial propagation were identified. Detailed turbulence statistics including the wall shear stress response during the acceleration were examined. The results reveal the changes in the near-wall structures during the acceleration. A self-sustaining mechanism of turbulence is proposed to explain the response of the turbulent flow to the temporal acceleration. Although the overall flow characteristics are similar between the channel and pipe flows, some differences were observed between the two flows.     

Effect of groove shape on flow characteristics around a circular cylinder with grooves

In the flow around a circular cylinder, a sudden decrease in the drag force occurs at a high Reynolds number, but the same phenomenon occurs at a lower Reynolds number in the case where there exist grooves or roughness on the circular cylinder surface. In this paper, in order to make clear the flow characteristics around a circular cylinder in the case of changing the shapes of grooves, the drag coefficient, pressure distribution, velocity distribution and turbulent distribution were measured. Moreover the flow around the cylinder was analyzed by applying the RNGk · ∈ turbulent model, and the surface flow pattern was investigated using the oil-film technique. From these results, it is clear that the drag coefficient of a circular cylinder with triangular grooves decreases by about 15% compared with that of a circular cylinder with arc grooves.     

Effect of Microbubble Gas Saturation on Near-Wall Turbulence and Drag Reduction

The article presents results of an experimental study of the effect of gravitational orientation of the flow along its lower/upper solid boundaries on reduction of turbulent drag and void fraction profiles with injection of gas through a porous channel wall. The shear stress on the wall was measured in the Reynolds number range Re_x = (0.23–1.1) × 10⁷ by floating element transducers; the void fraction profile was determined using a fiber-optic sensor. The void fraction in the inner (near-wall) region of the boundary layer was shown to be a key parameter for turbulent drag reduction. The size of the inner region depends on the gas flow rate, the fluid velocity, the distance downstream of the gas generator, and the gravitational orientation of the wall.     

A comparison of vortex and pseudo-spectral methods for the simulation of periodic vortical flows at high Reynolds numbers

We present a validation study for the hybrid particle-mesh vortex method against a pseudo-spectral method for the Taylor–Green vortex at Re_Γ = 1600 as well as in the collision of two antiparallel vortex tubes at Re_Γ = 10,000. In this study we present diagnostics such as energy spectra and enstrophy as computed by both methods as well as point-wise comparisons of the vorticity field. Using a fourth order accurate kernel for interpolation between the particles and the mesh, the results of the hybrid vortex method and of the pseudo-spectral method agree well in both flow cases. For the Taylor–Green vortex, the vorticity contours computed by both methods around the time of the energy dissipation peak overlap. The energy spectrum shows that only the smallest length scales in the flow are not captured by the vortex method.In the second flow case, where we compute the collision of two anti-parallel vortex tubes at Reynolds number 10,000, the vortex method results and the pseudo-spectral method results are in very good agreement up to and including the first reconnection of the tubes. The maximum error in the effective viscosity is about 2.5% for the vortex method and about 1% for the pseudo-spectral method. At later times the flows computed with the different methods show the same qualitative features, but the quantitative agreement on vortical structures is lost.     

Computational study of influences of a seam line of a ball for baseball on flows

Flows around a ball used in baseball games are calculated using third-order upwinddifference method with various seam positions determined by two rotation angles. Those are four-seam rotation with an angle: a and two-seam rotation with an angle: b. The computed results of the four-seam rotation are compared with experimental data measured in a wind tunnel and computed drag coefficients qualitatively agree well with experiments. However, lift coefficients do not agree well. The computed results and geometrical symmetry suggest that a supporting rod in the wind tunnel would have strong influence on the accuracy of the measurement. Flow changes in two-seam rotation are also simulated. It is found that the lowest drag force is observed atb=90 and that the value is less than half of the largest drag force ata=30 and 60 degrees. The largest lift force is observed atb=20 degree. In this case, a projection of the seam line on the top causes a large separation while smooth surface without the seam at the bottom dose not separate the flow. A pair of longitudinal vortices are found in the wake, which make wake slant and generate large lift force.     

Two-level discretizations of nonlinear closure models for proper orthogonal decomposition

Proper orthogonal decomposition has been successfully used in the reduced-order modeling of complex systems. Its original promise of computationally efficient, yet accurate approximation of coherent structures in high Reynolds number turbulent flows, however, still remains to be fulfilled. To balance the low computational cost required by reduced-order modeling and the complexity of the targeted flows, appropriate closure modeling strategies need to be employed. Since modern closure models for turbulent flows are generally nonlinear, their efficient numerical discretization within a proper orthogonal decomposition framework is challenging. This paper proposes a two-level method for an efficient and accurate numerical discretization of general nonlinear closure models for proper orthogonal decomposition reduced-order models. The two-level method computes the nonlinear terms of the reduced-order model on a coarse mesh. Compared with a brute force computational approach in which the nonlinear terms are evaluated on the fine mesh at each time step, the two-level method attains the same level of accuracy while dramatically reducing the computational cost. We numerically illustrate these improvements in the two-level method by using it in three settings: the one-dimensional Burgers equation with a small diffusion parameter ν = 10^?3, the two-dimensional flow past a cylinder at Reynolds number Re = 200, and the three-dimensional flow past a cylinder at Reynolds number Re = 1000.     

On mean flow universality of turbulent wall flows. II. Asymptotic flow analysis

Understanding of the structure of turbulent flows at extreme Reynolds numbers (Re) is relevant because of several reasons: almost all turbulence theories are only valid in the high Re limit, and most turbulent flows of practical relevance are characterized by very high Re. Specific questions about wall-bounded turbulent flows at extreme Re concern the asymptotic laws of the mean velocity and turbulence statistics, their universality, the convergence of statistics towards their asymptotic profiles, and the overall physical flow organization. In extension of recent studies focusing on the mean flow at moderate and relatively high Re, the latter questions are addressed with respect to three canonical wall-bounded flows (channel flow, pipe flow, and the zero-pressure gradient turbulent boundary layer). Main results reported here are the asymptotic logarithmic law for the mean velocity and corresponding scale-separation laws for bulk flow properties, the Reynolds shear stress, the turbulence production and turbulent viscosity. A scaling analysis indicates that the establishment of a self-similar turbulence state is the condition for the development of a strict logarithmic velocity profile. The resulting overall physical flow structure at extreme Re is discussed.     

Study on the optimum mixing throat length for drive nozzle position of the central jet pump

The present paper describes a numerical prediction of optimum mixing throat length for various drive nozzle positions of the central jet pump. The flow pattern and pressure distribution in the pump for various positions of the drive nozzle are investigated by three-dimensional numerical analysis using the RNG k-ε turbulent flow model. Numerical analysis was carried out for values of the nozzle throat ratiod/D of the jet pump of 0.5, 0.6 and 0.7, respectively. The static pressure in the flow field of the jet pump is calculated for the following conditions: (1) drive nozzle position from the entrance of the throatl /D=0 ∼ 2.0, (2) flow rate ratioM=0∼ 1.2, and (3) Reynolds numberRe=3.6×10⁵. These calculations revealed that (1) the optimum length of the mixing throat forl/D=0∼ 1.0 isLm/D=2.0 ∼ 3.5, (2) the length of the mixing throat forl/D=0 andM=0 (suction flow rate ratio=0) is approximatelyLm/D=3.5, and (3) the maximum efficiency is obtained ford/D=0.6 atl/D=0.5. Moreover, the flow pattern in the mixing throat is clarified through a spark tracing experiment. The results obtained in the visualization experiment and the numerically obtained mixing length agreed well.     

Vortex dynamics of a trapezoidal bluff body placed inside a circular pipe

The influence of Reynolds number and blockage ratio on the vortex dynamics of a trapezoidal bluff body placed inside a circular pipe is studied experimentally and numerically. Low aspect ratio, high blockage ratio, curved end conditions (junction of pipe and bluff body), axisymmetric upstream flow with shear and turbulence are some of the intrinsic features of this class of bluff body flows which have been scarcely addressed in the literature. A large range (200:200,000) of Reynolds number (Re_D) is covered in this study, encompassing all the three pipe flow regimes (laminar, transition and turbulent). Four different flow regimes are defined based on the distinct features of Strouhal number (St)–Re_D relation: steady, laminar irregular, transition and turbulent. The wake in the steady regime is stationary with no oscillations in the shear layer. The laminar regime is termed as irregular owing to irregular vortex shedding. The vortex shedding in this regime is observed to be symmetric. The emergence of separation bubble downstream of the bluff body on either side is another interesting feature of this regime, which is further observed to be symmetric. Two pairs of mean streamwise vortices are noticed in the near-wake regime, which are termed as reverse dipole-type wake topology. Beyond the irregular laminar regime, the Strouhal number falls gradually and vortex shedding becomes more periodic. This regime is named transition and occurs close to the Reynolds number at which transition to turbulence takes place in a fully developed pipe. The turbulent regime is characterised by a nearly constant Strouhal number. Typical Karman-type vortex shedding is noticed in this regime. The convection velocity, wake width formation length and irrecoverable pressure loss are quantified to highlight the influence of blockage ratio. These results will be useful to develop basic understanding of vortex dynamics of confined bluff body flow for several practical applications.