Numerical simulation of nonlinear propagation of sound waves in a finite horn

Nonlinear acoustic propagation generated by a piston in a finite horn is numerically studied.A quasi-one-dimensional nonlinear model with varying cross-section uses high-order low-dispersion numerical schemes to solve the governing equation.Because of the nonlinear wave distortion and reflected sound waves at the mouth,broadband time-domain impedance boundary conditions are employed.The impedance approximation can be optimized to identify the complex-conjugate pole-residue pairs of the impedance functions,which can be calculated by fast and efficient recursive convolution.The numerical results agree very well with experimental data in the situations of weak nonlinear wave propagation in an exponential horn,it is shown that the model can describe the broadband characteristics caused by nonlinear distortion.Moreover,finite-amplitude acoustic propagation in types of horns is simulated,including hyperbolic,conical,exponential and sinusoidal horns.It is found that sound pressure levels at the horn mouth are strongly affected by the horn profiles,the driving velocity and frequency of the piston.The paper also discusses the influence of the horn geometry on nonlinear waveform distortion.     

有限长号筒中非线性声传播的数值模拟研究

研究了大振幅活塞声源经有限长号筒向外辐射声波的非线性声学问题。采用具有频散保持特性的高精度计算格式求解了适用于变截面管道的一维非线性声场模型,并考虑到声波的非线性畸变及管口处的声反射,加入了宽频时域声阻抗边界条件。宽频阻抗模型的共轭复数系数采用优化拟合方法近似求解,并采用递推卷积算法快速求解出时域声阻抗。在弱非线性条件下模拟指数形号筒中的声传播取得了与已有实验相一致的结果,表明模型能够描述声波非线性畸变带来的宽频特性。进而采用本模型数值模拟了大振幅活塞声源在双曲形、锥形、指数形和正弦形号筒中的非线性声传播问题,结果表明号筒出口声压级受活塞振动速度、频率以及号筒形状的影响,并分析讨论了波形畸变与号筒几何形状之间的关系。      

一维非线性声波传播特性

针对一维非线性声波的传播问题进行了有限元仿真和实验研究. 首先推导了一维非线性声波方程的有限元形式, 含有高阶矩阵的非线性项导致声波具有波形畸变、谐波滋生、基频信号能量向高次谐波传递等非线性特性. 编制有限元程序对一维非线性声波进行了计算并对仿真得到的畸变非线性声波信号进行处理, 分析其传播性质和物理意义. 为验证有限元计算结果, 开展了水中的非线性声波传播的实验研究, 得到了不同输入信号幅度激励下和不同传播距离的畸变非线性声波信号. 然后对基波和二次谐波的传播性质进行详细讨论, 分析了二次谐波幅度与传播距离和输入信号幅度的变化关系及其意义, 拟合出二次谐波幅度随传播距离变化的方程并阐述了拟合方程的物理意义. 结果表明, 数值仿真信号及其频谱均与实验结果有较好的一致性, 证实计算方法和结果的正确性, 并提出了具有一定物理意义的二次谐波随传播距离变化的简单数学关系. 最后还对固体中的非线性声波传播性质进行了初步探讨. 本研究工作可为流体介质中的非线性声传播问题提供理论和实验依据.     

Estimation methodology for amplitude of second-order harmonic in nonlinear surface acoustic wave measured by optical diffraction method

An estimation methodology for amplitude from the fundamental to fourth-order harmonic composed of nonlinear surface acoustic waves on a LiNbO₃ substrate with a 128° Y-X propagation axis was demonstrated. First, equations for estimation of amplitudes from diffracted light intensity were shown and then numerical simulations of main factors were carried out. Normalized diffraction light intensities due to components of nonlinear surface acoustic waves were measured and their amplitudes were obtained from measured values. Finally, the amplitude of the second harmonic propagating on a gold thin film of about 50 nm thickness was estimated. It was found that the amplitude of the second harmonic increased continuously from 0.8 to 2.0 Å with propagation distance in the thin film.     

PARAMETRIC SIMULTONS IN ONE-DIMENSIONAL NONLINEAR LATTICES

Parametric simultaneous solitary wave (simulton) excitations are shown to be possible in nonlinear lattices. Taking a one-dimensional diatomic lattice with a cubic potential as an example, we consider the nonlinear coupling between the upper cut-off mode of acoustic branch (as a fundamental wave) and the upper cut-off mode of optical branch (as a second harmonic wave). Based on a quasi-discreteness approach the Karamzin－Sukhorukov equations for two slowly varying amplitudes of the fundamental and the second harmonic waves in the lattice are derived when the condition of second harmonic generation is satisfied. The lattice simulton solutions are given explicitly and the results show that these lattice simultons can be nonpropagating when the wave vectors of the fundamental wave and the second harmonic waves are exactly at π/a (where a is the lattice constant) and zero, respectively.     

Propagation of acoustic wave in viscoelastic medium permeated with air bubbles

Based on the modification of the radial pulsation equation of an individual bubble, an effective medium method (EMM) is presented for studying propagation of linear and nonlinear longitudinal acoustic waves in viscoelastic medium permeated with air bubbles. A classical theory developed previously by Gaunaurd (Gaunaurd GC and \"{U}berall H, {\em J. Acoust. Soc. Am}., 1978; 63: 1699--1711) is employed to verify the EMM under linear approximation by comparing the dynamic (i.e. frequency-dependent) effective parameters, and an excellent agreement is obtained. The propagation of longitudinal waves is hereby studied in detail. The results illustrate that the nonlinear pulsation of bubbles serves as the source of second harmonic wave and the sound energy has the tendency to be transferred to second harmonic wave. Therefore the sound attenuation and acoustic nonlinearity of the viscoelastic matrix are remarkably enhanced due to the system's resonance induced by the existence of bubbles.     

Acoustic particle velocity horns

The paper considers receiving acoustic horns designed for particle velocity amplification and suitable for use in vector sensing applications. Unlike conventional horns, designed for acoustic pressure amplification, acoustic velocity horns (AVHs) deliver significant velocity amplification even when the overall size of the horn is much less than an acoustic wavelength. An AVH requires an open-ended configuration, as compared to pressure horns which are terminated at the throat. The appropriate formulation, based on Webster's one-dimensional horn equation, is derived and analyzed for single conical and exponential horns as well as for double-horn configurations. Predicted horn amplification factors (ratio of mouth-to-throat radii) were verified using numerical modeling. It is shown that three independent geometrical parameters principally control a horn's performance: length l, throat radius R(1), and flare rate. Below a predicted resonance region, velocity amplification is practically independent of frequency. Acoustic velocity horns are naturally directional, providing maximum velocity amplification along the boresight.     

On reconstructing photon statistics by on/off detectors: Toward the multi-partite case

The propagation of pulses of the pump and its second harmonic in a quadratically nonlinear medium whose linear properties are characterized by a negative refractive index at the pump frequency and by a positive refractive index at the harmonic frequency is considered theoretically. In the case of a low intensity of the interacting waves, the pump and second-harmonic pulses propagate in the opposite directions, but sufficiently powerful pulses can form a simulton—a solitary two-frequency wave propagating in a certain direction as a single whole. Solutions to a set of equations are found which describe the steady-state propagation of a solitary wave and of a nonlinear periodic (cnoidal) wave.     

Computing for second harmonic generation in one-dimensional nonlinear photonic crystals

A numerical study of second harmonic generation (SHG) in one-dimensional nonlinear photonic crystals based on full nonlinear system of equations, implemented by a combination of the method of finite elements and fixed-point iterations, is reported. This model is derived from a nonlinear system of Maxwell’s equations, which partly overcomes the known shortcoming of some existing models relied on the undepleted-pump approximation. We derive a general solution of SHG in one-dimensional nonlinear photonic crystals structures. The convergence of our method is fast. Numerical simulations also show the conversion efficiency of SHG can be significantly enhanced when the frequencies of the fundamental wave are located at the photonic band edges or are assigned to the designed defect states.     

Particle Simulation of Electrostatic Waves Driven by an Electron Beam

A one-dimensional electrostatic particle-in-cell simulation is performed to study electrostatic wave excitation due to an electron beam in a plasma system. The excited fundamental and harmonic waves are analyzed with the fast Fourier transformation and the wavelet transformation. The second harmonic is suggested to be generated by wave-wave coupling during the nonlinear evolution, which involves forward propagating and backward propagating Langmuir waves. Furthermore, the background electrons may be heated and accelerated by the electrostatic waves.     

Reflection and refraction of surface acoustic waves by a periodic domain structure

Reflection and refraction of surface acoustic waves by a periodic domain structure formed in lithium niobate is studied. A second harmonic generation is observed. A mechanism underlying the linear and nonlinear interactions of acoustic waves with a periodic domain structure is proposed.     

Methods of enhancement of the efficiency of generation of the second optical harmonic of spatially limited laser beams

Frequency doubling of spatially limited laser beams in nonlinear optical crystals is analyzed. Conditions are found under which the diaphragm aperture effect limiting the efficiency of the process of second optical harmonic generation can be reduced. It is shown that a spatial shift of mutually orthogonally polarized beams of fundamental radiation in the direction perpendicular to the direction of propagation of the waves, as well as the optimization of the ratio of the spatial dimensions and the intensity at the entrance of the nonlinear medium, gives rise to an increase in the efficiency of second harmonic generation.     

Nonlinear imaging of isoechogenic phantoms using phase conjugation of the second acoustic harmonic

Experiments on nonlinear C-scan acoustic imaging were performed using phase conjugation of the second harmonic of a focused incident beam. This method made it possible to obtain images of objects with close linear properties but with different nonlinear parameters. Nonlinear backward propagation of the conjugate wave was also considered and improvement of the image contrast was demonstrated for this case. The results obtained were generally in agreement with the theoretical predictions.     

On the Evolution Equations of Nonlinearly Permissible,Coherent Hole Structures Propagating Persistently in Collisionless Plasmas

A fundamental study describing nonlinear plasma wave propagation is presented. Elementary linear wave theory describes small-amplitude random waves, but lacks information about coherent structures. This improved wave model arises from the fact that structure formation is inevitably associated with particle trapping, which can only be properly addressed by the pseudo-potential method instead of Bernstein, Greene, and Kruskal (BGK) - likemethods. Only by using this method can legitimate nonlinear dispersion relations be obtained and reconciled with trapping scenarios. This privilege is used to derive evolution equations for five structures, the derivation being simplified by the acoustic nature of the permitted modes. The focus is on a special structure, the solitary electron hole of negative polarity, with which it can explain a spacecraft observation for the first time. Furthermore, it is shown that an intrinsically nonlinear structure can become macroscopically linear and thus harmonic by suitably adjusting the trapping scenario. An example is the monochromatic ion acoustic wave that propagates at ion sound velocity without dispersion. In this literature research, it also takes a critical look at a recently awarded work.     

Acoustic wave propagation in a binary mixture of inviscid fluids

Linearized equations governing the thermo-mechanical behaviour of a binary mixture of inviscid fluids are derived. Restrictions which are sufficient for the equations to have a unique solution are imposed on some of the material constants. The propagation of plane harmonic waves of small amplitude in the mixture is examined and the inequalities are shown to ensure a physically reasonable response. As an application of the theory properties of acoustic waves in a binary mixture of ideal gases are evaluated numerically.     

谐振管内非线性驻波的有限体积数值算法

为了扩展谐振管内非线性驻波在工程中的应用, 以及克服现有数值计算方法仅局限于求解直圆柱形和指数形谐振管内非线性驻波的问题. 根据变截面的非稳态可压缩热黏性流体Navier-Stokes方程和空间守恒方程, 并基于求解压力速度耦合方程的半隐式算法和交错网格技术, 构建一种能够计算任意形状轴对称谐振管受活塞驱动时内部非线性驻波的有限体积算法. 分别对圆柱形、指数形和圆锥形谐振管内的非线性驻波进行仿真计算. 通过与现有试验结果以及数值仿真结果的对比, 验证了该方法的正确性.并获得除驻波声压之外的另外一些新的物理结果, 包括速度、密度、温度的瞬时变化.在直圆柱形谐振管内产生冲击声压波, 速度波形中出现钉状结构.而在指数形和圆锥形谐振管内产生高声压幅值的驻波, 没有出现冲击波, 速度波形中均未发现钉状结构. 计算结果表明谐振管内非线性驻波的物理属性与谐振管形状之间有密切关系.     

非线性声波方程二次谐波高阶近似解及实验验证

非线性声波方程的传统二次谐波二阶近似解不够精确,对测量材料的非线性系数造成较大误差。利用摄动法展开非线性声波方程得到一系列非齐次偏微分方程,根据低阶解的性质拟定高阶特解的形式,通过符号计算求出高次谐波特解,对所有二次谐波成分求和最终获得二次谐波的高阶近似解。在水中开展非线性声学实验验证高阶近似解,结果表明:二次谐波相对幅度随传播距离和输入信号均呈现先增加后减小的趋势,高阶近似解即使在传播距离不太短和输入信号不太弱的宽泛条件下也能更好地吻合实验结果。得到的高阶近似解不仅弥补了现有二次谐波幅度理论公式的不足,进行非线性声学实验时还可扩展测量范围,充分利用实验数据也能提高非线性系数的测量精度。      

高频聚焦超声声场和温度场的仿真研究

为探究临床常用的7 MHz高频聚焦超声在多层生物组织中的声传播以及毫秒级时间内的生物传热规律问题,基于Westervelt方程和Pennes传热方程,使用有限元方法建立高频聚焦超声辐照多层组织的非线性热黏性声传播及传热模型。首先分析了线性模型和非线性模型之间的差异,然后在非线性模型下探究换能器的参数对声场和温度场的影响。仿真结果显示：在7 MHz频率下,当换能器输出声功率超过5 W时,声波传播的非线性效应不可忽视(p <0.05);当声功率从5 W增大到15 W时,非线性模型与线性模型预测的温度偏差从20%增加到34.703%;高频聚焦超声波的非线性行为比低频更加显著,基频能量向高次谐波转移的程度增大,声功率为10 W和15 W时4次谐波与基波之比分别达到7.33%和12.12%;高频换能器参数的改变对组织中声场和温度场分布的影响较大,换能器焦距从12 mm减小到11.2 mm,焦点处最高温度增加了77%。结果表明,7 MHz聚焦超声的非线性声传播需要考虑到4次谐波的影响。该文提出的多层组织非线性仿真模型可为高频聚焦超声换能器参数优化及制定安全、有效的术前治疗方案提供理论参考。     

Improved contrast to tissue ratio at higher harmonics

The challenge in ultrasound contrast imaging is a better discrimination between the perfused tissue and the contrast bubbles, which is usually expressed by contrast to tissue ratio (CTR). Imaging based on the second harmonic frequency showed a higher CTR than imaging at the fundamental frequency. However, because of nonlinear propagation of ultrasound waves, harmonic frequencies are generated. These harmonic frequencies will be linearly reflected by the tissue and therefore limit the CTR at the second harmonic frequency. In order to reduce the scattering of tissue at harmonic frequencies and by that increase the CTR, nonlinear distortion has to be reduced. We demonstrate in this study that the CTR increases with the harmonic number. The increase is substantial when transmitting at lower frequencies. To take advantage of the higher harmonics (third, fourth, fifth and the ultraharmonics and termed here super harmonics), we have developed a new phased array transducer with a wide frequency band. In-vitro measurements using the new probe show an increase of 40 dB of the CTR for super harmonic components over the conventional second harmonic system. The increase in CTR is in agreement with the calculations using existing models for the response of encapsulated bubbles and known theory of nonlinear propagation.