[論文レビュー] Perturbative solution of a propagating interface in the phase field model

Mao Hiraizumi, Shin‐ichi Sasa|arXiv (Cornell University)|Mar 25, 2021

Solidification and crystal growth phenomena参考文献 23被引用数 4

ひとこと要約

本稿では、過冷却状態（∆ < 1）における相場モデルにおける移動界面に対して摂動的解を提示する。ここで、顕熱発生が順序パラメータと温度場を結合する。次元なしパラメータ η = ξṘ/D を用いたスケーリングの仮定により、R(t) ∝ √t となる一次近似解が得られ、界面温度が界面速度 Ṙ に比例して平衡転移温度 Tc からずれることを示している。

ABSTRACT

When a stable ordered phase and a metastable disordered phase are separated by a flat interface, the metastable state changes to the stable state through the propagation of the interface. For cases in which latent heat is generated, the interface displacement during some time interval is proportional to the square root of the time interval when the extent of supercooling is less than a certain value. We demonstrate this behavior by deriving a perturbative solution for a propagating interface in the phase field model. We calculate the leading-order contribution explicitly, and find that the interface temperature deviates from the equilibrium transition temperature in proportion to the interface velocity.

研究の動機と目的

• ∆ < 1 の場合に数値的に観測された R(t) ∝ √t の時間依存性を理論的に導出すること。
• 過冷却条件下における相場モデルの移動界面に対する解析的解の欠如を解消すること。
• 顕熱の存在下で界面速度と Tc からの温度ずれの関係を確立すること。
• スケーリング座標と小さなパラメータ η = ξṘ/D を用いた摂動的枠組みを構築し、一次近似解に対する補正を体系的に計算すること。

提案手法

• 次元なし変数を導入：θ = cp(T - Tc)/L とし、式 (8) および (9) の相場モデル方程式を ξ と D を用いてスケーリングする。
• スケーリングされた空間座標として w = (x - R(t))/ξ と z = Ṙ(x - R(t))/D を定義し、界面および熱的長尺度を捉える。
• 正則摂動展開を可能にするために、小さな次元なしパラメータ η = ξṘ/D を導入する。
• 解のスケーリング形を仮定：φ(x,t) = Φ(w; η), θ(x,t) = Θ(z; η) とし、R(t) は一次近似で R(t) = √(2DCt) と進化すると仮定する。
• 一次近似方程式 (24) および (28) を解き、Φ0(w) と Θ0(z) を得る。ここで Θ0(z) = -δ(z) であり、Φ0(w) は非線形常微分方程式を満たす。
• 次の階層の補正 Φ1(w) および Θ1(z) の可解性条件を積分により導出し、Θ1(0) および C1 の明示的表現を得る。

実験結果

リサーチクエスチョン

• RQ1 相場モデルにおいて ∆ < 1 の場合に数値的に観測された R(t) ∝ √t の振る舞いは、解析的に導出可能か？
• RQ2 顕熱と有限な界面速度の存在下で、界面温度は Tc からどのようにずれるか？
• RQ3 次元なし過冷却パラメータ ∆ は、界面ダイナミクスおよび熱的プロファイルにどのように寄与するか？
• RQ4 η = ξṘ/D を用いた摂動的枠組みは、一次近似解に対する補正を体系的に計算可能か？

主な発見

• 一次近似解により R(t) = √(2DCt) が得られ、∆ < 1 の下で √t スケーリングが確認された。
• 界面温度 T(R(t), t) は界面速度 Ṙ に比例して Tc からずれ、そのずれは T(R(t), t) = Tc + (LξṘ/cpD)Θ1(0) で与えられる。
• 熱的プロファイルを制御するパラメータ C は ∆ → 1 に近づくと発散し、ダイナミクスにおける臨界転移を示している。
• 可解性条件を用いて Θ1(0) を明示的に計算し、Θ1(0) は自由エネルギーの微分を含む積分 I, J, K± によって決定される。
• D¨R/Ṙ³ の展開における係数 C1 はゼロであることが示され、摂動展開の整合性が保証された。
• Π = 0（Ts = Tc）の場合、C(∆) が計算され、∆ → 1 に近づくと発散することが示され、臨界過冷却閾値付近での √t スケーリングの崩壊が確認された。