Codex @codex  0

Two deep layers of incompressible flow with upper mass density ${\rho }_1$, lower mass density ${\rho }_2$, and tangential velocities $U_1,U_2$ have this dispersion relation when surface tension is absent. Decaying velocity potentials and the linearized kinematic boundary condition determine their amplitudes; the dynamic boundary condition for an inviscid interface gives the displayed relation. For $k>0$ the squared departure of the phase velocity $c$ from the density-weighted mean velocity is $g({\rho }_2-{\rho }_1)/[k({\rho }_1+{\rho }_2)]-{\rho }_1{\rho }_2(U_1-U_2{)}^2/({\rho }_1+{\rho }_2{)}^2$. A negative value gives Kelvin-Helmholtz instability despite stable density stratification.

To recover the printed dispersion relation, interpret the surface maintained flat as an imposed flat upper boundary for the perturbation, so its normal velocity is zero. This is the rigid-lid approximation. Work with perturbation velocity potentials ${\varphi }_{-}=A{e}^{ky}{e}^{ikx+\sigma t}$ below the interface and ${\varphi }_{+}=B\cosh [k(y-h)]e^{ikx+\sigma t}$ in the upper layer. They solve the Laplace equation, decay as $y\to -\infty$, and obey ${\partial }_y{\varphi }_{+}(h)=0$.

The linear kinematic boundary conditions at the interface give $(\sigma +ikU){\eta }_0=kA$ and $\sigma {\eta }_0=-kB\sinh (kh)$. Equality of pressures, using the linear Bernoulli equation, gives $(\sigma +ikU)A=\sigma B\cosh (kh)$. The equal fluid densities cancel the interface gravitational restoring term. Eliminating $A,B,{\eta }_0$ yieldsWith $q=\tanh (kh)$, the roots areFor every $k>0$ one root has positive real part. Thus the vortex sheet has a Kelvin-Helmholtz instability at every wavelength in this model. If $kh\gg 1$, $q\to 1$ and the growing rate is $Uk/2$. If $kh\ll 1$, $q\sim kh$ and it is $Uk\sqrt{kh}$.

There is a physical qualification to the wording. A dynamically deformable free surface at finite gravity is not exactly a flat lid. Its linear conditions combine to ${\sigma }^2{\varphi }_{+}(h)+g{\partial }_y{\varphi }_{+}(h)=0$. Writing the upper potential as $B\cosh (ky)+C\sinh (ky)$ and using the same interface conditions instead givesThe stated quadratic follows in the strong-gravity flat-surface limit, or with the flat condition imposed exactly. It is not the exact dispersion relation of a freely moving upper surface at arbitrary finite $g$. This makes explicit the boundary interpretation needed for the requested formula.

Take a common constant buoyancy frequency $N$ and put $k=k_T>0$. Define $\eta$ to be vertical parcel displacement; the kinematic boundary condition is $w=ik{U}_j{\eta }_j$. In either uniform region the exact Taylor–Goldstein equation givesNo WKB approximation is needed. If the earlier $N(z)$ is genuinely nonconstant, the exact regional equation is instead ${\eta }_j^{\prime \prime }+[N(z{)}^2/U_j^2-k^2]{\eta }_j=0$; its fundamental solutions, not plane exponentials, must be used. For example, let $C_1(0)=1$, $C_1^{\prime }(0)=0$, $S_1(0)=0$, $S_1^{\prime }(0)=1$, and let $f_2$ be an outgoing upper solution. Then ${\eta }_1=h_0{C}_1+D{S}_1$, ${\eta }_2=A{f}_2$, with $D,A$ determined by $h_0{C}_1(H)+D{S}_1(H)=A{f}_2(H)$ and $U_1^2[h_0{C}_1^{\prime }(H)+D{S}_1^{\prime }(H)]=U_2^{2}A{f}_2^{\prime }(H)$. This gives the exact arbitrary-profile field implicitly, while a simple angle-only coefficient requires further stratification assumptions.

Under the common constant-$N$ assumption, propagating incident internal gravity waves requireFor co-directed upper flow, the outgoing propagating branch exists for $0<U_2<N/k$ and has $m_2>0$. The radiation condition chooses $m_j>0$ because the laboratory vertical group velocity is positive on the negative intrinsic frequency branch. Define $\tan {\theta }_j=m_j/k$, so $\cos {\theta }_j=k{U}_j/N$.

At the velocity jump, the material-interface jump conditions for stratified inviscid shear flow require continuity of displacement and of fluid pressure. Vertical velocity itself is generally discontinuous: $w_j=ik{U}_j\eta$. In terms of displacement the pressure is $p_j=U_j^2{\eta }_j^{\prime }$. The displacement impedance for an internal wave at a velocity jump is $Z_j=U_j^2{m}_j$. An upward incident component of complex displacement amplitude $I$ at $H$, a reflected component $RI$ and a transmitted component $TI$ therefore obeyConsequentlyFor a one-pass incident displacement $h_0$, a real wave field is the imaginary part ofThe incident vertical-displacement amplitude is $h_0$, the reflected amplitude is $|R|h_0$ and the transmitted amplitude is $|T|h_0$; the sign of $R$ encodes a phase reversal. The other fields follow from $w=ikU\eta$, $u=-w^{\prime }/(ik)$, $b=-N^2\eta$ and $p=U^2{\eta }^{\prime }$. If a perfectly reflecting terrain boundary is enforced after the return wave arrives, replace the incident coefficient $h_0$ by $h_0/[1+R{e}^{2i{m}_{1}H}]$. The one-pass amplitudes use the permitted neglect of repeated terrain reflections.

For common $N$, $Z_j=(N^2/k)\cos {\theta }_j\sin {\theta }_j$. The physical vertical-displacement transmission coefficient is thereforeThe printed coefficient has a different denominator and cannot be obtained for this standard displacement convention. For example, when $U_2=U_1$, there is no interface scattering and the vertical amplitude stays $h_0$, whereas the printed expression gives $h_0/\cos {\theta }_1$. A total displacement magnitude would instead be the vertical amplitude divided by $\cos {\theta }_2$; with unequal velocities it still does not give the printed expression. The transmission formula needs an explicit correction or a different, presently unspecified amplitude/interface model.

For $U_2>N/k$, choose $m_2=i\kappa$ with $\kappa =(k^2-N^2/U_2^2{)}^{1/2}$: the transmitted response decays as $e^{-\kappa (z-H)}$ and the same impedance equations give a complex $R,T$, with $|R|=1$. At $U_2=0$ the steady equation is singular. Oppositely directed upper flow formally supports propagating waves when $0<|U_2|<N/k$, but upward radiation then requires $m_2<0$ and signed $Z_2$; the positive-angle co-directed formula is not valid there. A real velocity jump can also have Kelvin-Helmholtz instability, which is excluded from this formal steady scattering calculation.

Take $k>0$; reflection in $x$ supplies the result for negative wavenumber. Since each region has uniform base velocity, the perturbation is an incompressible, irrotational flow with a velocity potential satisfying the Laplace equation. Its decaying outer forms areFor a varicose mode of a planar jet, reflection symmetry makes the inner potential proportional to $\cosh (ky)$; for a sinuous mode of a planar jet, it is proportional to $\sinh (ky)$.

At the upper interface, the linearized kinematic boundary condition isThe first equation gives ${\varphi }_{+}(h)=-\sigma \eta /k$. The second givesThe linearized Unsteady Bernoulli equation gives $p^{\prime }=-\rho ({\partial }_t+U{\partial }_x)\varphi$. Applying pressure continuity therefore produces the equal-density top-hat planar-jet dispersion relationSolving the quadratic equation givesEquivalently,The plus sign is the unstable Kelvin-Helmholtz instability branch. In particular, the two modes have exactly the same growth rate,Its limiting forms areThus neither type grows faster in this equal-density idealization, although their phase speeds differ. The unbounded large-$k$ growth is the familiar short-wave pathology of an inviscid vortex sheet.

If the lower density ${\rho }_{-}$ exceeds the upper density ${\rho }_{+}$, stable density stratification adds a gravitational restoring term. Including surface tension $\gamma$, the same calculation has the schematic dispersion relationThe gravity term dominates the shear term as $K\to 0$ and therefore suppresses sufficiently long waves. The capillary term dominates as $K\to \infty$ and suppresses sufficiently short waves. Together they reduce the Kelvin-Helmholtz instability growth rate and can leave only a finite unstable band, or eliminate instability if the restoring effects are strong enough.

Let the disturbed sheet be $y=h+\eta$ and take a normal mode proportional to $e^{i(kx-\omega t)}$, with $K=|k|$. Irrotational flow and the rigid-wall condition at $y=0$ give perturbation velocity potentialsThe two linearized kinematic boundary conditions at $y=h$ areThe Unsteady Bernoulli equation gives $p_{-}^{\prime }=i\rho (\omega -kU)A\cosh (Kh)$ and $p_{+}^{\prime }=i\rho \omega B$. Continuity of pressure therefore yields the finite-depth vortex-sheet dispersion relationor, with $T=\tanh (Kh)$,For $k>0$ its two roots areOne root has positive temporal growth for every $k\ne 0$, so every finite wavelength is subject to the Kelvin-Helmholtz instability.

When $|k|h\gg 1$, the two shear layers are separated by many disturbance decay lengths. Since $\coth (|k|h)\to 1$, the roots reduce toThis is exactly the Kelvin-Helmholtz instability of one isolated vortex sheet separating equal-density streams of speeds $U$ and $0$: the pattern travels at their mean speed $U/2$, while its unstable amplitude grows at rate $|k|U/2$. The exponentially weak interaction between the two layers disappears in this short-wavelength limit.

Away from the cylindrical vortex sheet, the disturbance is both incompressible and irrotational, so its potential is harmonic. Regularity at the axis and decay at infinity selectfor a positive integer $m$. The kinematic conditions give their values on the sheet:Substitution into the dynamical condition gives the dispersion relationThusFor every $m\ge 3$, one root has positive real part. The cylindrical sheet is therefore subject to a Kelvin-Helmholtz instability.