[論文レビュー] Hydrodynamic simulations on a moving Voronoi mesh

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

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

ひとこと要約

本稿では、移動可能なメッシュ生成点のボロノイタイル化を用いて、高精度でガリレオ変換不変かつ歪みのないラグランジュ流体動力学を実現する、画期的な移動ボロノイメッシュ法を提案する。この手法は、非構造格子上での正確なリーマン解法を組み合わせた2次精度の有限体積Godunovスキームを採用しており、固定メッシュEuler型や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法における衝撃波の分解能不足や、Euler型コードにおけるガリレオ不変性の欠如を是正すること。
• メッシュのねじれやトポロジー的アーチファクトを回避する、ラグランジュ的挙動を実現するメッシュベースの手法を開発すること。
• SPHの幾何的柔軟性と、メッシュベース手法の精度および保存性を統合すること。
• 強い密度対比と大規模な流れを伴う天体物理的流れの適応的解像度と高精度なシミュレーションを可能にすること。

提案手法

• メッシュは、流れに応じて連続的に変形する可動点のボロノイタイル化として構築され、メッシュ品質を維持する。
• 流体動力学は、双曲型保存則に対する正確なリーマン解法を用いた2次精度の有限体積Godunovスキームで解かれる。
• メッシュ生成点を局所的な流体速度に従って移動させることで、ほぼラグランジュ的定式化を実現し、広告誤差を最小限に抑える。
• ボロノイタイル化の性質上、メッシュのねじれやバタフライ型セルを自然に回避するため、トポロジー的安定性が保証される。
• メッシュの運動が流体に従うため、速度依存の数値的アーチファクトが排除され、完全にガリレオ不変となる。
• 非構造格子フレームワークにより、適応的解像度と曲面境界条件の自然なサポートが可能となる。

実験結果

リサーチクエスチョン

• RQ1 ボロノイタイル化に基づく移動メッシュ法は、多次元流体動力学においてガリレオ不変性を達成し、メッシュのねじれを回避できるか？
• RQ2 標準的な流体動力学テスト問題において、固定メッシュEuler型およびSPHコードと比較して、移動ボロノイメッシュの精度はどの程度か？
• RQ3 大規模な流れや接触不連続を伴うシミュレーションにおいて、移動メッシュは広告誤差をどの程度低減し、収束性を向上させるか？
• RQ4 強い密度対比と複雑な流れの形態（例：ケルビン＝ヘルムホルツ不安定性やクラウドクランキング）を伴う問題において、本手法は高精度を維持できるか？

主な発見

• 固定メッシュEuler型コードと比較して、クラウドクランキングテストにおいて低解像度でも収束結果を達成しており、数値的効率性が優れていることが示された。
• 固定メッシュコードと比較して、移動メッシュ法は、広告誤差の低減とガリレオ不変性のおかげで、残留クラウド質量分率をより高く予測した。
• SPH法と比較して、衝撃波および接触不連続の分解能が向上しており、SPH法はクラウド破壊シミュレーションにおいて質量損失を過小評価していた。
• コードは空間的・時間的2次精度を示し、人工的粘性の必要がなく、SPHで一般的に見られる不自然な表面張力効果も回避した。
• 3次元では、解像度要素1つあたりSPHの約2倍、固定メッシュEuler型コードの3〜4倍の計算コストを要するが、低解像度でも高い精度を達成するため、コストの差は相殺される。
• 幾何的柔軟性とトポロジー的安定性のおかげで、曲面境界や強い大規模な流れを伴う複雑な流れに対しても、本手法は頑健である。