Codex @codex  0

For the rapidly varying wave phase, define $k_i={\partial }_i\theta$ and $\omega =-{\partial }_t\theta$, with the common scaling by $\epsilon$ absorbed into these local variables. A packet moves with group velocity $c_{g,i}=\partial \Omega /\partial k_i$. Along that trajectory, $d/dt={\partial }_t+c_g\cdot {\nabla }_x$. Compatibility gives ${\partial }_t{k}_i=-{\partial }_i\omega$ and ${\partial }_i{k}_j={\partial }_j{k}_i$. Differentiating the local dispersion relation then cancels its implicit $k$ derivatives, yieldingThe last equality follows also by differentiating $\omega =\Omega (k,x,t)$ along the ray and cancelling the two Hamiltonian cross terms. The partial derivatives of $\Omega$ hold its other arguments fixed.

A mean flow adds the Doppler shift: the intrinsic frequency is $\hat{\omega }=\omega -kU(z)$, soon the positive-intrinsic-frequency branch launched here. The medium is stationary and independent of $x$, so $\omega$ and $k$ are constant. Initially $m=0,U=0$, hence $\omega =N$. Moreover $\dot{m}=-k{U}^{\prime }<0$, and then $\dot{z}=-Nkm/(k^2+m^2{)}^{3/2}>0$. The initial vertical group velocity is zero, but $m$ immediately becomes negative and the ray accelerates upwards.

If $U$ reaches $N/k$ at a finite height $z_c$, this is the critical level of an internal gravity wave, with $U(z_c)=\omega /k$. Below it,Assuming $0<U^{\prime }(z_c)<\infty$, put $\delta =z_c-z$. Then $N-kU\sim k{U}^{\prime }(z_c)\delta$, givingThus the ray cannot cross this level in finite time within the geometrical-optics approximation. The hypothesis $U^{\prime }>0,U(0)=0$ alone does not guarantee a finite $z_c$: for example $U(z)=\frac{N}{2k}(1-e^{-z})$ on $z\ge 0$ never reaches $N/k$. Existence of a finite critical level is an additional assumption. Very close to one, diverging wavenumber can invalidate the inviscid, linear, slowly varying approximation; the result is the formal ray prediction.

Put $D={\partial }_t+U{\partial }_x$, write $p$ for perturbation pressure divided by ${\rho }_0$, and set $b=-g{\rho }^{\prime }/{\rho }_0$. The linearized Boussinesq approximation, incompressibility and material derivative giveTaking the curl eliminates pressure; differentiating the resulting vorticity equation in $x$ and using incompressibility gives $D{\nabla }^{2}w-U^{\prime \prime }{w}_x=b_{xx}$. A further application of $D$ therefore yieldsFor a stationary nonzero horizontal Fourier mode, $D=U{\partial }_x$, so cancellation of ${\partial }_x^2$ gives the stationary Taylor–Goldstein equationThis division requires $U\ne 0$ on the interval considered: a zero of $U$ is a critical level of an internal gravity wave. The mean profiles must be sufficiently smooth for the displayed derivatives, with background hydrostatic pressure and stable density stratification, $N^2>0$, for propagating internal gravity waves. The horizontally uniform component is excluded from the cancellation. For a horizontal wavenumber $k$, local vertical propagation additionally requires $N^2/U^2-U^{\prime \prime }/U>k^2$; this is distinct from the smoothness restrictions. If stability of the background against other disturbances is needed, the Miles–Howard theorem supplies the sufficient condition $N^2\ge (U^{\prime }{)}^2/4$, rather than a necessary condition for deriving the equation.

Take $U_0,N_0,k>0$ and assume an initially propagating mode, $k<N_0/U_0$. Put $s=N_0/U_0$ and $\delta (z)=1+\gamma z$. Neglecting $U^{\prime \prime }/U$ as stipulated, the Scorer parameter and stationary internal-wave WKB solution areThe positive root $m>0$ selects upward group velocity for the negative intrinsic frequency branch $\hat{\omega }=-kU$. The WKB approximation needs a slowly varying background relative to the local vertical wavelength, in particular $|m^{\prime }|\ll m^2$ and $|m^{\prime \prime }|\ll |m{|}^3$, away from a turning level or singular mean profile.

The local dispersion relation is $\omega =kU-Nk/(k^2+m^2{)}^{1/2}=0$. Its intrinsic and laboratory group velocities are, with $K=(k^2+m^2{)}^{1/2}$,Thus intrinsic energy propagation is upstream and upward, while laboratory energy propagation is downstream and upward, parallel to $(k,m)$. Constant-phase lines of an internal gravity wave have slope $dz/dx=-k/m$ and normal spacing $2\pi /K$; the horizontal spacing remains $2\pi /k$.

If $\gamma <0$, $m$ grows with height. The wavevector and laboratory energy ray become more vertical; phase lines become flatter and closer together. For finite nonzero $N$, the heightis a critical level of an internal gravity wave: $U\to 0$ and the vertical wavenumber diverges. The laboratory group velocity tends to zero. Importantly, $|m^{\prime }|/m^2\to |\gamma |/s$, so the formal WKB condition need not deteriorate before this level if that ratio is small. Nevertheless the stationary equation is singular at the level, the phase has no finite limit, and the horizontal perturbation velocity grows as $m^{1/2}$; the inviscid linear wave cannot be continued uniformly through it without additional physics.

If $\gamma >0$, $m$ decreases to zero atThe wavevector and laboratory ray become horizontal, the phase lines become vertical, and the normal spacing increases. Above $z_t$, $m^2<0$ and the disturbance is an evanescent wave. Here $|m^{\prime }|/m^2$ diverges, so the WKB approximation fails in a turning region and must be replaced by a local connection solution.

Use a common factor $e^{i(kx-\omega t)}$, with $k=k_0>0$ and $\omega ={\omega }_0$. For the upward propagating branches with $0<{\hat{\omega }}_2=\omega -k{U}_2<N_2$, defineThese are the inclinations of upward-sloping constant-phase lines of an internal gravity wave to the horizontal. Write the incident, reflected and transmitted complex velocity amplitudes as $A_0,A_r,A_t$. Incompressibility gives $(u_i,w_i)=A_0(\cos {\theta }_1,\sin {\theta }_1)$, $(u_r,w_r)=A_r(\cos {\theta }_1,-\sin {\theta }_1)$, and $(u_t,w_t)=A_t(\cos {\theta }_2,\sin {\theta }_2)$.

The jump conditions for stratified inviscid shear flow require a common interface displacement, not a common vertical velocity: $w=-i\hat{\omega }\eta$ on each side. Since the background mass density is continuous, pressure is continuous without a hydrostatic jump. The horizontal momentum equation gives $p/{\rho }_0=(\hat{\omega }/k)u$. HenceThese two matching equations determine internal-wave transmission across a velocity jump:The amplitude formula printed in the PDF is inconsistent with these material-interface matching conditions. In particular, identical layers have $A_t/A_0=1$ and $A_r/A_0=0$, whereas its expression labelled “reflected” equals one. The displayed results above distinguish reflection from transmission and retain the intrinsic-frequency factors required by the kinematic boundary condition and pressure balance.

For equal buoyancy frequencies and the specified opposing current, ${\theta }_1=\pi /6$, ${\theta }_2=\arcsin (3/4)$, andThe incident phase lines have slope $1/\sqrt{3}$, the reflected lines slope $-1/\sqrt{3}$, and the transmitted lines slope $3/\sqrt{7}$: they steepen above the interface. For the other specified current, ${\hat{\omega }}_2=0$. The transmitted vertical wavenumber then diverges and there is no regular propagating upper-layer wave of that frequency: this is the critical level of an internal gravity wave limit. Approaching it from ${\hat{\omega }}_2>0$ gives $A_r/A_0\to -1$ and vanishing transmitted vertical energy flux; setting the intrinsic frequency to zero directly is outside the regular matching calculation.