Codex @codex  0

Assume again positive diffusion $\gamma >0$. Remove the drift with $A(x,t)=e^{Ux/(2\gamma )-i\omega t}f(x)$. Differentiating explicitly cancels the first derivative and givesBalance the second derivative with the linear confinement by taking $x=O({ϵ}^{-1/3})$, so the requested exponent is $\sigma =1/3$. More precisely, putThe eigenvalue problem reduces to $f_{\xi \xi }=(\xi -c)f$. Its two solutions are Airy functions. The growing $\mathrm{Bi}(\xi -c)$ is excluded by the condition at infinity: its growth proportional to $e^{2{\xi }^{3/2}/3}$ dominates any linear-in-$x$ drift factor. The decaying $\mathrm{Ai}(\xi -c)$ remains admissible even after multiplication by $e^{Ux/(2\gamma )}$. The wall imposes $\mathrm{Ai}(-c)=0$, so $-c=z_n$, giving the Airy global modes of a linearly confined Ginzburg-Landau equation:The least negative zero, $z_1\simeq -2.338107$, gives the largest temporal growth rate. HenceAt equality the leading global eigenmode is neutral, and below it all these eigenmodes decay. In contrast, freezing the coefficients locally gives absolute hydrodynamic instability wherever ${\mu }_0-ϵ\lambda x>U^2/(4\gamma )$. Such a pocket exists iff ${\mu }_0>U^2/(4\gamma )$ and extends to $x_a=[{\mu }_0-U^2/(4\gamma )]/(ϵ\lambda )$. Thus global instability implies a locally absolutely unstable pocket, while the converse fails: a pocket shorter than $|z_1|(\gamma /(ϵ\lambda ){)}^{1/3}$ does not overcome the wall and confinement losses. The additional threshold is order ${ϵ}^{2/3}$, precisely the loss produced by the Airy function spatial scale.

Away from the turning points, the Liouville–Green approximation gives in the allowed regionThe bound-state condition discards the growing exponential in each forbidden region, leavingThese approximations require variation slow compared with the local wavelength, and must be replaced near $a,b$ by Airy approximations.

Near $a$, set $z=-{\lambda }^{2/3}{q}^{\prime }(a{)}^{1/3}(x-a)$. The local differential equation becomes ${\psi }_{zz}-z\psi =0$. Decay on the left selects Airy function $\mathrm{Ai}(z)$ rather than its growing companion. Comparing the supplied asymptotics across the turning point yieldsAt $b$, use $z={\lambda }^{2/3}[-q^{\prime }(b){]}^{1/3}(x-b)$; decay on the right similarly yields $2C_R{q}^{-1/4}\cos (S_b(x)-\pi /4)$, where $S_b=\lambda {\int }_x^b\sqrt{q}ds$. Let $J=S_a+S_b=\lambda {\int }_a^b\sqrt{q}ds$. The left expression has equal cosine and sine coefficients in the variable $S_a$; the right expression has those coefficients proportional to $\cos (J-\pi /4)$ and $\sin (J-\pi /4)$. They can describe the same nonzero wave only when $J-\pi /4=\pi /4+n\pi$. Therefore the WKB quantization condition isThe two simple turning points contribute the two quarter-phase shifts, giving the half-integer correction.

The lowest Liouville–Green approximation comes from $f=e^S$, for which $S^{\prime \prime }+(S^{\prime }{)}^2=Q$. First take $S^{\prime }=\pm \sqrt{Q}$; the next transport term is $-Q^{\prime }/(4Q)$. Integrating givesMultiplicative constants and the lower limit are arbitrary. The negative sign selects the recessive approximation in the usual WKB approximation regime, where the exponent dominates the amplitude and the derivative corrections are small. Positivity of $Q$ by itself does not prove that $f_{-}$ tends to zero or that the approximation is valid: for example $Q=e^{-x}$ makes this expression grow like $e^{x/4}$ with a bounded integral in its exponent. For $Q=x$, all the required large-$x$ conditions hold and $f_{-}=x^{-1/4}{e}^{-2x^{3/2}/3}\to 0$.

Put $a=f_{-}$. Direct differentiation givesSubstituting $\mathrm{Ai}=wa$ into the Airy ordinary differential equation ${\mathrm{Ai}}^{\prime \prime }=x\mathrm{Ai}$ and dividing by $a$ gives $w^{\prime \prime }+2(a^{\prime }/a)w^{\prime }+5w/(16x^2)=0$. Multiplication by $x^2$ provesThe recessive branch has a constant leading amplitude $c$. To find its first correction write $w=c(1+b{x}^{-\delta }+\cdots )$. Balancing $-2x^{5/2}{w}^{\prime }$ against the constant term $5c/16$ forces $3/2-\delta =0$, so $\delta =3/2$. The constant coefficients then give $3b+5/16=0$. The remaining terms are of order $x^{-3/2}$ and are cancelled by a next term of order $x^{-3}$. ThusThe multiplicative constant depends on the chosen normalization of the Airy function and of $f_{-}$. The calculation selects the recessive solution; a growing Airy solution would not have this constant-amplitude expansion.

For a stationary component, the stationary internal-wave WKB solution has local vertical wavenumber and vertical-velocity amplitudewith the upward branch $m>0$. The WKB approximation requires $|m^{\prime }|\ll m^2$ and sufficient smoothness, with $|m^{\prime \prime }|\ll |m{|}^3$ a useful additional criterion. The background changes little over a vertical wavelength; this excludes a turning level of a stationary internal wave.

For increasing $N=N_0(1+z/H)$, the vertical wavenumber increases. The crest slope $-k/m$ tends toward zero, and successive crests get closer vertically. The laboratory energy ray obeys $dx/dz=k/m$, so it becomes more nearly vertical as it rises. There is no upper turning level in this ideal increasing-stratification profile.

For decreasing $N=N_0(1-z/H)$, the vertical wavenumber decreases to zero atThis assumes propagation at the bottom, $k{U}_0<N_0$. The crests become more nearly vertical and the incident energy ray becomes horizontal at $z_t$. Above it the wave is evanescent, and the incident wave is reflected downward. The divergent WKB amplitude at $m=0$ is an approximation failure, resolved by an Airy function transition, rather than an infinite physical wave amplitude. The decreasing formula for $N$ is only a positive-frequency background below $z=H$; the relevant turning height lies below that level.

Take $\alpha >0$; a real disturbance with negative wavenumber is recovered by complex conjugation. To obtain the pressure equation for a shear-flow normal mode, apply $i\alpha$ to the streamwise momentum equation and differentiate the normal momentum equation. The mass conservation equation cancels the viscous divergence and the terms proportional to $U-c$, leaving $2i\alpha U^{\prime }v={\alpha }^{2}p-p^{\prime \prime }$. Thus the fluid pressure equation and wall value areThe wall value follows from streamwise momentum with no-slip boundary condition. Normal momentum also gives $p^{\prime }(0)=v^{\prime \prime }(0)/R$.

In the lower viscous boundary layer, mass conservation gives $v=O(\alpha \delta )$ if $u=O(1)$. Since $U\sim y$, convection and shear have size $\alpha \delta$ when $c=O(\delta )$; matching these with transverse viscous diffusion gives $R^{-1}{\delta }^{-2}\sim \alpha \delta$. The pressure gradient balance then gives $p=O(\delta )$. Thus the lower scalings areThe neglected streamwise viscous term is smaller by $(\alpha \delta {)}^2$. In these variables the leading equations areDifferentiating the second and using the first yields the Airy reduction $\boxed{u_{YYY}=i(Y-C)u_Y}$.

Put $q=e^{i\pi /6}$, so $q^3=i$. The bounded matching branch of the Airy function gives $u_Y=A\mathrm{Ai}(q(Y-C))$; the other independent Airy solution grows along the ray and is excluded. Integrating with the two wall conditions giveswhile the wall momentum equation gives the constant fluid pressureDefine the convergent Airy displacement integralsThe Airy function equation and its decaying derivative imply $I_1-C{I}_0=iq{\mathrm{Ai}}^{\prime }(-qC)$. Consequently, as $Y\to \infty$,Both differentiated lower momentum and the constant of integration are satisfied; retaining that constant is essential for the final matching.

For the middle region introduce a perturbation stream function $\Phi$ with $u={\Phi }^{\prime }$ and $v=-i\alpha \Phi$. The leading inviscid flow streamwise equation isFor a shear flow, this is the integrated form of the Rayleigh equation for inviscid shear flow. The leading matched solution for fixed $y=O(1)$ is $\Phi \sim BU$, so $u\sim B{U}^{\prime }$ and $v\sim -i\alpha BU$. In the overlap with the lower viscous boundary layer, the first constant correction iswhich agrees with the integrated lower solution. For a regular smooth profile, the normal momentum equation gives $p^{\prime }\sim -{\alpha }^2(U-c)\Phi$. Thus fluid pressure variation across a fixed-width middle region is $O({\alpha }^{2}B)$; in the final distinguished scaling this is smaller than its $O(\delta )$ value, so $p\sim \delta P_0$ throughout that region. Since $U\to 1$, the middle velocity tends to $\boxed{v\sim -i\alpha B}$, which does not decay and requires another region.

The fluid pressure equation becomes $p^{\prime \prime }-{\alpha }^{2}p=0$ when $U\sim 1$. Therefore the upper scale is $\Delta =1/\alpha$, with $z=\alpha y$. The inviscid decaying upper solution isMatching velocity as $z\downarrow 0$ gives $P_{\infty }\sim \alpha B$; matching fluid pressure to the middle region gives $\delta P_0\sim \alpha B$. This is three-layer long-wave shear-flow matching, expressed as a matched asymptotic expansion.

Under $\alpha =k{R}^{-1/4}$ with fixed positive $k$, $\delta =k^{-1/3}{R}^{-1/4}$, and $\alpha /\delta =k^{4/3}$. Eliminating the nonzero amplitude $A$ gives the Airy wall-layer matching determinant and dispersion relationEquivalently, with $\mathcal{T}(C)={\int }_{-qC}^{\infty e^{i\pi /6}}\mathrm{Ai}(\zeta )d\zeta =q{I}_0(C)$,The determinant form does not divide by a possibly vanishing $I_0$. The wavespeed is $c=\delta C$; temporal growth has sign $\mathrm{Im}C$ for positive $\alpha$.

There is a genuine conflict in the printed assumptions: fixed-$k$ $\alpha =k{R}^{-1/4}$ cannot satisfy $R^{-1/7}\ll \alpha$ as $R\to \infty$. For order-one lower matching amplitudes, constant leading middle fluid pressure requires ${\alpha }^2\ll \delta$, equivalently $\alpha \ll R^{-1/7}$; this is the reverse inequality and is satisfied by the final distinguished scaling. Thinness also requires $\alpha R\gg 1$. The dispersion relation above is therefore derived under the final distinguished limit, for which the complete consistent window is $R^{-1}\ll \alpha \ll R^{-1/7}$. Both printed limits cannot be imposed simultaneously.

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.

Use the geometrical optics conventionBecause ${\partial }_t=\epsilon {\partial }_T$ and ${\partial }_x=\epsilon {\partial }_X$, differentiating $e^{i\theta /\epsilon }$ gives an order-one physical angular frequency and order-one physical wavenumbers. Absorb a constant leading amplitude phase into $\theta$ and take $A_0$ real; for a complex amplitude the corresponding quadratic density is $|A_0{|}^2$.

Substitution into the wave equation gives the eikonal equation for a variable-speed wave at leading order:Choose the positive-frequency branch $\omega =\Omega =cK$, $K=\sqrt{k^2+{\ell }^2}$; the negative branch can be treated in exactly the same way. The imaginary first-order terms giveMultiplying by $A_0$ puts this wave-action transport for a variable-speed wave in conservation form:Its density is $\omega A_0^2$ and its flux is that density times the group velocity. Independently, equality of mixed wave phase derivatives givesThis phase compatibility in geometrical optics is essential when converting the differentiated dispersion relation into equations along rays.

For completeness let $\omega =\Omega (\mathbf{k};\mathbf{X},T)$ be a general local dispersion relation. Define ray tracing by $d\mathbf{X}/dT={\mathbf{c}}_g={\nabla }_{\mathbf{k}}\Omega$. Using the chain rule in ${\nabla }_{\mathbf{X}}\omega$ and the symmetry ${\partial }_{X_i}{k}_j={\partial }_{X_j}{k}_i$ givesLikewise ${\omega }_T={\Omega }_T+{\mathbf{c}}_g\cdot {\mathbf{k}}_T$ cancels its advective term because ${\mathbf{k}}_T=-\nabla \omega$. Hence the Hamiltonian ray equations for a local dispersion relation areThe derivatives on the right hold the wavevector fixed. For this wave equation they reduce to

When $c=c(Y)$, the ray invariants are $k$ and $\omega$, while $d\ell /dT=-K{c}_Y<0$. For a monochromatic component with these fixed invariants write $\theta =kX-\omega T+S(Y)$, soUsing the stated independence of $A_0$ from $X,T$, the wave-action transport for a variable-speed wave becomes $(c^2\ell A_0^2{)}_Y=0$. ThereforeThis is the stationary wave amplitude in a stratified wavespeed; $\mu$ is the constant for that component or ray family.

At the proposed simple turning point of a variable-speed wave, take $c_s=c(Y_s)=\omega /k$ and $c_s^{\prime }=c_Y(Y_s)>0$. On the chosen positive branch this entails $k>0$. DefineAs the ray approaches from below, $\ell \sim [q(Y_s-Y){]}^{1/2}$ andThe amplitude diverges and the vertical group velocity $c^2\ell /\omega$ vanishes. The WKB approximation also fails its slow-variation condition $\epsilon |{\ell }_Y|/{\ell }^2\ll 1$ when $Y_s-Y=O({\epsilon }^{2/3})$. The turning-point enhancement of wave amplitude shows that this divergence is a failure of the outer ray approximation, rather than an actual infinite solution.

To derive the inner wavefield directly, separate the conserved horizontal wave phase and frequency in the original differential equation:Dividing by $c^2$ gives ${\epsilon }^2{\Psi }_{YY}+2{\epsilon }^2(c_Y/c){\Psi }_Y+Q(Y)\Psi =0$. Set $Y-Y_s=\delta \eta$. The second derivative is of size ${\epsilon }^2/{\delta }^2$, whereas $Q\sim -q\delta \eta$, so the Airy scaling at a variable-speed wave turning point isThe first-derivative term is smaller by $O(\delta )$, and higher coefficients in the expansion of $Q$ give the same relative order. Thus the general leading inner profile is a linear combination of the two Airy functions $\mathrm{Ai}(\eta )$ and $\mathrm{Bi}(\eta )$. For decaying Airy continuation and ray reflection, if the physical continuation is bounded and decays into $Y>Y_s$, its growing $\mathrm{Bi}$ coefficient is zero. In that usual evanescent continuation the wavefield has the formWithout a condition in the evanescent region the corresponding $\mathrm{Bi}$ term remains admissible locally; the given existence of a wavefield below the turning point alone does not select its coefficient.

Finally, the outer amplitude at distance $\delta$ is $O({\delta }^{-1/4})=O({\epsilon }^{-1/6})$ for $\mu =O(1)$ and fixed nonzero $q$. The oscillatory Airy turning-point connection formula has exactly the same $(-\eta {)}^{-1/4}$ envelope on the allowed side. Therefore the inner wavefield has magnitude $O({\epsilon }^{-1/6})$. The decaying Airy continuation matches an incident and reflected pair of rays, so a single upward $\ell >0$ branch must be supplemented by its reflected branch near the turning point; no formal coefficient matching is needed to identify this magnitude.

Write $W(s)=a_1\mathrm{Ai}(Z(s))+a_2\mathrm{Bi}(Z(s))$, using the independent Airy functions. To solve $(D^2-1)\hat{v}=W$, the variation of parameters formula, equivalently a direct convolution, gives the velocity reconstruction from the Orr-Sommerfeld vorticity:The reference point $y_0$ is arbitrary. Twice differentiating the integral produces $W(y)$ because $\sinh 0=0$ and its derivative at zero is one, verifying the formula. The four constants are independent: applying $D^2-1$ to an identically zero combination first forces $a_1=a_2=0$, then the complementary solutions force $C_{+}=C_{-}=0$. The boundary conditions would be $\hat{v}(y_j)=D\hat{v}(y_j)=0$, $j=1,2$. Their compatibility selects the eigenvalues; they need not be imposed here.

The hard wall imposes $\psi (0)=0$, while the increasing linear potential requires decay as $x\to \infty$. SetThe stationary Schrödinger equation becomes $d^2\psi /d{\xi }^2-\xi \psi =0$. The growing Airy function solution is excluded, leavingThis is the quantum bouncer. For the first four states, $z_n\simeq 2.3381,4.0879,5.5206,6.7867$. Their interior node counts are respectively zero, one, two and three. Each wavefunction starts at zero at the wall, oscillates in the classically allowed interval, and crosses smoothly into an exponentially decaying tail beyond the turning point $x_n=E_n/\beta =z_n\ell$. The turning point is not an extra hard wall: the wavefunction generally remains nonzero there and beyond it.

A high-Reynolds number long-wave mode near a wall has a viscous layer of thickness $\delta =(\alpha R{)}^{-1/3}$, an order-one region following the base profile, and a pressure-decay region of thickness $1/\alpha$. The lower Airy function solution fixes velocity displacement and fluid pressure. The middle Rayleigh equation for inviscid shear flow transports these to the outer constant-velocity region. Matching its decaying pressure produces an Airy wall-layer matching determinant. For order-one amplitudes, constant leading middle pressure requires ${\alpha }^2\ll \delta$; the distinguished balance $\alpha \sim \delta$ gives $\alpha \sim R^{-1/4}$.

For constant buoyancy frequency $N$ and $U=U_0+Sz$ with $S>0$, a stationary component of horizontal wavenumber $k$ turns at $z_t=(N/k-U_0)/S$ if it propagates at the source. Its vertical wavenumber vanishes there; WKB approximation fails and an Airy function transition connects reflection to an upper evanescent wave.

A vertically propagating stationary topographic internal gravity wave meets a turning level when its squared vertical wavenumber changes from positive to negative. The stationary internal-wave WKB solution then fails locally. A simple zero is resolved by an Airy function transition connecting incident and reflected waves below to a decaying evanescent wave above. For negligible wind curvature the condition is $|kU(z_t)|=N(z_t)$, distinct from the zero-intrinsic-frequency critical layer in a shear flow.

A simple ray turning point has outer amplitude proportional to $|Y-Y_s{|}^{-1/4}$. At the Airy scaling at a variable-speed wave turning point, $|Y-Y_s|=O({\epsilon }^{2/3})$, this gives an inner wavefield of size $O({\epsilon }^{-1/6})$ when the incident action normalization is order one. The enhancement is finite for nonzero $\epsilon$ and is resolved by the Airy function.