基于单层壁面粒子MPS方法的多相流及流固耦合算法研究.pdfVIP

摘要

以自由表面和界面大变形为特征的流动现象广泛存在于船舶及海洋工程中，如甲

板上浪现象、入水砰击现象以及破损船舶下沉入水时的涌流现象。这些流动现象有的

会对船舶及海洋工程结构造成损伤，有的则可能导致船舶沉没。因此，具有重要的研

究意义与价值。这些现象涉及不同流体间的相互影响以及流体与流场中结构的耦合作

用，具有强非线性及自由表面/流固界面的易变形性等特征，使得基于网格的方法需要

采用较复杂的数值格式才能实现有效的数值模拟。作为基于拉格朗日描述的一种粒子

类方法，移动粒子半隐式（MovingParticleSemi-Implicit，MPS）方法在处理大变形问

题时具有天然的优势。本文针对MPS求解压力时存在的震荡问题以及薄壁类结构适用

性差的问题进行了改进，提出了单层壁面粒子MPS（Single-layerWall-boundaryParticle

MPS，SWP-MPS）方法，并将之拓展至多相流和流固耦合应用，即开发了多相SWP-

MPS模型和基于多相SWP-MPS模型的流固耦合算法。本文主要研究内容分为以下三

部分：

（1）基于MPS方法框架，通过改进源项、边界条件以及粒子搜索方式，提出了

SWP-MPS方法，用于自由表面大变形问题的模拟。通过在自适应混合源项中添加压缩

项提出了微可压自适应混合源项（Slightly-compressibleAdaptiveMixedSourceterm，

SAMS），减少了数值噪声的累积，同时在自由表面采用概念粒子自由表面（Free

SurfacewithConceptualParticle，FSCP）边界条件，降低了表面粒子误判的影响，改善

了压力求解的稳定性；通过压力Neumann边界条件和Laplace算子补偿模型去除了固体

边界外的虚粒子，同时提出近壁面粒子搜索技术，解决了薄壁类结构粒子误搜索的问

题，克服了壁面边界条件导致的应用局限性。最后，对提出的模型及SWP-MPS方法进

行了验证，并采用该方法模拟了静水条件下二维破损舱体下沉入水的过程，探讨了破

损状态与舱内障碍物布置形式对舱体运动及舱室内部流动特性的影响。

（2）通过改善相界面密度及加速度的连续性，将SWP-MPS方法拓展至多相流问

题研究，提出了多相SWP-MPS模型，用于以自由表面/界面大变形为特征的多相流模

拟。提出了粒子数密度补偿（ParticleNumberDensityCompensation，PNDC）模型，以

去除固体边界外虚粒子的同时保证固体边界附近粒子数密度的计算准确性；通过一阶

精度密度光滑格式和多相密度模型提高了密度场的连续性，缓解了相界面处密度不连

续的影响；同时，采用多相粘度模型和连续压力梯度模型改善了加速度场的连续性。

最后，通过Rayleigh-Taylor不稳定性、液舱晃荡及溃坝等算例对所提多相SWP-MPS模

型进行了验证，该方法能够稳定、精确的模拟以界面大变形为特征的多相流问题。

（3）通过耦合刚体-模态叠加（CoupledRigid-bodyandModalSuperposition，

CRMS）模型，将SWP-MPS方法拓展至流固耦合问题研究，实现了涉及气、液、固三

相的强非线性流固耦合问题的有效求解。该算法采用多相SWP-MPS模型求解流体域，

采用CRMS模型求解固体域，实现了算法的全粒子离散；在流固界面采用了点对点信

息交互技术，简化了耦合条件的插值过程，提高了数值模拟过程中界面信息的传递效

率；同时，采用了Aitken松弛技术，以改善迭代过程中的收敛性。最后，通过破舱入

水冲击弹性舱壁、两相溃坝冲击弹性挡板等算例对所提流固耦合算法进行了验证，并

通过带弹性板两相液舱晃荡的模拟对晃荡过程中不同刚度挡板的变形量及其对晃荡特

性的影响进行了讨论。

关键词：移动粒子半隐式法；多相流；流固耦合；强非线性；大变形

ABSTRACT

Theflowphenomenacharacterizedbylargedeformationoffreesurfaceandinterfaceare

widelypresentinshipandoce