介观尺度流体绕流球体的耗散粒子动力学模拟

采用耗散粒子动力学(dissipative particle dynamics, DPD)方法, 对两平行平板间流体绕流三维球体进行了计算. 球体和平行平板由达到平衡状态的冻结DPD粒子组成, 流体在不同无量纲外力驱动下流动, 球体受力由组成球体的所有冻结DPD粒子求和得到. 流动达到充分发展后, 输出球体在流动方向的受力, 并计算球体的阻力系数, 与文献中的关联式进行了对比. 结果表明, 在Re≤qslant 100的范围内, DPD方法能较准确地计算出阻力系数, 在较大雷诺数时, 由于流     

耗散粒子动力学处理复杂固体壁面的一种有效方法

耗散粒子动力学(dissipative particle dynamics,DPD)作为一种介观尺度拉格朗日型粒子方法,已经成功地应用于微纳米流动和生化科技的研究中.复杂固体壁面的处理和壁面边界条件的实施一直是DPD方法发展及应用的一个障碍.提出了处理复杂固体壁面的一种新的方法.复杂固体区域通过冻结随机分布并且达到平衡状态的DPD粒子代表;所冻结的DPD粒子位于临近流动区域的一个截距内;在靠近固体壁面的流动区域中设置流动反弹层,当流动DPD粒子进入此流动层后反弹回流动区域.应用这种固体壁面处理方法对简单流动区域的Poiseuille流动和复杂多孔介质内的流动进行了分析.研究表明,这种新的固体壁面处理方法能够有效模拟复杂固体区域,准确实施壁面边界条件.     

微液滴动力学特性的耗散粒子动力学模拟

对传统的耗散粒子动力学方法进行了改进.改进的耗散粒子动力学方法采用了包含远程吸引力和近距排斥力的保守力势函数,从而使得用耗散粒子动力学方法模拟多相流动成为可能.应用改进的耗散粒子动力学方法,对微尺度下液滴的形成及液滴在微重力下的大幅度振荡变形进行了数值模拟.计算结果表明,改进的耗散粒子动力学(DPD)方法能够有效地描述微尺度下液滴的动力学特性,对研究复杂流体多相流动有着重要的意义. 关键词： 多相流 微液滴 耗散粒子动力学(DPD)方法 保守力势函数     

简单剪切流条件下液滴破碎和碰撞模拟

利用耗散粒子动力学(DPD)对简单剪切流条件下液滴的破碎和碰撞分离过程进行了模拟.液滴由多颗简单的DPD粒子构成,液滴与周围流体的不互溶性由提高它们之间的保守力系数获得.通道的上下壁面沿着相反方向运动从而使流体产生简单剪切流动.计算结果表明,在简单剪切流条件下,当界面张力数超过临界界面张力数后,液滴不稳定;液滴将被拉长并破碎成数个较小的稳定的液滴.当两个液滴相互碰撞并分离的现象发生时,两个液滴的质心在速度梯度方向上的距离与整个碰撞过程相关,碰撞之后的δz值比碰撞前的δz大.     

微通道疏水表面滑移的耗散粒子动力学研究

了解疏水表面的滑移规律对其在流动减阻方面的应用至关重要.利用耗散粒子动力学(dissipative particle dynamics, DPD)方法研究了微通道疏水表面的滑移现象.采用固定住的粒子并配合修正的向前反弹机制,构建了DPD固体壁面边界模型,利用该边界模型模拟了平板间的Couette流动.研究结果表明,通过调整壁面与流体间排斥作用强度,壁面能实现从无滑移到滑移的转变,壁面与流体间排斥作用越强,即疏水性越强,壁面滑移越明显,并且滑移长度与接触角之间存在近似的二次函数关系.无滑移时壁面附近密度分布均匀,有滑移时壁面附近存在低密度区域,低密度区域阻碍了动量传递,致使壁面产生滑移.     

适用复杂几何壁面的耗散粒子动力学边界条件

耗散粒子动力学(DPD)是一种针对介观流体的高效的粒子模拟方法,经过二十多年发展已经在诸如聚合物、红细胞、液滴浸润性等方面有了很多研究应用.但是因为其边界处理手段的不完善,耗散粒子动力学模拟仍局限于相对简单的几何边界问题中.本文提出一种能自适应各种复杂几何边界的处理方法,并能同时满足三大边界要求:流体粒子不穿透壁面、边界处速度无滑移、边界处密度和温度波动小.具体地,通过给每个壁面粒子赋予一个新的矢量属性—局部壁面法向量,该属性通过加权计算周围壁面粒子的位置得到;然后通过定义周围固体占比概念,仅提取固体壁面的表层粒子参与模拟计算,减少了模拟中无效的粒子;最后在运行中,实时计算每个流体粒子周围固体粒子占比,判断是否进入固体壁面内,如果进入则修正速度和位置.我们将这种方法应用于Poiseuille流动,验证了该方法符合各项要求,随后还在复杂血管网络和结构化固体壁面上展示了该边界处理方法的应用.这种方法使得DPD模拟不再局限于简单函数描述的壁面曲线,而是可以直接从各种设计图纸和实验扫描影像中提取壁面,极大地拓展了DPD的应用范围.     

耗散粒子动力学中施加无滑移边界条件的新方法

在耗散粒子动力学(DPD)中施加无滑移边界条件是较困难的,本文提出了一种新的处理方法来实现无滑移边界条件.通过选取合适的壁面与流体粒子之间的保守力系数可以控制近壁面区的粒子数,从而可以使密度波动降低至极小值,实现无滑移边界条件.通过该方法模拟泊肃叶流动求得的流体性质与通过使用Lees-Edwards边界条件模拟剪切流动所求得的流体性质一致,证明了该方法的可靠性.     

基于晶格-Boltzmann方法的纳米流体流动和传热模型

纳米流体是由流体与纳米粒子组成的悬浮体.悬浮在流体中的纳米粒子会受到运动阻力、布朗力、扩散力和重力等作用力的影响,因而其运动规律极其复杂.本文运用晶格-Boltzmann(LB)方法,建立纳米流体流动和传热模型模型,对纳米流体的流型结构和温度场进行了模拟和分析.     

小通道扁管内纳米流体流动与传热特性

建立了测量小通道扁管内纳米流体流动与对流换热性能的实验系统,测量了不同粒子体积份额的水-Cu纳米 流体的管内对流换热系数和摩擦阻力系数,实验结果表明,在相同雷诺数条件下,小通道扁管内纳米流体的对流换热系数 大于纯液体,且随粒子的体积份额的增加而增大,而纳米流体的阻力系数并未明显增大。     

SPH方法进出口数值边界条件研究

基于虚粒子概念,提出进出口数值边界条件处理方法.在进口边界外设置进口区域,并赋予进口虚粒子相应的物理量;出口边界外设置出口和缓冲区域,出口虚粒子的物理量由计算获得,缓冲区的虚粒子初始物理量是指定的;根据每个时间步的流动情况,改变粒子的区域属性,添加/删除相应的粒子.利用拉格朗日形式的SPH方法,通过管内流动问题验证数值进出口边界方法的适用性,研究进出口边界条件在激波管、绕流问题中的应用.提出的进出口边界处理方法,避免了边界附近流体粒子积分截断问题,保证流体粒子能够流出边界,激波能够透射边界.     

Drag and Lift Forces Acting on a Sphere in Shear Flow of Power-Law Fluid

Laminar flow of a power-law fluid over a sphere is considered for unbounded shear flow. The Navier–Stokes equations with power-law viscosity are solved numerically using an in-house developed CFD package. Vorticities structures downstream of particle are suppressed for powerlaw fluid. The shear rate influence on drag force is negligible for power index close to unit, and the drag force appreciably decreases with falling power index. For small Reynolds numbers, the lift force coefficient monotonically decreases against the power index and exhibits an opposite behavior for moderate values of Reynolds numbers. The results of the parametric studies are used to derive correlations for the drag force and to detect the hydrodynamic differences from uniform flow. The investigation parameters varied within the following ranges: power-law index 0.3 ≤ n ≤ 1, Reynolds number 0 < Re ≤ 150, and dimensionless shear rate 0.05 ≤ s ≤ 0.4.     

A Particle Resistance Model for Flow through Porous Media

A particle model for resistance of flow in isotropic porous media is developed based on the fractal geometry theory and on the drag force flowing around sphere. The proposed model is expressed as a function of porosity, fluid property, particle size, fluid velocity （or Reynolds number） and fractal characters D f of particles in porous media. The model predictions are in good agreement with the experimental data. The validity of the proposed model is thus verified.     

Calculations of the flow resistance and heat emission of a sphere in the laminar and high-turbulent gas flows

An early drag crisis can occur at high turbulence of incoming gas flow to a sphere. To study the influence of a crisis on heat transfer from a sphere to gas, a numerical experiment was carried out in which the free gas flow around a sphere with a temperature lower than the sphere temperature was simulated for two cases. The flow was laminar in the first case and highly turbulent in the second case. To take into account turbulence, the kinematic coefficient of turbulent viscosity with a value, which is much higher (up to 2000 times) than that for physical viscosity, was introduced. The results of calculations show that the early drag crisis occurs at Reynolds numbers of about 100 and results in considerable (by four to seven times) decrease in the hydrodynamic force and sphere drag coefficient C _d . The early drag crisis is also accompanied by the crisis of heat transfer from a sphere to gas with a decrease in Nusselt numbers Nu by three to six times.     

Experimental verification of the early crisis of drag using a single sphere

The recently discovered phenomenon of drag anomaly (early crisis) observed for droplets in a two-phase turbulent flow at transition Reynolds numbers (Re = 10–100) is experimentally supported in the case of a gas flow about a single hard sphere as well. Using a specially designed torsion balance, the drag force of a foam plastic sphere 3 mm in diameter placed in a blower-produced turbulent jet is measured. It is found that the turbulence of the flow about the sphere must be high for the drag anomaly to arise.     

Calculation of the flow about a sphere and the drag of the sphere under laminar and strongly turbulent conditions

To investigate the influence of a strongly turbulent incoming flow on the hydrodynamic drag of a body and occurrence of the early crisis of drag, a numerical experiment is conducted in which a free gas flow about a sphere is simulated for two cases, namely, for a laminar flow and for a strongly turbulent flow. Turbulence is simulated by assuming a high kinematic coefficient of turbulent viscosity. Calculation data lead us to conclude that the early crisis of drag at Reynolds numbers near 100, which shows up as a considerable (four-to sevenfold) decrease in the hydrodynamic force and the drag coefficient of the body, can be explained by the strong turbulence of the incoming flow.     

Mean force on a small sphere in a sound field in a viscous fluid

A mean force exerted on a small rigid sphere by a sound wave in a viscous fluid is calculated. The force is expressed as a sum of drag force coming from the external steady flow existing in the absence of the sphere and contributions that are cross products of velocity and velocity derivatives of the incident field. Because of the drag force and an acoustic streaming generated near the sphere, the mean force does not coincide with the acoustic radiation pressure, i.e., the mean momentum flux carried by the sound field through any surface enclosing the sphere. If the sphere radius R is considerably smaller than the viscous wave penetration depth delta, the drag force can give the leading-order contribution (in powers of delta/R) to the mean force and the latter can then be directed against the radiation pressure. In another limit, delta< or =R, the drag force and acoustic streaming play a minor role, and the mean force reduces to the radiation pressure, which can be expressed through source strengths of the scattered sound field. The effect of viscosity can then be significant only if the incident wave is locally plane traveling.     

Capillary force on colloidal particles in a freely suspended liquid thin film

The dynamic behavior of micron-sized polystyrene latex particles confined in a free-standing liquid film is experimentally studied. When the thickness of the film is less than the particle diameter and varies depending on position, the particles are accelerated toward the thicker region. Using a simple geometrical model and hydrodynamic theory, we calculate the capillary force on the particle. The drag coefficient of the particle is found to depend on the thickness of the film from the value of near zero to the Stokes' drag coefficient.     

Calculation of the drag and heat transfer from a sphere in the gas flow in a cylindrical channel

A numerical experiment on the simulation of heat transfer from a sphere to a gas flow in a cylindrical channel in the Stokes and transient flow regimes has been described. Radial and axial profiles of the gas temperature and the dependences of drag coefficient C_d of the body and Nusselt number Nu on Reynolds number Re have been calculated and analyzed. The problem of the influence of the early drag crisis for a sphere on its heat transfer to the gas flow has been considered. The estimation of this phenomenon has shown that the early drag crisis of the sphere in a strongly turbulent flow causes a reduction in heat transfer from the sphere to the gas by three to six times (in approximately the same proportion as for its drag coefficient).