Intensity and Phase Reinitialization
At a material discontinuity, the intensity of the reflected and refracted rays can be computed using the material properties of the media in the adjacent domains:
(8-2)
where ρ is the density of the medium and cc is the complex speed of sound,
where
c (SI unit: m/s) is the real speed of sound in the medium,
ω (SI unit: rad/s) is the angular frequency of the acoustic ray, and
α (SI unit: 1/m) is the attenuation coefficient.
The subscripts 1 and 2 denote the sides of the boundary into which the reflected and refracted ray propagate, respectively; see Ref. 4.
At a Wall, the reflection coefficient R can be defined explicitly or in terms of the absorption coefficient α or characteristic impedance of absorber Z1. In terms of the absorption coefficient,
Where the principal square root is used for purposes of phase reinitialization. In terms of the characteristic impedance of absorber,
The Wall feature also includes built-in options to compute the reflection coefficient for fluid-fluid interfaces, fluid-solid interfaces, and a fluid layer adjacent to a semi-infinite fluid domain. When treating the Wall boundary as a Fluid-fluid interface, the reflection coefficient is computed using Equation 8-2 as if the wall were a material discontinuity where the fluid properties of the adjacent medium are specified, except that no refracted ray is produced.
If the Wall boundary is modeled as a Fluid-solid interface, the reflection coefficient is instead defined as
where
where the subscripts p and s refer to the propagation of compressional and shear waves in the adjacent solid domain, respectively. For example, θp,t is the refraction angle computed using the compressional complex speed of sound in the adjacent solid domain.
If the Wall boundary is modeled as a Layered fluid half-space, the boundary is treated as a thin layer of one fluid backed by a semi-infinite domain of a second fluid. The reflection coefficient is
where the subscripts 1 and 2 refer to the properties of the thin layer and the semi-infinite fluid domain, respectively. For example, θt,1 is the angle of refraction in the thin layer. The angle is the phase delay in the thin layer,
where h1 (SI unit: m) is the layer thickness.
If the Rayleigh roughness model is used, the reflected intensity is multiplied by an additional factor to account for surface roughness:
where k (SI unit: rad/m) is the wave vector magnitude of the ray and σ (SI unit: m) is the RMS roughness of the surface.
If the phase of the ray is also computed, the phase shift of the reflected ray at a boundary is
At material discontinuities, the phase of the refracted ray is unchanged.