Eddy Structures in a Simulated Low Reynolds Number Turbulent Boundary Layer

Coherent structures in a tripped turbulent boundary layer havebeen analysed by applying a new conditional sampling algorithm tolarge-eddy simulation (LES) data. The space-time development of theevents and structural characteristics were examined. The new conditionalsampling scheme is shown to be very effective in the eduction of thecoherent eddy structures, allowing the simultaneous detection, trackingand averaging of several three-dimensional (3-D) structures. Alignmentof the triggering events in all spatial directions and flip-averaging(spatial reflection of the samples) enhance significantly the extractionof detailed features of the structures. The detection method minimisesthe smearing of the spatial details and avoids imposing artificialsymmetry, which is often intrinsic to other conditional samplingschemes.The results show the existence of cigar-shapedstreamwise vortices which are directly associated with negative pressurefluctuation peaks (positive source term of the pressure Poissonequation). They are inclined at 12^° to the wall andtilted laterally at an angle of 7^° to the streamwisedirection. The streamwise vortices induce ejections and sweeps throughan advection mechanism due to the tilting or inclination of thevortices. There is no evidence of hairpin vortices in either theconditional averages or instantaneous flow fields. Near-wall shearlayers are found to be related to the positive pressure fluctuationpeaks as a result of complex interactions between ejection and sweepevents. The shear layer structure has an inclination of10^°, being located between two tilted, streamwisevortices of opposite direction of rotation.The presentresults are very close to other Direct Numerical Simulation studiesusing different conditional sampling schemes. Conceptual models for thestreamwise vortices and shear layer structures are proposed to accountfor the results.     

Reynolds Number Dependence of Velocity Structure Functions in a Turbulent Pipe Flow

Low-order moments of the increments δu andδv where u and v are the axial and radial velocity fluctuations respectively, have been obtained using single and X-hot wires mainly on the axis of a fully developed pipe flow for different values of the Taylor microscale Reynolds numberR _λ. The mean energy dissipation rate〉ε〈 was inferred from the uspectrum after the latter was corrected for the spatial resolution of the hot-wire probes. The corrected Kolmogorov-normalized second-order structure functions show a continuous evolution withR _λ. In particular, the scaling exponentζ _v , corresponding to the v structure function, continues to increase with R _λ in contrast to the nearly unchanged value of ζ _u . The Kolmogorov constant for δu shows a smaller rate of increase with R _λ than that forδv. The level of agreement with local isotropy is examined in the context of the competing influences ofR _λ and the mean shear. There is close but not perfect agreement between the present results on the pipe axis and those on the centreline of a fully developed channel flow. This revised version was published online in July 2006 with corrections to the Cover Date.     

Large-Scale Streamwise Vortices in Turbulent Channel Flow Induced by Active Wall Actuations

Direct numerical simulations of turbulent flow in a plane channel using spanwise alternatively distributed strips (SADS) are performed to investigate the characteristics of large-scale streamwise vortices (LSSVs) induced by small-scale active wall actuations, and their role in suppressing flow separation. SADS control is obtained by alternatively applying out-of-phase control (OPC) and in-phase control (IPC) to the wall-normal velocity component of the lower channel wall, in the spanwise direction. Besides the non-controlled channel flow simulated as a reference, four controlled cases with 1, 2, 3 and 4 pairs of OPC/IPC strips are studied at M =?0.2 and R e =?6,000, based on the bulk velocity and the channel half height. The case with 2 pairs of strips, whose width is Δz ⁺ =?264 based on the friction velocity of the non-controlled case, is the most effective in terms of generating large-scale motions. It is also found that the OPC (resp. IPC) strips suppress (resp. enhance) the coherent structures and that leads to the creation of a vertical shear layer, which is responsible for the LSSVs presence. They are in a statistically steady state and their cores are located between two neighbouring OPC and IPC strips. These motions contribute significantly to the momentum transport in the wall-normal and spanwise directions showing potential for flow separation suppression.     

Statistical Structure at the Wall of the High Reynolds Number Turbulent Boundary Layer

The one and two-point statistical structure of very high Reynolds number turbulence in the surface layer near a rigid `wall' is analysed. The essential mechanisms for turbulent eddies impinging on the wall are studied using linearised rapid distortion theory, which show how the mean shear and blocking actions of the surface act first independently and then, over the life time of the eddy, interactively. Previous analytical results are reinterpreted and some new results are derived to show how the integral length scales, cross correlations and spectra of the different components of the turbulence are distorted depending on the form of the spectra of eddies above the surface layer and how they are related to motions of characteristic eddy structures near the surface. These results are applied to derive some quantitative and qualitative predictions in the surface layers (SL), where the eddies are affected by local shear dynamics, and in the `eddy surface layer' (ESL) where quasi independents loping elongated eddies interact directly with the wall, and where there is a large range of wave number within which the spectra of the horizontal velocity components are proportional to k ⁻¹. The longest eddies in the boundary layer occur near the wall. Field experiments agree with the theoretical model predictions for the quite different forms for the spectra, cospectra and cross correlations for the vertical and horizontal components of the velocity field. By showing that in SL the energy exchange between the large and small scale eddies is local(`staircase') energy cascade, whereas in ESL there is a direct nonlocal (`elevator-like')energy transfer to the small scales, it is shown why the thickness of the ESL increases over rougher surfaces and as the Reynolds number decreases. This revised version was published online in August 2006 with corrections to the Cover Date.     

The Relaminarization Mechanisms of Turbulent Channel Flow at Low Reynolds Numbers

The mechanisms of laminarization in wall-bounded flows have been investigated by performing direct numerical simulations (DNS) of turbulent channel flows. By decreasing Reynolds numbers systematically, the effects of the low Reynolds number are studied in connection with the near-wall turbulent structure and turbulent statistics. At approximately the critical Reynolds number, the turbulent skin friction is reduced, and the turbulent structure changes qualitatively in the very near-wall region. Instantaneous turbulent structures reveal that streamwise vortices, the cores of which are at y+ 10, disappear, although low speed streaks and Reynolds shear stress are still produced by larger streamwise vortices located in the buffer region y+ > 10. Sweep motions induced by these vortical structures are shifted toward the center of a channel and also significantly deterred, which may heighten the effects of the viscous sublayer over most of the channel section and suppress the regeneration mechanisms of new streamwise vortices in the very near-wall region. To investigate the details of how large-scale coherent vortices affect the viscous sublayer and the relevant small-scale streamwise vortices, a body force is virtually imposed in the wall-normal direction to enhance the large streamwise vortices. As a result, it is found that when they are sufficiently enhanced, the small-scale vortices reappear, and the sweep events are again dominant in the viscous sublayer.     

Approximate Wall Boundary Conditions in the Large-Eddy Simulation of High Reynolds Number Flow

The near-wall regions of high Reynolds numbers turbulent flows must be modelled to treat many practical engineering and aeronautical applications. In this review we examine results from simulations of both attached and separated flows on coarse grids in which the near-wall regions are not resolved and are instead represented by approximate wall boundary conditions. The simulations use the dynamic Smagorinsky subgrid-scale model and a second-order finite-difference method. Typical results are found to be mixed, with acceptable results found in many cases in the core of the flow far from the walls, provided there is adequate numerical resolution, but with poorer results generally found near the wall. Deficiencies in this approach are caused in part by both inaccuracies in subgrid-scale modelling and numerical errors in the low-order finite-difference method on coarse near-wall grids, which should be taken into account when constructing models and performing large-eddy simulation on coarse grids. A promising new method for developing wall models from optimal control theory is also discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Hydrodynamic Interactions at Low Reynolds Number

We consider the hydrodynamic interactions of low Reynolds number microswimmers, presenting a review of recent work based upon models of linked sphere swimmers. Particular attention is paid to those aspects that are generic, applicable to all microswimmers and not only to the simple models considered. The importance of the relative phase in swimmer–swimmer interactions is emphasised, as is the role of simple symmetry arguments in both understanding and constraining the hydrodynamic properties of microswimmers.     

On the Lower Critical Reynolds Number for Flow in a Circular Pipe

Flow in a circular pipe is investigated experimentally at Reynolds numbers higher than that at which the resistance coefficients calculated from the Blasius formula for laminar flow and from the Prandtl formula for turbulent flow are equal. The corresponding Reynolds number based on the mean-flow velocity and the pipe diameter is about 1000. The experiments were performed at a high level of inlet pulsations produced by feeding gas into the pipe through a hole with a diameter several times smaller than the pipe diameter. In our experiments the critical Reynolds number was determined as the value, independent of the distance from the inlet, at which the ratio of the axial to the mean-flow velocity as a function of the Reynolds number deviated from 2. At the maximum ratio of the pipe cross-sectional area to the area of the hole through which the gas entered the pipe, equal to 26, the critical Reynolds number was about 2300. After a fivefold increase in the hole area the critical Reynolds number increased by approximately 4%.At Reynolds numbers below 2000, after at a high level of the inlet pulsations an almost laminar flow had developed in the pipe, a perturbation was introduced by inserting a diametrically oriented cylindrical rod with a diameter 10–20 times smaller than the pipe diameter. In these experiments, at Reynolds numbers higher than 1000, at a distance from the rod equal to 50 pipe diameters the axial to mean-flow velocity ratio was less than 2, approaching this value again at large distances from the rod. The insertion of the rod led to a decrease in the critical Reynolds number by approximately 12%.     

Flow Rate Prediction for a Semi-permeable Membrane at Low Reynolds Number in a Circular Pipe

Transport in Porous Media - A concise and accurate prediction method is required for membrane permeability in chemical engineering and biological fields. As a preliminary study on this topic, we...     

Application of Differing Forcing Function Models on Simulated Flow Past an Oscillating Cylinder in a Uniform Low Reynolds Number Flow

A Computational Fluid Dynamics (CFD) model is presented for the uniform viscous two dimensional flow past an oscillating cylinder at low Reynolds number. Numerical simulations are made to study the effect of differing forced induced oscillation mechanisms with a large range of cylinder forcing frequencies. In the first case sinusoidal velocity slip boundary conditions are adopted for the cylinder surface to simulate cylinder oscillation. The implication suggests that no modification or additional term need to be added to the Navier-Stokes equations. In the second case this time extra body force terms which are assumed to account for velocity effects due to cylinder movement are included in the Navier-Stokes equations with the imposition of same boundary conditions. Drag and lift coefficients are extracted from present numerical results and other detailed computations of these coefficients are made at a Reynolds number of 80 and an amplitude-to diameter ratio 0.14. The results are found to be in agreement with each other at low force driving frequencies below and near lock-in. However, differences are found at higher frequencies above lock-in. Agreement are also found with experimental results at some frequency ranges.     

Turbulent Duct Flow Controlled with Spanwise Wall Oscillations

The spanwise oscillation of channel walls is known to substantially reduce the skin-friction drag in turbulent channel flows. In order to understand the limitations of this flow control approach when applied in ducts, direct numerical simulations of controlled turbulent duct flows with an aspect ratio of A R = 3 are performed. In contrast to channel flows, the spanwise extension of the duct is limited. Therefore, the spanwise wall oscillation either directly interacts with the duct side walls or its spatial extent is limited to a certain region of the duct. The present results show that this spanwise limitation of the oscillating region strongly diminishes the drag reduction potential of the control technique. We propose a simple model that allows estimating the achievable drag reduction rates in duct flows as a function of the width of the duct and the spanwise extent of the controlled region.     

雷诺数对三角翼绕流的影响

应用染色液流动显示技术研究了雷诺数对６０°尖前缘三角翼前缘涡破裂位置、背风面流动结构等的影响，并详细分析了背风面流动随攻角的变化．     

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.     

Reynolds Number Effects on Statistics and Structure of an Isothermal Reacting Turbulent Wall-Jet

Three-dimensional direct numerical simulation (DNS) is used to investigate the effects of changing the Reynolds number on dynamics of a reacting turbulent wall-jet. The flow is compressible and a single-step isothermal global reaction is considered. At the inlet, fuel and oxidizer enter the domain separately in a non-premixed manner. In this study, the bulk Reynolds number of the flow, in terms of the inlet quantities, varies from Re = 2000 to Re = 6000, which results in a comparable change in friction Reynolds numbers. The DNS database in Pouransari et al. (Phys. Fluids 23(085104), 2011) is used for the lower Reynolds number case and for the higher Reynolds number case, a new DNS is performed. One of the main objectives of this study is to compare the influences of changing the Reynolds number of the isothermal flow with the heat-release effects caused by the chemical reaction, that we studied earlier in Pouransari et al. (Int. J. Heat Fluid Flows 40, 65–80, 2013). While, both turbulent and flame structures become finer at the higher Reynolds number, the effect of decreasing the Reynolds number and adding the combustion heat release are compared with each other and found to be similar for some aspects of the flow, but are not always the same.     

Dissipation in Turbulent Shear Flows: A Wall Function Approach for High Reynolds Numbers

Dissipation effects that are neglected in almost all numerical codes when it comes to turbulent flows are analysed systematically. It is shown how the direct and the indirect dissipation can be implemented into the codes when asymptotic considerations about the flow behaviour close to a wall are taken into account. As an example, adiabatic pipe flow is analysed with special emphasis on the distribution of energy between the mechanical and thermal part of the total energy.     

Non-Local Interactions and Feedback Instability in a High Reynolds Number Flow

The question of non-locality is considered for a model supersonic flow at high Reynolds number in a channel formed between two parallel plates of different length, using the channel length as a control parameter. Examples are given of time-periodic stable and unstable flows forced by a disturbance positioned in the middle of the channel. It is shown that in certain parameter ranges the flow in a channel of ever increasing length is not approximated by the solutions obtained for infinitely long channels. This is interpreted in terms of a feedback interaction between the flow near the channel ends and the disturbance source. Feedback is shown to result from a slow upstream decay of disturbances coupled with a relatively fast downstream growth of instability waves. For a free (non-forced) flow, the feedback is found to lead to a form of global or resonant instability. Examples of growth rate calculations for the feedback modes are given.     

Turbulent Channel Flow with an Artificial Two-Dimensional Wall Layer

Turbulent flow of an incompressible fluid in a plane channel with parallel walls is considered. The three-dimensional time-dependent Navier-Stokes equations are solved numerically using the spectral finite-difference method. An artificial force which completely suppresses lateral oscillations of the velocity is introduced in the near-wall zone (10 % of the channel half-width in the neighborhood of each wall). Thus, the three-dimensional flow zone, in which turbulent oscillations can develop, is separated from the wall by a fluid layer. It is found that the elimination of three-dimensionality in the neighborhood of the walls leads to a significant reduction in the drag. However, complete laminarization does not occur. The flow in the stream core remains turbulent and can be interpreted as a turbulent flow in a channel with walls located on the boundary of the two-dimensional layer and traveling at the local mean-flow velocity. The oscillations developing inside the two-dimensional layer, which have significant amplitude, distort the flow only in the adjacent zone. Beyond this zone the distributions of the mean characteristics and the structure of instantaneous fields completely correspond to ordinary turbulent flow in a channel with rigid walls. The results obtained confirm the hypothesis of the unimportance of the no-slip boundary conditions for the fluctuating velocity component in the mechanism of onset and self-maintenance of turbulence in wall flows.     

Direct Numerical Simulation of Turbulent Pipe Flow at Moderately High Reynolds Numbers

Fully resolved direct numerical simulations (DNSs) have been performed with a high-order spectral element method to study the flow of an incompressible viscous fluid in a smooth circular pipe of radius R and axial length 25R in the turbulent flow regime at four different friction Reynolds numbers Re _τ ?=?180, 360, 550 and $1\text{,}000$ . The new set of data is put into perspective with other simulation data sets, obtained in pipe, channel and boundary layer geometry. In particular, differences between different pipe DNS are highlighted. It turns out that the pressure is the variable which differs the most between pipes, channels and boundary layers, leading to significantly different mean and pressure fluctuations, potentially linked to a stronger wake region. In the buffer layer, the variation with Reynolds number of the inner peak of axial velocity fluctuation intensity is similar between channel and boundary layer flows, but lower for the pipe, while the inner peak of the pressure fluctuations show negligible differences between pipe and channel flows but is clearly lower than that for the boundary layer, which is the same behaviour as for the fluctuating wall shear stress. Finally, turbulent kinetic energy budgets are almost indistinguishable between the canonical flows close to the wall (up to y ^?+??≈?100), while substantial differences are observed in production and dissipation in the outer layer. A clear Reynolds number dependency is documented for the three flow configurations.