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

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

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

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

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

使用移动粒子半隐式(MPS)法对典型的自由表面流动问题--二维溃坝问题进行了模拟,研究了不同粒子分布对模拟结果的影响,发现对于过大和过小的粒子间距均无法得到合理的结果.将模拟结果与SPH、VOF法结果和实验结果进行了比较,发现MPS的模拟结果介于VOF和SPH的模拟结果之间.在没有引入黏性的情况下,MPS法得到的结果比VOF法的结果更接近于实验结果,预期引入黏性后可达到与SPH法同样合理的结果.     

WC-CMPS法对溃坝溢流问题的数值模拟

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

表面自由能模型的改进及液滴浸润数值研究

通过优化极值点位置,对表面自由能模型进行了研究和改进,改进后的模型解决了原模型中表面粒子聚集和液滴形状畸变的问题.采用移动粒子半隐式法(MPS)结合改进后的表面自由能模型,对受表面张力作用的自由面运动问题进行了研究和分析.研究了包括气液、固液和不互溶液液界面在内的表面张力收缩和液滴浸润过程,并提取了接触角及浸润面积等特征量进行分析.数值计算结果与理论分析结果吻合较好,从结果中可以明显观察到固液浸润过程的3个发展阶段,浸润速度随时间递减且与浸润接触角成反比.研究结果验证了MPS方法结合改进后的表面自由能模型可以更加有效、准确地模拟多界面张力的作用.     

负压梯度下自由面流动模拟的无网格粒子法Dirichlet边界条件

为改善无网格粒子法自由表面附近计算精度及稳定性,尤其在负压梯度条件下较难收敛的状况,提出了一种基于虚拟边界的狄利克雷边界条件赋值方法。以移动粒子半隐式法为例,采用最小二乘法对自由表面粒子进行多项式拟合,以拟合曲线作为虚拟边界,用自由表面粒子到虚拟边界的偏差矢量对自由表面粒子赋值进行修正。通过方形液滴旋转算例验证了算法的有效性与准确性,获得了自由表面平滑的模拟结果。进一步研究了关键参数的选择对边界条件的表达和计算精度的影响,发现当拟合搜索半径为4倍初始粒子间距时,可减小拟合误差,防止边界畸变。提出的改进方法提高了自由表面核截断区域的计算精度与稳定性,可显著改善负压梯度下大变形自由曲面外缘粒子因梯度力偏差而导致的粒子飘散情况,为改善大变形复杂界面流动计算提供了新的思路。     

用无网格粒子法直接模拟多相多组分界面流

模拟了以下两个多相多组分含复杂界面变形的流动:在脉冲加速下的液滴破裂和突然减压下的闪蒸射流。采用了移动粒子半隐式法(MPS)直接数值模拟方法。MPS法是一种适用于不可压缩流体流动计算的无网格粒子法。结果表明液滴破裂的临界Weber数等于13,该值与Pilch和Erdman的实验值相吻合。对于闪蒸射流,射流由于表面闪蒸被剥离为锥状,且射流长度随着压力衰减比的增加而减少。这些特征与实验观察是一致的。这两种模拟的结果表明了MPS法在处理复杂界面流动计算中的优势以及其鲁棒性。这些结果对反应堆事故分析有意义。     

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

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

一种改进的RBF神经网络混合学习算法

提出一种基于粒子群优化算法、K-means算法及减聚类算法的径向基函数(RBF)神经网络混合学习算法.该算法使用减聚类方法确定隐层节点数,具有自适应确定隐层节点的能力,避免了调整隐层节点的人为干预.通过K-means算法形成粒子群优化(PSO)算法初始粒子群,避免了初始粒子群的随机性,提高了粒子群优化算法的优选能力;采用PSO算法训练RBF神经网络中的所有参数.数值结果表明,改进的混合算法具有更高的分类准确率。     

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

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

MPS方法在三维溃坝问题中的应用

MPS(Moving Particle Semi-Implicit)方法是一种拉格朗日粒子法,在处理大变形流动问题中具有很大的优势.将MPS方法应用到三维溃坝流动问题中,并对粒子的移动方式进行了修改,改进后的方法(XMPS)可以使粒子移动更加有序,避免粒子间的相互穿透.数值计算结果表明:MPS方法(包括XMPS)在处理复杂的自由面流动问题中具有很好的灵活性和可靠性,数值模拟能够比较准确地预测出自由面的形状和位置,即使在水面出现翻卷、入水和水柱撞击到障碍物时,数值结果仍然能够与实验相吻合.从细节上看,XMPS可以给出更为清晰、光滑的自由面形状,并能更好地捕捉到流动的细节特征.     

颗粒群组方法中关键参数的选取

为了深入了解颗粒群组方法的特点,提出了修正最大群组数方法,研究了颗粒群组方法中的两个关键参数(体积增量和最大群组数)对数量密度和体积分数的影响.结果表明:传统最大群组数在计算颗粒数量密度时比较准确,修正最大群组数在计算颗粒体积分数时比较准确.随着颗粒体积增量的增加,传统最大群组数得到的颗粒数量密度的预测值降低,而颗粒体积分数的预测值增加.随着颗粒体积增量的增加,修正最大群组数得到的颗粒数量密度的预测值下降,而颗粒体积分数的预测值呈振荡趋势.     

基于SPH方法的液体晃动数值模拟

采用光滑粒子流体动力学(SPH)方法模拟研究贮液容器中带有自由表面的液体晃动问题.应用密度对比法捕捉自由表面粒子,并采用边界排斥力法与镜像粒子相结合的方法处理固体边界.计算了贮液容器中粘性流体的自由晃动以及外部激励下的一阶模态晃动,得到了自由液面的波动和监测断面上的压强分布.将计算结果与理论值以及前人的实验数据对比表明,光滑粒子动力学方法作为一种无网格法,能较好地模拟带有自由液面的液体晃动问题.     

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

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

用差分和Monte Carlo求解相结合的方法计算有化学反应的活塞流

用差分和Monte Carlo求解相结合的方法对有化学反应的标量几率密度函数方程进行了求解,讨论了随机颗粒数和时间步长对计算结果的影响。计算结果表明,当化学反应速率是标量的线性函数时,结果与精确解很吻合;当反应速率满足Arrhenius定律时,结果也是合理的。     

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

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

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

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