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

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

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

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

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

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

TVD格式在超音速喷管三维粘性流动求解中的应用

详细给出了任意三维曲线坐标系中Ｎｏｖｉｅｒ－Ｓｔｏｋｅｓ方程的对流项ＴＶＤ格式的构造过程，建立了数值求解三维粘性流动的计算方法，应用该方法对三维超音速喷管中有激波及无激波情况下的两种工况的层流流场进行了数值求解，并与实验做了对比。结果表明本文建立的计算方法具有较高的精度，同时也证明ＴＶＤ格式具有分辩率高，稳定收敛等优点，为进一步开展叶栅流场及紊流的研究打下了基础。     

双颗粒在溶解条件下沉降的多相流动特性

对牛顿流体内溶解双颗粒在垂直管道中的沉降运动进行了直接数值模拟. 流体运动由守恒方程计算, 密度和黏性的变化考虑流场温度变化的影响, 通过积分粘性应力和压力获得颗粒的受力跟踪颗粒运动, 溶解引起的相变及其形状的变化由溶解潜热、溶解质量与分散相边界处的温度梯度的关系建立的方程决定. 通过颗粒和流体间相互的作用力和力矩及边界条件的施加实现相间耦合. 对双颗粒在等温流体无溶解条件和非等温流体溶解条件下的沉降过程进行了计算. 结果表明, 在一定雷诺数内, 热对流产生的颗粒尾迹处涡的脱落以及溶解引起的颗粒质量、颗粒表面形态的变化引起了颗粒的横向摆动, 并使颗粒沉降速度发生了变化.     

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

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

直接模拟中不同边界条件的实施及对沉降规律的影响

在牛顿流体中, 对颗粒在4种不同边界的垂直通道中的沉降运动进行了直接数值模拟. 计算结果表明:通过计算区域随颗粒运动而移动构建的无限长通道能准确模拟颗粒自由下落到稳定沉降的发展过程; 周期性边界条件由于流场变化, 对颗粒沉降产生了影响, 不能模拟颗粒的自由沉降过程; 底部封闭边界适合模拟封闭容器内颗粒与固壁的相互作用过程, 若颗粒达到稳定沉降, 也能模拟无限长通道内的沉降过程; 流化边界适合模拟流化床内气固两相流动. 计算结果有助于更好地理解和使用不同边界条件. 关键词： 直接数值模拟 边界条件 沉降 任意拉格朗日-欧拉方法     

溃坝问题的粒子方法模拟

溃坝问题包括水体自由表面变形、翻卷、破碎、复杂的紊流和漩涡,涉及非线性复杂水动力学问题,含有丰富而复杂的流动机理。本文改进了传统的物质点方法,将物质点法从模拟固体问题拓展到了流体问题。采用改进的光滑粒子流体动力学和物质点方法,对两种长宽比的水坝坍塌过程进行了数值模拟,得到了波浪前沿及自由表面顶面随时间的演化过程。计算结果表明,两种粒子方法模拟的水体发展过程与实验结果吻合较好,较VOF更接近实验值。对涌波与竖直壁面的作用进行了计算,结果表明,SPH方法与其它方法相比能更好地捕捉水腔形态。最后分析了两种粒子方法在处理流体大变形问题时的特点以及计算的准确性和精度。