AeroSPEED: a high order acoustic solver for aeroacoustic applications

Alberto Artoni,Paola F. Antonietti,Roberto Corradi,Ilario Mazzieri,Nicola Parolini,Daniele Rocchi,Paolo Schito,Francesco F. Semeraro

2023-03-18

Abstract:We propose AeroSPEED, a solver based on the Spectral Element Method (SEM) that solves the aeroacoustic Lighthill's wave equation. First, the fluid solution is computed employing a cell centered Finite Volume method. Then, AeroSPEED maps the sound source coming from the flow solution onto the acoustic grid, where finally the Lighthill's wave equation is solved. An ad-hoc projection strategy is adopted to apply the flow source term in the acoustic solver. A model problem with a manufactured solution and the Noise Box test case are used as benchmark for the acoustic problem. We studied the noise generated by the complex flow field around tandem cylinders as a relevant aeroacoustic application. AeroSPEED is an effective and accurate solver for both acoustics and aeroacoustic problems.

Numerical Analysis,Fluid Dynamics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop an efficient and accurate solver for simulating and analyzing aero - acoustic problems. Specifically, the AeroSPEED solver aims to solve the following problems:

1. **Development of a high - order acoustic solver**: The paper proposes a high - order acoustic solver AeroSPEED based on the Spectral Element Method (SEM) for solving the Lighthill wave equation, which is the fundamental equation for describing the noise generated by fluid flow.

2. **Handling of multi - scale problems**: Since flow and acoustic problems involve different scales, a separated approach is usually adopted to handle these problems. First, the fluid problem is solved using the Finite Volume Method (FVM), then the flow source terms are mapped onto the acoustic grid, and finally the Lighthill wave equation is solved.

3. **Noise generation in complex flow fields**: The noise generated in complex flow fields (such as turbulence around tandem cylinders) is studied. Such flow fields are highly unsteady and have complex structures, posing challenges to numerical simulations.

4. **Verification and comparison**: The accuracy and efficiency of AeroSPEED are verified through model problems and actual test cases (such as the Noise Box test of a simplified automobile cockpit), and compared with the commercial software COMSOL.

5. **Application examples**: AeroSPEED is applied to study the noise generated by complex flow fields around tandem cylinders, which is an important application example in the field of aero - acoustics.

### Formula summary

- **Lighthill tensor**:
\[
T=\rho_0u\otimes u
\]
where \(\rho_0\) is the fluid density and \(u\) is the velocity field.

- **Inhomogeneous acoustic wave equation**:
\[
\frac{\partial^2\rho}{\partial t^2}-c_0^2\Delta\rho = f,\quad\text{in }\Omega_A\times(0,T)
\]
where \(c_0\) is the speed of sound and \(f = \nabla\cdot\nabla\cdot T\) is the divergence of the Lighthill tensor.

- **Boundary conditions**:
\[
c_0^2\frac{\partial\rho}{\partial n}=0,\quad\text{on }\Gamma_B\times(0,T)
\]
\[
\frac{1}{\rho_0}\frac{\partial\rho}{\partial n}=-\frac{1}{Z}\frac{\partial\rho}{\partial t},\quad\text{on }\Gamma_Z\times(0,T)
\]
where \(Z\) is the impedance of the outer wall.

- **Time discretization**:
The Newmark method is used for time discretization:
\[
\rho_h^{k + 1}=\rho_h^k+\Delta t v_h^k+\Delta t^2\left(\beta_N A^{k + 1}+\left(\frac{1}{2}-\beta_N\right)A^k\right)
\]
\[
v_h^{k + 1}=v_h^k+\Delta t\left(\gamma_N A^{k + 1}+(1-\gamma_N)A^k\right)
\]
where \(A(\rho_h,t)=M^{-1}(f - c_0^2K\rho_h)\), and \(M\) and \(K\) are the mass respectively.