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

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

微液滴振荡过程的光滑粒子动力学方法数值模拟

本文对传统的光滑粒子动力学方法进行了改进, 改进的光滑粒子动力学方法对传统粒子方法中的核梯度进行了修正, 采用了一种新型的核函数和离散格式, 应用改进的光滑粒子动力学方法对微液滴振荡过程进行了数值研究. 研究了不同纵横比和雷诺数(Re)下振荡阻尼与振荡的周期、振幅与Re数的关系. 研究表明: 对于纵横比λ≤ 4时的微液滴振荡过程, 其他参数恒定不变的前提下, Re数越大, 液滴形状变化越剧烈, 波的阻尼作用越弱, 液滴振荡周期变长; 在Re数一定的前提下, 随着液滴初始的纵横比的增大, 液滴振动的振幅增大, 液滴振荡的周期变长.     

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

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

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

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

液滴在气固交界面变形移动问题的光滑粒子流体动力学模拟

为准确模拟液滴在气固交界面变形移动问题,对基于连续表面张力模型的表面张力光滑粒子流体动力学方法进行了改进.改进方法采用新的边界处理方式和界面法向修正方法,即将固体边界虚粒子色函数值根据液面的位置进行相应设定以保证气-液-固三相交界处流体粒子的界面法向沿接触线法线方向,引入Brackbill提出的壁面附着力边界条件处理方法,对在气-液-固三相交界处的流体粒子及部分固体边界虚粒子的界面法向进行修正,修正前后保持法向模值不变,得到了含壁面附着力边界条件的表面张力算法.模拟了受壁面附着力影响的水槽中液面的变化过程、液滴润湿壁面过程和剪切气流驱动液滴在固体表面变形脱落过程,并与流体体积函数方法进行了对比.结果表明,该方法在处理壁面附着力问题时精度较高,稳定性较好,适合处理工程中液滴在气固交界面变形移动问题.     

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

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

基于范德瓦尔斯表面张力模式液滴撞击疏水壁面过程的研究

液滴撞击疏水壁面过程的研究在介观流体力学和微流体作用材料科学的研究中具有重要的理论意义和工程价值. 论文在SPH方法中引入范德瓦尔斯状态方程处理液滴表面张力, 考虑流体粒子之间远程吸引, 近程排斥的内部作用力, 提出了流体粒子与疏水壁面粒子间势能函数与表面张力相结合的作用模式. 通过模拟真空条件下两个静止的等体积液滴相互融合的过程, 验证了计算模式在模拟液滴的表面张力中的有效性. 采用该模式模拟的液滴撞击疏水壁面过程, 不仅能够有效地模拟液滴撞击壁面后的变形过程, 而且清晰地模拟出液滴的回弹、腾空以及二次撞壁现象的完整过程. 模拟结果与液滴撞击疏水壁面的实验结果以及VOF模拟结果符合较好, 表明本文所提出的表面张力和疏水壁面作用力处理模式对模拟液滴撞壁过程具有实际应用价值.     

基于不可压缩光滑粒子动力学的黏性液滴变形过程仿真

应用基于投影算法的不可压缩光滑粒子动力学(incompressible smoothed particle hydrodynamics, ISPH)法对黏性液滴变形过程进行了数值仿真. 对于张力失稳导致的粒子非物理簇集问题, 采用粒子移位技术加以解决. 为了验证本文ISPH 算法的精度和稳定性, 分别模拟了圆形黏性液滴的拉伸变形过程以及方形液滴的旋转变形过程, 得到了不同时刻液滴内部的压力变化特征, 准确地捕捉了液滴自由面演化过程, 并将数值计算结果与文献中的解析解进行了比较.分析结果表明, 基于投影算法的不可压缩光滑粒子动力学方法结合粒子移位技术, 能够有效地模拟黏性液滴变形过程, 可以得到精确和稳定的结果. 关键词： 不可压缩光滑粒子动力学 黏性液滴 自由面流动 数值仿真     

一种新型光滑粒子动力学固壁边界施加模型

由于Lagrange粒子法的本质, 固壁边界条件的施加一直是光滑粒子动力学方法的难点之一. 本文从固壁边界的物理原理出发, 应用多层虚粒子表征固壁边界, 提出了一种新型固壁边界施加模型. 将虚粒子看作流体的扩展, 计算中虚粒子密度保持不变, 压力、速度等参数通过对流体粒子的插值获得, 虚粒子有条件的参与控制方程的计算, 对流体的密度/压力产生影响, 通过压力梯度隐式地表征壁面与流体之间的作用强度并对流体粒子施加沿壁面法线方向的斥力作用, 防止流体粒子对壁面的穿透. 数值算例测试结果表明, 与现有固壁边界施加方法相比, 本文方法更加符合流体与固壁边界作用的物理原理, 可以简单、有效地施加固壁边界条件, 方便地应用于具有复杂几何边界的问题, 获得稳定的流场形态、规则的粒子秩序及良好的速度、压力等参量的分布.     

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

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

自由上浮气泡运动特性的光滑粒子流体动力学模拟

基于虚功原理, 在Hu X Y等和Grenier N等的研究结果基础上推导了多相流光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)控制方程, 采用精度较高的黏性力和表面张力模型, 发展了一套适用于具有大密度比和大黏性比界面的多相流SPH方法. 首先, 通过施加人工位移修正, 适当背景压力和异相界面力, 使得计算全程粒子分布相对均匀, 改善了界面处的失稳现象, 防止了异相界面处粒子的非物理性穿透; 在此基础上, 利用方形流体团振荡模型对表面张力模型进行了验证, 数值结果与解析解甚为吻合; 然后采用上浮气泡经典数值算例对比研究了不同黏性力计算方法、不同核函数的适用性以及人工位移修正的效果; 最后, 对单个气泡的上浮、变形、撕裂以及垂向两个气泡的追赶、融合等现象进行了模拟, 初步揭示了气泡上浮过程中各种有趣物理现象的细节过程和动力学机理.     

光滑粒子流体动力学方法固壁处理的一种新型排斥力模型

与传统网格法相比, 光滑粒子流体动力学方法不能直接施加壁面边界条件, 这就限制了该方法在工程中的应用.为此, 本文基于Galerkin加权余量法并结合传统排斥力方法, 推导出一种新的排斥力公式来施加壁面边界条件.该方法不含未知参数, 能在不减小边界粒子尺寸的情形下有效地防止流体粒子穿透壁面, 同时可避免邻近边界的流体粒子的速度及压力振荡. 分别通过静止液柱算例、液柱坍塌算例、容器中液体静止算例及溃坝算 例来验证本文方法的有效性, 并与传统边界处理方法进行对比, 结果表明: 本文方法克服了传统方法存在的缺陷, 是一种有效的固壁边界处理方法. 关键词： 光滑粒子流体动力学法 固壁边界 排斥力 加权余量法     

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

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