缺陷岩体中弹性波传播的多次散射效应分析

弹性波在岩体中传播时与岩体缺陷相互作用形成复杂的传播图案。为研究缺陷对弹性波多次散射作用的影响,建立了双椭圆缺陷模型,基于Green函数基本解,采用边界积分的计算方法,得到了反映缺陷界面条件的刚度矩阵,分析了弹性波在双椭圆缺陷间的多次散射效应。结果表明：与单椭圆缺陷模型相比,双缺陷的相互作用使得弹性波频散和衰减效应增强,定量给出了缺陷的影响区域,从而明确了多次散射效应的尺度界限。进一步探讨了弹性波传播的多尺度效应,结果表明频散的Rayleigh峰、Mie峰和衰减的峰值频率同椭圆长轴和入射波波长两个尺度密切相关,存在明确的定量关系。相应的数值模拟结果表明,弹性波和缺陷相互作用在缺陷界面上诱发界面波,该界面波也存在频率相关性,影响了弹性波宏观传播的频散和衰减特征。     

Acoustical scattering identification with local impedance through a spectral approach

Acoustical scattering in waveguides is studied in this paper. The Wave Finite Element (WFE) approach is mainly used, since it allows the reduction of problems dealing with periodic waveguides. The paper deals with guided acoustical propagation, that is, propagation in a main direction is privileged. The scattering by a locally reacting lining is first studied. The liner can be characterised by its local impedance in this case. The equivalent surface impedance is therefore calculated. Then, scattering by a porous layer is considered. A full three-dimensional modelling of the lining is preferred since porous materials are bulk reacting. The scattering matrix of the lined part is computed, and acoustical scattering of high-order modes and conversion between modes are highlighted. The acoustic power attenuation is further evaluated. The response of ducts subjected to constraining boundary conditions is also calculated. Numerical results are presented and compared to those obtained with conventional approaches.     

Resonance scattering characteristics of double-layer spherical particles

Based on the principle of ultrasonic resonance scattering, sound-scattering characteristics of double-layer spherical particles in water were numerically studied in this paper. By solving the equations of the scattering matrix, the scattering coefficient determined by the boundary conditions can be obtained, thus the expression for the sound-scattering function of a single double-layer spherical particle can be derived. To describe the resonance scattering characteristics of a single particle, the reduced scattering cross section and reduced extinction cross section curves were found through numerical calculation. Similarly, the numerically calculated sound attenuation coefficient curves were used to depict the resonance scattering characteristics of monodisperse and polydisperse particles. The results of numerical calculation showed that, for monodisperse particles, the strength of the resonance was mainly related to the particle size and the total number of particles; while for polydisperse particles, it was primarily affected by the particle size, the coverage of the particle size distribution and the particle concentration.     

基于积分方程方法的弹性波多重散射计算

运用积分方程方法计算了含多个随机分布椭圆柱型孔洞的随机非均匀介质中相干波的速度和衰减系数,分析了这种介质的频散特性。首先,建立了散射位移场满足的积分方程,推导了单个椭圆柱孔洞的散射截面计算公式。接着分析了在含多个随机分布椭圆柱型孔洞的随机非均匀介质中弹性波的多重散射,给出在统计平均意义下的相干波的波速和衰减系数计算公式。然后用Matlab进行了编程,给出了一个数值算例,并将计算结果与波函数展开法进行了比较,分析了随机空隙介质的频散特征及其孔洞椭圆偏心率和材料空隙率的影响。     

Wave dispersion and attenuation in fresh mortar: theoretical predictions vs. experimental results

In the present paper, the dispersive and attenuative behavior of fresh cementitious material is examined through a series of ultrasonic, through-transmission measurements. The sand size and content dominate attenuation behavior, while the effect of entrapped air bubbles is mostly obvious at lower frequencies elevating phase velocity to values much higher than that of water. Theoretical investigation seems to explain the observed dispersion and attenuation mainly through two scattering interactions: sand embedded in paste and air bubbles in mortar. The predictions made by scattering theory follow closely the experimental data. The possibility of material characterization is discussed.     

Mean field attenuation and amplitude attenuation due to wave scattering

We point out the inadequacy of two widely used approaches of formulating the amplitude attenuation of seismic waves, the formulation of mean-field attenuation and that of scattering coefficient under the single scattering approximation. Using a one-dimensional layered slab, we show that the attenuation of the mean field is merely a statistical effect caused by phase interference among different realizations of the random wave ensemble, and does not represent the amplitude attenuation. We will call the attenuation coefficient of the mean field as the randomization coefficient in order to distinguish it from the amplitude attenuation coefficient.We also show that the scattering coefficient method leads to the same result as the first order approximation to the renormalized perturbation series (or the bilocal approximation to Dyson's equation) of the mean field. Therefore these two approaches are equivalent to a certain degree. This is also shown by using the one-dimensional layered slab model as an example.After pointing out the incorrectness of comparing experiments on amplitude attenuation with the mean field formulation, we suggest and discuss some methods of obtaining the mean field in experiments. For one-dimensional problems, the samples must be taken along the whole propagation path in order to use a spatial average to substitute for ensemble average. For a three-dimensional wave field, measurements over a large seismic array can be used to obtain the mean field. The data from Lasa measured by Aki are used to compare with theory; the agreement between them is good. Finally we compare the mean-field attenuation (randomization) and the amplitude attenuation using the back-halfspace-integration approximation introduced by Wu, and compare them with the measured data by Aki. The comparison shows further the inability of the mean-field formulation in dealing with the problem of amplitude attenuation.     

弹性波绕任意形状界面孔的散射

求解了弹性波绕任意形状界面孔的散射问题．通过入射波、反射波或折射波及孔的散射波场的叠加,得到了界面孔在SH波绕射下的总波场．总波场波函数的级数项待定系数可采用边界配点法来确定,该法不受边界正交性的限制,能够适用于任意形状的边界．最后,对界面椭圆孔进行了实例计算,得到了椭圆孔边的动应力集中系数．     

Wave propagation in gaseous small-scale channel flows

The propagation and attenuation of an initial shock wave through a mm-scale channel of circular cross-section over lengths up to 2,000 diameters is examined as a model problem for the scaling of viscous effects in compressible flows. Experimental wave velocity measurements and pressure profiles are compared with existing data and theoretical predictions for shock attenuation at large scales and low pressures. Significantly more attenuation is observed than predicted based on streamtube divergence. Simulations of the experiment show that viscous effects need to be included, and the boundary layer behavior is important. A numerical model including boundary layer and channel entrance effects reproduces the wave front velocity measurements, provided a boundary layer transition model is included. A significant late-time pressure rise is observed in experiments and in the simulations.     

Propagation and scattering of ultrasonic waves in polycrystals with arbitrary crystallite and macroscopic texture symmetries

A general ultrasonic attenuation model for a polycrystal with arbitrary macroscopic texture and triclinic ellipsoidal grains is described with proper accounting for the anisotropic Green’s function for the reference medium. The texture and the ellipsoidal grain frames in the model are independent and the wave propagation direction is arbitrary. The attenuation coefficients are obtained in the Born approximation accompanied by the Rayleigh and stochastic asymptotes. The scattering model displays statistical anisotropy due to two independent factors: (1) shape of the oriented grains and (2) preferred crystallographic orientation of the grains leading to macroscopic anisotropy of the homogenized reference medium. The model is applicable to most single phase polycrystalline materials that may occur as a result of thermomechanical manufacturing processes leading to different macrotextures and elongated-shaped grains. It predicts the strength of ultrasonic scattering and its dependence on frequency and propagation direction as a function of grain shape, grain crystallographic symmetry and macroscopic texture parameters and provides the texture-induced dependence of macroscopic ultrasonic velocity on propagation angle. It considers proper wave polarizations due to macroscopic anisotropy and scattering-induced transformations of waves with different polarizations. Competing effects of grain shape and texture on the attenuation are observed. In contrast to the macroscopically isotropic case, where in the stochastic regime the attenuation is highest in the direction of the longest ellipsoidal axis of the grain, the wave attenuation in the elongation direction may be suppressed or amplified by the texture with different effects on the quasilongitudinal and quasitransverse waves. The frequency behavior is also interestingly affected by texture: a hump in the total attenuation coefficient is found for the fast quasitransverse wave which is purely the result of macroscopic anisotropy and the existence of two quasitransverse waves; this hump is not observed in the macroscopically isotropic case. Striking differences of the texture effect on the directional dependences of the attenuation coefficients are found at low versus high frequencies.     

Scattering attenuation of elastic waves due to low-contrast inclusions

Seismic scattering attenuation due to random lithospheric heterogeneity has been theoretically modeled using two approaches. One approach is the Born approximation theory (BAT), which is primarily used to treat weak continuous heterogeneity, and the other approach is the Foldy approximation theory (FAT), which deals with sparsely distributed discrete inclusions. We apply the BAT to elastic wave scattering due to inclusions having low contrast with the matrix, and compare the results with those predicted by the FAT. We thus investigate the valid wavenumber range of the BAT based on a reasonable assumption that the inclusions are distributed so sparsely that the FAT is effectively correct for any wavenumber. For simplicity, we consider a specific type of round inclusion, which is either two- or three-dimensional and has a two-valued wave velocity and/or mass density. Both theories are confirmed to yield essentially equivalent results below a certain wavenumber limit, depending on the contrast. This is known as the Rayleigh-Gans scattering regime. Beyond the wavenumber limit, the BAT overestimates the attenuation for common-mode scattering due to wave-velocity contrast, but remains valid with respect to the attenuation for scattering due to mass-density contrast and/or conversion scattering. These conclusions are independent of the spatial dimensions of the media as well as the modes of the elastic waves (P or S). Some advantages of the BAT over the FAT for application to low-contrast inclusions are discussed.     

Scattered noise prediction using acoustic velocity formulations V1A and KV1A

Aeroacoustic scattering prediction generally relies on boundary integral methods which require evaluation of the impermeability condition on the scattering surface. The boundary condition implies zero normal velocity relative to the scattering surface. This condition has been expressed by relating the acoustic velocity to the acoustic pressure gradient, allowing indirect evaluation of the boundary condition by existent acoustic pressure gradient formulations. In the present paper, a direct evaluation of the hardwall boundary condition in scattering problems is demonstrated by time-domain analytic acoustic velocity formulae. Acoustic velocity formulations V1A and KV1A are implemented for acoustic scattering prediction, by hybrid approaches based on the FW–H equation and the Kirchhoff method These formulations can be coupled to any scattering solver, allowing time-domain prediction of the incident acoustic field when broadband noise generation is concerned. Formulation V1A offers mathematical simplicity and computational efficiency, which can be advantageous for realistic scattering applications. Implementation of formula KV1A enables acoustic scattering prediction by existing solvers based on the Kirchhoff method. The validity of the suggested methodology is assessed through the analytical test case of harmonic sound scattered by a rigid sphere. Sound propagation and scattering effects are analyzed by examination of the acoustic velocity field characteristics.     

SCATTERING OF CIRCULAR CAVITY IN RIGHT-ANGLE PLANAR SPACE TO STEADY SH-WAVE

Complex function method and multi-polar coordinate transformation technology are used here to study scattering of circular cavity in right-angle planar space to SH-wave with out-of-plane loading on the horizontal straight boundary. At first, Green function of right-angle planar space which has no circular cavity is constructed; then the scattering solution which satisfies the free stress conditions of the two right-angle boundaries with the circular cavity existing in the space is formulated. Therefore, the total displacement field can be constructed using overlapping principle. An infinite algebraic equations of unknown coefficients existing in the scattering solution field can be gained using multi-polar coordinate and the free stress condition at the boundary of the circular cavity. It can be solved by using limit items in the infinite series which can give a high computation precision. An example is given to illustrate the variations of the tangential stress at the boundary of the circular cavity due to different dimensionless wave numbers, the location of the circular cavity, the loading center and the distributing range of the out-of-plane loading. The results show the efficiency and effectiveness of the method introduced here.     

Dispersion and attenuation in thermoelastic multisize particulate composites

This paper deals with the problem of multiple scattering by a random distribution of spherical solid particles in a solid. The material properties of both media are taken as thermoelastic. The radii of the inclusions may be different. The self-consistent method in its variant of the effective medium is used to find the dispersion and attenuation of quasi-elastic, quasi-thermal and shear waves. The single scattering problem required by this technique is solved approximately by means of the Galerkin method applied to an integral equation using the Green function. Numerical results display a characteristic resonance phenomena which appears in the interval where the results are approximately valid, that is, for very long waves down to wavelengths about twice the largest diameter of the spheres. Examples are shown, for composites with two sets of inclusions, which have either a very similar or dissimilar size. Comparisons are made with the elastic counterpart. Among the material properties, the mass density ratio, inclusion to matrix, seems to play an important and simple role. Frequency intervals are distinguished and shown to depend on that ratio, where the attenuation and dispersion of quasi-elastic and P-waves are either very close to each other or not at all. The same applies to shear waves in either composite. The mass density ratio also displays a simple monotonic decreasing behaviour as a function of the frequency at the first attenuation maximum and velocity minimum. These results may be of interest for the nondestructive testing characterization of particulate composites.     

弹性波散射与多重散射的T矩阵方法

研究弹性波散射与多重散射的T矩阵方法。首先,基于Helmholtz体内和体外公式推导了对应于圆柱型散射体的T矩阵元素的具体表达式;接着分析了在含多个随机分布圆柱型散射体的随机非均匀介质中弹性波的多重散射并给出在统计平均意义下的相干波的定义以及波速和衰减系数计算公式;最后,针对Ge/Al、Sic/Al复合材料用Matlab进行了编程和数值计算;计算单个柱型散射体的散射截面以及随机非均匀介质中相干波的速度和衰减系数,分析了这种介质的频散特性。     

大延伸非均匀介质中地震波全弹性散射理论II-弹性波多次散射理论

推出大延伸非均匀连续介质的弹性波能量传输表达式,并以此为基础建立并发展了相应的弹性动力学能量传递理论.构造了一个基于非均匀薄层或非均匀相屏单次散射迭代法的多次散射模型.该模型既适用于弱散射,也适用于强散驰 既适用于普通散射,也适用于转换散射;尽管高频情况下只考虑普通散射及前向散射.应用该模型计算了弹性多次散射的能通量,处理了散射衰减问题。数值实验的结果表明,短周期地震图上的尾波主要来源于S波散射。     

Experimental study of transition to turbulence of a round thermal plume by ultrasound scattering

In this study, we carried out the characterization of the transition to turbulence of a thermal pure plume by using ultrasound scattering. For this, the position, amplitude and broadening of the scattering peak are analyzed. The technique is based upon the scattering of an ultrasound wave coupling with an unstable flow. The coupling between the acoustic mode with both vorticity and entropy modes is derived from non-linear terms of Navier–Stokes and energy equations. When the scattering mechanism occurs, the characteristic length scale of the flow structure under observation is comparable with the wavelength of incoming sound. Thus, the flow can be probed at different length scales by only changing the frequency of incoming sound. The thermal plume rises from a heated disk immersed into a quiescent medium and can reach transition and fully turbulent regimes. Criteria allowing the identification of both the beginning and the end of transition are derived from the results. The characteristics of the scattering process show evidence that allows us to discern the beginning of transition. The analysis of the amplitude of the scattering peak revealed a homogeneous behavior and led us to think of a possible principle of similarity. The evolution of both thermal and velocity fluctuations has made it possible to establish the limits of both the beginning and the end of transition, in terms of local Grashof number Gr_z and position of the measurement zone z/D. The limits for transition reported in this work are comparable in its magnitude order with those of the literature. It was verified that thermal and velocity transition are phenomena that begin and finish almost simultaneously.     

Particle size characterization by ultrasonic attenuation spectra

This paper contributes to extracting information from signals of broadband ultrasonic attenuation spectrum for effective utilization in particle size characterization. The single particle scattering model and the coupled-phase model are formulated simultaneously, the relationship between particle size distribution and ultrasonic spectrum is established, and a convergence criterion for calculation is quantified. Demonsa'ation inversion by the optimum regularization factor method is carded out to yield typical numerical results for discussion. With the experimental set-up developed by the Institute of Particle and Two-Phase Flow Measurement (IPTFM) at the University of Shanghai for Science and Technology, sand sediment particle size is measured by attenuation spectrum and analyzed using the above inversion algorithm and theoretical models. To validate the proposed ultrasonic spectrum particle sizing method, results are compared with those obtained by microscopy.     

SH波作用下地表软覆盖层中圆形夹杂的动应力分析

利用复变函数法和波函数展开法, 对地表软覆盖层中浅埋圆形夹杂在稳态SH波作用下的动应力集中问题进行研究并给出了解析解。根据SH波散射时的衰减特性, 采用了大圆弧假定的方法, 将半空间覆盖层直线边界问题转化为曲面边界问题。通过算例分析了SH波垂直入射时, 不同入射波波数和圆夹杂与半空间的波数比对圆形夹杂周边动应力集中因子的分布和动应力集中因子最大值变化的影响。算例表明, 圆形夹杂越“软”, 其波数越大, 夹杂周边的动应力集中因子越大; 入射波波数约0.35时, 夹杂周边的最大动应力集中因子达到最大值。     

Ultrasound scattering from a turbulent round thermal pure plume

This research work brings about additional contribution to validate the ultrasound scattering technique as a nonintrusive probe in the Fourier space for measurements performed in unsteady flows. In particular, this work reports experimental evidence of scattering from a turbulent thermal plume utilized as a testing flow. This technique is based upon the scattering of an ultrasound wave hitting and interacting with an unstable flow. The coupling among the acoustic mode with vorticity and entropy modes is derived from nonlinear terms of Navier–Stokes and energy equations. Scattering mechanism occurs when characteristic length scales of flows are comparable with wavelength of sound. Thus, it is possible to probe the flow at different length scales by changing the incoming frequency. The results allow verifying some theoretical predictions, such as the existence of a nonscattering angle. It was also observed, that both the phase and the Doppler shift of the Fourier's signal are linear, respectively, with respect to the time and the frequency of the incident wave. The Doppler shift allowed us to determine the advection velocity and has proved to be sensitive to the direction of the wave vector, to the scattering angle and also, we show that it is possible to have both positive and negative angles. The advection velocity increases with temperature and its values are coherent with those obtained with traditional techniques. Broadening and Doppler shift of the scattering signal allowed us to define the turbulence intensity, whose values are in agreement with those found in thermal plumes, where well-known techniques are currently used. This study has shown that the turbulence intensity increases weakly with temperature, nevertheless it seems more sensitive to the size of the structure under observation.     

二维直角平面内固定圆形夹杂对稳态入射反平面剪切波的散射

利用复变函数法、多极坐标及傅立叶级数展开技术求解了二维直角平面内固定圆形夹杂对稳态入射反平面剪切（shearing horizontal, SH）波的散射问题。首先构造出介质内不存在夹杂时的入射波场和反射波场,然后建立介质内存在夹杂时由夹杂边界产生的能够自动满足直角边应力自由条件的散射波解，从而利用叠加原理写出介质内的总波场。利用夹杂边界处位移条件和傅立叶级数展开方法列出求解散射波中未知系数的无穷代数方程组，在满足计算精度的前提下通过有限项截断，得到相应有限代数方程组的解，最后通过算例具体讨论了二维直角平面水平边界点的位移幅度比和相位随量纲一波数、入射波入射角及夹杂位置的不同而变化的情况，结果表明了算法的有效实用性。