利用海洋环境噪声互相关实现阵形校准

阵列信号处理需要准确的阵形,然而海底水平阵的阵形很难精确测量。利用海洋环境噪声可以提取两点间的格林函数,这为无源校准阵形提供了新的可能。通过海洋环境噪声互相关叠加可以获取两个阵元之间的距离信息,一种方法是可以通过该距离与位置之间的关系构建非线性方程组并利用最小二乘法(LS)求解阵形。本文利用多维尺度变换法(MDS)求解,通过对度量矩阵的特征值分解获得阵形。数值仿真和实验数据证明了该方法的可行性,并且MDS相比于LS有更稳健的性能。      

利用海浪噪声自相关实现散射体无源探测

提出了一种利用海浪噪声自相关实现散射体无源探测的新方法.将各接收器记录噪声信号的自相关减去所有接收器记录噪声信号自相关的平均值,得到散射信号的到达结构,然后结合基尔霍夫移位算法实现对散射体的探测.与利用背景噪声互相关提取格林函数从而实现散射体探测的方法不同,自相关无需考虑各个接收器之间的大量数据传输及时间同步问题,这为相距较远的多接收器和移动平台目标探测提供了极大的方便.将所提出的方法应用于实验数据中,最终探测结果与实际测量结果相比差别不大,验证了方法的有效性.     

交替投影盲反卷积的阵不变量被动定位*

本文针对低信噪比下传统盲反卷积技术提取格林函数过程中存在的图像模糊问题,提出了一种基于交替投影法的阵不变量浅海被动定位方法。该种方法利用交替投影的方法对波束信号进行分离,并利用分离出的波束分别提取格林函数。最后通过补偿时延的方式将多组波束提取的结果相干组合起来,实现格林函数的多波束联合估计。通过仿真实验证明了该方法有效增强了提取格林函数的强度,使得在低信噪比下利用阵不变量的连续声源定位方法更加稳健。     

扩散场重建格林函数检测钢轨近表面缺陷*

利用超声相控阵全矩阵捕获功能获取后期时间的扩散场信息,实现钢轨近表面缺陷成像。对扩散场信号进行互相关,重建阵元之间的格林函数,获取阵元之间未延时的响应,恢复被噪声湮没的早期缺陷信息,将波数成像方法用于近表面缺陷的快速成像。采用两种超声相控阵探头在钢轨表面采集实验数据,说明格林函数恢复效果与相控阵阵元的个数和激励频率有关。实验结果清晰地显示了距钢轨表面5～10 mm处的缺陷,证实了互相关方法重建格林函数对近表面缺陷检测的有效性。     

采用子带二值加权累积的经验格林函数提取

提出一种采用子带二值加权累积的海洋环境噪声互相关函数(NCF)提取经验格林函数(EGF)方法.首先将每一快拍NCF在频域划分为多个子带,每个子带内根据"累积后提取的EGF信噪比增加"的准则,确定各快拍NCF的加权系数为0或1,将各子带的加权累积结果谱白化后在频域拼接,再反傅里叶变换得到时域EGF.子带二值加权累积方法实...     

单边噪声源环境下的格林函数提取

对单边噪声源环境下空间中存在散射体时的格林函数到达时间结构提取展开了实验性研究。通过对位于沙滩上的传声器上记录的海浪噪声进行互相关处理,成功提取出两传声器间的格林函数到达时间结构。到达时间结构中观察到伪散射路径,并且伪散射路径在幅值上远高于散射路径。针对上述现象给出了理论及仿真解释:当传声器对受噪声源单边照射时,广义光学定理失效,伪散射路径将会出现;同时,由于伪散射路径由所有方向噪声源贡献,散射路径仅有稳相点附近声源贡献,因此伪散射路径在幅值上可能高于散射路径。结果表明,在单边噪声源环境下,基于互相关技术可以提取伪散射路径。伪散射路径的短时稳定提取可以为散射体无源成像提供更多信息,在海洋声学、地震学等领域具有潜在应用价值。      

利用矢量海洋环境噪声提取声场格林函数

考虑到矢量水听器在垂直方向上具有8字形指向性,能够有效抑制远方非平稳噪声源的干扰,提出了一种矢量环境噪声相关函数(NCF)提取声场时域格林函数(TDGF)的方法。基于简正波理论建立了声压和垂直振速垂直相关性模型。在此基础上,给出了声压和垂直振速相关函数提取声场纵向格林函数的过程.数值仿真对比和实验数据分析表明,相对于声压提取方法,垂直振速提取方法能够有效消除直达波前出现的亮纹与亮区干扰。此外,对于同等时间长度噪声序列,声压提取方法只提取到直达波路径,而垂直振速提取方法还提取到了我们更为关心的海底反射路径。利用直达波与海底反射波到达时延差估计的海深与实测海深吻合较好。      

基于子带二值加权累积的经验格林函数提取*

提出一种采用子带二值加权累积的海洋环境噪声互相关函数(NCF)提取经验格林函数(EGF)方法。首先将每一快拍NCF在频域划分为多个子带,每个子带内根据“累积后提取的EGF信噪比增加”的准则,确定各快拍NCF的加权系数为0或1,将各子带的加权累积结果谱白化后在频域拼接,再反傅里叶变换得到时域EGF。子带二值加权累积方法实现了每一快拍NCF累积前的频率自动选择,相比于已有的原始累积方法与时域二值加权累积算法,可以有效提高从较宽频率带宽的海洋环境噪声中提取EGF的信噪比。海试实验数据验证了该方法的有效性与优越性。     

利用环境噪声互相关实现散射体无源成像

最近研究表明利用环境噪声的互相关可以恢复两点之间的时域格林函数(声脉冲响应),这一原理在文献中被称为格林函数恢复。基于此原理,通过对多个传声器所接收的环境噪声进行互相关处理,获取与散射体相关联的散射波的到达时延信息,结合最小二乘反演算法和改进克希霍夫移位算法,分别获得道路交通噪声场中石柱以及海浪噪声场中塑料桶的空间位置,且其估计结果与实际测量相一致。实验结果表明将环境噪声作为探测信号进行散射体无源成像是可行的。这为设计室内无源声监测系统以及通过海洋环境噪声实现对水中静默目标成像提供了新思路和有益参考。      

海洋噪声场的时空相关函数

本文借助于方向性点源的平面波展开式导出了海洋中噪声场的时空相关函数的普遍表式。用演绎法将该表式用于分层不均匀浅海与分层不均匀深海,建立了噪声场的时空相关函数与海洋环境参数的有机联系。文中利用接收器附近点(噪声)源的奇性条件简化了浅海中噪声场时空相关函数的计算。接着用简正波理论计算了均匀浅海中噪声场的时空相关函数。     

The prediction of jet noise ground effects using an acoustic analogy and a tailored Green's function

An assessment of an acoustic analogy for the mixing noise component of jet noise in the presence of an infinite surface is presented. The reflection of jet noise by the ground changes the distribution of acoustic energy and is characterized by constructive and destructive interference patterns. The equivalent sources are modeled based on the two-point cross-correlation of the turbulent velocity fluctuations and a steady Reynolds-Averaged Navier–Stokes (RANS) solution. Propagation effects, due to reflection by the surface and refraction by the jet shear layer, are taken into account by calculating the vector Green's function of the linearized Euler equations (LEE). The vector Green's function of the LEE is written in relation to that of Lilley's equation; that is, it is approximated with matched asymptotic solutions and Green's function of the convective Helmholtz equation. The Green's function of the convective Helmholtz equation in the presence of an infinite flat plane with impedance is the Weyl–van der Pol equation. Predictions are compared with measurements from an unheated Mach 0.95 jet. Microphones are placed at various heights and distances from the nozzle exit in the peak jet noise direction above an acoustically hard and an asphalt surface. The predictions are shown to accurately capture jet noise ground effects that are characterized by constructive and destructive interference patterns in the mid- and far-field and capture overall trends in the near-field.     

Passive acoustic tomography of the ocean using arrays of unknown shape

It was recently established that ocean acoustic tomography based on an inversion of ray travel times can be implemented without use of any dedicated sound sources by cross-correlating the ambient noise recorded on two line arrays, the shapes of which are known. In contrast to active tomography, the amount of useful information from noise interferometry is proportional to the product of the numbers of receivers in the two arrays. In our study based on the 2D and 3D numerical experiments, we examine a hypothesis concerning the feasibility of simultaneous performance of a passive ray tomography and passive positioning of arrays through cross-correlation of ambient or shipping noise. The numerical experiments are conducted under conditions close to those of a field experiment on passive ocean tomography. It is demonstrated that, when using arrays of 20–40 hydrophones, the sound velocity profile and the array shape can be found from noise correlation to an accuracy adequate for oceanological and acoustic applications.     

Many‐body Green's function theory for thin ferromagnetic films: exact treatment of the single‐ion anisotropy

A theory for the magnetization of ferromagnetic films is formulated within the framework of many‐body Green's function theory which considers all components of the magnetization. The model Hamiltonian includes a Heisenberg term, an external magnetic field, a second‐ and fourth‐order uniaxial single‐ion anisotropy, and the magnetic dipole‐dipole coupling. The single‐ion anisotropy terms can be treated exactlyby introducing higher‐order Green's functions and subsequently taking advantage of relations between products of spin operators which leads to an automatic closure of the hierarchy of the equations of motion for the Green's functions with respect to the anisotropy terms. This is an improvement on the method of our previous work, which treated the corresponding terms only approximately by decoupling them at the level of the lowest‐order Green's functions. RPA‐like approximations are used to decouple the exchange interaction terms in both the low‐order and higher‐order Green's functions. As a first numerical example we apply the theory to a monolayer for spin S = 1 in order to demonstrate the superiority of the present treatment of the anisotropy terms over the previous approximate decouplings.     

Two-particle Green's functions in non-equilibrium matter

The method of the derivation of two-particle Green's functions in non-equilibrium matter is developed. The closed set of equations for the vertex functions and also for the two-particle Green's functions is obtained by means of the summation of the series of irreducible diagrams. The solution of such equations completely defines the two-particle Green's functions in matter.     

Exact Green's function for rectangular potentials and its application to quasi-bound states

In this work we calculate the exact Green's function for arbitrary rectangular potentials. Specifically we focus on Green's function for rectangular quantum wells enlarging the knowledge of exact solutions for Green's functions and also generalizing and resuming results in the literature. The exact formula has the form of a sum over paths and always can be cast into a closed analytic expression. From the poles and residues of the Green's function the bound states eigenenergies and eigenfunctions with the correct normalization constant are obtained. In order to show the versatility of the method, an application of the Green's function approach to extract information of quasi-bound states in rectangular barriers, where the standard analysis of quantum amplitudes fail, is presented.     

Visualization of aerodynamic noise source around a rotating fan blade

The purpose of this study is to understand the aerodynamic noise source distribution around a rotating fan blade by measuring the noise signal and velocity field around the blade. The local noise-level distribution over the fan blade is measured by microphone arrays, and the flow field is visualized by smoke and phase-averaged PIV measurement. The noise source distribution is examined by cross-correlation analysis between noise signal and velocity fluctuation. It is found that the noise source is located near the rotating fan blade, especially around leading and trailing edges. The separation and reattachment of flow are observed near the leading edge, and the tip vortices and vortex shedding are found near the trailing edge. The cross-correlation distribution of the noise signal and the radial velocity fluctuation shows large magnitude in the correlated regions, which indicates the noise generation by the formation of vortex structure around the blade.     

Electromagnetic volume scattering in the half-space: a moment-method analysis without Sommerfeld integrals or their approximations

Abstract

Traditionally, in moment-method analyses of electromagnetic scattering, the elements of the impedance matrix are calculated as convolutions of the basis elements with the appropriate dyadic Green's function. However, for scattering in the half-space, the vertical and azimuthal copolar terms of the Green's function require evaluation of Sommerfeld integrals which are computationally burdensome. In this paper, it is shown that, in populating the impedance matrix for the half-space problem, evaluation of Sommerfeld integrals is, in fact, not necessary. For monochromatic excitation, the plane-wave expansion of the scattered field constitutes a Fourier transform, in the horizontal plane, of a vector spectral function. This vector function results from the convolution, in the vertical dimension, of the respective angular spectra of the Green's function and the equivalent current. On application of the moment method, through the Weyl identity, the impedance-matrix elements corresponding to the singular terms of the Green's function are convolutions in the horizontal plane of spherical potentials, and Fourier transforms of scalar spectral functions. These scalar functions are derived from the basis elements and, with a judicious choice of basis, they are well behaved and of compact support, and consequently their Fourier transforms can be computed as FFTs.