Fluid-loaded membrane edge scattering Created 2026-10-06 Updated 2026-10-07
A pinned semi-infinite membrane scatters an incoming guided mode into reflected guided modes and outgoing sound. For unit incident displacement the scattered endpoint value is . Its Half-range Fourier transform therefore has nonzero endpoint terms even though the total displacement vanishes. The Wiener-Hopf equation kernel is with the outgoing branch .
Fluid-loaded plate. Use complex amplitudes with time dependence and take the velocity potential in the fluid above the elastic plate to be . Let and be the undisturbed mass density and speed of sound. The linearized Euler equations give ; the kinematic boundary condition is . 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 by
The Sommerfeld radiation condition selects , with and . Here and in the diffraction calculation below, the outgoing square root is specified by first taking . Draw a vertical cut upwards from and another downwards from , and put
Thus is analytic in the upper half-plane and in the lower half-plane. In the real-frequency limit,
The Fourier transform contour 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 , so . Eliminating the fluid pressure from the bending balance gives the line-forced fluid-loaded bending plate kernel
For harmonic line force amplitude at the origin, the plate displacement and fluid velocity potential are therefore
An arbitrary localized loading replaces by its Fourier transform. This representation both calculates the response and identifies its different wave components through poles and branch points.
For an observer , , with , the method of steepest descent uses the saddle point . At fixed angle away from grazing, the acoustic far field is
The spreading is that of cylindrical acoustic waves. The phase velocity parallel to the plate is : components with are supersonic in that sense and radiate, whereas those with 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 . Indeed, , , and
The residue theorem then gives an outgoing displacement contribution . It transports energy along the plate while its fluid field decays as . 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 . Its upward normal velocity is and its pressure is . 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 with prefactor . Integrating the radial acoustic intensity of the far field over a semicircle gives the same expression. The input power also feeds guided waves along the plate, so it must not be identified with alone.
Rigid half-plane diffraction. Suppress the time factor and assume . Write , , and denote the incident spatial field by . 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 be this jump, supported on . The transformed scattered field is
The 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 .
On the negative axis this derivative must cancel . Since the Half-range Fourier transform of is , the Wiener-Hopf equation has the form
where is the transform of an unknown function supported on and analytic in the lower half-plane. The incident pole lies below the Fourier contour. Divide by 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 , and the decay of the transforms at infinity force that function to vanish by Liouville's theorem. Consequently
This 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 . At that pole ; its residue produces
Thus the total geometrical optics field is plus a reflected wave in the upper reflection region; in the lower shadow region the scattered residue cancels , 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 , , , the saddle contribution uses and . The diffracted wave is
The total field is to this order. This far-field expression applies away from the shadow/reflection boundaries . 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 , put and , and suppress the common factor . The upper-half-plane incident acoustic plane wave is . A rigid zero-thickness screen imposes the same normal-derivative cancellation on both faces; the scattered velocity potential is consequently odd in . Across the open part 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 , with inverse exponential . The outgoing Helmholtz equation solution is for , where has on the propagating interval and positive real part on evanescent components. The limiting absorption principle supplies the contour prescription. The transform of the trace supported on is analytic below the contour. Let be the analytic upper transform of the unknown derivative on . The prescribed derivative on has transform
Since the upper normal derivative is , the Wiener-Hopf equation is
For the Wiener-Hopf factorization, choose analytic above and analytic below. Their branch cuts go into the opposite half-planes, with . Dividing by and performing pole subtraction in a Wiener-Hopf equation gives
The 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 , so . The normal derivative on the open side can be , giving . 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 therefore
Its forcing-pole residue at is . Deforming to the outgoing steepest descent contour crosses this pole precisely when , where , and . The residue contributes . Thus, away from transition boundaries, the scattered geometrical optics field is
Adding 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 : the reflected ray hits the screen at , and the transmitted incident ray crosses the open part at .
For the diffracted contribution, the outgoing phase has saddle , value , and second derivative . Stationary phase therefore supplies . Substituting , and using , gives
The full asymptotic scattered field is . 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 , the outgoing saddle is , not . The consistent inverse convention used above gives the requested outgoing field without that ambiguity.
Suppress the common time factor and invert the stated Fourier transform with . Use limiting absorption principle , and choose with positive real part for the evanescent spectrum. For its outgoing continuation is . The pole lies above the inversion contour and below it. These prescriptions fix both radiation and pole signs.
The even interior scattered transform solves , so it is . The exterior transform is . Let be the transform of their common upward normal derivative on . On the rigid plates set both derivatives to zero; on continuity requires them to agree. Hence this derivative is supported on the negative half-line and is analytic below the transform contour. It follows that
The difference between exterior and interior scattered traces on is , where . At the open negative half-line the total traces must match, so this difference equals the incident wave . Its negative-half transform is . The positive-half trace difference has an unknown upper-half-plane transform , giving the Wiener-Hopf equation
Factor , divide by , and subtract the forcing pole using :
The left side is analytic below and the right side above; the pole in the last quotient is removable. Analytic continuation makes their common value entire. The finite-energy edge condition and algebraic factor growth make this entire remainder vanish: the allowed normal-velocity edge singularity is at most inverse-square-root, so its half-line transform is , while . The same edge condition excludes a growing polynomial remainder. This pole subtraction in a Wiener-Hopf equation yields
The corresponding exterior result for the Wiener-Hopf solution for an open parallel-plate waveguide is
Both forms are useful: the first makes the common normal derivative transparent, and the second removes the apparent factor singularity at a grazing saddle.
The final unheaded request has a sign error in its printed exponent. To see this from the derived solution rather than an outgoing-decay assertion alone, simplify the interior transform to
Near , and for fixed . Thus
For , inversion of the first term gives , cancelling the continued incident potential. The remaining branch-point term uses on the specified contour, so it is . Its inverse transform at is proportional to . Therefore the upstream decay of an open parallel-plate waveguide is
The amplitude decays as , not . The positive exponent is present in the original PDF, so it cannot be proved under the stated outgoing problem. The incident pole cancellation and the explicit branch expansion give the corrected conclusion.
Use a Fourier transform parallel to the boundary. Splitting its traces into right- and left-supported parts gives Half-range Fourier transforms and , analytic respectively above and below a common inversion contour. A small positive imaginary part of the wavenumber implements the limiting absorption principle and separates the incident poles and outgoing branch points. Solve the transformed differential equation normal to the boundary using its outgoing or decaying solution. The boundary conditions then yield a Wiener-Hopf equation of the form
Factor its Wiener-Hopf kernel as , with each factor analytic and nonzero in its designated half-plane. Divide by and split , allocating poles according to the incident-wave prescription. Then
The two sides continue to a common entire function. Their growth at infinity, inherited from the field's edge behavior, determines that function, usually a constant or a polynomial; in a decaying case it is zero. Solve for the unknown transforms and apply Fourier inversion, preserving the outgoing contour prescription.
The principal analytic results used are as follows. The Paley–Wiener theorem and its half-line Laplace-transform version relate support on a half-line to holomorphy in a half-plane, with growth or square-integrability bounds. The Cauchy integral formula and Sokhotski–Plemelj theorem provide additive splitting: a Cauchy integral of a sufficiently decaying jump has upper and lower boundary values with that prescribed difference. Applying the same construction to a continuous logarithm of a nonvanishing scalar Wiener-Hopf kernel gives Wiener-Hopf factorization when its winding number is zero; a nonzero index requires an explicit index factor. The identity theorem gives analytic continuation across the common strip. Finally Liouville's theorem makes a bounded entire function constant; its polynomial-growth extension bounds the possible degree. The physical edge and radiation conditions select the remaining constants. Zeros, branch cuts, nonzero index and insufficient decay must be accounted for rather than discarded during splitting.
Suppress the common time factor and use the Fourier transform convention
A right Half-range Fourier transform is analytic above its convergence line and a left transform below it. For the limiting absorption principle with this time convention take , , , before taking . Set
This is the outgoing acoustic square-root branch: at a propagating real wavenumber, . Put and . Across the fluid part of , the field and its normal derivative are continuous; on the sheet both sides have the same normal velocity. The outgoing Helmholtz equation solutions therefore have equal normal derivatives and opposite scattered traces:
The jump is supported on , and so is a minus transform.
The linear fluid equation and the sheet velocity imply on either face. The incident contribution to the left derivative transform is . Writing , it follows that
The upward net pressure is lower minus upper, hence in transform space. If , the pinned endpoint gives
Combining these identities gives the Wiener-Hopf equation
This is the plane-wave forcing of a pinned elastic half-sheet.
Let and assume initially that the distinguished points are distinct. Besides the branch points , has simple poles at . Its physical-sheet zeros solve the fluid-loaded dispersion relation
For real positive frequency there is exactly one positive real root : the left side increases strictly from zero to infinity on that interval. Equivalently, with , . Its other roots have negative real part and do not belong to this outgoing sheet. Continue with small absorption; lie below the real line and their negative partners above it. The Wiener-Hopf factorization assigns the lower branch cut and pole to , and the upper branch cut and pole to . The continued has a zero at , and at ; these are poles of their reciprocals, rather than poles of the factors. Branch cuts can be drawn away from the real line into the corresponding half-planes.
For this standard assignment a common strip of analyticity and nonvanishing is
The strip shrinks to the real contour as absorption is removed. The factors continue meromorphically beyond their designated half-planes; the allocation of the poles and branch points is part of the radiation condition. The incident pole belongs above the inversion contour. For incidence with the real contour is a direct choice. Other nongrazing angles are obtained by analytic continuation with the incident pole kept on its prescribed side, deforming the contour locally if necessary. Coincident distinguished points and grazing incidence are treated by the corresponding limits, not by dividing by a coincident pole.
Divide the Wiener-Hopf equation by and use pole subtraction in a Wiener-Hopf equation. With , the minus part of is
Indeed these are exactly the incident-pole and upper bare-sheet-pole residues. At the lower bare-sheet pole , the forcing pole is canceled by the pole of before splitting. Thus is analytic on the plus side of the contour. The separated equation is
With the stipulated ,
The scattered field is consequently
The full time-dependent density perturbation is the incident field plus . The two traces at follow by limiting from above and below. A constant rescaling of the two factors does not change this field.
To determine the endpoint slope, reconstruct the sheet transform using the fluid condition:
At the continued has residue . Therefore the derivative transform generically has the lower-plane pole
This pole produces a forbidden incoming bare-sheet wave in the reconstructed displacement: travels towards the endpoint from . A minus transform of the specified outgoing solution must be analytic there. The incoming-pole cancellation at a pinned membrane edge is precisely
For generic parameters this determines
The lower pole is removed; the upper zero of can still produce the legitimate outgoing evanescent acoustic surface wave, whose phase travels to negative . This distinction separates the spurious bare-sheet wave from the coupled fluid-loaded mode. No explicit factorization has been used.
Suppress the common time factor , take , and fix the spatial Fourier transform convention
A plus transform is supported on and analytic above ; a minus transform is supported on and analytic below it. For the outgoing radiation condition with this time convention, initially take , , and pass to the limit at the end. Thus lies below the contour and above it.
Choose to have positive real part on the real transform line. Its branch cuts run from into the lower half-plane and from into the upper half-plane, without crossing ; downward and upward vertical rays are suitable. In the zero-absorption limit, for real and between the branch points. This ensures that represents decay or outgoing radiation, rather than an incoming exterior field.
Evenness in and the Helmholtz equation give the transformed fields
At , the one-sided normal derivatives of the scattered field agree for , because that part of the interface is open. For they both vanish by rigidity. Their common trace is therefore a minus function, denoted . Consequently
The jump in the total scattered transform is
For , continuity of the total mass density requires the scattered jump to cancel the incident jump, so . Its minus transform is . The unknown plate-side jump is the plus function . Thus , and the Wiener-Hopf equation follows:
Use the Wiener-Hopf factorization , with factors analytic and nonzero in their designated half-planes. Their analytic continuations allocate outgoing modal zeros to below the contour, and the opposite zeros to above it. Set . Multiplication by and pole subtraction give
The pole in the upper expression is removable by its numerator. The two expressions analytically continue to the common entire function, which is zero under the stipulated edge/growth assumption. Hence
and the transformed fields are
A constant reciprocal rescaling of the factors does not alter or the physical field.
Inside the guide, and are even entire functions of , while is analytic in the lower half-plane. Thus continuation across the lower branch cut changes none of the interior transform: that cut is removable. The exterior expression retains the cut, corresponding to radiation into the open exterior.
For , the inverse-transform contour closes downwards, clockwise. Away from modal cutoffs, its enclosed singularities are the outgoing simple poles
with positive real part for propagating modes and negative imaginary part for decaying modes; is the root. The opposite roots lie above the contour or cancel against zeros of . Symmetry excludes odd transverse modes. At , , and differentiation of gives
Each inverse-transform contribution is times its residue. Therefore the outgoing modes of an open rigid acoustic waveguide are
The original time factor multiplies this expression. Thus every cut-on mode propagates in the positive direction; cut-off modes are evanescent waves decaying into the guide, not additional backward waves. Strictly, a complete mode sum includes these evanescent modes as well as propagating ones. The displayed simple-pole amplitudes apply away from exact cutoffs; cutoff values use the outgoing limiting-absorption continuation before taking the limit. The reflected plane-wave amplitude relative to the unit incident wave is above. It has the expected negative sign in the long-wavelength leading kernel approximation .
Linearity permits subtraction of the incident acoustic membrane wave. The scattered pressure solves the homogeneous Helmholtz equation, and the kinematic relation defines its line displacement on both halves. On , subtracting the incident elastic membrane equation gives the scattered dynamic equation. On , there is no elastic membrane, so the total pressure jump must vanish, giving . The pinned endpoint has , hence , rather than zero.
Use full and Half-range Fourier transforms with the same convention as part (a), and put
Write for the left transform of the upper scattered pressure, for its right transform, and for the corresponding transforms of . Outgoing antisymmetry in implies
The second relation comes from the prescribed pressure cancellation on .
The endpoint terms in the left Fourier transform of are crucial:
Thus . Eliminating gives the Wiener-Hopf equation
There is a defect in the printed right side: with the stated incident exponential, pressure jump, and pinned total displacement, it omits the incident endpoint term and has the opposite sign on the incident-pressure term. The displayed corrected equation follows directly from both boundary traces. The given right side cannot be derived with these definitions. Using , an especially useful equivalent form is
The correction is necessary for the reconstructed pressure to satisfy the original physical boundary conditions.
For Wiener-Hopf factorization, take with the plus factor analytic and nonzero in the upper half-plane and the minus factor analytic and nonzero in the lower half-plane. Singularities listed next refer to continuation out of each factor's own analytic half-plane. Let , so and lies below the real axis. Generically inherits the pole at and the square-root branch point at , while inherits the pole at and the branch point at . The kernel is finite but nonanalytic at the acoustic branch points: its local behavior is a constant plus a square-root term, not an inverse-square-root divergence. Coincident branch points and poles require a limiting treatment.
Zeros of are the fluid-loaded dispersion relation roots. Lower-half-plane roots, including the incoming guided root , belong to the continued ; upper-half-plane roots, including the reflected root for the even kernel, belong to the continued . The roots represent membrane-guided modes; the branch cuts represent radiated acoustic waves. Only roots on the selected outgoing sheet are included, not spurious roots created by squaring the dispersion relation. No explicit factorization is required.
Here is an explicit solution in terms of those factors. Define
The pole of at has residue ; the pole at is canceled by the pole of . Thus have the required respective analyticity. Divide the corrected Wiener-Hopf equation by and split its right side as , where
The difference quotient in is removable at and is analytic below. Moving the plus and minus terms to opposite sides yields the common entire function. With the assumed , the solution is
Combining with the known gives the full scattered pressure transform
The requested pressure integral, for , is consequently
The real contour uses the causal continuation ; its undamped limit retains the induced pole and branch-cut prescriptions. The residue theorem shows that upper-half-plane zeros of yield left-going scattered elastic membrane modes. On the right, the lower incident-pole residue of the scattered integral cancels at the open line, as required.
Finally reconstruct the left scattered displacement:
For a generic unspecified , this has a lower-half-plane pole at the bare elastic membrane wavenumber . Such a pole is incompatible with analyticity of an outgoing left-supported scattered displacement: it would represent an additional right-going incoming elastic membrane contribution. It must be removed. Its residue is , so the incoming-pole cancellation at a pinned membrane edge condition is
Since is linear in , this fixes the endpoint slope generically. Explicitly it is
The physical lower incident pole of the total displacement is already prescribed by ; it must not be confused with this removable spurious pole of the scattered field.
For identical unit-density, unit-sound-speed fluids above and below a sheet on , use time and transform . The outgoing acoustic square-root branch obeys . If is the pressure jump upper minus lower, and the right Half-range Fourier transform of the scattered normal derivative, the Wiener-Hopf equation is
The total displacement is pinned, so its transformed second derivative is . The bare-sheet pole in reconstructed displacement must cancel: . This is an incoming-pole cancellation at a pinned membrane edge condition, distinct from the outgoing coupled-mode pole at the zero of .
Wiener-Hopf kernel 2026-10-14
The scalar coefficient multiplying an unknown half-line transform in a Wiener-Hopf equation is its kernel. This is a transform-domain function, distinct from the null space of a linear map. Its analytic zeros, poles and winding number govern Wiener-Hopf factorization. For mixed scalar Helmholtz half-plane diffraction a typical kernel is the outgoing vertical wavenumber .