Spherical angular spectrum and the fractional order Fourier transform

The notion of a spherical angular spectrum leads to the decomposition of the field amplitude on a spherical emitter into a sum of spherical waves that converge onto the Fourier sphere of the emitter. Unlike the usual angular spectrum, the spherical angular spectrum is propagated as the field amplitude, in a way that can be expressed by a fractional order Fourier transform.     

The far field asymptotics in the problem of diffraction of an acoustic plane wave by an impedance cone

The work deals with the far field asymptotics of the classical solution for the problem of diffraction by an impedance cone. The incident acoustic plane wave completely illuminates the semi-infinite conical surface. The scattered field contains different components in the asymptotics, namely, the spherical wave from the vertex of the cone, the reflected waves, and, under some conditions, also the surface waves of Rayleigh type. We give integral representations for the scattering diagram of the spherical wave. The uniform (with respect to the observation direction) asymptotic expression for the wave field is also addressed and described by the parabolic cylinder ansatz. Dedicated to the memory of Vladimir Borovikov     

Inversion of synthetic and experimental acoustical scattering data for the comparison of two reconstruction methods employing the Born approximation

This work is concerned with the reconstruction, from measured (synthetic or real) data, of a 2D penetrable fluid-like object of arbitrary cross-section embedded in a fluid of infinite extent and insonified by a plane acoustic wave. Green's theorem is used to provide a domain integral representation of the scattered field. The introduction therein of the Born approximation gives rise to a linearized form of the inverse problem. The actual inversion is carried out by two methods. The first diffraction tomography (DT), exhibits the contrast function very conveniently and explicitly in the form of a wave number/incident angle Fourier transform of the far backscattered field and thus requires measurements of this field for incident waves all around the object and at all frequencies. The second discretized domain integral equation with Born approximation method, is numerically more intensive, but enables a wider choice of configurations and requires less measurements (one or several frequencies, one or several incident waves, choice of measurement points) than the DT method. A comparison of the two methods is carried out by inversion of both simulated and experimental scattered field data.     

非均匀流体介质内部散射声场重建的逐层计算方法*

入射声波激励下非均匀流体介质内部散射声场的重建方法对超声层析成像具有重要意义。以往采用矩量法求解,但该方法全域离散形成的复数满秩矩阵规模随着分辨率与计算精度的提高而急剧增大,对算力具有很高的要求,一定程度上限制了其在实际中的应用。为克服上述缺陷,本文以逐层离散、逐层计算为核心思想,以声散射基本公式与近场声全息理论为基础,推导出逐层计算非均匀流体介质内部散射声场的理论公式并给出对应的几何离散模型。为验证该方法的可行性,以矩量法为参照,对同样的介质模型进行介质内部声场重构仿真。结果表明,逐层算法不仅可以有效地重建非均匀流体介质内部散射声场,且大幅度减小了求解规模。     

DIFFRACTION BY A HARD HALF-PLANE: USEFUL APPROXIMATIONS TO AN EXACT FORMULATION

In this paper, the problem of diffraction of a spherical wave by a hard half-plane is considered. The starting point is the Biot-Tolstoy theory of diffraction of a spherical wave by a fluid wedge with hard boundaries. In this theory, the field at a point in the fluid is composed eventually of a geometrical part: i.e., a direct component, one or two components due to the reflections on the sides of the hard wedge, and a diffracted component due exclusively to the presence of the edge of the wedge. The mathematical expression of this latter component has originally been given in an explicit closed form for the case of a unit momentum wave incidence, but Medwin has further developed its expression for the more useful case of a Dirac delta point excitation. The expression of this form is given in the time domain, but it is quite difficult to find exactly its Fourier transform for studying the frequency behaviour of the diffracted field. It is thus the aim of this paper to present various useful approximations of the exact expression. Among the approximations treated, three are most accurate for engineering purposes, and one of them is proposed, for its simplicity, as appropriate for most occurring practical situations.     

Reflexion of a spherical wave by the plane interface between a perfect fluid and a porous medium

By using the Fourier transform of the system of equations and continuity conditions, it is easy to obtain the Fourier transform of the solution. This last function is decomposed into several terms which are identified as Fourier images of known functions. An exact representation of the scattered sound pressure field is obtained as a combination of the radiation of the image source and layer potentials. Approximations are given when the point spherical source is located on the interface.     

Regularization methods for near-field acoustical holography.

The reconstruction of the pressure and normal surface velocity provided by near-field acoustical holography (NAH) from pressure measurements made near a vibrating structure is a linear, ill-posed inverse problem due to the existence of strongly decaying, evanescentlike waves. Regularization provides a technique of overcoming the ill-posedness and generates a solution to the linear problem in an automated way. We present four robust methods for regularization; the standard Tikhonov procedure along with a novel improved version, Landweber iteration, and the conjugate gradient approach. Each of these approaches can be applied to all forms of interior or exterior NAH problems; planar, cylindrical, spherical, and conformal. We also study two parameter selection procedures, the Morozov discrepancy principle and the generalized cross validation, which are crucial to any regularization theory. In particular, we concentrate here on planar and cylindrical holography. These forms of NAH which rely on the discrete Fourier transform are important due to their popularity and to their tremendous computational speed. In order to use regularization theory for the separable geometry problems we reformulate the equations of planar, cylindrical, and spherical NAH into an eigenvalue problem. The resulting eigenvalues and eigenvectors couple easily to regularization theory, which can be incorporated into the NAH software with little sacrifice in computational speed. The resulting complete automation of the NAH algorithm for both separable and nonseparable geometries overcomes the last significant hurdle for NAH.     

Fast prediction of scattered sound field based on Fourier diffraction theory under second-order born approximation

The directional pattern of sound waves scattered from an object insonified by a plane wave can be efficiently predicted using the Fourier diffraction theorem (FDT). This is achieved by sampling a circle in the discrete Fourier transform of the object/medium distribution. However, the FDT-based approach under the first-order Born approximation is only applicable to weak scattering. To improve the prediction accuracy and expand the method’s scope of applications, we introduce a second-order correction term to the solution, which is obtained by taking the first-order scattered waves as secondary incident sources, and calculate the “scattering” in the same way as in the first-order FDT-based approach. Adding the resulting correction term to the directional pattern based on the first-order Born approximation, the second-order prediction is obtained. Numerical results show that the proposed method can provide improved directional patterns of the scattered waves, and the range of applicability is significantly expanded.     

一种基于非衍射波的高帧率超声成像发射系统的研究

基于非衍射波的高帧率(high frame rate,HFR)超声成像的前提是发射阵列波(array beam)声场,但由于阵列波的正弦函数特征使对超声发射信号加权处理变得比较困难,导致了发射系统复杂.针对这一问题,本文提出一种有效的解决方案.基于傅里叶变换理论,从方波的能量成分主要集中在基波这一特性出发,采用单值方波对超声传感器发射信号进行加权,用一个功率信号源和简单的电子开关网络实现了HFR系统中所需要的发射电路.实验表明,这种方案对HFR的成像质量几乎没有什么影响,但大大简化了HFR成像系统,这为HFR超声成像方法在现有B超声系统上的实现提供了一个可行的方法.     

Scattering of waves by irregularities in periodic discrete lattice spaces. I. Reduction of problem to quadratures on a discrete model of the Schrödinger equation

The scattering of plane waves and of point source pulses by irregularities in a discrete lattice model of the Schrödinger equation is considered. Closed form expressions are derived for the scattered wave function in terms of lattice Green's functions in the case that a finite number of lattice points or “bonds” are defective. The scattered wave function appears in the form of the ratio of two determinants. While in continuum scattering theory the scatterer must have some symmetry, perhaps spherical, cylindrical or elliptical, in order to allow separation of variables in the basic scattering differential equation, such symmetries are not necessary for the construction of scattered wave functions on discrete lattices. When the number of irregularities becomes large, the determinants in the solution of the scattering problem become large.     

3D-solitons in a strongly magnetized plasma: Transition from flat to spherical waves of the Zakharov-Kuznetsov equation

The Zakharov-Kuznetsov equation has been used to describe ion-acoustic wave propagation in a strong magnetic plasma. An initial-value problem has been solved for this equation on the basis of a numerical method that uses the fast Fourier transform technique for calculating space derivatives and a fourth order Runge-Kutta method for the time scheme. Numerical simulations have shown that the disturbed flat solitary waves can break up into spherical ones.     

Surface wave scattering by a nano-object situated on a surface

The scattering of the surface electromagnetic waves by a nano-defect (object) on a surface was calculated. The scattered field has been considered as a field caused by the current generated by the self-consistent local field inside the defect. In turn, the self-consistent local field has been determined as a result of solution of the integral Lippmann-Schwinger equation. The effective susceptibility of the object has been calculated using a self-consistent procedure. The corrections of self-energy part due to direct and indirect electromagnetic interactions, as well as due to interaction with surface wave field are taken into account. The self-energy part is calculated analytically within the framework of the near-field approximation. The scattering indicatrisses in reciprocal space have been computed for different shapes of the scatterer. Strong dependence of the scattered field on geometry of the scatterer has been found and explained.     

The role of surface tension in elastic wave scattering in an inhomogeneous poroelastic medium

It is well known that many porous media such as rocks have heterogeneities at nearly all scales. We applied Biot's poroelastic theory to study the propagation of elastic waves in isotropic porous matrix with spherical inclusions. It is assumed that the heterogeneity dimension exceeds significantly the pore size. Modified boundary conditions on poroelastic interface are used to take into account the surface tension effects. The effective wavenumber is calculated using the Waterman and Truell multiple scattering theory, which relates the effective wave number to the amplitude of the wave field scattered by a single inclusion. The calculations were performed for a medium containing fluid-filled cavities or porous inclusions contrasting in saturating fluid elastic properties. The results obtained show that when we consider elastic wave propagation in poroelastic medium containing soft inclusions, it is necessary to take into account the capillary pressure. The influence of the surface tension depends on the diffraction parameter and it is a maximum in the low frequency range.     

The role of surface tension in elastic wave scattering in an inhomogeneous poroelastic medium

It is well known that many porous media such as rocks have heterogeneities at nearly all scales. We applied Biot's poroelastic theory to study the propagation of elastic waves in isotropic porous matrix with spherical inclusions. It is assumed that the heterogeneity dimension exceeds significantly the pore size. Modified boundary conditions on poroelastic interface are used to take into account the surface tension effects. The effective wavenumber is calculated using the Waterman and Truell multiple scattering theory, which relates the effective wave number to the amplitude of the wave field scattered by a single inclusion. The calculations were performed for a medium containing fluid-filled cavities or porous inclusions contrasting in saturating fluid elastic properties. The results obtained show that when we consider elastic wave propagation in poroelastic medium containing soft inclusions, it is necessary to take into account the capillary pressure. The influence of the surface tension depends on the diffraction parameter and it is a maximum in the low frequency range.     

An efficient method for computing backscattering from Born objects of arbitrary shape

A method is presented for efficiently computing the propagating pressure field backscattered by an arbitrarily shaped, weakly scattering, three-dimensional object. This is accomplished by drawing upon a previously reported relationship between the boundary condition on a two-dimensional radiating aperture and the pressure propagating along an axis normal to the aperture, and the fundamental theorem of diffraction tomography, which relates the Fourier transform of an object function to its scattered pressure field. Together, these two results are used to derive an integral formula that expresses the pressure field backscattered from an object as a one-dimensional Fourier transform of its scattering amplitude. This formula is then utilized to compute the backscattered pressure field from a uniform fluid sphere in the first Born approximation; the results of which are compared to the rigorous partial wave expansion.     

非傍轴高斯光束传输方程解的探讨

由洛伦兹规范下光场满足的矢势和标势方程出发, 探讨了初始光场具有任意对称性的一般非傍轴光束在空间的传输情况。利用角频谱并结合傅里叶变换处理得到了非傍轴光束的传输方程。利用扰动方法并结合新近似取法得到了非傍轴横向光场的扰动形式解, 该解对于近场弱非傍轴情况下的效果比强非傍轴条件下要好。给出了光束非傍轴传输高阶(八阶)修正解。对于零阶情况, 非傍轴光束退化为傍轴光束, 传输方程也相应退化为傍轴传输的模型方程。高斯解为扰动参量的零阶修正解。     

Evanescent Waves and Energy Conservation

ＥｖａｎｅｓｃｅｎｔＷａｖｅｓａｎｄＥｎｅｒｇｙＣｏｎｓｅｒｖａｔｉｏｎＣＡＯＱｉｎｇ；ＤＥＮＧＸｉｍｉｎｇ（ＮａｔｉｏｎａｌＬａｂｏｒａｔｏｒｙｏｎＨｉｇｈＰｏｗｅｒＬａｓｅｒａｎｄＰｈｙｓｉｃｓ，ＳｈａｎｇｈａｉＩｎｓｔｉｔｕｔｅｏｆＯｐｔｉｃｓ...     

Regularization for improving the deconvolution in real-time near-field acoustic holography

Near-field acoustic holography is a measuring process for locating and characterizing stationary sound sources from measurements made by a microphone array in the near-field of the acoustic source plane. A technique called real-time near-field acoustic holography (RT-NAH) has been introduced to extend this method in the case of nonstationary sources. This technique is based on a formulation which describes the propagation of time-dependent sound pressure signals on a forward plane using a convolution product with an impulse response in the time-wavenumber domain. Thus the backward propagation of the pressure field is obtained by deconvolution. Taking the evanescent waves into account in RT-NAH improves the spatial resolution of the solution but makes the deconvolution problem "ill-posed" and often yields inappropriate solutions. The purpose of this paper is to focus on solving this deconvolution problem. Two deconvolution methods are compared: one uses a singular value decomposition and a standard Tikhonov regularization and the other one is based on optimum Wiener filtering. A simulation involving monopoles driven by nonstationary signals demonstrates, by means of objective indicators, the accuracy of the time-dependent reconstructed sound field. The results highlight the advantage of using regularization and particularly in the presence of measurement noise.     

Spherical Wave Diffraction by a Slit in an Impedance Screen

An exact solution is developed for the problem of diffraction of an E-polarized and an H-polarized spherical waves by a slit in an impedance screen. This consideration is important in the sense that point sources are regarded as better substitutes for real sources than line sources/plane waves. The two independent problems are solved using the Fourier transform, the Wiener-Hopf technique and asymptotic approximations.     

A new formalism for time-dependent wave scattering from a bounded obstacle

A time-dependent three-dimensional acoustic scattering problem is considered. An incoming wave packet is scattered by a bounded, simply connected obstacle with locally Lipschitz boundary. The obstacle is assumed to have a constant boundary acoustic impedance. The limit cases of acoustically soft and acoustically hard obstacles are considered. The scattered acoustic field is the solution of an exterior problem for the wave equation. A new numerical method to compute the scattered acoustic field is proposed. This numerical method obtains the time-dependent scattered field as a superposition of time-harmonic acoustic waves and computes the time-harmonic acoustic waves by a new "operator expansion method." That is, the time-harmonic acoustic waves are solutions of an exterior boundary value problem for the Helmholtz equation. The method used to compute the time-harmonic waves improves on the method proposed by Misici, Pacelli, and Zirilli [J. Acoust. Soc. Am. 103, 106-113 (1998)] and is based on a "perturbative series" of the type of the one proposed in the operator expansion method by Milder [J. Acoust. Soc. Am. 89, 529-541 (1991)]. Computationally, the method is highly parallelizable with respect to time and space variables. Some numerical experiments on test problems obtained with a parallel implementation of the numerical method proposed are shown and discussed from the numerical and the physical point of view. The website: http://www.econ.unian.it/recchioni/w1 shows four animations relative to the numerical experiments.