改进的黎曼SPH方法及其在大密度比多相流中的应用研究.pdfVIP

改进的黎曼SPH方法及其在大密度比多相流中的应用研究

摘要

多相流广泛存在于船舶与海洋工程领域，如气泡上浮、结构入水砰击、水下爆炸等，

这些现象常涉及流体大变形、多相界面掺混、非线性流固耦合作用等问题。作为目前应

用最为广泛的无网格数值方法，光滑粒子流体动力学（SPH）方法得益于其天然的拉格

朗日特性，避免了网格畸变和界面捕捉等难题，在处理流体飞溅、自由面大变形以及多

相掺混问题中具有明显优势，因此是求解多相流问题的优势方法。近年来，SPH方法在

计算精度和效率方面不断完善，并在船海领域多相流问题应用中取得了丰硕的研究成果，

但仍存在一些问题有待解决。本文针对船海领域多相流问题，开发了适用于大密度比多

相流的黎曼SPH模型和算法，重点解决了无法考虑气相真实可压性、强可压缩流场计算

精度低、缺乏高精度的热传导模型等难题，实现了对计及气垫效应的结构入水砰击、气

泡与船体耦合作用、双水下高压气泡与边界耦合作用、以及气液多相流热传导等问题的

高精度模拟，为船海多相流问题研究提供了基础性技术支撑。

本文首先介绍了SPH方法的基本理论，包括核近似和粒子近似思想，以及基于这些

思想推导得到的SPH控制方程。然后，在SPH方法中引入黎曼思想，认为每对粒子中

点存在间断，粒子对组成激波管问题，那么粒子之间的相互作用可以采用黎曼求解器得

到。最后介绍了本文中采用的一些SPH数值处理技术，包括边界条件、多相流粒子位移

修正技术、多相流密度过滤技术以及SPH前后处理方法等，为后续的多相流问题研究奠

定了基础。

针对传统多相流SPH模型无法考虑气体真实可压性，对于存在“气垫效应”的入水

砰击问题载荷求解精度较差的难题，基于PVRS黎曼解，提出了一种低耗散SPH限制

器，建立了适用于大密度比多相流的黎曼SPH数值模型，克服了传统SPH模型无法考

虑真实气体压缩性的难题，实现了结构小攻角入水“气垫效应”的高精度重现，显著提

高了多相流砰击载荷的求解精度。通过模拟平板砰击试验验证了所开发的低耗散黎曼

SPH多相流数值模型的精度和收敛性，并讨论了有无空气存在以及不同板长对平板砰击

载荷的影响。随后计及砰击中的气垫效应，模拟了LNG储罐保温板的砰击问题，并分

析了砰击速度和砰击角度对该保温板砰击载荷特性的影响。最后，在三维实尺度下，考

虑空气影响，对水上飞机的迫降过程进行模拟，并分析了水上飞机的运动轨迹以及迫降

引起的不同位置砰击载荷特性。

针对黎曼多相流SPH模型无法考虑流体粘性的问题，基于建立的低耗散黎曼SPH

I

哈尔滨工程大学博士学位论文

模型，推导了粘性耗散限制器参数和雷诺数之间的关系，开发了一种稳健的适用于常压

大尺度上浮气泡的多相粘性流动黎曼SPH数值模型。通过模拟经典的气泡上浮过程，并

与传统多相流SPH方法结果进行对比，验证了新模型能够有效地提升计算精度。在此基

础上，考虑船体结构，建立单个和多个上浮气泡及其与船体的耦合动力学模型，对大尺

度上浮气泡及气泡群与船体的耦合作用进行了模拟，并探究了气泡大小和深度、气泡与

船体的相对位置、气泡群密度等参数对气泡和船体动力学特性的影响，为气泡沉船研究

提供了基础性技术支撑。

针对多相流黎曼SPH方法求解强可压缩流动问题时仍存在粒子体积剧烈变化导致

精度差，以及长时间模拟效率低等难题，在轴对称SPH框架下引入黎曼求解器，并提出

了高精度多相流SPH压力梯度算法，结合自适应的粒子分割与合并算法，建立了适用于

水下高压气泡脉动与射流的强可压缩多相SPH数值模型。同时为了提高计算效率，采用

了渐变粒子分布，并针对该分布提出了高效链表搜索算法。首先通过几个经典算例，包

括激波管问题、入水问题和水下爆炸气泡问题，验证了当前建立的强可压缩SPH数值模

型的准确性。进而，模拟了近壁面双大尺度水下爆炸气泡的脉动和射流过程，获得了双

气泡与固壁间的耦合动力学特性，揭示了距离和相位参数对流场载荷的影响规律。

针对传统SPH方法模拟气液多相流热力学问题精度较低的难题，计及热扩散比对

热传导效应的影响，引入黎曼思想处理热力学场间断难题，提出了适用于求解气液多相

流问题的黎曼SPH热传导模型，提高了SPH热传导问题的求解精度。通过模拟几个典

型的一维和二维热传导问题，并与传统的SPH热传