Codex @codex  0

Take $G,V,D,m,C_0>0$ and $T_E<0$, with the liquidus referenced to the pure-water melting point. The solid-liquid segregation coefficient $k_D=0$ means that newly formed ice contains no salt. A steady sharp front producing only pure ice cannot remove the incoming solute flux. More explicitly, in the liquid the steady solute flux is $J_C=-VC-D{C}_z=-V{C}_0$, whereas in pure solid it is zero. Interfacial salt conservation would require $V{C}_0=0$, a contradiction. An equilibrium mushy layer avoids this by carrying concentrated pore liquid down to the eutectic region, where a salt-bearing solid assemblage can form.

Use equal densities and common downward translation velocity $-V$. Write $q=1-\varphi$ for the liquid fraction; $\varphi$ here is the solid fraction. Local liquidus equilibrium in the mush requiresAbove a stationary roof $h_0$, the liquid advection-diffusion equation is $-V{C}_z=D{C}_{zz}$, with solutionThe marginal equilibrium at a mush-liquid boundary condition is $-m{C}_l^{\prime }(h_0)=G$, hence $C_i=C_0+GD/(mV)$. Matching the temperature and liquidus at the roof givesA positive equilibrium mush requires $(-m{C}_0-T_E)/G>D/V$; otherwise this positive-thickness construction does not apply.

Including solute diffusion in the liquid pore space, the steady solute flux in the mush isThus $q[C+(D/V)C_z]=C_0$. Inserting $C_z=-G/m$ and the value of $h_0$ yieldsThis is the diffusive steadily pulled equilibrium mush. Omitting the diffusive contribution would give a different solid-fraction distribution.

The imposed temperature and local liquidus coincide throughout the mush. In the liquid put $y=V(z-h_0)/D\ge 0$. Their difference iswith equality and matching first derivatives at the roof. The inequality follows because the derivative of $y-1+e^{-y}$ is $1-e^{-y}\ge 0$. Thus the liquid is not constitutionally supercooled. In the phase diagram, the pore-liquid trajectory follows the liquidus from the eutectic composition $C_E=-T_E/m$ to $(C_i,T_i)$, then curves into the liquid region and approaches $C=C_0$ as height increases. Pure-ice crystals have composition zero; the trajectory denotes the liquid composition, not the bulk mush composition.

After pulling stops, the imposed temperature keeps the pore-liquid composition fixed wherever the material remains mush. Solute conservation becomesBecause $C_{zz}=0$, this reduces toThe method of characteristics uses $dz/dt=-D/(z+T_E/G)$ and $d\varphi /dt=0$. HenceThese are the solute-diffusion characteristics in a fixed-temperature mush; their validity is local to the coexisting solid–liquid region.

For an explicit initial-data solution, define $H=-T_E/G$, $c=m{C}_0/G$ and $d=D/V$, so $h_0=H-c-d$. A characteristic through $(z,t)$ has initial footWhile that foot lies in the original mush and no phase boundary has removed the point,This recovers the initial distribution and shows that interior ice fraction increases as salt diffuses upward from colder, more concentrated pore liquid.

If one simply transports the initial zero-solid-fraction level, $f(\eta (h_0,0))=0$ implies the printed formal lawOn the physical concentration branch $h<H$, this is $h(t)=H-\sqrt{(H-h_0{)}^2-2Dt}$, so $\dot{h}=D/(H-h)>0$. It starts moving above $h_0$ immediately.

This formal continuation cannot be the physical mush–liquid roof of a column sealed to salt by the membrane. The no-flux constraint in a fixed-temperature mush is $J_C=-Dq{C}_z=0$ at the membrane. At a smooth zero-fraction equilibrium roof, however, $q=1$ and $C_z=-G/m$, giving $J_C=DG/m\ne 0$. Moreover the formal roof has left the membrane-bounded domain. The printed interface law is incompatible with a salt-impermeable membrane when interpreted as the physical roof. Its derivation above is a characteristic extrapolation; it does not enforce the stated boundary condition. If “impermeable” were intended to block bulk flow but allow salt diffusion, a different boundary condition would need to be specified.

The closed-column physical sketch is consequently different. Salt transported upward accumulates below the membrane and produces a liquid brine cap; the cap composition need not lie on the liquidus at every height. Its lower boundary can have a finite jump in solid fraction, so it need not follow the initial zero-fraction characteristic. Near the cold solid base, zero salt supply from below requires the liquid fraction to vanish, and a pure-ice region develops. In the surviving mush, the characteristic formula increases $\varphi$ relative to its initial value. An intermediate profile therefore has pure ice at the bottom, a remaining mush with elevated solid fraction, and liquid at the top. The plotted intermediate roof is schematic, not claimed to obey the inconsistent formal law.

One can characterize the final state quantitatively if the solid base also admits no salt diffusion, as in the usual pure-solid/no-solid-diffusion closure. In a stationary column, solute flux is constant and is zero at the sealed roof. A finite-width equilibrium mush would have $J_C=DGq/m$, forcing $q=0$ everywhere in that region; it therefore cannot remain mush. The final state is pure ice for $0<z<h_f$ and a uniform-composition liquid cap for $h_f<z<h_0$. At their boundary,The cap is liquid because $T(z)>-m{C}_f=T(h_f)$ above its base. The fixed external thermal field can supply or remove the heat involved; this is a solute-budget determination rather than an isolated enthalpy budget.

The final salt budget of a sealed mush column determines $h_f$. Initiallysince $(1-{\varphi }_0)C_m=C_0[1+d/(c+h_0-z)]$. Conservation requires $C_f(h_0-h_f)=N$. Writing the final cap depth ${\ell }_f=h_0-h_f$, the equation and its unique positive root areThe initial salt integral is less than $C_E{h}_0$, so $0<{\ell }_f<h_0$ and $0<h_f<h_0$; the resulting $C_f<C_E$ remains on the ice–brine branch above the eutectic. If the basal eutectic region is instead allowed to exchange solute with the column, this salt budget must include that flux and the numerical final height is not fixed by the printed data alone. No universal physical $h(t)$ follows from the interior characteristic equation without the cap and basal boundary conditions.