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.

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

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

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.

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.

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.

Full-waveform Velocity Inversion Based on the Acoustic Wave Equation

Full-waveform velocity inversion based on the acoustic wave equation in the time domain is investigated in this paper. The inversion is the iterative minimization of the misfit between observed data and synthetic data obtained by a numerical solution of the wave equation. Two inversion algorithms in combination with the CG method and the BFGS method are described respectively. Numerical computations for two models including the benchmark Marmousi model with complex structure are implemented. The inversion results show that the BFGS-based algorithm behaves better in inversion than the CG-based algorithm does. Moreover, the good inversion result for Marmousi model with the BFGS-based algorithm suggests the quasi-Newton methods can provide an important tool for large-scale velocity inversion. More computations demonstrate the correctness and effectives of our inversion algorithms and code.

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.

Space time fractional KdV Burgers equation for dust acoustic shock waves in dusty plasma with non-thermal ions

The KdV-Burgers equation for dust acoustic waves in unmagnetized plasma having electrons, singly charged non- thermal ions, and hot and cold dust species is derived using the reductive perturbation method. The Boltzmann distribution is used for electrons in the presence of the cold （hot） dust viscosity coefficients. The semi-inverse method and Agrawal variational technique are applied to formulate the space-time fractional KdV-Burgers equation which is solved using the fractional sub-equation method. The effect of the fractional parameter on the behavior of the dust acoustic shock waves in the dusty plasma is investigated.

Efficient Computational Approach for Predicting the 3D Acoustic Radiation of the Elastic Structure in Pekeris Waveguides

Ocean boundaries present a significant effect on the vibroacoustic characteristics and sound propagation of an elastic structure in practice.In this study,an efficient finite element/wave superposition method(FE/WSM)for predicting the three-dimen-sional acoustic radiation from an arbitrary-shaped radiator in Pekeris waveguides with a lossy seabed is proposed.The method is based on the FE method(FEM),WSM,and sound propagation models.First,a near-field vibroacoustic model is established by the FEM to obtain vibration information on a radiator surface.Then,the WSM based on the Helmholtz boundary integral is used to pre-dict the far-field acoustic radiation and propagation.Furthermore,the rigorous image source method and complex normal mode are employed to obtain the near-and far-field Green’s function(GF),respectively.The former,which is based on the spherical wave decomposition,is adopted to accurately solve the near-field source strength,and the far-field acoustic radiation is calculated by the latter and perturbation theory.The simulations of both models are compared to theoretical wavenumber integration solutions.Finally,numerical experiments on elastic spherical and cylindrical shells in Pekeris waveguides are presented to validate the accuracy and efficiency of the proposed method.The results show that the FE/WSM is adaptable to complex radiators and ocean-acoustic envi-ronments,and are easy to implement and computationally efficient in calculating the structural vibration,acoustic radiation,and sound propagation of arbitrarily shaped radiators in practical ocean environments.

Oblique collisional effects of dust acoustic waves in unmagnetized dusty plasma

Effects of oblique collisions of the dust acoustic(DA)waves in dusty plasma are studied by considering unmagnetized fully ionized plasma.The plasma consists of inertial warm negatively charged massive dusts,positively charged dusts,superthermal kappa distributed electrons,and isothermal ions.The extended Poincaré–Lighthill–Kuo(e PLK)method is employed for the drivation of two-sided Korteweg–de Vries(KdV)equations(KdVEs).The Kd V soliton solutions are derived by using the hyperbolic secant method.The effects of superthermality index of electrons,temperature ratio of isothermal ion to electron,and the density ratio of isothermal ions to negatively charged massive dusts on nonlinear coefficients are investigated.The effects of oblique collision on amplitude,phase shift,and potential profile of right traveling solitons of DA waves are also studied.The study reveals that the new nonlinear wave structures are produced in the colliding region due to head-on collision of the two counter propagating DA waves.The nonlinearity is found to decrease with the increasing density ratio of ion to negative dust in the critical region.The phase shifts decrease(increase)with increasing the temperature ratio of ion to electron(κe).The hump(compressive,κe<κec)and dipshaped(rarefactive,κe>κec)solitons are produced depending on the angle(θ)of oblique collision between the two waves.

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.

（2＋1）-Dimensional Electron Acoustic Solitary Waves in an Unmagnetized Collisionless Plasma with Vortex-Like Hot Electrons

In this paper, （2＋1）-dimensional electron acoustic waves （EAW） in an unmagnetized collisionless plasma have been studied by the linearized method and the reductive perturbation technique, respectively. The dispersion relation and a modified Kadomtsev-Petviashvili （KP） equation have been obtained for the EAW in the plasma considering a cold electron fluid and a vortex-like hot electrons. It is found from some numerical results that the parameter β（the ratio of the free hot electron temperature to the hot trapped electron temperature） effects on the amplitude and the Width of the electron acoustic solitary waves （EASW）. It can be indicated that the free hot electron temperature and the hot trapped electron temperature have very important effect on the characters of the propagation for the EASW.

Efficient solution of large-scale matrix of acoustic wave equations in 3D frequency domain

In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain finite-difference forward simulation algorithm based on the acoustic wave equation establishes a large bandwidth complex matrix according to the discretized acoustic wave equation,and then the frequency domain wave field value is obtained by solving the matrix equation.In this study,the predecessor's optimized five-point method is extended to a 3D seven-point finite-difference scheme,and then a perfectly matched layer absorbing boundary condition(PML)is added to establish the corresponding matrix equation.In order to solve the complex matrix,we transform it to the equivalent real number domain to expand the solvable range of the matrix,and establish two objective functions to transform the matrix solving problem into an optimization problem that can be solved using gradient methods,and then use conjugate gradient algorithm to solve the problem.Previous studies have shown that in the conjugate gradient algorithm,the product of the matrix and the vector is the main factor that affects the calculation efficiency.Therefore,this study proposes a method that transform bandwidth matrix and vector product problem into some equivalent vector and vector product algorithm,thereby reducing the amount of calculation and storage.

Nonlinear Acoustic Shadow Method to Reduce Reverberation Artifact

A novel technique for reducing reverberation artifact in acoustic shadow imaging using nonlinear ultrasound interaction, called nonlinear acoustic shadow method, has been developed and experimentally studied. In this technique, the conventional acoustic shadow method is modified by using the secondary wave generated by nonlinear interaction of two primary sound waves emitted from parametric array. Either conventional or nonlinear acoustic shadow imaging is carried out for aluminum square cylinder and the size of the shadow is compared. The result shows that the nonlinear acoustic shadow method reduces reverberation artifact inside the square cylinder and has better accuracy in the size measurement than conventional acoustic shadow method.

Nonlinear ion-acoustic solitary waves in an electron-positron-ion plasma with relativistic positron beam

The propagation characteristics of nonlinear ion–acoustic（IA） solitary waves（SWs） are studied in thermal electron–positron–ion plasma considering the effect of relativistic positron beam. Starting from a set of fluid equations and using the reductive perturbation technique, we derive a Korteweg–de Vries（KdV） equation which governs the evolution of weakly nonlinear IA SWs in relativistic beam driven plasmas. The properties of the IA soliton are studied, and it is shown that the presence of relativistic positron beam significantly modifies the characteristics of IA solitons.

Numerical Simulation of Acoustic-Gravity Waves Propagation in a Heterogeneous Earth-Atmosphere Model with Wind in the Atmosphere

A numerical-analytical solution for seismic and acoustic-gravity waves propagation is applied to a heterogeneous “Earth-Atmosphere” model. Seismic wave propagation in an elastic half-space is described by a system of first order dynamic equations of elasticity theory. Propagation of acoustic-gravity waves in the atmosphere is described by the linearized Navier-Stokes equations with the wind. The algorithm proposed is based on the integral Laguerre transform with respect to time, the finite integral Fourier transform along the spatial coordinate with the finite difference solution of the reduced problem.

Ion-acoustic waves in plasma of warm ions and isothermal electrons using time-fractional KdV equation

The ion-acoustic solitary wave in collisionless unmagnetized plasma consisting of warm ions-fluid and isothermal electrons is studied using the time fractional KdV equation. The reductive perturbation method has been employed to derive the Korteweg-de Vries equation for small but finite amplitude ion-acoustic wave in warm plasma. The Lagrangian of the time fractional KdV equation is used in a similar form to the Lagrangian of the regular KdV equation with fractional derivative for the time differentiation. The variation of the functional of this Lagrangian leads to the Euler-Lagrange equation that gives the time fractional KdV equation. The variational-iteration method is used to solve the derived time fractional KdV equation. The calculations of the solution are carried out for different values of the time fractional order. These calculations show that the time fractional can be used to modulate the electrostatic potential wave instead of adding a higher order dissipation term to the KdV equation. The results of the present investigation may be applicable to some plasma environments, such as the ionosphere plasma.

Propagations of Rayleigh and Love waves in ZnO films/glass substrates analyzed by three-dimensional finite element method

Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass substrates, the simulation results confirm that the Rayleigh waves along the [0001] direction and Love waves along the [1ˉ100] direction are successfully excited in the multilayered structures. Next, the crystal orientations of the ZnO films are rotated, and the influences of ZnO films with different crystal orientations on SAW characterizations, including the phase velocity, electromechanical coupling coefficient, and temperature coefficient of frequency, are investigated. The results show that at appropriate h/λ, Rayleigh wave has a maximum k^2 of 2.4% in（90°, 56.5°, 0°） ZnO film/glass substrate structure; Love wave has a maximum k^2 of 3.81% in（56°, 90°, 0°） ZnO film/glass substrate structure. Meantime, for Rayleigh wave and Love wave devices, zero temperature coefficient of frequency（TCF） can be achieved at appropriate ratio of film thickness to SAW wavelength. These results show that SAW devices with higher k^2 or lower TCF can be fabricated by flexibly selecting the crystal orientations of ZnO films on glass substrates.

Structural and Acoustic Response of A Finite Stiffened Submarine Hull

After borrowing the idea of precise integration method, a precise integration transfer matrix method （PITMM） is proposed by modifying traditional transfer matrix method. The submarine hull can be modeled as joined conical- cylindrical-spherical shells. By considering the effect of the ring-stiffeners, the field transfer matrixes of shells of revolution are obtained accurately by PITMM. After assembling the field transfer matrixes into an entire matrix, the dynamic model is established to solve the dynamic responses of the joined shell. By describing the sound pressure in fluid by modified wave superposition method （MWSM） and collocating points along the meridian line of the joined shell, finally the structural and acoustic responses of a finite stiffened submarine hull can be predicted by coupled PITMM and MWSM. The effectiveness of the present method has been verified by comparing the structural and acoustic responses of the spherical shell with existing results. Furthermore, the effects of the model truncation, stiffness and thickness on the structural and acoustic responses of the submarine hull are studied.

Acoustic field simulations of logging while drilling by cylindrical finite difference with variable grids

This study proposes an elastic finite difference(FD)time domain method with variable grids in three-dimensional cylindrical coordinates.The calculations will diverge and become less accurate by conventional cylindrical FD as the grid size gradually becomes more extensive with the increasing radius.To prevent grids from being too coarse in far fields,we compensate for the grid cell infl ation by refi ning the grid step in the azimuthal direction.The variable grid FD in the cylindrical coordinate systems has a higher effi ciency in solving acoustic logging while drilling(LWD)problems because the grid boundaries are consistent with those of the drill collar and the borehole.The proposed algorithm saves approximately 94%of the FD grids,80%of the computation time,and memory with a higher calculation accuracy than the FD on rectangular grids for the same models.We also calculate the acoustic LWD responses of the fl uid-fi lled borehole intersecting with fractures.Refl ections are generated at the fractures,which can be equivalent to an additional scattering source.The mode conversions between the collar and the Stoneley waves are revealed.The Stoneley spectra are more sensitive to the fracture.Finally,the logs in a heterogeneous formation with two refl ectors far from the borehole are modeled,and a means of estimating the azimuth of geological interfaces from refl ections is proposed.