# Approximate computation of acoustic reflection and shadow effects using the Kirchhoff method

Gilles Rahier, Christophe Peyret

<figure><img src="image_1.png" /></figure>

## To cite this version:

Gilles Rahier, Christophe Peyret. Approximate computation of acoustic reflection and shadow effects using the Kirchhoff method. Acta Acustica, 2024, 8 (58), 10.1051/aacus/2024041. hal-04799028

HAL Id: hal-04799028

https://hal.science/hal-04799028v1

Submitted on 22 Nov 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
The idea is to be able to evaluate the main reflection and shadow effects of a body, by a fast “first-order” method, without weighing the aerodynamic computation down. This approximate calculation can also show, depending on the results, whether a more precise but more expensive calculation is necessary before performing it. For example, a more precise calculation may be unnecessary if the reflected levels are deemed negligible. Conversely, it may be necessary if the levels found show a strong gradient in a critical region of interest. The simplified method can also speed up the optimization of the position of a reflective surface with the aim of reducing noise levels in a given direction, and save an elaborate calculation for the most promising configuration, to ensure an accurate level prediction. However, the method should not be limited to academic geometries or noise sources. It must work for bodies of any geometry (fuselage, wing, etc.) and noise sources whose spatial and spectral extent can be relatively large (jet noise, for example).

The solution proposed in this study is based on the Kirchhoff method in the time domain. It is therefore directly applicable to the flow disturbances provided by unsteady numerical simulations, without prior spectral analysis. It consists of a first calculation of free-field radiation to determine the incident acoustic pressure on a fictitious surface representing the skin of the reflecting body, then in a second radiation calculation to deduce the effects of this reflecting body on the complete acoustic field. This approach relies on two simplifying assumptions: a uniform flow around the body and short wavelengths compared to the local radius of curvature of the body surface. It can be noted that the Kirchhoff method has already been used for simplified reflection calculations [16–20], generally in the frequency domain, but the approach presented in this article is original in the simplicity and efficiency of the calculation algorithm that does not require any information on the origin of the incident acoustic field.

The article is organized as follows. In Section 2, the expression of the Kirchhoff surface integral is first reviewed and the method for calculating reflection and shadow effects is described. In Section 3, the consequences of the assumption of short wavelength compared to the local radius of curvature of the reflecting body are quantified by comparison with analytical solutions. An application example is then given for a wing in a Mach 0.5 flow, in Section 4. The calculation cost and a first assessment of the method are provided in the conclusion, together with leads for possible improvements.

## 2 Formulation

Let us define a reference frame $\mathbf{Y}$ in which the unperturbed fluid is in translation at a uniform velocity $U_0\vec{Y}_1$ and $\mathbf{S}$, the solid surface of a body fixed in this reference frame. Let $x_i$ and $y_i$ respectively be the space coordinates of an observation point $\vec{x}$ and a point $\vec{y}$ of surface $\mathbf{S}$ (Fig. 1). Let $p_{inc}$ be the pressure radiated on $\mathbf{S}$ by noise sources (that do not need to be determined for the reflection calculation).

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Schematic of acoustic waves towards the solid body and reflection towards the observation point.</figcaption></figure>

The radiation of the pressure $p(\vec{x}, t)$ reflected by $\mathbf{S}$ at the observation point $\vec{x}$ can then be calculated by a Kirchhoff integral starting from the reflected pressure $p_{ref}$ on $\mathbf{S}$ and its gradient $\nabla p_{ref}$, that can be expressed as (see Appendix A):

$$p(\vec{x}, t) = \int_{\tau} \int_{S} F_{ref}(\vec{x}, \vec{y}, \tau) G dS d\tau$$

$$\text{with } F_{ref}(\vec{x}, \vec{y}, \tau) = p_{ref} \frac{\beta^2}{d^2} (\vec{x} - \vec{y}) \cdot \vec{n} - \frac{\partial p_{ref}}{\partial n} + M_0^2 \frac{\partial p_{ref}}{\partial y_1} n_1 + \frac{1}{a_0} \left( \frac{(\vec{x} - \vec{y}) \cdot \vec{n}}{d} + M_0 n_1 \right) \frac{\partial p_{ref}}{\partial \tau} \quad (1)$$

where $G = \delta(g)/4\pi d$ the free-space Green's function in $R_3$ for the convective wave equation, $\delta(f)$ the Dirac distribution, $g = \tau - t + \sigma/a_0$, $\sigma = (d - M_0(x_1 - y_1))/\beta^2$, $d = \sqrt{C_i(x_i - y_i)^2}$, $C_i = (1 - \beta^2)\delta_{1i} + \beta^2$, $\delta_{ij}$ the Kronecker delta symbol (1 if $i = j$, zero otherwise), $\beta^2 = 1 - M_0^2$, $M_0 = U_0/a_0$, $a_0$ the speed of sound in the unperturbed medium, $\vec{n}$ the normal unit vector on the surface $\mathbf{S}$, pointing outward the solid body, $\frac{\partial}{\partial n} = \frac{\partial}{\partial y_i} n_i$, $\tau$ an emission time and $t$ a reception time.

The surface of the reflecting body is assumed to be rigid and parallel to the mean flow. The latter assumption should have limited consequences for aeronautical applications, where the mean flow approaches most aircraft surfaces with a low angle of incidence.

If the body's local radius of curvature is large enough compared to the wavelength, $\mathbf{S}$ can be locally considered as a plane surface from the acoustic point of view. The reflected pressure $p_{ref}$ on $\mathbf{S}$ and its gradient $\nabla p_{ref}$ can then be simply deduced from the incident pressure $p_{inc}$ and its gradient $\nabla p_{inc}$. The locally plane reflection hypothesis leads to the following relations:

$$\begin{gathered} P_{ref} = P_{inc} \\ \nabla p_{ref} = \nabla p_{inc} - 2 \frac{\partial p_{inc}}{\partial n} \vec{n} \end{gathered}$$

which lead to:

$$\frac{\partial p_{ref}}{\partial \tau} = \frac{\partial p_{inc}}{\partial \tau}$$
<table><caption>**Table 1.** Cases “Reflection side”, “Shadow side”, “Outgoing wave” and values of the coefficients $\kappa_1$ and $\kappa_2$ depending on the case.</caption><thead><tr><th></th><th>(a)<br>Reflection side<br>p<sub>inc</sub><br>Wall<br>Obs.</th><th>(b)<br>Shadow side<br>p<sub>inc</sub><br>Wall<br>Obs.</th><th>(c)<br>Outgoing wave (1)<br>Wall<br>p<sub>inc</sub><br>Obs.</th><th>(d)<br>Outgoing wave (2)<br>Wall<br>Obs.</th></tr></thead><tbody><tr><td>(x̄ - ȳ) ⋅ n̄</td><td>&gt;0</td><td>&lt;0</td><td>&gt;0</td><td>&lt;0</td></tr><tr><td>(∂p<sub>inc</sub>/∂n)(∂p<sub>inc</sub>/∂τ)</td><td>&gt;0</td><td>&gt;0</td><td>&lt;0</td><td>&lt;0</td></tr><tr><td>Values of κ<sub>1</sub> and κ<sub>2</sub></td><td>κ<sub>1</sub> = 1<br>κ<sub>2</sub> = 1</td><td>κ<sub>1</sub> = 1<br>κ<sub>2</sub> = 0</td><td>κ<sub>1</sub> = 0</td><td>κ<sub>1</sub> = 0</td></tr></tbody></table>

$$ \frac{\partial p_{ref}}{\partial n} = \frac{\partial p_{inc}}{\partial n} - 2 \frac{\partial p_{inc}}{\partial n} $$

$$ \frac{\partial p_{ref}}{\partial y_1} n_1 = \frac{\partial p_{inc}}{\partial y_1} n_1 - 2 \frac{\partial p_{inc}}{\partial n} n_1^2 $$

Expressing $p_{ref}$, $\frac{\partial p_{ref}}{\partial n}$ and $\frac{\partial p_{ref}}{\partial y_1} n_1$ in formula (1) by their expressions above and rearranging leads to:

$$ F_{ref}(\vec{x}, \vec{y}, \tau) = p_{inc} \frac{\beta^2}{d^2} (\vec{x} - \vec{y}) \cdot \vec{n} - \frac{\partial p_{inc}}{\partial n} + M_0^2 \frac{\partial p_{inc}}{\partial y_1} n_1 + 2(1 - M_0^2 n_1^2) \frac{\partial p_{inc}}{\partial n} + \frac{1}{a_0} \left( \frac{(\vec{x} - \vec{y}) \cdot \vec{n}}{d} + M_0 n_1 \right) \frac{\partial p_{inc}}{\partial \tau} \quad (2) $$

The source term $F_{ref}(\vec{x}, \vec{y}, \tau)$ of the reflected pressure can thus be expressed entirely as a function of the incident pressure and its gradient, which are input data assumed to be known.

In practice, in order to best model the reflection and shadow effects of the solid body, four cases must be distinguished for each surface element $\Delta S$ of its surface $S$. These four cases are illustrated in Table 1:

Case (a): the incident wave is incoming with respect to the surface element, and the observation point is on the same side as the local normal;

Case (b): same configuration for the incident wave, but the observation point is on the opposite side to the local normal;

Case (c) and case (d): the incident wave is outgoing with respect to the surface element.

Formula (2) only applies to an observer on the reflection side (case (a)). For an observer on the other side (shadow side, case (b)), the shadow effect of the body can be modeled in the calculation by a term creating a radiation opposite to the incident wave, so that the sum with the direct (free-field) pressure is zero. By taking into account the orientation of the local normal $\vec{n}$, the corresponding source term is simply formula (1) applied to $p^{inc}$:

$$ F_{shadow}(\vec{x}, \vec{y}, \tau) = p_{inc} \frac{\beta^2}{d^2} (\vec{x} - \vec{y}) \cdot \vec{n} - \frac{\partial p_{inc}}{\partial n} + M_0^2 \frac{\partial p_{inc}}{\partial y_1} n_1 + \frac{1}{a_0} \left( \frac{(\vec{x} - \vec{y}) \cdot \vec{n}}{d} + M_0 n_1 \right) \frac{\partial p_{inc}}{\partial \tau} \quad (3) $$

Finally, there is no reflection or shadow effect to calculate in the case of the incident wave coming out of the body surface (outgoing wave, cases (c) and (d)). The source term corresponding to these cases is therefore zero.

For each point on the surface $S$ and for each observation point, the direction of propagation of the incident wave and the relative position of the observation point with respect to the reflecting element can be determined by the sign of $(\partial p/\partial n)(\partial p/\partial \tau)$ and of $(\vec{x} - \vec{y}) \cdot \vec{n}$ as shown in Table 1.

All cases can therefore be easily managed by two coefficients $\kappa_1$ and $\kappa_2$ in the following expression of the source term:

$$ F(\vec{x}, \vec{y}, \tau) = \kappa_1 \begin{pmatrix} \frac{\beta^2}{d^2} (\vec{x} - \vec{y}) \cdot \vec{n} \\ -\frac{\partial p_{inc}}{\partial n} + M_0^2 \frac{\partial p_{inc}}{\partial y_1} n_1 \\ +\frac{1}{a_0} \left( \frac{(\vec{x} - \vec{y}) \cdot \vec{n}}{d} + M_0 n_1 \right) \frac{\partial p_{inc}}{\partial \tau} \\ +2\kappa_2 (1 - M_0^2 n_1^2) \frac{\partial p_{inc}}{\partial n} \end{pmatrix} \quad (4) $$

In the general case, the direction of the incident wave is not constant, and the coefficients $\kappa_1$ and $\kappa_2$ must obviously be determined for each time step of the reflection calculation. The discontinuities created by this coefficient switching during a time-domain calculation have no major impact on the results, mainly because the influence of these discontinuities is linked to the size of the elementary surface $\Delta S$ of the reflecting surface $S$, and the contribution decreases on finer meshes. The space-time shifting of these discontinuities on the surface $S$, during the simulation, also tends to spread their effects over time at the observation point.
In summary, a reflecting body can be taken into account in an approximate way in a radiation calculation in free field starting from the knowledge of the incident pressure, its time derivative and its gradient on the body surface. The incident pressure on the body surface is provided by the radiation calculation in free field and its time derivative can be easily deduced by finite differences. So, only the pressure gradient remains to be determined. In the general case, this can also be calculated by finite differences of the pressure on a double-layered mesh on either side of the body skin. This approach is used in the calculations that follow. It has the disadvantage of requiring to define the distance between the two layers so that the calculation of the pressure gradient is precise enough for the frequency band of interest (thickness small enough compared to the wavelengths while remaining compatible with the precision of the incident radiation calculation). A more accurate and safer way to proceed is possible. It consists in calculating the pressure gradient directly by the same integral formulation as that used for the pressure. Expressions of this gradient are provided in Appendix B, for fixed integration surfaces, for both Kirchhoff and FW-H integral formulations.

### *Remarks:*

* For a correct summation of direct and reflected pressures (Fig. 2), it is of course necessary for both pressures to be calculated with the same time reference point, so that the causal link between both pressures at the observer point is maintained. The calculation of the reflection of each elementary surface $\Delta S$ of $S$ must therefore be carried out from the reception times of the direct radiation on this elementary surface, so that the final reception times at the observer point include all propagation times.

* In the case of non-periodic acoustic fields (for which the signals cannot be folded back on a single period), the different propagation times must be taken into account to determine the part of the total signal that is valid, i.e. the part that has received both the direct signal and the entire contribution of the reflecting surface, as shown in Figure 2.

* The reflection calculation can include several reflective surfaces, but, for each surface, it only takes into account the direct radiation, not the radiation of other reflective surfaces. The method is therefore not intended for acoustic propagation in ducts or for the prediction of stationary waves that can form a standing wave pattern between two reflective surfaces.

## 3 Comparison with analytical solution

In this section, reflection calculation using the Kirchhoff surface integral is applied to the prediction of the scattering of a point source by a rigid sphere, and the results are compared with the analytical solution. The configuration is as follows: a sphere of radius $R = 0.17$ m in a fluid at rest ($M_0 = 0$), a monopole point source at frequency $f$ placed at $x_S = 0.34$ m from the center of the sphere in its horizontal median plane $z = 0$ (Fig. 3), 180 observation points evenly distributed on a circle of radius $r = 0.6$ m around the

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Illustration of the effect of propagation times on the duration of the valid total signal (case of non-periodic signals).</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Calculation of the scattering of a point source by a rigid sphere.</figcaption></figure>

sphere, in the horizontal plane. The calculations are performed for a speed of sound $a_0 = 340$ m/s and for four frequencies: $f = 1$ kHz ($\lambda/R = 2$), $f = 2.5$ kHz ($\lambda/R = 0.8$), $f = 5$ kHz ($\lambda/R = 0.4$) and $f = 7.5$ kHz ($\lambda/R = 0.27$).

The analytical solution of the scattered pressure in complex notations, at the observation points $(r, \theta)$ is [21]:
$$P_T(r, \theta) = \frac{ik}{4\pi} \sum_{m=0}^{\infty} (2m+1)P_m(\cos \theta) \\ \times \left[ j_m(kr)h_m(kx_S) - \frac{j'_m(kR)}{h'_m(kR)}h_m(kx_S)h_m(kr) \right] \quad r \le x_S \\ P_T(r, \theta) = \frac{ik}{4\pi} \sum_{m=0}^{\infty} (2m+1)P_m(\cos \theta) \\ \times \left[ j_m(kx_S)h_m(kr) - \frac{j'_m(kR)}{h'_m(kR)}h_m(kx_S)h_m(kr) \right] \quad r > x_S$$

with $i^2 = -1$, $k = 2\pi/\lambda$ and where $P_m$ is the Legendre polynomial of degree $m$, $j_m$ the first kind spherical Bessel function of order $m$, $j'_m$ its derivative, $h_m = j_m + iy_m$ the Hankel function of the first kind and order $m$, $h'_m$ its derivative and $y_m$ the second kind spherical Bessel function of order $m$. Below, these series are calculated for a number of terms that ensure an error of less than $10^{-10}$.

For the Kirchhof method, the surface of the sphere is represented by a uniform unstructured mesh composed of triangles. This mesh is duplicated on either side of the sphere surface (radial distance $\pm 0.00034$ m) for the calculation of the pressure gradient. The complex pressure radiated by the monopole on a node $\vec{x}$ of the double-layer mesh and at an instant $t$, is [22, 23]:

$$\hat{p}'(\vec{x}, t) = \rho_0 \left( \frac{\partial \phi}{\partial t} + U_0 \frac{\partial \phi}{\partial x_1} \right) \quad \text{where} \quad \phi(\vec{x}, \vec{y}, t) = \frac{A}{4\pi d} \exp^{i(2\pi f)(t-\sigma/a_0)} \text{ is the complex potential of the monopole,}$$

with $A$ its strength and $\vec{y} = (0.34, 0, 0)$. The input to Kirchhoff's calculation is the real part of this pressure:

$$p'(\vec{x}, t) = \frac{A\rho_0}{4\pi d} \left[ 2\pi f \left( 1 - \frac{U_0}{a_0} \frac{\partial \sigma}{\partial x_1} \right) \sin \left( 2\pi f \left( t - \frac{\sigma}{a_0} \right) \right) + \frac{U_0}{d} \frac{\partial d}{\partial x_1} \cos(2\pi f(t - \sigma/a_0)) \right].$$

The computation time step is equal to 1/100 of the period $T$ of the monopole.

As a first step, calculations are carried out for the highest frequency ($f = 7.5$ kHz), for three meshes of different densities, to check the stability of the results. The number of cells, the average distance $L$ between two nodes of a cell, and the $\lambda/L$ ratio are given in Table 2, for the three meshes and for the four frequencies. These values are given for guidance only. Since the angle of incidence of the acoustic waves on the sphere varies from 0° to 90°, choosing a uniform mesh based on the $\lambda/L$ criterion alone would be unnecessarily penalizing. The criterion considered is simply the minimum mesh density at which the results no longer change significantly.

The results for the three meshes are compared in Figure 4. The quantity plotted is the attenuation factor $\text{Gamma Total (GT)}$ defined as $\gamma_T(f) = \langle p_T(f) \rangle / \langle p_I(f) \rangle$, with $\langle \cdot \rangle$ indicating the ensemble average operator, and

<table><caption>Table 2. Number of cells, average distance $L$ and ratio $\lambda/L$.</caption><thead><tr><th>Mesh</th><th>Coarse</th><th>Medium</th><th>Fine</th></tr></thead><tbody><tr><td>Number of cells</td><td>3214</td><td>11,968</td><td>19,594</td></tr><tr><td>L (m)</td><td>0.0162</td><td>0.0084</td><td>0.0065</td></tr><tr><td>λ/L (f = 1 kHz)</td><td>21.0</td><td>40.5</td><td>52.3</td></tr><tr><td>λ/L (f = 2.5 kHz)</td><td>8.4</td><td>16.2</td><td>20.9</td></tr><tr><td>λ/L (f = 5 kHz)</td><td>4.2</td><td>8.1</td><td>10.5</td></tr><tr><td>λ/L (f = 7.5 kHz)</td><td>2.8</td><td>5.4</td><td>7.0</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>Figure 4. Influence of the mesh density on the Kirchhoff results ($f = 7.5$ kHz).</figcaption></figure>

$p_I$ and $p_T$, the incident and the total pressure fluctuations, respectively.

There are no significant differences between the directivities obtained with the medium and fine meshes. On the other hand, the coarse mesh leads to significantly different levels for observation angles between 45° and 315°. However, the differences remain moderate, indicating a certain robustness to the method. Given these results for the highest frequency, comparisons with analytical solutions are made with the medium mesh.

For all four frequencies, the Kirchhoff method provides realistic acoustic directivities with the correct number of lobes (Fig. 5). For the lowest frequency ($f = 1$ kHz), the wavelength is not small compared to the surface's radius of curvature ($\lambda/R = 2$) and the locally plane reflection assumption distorts the levels. In particular, the shadow effect in the half-space $x < 0$ is overestimated due to the incomplete modeling of refraction. The levels are correct only in a small angular sector in front of the source ($-10° < \theta < 10°$). The accuracy of calculations by the Kirchhoff method logically improves with increasing frequency, i.e. with decreasing the ratio $\lambda/R$ (Figs. 5b–5d). For the highest frequency ($f = 7.5$ kHz, $\lambda/R = 0.27$), the levels match the analytical solution throughout all the half-space $x > 0$ on the source side.
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Directivity of the scattering of a point source by a sphere. (a) $f = 1$ kHz, $\lambda/R = 2$; (b) $f = 2.5$ kHz, $\lambda/R = 0.8$; (c) $f = 5$ kHz, $\lambda/R = 0.4$ and (d) $f = 7.5$ kHz, $\lambda/R = 0.27$.</figcaption></figure>

Each of these calculations carried out for 11,968 reflective cells, 100 time steps, and 180 observation points, takes 12.3 s on a processor running at 3 GHz.

## 4 Comparisons with Euler calculations

In this section, reflection calculations by the Kirchhoff surface integral are carried out for an airfoil in a flow, as an application example. The results are compared with the numerical solution of the linearized Euler equations, in the presence of a mean flow.

The reference solution is a 2D Euler calculation in the frequency domain, based on a high-order discontinuous Galerkin (DG) method [24, 25]. The reflecting surface is a NACA 0012 profile of chord $C = 1$ m, without incidence in a flow at Mach number 0.5. The computational domain is a disk of radius 8 m around the profile. A noise source near the leading edge of the profile is modeled by a sinusoidal pressure disturbance at frequency $f = 2$ kHz, with a Gaussian distribution of width 0.02 m at the midpoint of the Gaussian. The center of the source is 0.1 m upstream

of the leading edge of the profile, 0.2 m below. Two unstructured meshes made of triangles are designed, one without the profile, the other with it, to highlight the effect of the profile on the pressure field radiated by the acoustic source (Fig. 6). Particular care is taken to ensure that these two meshes are as similar as possible in the common part of the fluid domain. The number of cells is respectively 10,324 for the mesh without profile and 8580 for the mesh with profile. The mean flow is uniform in the case without profile. It is determined by solving Euler's equations in the case with the profile. DG acoustic calculations are performed at order 8 in both cases. The harmonic solution of the density, velocity and pressure fields is projected onto 90 points distributed in a circle of radius 0.1 m around the source for FW-H free-field radiation calculations (the source term $F$ of this formulation is recalled in Appendix B).

This 2D configuration is best modeled in the 3D FW-H and Kirchhoff calculations presented below, with a wing and a cylindrical “source” surface obtained by duplicating the profile and the “source” circle in the third direction. This duplication is carried out over 50 m on either side of the DG calculation plane, with a step of 0.1 m, i.e. a total of
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Zoom in the unstructured meshes used for the 2D DG calculations in free field (left) and with the NACA profile (right). Circle of input data for FW-H calculations (in red).</figcaption></figure>

1001 vertical planes (Fig. 7). For FW-H calculations of the source free-field radiation, the “source” surface is therefore described by a structured grid of 1001 × 90 nodes. On each of the 90 lines of this grid, the harmonic solution (amplitude and phase) provided by the DG calculation on the “source” circle, is converted into the time domain with a time step equal to 1/50 of the period *T* of the noise source. For reflection calculations using the Kirchhoff method, the wing is described by a structured double-layer grid, with 1001 points in span and 99 points in chord, many of which concentrated at the leading edge for a precise description of its geometry. The thickness *e* of the double layer on either side of the wing skin is 0.001 m (λ/e = 170) for an accurate calculation of the pressure gradients at its leading edge. The calculations are carried out for a speed of sound $a_0 = 340$ m/s.

The calculation procedure is as follows. First, FW-H calculations of free-field radiation are performed for 109 × 90 observation points describing a 4 m long and 3 m high rectangle around the wing, in its median vertical plane (Fig. 7) and for the nodes of the wing’s double-layer mesh. Then, using the results on the wing, the Kirchhoff calculations of the pressure reflected by the wing is carried out for the observation points. Finally, the results of both calculations in the observation plane are summed. It should be noted that this procedure, combining FW-H calculations (with an integration surface surrounding the noise sources) and Kirchhoff calculations (to take into account the reflections of solid bodies outside the FW-H control surface), can be applied to most numerical aerodynamic simulations.

The results of the FW-H calculations are compared with the reference solution in Figures 8 and 9. The discrepancies between the two calculations are not discernible for instantaneous fields (Fig. 8). The agreement of the integrated levels (Fig. 9) confirms both the validity of the FW-H calculations starting from the DG fields on the “source” surface, and the sufficient span of the wing to find the 2D solution in the median plane. The slight loss of accuracy of the DG calculations in the bottom-left corner is due to the slightly less dense mesh in this area (Fig. 6). The small level oscillations in the FW-H results are due both to the rather coarse discretization of the “source” cylinder in span ($Δz = 0.1$ m) and to its finite length. Overall, these FW-H

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Configuration of the reflection calculation by Kirchhoff method: control surface for free-field acoustic radiation using FW-H formulation (in red), reflecting surface (in green), and observation plane (in blue).</figcaption></figure>

calculations appear to provide sufficiently accurate input data to the Kirchhoff reflection calculation to enable valid comparisons with the reference solution.

The results of the Kirchhoff reflection calculation are shown in Figure 10. It can be seen that the profile radiates reflected waves on the source side and waves in phase opposition with the incident pressure on the opposite side. Naturally, the integrated levels do not allow us to differentiate the constructive contributions of this radiation from the destructive ones.

The total pressure, resulting from the sum of the direct and reflected pressures, is plotted in Figures 11 and 12. For the frequency studied, they show that the profile strongly modifies the acoustic field of the noise source, with fairly marked areas of amplification or silence. The integrated levels predicted by the simplified method are close to those obtained by Euler’s calculations. The two main deviations from the reference solution are: less pronounced directivity lobes upstream of the profile, and a less destructive interference at 45° below the profile. The shadow zone above the profile is close to the reference solution, despite the assumption of a large local radius of curvature compared to the wavelength. This is a very satisfactory result, given the small radius of curvature of the profile trailing edge. The deviations in the bottom-left corner, due to slight
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Instantaneous pressure radiated by the noise source in free-field (without the profile, direct pressure). DG calculation (left) and FW-H calculation (right).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 9.** Same comparison as above for the integrated pressure levels.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 10.** Reflected pressure. Kirchhoff calculation. Instantaneous pressure (left) and integrated levels (right).</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 11.** Instantaneous pressure with the profile (total pressure). DG calculation (left) and FW-H calculation (right).</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>**Figure 12.** Same comparison as above for the integrated pressure levels.</figcaption></figure>

numerical dissipation in the DG calculation, show, conversely, the interest of Kirchhoff calculations for far-field acoustic predictions, especially when the mesh of aerodynamic simulation cannot be designed for acoustic propagation over long distances for reasons of computational cost.

The reflection calculation carried out with the Kirchhoff method for 99,000 reflective cells, 50 time steps, and 9810 observation points, take about 40 min on a processor running at 3 GHz.

## 5 Conclusion

Assuming locally plane reflections and a uniform mean flow parallel to the reflecting surface, a calculation of the reflection and shadow effects of a solid body (wing, fuselage, building, ground relief) in an acoustic field was obtained rather simply with a Kirchhoff surface integral. The calculation time of the reflected pressure field is about $5 \cdot 10^{-8}$ s/reflective cell/time step/obs on a processor running at 3 GHz. For example, for a reflective surface described by 20,000 cells and signals described by 5000 instants, the calculation time is about 5 s per observation point. In practice, a few hundred observation points are enough to characterize acoustic radiation in the far field. Calculating the reflected pressure field is therefore inexpensive in itself. In fact, what can be most costly is to determine the incident pressures on the body surface, but this input data is necessary whatever method is used.

Test cases show that the proposed model can provide realistic and useful results for a first analysis of acoustic scattering, for example as part of an optimization process. The results are all the more accurate as the assumption of short wavelength compared to the local radius of curvature of the reflective surface is satisfied. In the present study, the incident pressure gradient on the reflective surface was calculated by finite differences on two layers on either side of the surface, spaced by a prescribed value. When the incident pressure is calculated using the Kirchhoff or Ffowcs Williams and Hawkings formulation, it is preferable to determine this pressure gradient directly using the same integral formulation, with the corresponding source terms. This prevents the inaccuracies of the calculation by finite differences, especially for broadband noise sources, where

the distance between the two layers cannot logically be optimal for all frequencies.

Apart from its computational cost, another strength of the method is that it can be applied to unsteady aerodynamic numerical simulations without any sound source model. It works directly starting the incident pressure and its gradient on the reflective surfaces and requires no information on the origin of the incident acoustic field. Furthermore, time-domain calculation offers a significant advantage for periodic signals rich in harmonics or broadband signals. This method is therefore particularly attractive for BVI or HSI noise sources of rotors or propellers, and for noise sources of wakes or jets (recent example of application to a jet in [26]). For all these applications, BEM methods in the frequency domain remain considerably more expensive (when applicable).

The drawbacks of the present approach are mainly its lack of accuracy for high wavelengths (with respect to the body's radius of curvature), and its inability to handle sharp angles correctly. A possible improvement could be to use this high-frequency approximation in conjunction with a BEM method for the low frequencies. Once the spectral analysis of the results of the Kirchhoff calculation has been carried out, this would consist of replacing the imprecise levels at low frequencies with those obtained with the BEM calculations. It should also be possible to enrich the method to deal with sharp angles, although, from our point of view, at the cost of a level of complexity that would probably preclude interest of this approach for industrial applications.

## Conflicts of interest

The authors declare no conflicts of interest.

## Data availability statement

The data are available from the corresponding author on request.

## References

1. S.K. Lele, J.W. Nichols: A second golden age of aeroacoustics? Philosophical Transactions of the Royal Society A372 (2014) 20130321.
2. S. Schoder, M. Kaltenbacher: Hybrid aeroacoustic computations: state of art and new achievements. Journal of Theoretical and Computational Acoustics 27, 4 (2019) 1950020.

3. M.E. Golstein: Aeroacoustics. McGraw Hill International Book Company, New York City, New York, USA, 1976, pp. 59–61.

4. J.E. Ffowcs Williams, D.L. Hawkins: Sound generation by turbulence and surfaces in arbitrary motion. Philosophical Transactions of the Royal Society of London 264, A1151 (1969) 321–342.

5. J. Utzmann, C.-D. Munz, M. Dumbser, E. Sonnendrücker, S. Salmon, S. Jund, E. Frénod: Fluid-acoustic coupling and wave propagation. Numerical Simulation of Turbulent Flows and Noise Generation 104 (2009) 47–74.

6. G. Djambazov, C.-H. Lai, K. Pericleous: On the coupling of Navier-Stokes and linearised Euler equations for aeroacoustic simulation. Computing and Visualization in Science 3 (2000) 9–12.

7. R. Harris, E. Collins, E. Luke, A. Sescu: Coupled overset unstructured discontinuous Galerkin method for launch environment acoustics prediction. AIAA Journal 54, 6 (2016) 1932–1952.

8. A. Langenais, F. Vuillot, Ch. Peyret, G. Chaineray, Ch. Bailly: Assessment of a two-way coupling methodology between a flow and a high-order nonlinear acoustic unstructured solvers. Flow, Turbulence and Combustion 101 (2018) 681–703.

9. A. Tosh, P. Liever, F. Owens, Y. Liu: A high-fidelity CFD/BEM methodology for launch pad acoustic environment prediction. In: 18th AIAA/CEAS Aeroacoustics Conference (33rd AIAA Aeroacoustics Conference), 2012, p. 2107.

10. L. Dürrwächter, M. Kebler, E. Krämer: Numerical assessment of open-rotor noise shielding with a coupled approach. AIAA Journal 57, 5 (2019) 1930–1940.

11. D. Casalino, M. Barbarino, A. Visingardi: Simulation of helicopter community noise in complex urban geometry. AIAA Journal 49, 8 (2011) 1614–1624.

12. G.V. Groves: Geometrical theory of sound in the atmosphere. Journal of Atmospheric and Terrestrial Physics 7 (1955) 113–127.

13. R.J. Thompson: Ray theory for an inhomogeneous moving medium. JASA 51 (1965) 1675–1682.

14. S.M. Candel: Numerical solution of conservation equations arising in linear wave theory: application to aeroacoustics. Journal of Fluid Mechanics 83 (1977) 465–493.

15. Y.P. Guo, R.H. Thomas: Geometric acoustics for aircraft noise scattering. In: 28th AIAA/CEAS Aeroacoustics Conference, AIAA 2022-3077, June 2022.

16. T.V. Cox, Y.W. Lam: Evaluation of methods for predicting the scattering from simple rigid panels. Applied Acoustics 40, 2 (1993) 123–140.

17. K.G. Foote, D.T.I. Francis: Comparing Kirchhoff-approximation and boundary-element models for computing gadoid target strengths. Journal of the Acoustical Society of America 111, 4 (2002) 1644–1654.

18. F.-E. Aballéa, J. Defrance: Single and multiple reflections in plane obstacle using the parabolic equation method with a complementary Kirchhoff approximation. Acta Acustica 93 (2007) 22–30.

19. N.J. Pignier, C.J. O'Reilly, S. Boij: A Kirchhoff approximation-based numerical method to compute multiple acoustic scattering of a moving source. Applied Acoustics 96 (2015) 108–117.

20. Y.P. Guo, D.S. Pope, C.L. Burley, R.H. Thomas: Aircraft system noise shielding prediction with a Kirchhoff integral method. In: 23rd AIAA/CEAS Aeroacoustics Conference, AIAA 2017-3196, June 2017.

21. J.J. Bowman, T.B.A. Senior, P.L.E. Uslenghi: Electromagnetic and acoustics scattering by simple shapes. North Holland Publishing Company, Amsterdam, 1969.

22. D. Lockard: A comparison of Ffowcs Williams-Hawkings solvers for airframe noise applications. In: AIAA paper 2002-2580, 8th AIAA/CEAS Aeroacoustics Conference & Exhibit, Breckenridge, Colorado, June 2002.

23. A.P. Dowling, J.E. Ffowcs Williams: Sound and sources of sound, chap. 9. Horwood Publishing, Westergate, 1983, pp. 207–208.

24. Ph. Delorme, P. Mazet, Ch. Peyret, Y. Ventribout: Computational aeroacoustics applications based on a discontinuous Galerkin method. Comptes Rendus Mécanique 333 (2005) 676–682.

25. Ch. Peyret: A full high order method for computational aeroacoustics. In: 23rd AIAA/CEAS Aeroacoustics Conference, Denver, Colorado, June 2017.

26. M. Huet, F. Gand, G. Rahier: Simulation of isolated and installed jet noise at Mach=0.9: influence of numerical mesh and physical insights. Flow, Turbulence and Combustion, published on line: 31 August 2023.

27. J. Prieur, G. Rahier: Aeroacoustic integral methods, formulation and efficient numerical implementation. Aerospace Science and Technology 5, 7 (2001) 457–468.

## Appendix A: Expression of the Kirchhoff surface integral for a fixed integration surface and fixed observation points

The starting point of the formulation is the generalized Green's formula derived by Goldstein ([3], pp. 59–61):

$$p(\vec{x}, t) = \int_{\tau} \int_{S} \left[ \left( p \frac{\partial G}{\partial n} - G \frac{\partial p}{\partial n} \right) + \frac{M_0 n_1}{a_0} \left( G \frac{dp}{d\tau} - p \frac{dG}{d\tau} \right) \right] dS d\tau \\ \text{with } \frac{d}{d\tau} = \frac{\partial}{\partial \tau} + U_0 \frac{\partial}{\partial y_1} \text{ and the notations given in Section 2.}$$

Developing the derivatives $d/d\tau$ gives,

$$p(\vec{x}, t) = \int_{\tau} \int_{S} \left[ p \left( n_i \frac{\partial G}{\partial y_i} - \frac{M_0 n_1}{a_0} \frac{\partial G}{\partial \tau} - M_0^2 n_1 \frac{\partial G}{\partial y_1} \right) - G \left( \frac{\partial p}{\partial n} - \frac{M_0 n_1}{a_0} \frac{\partial p}{\partial \tau} - M_0^2 n_1 \frac{\partial p}{\partial y_1} \right) \right] dS d\tau$$

Knowing that, for a fixed surface $S$ and a fixed observation point $\vec{x}$, $\partial G/\partial \tau = 1$ and $\partial d/\partial \tau = 0$,

$$\int_{\tau} p \frac{\partial G}{\partial \tau} d\tau = \int_{g} p \frac{\partial G}{\partial g} dg = \int_{g} \frac{p}{4\pi d} \frac{\partial(\delta(g))}{\partial g} dg \\ = - \int_{g} \frac{\partial p}{\partial g} \frac{\delta(g)}{4\pi d} dg = \int_{\tau} \frac{\partial p}{\partial \tau} G d\tau$$

Developing $\frac{\partial G}{\partial y_i} = \frac{1}{4\pi d} \left( \frac{\partial(\delta(g))}{\partial g} \frac{\partial g}{\partial y_i} - \frac{\delta(g)}{d} \frac{\partial d}{\partial y_i} \right)$ gives:

$$\int_{\tau} p n_i \frac{\partial G}{\partial y_i} d\tau = \int_{g} \frac{p}{4\pi d} \frac{\partial G}{\partial y_i} n_i \frac{\partial(\delta(g))}{\partial g} dg - \int_{\tau} \frac{p}{d} \frac{\partial d}{\partial y_i} n_i \frac{\delta(g)}{4\pi d} d\tau$$
$$= - \int_{\tau} \frac{\partial}{\partial \tau} \left( p \frac{\partial d}{\partial y_i} n_i \right) \frac{\delta(g)}{4\pi d} d\tau - \int_{\tau} \frac{p}{d} \frac{\partial d}{\partial y_i} n_i \frac{\delta(g)}{4\pi d} d\tau$$

$$= - \int_{\tau} \frac{\partial p}{\partial \tau} \frac{\partial g}{\partial y_i} n_i G d\tau - \int_{\tau} \frac{p}{d} \frac{\partial d}{\partial y_i} n_i G d\tau$$

$$p(\vec{x}, t) \text{ can therefore write: } p(\vec{x}, t) = \int_{\tau} \int_{S} F(\vec{x}, \vec{y}, \tau) G dS d\tau$$

$$\text{with } F(\vec{x}, \vec{y}, \tau) = - \frac{\partial g}{\partial y_i} n_i \frac{\partial p}{\partial \tau} - \frac{\partial d}{\partial y_i} n_i \frac{p}{d} + \frac{M_0 n_1}{a_0} \frac{\partial p}{\partial \tau}$$

$$+ M_0^2 \frac{\partial g}{\partial y_1} n_1 \frac{\partial p}{\partial \tau} + M_0^2 \frac{\partial d}{\partial y_1} n_1 \frac{p}{d}$$

$$- \frac{\partial p}{\partial n} + \frac{M_0 n_1}{a_0} \frac{\partial p}{\partial \tau} + M_0^2 n_1 \frac{\partial p}{\partial y_1}$$

$$\text{Developing } \frac{\partial d}{\partial y_i} = \frac{-C_i(x_i - y_i)}{d} \quad \text{and} \quad \frac{\partial g}{\partial y_i} = \frac{1}{a_0 \beta^2} \left( \frac{-C_i(x_i - y_i)}{d} + M_0 \delta_{1i} \right) \text{ and rearranging terms eventually leads to:}$$

$$F(\vec{x}, \vec{y}, \tau) = p \frac{\beta^2}{d^2} (x_i - y_i) n_i - \frac{\partial p}{\partial y_i} n_i + M_0^2 \frac{\partial p}{\partial y_1} n_1$$

$$+ \frac{1}{a_0} \left( \frac{(x_i - y_i) n_i}{d} + M_0 n_1 \right) \frac{\partial p}{\partial \tau}$$

## Appendix B: Expression of pressure gradient

The integration surface and the observation points are assumed to be fixed in the calculation reference frame. The flow is uniform at velocity $U_0 = a_0 M_0$ in the direction of the x axis.

### General expression (see notations in Sect. 2)

Starting from the expression of the pressure for the Kirchhoff or FW-H formulation in the form [27]:

$$p(\vec{x}, t) \int_{\tau} \int_{S} F(\vec{x}, \vec{y}, \tau) G dS d\tau \quad \text{with} \quad G = \frac{\delta(g)}{4\pi d}$$

the gradient can write:

$$\frac{\partial p(\vec{x}, t)}{\partial x_j} = \int_{\tau} \int_{S} \frac{\partial F}{\partial x_j} G dS d\tau + \int_{\tau} \int_{S} F \frac{\partial G}{\partial x_j} dS d\tau$$

$$\text{and developing } \frac{\partial G}{\partial x_j} \text{ gives: } \frac{\partial G}{\partial x_j} = \frac{1}{4\pi d} \left( \frac{\partial(\delta(g))}{\partial g} \frac{\partial}{\partial x_j} - \frac{\delta(g)}{d} \frac{\partial}{\partial x_j} \right)$$

For fixed integration surfaces and observation points $\partial g / \partial \tau = 1$ and thus:

$$\int_{\tau} \frac{F}{4\pi d} \frac{\partial g}{\partial x_j} \frac{\partial(\delta(g))}{\partial g} d\tau = \int_{g} \frac{F}{4\pi d} \frac{\partial g}{\partial x_j} \frac{\partial(\delta(g))}{\partial g} dg$$

$$= - \int_{g} \frac{\partial}{\partial g} \left( \frac{\partial g}{\partial x_j} \frac{F}{4\pi d} \right) \delta(g) dg = - \int_{\tau} \frac{\partial g}{\partial x_j} \frac{\partial F}{\partial \tau} \frac{\delta(g)}{4\pi d} d\tau$$

Therefore, the pressure gradient writes:

$$\frac{\partial p(\vec{x}, t)}{\partial x_j} = \int_{\tau} \int_{S} \left( \frac{\partial F}{\partial x_j} - \frac{\partial g}{\partial x_j} \frac{\partial F}{\partial \tau} - \frac{\partial d}{\partial x_j} \frac{F}{d} \right) G dS d\tau$$

Writing $\hat{r} = \frac{x_i - y_i}{d}$:

$$\frac{\partial d}{\partial x_j} = C_j \hat{r}_j \quad \text{and} \quad \frac{\partial g}{\partial x_j} = \left( \frac{C_j \hat{r}_j - M_0 \delta_{1j}}{a_0 \beta^2} \right).$$

### Application to Kirchhoff formulation

The source term $F$ is: $F = M_0^2 n_1 \frac{\partial p}{\partial y_1} - \frac{\partial p}{\partial n} + \frac{1}{a_0} \left( M_0 n_1 + \frac{(x_i - y_i) n_i}{d} \right) \frac{\partial p}{\partial \tau} + \beta^2 \frac{(x_i - y_i) n_i}{d^2} p$

As the pressure $p$ does not depend on $x_j$:

$$\frac{\partial F}{\partial x_j} = \frac{1}{a_0} \left( \frac{n_j}{d} - \frac{(x_i - y_i) n_i}{d^2} \frac{\partial d}{\partial x_j} \right) \frac{\partial p}{\partial \tau} + \beta^2 \left( \frac{n_j}{d^2} - \frac{2(x_i - y_i) n_i}{d^3} \frac{\partial d}{\partial x_j} \right) p$$

Writing $\hat{r}_n = \hat{r}_i n_i$:

$$\frac{\partial F}{\partial x_j} = \frac{1}{a_0 d} (n_j - \hat{r}_n C_j \hat{r}_j) \frac{\partial p}{\partial \tau} + \frac{\beta^2}{d^2} (n_j - 2 \hat{r}_n C_j \hat{r}_j) p$$

$$\frac{\partial g}{\partial x_j} \frac{\partial F}{\partial \tau} = \frac{(C_j \hat{r}_j - M_0 \delta_{1j})}{a_0 \beta^2}$$

$$\times \left[ M_0^2 n_1 \frac{\partial^2 p}{\partial \tau \partial y_1} - \frac{\partial^2 p}{\partial \tau \partial n} + \frac{1}{a_0} (M_0 n_1 + \hat{r}_n) \frac{\partial^2 p}{\partial \tau^2} + \beta^2 \frac{\hat{r}_n}{d} \frac{\partial p}{\partial \tau} \right]$$

$$\frac{\partial d}{\partial x_j} \frac{F}{d} = \frac{C_j \hat{r}_j}{d} \left[ M_0^2 n_1 \frac{\partial p}{\partial y_1} - \frac{\partial p}{\partial n} + \frac{1}{a_0} (M_0 n_1 + \hat{r}_n) \frac{\partial p}{\partial \tau} + \beta^2 \frac{\hat{r}_n}{d} p \right]$$

Grouping terms finally leads to the following expression of the pressure gradient:

$$\frac{\partial p(\vec{x}, t)}{\partial x_j} = \int_{\tau} \int_{S} \left( \alpha_{1j} p + \alpha_{2j} \frac{\partial p}{\partial y_1} + \alpha_{3j} \frac{\partial p}{\partial n} + \alpha_{4j} \frac{\partial p}{\partial \tau} \right. + \alpha_{5j} \frac{\partial^2 p}{\partial y_1 \partial \tau} + \alpha_{6j} \frac{\partial^2 p}{\partial n \partial \tau} + \alpha_{7j} \frac{\partial^2 p}{\partial \tau^2} \left) G dS d\tau$$

with:

$$\alpha_{1j} = \frac{\beta^2}{d^2} (n_j - 3 \hat{r}_n C_j \hat{r}_j)$$

$$\alpha_{2j} = -M_0^2 n_1 \frac{C_j \hat{r}_j}{d}$$

$$\alpha_{3j} = \frac{C_j \hat{r}_j}{d}$$

$$\alpha_{4j} = \frac{1}{a_0 d} (n_j - (M_0 n_1 + 3 \hat{r}_n) C_j \hat{r}_j + M_0 \hat{r}_n \delta_{1j})$$

$$\alpha_{5j} = - \frac{M_0^2 n_1}{a_0 \beta^2} (C_j \hat{r}_j - M_0 \delta_{1j})$$
$$ \alpha_{6j} = \frac{1}{a_0 \beta^2} (C_j \hat{r}_j - M_0 \delta_{1j}) $$

$$ \alpha_{7j} = -\frac{1}{a_0^2 \beta^2} (C_j \hat{r}_j - M_0 \delta_{1j}) (M_0 n_1 + \hat{r}_n) $$

In the far field, the integration can be limited to the last three terms.

## Application to Ffowcs Williams and Hawking's formulation

$$ \text{The source term } F \text{ is: } F = \frac{\beta^2}{d^2} \Sigma_1 + \frac{\partial}{\partial \tau} \left( \frac{\Sigma_1}{a_0 d} + \frac{\Sigma_2}{\beta^2} \right) \\ \Sigma_1 = \left( A_i(x_i - y_i) + \frac{(M_0^2 A_1 - U_0 B)}{\beta^2} (x_1 - y_1) \right) \quad \text{and} \\ \Sigma_2 = B - \frac{M_0 A_1}{a_0} $$

$$ A_i = pn_i + pu_i(u_n + U_0 n_1) \quad \text{and} $$

$$ B = -\rho_0 U_0 n_1 + \rho(u_n + U_0 n_1) $$

where $\vec{u}$ is the velocity of the fluid in a reference frame in which the unperturbed fluid is at rest, and $u_n = \vec{u} \cdot \vec{n}$

$$ \frac{\partial F}{\partial x_j} = \frac{\beta^2}{d^2} \left( \frac{\partial \Sigma_1}{\partial x_j} - \frac{2\Sigma_1}{d} C_j \hat{r}_j \right) + \frac{\partial}{\partial \tau} \left( \frac{1}{a_0 d} \frac{\partial \Sigma_1}{\partial x_j} - \frac{\Sigma_1}{a_0 d^2} C_j \hat{r}_j + \frac{1}{\beta^2} \frac{\partial \Sigma_2}{\partial x_j} \right) $$

$$ \frac{\partial g}{\partial x_j} \frac{\partial F}{\partial \tau} = \left( \frac{C_j \hat{r}_j - M_0 \delta_{1j}}{a_0 \beta^2} \right) \left[ \frac{\beta^2}{d^2} \frac{\partial \Sigma_1}{\partial \tau} + \frac{\partial^2}{\partial \tau^2} \left( \frac{\Sigma_1}{a_0 d} + \frac{\Sigma_2}{\beta^2} \right) \right] $$

$$ \frac{\partial d}{\partial x_j} \frac{F}{d} = C_j \hat{r}_j \left[ \frac{\beta^2}{d^3} \Sigma_1 + \frac{\partial}{\partial \tau} \left( \frac{\Sigma_1}{a_0 d^2} + \frac{\Sigma_2}{\beta^2 d} \right) \right] $$

As $A_i$ and $B$ do not depend on $x_j$, $\frac{\partial \Sigma_2}{\partial x_j} = 0$.

Grouping terms and writing $\Sigma'_j = \frac{\partial \Sigma_1}{\partial x_j} - 3C_j \hat{r}_j \frac{\Sigma_1}{d} = A_j + \frac{(M_0^2 A_1 - U_0 B)}{\beta^2} \delta_{1j} - 3C_j \hat{r}_j \frac{\Sigma_1}{d}$ finally leads to the following expression of the pressure gradient:

$$ \frac{\partial p'(\vec{x}, t)}{\partial x_j} = \int_{\tau} \int_{S} \left( F_{1j} + \frac{\partial F_{2j}}{\partial \tau} + \frac{\partial^2 F_{3j}}{\partial \tau^2} \right) G dS d\tau $$

with:

$$ F_{1j} = \beta^2 \frac{\Sigma'_j}{d^2} $$

$$ F_{2j} = \frac{\Sigma'_j}{a_0 d} + M_0 \delta_{1j} \frac{\Sigma_1}{a_0 d^2} - \frac{C_j \hat{r}_j \Sigma_2}{\beta^2 d} $$

$$ F_{3j} = - \frac{(C_j \hat{r}_j - M_0 \delta_{1j})}{\beta^2} \left( \frac{\Sigma_1}{a_0^2 d} + \frac{\Sigma_2}{a_0 \beta^2} \right) $$

In the far field, the integration can be limited to the third term.

Cite this article as: Rahier G. & Peyret C. 2024. Approximate computation of acoustic reflection and shadow effects using the Kirchhoff method. Acta Acustica, **8**, 58.