Transient propagation in media with classical or power-law loss

This paper addresses the problem of small-signal transient wave propagation in media whose absorption coefficient obeys power-law frequency dependence, i.e., alpha infinity omega n. Our approach makes use of previously derived relations between the absorption and dispersion based on the Kramers-Kronig relations. This, combined with a recently obtained solution to a causal convolution wave equation enable expressions to be obtained for one-dimensional transient propagation when n is in the range 0 < n < 3. For n = 2, corresponding to no dispersion, straightforward analytical expressions are obtained for a delta-function and a sinusoidal step function sources and these are shown to correspond to relations previously derived. For other values of n, the effects of dispersion are accounted for by using Fourier transforms. Examples are used to illustrate the results for normal and anomalous dispersive media and to examine the question as to the conditions under which the effects of dispersion should be accounted for, especially for wideband ultrasound pulses of the type used in B-mode tissue imaging. It is shown that the product of the attenuation and total propagation path can be used as a criterion for judging whether dispersion needs to be accounted for.     

Stationary and quasi-stationary waves in even-order dissipative systems

A new equation was recently suggested by Rudenko and Robsman [1] for describing the nonlinear wave propagation in scattering media that are characterized by weak sound signal attenuation proportional to the fourth power of frequency. General self-similar properties of the solutions to this equation were studied. It was shown that stationary solutions to this equation in the form of a shock wave exhibit unusual oscillations around the shock front, as distinct from the classical Burgers equation. Here, similar solutions are studied in detail for nonlinear waves in even-order dissipative media; namely, the solutions are compared for the media with absorption proportional to the second, fourth, and sixth powers of frequency. Based on the numerical results and the self-similar properties of the solutions, the fine structure of the shock front of stationary waves is studied for different absorption laws and magnitudes. It is shown that the amplitude and number of oscillations appearing in the stationary wave profile increase with increasing power of the frequency-dependent absorption term. For initial disturbances in the form of a harmonic wave and a pulse, quasi-stationary solutions are obtained at the stage of fully developed discontinuities and the evolution of the profile and width of the shock wave front is studied. It is shown that the smoothening of the shock front in the course of wave propagation is more pronounced when the absorption law is quadratic in frequency.     

Modeling power law absorption and dispersion for acoustic propagation using the fractional Laplacian

The efficient simulation of wave propagation through lossy media in which the absorption follows a frequency power law has many important applications in biomedical ultrasonics. Previous wave equations which use time-domain fractional operators require the storage of the complete pressure field at previous time steps (such operators are convolution based). This makes them unsuitable for many three-dimensional problems of interest. Here, a wave equation that utilizes two lossy derivative operators based on the fractional Laplacian is derived. These operators account separately for the required power law absorption and dispersion and can be efficiently incorporated into Fourier based pseudospectral and k-space methods without the increase in memory required by their time-domain fractional counterparts. A framework for encoding the developed wave equation using three coupled first-order constitutive equations is discussed, and the model is demonstrated through several one-, two-, and three-dimensional simulations.     

Full wave modeling of therapeutic ultrasound: efficient time-domain implementation of the frequency power-law attenuation

For the simulation of therapeutic ultrasound applications, a method including frequency-dependent attenuation effects directly in the time domain is highly desirable. This paper describes an efficient numerical time-domain implementation of the power-law attenuation model presented by Szabo [Szabo, J. Acoust. Soc. Am. 96, 491-500 (1994)]. Simulations of therapeutic ultrasound applications are feasible in conjunction with a previously presented finite differences time-domain (FDTD) algorithm for nonlinear ultrasound propagation [Ginter et al., J. Acoust. Soc. Am. 111, 2049-2059 (2002)]. Szabo implemented the empirical frequency power-law attenuation using a causal convolutional operator directly in the time-domain equation. Though a variety of time-domain models has been published in recent years, no efficient numerical implementation has been presented so far for frequency power-law attenuation models. Solving a convolutional integral with standard time-domain techniques requires enormous computational effort and therefore often limits the application of such models to 1D problems. In contrast, the presented method is based on a recursive algorithm and requires only three time levels and a few auxiliary data to approximate the convolutional integral with high accuracy. The simulation results are validated by comparison with analytical solutions and measurements.     

Including dispersion and attenuation directly in the time domain for wave propagation in isotropic media

When sound propagates in a lossy fluid, causality dictates that in most cases the presence of attenuation is accompanied by dispersion. The ability to incorporate attenuation and its causal companion, dispersion, directly in the time domain has received little attention. Szabo [J. Acoust. Soc. Am. 96, 491-500 (1994)] showed that attenuation and dispersion in a linear medium can be accounted for in the linear wave equation by the inclusion of a causal convolutional propagation operator that includes both phenomena. Szabo's work was restricted to media with a power-law attenuation. Waters et al. [J. Acoust. Soc. Am. 108, 2114-2119 (2000)] showed that Szabo's approach could be used in a broader class of media. Direct application of Szabo's formalism is still lacking. To evaluate the concept of the causal convolutional propagation operator as introduced by Szabo, the operator is applied to pulse propagation in an isotropic lossy medium directly in the time domain. The generalized linear wave equation containing the operator is solved via a finite-difference-time-domain scheme. Two functional forms for the attenuation often encountered in acoustics are examined. It is shown that the presence of the operator correctly incorporates both, attenuation and dispersion.     

Fractional Laplacian time-space models for linear and nonlinear lossy media exhibiting arbitrary frequency power-law dependency

Frequency-dependent attenuation typically obeys an empirical power law with an exponent ranging from 0 to 2. The standard time-domain partial differential equation models can describe merely two extreme cases of frequency-independent and frequency-squared dependent attenuations. The otherwise nonzero and nonsquare frequency dependency occurring in many cases of practical interest is thus often called the anomalous attenuation. In this study, a linear integro-differential equation wave model was developed for the anomalous attenuation by using the space-fractional Laplacian operation, and the strategy is then extended to the nonlinear Burgers equation. A new definition of the fractional Laplacian is also introduced which naturally includes the boundary conditions and has inherent regularization to ease the hypersingularity in the conventional fractional Laplacian. Under the Szabo's smallness approximation, where attenuation is assumed to be much smaller than the wave number, the linear model is found consistent with arbitrary frequency power-law dependency.     

毫米波与太赫兹波在等离子体中传输特性

针对黑障区高速飞行器测控通信研究需求,在激波管设备上开展了Ka波段毫米波与太赫兹波在等离子体中传输特性的实验研究,获得了不同电子密度和碰撞频率的等离子体中传输衰减特性实验数据。采用辅助差分方程的时域有限差分（ADE-FDTD）方法对实验进行了数值模拟,数值模拟结果和实验结果有很好的一致性。实验结果和数值模拟结果均表明:太赫兹波信号在相同的实验等离子体中传输衰减比毫米波信号小得多,具有更强的穿透等离子体能力;毫米波与太赫兹波信号在等离子体中传输衰减量随着等离子体电子密度的增加而增加,二者的差值也增加;太赫兹波能够作为解决高密度等离子体中电磁波传输的有效技术手段。     

热粘弹波在变温非均匀合金熔体中的传播

机械波在金属凝固过程中传播的定量计算一直是一个难题,主要原因就是在这个过程中的熔体结构非常复杂.本研究考虑到熔体的变温、非均匀和粘弹性的特点,采用Kelvin粘弹性介质模型,建立了具有粘热损失特性的热粘弹性波动方程,通过隐式有限差分方法对波动方程进行求解,并以ZL203A合金熔体为研究对象,探究了热粘弹波在变温非均匀介质中的传播规律.结果表明:热粘弹波从合金熔体的低温区向高温区传播时,非均匀的温度场对波的传播有较大影响;相反,当波从合金熔体的高温区向低温区传播时,非均匀的温度场对波的传播几乎没有影响.热粘弹波在合金熔体中的衰减系数随频率的增大呈线性增大,而随温度的升高先增大后减小,在熔体的枝晶搭接温度附近达到最大值.     

曲线坐标系下孔隙介质圆柱纵向表面波的传播特性

为了研究孔隙介质圆柱纵向表面波的传播规律,分析其频散和衰减特性,在正交曲线坐标系下建立了表面波的频散方程,通过数值计算得到频散曲线,将纵向导波最低模态与表面波进行对比,并分析了曲率半径及孔隙参数对表面波频散和衰减的影响。结果表明,当频率足够大时,导波最低模态的频散曲线与表面波近似;在同一频率下,表面波的相速度随曲率半径的增大而增大,随孔隙度的增大而减小;表面波的衰减随孔隙度的增大而增大。研究结果为开展孔隙介质圆柱结构的超声无损评价提供了一定的理论参考。     

Finite element analysis of broadband acoustic pulses through inhomogenous media with power law attenuation

Acoustic waves in tissues and weakly attenuative fluids often have an attenuation parameter, alpha(omega), satisfying alpha(omega)= alpha0omegay in which alpha0 is a constant, omega is the frequency, and y is between 1 and 2. This power law attenuation is not predicted by the classical thermoviscous wave equation and researchers have proposed different modified viscous wave equations in which the loss term is a convolution operator or a fractional spatial or temporal derivative. In this paper, acoustic waves undergoing power law attenuation are modeled by a modification to the thermoviscous wave equation in which the time derivative of the viscous term is replaced by a fractional time derivative. An explicit time domain, finite element formulation leads to a stable algorithm capable of simulating axisymmetric, broadband acoustic pulses propagating through attenuative and dispersive media. The algorithm does not depend on the Born approximation, long wavelength limit, or plane wave assumptions. The algorithm is validated for planar and focused transducers and results include radiation patterns from a viscous scatterer in a lossless background and signals reflected from a viscous layer. The program can be used to determine scattering parameters for large, strong, possibly viscous scatterers, in either a lossless or viscous background, for which analytic results are scarce.     

非线性对大气介质中阵列聚焦声场分布影响的研究

基于时域有限差分算法将大气中近似到二阶微小项的非线性声波波动方程进行离散化,得到了模拟采用的差分波动方程.在此基础上,数值模拟了初始声压强弱不同的5个点声源组成的线阵列垂直或斜向辐射的连续正弦波在大气中传播时二维声场的分布情况.将线性条件下的模拟结果与非线性条件下的模拟结果进行比较后发现:弱非线性会对声场的分布和阵列聚焦增益产生一定的影响,使声场分布波形比线性条件下的声场分布波形更加靠近阵列,聚焦效果变差;强非线性会使波形发生更严重畸变,这是由于产生了基频以外的其他频率声波引起的;非线性对斜向传播时声场分布的影响与垂直传播时的影响效果基本相同,但由于斜向辐射时的声波几何扩展造成的轴向声压衰减要大于垂直辐射时的轴向声压衰减,因此聚焦增益和强非线性的影响都将小于垂直辐射时的情况.     

毫米波在等离子体中的衰减特性研究

为了研究毫米波与等离子体的相互作用,采用时域积分法对毫米波在等离子体中传输时的衰减特性进行了仿真,初步研究了毫米波在等离子体中传输时的衰减值与等离子体密度、碰撞频率及其工作频率之间的关系.通过理论结果和仿真结果的分析比较,验证了仿真可以为理论研究提供数据支撑的可行性. 关键词： 毫米波 等离子体 时域积分法 衰减     

核爆聚变电站发电时的冲击波超压传播规律

 核爆聚变电站在发电时会产生冲击波,其超压传播规律是核爆聚变电站洞壁设计的重要基础。现有公式在描述冲击波峰值超压时存在一定的不足（如相对误差大等）。通过理论分析得到新的冲击波峰值超压公式、正压作用时间τ₊和冲击波传播系数ω,最终得到冲击波超压的传播规律。通过对比分析表明,新公式比现有公式更加适合求解核爆聚变电站发电时所产生的冲击波峰值超压。     

Analysis of plasma sheath propagation attenuation based on ray tracing

Based on the complex Snell's law of lossy media, a ray tracing method is proposed to study the propagation attenuation characteristics of electromagnetic (EM) waves in plasma sheaths. The plasma sheath is modelled as layered media. This method considers the complex ray characteristics of inhomogeneous plane EM waves, tracks the propagation rays of EM waves in each layer of media, and calculates the propagation attenuation of EM waves in each layer of media according to the propagation direction of the complex rays. The attenuation during numerical cumulative propagation is the total attenuation of EM waves through the plasma sheath. By comparing the results with that obtained from the WKB method, the accuracy of the ray tracing algorithm is proved. The results of the propagation attenuation of a blunt cone model are calculated by the proposed method, and the effects of different parameters on the EM wave propagation attenuation in the plasma sheath are analysed at different heights, velocities, incident angles, and incident positions. Studying the propagation characteristics of EM waves in the plasma sheath is of importance in application for radar target tracking, blackout communication, and other issues.     

弛豫媒质中有限束超声波的非线性传播畸变理论

文中提出弛豫媒质中有限束非线性声波方程,并采用微扰法求得由非线性传播畸变产生的高次谐波的一般解.研究表明,对高斯型声波,其谐波畸变解可以解析给出,而且其径向分布始终维持高斯函数.虽然其频散量大小会影响各次谐波的振幅,但其相速的变化却仍与对在频率的小振幅波相同.文中还用Blackstock算子将所得的结果应用于任何吸收-频散媒质,包括只能从经验得到其吸收与频率关系的一些生物媒质.     

Nonlinear acoustic wave equations with fractional loss operators

Fractional derivatives are well suited to describe wave propagation in complex media. When introduced in classical wave equations, they allow a modeling of attenuation and dispersion that better describes sound propagation in biological tissues. Traditional constitutive equations from solid mechanics and heat conduction are modified using fractional derivatives. They are used to derive a nonlinear wave equation which describes attenuation and dispersion laws that match observations. This wave equation is a generalization of the Westervelt equation, and also leads to a fractional version of the Khokhlov-Zabolotskaya-Kuznetsov and Burgers' equations.     

Transient acoustic wave propagation in rigid porous media: a time-domain approach

Wave propagation of acoustic waves in porous media is considered. The medium is assumed to have a rigid frame, so that the propagation takes place in the air which fills the material. The Euler equation and the constitutive relation are generalized to take into account the dispersive nature of these media. It is shown that the connection between the fractional calculus and the behavior of materials with memory allows time-domain wave equations, the coefficients of which are no longer frequency dependent, to be worked out. These equations are suited for direct and inverse scattering problems, and lead to the complete determination of the porous medium parameters.     

Non-Steady Wall-Bounded Flows of Viscoelastic Fluids Under Periodic Forcing

The problem of oscillating flows inside pipes under periodic forcing of viscoelastic fluids is addressed here. Starting from the linear Oldroyd-B model, a generalized Darcy’s law is obtained in the frequency domain and an explicit expression for the dependence of the dynamic permeability on the fluid parameters and forcing frequency is derived. Previous results in both viscoelastic and Newtonian fluids are here shown to be particular cases of our results. On the basis of our calculations, a possible explanation for the observed damping of local dynamic response as the forcing frequency increases is given. Good fitting with recent experimental studies of wave propagation in viscoelastic media is here exhibited. Sound wave propagation in viscoelastic media flowing inside straight pipes is investigated. In particular, we obtain the local dynamic response for weakly compressible flows.     

基于弹性模量检测骨疲劳的超声导波方法研究

研究早期诊断骨疲劳的方法是当前骨质评价方面的研究热点之一. 本文对不同弹性模量下长骨中超声导波的传播特性进行了理论分析和仿真研究.首先, 通过数值计算得到导波在管状长骨中的理论解析解.然后对管状长骨进行了时域有限差分(FDTD) 仿真, 并验证了它与理论解析解的一致性, 同时得到长骨中不同模式导波群速度、 中心频率和衰减与弹性模量的关系.研究结果表明各个导波模式的群速度和中心频率均随弹性模量的增加而增加, 而衰减随弹性模量的增加而减小.说明超声导波的传播特性参量可以反映长骨弹性模量的变化, 从而为长骨的早期疲劳诊断提供理论依据.