Direct numerical simulation of turbulent pipe flow

Fully developed turbulent pipe flow at low Re-number is studied by means of direct numerical simulation (DNS). In contrast to many previous DNS's of turbulent flows in rectangular geometries, the present DNS code, developed for a cylindrical geometry, is based on the finite volume technique rather than being based on a spectral method. The statistical results are compared with experimental data obtained with two different experimental techniques. The agreement between numerical and experimental results is found to be good which indicates that the present DNS code is suited for this kind of numerical simulations.     

Theoretical research on rotational dispersion coefficient of fiber in turbulent shear flow of fiber suspension

The rotational dispersion coefficient of the fiber in the turbulent shear flow of fiber suspension was studied theoretically. The function of correlation moment between the different fluctuating velocity gradients of the flow was built firstly. Then the expression, dependent on the characteristic length, time, velocity and a dimensionless parameter related to the effect of wall, of rotational dispersion coefficient is derived. The derived expression of rotational dispersion coefficient can be employed to the inhomogeneous and non-isotropic turbulent flows. Furthermore it can be expanded to three-dimensional turbulent flows and serves the theoretical basis for solving the turbulent flow of fiber suspension.     

research on rotational dispersion coefficient of fiber in turbulent shear flow of fiber suspension

The rotational dispersion coefficient of the fiber in the turbulent shear flow of fiber suspension was studied theoretically. The function of correlation moment between the different fluctuating velocity gradients of the flow was built firstly. Then the expres- sion, dependent on the characteristic length, time, velocity and a dimensionless parameter related to the effect of wall, of rotational dispersion coefficient is derived. The derived expression of rotational dispersion coefficient can be employed to the inhomogeneous and non-isotropic turbulent flows. Furthermore it can be expanded to three-dimensional turbulent flows and serves the theoretical basis for solving the turbulent flow of fiber suspension.     

Dynamics of short fiber suspensions in bounded shear flow

We have studied the dynamics of non-colloidal short fiber suspensions in bounded shear flow using the Stokesian dynamics simulation. Such particles make up the microstructure of many suspensions for which the macroscopic dynamics are not well understood. The effect of wall on the fiber dynamics is the main focus of this work. For a single fiber undergoing simple shear flow between plane parallel walls the period of rotation was compared with the Jeffrey’s orbit. A fiber placed close to the wall shows significant deviation from Jeffrey’s orbit. The fiber moving near a solid wall in bounded shear flow follows a pole-vaulting motion, and its centroid location from the wall is also periodic. Simulations were also carried out to study the effect of fiber–fiber interactions on the viscosity of concentrated suspensions.     

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.     

Numerical prediction of fiber motion in a branching channel flow of fiber suspensions

Fiber orientation and dispersion in the dilute fiber suspension that flows through a T-shaped branching channel are simulated numerically based on the slender-body theory. The simulated results are consistent qualitatively with the experimental data available in the literature. The results show that the spatial distribution of fibers is dependent on the fiber aspect ratio, but has no relation with the volume fraction of fiber. The content ratio of fibers near the upper wall increases monotonically with an increasing Re number, and the situation is reverse for the region near the bottom wall. The orientation of fibers depends on Re number, however, the function of fiber volume fraction and aspect ratio is negligible. The fibers near the wall and in the central region of the channel align along the flow direction at all times, but the fibers in the other parts of the channel tend to align along the flow direction only in the downstream region.The project supported by the National Natural Science Foundation of China (10372090) and Doctoral Program of Higher Education in China (20030335001)The English text was polished by Ron Marshall     

EFFECTS OF TENSOR CLOSURE MODELS AND 3-D ORIENTATION ON THE STABILITY OF FIBER SUSPENSIONS IN A CHANNEL FLOW

IntroductionFlowoffibresuspensionshasbeenveryfamiliarinmanyindustrialfields.Fibreadditivesplayanimportantroleindragreductioninmanytypesofflow[1- 3].Inthesuspensions,somebehavioroftheflowmaybealteredbythefibres.Oneoftheimportantexamplesisthehydrodynamicsta…     

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.     

Degradation effects of dilute polymer solutions on turbulent drag reduction in pipe flows

An important practical problem in the application and study of drag reduction by polymer additives is the degradation of the polymer, for instance due to intense shearing, especially in recirculatory flow systems. Such degradation leads to a marked loss of the drag-reducing capability of the polymer.Three different polymer types were tested on degradation effects in a closed pipe flow system. The polymers used were Polyox WSR-301, Separan AP-273 and Superfloc A-110, dissolved in water in concentrations of 20 wppm each. The flow system consisted of a 16.3 mm pipe of 4.25 m length. Two different pumps were used: a centrifugal pump and a disc pump. Different solution-preparation procedures were tried and the experiments were performed at different flow rates.Superfloc A-110 proved to be both the most effective drag reducer and most resistant to degradation. Because of very fast degradation, Polyox WSR-301 was found to be unsuitable for being used as a drag reducer in re-circulatory systems. The disc pump proved to be much better suited for pumping the polymer solutions than the centrifugal pump. The degradation curve of the combination Superfloc/disc pump showed a plateau-like region with reasonable drag reduction, which makes it possible to perform (laser Doppler) measurements under nearly constant circumstances during a sufficient time.     

内管自转且公转时环管幂律流的受力分析

数值模拟了环管中内管偏心自转且公转时由轴向压力所驱动的幂律流体充分发展层流,分析了内管上的流体作用力。结果表明,内管偏心自转时流体作用力具有推动内管作和自转同向公转的效果。当只有外力矩驱动内管自转时,由于流体的作用,随内管线密度的不同,内管能达到的受力平衡态也不同:线密度较小时内管仅能在同心自转时达到受力平衡;线密度较大时内管能在作具有不变角速度和偏心率公转时达到受力平衡,且内管线密度越大,对应的受力平衡的公转的偏心率也越大。     

NUMERICAL SIMULATION OF COMPLICATED PIPE FLOW WITH k-εMODEL

In this paper, turbulence in a complicated pipe is simulated by using the k-ε model. The ladder-like mesh approximation is used to solve the problem of complicated boundary with the result of numerical simulation favorable. Two computational examples are given to validate the strong adaptability and stability of k-ε model.     

湍流横掠固定颗粒的全尺度直接数值模拟

气固两相流模拟中,当固相尺度接近或大于Kolmogorov尺度时,普通的点源模型将不再适用,固体相的体积效应和表面效应将对流体相产生显著的影响。通过采用直接数值模拟方法,结合内嵌边界方法对湍流中不同湍流强度流体横掠大于Kolmogorov尺度的固相颗粒进行了全尺度模拟,讨论分析了在两种湍流度下方形颗粒对湍流的调制影响以及颗粒的受力情况。     

Interaction of a gas-droplet turbulent jet with a cocurrent high-velocity high-temperature gas flow

A mathematical model and a method for calculating a gas-droplet turbulent jet with allowance for velocity nonequilibrium and virtual mass of the condensed phase during turbulent fluctuations and also heat and mass transfer within the three-temperature scheme are developed. Methodical calculations are performed. The results of these calculations are in reasonable agreement with available experimental data. The structure of the gas-droplet jet in a cocurrent high-velocity high-temperature gas flow is studied by numerical methods. The ratio of intensities of heat and mass transfer between the phases and turbulent diffusion transfers of substances is found to be different at the initial, transitional, and basic segments of the jet. This difference is responsible for the nonmonotonic axial distribution of vapor density and the lines of the halved mass flow of the condensed phase. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 85–94, May–June, 2008.     

Analysis of pipe flow transition. Part II. Energy transfer

The present article describes the results from a study of nonlinear mechanisms at work during the process of transition to turbulence in pipe flows. Using an accurate hybrid finite-difference code for the simulation of unsteady incompressible pipe flow, we have performed a direct numerical simulation designed to model experiments performed by Han, Tumin and Wygnanski [12]. Based on these numerical data, we have conducted a meticulous investigation of the dynamic interactions of the structures and flow modes that can be observed during this process. Based on this study, we can paint a detailed picture of the dynamical interactions of flow structures during both the linear and nonlinear stages of pipe flow transition. While this picture does have some similarities to earlier proposed mechanisms, we find that even for the simple cases considered here the structure of the pertinent interactions is much richer than suggested by these earlier models.     

Orientation distribution of fibers and rheological property in fiber suspensions flowing in a turbulent boundary layer

A model relating the translational and rotational transport of orientation distribution function （ODF） of fibers to the gradient of mean ODF and the dispersion coefficients is proposed to derive the mean equation for the ODE Then the ODF of fibers is predicted by numerically solving the mean equation for the ODF together with the equations of turbulent boundary layer flow. Finally the shear stress and first normal stress difference of fiber suspensions are obtained. The results, some of which agree with the available relevant experimental data, show that the most fibers tend to orient to the flow direction. The fiber aspect ratio and Reynolds number have significant and negligible effects on the orientation dis- tribution of fibers, respectively. The additional normal stress due to the presence of fibers is anisotropic. The shear stress of fiber suspension is larger than that of Newtonian solvent, and the first normal stress difference is much less than the shear stress. Both the additional shear stress and the first normal stress difference increase with increasing the fiber concentration and decreasing fiber aspect ratio.     

Research on the effect of cylinder particles on the turbulent properties in particulate flows

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｔｗｏ＿ｐｈａｓｅｆｌｏｗｗｉｔｈｃｙｌｉｎｄｅｒｐａｒｔｉｃｌｅｓｉｓａｄｉｆｆｉｃｕｌｔｍａｔｔｅｒｉｎｔｈｅｏｒｅｔｉｃａｌｒｅｓｅａｒｃｈｏｆｍｕｌｔｉｐｈａｓｅｆｌｏｗａｎｄｎｏｎ＿Ｎｅｗｔｏｎｉａｎｆｌｏｗ ,ｓｏｉｔｉｓａｃａｄｅｍｉｃａｌｌｙｖａｌｕａｂｌｅｔｏｓｔｕｄｙｉｔ.Ｍｅａｎｗｈｉｌｅ ,ｓｕｃｈｆｌｏｗｓｈａｖｅａｂｒｏａｄｂａｃｋｇｒｏｕｎｄｉｎｉｎｄｕｓｔｒｙ .Ｉｎｔｈｅｍａｔｅｒｉａｌｓｃｉｅｎｃｅ,ｔｈｅｍｏｌｄｉｎｇａｎｄｐｒｏｃｅｓ…     

Laser-Doppler measurements of laminar and turbulent flow in a pipe bend

Laser-Doppler measurements are reported for laminar and turbulent flow through a 90° bend of circular cross-section with mean radius of curvature equal to 2.8 times the diameter. The measurements were made in cross-stream planes 0.58 diameters upstream of the bend inlet plane, in 30, 60 and 75° planes in the bend and in planes one and six diameters downstream of the exit plane. Three sets of data were obtained: for laminar flow at Reynolds numbers of 500 and 1093 and for turbulent flow at the maximum obtainable Reynolds number of 43 000. The results show the development of strong pressure-driven secondary flows in the form of a pair of counter-rotating vortices in the streamwise direction. The strength and character of the secondary flows were found to depend on the thickness and nature of the inlet boundary layers, inlet conditions which could not be varied independently of Reynolds number. The quantitative anemometer measurements are supported by flow visualization studies. Refractive index matching at the fluid-wall interface was not used; the measurements consist, therefore, of streamwise components of mean and fluctuating velocities only, supplemented by wall pressure measurements for the turbulent flow. The displacement of the laser measurement volume due to refraction is allowed for in simple geometrical calculations. The results are intenden for use as benchmark data for calibrating flow calculation methods.