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 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.

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.

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.

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.

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.

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.

AN EFFICIENT PREDICTION METHOD FOR THE TWO-DIMENSIONAL COUPLED STRUCTURAL-ACOUSTIC ANALYSIS

A wave number method （WNM） is proposed to deal with the two-dimensional coupled structural-acoustic problem. Based on an indirect Trefftz approach, the displacement and the pressure response are approximated respectively by a set of wave functions, which exactly satisfy the governing equations and are independent of the size of the coupled system. The wave functions comprise the exact solutions of the homogeneous part of the governing equations and some particular solution functions, which arise from the external excitation. The weighting coefficients of the wave functions can be obtained by enforcing the pressure approximation to satisfy the boundary conditions and it is performed by applying the weighted residual formulation. The example is computed by the WNM and the BEM. The results show that, the WNM can attain the same accuracy and convergence as the BEM with less degrees of freedom.

A Differential Continuation－Regularization Method for Solving Inverse Problems of Acoustic Wave Equations

ＡＤｉｆｆｅｒｅｎｔｉａｌＣｏｎｔｉｎｕａｔｉｏｎ－ＲｅｇｕｌａｒｉｚａｔｉｏｎＭｅｔｈｏｄｆｏｒＳｏｌｖｉｎｇＩｎｖｅｒｓｅＰｒｏｂｌｅｍｓｏｆＡｃｏｕｓｔｉｃＷａｖｅＥｑｕａｔｉｏｎｓ￥（张大力）（韩波）（刘家琦）（姜立功）ＺＨＡＮＧＤａｌｉ；Ｈ...

Modeling and RTM for arbitrary-order pure acoustic wave equation in VTI media using normalized pseudo-analytical method

The conventional pseudo-acoustic wave equations（PWEs） in vertical transversely isotropic（VTI）media may generate SV-wave artifacts and propagation instabilities when anisotropy parameters cannot satisfy the pseudo-acoustic assumption. One solution to these issues is to use pure acoustic anisotropic wave equations, which can produce stable and pure P-wave responses without any SVwave pollutions. The commonly used pure acoustic wave equations（PAWEs） in VTI media are mainly derived from the decoupled P-SV dispersion relation based on first-order Taylor-series expansion（TE）, thus they will suffer from accuracy loss in strongly anisotropic media. In this paper, we adopt arbitrary-order TE to expand the square root term in Alkhalifah＇s accurate acoustic VTI dispersion relation and solve the corresponding PAWE using the normalized pseudoanalytical method（NPAM） based on optimized pseudodifferential operator. Our analysis of phase velocity errors indicates that the accuracy of our new expression is perfectly acceptable for majority anisotropy parameters. The effectiveness of our proposed scheme also can be demonstrated by several numerical examples and reverse-time migration（RTM） result.

Maximization of Acoustic Levitating Force for a Single-Axis Acoustic Levitation System Using the Finite Element Method

We investigate single-axis acoustic levitation using standing waves to levitate particles freely in a medium bounded by a driver and a reflector. The acoustic pressure at the pressure antinode of the standing wave counteracts the downward gravitational force of the levitating object. The optimal relationship between the air gap and the driving frequency leads to resonance and hence maximization of the levitating force. Slight deviation from the exact resonance condition causes a reduction in acoustic pressure at the pressure antinodes. This results in a significant reduction of the levitating force. The driving frequency is kept constant while the air gap is varied for different conditions. The optimal air gap for maximizing the levitation force is studied for first three resonance modes. Furthermore, a levitating particle is introduced between the driver and the reflector. The dependence of the resonance condition on the size of the levitating particle as well as the position of the particle between the driver and the reflector has also been studied. As the size of the levitating particle increases, the resonance condition also gets modified. Finite element results show a good agreement with the validated results available in the literature. Furthermore, the finite element approach is also used to study the variation of acoustic pressure at the pressure antinode with respect to the size of the reflector. The optimum diameter of the reflector is calculated for maximizing the levitating force for three resonance modes.

Wave Number Method for Three-Dimensional Steady-State Acoustic Problems

Based on the indirect Trefftz approach, a wave number method （WNM） is proposed to deal with three-dimensional steady-state acoustic problems. In the WNM, the dynamic pressure response variable is approximated by a set of wave functions, which exactly satisfy the Helmholtz equation. The set of wave functions comprise the exact solutions of the homogeneous part of the governing equations and some particular solution functions. The unknown coefficients of the wave functions can be obtained by enforcing the pressure approximation to satisfy the boundary conditions. Compared with the boundary element method （BEM）, the WNM have a smaller system matrix, and is applicable to the radiation problems since the wave functions are independent of the domain size. A 3D acoustic cavity is exemplified to show the properties of the method. The results show that the wave number method is more efficient than the BEM, and it is fairly accurate.

An Efficient and Accurate Spectral Method for Acoustic Scattering in Elliptic Domains

An efficient and accurate method for solving the two-dimensional Helmholtz equation in domains exterior to elongated obstacles is developed in this paper.The method is based on the so called transformed field expansion(TFE) coupled with a spectral-Galerkin solver for elliptical domain using Mathieu functions.Numerical results are presented to show the accuracy and stability of the proposed method.

Finite Element Method Computations of the Acoustics of the Human Head Based on the Projection Based Interpolation Data

In this paper we present the Projection Based Interpolation (PBI) technique for construction of continuous approximation of MRI scan data of the human head. We utilize the result of the PBI algorithm to perform three dimensional (3D) Finite Element Method (FEM) simulations of the acoustics of the human head. The computational problem is a multi-physics problem modeled as acoustics coupled with linear elasticity. The computational grid contains tetrahedral finite elements with the number of equations and polynomial orders of approximation varying locally on finite element edges, faces, and interiors. We utilize our own out-of-core parallel direct solver for the solution of this multi-physics problem. The solver minimizes the memory usage by dumping out all local systems from all nodes of the entire elimination tree during the elimination phase.

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.

Suppress numerical dispersion in reversetime migration of acoustic wave equation using optimal nearly analytic discrete method

Using staggered-grid finite difference method to solve seismic wave equation,large spatial grid and high dominant frequency of source cause numerical dispersion,staggeredgrid finite difference method,which can reduce the step spatial size and increase the order of difference,will multiply the calculation amount and reduce the efficiency of solving wave equation.The optimal nearly analytic discrete(ONAD)method can accurately solve the wave equation by using the combination of displacement and gradient of spatial nodes to approach the spatial partial derivative under rough grid and high-frequency condition.In this study,the ONAD method is introduced into the field of reverse-time migration(RTM)for performing forward-and reverse-time extrapolation of a two-dimensional acoustic equation,and the RTM based on ONAD method is realized via normalized cross-correlation imaging condition,effectively suppressed the numerical dispersion and improved the imaging accuracy.Using ONAD method to image the groove model and SEG/EAGE salt dome model by RTM,and comparing with the migration sections obtained by staggered-grid finite difference method with the same time order 2 nd and space order 4 th,results show that the RTM based on ONAD method can effectively suppress numerical dispersion caused by the high frequency components in source and shot records,and archive accurate imaging of complex geological structures especially the fine structure,and the migration sections of the measured data show that ONAD method has practical application value.