[論文レビュー] The Cahn-Hilliard equation on an evolving surface

David O'Connor, Björn Stinner|arXiv (Cornell University)|Jul 19, 2016

Solidification and crystal growth phenomena参考文献 22被引用数 5

ひとこと要約

本稿は、物質時間微分を用いて領域運動を扱うことで、変動する表面におけるCahn-Hilliard方程式の鋭い界面極限を、形式的マッチド漸近展開法を拡張して導出する。移動度スケーリングに応じて二つの異なる自由境界問題が特定される：定常移動度ではMullins-Sekerka型問題に輸送項が加わる一方、特異移動度では表面拡散方程式に曲率項と輸送項が寄与する。2次元および3次元設定において数値的に検証された。

ABSTRACT

We study the asymptotic limit of the Cahn-Hilliard equation on an evolving surface with prescribed velocity. The method of formally matched asymptotic expansions is extended to account for the movement of the domain. We consider various forms for the mobility and potential functions, in particular, with regards to the scaling of the mobility with the interface thickness parameter. Mullins-Sekerka, but also surface diffusion type problems, can be derived featuring additional terms which are due to the domain evolution. The asymptotic behaviour is supported and further explored with some numerical simulations.

研究の動機と目的

• Cahn-Hilliard方程式に対して、領域の幾何学的変化を考慮した移動表面における漸近解析手法の拡張を目的とする。
• 表面の変形（速度場および曲率を介して）が、相分離ダイナミクスの鋭い界面極限に与える影響を調査すること。
• 移動度スケーリング（定常移動度対特異移動度）に応じて、得られる自由境界問題を分類すること。
• 特に非定常状態において、表面変形が位相界面に引き起こすトポロジー的変化（自己交差や分裂など）を探索すること。
• 回転対称性を持つトーラスおよび球面の上での数値シミュレーションにより、理論的予測の妥当性を検証すること。

提案手法

• ε → 0 の極限において、形式的マッチド漸近展開を適用する。ここで ε は界面厚さパラメータである。
• 移動表面における内側（境界層）および外側（バルク）展開の一致を通じて、位相界面の時間発展を導出する。
• 支配的偏微分方程式における領域運動を考慮するため、物質時間微分 ∂•_t = ∂_t + v·∇ を組み込む。
• 表面微分作用素を用いる：∇_Γ(t)（表面勾配）、Δ_Γ(t)（Laplace-Beltrami作用素）、および表面発散。
• 二つの移動度領域を検討する：M(u) ∼ ε⁰（定常移動度）および M(u) ∼ ε⁻¹（特異移動度）、それぞれ異なる極限方程式を導出する。
• 表面速度に起因する化学ポテンシャル、平均曲率、および輸送項を含む、極限における鋭い界面問題を導出する。

実験結果

リサーチクエスチョン

• RQ1 領域の変動が、移動表面におけるCahn-Hilliard方程式の鋭い界面極限にどのように影響を与えるか？
• RQ2 移動度関数の異なるスケーリングのもとで、どのような異なる自由境界問題が出現するか？
• RQ3 漸近的手法を、表面速度および曲率を伴う移動領域に拡張可能か？
• RQ4 表面変形は、自己交差や分裂を含む位相界面のトポロジー的変化をどのように引き起こすか？
• RQ5 対数ポテンシャルのディープクエンチ極限において、速度場の表面発散が果たす役割は何か？

主な発見

• 定常移動度（M ∼ ε⁰）の場合、鋭い界面極限は、表面運動に起因する追加の輸送項を有するMullins-Sekerka型問題に帰着する。
• 特異移動度（M ∼ ε⁻¹）の場合、極限は表面拡散方程式に簡略化され、測地線曲率および領域変動に起因する輸送項を含む。
• 変動するトーラス上での数値シミュレーションでは、表面変形が単一の界面を自己交差させ、二つに分裂させることが確認され、静止表面では実現不可能なトポロジー的変化が生じる。
• 変動系におけるエネルギーの時間発展は、トポロジー的遷移の間には一時的に増加し、その後静止系よりも高いエネルギーに安定化することが示され、局所的最小エネルギー状態であることが示唆される。
• 対数ポテンシャルのディープクエンチ極限は、速度場が表面発散がゼロ（表面発散自由）である場合に限り、二重障害モデルに収束する。
• 回転対称性のある設定において、漸近的枠組みは正しいダイナミクスを正しく捉えており、数値結果は理論的予測を支持している。