超疏水小球低速入水空泡研究

物体入水问题是一类复杂的流固耦合问题,具有广泛的工程应用背景.物体在跨越自由液面入水的过程中,在一定的条件下,会向水中卷入空气形成空泡,空泡的运动还可能形成指向物体的射流,从而对物体的受力及其运动过程产生影响.超疏水表面能够在物体入水过程中形成多尺度流固耦合作用,进而影响物体的运动和宏观流动现象.而对于小尺度的小球低速入水问题,表面和界面力往往起主导作用.为了在更广的参数空间获得超疏水小球入水空泡类型和小球的运动特性,采用高速摄影实验方法,研究了半径0.175～10 mm的超疏水小球低速入水及空泡动力学行为,获得了小球漂浮振荡、准静态空泡、浅闭合空泡、深闭合空泡和表面闭合空泡5种类型的动力学行为,探讨了这些运动行为与韦伯数We和邦德数Bo之间的关系,并推导了小球漂浮振荡与下沉现象的无量纲关系.研究结果表明:超疏水小球的入水及空泡动力学行为主要与韦伯数We和邦德数Bo有关.在邦德数Bo O (10~(-1))范围内,表面张力对流动的影响显著,随着韦伯数We的增大,小球入水及空泡动力学行为依次经历漂浮振荡、准静态闭合、浅闭合、深闭合和表面闭合;在邦德数O(10~(-1))Bo O(1)范围内,漂浮振荡现象不再发生;当邦德数Bo O(1)后,浅闭合现象也不再发生;小球漂浮振荡与下沉现象的临界关系可以用相似律关系描述.     

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

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

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

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

基于半解析VOF-DEM的激光直接沉积多尺度过程模拟

与传统铸造技术相比, 基于金属粉末的增材制造技术因其生产周期短、可操作性强而在航空航天、生物医学等领域具有很好的优越性. 尤其是激光直接沉积技术, 因其自由度高, 在复杂构件制造、部件修复中有着广泛的运用. 但是该激光直接沉积过程涉及多物理场、跨尺度、极端高温高压环境和相变问题, 仅靠实验不能很好地研究其中的机理. 已有数值模拟技术一般通过预设或者射入拉格朗日点作为颗粒输入, 不能做到同时考虑环境气体、颗粒碰撞和相变过程. 本文在近期发展的基于核函数近似背景流场的半解析CFD-DEM耦合方法中引入了流体体积分数法(VOF), 发展了可以同时模拟含热、刚体颗粒、相变和自由液面及相变界面的半解析VOF-DEM (或半解析CFD-DEM-VOF)方法, 从而首次实现了真实物理环境下激光直接沉积技术的数值模拟. 其中, VOF中的气相为环境气体, 液相为熔融和凝固的金属相, 界面通过iso-Advector重构, DEM为未熔化的金属粉末, 且流体网格可解析离散元颗粒形状. 这一模拟框架可以有效复现颗粒之间的碰撞、粘结、熔化、融合, 以及熔池熔道的形成, 为激光直接沉积技术的数值模拟提供了开拓性的范式, 并可以被应用到其他带相变的颗粒系统中.      

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

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

光滑粒子动力学方法中粒子分布与数值稳定性分析

光滑粒子动力学(SPH)作为一种拉格朗日型无网格粒子方法,已经成功地应用于包括含多相流动界面以及移动边界的可压缩和不可压缩流体运动的研究中.通过对Poiseuille流动的深入研究,探索了SPH方法中粒子分布对计算精度的影响,揭示了一种因为粒子不规则分布而导致的数值不稳定现象.研究显示,这种数值不稳定性起源于SPH方法粒子近似过程中的不连续性.使用了一种新的粒子近似格式以确保SPH方法中粒子近似的连续性.计算结果表明,这种新的粒子近似格式对于规则和不规则的粒子分布都能得到稳定精度的结果.     

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

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

光滑粒子动力学方法的发展与应用

光滑粒子动力学(smoothed particle hydrodynamics,SPH)是一种拉格朗日型无网格粒子方法,已经成功地应用到了工程和科学的众多领域.SPH使用粒子离散及代表所模拟的介质,并且基于粒子体系估算和近似介质运动的控制方程.本文分析和综述了SPH模拟方法的发展历程、数值方法与应用进展.介绍了SPH方法的基本思想;从连续性、边界处理、稳定性和计算效率4个方面阐述了SPH方法的研究现状;介绍了SPH方法近年来在可压缩流动、不可压缩流动以及弹塑性材料高速变形与失效方面的一些典型应用;并对SPH方法的发展与应用进行了预测与展望.      

介观尺度通道内多相流动的耗散粒子动力学模拟

采用四次方光滑函数构造了包含远程吸引近距排斥的保守力势函数, 应用该势函数, 采用耗散粒子动力学方法, 对十字型微通道内的多相流动进行了计算. 计算结果表明, 该势函数能有效模拟十字型微通道内的流动过程及模式. 关键词： 多相流 耗散粒子动力学 保守力势函数     

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

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