[论文解读] Toward a 2D SPH Multiphysic Code with Solid-Solid & Fluid Interactions for Industrial Related Problems

David Cébron, Jean‐François Sigrist|arXiv (Cornell University)|Sep 30, 2008

Fluid Dynamics Simulations and Interactions参考文献 15被引用 2

一句话总结

本文提出 SPHYSICS2，一种增强的开源二维光滑粒子流体动力学（SPH）代码，扩展用于涉及多相流与固体相互作用的工业级模拟。通过集成固-固接触检测与流-固耦合，该代码能够准确模拟复杂的碰撞与滑动动力学，其结果在堤坝溃决与圆柱体入水案例中得到验证，尽管存在轻微的数值阻尼，但能量保持一致。

ABSTRACT

In the present study, applications of the SPH method to industrial related issues are considered by starting from an existing open source 2D SPH code, namely the SPHYSICS code, which offers an effective ground for numerical developments, which are performed in order to bring an answer to industrial problems, such as simulations of solid/fluid coupling in a free surface flow context. The purpose of the present paper is therefore to expose the numerical developments which yield an enhanced version (referred to as "SPHYSIC2") of the initial code. Firstly, the different features added to obtain the operational code needed for engineering applications are described, and so are the problems raised on this way, offering a kind of review of SPH methods for engineers. Secondly, the validation of the proposed code is partially presented with two well known but difficult test cases, namely the classical "dam break" and "wedge entry" problems. Thirdly, principles of a method to solve solid/solid contacts, frequently present in realistic configuration, are exposed and applied to achieve more complex simulations. Finally, perspectives for new features of the SPHYSIC2 code are exposed and discussed.

研究动机与目标

• 将开源 SPHYSICS 代码扩展，以支持多相流与多固体相互作用，适用于工业预设计应用。
• 实现一种鲁棒的固-固接触检测算法，计算成本低且运行高效。
• 在涉及流-结构与固-固相互作用的基准问题上验证增强后的代码。
• 展示复杂、真实场景模拟的可行性，例如圆柱体撞击并滑动于水箱中的水体。
• 为未来在安全评估与动态冲击研究中实现 SPH 在工业领域的应用奠定基础。

提出的方法

• 在 SPHYSICS 中增强接触检测算法，基于粒子邻近性，通过粒子与固体边界之间的简单距离检测实现。
• 实现一种接触力模型，施加法向与切向力，以模拟非弹性碰撞与滑动行为。
• 修改 SPH 动量方程，引入人工粘性与扩散项以提升稳定性和物理真实性。
• 采用具有紧支撑与归一化特性的核函数，确保粒子插值的准确性与收敛性。
• 采用弱可压缩 SPH 公式，结合立方样条核函数模拟流体动力学。
• 通过壁面粒子方法施加重力与边界条件，以模拟固体边界。

实验结果

研究问题

• RQ1 SPHYSICS 是否能有效扩展以模拟具有工业相关性的固-固与流-固相互作用？
• RQ2 增强后的 SPHYSICS2 模型在圆柱体撞击水面过程中的运动与能量耗散方面，其重现精度如何？
• RQ3 接触算法引入的数值阻尼程度如何，其对能量守恒的影响为何？
• RQ4 该代码能否在自由下落、流体撞击与边界滑动等复杂过程中保持动力学的一致性？
• RQ5 SPH 方法在船舶与机械工程的预设计分析中，其可行性在多大程度上仍成立？

主要发现

• SPHYSICS2 代码成功模拟了圆柱体落入装满水的水箱过程，准确捕捉了撞击、减速以及在水箱壁面滑动等关键动力学行为。
• 自由下落阶段的速度历史与理论预测高度一致，仅观察到轻微的数值阻尼。
• 自由下落阶段能量守恒良好，撞击过程中的动能损失主要源于流体阻力，而非数值耗散。
• 在短时接触事件中，接触检测算法引入的数值阻尼可忽略不计，证实其适用于非弹性碰撞建模。
• 圆柱体在水箱壁面滑动及在水面反弹的模拟结果，在定性与定量上均与预期物理行为一致。
• 该代码在复杂多物理场工业模拟方面展现出良好可行性，尽管仍需进一步验证。