声速不均匀介质的光声成像重建算法

为提高声速不均匀介质的光声成像精度,提出一种基于反卷积的光声成像重建算法.本算法不需预先知道介质的声速分布.先通过探测到的光声信号构造出一个新函数,并用不同位置探测到的光声信号间的相关性来估计空间两点间的声波传播时间,以补偿声速的不均匀性.然后基于反卷积方法由该函数解出待测组织内的电磁波吸收分布.仿真研究结果表明,当介质声速差异在10%以内时,重建图像能正确反映待测目标的大小、位置和电磁波吸收系数,算法具有良好的抗噪性能.由于生物软组织内的声速差异一般小于10%,因此本算法是一种有效的光声成像重建算法.     

人体组织相位畸变模型仿真研究

为研究医学超声成像中的相位畸变模型,通过Westervelt方程的时域有限差分方法,仿真比较了近场相位屏模型、全局声速误差模型以及分层介质模型在聚焦发射过程和回波信号接收过程的特点,并参考理想均匀介质与人体非均匀介质的结果,发现几种相位畸变模型都难以准确模拟人体腹壁对发射声场的影响,同时近场相位屏模型相对其他两种模型虽然能够较好地描述出非均匀组织的相位畸变特点,但无法表现出全局声速误差模型及分层介质模型中相位误差跟随目标方位变化的关系。     

深海上层声速剖面水平非均匀分布对声场会聚区特性预报的影响

针对深海上层声速剖面水平方向弱变化对声场会聚区特征的影响,通过对0°出射声线传播轨迹的分析,提出了只考虑声源及会聚区位置附近声速剖面水平变化实现深海声场会聚区特征准确预报的方法,该方法有助于降低深海声传播调查时声速剖面水平非均匀分布的测量工作量。通过模拟小振幅内波引起深海上层声速剖面水平非均匀变化,对该方法进行了仿真验证。南海声传播试验数据分析表明,当不考虑上层声速剖面的水平非均匀分布时,模型预报的声场会聚区位置及会聚区能量分布与实测结果有一定差异,该差异具有距离累积效应,到第二会聚区时位置累计误差可达3 km,当考虑会聚区宽度内声速剖面水平变化时,模型预报的声场会聚区位置与实测结果非常吻合,第二会聚区位置预报误差小于1 km。     

非均匀声速场中高精度逐次逼近声源定位方法

针对非均匀声速场中传统飞行时间法无法直接计算声源到各声接收点的精确距离从而影响定位精度的问题,提出了一种基于逐次逼近方法的高精度声源定位方法。建立二维平面模型,推导证明了均匀声速场和固定声接收阵型中逐次逼近定位理论的正确性。然后,在任意分布的非均匀声速场和接收阵型条件下,仿真证明逐次逼近法的有效性,以及迭代次数和接收点个数与定位精度的关系。结果表明,该算法可用于非均匀声速场中任意位置声源的声学定位,在125 m×125 m水平区域中,声速起伏范围为±5 m/s,当迭代次数大于26次时,定位误差可小于1 mm,且定位精度随着迭代次数和接收点个数的增加而提高。     

含时Schrdinger方程的高阶辛FDTD算法研究

提出了一种新的算法—高阶辛时域有限差分法(SFDTD(3,4):symplectic finite-difference time-domain)求解含时薛定谔方程.在时间上采用三阶辛积分格式离散,空间上采用四阶精度的同位差分格式离散,建立了求解含时薛定谔方程的高阶离散辛框架;探讨了高阶辛算法的稳定性及数值色散性.通过理论上的分析及数值算例表明:当空间采用高阶同位差分格式时,辛积分可提高算法的稳定度;SFDTD(3,4)法和FDTD(2,4)法较传统的FDTD(2,2)法数值色散性明显改善.对二维量子阱和谐振子的仿真结果表明:SFDTD(3,4)法较传统的FDTD(2,2)法及高阶FDTD(2,4)法有着更好的计算精度和收敛性,且SFDTD(3,4)法能够保持量子系统的能量守恒,适用于长时间仿真.     

含时Schroedinger方程的高阶辛FDTD算法研究

提出了一种新的算法一高阶辛时域有限差分法（SFDTD（3,4）：symplectic finite—difference time-domain）求解含时薛定谔方程．在时间上采用三阶辛积分格式离散,空间上采用四阶精度的同位差分格式离散,建立了求解含时薛定谔方程的高阶离散辛框架;探讨了高阶辛算法的稳定性及数值色散性．通过理论上的分析及数值算例表明：当空间采用高阶同位差分格式时,辛积分可提高算法的稳定度;SFDTD（3,4）法和FDTD（2,4）法较传统的FDTD（2,2）法数值色散性明显改善．对二维量子阱和谐振子的仿真结果表明：SFDTD（3,4）法较传统的FDTD（2,2）法及高阶FDTD（2,4）法有着更好的计算精度和收敛性,且SFDTD（3,4）法能够保持量子系统的能量守恒,适用于长时间仿真．     

非均匀组织医学超声发射声场仿真

为了分析医学超声在非均匀组织中的分布特性,建立了超声发射声场的计算模型。采用交错网格有限差分法对耦合超声非线性方程进行数值求解,获得了声速、密度及非线性参数非均匀分布情况下人体腹壁组织内的超声声场分布数据。同均匀介质相比:当声速均匀而密度非均匀时,声束仍聚焦良好,焦点处声能下降了1.8 dB;当密度均匀而声速非均匀时,声束发散严重,焦点处声能下降了3.8 dB,下降程度与非均匀组织接近。组织声速在空间分布的非均匀性是导致聚焦声束能量分布畸变的主要原因。     

含时Schrödinger方程的高阶辛FDTD算法研究

提出了一种新的算法——高阶辛时域有限差分法(SFDTD(3, 4): symplectic finite-difference time-domain)求解含时薛定谔方程.在时间上采用三阶辛积分格式离散, 空间上采用四阶精度的同位差分格式离散, 建立了求解含时薛定谔方程的高阶离散辛框架;探讨了高阶辛算法的稳定性及数值色散性.通过理论上的分析及数值算例表明:当空间采用高阶同位差分格式时, 辛积分可提高算法的稳定度;SFDTD(3, 4)法和FDTD(2, 4)法较传统的FDTD(2, 2)法数值色散性明显改善.对二维量子阱和谐振子的仿真结果表明: SFDTD(3, 4)法较传统的FDTD(2, 2)法及高阶FDTD(2, 4)法有着更好的计算精度和收敛性, 且SFDTD(3, 4)法能够保持量子系统的能量守恒, 适用于长时间仿真.     

海洋声学中三维抛物方程非均匀网格模型

在海洋声学中,三维抛物方程模型可以有效考虑三维空间的声传播效应。然而,采用三维抛物方程模型分析三维空间内的声传播问题时,计算时间较长,并且需要消耗较大的计算机内存,因此给远距离声场的快速精确计算带来了很大困难。为此,将非均匀网格Galerkin离散化方法用于三维直角坐标系下的水声抛物方程模型中,深度算子和水平算子Galerkin离散方式由均匀网格变为非均匀网格。仿真结果表明,三维直角坐标系下非均匀网格离散的抛物方程模型,在保持计算精度、提高计算速度的同时,可以实现远距离声场的快速预报。另外,针对远距离局部海底地形与距离有关的三维声传播问题,给出了声场快速计算方法;在海底保持水平的区域,采用经典Kraken模型,重构抛物方程算法的初始场,随后依次递推求解地形与距离有关海底下的三维声场。采用改进模型,证明了远距离楔形波导声强增强效应。     

一种计算分层半空间上方导线瞬态响应的新方法

提出了一种将时域积分方程(time domain integral equation, TDIE)方法和时域有限差分(finite differnce time domain, FDTD)方法相结合计算分层有耗半空间上方导线瞬态电磁响应的新方法.其中,一维FDTD方法用于计算入射电磁波经分层半空间反射的时域波形.TDIE用于求解细导线在加入两个激励源(直接入射电磁波和经分层半空间反射的电磁波)时的瞬态响应.相关计算理论和数值模拟结果说明了本文方法是一种解决了分层有耗介质上方水平放置导线瞬态响应的高效解决方案. 关键词： 时域积分方程 时域有限差分 细导线 分层半空间     

Sound field of thermoacoustic tomography based on a modified finite-difference time-domain method

A modified finite-difference time-domain （FDTD） method is proposed for the sound field simulation of the thermoacoustic tomography （TAT） in the acoustic speed inhomogeneous medium. First, the basic equations of the TAT are discretized to difference ones by the FDTD. Then the electromagnetic pulse, the excitation source of the TAT, is modified twice to eliminate the error introduced by high frequency electromagnetic waves. Computer simulations are carried out to validate this method. It is shown that the FDTD method has a better accuracy than the commonly used time-of-flight （TOF） method in the TAT with the inhomogeneous acoustic speed. The error of the FDTD is ten times smaller than that of the TOF in the simulation for the acoustic speed difference larger than 50%. So this FDTD method is an efficient one for the sound field simulation of the TAT and can provide the theoretical basis for the study of reconstruction algorithms of the TAT in the acoustic heterogeneous medium.     

Equations for finite-difference, time-domain simulation of sound propagation in moving inhomogeneous media and numerical implementation

Finite-difference, time-domain (FDTD) calculations are typically performed with partial differential equations that are first order in time. Equation sets appropriate for FDTD calculations in a moving inhomogeneous medium (with an emphasis on the atmosphere) are derived and discussed in this paper. Two candidate equation sets, both derived from linearized equations of fluid dynamics, are proposed. The first, which contains three coupled equations for the sound pressure, vector acoustic velocity, and acoustic density, is obtained without any approximations. The second, which contains two coupled equations for the sound pressure and vector acoustic velocity, is derived by ignoring terms proportional to the divergence of the medium velocity and the gradient of the ambient pressure. It is shown that the second set has the same or a wider range of applicability than equations for the sound pressure that have been previously used for analytical and numerical studies of sound propagation in a moving atmosphere. Practical FDTD implementation of the second set of equations is discussed. Results show good agreement with theoretical predictions of the sound pressure due to a point monochromatic source in a uniform, high Mach number flow and with Fast Field Program calculations of sound propagation in a stratified moving atmosphere.     

声诱导电磁场的赫兹矢量表示与多极声电测井模拟

在假设声场不受电磁场影响的前提下，将Pride声电耦合方程组化为具有电流源的麦克斯韦方程组.与空间位置固定的电流源产生的电磁场不同，孔隙地层中声波诱导的电磁场是由空间波动的电流源产生的.通过引入赫兹矢量，将求解麦克斯韦方程组问题转化为求解关于赫兹矢量的非齐次矢量赫姆霍兹方程组.通过求解该方程组，得出电磁场表达式.利用此方法，针对声电效应测井，分别计算了由单极声源、偶极声源、四极声源激发的井内声场及其诱导电磁场的全波波形. 关键词： 孔隙介质 诱导电磁场 测井 多极声源     

Fast volumetric integral-equation solver for acoustic wave propagation through inhomogeneous media

Elements are described of a volumetric integral-equation-based algorithm applicable to accurate large-scale simulations of scattering and propagation of sound waves through inhomogeneous media. The considered algorithm makes possible simulations involving realistic geometries characterized by highly subwavelength details, large density contrasts, and described in terms of several million unknowns. The algorithm achieves its competitive performance, characterized by O(N log N) solution complexity and O(N) memory requirements, where N is the number of unknowns, through a fast and nonlossy fast Fourier transform based matrix compression technique, the adaptive integral method, previously developed for solving large-scale electromagnetic problems. Because of its ability of handling large problems with complex geometries, the developed solver may constitute an efficient and high fidelity numerical simulation tool for calculating acoustic field distributions in anatomically realistic models, e.g., in investigating acoustic energy transfer to the inner ear via nonairborne pathways in the human head. Examples of calculations of acoustic field distribution in a human head, which require solving linear systems of equations involving several million unknowns, are presented.     

Distributed sound source localization algorithm with sound velocity calibration in windy environments

A new sound source localization method with sound speed compensation is proposed to reduce the wind influence on the performance of conventional TDOA(Time Difference of Arrival) algorithms. First, the sound speed is described as a set of functions of the unknown source location, to approximate the acoustic velocity field distribution in the wind field. Then,they are introduced into the TDOA algorithm, to construct nonlinear equations. Finally, the particle swarm optimization algorithm is used to estimate the source location. The simulation results show that the proposed algorithm can significantly improve the localization accuracy for different wind velocities, source locations and test area sizes. The experimental results show that the proposed method can reduce localization errors to about 40% of the original error in a four nodes localization system.     

非均匀海洋环境噪声中的匹配场处理舰船辐射噪声级估计方法

在实际舰船辐射噪声测量过程中,受非均匀海洋环境噪声的影响,导致常规测量方法的性能急剧下降,基于此,提出了一种非均匀海洋环境噪声背景中的垂直阵舰船辐射噪声测量方法。根据水声信道传播理论,建立了由海面多个空时独立均匀分布噪声源构成的非均匀背景噪声场模型,推导了非均匀海洋环境噪声场中垂直阵舰船辐射噪声估计的理论公式。针对典型的浅海水声信道,进行计算机仿真实验,分析了该方法的测量性能并与常规匹配场测量方法进行对比,结果表明:(1)该方法能有效克服非均匀海洋环境噪声对测量结果的影响,测量误差较小;(2)相同测量条件下,该方法测量性能优于常规匹配场舰船辐射噪声级测量方法;(3)当信噪比满足一定要求时,测量得到的声源级与实际声源级相比,误差小于1 dB。     

FDTD simulation of finite-amplitude pressure and temperature fields for biomedical ultrasound.

Full wave simulations provide a valuable tool for studying the spatial and temporal nature of an acoustic field. One method for producing such simulations is the finite-difference time-domain (FDTD) method. This method uses discrete differences to approximate derivatives in the governing partial differential equations. We used the FDTD method to model the propagation of finite-amplitude sound in a homogeneous thermoviscous fluid. The calculated acoustic pressure field was then used to compute the transient temperature rise in the fluid; the heating results from absorption of acoustic energy by the fluid. As an example, the transient temperature field was calculated in biological tissue in response to a pulse of focused ultrasound. Enhanced heating of the tissue from finite-amplitude effects was observed. The excess heating was attributed to the nonlinear generation of higher-frequency harmonics which are absorbed more readily than the fundamental. The effect of nonlinear distortion on temperature rise in tissue was observed to range from negligible at 1 MPa source pressure to an 80% increase in temperature elevation at 10 MPa source pressure.     

风场环境中声速修正的分布式声源定位算法

为减小声速误差对定位精度的影响,提出了一种基于声速修正的分布式声源定位方法。首先,将声速表示为未知声源位置的函数,逼近风场中的声速场分布,然后将其代入TDOA (Time Differences of Arrival)算法中,构建非线性超定方程组,最后采用粒子群优化算法求解声源位置。对不同风速、不同声源位置及不同测试区域进行仿真,结果表明:修正后的定位精度比修正前有明显提高,尤其对于大范围并且声源靠近测试区域边缘位置的定位系统,改善更加明显;4个节点的定位系统实验结果表明,修正后的定位误差可降至修正前的4l%,该方法能更好的应用于风场中的定位系统。     

有耗介质层上多导体传输线的电磁耦合时域分析方法

目前,针对空间电磁场作用有耗介质层上传输线的电磁耦合,仍缺乏有效的数值分析方法.因此,本文提出一种高效的时域混合算法,很好地解决了有耗介质层上传输线电磁耦合建模难的问题.首先,对经典传输线方程进行改进,推导了适用于有耗介质层上多导体传输线电磁耦合分析的修正传输线方程.然后,结合时域有限差分方法和相应插值技术,求解修正传输线方程,获得多导线及其端接负载上的电压和电流响应,并实现空间电磁场辐射与多导线瞬态响应的同步计算.最后,通过相应计算实例的数值模拟,与CST软件的仿真结果进行对比,验证了时域混合算法的正确性和高效性.     

High-order FDTD methods for transverse electromagnetic systems in dispersive inhomogeneous media

This Letter introduces a novel finite-difference time-domain (FDTD) formulation for solving transverse electromagnetic systems in dispersive media. Based on the auxiliary differential equation approach, the Debye dispersion model is coupled with Maxwell's equations to derive a supplementary ordinary differential equation for describing the regularity changes in electromagnetic fields at the dispersive interface. The resulting time-dependent jump conditions are rigorously enforced in the FDTD discretization by means of the matched interface and boundary scheme. High-order convergences are numerically achieved for the first time in the literature in the FDTD simulations of dispersive inhomogeneous media.