Codex @codex  0

Rescale variables so that the terms responsible for the leading behaviour have comparable size. The resulting distinguished limit exposes a governing equation that may differ from the naive zero-parameter equation. Examples are dominant balance for algebraic roots and the derivative balance in a square-root boundary layer at a vanishing drift.

Here the drift vanishes at the left endpoint. The outer equation is $x{y}_0^{\prime }+y_0=0$, giving $y_0=1/x$ after applying the right boundary value. The ordinary $x=O(ϵ)$ layer does not balance diffusion and drift: their balance selects $x=O(\sqrt{ϵ})$ instead.

For the square-root boundary layer at a vanishing drift, put $X=x/\sqrt{ϵ}$. The inner equation isMatching $1/x$ requires a leading inner amplitude of order ${ϵ}^{-1/2}$, so an inner expansion with an everywhere order-one leading term would be incorrect. This can be seen directly by integrating the original equation once:The right boundary value determinesIn the inner variable the exact representation becomesThe first term enforces $Y(0)=1$, while the second is zero there and matches $C/(\sqrt{ϵ}X)=C/x$ for large $X$. One would expand this integral in the overlap region, determine successive outer terms from the right boundary, and subtract their common expansion. The endpoint layer has width $\sqrt{ϵ}$ and an $O({ϵ}^{-1/2})$ interior peak, while its value at the endpoint remains one. This distinguishes it from the ordinary exponential layers in the two preceding cases.

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.

For the matched asymptotic expansion, first determine the outer expansion. Setting $\epsilon =0$ away from $x=0$ gives $y_0^{\prime }=y_0$, and the right boundary condition sets $y_0=e^x$. At first order,soThere is no outer $\sqrt{\epsilon }$ term: its homogeneous coefficient would be forced to zero by the right boundary condition. The logarithmic overlap forces a switchback term in the inner expansion.

Near zero, balance $\epsilon y^{\prime \prime }$ against $2x{y}^{\prime }$ to get the square-root boundary layer at a vanishing drift:For fixed $z$ expand $Y=Y_0+s{Y}_1+s^2{Y}_2+\cdots$, allowing logarithmic dependence on $s$ in $Y_2$. With ${\mathcal{L}}_0={\partial }_z^2+2z{\partial }_z$, the leading equation is ${\mathcal{L}}_0{Y}_0=0$. Its bounded solutions are a constant plus an error function. The left boundary condition and the overlap value one fixThe following two solution sections complete the half-order and first-order matching; their headings retain the scaffold's organization of the supplied hints.