Codex @codex  0

For the positive-frequency internal gravity wave branch, $\omega =N{K}_h/K$. Differentiating with respect to the components of the wave vector gives the group velocity,The wavefront-normal phase velocity is ${\mathbf{c}}_p=\omega \mathbf{k}/K^2$. Directly,The negative-frequency branch changes both propagation signs but not their orthogonality. Another proof is that $\omega$ is homogeneous of degree zero in the wave vector, so Euler's homogeneous-function identity gives $\mathbf{k}\cdot {\nabla }_{\mathbf{k}}\omega =0$. Energy propagates along constant-phase surfaces, rather than normal to them. When $m=0$, $\omega =N$ and ${\mathbf{c}}_g=0$; the scalar-product result remains true, but there is no nonzero group-velocity direction to describe geometrically.

Use the nonrotating linear Boussinesq approximation with reference density ${\rho }_0$ and $N^2=-(g/{\rho }_0)d\bar{\rho }/dz>0$:Let $K_h^2=k^2+l^2>0$, $K^2=K_h^2+m^2$ and $E=e^{i(kx+ly+mz-\omega t)}$. The horizontal momentum equations give $u^{\prime }=k{p}^{\prime }/({\rho }_0\omega )$ and $v^{\prime }=l{p}^{\prime }/({\rho }_0\omega )$. Continuity requires $k{u}^{\prime }+l{v}^{\prime }+m{w}^{\prime }=0$. Thus the internal-wave polarization isPhysical perturbations are the real parts. The density is in temporal quadrature with the vertical velocity, because buoyancy responds to the vertical fluid displacement. Substitution into the vertical momentum equation determines the internal gravity wave dispersion relation,A nontrivial propagating wave requires a nonzero horizontal wave vector and a nonzero frequency. The formulas are not to be divided by $K_h^2=0$: incompressibility excludes a nonzero oscillatory vertical velocity for a purely vertical wave vector.

Use the inviscid, nonrotating Boussinesq approximation, with a stable background mass density $\bar{\rho }(z)$ and reference density ${\rho }_{\ast }$. Write $B=-g{\rho }^{\prime }/{\rho }_{\ast }$ for the perturbation buoyancy and $\pi =p^{\prime }/{\rho }_{\ast }$ for the kinematic pressure. The buoyancy frequency is $N_0^2=-(g/{\rho }_{\ast }){\bar{\rho }}_z>0$. Dropping products of perturbations in the Boussinesq equations gives the Linearized Boussinesq equationsThe last equation expresses incompressible flow; the second follows by advecting the background mass density gradient. Neglecting rotation and viscosity is part of this internal gravity wave model.

For a plane internal gravity wave proportional to $\exp [i({\mathbf{k}}_h\cdot {\mathbf{x}}_h+mz-\omega t)]$, put $\kappa =|{\mathbf{k}}_h|>0$ and $K^2={\kappa }^2+m^2$. Eliminating the horizontal velocity, pressure and buoyancy from the linear equations givesFor example, the horizontal momentum and continuity equations give $\pi =-\omega mw/{\kappa }^2$; substituting $B=-i{N}_0^{2}w/\omega$ into vertical momentum yields the displayed dispersion relation. Thus an IGW has $0<|\omega |\le N_0$ in this model.

On the positive-frequency branch, the phase velocity normal to a constant-phase plane and the group velocity areConsequently $({\mathbf{k}}_h,m)\cdot {\mathbf{c}}_g=0$, and phase and group velocity are perpendicular. Equivalently, the dispersion relation is homogeneous of degree zero in the wave vector, so differentiating with respect to its scale proves the same orthogonality. Energy travels with the group velocity, along the phase planes. At the degenerate limit $m=0$, $\omega =N_0$ and the group velocity vanishes; the orthogonality statement then has this limiting interpretation.

For a two-dimensional internal gravity wave with fixed horizontal wave number $k\ne 0$, the vertical velocity amplitude satisfiesThis follows from the Linearized Boussinesq equations just as in part (a), allowing the background buoyancy frequency to vary. The WKB approximation gives locallyIt requires slowly varying stratification compared with the local vertical wavelength, in particular $|m_z|\ll m^2$, and a small-amplitude wave in the inviscid, nonrotating regime. Propagation requires $N>|\omega |$; where $N<|\omega |$, the local vertical wave number is imaginary and the solution is evanescent. Turning levels and abrupt changes need a separate connection calculation. Finite viscosity, background flow, rotation or nonlinear effects require modifications of this model.

The ray tracing equations use the local dispersion relation as a Hamiltonian:The stationary, horizontally uniform background conserves the frequency and horizontal wave number. In two dimensions the group velocity gives $dx/dz=-m/k$. For the specified buoyancy frequency and $\omega =1$, $m^2=256k^2{z}^2$. Hence for $z>0$ the two geometric internal-wave rays through the specified point areEach curve has both directed propagation senses, giving four local directed rays. Their formal intersections with $z=0$ are $x=-8$ and $x=8$.

At $z=0$, $N=\omega$ and $m=0$: this is a double-zero internal-wave turning level, since $m^2$ touches zero without changing sign. The ray becomes vertical geometrically, but its group velocity tends to zero. Taking $k>0$ without loss of generality, the downward branch hasThe formal travel time diverges logarithmically as $z\to 0$. Meanwhile $|m_z|/m^2=1/(16|k|z^2)$ grows, so the WKB approximation fails before that limit; its predicted $|m{|}^{-1/2}$ amplitude cannot be extrapolated to infinity. The spatial ray reaches a finite limiting position, but finite-wavelength behaviour near the zero must be found beyond WKB. There is no evanescent half-space supplied by this particular $N(z)$, so ordinary simple-turning-point reflection cannot be assumed solely from the local ray construction.

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.

Use ${\varphi }_j=U_{j}x-\frac{1}{2}{U}_j^{2}t+{\varphi }_j^{\prime }$ and retain terms linear in the displacement and perturbation potentials. Evaluation at $z=\eta$ may then be replaced by evaluation at $z=0$: the correction to a perturbation is second order. The kinematic boundary conditions becomeThe linearized pressure continuity isThe decaying normal modes give $-k{A}_1=ik(U_1-c)B$ and $k{A}_2=ik(U_2-c)B$. Thus $A_1=-i(U_1-c)B$, $A_2=i(U_2-c)B$. Substitution into the dynamic condition yields the Kelvin-Helmholtz dispersion relation with gravity,Put $\rho ={\rho }_1+{\rho }_2$ and $\overline{U}=({\rho }_1{U}_1+{\rho }_2{U}_2)/\rho$. Completing the square givesFor $k>0$, the temporal growth rate is $k{c}_i$. A growing complex conjugate branch therefore exists precisely when the radicand of the square root is negative:At equality the two phase velocities coincide and there is no exponential growth. Below it the roots describe real internal gravity waves. The denser lower fluid gives stable density stratification, but the tangential shear can overcome it. The assumption $U_1\ne U_2$ is needed for the displayed threshold; equal streams have no such shear instability.

Let $p$ denote the fluid pressure perturbation divided by the reference mass density, and let $b=-g{\rho }^{\prime }/{\rho }_0$ be the buoyancy perturbation. The buoyancy frequency satisfies $N^2=-g{\hat{\rho }}_z/{\rho }_0$. Linearization requires small boundary slope $k|{\eta }_0|\ll 1$, small displacements compared with the vertical scales of the disturbance and background, and advection small compared with the oscillatory acceleration. For a propagating internal gravity wave this includes $|m{\eta }_0|\ll 1$; boundary slope alone is insufficient when $|m|\gg k$. The Boussinesq approximation also requires small relative mass density differences over the region of interest. These conditions must hold for the resulting disturbance, including any amplification by resonance.

The inviscid Linearized Boussinesq equations and linear kinematic boundary condition areWrite the velocity field as $(u,w)=\mathrm{Re}\{(U(z),W(z))e^{i(kx-\omega t)}\}$. Incompressibility and horizontal momentum balance giveThe vertical momentum balance therefore yieldsFor a real vertical wavenumber $m$, this gives the dispersion relationWhen $0<\omega <N$, define $\ell =k\sqrt{N^2/{\omega }^2-1}$. The radiation condition selects $m=-\ell$, because energy must leave the boundary upward. The boundary-forced internal gravity wave isIts wavevector is $\mathbf{K}=(k,m)$, its phase velocity is ${\mathbf{c}}_p=\omega \mathbf{K}/(k^2+m^2)$, and its group velocity isThus the phase velocity points down and right, while the group velocity points up and right. They are perpendicular. Internal-wave polarization makes the particle motion an oscillation along the group velocity direction, with zero first-order mean transport. The constant-phase lines of an internal gravity wave are also parallel to the group velocity. Their inclination $\vartheta$ above the horizontal satisfies $\sin \vartheta =\omega /N$.

When $\omega >N$, set $\kappa =k\sqrt{1-N^2/{\omega }^2}$. Boundedness at infinity selects the evanescent waveThe horizontal and vertical velocity components are in quadrature: fluid particles describe small ellipses, and the response decays over ${\kappa }^{-1}$. There is no upward time-averaged energy flux, because $P$ and $W$ are in quadrature. A real vertical group velocity is not defined for this evanescent wave. The pattern travels horizontally with phase velocity $\omega /k$.

For $N=0$ the response is the evanescent wave with $\kappa =k$; the ratios $\omega /N$ should not be used. At the cutoff $\omega =N>0$, the bounded harmonic solution has $W=-i\omega {\eta }_0$, $U=0$: it neither decays nor has nonzero upward group velocity. It is the limiting cutoff response, rather than a localized radiating disturbance. A uniform nonzero buoyancy frequency in an infinitely deep Boussinesq approximation is itself a local idealization of the background mass density.

The arrows distinguish the phase velocity, group velocity, and oscillatory particle motion of the internal gravity wave. The right panel shows the decay envelope and particle ellipses of the evanescent wave.

At the change in buoyancy frequency, continuity of vertical velocity and fluid pressure givesThe background mass density is continuous, so there is no interfacial density jump or separate surface restoring force. Since $P=i\omega W^{\prime }/k^2$, these internal-wave transmission across a stratification step conditions are continuity of $W$ and $W^{\prime }$. The buoyancy perturbation need not be continuous, because the background density gradient jumps.

Let $W_0=-i\omega {\eta }_0$ and $\ell =k\sqrt{N_2^2/{\omega }^2-1}$. The lower unstratified layer obeys $W^{\prime \prime }-k^{2}W=0$ and the upper radiation condition givesPropagating these matching data downward givesThe kinematic boundary condition at the moving lower boundary therefore givesThe upper vertical-displacement complex amplitude at $H$ isIn particular, $|{\eta }_H/{\eta }_0|=[{\cosh }^2(kH)+(\ell /k{)}^2{\sinh }^2(kH){]}^{-1/2}$. The unstratified layer attenuates and phase-shifts the transmitted internal gravity wave. The matching height is $H$ in the original PDF; the TeX transcription's $H^{\prime }$ is a typo.

The central stratified layer supports upward and downward internal gravity waves, whose counterpropagating components form a standing wave. The unstratified regions support evanescent waves, and the upper disturbance decays as $e^{-k(z-2H)}$. Consequently the central layer is a stratified internal-wave guide, rather than a source of propagating energy at infinity.

The largest response occurs at its trapped normal mode frequencies. They can be specified without solving the forced problem. Set the boundary forcing to zero to find a free normal mode. The lower solution is proportional to $\sinh (kz)$, so continuity of fluid pressure and vertical velocity givesIn the central layer writeThe upper Robin boundary condition then gives the exact trapped-mode conditionThe left side increases strictly from zero to infinity, so there is one positive root for each $n$; the frequencies accumulate at zero. This is constructive phase matching after reflection at both ends. The phase shifts from the evanescent waves matter: simply imposing integer half-wavelengths across the stratified layer is generally incorrect.

In the ideal inviscid model, exact resonant forcing has no bounded steady harmonic solution. The undamped normal mode grows secularly under sustained forcing. Weak viscosity or other losses would produce large finite peaks near the displayed frequencies; the linear approximation eventually fails if the disturbance becomes too large.

The illustrated free normal mode has an oscillatory central region and evanescent wave tails. Its lower tail reaches the fixed zero-displacement boundary; its upper tail decays to infinity.

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.

Let every perturbation be proportional to $\exp [i(kx+mz-\omega t)]$ and define the intrinsic frequencyEliminating pressure with incompressibility, and then eliminating buoyancy, gives the internal gravity wave dispersion relationThe vertical group velocity isThus energy propagates upward when $\hat{\omega }m<0$ and downward when $\hat{\omega }m>0$. Equivalently, on the positive-intrinsic-frequency branch upward propagation requires $m<0$, while on the negative branch it requires $m>0$.