有限振幅波的三维时域建模

通过对二维算法(Texas code)做扩展,建立了求解扩展版KZK(Khokhlov-Zabolotskaya-Kuznetsov)方程的三维时域有限差分算法。该算法先将KZK方程变换成TBE(Transformed Beam Equation),然后依次求解衍射(抛物线近似条件下)、热黏滞吸收、弛豫和非线性作用。数值仿真了圆形、矩形和方形阵列的非线性声场,与前人的结果相符,证明了算法的有效性,并分析了把弛豫吸收系数当成热黏滞吸收系数引入空气中的参量阵的模拟会带来的误差。     

Time-domain modeling of nonlinear distortion of pulsed finite amplitude sound beams

This work aims to validate a time domain numerical model for the nonlinear propagation of a short pulse of finite amplitude sound beam propagation in a tissue-mimicking liquid. The complete evolution equation is simply derived by a superposition of elementary operators corresponding to the 'one effect equation'. Diffraction LD, absorption and dispersion LAD, and nonlinear distortion LNL effects are treated independently using a first order operator-splitting algorithm. Using the method of fractional steps, the normal particle velocity and the acoustical pressure are calculated plane by plane, at each point of a two-dimensional spatial grid, from the surface of the plane circular transducer to a specified distance. The LA operator is a time convolution between the particle velocity and the causal attenuation filter built after the Kramers-Kroning relations. The LNL operator is a time-based transformation obtained by following an implicit Poisson analytic solution. The LD operator is the usual Rayleigh integral. We present a comparison between theoretical and experimental temporal pressure waveform and axial pressure curves for fundamental (2.25 MHz), second, third and fourth harmonics, obtained after spectral analysis.     

Model equation for strongly focused finite-amplitude sound beams

A model equation that describes the propagation of sound beams in a fluid is developed using the oblate spheroidal coordinate system. This spheroidal beam equation (SBE) is a parabolic equation and has a specific application to a theoretical prediction on focused, high-frequency beams from a circular aperture. The aperture angle does not have to be small. The theoretical background is basically along the same analytical lines as the composite method (CM) reported previously [B. Ystad and J. Berntsen, Acustica 82, 698-706 (1996)]. Numerical examples are displayed for the amplitudes of sound pressure along and across the beam axis when sinusoidal waves are radiated from the source with uniform amplitude distribution. The primitive approach to linear field analysis is readily extended to the case where harmonic generation in finite-amplitude sound beams becomes significant due to the inherent nonlinearity of the medium. The theory provides the propagation and beam pattern profiles that differ from the CM solution for each harmonic component.     

Modeling of pulsed finite-amplitude focused sound beams in time domain.

A time-domain numerical model is presented for simulating the finite-amplitude focused acoustic pulse propagation in a dissipative and nonlinear medium with a symmetrical source geometry. In this method, the main effects responsible in finite-amplitude wave propagation, i.e., diffraction, nonlinearity, and absorption, are taken into account. These effects are treated independently using the method of fractional steps with a second-order operator-splitting algorithm. In this method, the acoustic beam propagates, plane-by-plane, from the surface of a highly focused radiator up to its focus. The results of calculations in an ideal (linear and nondissipative) medium show the validity of the model for simulating the effect of diffraction in highly focused pulse propagation. For real media, very good agreement was obtained in the shape of the theoretical and experimental pressure-time waveforms. A discrepancy in the amplitudes was observed with a maximum of around 20%, which can be explained by existing sources of error in our measurements and on the assumptions inherent in our theoretical model. The model has certain advantages over other time-domain methods previously reported in that it: (1) allows for arbitrary absorption and dispersion, and (2) makes use of full diffraction formulation. The latter point is particularly important for studying intense sources with high focusing gains.     

柱坐标系中截断光束的广义M2因子

推导出截断圆对称光束二阶矩的传输方程,它们与无截断光束的二阶矩满足相似的传输定律.因此,基于广义截断二阶矩方法,将直角坐标系中截断二维光束的二阶矩特征参数推广到柱坐标系中的截断圆对称光束,类似方法定义的广义M2因子(M2G因子)是一个传输不变量.对理论的自洽性作了物理解释.推导出柱坐标系中截断超高斯光束的二阶矩参数和M2G因子.对一些有意义的特殊情况作了讨论,并以数值计算例作了说明.     

Propagation of Cartesian beams in nonlocal nonlinear media

We introduce a very general self-trapped beam solution of the Snyder-Mitchell linear model in Cartesian coordinates. We name such a field a self-trapped Cartesian beam (CB) which is characterized by two parameters. The complex amplitude of the self-trapped CBs is described by the product of the parabolic cylinder functions and the Gaussian function. The self-trapped standard, elegant, and generalized Hermite-Gaussian beams can be obtained by treating them as the special cases of the self-trapped CBs.     

Time-domain simulation of sound production of the sho

A physical model based on the sound production mechanism of the sho is proposed with intention of applying it to sound synthesis. Time-domain simulation was done using this model, and effects of the tube length and blowing pressure on the sounding frequency and sounds spectra were investigated. The reed vibration, pressure variation inside the tube, and threshold blowing pressure for oscillation were measured by artificially blowing air into the sho. The experimental results are in acceptable agreement with simulation results in terms of the relationships between tube length and threshold pressure and between tube length and the sounding frequency. In addition, recorded sound waveforms and simulated ones have a common feature in the sense that high-frequency components of their spectra increase with increasing blowing pressure. Further, it is concluded that a sho reed acts as an "outward-striking valve."     

Time-domain simulations of sound propagation through screen-induced turbulence

A flow model in combination with a statistical-dynamical turbulence generator and a linearised Euler time-domain model for sound waves were used to simulate the effect of screen-induced turbulence on the noise level in the acoustical shadow of a screen in wind. Instead of simulating a great number of different frozen turbulence realisations, the concept of transient turbulence was successfully tested and applied. This concept is adequate to the time-domain model and reduces the computational demands. Several two-dimensional simulations allowed to isolate the individual effects of wind and screen on the propagation of 500 Hz sound waves over a 4-m high noise barrier. At a distance of 250 m from the source (240 m behind the screen) the sheltering effect of the screen and the refraction effect of the wind are in the order of 6 and 4 dB, respectively. The screen-induced turbulence leads to fluctuations in the noise level with a standard deviation of 1.2 dB and a maximum amplitude of 3 dB. However, the time averaged effect turned out to be in the order of merely 0.2 dB. The effect of the screen-induced turbulence on the average noise level behind the screen is therefore negligible.     

Simultaneous Cartesian coordinate display of defocused optical transfer functions

We show that by using a binary spatial filter and a square-law detector we can display all the defocused optical transfer functions (OTF's) in a given study in a single picture. The resulting unique picture has as its horizontal coordinates the spatial frequency and as its vertical coordinates the amount of defocus. The gray-level variations are proportional to the modulus of the OTF, that is, the modulation transfer function. Numerical simulations are included.     

球坐标中三维各向同性谐振子的类经典态

把坐标平均值随时间的变化和在宏观条件下与经典解相同的量子态定义为类经典态(NCS), 并求解球坐标中三维各向同性谐振子的NCS问题, 有助于从波动力学角度理解量子到经典过渡的问题. 选与经典态相应的大量子数附近的矩形波包作为NCS, 得到与经典解一致的结果, 但NCS不是惟一的. 一个经典态可以有很多NCS与之对应, 就像一个热力学态可以有无数力学态与之对应一样, 从量子到经典的描述是一个粗粒化和信息丢失的过程.     

有限振幅声波间的非线性相互作用对声源远场指向性的影响

根据Fenlon理论推导得到了多频声源的辐射声压.基于单频声源谐波指向性的求解方法,得到了二阶近似下的双频声源辐射出的声波相互作用时的远场指向性.分别研究和讨论了在初始辐射声压和频率不同情况下,两列波的相互作用对其中一列波的一阶波和二阶波远场指向性的影响.结果表明,声波间的相互作用对声源远场指向性的影响根据各波之间的相对初始辐射声压和相对频率的不同而有所变化. 关键词： 有限振幅声波 非线性相互作用 远场指向性     

Measurement and modeling of three-dimensional sound intensity variations due to shallow-water internal waves

Broadband acoustic data (30-160 Hz) from the SWARM'95 experiment are analyzed to investigate acoustic signal variability in the presence of ocean internal waves. Temporal variations in the intensity of the received signals were observed over periods of 10 to 15 min. These fluctuations are synchronous in depth and are dependent upon the water column variability. They can be explained by significant horizontal refraction taking place when the orientation of the acoustic track is nearly parallel to the fronts of the internal waves. Analyses based on the equations of vertical modes and horizontal rays and on a parabolic equation in the horizontal plane are carried out and show interesting frequency-dependent behavior of the intensity. Good agreement is obtained between theoretical calculations and experimental data.     

Time-domain calculations of sound interactions with outdoor ground surfaces

A time-domain formulation for sound propagation in rigid-frame porous media, including waveform attenuation and dispersion, is developed. The new formulation is based on inversion of the relaxation functions from a previous model [Wilson DK, Ostashev VE, Collier SL. J Acoust Soc Am 2004;116:1889-92], thereby casting the convolution integrals in a form amenable to numerical implementation. Numerical techniques are developed that accurately implement the relaxational equations and transparently reduce to previous results in low- and high-frequency limits. The techniques are demonstrated on calculations of outdoor sound propagation involving hills, barriers, and ground surfaces with various material properties. We also compare the relaxation formulation to a widely applied phenomenological model developed by Zwikker and Kosten. The two models can be made equivalent if the resistance constant, structure constant, and compression modulus in the ZK model are allowed to be weakly frequency dependent. But if the ZK parameters are taken to be constant, as is typically the case, the relaxation model provides more accurate calculations of attenuation by acoustically soft porous materials such as snow, gravel, and forest litter.     

阻尼受击板的时域建模及声音合成

为满足声源辨识中对合成冲击声的迫切需求,建立了球-板撞击的时域模型,提出一种时域快速求解方法,并进行了实验验证.首先给出一种将时域有限差分法(FDTD)和模态展开法(MEM)相结合的时域混合方法,求解板的振动方程,并解决了混合方法中MEM的模态截断和初值问题,及两种方法中阻尼的一致性问题;随后,给出了简支矩形板的冲击声计算结果,通过与FDTD方法的运算量进行对比,验证了混合方法的高效性;最后,针对自由边界下的L形板进行了实验验证.结果表明,与传统FDTD方法相比,时域混合方法在保证合成冲击声精度的前提下可将计算效率提高100至1200倍。     

柱坐标系中截断光束的广义M2因子

推导出截断圆对称光束二阶矩的传输方程，它们与无截断光束的二阶矩满足相似的传输定律.因此，基于广义截断二阶矩方法，将直角坐标系中截断二维光束的二阶矩特征参数推广到柱坐标系中的截断圆对称光束，类似方法定义的广义M²因子（M²_G因子）是一个传输不变量.对理论的自洽性作了物理解释.推导出柱坐标系中截断超高斯光束的二阶矩参数和M²_G因子.对一些有意义的特殊情况作了讨论，并以数值计算例作了说明. 关键词： 光束描述 截断光束 圆对称性 广义截断二阶矩(GTSM) 2因子（M²_G因子）')" href="#">广义M²因子（M²_G因子）     

Numerical modeling of the exterior-to-interior transmission of impulsive sound through three-dimensional,thin-walled elastic structures

Exterior propagation of impulsive sound and its transmission through three-dimensional, thin-walled elastic structures, into enclosed cavities, are investigated numerically in the framework of linear dynamics. A model was developed in the time domain by combining two numerical tools: (i) exterior sound propagation and induced structural loading are computed using the image-source method for the reflected field (specular reflections) combined with an extension of the Biot–Tolstoy–Medwin method for the diffracted field, (ii) the fully coupled vibro-acoustic response of the interior fluid–structure system is computed using a truncated modal-decomposition approach. In the model for exterior sound propagation, it is assumed that all surfaces are acoustically rigid. Since coupling between the structure and the exterior fluid is not enforced, the model is applicable to the case of a light exterior fluid and arbitrary interior fluid(s). The structural modes are computed with the finite-element method using shell elements. Acoustic modes are computed analytically assuming acoustically rigid boundaries and rectangular geometries of the enclosed cavities. This model is verified against finite-element solutions for the cases of rectangular structures containing one and two cavities, respectively.     

Time-reversed sound beams of finite amplitude.

Numerical simulations based on the nonlinear parabolic wave equation are used to investigate time reversal of sound beams radiated by unfocused and focused sources. Emphasis is placed on nonlinear propagation distortion in the time-reversed beam, and specifically its effect on field reconstruction. Distortion of this kind, due to amplification during time reversal, has been observed in recent experiments [A. P. Brysev et al., Acoust. Phys. 44, 641-650 (1998)]. Effects of diffraction introduced by time-reversal mirrors with finite apertures are also considered. It is shown that even in the presence of shock formation, the ability of time reversal to retarget most of the energy on the source or focal region of the incident beam is quite robust.     

Diffraction of three-dimensional beams in uniaxial media

Within here we make an analysis of two and three-dimensional beams synthesized by extraordinary waves that propagate in uniaxial media. A relation between the geometric place of the first interference maxima and the energy flow direction is established. In addition, we determine that three-dimensional extraordinary symmetrical beams can be obtained from two two-dimensional ones that are not to be necessarily orthogonal.