Pipe flow of suspensions

This study shows that fully developed pipe flow of a particulate suspension is defined by four dimensionless parameters of particle-fluid interactions in addition to the Reynolds number. Effects accounted for include the Magnus effect due to fluid shear, electrostatic repulsion due to electric charges on the particles, and Brownian or turbulent diffusion. In the case of a laminar liquid-solid suspension, electrostatic effect is negligible but shear effect is prominent. Solution of the basic equations gives the density distribution of particles with a peak at the center (Einstein, Jeffery) or at other radii between the center and the pipe wall (Segré et al) depending on the magnitudes of the various flow parameters. In the case of a turbulent gas-solid suspension, the Magnus effect is significant only within the thickness of the laminar sublayer. However, charges induced on the particles by the impact of particles at the wall produce a higher density at the wall than at the center of the pipe. The velocity distribution of particles is characterized by a slip velocity at the wall and a lag in velocity in the core from the fluid phase. These results are verified by earlier measurements.     

Measurement of local size and velocity probability density distributions in two-phase suspension flows by Laser-Doppler technique

A technique is presented for the simultaneous measurement of the local number and velocity probability densities of a dilute two-phase suspension which has a distribution of particle sizes and a predominate direction of flow orientation such as in the cases of pipe and boundary-layer flows. It is shown that by a suitable scheme of discrimination on the amplitude as well as the residence time and frequency of the individual Laser-Doppler bursts, one can obtain the statistics on the size number density distribution and, for each size range, velocity distribution of the particulate phase together with the velocity probability distribution of the fluid phase.Results have been obtained for experiments conducted on a laminar uniform flow and a turbulent shear flow of a dilute glass particle-water suspension having a particle size distribution. Calibration needed for the scheme was accomplished by analyzing particle size and number density distribution data obtained from a Coulter particle sizing counter on a sample taken with an isokinetic probe.     

流体黏性对输流管道运动方程及临界流速的影响

考虑实际流体黏性引起的管内流速非均匀分布,针对层流和两种不同的湍流流态,对理想流体情况下输流管道运动方程中的离心力项进行了修正,得到的修正系数分别为1.333(圆管层流)、1.020(光滑管壁圆管湍流)和1.037～1.055(粗糙管壁圆管湍流).根据修正后的运动方程得到的上述3种情况下的发散失稳临界流速比理想流体流动情况下依次分别低13.4%,1.0%和1.8%～2.6%.流体黏性对输流管道运动方程及临界流速的影响只与流态有关,雷诺数决定流态,而黏性系数通过雷诺数间接起作用.     

Single-Fluid Model of a Mixture for Laminar Flows of Highly Concentrated Suspensions

A model of laminar flow of a highly concentrated suspension is proposed. The model includes the equation of motion for the mixture as a whole and the transport equation for the particle concentration, taking into account a phase slip velocity. The suspension is treated as a Newtonian fluid with an effective viscosity depending on the local particle concentration. The pressure of the solid phase induced by particle-particle interactions and the hydrodynamic drag force with account of the hindering effect are described using empirical formulas. The partial-slip boundary condition for the mixture velocity on the wall models the formation of a slip layer near the wall. The model is validated against experimental data for rotational Couette flow, a plane-channel flow with neutrally buoyant particles, and a fully developed flow with heavy particles in a horizontal pipe. Based on the comparison with the experimental data, it is shown that the model predicts well the dependence of the pressure difference on the mixture velocity and satisfactorily describes the dependence of the delivered particle concentration on the flow velocity.     

Rotation effects on a fully-developed turbulent pipe flow

A fully-developed turbulent pipe flow is allowed to pass through a rotating pipe section, whose axis of rotation coincides with the pipe axis. At the exit end of the rotating section, the flow passes into a stationary pipe. As a result of the relaxation of surface rotation, the turbulent flow near the pipe wall is affected by extra turbulence production created by the large circumferential shear strain set up by the rapid decrease of the rotational velocity to zero at the wall. However, the flow in the most part of the pipe is absent of this extra turbulence production because the circumferential strain is zero as a result of the solid-body rotation imparted to the flow by the rotating pipe section. The combined effect of these two phenomena on the flow is investigated in detail using hot-wire anemometry techniques. Both mean and turbulence fields are measured, together with the wall shear and the turbulent burst behavior at the wall. A number of experiments at different rotational speeds are carried out. Therefore, the effects of rotation on the behavior of wall shear, turbulent burst at the wall, turbulence production and the near-wall flow can be documented and analysed in detail.     

On the motion of particles in turbulent duct flows

Existing knowledge on particle deposition rates on walls from turbulent pipe and channel flows is summarized and it is shown that discrepancies exist between experimental and theoretical findings. To contribute to the existing experimental information, laser Doppler measurements are reported of the flow field of a glass particle-air two-phase flow. The results reveal certain seemingly peculiar behaviors of the particles which obviously defy the predictions of the conventional analyses of turbulent two-phase suspension flows.In an accompanying approximate, yet pragmatic theoretical approach, an attempt is made to find a rational basis for the explanation of these experimentally observed particle behaviors. It is shown for the particles in the present study, there exists a limiting size above which their response to the agitation of the fluctuating motion of the surrounding fluid could be treated as if the flow were laminar. On this rational basis, these experimentally observed particle behaviors can then be qualitatively explained by the existing theory of particle excursion in a laminar shear flow field.Reported also is a suggestion to extend the present analysis to a dispersion of particles of multiple sizes.     

输气管道壁面涂料减阻机理的实验研究

用IFA-300热线风速仪以高于对应最小湍流时间尺度的分辨率精细测量了风洞中不同壁面涂料的管道湍流边界层不同法向位置流向速度分量的时间序列信号,利用湍流边界层近壁区域对数律平均速度剖面与壁面摩擦速度、流体黏性系数等内尺度物理量的关系和壁面摩擦速度与壁面摩擦切应力的关系,在准确测量湍流边界层近壁区域对数律平均速度剖面的基础上,间接测量湍流边界层的壁面摩擦阻力．对不同壁面涂料的壁湍流脉动速度信号用子波分析进行多尺度分解,用子波系数的瞬时强度因子和平坦因子检测管道湍流边界层中的多尺度相干结构,提取不同尺度相干结构的条件相位平均波形,对比研究输气管道壁面涂料的减阻机理．     

暖季强降雨对多年冻土南界斜坡路基稳定性影响分析

为阐明表面活性剂水溶液的减阻作用,使用LDV对零压梯度的二维湍流平板边界层中的CTAB 表面活性剂水溶液的湍流特性进行了实验研究. 结果表明：与牛顿流体相比,CTAB水溶液边 界层的粘性底层增厚;主流时均速度分布有被层流化的趋势,对数分布域上移;主流方向速 度湍动强度峰值减小,且远离壁面,在靠近边界层中部,出现第2峰值;垂直于主流方向的 速度湍动强度受到了大幅度抑制,雷诺应力沿着边界层厚度方向几乎为零. 结果说明CTAB 水溶液具有减弱湍流湍动各个成分相关度的作用,从而能够使雷诺应力降低、湍流能量生成 项减小最终降低流体的输送动力.     

Particles migration in laminar boundary layer flow

The behavior of a suspension of neutrally buoyant solid particles in the moderate relative velocity region (close to the leading edge) in a laminar boundary layer along a flat plate, is investigated. The velocity fields of both the fluid and the solid particles and the density distribution of the solid particles are simultaneously obtained by the use of the Laser-Doppler technique.It is shown that in this longitudinal region of the two-phase suspension boundary layer, the solid particles which initially lead the fluid decelerate as they move downstream, until at a certain vertical distance they lag the fluid. Clearly defined one-way transverse migration of particles directed away from the plate is characteristic of this region.     

Particle drag in a dilute turbulent two-phase suspension flow

The drag between phases plays an important role in the study of a turbulent two-phase suspension flow and its physical understanding will greatly promote progress in theoretical treatments of a whole range of important industrial and technical problems involving such a flow. The conventional practice of using the results of measurements based on a single particle in a laminar stream for the case of a turbulent flow of a dilute suspension is questioned. An analysis of the results of local measurements of upward turbulent flows of a solid particle-air two-phase suspension leads to the determination of the realistic particle drag coefficient over a wide range of flow conditions. It is established that the particle drag can be described by the simple Stokes law, based on an apparent turbulent viscosity of the fluid for the particles in the suspension flow. A correlation is provided for this apparent turbulent viscosity in terms of the particle size and concentration in the suspension, the local flow turbulence Reynolds number and the particle-to-fluid density ratio.     

Theory on transverse migration of particles in a turbulent two-phase suspension flow due to turbulent diffusion—I

A new theoretical model has been developed to explain the behavior of transverse particle transport in turbulent flow of a dilute two-phase suspension due to turbulent diffusion. This model is based on the ability of a particle to respond to surrounding fluid motion and depends on particle size and density relative to the carrier fluid, the fractional variation in particle concentration in the transverse direction as well as the existing turbulence structure of the surrounding fluid. The model developed in this investigation has been formulated by dividing the transverse fluid velocity, as seen by a particular particle, into two superimposed components representing, respectively, the transverse turbulent fluid fluctuations and an apparent transverse local fluid drifting velocity due to the effect on the transverse oscillatory component of fluid motion by the transverse concentration distribution of particles. A subsequent paper will show that the theory (together with other new results on the concentration effects on particle drag and lift and fluid turbulence properties) can help to explain the phenomena measured previously.     

Direct numerical simulation of the turbulent energy balance and the shear stresses in power-law fluid flows in pipes

The results of direct numerical simulation of turbulent flows of non-Newtonian pseudoplastic fluids in a straight pipe are presented. The data on the distributions of the turbulent stress tensor components and the shear stress and turbulent kinetic energy balances are obtained for steady turbulent flows at the Reynolds numbers of 10⁴ and 2×10⁴. As distinct from Newtonian fluid flows, the viscous shear stresses turn out to be significant even far from the wall. In power-law fluid flows the mechanism of the energy transport from axial to transverse component fluctuations is suppressed. It is shown that with decrease in the fluid index the turbulent transfer of the momentum and the velocity fluctuations between the wall layer and the flow core reduces, while the turbulent energy flux toward the wall increases. The earlier-proposed models for the average viscosity and the non-Newtonian one-point correlations are in good agreement with the data of direct numerical simulation.     

Investigation of turbulence effects on forced convection in a composite porous/fluid duct: Constant wall flux and constant wall temperature cases

Forced convection heat transfer in composite porous/fluid domains is of great practical and theoretical significance. However, research in this area traditionally addressed only the laminar flow case in both homogeneous fluid and porous regions of the domain. This paper investigates the interaction between turbulent flow in the center of a circular tube filled with a homogeneous fluid and laminar flow in the porous layer adjacent to the tube wall. A two-layer algebraic turbulence model suggested by Cebeci and Smith is utilized for the flow in the central region of the tube. The effects of turbulence in the central region on velocity and temperature distributions as well as on the Nusselt number are analyzed.     

Pathline analysis of traveling wavy blowing and suction control in turbulent pipe flow for drag reduction

Skin friction drag is much greater in turbulent flows as compared with that in laminar flows. It is well known that traveling wave control can be used to achieve a large drag reduction. In the present study, a direct numerical simulation of a turbulent pipe flow was performed to clarify the mechanism of the drag reduction caused by the traveling wave control. The flow induced by the control was evaluated using pathline analysis. Near the wall, a “closed flow” was formed, wherein the injected particles return to the wall owing to the suction flow. The random component of Reynolds shear stress was perfectly suppressed in the closed flow, which suggests that there was no turbulence. The controlled flow was categorized into four patterns, and each flow characteristic and drag reduction effect was discussed. When the closing rate is high, the drag decreases, while when the closing rate is low, i.e., when the injected particles are released into the main flow, the turbulence is maintained. If the thickness of the layer suppressing turbulence is insufficient, a significant effect in terms of the drag reduction cannot be expected. The large drag reduction owing to the traveling wave control can be attributed to the elimination of turbulence in the region near the wall.     

An experimental study of two flows through an axisymmetric sudden expansion

Two turbulent separated and reattaching flows produced by a sudden expansion in a pipe have been studied. The first was produced by a simple axisymmetric sudden enlargement from a nozzle of diameter 80 mm to a pipe of diameter 150 mm. The second was the flow at the same enlargement with the addition of a centerbody 90 mm downstream of the nozzle exit. Detailed measurements of velocity and skin friction (made primarily using pulsed wires) and of wall static pressure are presented. Without the centerbody the flow structure is similar to that observed in other sudden pipe expansions and over backward-facing steps. A turbulent free shear layer, bearing some similarity to that of a round jet, grows from separation and then reattaches to the pipe wall downstream. Reattachment is a comparatively gradual process, the shear layer approaching the wall at a glancing angle. The introduction of the centerbody causes the shear layer to curve towards the wall and reattach at a much steeper angle. Reattachment is much more rapid; gradients of skin friction and pressure along the wall are many times those without the centerbody. The high curvature of the shear layer strongly influences its turbulent structure, locally suppressing turbulence levels and reducing its growth rate.     

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.     

Propagation of a pre-existing turbulent fluid fracture

The propagation of rough and smooth wall pre-existing turbulent fluid fractures is investigated. The laminar fluid fracture is included as a special case for comparison. Lubrication theory is assumed to apply in the fracture and turbulence is introduced through the wall shear stress. The Perkins–Kern–Nordgren approximation is made in which the fluid pressure is proportional to the half-width of the fracture. The fracture half-width satisfies a non-linear diffusion equation. By using a linear combination of the Lie point symmetries of the non-linear diffusion equation a group invariant solution for the fracture length, volume and half-width is derived. The evolution of the length, half-width and mean flow velocity is analysed for a range of working conditions at the fracture entry. It is found that the mean flow velocity increases approximately linearly along the fracture.     

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.     

Direct numerical simulation of particle-laden plane turbulent wall jet and the influence of Stokes number

The distribution and motion of inertial particles in plane turbulent wall jet are investigated using direct numerical simulation, under the assumption of one-way coupling. To our knowledge, this appears to be the first direct numerical simulation of a particle-laden plane turbulent wall jet. It is shown that, in outer part of the wall jet, the behaviour of particles closely resembles that of a free plane jet. Due to the streamwise decay of particle Stokes number, the particle streaks formed in the near wall region of the wall jet are characterized by their intensity variation, which differs significantly from those in the channel flow. The streamwise growth of the particle velocity half-width is approximately equal to that of the fluid velocity half-width and the maximum velocity of particles decays slower than that of fluid due to inertia. The outer scaling can collapse the mean particle velocity in both the inner and outer region for heavier particles. In the buffer region, the particle–fluid velocity difference can be negative or positive depending on the Stokes number since there are two competing effects, namely the memory effect and turbophoresis. In the viscous region, the larger particles are on average faster than fluid and the velocity difference is found to be self-similar depending on outer Stokes number. The near-wall distribution of velocity difference is significantly correlated with the presence of high-momentum particles which are entrained by vortical structures generated in the outer region of the wall jet. These results are useful for environmental and engineering applications.     

Direct numerical simulation of supersonic pipe flow at moderate Reynolds number

We study compressible turbulent flow in a circular pipe at computationally high Reynolds number. Classical related issues are addressed and discussed in light of the DNS data, including validity of compressibility transformations, velocity/temperature relations, passive scalar statistics, and size of turbulent eddies. Regarding velocity statistics, we find that Huang’s transformation yields excellent universality of the scaled Reynolds stresses distributions, whereas the transformation proposed by Trettel and Larsson (2016) yields better representation of the effects of strong variation of density and viscosity occurring in the buffer layer on the mean velocity distribution. A clear logarithmic layer is recovered in terms of transformed velocity and wall distance coordinates at the higher Reynolds number under scrutiny (Re_τ ≈ 1000), whereas the core part of the flow is found to be characterized by a universal parabolic velocity profile. Based on formal similarity between the streamwise velocity and the passive scalar transport equations, we further propose an extension of the above compressibility transformations to also achieve universality of passive scalar statistics. Analysis of the velocity/temperature relationship provides evidence for quadratic dependence which is very well approximated by the thermal analogy proposed by Zhang et al. (2014). The azimuthal velocity and scalar spectra show an organization very similar to canonical incompressible flow, with a bump-shaped distribution across the flow scales, whose peak increases with the wall distance. We find that the size growth effect is well accounted for through an effective length scale accounting for the local friction velocity and for the local mean shear.