[論文レビュー] Explicit temperature coupling in phase-field crystal models of solidification

Maik Punke, Steven M. Wise+2 more|arXiv (Cornell University)|May 30, 2022

Solidification and crystal growth phenomena参考文献 54被引用数 10

ひとこと要約

本稿では、温度依存的な格子定数を用いて温度の時間発展と熱膨張を明示的に結合する、相場結晶（PFC）モデルを提示する。このモデルにより、凝固過程における熱応力と形態変化のシミュレーションが可能になる。古典的PFCモデルを拡張し、熱フラックスと温度依存周期性を組み込み、Eshelbyインクルージョンおよび枝晶成長シミュレーションにより妥当性を検証した。熱応力の予測精度と、温度変化に伴う格子間隔の変調が正確に再現された。

ABSTRACT

We present a phase-field crystal (PFC) model for solidification that accounts for thermal transport and a temperature-dependent lattice parameter. Elasticity effects are characterized through the continuous elastic field computed from the microscopic density field. We showcase the model capabilities via selected numerical investigations which focus on the prototypical growth of two-dimensional crystals from the melt, resulting in faceted shapes and dendrites. This work sets the grounds for a comprehensive mesoscale model of solidification including thermal expansion.

研究の動機と目的

• ミクロスケールのPFCモデルを構築し、熱輸送と結晶成長ダイナミクスを明示的に結合する。
• PFC自由エネルギー関数に温度依存格子定数を組み込み、熱膨張効果を表現する。
• 拡散的時間スケールにおける凝固結晶内の熱応力と弾性緩和をシミュレート可能にする。
• 既知の解析解（例：Eshelbyインクルージョン）との比較による妥当性評価を行い、枝晶成長および多面体成長の予測能力を示す。

提案手法

• 波数 q0 を β(T, T0)q0 に置き換えた修正されたPFC自由エネルギー関数を用いる。β は εth(T, T0) = αth(T − T0) を通じて熱膨張を記述する。
• 連成された時間発展方程式には、密度場 ψ に対する保存的 H−1 勾配流れと、温度場 T に対する拡散反応方程式が含まれる。熱フラックスは熱拡散率 M を用いてモデル化される。
• 温度場 T は無次元量であり、融点で T = 1 に正規化される。ψ と T 間の明示的結合は γ∂tT − ϑ∂tψ 項によって実現される。
• 数値解法には、スペクトル空間離散化と後退オイラー時刻積分を用いたフーリエ擬スペクトル法が用いられる。
• 連続的弾性場は、密度場 ψ からひずみテンソルを介して導出され、機械的応力成分 σij の計算が可能になる。
• シミュレーションは、空間的に一様な T = T0 と、液体（A）または固体（A, φ±）相に設定された初期条件で開始され、ハニカムおよび三角格子構造を再現するようにパラメータが調整される。

実験結果

リサーチクエスチョン

• RQ1 温度依存的格子間隔は、液体からの二次元結晶成長過程における形態と応力の進化にどのように影響するか？
• RQ2 PFCモデルは、固体相における熱膨張に起因する熱応力をどの程度正確に再現できるか？
• RQ3 熱フラックスと温度発展を組み込むことで、等温PFCモデルと比較して枝晶成長および多面体成長パターンにどのような影響を与えるか？
• RQ4 熱膨張および温度勾配が存在する条件下で、Eshelbyインクルージョン問題を正確に捉えることができるか？

主な発見

• 本モデルはEshelbyインクルージョン問題を正確に再現した。応力場の解析解と数値解の間で良好な一致が得られた。
• 熱応力 σα_ij が正しく予測された：T > T0 の場合正、T < T0 の場合負。A1（T < T0）とA2（T > T0）で初期化した結晶では符号が逆転し、物理的整合性が確認された。
• 格子間隔は温度に依存し、σα_ij の大きさは約 10⁻⁵ に達する。これにより、熱膨張が界面間隔および結晶融合ダイナミクスに与える影響が明確に示された。
• 温度勾配下では、自然に枝晶および多面体成長形態が出現し、固体-液体界面で明確な応力パターンが観察された。
• 液体相では応力場 σij が消滅し、固体相では空間的に強く変動する。σxx(A1) と σxx(A2) はそれぞれ σ0_ij と σα_ij 成分に分解され、対称的および反対称的であることが明らかになった。
• 本モデルにより、弾性緩和と拡散的ダイナミクスの分離が可能となり、標準的PFCと同様に緩和が遅い挙動を示した。これにより、応力寄与の解釈が大幅な手法変更なしに行えるようになった。