[论文解读] Hydrodynamic simulations on a moving Voronoi mesh

Volker Springel|arXiv (Cornell University)|Sep 10, 2011

Fluid Dynamics Simulations and InteractionsEngineering参考文献 2被引用 21

一句话总结

本文提出了一种用于流体动力学模拟的新型移动Voronoi网格方法，通过可移动网格生成点的Voronoi剖分实现高精度、伽利略不变性且无畸变的拉格朗日流体动力学。该方法结合了非结构化网格上的二阶精度有限体积Godunov格式与精确黎曼求解器，相较于固定欧拉网格或SPH方法，在较低分辨率下展现出更优的精度与收敛性，尤其在涉及激波、接触间断和大尺度流动的问题中表现更佳。

ABSTRACT

At the heart of any method for computational fluid dynamics lies the question of how the simulated fluid should be discretized. Traditionally, a fixed Eulerian mesh is often employed for this purpose, which in modern schemes may also be adaptively refined during a calculation. Particle-based methods on the other hand discretize the mass instead of the volume, yielding an approximately Lagrangian approach. It is also possible to achieve Lagrangian behavior in mesh-based methods if the mesh is allowed to move with the flow. However, such approaches have often been fraught with substantial problems related to the development of irregularity in the mesh topology. Here we describe a novel scheme that eliminates these weaknesses. It is based on a moving unstructured mesh defined by the Voronoi tessellation of a set of discrete points. The mesh is used to solve the hyperbolic conservation laws of ideal hydrodynamics with a finite volume approach, based on a second-order Godunov scheme with an exact Riemann solver. A particularly powerful feature of the approach is that the mesh-generating points can in principle be moved arbitrarily. If they are given the velocity of the local flow, a highly accurate Lagrangian formulation of continuum hydrodynamics is obtained that is free of mesh distortion problems, while it is at the same time fully Galilean-invariant, unlike ordinary Eulerian codes. We describe the formulation and implementation of our new Voronoi-based hydrodynamics, and we discuss a number of illustrative test problems that highlight its performance in practical applications.

研究动机与目标

• 为解决现有流体动力学方法的根本局限性，例如SPH方法在激波解析上的不足以及欧拉格式缺乏伽利略不变性。
• 开发一种基于网格的方法，实现类似拉格朗日的行为，同时避免网格纠缠或拓扑伪影。
• 结合SPH的几何灵活性与基于网格方法的精度和守恒特性。
• 实现自适应分辨率，并准确模拟具有强密度对比和大尺度运动的天体物理流动。

提出的方法

• 网格通过一组可移动点的Voronoi剖分构建，这些点随流场连续变形以维持网格质量。
• 流体动力学方程采用二阶精度的有限体积Godunov格式求解，结合双曲守恒律的精确黎曼求解器。
• 网格生成点随局部流体速度移动，以实现近乎拉格朗日的格式，最大限度减少对流误差。
• 由于Voronoi剖分天然避免网格纠缠或蝴蝶结形单元，该方法确保了拓扑鲁棒性。
• 由于网格运动跟随流体，该方法完全具备伽利略不变性，消除了与速度相关的数值伪影。
• 该格式通过非结构化网格框架自然支持自适应分辨率和曲面边界条件。

实验结果

研究问题

• RQ1 基于Voronoi剖分的移动网格方法是否能在多维流体动力学中实现伽利略不变性并避免网格纠缠？
• RQ2 在标准流体动力学测试问题中，移动Voronoi网格方法的精度与固定欧拉网格和SPH代码相比如何？
• RQ3 在具有大尺度流动或接触间断的模拟中，该移动网格在多大程度上减少了对流误差并改善了收敛性？
• RQ4 在具有强密度对比和复杂流动形态（如Kelvin-Helmholtz不稳定性或云体粉碎）的问题中，该方法能否保持高精度？

主要发现

• 在云体粉碎测试中，移动Voronoi网格在比固定网格欧拉格式更低的分辨率下即可实现收敛，表明其具有更高的数值效率。
• 由于减少了对流误差和具备伽利略不变性，该移动网格方法始终预测出比固定网格格式更高的残余云体质量分数。
• 与SPH相比，该方法在激波和接触间断的解析上表现更优，SPH在云体破坏模拟中会低估质量损失。
• 该代码在空间和时间上均表现出二阶精度，无需人工粘性，并避免了SPH中常见的虚假表面张力效应。
• 在三维模拟中，每分辨率单元的计算速度约为SPH的2倍，固定网格欧拉格式的3–4倍，但其在较低分辨率下的更高精度可弥补此性能代价。
• 由于其几何灵活性与拓扑稳定性，该方法在复杂流动（包括曲面边界和强大尺度运动）中表现出良好的鲁棒性。