等直径微液滴碰撞过程的改进光滑粒子动力学模拟

运用一种改进光滑粒子动力学(SPH)方法模拟了相溶和不相溶两种情况下的等直径微液滴碰撞动力学过程. 为提高传统SPH方法的数值精度和稳定性, 采用一种不涉及核导数计算的核梯度改进形式; 为处理液滴界面张力采用修正的van der Waals表面张力模型. 通过模拟牛顿液滴碰撞聚并变形过程并与相关文献或试验结果进行对比, 验证了改进SPH 方法模拟微液滴碰撞过程的可靠性. 随后, 研究了基于van der Waals模型相溶聚合物微液滴碰撞聚并变形过程及不相溶微液滴碰撞后的反弹、分离过程, 讨论了碰撞过程中碰撞速度、碰撞角度、密度比等参数对碰撞变形过程的影响, 分析了流体桥、旋转角度等因素的变化情况. 关键词： 光滑粒子动力学 微液滴 聚合物液滴 碰撞     

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

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

黏性液滴变形过程的核梯度修正光滑粒子动力学模拟

本文提出了一种核梯度改进光滑粒子动力学(KGC-SPH)方法,模拟了黏性液滴形变自由表面问题.首先,通过模拟等温黏性液滴拉伸和旋转变形,验证了KGC-SPH法较SPH法具有较高精度和更好稳定性,且能很好地保持总角动量守恒.其次,基于非等温van der Waals模型对平衡态圆形液滴的形成过程进行数值研究,观察到小幅度振荡现象,并给出了一种新的克服张力不稳定性的方法和一种适合KGC-SPH方法的新的表面张力处理技术.最后,研究了van der Waals液滴的周期性振荡现象,讨论了初始椭圆形液滴长短半轴比 关键词： 光滑粒子动力学 黏性液滴 van der Waals模型 表面张力     

基于光滑粒子动力学方法的液滴冲击固壁面问题数值模拟

采用改进的光滑粒子动力学(SPH)方法对液滴冲击固壁面问题进行了数值模拟. 为了提高传统SPH方法的计算精度和数值稳定性, 在传统的SPH方法的基础上对粒子方法中的密度和核梯度进行了修正, 采用了考虑黎曼解法的SPH流体控制方程, 构造了一种新型的粒子间相互作用力(IIF)模型来模拟表面张力的影响. 应用改进的SPH方法对液滴冲击固壁面问题进行了数值模拟. 计算结果表明：新型的IIF 模型能够较好地模拟表面张力的影响, 改进的SPH方法能够精细地描述液滴与固壁面相互作用过程中液滴的内部压力场演变和自由面形态变化, 液滴的铺展因子随初始韦伯数的增大而增大, 数值模拟结果与实验得到的结果基本一致. 关键词： 液滴 固壁面 光滑粒子动力学 表面张力     

液滴溅落问题的光滑粒子动力学模拟

对传统的光滑粒子动力学方法进行了改进, 改进的光滑粒子动力学方法对传统粒子方法中的核近似式和粒子近似式进行了修正, 采用Riemann 算法求解光滑粒子动力学流体控制方程, 添加了表面张力的计算程序, 考虑了表面张力对液滴溅落的影响. 应用改进的光滑粒子动力学方法对液滴静止状态下冲击液面的飞溅过程进行了数值模拟. 计算结果表明, 改进的光滑粒子动力学方法能够有效地描述液滴溅落液面的动力学特性和自由表面变化特征, 能够得到稳定精度的结果.     

液滴冲击液膜问题的光滑粒子动力学模拟

本文对传统的光滑粒子动力学方法进行了改进．改进的光滑粒子动力学方法对传统粒子方法中的核梯度进行了修正,采用了一种新型的耦合边界条件,添加了表面张力和人工应力的计算程序．应用改进的光滑粒子动力学方法对液滴冲击液膜问题进行了数值模拟．得到了不同时刻液滴内部的压力变化特征,精细地捕捉了不同时刻的自由面,从机理上分析了液滴产生飞溅的条件,探讨了韦伯数,表面张力对液滴冲击液膜问题的影响．计算结果表明,改进光滑粒子动力学方法能够有效地描述液滴冲击液膜的动力学特性和自由表面变化特征,能够得到稳定精度的结果．     

二维管道充填过程的修正SPH模拟

通过施加一种密度初始化方法对传统光滑粒子动力学(SPH)方法进行修正,提出一种修正SPH方法.同时,为了提高边界上数值计算的准确性,提出一种新的固壁边界处理方法.通过修正SPH方法模拟液滴拉伸问题和溃坝问题,验证修正SPH方法的准确性和可靠性.随后,对研究很少的管道充填过程进行修正SPH模拟,并讨论Re对流场及涡的影响.数值结果表明,修正SPH方法能够准确模拟牛顿流体管道充填过程,且流动受Re的影响较大.     

基于MRT-LB伪势模型的液滴撞击液膜数值研究

运用改进的格子玻尔兹曼（LB）伪势多松弛多相模型,研究单/双液滴撞击液膜的流动特性.考察单液滴在不同气液相密度比时撞击液膜的发展.随着密度比的减小,冠状水花顶端开始向内弯曲,且底部半径显著减小.在大密度比情况下研究双液滴撞击液膜.结果表明：双液滴撞击液膜有中心射流的形成;液滴水平间距的增大,延缓中心射流的出现,并降低初始中心射流的高度;随Re数的增加,中心射流的高度明显增大.     

微通道内椭球颗粒惯性聚焦行为数值研究

基于浸没边界-格子Boltzmann方法,对方形截面微通道内椭球颗粒的惯性迁移与旋转动力学行为进行数值研究,发现微通道内椭球颗粒的惯性迁移存在两种主要的运动状态：①翻转状态,即椭球颗粒前进过程中长轴始终在中心对称平面内;②滚动状态,即椭球颗粒前进过程中长轴始终垂直于中心对称平面.研究表明：在低Re数（Re=10）下颗粒以两种状态随流体迁移至平衡位置;在较大Re数（50≤Re≤200）下最终均以翻转状态随流体迁移,随Re数增加,平衡位置先逼近壁面后远离壁面.通过对不同运动状态下椭球颗粒周围的微观流场进行分析,提示该微观流动在颗粒惯性聚焦行为特征中有重要影响,并从流体和颗粒的惯性角度对颗粒不同运动状态的转换机理给出解释.     

eXtended Pom-Pom黏弹性流体的改进光滑粒子动力学模拟

黏弹性流体广泛存在于自然界和工业生产中,对其复杂流变特性的研究具有重要的学术价值和应用意义.本文提出一种改进的光滑粒子动力学方法,对基于eXtended Pom-Pom模型的黏弹性流动进行了数值模拟.为了提高计算精度,采用一种不含核导数计算的核梯度修正离散格式.为了防止粒子穿透固壁,提出一种增强型的边界处理技术.为了消除张力不稳定性,将人工应力耦合到动量守恒方程中.运用改进光滑粒子动力学方法数值模拟了基于eXtended Pom-Pom模型的黏弹性Poiseuille流和黏弹性液滴撞击固壁问题,通过与解析解或有限差分方法解的比较以及对数值收敛性的评价,验证了改进光滑粒子动力学方法的有效性和优势,并在此基础上,深入分析了Reyonlds数、Weissenberg数、溶剂黏度比、各向异性参数、松弛时间比和分子链臂数等流变参数对流动过程的影响.     

修正表面张力算法的SPH方法及其实现

在Morris提出的表面张力SPH方法基础上,通过引入CSPM方法对边界法向的计算和曲率的计算进行修正,得到表面张力修正方程组;通过半圆形算例测试方法和Morris方法在边界定位、法向计算和曲率计算等影响表面张力关键因素的求解精度,研究曲率计算中应采用的光滑长度值.模拟初始方形液滴在表面张力作用下的自然变化过程,并与M...     

非等温非牛顿黏性流体流动问题的修正光滑粒子动力学方法模拟

为准确、有效地模拟非等温非牛顿黏性流体的流动问题,本文基于一种不含核导数计算的核梯度修正格式和不可压缩条件给出了一种改进光滑粒子动力学(SPH)离散格式,它较传统SPH离散格式具有较高精度和较好稳定性.同时,为准确地描述温度场的演化过程,建立了非牛顿黏性的SPH温度离散模型.通过对等温Poiseuille流、喷射流和非等温Couette流、4:1收缩流进行模拟,并与其他数值结果作对比,分别验证了改进SPH方法模拟非牛顿黏性流动问题的可靠性和提出的SPH温度离散模型求解非等温流动问题的有效性和准确性.随后,运用改进SPH方法结合SPH温度离散模型对环形腔和C形腔内非等温非牛顿黏性流体的充模过程进行了试探性模拟研究,分析了数值模拟的收敛性,讨论了不同位置处热流参数对温度和流动的影响.     

SPH simulation of fuel drop impact on heated surfaces

The interaction of liquid drops and heated surfaces is of great importance in many applications. This paper describes a numerical method, based on smoothed particle hydrodynamics (SPH), for simulating n-heptane drop impact on a heated surface. The SPH method uses numerical Lagrangian particles, which obey the laws of fluid dynamics, to describe the fluid flows. By incorporating the Peng–Robinson equation of state, the present SPH method can directly simulate both the liquid and vapor phases and the phase change process between them. The numerical method was validated by two experiments on drop impact on heated surfaces at low impact velocities. The numerical method was then used to predict drop-wall interactions at various temperatures and velocities. The model was able to predict the different outcomes, such as rebound, spread, splash, breakup, and the Leidenfrost phenomenon, consistent with the physical understanding.     

重构核粒子法对光滑粒子法的改进效果

光滑粒子法通过核函数进行近似估计,在计算域边界附近,核估计的精度明显下降.重构核粒子法通过校正函数对核函数进行重新构造,提高核估计方法在边界点和内点上对函数的估计精度以及计算的稳定性.研究发现,虽然校正函数的构造立足于对函数的精确估计,但这个优势同样能在对导数的估计中继续保持.通过理论研究及数学、物理模型的模拟,展示重构核粒子法的改进效果,揭示其能够提高精度的原因.     

瑞利-泰勒不稳定问题的光滑粒子法模拟研究

提出了一种适用于模拟多相流的光滑粒子法,该方法对密度方程在交界面处的离散格式进行了修正以适应多相流所涉及的大密度比问题,在不同相粒子之间施加了很小的排斥力以防止粒子穿透交界面,并采用了最新发展的双曲型光滑函数以消除应力不稳定问题.应用该多相流光滑粒子法模拟研究了单模态和多模态瑞利-泰勒不稳定问题.通过与文献中结果的对比研究表明:在模拟瑞利-泰勒不稳定问题时,本文方法的结果明显优于文献中的大部分光滑粒子法模拟结果,与Grenier等(2009 J.Comput.Phys.228 8380)的结果相当,但本文方法比Grenier等的方法简单方便.对于单模态瑞利-泰勒不稳定问题,研究了交界面的形态,涡结构的演化过程以及贯穿深度随时间的变化关系.对于多模态瑞利-泰勒不稳定问题,研究了交界面演化过程中小尺度结构合并成大尺度结构的过程,水平方向的平均密度随高度的变化关系,以及贯穿深度随时间的变化关系.     

基于弱可压与不可压光滑粒子动力学方法的封闭方腔自然对流数值模拟及算法对比

为了对比研究弱可压光滑粒子动力学(WCSPH)方法和不可压光滑粒子动 力学(ISPH)方法在模拟封闭方腔自然对流问题时的特性, 采用粒子位移技术有效地解决了高瑞利数条件下, 拉格朗日型SPH方法模拟封闭方腔自然对流时流体域内的粒子聚集和空穴问题, 将拉格朗日型SPH 方法求解封闭方腔自然对流问题的最高瑞利数提高到了10⁶; 进而通过对比瑞利数分别为10⁴, 10⁵, 10⁶的条件下, 采用拉格朗日型WCSPH、 拉格朗日型ISPH、欧拉型ISPH三种SPH方法模拟得到的封闭方腔速度分布云图、 温度分布云图、壁面努赛尔特数分布曲线和平均努塞尔特数, 分析了三种SPH方法在模拟封闭方腔自然对流时的精度、稳定性和计算效率. 结果表明: 在低瑞利数条件下, 以上三种SPH方法都可以较好地模拟此问题, 在高瑞利数条件下, 欧拉型ISPH方法的模拟结果最为精确; 拉格朗日型WCSPH方法模拟所得结果比拉格朗日型ISPH方法模拟所得结果稍好些. 关键词： 光滑粒子动力学 不可压光滑粒子动力学 粒子位移技术 自然对流