Codex @codex  0

Fluid-loaded plate. Use complex amplitudes with time dependence $e^{-i\omega t}$ and take the velocity potential in the fluid above the elastic plate to be $\Phi$. Let ${\rho }_0$ and $c_0$ be the undisturbed mass density and speed of sound. The linearized Euler equations give $p=i{\rho }_0\omega \Phi$; the kinematic boundary condition is ${\Phi }_y(x,0)=-i\omega \eta (x)$. We model loading from one fluid half-space, with no fluid loading on the other face. Loading by identical fluids on both faces instead doubles the fluid term below.

Define the Fourier transform and its inverse byThe Sommerfeld radiation condition selects $\hat{\Phi }(k,y)=\hat{\Phi }(k,0)e^{-\gamma y}$, with ${\gamma }^2=k^2-k_0^2$ and $k_0=\omega /c_0$. Here and in the diffraction calculation below, the outgoing square root is specified by first taking $\mathrm{Im}\omega >0$. Draw a vertical cut upwards from $+k_0$ and another downwards from $-k_0$, and putThus ${\gamma }_u$ is analytic in the upper half-plane and ${\gamma }_l$ in the lower half-plane. In the real-frequency limit,The Fourier transform contour $C$ follows the real axis with the limiting-absorption pole prescriptions. In particular, an outgoing propagating component has positive vertical wavenumber.

The kinematic boundary condition gives $\hat{\Phi }(0)=i\omega \hat{\eta }/\gamma$, so $\hat{p}(0)=-{\rho }_0{\omega }^2\hat{\eta }/\gamma$. Eliminating the fluid pressure from the bending balance gives the line-forced fluid-loaded bending plate kernelFor harmonic line force amplitude $F_0$ at the origin, the plate displacement and fluid velocity potential are thereforeAn arbitrary localized loading replaces $F_0$ by its Fourier transform. This representation both calculates the response and identifies its different wave components through poles and branch points.

For an observer $x=r\cos \theta$, $y=r\sin \theta$, with $0<\theta <\pi$, the method of steepest descent uses the saddle point $k_s=k_0\cos \theta$. At fixed angle away from grazing, the acoustic far field isThe $r^{-1/2}$ spreading is that of cylindrical acoustic waves. The phase velocity parallel to the plate is $\omega /|k|$: components with $|k|<k_0$ are supersonic in that sense and radiate, whereas those with $|k|>k_0$ are subsonic and exponentially decay normally to the plate. The latter have reactive fluid loading and zero mean normal acoustic energy flux.

There is also a real guided root $\kappa >k_0$. Indeed, $K(k_0+)=-\infty$, $K(+\infty )=+\infty$, andThe residue theorem then gives an outgoing displacement contribution $i{F}_0{e}^{i\kappa |x|}/K^{\prime }(\kappa )$. It transports energy along the plate while its fluid field decays as $e^{-\gamma (\kappa )y}$. Other zeros on continued sheets can describe leaky or resonant waves; branch point contributions describe near and grazing fields. A fixed-angle acoustic far field expansion is not uniform near the plate, where these components can matter.

For a propagating spectral component, write $\ell =(k_0^2-k^2{)}^{1/2}$. Its upward normal velocity is $i\ell \hat{\Phi }$ and its pressure is $i{\rho }_0\omega \hat{\Phi }$. Applying Parseval's identity to the time-averaged acoustic energy flux gives the upward power per unit length of the forcing line:Equivalently, the integrand is $|\hat{\eta }{|}^2/\ell$ with prefactor ${\rho }_0{\omega }^3/(4\pi )$. Integrating the radial acoustic intensity of the far field over a semicircle gives the same expression. The input power $\omega \mathrm{Im}(F_0^{\ast }\eta (0))/2$ also feeds guided waves along the plate, so it must not be identified with $P_y$ alone.

Rigid half-plane diffraction. Suppress the time factor and assume $0<{\theta }_0<\pi$. Write $a=k_0\cos {\theta }_0$, $b=k_0\sin {\theta }_0$, and denote the incident spatial field by $I=e^{-iax-iby}$. The Helmholtz equation, Sommerfeld radiation condition, and square-root branches are as above. Across the open positive axis the scattered field and its normal derivative are continuous. Across the rigid negative axis its normal derivative is also continuous, while its potential can jump. Let $J(x)$ be this jump, supported on $x<0$. The transformed scattered field isThe jump transform is analytic in the upper half-plane. The absence of an even scattered component follows from derivative continuity and the outgoing condition. Its normal derivative on either face is $-\gamma {\hat{J}}_u/2$.

On the negative axis this derivative must cancel $-ib{e}^{-iax}$. Since the Half-range Fourier transform of $H(-x)e^{-iax}$ is $i/(k+a)$, the Wiener-Hopf equation has the formwhere $G_l$ is the transform of an unknown function supported on $x>0$ and analytic in the lower half-plane. The incident pole $k=-a$ lies below the Fourier contour. Divide by ${\gamma }_l$ and perform pole subtraction in a Wiener-Hopf equation:The first term is upper analytic and the second lower analytic, since the latter's apparent pole cancels. Equality of the two continued sides gives an entire function. The finite-energy edge condition, with $J(x)=O(\sqrt{-x})$, and the decay of the transforms at infinity force that function to vanish by Liouville's theorem. ConsequentlyThis is the Wiener-Hopf solution of rigid half-plane diffraction with the present transform and time conventions.

Deforming to a steepest descent contour gives a residue from the incident pole precisely when $x+|y|\cot {\theta }_0<0$. At that pole $\gamma (-a)=-ib$; its residue producesThus the total geometrical optics field is $I$ plus a reflected wave in the upper reflection region; in the lower shadow region the scattered residue cancels $I$, and outside that shadow the direct incident wave remains. This gives both the amplitude and the spatial support of the geometrical optics fields rather than merely naming a reflected wave.

For $x=r\cos \theta$, $y=r\sin \theta$, $-\pi <\theta <\pi$, the saddle contribution uses ${\gamma }_l(-a)=-i\sqrt{2k_0}\cos ({\theta }_0/2)$ and ${\gamma }_u(k_s)=\sqrt{2k_0}\cos (\theta /2)$. The diffracted wave isThe total field is $I+{\varphi }_{\mathrm{GO}}+{\varphi }_d$ to this order. This far-field expression applies away from the shadow/reflection boundaries $|\theta |=\pi -{\theta }_0$. There the pole and saddle coalesce, and a uniform Fresnel integral approximation replaces the divergent separate terms by a smooth transition. The exact contour integral remains well defined.

Take $0<{\theta }_0<\pi$, put $a=k_0\cos {\theta }_0$ and $b=k_0\sin {\theta }_0>0$, and suppress the common factor $e^{i\omega t}$. The upper-half-plane incident acoustic plane wave is $e^{iax+iby}$. A rigid zero-thickness screen imposes the same normal-derivative cancellation on both faces; the scattered velocity potential is consequently odd in $y$. Across the open part $x>0$ it is continuous, so its trace there is zero. These complementary half-line data are what make the Wiener-Hopf method applicable.

Use the Fourier transform $\hat{\varphi }(\kappa ,y)={\int }_{\mathbb{R}}\varphi (x,y)e^{i\kappa x}dx$, with inverse exponential $e^{-i\kappa x}$. The outgoing Helmholtz equation solution is $\hat{\varphi }=F^{-}(\kappa )e^{-\gamma y}$ for $y>0$, where $\gamma =\sqrt{{\kappa }^2-k_0^2}$ has $\gamma =i\sqrt{k_0^2-{\kappa }^2}$ on the propagating interval and positive real part on evanescent components. The limiting absorption principle supplies the contour prescription. The transform $F^{-}$ of the trace supported on $x<0$ is analytic below the contour. Let $G^{+}$ be the analytic upper transform of the unknown derivative on $x>0$. The prescribed derivative on $x<0$ has transformSince the upper normal derivative is $-\gamma F^{-}$, the Wiener-Hopf equation isFor the Wiener-Hopf factorization, choose ${\gamma }_{+}=(\kappa -k_0{)}^{1/2}$ analytic above and ${\gamma }_{-}=(\kappa +k_0{)}^{1/2}$ analytic below. Their branch cuts go into the opposite half-planes, with ${\gamma }_{+}(-a)=i\sqrt{k_0+a}$. Dividing by ${\gamma }_{+}$ and performing pole subtraction in a Wiener-Hopf equation givesThe forcing pole is allocated to the lower analytic part; the second quotient has a removable pole and is upper analytic. Subtracting these parts leaves a common entire function. The finite-energy edge condition excludes an inverse-square-root singularity in the potential: its odd scattered trace is $O(\sqrt{|x|})$, so $F^{-}=O(|\kappa {|}^{-3/2})$. The normal derivative on the open side can be $O(x^{-1/2})$, giving $G^{+}=O(|\kappa {|}^{-1/2})$. Both sides of the entire remainder tend to zero. The Liouville theorem sets that remainder to zero, giving the explicit Wiener-Hopf solution of rigid half-plane diffraction:The scattered field on either side of the screen is thereforeIts forcing-pole residue at $\kappa =-a$ is $-i$. Deforming to the outgoing steepest descent contour crosses this pole precisely when $|\theta |>\pi -{\theta }_0$, where $x=r\cos \theta$, $y=r\sin \theta$ and $-\pi <\theta <\pi$. The residue contributes $\mathrm{sgn}(y)e^{iax-ib|y|}$. Thus, away from transition boundaries, the scattered geometrical optics field isAdding the incident acoustic plane wave gives incident plus unit specular reflection in the upper reflected sector, just the incident wave in the remaining illuminated sectors, and zero in the lower shadow sector. The geometrical shadow boundaries also follow directly by tracing a straight ray to $y=0$: the reflected ray hits the screen at $x+y\cot {\theta }_0<0$, and the transmitted incident ray crosses the open part at $x-y\cot {\theta }_0>0$.

For the diffracted contribution, the outgoing phase $-i\kappa r\cos \theta -\gamma r|\sin \theta |$ has saddle ${\kappa }_s=k_0\cos \theta$, value $-i{k}_{0}r$, and second derivative $ir/(k_0{\sin }^2\theta )$. Stationary phase therefore supplies $\sqrt{2\pi k_0/r}|\sin \theta |e^{-i{k}_{0}r+i\pi /4}$. Substituting $F^{-}({\kappa }_s)$, and using $\mathrm{sgn}(y)|\sin \theta |/\sqrt{1+\cos \theta }=\sqrt{2}\sin (\theta /2)$, givesThe full asymptotic scattered field is ${\varphi }_{\mathrm{GO}}+{\varphi }_d$. The angular edge-wave expression is not uniform on the geometrical boundaries. There is also a sign inconsistency in the supplied saddle formula: with its printed inverse exponential $e^{+i\kappa x}$, the outgoing saddle is $-k_0\cos \theta$, not $+k_0\cos \theta$. The consistent inverse convention used above gives the requested outgoing field without that ambiguity.

Take a cut-on mode with $k>0$, put $\lambda =n\pi /h$ and $\alpha ={\tan }^{-1}(\lambda /k)$, so $k=k_0\cos \alpha$, $\lambda =k_0\sin \alpha$, and initially take $n>0$. The acoustic waveguide mode splits into two acoustic plane waves:At the upper edge $y=h$, the upward-going constituent is incident from below. Use the reflected local coordinate $Y_{+}=h-y$; its incident potential becomes $(-1{)}^n{e}^{-ikx+i\lambda Y_{+}}/2$, with local incidence angle ${\theta }_0=\pi -\alpha$. At the lower edge use $Y_{-}=y+h$ and the downward-going constituent; it has the same local incident amplitude and incidence angle. At each wall the other constituent is the geometrical optics reflection of this incident half-amplitude, so it must not be counted again as an independent incident wave at the same edge.

Let $D(r)=\sqrt{2/(\pi k_{0}r)}{e}^{-i{k}_{0}r-i\pi /4}$. For a far observer, the upper and lower local observation angles are $-\theta$ and $\theta$, and their distances are $r_{+}=r-h\sin \theta +O(h^2/r)$ and $r_{-}=r+h\sin \theta +O(h^2/r)$. Applying the Wiener-Hopf solution of rigid half-plane diffraction separately givesThe sign difference is the reversal of the upper-edge transverse coordinate. Their sum is the leading independent-edge radiation from an open planar duct:Restore $e^{i\omega t}$ if a time-dependent potential is desired. This is the leading large-$k_{0}h$ approximation, neglecting further diffraction between edges, and assumes a far-field distance $r\gg k_0{h}^2$ as well as $r\gg h$. It is not a uniform approximation for a mode approaching cutoff. The zero mode follows as the grazing-incidence limit $\alpha =0$; its beam-direction limit is treated separately below.