In the method of far-field acoustic holographic noise identification, the Kirchhoff formula has been widely used. Achievement of reconstruction is the outstanding advantage of acoustical holography over the widely-used beam-forming method. In this paper, the mechanism of acoustical holography based on finite Kirchhoff diffraction is established, the relationship between holographic aperture angle and the quantitative accuracy is revealed. With simulations for known sound source, the quantitative mechanism is validated successfully.
Quantitative accuracy and boundary conditions of Kirchhoff diffraction in optics are studied and verified ; in which infinite Kirchhoff diffraction model was established and Fresnel wave integral was employed to calculate the quantitative results. The integral results were correlated with the inclination factor and wave front which was determined by integral domain as integration of the surface instead of the aperture plane. Inspired by this work, in our recent study, quantitative accuracy of Kirchhoff diffraction Acoustical Holography was migrated and studied in sound source identification and reconstruction. The relationship between quantitative accuracy of it and the aperture angle (inclination factor) when applied in sound source reconstruction is established.
2. Quantitative mechanism of Kirchhoff diffraction acoustical holography
Kirchhoff diffraction theory is traditional and widely used in optics and acoustics. The classic infinite Kirchhoff diffraction formula is as Eq. (1):
Based on Eq. (1), Kirchhoff diffraction acoustical holography method was established and originally introduced in references [2-4]. According to Kirchhoff diffraction Acoustical Holography principle, the sound pressure distribution in the holographic plane can be reconstructed with Eq. (2):
where is the frequency domain sound pressure of the point in holographic plane, is the wave number, is the holography constant, is the distance between point in measurement plane and . is the conjugate sound pressure of point . From Eq. (1) and Eq. (2), it can be seen that without considering Doppler-effect and discretization error, the main quantitative error theoretically comes from the finite integral domain. In another word, the integral domain error stands for the quantitative accuracy of Kirchhoff diffraction Acoustical Holography.
The Kirchhoff diffraction model of point source is shown in Fig. 1. For the ideal point source model in the free space, the parameters in Eq. (1) are as shown in Fig. 1(a). The strict Kirchhoff diffraction assumes that point source is located in which radiates the sound wave conjugated with the wave radiated from the point source in position . If the sound field distribution in the closed spherical surfaceis obtained, Kirchhoff diffraction can be applied to reconstruct the accurate sound source with Eq. (1).
Fig. 1Ideal point source model in free space without diffraction: a) complete spherical closed surface, b) finite plane
In the practical measurement of two-dimensional microphones-array, the surface S is a finite and discrete plane. For finite Kirchhoff diffraction conditions, Huygens-Fresnel principle and Fresnel diffraction schematic diagram are shown in Fig. 1 (b). Spherical waves are radiated by any secondary point source in the plane . Located in space , the coordinates of the point source is . is the sound pressure. is the displacement vector from to . For point sources, Eq. (3) and Eq. (4) From Eq. (1) as follows:
where is the sound source amplitude and is the angle between sound source and holographic plane. Then the closed spherical surface and finite planar Kirchhoff diffraction models are separately shown in Fig. 2.
Fig. 2Fresnel’s zone construction algorithm schematic for spherical and planar model: a) complete spherical enclosed surface, b) finite plane
With Eq. (3) and Eq. (4), the original Kirchhoff diffraction formula can be converted into Eq. (5):
where is a dimensionless inclination factor which changes with the angle . In Eq. (5), 1, , . In order to calculate the integral in Eq. (5), Fresnel’s zone method is applied as shown in Fig. 2. Take point as the centre and , , , , ..., as radius to draw a series of spherical surface , , ..., , …, are deemed same between neighbouring under the condition that the zone is narrow enough. In plane , the differential surface element is calculated as Eq. (6):
where is azimuthal angle round axis. This formula converts into three steps as Eq. (7) shows:
So, the contribution of the th wave zone to is:
Calculate the integration, can be expressed as Eq. (9):
The factor can be neglected in the far field condition, which means that the measurement distance is more than 20. Under this approximation, the integral result can be converted to Eq. (10):
Eq. (10) shows the relative value between the reconstructed holographic result and the theoretical reference value (). Fig. 3 gives the curve of changing with the aperture angle from 0 to .
Fig. 3 shows the calculated quantitative accuracy against the aperture angle.
Fig. 3Curve of ratio between Kirchhoff diffraction acoustical holography value and the reference value changing with aperture angle
The far-field sound source identification simulation of finite Kirchhoff diffraction is carried out by Huygens-Fresnel theory, and the limit of accurate prediction is studied. A near-uniform grid microphone array measures a 1000 Hz point source at a distance of 2 meters. Table 1 and Fig. 4 show the relationship between the reconstruction result and the inclination coefficient considering the aperture angle.
Table 1Results comparing with the inclination factor.
Scale factor ()
As can be clearly seen from Fig. 4, when the aperture angle increases, the reconstruction results gradually reach the actual value. For tilt factor curves and curves, the fitting of the two curves proves that the critical holographic aperture established in our work is effective for determining the boundary conditions of Kirchhoff diffraction acoustic holography, and the fitting curves are in good agreement with the theory value.
The results can provide the basis for improving the quantitative accuracy of acoustic holography. The reconstruction with a fixed coefficient which can be obtained based on the known sound source.
From the theoretical derivation, the limiting conditions are clear that the sound source must be high-frequency to ensure the Fresnel’s zone is narrow enough for approximation. The microphone array should be extremely close-set to ignore the discretization error, although this will affect the result of flow filed computation greatly.
Fig. 4Comparison of theoretical and experimental values
In this paper, the quantitative mechanism of far-field Kirchhoff acoustical holography has been elaborated systematically. Holographic aperture angle is a key factor influencing the quantitative accuracy of acoustical holography.
The results were shown that the spherical wave incident result produced by the wave-front preparation is equivalent to the result of the aperture plane formula at the incident wave plane. The simulation results were shown that the fitting curve of Kirchhoff acoustic holography agrees well with the theoretical value. However, according to the assumptions adopted by our approach, there are still some academic research will be done in the future.
This work is based on the research of Project 51805229 and Project 51705221, supported by the National Natural Science Foundation of China.
Wang Z., Wu S. F. Helmholtz equation least-squares method for reconstructing the acoustic pressure field. The Journal of the Acoustical Society of America, Vol. 102, Issue 4, 1997, p. 2020-2032.
Wu S. F., Zhao X. Combined Helmholtz equation-least method for reconstructing the acoustic radiation from arbitrary shaped objects. The Journal of the Acoustical Society of America, Vol. 112, Issue 1, 2002, p. 179-188.
Earl G. Williams, HoustonBrian H. Fast Fourier transform and singular value decomposition formulations for patch nearfield acoustical holography. The Journal of the Acoustical Society of America, Vol. 114, Issue 3, 2003, p. 1322-1333.
Song L., Koopmann G. H., Fahnline J. Numerical errors associated with the method of super-position for computing acoustic fields. Journal of the Acoustical Society of America, Vol. 89, Issue 6, 1991, p. 2625-2633.