Pulse propagation of acoustic waves scattered in a three-dimensional channel

Abstract

In a previous paper (Whitman et al 1999 Waves Random Media 9 1–11) we discussed the scattering of acoustic waves by random sound-speed fluctuations in a two-dimensional channel and presented an asymptotic form for an acoustic pulse propagating in the channel. Here we include the three-dimensional effect of transverse scattering. We find an asymptotic solution in which initially the two-dimensional mode-transfer effect is more important than the transverse scattering effect. However, for large enough propagation distances the transverse scattering effect dominates the pulse spread. In this paper we shall show the form of the pulse shape in both propagation ranges as well as in the transition regime. We shall begin with a discussion of the physics of the problem and then present a mathematical discussion.     

Propagation of light through an optically thick flat homogeneous cloud layer

Abstract

Monte Carlo simulations of light pulse transfer through a flat homogeneous cloud layer are used to compute zero, first and second moments of the light power temporal distribution at the input of a receiver placed at the cloud bottom or removed from it to a distance H in a pure atmosphere. The upper boundary of a cloud layer is considered to be illuminated by an infinitely extended δ-pulsed light beam. Transmittance of the cloud layer T and radiance angular distribution B(θ) at its bottom (these data define the zero moment or the pulse-response energy under steady illumination W), as well as mean propagation time t¯ and variance (width) (Δt)² of the pulse-response are calculated for cloud optical thicknesses τ=0, …, 50, incidence angle θ₀=0, …, 60°, single scattering albedo ω₀=0.995, …, 1.0 and the receiver view angle 2 β=40° and 90°. Monte Carlo simulation results are compared with well known Rosenberg asymptotic formulae for T and B(θ) as well as with a simple model developed for calculating t¯ and Δt. The comparison reveals high accuracy of all formulae considered for the optical thickness τ>10, …, 15. A simple procedure for computing W, t¯ and Δt parameters for a three-layer medium (cloud/pure atmosphere/sea water) is outlined.     

Propagation properties of partially coherent Laguerre-Gaussian beams in turbulent atmosphere

To study the propagation properties of partially coherent Laguerre-Gaussian (PLG) beams through turbulent atmosphere, the analytical formulas are derived for the angular width and the beam-propagation factor (M²-factor) of PLG beams by using the extended Huygens-Fresnel principle and the second-order moments of the Wigner distribution function (WDF). The corresponding numerical results are also calculated. When propagation distance increases, the angular width is found to spread faster for PLG beams with higher beam order, smaller correlation length and bigger structure constant The angular width of PLG beams decreases with increase in waist width (w₀).The M²-factor of PLG beams with higher beam order and smaller correlation length is less affected by turbulence with increase in propagation distance. The propagation properties of the M²-factor for PLG beams with the smaller structure constant are better than that with bigger structure constant . The M²-factor of PLG beams decreases with increase in the wavelength λ, and it is also less affected by turbulence for beams with higher order and smaller correlation length. Furthermore, for the PLG beams with the same beam order, the angular width and the M²-factor keep invariable in free space.     

A high-order discontinuous Galerkin method for wave propagation through coupled elastic–acoustic media

We introduce a high-order discontinuous Galerkin (dG) scheme for the numerical solution of three-dimensional (3D) wave propagation problems in coupled elastic–acoustic media. A velocity–strain formulation is used, which allows for the solution of the acoustic and elastic wave equations within the same unified framework. Careful attention is directed at the derivation of a numerical flux that preserves high-order accuracy in the presence of material discontinuities, including elastic–acoustic interfaces. Explicit expressions for the 3D upwind numerical flux, derived as an exact solution for the relevant Riemann problem, are provided. The method supports h-non-conforming meshes, which are particularly effective at allowing local adaptation of the mesh size to resolve strong contrasts in the local wavelength, as well as dynamic adaptivity to track solution features. The use of high-order elements controls numerical dispersion, enabling propagation over many wave periods. We prove consistency and stability of the proposed dG scheme. To study the numerical accuracy and convergence of the proposed method, we compare against analytical solutions for wave propagation problems with interfaces, including Rayleigh, Lamb, Scholte, and Stoneley waves as well as plane waves impinging on an elastic–acoustic interface. Spectral rates of convergence are demonstrated for these problems, which include a non-conforming mesh case. Finally, we present scalability results for a parallel implementation of the proposed high-order dG scheme for large-scale seismic wave propagation in a simplified earth model, demonstrating high parallel efficiency for strong scaling to the full size of the Jaguar Cray XT5 supercomputer.     

Uniqueness for a Class of Spatially Homogeneous Boltzmann Equations Without Angular Cutoff

We consider the 3-dimensional spatially homogeneous Boltzmann equation, which describes the evolution in time of the velocity distribution in a gas, where particles are assumed to undergo binary elastic collisions. We consider a cross section bounded in the relative velocity variable, without angular cutoff, but with a moderate angular singularity. We show that there exists at most one weak solution with finite mass and momentum. We use a Wasserstein distance. Although our result is far from applying to physical cross sections, it seems to be the first one which deals with cross sections without cutoff for non Maxwellian molecules. MSC 2000 : 82C40.     

Natural beam focusing of non-axisymmetric guided waves in large-diameter pipes

The propagation of non-axisymmetric guided waves in larger diameter pipes is studied in this paper by treating the guided waves as corresponding Lamb waves in an unwrapped plate. This approximation leads to a simpler method for calculating the phase velocities of hollow cylinder guided waves, which reveals a beam focusing nature of non-axisymmetric guided waves generated by a partial source loading. The acoustic fields in a pipe generated by a partial-loading source includes axisymmetric longitudinal modes as well as non-axisymmetric flexural modes. The circumferential distribution of the total acoustic field, also referred as an angular profile, diverges circumferentially while guided waves propagate with dependence on such factors as mode, frequency, cylinder size, propagation distance, etc. Exact prediction of the angular profile of the total field can only be realized by numerical calculations. In particular cases, however, when the wall thickness is far less than the cylinder diameter and the wavelength is smaller than or comparable to the pipe wall thickness, the acoustic field can be analyzed based on the characteristics of Lamb waves that travel along a periodic unwrapped plate. Based on this assumption, a simplified model is derived to calculate the phase velocities of non-axisymmetric flexural mode guided waves. The model is then applied to discussions on some particular characteristics of guided-wave angular profiles generated by a source loading. Some features of flexural modes, such as cutoff frequency values are predicted with the simpler model. The relationship between the angular profiles and other factors such as frequency, propagation distance, and cylinder size is obtained and presented in simple equations. The angular profile rate of change with respect to propagation distance is investigated. In particular, our simplified model for non-axisymmetric guided waves predicts that the wave beam will converge to its original circumferential shape after the wave propagates for a certain distance. A concept of "natural focal point" is introduced and a simple equation is derived to compute the 1st natural focal distance of non-axisymmetric guided waves. The applicable range of the simplified equation is provided. Industrial pipes meet the requirement of wall thickness being far less than the pipe diameter. The approximate analytical algorithms presented in this paper provides a convenient method enabling quick acoustic field analysis on large-diameter industrial pipes for NDE applications.     

On the energy distribution and angular distribution of sputtered particles

Abstract

The influence of the anisotropy of the momentum recoil density of collision cascades in solids on the energy and angular distribution of sputtered particles is considered. Various sputtering mechanisms are discussed. For some of those mechanisms this anisotropy leads to deviations from a ε^?2 asymptotic form in the energy distribution and to an asymmetry in the angular distribution as a function of the angle of incidence of the projectile.     

Characteristics of the reverberation signal at a vertical array with tonal insonification of a shallow-water basin

The angular and spectral characteristics of reverberation signals on a fixed path are studied by using vertical receiving arrays with a tonal insonification of the basin. The choice of the frequency range F ～ 200–300 Hz is caused, first, by the low propagation loss and the availability of high-power acoustic transducers for this range and, second, by the possibility to study the phenomena of sound scattering by the sea surface for comparable wavelengths of surface and acoustic waves. The data of experiments on two paths are considered: the first path is in the Barents Sea with propagation conditions corresponding to a uniform waveguide with sea depths H ～ 110–120 m, whereas the second path is in the White Sea with propagation conditions corresponding to a nonuniform waveguide where the sea depth varies from 30 to 250 m.     

Pulse propagation of acoustic waves scattered in a channel

When acoustic waves are scattered by random sound-speed fluctuations in a two-dimensional channel the energy is continually transferred between the propagating modes. In the multiple- scattering region the energy flux assumes an asymptotic form in which there is equal energy flux propagating in each mode. Here we shall make use of this well known result to show how to obtain an asymptotic form for a pulse of acoustic energy propagating in the channel. In the multiple-scattering region the speed of the acoustic waves in the pulse continually changes as the energy is transferred between the modes. The process is basically a diffusion process around the mean speed of propagation. We shall first show, using physical arguments, that the diffusion coefficient is proportional to the square root of the propagation distance times the mean free path of scattering. The theory governing the acoustic propagation in the channel is formulated in terms of modal coherence equations and we shall next give a brief review of the definitions of the coherence functions and a discussion of how the equations governing the propagation of the modal coherence functions are derived. We shall then show how the pulse shape and the relevant parameters may be obtained by solving the basic modal coherence equations at large propagation distances.     

楔形海域中三维效应对声源距离估计的影响分析*

楔形海域中声传播会出现显著的三维效应,导致水平方向上特定接收位置的声信号出现模态多次到达或者模态影区。利用射线模型,首先对楔形波导中的传播损失进行计算,分析三维效应的存在对于声波传播的影响, 结果表明,随着距离的增加,水平折射效应愈加明显,所得的三维结果与利用二维模型计算所得相差逐渐增大。然后对比二维声场和三维声场中的接收波形,并对其进行模态分离,证明在某些接收位置将会出现明显的模态多次到达或模态影区。最后通过匹配全波解以及分离出的各号简正波所对应的脉冲波形进行声源距离估计研究,并分析三维效应对于声源距离估计结果的影响。结果表明,由于三维效应的存在,使用三维模型能够得到准确估计声源距离,而利用二维模型进行估计得到的结果与实际距离相差较多,总体上偏大。     

Evaluation of a wave-vector-frequency-domain method for nonlinear wave propagation

A wave-vector-frequency-domain method is presented to describe one-directional forward or backward acoustic wave propagation in a nonlinear homogeneous medium. Starting from a frequency-domain representation of the second-order nonlinear acoustic wave equation, an implicit solution for the nonlinear term is proposed by employing the Green's function. Its approximation, which is more suitable for numerical implementation, is used. An error study is carried out to test the efficiency of the model by comparing the results with the Fubini solution. It is shown that the error grows as the propagation distance and step-size increase. However, for the specific case tested, even at a step size as large as one wavelength, sufficient accuracy for plane-wave propagation is observed. A two-dimensional steered transducer problem is explored to verify the nonlinear acoustic field directional independence of the model. A three-dimensional single-element transducer problem is solved to verify the forward model by comparing it with an existing nonlinear wave propagation code. Finally, backward-projection behavior is examined. The sound field over a plane in an absorptive medium is backward projected to the source and compared with the initial field, where good agreement is observed.     

海底爆破辐射声场的理论及数值研究

对固体基底浅海中声波传播的色散特性作了数值模拟与分析. 利用液体及固体中声传播理论导出了固体基底内爆破源在浅海中所激发声场的解析解, 在此基础上, 通过数值计算, 得出了浅海中不同空间点处的瞬态声压时域波形, 并在瞬态声压时域波形的平滑伪Wigner-Ville分布基础上, 对固体基底内爆破源在浅海中所激发瞬态声场的特性作了详细分析. 关键词： 海底爆破 声场 色散 瞬态声压     

Wave propagation in a waveguide with continuous right-angled corners: Numerical simulations and experiment measurements

This paper studied the acoustic wave propagation in a waveguide with continuous right-angled corners, with emphasis on the effect brought by the distance between the corners. The numerical analyses showed that at middle to high frequencies, the transmission loss (TL) of a multi-cornered waveguide was 2–5 dB higher than that of single-cornered and varied with frequency. To explain the performance at peaks and dips in the TL curve, analyses on eigenmodes and sound intensity distribution were conducted. The performance of multi-cornered waveguides was experimentally investigated, which fit well with the numerical results. The present study indicates that, for a waveguide with continuous corners, its acoustic performance is not simply a “summation” of two individual single-cornered ones. Both the standing wave modes and the evanescent modes between the corners lead to its complicated frequency performance.     

Angular spectrum approach for the computation of group and phase velocity surfaces of acoustic waves in anisotropic materials

The decomposition of an acoustic wave into its angular spectrum representation creates an effective base for the calculation of wave propagation effects in anisotropic media. In this method, the distribution of acoustic fields is calculated in arbitrary planes from the superposition of the planar components with proper phase shifts. These phase shifts depend on the ratio of the distance between the planes to the normal component of the phase slowness vector. In anisotropic media, the phase shifts depend additionally on the changes of the slowness with respect to the direction of the propagation vector and the polarization. Those relations are obtained from the Christoffel equation. The method employing the fast Fourier transformation algorithm is especially suited for volume imaging in anisotropic media, based on holographic detection in transmission of acoustic waves generated by a point source. This technique is compared with measurements on crystals performed by phase-sensitive scanning acoustic microscopy.     

Numerical simulations of three-dimensional nonlinear acoustic waves in bubbly liquids

This paper presents three-dimensional simulations of nonlinear propagation of ultrasonic waves through bubbly liquids, which represent the continuity of our previous works included in the numerical tool SNOW-BL. The behavior of three-dimensional nonlinear acoustic waves in bubbly liquids is analyzed by means of numerical predictions. Nonlinearity, attenuation, and dispersion due to the presence of bubbles in the liquid are taken into account. The numerical solution to the differential problem is obtained by means of a finite-difference scheme. The simulations we present here consider a homogeneous distribution of bubbles in the liquid. Results compare high and low-amplitude waves to detect the nonlinear effects of the bubbles. Results are shown for radiation and enclosure problems.     

Refinement of a formula for decay after weak coherent excitation of a sphere

A new formula has been given recently by A.A. Svidzinsky and M.O. Scully to describe the temporal evolution of the excitation function in a large sphere satisfying the Markov condition after excitation by a single photon. This formula is based on a physically reasonable Ansatz from which differential equations are inferred for the undetermined radial functions in the Ansatz. The solution to these differential equations leads to the formula for β. Numerical calculations from this formula yield a value ∼10% for the maximum probability of occupancy of secondary excited states.In this paper, we refine the formula of Svidzinsky and Scully by allowing the radial functions in the Ansatz to depend on the angular index l of the spherical Bessel functions. By using the Debye formula for the asymptotic behavior of j_l(u) for large l as well as u, we obtain differential equations in each angular sector, similar to theirs but with a dependence on l. The solution to these equations yields our improved formula, from which we calculate 17.1% for the maximum probability of secondary excited states.     

Pressure dependence of the refractive index and of the sound velocity of methanol-ethanol solution: A Brillouin scattering study

Abstract

The pressure dependence of the refractive index n(P) and of the longitudinal acoustic velocities of 4.1 methanol-ethanol solution have been measured in a diamond anvil cell up to 8.1 GPa. We utilize Brillouin scattering methods for this determination and detect, in the same back scattering configuration, acoustic wave propagation parallel to the faces of the diamond anvils as well as along the incident laser radiation direction. We also deduce that the polarizability of the fluid is reduced by about 15% over the pressure range studied.     

On the Application of Asymptotic Formulae Based on the Modified Maslov Canonical Operator to the Modeling of Acoustic Pulses Propagation in Three-Dimensional Shallow-Water Waveguides

Acoustical Physics - In this study a technique for the modeling of propagation of acoustic pulses in shallow-water waveguides with three-dimensional bottom inhomogeneities is described. The...