Simulation analysis of underwater structure based on local cylindrical acoustic holography

Gaoquan Gu1 , Jing Zhang2 , Yingchao Shi3

1, 3College of Electronic Engineering, Naval University of Engineering, Wuhan, China

2China Ship Development and Design Center, Wuhan, China

1Corresponding author

Vibroengineering PROCEDIA, Vol. 23, 2019, p. 133-137.
Received 17 March 2019; accepted 29 March 2019; published 25 April 2019

Copyright © 2019 Gaoquan Gu, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License

The local cylindrical acoustic holography technology is introduced, and its application in identifying and locating the column is introduced. The implementation algorithm of local cylindrical holographic reconstruction is given. In engineering operations, columnar sound sources such as electric motors, water pipes, and underwater structures are common noise sources in engineering. Therefore, this paper is of great significance for the location identification and control of such sound sources. The simulation results show that the proposed local cylindrical acoustic holography algorithm is useful to identify the noise source and reconstruct the sound field correctly.

Keywords: local cylindrical acoustic holography, radiation sound field, underwater structure, space splice.

1. Instruction

With the development of underwater acoustics and signal processing technology, near-field acoustic holography technology has made great progress. By processing and inversing the information of amplitude, phase and pressure of sound wave, the directivity information of sound field and noise source can be obtained [1]. In theory, near-field acoustic holography algorithms include equivalent source method [2], boundary element method [3], and acoustic field space transformation method [4], among them. The space transformation method of sound field has high computational efficiency and resolution, and its practicability is strong. However, this method requires that the measurement surface and the sound source surface are of the same shape. Domestic scholars have done a lot of theoretical and Experimental Research on this method [5-6], but they all focus on the plane. If the sound source is cylindrical, it will cause greater errors.

For the non-planar sound source, the planar near-field acoustic holography can only obtain the sound field distribution of the equivalent planar sound source on the “virtual source surface” on the reconstructed surface, and can't know the true sound field distribution of the non-planar sound source surface, thus affecting accurate identification and analysis of non-planar sound sources, columnar sound sources or columnar sound sources are sound source structures often encountered in engineering, such as aero engines, electric motors, water pumps and various pipes. For columnar sound source structures, cylindrical near-field acoustic holography must be used to obtain the true sound field distribution of the sound source surface.

According to the above problem, this paper proposes a method of partial cylindrical acoustic holography.

2. The realization process of cylindrical acoustic holography

2.1. Deduction process

The measurement of the cylindrical near-field acoustic hologram is shown in Fig. 1: the cylinder with the radius of rH is the holographic measuring surface, the length of the measuring surface along the axial direction is 2L, and the measuring surface is expanded as shown in the right side of the figure; The cylinder with radius rS is the reconstruction surface.

Fig. 1. Schematic diagram of measurement and reconstruction of cylindrical acoustic holography

Schematic diagram of measurement and reconstruction of cylindrical acoustic holography

As shown in Fig. 1, the axial measurement interval on the holographic surface is z=2L/(2m+1); The circumferential angle of the measurement is ϕ=2π/(2N+1), where M and N are positive integers, the (2M+1)×(2N+1) point discrete sound pressure is obtained on the holographic surface. And can be expressed as -NnN,-MmM. Then, pn(rH,kz) discrete space Fourier transform that can approximate the holographic surface sound pressure of p(rH,nϕ,mz):

p n ( r H , i k z ) = z ϕ n = - N N m = - M M p ( r H , n ϕ , m z ) w 2 M + 1 i m w 2 N + 1 j n ,

where m and j are integers from one M to M; n and j are integers from one N to N; kz=π/L; w2M+1 and w2N+1 are respectively e-2πi/(2M+1) and e-2πi/(2N+1). From the discretized pn(rH,ikz) and Eq. (4):

p n r S , k z H n 1 k r r s H n 1 k r r H p n r H , k z ,

the sound pressure cylinder spectrum on the source surface can be approximated as:

p n r s , m k z = H n 1 k   2 - ( m k z ) 2 r s H n 1 k   2 - ( m k z ) 2 r H p n r H , m k z .

Therefore, the sound pressure at the discrete points on the reconstructed surface can be obtained by the discrete space Fourier inverse transform:

p n ( r s , n ϕ ) = 1 4 π L j = - N N i = - M M p n ( r s , i k z ) w 2 M + 1 - i m w 2 N + 1 - j m .

From the discretized implementation Eqs. (3), (5) and (6), the cylindrical near-field acoustic holography can be transformed by the computer to obtain the sound pressure field on the reconstructed or predicted surface.

2.2. Cylindrical spectral filtering

Wavenumber domain windowing filtering can be used to suppress the amplification process of high spatial frequency component errors. The window function used is similar to the planar near-field acoustic holography, just replace kx and ky in kr with the axial wave number kx and the circumferential wave number n/r:

I I ( n , k z ) = 1 - 1 2 exp k r k c - 1 a ,           k r k c , 1 2 exp 1 - k r k c a ,         k r > k c .

In the formula, kr=kz2+(n/r)2; a is the steepness coefficient of the window function; kr is the cutoff wave number.

This is a circularly symmetric exponential filter. The choice of kc is very important and directly affects the accuracy of the reconstruction results. A method for selecting the holographic surface and the radius of the reconstructed surface and the signal to noise ratio of the holographic measuring system is given below.

Let the sound pressure signal-to-noise ratio of each point of the holographic surface be SNR, pnt(rs,kz) is the cylindrical spectrum of the theoretical sound pressure pt(rH,kz) of the reconstructed surface with radius rs, known by Eq. (8):

k n ( k r 1 r ) k n ( k r 1 a ) a r e - k r ' r - a .

The theoretical sound pressure of the holographic surface with radius rH is the relationship between pnt(rH,kz) and pnt(rs,kz):

p n t r H , k z p n t r S , k z r s r H e - k r r H - r s .

Let pnt(rH,kz) be the spiral spectrum of the noise signal component pt(rH,ϕ,z) contained in the holographic surface sound pressure, then the cylindrical spectrum p(rH,kz) of the actual sound pressure on the holographic surface can be expressed as the sum of the theoretical sound pressure cylinder spectrum and the noise cylinder spectrum:

p r H , k z = p n t r S , k z r s r H e - k r r H - r s + p n e r H , k z .

According to Eqs. (8), the reconstructed surface acoustic pressure spiral spectrum is:

p r S , k z = p n r H , k z r H r S e - k r r S - r H .

Available from Eqs. (10) and (11):

p r S , k z = p n t r H , k z r H r S p n t e - k r r S - r H .

At this point, the reconstructed surface noise spiral spectrum pnt(rS,kz) can be expressed as:

p n e r S , k z = p r S , k z - p n t r H , k z + r H r S p n t r H , k z e - k r r S - r H .

To ensure that the reconstruction results are meaningful, it must be ensured that the theoretical sound pressure cylinder spectrum on the reconstructed surface is not overwhelmed by noise, that is, it should be satisfied:

p n e ( r S , k z ) p n t r S , k z .

According to Eq. (13), it needs to be satisfied:

r H r S p n e ( r H , k z ) e - k r ( r s - r H ) p n t r S , k z .

According to the previous analysis, in the process of sound wave propagation from reconstruction to holographic surface, the amplitude of the propagation wave is attenuated by the square root of the distance, and the amplitude of the evanescent wave component is attenuated exponentially with the propagation distance. Therefore, the following formula is established:

k z 2 + k o 2 k 2 + 1 ( r H - r S ) 2 S N R 1 n 10 20 + 1 n r S r H 2 = k q .

Eq. (16) shows that as long as the wavenumber component of kz2+ko2>kq in the cyber surface acoustic pressure cylinder spectrum is filtered out, it can be ensured that any cylindrical wave component obtained by reconstruction is not overwhelmed by noise. Thus, according to the Eq. (16), it can be determined that the highest wave value used in the reconstruction process is kq, and this value is used as the cutoff wave number of the wavenumber domain filtering in the reconstruction process.

3. Simulation analysis

3.1. Simulation conditions and steps

A group of wideband signals is simulated, and the sound field of each frequency point is reconstructed by local aperture acoustic holography algorithm. The sound field reconstruction accuracy of local aperture acoustic holography algorithm is analyzed, and the influence on the calculation of relative local radiated sound power level is evaluated. Sound source: point source, coordinates (0.75, 075); simulation size: S1= 28×16 = 13.5 m×7.5 m; array spacing: detal = 0.5 m; the distance of the holographic surface from the sound source: Zh= 0.5 m; reconstruction distance from the sound source: Zs= 0.4 m; simulation frequency: f= 20 Hz-1 kHz.

3.2. Calculation steps

1. Simulate a set of Gaussian noise, through low-pass filtering, the upper limit of the low-pass filter cutoff frequency is 1 kHz.

2. Setting parameters and bringing into the broadband near-field acoustic holography algorithm module.

3.Calculate the reconstructed sound field and relative local radiated sound power level for each frequency point, and the total relative local radiated sound power level of the reconstructed sound field and the theoretical sound field within the bandwidth.

4. Conclusion

It can be concluded from the simulation that the local cylindrical acoustic holography can accurately reconstruct the noise source position and the relative local radiated sound power level.

Fig. 2. Spectrogram of the first channel

 Spectrogram of the first channel

Fig. 3. Theoretical and reconstructed values of the sound field at 0.4 m

 Theoretical and reconstructed values of the sound field at 0.4 m


  1. Maynard J. D., Williams E. G., Lee Y. Near field acoustical holography I: theory of generalized holography and the development of NAH. Journal of the Acoustical Society of America, Vol. 78, Issue 4, 1985, p. 1395-1413. [Publisher]
  2. Zhang H. B., Jiang W. K., Wan Q. Theoretical study of cyclically stationary near-field acoustic holography by boundary element method. Journal of Acoustics of China, Vol. 33, Issue 3, 2008, p. 231-237. [CrossRef]
  3. Sun C., He Y. A., Liu Y. C. Wave superposition method for reconstructing radiated acoustic field of cylindrical shells. Noise and Vibration Control, Vol. 10, Issue 5, 2012, p. 31-36. [CrossRef]
  4. Lee L. Z., Lee J., Lu B. W. Determination of regularization parameters in planar near-field acoustic holography. Journal of Acoustics of China, Vol. 28, Issue 6, 2003, p. 555-560. [CrossRef]
  5. Xu L., Bi C. X., Chen X. Z. Application of Papoulis-Gerchberg algorithm in Patch near-field acoustic holography. Inter-Noise, Shanghai, China, 2008, p. 2330-2337. [CrossRef]
  6. Lee Q., Bi C. X., Xu L. A near-field acoustic holographic noise source recognition system based on Labview. Noise and Vibration Control, Vol. 12, Issue 6, 2010, p. 140-144. [CrossRef]