一维非线性声波传播特性

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

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.     

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

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

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.     

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.     

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.     

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.     

Nonlinear elastic properties of a model one-dimensional granular unconsolidated structure

The characteristic features of the propagation of finite-amplitude elastic waves in a model one-dimensional unconsolidated granular medium are investigated. The model medium is represented by a linear chain of 80 steel balls with a diameter of 6.5 mm each, this chain being preliminarily loaded with an external static force F. The elastic properties of the model are analyzed. The theoretical dependences of the coefficients of elasticity of the second, third, and fourth orders on the force F are obtained. The experimental setup is described. The results of studying the nonlinear effects, namely, the higher harmonic generation and the wave generation at combination frequencies, which accompany the acoustic wave propagation in the chain, are presented. For the chain of balls under study, a structural phase transition from the 1D structure to a 2D one is observed with an increase in the external compression force F applied to the balls. The results of the study are analyzed using the Hertz theory of contact interactions.     

分层半空间表面非线性瑞利波的激发*

针对非线性瑞利波在均匀分层半空间结构中的激发和传播规律进行研究。根据摄动理论和模态分解将分层半空间结构中瑞利波的二次谐波声场表示为二倍频瑞利波模式的线性组合,经由互易关系得到各模式的展开系数表达式。对不同分层半空间结构中瑞利波二次谐波的激发和传播特性进行讨论,结果表明相速度匹配的瑞利波模式其二次谐波分量随传播距离线性增长,非匹配模式的二次谐波分量则沿传播方向周期震荡传播。此外,文中定义非线性参数表征瑞利波模式产生的非线性程度,这有利于选择出具有明显非线性效应的匹配点,为后续检测工作提供理论依据,具有指导意义。     

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.     

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

为探究临床常用的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次谐波的影响。该文提出的多层组织非线性仿真模型可为高频聚焦超声换能器参数优化及制定安全、有效的术前治疗方案提供理论参考。     

非线性有限束光声效应理论

当一束受正弦调制的激光入射于光声腔中的固体样品上时,由于非线性的光声效应,在光声腔中不仅能接收到基频成份的声信号,还能接收到其二次谐波成份。本文提出一个非线性热波束方程及其相应的非线性边界条件。借助于逐步近似法在光源为高斯径向分布的情形下,求解这一方程。利用Hankel变换获得这一方程的一级与二级近似解。解析结果表明,二次谐波的热波束仍然维持高斯径向分布,而其高斯半径比基频成份小。分析还表明,二次谐波的振幅不仅与线性热参数,而且也与非线性热参数有关。后者或许能提供样品的更多的有意义的信息。综合其各种特点, 关键词：     

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

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

非线性声波方程的传统二次谐波二阶近似解不够精确,对测量材料的非线性系数造成较大误差.利用摄动法展开非线性声波方程得到一系列非齐次偏微分方程,根据低阶解的性质拟定高阶特解的形式,通过符号计算求出高次谐波特解,对所有二次谐波成分求和最终获得二次谐波的高阶近似解.在水中开展非线性声学实验验证高阶近似解,结果表明:二次谐波相对...     

Finite-difference and finite-volume methods for nonlinear standing ultrasonic waves in fluid media

In the framework of the application of high-power ultrasonics in industrial processing in fluid media, the mathematical prediction of the acoustical parameters inside resonators should improve the development of practical systems. This can be achieved by the use of numerical tools able to treat the nonlinear acoustics involved in these phenomena. In particular, effects like nonlinear distortion and nonlinear attenuation are fundamental in applications. In this paper, three one-dimensional numerical models in the time domain for calculating the nonlinear acoustic field inside a one-dimensional resonant cavity are presented and compared. They are based on the finite-difference and the finite-volume methods. These different algorithms solve the differential equations, from the linear up to the strongly nonlinear case (including weak shock). Some physical results obtained from the modelling of ultrasonic waves and a comparison of the efficiency of the different algorithms are presented.     

Cumulative solutions of nonlinear longitudinal vibration in isotropic solid bars

Based on the strain invariant relationship and taking the high-order elastic energy into account, a nonlinear wave equation is derived, in which the excitation, linear damping, and the other nonlinear terms are regarded as the first-order correction to the linear wave equation. To solve the equation, the biggest challenge is that the secular terms exist not only in the fundamental wave equation but also in the harmonic wave equation （unlike the Duffing oscillator, where they exist only in the fundamental wave equation）. In order to overcome this difficulty and to obtain a steady periodic solution by the perturbation technique, the following procedures are taken： （i） for the fundamental wave equation, the secular term is eliminated and therefore a frequency response equation is obtained; （ii） for the harmonics, the cumulative solutions are sought by the Lagrange variation parameter method. It is shown by the results obtained that the second- and higher-order harmonic waves exist in a vibrating bar, of which the amplitude increases linearly with the distance from the source when its length is much more than the wavelength; the shift of the resonant peak and the amplitudes of the harmonic waves depend closely on nonlinear coefficients; there are similarities to a certain extent among the amplitudes of the odd- （or even-） order harmonics, based on which the nonlinear coefficients can be determined by varying the strain and measuring the amplitudes of the harmonic waves in different locations.     

Nonlinear propagation and control of acoustic waves in phononic superlattices

The propagation of intense acoustic waves in a one-dimensional phononic crystal is studied. The medium consists in a structured fluid, formed by a periodic array of fluid layers with alternating linear acoustic properties and quadratic nonlinearity coefficient. The spacing between layers is of the order of the wavelength, therefore Bragg effects such as band gaps appear. We show that the interplay between strong dispersion and nonlinearity leads to new scenarios of wave propagation. The classical waveform distortion process typical of intense acoustic waves in homogeneous media can be strongly altered when nonlinearly generated harmonics lie inside or close to band gaps. This allows the possibility of engineer a medium in order to get a particular waveform. Examples of this include the design of media with effective (e.g., cubic) nonlinearities, or extremely linear media (where distortion can be canceled). The presented ideas open a way towards the control of acoustic wave propagation in nonlinear regime.     

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.