An Extension of the Virtual Rotating Array Method Using Arbitrary Microphone Configurations for the Localization of Rotating Sound Sources
An Extension of the Virtual Rotating Array Method Using Arbitrary Microphone Configurations for...
Jekosch, Simon;Sarradj, Ennes
acoustics Article An Extension of the Virtual Rotating Array Method Using Arbitrary Microphone Conﬁgurations for the Localization of Rotating Sound Sources Simon Jekosch * and Ennes Sarradj Institute of Fluid Mechanics and Engineering Acoustics, Technical University Berlin, Einsteinufer 25, 10587 Berlin, Germany; email@example.com * Correspondence: firstname.lastname@example.org; Tel.: +49-30-314-25211 Received: 17 April 2020; Accepted: 12 May 2020; Published: 15 May 2020 Abstract: The characterization of rotating aeroacoustic sources using microphone array methods has been proven to be a useful tool. One technique to identify rotating sources is the virtual rotating array method. The method interpolates the pressure time data signals between the microphones in a stationary array to compensate the motion of the rotating sources. One major drawback of the method is the requirement of ring array geometries that are centred around the rotating axis. This contribution extends the virtual rotating array method to arbitrary microphone conﬁgurations. Two different ways to interpolate the time signals between the microphone locations are proposed. The ﬁrst method constructs a mesh between the microphone positions using Delaunay-triangulation and interpolates over the mesh faces using piecewise linear functions. The second one is a meshless technique which is based on radial basis function interpolation. The methods are tested on synthetic array data from a benchmark test case as well as on experimental data obtained with a spiral array and a ﬁve-bladed fan. Keywords: rotating sources; virtual rotating microphone array; beamforming; fan noise 1. Introduction Beamforming and other array methods for acoustic source localisation have been extensively investigated for multichannel microphone measurements. The characterisation of rotating aeroacoustic sources using microphone array methods has been proven to be a useful tool. Various methods for identifying rotating noise source have been proposed. They can be classiﬁed as methods in time domain or in frequency domain. Rotating beamforming in time domain with a rotating focus point has been applied using the Rotating Source Identiﬁer (ROSI)  for wind turbine and helicopter blades. The method was later combined with deconvolution algorithms as a hybrid time and frequency domain method  and applied to rotating sources . An improved way to calculate the point spread function for the deconvolution of rotating sources was derived by Debrouwere et al. . Rotating beamforming in the frequency domain was ﬁst applied for induct measurements. The method solves the Green’s function for a rotating acoustic source in a duct [5,6]. The free ﬁeld solution of the Green’s function for a rotating monopole in spherical coordinates was deduced by Pannert and Maier . This method is known as modal decomposition (MD). Another method in the frequency domain to determine the location of acoustic sources is the virtual rotating array (VRA) method introduced by Herold and Sarradj . It compensates the motion of the source by virtually rotating the array microphone positions. Since both methods work directly in the frequency domain the construction of the cross-spectral matrix and the application of subsequent beamforming, deconvolution and inverse techniques are easily possible. Acoustics 2020, 2, 330–342; doi:10.3390/acoustics2020019 www.mdpi.com/journal/acoustics Acoustics 2020, 2 331 Various comparisons of the three methods have been conducted [9–11]. The resulting delay and sum beamforming output for source localisation have been found to be almost equal for all the methods. Though, the computational cost of ROSI was found to be two to three orders of magnitude higher than that of the VRA method. The major drawback of the VRA method is the restrictive microphone arrangement. It relies on evenly spaced microphone placements in one or multiple concentric rings around the center of the rotation. Geometries with uneven microphone placements been studied, but still rely on positioning on a ring . To overcome this drawback, a new interpolation method for the virtual rotating array for arbitrary two dimensional and three dimensional microphone placements is proposed here. The remainder of the paper is organised as follows. In Section 2, the beamforming methods and the implementation of the virtual rotation are described. Application of the methods on simulated benchmark data is presented in Section 3. The test of the methods with experimental data is discussed in Section 4. Section 5 concludes the paper. 2. Methods 2.1. Beamforming and Deconvolution Methods The virtual rotating array method is based on beamforming in the frequency domain. The beamforming output (the squared sound pressures at one of the N focus points x ) in the frequency domain is calculated via b(x ) = h (x ) C h(x ), s = 1, ..., N (1) s s s with the cross-spectral matrix C of the microphone signals, the steering vector h and the hermitian of the steering vector h . Here, the latter is chosen according to formulation III from Sarradj : jk(r r ) s,m s,0 h = e , m = 1...M. (2) r r r s,0 s,m å s,l l=1 with the number of microphones M, the wavenumber k, r and r indicating the distances between s,m s,l the sources implicitly assumed stationary and the microphone locations and r between the the s,0 sources and the array centre location. This formulation has been shown theoretically to yield the correct source level . In this contribution the (DAMAS)  deconvolution is used. The DAMAS algorithm assumes that the beamforming output b consists of the actual source distribution being convoluted by the point spread function from each focus point to all focus points P. The correct source distribution b can be revealed by solving the system of equations: b = Pb (3) The companding convex optimization problem can be written as: 0 0 arg min kb Pbk , b 0 (4) with kk denoting the L2-Norm and arg min being the arguments of the minima function which 0 0 returns the set of points b for which kb Pbk attains the smallest value. The problem can be solved using a non negative least squares (NNLS) algorithm . Since beamforming and deconvolution in this form only works for stationary sound sources, the motion of the rotating source has to be compensated. Acoustics 2020, 2 332 2.2. Virtual Rotating Array The motion compensation with a virtual rotating array (VRA)  tracks the angular position of the rotating source and rotates the virtual microphone positions at the same rate as the source. The signal at each virtual microphone positions is interpolated linearly from the two neighbouring microphones. The VRA method with a ring array uses only the two neighbouring microphones to interpolate the signal at the virtual microphone positions. Since every point on the ring has exactly 2 neighbouring microphones this approach is suitable. When using arbitrary microphone geometries the neighbouring microphones cannot be deduced from the virtual microphone position, which makes this interpolation method insufﬁcient. For the extension of this algorithm to 2D and 3D array geometries an interpolation process that is based on the microphone positions has to be applied. Two methods are proposed here for the interpolation over irregularly distributed microphone positions: mesh based triangulation and radial basis function approximation. 2.3. Mesh Based Triangulation For the ﬁrst method, a two- or three-dimensional mesh needs to be constructed between all neighbouring microphones. A common way to derive such a mesh is the Delaunay triangulation . The triangulation for the microphone positions x is constructed in such a way that all circumcircles of all triangles have empty interiors, i.e., no other microphone position lies inside the circle. For this contribution all Delaunay triangulations have been calculated with the qhull package which is based on the quickhull algorithm by Barber et al. . Figure 1 shows the Delaunay triangulation of a 2-D spiral microphone geometry. The array consists of 64 microphones which are ordered in a single sunﬂower spiral . The spiral parameter H which deﬁnes the radial distribution of the microphones is chosen to be 1.0 and the parameter V deﬁning the azimuthal distribution of points is set to 5.0. The parameters are chosen according to Sarradj . The spiral array has an aperture of d = 1 m. ring Since the virtual rotation is carried out in the angular domain the triangulation is also transferred in polar coordinates. The transferred Delauney triangulation is shown in Figure 1b. The angle of interest is set in the domain f 2 [ p, p]. To achieve angular periodicity of the interpolation, the mesh is extended to both sides with the repeating microphone position points in angular direction. 0.6 0.6 0.4 0.4 0.2 0.2 -0.2 -0.4 -0.6 -π -π/ 2 0 π/ 2 π -0.6 -0.4 -0.2 0 0.2 0.4 0.6 x/m φ/rad (a) (b) Figure 1. Delaunay triangulation of the geometry in Cartesian coordinates (a) and in cylinder coordinates (b) with repeating microphone position points in angular direction. The black dots denote the microphone positions, the grey lines indicate the edges of each mesh triangle. The red and blue points indicate transferred microphone positions. y/m r/m Acoustics 2020, 2 333 The interpolation between the points in the Delaunay triangulation is achieved by barycentric interpolation between neighbouring points in the cylindrical-coordinate space. For a 2D Microphone geometry each point in the plane has three vertices of the triangle it is part of for the interpolation, for 3D geometries each point inside the triangulation is interpolated by the four tetrahedron vertices. This algorithm is limited to interpolation inside the convex hull of the original array. Therefore the aperture of the virtual array has to be the same or smaller than the microphone array. This strict limitation can be overcome with the meshless algorithm which interpolates for the whole domain. 2.4. Radial Basis Functions The meshless method for the virtual rotation uses radial basis function interpolation. A radial basis function (RBF) is an n-dimensional radially symmetric function F : R ! R which depends only on the Euclidean distance between the center of the RBF x and the point of evaluation x. The Euclidean distance is deﬁned as r = (kx x k ). There are various RBF’s available in the literature that can be used for interpolation . The most commonly used RBF’s are listed in Table 1, with e being a scaling factor. Table 1. Radial basis functions. Multiquadric (r e) + 1 Inverse 1.0/ (r e) + 1 (re) Gaussian e Cubic r Thin plate r log(r) The pressure time signal is reconstructed at any point x by a linear combination of the RBF’s centred at the original microphone positions x multiplied by the weight coefﬁcients w : i i p(x) = w F(kx x k ) (5) å i i 2 i=1 and the weight coefﬁcients are determined by solving the linear system 2 32 3 2 3 F(kx x k ) F(kx x k ) F(kx x k ) w p 1 1 2 2 1 2 M 1 2 1 1 6 76 7 6 7 F(kx x k ) F(kx x k ) F(kx x k ) w p 6 76 7 6 7 1 2 2 2 2 2 M 2 2 2 2 6 76 7 6 7 = (6) . . . . . 6 76 7 6 7 . . . . . . 4 . . . 54 . 5 4 . 5 F(kx x k ) F(kx x k ) F(kx x k ) w p 1 M 2 2 M 2 M M 2 M M with p p being the discrete pressure values at the microphone positions. This linear system 0 M can be solved iteratively and results in new time signals p(x) at any virtual microphone position. This offers the possibility to extrapolate outside of the array aperture. Although, this only works for points close to the microphone positions. For points further away from the array, the extrapolation is not valid and strongly dependent on the chosen basis function. 3. Simulation 3.1. Setup The data used for evaluation are a part of a benchmark series for microphone array measurements. The benchmarks were put together to provide an evaluation standard for array method algorithms and their individual implementations . The benchmark that is used in this contribution is the B11 benchmark . The data of this benchmark comprise two sub-cases, each featuring simulated rotating monopole sources. The calculations of this contribution focus on the sub-case B. This test case consists Acoustics 2020, 2 334 of three rotating monopoles with a rotational rate that slightly varies over time. Source 1 and 2 are on a radius of 0.25 m with an angle difference Dj = 40 . The source power of the second source is 3 dB below the power of the ﬁrst source. The third source rotates on a radius of 0.125 m and has a 3 dB lower strength than the second source. The rotational rate is tracked by a tacho signal which provides one trigger pulse per revolution. The rotational angle is calculated from the tacho signal. In order to evaluate the new virtual rotating array method the experimental setup in the Benchmark B11 has been extended. The microphone geometry was changed from an equidistant ring array to a planar spiral array. The number of microphones is kept at 64. The parameters of the spiral are the same as described in Section 2.3. Figure 2 shows the B11 benchmark setup as well as the extended case with the spiral array geometry and Table 2 shows the evaluation parameters for the calculation of the beamforming results. 0.4 0.4 0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 x/m x/m (a) (b) Figure 2. Experimental setup of the benchmark B11b (a) and the extended benchmark with a spiral microphone array (b). Blue dots mark the microphone positions and red dots mark the rotating source positions. The black circles show the trajectory of the rotating sources. Table 2. Simulation data processing parameters. Number of microphones 64 Array aperture 1 m Sampling rate 48 kHz FFT block size 1024 FFT window/overlap von Hann/50% Focus grid distance 0.5 m Focus grid resolution 0.01 m Focus grid area 1.00 m 1.00 m 3.2. Results Two different interpolation methods are shown in this contribution. The modiﬁed test case with the spiral array geometry with Delaunay triangulation and barycentric interpolation as well as radial basis function interpolation with cubic basis functions. They are compared to the beamforming results without motion compensation and the VRA method using a ring array and linear interpolation. The Figures 3–6 show the octave band source maps for three different center frequencies. The beamforming output was calculated using conventional beamforming in the frequency domain. For the beamforming result without motion compensation the sources are ‘smeared’ over their trajectory along the ring. The octave band around 2500 Hz shows only one ring for all sources, while to two higher frequency bands show the separate third source at the 0.125 m radius. At 10 kHz, the effect of the beamforming sidelobes is clearly visible. y/m y/m Acoustics 2020, 2 335 Figure 4 shows the beamforming result for the B11 benchmark setup with a ring array and linear interpolation. The Beamforming results with a spiral array are shown in Figure 5 using motion compensation with Delaunay triangulation and in Figure 6 using motion compensation with radial basis function interpolation. For all three cases the sources appear at the positions given in the benchmark B11, which indicates that the virtual rotation of the microphone position is feasible for reconstructing rotating source positions. The methods using the spiral array are able to interpolate the pressure values at the virtual microphone positions in a similar way. The width of the beamforming lobes depends strongly on the frequency. For the lowest octave at 2500 Hz the source positions are correct but the source separation of the three single sources is not possible because at this frequency the distance of the sources is below the Rayleigh criterion . The source maps for the 5 kHz octave band are almost equal for the Delaunay and RBF interpolation. The three sources are at the correct positions and distinguishable. For very high frequencies around the 10 kHz band the sources are also displayed correctly but the RBF interpolation produces slightly more artefacts than the Delauney triangulation outside the source regions. A quantitative comparison of the side lobe levels and the beam width for source 1 is given in Table 3. The beam with is deﬁned by the diameter of the 3 dB level difference distance from source 1’s maximum level. The dynamic is the level difference between source 1’s maximum level and the level of the highest side lobe. In comparison the interpolation method using Delauney triangulation leads to a slightly better performances than the radial basis function interpolation in beam width as well as dynamic for this test case. The resulting beamforming output using the ring array produces less artefacts at higher frequencies than the spiral array geometry, but produce higher side lobes levels at all octave bands. Those differences are related to the array geometries. The ring array uses the full array aperture and therefore produces a smaller main lobe than the spiral array. The spiral array has a higher dynamic range than the ring array and produces less side lobes. The ring array is less sensitive to artefacts at high frequencies due to the smaller distances between the microphone positions compared to the spiral array. 2500 Hz 5000 Hz 10000 Hz Max 0.4 0.4 0.4 -3 0.2 0.2 0.2 -6 0.0 0.0 0.0 -9 − 0.2 − 0.2 − 0.2 -12 − 0.4 − 0.4 − 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 x/m Figure 3. Octave band source maps for 2.5, 5 and 10 kHz obtained with the spiral array and no virtual rotation of the microphones. Table 3. Comparison of array geometry and interpolation method regarding side lobe dynamic and beam width of the main lobe. All results are calculated for source 1. Triangulation, Spiral RBF, Spiral Linear, Ring Frequency 2500 Hz 5000 Hz 10,000 Hz 2500 Hz 5000 Hz 10,000 Hz 2500 Hz 5000 Hz 10,000 Hz Beam width 0.125 m 0.078 m 0.055 m 0.129 m 0.085 m 0.058 m 0.071 m 0.033 m 0.019 m Dynamic 19.95 dB 12.64 dB 9.71 dB 14.45 dB 10.12 dB 9.67 dB 9.53 dB 8.81 dB 6.18 dB y/m L /dB p Acoustics 2020, 2 336 2500 Hz 5000 Hz 10000 Hz Max 0.4 0.4 0.4 -3 0.2 0.2 0.2 -6 0.0 0.0 0.0 -9 − 0.2 − 0.2 − 0.2 -12 − 0.4 − 0.4 − 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 x/m Figure 4. Octave band source maps for 2.5, 5 and 10 kHz obtained with the ring array and linear interpolation between the microphones. 2500 Hz 5000 Hz 10000 Hz Max 0.4 0.4 0.4 -3 0.2 0.2 0.2 -6 0.0 0.0 0.0 -9 − 0.2 − 0.2 − 0.2 -12 − 0.4 − 0.4 − 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 x/m Figure 5. Octave band source maps for 2.5, 5 and 10 kHz obtained with the spiral array and barycentric interpolation between the microphones. 2500 Hz 5000 Hz 10000 Hz Max 0.4 0.4 0.4 -3 0.2 0.2 0.2 -6 0.0 0.0 0.0 -9 − 0.2 − 0.2 − 0.2 -12 − 0.4 − 0.4 − 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 x/m Figure 6. Octave band source maps for 2.5, 5 and 10 kHz obtained with the spiral array and RBF interpolation between the microphones. 3.3. Computation Time The computation of the beamforming results was done on a Intel I5-7200U CPU. All results have been calculated with the recent version of the open-source software Acoular . The calculation time for the beamforming without motion compensation was 90 s. Using the motion compensation with the ring array and linear interpolation increased the calculation time to 126 s. For the motion compensation with the Delaunay triangulation the calculation time was 190 s and for the radial basis function interpolation 545 s. While the delaunay triangulation is only 100 s slower than the beamforming without motion compensation, the RBF interpolation takes about ﬁve times longer. This is mainly because the triangulation is calculated only once but the linear system from the RBF interpolation has to be solved for every time step. y/m y/m y/m L /dB L /dB L /dB p p p Acoustics 2020, 2 337 4. Experimental Validation 4.1. Setup In order to test the different interpolation techniques acoustic measurements of a ﬁve-bladed fan have been conducted. The axial fan has a diameter of 800 mm and a rated rotational speed of 1020 revolutions per minute. The array parameters of the spiral array that is used for the measurement are the same as described in Section 2.3. The array aperture of the array is d = 1.5 m. The type of microphones used for spiral both arrays is GRAS 40PL-1 Short CCP and the distance between the array center and the fan was 0.74 m. In addition to the collected sound pressures, the rotational speed of the fan was tracked using a laser trigger. The trigger-per-revolution signal was synchronised to the pressure measurement to take into account small variations in the rotational speed of the fan. The measurement time was 40 s with a sampling frequency of 51,200 Hz. The beamforming output was calculated using conventional beamforming in the frequency domain as well as DAMAS deconvolution. Diagonal removal was applied to the cross spectral matrix. Diagonal removal removes the main diagonal elements of the cross spectral matrix to neglect the contribution of noise contamination which is incoherent for all the array microphones, e.g., self-induced ﬂow noise . The measurement was conducted in the large anechoic chamber at TU Berlin. Figure 7 shows a photograph of the experimental setup and Figure 8 shows a sketch of the setups side view as well as the front view of the fan. The array measurement and data processing parameters are summarised in Table 4. Table 4. Experimental data acquisition and processing parameters. Number of microphones 63 Spiral array aperture 1.5 m Sampling rate 51.2 kHz FFT block size 1024 FFT window/overlap von Hann/50% Focus grid distance 0.74 m Focus grid points 61 61 Figure 7. Photograph of the fan measurement setup in the anechoic chamber at TU Berlin. Acoustics 2020, 2 338 rotational axis ∅800 mm 740 mm microphone fan plane array plane Ø186 mm (a) (b) Figure 8. Sketch of the fan measurement setup. (a) shows the side view of the setup and (b) shows the front view of the fan. 4.2. Results The Figures 9 and 10 show the source distribution in terms of the sound pressure level at the array center for three different one-third octave bands that are calculated with conventional beamforming. The displayed dynamic range of the sound pressure level is set to 15 dB. For all three frequency bands the source positions at the tip of the rotor blade are the same for both interpolation algorithms. Previous studies on fan noise localisation have already shown that these source locations are plausible [26,27]. At high ﬂow rate coefﬁcients, as in this setup without static pressure differences between suction and the pressure side of the fan, the main noise sources are at the tips of the rotor blade. At 2.5 and 5 kHz both methods produce a similar beamforming output. The dynamic at the center of the grid is higher with RBF interpolation for 2.5 kHz. At the highest frequency band, sources at the center of the grid are visible. The spatial resolution of the source maps calculated with conventional beamforming is degraded by the inﬂuence of the point spread function of the microphone array. The absolute values of the conventional beamforming are overestimated because of the convoluted map. DAMAS deconvolution is applied to predict a correct source strength and enhance the localization. The source maps obtained with DAMAS deconvolution are shown in the Figures 11 and 12. The dynamic range of the maps is increased to 30 dB. For the 2.5 kHz band the distribution of the sources along the 5 blades is visible for both methods. The accuracy of the source positions at the blade tips seems more plausible using the barycentric interpolation. At 5 kHz the sources at the blade tips are reconstructed equally well with both interpolation methods. At the 10 kHz band sources in the center of the focus grid remain visible after deconvolution. The effect is stronger for the radial basis function interpolation. The localisation of the sources at the blade tips also is in better accordance with previous results using barycentric interpolation at this frequency. In order to compare the quantitative results of the interpolation methods one-third octave band spectra were computed by integrating the sound pressure over the blade areas as well as over the whole focus area. The spectra are shown in Figure 13a using barycentric and in Figure 13b using radial basis function interpolation. Both spectra are based on the DAMAS deconvolution results. The contributions of the individual blades are shown in blue. The sum of the 5 blades is shown in green. The sound pressure level over the whole focus area is shown in black. For both interpolation algorithms the individual blades equally contribute to the overall sound pressure level. Moreover, the blades act as the main noise source, except for very low and very high frequencies. For high frequencies this stems 5x 72˚ Ø786 mm Acoustics 2020, 2 339 from the sources at the center of the calculation area. For the lower frequencies the differences in the sound pressure level occur due to the insufﬁcient spatial resolution of the source localisation. 43 26 2500 Hz 5000 Hz 10000 Hz 0.4 0.4 0.4 0.2 0.2 0.2 37 20 0.0 0.0 0.0 − 0.2 − 0.2 − 0.2 31 14 − 0.4 − 0.4 − 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 x/m Figure 9. One-third octave band source maps for 2.5, 5 and 10 kHz obtained with the spiral array and barycentric interpolation between the microphones using conventional beamforming. 2500 Hz 5000 Hz 10000 Hz 0.4 0.4 0.4 41 25 0.2 0.2 0.2 38 22 0.0 0.0 0.0 − 0.2 − 0.2 − 0.2 32 16 − 0.4 − 0.4 − 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 x/m Figure 10. One-third octave band source maps for 2.5, 5 and 10 kHz obtained with the spiral array and cubic radial basis function interpolation between the microphones using conventional beamforming. 37 30 21 2500 Hz 5000 Hz 10000 Hz 34 27 18 0.4 0.4 0.4 31 24 15 28 21 12 0.2 0.2 0.2 25 18 9 0.0 22 0.0 15 0.0 6 19 12 3 − 0.2 − 0.2 − 0.2 16 9 0 13 6 − 3 − 0.4 − 0.4 − 0.4 10 3 − 6 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 x/m Figure 11. One-third octave band source maps for 2.5, 5 and 10 kHz obtained with the spiral array and barycentric interpolation between the microphones using DAMAS. y/m y/m y/m L /dB L /dB L /dB p p L /dB Lp/dB Lp/dB L /dB L /dB L /dB p p p Acoustics 2020, 2 340 37 30 21 2500 Hz 5000 Hz 10000 Hz 34 27 18 0.4 0.4 0.4 31 24 15 28 21 12 0.2 0.2 0.2 25 18 9 0.0 22 0.0 15 0.0 6 19 12 3 − 0.2 − 0.2 − 0.2 16 9 0 13 6 − 3 − 0.4 − 0.4 − 0.4 10 3 − 6 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 − 0.4 − 0.2 0.0 0.2 0.4 x/m Figure 12. One-third octave band source maps for 2.5, 5 and 10 kHz obtained with the spiral array and cubic radial basis function interpolation between the microphones using DAMAS. 70 70 focus area individual blades focus area individual blades all blades all blades 60 60 50 50 40 40 30 30 20 20 1000 2000 4000 8000 16000 1000 2000 4000 8000 16000 f/Hz f/Hz (a) (b) Figure 13. One-third octave band spectra for all blades (green), the individual blades (blue) and the whole focus area (black) obtained with the spiral array and DAMAS deconvolution. (a) Using barycentric interpolation and (b) using radial basis function interpolation. 5. Summary An extension to the original virtual rotating array method has been proposed to overcome the limitation to use only ring arrays. The developed method is able to virtually rotate microphone array geometries with arbitrary microphone arrangements. Two different interpolation methods have been shown. The ﬁrst one is based on Delaunay triangulation in combination with barycentric interpolation and the second one on radial basis function interpolation. The methods have been tested on simulated benchmark data and have been compared with the original VRA method using a ring array. The location of the simulated sources could be recovered with both methods, but the interpolation bases on Delauney triangulation archived favourable results. Measurements with a spiral array have been performed on a fan. The experimental data were evaluated with the two interpolation methods. Beamforming in frequency domain as well as DAMAS deconvolution has been applied to generate source maps and spectra. This study has shown that both methods have been able to compensate the motion of rotating sources and reconstruct the original source positions and levels with similar quality. Author Contributions: Funding acquisition, E.S.; Methodology, S.J. and E.S.; Software, S.J.; Supervision, E.S.; Writing—original draft, S.J.; Writing—review and editing, E.S. All authors have read and agreed to the published version of the manuscript. Acknowledgments: The authors thankfully acknowledge the support of this research by the German Research Foundation. Conﬂicts of Interest: The authors declare no conﬂict of interest. y/m Lp/dB L /dB L /dB L /dB L /dB p Acoustics 2020, 2 341 Abbreviations The following abbreviations are used in this manuscript: ROSI Rotating Source Identiﬁer MD modal decomposition VRA virtual rotating array DAMAS Deconvolution Approach for the Mapping of Acoustic Sources NNLS non negative least squares RBF radial basis function References 1. Sijtsma, P.; Oerlemans, S.; Holthusen, H. Location of rotating sources by phased array measurements. In Proceedings of the 7th AIAA/CEAS Aeroacoustics Conference and Exhibit, Maastricht, The Netherlands, 28–30 May 2001. [CrossRef] 2. Gomes, J. Noise Source Identiﬁcation with Blade Tracking on a Wind Turbine. Internoise 2012, 2012, 5115–5126. 3. Zhang, X.; Chu, Z.; Yang, Y.; Zhao, S.; Yang, Y. An Alternative Hybrid Time-Frequency Domain Approach Based on Fast Iterative Shrinkage-Thresholding Algorithm for Rotating Acoustic Source Identiﬁcation. IEEE Access 2019, 7, 59797–59805. [CrossRef] 4. Debrouwere, M.; Angland, D. Airy pattern approximation of a phased microphone array response to a rotating point source. J. Acoust. Soc. Am. 2017, 141, 1009–1018. [CrossRef] [PubMed] 5. Dougherty, R.P.; Walker, B. Virtual Rotating Microphone Imaging of Broadband Fan Noise. In Proceedings of the 15th AIAA/CEAS Aeroacoustics Conference, Miami, FL, USA, 11–13 May 2009. 6. Lowis, C.; Joseph, P. Determining the strength of rotating broadband sources in ducts by inverse methods. J. Sound Vib. 2006, 295, 614–632. [CrossRef] 7. Pannert, W.; Maier, C. Rotating beamforming - Motion-compensation in the frequency domain and application of high-resolution beamforming algorithms. J. Sound Vib. 2014, 333, 1899–1912. [CrossRef] 8. Herold, G.; Sarradj, E. Microphone array method for the characterization of rotating sound sources in axial fans. Noise Control. Eng. J. 2015, 63, 546–551. [CrossRef] 9. Kotán, G.; Tóth, B.; Vad, J. Comparison of the rotating source identiﬁer and the virtual rotating array method. Period. Polytech. Mech. Eng. 2018, 62, 261–268. [CrossRef] 10. Wei, M.; Huan, B.; Zhang, C.; Liu, X. Beamforming of Phased Microphone Array for Rotating. J. Sound Vib. 2019, 467. [CrossRef] 11. Herold, G.; Ocker, C.; Sarradj, E.; Pannert, W. A comparison of microphone array methods for the characterization of rotating sound sources. In Proceedings of the 7th Berlin Beamforming Conference, Berlin, Germany, 5–6 March 2018; pp. 1–12. 12. Tóth, B.; Kalmár-Nagy, T.; Vad, J. Rotating beamforming with uneven microphone placements. In Proceedings of the 7th Berlin Beamforming Conference, Berlin, Germany, 5–6 March 2018; pp. 1–13. 13. Sarradj, E. Three-dimensional acoustic source mapping with different beamforming steering vector formulations. Adv. Acoust. Vib. 2012, doi:10.1155/2012/292695. [CrossRef] 14. Brooks, T.F.; Humphreys, W.M. A deconvolution approach for the mapping of acoustic sources (DAMAS) determined from phased microphone arrays. J. Sound Vib. 2006, 294, 856–879. [CrossRef] 15. Lawson, C.L. Solving Least Squares Problems, 1974th ed.; SIAM: Philadelphia, PA, USA, 1995. 16. Delaunay, B. Sur la sphère vide. Cl. Des Sci. Math. 1934, 6, 793–800. 17. Barber, C.B.; Dobkin, D.P.; Huhdanpaa, H. The Quickhull Algorithm for Convex Hulls. ACM Trans. Math. Softw. 1996, 22, 469–483. [CrossRef] 18. Vogel, H. A Better Way to Construct the Sunﬂower. Math. Biosci. 1979, 44, 179–189. [CrossRef] 19. Sarradj, E. A Generic Approach To Synthesize Optimal Array Microphone Arrangements. In Proceedings of the 6th Berlin Beamforming Conference, Berlin, Germany, 29 February–1 March 2016; pp. 1–12. 20. Buhmann, M.D. Radial Basis Functions: Theory and Implementations; Cambridge University Press: Cambridge, UK, 2003; Volume 12. Acoustics 2020, 2 342 21. Sarradj, E.; Herold, G.; Sijtsma, P.; Merino Martinez, R.; Geyer, T.F.; Bahr, C.J.; Porteous, R.; Moreau, D.; Doolan, C.J. A Microphone Array Method Benchmarking Exercise using Synthesized Input Data. In Proceedings of the 23rd AIAA/CEAS Aeroacoustics Conference, Denver, CO, USA, 5–9 June 2017. [CrossRef] 22. Herold, G. Microphone Array Benchmark b11: Rotating Point Sources. 2017. Available online: https: //depositonce.tu-berlin.de/handle/11303/9404 (accessed on 2 March 2020). 23. Johnson, D.H.; Dudgeon, D.E. Array Signal Processing: Concepts and Techniques; P T R Prentice Hall: Upper Saddle River, NJ, USA, 1993. 24. Sarradj, E.; Herold, G. A Python framework for microphone array data processing. Appl. Acoust. 2017, 116, 50–58. [CrossRef] 25. Hald, J. Removal of incoherent noise from an averaged cross-spectral matrix. J. Acoust. Soc. Am. 2017, 142, 846–854. [CrossRef] [PubMed] 26. Zenger, F.J.; Herold, G.; Becker, S.; Sarradj, E. Sound source localization on an axial fan at different operating points. Exp. Fluids 2016, 57, 1–10. [CrossRef] 27. Herold, G.; Zenger, F.; Sarradj, E. Inﬂuence of blade skew on axial fan component noise. Int. J. Aeroacoustics 2017, 16, 418–430. [CrossRef] c 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.pngAcousticsMultidisciplinary Digital Publishing Institutehttp://www.deepdyve.com/lp/multidisciplinary-digital-publishing-institute/an-extension-of-the-virtual-rotating-array-method-using-arbitrary-jFuwTEbnau
An Extension of the Virtual Rotating Array Method Using Arbitrary Microphone Configurations for the Localization of Rotating Sound Sources