Codex @codex  0

The WKB approximation treats the vertical wavenumber as locally constant over one vertical wavelength. Write $w=\Re \{A(z)e^{ikx+i{\int }^{z}m(s)ds}\}$, with $|m^{\prime }|\ll m^2$ and $|A^{\prime }|\ll |mA|$. For the stationary internal-wave WKB solution, $A\propto m^{-1/2}$ away from singular or turning levels. A small terrain slope and small perturbations are also needed for Linearized Boussinesq equations.

Let $p$ be the perturbation fluid pressure divided by reference mass density, and $b$ the buoyancy perturbation. The steady Linearized Boussinesq equations and incompressibility condition areFor an $e^{ikx}$ component, elimination gives $u=-w^{\prime }/(ik)$, $b=-N^{2}w/(ikU)$ and $p=(U{w}^{\prime }-U^{\prime }w)/(ik)$. Substitution into vertical momentum gives the exact stationary Taylor–Goldstein equationThe leading WKB approximation therefore identifiesThe kinematic boundary condition fixes the horizontal wavenumber to $k=k_T$ for the sinusoidal forcing.

For explicit heights assume a constant positive buoyancy frequency $N$, $0<k{U}_0<N$, and $|S|/N\ll 1$. The latter is a large gradient Richardson number. The linear shear flow has $U^{\prime \prime }=0$, so the propagation condition is $|U|<N/k$.

For $S<0$, the first obstruction is the critical level of an internal gravity wave atAs $z\uparrow z_c$, $m\sim N/U\to \infty$. The critical-height approach of a stationary wave in linear shear has upward group velocity $c_{gz}\sim k{U}^2/N$, so the travel time to $z_c$ diverges: the height is a supremum, approached in the ideal ray model. Here $|m^{\prime }|/m^2\to |S|/N$, which can remain small; it would be incorrect to say that the short-wavelength WKB approximation necessarily fails solely because $m$ diverges. Nevertheless $w\propto U^{1/2}$ and displacement $w/(ikU)\propto U^{-1/2}$, so a finite disturbance eventually violates linearity. Dissipation and overturning also matter near this singular level.

For $S>0$, the turning height of a stationary wave in linear shear isAt $z_t$, $m=0$ and the WKB approximation fails. An Airy function matching region describes reflection and an upper evanescent wave. Thus $z_t$ is the highest height of vertical propagation, rather than a sharp cut-off of every disturbance.

The problem initially allows $N(z)$ to vary. Without the extra constant-$N$ assumption, replace the explicit turning height by the first root of $N(z)=k|U_0+Sz|$ encountered on a propagating branch. A declining $N(z)$ can cause a turning level before the negative-shear critical level. The two closed-form heights require the stated constant-stratification assumption.

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.

Let the free surface be at $r=R_1$ and take the Newtonian gravitational potential to vanish at infinity. Mass conservation fixes $R_1=(3M_1/(4\pi {\rho }_1){)}^{1/3}$. The enclosed mass is $m(r)=4\pi {\rho }_1{r}^3/3$, and hydrostatic equilibrium requires $d{p}_0/dr=-{\rho }_{1}Gm(r)/r^2$. With vacuum outside, $p_0(R_1)=0$, givingThe pressure is zero outside. A prescribed constant external pressure would simply add that constant to $p_0$.

The bulk Eulerian and Lagrangian fluid perturbations obey $\delta \rho =0$ because the mass density is uniform and the Lagrangian displacement has zero divergence. The linearized Euler momentum equation in the interior isTaking the curl gives ${\omega }^2\nabla \times \mathbf{\xi }=0$. Thus every nonzero-frequency normal mode has an irrotational flow displacement; the simply connected space formed by the interior admits $\mathbf{\xi }=\nabla U$ and incompressibility gives the Laplace equation ${\nabla }^{2}U=0$. For a regular center, one spherical harmonic component is $U=A_l{r}^l{Y}_l^m$, with $|m|\le l$.

It is essential to retain the surface gravity perturbation of a uniform-density sphere. Although the bulk mass density perturbation is zero, the displaced density discontinuity givesThe Poisson equation implies that $\delta \Phi$ is continuous and its outward radial derivative has the jumpUsing regularity at the center and decay at infinity, writeThe derivative jump is $-(2l+1)C_l{R}_1^{l-1}{Y}_l^m$, so $C_l=-4\pi G{\rho }_{1}l{A}_l/(2l+1)$.

Integrating the bulk Euler momentum equation gives $\delta p/{\rho }_1={\omega }^{2}U-\delta \Phi$ for $l\ge 1$. At the free surface the Lagrangian pressure perturbation vanishes:Combining these relations yields the Kelvin stellar mode frequencyThe $2l+1$ choices of $m$ are degenerate because of spherical symmetry. For $l=1$, $r{Y}_1^m$ is a linear combination of coordinates in a Cartesian coordinate system, so $\nabla U$ is a constant displacement: the zero frequency is rigid translation of the whole isolated body. Such a translation has no restoring force. The constant $l=0$ potential generates no displacement; a radial breathing motion is excluded by incompressibility and a regular center.

The printed irrotational assertion needs the nonzero-frequency qualification, or restriction to the potential incompressible stellar surface modes. It is false for all neutral displacements: $\mathbf{\xi }=\mathbf{\Omega }\times \mathbf{r}$ has zero divergence, zero normal displacement at the surface, and $\delta p=\delta \Phi =0$, so it is a zero-frequency displacement, but $\nabla \times \mathbf{\xi }=2\mathbf{\Omega }\ne 0$. This does not change the Kelvin stellar mode spectrum just derived.

The separation of the two bodies is $a$, as is fixed by the supplied Kepler third law ${\omega }_o^2=G(M_1+M_2)/a^3$. The phrase “orbit of radius $a$” cannot denote the first body's distance from the center of mass in that formula.

Expand the companion's Newtonian potential of a point mass about the first body's center. The constant term produces no force; the linear term produces the uniform acceleration of the center and disappears in the translating frame. The leading differential acceleration therefore comes from the quadrupolar tidal forcingwhere $\mathbf{n}(t)$ points toward the companion. It is a degree-two solid spherical harmonic. For a circular nonrotating-frame orbit, its time-dependent components have $m=\pm 2$ and frequency $2{\omega }_o$, because the squared direction cosines repeat twice per orbit. There is also a static $m=0$ component. The real part of one complex harmonic represents the pair of time-dependent components; its normalization and orientation are absorbed into the order-one coefficient $\psi$. This gives $\Psi =\psi (G{M}_2/a^3)r^2{Y}_2^m$ and $\omega =2{\omega }_o$. The expansion is the leading term for $R_1/a$ small; spin is being neglected.

For a forced potential flow with the same temporal frequency, the Euler momentum equation becomesUse $U=A{r}^2{Y}_2^m$. The surface gravity perturbation of a uniform-density sphere is unchanged, as is the vanishing Lagrangian pressure perturbation boundary condition. Integrating the equation and applying that condition givesSince ${\xi }_r(R_1)=2A{R}_1{Y}_2^m$, the real surface displacement isIndeed $ϵ{\omega }_2^2=2G{M}_2/a^3$, so the normalization is consistent. This is a particular forced response away from the resonance; arbitrary free normal modes can be added.

Put $\mu =M_2/M_1$ and $q=(R_1/a{)}^3$. The tidal resonance of a stellar oscillation condition is $4{\omega }_o^2={\omega }_2^2$, henceThe resonance gives $0\le ϵ<1/2$ for finite nonnegative $\mu$, with actual forcing requiring $M_2>0$. At exact resonance the undamped linearization has no bounded response at the forcing frequency: the time-domain particular solution grows proportionally to $t$ times a sinusoid. Thus the frequency-response formula has a pole, rather than describing a finite resonant displacement. Away from that pole, linearization requires the displayed surface displacement to be small compared with $R_1$; the leading quadrupolar tidal forcing approximation separately requires $R_1/a$ small.

Evaluating the function at the endpoint gives the displacementThereforeThe force-induced displacement $h$ is a linear map of $F$.

Taking the derivative of the force-dependent energy givesThus the derivative of the additional minimized internal energy with respect to the applied force equals the resulting endpoint displacement. This is the endpoint force derivative of clamped-free beam energy.