[論文レビュー] Calculation of interface curvature with the level-set method

Karl Yngve Lervåg|arXiv (Cornell University)|Jan 1, 2011

Fluid Dynamics and Heat Transfer参考文献 7被引用数 4

ひとこと要約

本稿では、ドロップレットの融合のようなトポロジー的変化を伴う領域においても、界面の曲率および法線ベクトルを計算するための頑健で実装が簡単な離散化スキームを提案する。幾何学的特徴を考慮した曲率推定と、折りたたみ部における方向差分の適応により、標準的な有限差分法では生じる不自然な曲率スパイクや誤った圧力場を解消し、トポロジー変化を伴う二相流の正確なシミュレーションを可能にする。

ABSTRACT

The level-set method is a popular method for interface capturing. One of the advantages of the level-set method is that the curvature and the normal vector of the interface can be readily calculated from the level-set function. However, in cases where the level-set method is used to capture topological changes, the standard discretization techniques for the curvature and the normal vector do not work properly. This is because they are affected by the discontinuities of the signed-distance function half-way between two interfaces. This article addresses the calculation of normal vectors and curvatures with the level-set method for such cases. It presents a discretization scheme that is relatively easy to implement in to an existing code. The improved discretization scheme is compared with a standard discretization scheme, first for a case with no flow, then for a case where two drops collide in a shear flow. The results show that the improved discretization yields more robust calculations in areas where topological changes are imminent.

研究の動機と目的

• トポロジー的変化が生じる場合に、標準的な有限差分スキームを用いたレベルセット法における曲率および法線ベクトル計算の不正確さを是正すること。
• 近接する界面間の折りたたみ部における符号付き距離関数の不連続性に起因する曲率スパイクの問題を解決すること。
• 既存のレベルセットコードベースへの大規模なアーキテクチャ変更なしに、容易に統合可能な、頑健で正確な曲率離散化を開発すること。
• 静的および動的二相流れケース、特に融合付近において、標準的手法と比較して改良されたスキームの有効性を検証すること。
• 改良された手法が、曲率誤差に起因する誤った圧力場を防ぐことで、物理的融合を可能にすることを示すこと。

提案手法

• マククリン＆ルーウェンブルグ（2011）の方向差分法を用いて法線ベクトルを計算し、標準的な勾配法よりも折りたたみ部をより頑健に扱う。
• マククリンら（2012）のインスパイアされた幾何学的特徴を考慮した曲率離散化を採用し、レベルセット関数の2階微分に依存せずに、局所的な界面幾何をフィットして曲率を推定する。
• 固定グリッド上の有限差分フレームワークに改良された曲率スキームを適用し、既存のレベルセット実装への変更を最小限に抑える。
• 連続表面力（CSF）法を用いて界面張力力をモデル化した、2次元非圧縮性二相流れソルバにこの手法を実装する。
• 符号付き距離レベルセット関数を用い、再初期化を適用して距離関数の性質を維持し、数値誤差を低減する。
• ナビエ＝ストークス方程式およびレベルセットの対流を固定グリッド有限差分スキームで解き、各タイムステップで曲率および法線ベクトルを計算する。

実験結果

リサーチクエスチョン

• RQ1 改良された曲率離散化は、2つの近接した界面間の折りたたみ部における不自然な曲率スパイクを排除できるか？
• RQ2 近接する2つのドロップレット間の薄膜領域における圧力場は、改良された手法によってどのように変化するか？
• RQ3 標準的手法が失敗するせん断流れシミュレーションにおいて、改良された離散化は物理的ドロップレット融合を可能にするか？
• RQ4 曲率の精度と安定性の観点から、改良されたスキームは標準的な有限差分法と比較してどの程度優れているか？
• RQ5 改良された曲率離散化は、大規模な変更なしに、既存のレベルセットコードに頑健かつ実用的に実装可能か？

主な発見

• 標準的な有限差分曲率離散化は、近接する2つのドロップレットの折りたたみ部で顕著な曲率スパイクを生じさせ、物理的に不自然な圧力歪みを引き起こす。
• 改良された曲率スキームは、折りたたみ部における曲率スパイクを効果的に排除し、滑らかで物理的に整合性のある曲率分布を実現する。
• 標準手法では、曲率スパイクに起因する誤った圧力場が融合を抑制し、せん断流れケースでドロップレットの融合を阻止する。
• 改良された離散化により、薄膜領域における圧力場は物理的に正確な分布を示す——中心部が高く、縁部が低く、外向きの流れを誘発し、融合を促進する。
• 改良された手法のシミュレーションでは融合が成功するが、標準手法では曲率誤差に起因する不自然な力により融合が阻止される。
• 改良された曲率スキームは最小限のコード変更で実装可能であり、二相流れシミュレーションにおける既存のレベルセット実装への実用的で頑健な改善手段である。