The viscous strip approach to simplify the calculation of the surface acoustic wave generated streaming

In recent decades,the importance of surface acoustic waves,as a biocompatible tool to integrate with microfluidics,has been proven in various medical and biological applications.The numerical modeling of acoustic streaming caused by surface acoustic waves in microchannels requires the effect of viscosity to be considered in the equations which complicates the solution.In this paper,it is shown that the major contribution of viscosity and the horizontal component of actuation is concentrated in a narrow region alongside the actuation boundary.Since the inviscid equations are considerably easier to solve,a division into the viscous and inviscid domains would alleviate the computational load significantly.The particles'traces calculated by this approximation are excellently alongside their counterparts from the completely viscous model.It is also shown that the optimum thickness for the viscous strip is about 9-fold the acoustic boundary layer thickness for various flow patterns and amplitudes of actuation.

Ion acoustic solitary waves in an adiabatic dusty plasma:Roles of superthermal electrons,ion loss and ionization

We investigate propagation of dust ion acoustic solitary wave(DIASW)in a multicomponent dusty plasma with adiabatic ions,superthermal electrons,and stationary dust.The reductive perturbation method is employed to derive the damped Korteweg-de Vries(DKdV)equation which describes DIASW.The result reveals that the adiabaticity of ions significantly modifies the basic features of the DIASW.The ionization effect makes the solitary wave grow,while collisions reduce the growth rate and even lead to the damping.With the increases in ionization cross sectionΔσ/σ_(0),ion-to-electron density ratioδ_(ie)and superthermal electrons parameterκ,the effect of ionization on DIASW enhances.

Numerical Simulation and Parallel Computing of Acoustic Wave Equation in Isotropic-Heterogeneous Media

In this paper,we consider the numerical implementation of the 2D wave equation in isotropic-heterogeneous media.The stability analysis of the scheme using the von Neumann stability method has been studied.We conducted a study on modeling the propagation of acoustic waves in a heterogeneous medium and performed numerical simulations in various heterogeneous media at different time steps.Developed parallel code using Compute Unified Device Architecture(CUDA)technology and tested on domains of various sizes.Performance analysis showed that our parallel approach showed significant speedup compared to sequential code on the Central Processing Unit(CPU).The proposed parallel visualization simulator can be an important tool for numerous wave control systems in engineering practice.

3D near-surface P-wave velocity structure imaging with Distributed Acoustic Sensing and electric hammer source

Distributed Acoustic Sensing(DAS) is an emerging technique for ultra-dense seismic observation, which provides a new method for high-resolution sub-surface seismic imaging. Recently a large number of linear DAS arrays have been used for two-dimensional S-wave near-surface imaging in urban areas. In order to explore the feasibility of three-dimensional(3D) structure imaging using a DAS array, we carried out an active source experiment at the Beijing National Earth Observatory. We deployed a 1 km optical cable in a rectangular shape, and the optical cable was recast into 250 sensors with a channel spacing of 4 m. The DAS array clearly recorded the P, S and surface waves generated by a hammer source. The first-arrival P wave travel times were first picked with a ShortTerm Average/Long-Term Average(STA/LTA) method and further manually checked. The P-wave signals recorded by the DAS are consistent with those recorded by the horizontal components of short-period seismometers. At shorter source-receiver distances, the picked P-wave arrivals from the DAS recording are consistent with vertical component recordings of seismometers, but they clearly lag behind the latter at greater distances.This is likely due to a combination of the signal-to-noise ratio and the polarization of the incoming wave. Then,we used the Tomo DD software to invert the 3D P-wave velocity structure for the uppermost 50 m with a resolution of 10 m. The inverted P-wave velocity structures agree well with the S-wave velocity structure previously obtained through ambient noise tomography. Our study indicates the feasibility of 3D near-surface imaging with the active source and DAS array. However, the inverted absolute velocity values at large depths may be biased due to potential time shifts between the DAS recording and seismometer at large source-receiver distances.

A stable staggered-grid finite-difference scheme for acoustic modeling beyond conventional stability limit

Staggered-grid finite-difference(SGFD)schemes have been widely used in acoustic wave modeling for geophysical problems.Many improved methods are proposed to enhance the accuracy of numerical modeling.However,these methods are inevitably limited by the maximum Courant-Friedrichs-Lewy(CFL)numbers,making them unstable when modeling with large time sampling intervals or small grid spacings.To solve this problem,we extend a stable SGFD scheme by controlling SGFD dispersion relations and maximizing the maximum CFL numbers.First,to improve modeling stability,we minimize the error between the FD dispersion relation and the exact relation in the given wave-number region,and make the FD dispersion approach a given function outside the given wave-number area,thus breaking the conventional limits of the maximum CFL number.Second,to obtain high modeling accuracy,we use the SGFD scheme based on the Remez algorithm to compute the FD coefficients.In addition,the hybrid absorbing boundary condition is adopted to suppress boundary reflections and we find a suitable weighting coefficient for the proposed scheme.Theoretical derivation and numerical modeling demonstrate that the proposed scheme can maintain high accuracy in the modeling process and the value of the maximum CFL number of the proposed scheme can exceed that of the conventional SGFD scheme when adopting a small maximum effective wavenumber,indicating that the proposed scheme improves stability during the modeling.

Acoustic propagation uncertainty in internal wave environments using an ocean-acoustic joint model

An ocean-acoustic joint model is developed for research of acoustic propagation uncertainty in internal wave environments.The internal waves are numerically produced by tidal forcing over a continental slope using an ocean model.Three parameters(i.e.,internal wave,source depth,and water depth)contribute to the dynamic waveguide environments,and result in stochastic sound fields.The sensitivity of the transmission loss(TL)to environment parameters,statistical characteristics of the TL variation,and the associated physical mechanisms are investigated by the Sobol sensitivity analysis method,the Monte Carlo sampling,and the coupled normal mode theory,respectively.The results show that the TL is most sensitive to the source depth in the near field,resulted from the initial amplitudes of higher-order modes;while in middle and far fields,the internal waves are responsible for more than 80%of the total acoustic propagation contribution.In addition,the standard deviation of the TL in the near field and the shallow layer is smaller than those in the middle and far fields and the deep layer.

Response characteristics of drill-string guided wave in downhole acoustic telemetry

Modeling of a drill-string acoustic channel has been an important topic in downhole telemetry for a long time.The propagation of drill-string guided waves in the borehole contains excitation,attenuation,and mode conversion issues that have not been considered by existing modeling methods.In this article,we formulate a hybrid modeling method to investigate the response characteristics of a fundamental-mode drill-string wave in various borehole environments.This hybrid method provides channel functions,including transmitting and receiving deployments,periodicity of the structure,and formation property changes.The essential physics of the drill-string wave propagation is captured with a one-dimensional model.The analytical solutions of the wavefield in multilayered cylindrical structures are introduced into a propagation matrix to express drill-string-wave interactions with the borehole environments.The effectiveness of the proposed method is confirmed through comparison with the finite-difference method.In addition,by designing numerical models,we investigate the conversion effect of the drill-string wave at the tool joint.We demonstrate that the conversion intensity of the drill-string wave is positively correlated not only with the cross-sectional area of the tool joint but also with the wave impedance of the outer formation.Hard formation outside the borehole reduces the energy leakage while intensifying the conversion of drill-string waves to Stoneley waves,and the opposite is true for the drill string in an infinite fluid.The converted Stoneley waves interfere with the drill-string waves,resulting in variations of bandgap distribution,which challenges the reliability of the data transmission.

The(1+1)-dimensional nonlinear ion acoustic waves in multicomponent plasma containing kappa electrons

Large amplitude (1+1)-dimensional nonlinear ion acoustic waves are theoretically studied in multicomponent plasma consisting of positively charged ions and negatively charged ions, ion beam, kappa-distributed electrons, and dust grains,respectively. By using the Sagdeev potential method, the dynamical system and the Sagdeev potential function are obtained.The important influences of system parameters on the phase diagram of this system are investigated. It is found that the linear waves, the nonlinear waves and the solitary waves are coexistent in the multicomponent plasma system. Meanwhile,the variations of Sagdeev potential with parameter can also be obtained. Finally, it seems that the propagating characteristics of (1+1)-dimensional nonlinear ion acoustic solitary waves and ion acoustic nonlinear shock wave can be influenced by different parameters of this system.

Three-dimensional acoustic wave equation modeling based on the optimal finite-difference scheme

Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.

Implications for identification of principal stress directions from acoustic emission characteristics of granite under biaxial compression experiments

The rock fracture characteristics and principal stress directions are crucial for prevention of geological disasters.In this study,we carried out biaxial compression tests on cubic granite samples of 100 mm in side length with different intermediate principal stress gradients in combination with acoustic emission(AE)technique.Results show that the fracture characteristics of granite samples change from‘sudden and aggregated’to‘continuous and dispersed’with the increase of the intermediate principal stress.The effect of increasing intermediate principal stress on AE amplitude is not significant,but it increases the proportions of high-frequency AE signals and shear cracks,which in turn increases the possibility of unstable rock failure.The difference of stress in different directions causes the anisotropy of rock fracture and thus leads to the obvious anisotropic characteristics of wave velocity variations.The anisotropy of wave velocity variations with stress difference is probable to identify the principal stress directions.The AE characteristics and the anisotropy of wave velocity variations of granite under two-dimensional stress are not only beneficial complements for rock fracture characteristic and principal stress direction identification,but also can provide a new analysis method for stability monitoring in practical rock engineering.

Excitation and Detection of Shear Horizontal Waves with Electromagnetic Acoustic Transducers for Nondestructive Testing of Plates

The fundamental shear horizontal（SH0） wave has several unique features that are attractive for long-range nondestructive testing（NDT）. By a careful design of the geometric configuration, electromagnetic acoustic transducers（EMATs） have the capability to generate a wide range of guided wave modes, such as Lamb waves and shear-horizontal（SH） waves in plates. However, the performance of EMATs is influenced by their parameters. To evaluate the performance of periodic permanent magnet（PPM） EMATs, a distributed-line-source model is developed to calculate the angular acoustic field cross-section in the far-field. Numerical analysis is conducted to investigate the performance of such EMATs with different geometric parameters, such as period and number of magnet arrays, and inner and outer coil widths. Such parameters have a great influence on the directivity of the generated SH0 waves that arises mainly in the amplitude and width of both main and side lobes. According to the numerical analysis, these parameters are optimized to obtain better directivity. Optimized PPM EMATs are designed and used for NDT of strip plates. Experimental results show that the lateral boundary of the strip plate has no perceivable influence on SHO-wave propagation, thus validating their used in NDT. The proposed model predicts the radiation pattern ofPPM EMATs, and can be used for their parameter optimization.

Fabrication and characterization of one-port surface acoustic wave resonators on semi-insulating GaN substrates

One-port surface acoustic wave resonators(SAWRs) are fabricated on semi-insulating high-quality bulk GaN and GaN film substrates, respectively. The semi-insulating GaN substrates are grown by hydride vapor phase epitaxy(HVPE)and doped with Fe. The anisotropy of Rayleigh propagation and the electromechanical coupling coefficient in Fe-doped GaN are investigated. The difference in resonance frequency between the SAWs between [1120] GaN and [1100] GaN is about 0.25% for the Rayleigh propagation mode, which is smaller than that of non-intentionally doped GaN film(~1%)reported in the literature. The electromechanical coupling coefficient of Fe-doped GaN is about 3.03%, which is higher than that of non-intentionally doped GaN film. The one-port SAWR fabricated on an 8-μm Fe-doped GaN/sapphire substrate has a quality factor of 2050, and that fabricated on Fe-doped bulk GaN has a quality factor as high as 3650. All of our results indicate that high-quality bulk GaN is a very promising material for application in SAW devices.

Imaging shallow structure with active-source surface wave signal recorded by distributed acoustic sensing arrays

Distributed acoustic sensing(DAS) is one recently developed seismic acquisition technique that is based on fiber-optic sensing. DAS provides dense spatial spacing that is useful to image shallow structure with surface waves.To test the feasibility of DAS in shallow structure imaging,the PoroTomo team conducted a DAS experiment with the vibroseis truck T-Rex in Brady’s Hot Springs, Nevada, USA.The Rayleigh waves excited by the vertical mode of the vibroseis truck were analyzed with the Multichannel Analysis of Surface Waves(MASW) method. Phase velocities between5 and 20 Hz were successfully extracted for one segment of cable and were employed to build a shear-wave velocity model for the top 50 meters. The dispersion curves obtained with DAS agree well with the ones extracted from co-located geophones data and from the passive source Noise Correlation Functions(NCF). Comparing to the co-located geophone array, the higher sensor density that DAS arrays provides help reducing aliasing in dispersion analysis, and separating different surface wave modes. This study demonstrates the feasibility and advantage of DAS in imaging shallow structure with surface waves.

Acoustic wave detection of laser shock peening

In order to overcome the existing disadvantages of offline laser shock peening detection methods, an online detection method based on acoustic wave signals energy is provided. During the laser shock peening, an acoustic emission sen- sor at a defined position is used to collect the acoustic wave signals that propagate in the air. The acoustic wave signal is sampled, stored, digitally filtered and analyzed by the online laser shock peening detection system. Then the system gets the acoustic wave signal energy to measure the quality of the laser shock peening by establishing the correspondence between the acoustic wave signal energy and the laser pulse energy. The surface residual stresses of the samples are measured by X-ray stress analysis instrument to verify the reliability. The results show that both the surface residual stress and acoustic wave signal energy are increased with the laser pulse energy, and their growth trends are consistent. Finally, the empirical formula between the surface residual stress and the acoustic wave signal energy is established by the cubic equation fitting, which will provide a theoretical basis for the real-time online detection of laser shock peening.

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 UEberall H, 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.

PROPAGATION OF SURFACE ACOUSTIC WAVES IN PRESTRESSED ANISOTROPIC LAYERED PIEZOELECTRIC STRUCTURES

The propagation of surface acoustic waves in layered piezoelectric structureswith initial stresses is investigated. The phase velocity equations are obtained for electricallyfree and shorted cases, respectively. Effects of the initial stresses on the phase velocity and theelectromechanical coupling coefficient for the fundamental mode of the layered piezoelectricstructures are discussed. Numerical results for the c-axis oriented film of LiNbO_3 on a sapphiresubstrate are given. It is found that the fractional change in phase velocity is a linear functionwith the initial stresses, and the electromechanical coupling factor increases with an increase ofthe absolute values of the compressive initial stresses. The results are useful for the design ofsurface acoustic wave devices.

Lamb wave signal selective enhancement by an improved design of meander-coil electromagnetic acoustic transducer

In this paper, we investigate a method of selectively enhancing the single mode signal of a Lamb wave by using a meander-coil electromagnetic acoustic transducer （EMAT） with a new magnetic configuration. We use the Lamb antisym- metric （A0） mode and symmetric （SO） mode as an example for analysis. The analytical expression of the magnitude of the spatial Fourier transform of the Lorentz force generated by different meander coils is used to determine the optimal driving frequency for single mode generation. The numerical calculation is used to characterize the new magnetic configuration and the conventional EMAT magnet. Experimental examinations of each meander coil in combination with the conventional and new magnetic configuration show that the Lamb wave signal can be selectively enhanced by choosing the appropriate driving frequency and coil parameters through using the improved meander-coil EMAT.

3-D acoustic wave equation forward modeling with topography

In order to model the seismic wave field with surface topography, we present a method of transforming curved grids into rectangular grids in two different coordinate systems. Then the 3D wave equation in the transformed coordinate system is derived. The wave field is modeled using the finite-difference method in the transformed coordinate system. The model calculation shows that this method is able to model the seismic wave field with fluctuating surface topography and achieve good results. Finally, the energy curves of the direct and reflected waves are analyzed to show that surface topography has a great influence on the seismic wave＇s dynamic properties.

Controlling an acoustic wave with a cylindrically-symmetric gradient-index system

We present a detailed theoretical description of wave propagation in an acoustic gradient-index system with cylindrical symmetry and demonstrate its potential to numerically control acoustic waves in different ways.The trajectory of an acoustic wave within the system is derived by employing the theory of geometric acoustics,and the validity of the theoretical descriptions is verified numerically by using the finite element method simulation.The results show that by tailoring the distribution function of the refractive index,the proposed system can yield a tunable manipulation of acoustic waves,such as acoustic bending,trapping,and absorbing.

WAVE SUPERPOSITION METHOD BASED ON VIRTUAL SOURCE BOUNDARY WITH COMPLEX RADIUS VECTOR FOR SOLVING ACOUSTIC RADIATION PROBLEMS

By virtue of the comparability between the wave superposition method and the dynamic analysis of structures, a general format for overcoming the non-uniqueness of solution induced by the wave superposition method at the eigenfrequencies of the corresponding interior problems is proposed. By adding appropriate damp to the virtual source system of the wave superposition method, the unique solutions for all wave numbers can be ensured. Based on this thought, a novel method-wave superposition method with complex radius vector is constructed. Not only is the computational time of this method approximately equal to that of the standard wave superposition method, but also the accuracy is much higher compared with other correlative methods. Finally, by taking the pulsating sphere and oscillating sphere as examples, the results of calculation show that the present method can effectively overcome the non-uniqueness problem.