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.     

Direct Numerical Simulation of an Accelerating Channel Flow

Direct Numerical Simulation of a linearly accelerating channel flow starting from an initially statistically steady turbulent flow has been performed. It is shown that the response of the accelerating flow is fundamentally the same as that of the step-change transient flow described in He and Seddighi (J Fluid Mech 715:60–102, 2013). The flow structure again behaves like a boundary layer bypass transition undergoing three distinct phases, namely, (i) initially (pre-transition), the flow is laminar-like and the pre-existing turbulent structures are modulated resulting in elongated streaks leading to a strong and continuous increase in the streamwise fluctuating velocity but little changes in the other two components; (ii) it then undergoes transition when isolated turbulent spots are generated which spread and merge with each other, and (iii) they eventually cover the entire surface of the wall when the flow is fully turbulent. The similarity between the turbulence responses in the two flows is significant noting the contrasting features of the two types of mean flow unsteadiness: in the step-change flow, a sharp boundary layer is resulted in nearly instantly on the wall which closely resembles the spatially developing boundary layer, whereas the linear flow acceleration causes a continuing change of velocity gradient adjacent to the wall which propagates into the flow field with time, resulting in a gradually-developing boundary layer. There are, however, quantitative differences in the detailed behavior of the two flows and especially the transition is much delayed in the accelerating flow. It is also shown that the late pre-transition and early transition stages in both flows are characterised by significantly increased inwards sweep events in the wall region and ejection events in the outer layer. The flatness of the wall-normal velocity increases markedly near the wall around the time of onset of transition as a consequence of the huge intermittency of the velocity fluctuations. That is, there are long periods of quiescent flow coupled with occasional turbulent bursts.     

The turbulent boundary layer of dissociated gas in the initial section of a tube

Many theoretical and experimental papers [1–4] have been devoted to investigating the turbulent boundary layer in the initial section of a channel. For the most part, however, the flow of an incompressible fluid with constant parameters is considered. There are many practical cases in which it is of interest to treat the development of the turbulent boundary layer of gas in the initial section of a pipe when conditions are strongly nonisothermal. A solution of a problem of this type, based on the theory of limit laws, is given in paper [1]. The present article extends this solution to the case of the flow of a high-enthalpy gas when the effect of gas dissociation on the turbulent boundary layer characteristics must be taken into account. We shall consider the flow of a mixture of i gases which is in a frozen state inside the boundary layer, and in an equilibrium state on its boundaries. Formulas are derived for the laws of friction and heat exchange, and a solution is given for the turbulent boundary layer equations in the initial section of the pipe when the wall temperature is constant and the gas flows at a subsonic velocity.Finally the authors are grateful to S. S. Kutateladze for discussing the paper.     

Nonlinear dynamics of large-scale vortex structures in a turbulent Ekman layer

Nonlinear regimes of roll structure development in the Ekman atmospheric boundary layer are considered. Within the framework of a turbulent boundary layer theory that takes the inherent helicity into account, the helicity effect on the stability and development of secondary vortices is studied. The stabilizing role of helicity in the Ekman layer dynamics is detected. With increasing nonlinearity the rolls become time-dependent self-oscillatory structures exchanging energy and helicity with the main flow. Multi-scale structure regimes are obtained, similar to those observable in the atmospheric boundary layer and the solitary vortex structures previously discovered in bench and numerical experiments for Ekman flows between disks rotating in opposite directions.     

Turbulent particle dispersion in arbitrary wall-bounded geometries: A coupled CFD-Langevin-equation based approach

A Lagrangian continuous random walk (CRW) model is developed to predict turbulent particle dispersion in arbitrary wall-bounded flows with prevailing anisotropic, inhomogeneous turbulence. The particle tracking model uses 3D mean flow data obtained from the Fluent CFD code, as well as Eulerian statistics of instantaneous quantities computed from DNS databases. The turbulent fluid velocities at the current time step are related to those of the previous time step through a Markov chain based on the normalized Langevin equation which takes into account turbulence inhomogeneities. The model includes a drift velocity correction that considerably reduces unphysical features common in random walk models. It is shown that the model satisfies the well-mixed criterion such that tracer particles retain approximately uniform concentrations when introduced uniformly in the domain, while their deposition velocity is vanishingly small, as it should be. To handle arbitrary geometries, it is assumed that the velocity rms values in the boundary layer can locally be approximated by the DNS data of fully developed channel flows. Benchmarks of the model are performed against particle deposition data in turbulent pipe flows, 90° bends, as well as more complex 3D flows inside a mouth-throat geometry. Good agreement with the data is obtained across the range of particle inertia.     

A ghost-cell immersed boundary method for large-eddy simulations of compressible turbulent flows

A methodology to perform a ghost-cell-based immersed boundary method (GCIBM) is presented for simulating compressible turbulent flows around complex geometries. In this method, the boundary condition on the immersed boundary is enforced through the use of ‘ghost cells’ that are located inside the solid body. The computations of variables on these ghost cells are achieved using linear interpolation schemes. The validity and applicability of the proposed method is verified using a three-dimensional (3D) flow over a circular cylinder, and a large-eddy simulation of fully developed 3D turbulent flow in a channel with a wavy surface. The results agree well with the previous numerical and experimental results, given that the grid resolution is reasonably fine. To demonstrate the capability of the method for higher Mach numbers, supersonic turbulent flow over a circular cylinder is presented. While more work still needs to be done to demonstrate higher robustness and accuracy, the present work provides interesting insights using the GCIBM for the compressible flows.     

Low Reynolds number axisymmetric turbulent boundary layer on a cylinder

In the present study, an axisymmetric turbulent boundary layer growing on a cylinder is investigated experimentally using hot wire anemometry. The combined effects of transverse curvature as well as low Reynolds number on the mean and turbulent flow quantities are studied. The measurements include the mean velocity, turbulence intensity, skewness and flatness factors in addition to wall shear stress. The results are presented separately for the near wall region and the outer region using dimensionless parameters suitable for each case. They are also compared with the results available in the open literature.The present investigation revealed that the mean velocity in near wall region is similar to other simple turbulent flows (flat plate boundary layer, pipe and channel flows); but it differs in the logarithmic and outer regions. Further, for dimensionless moments of higher orders, such as skewness and flatness factors, the main effects of the low Reynolds number and the transverse curvature are present in the near wall region as well as the outer region.     

Effects of free‐stream turbulence on large‐scale coherent structures of separated boundary layer transition

It has been well established that large‐scale structures, usually called coherent structures, exist in many transitional and turbulent flows. The topology and range of scales of those large‐scale structures vary from flow to flow such as counter‐rotating vortices in wake flows, streaks and hairpin vortices in turbulent boundary layer. There has been relatively little study of large‐scale structures in separated and reattached transitional flows. Large‐eddy simulation (LES) is employed in the current study to investigate a separated boundary layer transition under 2% free‐stream turbulence on a flat plate with a blunt leading edge. The Reynolds number based on the inlet free stream velocity and the plate thickness is 6500. A dynamic subgrid‐scale model is employed to compute the subgrid‐scale stresses more accurately in the current transitional flow case. Flow visualization has shown that the Kelvin–Helmholtz rolls, which have been so clearly visible under no free‐stream turbulence (NFST) are not as apparent in the present study. The Lambda‐shaped vortical structures which can be clearly seen in the NFST case can hardly be identified in the free‐stream turbulence (FST) case. Generally speaking, the effects of free‐stream turbulence have led to an early breakdown of the boundary layer, and hence increased the randomization in the vortical structures, degraded the spanwise coherence of those large‐scale structures. Copyright © 2005 John Wiley & Sons, Ltd.     

Combined laser-doppler and cold wire anemometry for turbulent heat flux measurement

 combined laser-doppler and cold wire anemometry technique for determining turbulent heat flux is described. The system can be used in flows of arbitrarily high turbulent intensity and large temperature variations. Its potential is demonstrated via measurements in a simulated stable atmospheric boundary layer, for which the Monin-Obukhov length scale was about 70% of the boundary layer depth. Mean and turbulence properties were obtained throughout the boundary layer and the results are shown to be both internally consistent and similar to corresponding field data. Measurements in the highly turbulent, separated flow behind a bluff body mounted in the stable boundary layer are also presented. Received: 9 May 1997 / Accepted: 2 September 1997     

Wave regimes of ascending flow of a thin layer of a viscous liquid in contact with a gas

The flow of a liquid in thin layers is one of the hydrodynamic problems of chemistry and heat engineering. The large surface area of films and their small thickness make it possible to accelerate thermal, diffusive, and chemical processes at the gas-liquid boundary.Theoretical studies of liquid flow in a vertical descending thin layer are presented in [1–4]. In this paper we study ascending wave flows of a liquid in a thin vertical layer in contact with a gas, i.e., flows in the direction opposite the action of the force due to gravity, with account for the action of the gas on the liquid surface. Such motions are encountered when oil is extracted from strata that are saturated with gas. At some distance from the stratum the oil and gas separate: the gas travels at high velocity inside the pipe, occupying a considerable portion of the pipe, and the liquid is displaced toward the pipe walls, forming a thin film. In certain cases a wave-like interface develops between the oil and gas that travels with a velocity greater than that of the liquid but less than the average gas velocity. Similar phenomena are observed in high velocity mass exchangers.We examine the effect of the gas for both laminar and turbulent flow.Studies that neglect the effect of the gas flow on the liquid show that for waves on the film surface whose lengths are considerably longer than the average thickness of the layer, the liquid motion in the film is described by boundary layer equations in which account is taken of the mass force, i.e., the force due to gravity. With some approximation, we can assume that in accounting for the effect of the gas on the liquid the liquid flow is described by these same equations.     

A vorticity stretching diagnostic for turbulent and transitional flows

Vorticity stretching in wall-bounded turbulent and transitional flows has been investigated by means of a new diagnostic measure, denoted by Γ, designed to pick up regions with large amounts of vorticity stretching. It is based on the maximum vorticity stretching component in every spatial point, thus yielding a three-dimensional scalar field. The measure was applied in four different flows with increasing complexity: (a) the near-wall cycle in an asymptotic suction boundary layer (ASBL), (b) K-type transition in a plane channel flow, (c) fully turbulent channel flow at Re _τ = 180 and (d) a complex turbulent three-dimensional separated flow. Instantaneous data show that the coherent structures associated with intense vorticity stretching in all four cases have the shape of flat ‘pancake’ structures in the vicinity of high-speed streaks, here denoted ‘h-type’ events. The other event found is of ‘l-type’, present on top of an unstable low-speed streak. These events (l-type) are further thought to be associated with the exponential growth of streamwise vorticity in the turbulent near-wall cycle. It was found that the largest occurrence of vorticity stretching in the fully turbulent wall-bounded flows is present at a wall-normal distance of y ⁺?=?6.5, i.e. in the transition between the viscous sublayer and buffer layer. The associated structures have a streamwise length of ~200–300 wall units. In K-type transition, the Γ-measure accurately locates the regions of interest, in particular the formation of high-speed streaks near the wall (h-type) and the appearance of the hairpin vortex (l-type). In the turbulent separated flow, the structures containing large amounts of vorticity stretching increase in size and magnitude in the shear layer upstream of the separation bubble but vanish in the backflow region itself. Overall, the measure proved to be useful in showing growing instabilities before they develop into structures, highlighting the mechanisms creating high shear region on a wall and showing turbulence creation associated with instantaneous separations.     

Higher-order and length-scale statistics of velocity and temperature fluctuations in turbulent boundary layer along a heated vertical flat plate

Time-developing direct numerical simulation (DNS) was performed to clarify the higher-order turbulent behaviors in the thermally-driven boundary layers both in air and water along a heated vertical flat plate. The predicted statistics of the heat transfer rates and the higher-order turbulent behaviors such as skewness factors, flatness factors and spatial correlation coefficients of the velocity and temperature fluctuations in the natural-convection boundary layer correspond well with those obtained from experiments for space-developing flows. The numerical results reveal that the turbulent structures of the buoyancy-driven boundary layers are mainly controlled by the fluid motions in the outer region of the boundary layer, and these large-scale structures are strongly connected with the generation of turbulence in the thermally-driven boundary layers, in accordance with the actual observations for space-developing flows. Moreover, to specify the turbulence structures of the boundary layers, the cross-correlation coefficients and the characteristic length scales are examined for the velocity and thermal fields. Consequently, it is found that with a slight increase in freestream velocity, the cross-correlation coefficient for the Reynolds shear stress and turbulent heat flux increases for opposing flow and decreases for aiding flow, and the integral scales for the velocity and temperature fields become larger for opposing flow and smaller for aiding flow compared with those for the pure natural-convection boundary layer.     

Quantitative measurement of spatio-temporal heat transfer to a turbulent water pipe flow

A technique using high-speed infrared thermography was applied to measure the spatio-temporal heat transfer to a turbulent water flow in a horizontal circular pipe. The instantaneous distribution of the heat transfer coefficient and its temporal fluctuation was evaluated by solving inverse heat conduction equation of the heated thin-test-surface. As a result, it was demonstrated that the quantitative measurement, not only the time-averaged heat transfer, but also the statistics of the spatio-temporal fluctuation, was possible using this technique. In addition, a unique feature of the spatio-temporal heat transfer was clearly visualized for the turbulent pipe flow, which was dominated by the streaky structure similar to that for the turbulent boundary layer and the turbulent channel flow.     

明渠湍流边界层中颗粒的运动与分布

利用直接数值模拟、点球浸入边界法和颗粒离散元法相结合的方法, 模拟了颗粒在明渠湍流边界层中的运动, 并对颗粒的瞬时位置进行了Voronoi 分析, 定量研究了颗粒在湍流边界层中的运动和分布规律. 研究发现:颗粒的输运对湍流的统计特征有影响, 其运动与近壁区湍流拟序结构密切相关, 在"喷发"结构作用下被带离壁面, 在"扫掠" 结构和自身重力作用下回到壁面; 在湍流边界层中, 颗粒倾向于聚集在低流速带, 呈条带状分布;颗粒在大部分时间处于"簇"状态, 偶尔跳跃到"空" 状态, 但能够很快回到邻近低速区域.      

Two- and Four-Way Coupled Euler–Lagrangian Large-Eddy Simulation of Turbulent Particle-Laden Channel Flow

Large-eddy simulations (LES) of a vertical turbulent channel flow laden with a very large number of solid particles are performed. The motivation for this research is to get insight into fundamental aspects of co-current turbulent gas-particle flows, as encountered in riser reactors. The particle volume fraction equals about 1.3%, which is relatively high in the context of modern LES of two-phase flows. The channel flow simulations are based on large-eddy approximations of the compressible Navier–Stokes equations in a porous medium. The Euler–Lagrangian method is adopted, which means that for each individual particle an equation of motion is solved. The method incorporates four-way coupling, i.e., both the particle-fluid and particle–particle interactions are taken into account. The results are compared to single-phase channel flow in order to investigate the effect of the particles on turbulent statistics. The present results show that due to particle–fluid interactions the mean fluid profile is flattened and the boundary layer is thinner. Compared to single-phase turbulent flow, the streamwise turbulence intensity of the gas phase is increased, while the normal and spanwise turbulence intensities are reduced. This finding is generally consistent with existing experimental data. The four-way coupled simulations are also compared with two-way coupled simulations, in which the inelastic collisions between particles are neglected. The latter comparison clearly demonstrates that the collisions have a large influence on the main statistics of both phases. In addition, the four-way coupled simulations contain stronger coherent particle structures. It is thus essential to include the particle–particle interactions in numerical simulations of two-phase flow with volume fractions around one percent.     

Calculation of the initial segment and the stabilized flow segment in two-dimensional separation-free diffusers

The aerodynamic characteristics of a series of plane diffusere with straight walls are calculated for a broad range of divergence angles, Reynolds numbers, and the parameter which characterizes the initial flow nonuniformity.The fluid is assumed to be incompressible. The calculations are made for the case in which the boundary layer is fully turbulent, i. e., there is no laminar flow segment near the entrance section.The calculation of separation-free flow in diffuser channels is based on the use of boundary layer theory [1]. It has now become possible to carry out large-scale calculations for diffusers whose geometric and aerodynamic parameters vary over rather wide limits. This is the result both of the use of computers and of the fact that the modern approximate methods for calculating the turbulent boundary layer have been reduced to comparatively simple interpolation formulas [2].Usually, in the calculation of diffusers we examine only the initial flow segment, within the limits of which the boundary layers which form on the walls do not come together, i. e., there is a potential core. The laws governing diffuser flow in the absence of the potential core have been studied relatively little; the only known solutions are those of [3], which are valid at a very great distance from the entrance section.In this study we examine three characteristic flow zones: the initial segment, extending from the entrance section to the section at which the boundary layers come together; the stabilized-flow zone with closed boundary layers comprised of two characteristic segments-the transitional segment extending from the plane where the boundary layers join to the beginning of the radial flow segment; and, finally, the radial (self-similar) flow segment, characterized by constancy of all the dimensionless boundary layer characteristics along the flow. It is obvious that this division into characteristic zones is arbitrary: a consequence of the adopted flow idealization is a break in the curves expressing the aerodynamic characteristics as a function of the axial coordinate at the junction of the initial segment and the stabilized-flow segment. It is well known that a similar phenomenon occurs in the calculation of free turbulent jets based on arbitrary division of the jet into two segments-initial and primary segments [4].The computer calculation of the initial segments was performed by A. N. Smol'yaninova, and the stabilized-flow segments were calculated by I. N. Podol'nyi.     

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.     

Energy redistribution dynamics in coupled Couette–Poiseuille flows using large-Eddy simulation

The problem of turbulent Couette flow driven by a statistically steady external wind is studied in the framework of spatially filtered Navier–Stokes equations. The phenomenon of wind-driven flow of water is represented by a layer of air modeled as Poiseuille flow (air sub-domain), coupled to a layer of water modeled as Couette flow (water sub-domain). We focus on changes in the statistics in either the air or the water sub-domain, due to the coupling with the other sub-domain. We also highlight dynamic flow structures forming near the air-water interface. Simulations based on different Reynolds numbers in the air and the water sub-domains are compared to computationally less demanding simulations with equal Reynolds numbers. Results of these simulations indicate strong similarities, i.e., the flow is well approximated by simulating air and water at the same Reynolds numbers. Further analysis shows that the flow in the water domain shares important features with classical Couette flows. The horizontal turbulent mixing renders a thinner boundary layer in the water sub-domain. Moreover, an increased intermittency in the flow velocities is observed, which may be linked to so-called splat events near the air-water interface. These splats characterize the interaction of coherent structures across the interface, being stronger in the water phase. An analysis of the pressure-strain correlation near the air-water interface on the water side shows that such splats are responsible for redistributing energy from the streamwise and spanwise directions, to the vertical direction. This behavior, although qualitatively similar to wall-bounded flows, differ mainly on the fact that most of the energy drained comes from the streamwise direction: in wall-bounded the main contributor is the spanwise direction. The boundary layers near the air-water interface show inclined vortical structures. Unlike in coupled Couette–Couette flow, the peak in the Reynolds stress is displaced from the channel’s center into the buffer region of the water sub-domain.     

Skin Friction following BLADE manipulation in a turbulent pipe flow

This study attempts to analyze the measured wall shear stress distribution downstream of single and tandem BLADEs in fully developed pipe flow. Previous works have indicated the adverse effect of overall drag increase with the single BLADE in both channel and pipe flows, and an even larger drag increase with the tandem BLADES. This is contrary to that observed for external boundary layer flow. Extensive comparisons are then made to the wall shear stress distribution following BLADEs in boundary layer flow, leading to the conclusion of little or no potential in the application of BLADEs alone to pipe flow.     

A time-resolved hot-wire shear stress probe for turbulent flow: use of laminar flow calibration

A specially-designed rotating rig for producing near Couette flow was used in the calibration of a marginally elevated hot-wire shear stress probe. The probe was then used for measurements in both the turbulent boundary layer and pipe flows. Results showed that the mean wall shear stress can be accurately predicted and the near wall statistical quantities of intensity, skewness and flatness of shear stress fluctuations concurred well with previous works, thereby supporting the notion of a time-resolved shear stress probe for turbulent flows.