[論文レビュー] Linking Phase-Field and Atomistic Simulations to Model Dendritic Solidification in Highly Undercooled Melts

J. Bragard, Alain Karma+2 more|arXiv (Cornell University)|Dec 10, 2001

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

ひとこと要約

本研究は、分子動力学から導かれた界面の異方性を用いて、過冷却ニッケル溶融状態における枝晶凝固を、位相場法と原子論的シミュレーションを統合してモデル化する。運動的異方性が実験的枝晶先端速度を定量的に再現するために不可欠であることが示され、標準的な拡散モデルでは観察されたV–ΔT曲線の速度折り返しを再現できないことが判明し、不純物や流れの効果が関係している可能性を示唆する。

ABSTRACT

Even though our theoretical understanding of dendritic solidification is relatively well developed, our current ability to model this process quantitatively remains extremely limited. This is due to the fact that the morphological development of dendrites depends sensitively on the degree of anisotropy of capillary and/or kinetic properties of the solid-liquid interface, which is not precisely known for materials of metallurgical interest. Here we simulate the crystallization of highly undercooled nickel melts using a computationally efficient phase-field model together with anisotropic properties recently predicted by molecular dynamics simulations. The results are compared to experimental data and to the predictions of a linearized solvability theory that includes both capillary and kinetic effects at the interface.

研究の動機と目的

• 未知の界面異方性に敏感に依存する形態的進化を示す、過冷却状態における枝晶凝固を定量的にモデル化するという長年の課題に対処すること。
• 原子論的シミュレーションから得た異方性界面特性を位相場シミュレーションに組み込むことで、理論的モデルと実験的観察のギャップを埋めること。
• ニッケルのV–ΔT曲線で広く観察される速度折り返し現象が、標準的な拡散輸送モデルによって説明可能かどうかを検証すること。
• 分子動力学から得た正確な界面パラメータを用いて、実験データおよび可解性理論と照合するための位相場モデルを検証すること。

提案手法

• 位相場モデルは、分子動力学シミュレーションから導かれた異方性キャピラリー係数および運動的係数を用いて、過冷却ニッケル溶融状態における枝晶凝固を模擬する。
• 界面エネルギー異方性εc = 0.018および運動的異方性εk = 0.13を原子論的研究から得て、曲率および運動的異方性を含むギブズ＝トムソン条件をモデルに組み込む。
• 過冷却度（ΔT/(L/cp) = 0.2 から 0.8）および界面厚さ（λ = 2 から 30）の範囲でシミュレーションを実施し、定常状態における速度および形態を解析する。
• LumらおよびWillneckerらの実験データ、およびキャピラリーおよび運動的効果を両方含む線形化された可解性理論の予測と結果を比較する。
• 運動的異方性を変化させることで、枝晶構造と密な分岐構造との間の遷移を検討し、安定な枝晶成長のための臨界閾値を特定する。
• 3次元シミュレーションで界面の時間発展および速度を追跡し、スナップショットおよび速度プロファイルを用いて形態的安定性および選択メカニズムを解析する。

実験結果

リサーチクエスチョン

• RQ1 原子論的シミュレーションから導かれた異方性界面パラメータを用いた位相場シミュレーションは、過冷却ニッケルにおける実験的枝晶先端速度を定量的に再現できるか？
• RQ2 標準的な位相場モデルが純粋な拡散輸送を仮定する場合、実験的に観察されるV–ΔT曲線の速度折り返し現象がなぜ再現されないのか？
• RQ3 運動的異方性は、特に密な分岐構造への遷移付近における枝晶成長の安定性および形態的特性にどのように影響を与えるか？
• RQ4 キャピラリーおよび運動的異方性を両方含む場合、可解性理論の予測と位相場シミュレーションの結果がどの程度一致するか？
• RQ5 ニッケル凝固における観察された速度折り返し現象は、内在的な拡散的物理法則に起因するのか、それとも不純物や流体力学的効果といった外部要因に起因するのか？

主な発見

• 原子論的シミュレーションから導かれた異方性パラメータ（εc = 0.018, εk = 0.13）を用いた位相場シミュレーションは、さまざまな過冷却度において実験的枝晶先端速度を定量的に再現した。
• 実験的に観察されるが、標準的な拡散モデルでは再現されないV–ΔT曲線の速度折り返し現象は、現在の位相場モデルでも再現されないため、物理的メカニズムが欠落していることを示唆する。
• 実験データと一致させるためには、運動的異方性εkが大きく（≥0.13）なる必要があることが確認され、分子動力学シミュレーションからの予測を裏付けた。
• 運動的異方性が低い（εk = 0）場合、枝晶構造から密な分岐構造への遷移が発生し、臨界過冷却度で鞍点分岐を介して発生する。
• キャピラリーおよび運動的異方性を両方含む場合、定常状態における枝晶成長の可解性理論の予測は、位相場シミュレーション結果と定量的に一致した。
• シミュレーションで速度折り返しが再現されないことは、実験的挙動の背後にある要因が、内在的な拡散的物理法則ではなく、不純物や流れの効果である可能性を示唆する。