Sound speed in air-filled sand and 03ine shallow-layer sandy sediments

Based on the author's theory for acoustic propagation in granular media and by employing the extinction theorem, the sound speed formulae in these media were derived. The numerical computations of sound speed in 03ine sediments and air-filled sand were carried out, and the results demonstrated that under the normal atmosphere the sound speed in air- filled sand is lower than the sound speed in the air. The numerical results also indicated that the influence of scattering interaction between the grains upon the sound speed in the 03ine shallow-layer sandy sediments has to be taken into account; however, the influence of the viscous-wave interaction can be neglected. The theoretical results obtained from the rigid-granular model seem to match the measured data better than from the elastic-granular model, even the latter model fits better for the real situation, indicating further measurements are necessary in order to gain an insight into this problem thoroughly. Through an analysis of experimental data published in journals a conclusion can be drawn that the theory of granular media is suited to deal with the problems of the sound propagation in air-filled sand better than the theory of porous media.     

Prediction-step staggered-in-time FDTD: An efficient numerical scheme to solve the linearised equations of fluid dynamics in outdoor sound propagation

The finite-difference time-domain (FDTD) method to solve the linearised equations of fluid dynamics has shown to be very powerful and useful in outdoor sound propagation. Practical applications are however limited due to the large need for computational resources. The numerical discretisation influences computational efficiency to an important degree. In this paper, some possible ways to discretise temporal derivatives are studied. Two obvious ways of time-discretisation namely staggered-in-time (SIT) and a simple collocated-in-time (CIT) scheme are compared to the prediction-step staggered-in-time (PSIT) scheme. The latter is intended to be used for the calculation of sound propagation in the typical low wind speeds encountered in the outdoor environment at low heights above the earth’s surface. It was shown that the PSIT scheme is more stable than the SIT scheme, so practical calculations are possible. Computational efficiency is increased to an important degree compared to the CIT scheme. The numerical accuracy (more precisely the amplitude error) of the PSIT scheme is an important improvement upon SIT. The CIT scheme on the other hand conserves amplitude better. The amplitude error becomes larger with increasing wind speed because of some simplifications during the numerical discretisation. In low wind speeds, the PSIT algorithm can serve as an interesting compromise between numerical accuracy and the required amount of computing power.     

Tracking of time-evolving sound speed profiles with an auto-regressive state-space model

An approach for time-evolving sound speed profiles tracking in shallow water is discussed.The inversion of time-evolving sound speed profiles is modeled as a state-space estimation problem,which includes a state equation for predicting the time-evolving sound speed profile and a measurement equation for incorporating local acoustic measurements.In the paper,auto-regression(AR) method is introduced to obtain a high-order AR evolution model of the sound speed field time variations,and the ensemble Kalman filter is utilized to track the sound speed field.To validate the approach,the accuracy in sound speed estimation is analyzed via a numerical implementation using the ASIAEX experimental environment and the sound velocity measurement data.Compared with traditional approaches based on the state evolution represented as a random walk,simulation results show the proposed AR method can effectively reduce the tracking errors of sound speed,and still keep good tracking performance at low signal-to-noise ratios.     

Application of the Transmission Line Matrix method for outdoor sound propagation modelling – Part 2: Experimental validation using meteorological data derived from the meso-scale model Meso-NH

Outdoor sound prediction is both a societal concern and a scientific issue. This paper deals with numerical simulations of micrometeorological (temperature and wind) fields for environmental acoustics. These simulations are carried out using the reference meso-scale meteorological model at the Meteo-France weather agency (Meso-NH). Meso-NH predictions at very fine scales (up to 3 m), including new developments (drag force approach), are validated both numerically and experimentally under stable, unstable and neutral conditions. Then, this information can be used as input data for the acoustic propagation model. The time-domain acoustic model is based on the Transmission Line Matrix method. Its development has also been promoted for application to outdoor sound propagation, i.e. to take into account topography, ground impedance, meteorological conditions, etc. In part 1, the presentation and evaluation of the Transmission Line Matrix method showed the relevance of this method’s use in the context of environmental acoustics. Finally, simulated noise levels under different propagation conditions were compared to in situ measurements. Satisfactory results were obtained regarding the variability of the observed phenomena.     

浅海环境下数据同化声层析方法研究

讨论了适用于浅海的基于声速局部测量模型、声传播模型及声压场局部测量模型的声学数据同化方法,并给出了具体的执行算法。该算法根据最优化准则把局部声速和声压场的先验测量信息,以及声传播模型进行了有效地融合,提高了海洋环境参数估计的精度,为浅海声层析提供了新思路。同时,利用实验浅海声速测量数据,通过经验正交函数对实际海洋声速剖面进行了估计,并分析了各类噪声对水体声速场及海底声学参数估计精度的影响,验证了该执行算法的有效性。      

自回归状态空间模型下时变声速剖面跟踪方法

讨论了一种适用于浅海的时变声速剖面跟踪方法。该方法将时变水体声速剖面的反演问题建模为由描述声速剖面时变特性的状态方程与包含声压场局部测量信息的测量方程组成的状态-空间模型,提出利用自回归分析拟合方法将声速场扰动建模为高阶自回归演化模型,并通过集合卡尔曼滤波序贯地估计时间演化的海洋声速场。利用2001年亚洲海实验环境与声速测量数据,仿真分析了基于高阶自回归演化模型的时变声速剖面集合卡尔曼滤波估计方法。结果表明,相比于利用传统随机游走状态演化模型的估计方法,该改进方法可有效降低声速的跟踪误差,并且在较低信噪比条件下仍具有较好的跟踪性能。      

Simulation of the acoustic field emitted from medical linear transducer in a heterogeneous tissue

A numerical model is developed to simulate the acoustic field in heterogeneous tissue from a medical linear transducer.The coupled full-wave equation for nonlinear ultrasound is solved using a staggered-grid finite difference time domain method.The distribution of acoustic pressure and power in human abdominal wall with heterogeneities in sound speed,density,and nonlinear parameter are obtained.Compared with homogeneous medium,when sound speed in tissue is uniform and density unchanged,the acoustic energy decreases only1.8 dB in the focal region;when density in tissue is uniform and sound speed unchanged,the energy decreases 3.8 dB in the focal region,which is almost the same as heterogeneous tissue.Thus,the primary factor of the aberration of focused beam is the heterogeneous distribution of the tissue sound speed.     

Analysis and algorithms for a regularized cauchy problem arising from a non-linear elliptic PDE for seismic velocity estimation

In the present work we derive and study a non-linear elliptic PDE coming from the problem of estimation of sound speed inside the Earth. The physical setting of the PDE allows us to pose only a Cauchy problem, and hence is ill-posed. However, we are still able to solve it numerically on a long enough time interval to be of practical use. We used two approaches. The first approach is a finite difference time-marching numerical scheme inspired by the Lax–Friedrichs method. The key features of this scheme is the Lax–Friedrichs averaging and the wide stencil in space. The second approach is a spectral Chebyshev method with truncated series. We show that our schemes work because of (i) the special input corresponding to a positive finite seismic velocity, (ii) special initial conditions corresponding to the image rays, (iii) the fact that our finite-difference scheme contains small error terms which damp the high harmonics; truncation of the Chebyshev series, and (iv) the need to compute the solution only for a short interval of time. We test our numerical schemes on a collection of analytic examples and demonstrate a dramatic improvement in accuracy in the estimation of the sound speed inside the Earth in comparison with the conventional Dix inversion. Our test on the Marmousi example confirms the effectiveness of the proposed approach.     

Nonlinear pulsed ultrasound beams radiated by rectangular focused diagnostic transducers

A numerical model for simulating nonlinear pulsed beams radiated by rectangular focused transducers, which are typical of diagnostic ultrasound systems, is presented. The model is based on a KZK-type nonlinear evolution equation generalized to an arbitrary frequency-dependent absorption. The method of fractional steps with an operator-splitting procedure is employed in the combined frequency-time domain algorithm. The diffraction is described using the implicit backward finite-difference scheme and the alternate direction implicit method. An analytic solution in the time domain is employed for the nonlinearity operator. The absorption and dispersion of the sound speed are also described using an analytic solution but in the frequency domain. Numerical solutions are obtained for the nonlinear acoustic field in a homogeneous tissue-like medium obeying a linear frequency law of absorption and in a thermoviscous fluid with a quadratic frequency law of absorption. The model is applied to study the spatial distributions of the fundamental and second harmonics for a typical diagnostic ultrasound source. The nonlinear distortion of pulses and their spectra due to the propagation in tissues are presented. A better understanding of nonlinear propagation in tissue may lead to improvements in nonlinear imaging and in specific tissue harmonic imaging. Published in Russian in Akusticheskiĭ Zhurnal, 2006, Vol. 52, No. 4, pp. 560–570. This article was translated by the authors.     

The wave equation with viscoelastic attenuation and its application in problems of shallow-sea acoustics

A suitable tool for the simulation of low frequency acoustic pulse signals propagating in a shallow sea is the numerical integration of the nonstationary wave equation. The main feature of such simulation problems is that in this case the sound waves propagate in the geoacoustic waveguide formed by the upper layers of the bottom and the water column. By this reason, the correct dependence of the attenuation of sound waves in the bottom on their frequency must be taken into account. In this paper we obtain an integro-differential equation for the sound waves in the viscoelastic fluid, which allows to simulate the arbitrary dependence of acoustic wave attenuation on frequency in the time domain computations. The procedure of numerical solution of this equation based on its approximation by a system of differential equations is then considered and the methods of artificial limitation of computational domain are described. We also construct a simple finite-difference scheme for the proposed equation suitable for the numerical solution of nonstationary problems arising in the shallow-sea acoustics.     

风扇/压气机数值稳定性模型

本文给出了一种用于预测轴流压气机旋转失速发生的数值稳定性模型,并描述了模型的理论基础以及求解模型方程的方法。模型的建立基于线性稳定性理论的基础之上,通过特征值的虚部来判断系统的稳定性。一种全局性方法被用于求解离散系统的所有特征值。离散的方法包括二阶精度的有限差分方法以及切比雪夫谱配置的方法。计算结果表明,用这两种方法均可以得到准确的特征值,但是收敛的速度有所不同。另外,将文中的分析与实验进行了比较,结果表明实验中的失速点可以用这种模型进行合理的预测。     

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

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

Multiple-order derivatives of a waveguide acoustic field with respect to sound speed, density, and frequency

An adjoint perturbative method is used to derive expressions for the first- through third-order derivatives of a pressure field with respect to sound speed, density, and frequency, for the restricted case of a laterally homogenous waveguide in which environmental parameters are only a function of depth. By using a normal-mode Green's function, the three-dimensional spatial correlation required by the standard acoustic adjoint equation can be reduced to a set of one-dimensional depth integrals. The resulting expressions for the first-order derivative are similar to those obtained by previous perturbative approaches based on the depth-separated wave equation, but the approach followed here permits straightforward extension to higher-order derivatives. Explicit evaluations of the expressions for a representative shallow-water waveguide model are in excellent agreement with numerical finite-difference computations. An analysis of the expressions as a function of source-receiver range finds the contributions to the mode amplitude derivatives to be non-negligible at ranges less than a few modal interference lengths, for parameters associated with the ocean bottom. Therefore, linear perturbative inversion methods that perturb only horizontal wavenumbers and not mode amplitudes should either be used with caution or modified to incorporate the expressions presented here.     

复杂结构冲击声合成及实验验证

以复杂结构受击振动响应的时域计算为目的, 讨论了结构阻尼的计算方法, 给出一种用于冲击声合成的综合数值方法, 并进行了实验验证. 首先, 考虑到阻尼是影响瞬态振动时变特性的重要因素, 详细讨论了两种模态阻尼的计算方法; 其次, 对阻尼板的受击振动和声辐射进行了时域仿真, 并与时域有限差分法的计算结果进行对比, 显示出两种声音合成方法的计算结果具有高度的一致性; 最后, 针对有限长圆柱壳的受击振动, 将合成声与实验录音进行了对比研究. 结果表明, 合成声与实际录音的时域包络、频谱结构以及衰减趋势基本一致, 证明了采用数值方法进行冲击声合成的有效性. 关键词： 声音合成 模态阻尼 冲击声 数值方法     

Evaluation of the vibrational modes of the human skull as it relates to bone-conducted sound

Analytic and numerical models are used to study bone-conducted sound and how it relates to the vibrational modes of the human skull. The analytic model is based on the solution to the acoustic and elastic wave equations and the constraining boundary conditions for a fluid-filled elastic sphere. Both models predict that most of the acoustic energy of bone-conducted sound exists in the form of surface wave vibrations at the interface between two acoustic media rather than in the bone or cranial chamber. These surface waves have phase speeds much slower than the bulk sound speed for bone. The analytic model, based on spherical elastic shells, predicts a phase speed of 775 m/s and the first resonance frequency at 1500 Hz while the numerical solution yields approximate phase speeds of 450 m/s and provides a visual display of the surface waves and diffraction effects.     

Determination of the inhomogeneous structure of a medium from its plane wave reflection response,part II: A numerical approximation

This paper is concerned with a standard one dimensional inverse scattering problem: given the reflection response of an unknown inhomogeneous medium for plane waves under normal or oblique incidence, determine its sound speed and density structures. The problem is solved by means of a simple numerical technique which involves only fast Fourier transform operations and numerical integration of ordinary differential equations. Three cases are specifically considered: (a) sound speed is unknown, density is known; (b) sound speed is known, density is unknown; (c) sound speed and density are to be determined simultaneously. Numerical simulations performed on reflection coefficients computed in Part I for a limited band of frequencies lead to accurate reconstructions of the original structures of various media.     

A higher-order split-step Fourier parabolic-equation sound propagation solution scheme

A three-dimensional Cartesian parabolic-equation model with a higher-order approximation to the square-root Helmholtz operator is presented for simulating underwater sound propagation in ocean waveguides. The higher-order approximation includes cross terms with the free-space square-root Helmholtz operator and the medium phase speed anomaly. It can be implemented with a split-step Fourier algorithm to solve for sound pressure in the model. Two idealized ocean waveguide examples are presented to demonstrate the performance of this numerical technique.     

An Eulerian hybrid WENO centered-difference solver for elastic–plastic solids

We present a finite-difference based solver for hyper-elastic and viscoplastic systems using a hybrid of the weighted essentially non-oscillatory (WENO) schemes combined with explicit centered difference to solve the equations of motion expressed in an Eulerian formulation. By construction our approach minimizes both numerical dissipation errors and the creation of curl-constraint violating errors away from discontinuities while avoiding the calculation of hyperbolic characteristics often needed in general finite-volume schemes. As a result of the latter feature, the formulation allows for a wide range of constitutive relations and only an upper-bound on the speed of sound at each time is required to ensure a stable timestep is chosen. Several one- and two-dimensional examples are presented using a range of constitutive laws with and without additional plastic modeling. In addition we extend the reflection technique combined with ghost-cells to enforce fixed boundaries with a zero tangential stress condition (i.e. free-slip).     

Sound propagation over layered poro-elastic ground using a finite-difference model

This article presents an axisymmetric pressure-velocity finite-difference formulation (PV-FD) based on Biot's poro-elastic theory for modeling sound propagation in a homogeneous atmosphere over layered poro-elastic ground. The formulation is coded in a computer program and a simulation of actual measurements from airblast tests is carried out. The article presents typical results of simulation comprising synthetic time histories of overpressure in the atmosphere and ground vibration as well as snapshots of the response of the atmosphere-ground system at selected times. Comparisons with the measurements during airblast tests performed in Haslemoen, Norway, as well as the simulations by a frequency-wave number FFP formulation are presented to confirm the soundness of the proposed model. In particular, the generation of Mach surfaces in the ground motion, which is the result of the sound speed being greater than the Rayleigh wave velocity in the ground, is demonstrated with the help of snapshot plots.     

Effects of sediment acoustic characteristics on low-frequency acoustic field distribution in shallow water

Most of the published literatures on low-frequency underwater sound propagation are focused on the sound propagation features in the water column,while studies on sound propagation features in the sediment layer or the semi-infinite basement are rare.In this paper,based on the wave equation,a computational model for sound energy flux in the sediment layer and the basement as well as in the water column is proposed under a cylindrical coordinate system.On this basis,the effects of various sediment acoustic parameters on the sound energy distribution and the corresponding mechanisms are elaborated through numerical examples and acoustic theory.Simulation results reveal that,in a situation where sediment P-wave speed>water sound speed> sediment S-wave speed,the greater the values of density and P-wave speed in sediment,the more likely it is that the sound energy remains in the water column without leaking to the sea floor.Conversely,the influence of the variation of S-wave speed is reversed.Basement influence on the sound propagation in the fluid layer is approximately negligible if the sediment layer is sufficiently thick.