基于摩擦-动力应力模型模拟喷动床气固流动

建立同时考虑颗粒和颗粒之间瞬间碰撞作用(瞬间接触)和非流化下颗粒滑动和滚动运动产生的颗粒半接触作用的颗粒碰撞-摩擦应力模型。推导具有普适性的固相动量守恒方程,建立高颗粒浓度气固两相流动模型,对喷动床内气体颗粒流动过程进行了数值模拟。模拟计算预测喷动床内喷射区、环隙区和喷泉区颗粒流动特性。预测喷动床内颗粒浓度和速度分布与他人实验结果相吻合。     

喷动床内纳米颗粒气固两相流动的数值模拟

采用欧拉-欧拉双流体模型,颗粒动理学理论模拟颗粒相流动,采用周涛和李洪钟(1999)的力平衡模型预测纳米颗粒聚团尺寸,对喷动床内纳米颗粒聚团流化过程进行了数值模拟,得到了喷动床内纳米颗粒聚团的流化过程,获得喷射区和环隙区内颗粒相速度和浓度分布。分析了喷动床结构和进口气体速度等对纳米颗粒聚团流化特性的影响。由于纳米颗粒的特殊性质,不易形成喷泉区。适当的喷动床结构和进口气体速度有助于形成稳定喷动。     

喷动流化床的数值模拟研究

本文借鉴 Ｃｕｎｄａｌｌ＆ Ｓｔｒａｃｋ提出的碰撞模型,将欧拉方法和拉格朗日法结合起来,建立了喷动流化床气固两相流动的数学模型,在计算机上对试验台进行了模拟,探讨了床内颗粒的运动特征。按１：１的比例建立了喷动流化床的模型试验台并进行了试验,对模拟结果和试验结果进行了比较与分析。颗粒在喷动床内的运动特征,模拟结果都能够再现出来,模拟结果与试验结果取得了很好的一致性。     

数值模拟颗粒旋转对流化床内气固两相流动特性的影响

流化床内颗粒自旋转将影响颗粒相的流动特性.本文运用基于颗粒动理学理论的欧拉-欧拉气固多相流模型,考虑颗粒自旋转流动对颗粒碰撞能量交换和耗散的影响,数值模拟流化床内气体颗粒两相流动特性.计算结果表明颗粒的自旋转使得床内更容易形成气泡,颗粒浓度分布变化增大.颗粒自旋转运动将导致床内非均匀结构更明显.     

湍动流化床内气固两相流动特性的数值模拟

采用欧拉-欧拉双流体模型,颗粒动理学方法模拟颗粒脉动流动和κ-ε双方程模型模拟气相湍流流动,考虑气固两相间耦合作用,数值模拟湍动流化床内气固两相流动行为,获得颗粒浓度和颗粒速度分布.计算结果表明湍动流化床呈现下部密相区、上部稀相区的颗粒分布特性.在密相区,沿床径向方向颗粒浓度在床中心处低、壁面逐渐增高;在稀相区颗粒浓度分布较均匀.沿轴向方向颗粒浓度呈底部浓度高、顶部浓度低的"S"型分布.     

颗粒数方程的求解及流化床数值模拟

流化床内颗粒聚并和破碎将影响颗粒相的流动特性.本文运用基于颗粒动理学理论的欧拉一欧拉气固多相流模型,利用直接矩积分方法求解颗粒数平衡方程,建立颗粒数密度与连续性方程、动量方程之间的关系,数值模拟流化床内两种不同直径颗粒发生聚并时气固两相流动特性。计算结果表明,颗粒聚并伴随着床内颗粒直径逐渐增大,床内颗粒流化状态逐渐变为固定床状态,两种颗粒直径均增加,且小颗粒的体积分数逐渐减小、大颗粒的体积分数增加。当仅考虑聚并过程时增加流化速度将导致床内颗粒体积平均直径变大。随着颗粒密度减小,床内体积平均直径增加。     

过热蒸汽/热载体/生物质三相混合数值模拟

建立了工业尺度数值试验的过热蒸汽/热载体/生物质三相欧拉数理模型。水蒸气设为主连续相,采用k-ε方程描述;热载体和生物质设为两种不同的拟流体相,其黏度和压力采用颗粒动理论方程进行封闭;气固相互作用力采用Symbol-O'Brein曳力模型描述,固固之间作用力采用morsi-alexander模型描述;气固之间传热采用Gunn模型描述。基于该模型的数值试验,实现了喷动流化床混合器的优化设计使热载体与生物质的混合程度达70%以上;并研究了物料性质对生物质终温的影响,生物质终温受生物质粒径的影响最为显著。     

循环流化床内颗粒团运动的数值模拟

本文引入聚合力来表征稠密气固两相流动中颗粒所受到的团聚效应,井将颗粒团视为床内的离散相,建立了颗粒团的碰撞、形成及破碎模型。用此模型对循环流化床的两相流动进行数值模拟,得到了较详细的床内流动特性。计算结果表明,用本文的模型和算法模拟工程意义上的循环流化床内的两相流动是可行的。     

喷动床内颗粒浓度的实验测量和理论模拟

颗粒浓度对于稠相流动体系中颗粒应力本构方程建立具有决定性的作用。本文采用双流体模型,并结合颗粒流动力学理论来描述固相应力的基础上实现了对喷动床的数值模拟,探讨了从喷射区到环隙区颗粒浓度的理沦变化趋势。同时采用自主开发的电容探针技术,不仅颗粒浓度测量值和理论值吻合较好,而且也准确地验证了颗粒非均匀夹带这一现象。电容探针还可灵敏地分辨出环隙区内不同床高处以及接近壁面处的颗粒浓度的细微变化。     

四角喷燃炉膛内湍流三维气相燃烧和气固两相流动的数值模拟

目前,已建立的绝大多数炉内三维数值模拟都以颗粒轨道模型为基础。颗粒相连续介质模型已成功地应用于模拟气固两相射流及煤粉燃烧室内的三维气固两相流动。本文将它应用于四角喷燃炉膛内流动的模拟,作为炉内粉煤燃烧全过程数值模拟的第一步。     

静止水中单个上升气泡的直接数值模拟

本文发展了基于Front Tracking的直接数值模拟方法研究气液两相界面的迁移特性,该方法对气液两相采用半隐式的分步法直接求解N-S方程,耦合Front Tracking Method获得两相界面的三维变形。针对无边界以及垂直壁面附近静止水中的单个气泡上升过程进行模拟,研究气泡运动的机理以及气泡与壁面的相互作用。数值模拟准确再现了气泡的上升过程和变形,不同Re数下气泡的上升速度计算结果同经验关联式非常吻合,验证了该方法的有效性。随后分析了气泡周围流场的结构,发现壁面对气泡周围流场的抑制是壁面对气泡作用力的主要原因,将导致气泡逐渐偏离垂直壁面。     

Hydrodynamic binary coalescence of droplets under air flow in a hydrophobic microchannel

Based on the volume of fluid(VOF) method, we conduct a numerical simulation to study the hydrodynamic binary coalescence of droplets under air flow in a hydrophobic rectangular microchannel. Two distinct regimes, coalescence followed by sliding motion and that followed by detaching motion, are identified and discussed. Additionally, the detailed hydrodynamic information behind the binary coalescence is provided, based on which a dynamic mechanical analysis is conducted to reveal the hydrodynamic mechanisms underlying these two regimes. The simulation results indicate that the sliding motion of droplets is driven by the drag force and restrained by the adhesion force induced by the interfacial tension along the main flow direction. The detachment(i.e., upward motion) of the droplet is driven by the lift force associated with an aerodynamic lifting pressure difference imposed on the coalescent droplet, and also restrained by the adhesion force perpendicular to the main flow direction. Especially, the lift force is mainly induced by an aerodynamic lifting pressure difference imposed on the coalescent droplet. Two typical regimes can be quantitatively recognized by a regime diagram depending on Re and We. The higher Re and We respectively lead to relatively larger lift forces and smaller adhesion forces acting on the droplet, both of which are helpful to detachment of the coalesced droplet.     

LES-DSMC方法研究超细颗粒气固两相流动过程

采用考虑颗粒脉动流动对气相湍流流动影响的大涡模拟(LES)研究气相湍流,采用直接模拟蒙特卡罗方法(DSMC)模拟颗粒间的碰撞。单颗粒运动满足牛顿第二定律,颗粒相和气相相间作用的双向耦合由牛顿第三定律确定,考虑超细颗粒间的van der Waals作用力。数值模拟垂直管内超细颗粒气固两相流动,对颗粒相速度、浓度以及团聚物流动过程进行分析。     

Rotational intermittency and turbulence induced lift experienced by large particles in a turbulent flow

The motion of a large, neutrally buoyant, particle freely advected by a turbulent flow is determined experimentally. We demonstrate that both the translational and angular accelerations exhibit very wide probability distributions, a manifestation of intermittency. The orientation of the angular velocity with respect to the trajectory, as well as the translational acceleration conditioned on the spinning velocity, provides evidence of a lift force acting on the particle.     

On the diffusion of a fast molecule

We consider the motion of a heavy particle in interaction with an infinite ideal gas of slow atoms. We prove that the velocity of the heavy particle is, in a suitable limit, modeled by a deterministic process. We also treat the process of rescaled velocity fluctuations around a certain deterministic motion and show that this is appropriately modeled by a nonhomogeneous diffusion process.Supported in part by NSF grants PHY-8201708 and DMR81-14726.Supported in part by NSF grant PHY-8003298 and the Seton Hall University, Research and Faculty Development Council.     

Inertial microfluidics with multi-particle collision dynamics

Using the method of multi-particle collision dynamics (MPCD), we investigate inertial focussing in microfluidic channels that gives rise to the Segré-Silberberg effect. At intermediate Reynolds numbers, we model the motion of a spherical colloid in a circular microchannel under pressure-driven flow. We determine the radial distribution function and show how its width and the location of its maximum are strongly influenced by the colloid size and the Reynolds number of the Poiseuille flow. We demonstrate that MPCD is well suited for calculating mean values for the lift force acting on the colloid in the cross-sectional plane and for its mean axial velocity. We introduce a Langevin equation for the cross-sectional motion whose steady state is the Boltzmann distribution that contains the integrated lift force as potential energy. It perfectly coincides with the simulated radial distribution function.     

Finite element simulation of nanostructuring burnishing

Dynamics simulation of burnishing of a thin steel layer beneath an indenter applied with a constant force and then moving with a constant velocity was performed by the finite element method in the plane strain approximation. The indenter was modeled by a perfectly rigid body, and the steel was modeled by an elastoplastic body with isotropic hardening according to an experimentally defined law. The regularities of changes in the stressstrain state of the material near the treated surface were studied and mechanisms of the formation of a nanostructured layer were disclosed. The effect of the friction coefficient and the burnishing force on the height of a bulge of edged material was analyzed. The results of studies agree well with experimental data.     

Drag induced lift in granular media

Laboratory experiments and numerical simulation reveal that a submerged intruder dragged horizontally at a constant velocity within a granular medium experiences a lift force whose sign and magnitude depend on the intruder shape. Comparing the stress on a flat plate at varied inclination angle with the local surface stress on the intruders at regions with the same orientation demonstrates that intruder lift forces are well approximated as the sum of contributions from flat-plate elements. The plate stress is deduced from the force balance on the flowing media near the plate.     

气温和颗粒密度对声场中颗粒动力学影响的数值模拟

综合考虑黏性夹带力、Basset力、虚拟质量力和压力梯度力,建立颗粒在声场中的动力学模型,利用变步长四阶RungeKutta算法和二阶隐式Adams插值算法对颗粒的受力和运动进行数值模拟。将模拟和实验得到的颗粒运动特性进行对比,验证数值模拟的正确性。在此基础上,研究气温和颗粒密度对颗粒动力学的影响规律。结果表明,黏性夹带力对颗粒运动起主导作用;气温升高,压力梯度力与黏性夹带力之间的相位差减小,Basset力、虚拟质量力与黏性夹带力之间的相位差增大。研究还发现,气温较低时,颗粒密度对颗粒运动有重要影响,夹带系数随着密度的增加而迅速下降;气温较高时,颗粒密度对颗粒运动的影响较小,颗粒位移振幅和夹带系数相对低温时明显增加。      

II. Recirculation flow induced by a bubble stream rising in a viscous liquid

The recirculation flow induced by the rising motion of a bubble stream in a viscous fluid within an open-top rectangular enclosure is studied. The three-dimensional volume averaged conservation equations are solved by a control-volume method using a hybrid finite differencing scheme to describe the liquid phase hydrodynamics. The momentum exhange between the bubbles and the liquid phase is modeled with a source term equals to the volumetric buoyancy force acting on the gas in the bubble stream. The volumetric buoyancy force accounts for in line interactions between bubbles through the average gas volume fraction in the gas liquid column which depends on the size and the rising velocity of bubbles. The fluid flow within an open-top rectangular enclosure is further investigated by particle image velocimetry for a bubble stream rising in a water-glycerol solution. The measured fluid velocities in a vertical plane are compared with the predictions of the numerical model over a wide range of fluid viscosity (43 mPa s-800 mPa s) and gas flow rates. Finally, the recirculation flows resulting from the interaction of two neighbouring vertical bubble streams are studied. Received: 23 July 1997 / Revised: 19 December 1997 / Accepted: 11 May 1998