[论文解读] Perturbative solution of a propagating interface in the phase field model

Mao Hiraizumi, Shin‐ichi Sasa|arXiv (Cornell University)|Mar 25, 2021

Solidification and crystal growth phenomena参考文献 23被引用 4

一句话总结

本文针对超冷条件（∆ < 1）下相场模型中传播界面提出了微扰解法，其中潜热生成使序参量场与温度场相互耦合。通过引入以小无量纲参数 η = ξṘ/D 为尺度的标度假设，作者推导出主导阶解，表明 R(t) ∝ √t，并证明界面温度偏离平衡转变温度 Tc 的程度与界面速度 Ṙ 成正比。

ABSTRACT

When a stable ordered phase and a metastable disordered phase are separated by a flat interface, the metastable state changes to the stable state through the propagation of the interface. For cases in which latent heat is generated, the interface displacement during some time interval is proportional to the square root of the time interval when the extent of supercooling is less than a certain value. We demonstrate this behavior by deriving a perturbative solution for a propagating interface in the phase field model. We calculate the leading-order contribution explicitly, and find that the interface temperature deviates from the equilibrium transition temperature in proportion to the interface velocity.

研究动机与目标

• 理论推导相场模型在 ∆ < 1 条件下数值模拟中观察到的 R(t) ∝ √t 时间依赖关系。
• 解决相场模型在超冷条件下传播界面缺乏解析解的问题。
• 建立在潜热存在下，界面速度与偏离 Tc 的温度偏差之间的关系。
• 提出一种基于标度坐标和小参数 η = ξṘ/D 的微扰框架，以系统计算主导阶解的修正项。

提出的方法

• 引入无量纲变量：θ = cp(T - Tc)/L，并利用 ξ 和 D 对相场模型方程 (8) 和 (9) 进行无量纲化。
• 定义标度空间坐标 w = (x - R(t))/ξ 和 z = Ṙ(x - R(t))/D，以捕捉界面和热扩散长度尺度。
• 引入一个小的无量纲参数 η = ξṘ/D，以实现关于 η 的正规微扰展开。
• 假设解具有如下标度形式：φ(x,t) = Φ(w; η)，θ(x,t) = Θ(z; η)，其中 R(t) 在主导阶下演化为 R(t) = √(2DCt)。
• 求解主导阶方程 (24) 和 (28)，得到 Φ0(w) 和 Θ0(z)，其中 Θ0(z) = -δ(z)，且 Φ0(w) 满足一个非线性常微分方程。
• 通过积分推导下一阶修正项 Φ1(w) 和 Θ1(z) 的可解性条件，获得 Θ1(0) 和 C1 的显式表达式。

实验结果

研究问题

• RQ1 能否在 ∆ < 1 条件下，从理论上推导出相场模型中数值模拟观察到的 R(t) ∝ √t 行为？
• RQ2 在潜热和有限界面速度存在下，界面温度如何偏离 Tc？
• RQ3 无量纲过冷参数 ∆ 在决定界面动力学和温度分布中起什么作用？
• RQ4 基于 η = ξṘ/D 的微扰框架如何实现对主导阶解修正项的系统计算？

主要发现

• 主导阶解给出 R(t) = √(2DCt)，确认在 ∆ < 1 条件下界面位置的 √t 标度关系。
• 界面温度 T(R(t), t) 偏离 Tc 的量与界面速度 Ṙ 成正比，其偏差表达式为 T(R(t), t) = Tc + (LξṘ/cpD)Θ1(0)。
• 控制温度分布的参数 C 在 ∆ → 1 时从下方发散，表明动力学中存在临界转变。
• 通过可解性条件显式计算了下一阶修正项 Θ1(0)，其值由涉及自由能导数的积分 I, J, K± 决定。
• 在 D¨R/Ṙ³ 展开中，系数 C1 被确定为零，确保了微扰展开的一致性。
• 对于 Π = 0（即 Ts = Tc）的情形，C(∆) 被计算并显示当 ∆ → 1 时发散，证实了在临界过冷阈值附近 √t 标度关系的失效。