Folch R, Plapp M
Laboratoire de Physique de la Matière Condensée, CNRS/Ecole Polytechnique, Palaiseau, France.
Phys Rev E Stat Nonlin Soft Matter Phys. 2005 Jul;72(1 Pt 1):011602. doi: 10.1103/PhysRevE.72.011602. Epub 2005 Jul 7.
A phase-field model that allows for quantitative simulations of low-speed eutectic and peritectic solidification under typical experimental conditions is developed. Its cornerstone is a smooth free-energy functional, specifically designed so that the stable solutions that connect any two phases are completely free of the third phase. For the simplest choice for this functional, the equations of motion for each of the two solid-liquid interfaces can be mapped to the standard phase-field model of single-phase solidification with its quartic double-well potential. By applying the thin-interface asymptotics and by extending the antitrapping current previously developed for this model, all spurious corrections to the dynamics of the solid-liquid interfaces linear in the interface thickness W can be eliminated. This means that, for small enough values of W, simulation results become independent of it. As a consequence, accurate results can be obtained using values of W much larger than the physical interface thickness, which yields a tremendous gain in computational power and makes simulations for realistic experimental parameters feasible. Convergence of the simulation outcome with decreasing W is explicitly demonstrated. Furthermore, the results are compared to a boundary-integral formulation of the corresponding free-boundary problem. Excellent agreement is found, except in the immediate vicinity of bifurcation points, a very sensitive situation where noticeable differences arise. These differences reveal that, in contrast to the standard assumptions of the free-boundary problem, out of equilibrium the diffuse trijunction region of the phase-field model can (i) slightly deviate from Young's law for the contact angles, and (ii) advance in a direction that forms a finite angle with the solid-solid interface at each instant. While the deviation (i) extrapolates to zero in the limit of vanishing interface thickness, the small angle in (ii) remains roughly constant, which indicates that it might be a genuine physical effect, present even for an atomic-scale interface thickness.
开发了一种相场模型,该模型能够在典型实验条件下对低速共晶和包晶凝固进行定量模拟。其核心是一个光滑的自由能泛函,该泛函经过专门设计,使得连接任意两相的稳定解完全不含第三相。对于该泛函的最简单选择,两个固液界面各自的运动方程可以映射到具有四次双阱势的单相凝固标准相场模型。通过应用薄界面渐近法并扩展先前为该模型开发的反俘获电流,可以消除对固液界面动力学的所有与界面厚度W成线性关系的虚假修正。这意味着,对于足够小的W值,模拟结果变得与W无关。因此,使用比物理界面厚度大得多的W值可以获得准确结果,这在计算能力上有巨大提升,并使针对实际实验参数的模拟变得可行。明确展示了模拟结果随W减小的收敛情况。此外,将结果与相应自由边界问题的边界积分公式进行了比较。除了在分岔点附近这一非常敏感的情况下出现明显差异外,发现两者具有极好的一致性。这些差异表明,与自由边界问题的标准假设相反,如果处于非平衡状态,相场模型的扩散三叉结区域可以(i)与接触角的杨氏定律略有偏差,并且(ii)在每个瞬间朝着与固 - 固界面形成有限角度的方向推进。虽然偏差(i)在界面厚度趋于零的极限情况下外推至零,但(ii)中的小角度大致保持不变,这表明它可能是一种真正的物理效应,即使对于原子尺度的界面厚度也存在。