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.

Acoustic Simulation with Dynamic Mechanisms in Virtual Reality

Although most invesigators have realized the importance of acoustic simulation in sophisticated VR systems, large computational load involved in this process often contradicts the requirements of real-time interaction, which in return brings on applying the expensive hardware or VR-specific workstations to this area. In order to reduce the computational cost and try to realize the real-time acoustic simulation in software with (or even without) some low-cost hardware, this paper proposes some dynamic mechanisms which can be used as possible strategies embedded into acoustic simulation in VR. Preliminary implementation of those mechanisms has proved to be fairly effective.

SYLLABLE-BASED ACOUSTIC SIMULATION PROCESS IN SECOND LANGUAGE PRODUCTION

This paper discusses how Chinese speakers produce English sonorant consonants embedded in syllable structures novel to them.It shows that speakers may simulate target second language(L2)sounds based on acoustically similar first language(L1)sounds and it is a rather complex process in which sonorant type,vowel context,and articulatory constraints all come into play.Particularly,simulation of L2 speech based on acoustic cues is found to occur not only at the segmental but also at the syllabic level in the English production of Chinese speakers.To explain this finding,the present study proposes an acoustic-articulatory model of L2 syllable production and assumes that the acoustic simulation process is syllable-based and the output form is guided by salient acoustic cues and modulated by both language general and specific coarticulatory mechanisms.

Method of fuzzy uniformity analysis for simulation underwater acoustic signals

Since the simulation underwater acoustic signal is used in the semi-object simulation experiment of underwater weapons, it has great impression upon simulation fidelity. It is asked that whether simulation signals can replace the real signal effectually. Considering the randomness of signals, the interval estimation of feature parameters of simulation signals is made. By comparing the obtained confidence interval with the corresponding accept interval, the concept of similarity coefficient of simulation signals is given. By making a statistical analysis for similarity coefficient, the uniformity information of simulation signals is extracted, and the fuzzy number which expresses the fuzzy uniformity level of simu- lation signals is obtained. The analysis method on fuzzy uniformity of simulation underwater acoustic signals is presented. It is indi- cated by the application in simulation of target radiated-noises that the method is suitable and effectual for the simulation research on underwater acoustic signals, and the analysis result may provide support for decision-making relative to perfecting simulation sys- tems and applying simulation signals.

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.

Simulation of phased array S-scan acoustic field in FSW joint of aluminum alloy extrudate with complex shape

The shape of aluminum alloy extrudate used in high-speed train is complex, structural noises from the surfaces of the extrudate will be received when using ultrasonic phased array to detect the flaws in FSW. To solve this problem, ultrasonic phased array acoustic field model and propagation simulation of acoustic waves were introduced to simulate the acoustic pressure distribution and the propagation of the acoustic waves. With the methods above, the detection parameters can be optimized and as a result, the experimental process can be simplified and the detection efficiency can be improved. Meanwhile, the echoes in the S-scan images can be predicted, which can help analyze the detection results and judge the defects.

Acoustic Prediction and Risk Evaluation of Shallow Gas in Deep-Water Areas

Shallow gas is a potential risk in deep-water drilling that must not be ignored,as it may cause major safety problems,such as well kicks and blowouts.Thus,the pre-drilling prediction of shallow gas is important.For this reason,this paper conducted deep-water shallow gas acoustic simulation experiments based on the characteristics of deep-water shallow soil properties and the theory of sound wave speed propagation.The results indicate that the propagation speed of sound waves in shallow gas increases with an in-crease in pressure and decreases with increasing porosity.Pressure and sound wave speed are basically functions of the power expo-nent.Combined with the theory of sound wave propagation in a saturated medium,this paper establishes a multivariate functional relationship between sound wave speed and formation pressure and porosity.The numerical simulation method is adopted to simulate shal-low gas eruptions under different pressure conditions.Shallow gas pressure coefficients that fall within the ranges of 1.0-1.1,1.1-1.2,and exceeding 1.2 are defined as low-,medium-,and high-risk,respectively,based on actual operations.This risk assessment me-thod has been successfully applied to more than 20 deep-water wells in the South China Sea,with a prediction accuracy of over 90%.

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.

Development of acoustic computer simulation for performance spaces: A systematic review and meta-analysis

This article aims to review the development of acoustic computer simulation for performance spaces. The databases of Web of Science and Scopus were searched for peer-reviewed journal articles published in English between 1960 and 2021, using the keywords for “simulation”, “acoustic”, “performance space”, “measure”, and their synonyms. The inclusion criteria were as follows: (1) the searched article should be focused on the field of room acoustics (reviews were excluded);(2) a computer simulation algorithm should be used;(3) it should be clearly stated that the simulated object is a performance space;and (4) acoustic measurements should be used for comparison with the simulation. Finally, twenty studies were included. A standardised data extraction form was used to collect the modelling information, software/algorithm, indicators for comparison, and other information. The results revealed that the most used acoustic indicators were early decay time (EDT), reverberation time (T30), strength (G), and definition (D50). The accuracy of these indicators differed greatly. For non-iterative simulation, the simulation accuracies of most indicators were outside their respective just noticeable differences. Although a larger sample size was required for further validation, simulations of T30, EDT, and D50 all showed an increase in accuracy with increasing time from 1979 to 2020, except for G. In terms of frequency, the simulation was generally less accurate at lower frequencies, which occurred at T30, G, D50 and T20. However, EDT accuracy did not exhibit significant frequency sensitivity. The prediction accuracy of inter-aural cross-correlation coefficients (IACC) was even higher at low frequencies than it was at high frequencies. The average value of most indicators showed a clear systematic deviation from zero, providing hints for future algorithm improvements. Limitations and the risks of bias in this review were discussed. Finally, various types of benchmark tests were suggested for various comparison goals.

Laboratory simulation on acoustic well-logging with phased array transmitter

Two small scale acoustic phased arrays with 4 elements have been designed and assembled in the laboratory. Experiments have been carried out with them. It is found that both directivity and radiation lobe width of the phased array can be regulated by changing the time delay between the input signals on neighboring elements. Results measured are in good agreement with those calculated. By using the phased array as an acoustic transmitter and hydrophone as a receiver, small scale acoustic well-logging simulations have been carried out both on an aluminum modei well and on a concrete one. Experimental results show that, by increasing the time delay of the input signals on neighboring elements, the steered radiation angle of the phased array becomes larger and larger, and generation conditions of the refracted compressional wave and the refracted shear wave are reached successively, and the refracted compressional wave, the refracted shear wave and the Stoneley wave are strengthened, respec-tively. Therefore, by choosing element spacing of a phased array and acoustic wave frequency appropriately, the main radiation lobe of the phased array can be widened to cover the first critical angle of all kinds of formations, which makes it possible to apply phased array acoustic well-logging in any formation continuously without regulating directivity of the phased array.

Numerical study on static and dynamic fracture evolution around rock cavities

In this paper,a numerical code,RFPA2D（rock failure process analysis）,was used to simulate the initiation and propagation of fractures around a pre-existing single cavity and multiple cavities in brittle rocks.Both static and dynamic loads were applied to the rock specimens to investigate the mechanism of fracture evolution around the cavities for different lateral pressure coefficients.In addition,characteristics of acoustic emission（AE） associated with fracture evolution were simulated.Finally,the evolution and interaction of fractures between multiple cavities were investigated with consideration of stress redistribution and transference in compressive and tensile stress fields.The numerically simulated results reproduced primary tensile,remote,and shear crack fractures,which are in agreement with the experimental results.Moreover,numerical results suggested that both compressive and tensile waves could influence the propagation of tensile cracks;in particular,the reflected tensile wave accelerated the propagation of tensile cracks.

Simulation of the acoustic field emitted from medical linear transducer in a heterogeneous tissue

A numerical model is developed to simulate the acoustic field in heterogeneous tissue from a medical linear transducer.The coupled full-wave equation for nonlinear ultrasound is solved using a staggered-grid finite difference time domain method.The distribution of acoustic pressure and power in human abdominal wall with heterogeneities in sound speed,density,and nonlinear parameter are obtained.Compared with homogeneous medium,when sound speed in tissue is uniform and density unchanged,the acoustic energy decreases only1.8 dB in the focal region;when density in tissue is uniform and sound speed unchanged,the energy decreases 3.8 dB in the focal region,which is almost the same as heterogeneous tissue.Thus,the primary factor of the aberration of focused beam is the heterogeneous distribution of the tissue sound speed.

Correlation equalizer and its simulation for underwater acoustic channels

A robust correlation equalizer is presented to resolve the equalization divergence,which is caused by the distribution of channel response function's zero points outside the unit circle. The equalizer is robust except that the channel response function has zero points just on the wht circle. Noncausal two-side prediction is also proposed. Based on the prediction and FFT the transversal form of the correlation equdrier and the fast algorithIn are developed. The simulation results of the correlation equallzer of transversal form are given.

The noise spectral characteristics and noise reduction schemes of screw air-source heat pump:A case study

The screw air-source heat pump can cause incessant high noise levels during operation,which might hinder adoption of this energy-efficient heat pump.First,acoustic measurements and comparison testing were performed in this research.The measurements revealed that the compressor is the main noise source of the heat pump,and it shows a multipeak frequency distribution and a wide frequency spectrum under different work conditions,with multiple peaks at 63,250,and 1000 Hz.Then,a compressor sound insulation cover with broadband absorption was proposed,and it was experimentally proven that the insulation cover can reduce the maximum sound pressure level of one unit from 89.8 dBA to 79.1 dBA.Third,we proposed several noise reduction strategies and compared their noise reduction effects using computer simulation.The results showed that the noise problem can be effectively improved through the rational design of the sound barrier and the layout and opening options of heat pump.The distance between the sound barrier and heat pump and the sound attenuation due to diffraction ALa exhibit a U-shaped relation.For buildings of different heights,the optimal heights of noise barrier are proposed.The 5.5-meter is the optimal height of the sound barrier for single-story buildings.The conclusions can be applied to other building projects for heat pump noise reduction.

Distinguishing Oil and Water Layers in a Porous Cracked Medium by Interpreting Acoustic Logging Data on the Basis of Hudson Theory

During surveys, water layers may interfere with the detection of oil layers. In order to distinguish between oil and water layers in a porous cracked medium, research on the properties of cracks and oil and water layers and their relation to acoustic logging rules is essential. On the basis of Hudson＇s crack theory, we simulated oil and water layers in crack-porous medium with different crack parameters corresponding to the well-field response. We found that in a cracked medium with high crack angle or low number density of cracks, compressional and shear wave velocities are sensitive to crack characteristics; further, these velocities are more sensitive to crack characteristics when the waves propagate through the water layer than when they propagate through the oil layer. Compressional and shear wave velocities increase with an increase in crack angle： in the water layer, the increase is approximately linear. On comparing the full waveforms observed in the oil and water layers, we find that the amplitudes of most waves are higher in the water layer. Among the considered waves, the Stoneley wave suffers maximum amplitude attenuation in the oil layer. The maximum excitation intensity for oil layer is greater than that for the water layer. These results can guide further cracked media logging field exploration work.

Simplified Method of Simulating Double-Layer Micro-Perforated Panel Structure

The micro-perforated panel(MPP)structure has been widely used in various noise control applications,and thus its acoustic performance prediction has been receiving increasing attention.The acoustic performance of simple MPP structures,such as a MPPsound absorber,has been predicted using an analytical calculation method.However,this is not a suitable approach toward predicting the acoustic performance of complicated MPP structures,owing to the structural complexity of these structures.Moreover,the many perforations of submillimeter scale diameter render the MPP structures very difficult to analyze using numerical simulation.Thus,this study focused on two different simplified MPP simulation methods:the transfer admittance method and the equivalent fluid method,and their application on double-layer MPP structures.Based on the two simplified MPP simulation methods,the transmission loss value of the double-layer MPP mufflers with two sets of different structural parameters was calculated,respectively.The predicted results were compared with the impedance tube measurements.The results revealed that the two simplified MPP simulation methods could effectively predict the acoustic performance of doublelayer MPP structures.Moreover,the prediction based on the transfer admittance method can outperform the two simplified simulation methods.