【计算题】推导成分过冷的判别式并求出成分过冷的过冷度。

【计算题】推导成分过冷的判别式并求出成分过冷的过冷度。

用户现在需要推导成分过冷的判别式并求出成分过冷的过冷度。首先，我得回忆一下成分过冷的基本概念。成分过冷是凝固过程中，由于固液界面前沿溶质再分配导致熔体中溶质浓度变化，进而引起熔点变化，使得熔体实际温度低于其熔点而产生的过冷现象，区别于纯粹的热过冷。

首先，推导判别式需要考虑几个关键因素：固液界面前沿的温度梯度（G）、晶体生长速度（R）、溶质扩散系数（D）、平衡分配系数（k）、合金原始浓度（C0）以及液相线斜率（m）。这些都是合金凝固中的基本参数，需要明确它们的定义和关系。

先从平衡分配系数k开始，k = Cs/Cl，其中Cs是固相溶质浓度，Cl是液相溶质浓度，在平衡凝固时成立。当凝固以有限速度R进行时，固液界面前沿会形成溶质边界层，浓度分布不再均匀。根据稳态溶质扩散方程，可以推导出界面前沿液相中的溶质浓度分布。

稳态下，溶质扩散方程应该是：D*(d²C/dx²) = R*(dC/dx)，这个方程的解需要边界条件。在界面处（x=0），根据质量守恒，固相排出的溶质必须通过液相扩散走，所以有：R*(Cl - Cs) = -D*(dC/dx)ₓ=0。结合k = Cs/Cl，所以Cs = kCl，代入得R(Cl - kCl) = -D(dC/dx)ₓ=0，即RCl(1 - k) = -D*(dC/dx)ₓ=0，整理得(dC/dx)ₓ=0 = -RCl(1 - k)/D。这里的Cl应该是界面处的液相浓度，记为C₁*，所以(dC/dx)ₓ=0 = -RC₁(1 - k)/D。

然后解扩散方程，D*(d²C/dx²) = R*(dC/dx)，令u = dC/dx，则方程变为D*(du/dx) = Ru，分离变量得du/u = (R/D)dx，积分得lnu = (R/D)x + C，所以u = dC/dx = Aexp(Rx/D)。应用边界条件x=0时，dC/dx = -RC₁(1 - k)/D，所以A = -RC₁(1 - k)/D，因此dC/dx = -[RC₁(1 - k)/D]exp(Rx/D)。

再积分一次求C(x)：C(x) = -[RC₁(1 - k)/D](D/R)exp(Rx/D) + B = -C₁(1 - k)exp(Rx/D) + B。当x→∞时，液相浓度应该恢复到原始浓度C0，所以C(∞) = B = C0，因此C(x) = C0 - C₁*(1 - k)exp(Rx/D)。在x=0处，C(0) = C₁，代入得C₁ = C0 - C₁*(1 - k)1，整理得C₁ = C0 / [1 - (1 - k)exp(0)]？不对，x=0时exp(0)=1，所以C₁ = C0 - C₁(1 - k)，移项C₁ + C₁*(1 - k) = C0，C₁*[1 + 1 - k] = C0？不对，应该是C₁ = C0 - C₁*(1 - k) → C₁ + C₁*(1 - k) = C0 → C₁*[1 + (1 - k)] = C0？不，等号左边是C₁ = C0 - C₁*(1 - k)，把右边的C₁*(1 - k)移到左边：C₁ + C₁*(1 - k) = C0 → C₁*[1 + (1 - k)] = C0？那就是C₁*(2 - k) = C0？这不对，我肯定哪里错了。

哦，对了，平衡分配系数k的定义是k = Cs*/Cl*，其中Cs是界面处固相浓度，Cl是界面处液相浓度。在稳态凝固时，固相中的溶质浓度Cs* = kCl，而原始合金浓度C0应该等于固相浓度乘以固相体积分数加上液相浓度乘以液相体积分数，但在稳态下，沿着生长方向，固液界面以速度R移动，所以在单位时间内，通过单位面积的溶质守恒应该是：RC0 = RCs* + D*(dCl/dx)ₓ=0。这个才是正确的溶质守恒条件。因为左边是流入的溶质（原始合金浓度C0以速度R向前移动），右边是固相带走的溶质（Cs以速度R）加上液相中扩散走的溶质（扩散通量D(dCl/dx)ₓ=0，注意符号，浓度梯度如果是负的，即dCl/dx为负，那么扩散通量是正方向，即离开界面）。

对，这个才是正确的边界条件。之前的那个可能记错了。正确的质量守恒方程应该是：在固液界面处，单位时间内从液相区进入界面的溶质量等于固相带走的溶质量加上通过液相扩散离开界面的溶质量。所以：RC0 = RCs* + D*(dCl/dx)ₓ=0。因为Cs* = kCl，代入得：RC0 = RkCl + D*(dCl/dx)ₓ=0，所以D*(dCl/dx)ₓ=0 = R*(C0 - kCl)，因此(dCl/dx)ₓ=0 = R*(C0 - kCl)/D。这才是正确的边界条件，之前的那个错了，应该用这个。

然后，液相中的溶质扩散方程仍然是稳态的，即∂Cl/∂t = 0（稳态），所以对流扩散方程中，只有扩散项和迁移项：D*(d²Cl/dx²) = R*(dCl/dx)，这个是对的，因为在随界面移动的坐标系中，对流项是-R*(dCl/dx)，所以扩散方程∂Cl/∂t + v·∇Cl = D∇²Cl，稳态下∂Cl/∂t=0，v=R（沿x方向），所以R*(dCl/dx) = D*(d²Cl/dx²)，即D*(d²Cl/dx²) - R*(dCl/dx) = 0，这是二阶线性常微分方程，通解为Cl(x) = A + B*exp(Rx/D)，其中x是距离界面的距离，x=0为界面位置，x>0为液相区。

现在应用边界条件：1. 当x→∞时，Cl(x) = C0（远离界面处液相浓度恢复到原始浓度）；2. 在x=0处，由质量守恒得到的(dCl/dx)ₓ=0 = R*(C0 - kCl)/D，其中Cl*是界面处液相浓度（x=0处的Cl）。

由第一个边界条件，当x→∞时，exp(Rx/D)→∞，如果B是正的，Cl(x)会趋于无穷，这不可能，所以B必须为0？不对，x是从界面向外的距离，应该是x=0在界面，x>0指向液相内部，当x→∞时，Cl(x)=C0，所以Cl(x) = A + Bexp(Rx/D)，当x→∞，exp(Rx/D)→∞，除非B=0，否则Cl(x)→∞，所以B=0的话Cl(x)=A=C0，这不对，那说明我的坐标系可能反了，应该取x=0为界面，x<0指向液相内部，这样当x→-∞时，Cl(x)=C0，此时exp(Rx/D)当x→-∞时→0，所以Cl(x)=A + Bexp(Rx/D)，当x→-∞，exp项→0，所以Cl(-∞)=A=C0，正确。然后在x=0处，Cl(0)=Cl*=A + Bexp(0)=C0 + B，所以B=Cl - C0。

然后对Cl(x)求导：dCl/dx = B*(R/D)exp(Rx/D) = (Cl - C0)(R/D)exp(Rx/D)。在x=0处，dCl/dx|ₓ=0 = (Cl - C0)(R/D)。根据之前的质量守恒方程，dCl/dx|ₓ=0 = R*(C0 - kCl)/D。所以联立：

(Cl* - C0)(R/D) = R(C0 - kCl)/D

两边消去R/D：Cl* - C0 = C0 - kCl → Cl* + kCl = 2C0？不对，Cl* - C0 = C0 - kCl → Cl* + kCl = C0 + C0 → Cl*(1 + k) = 2C0？这也不对，我肯定哪里又错了。哦，符号问题！扩散通量的方向，如果浓度梯度是正的（dCl/dx>0），则扩散方向是沿x负方向（指向界面），而之前的质量守恒方程中的D*(dCl/dx)应该是扩散通量，其方向应该是离开界面还是进入界面？

正确的溶质守恒：在界面处，固相以速度R生长，固相中的溶质浓度是Cs*，所以单位时间单位面积固相带走的溶质是RCs。液相中的溶质以浓度Cl到达界面，所以单位时间单位面积液相带来的溶质是RCl*。两者的差值就是必须通过扩散离开界面的溶质量，即R*(Cl* - Cs*) = -D*(dCl/dx)ₓ=0。这里负号是因为如果dCl/dx是负的（浓度随x增加而减小，即界面处浓度高于远处），则扩散通量是正方向（离开界面）。Cs*=kCl，所以左边=R*(Cl* - kCl)=RCl(1 - k)。右边=-D*(dCl/dx