[论文解读] A Vortex Method for Bi-phasic Fluids Interacting with Rigid Bodies

Mathieu Coquerelle, Jérémie Allard+2 more|arXiv (Cornell University)|Jul 24, 2006

Computer Graphics and Visualization Techniques参考文献 30被引用 6

一句话总结

本文提出了一种新颖的涡粒子方法，用于基于涡度的公式模拟双相流体（如液态和气态）与刚性体之间的相互作用。通过结合涡粒子、水平集方法和浸入边界技术，该方法实现了复杂流体动力学（包括飞溅、湍流和气泡）的稳定、高精度模拟，并实现了高效的双向流体-刚性体耦合，经基准实验验证，其效率和视觉保真度优于先前的方法。

ABSTRACT

We present an accurate Lagrangian method based on vortex particles, level-sets, and immersed boundary methods, for animating the interplay between two fluids and rigid solids. We show that a vortex method is a good choice for simulating bi-phase flow, such as liquid and gas, with a good level of realism. Vortex particles are localized at the interfaces between the two fluids and within the regions of high turbulence. We gain local precision and efficiency from the stable advection permitted by the vorticity formulation. Moreover, our numerical method straightforwardly solves the two-way coupling problem between the fluids and animated rigid solids. This new approach is validated through numerical comparisons with reference experiments from the computational fluid community. We also show that the visually appealing results obtained in the CG community can be reproduced with increased efficiency and an easier implementation.

研究动机与目标

• 开发一种拉格朗日方法，以准确模拟涉及液态和气态的双相流体与刚性体之间的相互作用。
• 解决在无需复杂边界条件的情况下实现双向流体-结构耦合的挑战。
• 实现对飞溅、湍流、气泡和破浪等复杂现象的真实感模拟。
• 证明涡粒子方法在多相流模拟中既高效又准确，克服欧拉方法和仅粒子方法的局限性。
• 通过与既定的计算流体动力学基准和计算机图形学社区标准进行对比，验证该方法的有效性。

提出的方法

• 该方法采用纳维-斯托克斯方程的涡度公式，将计算资源集中于高涡度和湍流区域。
• 使用涡粒子表示流体涡度，将其定位在流体界面和湍流区域，以实现局部高精度。
• 水平集函数用于追踪两种流体之间以及流体与刚性体之间的界面，每时间步均执行迎风和重初始化操作。
• 浸入边界方法通过将流体力分布到固体表面，实现对刚性体上无滑移和无穿透过条件的强制满足。
• 采用三维结构化笛卡尔网格通过求解泊松方程计算速度场，并实现粒子重分布，以确保稳定性和效率。
• 通过在流体与固体之间强制保持速度连续性，计算刚性体运动，从而实现精确的双向耦合。

实验结果

研究问题

• RQ1 涡粒子方法能否有效模拟涉及液态和气态的双相流体流动，实现高保真度与高效率？
• RQ2 如何在无需复杂边界条件强制的情况下，稳健地建模双向流体-结构相互作用？
• RQ3 在多相流中，涡方法是否能比欧拉方法或SPH方法更好地保持涡度和动态特征？
• RQ4 水平集与涡粒子技术在多大程度上可结合使用，以模拟复杂界面和流体行为？
• RQ5 所提出的方法在计算效率上是否足够高，能够超越现有粒子或网格基方法，在复杂流体-刚性体相互作用中表现更优？

主要发现

• 该方法成功再现了波浪破碎、飞溅和气泡动力学等复杂流体现象，具有极高的视觉保真度。
• 在模拟两个下沉小球的场景中，其动力学行为与先前方法[5]相当，但性能提升了58倍，且使用更小的时间步长和更少的网格点。
• 性能分析显示，42.2%的计算时间用于求解泊松方程，而粒子推进和重分布仅占10.32%。
• 水平集操作占总计算时间的21.2%，表明界面处理是主要计算开销，但仍在可管理范围内。
• 该方法实现了稳定、准确的流体-结构耦合，无流体渗透现象，并在固体边界处保持了连续的速度场。
• 该方法能够高效模拟复杂场景，如金字塔撞击水墙以及不同姿态的杯子进入流体，充分展示了其鲁棒性和真实感。