[论文解读] Review on Smoothed Particle Hydrodynamics: Methodology development and recent achievement

Chi Zhang, Yujie Zhu+2 more|arXiv (Cornell University)|May 6, 2022

Fluid Dynamics Simulations and Interactions参考文献 249被引用 5

一句话总结

本文综述了光滑粒子流体动力学（SPH）的最新进展，重点包括基于黎曼求解器的数值方法（含耗散控制）与高阶重构（MUSCL、WENO、MOOD）、高效的粒子邻近搜索与时间推进方法、带有沙漏模态控制的稳定总拉格朗日格式、基于多分辨率格式的流固耦合模拟，以及粒子松弛技术在网格与粒子生成中的创新应用。主要贡献在于对前沿SPH方法的全面整合，并通过严格的基准测试验证了其在精度、收敛性与鲁棒性方面的表现。

ABSTRACT

Since its inception, the full Lagrangian meshless smoothed particle hydrodynamics (SPH) method has experienced a tremendous enhancement in methodology and impacted a range of multi-physics applications in science and engineering. The paper presents a concise review on latest developments and achievements of the SPH method, including (1) brief review of theory and fundamental with kernel corrections, (2) the Riemann-based SPH method with dissipation limiting and high-order data reconstruction by using MUSCL, WENO and MOOD schemes, (3) particle neighbor searching with particle sorting and efficient dual-criteria time stepping schemes, (4) total Lagrangian formulation with stablized, dynamics relaxation and hourglass control schemes, (5) fluid-structure interaction scheme with interface treatments and multi-resolution discretizations, (6) novel applications of particle relaxation for mesh and particle generations. Last but not least, benchmark tests for validating computational accuracy, convergence, robustness and efficiency are also supplied accordingly.

研究动机与目标

• 提供计算流体力学与固体力学中SPH方法发展的简明、最新综述。
• 解决SPH中的关键挑战，包括一致性、稳定性、边界处理与收敛性，特别是在复杂多物理场场景中。
• 突出基于黎曼求解器的SPH在高阶重构与耗散限制方面的最新创新，以提升精度。
• 展示粒子松弛技术在复杂几何中生成高质量无结构网格与粒子分布的新应用。
• 通过全面的基准测试验证所提方法，评估其计算精度、收敛性、鲁棒性与效率。

提出的方法

• 采用基于黎曼求解器的SPH格式，利用近似黎曼求解器计算粒子界面处的通量，实现对激波与间断的精确捕捉。
• 在SPH框架内应用高阶数据重构技术——MUSCL、WENO与MOOD，以提升空间精度并减少数值耗散。
• 实施粒子排序与双准则时间推进策略，加速邻近粒子搜索，提升大规模模拟的计算效率。
• 采用总拉格朗日格式，结合稳定化方案、动态松弛与沙漏模态控制，以抑制虚假模态并提升固体力学中的稳定性。
• 提出一种基于水平集插值的表面边界方法，通过将粒子重新定位至φ = -0.5Δx的界面，实现对复杂几何边界的精确边界条件施加。
• 在无网格框架中利用粒子松弛技术，生成高质量、自适应、各向同性且全局优化的无结构网格，以及适用于复杂三维体的分布规则的粒子分布。

实验结果

研究问题

• RQ1 基于黎曼求解器的SPH结合高阶重构（MUSCL、WENO、MOOD）如何提升复杂含激波与间断流场模拟的精度与稳定性？
• RQ2 在大规模SPH模拟中，最有效的粒子邻近搜索与时间推进策略是什么？
• RQ3 如何通过稳定化手段抑制沙漏模态，确保总拉格朗日SPH格式在固体力学与大变形问题中的鲁棒性？
• RQ4 SPH中流固耦合的最优界面处理方法与多分辨率离散化格式是什么？
• RQ5 粒子松弛技术能否被有效重用于生成复杂几何中高质量的无结构网格与粒子分布？

主要发现

• 基于黎曼求解器、耗散控制与高阶重构（MUSCL、WENO、MOOD）的SPH方法显著提升了对激波与接触间断的分辨率，同时保持了良好的数值稳定性。
• 粒子排序与双准则时间推进策略显著提升了计算效率与可扩展性，尤其在粒子数量庞大的模拟中表现突出。
• 结合动态松弛与沙漏模态控制的总拉格朗日SPH格式能有效抑制虚假模态，实现对大变形与断裂问题的稳定模拟。
• 基于水平集插值的表面边界方法通过将粒子重定位至φ = -0.5Δx，成功在复杂几何上施加边界条件，实现贴体粒子分布。
• 粒子松弛技术可生成高质量、全局优化、自适应且各向同性的无结构网格，以及适用于高几何复杂度三维体的分布规则的粒子分布。
• 全面的基准测试结果表明，所提出的SPH方法在多种多物理场应用中均实现了高计算精度、最优收敛率、强鲁棒性与高效性。