含悬浮固粒的旋转射流剪切层稳定性研究

本文在理想不可压旋转圆射流的运动方程中添加了固粒作用项，由此推得了时间增长率的表达式，进而得关于含悬浮固粒放置射流稳定性的修正瑞利稳定性准则，求出了不同固粒质量密度固－气脉动速度比值，固气脉动速度相位差及Ｓｔｏｋｅｓ数情况下旋转射流场的增长率与径向空间波数的关系曲线，在比较这些曲线的基础上，给出了关于固粒属性对旋转射流场稳定性影响的几个重要结论为控制旋转射流场和后续发展提供了依据。     

Research on stability of moving jet containing dense suspended solid particles

Ｎｏｍｅｎｃｌａｔｕｒｅαｊ:ｐｈａｓｅｊｖｏｌｕｍｅｆｒａｃｔｉｏｎ;ｕｊ:ｐｈａｓｅｊｖｅｌｏｃｉｔｙｖｅｃｔｏｒ;Ｕ :ｍｅａｎｖｅｌｏｃｉｔｙｖｅｃｔｏｒ;Ｕ :ｍｅａｎａｘｉａｌｖｅｌｏｃｉｔｙ;Ｕｏ:ｍｅａｎｘ＿ａｘｉａｌｖｅｌｏｃｉｔｙｏｆｔｈｅｅｘｉｔ;Ｕ∞ :ｖｅｌｏｃｉｔｙｏｆｊｅｔｔｉｎｇｄｅｖｉｃｅ;ΔＵ :ｊｕｍｐｉｎｇｖｅｌｏｃｉｔｙｏｆｍｏｖｉｎｇｊｅｔ,ｃｈａｒａｃｔｅｒｉｓｔｉｃｖｅｌｏｃｉｔｙｏｆｆｌｏｗ＿ｆｉｅｌｄ ,ΔＵ =Ｕｏ-Ｕ∞ ;ｐ:ｐｒｅｓｓｕｒｅ;Ｐ :ｅｉｇｅｎｐｒｅｓ…     

Numerical solution and sensitivity analysis of filter cake permeability and resistance to model parameters

The migration and capture of solid particles in porous media occur in fields as diverse as water and wastewater treatment, well drilling, and in various liquid-solid separation processes. Filter cakes are formed when a liquid containing solid particles is forced through a pervious surface which allows the liquid transport while retaining solid particles. Following a literature survey, a governing equation for the cake thickness is obtained by considering the instantaneous mass balance. Later, numerical solutions for the cake thickness, cake permeability, cake resistance, solid particle velocity (cake compression rate) and concentration of suspended particles are obtained and a sensitivity analysis is conducted. The sensitivity analysis shows that the cake permeability and cake resistance are more sensitive to the rate constant of cake erosion than they are to the rate constant of particle capture. However, the concentration of suspended solid particles, and the solid velocity are mostly sensitive to the slurry parameter and the rate constant of particle trapping. Moreover, cake permeability, compressibility, concentration of suspended particles, and the solid velocity are very sensitive to the concentration at the filter septum. Finally, as expected, with a thicker slurry, more particles are captured inside the cake, thus forming a thicker and more resistant cake. Also, as more particles are being filtered at the filter septum, a thinner cake is formed and a smaller effluent concentration is achieved.     

无黏流体与固体界面的漏瑞利波

漏瑞利波存在于半无限无黏性流体和半无限固体媒质的界面处. 首先推导流固无限各向同性介质界面处漏瑞利波的特征方程和位移及应力的解析计算公式. 然后结合典型结构通过数值计算研究了漏瑞利波特性以及位移和应力在流体和固体中的分布规律. 数值计算结果表明漏瑞利波的相速度和衰减随流固密度比的增大而增大, 在流固界面上法向位移连续而切向位移不连续. 流固密度比对固体媒质中沿垂直于漏瑞利波的传播方向的位移、正应力和剪应力有比较大的影响,而对沿漏瑞利波的传播方向的正应力几乎没影响. 为利用漏瑞利波的无损检测与评价提供了理论基础.      

Research on the flow stability in a cylindrical particle two-phase boundary layer

IntroductionThecylindricalparticletwo_phaseflowsareofparticularinterestintheprocessingofcompositematerials ,textileindustry ,papermaking ,chemicalengineering ,foodprocessing[1].Thecylindricalparticlesinaflowcanmakethereinforcementofmaterials,thechangeofphysicalpropertyformaterialsandthereductionofdrag .Arranaga[2 ]reportedthatdragreductioneffectsareupto 60 %inpipeflowsbyaddingcylindricalparticlestoflow .Thecylindricalparticleshavealsoeffectsonthemechanismsofflowstability .Theeffectofcylindric…     

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.     

紊动流场中悬浮颗粒分布的随机理论

通过分析固体颗粒在紊动流场中的随机运动,建立了二维流场中垂直于时均流动的方向上颗粒随机位移的概率密度分布函数所满足的方程。由该方程解出的分布函数在一定条件下即相当于颗粒浓度分布函数。运用这一方法研究了[1]、[2]中报道的壁面附近颗粒浓度降低的现象。     

Characterization of the Mixing Layer Resulting from the Detonation of Heterogeneous Explosive Charges

A dense, two-phase numerical methodology is used to study the mixing layer developing behind the detonation of a heterogeneous explosive charge, i.e., a charge comprising of a high explosive with metal particles. The filtered Navier–Stokes equations are solved in addition to a sub-grid kinetic energy equation, along with a recently developed Eulerian–Lagrangian formulation to handle dense flow-fields. The mixing layer resulting from the post-detonation phase of the explosion of a nitromethane charge consisting of inert steel particles is of interest in this study. Significant mixing and turbulence effects are observed in the mixing layer, and the rms of the radial velocity component is found to be about 25% higher than that of the azimuthal and zenith velocity components due to the flow being primarily radial. The mean concentration profiles are self-similar in shape at different times, based on a scaling procedure used in the past for a homogeneous explosive charge. The peak rms of concentration profiles are 23–30% in intensity and decrease in magnitude with time. The behavior of concentration gradients in the mixing layer is investigated, and stretching along the radial direction is observed to decrease the concentration gradients along the azimuth and zenith directions faster than the radial direction. The mixing and turbulence effects in the mixing layer subsequent to the detonation of the heterogeneous explosive charge are superior to that of a homogeneous explosive charge containing the same amount of the high explosive, exemplifying the role played by the particles in perturbing the flow-field. The non-linear growth of the mixing layer width starts early for the heterogeneous explosive charge, and the rate is reduced during the implosion phase in comparison with the homogeneous charge. The turbulence intensities in the mixing layer for the heterogeneous explosive charge are found to be nearly independent of the particle size for two different sizes considered in the initial charge. Overall, this study has provided some useful insights on the mixing layer characteristics subsequent to the detonation of heterogeneous explosives, and has also demonstrated the efficacy of the dense, multiphase formulation for such applications.     

高浓度固-液两相流紊流的动理学模型

采用分子动理学方法,基于固-液两相流液相分子或颗粒相颗粒的Boltzmann方程,对Boltzmann方程分别取零矩和一次矩,则得到高浓度固-液两相流紊流的连续方程和动量方程,再和较成熟的低浓度两相流连续方程和动量方程比较,取低浓度两相流控制方程中较成熟合理的有关项和高浓度时由动理学方法推导出的颗粒间碰撞项,则得到高浓度固-液两相流紊流的最终控制方程:连续方程和动量方程.     

Numerische Berechnung des Partikeltransportes zu einer laminar angeströmten Kreisscheibe

Increasingly process steps become important, in which particles as product particles or contaminants are deposited on substrates out of the gas phase. In this paper the particles transport processes are investigated close to the surface of a circular plate surrounded by a laminar flow. The analogy between the governing equations of momentum, energy and mass is applied to the extended diffusion equation. In the nondimensional form the results of the numerical calculations give informations about velocity, temperature and particle concentration boundary layer thickness as well as their distributions. Especially the impact of external forces on particle concentration boundary layer thickness and profile is discussed. The transport of submicron particles to the surface due to convection, diffusion, gravity and thermophoretic forces acting independently is investigated. In the used normalized form the different forces are acting as one resulting force independently of their origin. Their resulting effect in comparison to the effect due to convective diffusive transport is important for particle deposition.     

饱和砂土中泥浆渗透的变形-渗流-扩散耦合计算模型

传统的泥浆渗透计算中没有考虑土体变形和浆液流速的影响.根据泥浆颗粒的质量守恒定律推导了耦合流速的浓度扩散方程,并通过在浓度方程中引入沉积系数进一步计算得到沉积颗粒的质量;同时,以沉积量作为耦合项对毕奥固结方程中的水量连续方程进行了修正,在此基础上建立了变形-渗流-扩散耦合的控制方程及其变分原理. 采用有限单元法求解基本方程,运用了时间增量法与直接迭代法,并利用一维试验验证计算方法的可靠性,并与赫齐格的经典模型的计算结果进行了比较,结果表明,本文建立的模型的计算结果可以较好地预测各组试验中颗粒的沉积规律,且吻合程度优于仅考虑颗粒对流和扩散的传统计算方法. 最后,将泥浆在槽壁中的渗透简化为二维问题并进行了计算,计算结果与工程认识相符合,泥浆的沉积填充效应随深度的增加而增大,施工时需要严格控制浅层作业段的机械垂直度;成槽机的下斗抓挖时机可以根据地层填充的致密程度进行计算,对现场施工具有一定的指导意义.      

Numerical study of concentrated fluid–particle suspension flow in a wavy channel

The particle migration effects and fluid–particle interactions occurring in the flow of highly concentrated fluid–particle suspension in a spatially modulated channel have been investigated numerically using a finite volume method. The mathematical model is based on the momentum and continuity equations for the suspension flow and a constitutive equation accounting for the effects of shear‐induced particle migration in concentrated suspensions. The model couples a Newtonian stress/shear rate relationship with a shear‐induced migration model of the suspended particles in which the local effective viscosity is dependent on the local volume fraction of solids. The numerical procedure employs finite volume method and the formulation is based on diffuse‐flux model. Semi‐implicit method for pressure linked equations has been used to solve the resulting governing equations along with appropriate boundary conditions. The numerical results are validated with the analytical expressions for concentrated suspension flow in a plane channel. The results demonstrate strong particle migration towards the centre of the channel and an increasing blunting of velocity profiles with increase in initial particle concentration. In the case of a stenosed channel, the particle concentration is lowest at the site of maximum constriction, whereas a strong accumulation of particles is observed in the recirculation zone downstream of the stenosis. The numerical procedure applied to investigate the effects of concentrated suspension flow in a wavy passage shows that the solid particles migrate from regions of high shear rate to low shear rate with low velocities and this phenomenon is strongly influenced by Reynolds numbers and initial particle concentration. Copyright © 2008 John Wiley & Sons, Ltd.     

Fully resolved simulations of viscoelastic suspensions by an efficient immersed boundary-lattice Boltzmann method

An efficient immersed boundary-lattice Boltzmann method (IB-LBM) is proposed for fully resolved simulations of suspended solid particles in viscoelastic flows. Stress LBM based on Giesekus and Oldroyd-B constitutive equation are used to model the viscoelastic stress tensor. A boundary thickening-based direct forcing IB method is adopted to solve the particle–fluid interactions with high accuracy for non-slip boundary conditions. A universal law is proposed to determine the diffusivity constant in a viscoelastic LBM model to balance the numerical accuracy and stability over a wide range of computational parameters. An asynchronous calculation strategy is adopted to further improve the computing efficiency. The method was firstly applicated to the simulation of sedimentation of a single particle and a pair of particles after good validations in cases of the flow past a fixed cylinder and particle migration in a Couette flow against FEM and FVM methods. The determination of the asynchronous calculation strategy and the effect of viscoelastic stress distribution on the settling behaviors of one and two particles are revealed. Subsequently, 504 particles settling in a closed cavity was simulated and the phenomenon that the viscoelastic stress stabilizing the Rayleigh–Taylor instabilities was observed. At last, simulations of a dense flow involving 11001 particles, the largest number of particles to date, were performed to investigate the instability behavior induced by elastic effect under hydrodynamic interactions in a viscoelastic fluid. The elasticity-induced ordering of the particle structures and fluid bubble structures in this dense flow is revealed for the first time. These simulations demonstrate the capability and prospects of the present method for aid in understanding the complex behaviors of viscoelastic particle suspensions.     

Rayleigh-Bénard convection heat transfer in nanoparticle suspensions

Natural convection heat transfer of nanofluids in horizontal enclosures heated from below is investigated theoretically. The main idea upon which the present work is based is that nanofluids behave more like a single-phase fluid rather than like a conventional solid-liquid mixture, which implies that all the convective heat transfer correlations available for single-phase flows can be extended to nanoparticle suspensions, provided that the thermophysical properties appearing in them are the nanofluid effective properties calculated at the reference temperature. In this connection, two empirical equations, based on a wide variety of experimental data reported in the literature, are developed for the evaluation of the nanofluid effective thermal conductivity and dynamic viscosity, whereas the other effective properties are evaluated by the traditional mixing theory. The heat transfer enhancement that derives from the dispersion of nano-sized solid particles into the base liquid is calculated for different operating conditions, nanoparticle diameters, and combinations of solid and liquid phases. One of the fundamental results is the existence of an optimal particle loading for maximum heat transfer across the bottom-heated enclosure. In particular, for any assigned combination of suspended nanoparticles and base liquid, it is found that the optimal volume fraction increases as the nanofluid average temperature increases, and may either increase or decrease with increasing the nanoparticle size according as the flow is laminar or turbulent. Moreover, the optimal volume fraction has a peak at a definite value of the Rayleigh number of the base fluid, that depends on both the average temperature of the nanofluid and the diameter of the suspended nanoparticles.     

Numerical simulations of dispersive mixing of viscoelastic suspensions in a four-roll mill

The dispersive mixing of particles suspended in Newtonian and viscoelastic fluids in a four-roll mill is studied by direct numerical simulations. A fictitious domain method is used to handle the particle motion. To quantify the mixing, a proper mixing distribution function is defined. The combined effect of fluid rheology and particle-particle/particle-wall hydrodynamic interactions is addressed. At variance with the Newtonian case where the particle distribution remains uniform, the viscoelasticity-induced migration leads to a significant segregation process. The effect of the Deborah number (the product of the fluid relaxation time and the roll angular velocity), shear-thinning, particle concentration, and size on the microstructure evolution is thoroughly investigated.     

Influences of the fluidizing and spouting pulsation on particle motion in spout-fluid beds

Effects of variable airflow on particle motion in spout-fluid beds are studied. Computational fluid dynamics using Navier–Stokes equations for the gas phase coupled with the discrete element method using Newton’s laws for the solid phase have been employed. Results indicate that increasing the fluidizing velocity diminishes dead zones and increases both the total height of the bed and the traversed distance by particles in the steady spout-fluid bed. In pulsed airflows, two configurations are investigated, namely, the spouted pulsed-fluidized bed with pulsed flow of the fluidizing velocity, and the pulsed-spouted fluidized bed with pulsed flow of the spouting velocity. The positive effect of pulsation on particle motion is shown and the effects of parameters, such as amplitude and frequency, on the dynamics of the bed are investigated in each configuration. An increase of up to 19% in traversed distance is found for the range studied, which suggests flow pulsation as a promising technique for increasing particle mixing in spout-fluid beds.     

Steady two-phase flow in a vertical tube subject to combined free and forced convection

Using the two-velocity, two-temperature model of a continuous medium, the viscousgravitational flow of a mixture of incompressible liquid and solid particles in a vertical round tube is considered. The free-convection equations are written down on the basis of the general equation of motion and the energy equation of a two-phase medium [1, 2]. Using a finite Hankel integral transformation, a solution is constructed for the case of a linear wall-temperature distribution along the tube. The results of some practical calculations of the velocity and temperature fields over the cross section of the tube are presented, together with the dimensionless heat-transfer coefficient expressed as a function of the Rayleigh number and phase concentration. Here it is assumed that the dynamic and thermal-interaction coefficients between the phases correspond to the Stokes mode of flow for each particle, as a result of which the velocity and thermal phase lag is very small [3].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 132–136, July–August, 1975.     

Dispersion characteristics of wave propagation in layered piezoelectric structures: Exact and simplified models

This article presents a study of the dispersion characteristics of wave propagation in layered piezoelectric structures under plane strain and open-loop conditions. The exact dispersion relation is first determined based on an electro-elastodynamic analysis. The dispersion equation is complicated and can be solved only by numerical methods. Since the piezoelectric layer is very thin and can be modeled as an electro-elastic film, a simplified model of the piezoelectric layer reduces this complex problem to a non-trivial solution of a series of quadratic equations of wave numbers. The model is simple, yet captures the main phenomena of wave propagation. This model determines the dispersion curves of PZT4-Aluminum layered structures and identifies the two lowest modes of waves: the generalized longitudinal mode and the generalized Rayleigh mode. The model is validated by comparing with exact solutions, indicating that the results are accurate when the thickness of the layer is smaller or comparable to the typical wavelength. The effect of the piezoelectricity is examined, showing a significant influence on the generalized longitudinal wave but a very limited effect on the generalized Rayleigh wave. Typical examples are provided to illustrate the wave modes and the effects of layer thickness in the simplified model and the effects of the material combinations.