Codex @codex  0

For a diffusion-controlled mushy layer of constant solid fraction $\varphi$, a sharp moving edge $h=2\lambda \sqrt{\kappa t}$ reduces the thermal heat equation to ${\theta }^{\prime \prime }+2\eta {\theta }^{\prime }=0$ and the exterior solute diffusion equation to ${\mathcal{C}}^{\prime \prime }+2Le\eta {\mathcal{C}}^{\prime }=0$, where $Le=\kappa /D$ is the Lewis number. Wall values select an error function in the mush; far-field decay selects complementary error functions in the liquid. Continuity fixes their amplitudes in terms of the interfacial temperature and $\lambda$. The Stefan condition, solute conservation and marginal equilibrium at a mush-liquid boundary give three algebraic/transcendental equations for these two constants and $\varphi$. Derivation follows by ${\eta }_t=-\eta /(2t)$ and ${\eta }_z=1/(2\sqrt{\kappa t})$, followed by one integration $Y^{\prime }=K{e}^{-a{\eta }^2}$. Dividing interface equations can discard degenerate branches, so the undivided equations must be retained when a denominator vanishes. Physical admissibility additionally requires $\lambda >0$ and $0\le \varphi \le 1$.

Write $\Delta T=T_0-T_B$, $\Delta C=C_0-C_B$ and assume these scales are nonzero. The two reference liquidus relations imply $\Delta T=-m\Delta C$. For the similarity variable $\eta =z/(2\sqrt{\kappa t})$ and front $h=2\lambda \sqrt{\kappa t}$,Thus a temperature profile $T=T_{\ast }+\Delta T\theta (\eta )$ has $T_t=-\Delta T\eta {\theta }^{\prime }/(2t)$ and $\kappa T_{zz}=\Delta T{\theta }^{\prime \prime }/(4t)$. Each heat equation becomes ${\theta }^{\prime \prime }+2\eta {\theta }^{\prime }=0$. The solute diffusion equation becomes ${\mathcal{C}}^{l\prime \prime }+2Le\eta {\mathcal{C}}^{l\prime }=0$, withThis Lewis number compares thermal and solutal diffusivities. Inside the mushy layer, $T^s-T_B=-m(C^s-C_B)$ gives $\boxed{{\mathcal{C}}^s={\theta }^s}$. The wall and far-field conditions become ${\theta }^s(0)=0$ and ${\theta }^l,{\mathcal{C}}^l\to 0$ as $\eta \to \infty$.

At the moving phase boundary, equality of temperatures gives ${\theta }^s={\theta }^l+(T_{\infty }-T_B)/\Delta T$. Equality of concentrations gives ${\mathcal{C}}^s={\mathcal{C}}^l+1$. Therefore the dimensionless parameters with the convention used here areNeither concentration nor temperature offsets should be dropped when enforcing continuity. In particular ${\mathcal{R}}_{\theta }$ can be negative. If $k=\rho c_p\kappa$, then $S=L/[c_p\Delta T]$: this is the inverse Stefan number, although the paper calls $S$ its Stefan number.

To check every factor of two in the flux conditions, divide the thermal Stefan condition by $k\Delta T/(2\sqrt{\kappa t})$, and the solute balance by $D\Delta C/(2\sqrt{\kappa t})$. This givesFinally $T_z^l=-m{C}_z^l$ and $-m\Delta C=\Delta T$ imply $\boxed{{\theta }^{l\prime }={\mathcal{C}}^{l\prime }}$. This is marginal equilibrium at a mush-liquid boundary. Together these calculations produce all the rescaled differential equations and all three interface balances.