[論文レビュー] A sharp-interface model and its numerical approximation for solid-state dewetting with axisymmetric geometry

Quan Zhao|arXiv (Cornell University)|Nov 7, 2017

Fluid Dynamics and Thin Films参考文献 16被引用数 4

ひとこと要約

本稿では、表面拡散および接触線運動を支配要因とする固体状態におけるドレーピング現象を、軸対称な鋭い界面モデルとして厳密に提示する。エネルギー変分原理を用いて導出されたこのモデルは、半implicitなパラメトリック有限要素法（PFEM）を用いて、リム形成、穴の生成、トロイダル膜の収縮といった複雑な形態的変化をシミュレートしており、空洞形成が異方性パラメータσおよび初期膜形状に強く依存することを示している。

ABSTRACT

Based on the thermodynamic variation, we rigorously derive the sharp-interface model for solid-state dewetting on a flat substrate in the form of cylindrical symmetry. The governing equations for the model belong to fourth-order geometric curve evolution partial differential equations, with proper boundary conditions such that the total volume of the system is conserved and the total energy is dissipative during the time evolution. We propose a variational formulation for the sharp-interface model and then apply the parametric finite element method for solving it efficiently. Extensive numerical simulation results are presented lastly to demonstrate the morphological characteristics for solid-state dewetting.

研究の動機と目的

• 熱力学的変分を用いて、軸対称幾何における固体状態ドレーピングの数学的に厳密な鋭い界面モデルを導出すること。
• 動的境界条件を伴う4階の幾何的偏微分方程式を安定かつ効率的に解く数値スキームを開発すること。
• 3次元固体状態ドレーピングにおける、縁の収縮、ピンチオフ、空洞形成といった複雑な形態的変化を調査すること。
• 表面エネルギーの異方性および初期膜形状が平衡状態および動的形態に与える影響を分析すること。
• 軸対称極限における主要なダイナミクスを検証することで、将来的な完全3次元シミュレーションの基盤を提供すること。

提案手法

• 全表面エネルギー汎関数の第一変分を、軸対称制約下で行い、軸断面のための4階の幾何的曲線発展偏微分方程式が得られる。
• 境界条件は体積保存とエネルギー散逸を強制し、接触線運動は緩和された接触角条件に従う。
• エネルギー安定性を保証し、効率的な数値的取り扱いを可能とするための変分的定式化が導入された。
• パラメトリック有限要素法（PFEM）が、接線速度制御によるメッシュ品質の維持を目的とした半implicitな混合有限要素フレームワークで適用された。
• 動的メッシュ再配分が暗黙的に処理され、再メッシュ化を回避し、明示的マーカー粒子法に比べて大きな時間ステップを可能にした。
• 数値的シミュレーションは、弧長パラメータ化された軸断面曲線を対象とし、表面拡散と接触線運動の完全な結合が実現された。

実験結果

リサーチクエスチョン

• RQ1 表面エネルギーの異方性は、軸対称的固体状態ドレーピングの形態的変化にどのように影響するか？
• RQ2 トロイダル薄膜の収縮過程で空洞が形成される条件は何か？
• RQ3 初期膜形状と異方性パラメータσは、最終的な平衡形状をどのように決定するか？
• RQ4 縁の収縮やピンチオフといった、軸対称的固体状態ドレーピングにおける主要な不安定性および遷移メカニズムは何か？
• RQ5 PFEMは、安定的かつ効率的に、膜の破断や自己融合といった複雑なトポロジカル変化を捉えることができるか？

主な発見

• モデルは、膜/気相界面の動的変化をH⁻¹勾配流れ、接触線運動をL²勾配流れとして正しく捉えており、エネルギー散逸と体積保存を保証している。
• 数値的シミュレーションにより、4つの明確な形態的変化タイプが同定された：(1) 安定な単一膜、(2) 空洞を伴うトロイダル膜、(3) 空洞を伴わないトロイダル膜、(4) 中央膜と外側トロイダスへの破断。
• 等方的および異方的表面エネルギーに対するフェーズダイアグラムから、空洞形成が異方性パラメータσおよび初期膜のアスペクト比に強く依存することが明らかになった。
• PFEMは、従来のマーカー粒子法に比べ、はるかに大きな時間ステップと優れたメッシュ品質を実現する安定なシミュレーションを可能にした。
• 軸対称モデルは、軸方向収縮、リム形成、および方位角方向のレイリー様の破壊といった主要な不安定性を明らかにしたが、完全な3次元効果は捉えていない。
• モデルは、脂肪なトロイダル膜のみが、方位角方向の破壊よりも前に中心へ移動して合体可能であると予測しており、実験的観察と整合的である。