基于SPH方法的瞬态粘弹性流体的数值模拟

运用SPH(Smoothed Particle Hydrodynamics)方法模拟基于Oldroyd-B模型的平面突然起动Couette流,通过数值解与解析解的比较,验证SPH方法模拟瞬态粘弹性流动的准确性;且对基于Oldroyd-B模型的方腔驱动流进行SPH模拟.采用一种新的固壁边界处理方法,有效地防止了粒子穿透,提高数值计算的准确性.用数值算例验证SPH方法对粘弹性流体模拟的有效性和稳定性.     

耦合方法在超高速碰撞数值模拟中的应用

 相比于传统的数值方法，无网格光滑粒子法（Smoothed Particle Hydrodynamics，SPH）更适合模拟超高速碰撞问题，它能够清晰地模拟材料界面，并且能够克服中网格大变形导致的数值不稳定。然而在SPH方法中，为了得到t时刻一个节点上的物理量，需要先找出与它相互作用的临近节点，导致计算速度很慢，计算所需的时间远远大于有限元方法所需的时间。采用SPH与Lagrange法相耦合的方法对3D超高速碰撞过程进行了模拟研究。计算结果表明：SPH与Lagrange法耦合方法的计算结果与SPH计算结果精度接近，并且能够大大地节省计算时间。     

圆锥边界附近激光空泡溃灭行为的研究

为了研究刚性圆锥边界锥角对激光空泡溃灭行为的影响,文章建立了虚拟平面边界模型, 同时采用阴影摄影术、光偏转法以及数值计算的手段对边界附近空泡溃灭过程进行了研究. 结果表明边界的锥角对空泡的形状、溃灭时间以及液体射流形成均有明显影响. 空泡形状偏离球形的程度和溃灭时间均随锥角的增大而增大,且增大锥角度可以促使射流的形成. 空泡溃灭时间的实验值同理论值具有较高的一致性,验证了虚拟平面边界假设及无量纲距离修正的有效性.     

激光空泡在文丘里管中运动的动力学特性

以水为工作介质,在不同文丘里管入口压力下,利用YAG激光器产生的激光轰击水中的金属靶材产生空泡,借助高速摄像系统记录激光空泡在文丘里管中的运动过程,并采用流体动力学模拟对文丘里管中的流场特性和空泡的溃灭特性进行分析.结果表明:激光空泡在文丘里管中的运动,其形状变化可分为产生阶段、挤压阶段、溃灭初始阶段和溃灭阶段等四个阶段.空泡的溃灭取决于流场状态,当流动为层流时,空泡不发生溃灭;当流动为湍流时,空泡发生溃灭,且湍流程度越剧烈,溃灭现象越显著.搭建的激光空泡生成与运动系统能够实现空泡的定点溃灭.     

聚焦换能器中心轴线上泡群的平移运动

空化中空泡群的速度计算具有显著的理论和应用意义。本文基于试验方法,在自行搭建的试验台上利用高速观察摄影设备获得空泡群平移运动的连续图像记录后,利用PIV(Particle Image Velocitry)中相关性的原理,对空泡群的平移运动进行了分析,得到了换能器中心轴线上空泡运动初步的量化值。在此基础上,对空泡受到声辐射压力作用下的加速运动进行了理论推导,很好地解释了计算中得到的速度结果。     

SPH方法对金属表面微射流的数值模拟

 光滑粒子流体动力学（Smoothed Particle Hydrodynamics）方法是一种无网格的拉氏计算方法,具有计算格式简单、易于计算大变形问题等优势。用二维SPH程序模拟了冲击加载下金属铝表面的沟槽状缺陷产生微射流的过程,得到了喷射物的总质量、最大速度和质量-速度曲线。还计算了相同深度、不同夹角的沟槽微射流,并与实验结果进行了比较。结果表明：计算得到的喷射物总质量和最大速度与实验结果符合较好,SPH方法对于估算金属表面的微射流提供了一种有效的数值模拟手段。     

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

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

完全变光滑长度SPH法及其实现

提出完全变光滑长度SPH法及其算法实现.方程组基于对称形式核函数近似,SPH密度演化方程与变光滑长度方程隐式关联;在Springel提出的全守恒SPH方程组基础上,通过将分散核近似形式改进为对称核近似形式得到SPH动量方程和能量方程.采用迭代求解密度演化方程和变光滑长度方程,显式求解SPH动量方程和能量方程,增加的计算量相对很少.给出三个lD激波管算例和2D Sedov算例验证方法的有效性.数值结果表明,算法保持动量和能量的守恒律,解决了传统SPH法中由于变光滑长度影响带来的计算误差,且在模拟2D Sedov问题时能得到比Springel方法更准确的压强峰值位置和中心压强值.特别适合于模拟爆炸与冲击、大变形大扭曲等密度梯度和光滑长度梯度剧烈变化的问题.     

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

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

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

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