# Optimization of Equivalent Source Configuration for an Independent-Equivalent Source Method in Half-Space Sound Field

Authors: Wen-Qian Jing; Huawei Wu; Jin-Quan Nie
Journal: Shock and Vibration
 (2020)

---

## Abstract

In the situation that vibrating objects are located above a reflecting plane, an independent-equivalent source method (I-ESM) regards the reflections due to the plane as being radiated by equivalent sources placed under the plane and then the half-space sound field is reconstructed by matching the measured pressure with the equivalent sources distributed within the vibrating object and those substituting for reflections. But, this method heavily depends on the equivalent source configuration and may obtain bad reconstruction results if the equivalent sources are arranged incorrectly. This paper deals with the optimization of the equivalent source configuration to ensure I-ESM always perform well. Through numerical simulations and experiments, the influence of equivalent source configurations on the reconstruction accuracy was studied and optimal configurations were acquired and confirmed.


---

## Body

## 1. Introduction

In practice, most vibrating objects are mounted on or located above a plane, resulting in a half-space sound field, and then the reflections from that plane should be considered when using the nearfield acoustic holography (NAH) [1–5] to realize the sound field reconstruction.By combining with the conventional NAH, the field separation techniques [6–8] can be used to reconstruct the half-space sound field. Also, these techniques can work in the situation that the surface impedance of the reflecting plane is unknown, but they require that the measurement surfaces should surround the vibrating object at a close range. In addition, the field separation techniques require both the pressure and particle velocity measurements if the measurements are made on one single-layer surface [6], or require double-layer measurements if only pressure [7] or particle velocity [8] is measured. As a result, high measurement costs would be needed.On the condition that the surface impedance is knowna prior, half-space Green’s function can be introduced into the conventional NAH [9–13] to account for the reflections. Both in [10, 11], the reflections were assumed to be as plane waves, which might not obtain a good result if the vibrating object was not far away from the reflecting plane. Based on this, half-space Green’s function [14, 15] that automatically satisfies the boundary condition of the reflecting plane was introduced into the equivalent source method- (ESM-) based NAH, and this method performed well [16]. But, this method depends on the surface impedance of the reflecting plane and is named dependence (D)-ESM.To develop a method for the situations that the surface impedance is unknown and to avoid high measurement costs, the reflections due to the reflecting plane can be regarded as being radiated by equivalent sources located under the plane and then both the direct sound and reflected sound can be expressed as the superposition of a series of equivalent sources [16]. Because this method is based on ESM and is independent of the surface impedance, it is named independent (I)-ESM. It is predictable that I-ESM may give bad reconstruction results and even fail in reconstructing if the equivalent sources are arranged improperly because the equivalent source configuration affects ESM heavily.In [16], two types of configurations of equivalent sources, inside the image source (IIS) and close to the reflecting plane (CRP), are considered to substituting for the reflections in I-ESM. If IIS configuration is used, the reflections are regarded as being radiated by a source with the same geometry of the vibrating object, and this source and the vibrating object are placed symmetrically with respect to the reflecting plane. If CRP configuration is used, the reflections are regarded as being radiated by a planar source coincident with the reflecting plane. The results show that the I-ESM with IIS configuration performs well, whereas the I-ESM with CRP configuration cannot provide stable reconstruction accuracy. This is probably because the equivalent sources substituting for the reflections are set improperly when using CRP configuration, as predicated above.This paper deals with the optimization of the equivalent source configuration to ensure I-ESM always performs well. Because the problem focuses on the reflected sound, the configuration of equivalent sources substituting for the direct sound is not considered here. The configuration mentioned in the following just considers the equivalent sources substituting for the reflections.This paper realizes the optimization of equivalent source configuration for I-ESM as follows. Firstly, analyze and discuss the influence of CRP configurations on the reconstruction accuracy to search for the appropriate CRP configurations. Then, study the relation between the appropriate CRP configuration and the IIS configuration and acquire the optimal configurations. Finally, confirm the optimal equivalent source configurations.


## 2. Outline of Theory

According to the idea of I-ESM [16], the sound field above the reflecting plane can be approximated by the superposition of the fields generated by the equivalent sources on a fictitious surface Γ within the vibrating object and those on a fictitious surface Ω substituting for the reflections, as shown in Figure 1. Given a field point r, the pressure can be expressed as(1)pr=iρω∑qiΓgpr,riΓ+∑qiΩgpr,riΩ,where i means the imaginary, ρ is the density of air, ω is the angular frequency, qiΓ and qjΩ are the strength of the ith and jth equivalent source on Γ and Ω, respectively, riΓ and rjΩ are the position of the ith and jth equivalent source, respectively, and Green’s function can be expressed as(2)gpr,riΓ=eikRiΓ4πRiΓ,gpr,rjΩ=eikRjΩ4πRjΩ,where k is the wave number, RiΓ is the distance between r and riΓ, and RjΩ is the distance between r and rjΩ.Figure 1
Positions of the spherical sourceS, the hologram surface H, and the fictitious surface Γ and Ω. (a) Top view. (b) Side view.
(a)
(b)Given a hologram surfaceH, equation (1) can be written in a matrix form as(3)PH=iρωGpΓHQΓ+GpΩHQΩ=iρωGpΣHQΣ,with(4)GpΣH=GpΓHGpΩH,QΣ=QΓQΩT,where GpΓHGpΩH are the transfer matrix relating the pressures on H to the equivalent sources on Γ and Ω, respectively, and QΓ and QΩ are the equivalent source strength vector on Γ and Ω, respectively. Subsequently, the regularized solution of QΣ can be obtained as(5)QΣ=1iρωGpΣHHGpΣH+εE−1GpΣHHPH,where the superscript “H” denotes the Hermitian transpose, “−1” denotes the inverse matrix, ε is the regularization parameter, and E is a unit matrix.By substitutingQΣ into equation (3), the pressure at any field point in half-space can be reconstructed. Also, the normal velocity on the surface of the vibrating object can also be reconstructed by using the following formula:(6)VS=GvΣSQΣ,where GvΣS is the transfer matrix relating the normal velocities on S to the equivalent sources on Γ and Ω, and it can be expressed as(7)GvΣS=GvΓSGvΩS,where GvΓS and GvΩS are, respectively, constituted of(8)gvrnS,riΓ=1−ikRjnΓSeikRinΓS4πRinΓS2cosϕinΓS,cosϕinΓS=rnS−riΓ⋅nSrnS−riΓ,(9)gvrnS,rjΩ=1−ikRjnΩSeikRjnΩS4πRjnΩS2cosϕjnΩS,cosϕjnΩS=rnS−rjΩ⋅nSrnS−rjΩ,where “⋅” denotes the dot product operator, nS is the unit normal vector on the surface of the vibrating object, RinΓS is the distance between the nth surface node rnS and riΓ, and RjnΩS is the distance between rnS and rjΩ.


## 3. Numerical Simulations

This section explains optimization of the configuration of equivalent sources substituting for reflections, i.e., the equivalent sources on the fictitious surface Ω, a planar surface coincident with the reflecting plane, as shown in Figure1.


### 3.1. Parameters in Numerical Simulations

A spherical source with a radius of 0.1 m was placed above a reflecting plane. The reflecting plane was assumed to be of the Delany and Bazley type [17] and to be infinitely thick, and then the surface impedance Z was calculated by using the following formula:(10)Z=ρc1+9.08fσ0−0.75+i×11.9fσ0−0.73,where f is the frequency and σ0 is the flow resistivity in cgs units. Note that this equation does not take into account the incidence effect; thus, the reflecting plane considered here is of local reaction.The center of the spherical source was located 0.5 m above the reflecting plane and the hologram planeH was located at y = 0.2 m, as shown in Figure 1. The dimensions of the hologram plane were 0.5 m × 0.5 m with measurement points distributed from −0.25 m to 0.25 m in the x direction and from 0.25 m to 0.75 m in the z direction, and the sampling interval was 0.05 m in both x and z directions.The nodes on the source surface were evenly distributed with intervals ofπ/4 in the azimuth angle direction and π/6 in the polar angle direction. The equivalent sources substituting for the direct sound were placed on a smaller concentric spherical surface Γ inside the spherical source S with a radius of 0.02 m, which was less than 0.4 times the radius of the spherical source [18]. The distribution of equivalent sources on Γ was the same as that of surface nodes.The normal velocity on the spherical source surface was reconstructed to optimize the equivalent source configuration for I-ESM. The theoretical normal velocity is(11)vS=v0zS−zara,where v0 is the uniform velocity in the radial direction, ra is the radius of the spherical source, and zS and za are the z coordinates of the source node and the center of the source sphere, respectively. The holographic pressure was calculated numerically by using the boundary element method codes [19]. Also, the surface impedance of the reflecting plane was set as 35 cgs units. Besides, the white noise with a signal-to-noise ratio of 30 dB was added to the holographic pressure.


### 3.2. Discussion on the Equivalent Source Configuration

When using CRP configuration, the retreat distance of Ω from the reflecting plane,hz shown in Figure 1(b), should be considered. Here, hz is set as −0.001 m, −0.2 m, −0.5 m, and −1 m, respectively. And the dimension of Ω and the distribution interval of equivalent sources should be considered. The dimension in the x direction is indicated by the minimum and maximum x coordinates, x1 and x2, as shown in Figure 1(a). Similarly, y1 and y2 indicate the dimension in the y direction. The distribution intervals are indicated by dx and dy. Here, 22 groups of parameters x1,x2,y1,y2,dx,dy are considered and they are given in Table 1.Table 1
Parameters groupsx1,x2,y1,y2,dx,dy in numerical simulations.
Number
x
1
,
x
2 (m)
y
1
,
y
2 (m)
d
x (m)
d
y (m)
1
(−0.2, 0.2)
(−0.2, 0.2)
0.050
0.050
2
(−0.2, 0.2)
(−0.2, 0.2)
0.025
0.025
3
(−0.1, 0.1)
(−0.1, 0.1)
0.050
0.050
4
(−0.1, 0.1)
(−0.1, 0.1)
0.025
0.025
5
(−0.1, 0.1)
(−0.1, 0.1)
0.010
0.010
6
(−0.05, 0.05)
(−0.05, 0.05)
0.050
0.050
7
(−0.05, 0.05)
(−0.05, 0.05)
0.025
0.025
8
(−0.05, 0.05)
(−0.05, 0.05)
0.010
0.010
9
(−0.25, 0.25)
(0.1, 0.2)
0.050
0.050
10
(−0.25, 0.25)
(0.1, 0.2)
0.050
0.025
11
(−0.25, 0.25)
(0, 0.2)
0.050
0.050
12
(−0.25, 0.25)
(0, 0.2)
0.050
0.025
13
(−0.25, 0.25)
(−0.1, 0.2)
0.050
0.050
14
(−0.25, 0.25)
(−0.1, 0.2)
0.050
0.025
15
(−0.25, 0.25)
(−0.2, 0.2)
0.050
0.050
16
(−0.25, 0.25)
(−0.2, 0.2)
0.050
0.025
17
(−0.25, 0.25)
(−0.1, 0.1)
0.050
0.050
18
(−0.25, 0.25)
(−0.1, 0.1)
0.050
0.025
19
(−0.25, 0.25)
(−0.1, 0.1)
0.050
0.010
20
(−0.25, 0.25)
(−0.05, 0.05)
0.050
0.050
21
(−0.25, 0.25)
(−0.05, 0.05)
0.050
0.025
22
(−0.25, 0.25)
(−0.05, 0.05)
0.050
0.010With different CRP configurations, the reconstruction errors of normal velocity on the spherical source surface at 500 Hz were obtained by using I-ESM and they are given in Figure2.(1)
Whenhz takes −0.001 m, −0.2 m, and −0.5 m successively, it is obvious that the reconstruction error decreases with the decrease of hz. However, the error does not get smaller when hz takes smaller value −1 m and even increases when taking no. 3–12 groups of parameters. Thus, −0.5 m is a suitable choice as the retreat distance. It is interesting that the spherical source was located right 0.5 m above the reflecting plane. So, it can be concluded that the retreat distance below the plane should approximate with the distance of the spherical source above the plane, i.e., Ω should be placed at the distance of the source image.(2)
In the casehz<−0.5m, i.e., hz=−1m, Ω is beyond the distance of the source image. It can be seen that the reconstruction error changes little versus the group of parameters and generally can be acceptable.(3)
In the casehz>−0.5m, i.e., hz=−0.001m and hz=−0.2m, Ω is within the distance of the source image. It can be seen that the two error curves are consistent, showing the same turning points. From Table 1, it can be found that the turning points are just the cases that the dimension x1,x2,y1,y2 changes. But, the error changes little with the distribution interval dx,dy. This shows that the dimension has heavy influence on the error, while the distribution interval affects the error little.Figure 2
The reconstruction errors of normal velocity at 500 Hz when using different CRP configurations.Besides, it can be seen that the error was the smallest when taking No. 9 group of parameters and was the largest when taking No. 15 group of parameters. Figure3 gives Ω9 and Ω15, which are, respectively, drawn by No. 9 and No. 15 group of parameters. It can be seen that Ω9 just covers the area between the source and hologram plane and Ω15 is the largest surface covering Ω9. This indicates that the equivalent sources should be arranged just in the area between the source and hologram plane when Ω is within the distance of the source image.(4)
In the casehz=−0.5m, Ω is right at the distance of the source image. It can be found that the reconstruction error is quite small when taking No. 3, No. 6, No. 9, No. 11, or No. 20 group of parameters. According to Figure 3, Ω3 and Ω20 just cover the spherical source and Ω6 gets smaller covering the central part of the spherical source. Also, Ω9 just covers the area between the source and the hologram plane and Ω11 also covers this area though with a little larger dimension.Figure 3
Geometric sketch of the fictitious surface Ω.In summary, the equivalent sources substituting for the reflections should be placed at either the position of the source image or the area between the source and hologram plane. How interesting it is that the equivalent sources distribute just within the source image when using IIS configuration. Thus, it is expectable that the I-ESM with IIS configuration can perform well.


### 3.3. Confirmation of the Optimal Equivalent Source Configuration

To confirm the optimal equivalent source configuration, Table2 gives three sets of configurations for further investigation. The first set (hz=−0.001m, No. 15) represents the inappropriate configuration, while the third set (hz=−0.5m, No. 6) represents the optimal configuration. Also, the second set (hz=−0.2m, No. 1) represents an acceptable configuration.Table 2
Equivalent source configurations for confirmation.
Retreat distancehz (m)
Number of parameter groupx1,x2,y1,y2,dx,dy
−0.001
15
−0.2
1
−0.5
6With the configurations given in Table2, the reconstruction errors of normal velocity were obtained by using I-ESM. Because the D-ESM performs well, it was also used to obtain the reconstruction error, and the error acts as the reference.Figure4 gives the reconstruction error of normal velocity versus the frequency, and Figure 5 gives the reconstruction error versus the flow resistivity at 500 Hz. Both figures illustrate that the errors obtained by using I-ESM with hz=−0.5m are comparable to the reference errors, while the errors obtained by using I-ESM with hz=−0.001m are much larger than those. This confirms that the optimal configuration, placing the equivalent sources at the position of source image, is indeed a good choice.Figure 4
The reconstruction error of normal velocity versus the frequency.Figure 5
The reconstruction error of normal velocity versus the flow resistivity.


## 3.1. Parameters in Numerical Simulations

A spherical source with a radius of 0.1 m was placed above a reflecting plane. The reflecting plane was assumed to be of the Delany and Bazley type [17] and to be infinitely thick, and then the surface impedance Z was calculated by using the following formula:(10)Z=ρc1+9.08fσ0−0.75+i×11.9fσ0−0.73,where f is the frequency and σ0 is the flow resistivity in cgs units. Note that this equation does not take into account the incidence effect; thus, the reflecting plane considered here is of local reaction.The center of the spherical source was located 0.5 m above the reflecting plane and the hologram planeH was located at y = 0.2 m, as shown in Figure 1. The dimensions of the hologram plane were 0.5 m × 0.5 m with measurement points distributed from −0.25 m to 0.25 m in the x direction and from 0.25 m to 0.75 m in the z direction, and the sampling interval was 0.05 m in both x and z directions.The nodes on the source surface were evenly distributed with intervals ofπ/4 in the azimuth angle direction and π/6 in the polar angle direction. The equivalent sources substituting for the direct sound were placed on a smaller concentric spherical surface Γ inside the spherical source S with a radius of 0.02 m, which was less than 0.4 times the radius of the spherical source [18]. The distribution of equivalent sources on Γ was the same as that of surface nodes.The normal velocity on the spherical source surface was reconstructed to optimize the equivalent source configuration for I-ESM. The theoretical normal velocity is(11)vS=v0zS−zara,where v0 is the uniform velocity in the radial direction, ra is the radius of the spherical source, and zS and za are the z coordinates of the source node and the center of the source sphere, respectively. The holographic pressure was calculated numerically by using the boundary element method codes [19]. Also, the surface impedance of the reflecting plane was set as 35 cgs units. Besides, the white noise with a signal-to-noise ratio of 30 dB was added to the holographic pressure.


## 3.2. Discussion on the Equivalent Source Configuration

When using CRP configuration, the retreat distance of Ω from the reflecting plane,hz shown in Figure 1(b), should be considered. Here, hz is set as −0.001 m, −0.2 m, −0.5 m, and −1 m, respectively. And the dimension of Ω and the distribution interval of equivalent sources should be considered. The dimension in the x direction is indicated by the minimum and maximum x coordinates, x1 and x2, as shown in Figure 1(a). Similarly, y1 and y2 indicate the dimension in the y direction. The distribution intervals are indicated by dx and dy. Here, 22 groups of parameters x1,x2,y1,y2,dx,dy are considered and they are given in Table 1.Table 1
Parameters groupsx1,x2,y1,y2,dx,dy in numerical simulations.
Number
x
1
,
x
2 (m)
y
1
,
y
2 (m)
d
x (m)
d
y (m)
1
(−0.2, 0.2)
(−0.2, 0.2)
0.050
0.050
2
(−0.2, 0.2)
(−0.2, 0.2)
0.025
0.025
3
(−0.1, 0.1)
(−0.1, 0.1)
0.050
0.050
4
(−0.1, 0.1)
(−0.1, 0.1)
0.025
0.025
5
(−0.1, 0.1)
(−0.1, 0.1)
0.010
0.010
6
(−0.05, 0.05)
(−0.05, 0.05)
0.050
0.050
7
(−0.05, 0.05)
(−0.05, 0.05)
0.025
0.025
8
(−0.05, 0.05)
(−0.05, 0.05)
0.010
0.010
9
(−0.25, 0.25)
(0.1, 0.2)
0.050
0.050
10
(−0.25, 0.25)
(0.1, 0.2)
0.050
0.025
11
(−0.25, 0.25)
(0, 0.2)
0.050
0.050
12
(−0.25, 0.25)
(0, 0.2)
0.050
0.025
13
(−0.25, 0.25)
(−0.1, 0.2)
0.050
0.050
14
(−0.25, 0.25)
(−0.1, 0.2)
0.050
0.025
15
(−0.25, 0.25)
(−0.2, 0.2)
0.050
0.050
16
(−0.25, 0.25)
(−0.2, 0.2)
0.050
0.025
17
(−0.25, 0.25)
(−0.1, 0.1)
0.050
0.050
18
(−0.25, 0.25)
(−0.1, 0.1)
0.050
0.025
19
(−0.25, 0.25)
(−0.1, 0.1)
0.050
0.010
20
(−0.25, 0.25)
(−0.05, 0.05)
0.050
0.050
21
(−0.25, 0.25)
(−0.05, 0.05)
0.050
0.025
22
(−0.25, 0.25)
(−0.05, 0.05)
0.050
0.010With different CRP configurations, the reconstruction errors of normal velocity on the spherical source surface at 500 Hz were obtained by using I-ESM and they are given in Figure2.(1)
Whenhz takes −0.001 m, −0.2 m, and −0.5 m successively, it is obvious that the reconstruction error decreases with the decrease of hz. However, the error does not get smaller when hz takes smaller value −1 m and even increases when taking no. 3–12 groups of parameters. Thus, −0.5 m is a suitable choice as the retreat distance. It is interesting that the spherical source was located right 0.5 m above the reflecting plane. So, it can be concluded that the retreat distance below the plane should approximate with the distance of the spherical source above the plane, i.e., Ω should be placed at the distance of the source image.(2)
In the casehz<−0.5m, i.e., hz=−1m, Ω is beyond the distance of the source image. It can be seen that the reconstruction error changes little versus the group of parameters and generally can be acceptable.(3)
In the casehz>−0.5m, i.e., hz=−0.001m and hz=−0.2m, Ω is within the distance of the source image. It can be seen that the two error curves are consistent, showing the same turning points. From Table 1, it can be found that the turning points are just the cases that the dimension x1,x2,y1,y2 changes. But, the error changes little with the distribution interval dx,dy. This shows that the dimension has heavy influence on the error, while the distribution interval affects the error little.Figure 2
The reconstruction errors of normal velocity at 500 Hz when using different CRP configurations.Besides, it can be seen that the error was the smallest when taking No. 9 group of parameters and was the largest when taking No. 15 group of parameters. Figure3 gives Ω9 and Ω15, which are, respectively, drawn by No. 9 and No. 15 group of parameters. It can be seen that Ω9 just covers the area between the source and hologram plane and Ω15 is the largest surface covering Ω9. This indicates that the equivalent sources should be arranged just in the area between the source and hologram plane when Ω is within the distance of the source image.(4)
In the casehz=−0.5m, Ω is right at the distance of the source image. It can be found that the reconstruction error is quite small when taking No. 3, No. 6, No. 9, No. 11, or No. 20 group of parameters. According to Figure 3, Ω3 and Ω20 just cover the spherical source and Ω6 gets smaller covering the central part of the spherical source. Also, Ω9 just covers the area between the source and the hologram plane and Ω11 also covers this area though with a little larger dimension.Figure 3
Geometric sketch of the fictitious surface Ω.In summary, the equivalent sources substituting for the reflections should be placed at either the position of the source image or the area between the source and hologram plane. How interesting it is that the equivalent sources distribute just within the source image when using IIS configuration. Thus, it is expectable that the I-ESM with IIS configuration can perform well.


## 3.3. Confirmation of the Optimal Equivalent Source Configuration

To confirm the optimal equivalent source configuration, Table2 gives three sets of configurations for further investigation. The first set (hz=−0.001m, No. 15) represents the inappropriate configuration, while the third set (hz=−0.5m, No. 6) represents the optimal configuration. Also, the second set (hz=−0.2m, No. 1) represents an acceptable configuration.Table 2
Equivalent source configurations for confirmation.
Retreat distancehz (m)
Number of parameter groupx1,x2,y1,y2,dx,dy
−0.001
15
−0.2
1
−0.5
6With the configurations given in Table2, the reconstruction errors of normal velocity were obtained by using I-ESM. Because the D-ESM performs well, it was also used to obtain the reconstruction error, and the error acts as the reference.Figure4 gives the reconstruction error of normal velocity versus the frequency, and Figure 5 gives the reconstruction error versus the flow resistivity at 500 Hz. Both figures illustrate that the errors obtained by using I-ESM with hz=−0.5m are comparable to the reference errors, while the errors obtained by using I-ESM with hz=−0.001m are much larger than those. This confirms that the optimal configuration, placing the equivalent sources at the position of source image, is indeed a good choice.Figure 4
The reconstruction error of normal velocity versus the frequency.Figure 5
The reconstruction error of normal velocity versus the flow resistivity.


## 4. Experiments

Two experiments were carried out in a semianechoic chamber to optimize the equivalent source configuration for I-ESM. Also, two kinds of reflecting materials, a marble floor and a piece of sponge, were, respectively, used as the reflecting plane.


### 4.1. Parameters in Experiments

A loudspeaker with the size of 0.60 m × 0.15 m × 0.14 m was used as the source. Also, the positions of source and measurement surfaces are shown in Figure6. It can be seen that the loudspeaker was placed 0.18 m above the reflecting plane and the centers of two cores of the loudspeaker were located at the points (−0.15 m, 0, 0.18 m) and (0.10 m, 0, 0.18 m). The measurement planes H1 and H2 were located 0.14 m and 0.20 m away from the source plane, respectively. The dimensions of the measurement surface were 0.7 m × 0.4 m with the interval of 0.05 m in both x and z directions. The fictitious surface Γ was located at y=−0.01m with the same dimension and interval as those of measurement surfaces.Figure 6
Location of the loudspeaker, the measurement surfaces H1 and H2, and the fictitious surface Ω. (a) Top view. (b) Side view.
(a)
(b)Using the measured pressure onH2 as the input, the pressure on H1 was reconstructed by using I-ESM, and the reconstructed results were compared with the measured values on H1. Then, the pressure reconstruction error was used to realize the optimization of equivalent source configuration for I-ESM.


### 4.2. Discussion on the Equivalent Source Configuration

Table3 gives 23 groups of parameters x1,x2,y1,y2,dx,dy to show the dimension and distribution interval of Ω. Also, the retreat distance hz is set as −0.001 m, −0.18 m, −0.5 m and −1 m, respectively.Table 3
Parameters groupsx1,x2,y1,y2,dx,dy in experiments.
Number
x
1
,
x
2 (m)
y
1
,
y
2 (m)
d
x (m)
d
y (m)
1
(−0.2, 0.2)
(−0.2, 0.2)
0.050
0.050
2
(−0.2, 0.2)
(−0.2, 0.2)
0.025
0.025
3
(−0.15, 0.15)
(−0.15, 0.15)
0.050
0.050
4
(−0.15, 0.15)
(−0.15, 0.15)
0.025
0.025
5
(−0.15, 0.15)
(−0.15, 0.15)
0.010
0.010
6
(−0.1, 0.1)
(−0.1, 0.1)
0.050
0.050
7
(−0.1, 0.1)
(−0.1, 0.1)
0.025
0.025
8
(−0.1, 0.1)
(−0.1, 0.1)
0.010
0.010
9
(−0.35, 0.35)
(0, 0.2)
0.050
0.050
10
(−0.35, 0.35)
(0, 0.2)
0.050
0.025
11
(−0.35, 0.35)
(−0.15, 0.2)
0.050
0.050
12
(−0.35, 0.35)
(−0.15, 0.2)
0.050
0.025
13
(−0.35, 0.35)
(−0.2, 0.2)
0.050
0.050
14
(−0.35, 0.35)
(−0.2, 0.2)
0.050
0.025
15
(−0.35, 0.35)
(0, 0.15)
0.050
0.050
16
(−0.35, 0.35)
(0, 0.15)
0.050
0.025
17
(−0.35, 0.35)
(0, 0.15)
0.050
0.010
18
(−0.35, 0.35)
(−0.15, 0.15)
0.050
0.050
19
(−0.35, 0.35)
(−0.15, 0.15)
0.050
0.025
20
(−0.35, 0.35)
(−0.15, 0.15)
0.050
0.010
21
(−0.35, 0.35)
(−0.1, 0.1)
0.050
0.050
22
(−0.35, 0.35)
(−0.1, 0.1)
0.050
0.025
23
(−0.35, 0.35)
(−0.1, 0.1)
0.050
0.010Figure7 gives the reconstruction errors of pressure at 700 Hz obtained by using I-ESM with different CRP configurations when a marble floor was used as the reflecting plane.(1)
It is obvious that the reconstruction error is the smallest whenhz=−0.18m, showing that −0.18 m is a nearly perfect choice as the retreat distance. Also, the loudspeaker was placed right 0.18 m above the reflecting plane. This illustrates that the retreat distance below the reflecting plane should be set as the distance of the source image, which is consistent with the conclusion obtained in numerical simulations.(2)
In the casehz<−0.18m, i.e., hz=−0.5m and hz=−1m, it can be seen that the reconstruction results are rather stable, but all of them are a little large.(3)
In the casehz>−0.18m, i.e., hz=−0.001m, it can be seen that the error has large fluctuations versus the group of parameters. Each time the dimension x1,x2,y1,y2 changes, the error changes significantly. Also, the distribution interval dx,dy also affects the error much. Thus, if the dimension and distribution interval are set properly, e.g. taking No. 15 or No. 18 group of parameters, the error can be quite small. From Table 3, it can be found that Ω15 just covers the area between the source and the reconstructed plane and Ω18 also covers this area though with a little larger dimensions.(4)
In the casehz=−0.18m, it is obvious that there is an abnormal point, i.e., when taking No. 6 group of parameters. From Table 3, it can be seen that Ω6 does not cover both cores of the loudspeaker. Thus, the abnormal point is understandable.Figure 7
The reconstruction errors of pressure when using different CRP configurations with a marble floor as the reflecting plane.Figure8 gives the reconstruction errors of pressure at 700 Hz when a piece of sponge was used as the reflecting plane.(1)
It is obvious that the reconstruction error fluctuates slightly versus the group of parameters whenhz≤−0.18m.(2)
In the casehz>−0.18m, i.e., hz=−0.001m, it can be seen that the error fluctuates significantly versus the group of parameters. Also, the error was the smallest when taking No. 6 or No. 21 group of parameters and the largest when taking No. 13 or No. 14 group of parameters. From Table 3, it can be found that both Ω6 and Ω21 are very small, just covering the central part of the source, while Ω13 and Ω14 are the largest dimensions. In addition, it can be found that the fewer the equivalent sources, the better the reconstruction results. This is probably because that the highly absorptive sponge contributes little for the half-space sound field and few equivalent sources are needed.(3)
In the casehz=−0.18m, it can be found that the error is the smallest when taking No. 6 or No. 21 group of parameters. This is consistent with the conclusion presented above that Ω should just cover the central part of the source and the equivalent sources should be as few as possible.Figure 8
The reconstruction errors of pressure when using different CRP configurations with a piece of sponge as the reflecting plane.


### 4.3. Confirmation of the Optimal Equivalent Source Configuration

According to the discussion presented in Section4.2, Table 4 gives six sets of equivalent source configurations for further investigations to confirm the optimal configurations. The two configurations when hz=−0.18m represent the optimal ones. Also, other four configurations, respectively, represent the worst and the best in the case hz=−0.001m.Table 4
Equivalent source configurations for confirmation in experiments.
Reflecting plane
Retreat distancehz (m)
Number of parameter groupx1,x2,y1,y2,dx,dy
A marble floor
−0.001
5
A marble floor
−0.001
18
A marble floor
−0.18
6
A piece of sponge
−0.001
13
A piece of sponge
−0.001
6
A piece of sponge
−0.18
21Figures9 and 10 give the reconstruction errors of pressures versus the frequency obtained by using I-ESM with the CRP configurations shown in Table 4. It can be seen that I-ESM provides large and unstable errors when hz=−0.001m with No. 5 or No. 13 group of parameters. But, the errors obtained when hz=−0.18m are smaller and more stable, confirming that the distance of the source image is suitable to be set as the retreat distance.Figure 9
The reconstruction error of pressure versus the frequency with a marble floor as the reflecting plane.Figure 10
The reconstruction error of pressure versus the frequency with a piece of sponge as the reflecting plane.Note that the I-ESM performs best when taking No. 6 group of parameters but fails to reconstruct the half-space sound field at 100 Hz and 500 Hz when taking No. 13 group of parameters in the case that a piece of sponge was used as the reflecting plane andhz=−0.001m. This confirms that placing a small amount of equivalent sources very close to the reflecting plane and just in the central part of source is a good choice when the reflecting plane is highly absorptive.


## 4.1. Parameters in Experiments

A loudspeaker with the size of 0.60 m × 0.15 m × 0.14 m was used as the source. Also, the positions of source and measurement surfaces are shown in Figure6. It can be seen that the loudspeaker was placed 0.18 m above the reflecting plane and the centers of two cores of the loudspeaker were located at the points (−0.15 m, 0, 0.18 m) and (0.10 m, 0, 0.18 m). The measurement planes H1 and H2 were located 0.14 m and 0.20 m away from the source plane, respectively. The dimensions of the measurement surface were 0.7 m × 0.4 m with the interval of 0.05 m in both x and z directions. The fictitious surface Γ was located at y=−0.01m with the same dimension and interval as those of measurement surfaces.Figure 6
Location of the loudspeaker, the measurement surfaces H1 and H2, and the fictitious surface Ω. (a) Top view. (b) Side view.
(a)
(b)Using the measured pressure onH2 as the input, the pressure on H1 was reconstructed by using I-ESM, and the reconstructed results were compared with the measured values on H1. Then, the pressure reconstruction error was used to realize the optimization of equivalent source configuration for I-ESM.


## 4.2. Discussion on the Equivalent Source Configuration

Table3 gives 23 groups of parameters x1,x2,y1,y2,dx,dy to show the dimension and distribution interval of Ω. Also, the retreat distance hz is set as −0.001 m, −0.18 m, −0.5 m and −1 m, respectively.Table 3
Parameters groupsx1,x2,y1,y2,dx,dy in experiments.
Number
x
1
,
x
2 (m)
y
1
,
y
2 (m)
d
x (m)
d
y (m)
1
(−0.2, 0.2)
(−0.2, 0.2)
0.050
0.050
2
(−0.2, 0.2)
(−0.2, 0.2)
0.025
0.025
3
(−0.15, 0.15)
(−0.15, 0.15)
0.050
0.050
4
(−0.15, 0.15)
(−0.15, 0.15)
0.025
0.025
5
(−0.15, 0.15)
(−0.15, 0.15)
0.010
0.010
6
(−0.1, 0.1)
(−0.1, 0.1)
0.050
0.050
7
(−0.1, 0.1)
(−0.1, 0.1)
0.025
0.025
8
(−0.1, 0.1)
(−0.1, 0.1)
0.010
0.010
9
(−0.35, 0.35)
(0, 0.2)
0.050
0.050
10
(−0.35, 0.35)
(0, 0.2)
0.050
0.025
11
(−0.35, 0.35)
(−0.15, 0.2)
0.050
0.050
12
(−0.35, 0.35)
(−0.15, 0.2)
0.050
0.025
13
(−0.35, 0.35)
(−0.2, 0.2)
0.050
0.050
14
(−0.35, 0.35)
(−0.2, 0.2)
0.050
0.025
15
(−0.35, 0.35)
(0, 0.15)
0.050
0.050
16
(−0.35, 0.35)
(0, 0.15)
0.050
0.025
17
(−0.35, 0.35)
(0, 0.15)
0.050
0.010
18
(−0.35, 0.35)
(−0.15, 0.15)
0.050
0.050
19
(−0.35, 0.35)
(−0.15, 0.15)
0.050
0.025
20
(−0.35, 0.35)
(−0.15, 0.15)
0.050
0.010
21
(−0.35, 0.35)
(−0.1, 0.1)
0.050
0.050
22
(−0.35, 0.35)
(−0.1, 0.1)
0.050
0.025
23
(−0.35, 0.35)
(−0.1, 0.1)
0.050
0.010Figure7 gives the reconstruction errors of pressure at 700 Hz obtained by using I-ESM with different CRP configurations when a marble floor was used as the reflecting plane.(1)
It is obvious that the reconstruction error is the smallest whenhz=−0.18m, showing that −0.18 m is a nearly perfect choice as the retreat distance. Also, the loudspeaker was placed right 0.18 m above the reflecting plane. This illustrates that the retreat distance below the reflecting plane should be set as the distance of the source image, which is consistent with the conclusion obtained in numerical simulations.(2)
In the casehz<−0.18m, i.e., hz=−0.5m and hz=−1m, it can be seen that the reconstruction results are rather stable, but all of them are a little large.(3)
In the casehz>−0.18m, i.e., hz=−0.001m, it can be seen that the error has large fluctuations versus the group of parameters. Each time the dimension x1,x2,y1,y2 changes, the error changes significantly. Also, the distribution interval dx,dy also affects the error much. Thus, if the dimension and distribution interval are set properly, e.g. taking No. 15 or No. 18 group of parameters, the error can be quite small. From Table 3, it can be found that Ω15 just covers the area between the source and the reconstructed plane and Ω18 also covers this area though with a little larger dimensions.(4)
In the casehz=−0.18m, it is obvious that there is an abnormal point, i.e., when taking No. 6 group of parameters. From Table 3, it can be seen that Ω6 does not cover both cores of the loudspeaker. Thus, the abnormal point is understandable.Figure 7
The reconstruction errors of pressure when using different CRP configurations with a marble floor as the reflecting plane.Figure8 gives the reconstruction errors of pressure at 700 Hz when a piece of sponge was used as the reflecting plane.(1)
It is obvious that the reconstruction error fluctuates slightly versus the group of parameters whenhz≤−0.18m.(2)
In the casehz>−0.18m, i.e., hz=−0.001m, it can be seen that the error fluctuates significantly versus the group of parameters. Also, the error was the smallest when taking No. 6 or No. 21 group of parameters and the largest when taking No. 13 or No. 14 group of parameters. From Table 3, it can be found that both Ω6 and Ω21 are very small, just covering the central part of the source, while Ω13 and Ω14 are the largest dimensions. In addition, it can be found that the fewer the equivalent sources, the better the reconstruction results. This is probably because that the highly absorptive sponge contributes little for the half-space sound field and few equivalent sources are needed.(3)
In the casehz=−0.18m, it can be found that the error is the smallest when taking No. 6 or No. 21 group of parameters. This is consistent with the conclusion presented above that Ω should just cover the central part of the source and the equivalent sources should be as few as possible.Figure 8
The reconstruction errors of pressure when using different CRP configurations with a piece of sponge as the reflecting plane.


## 4.3. Confirmation of the Optimal Equivalent Source Configuration

According to the discussion presented in Section4.2, Table 4 gives six sets of equivalent source configurations for further investigations to confirm the optimal configurations. The two configurations when hz=−0.18m represent the optimal ones. Also, other four configurations, respectively, represent the worst and the best in the case hz=−0.001m.Table 4
Equivalent source configurations for confirmation in experiments.
Reflecting plane
Retreat distancehz (m)
Number of parameter groupx1,x2,y1,y2,dx,dy
A marble floor
−0.001
5
A marble floor
−0.001
18
A marble floor
−0.18
6
A piece of sponge
−0.001
13
A piece of sponge
−0.001
6
A piece of sponge
−0.18
21Figures9 and 10 give the reconstruction errors of pressures versus the frequency obtained by using I-ESM with the CRP configurations shown in Table 4. It can be seen that I-ESM provides large and unstable errors when hz=−0.001m with No. 5 or No. 13 group of parameters. But, the errors obtained when hz=−0.18m are smaller and more stable, confirming that the distance of the source image is suitable to be set as the retreat distance.Figure 9
The reconstruction error of pressure versus the frequency with a marble floor as the reflecting plane.Figure 10
The reconstruction error of pressure versus the frequency with a piece of sponge as the reflecting plane.Note that the I-ESM performs best when taking No. 6 group of parameters but fails to reconstruct the half-space sound field at 100 Hz and 500 Hz when taking No. 13 group of parameters in the case that a piece of sponge was used as the reflecting plane andhz=−0.001m. This confirms that placing a small amount of equivalent sources very close to the reflecting plane and just in the central part of source is a good choice when the reflecting plane is highly absorptive.


## 5. Conclusions

This paper optimizes the configuration of equivalent sources substituting for the reflections for I-ESM. A spherical source was used to carry out numerical simulations, and a loudspeaker was used to carry out two experiments, in which a marble floor and a piece of sponge were, respectively, used as the reflecting plane. By analyzing and discussing the influence of the equivalent source configurations on the reconstruction accuracy, the appropriate CRP configuration is acquired. That is, the equivalent sources should be placed at the position of the image source, verifying that the IIS configuration is indeed an appropriate one. But, if the reflecting plane is consisted of highly absorptive material, such as sponge, it will be better to place few equivalent sources very closely to the plane and just at the central part of source projection.Moreover, both the numerical simulations and experiments indicate that the optimal equivalent source configuration can improve the performance of I-ESM.



---

*Source: 6029393-2020-03-11.xml*