Method for solving the problems of multiple scattering by several bodies in a homogeneous unbounded medium

A method is developed for solving problems of multiple scattering by an aggregate of bodies in a homogeneous unbounded medium. For this purpose, the problem on the multiple scattering produced by two bodies in the field of a plane wave is first considered under the assumption that the initial unperturbed scattering amplitudes of both scatterers are known. The solution is constructed by considering plane waves multiply rescattered by the scatterers. Integral equations are obtained that allow one to calculate the resulting scattering amplitude of each scatterer and the combined scattering amplitude of the system of two scatterers. It is shown that knowledge of the solution to this problem is sufficient to solve the problem on the scattering field of a system consisting of an arbitrary number of scatterers. Expressions for the scattering amplitude in the case of an arbitrary primary field are presented. The relationship between the integral equations describing the multiple scattering in a homogeneous space and the multiple scattering by a single scatterer located near an interface is demonstrated. Approximate expressions are given for calculating the scattering amplitude in the case of multiple scattering.     

Acoustical scattering by multilayer spherical elastic scatterer containing electrorheological layer

A computational procedure for analyzing acoustical scattering by multilayer concentric spherical scatterers having an arbitrary mixture of acoustic and elastic materials is proposed. The procedure is then used to analyze the scattering by a spherical scatterer consisting of a solid shell and a solid core encasing an electrorheological (ER) fluid layer, and the tunability in the scattering characteristics afforded by the ER layer is explored numerically. Tunable scatterers with two different ER fluids are analyzed. One, corn starch in peanut oil, shows that a significant increase in scattering cross-section is possible in moderate frequencies. Another, fine poly-methyl methacrylate (PMMA) beads in dodecane, shows only slight change in scattering cross-sections overall. But, when the shell is thin, a noticeable local resonance peak can appear near ka=1, and this resonance can be turned on or off by the external electric field.     

Multilevel fast multipole algorithm for acoustic wave scattering by truncated ground with trenches

The multilevel fast multipole algorithm (MLFMA) is extended to solve for acoustic wave scattering by very large objects with three-dimensional arbitrary shapes. Although the fast multipole method as the prototype of MLFMA was introduced to acoustics early, it has not been used to study acoustic problems with millions of unknowns. In this work, the MLFMA is applied to analyze the acoustic behavior for very large truncated ground with many trenches in order to investigate the approach for mitigating gun blast noise at proving grounds. The implementation of the MLFMA is based on the Nystrom method to create matrix equations for the acoustic boundary integral equation. As the Nystrom method has a simpler mechanism in the generation of far-interaction terms, which MLFMA acts on, the resulting scheme is more efficient than those based on the method of moments and the boundary element method (BEM). For near-interaction terms, the singular or near-singular integrals are evaluated using a robust technique, which differs from that in BEM. Due to the enhanced efficiency, the MLFMA can rapidly solve acoustic wave scattering problems with more than two million unknowns on workstations without involving parallel algorithms. Numerical examples are used to demonstrate the performance of the MLFMA with report of consumed CPU time and memory usage.     

中红外散射的基质折射率效应研究

本文利用Mie散射理论对中红外波段三种不同类型散射体的散射和消光进行了研究,发现基质在三种不同类型散射体的散射和消光中所起作用完全不同。对于无吸收的高折射率散射体,基质折射率越大面积散射比(或面积消光比)越低;对于反常色散散射体,基质折射率越大面积散射比和面积肖光比越高;而对于金属散射体,基质折射率的大小对于散射和消光影响极小。     

On acoustic scattering by a shell-covered seafloor

Acoustic scattering by the seafloor is sometimes influenced, if not dominated, by the presence of discrete volumetric objects such as shells. A series of measurements of target strength of a type of benthic shelled animal and associated scattering modeling have recently been completed (Stanton et al., "Acoustic scattering by benthic and planktonic shelled animals," J. Acoust. Soc. Am., this issue). The results of that study are used herein to estimate the scattering by the seafloor with a covering of shells at high acoustic frequencies. A simple formulation is derived that expresses the area scattering strength of the seafloor in terms of the average reduced target strength or material properties of the discrete scatterers and their packing factor (where the reduced target strength is the target strength normalized by the geometric cross section of the scatterers and the averaging is done over orientation and/or a narrow range of size or frequency). The formula shows that, to first order, the backscattering at high acoustic frequencies by a layer of shells (or other discrete bodies such as rocks) depends principally upon material properties of the objects and packing factor and is independent of size and acoustic frequency. Estimates of area scattering strength using this formula and measured values of the target strength of shelled bodies from Stanton et al. (this issue) are close to or consistent with observed area scattering strengths due to shell-covered seafloors published in other papers.     

A successive scattering methodology with application to two cylinders in the Fresnel region of each other

We present an efficient approach to compute the second-order scattering of an electromagnetic wave by two discrete scatterers in proximity to each other. Such a two-body system represents the simplest canonical arrangement to address near-field volume scattering phenomena in microwave remote sensing models of vegetation. Using an analytical wave-based approach, a successive scattering methodology is employed to derive the first interaction term in multiple scattering by two arbitrary scatterers in terms of their transition operators. The general formulation is applied to find the second-order bistatic scattering amplitude for a pair of finite length thin cylinders at arbitrary interaction distances using the exact Green's function. To improve computational efficiency, the solution is then specialized to the Fresnel region. These second-order bistatic scattering amplitude results are in agreement with the exact Green's function model when the scatterers are in the Fresnel region of each other. Additionally, it is demonstrated that using the far field approximation in the Fresnel region can yield significant deviations from the exact results. The Fresnel model, unlike the far field approximation, accurately predicts the scattering amplitude peak values and null locations, and is suited to fast solutions in realistic canopy simulations.     

Solution of problem of electromagnetic waves scattering on inhomogeneously layered scatterers using pattern equation method

In this paper, the pattern equation method (PEM) is used to solve the problems of electromagnetic waves scattering by inhomogeneously layered scatterers. The numerical algorithm based on PEM for multilayered scatterers of arbitrary geometry is developed. The numerical results are presented for magneto-dielectric objects coated with magneto-dielectric materials to demonstrate the stability, accuracy, and efficiency of the proposed method. A comparison of the results obtained by the present method and by the other methods are shown and the limits of applicability of method are outlined.     

Application of the mutual-interaction method to a class of two-scatterer systems: I. Two discrete scatterers

In a recent paper we developed a formalism that fully accommodates the mutual interactions among scatterers separable by parallel planes. The total fields propagating away from these planes are the unknowns of a system of difference equations. Each scatterer is characterized by a scattering function that expresses the scattered wave amplitude as a function of the incident and scattered wavevectors for a unit-amplitude plane wave scattered from the object in isolation. This function can be derived completely from the scattered far field with the help of analytic continuation. For a two-scatterer system the mutual-interaction equations reduce to a single Fredholm integral equation of the second kind. It turns out that analytic solutions are tractable for those scattering functions that are Dirac deltas or a sum of products of separable functions of the incident and scattered wavevectors. Scattering functions for planes and isotropic scatterers, as well as electric and magnetic dipoles all possess this property and are considered. The exact scattering functions agree with results obtained by analytic continuation. This paper consists of two parts. Part I derives analytic solutions for two discrete scatterers (isotropic scatterers. electric dipoles, magnetic dipoles). Part II is devoted to scattering from an object (isotropic or dipole scatterer) near an interface separating two semi-infinite uniforn-media. Because the results in this paper are exact, the effects of near-field interactions can be assessed. The forms of the scattering solutions can be adapted to objects that are both radiating and scattering.     

水下任意刚性散射体对Bessel波的散射特性分析

本文利用Bessel波的谐波展开式, 采用T矩阵方法的推导思路, 建立了水下任意刚性散射体在Bessel波照射下的声散射场计算公式. 以水下刚性椭球体和两端附连半球的刚性圆柱体为例, 计算了在不同波锥角β 下的反向散射形态函数, 同时, 依据镜反射波和绕行波的干涉物理模型, 给出了预报Bessel波照射下的反向散射形态函数峰峰间隔值的计算模型. 仿真结果表明本文提出的Bessel波照射下反向散射形态函数峰峰间隔值预报方法是准确有效的, 同时也说明, 本文建立的基于T矩阵法计算水下任意刚性散射体在Bessel波束下的声散射场方法是有效的, 这拓展了T矩阵法的应用领域.     

Application of the mutual-interaction method to a class of two-scatterer systems: I. Two discrete scatterers

Abstract

In a recent paper we developed a formalism that fully accommodates the mutual interactions among scatterers separable by parallel planes. The total fields propagating away from these planes are the unknowns of a system of difference equations. Each scatterer is characterized by a scattering function that expresses the scattered wave amplitude as a function of the incident and scattered wavevectors for a unit-amplitude plane wave scattered from the object in isolation. This function can be derived completely from the scattered far field with the help of analytic continuation. For a two-scatterer system the mutual-interaction equations reduce to a single Fredholm integral equation of the second kind. It turns out that analytic solutions are tractable for those scattering functions that are Dirac deltas or a sum of products of separable functions of the incident and scattered wavevectors. Scattering functions for planes and isotropic scatterers, as well as electric and magnetic dipoles all possess this property and are considered. The exact scattering functions agree with results obtained by analytic continuation. This paper consists of two parts. Part I derives analytic solutions for two discrete scatterers (isotropic scatterers. electric dipoles, magnetic dipoles). Part II is devoted to scattering from an object (isotropic or dipole scatterer) near an interface separating two semi-infinite uniforn-media. Because the results in this paper are exact, the effects of near-field interactions can be assessed. The forms of the scattering solutions can be adapted to objects that are both radiating and scattering.     

Time-reversal analysis for scatterer characterization

A new application of time-reversal processing of wave scattering data permits characterization of scatterers by analyzing the number and nature of the singular functions (or eigenfunctions) associated with individual scatterers when they have multiple contributions from monopole, dipole, and/or quadrupole scattering terms. We discuss acoustic, elastic, and electromagnetic scattering problems for low frequencies. Specific examples for electromagnetic scattering from one of a number of small conducting spheres show that each sphere can have up to six distinct time-reversal eigenfunctions associated with it.     

Focused antenna array in a strongly inhomogeneous medium

A focused acoustic antenna array is considered in a strongly inhomogeneous stationary medium. An opportunity is indicated to determine the coordinates of a number of objects by active location. It is assumed that, in insonifying the objects by a wave with an arbitrary wave front, they scatter spherical waves and are sufficiently separated in distance and angle to be resolved by the same array in a homogeneous medium. The procedure of determining the coordinates of the objects involves a wave front inversion for distinguishing between the signals from different objects. The coordinates are determined by estimating the parameters for each individual object. The parameter estimation procedure is shown to provide a high efficiency of extracting the argument of a complex signal. The results of the numerical modeling and solution of the problem are presented.     

Modeling of the exact solution of the inverse scattering problem by functional methods

Relations between different functional algorithms for solving the inverse scattering problem are analyzed. It is shown that the Rose algorithm does not provide a unique solution, but can be used as a means to improve the interference resistance in reconstruction algorithms that provide unique restoration of scatterer characteristics. The possibility of unique reconstruction of refractive-absorbing scatterers by the modified Rose algorithm, which includes the Sokhotsky equation, is illustrated. Results of numeric simulation of the Novikov-Grinevich-Manakov algorithm, which is efficient in reconstructing two-dimensional acoustic refractive-absorbing scatterers of actually arbitrary shape and strength, are presented. The algorithm rigorously allows for multiple scattering effects. It is promising for tomography-like application problems and features a sufficiently high interference resistance.     

Finite difference computation of acoustic scattering by small surface inhomogeneities and discontinuities

The use of finite difference schemes to compute the scattering of acoustic waves by surfaces made up of different materials with sharp surface discontinuities at the joints would, invariably, result in the generations of spurious reflected waves of numerical origin. Spurious scattered waves are produced even if a high-order scheme capable of resolving and supporting the propagation of the incident wave is used. This problem is of practical importance in jet engine duct acoustic computation. In this work, the basic reason for the generation of spurious numerical waves is first examined. It is known that when the governing partial differential equations of acoustics are discretized, one should only use the long waves of the computational scheme to represent or simulate the physical waves. The short waves of the computational scheme have entirely different propagation characteristics. They are the spurious numerical waves. A method by which high wave number components (short waves) in the wave scattering process is intentionally removed so as to minimize the scattering of spurious numerical waves is proposed. This method is implemented in several examples from computational aeroacoustics to illustrate its effectiveness, accuracy and efficiency. This method is also employed to compute the scattering of acoustic waves by scatterers, such as rigid wall acoustic liner splices, with width smaller than the computational mesh size. Good results are obtained when comparing with computed results using much smaller mesh size. The method is further extended for applications to computations of acoustic wave reflection and scattering by very small surface inhomogeneities with simple geometries.     

The scanning acoustic microprobe: I. Analysis and synthesis of a spherically symmetric point spread function

The scanning acoustic microprobe is a novel system which probes and characterizes, from a limited acoustic window, the fine-scale structural features of an object, point-by-point, using a multiplicity of acoustic pulses all aimed and focused at that point. Spherically symmetric, three-dimensional Gaussian pulses are synthesized to measure the backscatter diffraction pattern of the least-resolvable volume of scatterers centered at the point in question. The size and distribution of the scattering volume is forced to be constant, independent of frequency and angle. This method is analytically simple, compared with other pulse-echo techniques, and is applicable to scatterers ranging continuously in size from Rayleigh scatterers to specular reflectors. This is the first of a number of papers describing the development and application of systems based on these concepts. The analytical principles will be described herein for examination of one point at a time. In a companion paper appearing in this issue [F. E. Barber, J. Acoust. Soc. Am. 90, 11-19 (1991)], application to measurement and characterization of a discrete, flat, circular "piston" will be presented. Application to human tissue imaging and tissue characterization will be described in a subsequent third paper. The primary features detected experimentally are the strength of nondirective patterns, and the strength, orientation, and directivity of angle-dependent echo functions associated with planar or layered structures. Fine-scale structural features of a scattering center are obtained either by pattern recognition in k(data) space or inverse Fourier transformation. It is shown that when the bandwidth criteria are met to produce a spherically symmetric point spread function, scattering phenomena are completely described by only two parameters, namely the center frequency of the pulse-echo system and the characteristic diameter of the Gaussian point spread function.     

丛聚的含气泡水对线性声传播的影响

含气泡水内气泡的空间分布会对线性声传播产生影响,导致实验结论与理论预测存在较大偏差.为解决这一问题,将准晶体近似引入到自洽方法中,导出了考虑空间分布时多分散含气泡水的等效声波波数.考虑到含气泡水内,气泡间存在小范围的聚集趋势(简称丛聚现象),在此基础上引入Neyman-Scott点过程描述了含气泡水内气泡的丛聚现象.分析发现,丛聚时,声速、声衰减的峰值将受到抑制,并向低频偏移,且抑制和频偏现象会随丛聚加剧而变强;随频率远离峰值段,丛聚对声传播的影响逐渐减弱.此外,考虑到空间分布的统计信息提取对相关研究的精确与否起到重要作用,引入了一种比例无偏估计,通过该方法获得了仿真环境下丛聚含气泡水模型的相速度及衰减系数,该建模及统计方法也可为相关实验工作提供理论基础.     

Solution of the three-dimensional problem of wave diffraction by an ensemble of objects

The pattern equations method is extended to solving three-dimensional problems of wave diffraction by an ensemble of bodies. The method is based on the reduction of the initial problem to a system of N (N is the number of scatterers in the ensemble) integro-operator equations of the second kind for the scattering patterns of scatterers. With the use of the series expansions of the scattering patterns in angular spherical harmonics, the problem is reduced to an algebraic system of equations in the expansion coefficients. An explicit (asymptotic) solution to the problems is obtained in the case when the scattering bodies are separated by sufficiently long distances. It is shown that the method can be used to model the characteristics of wave scattering by complex-shaped bodies.     

3D analysis of in-filled trench as passive barriers for ground vibration isolation

The three-dimensional (3D) problem of the ground vibration isolation by an in-filled trench as a passive barrier is studied theoretically. Integral equations governing Rayleigh wave scattering are derived based on the Green’s solution of Lamb problem. The integral equations are solved accurately and efficiently with an iteration technique. They are used to evaluate the complicated Rayleigh wave field generated by irregular scatterers embedded in an elastic half-space solid. The passive isolation effectiveness of ground vibration by the in-filled trench for screening Rayleigh wave is further studied in detail. Effects of relevant parameters on the effectiveness of vibration isolation are investigated and presented. The results show that a trench filled with stiff backfill material gets a better isolation effect than a soft one, and increasing the depth or width of the in-filled trench also improves its screening effectiveness. The effectiveness and the area of the screened zone are surging with the increase in the length of the in-filled trench. Supported by the National Natural Science Foundation of China (Grant Nos. 50678128 and 50538010) and the Research Fund for PhD Student of Chinese College (Grant No. 20050247030)     

On the numerical approximation of high-frequency acoustic multiple scattering problems by circular cylinders

The aim of this paper is to propose a numerical strategy for computing the solution of two-dimensional time-harmonic acoustic multiple scattering problems at high-frequency. The scatterers are assumed to be circular, leading therefore to semi-analytical representation formulae of the scattered field through the solution of a large linear system of equations. Taking advantage of the special block Toeplitz structure of the matrix of the linear system, a fast iterative and preconditioned numerical method yielding large memory savings is proposed. Several numerical experiments for general configurations are presented to show the efficiency of the numerical method.     

开口/封闭薄壳体声辐射和散射的统一边界积分方程解法

推导了在二维和三维空间下开口和封闭薄壳体在任意阻抗边界条件下声辐射和散射的统一边界积分方程.相对于以前的求解方法,该方程求解声辐射和散射问题具有相同的影响矩阵,能够同时求解薄壳体气动和振动噪声的辐射和散射现象,以及分析壳体声阻抗对声波传播的影响.推导的方程可以应用于叶轮机械、管道等噪声和消声器消声性能的预测等方面.在此方程基础上,可以进一步考虑运动边界和运动介质对声辐射和散射的影响. 关键词： 薄壳体 声阻抗 积分方程 边界元方法