Some new results of the observation on the coherent structures of turbulent boundary layer in flow with adverse pressure gradient

By a suitable manipulation of hydrogen bubble generation, some new results were obtained: (1) The long-streaks are generated along the interfaces between low and high-speed streaks. The long-streaks are generally stretching and are moving faster than its neighboring high-speed streaks. The hydrogen bubbles in long-streaks have longer life. (2) The stream-wise vortices are also generated along the interfaces. The project is supported by the National Natural Science Foundation of China.     

Optimal transient growth in turbulent pipe flow

The optimal transient growth process of perturbations driven by the pressure gradient is studied in a turbulent pipe flow. A new computational method is proposed, based on the projection operators which project the governing equations onto the subspace spanned by the radial vorticity and radial velocity. The method is validated by comparing with the previous studies. Two peaks of the maximum transient growth amplification curve are found at different Reynolds numbers ranging from 20 000 to 250 000. The optimal flow structures are obtained and compared with the experiments and DNS results. The location of the outer peak is at the azimuthal wave number n=1, while the location of the inner peak is varying with the Reynolds number. It is observed that the velocity streaks in the buffer layer with a spacing of 100δ_v are the most amplified flow structures. Finally, we consider the optimal transient growth time and its dependence on the azimuthal wave length. It shows a self-similar behavior for perturbations of different scales in the optimal transient growth process.     

逆压梯度边界层未分离情况下的湍流模式适用性研究Turbulence model applicability for boundary layer flow with adverse pressure gradient in attached case

边界层逆压梯度作用下的流动是许多工程中的一个基础问题,由于逆压梯度作用,流动形态复杂,使得数值模拟有很大的难度。基于雷诺平均纳维‐斯托克斯RANS（Reynolds Averaged Navier‐Stokes）方程对二维平板逆压梯度边界层作数值计算研究,选取6种代表性的湍流模式,得到局部摩擦系数的数值解,与实验值比较,发现k‐ω模式具有很好的精度。基于该湍流模式,给出了湍动能分布,该结果有助于认识逆压梯度边界层流动的复杂特征。     

The turbulent/non-turbulent interface in an adverse pressure gradient turbulent boundary layer

The turbulent/non-turbulent interface (TNTI) in an adverse pressure gradient (APG, β = 1.45) turbulent boundary layer (TBL) is explored here by using direct numerical simulation (DNS) data; β is the Clauser pressure gradient parameter. For comparison, the DNS data for a zero pressure gradient (ZPG) TBL is included. The interface is extracted with an approach based on enstrophy criteria. Depending on the enstrophy, the outer boundary layer flow can be classified into the free stream, boundary layer wake, and intermittent flow regimes. The fractal dimension of the interface is obtained by using the box-counting algorithm, and was found to be constant over a long range of box sizes. The TNTI shows a monofractal behavior. The geometric complexity of a TNTI can be determined in terms of the genus, which is defined as the number of handles in a geometric object. We examine the volume and projection area of the genus of the TNTI to analyze the entrainment process. The geometric complexity of the APG TBL interface and the local entrainment are greater than those of the ZPG TBL, as is evident in the increases in the genus near the interface. The local entrainment velocity is dominantly affected by the viscous diffusion at the interface.     

Turbulent flow in a conical diffuser: A review

This is a review of experimental studies of turbulent flow in a conical diffuser by eight Ph.D. students, eleven M.Sc. students, one M.Eng. student, and myself in the past 29 years. During this time, two conical diffusers were constructed: the first was of cast aluminum construction, and the second was of plastic fabrication. These two diffusers were basically the same in geometry except that the pipe section was constructed as an integral part of the plastic diffuser to avoid the lip at the junction of the inlet pipe and the diffuser. The conical diffuser had a total divergence angle of 8°, an area ratio of 4:1, and an inlet diameter of 0.1016 m (4 in.).

The flow at the inlet of the diffuser was usually fully developed pipe flow, but sometimes it was boundary layer grown on the pipe wall. Hot-wire and pulse-wire anemometry together with computer facilities were used to obtain the results of complex flow present in the conical diffuser. Mean velocity profiles were obtained throughout the diffuser, which in turn were used to obtain strain rates and their principal direction. Turbulence moments up to fourth order were measured. The results were used to assess momentum, turbulent kinetic energy, and shear stress equations. Other features such as instantaneous flow reversals in the wall region, relative strength of large eddies, extra strain rate, and the production of kinetic energy also were investigated to find the dynamical picture in the diffuser flow.     

A robust numerical method for flow through a pipe driven by an oscillating pressure gradient

The problem of periodic flow of an incompressible fluid through a pipe, which is driven by an oscillating pressure gradient (e.g. a reciprocating piston), is investigated in the case of a large Reynolds number. This process is described by a singularly perturbed parabolic equation with a periodic right‐hand side, where the singular perturbation parameter is the viscosity ν. The periodic solution of this problem is a solution of the Navier–Stokes equations with cylindrical symmetry. We are interested in constructing a parameter‐robust numerical method for this problem, i.e. a numerical method generating numerical approximations that converge uniformly with respect to the parameter ν and require a bounded time, independent of the value of ν, for their computation. Our method comprises a standard monotone discretization of the problem on non‐standard piecewise uniform meshes condensing in a neighbourhood of the boundary layer. The transition point between segments of the mesh with different step sizes is chosen in accordance with the behaviour of the analytic solution in the boundary layer region. In this paper we construct the numerical method and discuss the results of extensive numerical experiments, which show experimentally that the method is parameter‐robust. Copyright © 2006 John Wiley & Sons, Ltd.     

具有壁面射流激励的圆管内气相流动大涡模拟

采用大涡模拟方法和单方程亚格子模式对小尺度量进行模拟。研究了不同强度壁面射流激励对圆管内气相流动的影响,模拟结果给出了射流对瞬态拟序结构发展、时平均流向速度分布的影响。随着射流强度的增加,射流入口附近流体的回流现象增强。射流强度足够大时可以减小管壁处的切应力值,同时会减小壁面附近流动速度,这种速度分布会导致气体夹带颗粒的能力下降,从而在实际两相流动中容易造成壁面附近的气粒返混现象。     

On the flows produced by sudden application of a constant pressure gradient or by impulsive motion of a boundary

Some properties of the time-dependent Navier-Stokes equations are discussed for flows impulsively started from rest by sudden application of a constant pressure gradient or by the impulsive motion of a boundary. Five illustrative examples are given. They are: unsteady flow in a circular cylinder moving parallel to its length, starting flow in a circular pipe, unsteady flow in a rotating cylinder, starting flow in a rectangular channel moving parallel to its length and unsteady flow in a channel of rectangular cross-section. It is found that the expressions of the quantities such as velocity, flux and skin friction are in series forms which may be rapidly convergent for large values of the time but slowly convergent for small values of the time or vice versa. It is shown that if their expressions can be found for one of large values of the time or small values of the time, these expressions can be used for the other.     

Heat transfer and skin-friction in a turbulent boundary layer under a non-equilibrium longitudinal adverse pressure gradient

The experimental data on the effect of weak and moderate non-equilibrium adverse pressure gradients (APG) on the parameters of dynamic and thermal boundary layers are presented. The Reynolds number based on the momentum thickness at the beginning of the APG region was Re^** = 5500. The APG region was a slot channel with upper wall expansion angles from 0 to 14°. The profiles of the mean and fluctuation velocity components were measured using a single-component hot-wire anemometer. The friction coefficients were determined using two methods, namely, the indirect Clauser method and the direct method of weighting the lower wall region on a single-component strain-gage balance. The heat transfer coefficients were determined by a transient method using an IR camera. It is noticed that in the pressure gradient range realized the universal logarithmic region in the boundary layer profile is conserved. The values of the relative (divided by the parameters in zero gradient flow at the same value of Re^**) friction and heat transfer coefficients, together with the Reynolds analogy factor, are determined as functions of the longitudinal pressure gradient. The values of the relative friction coefficient reduced to c_f/c_f0 = 0.7 and those of the heat transfer to St/St₀ = 0.9. A maximum value of the Reynolds analogy factor (St/St₀)/(c_f/c_f0) = 1.16 was reached for the pressure gradient parameter β = 2.9.     

Numerical study of an oscillatory turbulent flow over a flat plate

Oscillatory turbulent flow over a flat plate is studied using large eddy simulation (LES) and Reynolds-average Navier-Stokes (RANS) methods. A dynamic subgrid-scale model is employed in LES and Saffman's turbulence model is used in RANS. The flow behaviors are discussed for the accelerating and decelerating phases during the oscillating cycle. The friction force on the wall and its phase shift from laminar to turbulent regime are also investigated for different Reynolds numbers. The project supported by the Youngster Funding of Academia Sinica and by the National Natural Science Foundation of China     

TRPIV measurement of turbulent coherent structures in a drag-reducing flow bypolymers

Fully developed turbulentflow fields with and without polymer solution at the same Reynolds number were measured by time-resolved particle image velocimetry (TRPIV) in a water channel toinvestigate the mechanism of drag-reducingsolution from theview of coherent structures manipulation. The streamwise mean velocity and Reynolds stress profiles in thesolution werecompared with those in water. After adding the polymer solution, the Reynolds stress in the near-wall area decreases significantly. Theresult relates tightly to the decease of the coherent structures' bursting. The spatial topology of coherentstructures duringbursts has been extracted by the new mu-level criterion based on locally averaged velocitystructure function.The effect of polymers onturbulent coherentstructures mainly reflects in the intensity, not in the shape. In the solution, it is by suppressing thecoherent structuresthat the wall friction isreduced.     

Effect of an adverse pressure gradient on the streamwise Reynolds stress profile maxima in a turbulent boundary layer

It is widely accepted that in a turbulent boundary layer （TBL） with adverse pressure gradient （APG） an outer peak usually appears in the profile of streamwise Reynolds stress. However, the effect of APG on this outer peak is not clearly understood. In this paper, the effect of APG is analysed using the numerical and experimental results in the literature. Because the effect of upstream flow is inherent in the TBL, we first analyse this effect in TBLs with zero pressure gradient on flat plates. Under the individual effect of upstream flow, an outer peak already appears in the profile of streamwise Reynolds stress when the TBL continues developing in the streamwise direction. The APG accelerates the appearance of the outer peak, instead of being a trigger.     

Determination of turbulent pressure loss of non-newtonian oil flow in rough pipes

Summary A detailed theoretical analysis is carried out of the turbulent flow of non-newtonian fluids in rough pipes in which approximate assumptions are of negligible importance. A new friction-factor equation (BNS-equation) is obtained which is valid for turbulent flow of power-law type non-newtonian fluids in the transition region between smooth and wholly rough wall turbulence.The accuracy of the deduced formula was checked by in situ measurements on a crude oil pipeline having a diameter of 305 mm and a length of 161 km. For this pipeline the Reynolds numbers ranged between 10⁴ and 10⁵. The computed and the measured friction factors showed good agreement. The mean of the absolute values of the relative error was 4.1%, while the standard deviation was 0.81%, these results appear to give better fitting than other similar equations published earlier.The BNS-equation is the first analytically deduced relationship for non-newtonian crude oils which includes the effects of both Reynolds number and relative roughness of the pipe.With 6 figures and 2 tables     

Coherent contribution to turbulent mixing of a jet in cross flow

In this work, an experimental study of a jet issuing in cross flow is reported. Support is offered to a theory on the origin of Karman-like vortices, in the wake of jets issuing in boundary layers. Further, a coherent structure with strong contribution to the energetics of the flow field is identified here. This structure has not been reported earlier, according to the present authors' knowledge, and it is shown to be related to engulfment of external fluid at the bottom of the jet, and to interesting stochastic and spectral characteristics of the flow field. Three-dimensional plots of the coherent quantities, based on experimental data, reveal a double-helical morphology of the coherent structure. The same morphology has recently been proposed for the far field of jets issuing in stagnant fluid (i.e., without cross flow). The results of this study are expected to support theoretical and numerical work on jets issuing in cross flow.     

Heat transfer and pressure drop characteristics of turbulent flow in a circular tube fitted with regularly spaced twisted-tape elements

The results of an experimentalnvestigation of turbulent flow heat transfer and pressure drop characteristics in a circular tube fitted with regularly spaced twisted-tape elements connected by thin circular rods are reported. The characteristics are governed by Reynolds number, Prandtl number, twist ratio, space ratio, and rod-to-tube diameter ratio. Correlations for friction factor and Nusselt number are also reported. It is shown that on the basis of both constant pumping power and constant heat duty, regularly spaced twisted-tape elements do not perform better than full-length twisted tapes.     

Results of direct numerical simulation of turbulent flow in a pipe of elliptical cross-section

Within the framework of the complete Navier-Stokes equations the turbulent flow in a pipe of elliptical cross-section with semiaxis ratio b/a = 0.5 is directly calculated for the Reynolds number Re = 6000 determined from the mean-flow velocity and the hydraulic diameter. The distribution of the average and pulsatory flow characteristics over the pipe cross-section are obtained. In particular, the secondary flow in the cross-section plane, typical of turbulent flows in noncircular pipes, is calculated. The equation for the longitudinal vorticity which determines the shape and intensity of the secondary flow is analyzed. In the balance equation for the pulsation kinetic energy the behavior of all the terms that characterize energy production, dissipation and redistribution over the pipe cross-section is described.