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

Suppose a $C^1$ classical solution existed and set $w=u+iv$. The system is preciselywhich are the Cauchy-Riemann equations for $w$ as a function of the complex number $z=x+it$. Hence $w$ is holomorphic near $z=0$. Every holomorphic function is real analytic, so its restrictionto the real axis is real analytic near zero. Its real part and imaginary part show that both $f$ and $g$ must be real analytic, contradicting the hypothesis. Therefore no such $C^1$ solution exists.

For an unstable Taylor–Goldstein equation mode in a finite channel, take $a=1$ in the power-transformed Taylor–Goldstein energy identity. Its imaginary part givesThus the real phase velocity lies strictly between the extremes of a nonconstant smooth shear profile. Equality would force the regular eigenfunction to vanish on an interval, hence everywhere by uniqueness for its ordinary differential equation.