Porous absorber

Note: This page has not been verified by our editors.

Figure 1: Sketch of the geometry. The porous material is marked in blue, the acoustic source by \(q\).

Figure 2: Frequency response function.

Description

The Porous absorber benchmark models the sound pressure in a cavity excited by a single harmonic load. One side of the cavity is covered by a layer of poroelastic material, which adds dissipation to the system. The geometry of this model follows . Various projection-based model order reduction methods have been applied and compared using this example as a benchmark in .

The cavity has the dimensions \(0.75 \times 0.6 \times 0.4\,\mathrm{m}\) and one wall is covered by a \(0.05\,\mathrm{m}\) thick poroelastic layer acting as a sound absorber. The poroelastic material is described by the Biot theory and the system is excited by a point source located in a corner opposite of the porous layer. The material parameters for the acoustic fluid and the poroelastic material have been chosen according to. The transfer function measures the mean acoustic pressure inside the cavity.

Dimensions

System structure: \[ \begin{align} \left( K + \tilde{\gamma}(s) K_{p,1} + \tilde{\rho}_f(s) K_{p,2} + s^2 M + s^2 \tilde{\gamma}(s) M_{p,1} + s^2 \tilde{\rho}(s) M_{p,2} + \frac{s^2 \phi^2}{\tilde{R}(s)} M_{p,3} \right) x(s) &= B, \\ y(s) &= C x(s), \end{align} \] with the frequency dependent functions for the effective densities \(\tilde{\rho}(s), \tilde{\rho}_f(s)\), the parameter \(\tilde{\gamma}(s)\) relating the effective densities and the frequency dependent elasticity coefficients to the porosity, and the scaled effective bulk modulus \(\tilde{R}(s)\). For more details on the functions, see .

System dimensions\[K, K_{p,1}, K_{p,2}, M, M_{p,1}, M_{p,2}, M_{p,3} \in \mathbb{R}^{n \times n}\], \(B \in \mathbb{R}^{n \times 1}\), \(C \in \mathbb{R}^{1 \times n}\), with \(n=386\,076\).

Data

The data is available at Zenodo.

Remarks

• The numerical model resembles the results from in a frequency range from \(100\,\mathrm{Hz}\) to \(1000\,\mathrm{Hz}\). The frequency response in this range is also included in the dataset.
• The finite element discretization has been performed with Kratos Multiphysics.
• A comparison of different interpolation-based MOR methods using this benchmark example is available in.

Citation

To cite this benchmark, use the following references:

• For the benchmark itself and its data:

 @Misc{dataAum23,
   author =       {Aumann, Q.},
   title =        {Matrices for an acoustic cavity with poroelastic layer},
   howpublished = {hosted at {MORwiki} -- Model Order Reduction Wiki},
   year =         2023,
   doi =          {10.5281/zenodo.8087341}
 }

• For the background on the benchmark:

 @Article{AumW23,
   author =       {Aumann, Q. and Werner, S.~W.~R.},
   title =        {Structured model order reduction for vibro-acoustic problems using interpolation and balancing methods},
   journal =      {Journal of Sound and Vibration},
   volume =       543,
   year =         2023,
   pages =        {117363},
   doi =          {10.1016/j.jsv.2022.117363},
   publisher =    {Elsevier {BV}}
 }