[論文レビュー] Finite-Element Approximation of One-Sided Stefan Problems with Anisotropic, Approximately Crystalline, Gibbs--Thomson Law

John W. Barrett, Harald Garcke+1 more|arXiv (Cornell University)|Jan 9, 2012

nanoparticles nucleation surface interactions被引用数 8

ひとこと要約

本稿では、異方的で近似的に結晶的なギブズ＝トムソン則を伴う片側ステファンおよびムリンズ＝セケルカ問題をシミュレートするための新規な有限要素法を提示する。この手法は、パrametric境界の進化とバルク場の近似を結合し、滑らかさの処理を要せず安定性を確保するとともに、優れたメッシュ品質を維持する。これにより、斜面の破壊を再現できる物理的に厳密な3次元雪氷晶成長の初回のシミュレーションが可能になった。

ABSTRACT

We present a finite-element approximation for the one-sided Stefan problem and the one-sided Mullins--Sekerka problem, respectively. The problems feature a fully anisotropic Gibbs--Thomson law, as well as kinetic undercooling. Our approximation, which couples a parametric approximation of the moving boundary with a finite element approximation of the bulk quantities, can be shown to satisfy a stability bound, and it enjoys very good mesh properties which means that no mesh smoothing is necessary in practice. In our numerical computations we concentrate on the simulation of snow crystal growth. On choosing realistic physical parameters, we are able to produce several distinctive types of snow crystal morphologies. In particular, facet breaking in approximately crystalline evolutions can be observed.

研究の動機と目的

• 異方的で近似的に結晶的な界面条件を有する片側ステファンおよびムリンズ＝セケルカ問題に対して、安定でメッシュを保全する有限要素法の開発。
• 拡散および界面動力学に基づく熱力学的に整合性のあるモデルを用いて、3次元における現実的な雪氷晶成長のシミュレーション。
• 板状、柱状、デンドリテ、斜面の破壊を含む、従来の手法では困難な複雑な結晶形状の再現。
• 雪氷晶形成に関連する物理的パラメータを用いた数値実験を通じて、手法の精度と頑健性の検証。
• 移動境界アプローチを用いて、詳細な表面特徴を有する高構造的でファセットを持つ雪氷晶のシミュレーションの可能性の提示。

提案手法

• 移動結晶境界の進化にパラメトリック有限要素法を用い、頂点が均等に配置され、メッシュ品質が維持されるようにする。
• この手法は、気相または液相反で標準的な有限要素離散化を用いたバルク拡散方程式と結合する。
• 異方的ギブズ＝トムソン則は、Wulff形状に基づく表面エネルギー汎関数を介して組み込み、界面速度条件に運動的冷却を含める。
• 界面速度は、曲率、冷却度、および異方的表面エネルギーを含む法線速度条件から変分原理に基づいて導出される。
• 半陰的時間離散化を採用し、固定されたバルクメッシュと、精度を維持するための境界近似における適応的メッシュ細分化を実施する。
• この手法は離散的安定性バインドを満たしており、実際にはメッシュ品質の内因的保持のおかげで、メッシュの滑らかさ処理が不要であることが示された。

実験結果

リサーチクエスチョン

• RQ1 異方的で近似的に結晶的なギブズ＝トムソン則を伴う片側ステファンおよびムリンズ＝セケルカ問題を正確にシミュレートできる有限要素法を開発できるか？
• RQ2 メッシュの滑らかさ処理を要せず、長時間にわたる進化においても良好なメッシュ品質を維持できるか？
• RQ3 現実的な物理的パラメータ下で、実験的に観察された雪氷晶の形状、特に斜面の破壊を再現できるか？
• RQ4 ある値の異方性パラメータ γ_TB に対して、進化が自己相似的であるか、数値的に検証可能か？
• RQ5 本手法は、デンドリテ、キャップドカラム、板状のスクロールを含む、複雑な3次元雪氷晶構造を高い忠実度でシミュレートできるか？

主な発見

• 本手法は、板状、柱状、中空の柱、デンドリテ、キャップドカラム、板状のスクロールなど、広範な雪氷晶形状をシミュレートでき、実験的観察と一致した。
• ほぼ結晶的であるとされる異方性下で、斜面の破壊が数値的に観測され、理論的予測および現実の雪氷晶挙動と整合的であった。
• γ_TB ≈ 0.92 の場合、円柱状結晶の進化は自己相似的であることが示され、高さ対底面直径比が期待値に収束した。これは安定で対称的な成長モードを示している。
• 進化の全過程にわたり優れたメッシュ品質が維持され、メッシュの滑らかさ処理が不要であり、固定されたバルクメッシュで安定な時間積分が達成された。
• 本手法による数値結果は、バルクに拡散方程式を用いる移動境界アプローチを用いた3次元雪氷晶成長シミュレーションとして、これまでで最も複雑かつ高精度なものである。
• シミュレーションは、ステップやリッジなどの複雑な表面特徴を捉える能力を確認したが、計算コストが高いため、将来的には高解像度および確率的パrameter化による改良の余地がある。