[论文解读] 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等人提供的实验数据，以及同时包含比表面能和动力学效应的线性化可解性理论预测进行比较。
• 通过改变动力学各向异性，研究枝晶与致密分枝形貌之间的转变，揭示稳定枝晶生长的临界阈值。
• 在三维模拟中追踪界面演化和速度，利用快照和速度分布分析形貌稳定性和选择机制。

实验结果

研究问题

• RQ1 基于原子尺度推导的各向异性界面参数的相场模拟，能否在过冷镍中定量再现实验测得的枝晶尖端速度？
• RQ2 为何实验中观测到的V–ΔT曲线中的速度转折在仅采用纯扩散传输的标准相场模型中无法出现？
• RQ3 动力学各向异性如何影响枝晶生长的稳定性和形貌，特别是在向致密分枝结构过渡的区域？
• RQ4 当同时包含比表面能和动力学各向异性时，可解性理论预测与相场模拟结果的吻合程度如何？
• RQ5 镍凝固过程中观测到的速度转折是源于本征扩散物理机制，还是外部因素（如杂质或流体动力效应）所致？

主要发现

• 采用原子尺度推导的各向异性参数（εc = 0.018，εk = 0.13）的相场模拟，在多种过冷度下能定量匹配实验测得的枝晶尖端速度。
• 实验中观测到的V–ΔT曲线中的速度转折在当前相场模型中未被再现，表明存在缺失的物理机制。
• 为匹配实验数据，动力学各向异性εk必须足够大（≥0.13），与分子动力学模拟的预测一致。
• 当动力学各向异性εk = 0时，枝晶形态向致密分枝形态转变，该转变通过在临界过冷度下的鞍结分岔发生。
• 当同时包含比表面能和动力学各向异性时，可解性理论对稳态枝晶生长的预测与相场模拟结果定量吻合。
• 模拟中未出现速度转折，提示实验中观测到的行为可能由杂质或流动效应引起，而非本征扩散物理机制。