Transition of Stratified Flow in Horizontal or Slightly Inclined Pipe with Elliptical Cross Section

The study concerns the separated two-phase gas-liquid flow in a horizontal or slightly inclined pipe. The cross section of the cylindrical pipe is elliptical. After deriving the partial differential equations, the analytical prediction of the transition between stratified/non-stratified flows is obtained by using a linear stability analysis. The stratified flow, in a pipe having one of its cross sections axes disposed horizontally, is more stable when this axis is the great major.     

Start-up Flow of a Bingham Fluid in a Pipe

We present an analytical solution of axisymmetric motion for a Bingham fluid initially at rest subjected to a constant pressure gradient applied suddenly. Using the Laplace transform, we obtain expressions which allow the calculation of the instantaneous velocity, plug radius and rate of flow as a function of time. We also give a relation for the shear stress in the plug and in the region where the behaviour of the fluid is Newtonian.     

Transport and deposition of neutral particles in magnetohydrodynamic turbulent channel flows at low magnetic Reynolds numbers

The effect of Lorentz force on particle transport and deposition is studied by using direct numerical simulation of turbulent channel flow of electrically conducting fluids combined with discrete particle simulation of the trajectories of uncharged, spherical particles. The magnetohydrodynamic equations for fluid flows at low magnetic Reynolds numbers are adopted. The particle motion is determined by the drag, added mass, and pressure gradient forces. Results are obtained for flows with particle ensembles of various densities and diameters in the presence of streamwise, wall-normal or spanwise magnetic fields. It is found that the particle dispersion in the wall-normal and spanwise directions is decreased due to the changes of the underlying fluid turbulence by the Lorentz force, while it is increased in the streamwise direction. The particle accumulation in the near-wall region is diminished in the magnetohydrodynamic flows. In addition, the tendency of small inertia particles to concentrate preferentially in the low-speed streaks near the walls is strengthened with increasing Hartmann number. The particle transport by turbophoretic drift and turbulent diffusion is damped by the magnetic field and, consequently, particle deposition is reduced.     

Importance of unsteady contributions to force and heating for particles in compressible flows. Part 2: Application to particle dispersal by blast waves

Particle dispersal by blast waves is an interesting phenomenon. A model problem, i.e., a sudden release of a compressed gas–particle mixture contained in a spherical container, is employed to investigate the fundamental physics of particle dispersal. The problem is simulated by the multiphase flow models proposed in Part 1 of this article that include unsteady contributions in momentum and energy exchange between gas and particles. At early times, when particles are accelerated in the expansion fan, unsteady force and heating contributions are much larger than the corresponding quasi-steady contributions. Consequently, neglecting unsteady contributions leads to significant errors in particle front location (the boundary of the particle cloud). The complex wave interactions in the flow have strong influence on the particle motion. As a result, the particle motion is a non-monotonic function of particle density or diameter and the evolution of particle concentration is non-uniform and unsteady.     

Finite-length effects on dynamical behavior of rod-like particles in wall-bounded turbulent flow

Combined Particle Image Velocimetry (PIV) and Particle Tracking Velocimetry (PTV) measurements have been performed in dilute suspensions of rod-like particles in wall turbulence. PIV results for the turbulence field in the water table flow apparatus compared favorably with data from Direct Numerical Simulations (DNS) of channel flow turbulence and the universality of near-wall turbulence justified comparisons with DNS of fiber-laden channel flow. In order to examine any shape effects on the dynamical behavior of elongated particles in wall-bounded turbulent flow, fibers with three different lengths but the same diameter were used. In the logarithmic part of the wall-layer, the translational fiber velocity was practically unaffected by the fiber length l. In the buffer layer, however, the fiber dynamics turned out to be severely constrained by the distance z to the wall. The short fibers accumulated preferentially in low-speed areas and adhered to the local fluid speed. The longer fibers (l/z > 1) exhibited a bi-modal probability distribution for the fiber velocity, which reflected an almost equal likelihood for a long fiber to reside in an ejection or in a sweep. It was also observed that in the buffer region, high-speed long fibers were almost randomly oriented whereas for all size cases the slowly moving fibers preferentially oriented in the streamwise direction. These phenomena have not been observed in DNS studies of fiber suspension flows and suggested l/z to be an essential parameter in a new generation of wall-collision models to be used in numerical studies.     

Numerical study of short-circuiting flow and particles in a gas cyclone

Short-circuiting flow is an important secondary flow in gas cyclones, which has a negative impact on the separation performance. To improve the understanding of the short-circuiting flow and guide the optimization of gas cyclones, this paper presents a numerical study of a cyclone using computational fluid dynamics. Based on the steady flow field, three methods were adopted to investigate the formation mechanism and characteristics of the short-circuiting flow and particles. The temporal variation of the tracer species concentration distribution reveals that the formation mechanism of the short-circuiting flow is the squeeze between the airflows entering the annular space of the gas cyclone at different times. The short-circuiting flow region, distinguished through the spatial distribution of the moments of age, is characterized by a small mean age and a large coefficient of variation. The proportion of the short-circuiting particles increases with the increase of the inlet velocity only for small particles. But with the increase of particle size, the proportion of the short-circuiting particles decreases faster at higher inlet velocities, resulting in significant differences in collection efficiency curves.     

Experimental and CFD investigation of flow behavior and sand erosion pattern in a horizontal pipe bend under annular flow

The internal erosion of pipelines in oil and gas storage and transportation engineering is highly risky. High gas velocity of annular flow entrained sand will cause damage to the pipelines, and may further result in thinning of the wall. If this damage lasts for a long time, it may cause pipeline leakage and cause huge economic losses and environmental problems. In this research, an experimental device for studying multiphase flow erosion is designed, including an erosion loop and an experimental elbow that can test the erosion rate. The annular flow state and pipe wall erosion morphology can also be tested by the device. The computational fluid dynamics (CFD) method is combined with the experiment to further study the annular flow erosion mechanism in the pipeline. The relationship between gas-liquid-solid distribution and erosion profile was studied. The results show that the most eroded region occurs between 22.5° and 45° in the axial angle direction and between 90° and 135° in the circumferential angle direction of the elbow. The pits and deep scratches form on the surface of the sample after the sand collision.     

Stokes eigenfunctions and Galerkin projection of the disturbance equations in plane Poiseuille flow: a systematic analytical approach

In the present paper, the L ²-normalized Stokes eigenfunctions for plane Poiseuille flow, which form an orthonormal functional basis for the space of disturbances, are written in a general exponential form. Then, the evolution equations for the disturbances are Galerkin-projected on the considered basis functions, and all the terms of the resulting dynamical system are expressed systematically in analytical form. Finally, a numerical example is given in which the proposed basis functions are used for the simulation of the time evolution of the critical disturbance predicted by the energetic stability theory.     

高导无氧铜在大变形、不同温度和不同应变率下的流动应力和本构模型

为了理解高导无氧铜(OFHC Cu)的塑性流动行为,采用Instron液压试验机和分离式Hopkinson压杆,系统地对OFHC Cu进行了温度为77 ~1 000 K,应变率为0.001 ~7 000 s-1,以及真实应变超过80%的单轴压缩试验。结果表明：在0.001 s-1应变率下, OFHC Cu在约500 K呈现动态应变时效现象。随应变率增高,动态应变时效温度区域向更高温度移动,甚至动态应变时效现象消失。在高应变变形区域,相对温度来说,OFHC Cu塑性流动应力对应变率依赖更强。基于位错运动学和动力学概念,考虑位错在高温和高应变率的粘-曳阻力现象,结合试验结果,导出一个基于物理概念的本构模型。此模型可预测从低到高不同应变率不同温度下OFHC Cu的塑性流动应力。通过比较表明,本构模型预测结果与试验结果吻合较好。     

Dynamic coupling of fluid and structural mechanics for simulating particle motion and interaction in high speed compressible gas particle laden flow

In this work, structural finite element analyses of particles moving and interacting within high speed compressible flow are directly coupled to computational fluid dynamics and heat transfer analyses to provide more detailed and improved simulations of particle laden flow under these operating conditions. For a given solid material model, stresses and displacements throughout the solid body are determined with the particle–particle contact following an element to element local spring force model and local fluid induced forces directly calculated from the finite volume flow solution. Plasticity and particle deformation common in such a flow regime can be incorporated in a more rigorous manner than typical discrete element models where structural conditions are not directly modeled. Using the developed techniques, simulations of normal collisions between two 1 mm radius particles with initial particle velocities of 50–150 m/s are conducted with different levels of pressure driven gas flow moving normal to the initial particle motion for elastic and elastic–plastic with strain hardening based solid material models. In this manner, the relationships between the collision velocity, the material behavior models, and the fluid flow and the particle motion and deformation can be investigated. The elastic–plastic material behavior results in post collision velocities 16–50% of their pre-collision values while the elastic-based particle collisions nearly regained their initial velocity upon rebound. The elastic–plastic material models produce contact forces less than half of those for elastic collisions, longer contact times, and greater particle deformation. Fluid flow forces affect the particle motion even at high collision speeds regardless of the solid material behavior model. With the elastic models, the collision force varied little with the strength of the gas flow driver. For the elastic–plastic models, the larger particle deformation and the resulting increasingly asymmetric loading lead to growing differences in the collision force magnitudes and directions as the gas flow strength increased. The coupled finite volume flow and finite element structural analyses provide a capability to capture the interdependencies between the interaction of the particles, the particle deformation, the fluid flow and the particle motion.     

煤仓内煤散料流动状态与力学行为影响因素

针对煤仓内煤散料流动问题及其力学行为,采用三维颗粒流模拟程序PFC3D建立了某型号煤仓与某种煤散料的离散元模型,简述了其力学模型与求解步骤,模拟分析了煤仓内煤散料卸料流动状态。通过分析水平向侧压力、颗粒速度场和接触力场,重点讨论了煤仓下部锥体内壁面摩擦系数、锥仓倾角和卸料口径等对煤散料颗粒流动状态和力学行为的影响。结果显示,深仓卸料流动为整体流动与中心流动混合状态,煤仓内壁摩擦系数、锥体倾角和卸料孔开口半径均对煤散料流动和水平侧压力有较大影响。     

Experimental-computational study of fibrous particle transport and deposition in a bifurcating lung model

Experiments carried out using a lung model with a single horizontal bifurcation under different steady inhalation conditions explored the orientation of depositing carbon fibers, and particle deposition fractions. The orientations of deposited fibers were obtained from micrographs. Specifically, the effects of the sedimentation parameter (γ), fiber length, and flow rate on orientations were analyzed. Our results indicate that gravitational effect on deposition cannot be neglected for 0.0228 < γ < 0.247. The absolute orientation angle of depositing fibers decreased linearly with increasing γ for values 0.0228 < γ < 0.15. Correspondence between Stokes numbers and γ suggests these characteristics can be used to estimate fiber deposition in the lower airways. Computer simulations with sphere-equivalent diameter models for the fibers explored deposition efficiency vs. Stokes number. Using the volume-equivalent diameter model, our experimental data for the horizontal bifurcation were replicated. Results for particle deposition using a lung model with a vertical bifurcation indicate that body position also affects deposition.     

Dispersion de particules solides en mouvement de saltation dans un écoulement turbulentDispersion of solid saltating particles in a turbulent flow

The solid particle dispersion in saltating motion is studied in an homogeneous turbulence and in a turbulent boundary layer. The fluid velocity along the particle trajectory is estimated using a continuous stochastic differential equation in which the correlation integral time takes into account gravity and inertia effects. As far as the boundary layer is concerned, the aerodynamic entrainment of particles and the rebound are modelised as random variables with Gaussian probability density functions. Compared with experimental results, the numerical results show good agreement for dispersion, although velocity fluctuations are slightly under evaluated. To cite this article: C. Aguirre et al., C. R. Mecanique 332 (2004).     

An experimental study of drag on a single tube and on a tube in an array under two-phase cross flow

In this work, the drag coefficient and the void fraction around a tube subjected to two-phase cross flow were studied for a single tube and for a tube placed in an array. The drag coefficients were determined by measuring the pressure distribution around the perimeter of the tube. Single tube drag data were taken when the tube was held both rigidly and flexibly. The test tube was made of acrylic and was 2.2 cm in diameter and 20 cm in length. In the experiments, liquid Reynolds number ranged from 430 to 21,900 for the single tube and liquid gap Reynolds number ranged from 32,900 and 61,600 for the tube placed in a triangular array. Free stream void fraction was varied from 0 to 0.4. At low Reynolds numbers, the ratio of two-phase to single-phase drag coefficient is found to be a strong function of εGr/Re2. However, at high Reynolds numbers only void fraction is the important parameter. Empirical correlations have been developed for the ratio of two-phase drag on a single tube and on a tube placed in an array.     

Comparison of spectral and finite element methods applied to the study of the core‐annular flow in an undulating tube

A Galerkin/finite element and a pseudo‐spectral method, in conjunction with the primitive (velocity‐pressure) and streamfunction‐vorticity formulations, are tested for solving the two‐phase flow in a tube, which has a periodically varying, circular cross section. Two immiscible, incompressible, Newtonian fluids are arranged so that one of them is around the axis of the tube (core fluid) and the other one surrounds it (annular fluid). The physical and flow parameters are such that the interface between the two fluids remains continuous and single‐valued. This arrangement is usually referred to as Core‐Annular flow. A non‐orthogonal mapping is used to transform the uneven tube shape and the unknown, time dependent interface to fixed, cylindrical surfaces. With both methods and formulations, steady states are calculated first using the Newton–Raphson method. The most dangerous eigenvalues of the related linear stability problem are calculated using the Arnoldi method, and dynamic simulations are carried out using the implicit Euler method. It is shown that with a smooth tube shape the pseudo‐spectral method exhibits exponential convergence, whereas the finite element method exhibits algebraic convergence, albeit of higher order than expected from the relevant theory. Thus the former method, especially when coupled with the streamfunction‐vorticity formulation, is much more efficient. The finite element method becomes more advantageous when the tube shape contains a cusp, in which case the convergence rate of the pseudo‐spectral method deteriorates exhibiting algebraic convergence with the number of the axial spectral modes, whereas the convergence rate of the finite element method remains unaffected. Copyright © 2002 John Wiley & Sons, Ltd.     

Influence of drag laws on pressure and bed material recirculation rate in a cold flow model of an 8 MW dual fluidized bed system by means of CPFD

A cold flow model of an 8 MW dual fluidized bed (DFB) system is simulated using the commercial computational particle fluid dynamics (CPFD) software package Barracuda. The DFB system comprises a bubbling bed connected to a fast fluidized bed with the bed material circulating between them. As the hydrodynamics in hot DFB plants are complex because of high temperatures and many chemical reaction processes, cold flow models are used. Performing numerical simulations of cold flows enables a focus on the hydrodynamics as the chemistry and heat and mass transfer processes can be put aside. The drag law has a major influence on the hydrodynamics, and therefore its influence on pressure, particle distribution, and bed material recirculation rate is calculated using Barracuda and its results are compared with experimental results. The drag laws used were energy-minimization multiscale (EMMS), Ganser, Turton–Levenspiel, and a combination of Wen–Yu/Ergun. Eleven operating points were chosen for that study and each was calculated with the aforementioned drag laws. The EMMS drag law best predicted the pressure and distribution of the bed material in the different parts of the DFB system. For predicting the bed material recirculation rate, the Ganser drag law showed the best results. However, the drag laws often were not able to predict the experimentally found trends of the bed material recirculation rate. Indeed, the drag law significantly influences the hydrodynamic outcomes in a DFB system and must be chosen carefully to obtain meaningful simulation results. More research may enable recommendations as to which drag law is useful in simulations of a DFB system with CPFD.