[論文レビュー] Low artificial anisotropy cellular automaton model and its applications to the cell-to-dendrite transition in directional solidification

Lei Wei, Xin Lin+2 more|arXiv (Cornell University)|Oct 15, 2014

Solidification and crystal growth phenomena参考文献 1被引用数 3

ひとこと要約

本稿では、方向性凝固中の微細組織進化をシミュレートするため、限定的近隣固体分率（LNSF）捕獲則を用いた低人工的異方性セルオートマトン（CA）モデルを提示する。シミュレーションおよびSCN-aceおよびSCN-カムフラール系の実験データによる検証を通じて、溶質分配係数k₀が臨界閾値0.125を超えると、セルから枝晶への遷移速度（Vcd）がk₀に比例するスケーリング（Vcd ∝ Vc/k₀）から、定数値（Vcd = 8Vc）への遷移を示すことが明らかになった。

ABSTRACT

A low artificial anisotropy cellular automaton (CA) model is developed for the simulation of microstructure evolution in directional solidification. The CA model's capture rule was modified by a limited neighbor solid fraction (LNSF) method. Various interface curvature calculation methods have been compared. The simulated equilibrium shapes agree with the theoretical shapes, when the interface energy anisotropy coefficient is ε=0.01, ε=0.03 and ε=0.05, respectively.The low artificial anisotropy CA model is used in the numerical simulation of the cell-to-dendrite transition (CDT) in directional solidification. The influence of physical parameters (Γ, Dl, k0, ml) on CDT has been investigated. The main finding in this paper is the discovery of the changing behavior of the Vcd when the solute partition coefficient k0 is larger than a critical value. When k0 is less than 0.125, the Vcd follows the Kurz and Fisher criterion Vc/k0; while when k0>0.125, the Vcd equals to 8Vc. The experimental data of succinonitrile-acetone (SCN-ace, k0=0.1) and SCN-camphor (k0=0.33) support the conclusion from CA simulations.

研究の動機と目的

• 方向性凝固における微細組織の正確なシミュレーションを可能にする、最小限の人工的異方性を有するセルオートマトンモデルの開発。
• Γ、Dl、k₀、mlといった主要な物理パラメータがセルから枝晶への遷移（CDT）に与える影響の調査。
• 溶質分配係数k₀の関数としてのセルから枝晶への遷移速度（Vcd）の遷移挙動の解明。
• k₀値が異なるサリチルニトリル系系における実験データとシミュレーション結果の一致を確認すること。

提案手法

• 人工的異方性の低減を目的として、CA捕獲則を修正するための限定的近隣固体分率（LNSF）法を導入した。
• 枝晶成長の正確なシミュレーションを確保するため、複数の界面曲率計算手法を評価・比較した。
• 低異方性係数（ε = 0.01, 0.03, 0.05）における理論的平衡形状との一致を確認するため、モデルを理論的形状でキャリブレーションした。
• k₀、Dl、ml、Γを含む多様な物理パラメータの範囲で、セルから枝晶への遷移の数値シミュレーションを実施した。
• k₀ = 0.1（SCN-ace）およびk₀ = 0.33（SCN-カムフラール）の系における方向性凝固をモデルで再現し、実験観察と一致した。

実験結果

リサーチクエスチョン

• RQ1 溶質分配係数k₀は、方向性凝固におけるセルから枝晶への遷移速度（Vcd）にどのように影響を与えるか？
• RQ2 Vcdの挙動がVc/k₀スケーリングから定数値に遷移するk₀の臨界値は何か？
• RQ3 低人工的異方性CAモデルは、低異方性条件下で理論的平衡形状をどの程度正確に再現できるか？
• RQ4 SCN-aceおよびSCN-カムフラール系におけるシミュレートされたVcd値は、実験データとどの程度一致するか？
• RQ5 界面エネルギー異方性（ε）は、成長する枝晶の安定性および形態にどのような役割を果たすか？

主な発見

• 低異方性係数（ε = 0.01, 0.03, 0.05）における理論的平衡形状が、モデルにより正確に再現され、最小限の人工的異方性が確認された。
• k₀ < 0.125の場合、セルから枝晶への遷移速度（Vcd）は、KurzとFisherの基準に従い、Vc/k₀に比例して増加する。
• k₀ > 0.125の場合、Vcdは定数値8Vcに安定化し、成長ダイナミクスに顕著な遷移が生じることを示した。
• SCN-ace（k₀ = 0.1）およびSCN-カムフラール（k₀ = 0.33）系におけるシミュレートされたVcd値は、実験観察と良好に一致し、モデルの妥当性が裏付けられた。
• k₀ = 0.125における遷移は、溶質拡散制限から熱力学的制御への成長メカニズムの本質的転換を示している。