摘要:研究了三元过冷熔体中柱状晶体在非等温条件下受straining流作用的生长问题, 给出了柱状晶体生长形态的近似解析表达式. 发现流入的straining流加快了界面的生长速度, 而流出的straining流减缓了界面的生长速度, 即straining流使得柱状晶体的界面发生变形. 同时发现, 随着流动速度的增大, 界面变形也更为显著. 通过比较straining流对纯熔体、二元熔体、三元熔体中柱状晶体界面的影响, 发现相比于纯熔体, 柱状晶体在稀合金熔体中的界面形态受straining流的影响更大.
关键词: 柱状晶体/
三元熔体/
straining流/
界面形态

English Abstract

全文HTML

--> --> -->

在材料科学领域, 凝固过程可以看作是传热、传质的过程, 它决定着材料晶体微结构的形态及生长, 进而决定着材料的性能. 自Mullins和Sekerka^[1]开创性地研究了熔体中球形粒子的生长问题以来, 许多研究者开始关注这一理论, 并陆续推广到不同几何形状粒子的生长问题^[2-12]. 柱状晶体作为一种重要的微结构, 其生长技术特别是单个晶体的生长技术对于提升半导体、光学等相关产品的性能有着重要的作用. Coriell和Parker^[13]研究了溶质扩散、表面张力各向异性及界面动力学对柱状冰界面形态稳定性的影响, 得出了柱状晶体生长的临界半径. 陈亚军等^[14]通过实验说明了柱状晶体受晶粒淘汰机制和合并机制的共同作用. Du等^[15]利用相场法研究了柱状晶体的生长问题, 指出生长前沿浓度的不均匀是产生不规则柱状晶体的原因. 因为流动可以促使凝固前沿浓度的不均匀, 当然也就可以通过流动来控制柱状晶体的界面形态. 基于上述思想一些研究者开始研究各种流动对柱状晶体生长界面形态的影响. 例如, Murakami等^[16]通过实验研究发现柱状晶体在Cu-Al二元熔体中生长时向来流方向偏转, 且其偏转角是流动大小和熔体浓度共同作用的结果. Szajnar^[17]研究了Al-Cu合金的凝固过程, 指出温度的变化不足以使柱状晶向等轴晶转变, 其在实验中通过强加磁场使液态金属产生流动, 从而使凝固界面的浓度发生变化.
在许多实际应用中一些研究者研究了多组分熔体的凝固问题. 例如, Altieri和Davis^[18]研究了多组分熔体的凝固问题, 通过线性稳定性分析得出了随着组分的增加凝固前端变得更加不稳定的结论. 由此可见需要研究多组分熔体的凝固问题, 并深入研究多组分熔体中不同形状晶体微结构的生长问题. Colin^[19]研究了柱状晶体在三元熔体中弹性效应对界面形态影响的生长问题, 发现了一个由纯弹性效应引起的不稳定性. 由于熔体中不可避免地存在流动, 这就需要研究存在流动的多组分熔体中柱状晶体的生长问题.
本文研究了straining流对多组分柱状晶体界面形态的影响, 利用渐近分析方法^[20]给出了柱状晶体界面形态的近似表达式. 分析了流动对柱状晶体界面形态的影响, 并比较了straining流对多组分熔体与单组分熔体柱状晶体界面形态影响的异同.

考虑单个初始半径为a的无限长柱状晶体在straining流作用下的多组分过冷稀合金熔体中的生长问题. 建立以柱状晶体的中心轴为z轴, 中心轴的中点为坐标原点O, 垂直于中心轴的横截面为$Oxy$平面的柱坐标系$(r, \theta, z)$. 流场的速度U是由远离柱状晶体的straining流${{U}}={c_1}x{{i}} + {c_2}y{{j}}$引起的, 其中${{i}}, {{j}}$为直角坐标系下的单位向量, $x, y$为直角坐标, ${c_1}, {c_2}$为常数且满足${c_1} + {c_2} = 0$, 即连续性条件$\nabla \cdot {{U}} = 0$成立. 用${T_{\rm{L}}}$和${T_{\rm{S}}}$分别表示液相和固相的温度, ${T_{\rm{I}}}$表示界面温度, $C_{\rm{S}}^i$表示固相第个i组分的浓度, $C_{\rm{L}}^i$表示液相第个i组分的浓度, 并做如下假设: 1)扩散的相互作用、浮力效应和熔体压力忽略不计; 2)液相和固相中的密度相等; 3)由于实际中固体扩散系数比液体扩散系数小几个数量级, 故而假设溶质在液相中的扩散速率远大于固相, 即$D_{\rm{L}}^i \gg D_{\rm{S}}^i$; 4) 由于柱状晶体相对于整个熔体来说较小, 故而根据Mullins和Sekerka的观点, 认定远离柱状晶体几个波长的熔体为远场; 5)假设远场温度为${T_\infty }$, 远场浓度为$C_\infty ^1$和$C_\infty ^2$且都为常数, 并视熔体为不可压缩的牛顿流体.
为了突出考虑流动对柱状晶体生长的影响, 采用如下尺度: 初始半径a为长度尺度; 特征速度$V=\dfrac{{{k_{\rm{L}}}\Delta T}}{{a\Delta H}}$为速度尺度, 其中${k_{\rm{L}}}$表示液相的热传导系数, $\Delta T = {T_{\rm{e}}} - {T_\infty }$, ${T_{\rm e}}$表示平界面的液相平衡温度, ${T_\infty }$表示远场温度, $\Delta H$表示熔体单位体积的潜热; $C_{\rm{e}}^i - C_{\rm{S}}^i$为浓度尺度, 其中$C_{\rm{e}}^i$表示第i个组分的液相平衡浓度, $C_{\rm{S}}^i$表示第i个杂质的固相浓度; $\Delta H/({c_{\rm{p}}}{\rho _{\rm{L}}})$为温度尺度; ${{{k_{\rm{L}}}V}}/({{{c_{\rm{p}}}a}})$为压力尺度, 其中${\rho _{\rm{L}}}$为熔体的密度, ${c_{\rm{p}}}$为熔体的定压比热. 引入无量纲量

$\begin{split}&{\bar T_{\rm{L}}} = \frac{{{T_{\rm{L}}} - {T_{\rm{M}}}}}{{\Delta H/({c_{\rm{p}}}{\rho _{\rm{L}}})}},\;\;{\bar T_{\rm{S}}} = \frac{{{T_{\rm{S}}} - {T_{\rm{M}} }}}{{\Delta H/({c_{\rm{p}}}{\rho _{\rm{L}}})}},\\&\bar C_{\rm{L}}^1 = \frac{{C_{\rm{L}}^1 - C_{\rm{e}}^1}}{{C_{\rm{e}}^1 - C_{\rm{S}}^1}},\;\;\bar C_{\rm{L}}^2 = \frac{{C_{\rm{L}}^2 - C_{\rm{e}}^2}}{{C_{\rm{e}}^2 - C_{\rm{S}}^2}},\\& \bar P = \frac{P}{{{k_{\rm{L}}}V/({c_{\rm{p}}}a)}},\;\;\bar {{U}} = \frac{{{U}}}{V},\;\;\bar r = \frac{r}{a},\;\;\bar t = \frac{t}{{{a / V}}}.\end{split}$

为方便起见, 省略无量纲头上的“–”, 得到无量纲控制方程

$\varepsilon \left(\frac{{\partial {T_{\rm{L}}}}}{{\partial t}} + {{U}} \cdot \nabla {T_{\rm{L}}}\right) = {\nabla ^2}{T_{\rm{L}}},$

$\varepsilon {\lambda _{\rm{T}}}\frac{{\partial {T_{\rm{S}}}}}{{\partial t}} = {\nabla ^2}{T_{\rm{S}}},$

$\varepsilon \left(\frac{{\partial C_{\rm{L}}^1}}{{\partial t}} + {{U}} \cdot \frac{{\partial C_{\rm{L}}^1}}{{\partial r}}\right) = \lambda _{\rm{D}}^1{\nabla ^2}C_{\rm{L}}^1,$

$\varepsilon \left(\frac{{\partial C_{\rm{L}}^2}}{{\partial t}} + {{U}} \cdot \frac{{\partial C_{\rm{L}}^2}}{{\partial r}}\right) = \lambda _{\rm{D}}^2{\nabla ^2}C_{\rm{L}}^2,$

$\varepsilon ({{U}} \cdot \nabla ){{U}} + \nabla P = Pr{\nabla ^2}{{U}},$

其中

$\begin{split}&\varepsilon = \frac{{\Delta T}}{{\Delta H/{c_{\rm{p}}}{\rho _{\rm{L}}}}},\;{\lambda _{\rm{T}}} = \frac{{{\kappa _{\rm{L}}}}}{{{\kappa _{\rm{S}}}}},\;\lambda _{\rm{D}}^1 = \frac{{D_{\rm{L}}^1}}{{{\kappa _{\rm{L}}}}},\\&\lambda _{\rm{D}}^2 = \frac{{D_{\rm{L}}^2}}{{{\kappa _{\rm{L}}}}},\;{\kappa _{\rm{L}}} = \frac{{{k_{\rm{L}}}}}{{{c_{\rm{p}}}{\rho _{\rm{L}}}}},\;Pr = \frac{\upsilon }{{{\kappa _{\rm{L}}}}},\end{split}$

式中的${\kappa _{\rm{L}}}$表示液相的热扩散系数, ${\kappa _{\rm{S}}}$表示固相的热扩散系数, $D_{\rm{L}}^1$和$D_{\rm{L}}^2$分别表示液相中第一个组分和第二个组分的扩散系数, $Pr$是普朗特数, $\upsilon $表示熔体的黏度, $\nabla $是梯度算子, ${\nabla ^2}$表示拉普拉斯算子. $\varepsilon $是无量纲参数, 其值等于Stefan数的倒数. 无量纲参数${\lambda _{\rm{T}}}$表示液相和固相热扩散系数的比值.
晶体生长控制方程(2)—(6)在界面$r = R(\theta, t)$上满足温度连续性条件、Gibbs-Thomson条件、热量守恒条件、质量守恒条件、界面无滑移条件.

${T_{\rm{L}}} = {T_{\rm{S}}} = {T_{\rm{I}}},$

$\begin{split}{T_{\rm{I}}} =\; & \varepsilon \varGamma 2K - \varepsilon {E^{ - 1}}{M_{\rm{K}}}{{{U}}_{\rm{I}}} - M_{\rm{C}}^1C_{\rm{L}}^1 - \varepsilon M_{\rm{e}}^1 \\ &- M_{\rm{C}}^2C_{\rm{L}}^2 - \varepsilon M_{\rm{e}}^2,\end{split}$

$\varepsilon {{{U}}_{\rm{I}}} = ({K_{\rm{T}}}\nabla {T_{\rm{S}}} - \nabla {T_{\rm{L}}}) \cdot {{n}},$

$\varepsilon S_{\rm{p}}^1{{{U}}_{\rm{I}}} = - \nabla C_{\rm{L}}^1{{n}},$

$\varepsilon S_{\rm{p}}^2{{{U}}_{\rm{I}}} = - \nabla C_{\rm{L}}^{\rm{2}}{{n}},$

${{U}} \cdot {{n}} = 0,\;{{U}} \cdot \tau = 0,$

其中

$ \begin{split} & \varGamma = \frac{{\gamma \Delta {T_{\rm{M}} }}}{{a\Delta H\Delta T}},\;\;E = \frac{{\Delta T}}{{{T_{\rm{M}}}}},\;\;{M_{\rm{K}}} = \frac{V}{{\mu {T_{\rm{M}}}}},\\ &{K_{\rm{T}}} = \frac{{{k_{\rm{S}}}}}{{{k_{\rm{L}}}}},\;\;S_{\rm{p}}^1 = \frac{{({k_1} - 1)(1 - C_0^1)}}{{{k_1}\lambda _{\rm{D}} ^1}},\\ & S_{\rm{p}}^2 = \frac{{\left( {{k_2} - 1} \right)\left( {1 - C_0^2} \right)}}{{{k_2}\lambda _{\rm{D}}^2}},\;\;C_0^1 = \frac{{C_{\rm{e}}^1}}{{C_{\rm{e}}^1 - C_{\rm{S}}^1}},\\ &C_0^2 = \frac{{C_{\rm{e}}^2}}{{C_{\rm{e}}^2 - C_{\rm{S}}^2}},\;\;M_{\rm{C}}^1 = - \frac{{{m_1}(C_{\rm{e}}^1 - C_{\rm{S}}^1)}}{{\Delta H/{c_{\rm{p}}}{\rho _{\rm{L}}}}},\\ & M_{\rm{C}}^2 = - \frac{{{m_2}(C_{\rm{e}}^2 - C_{\rm{S}}^2)}}{{\Delta H/{c_{\rm{p}}}{\rho _{\rm{L}}}}},\;\;M_{\rm{e}}^1 = - \frac{{{m_1}(C_{\rm{e}}^1 - C_\infty ^1)}}{{\Delta T}},\\ &M_{\rm{e}}^2 = - \frac{{{m_2}(C_{\rm{e}}^2 - C_\infty ^2)}}{{\Delta T}}, \end{split} $

其中K为界面的局部平均曲率, ${{{U}}_{\rm{I}}}$表示界面生长速度, $\gamma $为表面自由能, $\mu $为界面动力学系数, ${k_{\rm{S}}}$和${k_{\rm{L}}}$分别表示固相和液相的热传导系数, n表示单位外法线向量, ${{\tau }}$表示单位切向向量.
晶体控制方程(2)—(6)还应满足下面的远场条件、初始条件、界面分离条件.
远场条件, 当$r \to \infty $时, 有

$\begin{split}& {T_{\rm{L}}} \to - \varepsilon,\;\;C_{\rm{L}}^1 \to - C_{{\rm{L}},\infty }^1\varepsilon,\\ & C_{\rm{L}}^2 \to - C_{{\rm{L}},\infty }^2\varepsilon,\;\;{{U}} \to Ax{{i}} + By{{j}},\end{split}$

其中

$C_{{\rm{L}},\infty }^1 = \frac{{(C_{\rm{e}}^1 - C_\infty ^1)\Delta H}}{{(C_{\rm{e}}^1 - C_{\rm{S}}^1){c_{\rm{p}}}{\rho _{\rm{L}}}\Delta T}},$

$C_{{\rm{L}},\infty }^2 = \frac{{(C_{\rm{e}}^2 - C_\infty ^2)\Delta H}}{{(C_{\rm{e}}^2 - C_{\rm{S}}^2){c_{\rm{p}}}{\rho _{\rm{L}}}\Delta T}},$

$A = {{{c_1}a}}/{V},~~ B = {{{c_2}a}}/{V},$

且有$A + B = 0$成立.
初始条件

$R(\theta,0) = 1.$

界面分离条件

$C_{\rm{S}}^i = {k_i}C_{\rm{L}}^i,$

其中${k_i}$表示液相第i个组分的分离系数.
由于界面附近的热传递、质量传递与远离界面处的情况有较大的区别, 因而界面附近的区域和远离界面的熔体区有不同的尺度. 故引入慢变量$\overline r = \varepsilon r,$ 把$\overline r,$ $r, $ $\theta, $ t视为独立的变量, 则(2)—(16)式变为

$\begin{split} &\varepsilon \frac{{\partial {T_{\rm{L}}}}}{{\partial t}}+\varepsilon ({{U}} \cdot \nabla ){T_{\rm{L}}} + {\varepsilon ^2}{u_r}\frac{{\partial {T_{\rm{L}}}}}{{\partial \bar r}} \\=\;& {\nabla ^2}{T_{\rm{L}}}+2\varepsilon \frac{{{\partial ^2}{T_{\rm{L}}}}}{{\partial r\partial \bar r}}+{\varepsilon ^2}\frac{{{\partial ^2}{T_{\rm{L}}}}}{{\partial {{\bar r}^2}}}+\frac{\varepsilon }{r}\frac{{\partial {T_{\rm{L}}}}}{{\partial \bar r}},\end{split}$

$\varepsilon {\lambda _{\rm{T}}}\frac{{\partial {T_{\rm{S}}}}}{{\partial t}} = {\nabla ^2}{T_{\rm{S}}}+2\varepsilon \frac{{{\partial ^2}{T_{\rm{S}}}}}{{\partial r\partial \bar r}}+{\varepsilon ^2}\frac{{{\partial ^2}{T_{\rm{S}}}}}{{\partial {{\bar r}^2}}}+\frac{\varepsilon }{r}\frac{{\partial {T_{\rm{S}}}}}{{\partial \bar r}},$

$\begin{split} &\varepsilon \frac{{\partial C_L^1}}{{\partial t}}+\varepsilon ({{U}} \cdot \nabla )C_{\rm{L}}^1 + {\varepsilon ^2}{u_r}\frac{{\partial C_{\rm{L}}^1}}{{\partial t}} \\ =\; & \lambda _{\rm{D}}^1\left({\nabla ^2}C_{\rm{L}}^1 \!+\! 2\varepsilon \frac{{{\partial ^2}C_{\rm{L}}^1}}{{\partial r\partial \bar r}}\!+\!{\varepsilon ^2}\frac{{{\partial ^2}C_{\rm{L}}^1}}{{\partial {{\bar r}^2}}}\!+\!\frac{\varepsilon }{r}\frac{{\partial C_{\rm{L}}^1}}{{\partial \bar r}}\right),\end{split}$

$\begin{split} & \varepsilon \frac{{\partial C_L^2}}{{\partial t}}+\varepsilon ({{U}} \cdot \nabla )C_{\rm{L}}^2 + {\varepsilon ^2}{u_r}\frac{{\partial C_{\rm{L}}^2}}{{\partial t}} \\=\, & \lambda _{\rm{D}}^2\left({\nabla ^2}C_{\rm{L}}^2+2\varepsilon \frac{{{\partial ^2}C_{\rm{L}}^2}}{{\partial r\partial \bar r}}\!+\! {\varepsilon ^2}\frac{{{\partial ^2}C_{\rm{L}}^2}}{{\partial {{\bar r}^2}}}\!+\! \frac{\varepsilon }{r}\frac{{\partial C_{\rm{L}}^2}}{{\partial \bar r}}\right),\end{split}$

$\nabla {{U}}+\varepsilon \frac{{\partial {u_r}}}{{\partial \bar r}} = 0,$

$\begin{split} &\varepsilon ({{U}} \cdot \nabla ){{U}}+{\varepsilon ^2}{u_r}\frac{{\partial {{U}}}}{{\partial \bar r}} \\=\; & \Pr \left({\nabla ^2}{{U}} + 2\varepsilon \frac{{{\partial ^2}{{U}}}}{{\partial r\partial \bar r}} + {\varepsilon ^2}\frac{{{\partial ^2}{{U}}}}{{\partial {{\bar r}^2}}} + \frac{\varepsilon }{r}\frac{{\partial {{U}}}}{{\partial \bar r}}\right),\end{split}$

其中${{U}} = ({u_r}, {u_\theta })$.
界面条件为

$\begin{split} {T_{\rm{I}}} =\;& {T_{\rm{L}}} = {T_{\rm{S}}} = 2\varepsilon \varGamma K - \varepsilon {E^{ - 1}}{M_{\rm{K}}}{{{U}}_{\rm{I}}} - M_{\rm{C}}^1C_{\rm{L}}^1 \\ &- M_{\rm{C}}^2C_{\rm{L}}^2 - \varepsilon M_{\rm{e}}^1 - \varepsilon M_{\rm{e}}^2,\end{split}$

$\varepsilon {{{U}}_{\rm{I}}} = ({K_{\rm{T}}}\nabla {T_{\rm{S}}} - \nabla {T_{\rm{L}}}) \cdot {{n}} + \varepsilon \frac{\partial }{{\partial \bar r}}({K_{\rm{T}}}{T_{\rm{S}}} - {T_{\rm{L}}}),$

$\varepsilon S_{\rm{p}}^1{{{U}}_I} = - \nabla C_{\rm{L}}^1 \cdot {{n}} - \varepsilon \frac{\partial }{{\partial \bar r}}C_{\rm{L}}^1,$

$\varepsilon S_{\rm{p}}^2{{{U}}_{\rm{I}}} = - \nabla C_{\rm{L}}^2 \cdot {{n}} - \varepsilon \frac{\partial }{{\partial \bar r}}C_{\rm{L}}^2.$

远场条件, 当$r \to \infty $时, 有

$\begin{split}& {T_{\rm{L}}} \to - \varepsilon,\;\;C_{\rm{L}}^1 \to - C_{{\rm{L}},\infty }^1\varepsilon,\;\;C_{\rm{L}}^2 \to - C_{{\rm{L}},\infty }^2\varepsilon,\\& {u_r} \to r(A{\cos ^2}\theta + B{\sin ^2}\theta ),\\&{u_\theta } \to r(B - A)\sin \theta \cos \theta.\end{split}$

初始条件为

$R(\theta,0) = 1.$

对于典型的金属, $\Delta H/({c_{\rm{p}}}{\rho _{\rm{L}}})$通常为几百K. 如 Cu的$\Delta H=1.830 \times {10^9}$ J·m^–3, ${c_{\rm{p}}}$ = 390 J·kg^–1·K^–1, ${\rho _{\rm{L}}}$ = 8930 kg·m^–3, $\Delta H/({c_{\rm{p}}}{\rho _{\rm{L}}})$ = 525 K. Al的$\Delta H=1.{\rm{0676}} \times {10^9}$ J·m^–3, ${c_{\rm{p}}}$ = 1084 J·kg^–1·K^–1, ${\rho _{\rm{L}}}$ = 2700 kg·m^–3, $\Delta H/({c_p}{\rho _L})$ = 365 K, 因此无量纲参数$\varepsilon $实际上是一个小量, 即

$\varepsilon =\frac{{\Delta T}}{{\Delta H/{C_{\rm{P}}}{\rho _{\rm{L}}}}} \ll 1.$

当$\varepsilon \to 0$时, 采用类似于Vogel和Cantor^[20]的方法, 求解柱状晶体生长系统(17)式—(28)式的渐近解. 在界面附近的区域, 做如下展开:

$\begin{split}& {T_{\rm{L}}} = \varepsilon {T_{{\rm{L}}0}} + {\varepsilon ^2}{T_{{\rm{L}}1}} + \cdots,\quad {T_{\rm{S}}} = \varepsilon {T_{{\rm{S}}0}} + {\varepsilon ^2}{T_{{\rm{S}}1}} + \cdots, \\& C_{\rm{L}}^1 = \varepsilon C_{{\rm{L}}0}^1 + {\varepsilon ^2}C_{{\rm{L}}1}^1 + \cdots,~~ C_{\rm{L}}^2 \!=\! \varepsilon C_{{\rm{L}}0}^2 \!+\! {\varepsilon ^2}C_{{\rm{L}}1}^2 + \cdots,\\& R = {R_0} + \varepsilon {R_1} + \cdots, \quad \bar R = {{\bar R}_0} + \varepsilon {{\bar R}_1} + \cdots, \\[-12pt]\end{split}$

其中$\bar R = \varepsilon R$, ${\bar R_0} = \varepsilon {R_0}$, ${\bar R_1} = \varepsilon {R_1}.$
界面平均曲率可以展开为

$K = - \frac{1}{{{R_0}}} + \frac{\varepsilon }{{R_0^2}}\left(\frac{{{\partial ^2}}}{{\partial {\theta ^2}}} + 1\right){R_1} + \cdots .$

假设${\partial }/{{\partial t}} = O(1)$, 把(29)式—(31)式代入(17)式—(28)式可得出首阶近似.
首阶解的控制方程为

$\begin{split}&{\nabla ^2}{T_{{\rm{L}}0}}=0,\;\;{\nabla ^2}{T_{{\rm{S}}0}}=0,\;\;{\nabla ^2}C_{{\rm{L}}0}^1=0,\\ &{\nabla ^2}C_{{\rm{L}}0}^2=0,\;\;\nabla {{{U}}_0}=0,\;\;{\nabla ^2}{{{U}}_0}=0.\end{split}$

在柱晶界面$R = {R_0}$上, 有

$\begin{split}{T_{{\rm{L}}0}} = {T_{{\rm{S0}}}} =\;& - \frac{{2\varGamma }}{{{R_0}}} - {E^{ - 1}}{M_{\rm{K}}}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}} - M_{\rm{C}}^1C_{L0}^1 \\ &- M_{\rm{C}}^2C_{{\rm{L}}0}^2 - M_{\rm{e}}^1 - M_{\rm{e}}^2,\end{split}$

$\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}} = {K_{\rm{T}}}\frac{{\partial {T_{{\rm{S}}0}}}}{{\partial r}} - \frac{{\partial {T_{{\rm{L}}0}}}}{{\partial r}},$

$S_{\rm{p}}^1\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}} = - \frac{{\partial C_{{\rm{L0}}}^1}}{{\partial r}},$

$S_{\rm{p}}^2\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}} = - \frac{{\partial C_{{\rm{L0}}}^2}}{{\partial r}},$

${{{U}}_0}=0,$

远场条件: 当$r \to \infty $时

$\begin{split}& {T_{{\rm{L}}0}} \to - 1,\;\;C_{{\rm{L}}0}^1 \to - C_{{\rm{L}},\infty }^1,\;\;C_{{\rm{L}}0}^2 \to - C_{{\rm{L}},\infty }^2,\\&{u_r}=r(A{\cos ^2}\theta + B{\sin ^2}\theta ),\\&{u_\theta } \to r(B - A)\sin \theta \cos \theta.\end{split}$

初始条件:

$R(\theta,0) = 1.$

因为${{U}} = ({u_r}, {u_\theta })$满足方程(32), 故有下列方程成立:

$\frac{{\partial {u_r}}}{{\partial r}} + \frac{{{u_r}}}{r} + \frac{1}{r}\frac{{\partial {u_\theta }}}{{\partial \theta }} = 0,$

$\frac{1}{r}\frac{{\partial {u_r}}}{{\partial r}} + \frac{{{\partial ^2}{u_r}}}{{\partial {r^2}}} + \frac{1}{{{r^2}}}\frac{{{\partial ^2}{u_r}}}{{\partial {\theta ^2}}} - \frac{2}{{{r^2}}}\frac{{\partial {u_\theta }}}{{\partial \theta }} - \frac{1}{{{r^2}}}{u_r} = 0,$

$\frac{1}{r}\frac{{\partial {u_\theta }}}{{\partial r}} + \frac{{{\partial ^2}{u_\theta }}}{{\partial {r^2}}} + \frac{1}{{{r^2}}}\frac{{{\partial ^2}{u_\theta }}}{{\partial {\theta ^2}}}+\frac{2}{{{r^2}}}\frac{{\partial {u_r}}}{{\partial \theta }} - \frac{1}{{{r^2}}}{u_\theta } = 0,$

由远场条件(38)式可知, 当$r \to \infty $时有

${u_r} \to Ar\cos (2\theta ),\;\;{u_\theta } \to - Ar\sin (2\theta ).$

由(40)式—(43)式解得

$\begin{split}& {u_r} = \left(r - \frac{{R_0^4}}{{{r^3}}}\right)A\cos (2\theta ), \\& {u_\theta } = \left( - r - \frac{{R_0^4}}{{{r^3}}}\right)A\sin (2\theta ),\end{split}$

类似于文献[21], 可求得温度场和浓度场的首阶近似解为:

${T_{{\rm{L}}0}}(t) \!=\! {R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{\rm{(}}{A_\lambda }+\ln {R_0} - \ln r){{\rm{e}}^{{{\overline R}_0} - \overline r}} \!-\! 1,$

${T_{{\rm{S0}}}}={A_\lambda }{R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}} - 1,$

$C_{{\rm{L0}}}^{\rm{1}}(t) = S_{\rm{p}}^{\rm{1}}{R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}({A_\lambda } + \ln {R_0} - \ln r){{\rm{e}}^{{{\bar R}_0} - \bar r}} - C_{{\rm{L}},\infty }^1,$

$C_{{\rm{L0}}}^2(t) = S_{\rm{p}}^2{R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}({A_\lambda } + \ln {R_0} - \ln r){{\rm{e}}^{{{\bar R}_0} - \bar r}} - C_{{\rm{L}},\infty }^2,$

其中${A_\lambda } = \ln {R_\infty } - \ln {R_0} = {N}/({{2{\lambda ^2}}}) > 0,$ $\lambda $是方程${\lambda ^2}\ln ({\nu ^2}{\lambda ^2}) \!+\! N \!=\! 0$的一个解. $N \!\ll\! 1,$ $\ln \nu^2 \!=\! 0.5772$.
由界面条件可得柱状晶体半径${R_0}$的首阶近似解满足下列微分方程:

$\begin{split} & {{{\rm{d}}{R_0}}}/{{{\rm{d}}t}}= \\ \; &\frac{{{R_0} - 2\varGamma }}{{{R_0}({A_\lambda }{R_0} + M_{\rm{C}}^1S_{\rm{p}}^1{R_0}{A_\lambda } + M_{\rm{C}}^2S_{\rm{p}}^2{R_0}{A_\lambda } + {E^{ - 1}}{M_{\rm{K}}})}},\\ &\qquad {R_0}(0) = 1. \\[-12pt]\end{split}$

因此有

$\begin{split}t =\;& [4{\varGamma ^2}{A_\lambda }(M_{\rm{C}}^1S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2 + 1) + 2\varGamma {M_{\rm{K}}}{E^{ - 1}}]\\& \times \ln [({{{R_0} - 2\varGamma }})/({{1 - 2\varGamma }})]\\&+ [2\varGamma {A_\lambda }(M_{\rm{C}}^1S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2 + 1) + {E^{ - 1}}{M_{\rm{K}}}]({R_0} - 1)\\ &+ \frac{{{A_\lambda }}}{2}(M_{\rm{C}}^1S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2 + 1)(R_0^2 - 1).\\[-15pt]\end{split}$

由(48)式可知, 当$\dfrac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}} < 0$时, 晶体收缩; 当$\dfrac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}} > 0$时, 晶体生长.
一阶渐近解的控制方程为

$\frac{{\partial {T_{{\rm{L}}0}}}}{{\partial t}} + ({{{U}}_0} \cdot \nabla ){T_{{\rm{L}}0}} = {\nabla ^2}{T_{{\rm{L}}1}} + 2\frac{{{\partial ^2}{T_{{\rm{L}}0}}}}{{\partial r\partial \bar r}} + \frac{1}{r}\frac{{\partial {T_{{\rm{L}}0}}}}{{\partial \bar r}},$

${\lambda _{\rm{S}}}\frac{{\partial {T_{{\rm{S}}0}}}}{{\partial t}} = {\nabla ^2}{T_{{\rm{S}}1}} + 2\frac{{{\partial ^2}{T_{{\rm{S}}0}}}}{{\partial r\partial \bar r}} + \frac{1}{r}\frac{{\partial {T_{{\rm{S}}0}}}}{{\partial \bar r}},$

$\begin{split}& \frac{{\partial C_{{\rm{L}}0}^1}}{{\partial t}} + ({{{U}}_0} \cdot \nabla )C_{{\rm{L}}0}^1 \\={}& \lambda _{\rm{D}}^1\left({\nabla ^2}C_{{\rm{L}}1}^1 + 2\frac{{{\partial ^2}C_{{\rm{L}}0}^1}}{{\partial r\partial \bar r}} + \frac{1}{r}\frac{{\partial C_{{\rm{L}}0}^1}}{{\partial \bar r}}\right),\end{split}$

$\begin{split}& \frac{{\partial C_{{\rm{L}}0}^2}}{{\partial t}} + ({{{U}}_0} \cdot \nabla )C_{{\rm{L}}0}^2 \\={}& \lambda _{\rm{D}}^2\left({\nabla ^2}C_{{\rm{L}}1}^2 + 2\frac{{{\partial ^2}C_{{\rm{L}}0}^2}}{{\partial r\partial \bar r}} + \frac{1}{r}\frac{{\partial C_{{\rm{L}}0}^2}}{{\partial \bar r}}\right),\end{split}$

其满足下面的界面条件:

${T_{{\rm{L1}}}} = {T_{{\rm{S1}}}}+\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{R_1}+{A_\lambda }{R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{\bar R_1},$

$ \begin{split} {T_{{\rm{S1}}}} =\;& \frac{{2\varGamma }}{{R_0^2}}\left(\frac{{{\partial ^2}}}{{\partial {\theta ^2}}} + 1\right){R_1} - {E^{ - 1}}{M_{\rm{K}}}\frac{{{\rm{d}}{R_1}}}{{{\rm{d}}t}} - M_{\rm{C}}^{\rm{1}}C_{{\rm{L1}}}^{\rm{1}} \\ & - M_{\rm{C}}^{\rm{2}}C_{{\rm{L}}1}^2 + {\rm{(}}M_{\rm{C}}^{\rm{1}}S_{\rm{p}}^{\rm{1}} + M_{\rm{C}}^2S_{\rm{p}}^{\rm{2}}{\rm{)}}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{R_1}\\ &+ (M_{\rm{C}}^1S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2){R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{\bar R_1},\\[-15pt] \end{split} $

$\begin{split}\frac{{{\rm{d}}{R_1}}}{{{\rm{d}}t}} =\; & {K_T}\frac{{\partial {T_{{\rm{S}}1}}}}{{\partial r}} - \frac{{\partial {T_{{\rm{L1}}}}}}{{\partial r}} - \frac{{{R_1}}}{{{R_0}}}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}} \\ & + {A_\lambda }{R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}} - {A_\lambda }{R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{\bar R_1},\end{split}$

$\begin{split} S_{\rm{p}}^1\frac{{{\rm{d}}{R_1}}}{{{\rm{d}}t}} =\; & - \frac{{\partial C_{{\rm{L1}}}^1}}{{\partial r}} - \frac{{S_{\rm{p}}^1}}{{{R_0}}}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{R_1}\\ & - S_{\rm{P}}^1\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{\bar R_1} + S_{\rm{p}}^1{R_0}{A_\lambda }\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}},\end{split}$

$\begin{split}S_{\rm{p}}^2\frac{{{\rm{d}}{R_1}}}{{{\rm{d}}t}} =\; & - \frac{{\partial C_{{\rm{L1}}}^2}}{{\partial r}} - \frac{{S_{\rm{p}}^2}}{{{R_0}}}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{R_1} \\ &- S_{\rm{p}}^2\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{\bar R_1} + S_{\rm{p}}^2{R_0}{A_\lambda }\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}.\end{split}$

远场条件: 当$r \to \infty $时, 有

${T_{{\rm{L1}}}} \to 0,\;\;C_{{\rm{L1}}}^{\rm{1}} \to 0,\;\;C_{{\rm{L1}}}^{\rm{2}} \to 0,$

初始条件

${R_1}(\theta,0) = 0.$

类似于文献[21], 可求得一阶近似解为

${T_{{\rm{L}}1}} = {C_{{\rm{L}}0}} + {C_{{\rm{L}}2}}{r^{ - 2}}\cos (2\theta ) + T_{{\rm{L}}10}^* + T_{{\rm{L}}12}^*\cos (2\theta ),$

${T_{{\rm{S}}1}} = {C_{{\rm{S0}}}} + {C_{{\rm{S2}}}}{r^2}\cos (2\theta ) + T_{{\rm{S1}}0}^*,$

$C_{{\rm{L}}1}^1 = C_{{\rm{C0}}}^1 + C_{{\rm{C2}}}^{\rm{1}}{r^{ - 2}}\cos (2\theta ) + C_{{\rm{C}}10}^{1*} + C_{{\rm{C12}}}^{{\rm{1*}}}\cos (2\theta ),$

$C_{{\rm{L}}1}^2 = C_{{\rm{C}}0}^2 + C_{{\rm{C}}2}^2{r^{ - 2}}\cos (2\theta ) + C_{{\rm{C10}}}^{{\rm{2*}}} + C_{{\rm{C}}12}^{2*}\cos (2\theta ),$

${R_1} = {g_0} + {g_2}\cos (2\theta ),$

其中参数${C_{{\rm{L}}0}}, $ ${C_{{\rm{L}}2}}, $ $T_{{\rm{L}}10}^*, $ $T_{{\rm{L}}12}^*, $ ${C_{{\rm{S}}0}}, $ ${C_{{\rm{S2}}}}{\rm{, }}$ $T_{{\rm{S10}}}^{\rm{*}}, $ $C_{{\rm{C}}0}^1, $ $C_{{\rm{C}}2}^1, $ $C_{{\rm{C}}10}^{1*}, $ $C_{{\rm{C}}12}^{1*}, $ $C_{{\rm{C}}0}^2, $ $C_{{\rm{C}}2}^2, $ $C_{{\rm{C}}10}^{2*}, $ $C_{{\rm{C}}12}^{2*}, $ $\dfrac{{{\rm{d}}{g_0}}}{{{\rm{d}}t}}, $ $\dfrac{{{\rm{d}}{g_2}}}{{{\rm{d}}t}}, $ $\dfrac{{{\rm{d}}{b_0}}}{{{\rm{d}}t}}$的取值见附录.
综上可得温度场、浓度场及界面半径的渐近解为:

$\begin{split}{T_{\rm{L}}} =\; & \varepsilon {T_{{\rm{L}}0}} + {\varepsilon ^2}[{C_{{\rm{L}}0}} + {C_{{\rm{L}}2}}{r^{ - 2}}\cos (2\theta ) \\ &+ T_{{\rm{L}}10}^* + T_{{\rm{L}}12}^*\cos (2\theta )] + \cdots,\end{split}$

${T_{\rm{S}}} = \varepsilon {T_{{\rm{S}}0}} + {\varepsilon ^2}[{C_{{\rm{S}}0}} + {C_{{\rm{S}}2}}{r^2}\cos (2\theta ) + T_{{\rm{S}}10}^*] + \cdots,$

$\begin{split} C_{\rm{L}}^1 =\; & \varepsilon C_{{\rm{L}}0}^1 + {\varepsilon ^2}[C_{{\rm{C}}0}^1 + C_{{\rm{C2}}}^{\rm{1}}{r^{ - 2}}\cos (2\theta ) \\ &+ C_{{\rm{C}}10}^{1*} + C_{{\rm{C}}12}^{1*}\cos (2\theta )],\end{split}$

$\begin{split} C_{\rm{L}}^2 =\; & \varepsilon C_{{\rm{L}}0}^2 + {\varepsilon ^2}[C_{{\rm{C}}0}^2 + C_{{\rm{C}}2}^2{r^{ - 2}}\cos (2\theta )\\ & + C_{{\rm{C}}10}^{2*} + C_{{\rm{C}}12}^{2*}\cos (2\theta )],\end{split}$

$R = {R_0} + \varepsilon [{g_0} + {g_2}\cos (2\theta )] + \cdots .$

界面生长速率是

${{{U}}_I} = \frac{{\partial R}}{{\partial t}} = \frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}} + \varepsilon \left[ {\frac{{{\rm{d}}{g_0}}}{{{\rm{d}}t}} + \frac{{{\rm{d}}{g_2}}}{{{\rm{d}}t}}\cos (2\theta )} \right] + \cdots .$

利用渐近方法研究了非等温条件下多组分稀合金熔体中柱状晶体生长的数学模型, 获得了柱状晶体生长模型界面的近似解析解. 由所得出的解析解表达式(71)式, 可以定量分析straining流对柱状晶体界面形态的影响. 由图1可以直观地发现在straining流的作用下, 柱状晶体的界面形态发生了明显的变化, 即柱状晶体生长的横截面不再保持圆形. 图2给出了柱状晶体在straining流下生长的界面形态. 图3给出了不同强度的straining流对柱状晶体界面形态的影响, 可以看出流入的straining流加强了界面的生长速度, 流出的straining流降低了界面的生长速度, 由此产生了不规则的柱状晶体. 这从物理含义方面理解实际上是由于流动改变了界面前沿的浓度, 导致晶体生长前沿浓度的不均匀, 从而产生了不规则的柱状晶体. 这一结论与文献[15]的结论是一致的, 但需要指出的是, 我们给出的解析解可以定量地描述这一问题, 从而为将来通过流动来精确控制柱状晶体的界面形态打下良好的基础. 图4给出了在固定强度straining流影响下, 柱状晶体界面形态随时间变化的生长过程, 这一发现的物理意义在于可以动态地预测某一时刻柱状晶体界面的形态. 图5比较了straining流对多组分熔体和纯熔体中柱状晶体界面形态影响, 发现多组分稀合金熔体中的柱状晶体界面形态受straining流的影响更大. 文献[18]在研究多组分熔体的凝固问题时指出, 随着组分的增加凝固前端变得更加不稳定, 这与我们的结论是一致的. 图6比较了straining流对二元熔体和三元熔体中柱状晶体界面形态影响, 发现熔体中组分的多少, 并不是柱状晶体界面形态变化的决定性因素, 即不能说明三元柱状晶体的界面形态比二元柱状晶体的界面形态更容易受straining流的影响. 由图7柱状晶体界面杂质浓度$C_{\rm{L}}^1$随$\theta $的变化情况, 可以看出柱状晶体的界面形态主要是由界面附近的溶质浓度决定的. 实际上通过(71)式可以发现, 相比于纯熔体中柱状晶体的生长, 在多组分熔体中生长的柱状晶体的界面形态表达式更为复杂, 其中参数$M_{\rm{C}}^1 S_{\rm{p}}^1 + M_{\rm{C}}^2 S_{\rm{p}}^2$和$\dfrac{1}{{\lambda _{\rm{D}}^1}}M_{\rm{C}}^1 S_{\rm{p}}^1 + \dfrac{1}{{\lambda _{\rm{D}}^2}}M_{\rm{C}}^2 S_{\rm{p}}^2$对界面形态有着重要的影响. 为了便于分析我们假设$D_{\rm{L}}^1 = D_{\rm{L}}^2$, 即液相中第一个组分和第二个组分的扩散系数相等, 容易得出$\lambda _{\rm{D}}^1 = \lambda _{\rm{D}}^2$. 这时只需分析$M_{\rm{C}}^1 S_{\rm{p}}^1 + M_{\rm{C}}^2 S_{\rm{p}}^2$对于界面形态的影响即可. 由$M_{\rm{C}}^1 S_{\rm{p}}^1 + M_{\rm{C}}^2 S_{\rm{p}}^2$容易看出相比于纯熔体, 多组分熔体中柱状晶体的界面形态还主要受液相线斜率${m_1}, {m_2}$以及液相中第i个杂质组分浓度$C_{\rm{L}}^i$的影响, 实际上这正是受到成分过冷(constitutional undercooling)影响的结果. 通过上面的分析可以得到以下结论, 稀合金熔体中柱状晶体的界面形态比纯熔体中的柱状晶体的界面形态受straining流的影响更为显著.
图 1 在$Oxy$平面上应变流对柱状晶体形态演化的影响, 其中$t = 396, $ $\varGamma = 0.25, $ $M_{\rm{C}}^1=0.01, $ $M_{\rm{C}}^2=0.02, $ $C_{{\rm{L}}, \infty }^1 = $1.0, $C_{{\rm{L}}, \infty }^2 = 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ ${m_2} = $ –2.33, $E = 0.3, $ ${M_{\rm{k}}} = 0.01, $ ${k_{\rm{T}}} = 1.23, $ ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^1 = $ 0.01, $\lambda _{\rm{D}}^2 = 0.02, $ $\varepsilon = 0.05$
Figure1. The morphology evolution of columnar crystal in a straining flow on the cross-section of $Oxy$plane at $t = 396, $ where $\varGamma = 0.25, $ $M_{\rm{C}}^1=0.01, $ $M_{\rm{C}}^2=0.02, $ $C_{{\rm{L}}, \infty }^1 = 1.0, $ $C_{{\rm{L}}, \infty }^2 = 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ ${m_2} = - 2.33, $ E = 0.3, ${M_{\rm{k}}} = 0.01, $ ${k_{\rm{T}}} = 1.23, $ ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^1 = 0.01, $ $\lambda _{\rm{D}}^2 = $ 0.02, $ \varepsilon = 0.05$

图 2 在$t = 396$时, 柱状晶体的界面形态. 其中$\varGamma = 0.25, $ $M_{\rm{C}}^{\rm{1}}=0.01, $ $M_{\rm{C}}^{\rm{2}}=0.02, $ $C_{{\rm{L}}, \infty }^1 = 1.0, $ $C_{{\rm{L}}, \infty }^2 = 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ ${m_2} = - 2.33, $ $E = 0.3, $ ${M_{\rm{k}}} = 0.01, $ ${k_{\rm{T}}} = 1.23, $ ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^1 = 0.01, $ $\lambda _{\rm{D}}^{\rm{2}} = 0.02, $ $\varepsilon = 0.05, $ $A = 0.9$
Figure2. The morphology evolution of columnar crystal in a straining flow at $t = 396, $ where $\varGamma = 0.25, $ $M_{\rm{C}}^{\rm{1}}=0.01, $ $M_{\rm{C}}^{\rm{2}}=0.02, $ $C_{{\rm{L}}, \infty }^1 = 1.0, $ $C_{{\rm{L}}, \infty }^2 = 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ ${m_2} = - 2.33, $ $E = 0.3, $ ${M_{\rm{k}}} = 0.01, $ ${k_{\rm{T}}} = 1.23, $ ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^1 = 0.01, $ $\lambda _{\rm{D}}^{\rm{2}} = 0.02, $ $\varepsilon = 0.05, $ $A = 0.9$

图 3 不同强度的应变流对柱状晶体界面形态的影响, 其中$t \!=\! 256, $ $\varGamma \!=\! 0.25, $ $M_{\rm{C}}^1 \!=\! 0.01, $ $M_{\rm{C}}^2\! =\! 0.02, $ $C_{{\rm{L}}, \infty }^1 = 1.0, $ $C_{{\rm{L}}, \infty }^2 = 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ ${m_2} = - 2.33, $ $E = 0.3, $ ${M_{\rm{k}}} = 0.01, $ ${k_{\rm{T}}} = 1.23, $ ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^1 = 0.01, $ $\lambda _{\rm{D}}^{\rm{2}} = 0.02, $ $\varepsilon = 0.05.$ A由左向右分别为0.9, 0.6, 0.3, 0
Figure3. Interface morphology of columnar crystals affected by different sizes of straining flow, where $t = 256, $ $\varGamma = 0.25, $ $M_{\rm{C}}^1=0.01, $ $M_{\rm{C}}^2=0.02, $ $C_{{\rm{L}}, \infty }^1 = 1.0, $ $C_{{\rm{L}}, \infty }^2 = 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ ${m_2} = - 2.33, $ $E = 0.3, $ ${M_{\rm{k}}} = 0.01, $ ${k_{\rm{T}}} = 1.23, $ ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^1 = 0.01, $ $\lambda _{\rm{D}}^{\rm{2}} = 0.02, $ $\varepsilon = 0.05.$ A is 0.9, 0.6, 0.3, 0 from left to right, respectively.

图 4 在$Oxy$平面上柱状晶体界面随时间的演化, 其中$\varGamma \!=\! 0.25, $ $M_{\rm{C}}^1 \!=\! 0.01, $ $M_{\rm{C}}^{\rm{2}} \!=\! 0.02, $ $C_{{\rm{L}}, \infty }^1 \!=\! 1.0, $ $C_{{\rm{L}}, \infty }^2 \!=\! 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ ${m_2} = - 2.33, $ E = 0.3, M_k = 0.01, k_T = 1.23, ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{1}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{2}} = 0.02, $ $\varepsilon = 0.05, $ $A = 0.9$
Figure4. Evolution of columnar crystal interface with time in the $Oxy$ plane, where $\varGamma = 0.25, $ $M_{\rm{C}}^1=0.01, $ $M_{\rm{C}}^{\rm{2}}=0.02, $ $C_{{\rm{L}}, \infty }^1 = 1.0, $ $C_{{\rm{L}}, \infty }^2 = 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ m₂ = –2.33, E = 0.3, M_k = 0.01, k_T = 1.23, ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{1}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{2}} = 0.02, $ $\varepsilon = 0.05, $ $A = 0.9$

图 5 应变流对不同杂质含量柱状晶体界面形态的影响, 其中$t = 256, $ $\varGamma = 0.25, $ $C_{{\rm{L}}, \infty }^1 = 1.0, $ $C_{{\rm{L}}, \infty }^2 = 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ ${m_2} = - 2.33, $ $E = 0.3, $ ${M_{\rm{k}}} = 0.01, $ ${k_{\rm{T}}} = 1.23, $ ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{1}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{2}} = 0.02, $ $\varepsilon = 0.05$
Figure5. Effect of straining flow on the interface morphology of columnar crystals in different impurity content, where $t = 256, $ $\varGamma = 0.25, $ $C_{{\rm{L}}, \infty }^1 = 1.0, $ $C_{{\rm{L}}, \infty }^2 = 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ ${m_2} = - 2.33, $ $E = 0.3, $ ${M_{\rm{k}}} = 0.01, $ ${k_{\rm{T}}} = 1.23, $ ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{1}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{2}} = 0.02, $ $\varepsilon = 0.05$

图 6 应变流对不同杂质含量柱状晶体界面形态的影响, 其中$t = 256, $ $\varGamma = 0.25, $ $C_{{\rm{L}}, \infty }^1 = 1.0, $ $C_{{\rm{L}}, \infty }^2 = 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ ${m_2} = - 2.33, $ $E = 0.3, $ ${M_{\rm{k}}} = 0.01, $ ${k_{\rm{T}}} = 1.23, $ ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{1}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{2}} = 0.02, $ $\varepsilon = 0.05$
Figure6. Effect of straining flow on the interface morphology of columnar crystals in different impurity content, where $t = 256, $ $\varGamma = 0.25, $ $C_{{\rm{L}}, \infty }^1 = 1.0, $ $C_{{\rm{L}}, \infty }^2 = 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ ${m_2} = - 2.33, $ $E = 0.3, $ ${M_{\rm{k}}} = 0.01, $ ${k_{\rm{T}}} = 1.23, $ ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{1}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{2}} = 0.02, $ $\varepsilon = 0.05$

图 7 在$R, \theta $平面上, 柱状晶体界面杂质浓度$C_{\rm{L}}^1$随$\theta $的变化情况, 其中$t = 256, $ $\varGamma = 0.25, $ $M_{\rm{C}}^{\rm{1}}=0.01, $ $M_{\rm{C}}^{\rm{2}}=0.02, $ $C_{{\rm{L}}, \infty }^1 = 1.0, $ $C_{{\rm{L}}, \infty }^2 = 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ ${m_2} = - 2.33, $ $E = 0.3, $ ${M_{\rm{k}}} = 0.01, $ ${k_{\rm{T}}} = 1.23, $ ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{1}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{2}}=0.02, $ $\varepsilon = 0.05, $ $A = 0.9$
Figure7. The change of impurity concentration at the interface of columnar crystal in the $R, \theta $ plane, where $t = 256, $ $\varGamma = 0.25, $ $M_{\rm{C}}^{\rm{1}}=0.01, $ $M_{\rm{C}}^{\rm{2}}=0.02, $ $C_{{\rm{L}}, \infty }^2 = 3.0, $ ${A_\lambda } = 3.3, $ ${m_1} = - 1.6, $ ${m_2} = $ –2.33 $E = 0.3, $ ${M_{\rm{k}}} = 0.01, $ ${k_{\rm{T}}} = 1.23, $ ${\lambda _{\rm{S}}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{1}} = 0.01, $ $\lambda _{\rm{D}}^{\rm{2}}=0.02, $ $\varepsilon = 0.05, $ $A = 0.9$

本文从理论上研究了straining流对三元过冷熔体中柱状晶体生长的影响. 通过计算给出了晶体生长界面的解析表达式, 由表达式可以看出straining流是产生不规则柱状晶体的重要原因. 在分析straining流对三元熔体中柱状晶体生长的影响时, 发现流入的straining流加快了界面的生长速度, 而流出的straining流减缓了界面的生长速度. 同时发现随着流动速度的增大, 界面变形也更为显著. 我们比较了纯熔体、二元熔体和三元熔体中柱状晶体的生长, 发现合金熔体中柱状晶体的生长受straining流的影响更为显著.

$\small\begin{split}{C_{{\rm{L}}0}} =\;& \frac{{2\varGamma }}{{{R_0}}}{g_0} - {E^{ - 1}}{M_{\rm{K}}}\frac{{{\rm{d}}{g_0}}}{{{\rm{d}}t}} + (M_{\rm{C}}^{\rm{1}}S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2)({g_0} + {R_0}{\bar g_0} - {A_\lambda }R_0^2 + {g_0})\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}\\ & {{ + (}}M_{\rm{C}}^1S_{\rm{p}}^1+M_{\rm{C}}^2S_{\rm{p}}^2)({R_0} + {R_0}{\bar g_0})\frac{{{\rm{d}}{g_0}}}{{{\rm{d}}t}} - \frac{1}{4}R_0^2{\lambda _{\rm{S}}}{A_\lambda }\frac{{{\rm{d}}{b_0}}}{{{\rm{d}}t}}+\frac{1}{4}({\lambda _{\rm{S}}}{A_\lambda } + {A_\lambda } + 1)R_0^2\frac{{{\rm{d}}{b_0}(t)}}{{{\rm{d}}t}} \\ & + \frac{1}{4}({A_\lambda } + 1)R_0^2{b_0}(t)\frac{{{\rm{d}}{{\bar R}_0}}}{{{\rm{d}}t}} + {A_\lambda }{R_0}{b_0}(t) +\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}({A_\lambda }{R_0}{\bar g_0}(t) - \frac{1}{4}{b_0}(t){R_0} + {g_0}(t)),\end{split}\tag{A1}$

$\small\begin{split}{C_{{\rm{L2}}}} =\;& - 6\varGamma {g_2} - {E^{ - 1}}{M_{\rm{K}}}R_0^2\frac{{{\rm{d}}{g_2}}}{{{\rm{d}}t}} + \left(\frac{1}{{\lambda _{\rm{D}}^1}}M_{\rm{C}}^1S_{\rm{p}}^1 + \frac{1}{{\lambda _{\rm{D}}^2}}M_{\rm{C}}^2S_{\rm{p}}^2\right)\frac{1}{4}{A_0}R_0^5\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}(1 + 2\ln {R_0})\\ &+ (M_{\rm{C}}^1S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2)\left(\frac{1}{2}R_0^2\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{g_2} + \frac{1}{2}R_0^3\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{\bar g_2} - \frac{1}{2}R_0^3\frac{{{\rm{d}}{g_2}}}{{{\rm{d}}t}}\right)\\ &+ R_0^2\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{g_2} + {A_\lambda }R_0^3\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{\bar g_2} + \frac{1}{2}{A_0}R_0^5\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}\ln {R_0},\end{split}\tag{A2}$

$\small\begin{split}T_{{\rm{L}}10}^* =\; &\frac{1}{4}\frac{{{b_0}(t)}}{{{R_0}}}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{r^2}{{\rm{e}}^{{{\bar R}_0} - {{\bar r}_0}}} + \frac{1}{4}{r^2}({A_\lambda } + \ln {R_0} + 1) \cdot \left(\frac{{{\rm{d}}{b_0}(t)}}{{{\rm{d}}t}} + {b_0}(t)\frac{{{\rm{d}}{{\bar R}_0}}}{{{\rm{d}}t}}\right){{\rm{e}}^{{{\bar R}_0} - {{\bar r}_0}}} - r\ln r{b_0}(t){{\rm{e}}^{{{\bar R}_0} - {{\bar r}_0}}}\\ & + ({A_\lambda } + \ln {R_0})r{b_0}(t){{\rm{e}}^{{{\bar R}_0} - {{\bar r}_0}}} + \frac{1}{4}{r^2}\ln r\left(\frac{{{\rm{d}}{b_0}(t)}}{{{\rm{d}}t}} + {b_0}(t)\frac{{{\rm{d}}{{\bar R}_0}}}{{{\rm{d}}t}}\right){{\rm{e}}^{{{\bar R}_0} - {{\bar r}_0}}},\end{split}\tag{A3}$

$\small T_{{\rm{L}}12}^* = - \frac{1}{4}{A_0}R_0^5\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{{\rm{e}}^{{{\bar R}_0} - {{\bar r}_0}}}{r^{ - 2}}\ln r - \frac{1}{4}{A_0}{R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{{\rm{e}}^{{{\bar R}_0} - {{\bar r}_0}}}{r^2}\ln r,\tag{A4}$

$\small\begin{split}{C_{{\rm{S0}}}}=\; &\frac{{2\varGamma }}{{R_0^2}}{g_0} - {E^{ - 1}}{M_{\rm{k}} }\frac{{{\rm{d}}{g_0}}}{{{\rm{d}}t}} + (M_{\rm{C}}^1S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2)({g_0} + {R_0}{\bar g_0} - {A_\lambda }R_0^2 + {g_0})\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}\\& + (M_{\rm{C}}^1S_{\rm{p}}^1+M_{\rm{C}}^2S_{\rm{p}}^2)({R_0} + {R_0}{\bar g_0})\frac{{{\rm{d}}{g_0}}}{{{\rm{d}}t}} - \frac{1}{4}R_0^2{\lambda _{\rm{S}}}{A_\lambda }\frac{{{\rm{d}}{b_0}}}{{{\rm{d}}t}},\end{split}\tag{A5}$

$\small \begin{split}{C_{{\rm{S2}}}} =\; & - \frac{{6\varGamma }}{{R_0^4}}{g_2} - {E^{ - 1}}{M_{\rm{K}}}\frac{1}{{R_0^2}}\frac{{{\rm{d}}{g_2}}}{{{\rm{d}}t}} + \left(\frac{1}{{\lambda _{\rm{D}}^1}}M_{\rm{C}}^1S_{\rm{p}}^1 + \frac{1}{{\lambda _{\rm{D}}^2}}M_{\rm{C}}^2S_{\rm{p}}^2\right)\frac{1}{4}{A_0}{R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}(1 + 2\ln {R_0})\\&+ (M_{\rm{C}}^1S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2)\left(\frac{1}{2}\frac{1}{{R_0^2}}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{g_2} + \frac{1}{2}\frac{1}{{{R_0}}}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{\bar g_2} - \frac{1}{2}\frac{1}{{{R_0}}}\frac{{{\rm{d}}{g_2}}}{{{\rm{d}}t}}\right),\end{split}\tag{A6}$

$\smallT_{{\rm{S10}}}^{\rm{*}} = \frac{1}{4}{r^2}{\lambda _{\rm{S}}}{A_\lambda }\frac{{{\rm{d}}{b_0}(t)}}{{{\rm{d}}t}},\;\;{\text{其中}}~{b_0}(t) = {R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}},\tag{A7}$

$\small\begin{split}C_{{\rm{C0}}}^1=\;& - S_{\rm{p}}^1({R_0} + {R_0}{\bar g_0})\frac{{{\rm{d}}{g_0}}}{{{\rm{d}}t}} + S_{\rm{p}}^1\left(2{A_\lambda }R_0^2 - {g_0} - \frac{{{R_0}}}{{4\lambda _{\rm{D}}^{\rm{1}}}}{b_0}\right)\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}\\ & - \frac{{R_0^2}}{4}\frac{1}{{\lambda _{\rm{D}}^1}}S_{\rm{p}}^1\frac{{{\rm{d}}{b_0}(t)}}{{{\rm{d}}t}}({A_\lambda } + 1) - \frac{{R_0^2}}{{4\lambda _{\rm{D}}^1}}S_{\rm{p}}^1{b_0}(t)\frac{{{\rm{d}}{{\bar g}_0}}}{{{\rm{d}}t}}({A_\lambda } + 1),\end{split}\tag{A8}$

$\small C_{{\rm{C2}}}^1 = - \frac{1}{4}\frac{1}{{\lambda _{\rm{D}}^1}}{A_0}S_{\rm{p}}^1R_0^5\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}} + \frac{1}{2}S_{\rm{p}}^1R_0^3\frac{{{\rm{d}}{g_2}}}{{{\rm{d}}t}} + \frac{1}{2}S_{\rm{p}}^1R_0^2\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{g_2} + \frac{1}{2}S_{\rm{p}}^1R_0^3\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{\bar g_2},\tag{A9}$

$\small\begin{split}C_{{\rm{C10}}}^{1*} =\; & \frac{1}{4}{r^2}\frac{1}{{\lambda _{\rm{D}}^1}}S_{\rm{p}}^1\left(\frac{{{\rm{d}}{b_0}(t)}}{{{\rm{d}}t}} + {b_0}(t)\frac{{{\rm{d}}{{\bar R}_0}}}{{{\rm{d}}t}}\right)({A_\lambda } + \ln {R_0} + 1 - \ln r){{\rm{e}}^{{{\bar R}_0} - {{\bar r}_0}}}\\ &+ rS_{\rm{p}}^1{R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}(\ln r - {A_\lambda } - \ln {R_0}){{\rm{e}}^{{{\bar R}_0} - {{\bar r}_0}}} + \frac{{{r^2}}}{4}\frac{1}{{\lambda _{\rm{D}}^1}}S_{\rm{p}}^1\frac{{{b_0}}}{{{R_0}}}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{{\rm{e}}^{{{\bar R}_0} - {{\bar r}_0}}},\end{split}\tag{A10}$

$\small C_{{\rm{C12}}}^{1*} = - \frac{1}{4}\frac{1}{{\lambda _{\rm{D}}^1}}{A_0}R_0^5S_{\rm{p}}^1\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{{\rm{e}}^{{{\bar R}_0} - \bar r}}{r^{ - 2}}\ln r - \frac{1}{4}\frac{1}{{\lambda _{\rm{D}}^1}}{A_0}{R_0}S_{\rm{p}}^1\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{{\rm{e}}^{{{\bar R}_0} - \bar r}}{r^2}\ln r,\tag{A11}$

$\small\begin{split}C_{{\rm{C0}}}^2=\; & - S_{\rm{p}}^2({R_0} + {R_0}{\bar g_0})\frac{{{\rm{d}}{g_0}}}{{{\rm{d}}t}} + S_{\rm{p}}^2\left(2{A_\lambda }R_0^2 - {g_0} - \frac{{{R_0}}}{{4\lambda _{\rm{D}}^2}}{b_0}\right)\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}\\ &- \frac{{R_0^2}}{4}\frac{1}{{\lambda _{\rm{D}}^2}}S_{\rm{p}}^2\frac{{{\rm{d}}{b_0}(t)}}{{{\rm{d}}t}}({A_\lambda } + 1) - \frac{{R_0^2}}{{4\lambda _{\rm{D}}^2}}S_{\rm{p}}^2{b_0}(t)\frac{{{\rm{d}}{{\bar g}_0}}}{{{\rm{d}}t}}({A_\lambda } + 1),\end{split}\tag{A12}$

$\small C_{{\rm{C2}}}^2 = - \frac{1}{4}\frac{1}{{\lambda _{\rm{D}}^2}}{A_0}S_{\rm{p}}^2R_0^5\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}} + \frac{1}{2}S_{\rm{p}}^2R_0^3\frac{{{\rm{d}}{g_2}}}{{{\rm{d}}t}} + \frac{1}{2}S_{\rm{p}}^2R_0^2\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{g_2} + \frac{1}{2}S_{\rm{p}}^2R_0^3\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{\bar g_2},\tag{A13}$

$\small\begin{split}C_{{\rm{C1}}0}^{2*} =\;& \frac{1}{4}{r^2}\frac{1}{{\lambda _{\rm{D}}^2}}S_{\rm{p}}^{\rm{2}}\left(\frac{{{\rm{d}}{b_0}(t)}}{{{\rm{d}}t}} + {b_0}(t)\frac{{{\rm{d}}{{\bar R}_0}}}{{{\rm{d}}t}}\right)({A_\lambda } + \ln {R_0} + 1 - \ln r){{\rm{e}}^{{{\bar R}_0} - {{\bar r}_0}}}\\ &+ rS_{\rm{p}}^2{R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}(\ln r - {A_\lambda } - \ln {R_0}){{\rm{e}}^{{{\bar R}_0} - {{\bar r}_0}}} + \frac{{{r^2}}}{4}\frac{1}{{\lambda _{\rm{D}}^2}}S_{\rm{p}}^2\frac{{{b_0}}}{{{R_0}}}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{{\rm{e}}^{{{\bar R}_0} - {{\bar r}_0}}},\end{split}\tag{A14}$

$\small C_{{\rm{C}}12}^{2*} = - \frac{1}{4}\frac{1}{{\lambda _{\rm{D}}^2}}{A_0}R_0^5S_{\rm{p}}^2\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{{\rm{e}}^{{{\bar R}_0} - \bar r}}{r^{ - 2}}\ln r - \frac{1}{4}\frac{1}{{\lambda _{\rm{D}}^2}}{A_0}{R_0}S_{\rm{p}}^2\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{{\rm{e}}^{{{\bar R}_0} - \bar r}}{r^2}\ln r,\tag{A15}$

$\small\begin{split}\frac{{{\rm{d}}{g_0}}}{{{\rm{d}}t}} =\; & - \frac{1}{{{R_0}}}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{g_0} + \frac{1}{2}{b_0}(t)\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}} + \frac{1}{2}\left({k_{\rm{T}}}{\lambda _{\rm{S}}}{A_\lambda } + {A_\lambda } + \frac{1}{2}\right){R_0}(t)\frac{{{\rm{d}}{b_0}}}{{{\rm{d}}t}}{{ + (2}}{A_\lambda } - 1){R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}\\ &- {A_\lambda }{R_0}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{\bar g_0} + \frac{1}{2}({A_\lambda } + \frac{1}{2})R_0^2\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}\frac{{{\rm{d}}{{\bar R}_0}}}{{{\rm{d}}t}},\end{split}\tag{A16}$

$\small\begin{split}\frac{{{\rm{d}}{g_2}}}{{{\rm{d}}t}} =\; & \frac{{ - 12\varGamma R_0^{ - 2}(1 + {K_{\rm{T}}})}}{{{R_0} + 2{E^{ - 1}}{M_{\rm{K}}}(1 + {K_{\rm{T}}}) + (M_{\rm{C}}^1S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2){R_0}(1 + {K_{\rm{T}}})}}{g_2}\\ &+ \frac{{(M_{\rm{C}}^1S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2)(1 + {K_{\rm{T}}}) + 1}}{{{R_0} + 2{E^{ - 1}}{M_{\rm{K}} }(1 + {K_{\rm{T}}}) + (M_{\rm{C}}^1S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2){R_0}(1 + {K_{\rm{T}}})}}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{g_2}\\ &+ \frac{{(M_{\rm{C}}^1S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2)(1 + {K_{\rm{T}}}){R_0} + {A_\lambda }{R_0}(2 - {R_0})}}{{{R_0} + 2{E^{ - 1}}{M_{\rm{K}}}(1 + {K_{\rm{T}}}) + (M_{\rm{C}}^1S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2){R_0}(1 + {K_{\rm{T}}})}}\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}{\bar g_2}\\ &+ \dfrac{{\Big(\dfrac{1}{{\lambda _{\rm{D}}^1}}M_{\rm{C}}^1S_{\rm{p}}^1 + \dfrac{1}{{\lambda _{\rm{D}}^2}}M_{\rm{C}}^2S_{\rm{p}}^2\Big)R_0^3(1 + 2\ln {R_0})(1 + {K_{\rm{T}}}) + R_0^3(1 + 2\ln {R_0})}}{{{R_0} + 2{E^{ - 1}}{M_{\rm{K}}}(1 + {K_{\rm{T}}}) + (M_{\rm{C}}^1S_{\rm{p}}^1 + M_{\rm{C}}^2S_{\rm{p}}^2){R_0}(1 + {K_{\rm{T}}})}}\frac{1}{2}A\frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}},\end{split}\tag{A17}$

$\small\frac{{{\rm{d}}{b_0}}}{{{\rm{d}}t}} = \frac{{{E^{ - 1}}{M_{\rm{K}}} + 2\varGamma {A_\lambda }(1 + M_{\rm{C}}^1S_{\rm{p}} ^1 + M_{\rm{C}}^2S_{\rm{p}} ^2)}}{{{{[{A_\lambda }{R_0}(1 + M_{\rm{C}}^1S_{\rm{p}} ^1 + M_{\rm{C}}^2S_{\rm{p}} ^2) + {E^{ - 1}}{M_{\rm{K}}}]}^2}}} \cdot \frac{{{\rm{d}}{R_0}}}{{{\rm{d}}t}}.\tag{A18}$