Sound propagation above an inhomogeneous plane: Boundary integral equation methods

A boundary integral equation method is used to compute the sound pressure emitted by a harmonic source above an inhomogeneous plane. First, the theoretical aspects of the problem (behaviour of the pressure around the discontinuities,…) are studied. Then, a comparison between theoretical levels and experimental levels obtained in an anechoic room is presented. It shows that the boundary integral equation (BIE) method is quite convenient for solving this kind of problem. Two interesting results are pointed out: (i) if only a prediction of maximum sound levels is needed, the attenuation is the same for a cylindrical source, a spherical source and N spherical sources, and so it is possible to transform some three-dimensional problems into two-dimensional ones; (ii) a numerical method of computation of the sound field above an inhomogeneous plane does not provide a correct prediction if each part of the plane is not accurately described by the boundary condition chosen.     

Sound propagation along an impedance plane

The sound field due to a point source above a plane boundary with a constant normal impedance is obtained by a double saddle point method of integration. Variations in previous studies by Ingard, by Lawhead and Rudnick and by Wenzel are clarified.     

Sound wave propagation in a vertically inhomogeneous ocean

Propagation of sound waves generated by a time dependent acoustic source in a vertically inhomogeneous ocean is considered. The effect of the solid bottom is included so that both the longitudinal and shear waves can be excited inside the bottom. The possibility of exciting a lateral shear wave by the acoustic source is also discussed. Although the results presented here are formal and general, physical interpretations have been offered whenever possible.     

The complex equivalent source method for sound propagation over an impedance plane

The sound field caused by a monopole source above an impedance plane can be calculated by using a superposition of equivalent point sources located along a line in the mirror space below the plane. Originally, such an approach for representing the half-space Green's function was described by Sommerfeld at the beginning of the last century, in order to treat half-space problems of heat conduction. However, the representation converges only for masslike impedances and cannot be used for the more important case of reflecting planes with springlike surface impedances. The singular part of the line integral can be transformed into a Hankel function, which shows that surface waves are contained in the whole solution. Unfortunately, this representation suffers from the lack of validity at certain receiver points and from restrictions on wave number and impedance range to ensure the necessary convergence. The main idea of the present method is to use also a superposition of equivalent point sources, but to allow that these sources can be located at complex source points. The corresponding form of the half-space Green's function is suitable for both masslike and springlike surface impedances, and can be used as a cornerstone for a boundary element method.     

Sound and light propagation in a weakly inhomogeneous atmosphere

Small-scale inhomogeneities caused by atmospheric turbulence have a considerable effect on sound and light propagation, producing the fluctuations of these wave fields. V A Krasilnikov [1, 2] performed experiments on phase and amplitude variations of a sound wave propagating through the atmosphere. Fluctuations of light-wave parameters occur, for example, in the well known phenomenon of star scintillation, apparently strongly connected with turbulent irregularities of the atmospheric temperature field [3, 4]‡.

Some calculations of phase (arrival angle) and amplitude fluctuations for a wave propagating through a turbulent medium are described by Krasilnikov [1, 3, 8]. All these and similar calculations are based on the geometrical optics (acoustics) approximation, which may be the reason for disagreement between calculation and experimental data in some cases. Thus, for example, amplitude fluctuations in the geometrical approximation turn out to be proportional to the distance of propagation through a turbulent medium to the power of 3/2. However, observations usually show much slower fluctuation growth.

This paper represents an attempt to consider the problem of amplitude and phase variations for a scalar wave field in terms of more general equations including some diffraction effects. Incidentally, the range of validity of geometrical optics theory becomes clear.     

Normal mode solution for low-frequency sound propagation in a downward refracting atmosphere above a complex impedance plane.

The development of the fast field and parabolic equation solutions to the wave equation has made it possible to solve for the combined effects of refraction in a layered atmosphere and the interaction of sound with a complex impedance ground surface. In many respects the numerical methods have advanced beyond our understanding of the basic phenomena. In an earlier study [J. Acoust. Soc. Am. 89, 107-114 (1991)], the residue series solution for upward refraction was investigated and provided insight into the nature of the interaction of refraction and ground reflection. In this paper results are presented of a similar normal mode solution for downward refraction above a complex impedance ground surface. This model is used to investigate when the surface wave is excited for downward refraction conditions and to develop criteria for the maximum range of cylindrical decay as a function of phase and magnitude of the ground impedance and the magnitude of the sound velocity gradient.     

Sound propagation in inhomogeneous waveguides with sound-speed profiles using the multimodal admittance method

The multimodal admittance method and its improvement are presented to deal with various aspects in underwater acoustics, mostly for the sound propagation in inhomogeneous waveguides with sound-speed profiles, arbitrary-shaped liquid-like scatterers, and range-dependent environments. In all cases, the propagation problem governed by the Helmholtz equation is transformed into initial value problems of two coupled first-order evolution equations with respect to the modal components of field quantities(sound pressure and its derivative), by projecting the Helmholtz equation on a constructed orthogonal and complete local basis. The admittance matrix, which is the modal representation of Direchlet-to-Neumann operator, is introduced to compute the first-order evolution equations with no numerical instability caused by evanescent modes. The fourth-order Magnus scheme is used for the numerical integration of differential equations in the numerical implementation. The numerical experiments of sound field in underwater inhomogeneous waveguides generated by point sources are performed. Besides, the numerical results computed by simulation software COMSOL Multiphysics are given to validate the correction of the multimodal admittance method. It is shown that the multimodal admittance method is an efficient and stable numerical method to solve the wave propagation problem in inhomogeneous underwater waveguides with sound-speed profiles, liquid-like scatterers, and range-dependent environments. The extension of the method to more complicated waveguides such as horizontally stratified waveguides is available.     

Sound propagation in an unstable atmospheric layer

The equation of sound propagation in an unstable medium produced by the presence of supersaturated water vapor in it, which can appear in a hurricane area, is derived. This equation takes into account the effects of sound velocity dispersion, amplification, damping, and nonlinear effects. Some solutions to this equation are presented, illustrating the specific features of sound propagation in an unstable medium. Published in Akusticheskiĭ Zhurnal, 2007, Vol. 53, No. 3, pp. 477–480.     

The transverse acoustic impedance of an inhomogeneous viscous fluid

The transverse acoustic impedance Z of a fluid can be determined from measurements of the complex reflection coefficient of a transverse ultrasonic wave incident on a plane solid/fluid interface. Inhomogeneities in the fluid close to the solid surface may have a significant effect on the measured values of Z. We have derived the Riccati equation which determines Z in an inhomogeneous viscous fluid, using transmission line theory. This equation was integrated numerically to obtain the impedance of viscous films and for inhomogeneities due to healing lengths and van der Waals forces near the solid/fluid interface. The results show that the measurement of both the real and imaginary parts of Z can be a powerful technique for investigating any inhomogeneities which occur on a length scale comparable with the viscous penetration depth in the fluid.     

Sound propagation above a porous road surface with extended reaction by boundary element method

Acoustic impedance of an absorbing interface is easily introduced in boundary element codes provided that a local reaction is assumed. But this assumption is not valid in the case of porous road surface. A two-domain approach was developed for the prediction of sound propagation above a porous layer that takes into account the sound propagation inside the porous material. The porous material is modeled by a homogeneous dissipative fluid medium. An alternative to this time consuming two-domain approach is proposed by using the grazing incidence approximate impedance in the traditional single-domain boundary element method (BEM). It can be checked that this value is numerically consistent with the surface impedance calculated at the interface from the pressure and surface velocity solutions of the two-domain approach. The single-domain BEM introducing this grazing incidence impedance is compared in terms of sound attenuation with analytical solutions and two-domain BEM. The comparison is also performed with the single-domain BEM using the normal incidence impedance, and reveals a much better accuracy for the prediction of sound propagation above a porous interface.