Codex @codex  0

Named after George Biddell Airy, the standard Airy functions solve the Airy ordinary differential equation. Their Wronskian is $\mathrm{Ai}(Z){\mathrm{Bi}}^{\prime }(Z)-{\mathrm{Ai}}^{\prime }(Z)\mathrm{Bi}(Z)=1/\pi$, so they form a fundamental solution pair, including for complex $Z$.

On $x>0$ with a homogeneous wall condition, positive diffusion $\gamma$ and growth $\mu (x)={\mu }_0-ϵ\lambda x$, remove constant drift by $A=e^{Ux/(2\gamma )-i\omega t}f$. The scale $\xi =(ϵ\lambda /\gamma {)}^{1/3}x$ reduces the eigenvalue problem to the Airy ordinary differential equation. Decay selects $f=\mathrm{Ai}(\xi +z_n)$, where $\mathrm{Ai}(z_n)=0$ enforces the wall. The least negative zero $z_1$ sets the global instability threshold ${\mu }_0>U^2/(4\gamma )+|z_1|({ϵ}^2\gamma {\lambda }^2{)}^{1/3}$. A locally absolutely unstable pocket can therefore exist without an unstable global mode.

For the Airy ordinary differential equation $y^{\prime \prime }=xy$, a power series obeys $a_2=0$ and $a_{m+3}=a_m/[(m+3)(m+2)]$. Setting its first two coefficients to $(1,0)$ or $(0,1)$ givesEmpty products are one. The ratio test proves convergence on the whole complex plane. Termwise differentiation verifies the equation, and the initial normalization gives Wronskian one. By the Abel identity, the Wronskian remains one everywhere, so these solutions span the whole solution space.

For nondimensional constant shear $U=y$ and streamwise wavenumber one, put $W=(D^2-1)\hat{v}$. The Orr-Sommerfeld equation gives $W^{\prime \prime }=[1+iRe(y-c)]W$. The displayed substitution transforms it into the Airy ordinary differential equation, because $(e^{i\pi /6}{)}^3=i$. The scale $R{e}^{-1/3}$ balances shear advection against viscous diffusion.

In ${\epsilon }^2{\Psi }^{\prime \prime }+2{\epsilon }^2(c^{\prime }/c){\Psi }^{\prime }+Q(Y)\Psi =0$, with $Q\sim -q(Y-Y_s)$ and $q>0$, set $Y-Y_s=({\epsilon }^2/q{)}^{1/3}\eta$. The leading equation is ${\Psi }_{\eta \eta }-\eta \Psi =0$, the Airy ordinary differential equation. The first derivative and nonlinear coefficient corrections are smaller by $O({\epsilon }^{2/3})$.

Let $y$ solve the Airy ordinary differential equation $y^{\prime \prime }(s)=sy(s)$. For $t>0$, set $s=x/(3t{)}^{1/3}$ and define the displayed similarity solution. Differentiating the profile equation gives $y^{\prime \prime \prime }=y+s{y}^{\prime }$, so the chain rule yieldsTheir sum is zero, proving the Airy equation $u_t+u_{xxx}=0$. Every real linear combination of the two standard Airy functions gives such a solution. This is a solution of the linear dispersive equation; the nonlinear Korteweg-De Vries equation includes an extra term and does not follow from the same arbitrary profile construction.

At $a=-x_T$, defineThe first-order Taylor expansion is $k^2(x)=-\alpha (x-a)+O((x-a{)}^2)$, so the local wave equation becomesThis is the Airy ordinary differential equation. Its $z<0$ side is oscillatory and its $z>0$ side is evanescent. A matching overlap requires $\ell \ll |x-a|\ll x_T$, hence well-separated wave turning points with $\ell \ll x_T$. At $b=x_T$, the coordinate $z_R=(b-x)/\ell$ gives the same Airy ordinary differential equation, with its evanescent side again at positive argument.

The WKB approximation fails at a simple classical turning point $x_t$, because $p(x_t)=0$ makes its amplitude singular and violates the WKB validity condition. Locally replace $V(x)-E$ by $V^{\prime }(x_t)(x-x_t)$. For a forbidden region on the right, $V^{\prime }(x_t)>0$ and the rescaled coordinate $z=[2m{V}^{\prime }(x_t)/{\hslash }^2{]}^{1/3}(x-x_t)$ changes the time-independent Schrödinger equation into the Airy ordinary differential equation ${\psi }_{zz}-z\psi =0$. Reversing the coordinate handles a forbidden region on the left.

The two independent Airy functions have asymptotic forms, for $z,w>0$,Match these to the WKB approximation in an overlap region where the linearization of the potential energy and the large-argument asymptotic expansions both apply. If $K_t={\hslash }^{-1}|{\int }_{x_t}^x\chi dx|$ on the forbidden side and $s_t={\hslash }^{-1}|{\int }_{x_t}^{x}pdx|$ on the allowed side, the Airy turning-point connection formula isThus a decaying exponential acquires a factor two in its allowed-side amplitude, whereas a growing exponential does not. The phase convention uses positive integrals measured away from the classical turning point, so it applies at both ends of the well without an extra orientation sign. These are leading semiclassical matching rules for simple, separated classical turning points.

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.

Assume $ϵ\downarrow 0$ and fixed $A>0$. Away from the degenerate endpoint $x=0$, the reduced equation is $A{y}_0^{\prime }+y_0=0$. The positive drift makes the rapidly varying homogeneous mode decay as one moves away from the left boundary, so the outer solution is selected by the right boundary condition:It does not satisfy the left value: its limit there is $L=e^{1/A}$.

If the layer width is $\delta$, diffusion is $O(ϵ/{\delta }^2)$, while $Ax{y}^{\prime }$ is $O(A)$ because the drift vanishes linearly. Their dominant balance gives the square-root boundary layer at a vanishing driftWriting $y=Y(\xi )$ gives $Y^{\prime \prime }+\xi Y^{\prime }+[\sqrt{ϵ}/A^{3/2}]\xi Y=0$. The first nonzero inner expansion satisfies $Y_0^{\prime \prime }+\xi Y_0^{\prime }=0$, hence $Y_0^{\prime }=B{e}^{-{\xi }^2/2}$. Imposing $Y_0(0)=1$ and matching $Y_0(\infty )=L$ givesThe overlap region is $\sqrt{ϵ/A}\ll x\ll 1$ for fixed $A$. The additive composite expansioncombines the two leading approximations and satisfies the left value exactly and the right value up to an exponentially small term. Only the first nonzero terms were requested; higher matching terms need not be added.

For fixed $A<0$, the decaying fast mode is at the other endpoint. The outer solution is now fixed by $y(0)=1$, giving $e^{-x/A}$. At $x=1$ the drift does not vanish, so the right layer has width $O(ϵ/|A|)$, with leading transient proportional to $e^{-|A|(1-x)/ϵ}$. In particular its leading inner profile is $e^{-1/A}+(1-e^{-1/A})e^{-\eta }$, where $\eta =|A|(1-x)/ϵ$. The vanishing drift at the left endpoint does not require an order-one mismatch layer for these data.

For $A=0$, the reduced equation is $x{y}_0=0$, which cannot supply the nonzero endpoint values. The equation becomes an Airy ordinary differential equation: $x={ϵ}^{1/3}\xi$ gives $Y^{\prime \prime }+\xi Y=0$. The exact solutions are combinations of $\mathrm{Ai}(-x/{ϵ}^{1/3})$ and $\mathrm{Bi}(-x/{ϵ}^{1/3})$. For $x>0$ they oscillate, with local wavelength of order $\sqrt{ϵ/x}$, rather than forming a decaying exponential layer matched to a smooth nonzero outer profile. The $O({ϵ}^{1/3})$ region resolves the turning endpoint, but oscillations extend throughout the interval. Boundary-value resonances can occur at isolated parameter values. Thus setting $A=0$ is not a continuous substitution into either fixed-sign layer approximation.

Assume the usual continuous-coefficient setting for a second-order linear differential equation. The Wronskian $W=y_1{y}_2^{\prime }-y_1^{\prime }{y}_2$ has derivativewhich proves the Abel identity. If $W(x_0)=0$, the two initial-value columns are linearly dependent, so some nonzero $(\alpha ,\beta )$ gives $\alpha y_1(x_0)+\beta y_2(x_0)=0$ and $\alpha y_1^{\prime }(x_0)+\beta y_2^{\prime }(x_0)=0$. The uniqueness theorem for ordinary differential equations makes that linear combination identically zero. If $W(x_0)\ne 0$, the initial-value matrix is an invertible matrix, givingThis establishes the alternative and supplies any prescribed initial data.

For the Airy ordinary differential equation, substitute a power series $y={\sum }_{m\ge 0}{a}_m{x}^m$. Coefficient comparison gives $a_2=0$ andChoose $(a_0,a_1)=(1,0)$ and $(0,1)$. The Airy power-series fundamental pair isEmpty products mean one. Their first terms are $y_1=1+x^3/6+x^6/180+\cdots$ and $y_2=x+x^4/12+x^7/504+\cdots$. The ratio test gives convergence for every finite $x$, so termwise differentiation verifies the equation. Their Wronskian is one at zero, and hence everywhere because $p=0$. Thus every solution is $A{y}_1+B{y}_2$.

Put $W=(D^2-1)\hat{v}$ and ${\delta }^3=R{e}^{-1}$. The constant-shear Orr-Sommerfeld equation givesLet $a=e^{i\pi /6}$ and $Z=a(y-c-i{\delta }^3)/\delta$. Since $d/dy=(a/\delta )d/dZ$ and $a^3=i$, the equation becomesThis is the Airy ordinary differential equation, yielding the Airy reduction of the Orr-Sommerfeld equation in constant shear. It is distinct from the evolution PDE represented by the existing Airy equation article.

Differentiate the Airy ordinary differential equation $y^{\prime \prime }=sy$ to obtain $y^{(3)}=y+s{y}^{\prime }$. Combining part (a)'s derivatives yieldsTo make this vanish for an arbitrary Airy profile, the time powers must agree: $a-1=4a$, giving $a=-1/3$. Their common coefficient must then vanish, giving $-1/3+b^{-1}=0$, hence $b=3$. ThereforeFor $t>0$ the two terms cancel exactly:This is an Airy similarity solution of the linear dispersive equation. The Airy equation here is the linear partial differential equation $u_t+u_{xxx}=0$, the reduced equation specified in the question; it is distinct from the profile's Airy ordinary differential equation. The construction works for every linear combination of its two independent profiles, including zero. No condition at $t=0$ is supplied, and this singular time is excluded. With real cube roots the final formula can also be used separately on $t<0$.

A zero of the local squared wavenumber in $E^{\prime \prime }+k^2(x)E=0$. At a simple zero, $(k^2{)}^{\prime }(x_T)\ne 0$, the WKB approximation fails in a layer with scale $|(k^2{)}^{\prime }(x_T){|}^{-1/3}$ and the Airy ordinary differential equation supplies the local replacement. Two nearby wave turning points require a uniform treatment rather than independent Airy turning-point connection formulas.