RESEARCH ON THE MOTION OF PARTICLES IN THE TURBULENT PIPE FLOW OF FIBER SUSPENSIONS

The motion of fibers in turbulent pipe flow was simulated by 3-D integral method based on the slender body theory and simplified model of turbulence.The orientation distribution of fibers in the computational area for different Re numbers was computed.The results which were consistent with the experimental ones show that the fluctuation velocity of turbulence cause fibers to orient randomly.The orientation distributions become broader as the Re number increases.Then the fluctuation velocity and angular velocity of fibers were obtained.Both are affected by the fluctuation velocity of turbulence.The fluctuation velocity intensity of fiber is stronger at longitudinal than at lateral,while it was opposite for the fluctuation angular velocity intensity of fibers.Finally,the spatial distribution of fiber was given.It is obvious that the fiber dispersion is strenghened with the increase of Re numbers.

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.

On the structural features of fiber suspensions in converging channel flow

The structural features of fiber suspensions are dependent on the fiber alignment in the flows. In this work the orientation distribution function and orientation tensors for semi-concentrated fiber suspensions in converging channel flow were calculated, and the evolutions of the fiber alignment and the bulk effective vis-cosity were analyzed. The results showed that the bulk stress and the effective viscosity were functions of therate-of-strain tensor and the fiber orientation state ; and that the fiber suspensions evolved to steady alignment and tended to concentrate to some preferred directions close to but not same as the directions of local stream-lines. The bulk effective viscosity depended on the product of Reynolds number and time. The decrease of ef-fective viscosity near the boundary benefited the increase of the rate of flow. Finally when the fiber alignment went into steady state, the structural features of fiber suspensions were not dependent on the Reynolds numberbut on the converging channel angle.

Numerical Research on the Fiber Suspensions in a Turbulent T-shaped Branching Channel Flow

The concentration and orientation of fiber in a turbulent T-shaped branching channel flow are investi-gated numerically. The Reynolds averaged Navier-Stokes equations together with the Reynolds stress turbulent model are solved for the mean flow field and the turbulent kinetic energy. The fluctuating velocities of the fluid are assumed as a random variable with Gaussian distribution whose variance is related to the turbulent kinetic energy. The slender-body theory is used to simulate the fiber motion based on the known mean and fluctuating velocities of the fluid. The results show that at low Reynolds number, fiber concentration is high in the flow separation regions, and fiber orientation throughout the channel is widely distributed with a slight preference of aligning along the horizontal axis. With increasing of Re, the high concentration region disappears, and fiber orientation becomes ho-mogeneous without any preferred direction. At high Reynolds number, fiber concentration increases gradually along the flow direction. The differences in the distribution of concentration and orientation between different fiber aspect ratio are evident only at low Re. Both Re and fiber aspect ratio have small effect on the variance of orientation angle.

RESEARCH ON DISTRIBUTION FUNCTION AND ORIENTATION TENSORS OF FIBER SUSPENSIONS IN WEDGE-SHAPED FLOW FIELD

The distribution function and orientation tensor of fiber suspensions in wedge shaped flow field were computed . The results indicate that with time increasing, the fiber orient themselves to flow direction gradually. At the locations with same pole radii, the angle between fiber orientation and centerline, which occurs with the most probability, decreases with the pole angle increasing . At the locations with same pole angle, the angle between fiber orientation and centerline increases with the pole radii decreasing. The second order tensors get steady more quickly at the points where the velocity g radients are larger. At the locations with same pole angle, the steady values of orientation tensors are identical. At the locations with same pole radii is, streamline becomes flatter as the pole angle decreases.

FIBER ORIENTATION DISTRIBUTIONS IN SLIT CHANNEL FLOWS WITH ABRUPT EXPANSION FOR FIBER SUSPENSIONS

The lattice Boltzmann method was used to investigate numerically the fiber orientation distributions in slit channel flows with abrupt expansion for fiber suspensions even in the concentrated regime. The channels have a thin slit geometry with 1：4 and 1：3 expansions. Both the interactions between fibers and that between fibers and channel walls were taken into consideration. Some of numerical results are qualitatively in agreement with the experiment data. It is found that most of fibers are aligned in the flow direction in all the suspensions in the entrance region of the expansion. Fiber orientation distributions, having different patterns in different regions of the flow, depend on the expansion ratio of the channel. The mechanical fiber-fiber interaction largely affects the fiber orientation in the downstream of the expansion and in the salient corner for the cases of concentrated suspensions. The hydrodynamic interaction plays an important role on the fiber orientation in the dilute suspension.

Research on the orientation distribution of fibers immersed in a pipe flow

The computed orientation distribution of fibers immersed in laminar pipe flows showed that the longitudinal distributions are wide for small Reynolds numbers and become narrower with increasing Re.For low Re number.the axial orientation distributions are broad with almost no preferred orientations.For high Re numer,the axial distribution becomes narrow.with sharp maxima..The mean values of the longitudinal orientation depend strongly on the Re number.The computed results are in qualitative agreemernt with relevant experimental results.

Research on the behavior of fiber orientation probability distribution function in the planar flows

The equation of two-dimensional fiber direction vector was solved theoretically to give the fiber orientation distri- bution in simple shear flow, flow with two direction shears, extensional flow and arbitrary planar incompressible flow. The Fok- ker-Planck equation was solved numerically to validify the theoretical solutions. The stable orientation and orientation period of fiber were obtained. The results showed that the fiber orientation distribution is dependent on the relative not absolute magnitude of the matrix rate-of-strain of flow. The effect of fiber aspect ratio on the orientation distribution of fiber is insignificant in most conditions except the simple shear case. It was proved that the results for a planar flow could be generalized to the case of 3-D fiber direction vector.

Science Letters:Research on 3D fiber orientation distribution in arbitrary planar flows

A non-stretchable fiber rotation in planar flows has been solved. The fiber will rotate periodically or run to the asymptotical direction decided by a discriminant defined in the paper involving the fiber aspect ratio and the flow characteristics. Subsequently the fiber orientation distribution is derived directly without the bother of solving the Fokker-Planck equation. The research clearly indicates the overall configuration of a fiber rotation movement in planar flows.

Orientation of the fiber suspending in the flow through a tube containing a sphere

Fiber suspensions flow through a tube containing a sphere in the dilute and concentrated regimes is simulated numerically with the Lattice Boltzmann Method （LBM）. The numerical results of fiber orientation distribution based on a statistical scheme are obtained and agree qualitatively with the experimental ones for the flow through a parallel plate channel comaining a cylinder. The results show that the sphere in the tube results in a change in the fiber orientation distribution downstream of the sphere along the flow and transverse directions. The influences of the sphere on the fiber orientation distribution are more significant for the concentrated suspensions than for the dilute one. The effect of the initial fiber orientations on the fiber orientation distribution is significant upstream of the sphere but small downstream of the sphere.

圆管纤维悬浮湍流场中粒子运动的研究

利用从细长体理论出发得到的三维分段积分法和湍流简化方法模拟了大量纤维粒子在圆管湍流内的运动· 统计了不同Re数下计算区域内的纤维的取向分布,计算结果与实验结果基本吻合,结果表明湍流的脉动速度导致纤维取向趋于无序,且随着Re数的增加,纤维取向的分布越来越趋于均匀· 其后又考虑了纤维速度和角速度的脉动,二者都充分体现了流体速度脉动的影响,且纤维速度的脉动在流向上的强度大于横向,而其角速度的脉动在流向上的强度小于横向· 最后统计了纤维在管道截面上的位置分布。

纤维悬浮槽道湍流中纤维取向分布和流变特性

为了解悬浮槽道湍流中纤维取向分布和流变特性,对纤维悬浮湍流进行了理论分析和数值模拟.推导了纤维取向平均概率分布函数方程,得到了脉动流场速度和脉动取向概率分布函数关联项,建立了关联项与平均概率分布函数梯度以及扩散系数之间的关系,数值求解了平均概率分布函数的方程.计算结果表明,随着长径比的增大,纤维的轴向更趋向于流动方向;悬浮流的剪切应力和第一法向应力差由壁面到槽道中心逐渐接近于0;剪切应力几乎不随长径比和体积分数变化,而第一法向应力差的绝对值随体积分数的增大和长径比的减小而增大.