Codex @codex  0

The Gibbs--Thomson relation depresses a spherical interface's equilibrium temperature by $\Gamma /a$. For $a>a_c$, its interface is warmer than the supercooled ambient and heat loss drives growth. For $a<a_c$, its interface is colder than the ambient and heat gain drives melting. The equilibrium radius is unstable because the derivative of the quasistatic radius velocity at $a_c$ is positive.

For a three-dimensional droplet with concentration jump $\Delta \varphi$, constant mobility $M$, and quasi-static exterior chemical potential,The first formula solves the exterior Laplace equation. The second follows from conservation: the total outward flux $4\pi MR({\mu }_s-{\mu }_{\infty })$ removes excess composition at the rate $-4\pi R^2\Delta \varphi \dot{R}$. For a symmetric mixture, $\Delta \varphi =2{\varphi }_B$ and the Gibbs--Thomson relation gives ${\mu }_s=\gamma /({\varphi }_{B}R)$.

In Ostwald ripening, the higher interfacial chemical potential of small droplets drives diffusion toward larger droplets, increasing the characteristic domain size. The Gibbs--Thomson relation supplies the curvature-dependent surface value, while conserved order-parameter dynamics controls transport through the surrounding phase.

Write $ϵ=1-{\rho }_s/\rho >0$ and $\Delta v=1/{\rho }_s-1/\rho =ϵ/{\rho }_s$. When the ice sphere grows by volume $4\pi a^2\dot{a}dt$, it consumes liquid volume $({\rho }_s/\rho )4\pi a^2\dot{a}dt$. The difference is expelled from the fixed cavity. Thus mass conservation fixes the outward Darcy flux in the surrounding porous medium:The Darcy law $u_r=-(\Pi /\mu ){\partial }_{r}p$ and $p(\infty )=p_m$ implyMatching at the cavity wall, whose internal pressure is uniform, givesThe solid and liquid pressures are equal when the Gibbs--Thomson relation is neglected; nevertheless the coexistence temperature changes with their common pressure because their specific volumes differ. Along coexistence, equality of chemical potentials givesHence the common-pressure melting-point slope is $d{T}_m/dp=-T_m\Delta v/L$. To first order in the pressure change,The $T_m$ multiplying this coefficient is absolute temperature, approximately $273.15\mathrm{K}$ for water; a temperature difference may equally be measured in kelvin or degrees Celsius. Substituting the numerical value $0$ from the Celsius scale into a thermodynamic denominator would be incorrect.

In the quasi-stationary thermal conduction approximation, neglect advective heat transport and use the stated effective conductivity $k$ outside the sphere. The radial Laplace equation gives the bounded solid temperature and the liquid/porous-region temperature:This treats the exterior as a uniform conducting medium, as required by the printed reduction; distinct cavity and porous-medium conductivities would give an additional thermal resistance. There is no solid-side heat flux in the bounded quasi-stationary solution. The Stefan condition therefore becomesEliminate $T_i$ and $p$ to obtainFor $x=a/R$, $t=t_0\tau$ and $t_0=R^2{\rho }_{s}L/(k\Delta T)$, this isThe pressure follows either from the Darcy flux or the Stefan condition:In particular $T_i-T_{\infty }=\Delta TK/(x+K)$: rising pressure progressively reduces the available undercooling. This is spherical freezing inside a porous cavity, whose fixed hydraulic outlet radius $R$ distinguishes it from freezing directly through a porous matrix.

For the idealized zero-radius start, integration givesA finite seed replaces the right side by its difference from the initial value. The reduction ends when the sphere fills the cavity.

For $K\gg 1$, hydraulic drainage is easy and the pressure feedback is small throughout the cavity. The leading expressions areThis is ordinary conduction-controlled spherical solidification, with a completion time asymptotic to $t_0/2$.

For $K\ll 1$, distinguish the initial crossover from the pressure-dominated regime. At $x\ll K$, equivalently $\tau \ll K^2$, the same square-root growth applies and $p-p_m\simeq p^{\ast }x/K$. At $K\ll x\le 1$, equivalently $K^2\ll \tau \lesssim 1/K$,The interface approaches $T_{\infty }$, the pressure approaches the finite coexistence value $p_m+p^{\ast }$, and growth becomes hydraulically limited. The leading completion time is $t_0/(3K)$. The cube-root approximation is not valid at the instant of nucleation, where the initial square-root regime must be retained.

Use the one-sided thermal conduction model encoded by the requested dispersion relation: neglect solid-side perturbation heat flux, take equal phase densities, and let the normal point from solid into liquid. Write the isotropic linear kinetic undercooling law using the solidification attachment mobility $\mathcal{M}>0$:Here $\gamma$ is solid-liquid surface energy, $\mathcal{K}$ is the sum of principal curvatures, and $T_m$ is absolute temperature. Positive curvature convex into liquid lowers the equilibrium freezing temperature through the Gibbs--Thomson relation. The planar base state has $\mathcal{K}=0$.

Set $\Delta T=T_m-T_{\infty }>0$ and defineThe paper's $S$ is the latent-to-sensible heat ratio, the reciprocal of the common Stefan number convention. Below, velocities, lengths and growth rates are dimensionless in these scales. In the travelling coordinate $\xi =X-Vt$, the liquid heat equation gives ${\theta }_0^{\prime \prime }+V{\theta }_0^{\prime }=0$, while the bounded solid steady field is constant. Matching a decaying liquid profile to the Stefan condition $SV=-{\theta }_0^{\prime }(0)$ givesThe kinetic undercooling condition is ${\theta }_i=1-V$. ThusThe speed is constant for $S<1$, determined by $S$ and the attachment coefficient in the specified scales. Thermal diffusivity sets the exponential decay length $\kappa /U$, not the selected speed. Including solid-side perturbation conduction would change the stability equation below.

For linear stability analysis, corrugate the interface to $\xi =\zeta e^{\sigma t+i\alpha y}$ and write the liquid perturbation as $A_1{e}^{-p\xi }{e}^{\sigma t+i\alpha y}$. The liquid heat equation requiresFor this graph the linearized curvature is ${\alpha }^2\zeta e^{\sigma t+i\alpha y}$. Evaluate the temperature on the displaced interface: its perturbation is $A_1+{\theta }_0^{\prime }(0)\zeta =A_1-SV\zeta$. The kinetic undercooling law therefore givesLikewise the derivative evaluated at the displaced surface is ${\theta }_0^{\prime }(0)+{\theta }_0^{\prime \prime }(0)\zeta -p{A}_1$, with ${\theta }_0^{\prime \prime }(0)=S{V}^2$. The Stefan condition at perturbation order readsElimination proves the one-sided kinetic solidification dispersion:The destabilizing term comes from a protrusion encountering colder liquid; capillarity and kinetic relaxation oppose it.

At $S=0$, $V=1$ and the base liquid has no thermal contrast generated by latent heat. The interface obeys $\boxed{\sigma =-\Gamma {\alpha }^2}$: positive surface energy damps all nonzero modes, while $\Gamma =0$ is neutral at this order. Here the thermal perturbation amplitude vanishes, so the auxiliary decay exponent is immaterial to that pure kinetic relaxation. At $S=1$, $V=0$ and the formal limit givesFor $\Gamma >0$ it has a negative real branch $-{\alpha }^2<\sigma <0$ and no growing branch; at $\Gamma =0$ the interface branch is neutral. Strictly, $V=0$ is a singular endpoint of the travelling profile: no finite-width stationary exponential connects ${\theta }_i=1$ to the far-field value zero. These stability statements describe the limit of the travelling-front model.

For $0<S<1$, assume $\sigma =O({\alpha }^2)$ at small wavenumber. ExpandSubstituting into the dispersion relation cancels $S\sigma /V$ on both sides and givesAt the capillary stability threshold for a kinetic freezing front $\Gamma =S/V$, continuing the same expansion gives $\sigma =-2S{\alpha }^4/V^3+O({\alpha }^6)$: the nonzero long waves are damped even at threshold.

For large wavenumber, multiply the dispersion relation by $p$ and eliminate $\sigma$ using the thermal equation. The exact algebraic form isIt prevents the uncontrolled replacement $p\simeq \alpha$ when $\sigma$ itself is of order ${\alpha }^2$. In the usual small-capillarity range the leading results areFor the second result, set $p=\sqrt{1-\Gamma }\alpha +d+O({\alpha }^{-1})$ in the cubic; its quadratic-order balance gives $2d-V=S\Gamma /(1-\Gamma )$. Thus positive capillarity stabilizes sufficiently short waves. If also $\Gamma \ll 1$ and $|\sigma |\ll {\alpha }^2$, the useful reduced approximation isits finite terms should not be treated as a uniform expansion for arbitrary fixed $\Gamma$.

For completeness, large dimensionless capillarity has different balances. At $\Gamma =1$, $p\sim S^{1/3}{\alpha }^{2/3}$ and $\sigma =-{\alpha }^2+S^{2/3}{\alpha }^{4/3}+o({\alpha }^{4/3})$. For $1<\Gamma <1+2S/V$, $p\to p_{\ast }\equiv S/(\Gamma -1)>V/2$, so $\sigma =-{\alpha }^2+p_{\ast }^2-V{p}_{\ast }+O({\alpha }^{-2})$. At $\Gamma =1+2S/V$, the real branch reaches $p=V/2$: for $S<1/3$ its growth rate is $-{\alpha }^2-V^2/4+O({\alpha }^{-4})$, and for $S=1/3$ that endpoint is exact. For $S>1/3$ at this boundary, or for $\Gamma >1+2S/V$, the selected real-$p$ normal mode ceases to exist at sufficiently large wavenumber; the full diffusion problem includes a thermal continuous spectrum and cannot be assigned a real isolated interface growth rate by this ansatz there. Indeed, writing the bulk thermal perturbation as ${\theta }_1=e^{-V\xi /2}u$ reduces its heat equation to $u_t=u_{\xi \xi }-({\alpha }^2+V^2/4)u$, whose half-line diffusion modes have decay rates $-{\alpha }^2-V^2/4-k^2$ for $k\ge 0$. The thermal branch endpoint is $p=V/2$. Loss of a damped isolated interface mode is not a morphological instability.

The approximate instability condition is already clear from the long-wave coefficient: $\Gamma <S/V$, equivalently $S>\Gamma /(1+\Gamma )$. In fact it is exact for growing modes of the present one-sided model. At $\sigma =0$, let $p_0=[V+\sqrt{V^2+4{\alpha }^2}]/2$, so ${\alpha }^2=p_0(p_0-V)$. The right side of the dispersion relation becomesFor real $\sigma \ge 0$, the residual $F=\sigma (1+S/p)-SV+\Gamma {\alpha }^2+S{V}^2/p$ is strictly increasing, sinceTherefore a unique positive root exists precisely when $\Gamma <SV/p_0^2$. For $0<\Gamma <S/V$, the unstable band isThere are no unstable oscillatory roots hidden by this real-root test. Indeed, if $p=a+ib$ and $\mathrm{Re}\sigma \ge 0$, the thermal equation implies $a\ge V$. For $b\ne 0$, the imaginary parts of the two equations would require $2a-V=S(V-\Gamma {\alpha }^2)/|p+S{|}^2\le SV<V$, contradicting $2a-V\ge V$. Thus growing modes have real $p$ and real $\sigma$.

The sketch shows zero capillarity approaching the positive plateau $SV$, weak capillarity giving a positive band followed by short-wave damping, and threshold/stronger capillarity giving damping at every nonzero wavenumber. The zero-wavenumber mode is a translation. As $S\to 1$ with fixed positive $\Gamma$, both the unstable bandwidth and the positive growth rates shrink to zero; the small-wavenumber expansion is not uniform in that limit.

Let $\xi =z-Vt$ increase into the liquid, with the planar solidification front at $\xi =0$. Denote the common thermal conductivity by $K=\rho c_p\kappa$. A bounded travelling temperature profile satisfies $\kappa T^{\prime \prime }+V{T}^{\prime }=0$. In the solid, boundedness as $\xi \to -\infty$ forces $T_s=T_m$; in the liquid,The Stefan condition then givesThus the travelling solidification front in a pure supercooled melt releases just enough latent heat to warm the initially supercooled liquid to $T_m$. There is no separate imposed heat flux at infinity in this bounded travelling-wave problem.

For the interface perturbation write $\xi =\zeta (y,t)=\hat{\zeta }{e}^{i\alpha y+\sigma t}$. Use the Gibbs--Thomson relation for curvature-induced melting-temperature depressionHere $\gamma$ is the solid-liquid surface energy, $T_m$ in this expression is absolute temperature in kelvin, and $\mathcal{K}$ is the sum of principal curvatures, positive for solid convex into liquid. For this graph, $\mathcal{K}=-{\zeta }_{yy}={\alpha }^2\zeta$. If mean curvature is defined as half this sum, $\mathcal{K}=2H$; this fixes the capillary factor unambiguously.

Under the small-wavelength thermal approximation for a solidification front, take harmonic perturbation temperatures that decay away from the interface:This approximates their temperature equation by Laplace's equation, neglecting perturbation time dependence and the moving-frame advection. Expanding the interfacial temperature at the displaced surface givesThe second expression includes displacement through the negative base liquid gradient. Linearizing the Stefan condition also requires evaluation of the base gradient at the displaced surface:Using $S=1$ gives the requested reduced dispersion relation:orThe destabilizing $V\alpha$ term describes a protrusion entering colder liquid. The cubic term expresses stabilization by the curvature-induced temperature decrease.

There is an important accuracy qualification. The harmonic approximation gives the leading short-wave thermal exponents, but the retained constant term is not a uniformly controlled next-order correction to the full moving-boundary problem. If perturbation diffusion is retained, its exponents areSubstitution into exactly the same linearized boundary balances gives the full thermal stability of a travelling solidification front:For example, with capillarity removed, $\sigma =V\alpha$ solves this relation exactly, whereas the harmonic reduction gives $V\alpha -V^2/\kappa$. The difference is smaller than the leading term when $\alpha \gg V/\kappa$, but demonstrates why the reduced formula should not be asserted as an exact dispersion law or a precise onset calculation. To neglect perturbation time dependence one also needs $|\sigma |\ll \kappa {\alpha }^2$; extremely large capillary damping is outside that approximation.

To analyse the displayed cubic, defineThen $\kappa \sigma /V^2=g(q)=-1+q-2q^3/R_c$. It is negative at zero and at sufficiently large $q$. Its positive stationary point is the unique maximum:Thus the capillary threshold of a reduced solidification dispersion iswith equality giving a neutral maximum and strict inequality giving a band of unstable positive wavenumbers. At equality, however, $q_{\ast }=3/2$, not $q_{\ast }\gg 1$. The number $27/2$ is therefore the threshold of the stated reduced cubic, not a controlled threshold of the full thermal problem under the stated short-wave ordering. Far above this threshold its fastest-growing scale does satisfy $q_{\ast }\gg 1$, making the short-wave mechanism meaningful.

The analogy for a snowflake is dendritic growth of a snowflake. A protrusion has better access to the driving supersaturation or undercooling and grows faster than a sheltered region; this amplifies tips and promotes branching into crystal dendrites. Interfacial capillarity suppresses arbitrarily fine tips. In an actual snowflake, the main growth is deposition from water vapour, so the diffusion field and its driving variable differ from this pure-melt thermal example. The hexagonal ice crystal anisotropy supplies six equivalent preferred arm directions; an isotropic planar stability calculation alone cannot predict six arms. Nonlinear tip selection, vapour transport and side branching are needed for the later detailed morphology. The analysis explains the competition between diffusion-driven protrusion growth and capillary smoothing, with this physically necessary qualification.

Melting raises permeability of a porous medium. An advancing protrusion therefore offers a less resistive path to hot liquid, attracting larger Darcy flux and receiving more heat. The resulting faster melting reinforces the protrusion: this is a reactive infiltration instability. A permeability decrease would reverse this feedback, while zero contrast gives no growth in this ideal model.

For positive contrast, the growth rate is proportional to transverse wavenumber, so the model has no finite fastest-growing wavelength. The latent-to-sensible heat ratio slows the front but supplies no short-wave cutoff. Finite thermal transport and a resolved melting zone can introduce a length scale; pore geometry eventually limits the continuum description. If interface curvature changes the equilibrium melting temperature through the Gibbs--Thomson relation, that can also oppose short-wave corrugations. The modified dispersion relation must include such effects before a preferred wavelength can be calculated; a length scale alone does not guarantee that every added mechanism selects one.

Apply force balance to one hemisphere. The pressure force is $\pi R^2({\Pi }_{\mathrm{in}}-{\Pi }_{\mathrm{out}})$, while the interfacial tension around its equator gives $2\pi R\gamma$. The Young–Laplace equation therefore yieldsTo leading order in curvature, let the two bulk compositions be ${\varphi }_{+}={\varphi }_B+{\delta }_{+}$ and ${\varphi }_{-}=-{\varphi }_B+{\delta }_{-}$. Their curvatures agree, $\alpha =f^{\prime \prime }(\pm {\varphi }_B)=-2a$, so equal chemical potentials require $\alpha {\delta }_{+}=\alpha {\delta }_{-}$ and hence a common shift $\delta$.

Since $d\Pi /d\varphi =\varphi f^{\prime \prime }(\varphi )$, expansion about the flat coexistence compositions gives ${\Pi }_{+}-{\Pi }_{-}=2{\varphi }_B\alpha \delta +O({\delta }^2)$. Matching the capillary pressure proves the Gibbs--Thomson relation in these variables:The pressure jump balances interfacial forces; it does not permit unequal chemical potentials. A common shift in both compositions supplies a common interfacial chemical potential while their pressures differ. The expansion assumes $R$ large compared with the phi-four diffuse interface width.

The Young–Laplace equation gives a pressure excess $2\sigma /R$ inside a spherical droplet. Near coexistence, the common chemical-potential shift is $\mu \simeq \alpha \delta$, where $\alpha =f^{\prime \prime }({\varphi }_B)$, while changing phase changes the composition by approximately $2{\varphi }_B$. Balancing the capillary pressure against this thermodynamic shift gives the Gibbs--Thomson relationHence the compositions immediately outside and inside arewith

Write the base-state gradients asThe base heat flux and Stefan condition relations are

Under the quasi-stationary approximation, and with advective heat transport neglected, each temperature perturbation satisfies the Laplace equation. The decaying normal-mode solutions areKeeping the displaced ice interface at $T_m$ omits the curvature correction described by the Gibbs--Thomson relation.
The fixed melting temperature at the displaced ice interface givesExpanding temperature and conductive heat flux continuity at the displaced outer interface gives

Put $r=k_a/k_w$. Since $G_a=G_w/r$, the temperature condition in the limits $r\ll 1$ and $\alpha h\ll 1$ givesAlthough $r$ is small, the product $rC$ is therefore order one and cannot be discarded. The flux condition then givesThis is why the stated asymptotic warning matters.

The perturbed Stefan condition at $x=0$ isUsing $B\sim {\eta }_2{G}_w$ and $-k_w{G}_w={\rho }_i{L}_{f}V$ yieldsSubstitution of the interface-amplitude ratio from part i produces the thin-film icicle-ripple instability dispersion relation

Set $b=\Gamma h{\alpha }^3$. Rationalizing this complex number givesso its growth rate and imaginary part areDifferentiating the growth rate with respect to the wavenumber shows that its positive stationary point satisfiesIt is the unique global maximum, and henceAt this wavenumber $b=1/\sqrt{5}$. With the convention $e^{i\alpha z+\sigma t}$, a constant phase travels with phase velocityThereforeBecause $z$ increases downward, the negative sign means that the icicle ripples migrate upward with speed $\sqrt{5}V/6$.

For $a<0$ and $b>0$, planar coexistence occurs at $\varphi =\pm {\varphi }_B$ with ${\varphi }_B=\sqrt{-a/b}$. Put $\alpha =f^{\prime \prime }({\varphi }_B)=a+3b{\varphi }_B^2=-2a>0$. A common small shift $\delta$ gives equal chemical potentials to first order,The thermodynamic pressure is $P=\varphi \mu -f$. Its difference between the positive interior and negative exterior is thereforeFor a sphere the Young–Laplace equation gives $P_{+}-P_{-}=2\sigma /R$. Consequently the Gibbs--Thomson relation is

Approximating perturbation temperatures by harmonic fields decaying as $e^{\pm \alpha \xi }$, while retaining displacement through the base gradient, produces this reduced dispersion relation. The Gibbs--Thomson relation supplies $\Gamma =\gamma T_m/(\rho L)$. This short-wave reduction suppresses perturbation advection and time dependence; its constant term and a threshold at $\alpha \kappa /V=O(1)$ are not controlled next-order results of the full diffusion problem.

The attachment mobility $\mathcal{M}>0$ is the coefficient converting interfacial undercooling into a solidification speed. In the linear isotropic law, $U=\mathcal{M}(T_{\mathrm{eq}}-T_i)$, and capillarity makes $T_{\mathrm{eq}}=T_m-{\gamma }_T\mathcal{K}$. Thus $T_i=T_m-U/\mathcal{M}-{\gamma }_T\mathcal{K}$. This kinetic coefficient refers to a moving phase boundary, rather than an order-parameter relaxation coefficient. The Gibbs--Thomson relation determines ${\gamma }_T=T_m\gamma /(\rho L)$ for isotropic solid-liquid surface energy $\gamma$.

The excess free energy per area of a solid-liquid phase boundary. Curvature changes the equilibrium temperature through the Gibbs--Thomson relation. Crystal anisotropy can make this energy and its associated surface stiffness depend on orientation.