[论文解读] A modified SPH approach for fluids with large density differences

F. Ott, Erik Schnetter|arXiv (Cornell University)|Mar 27, 2003

Fluid Dynamics Simulations and Interactions参考文献 9被引用 21

一句话总结

该论文提出了一种改进的光滑粒子流体动力学（SPH）方法，通过使用粒子密度而非质量密度进行离散化，以稳定模拟具有极大密度差异的多相流。通过重新定义密度重构以使用粒子密度，该方法消除了虚假压力振荡，并在激波管、声波传播和对流问题中提高了精度，尤其在粒子质量或空间分辨率差异显著时表现更优。

ABSTRACT

We introduce a modified SPH approach that is based on discretising the particle density instead of the mass density. This approach makes it possible to use SPH particles with very different masses to simulate multi-phase flows with large differences in mass density between the phases. We test our formulation with a simple advection problem, with sound waves encountering a density discontinuity, and with shock tubes containing a contact discontinuity between air and Diesel oil. For all examined problems where particles have different masses, the new formulation yields better results than standard SPH. This is also the case for problems in which different spatial resolutions are needed while the mass density does not change.

研究动机与目标

• 解决标准SPH在模拟具有大密度差异的多相流时出现的不稳定性和不准确性，特别是在流体界面处。
• 克服标准SPH中粒子间质量差异过大导致在不连续区域附近密度和压力计算不准确的根本性问题。
• 实现对质量不同、空间分辨率不同的粒子的鲁棒模拟，即使各相间质量密度保持恒定。
• 开发一种在具有刚性状态方程（如液-气或液-液界面）的可压缩流中保持稳定性和准确性的公式。
• 为涉及接触不连续、激波和高密度对比下的声波问题，提供一种数值稳定的标准SPH替代方案。

提出的方法

• 用基于粒子密度的公式替代标准SPH中基于质量密度的密度重构，其中粒子密度定义为粒子体积的倒数。
• 重新解释SPH核求和，以计算粒子密度而非质量密度，使用公式 $ \tilde{\rho}_{\text{part}}(\mathbf{x}_i) = \sum_j V_j W_{ij} $，其中 $ V_j $ 为粒子体积。
• 利用粒子质量 $ m_i $ 和粒子密度 $ \rho_{\text{part},i} $ 计算实际流体质量密度 $ \rho_i = m_i / V_i $，以确保与粒子运动的一致性。
• 将改进的公式应用于标准SPH运动方程，包括动量和能量演化，同时保持方法的拉格朗日特性。
• 采用人工粘性，$ \alpha = 0.5 $，以稳定激波前沿，作用范围约为8个光滑长度，且不会放大初始不连续性。
• 仅在必要时对初始数据进行平滑处理，但表明不连续的初始条件不会导致新公式的不稳定。

实验结果

研究问题

• RQ1 基于粒子密度而非质量密度的改进SPH公式是否能提升大密度对比下多相流的稳定性和准确性？
• RQ2 与标准SPH相比，新公式在模拟密度不连续对流时的表现如何？
• RQ3 改进的SPH方法能否准确捕捉高密度对比激波管（如空气与柴油油）中的激波和稀疏波？
• RQ4 该方法在粒子质量与空间分辨率不同时是否能避免虚假压力振荡？
• RQ5 当初始数据包含不连续性或高频模式时，新公式是否具有鲁棒性？

主要发现

• 改进的SPH公式成功模拟了一维密度不连续对流问题，未出现虚假压力振荡，且界面保持准确。
• 在声波测试中，新方法正确捕捉了透射与反射系数，而标准SPH则产生错误的波行为和压力尖峰。
• 在空气与柴油油（密度比约1:80）的激波管中，接触不连续面保持稳定且清晰分辨，且随时间推移无压力振荡增长。
• 激波速度（$ v_s = 1056 \pm 20 \, \text{m/s} $）、激波后速度（$ v_D = 5.94 \pm 0.02 \, \text{m/s} $）和压力升高（$ \Delta p_D = 4.94 \pm 0.02 \, \text{MPa} $）的模拟结果与解析解高度一致。
• 该方法与解析激波关系高度吻合，所有关键激波特性的相对误差低于2%，表现出高精度。
• 标准SPH在相同条件下无法产生稳定模拟，因出现非物理解压错误，其数值可超过实际压力达数个数量级，尤其在刚性流体中。