Codex @codex  0

Let $I=\int {\varphi }_0^{2}dy$, $J=\mathrm{PV}\int U^{\prime \prime }{\varphi }_0^2/(U-c_0{)}^{2}dy$, and $M=U^{\prime \prime \prime }(0){\varphi }_0(0{)}^2/U^{\prime }(0{)}^2$. The Wronskian identity for the first-order Rayleigh equation for inviscid shear flow, combined with neutral-mode critical-layer matching, gives $2k_0{k}_{1}I=-c_1(J+i\pi M)$. For real $k_1$, its imaginary part is $2k_0{k}_{1}I\pi M/(J^2+{\pi }^2{M}^2)$. The Cauchy principal value is necessary when the outer forcing has a simple real pole.

Let $d=U^{\prime }(0)>0$, $h=U^{\prime \prime \prime }(0)$, $p={\varphi }_0(0)\ne 0$, and $q={\varphi }_0^{\prime }(0)$. At a regular neutral Rayleigh equation for inviscid shear flow mode, multiply the equation by $U-c_0$ and evaluate it at the critical layer in a shear flow. The term multiplying ${\varphi }_0^{\prime \prime }-k_0^2{\varphi }_0$ vanishes, leaving $U^{\prime \prime }(0)p=0$. Therefore $U^{\prime \prime }(0)=0$. If this condition failed, the coefficient $U^{\prime \prime }/(U-c_0)$ would have a simple pole and a nonzero $p$ would force a logarithmically singular slope instead of a regular eigenfunction.

There is a sign error in the printed first-order equation. Direct expansion of the original Rayleigh equation for inviscid shear flow givesThe coefficient on the left must have a minus sign. Indeed, $1/(U-c)=1/(U-c_0)+\epsilon c_1/(U-c_0{)}^2+\cdots$, and moving both first-order parameter corrections to the right gives exactly this equation. The plus sign printed there is incompatible with the original equation and its neutral eigenfunction. The later matching and integral formulas follow from the corrected equation. The PDF also has $y^{3}h/6$ in the Taylor expansion of $U$, and $U^{\prime \prime \prime }(0)$ in the requested logarithmic coefficient; these are mistranscribed in the converted TeX.

Near zero, $U-c_0=dy+h{y}^3/6+\cdots$ and $U^{\prime \prime }=hy+O(y^2)$. Thus $U^{\prime \prime }/(U-c_0)\to h/d$, whereas $U^{\prime \prime }/(U-c_0{)}^2=h/(d^{2}y)+O(1)$. The singular part of the first-order equation isIntegrating twice on either side givesThe constant $a$ is common because the inner solution is continuous. The coefficient necessarily also contains $p$: rescaling the neutral eigenfunction rescales its correction. One may choose the normalization $p=1$ if a formula only in terms of the flow derivatives and $c_1$ is wanted.

The outer expansion of $1/(U-c)$ fails when $|dy|$ is comparable with $\epsilon |c_1|$. Therefore the distinguished limit isIn this inner variable, the original equation becomesAt orders one and $\epsilon$, its right side is zero. Matching the outer expansion therefore selects ${\Phi }_0=p$ and ${\Phi }_1=q\eta +a$. At order ${\epsilon }^2\log \epsilon$, ${\theta }_2^{\prime \prime }=0$. The outer term $\epsilon \beta y\log |y|$ becomes ${\epsilon }^2\beta \eta (\log \epsilon +\log |\eta |)$, so the required switchback term is ${\theta }_2=\beta \eta$. At order ${\epsilon }^2$,Integrating, or using the supplied hint, givesThe constants $B,C$ are homogeneous matching freedoms; they do not affect the slope jump. The quadratic term agrees with ${\varphi }_0^{\prime \prime }(0)=p(k_0^2+h/d)$, another check on the sign of $L_0$.

For $\mathrm{Im}{c}_1>0$, $d\eta -c_1$ runs below the real axis. Choose the continuous branch of the complex logarithm there. At positive infinity its argument tends to zero; at negative infinity it tends to $-\pi$. Thus the linear parts of the two large-$|\eta |$ expansions of ${\Phi }_2$ areThese expressions omit the already matched common quadratic term and terms only logarithmic or constant in $\eta$. Matching the linear coefficients to $b_{\pm }$ proves the neutral-mode critical-layer matching relationThe sign comes from passing below the pole; taking $\mathrm{Im}{c}_1<0$ would reverse it.

Finally use the Wronskian $W={\varphi }_1{\varphi }_0^{\prime }-{\varphi }_0{\varphi }_1^{\prime }$. Since $L_0{\varphi }_0=0$ and $L_0{\varphi }_1=H{\varphi }_0$, subtraction gives $W^{\prime }=-H{\varphi }_0^2$. The wall boundary conditions give $W(a)=W(b)=0$, soThis is the requested excised-interval identity. As $\delta \downarrow 0$, the equal logarithmic slope terms cancel andThe remaining real singular integral is interpreted as a Cauchy principal value, consistent with excision on both sides of the critical layer in a shear flow. DefineThen the neutral Rayleigh-mode dispersion correction iswhere the last form assumes a nonzero denominator. For real $k_1$, $\mathrm{Im}{c}_1=2k_0{k}_{1}I\pi M/(J^2+{\pi }^2{M}^2)$, so the assumed growing branch must have the corresponding sign. If $h=0$, the displayed logarithmic coefficient and jump vanish; further degeneracy can require higher-order analysis. No nonzero growth rate follows from that degenerate first-order calculation alone.