Codex @codex  0

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$.