自由表面流动的移动粒子半隐式模拟方法

引入可压缩因子，对移动粒子半隐式方法（MPS）进行了改进，解决了MPS方法在计算中容易出现的不稳定问题．使用改进后的方法对液柱倒塌算例进行了模拟计算，并将计算结果与实验及其他数值计算结果进行了比较．计算表明，引入可压缩因子大大提高了计算的稳定性，且方法简单易于实现，在液柱与挡板发生碰撞、流体发生大变形时依然能进行有效地模拟。     

一种适用于移动粒子半隐式法的自由表面识别方法

针对移动粒子半隐式 (MPS) 法中自由表面判定方法的不足,提出了一种基于邻居搜索的自由表面识别方法,其介于MPS法每时间步的显式和隐式修正之间,基于粒子周围的邻居分布,将周围区域分为8块,其中相邻区域无邻居时的粒子属于自由表面粒子,且在实际中,该方法与粒子数密度参数混合使用,在可能被误判的数密度参数区域使用邻居搜索.结果表明,该混合方法能够有效识别自由表面,避免非自由表面粒子被误判,且对算法的效率影响不大.     

弱可压缩移动粒子半隐式方法对溃坝溢流问题的数值模拟

无网格粒子法(拉格朗日法)是计算流体力学领域最新一代的计算方法。由于不需要生成网格,无网格粒子法能够模拟任意形状的界面及其变形和碎裂,尤其适用于计算如溃坝溢流及波浪破碎等具有复杂自由表面运动的流体力学问题。其中,光滑粒子流体动力学方法(SPH)和移动粒子半隐式法(MPS)是无网格计算方法的典型代表。采用基于MPS法的修正弱可压缩移动粒子半隐式法(weakly compressible corrected moving particle semi-implicit method)构建了用于求解自由表面流问题的数值计算模型,并用于模拟干床面与湿床面上的溃坝溢流问题。计算结果与实验结果的对比分析表明,采用基于WCCMPS无网格法建立的数值计算模型对干、湿床面上溃坝溢流的溢流形态、水柱高度以及波前速度等均给出了良好的计算结果。     

液滴斜向撞击液面的移动粒子半隐式法数值模拟

采用移动粒子半隐式法在三维情况下对液滴斜向撞击液面的过程进行了数值模拟,很好地解决了撞击过程中自由液面的大变形问题,捕捉了液滴变形、冲击凹坑和液柱凸起3个特征现象及其发展变形过程,分析了不同撞击角度对液柱凸起的高度、角度及运动前沿的影响.结果表明,液柱凸起最高点的运动在水平方向上基本匀速,液柱凸起与水平面的夹角随撞击角度减小而减小,当凸起达到最高点后,该夹角将在重力与黏性作用下骤减.     

移动粒子半隐式方法的计算效率优化

针对无网格法计算效率普遍不高的局限性,以移动粒子半隐式(MPS)方法为例,将计算耗时较多的程序模块作为优化对象,以降低计算等级和减少计算次数为优化手段,从搜索算法、数学运算和计算实现形式3个主要方面进行了改进.提出了用低级运算来代替高级运算、用方位比较来代替距离比较的位置比较消冗法来优化搜索算法,消除了大部分无效计算,大大降低了计算的复杂程度和计算量,使得该方法的计算效率得到显著提高,增强了MPS方法的应用范围和工程时效性.该算法具有非常好的普适性,容易在类似算法中移植,还有随着计算规模的增大优化效果更加明显的特点.同时,发现多种优化方式的组合使用并非效果叠加,需考虑它们的兼容关系.     

移动粒子半隐式法流固耦合模型及自由浮体数值研究

对流固耦合问题进行了系统的分类并在移动粒子半隐式法中建立了刚体与流体之间的单向及双向耦合模型。对二维溃坝问题和活塞式造波板的造波问题进行了数值模拟,结果与实验值或分析解符合得较好。对浮体在静水中的自由振荡问题进行了计算,浮体质心高度与理论值的最大偏差控制在0.5%以内,达到了工程计算的精度。在长板形自由浮体对正弦波浪的运动响应算例中,综合应用了单向和双向耦合模型,结果表明长度小于波长的浮体随波浪做正弦规律的起伏运动和摇摆运动,同时沿行波方向迂回前进,当浮体长度增加时,浮体运动幅度减小。     

用移动粒子半隐式法模拟液舱横摇晃荡现象

应用移动粒子半隐式法（MPS）模拟了半满粘性液舱的横摇晃荡运动.模拟中每个时间步分为显示和隐式两步,控制方程中利用粒子与周围粒子之间的关系取代了传统方法中的差分项,利用粒子数密度的变化来判别自由表面.MPS法作为一种Lagrange数值方法,存在着粒子穿透现象.文中使用了一种防穿透方法,取得了较好的结果.结果表明,MPS法在复杂曲面的布点及自由表面的表现上有一定的优势.     

一种适用于移动半隐式法的缓解压力振荡方法

针对无网格移动粒子半隐式（MPS）法压力振荡性过大,且难以处理压力的问题,提出一种新的压力缓解方法.该方法在MPS法压力计算方程中加入用以保证流体不可压缩的速度散度为零的条件,在MPS法计算完成后对压力进行适当的时间平均,得出最终计算压力,在一定程度上缓解压力的振荡.通过溃坝和船体横摇晃荡2个实例,将新方法和MPS法进行比较,验证新方法在不增加计算量的基础上对缓解压力振荡有很好的效果.     

基于移动粒子半隐式法的齿轮搅油损失分析与试验验证

将移动粒子半隐式法应用于齿轮搅油损失分析,以研究齿轮传动中搅油损失的内在特性.以单齿斜齿轮为对象分析其在不同浸油深度、齿轮转速和齿轮宽度条件下的搅油损失情况.分析结果表明,齿轮转速对搅油损失的影响最大;浸油深度对搅油损失的影响也较为显著;齿轮宽度对搅油损失的影响则不太明显.根据仿真分析结果建立了单齿斜齿轮搅油损失数学模型.搭建试验台架,通过试验与仿真数据的比较,验证了移动粒子半隐式法在计算搅油损失方面的可行性和准确性,从而为齿轮传动搅油损失的进一步研究提供了一种新方法.     

移动粒子半隐式法核函数特征对压力求解稳定性的影响

针对移动粒子半隐式法在求解特定问题时,压力求解可能会出现一定程度的波动,分析了移动粒子半隐式法中核函数曲线形状特征对压力求解稳定性的影响,构造了一种指数多项式型核函数。模拟了典型静压(静水压力问题)和动压(液体晃动问题)算例,并将模拟结果与理论解或实验值进行对比,研究结果表明:改进的核函数可有效抑制模拟过程中压力求解的振荡现象;核函数与对应粒子数密度比值曲线的形状特征可真实反映粒子间相互作用关系,在稳定性分析中起着至关重要的作用。当核函数是光滑单调递减非负函数且最大值为有限值、两粒子间距离与影响半径的比r/r_e在[0,1]区间时,曲线两端附近核函数数值变化平缓更有利于使粒子保持合理距离,压力求解更加稳定;在r/r_e为0.8的附近,核函数值过小时会影响系统的动力学性能。     

单液膜气泡聚并和连接过程的移动粒子半隐式法数值模拟

为了研究单液膜气泡单层液膜包裹气体的特殊结构形式和内外壁面均受气液表面张力的特殊力学形式,在拉格朗日框架下采用无网格移动粒子半隐式法并基于表面自由能表面张力模型,建立了单液膜双气液界面表面张力模型,从而实现了单液膜气泡振荡变形过程中的复杂界面计算和捕捉。在此基础上对2个单液膜气泡的聚并和连接过程进行了模拟分析,获得了典型的流动现象和液膜变形特征与规律,发现减小表面张力系数或增大黏性系数均会减弱气泡变形过程中表面张力项的变形主导作用。为此,提出了凹点切线法用于计算连接型气泡的液膜夹角,明晰了连接型气泡的形状。计算结果可为工业消泡技术提供一定的理论依据。     

RDSPH:一种适用于一维非连续条件的新SPH方法

基于恢复粒子一致性的光滑粒子流体动力学(RSPH)方法,同时借鉴DCSMP方法,在非连续区间上对未知函数分段泰勒展开,并保留零阶项和一阶项,构建了一套适合于模拟一维非连续物理现象的新SPH方法(RDSPH).比较不同SPH方程近似二次函数的结果,新方法不仅改进了传统SPH方法存在的边界缺陷问题,同时在非连续区域内也能更有效地修复由非连续引起的核近似截断积分和消除粒子不一致性所造成的误差.     

自由界面问题的拉格朗日粒子和流体体积耦合算法

针对自由界面问题,构建了拉格朗日粒子和流体体积(volume of fluid,VOF)耦合算法.拉格朗日粒子方法可以准确追踪运动界面,但是一般很难保证流体的质量守恒性.VOF方法可以保证很好的质量守恒性,但是不容易计算界面的几何信息.因此,本文构造了一种耦合算法,吸收两种方法的优点.耦合算法中还引入了四叉树自适应网格...     

二维自由表面流动的光滑粒子动力学方法数值模拟

在流体微可压缩的假定下,运用光滑粒子动力学(SPH)方法模拟了二维自由表面流动问题.选用正则化五次样条函数作为SPH核函数,由状态方程求得压力,采用显式三步法进行时间积分,利用链表法作为粒子搜索策略,并引入动态粒子模拟边界条件,提出了多项修正措施以弥补SPH方法的缺点.研究表明,文中计算结果与实验数据吻合,计算精度理想,且消除了流体粒子附着在边界上的不合理现象.     

Buoy Relay Method for Instantaneous Fluid Flow with Free Surface

ＩｎｔｒｏｄｕｃｔｉｏｎＦｒｅｅｂｏｕｎｄａｒｉｅｓａｒｅｓｕｒｆａｃｅｓｏｎｗｈｉｃｈｄｉｓｃｏｎｔｉｎｕｉｔｉｅｓｅｘｉｓｔｉｎｏｎｅｏｒｍｏｒｅｖａｒｉａｂｌｅｓ.Ｅｘａｍｐｌｅｓａｒｅｆｒｅｅｓｕｒｆａｃｅｓ(ｂｕｂｂｌｅｓｏｒｄｒｏｐｌｅｔｓ),ｍａｔｅｒｉａｌｉｎｔｅｒｆａｃｅｓｏｒｉｎｔｅｒｆａｃｅｓｂｅｔｗｅｅｎｆｌｕｉｄｓａｎｄｆｌｅｘｉｂｌｅｓｔｒｕｃｔｕｒｅｓ.Ｔｈｒｅｅｔｙｐｅｓｏｆｐｒｏｂｌｅｍｓａｒｉｓｅｉｎｔｈｅｎｕｍｅｒｉｃａｌｔｒｅａｔｍｅｎｔｏｆｆｒｅｅｂｏｕｎｄ…