Piecewise Acoustic Source Imaging with Unknown Speed of Sound Using a Level-Set Method

We investigate the following inverse problem:starting from the acoustic wave equation,reconstruct a piecewise constant passive acoustic source from a single boundary temporal measurement without knowing the speed of sound.When the amplitudes of the source are known a priori,we prove a unique determination result of the shape and propose a level set algorithm to reconstruct the singularities.When the singularities of the source are known a priori,we show unique determination of the source amplitudes and propose a least-squares fitting algorithm to recover the source amplitudes.The analysis bridges the low-frequency source inversion problem and the inverse problem of gravimetry.The proposed algorithms are validated and quantitatively evaluated with numerical experiments in 2D and 3D.

A joint absorbing boundary for the multiple-relaxation-time lattice Boltzmann method in seismic acoustic wavefield modeling

Conventional seismic wave forward simulation generally uses mathematical means to solve the macroscopic wave equation,and then obtains the corresponding seismic wavefield.Usually,when the subsurface structure is finely constructed and the continuity of media is poor,this strategy is difficult to meet the requirements of accurate wavefield calculation.This paper uses the multiple-relaxation-time lattice Boltzmann method(MRT-LBM)to conduct the seismic acoustic wavefield simulation and verify its computational accuracy.To cope with the problem of severe reflections at the truncated boundaries,we analogize the viscous absorbing boundary and perfectly matched layer(PML)absorbing boundary based on the single-relaxation-time lattice Boltzmann(SRT-LB)equation to the MRT-LB equation,and further,propose a joint absorbing boundary through comparative analysis.We give the specific forms of the modified MRT-LB equation loaded with the joint absorbing boundary in the two-dimensional(2D)and three-dimensional(3D)cases,respectively.Then,we verify the effects of this absorbing boundary scheme on a 2D homogeneous model,2D modified British Petroleum(BP)gas-cloud model,and 3D homogeneous model,respectively.The results reveal that by comparing with the viscous absorbing boundary and PML absorbing boundary,the joint absorbing boundary has the best absorption performance,although it is a little bit complicated.Therefore,this joint absorbing boundary better solves the problem of truncated boundary reflections of MRT-LBM in simulating seismic acoustic wavefields,which is pivotal to its wide application in the field of exploration seismology.

A Modified Formulation of Singular Boundary Method for Exterior Acoustics

This paper proposes amodified formulation of the singular boundarymethod(SBM)by introducing the combined Helmholtz integral equation formulation(CHIEF)and the self-regularization technique to exterior acoustics.In the SBM,the concept of the origin intensity factor(OIF)is introduced to avoid the singularities of the fundamental solutions.The SBM belongs to the meshless boundary collocation methods.The additional use of the CHIEF scheme and the self-regularization technique in the SBM guarantees the unique solution of the exterior acoustics accurately and efficiently.Consequently,by using the SBM coupled with the CHIEF scheme and the self-regularization technique,the accuracy of the numerical solution can be improved,especially near the corresponding internal characteristic frequencies.Several numerical examples of two-dimensional and threedimensional benchmark examples about exterior acoustics are used to verify the effectiveness and accuracy of the proposed method.The proposed numerical results are compared with the analytical solutions and the solutions obtained by the other numerical methods.

Topology Optimization of Sound-Absorbing Materials for Two-Dimensional Acoustic Problems Using Isogeometric Boundary Element Method

In this work,an acoustic topology optimizationmethod for structural surface design covered by porous materials is proposed.The analysis of acoustic problems is performed using the isogeometric boundary elementmethod.Taking the element density of porousmaterials as the design variable,the volume of porousmaterials as the constraint,and the minimum sound pressure or maximum scattered sound power as the design goal,the topology optimization is carried out by solid isotropic material with penalization(SIMP)method.To get a limpid 0–1 distribution,a smoothing Heaviside-like function is proposed.To obtain the gradient value of the objective function,a sensitivity analysis method based on the adjoint variable method(AVM)is proposed.To find the optimal solution,the optimization problems are solved by the method of moving asymptotes(MMA)based on gradient information.Numerical examples verify the effectiveness of the proposed topology optimization method in the optimization process of two-dimensional acoustic problems.Furthermore,the optimal distribution of sound-absorbingmaterials is highly frequency-dependent and usually needs to be performed within a frequency band.

A nonparametric spectrum estimation method for dispersion and attenuation analysis of borehole acoustic measurements

Dispersion and attenuation analysis can be used to determine formation anisotropy induced by fractures,or stresses.In this paper,we propose a nonparametric spectrum estimation method to get phase dispersion characteristics and attenuation coefficient.By designing an appropriate vector filter,phase velocity,attenuation coefficient and amplitude can be inverted from the waveform recorded by the receiver array.Performance analysis of this algorithm is compared with Extended Prony Method(EPM)and Forward and Backward Matrix Pencil(FBMP)method.Based on the analysis results,the proposed method is capable of achieving high resolution and precision as the parametric spectrum estimation methods.At the meantime,it also keeps high stability as the other nonparametric spectrum estimation methods.At last,applications to synthetic waveforms modeled using finite difference method and real data show its efficiency.The real data processing results show that the P-wave attenuation log is more sensitive to oil formation compared to S-wave;and the S-wave attenuation log is more sensitive to shale formation compared to P-wave.

Three-dimensional acoustic propagation model for shallow waters based on an indirect boundary element method

This work has a two-fold purpose.On the one hand,the theoretical formulation of a three-dimensional(3D)acoustic propagation model for shallow waters with a constant sound speed is presented,based on the boundary element method(BEM),which uses a half-space Green function instead of the more conventional free-space Green function.On the other hand,a numerical implementation is illustrated to explore the formulation in simple idealized cases,controlled by a few parameters,which provides necessary tests for the accuracy and performance of the model.The half-space Green's function,which has been previously used in scattering and diffraction,adds terms to the usual expressions of the integral operators without altering their continuity properties.Verifications against the wavenumber integration solution of the Pekeris waveguide suggest that the model allows an adequate prediction for the acoustic field.Likewise,numerical experiments in relation to the necessary mesh size for the description of the water-marine sediment interface lead to the conclusion that a transmission loss prediction with acceptable accuracy can be obtained with the use of a limited mesh around the desired evaluation region.

Three-dimensional acoustic radiation force of a eukaryotic cell arbitrarily positioned in a Gaussian beam

Expressions are derived for calculating the three-dimensional acoustic radiation force(ARF)on a multilayer microsphere positioned arbitrarily in a Gaussian beam.A theoretical model of a three-layer microsphere with a cell membrane,cytoplasm,and nucleus is established to study how particle geometry and position affect the three-dimensional ARF,and its results agree well with finite-element numerical results.The microsphere can be moved relative to the beam axis by changing its structure and position in the beam,and the axial ARF increases with increasing outer-shell thickness and core size.This study offers a theoretical foundation for selecting suitable parameters for manipulating a three-layer microsphere in a Gaussian beam.

A NOVEL SEMI-ANALYTICAL METHOD FOR SOLVING ACOUSTIC RADIATION FROM LONGITUDINALLY STIFFENED INFINITE NON-CIRCULAR CYLINDRICAL SHELLS IN WATER

Based on the extended homogeneous capacity high precision integration method and the spectrum method of virtual boundary with a complex radius vector, a novel semi-analytical method, which has satisfactory computation efectiveness and precision, is presented for solving the acoustic radiation from a submerged infnite non-circular cylindrical shell stifened by longitudinal ribs by means of the Fourier integral transformation and stationary phase method. In this work, besides the normal interacting force, which is commonly adopted by some researchers, the other interacting forces and moments between the longitudinal ribs and the non-circular cylindrical shell are considered at the same time. The efects of the number and the size of the cross-section of longitudinal ribs on the characteristics of acoustic radiation are investigated. Numerical results show that the method proposed is more efcient than the existing mixed FE-BE method.

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.

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.

Application of dual reciprocity boundary element method to predict acoustic attenuation characteristics of marine engine exhaust silencers

In marine engine exhaust silencing systems, the presence of exhaust gas flow influences the sound propagation inside the systems and the acoustic attenuation performance of silencers. In order to investigate the effects of three-dimensional gas flow and acoustic damping on the acoustic attenuation characteristics of marine engine exhaust silencers, a dual reciprocity boundary element method (DRBEM) was developed. The acoustic governing equation in three-dimensional potential flow was derived first, and then the DRBEM numerical procedure is given. Compared to the conventional boundary element method (CBEM), the DRBEM considers the second order terms of flow Mach number in the acoustic governing equation, so it is suitable for the cases with higher Mach number subsonic flow. For complex exhaust silencers, it is difficult to apply the single-domain boundary element method, so a substructure approach based on the dual reciprocity boundary element method is presented. The experiments for measuring transmission loss of silencers are conducted, and the experimental setup and measurements are explained. The transmission loss of a single expansion chamber silencer with extended inlet and outlet were predicted by DRBEM and compared with the measurements. The good agreements between predictions and measurements are observed, which demonstrated that the derived acoustic governing equation and the DRBEM numerical procedure in the present study are correct.

PREDICTION OF ACOUSTIC PERFORMANCE IN EXPANSION CHAMBER MUFFLERS WITH MEAN FLOW BY FINITE ELEMENT METHOD

The equation of wave propagation in a circular chamber with mean flow is obtained. Computational solution based on finite element method is employed to determine the transmission loss of expansive chamber. The effect of the mean flow and geometry (length of expansion chamber and expansion ratio)on acoustic attenuation performance is discussed, the predicted values of transmission loss of expansion chamber without and with mean flow are compared with those reported in the literature and they agree well. The accuracy of the prediction of transmission loss implies that finite element approximations are applicable to a lot of practical applications.

Diagonal form fast multipole boundary element method for 2D acoustic problems based on Burton-Miller boundary integral equation formulation and its applications

This paper describes formulation and implementation of the fast multipole boundary element method （FMBEM） for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the FMBEM are described in details. They are moment calculation, moment to moment translation, moment to local translation, and local to local translation. A data structure for the quad-tree construction is proposed which can facilitate implementation. An analytical moment expression is derived, which is more accurate, stable, and efficient than direct numerical computation. Numerical examples are presented to demonstrate the accuracy and efficiency of the FMBEM, and radiation of a 2D vibration rail mode is simulated using the FMBEM.

Circuit-field coupled finite element analysis method for an electromagnetic acoustic transducer under pulsed voltage excitation

This paper presents an analytical method for electromagnetic acoustic transducers （EMATs） under voltage excitation and considers the non-uniform distribution of the biased magnetic field. A complete model of EMATs including the non-uniform biased magnetic field, a pulsed eddy current field and the acoustic field is built up. The pulsed voltage excitation is transformed to the frequency domain by fast Fourier transformation （FFT）. In terms of the time harmonic field equations of the EMAT system, the impedances of the coils under different frequencies are calculated according to the circuit-field coupling method and Poynting＇s theorem. Then the currents under different frequencies are calculated according to Ohm＇s law and the pulsed current excitation is obtained by inverse fast Fourier transformation （IFFT）. Lastly, the sequentially coupled finite element method （FEM） is used to calculate the Lorentz force in the EMATs under the current excitation. An actual EMAT with a two-layer two-bundle printed circuit board （PCB） coil, a rectangular permanent magnet and an aluminium specimen is analysed. The coil impedances and the pulsed current are calculated and compared with the experimental results. Their agreement verified the validity of the proposed method. Furthermore, the influences of lift-off distances and the non-uniform static magnetic field on the Lorentz force under pulsed voltage excitation are studied.

Software for Acoustic Design

We present the СATEС software, which implements the solution to the problems of computational acoustics. The software is based on the use of the superelement method and finite element modeling algorithms, in-cluding hydrodynamic noise. The paper presents the main possibilities of software for solving acoustic design problems. .

Calculation and Analysis of Acoustic Characteristics of Straight-Through Perforated Pipe Muffler Based on Multilayer Sound Absorbing Material

Using the multi-physical field simulation software COMSOL,the acoustic characteristics of the multilayer sound absorbing material straight-through perforated pipe muffler are studied by the finite element method.The results show that the finite element calculation of the multilayer sound absorbing material straight-through the perforated pipe muffler agrees well with the experimental measurement results.The reliability of the finite element method for studying the acoustic performance of the straight-through perforated pipe muffler with multilayer sound absorbing materials is shown.Furthermore,the influence of some structural parameters of porous sound absorbing material and micro-perforated plate on the acoustic performance of the multilayer sound absorbing material straight-through perforated pipe muffler is analyzed.The muffler based on multilayer sound absorbing material is different from the traditional muffler.After applying the multilayer sound absorbing material to the straight-through perforated pipe muffler,the transmission loss value greatly increases,which provides new ideas and directions for future research on the muffler.

Mechanical responses and acoustic emission behaviors of coal under compressive differential cyclic loading(DCL):a numerical study via 3D heterogeneous particle model

The stability of coal walls(pillars)can be seriously undermined by diverse in-situ dynamic disturbances.Based on a 3D par-ticle model,this work strives to numerically replicate the major mechanical responses and acoustic emission(AE)behaviors of coal samples under multi-stage compressive cyclic loading with different loading and unloading rates,which is termed differential cyclic loading(DCL).A Weibull-distribution-based model with heterogeneous bond strengths is constructed by both considering the stress-strain relations and AE parameters.Six previously loaded samples were respectively grouped to indicate two DCL regimes,the damage mechanisms for the two groups are explicitly characterized via the time-stress-dependent variation of bond size multiplier,and it is found the two regimes correlate with distinct damage patterns,which involves the competition between stiffness hardening and softening.The numerical b-value is calculated based on the mag-nitudes of AE energy,the results show that both stress level and bond radius multiplier can impact the numerical b-value.The proposed numerical model succeeds in replicating the stress-strain relations of lab data as well as the elastic-after effect in DCL tests.The effect of damping on energy dissipation and phase shift in numerical model is summarized.

The Application of the Acoustic Method to the Exploration of Urban Active Faults and a Comparison with the Artificial Earthquake, and Radar Methods

The depth of upper fault point is the key data for ascertaining the active age of a buried fault on a plain. The difference of depth obtained from same fault may be dozens to several hundred meters when using different geophysical methods. It can result in the absolutely opposite conclusions when judging fault activity. Because of a lack of an artificial earthquake source with wide band and high central frequency, many kinds of methods have to be used together. The higher the frequency of the artificial earthquake wave, electromagnetic wave and sonic wave, the higher the resolution. However the attenuation is also very fast and the exploration depth is very shallow. The reverse is also true. The frequency of artificial seismic waves is in the tens of Hz. Its exploration depth is big and the resolution is poor. The frequency of radar electromagnetic waves is about a million Hz, indicating that the resolving power is better, but the exploration depth is very shallow. However, the acoustic frequency is thousands of Hz, its resolving power is better than that of the artificial earthquake method and the exploration depth is larger than that of the radar method. So it is suitable for extra shallow exploration in the thick deposit strata of the Quaternary. The preliminary results detected using the high frequency acoustic method in extra shallow layers indicates that previous inferences about some fault activity in the eastern part of the North China plain may need to be greatly corrected.

RLS and LMS blind adaptive multi-user detection method and comparison in acoustic communication

RLS and LMS blind adaptive multi-user detection algorithm and multi-user detector was proposed to solve the problem of multi-user signal detection problem encountered in underwater acoustic communication networks.In simulation analysis,RLS and the LMS blind adaptive multi-user detector were designed and tested for synchronous and asynchronous multi-user communication process.The results of SIR comparison and MMSE comparison show that,both of the two methods can realize blind adaptive detection when any user change in multi-user communication,during this process,the training communication sequences are not needed.The RLS algorithm has about 5 dB higher in SIR compared with LMS algorithm,and the convergence velocity of RLS algorithm is also higher than LMS algorithm when the communication users change.RLS algorithm has better ability in multi-user detection than that of LMS algorithm,and it has great attraction and guiding significance for solving the problem of multiple access interference(MAI) in multi-user communication.

Novel damage constitutive models and new quantitative identification method for stress thresholds of rocks under uniaxial compression

Four key stress thresholds exist in the compression process of rocks,i.e.,crack closure stress(σ_(cc)),crack initiation stress(σ_(ci)),crack damage stress(σ_(cd))and compressive strength(σ_(c)).The quantitative identifications of the first three stress thresholds are of great significance for characterizing the microcrack growth and damage evolution of rocks under compression.In this paper,a new method based on damage constitutive model is proposed to quantitatively measure the stress thresholds of rocks.Firstly,two different damage constitutive models were constructed based on acoustic emission(AE)counts and Weibull distribution function considering the compaction stages of the rock and the bearing capacity of the damage element.Then,the accumulative AE counts method(ACLM),AE count rate method(CRM)and constitutive model method(CMM)were introduced to determine the stress thresholds of rocks.Finally,the stress thresholds of 9 different rocks were identified by ACLM,CRM,and CMM.The results show that the theoretical stress−strain curves obtained from the two damage constitutive models are in good agreement with that of the experimental data,and the differences between the two damage constitutive models mainly come from the evolutionary differences of the damage variables.The results of the stress thresholds identified by the CMM are in good agreement with those identified by the AE methods,i.e.,ACLM and CRM.Therefore,the proposed CMM can be used to determine the stress thresholds of rocks.