Anonymous
×
Create a new article
Write your page title here:
We currently have 106 articles on MOR Wiki. Type your article name above or click on one of the titles below and start writing!



Porous absorber: Difference between revisions

Aumann (talk | contribs)
No edit summary
Aumann (talk | contribs)
No edit summary
Line 18: Line 18:
:<math>
:<math>
\begin{align}
\begin{align}
\left( A_1 + \tilde{\gamma}(s) A_2 + \tilde{\rho}_f(s) A_3 + s^2 A_4 + s^2 \tilde{\gamma}(s) A_5 + s^2 \tilde{\rho}(s) A_6 + \frac{s^2 \phi^2}{\tilde{R}(s)} A_7 \right) x(s) &= B, \\
\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),
y(s) &= C x(s),
\end{align}
\end{align}
Line 27: Line 27:
System dimensions:
System dimensions:


<math>A_i, \in \mathbb{R}^{n \times n}</math>, with <math>i=1, \ldots, 7 </math>,  
<math>K, K_{p,1}, K_{p,2}, M, M_{p,1}, M_{p,2}, M_{p,3} \in \mathbb{R}^{n \times n}</math>,  
<math>B \in \mathbb{R}^{n \times 1}</math>,
<math>B \in \mathbb{R}^{n \times 1}</math>,
<math>C \in \mathbb{R}^{1 \times n}</math>,
<math>C \in \mathbb{R}^{1 \times n}</math>,

Revision as of 15:25, 27 June 2023

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 [1]. Various projection-based model order reduction methods have been applied and compared using this example as a benchmark in [2].

The cavity has the dimensions 0.75×0.6×0.4m and one wall is covered by a 0.05m thick poroelastic layer acting as a sound absorber. The poroelastic material is described by the Biot theory[3] 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[1].

Dimensions

System structure:

(K+γ~(s)Kp,1+ρ~f(s)Kp,2+s2M+s2γ~(s)Mp,1+s2ρ~(s)Mp,2+s2ϕ2R~(s)Mp,3)x(s)=B,y(s)=Cx(s),

with the frequency dependent functions for the effective densities ρ~(s),ρ~f(s), the parameter γ~(s) relating the effective densities and the frequency dependent elasticity coefficients to the porosity, and the scaled effective bulk modulus R~(s). For more details on the functions, see [1].


System dimensions:

K,Kp,1,Kp,2,M,Mp,1,Mp,2,Mp,3n×n, Bn×1, C1×n, with n=386076.


Data

The data is available at Zenodo.


References

  1. 1.0 1.1 1.2 R. Rumpler, P. Göransson, J.-F. Deü. "A finite element approach combining a reduced-order system, Padé approximants, and an adaptive frequency windowing for fast multi-frequency solution of poro-acoustic problems", International Journal for Numerical Methods in Engineering, 97: 759-784, 2014.
  2. Q. Aumann, S. W. R. Werner. "Structured model order reduction for vibro-acoustic problems using interpolation and balancing methods", Journal of Sound and Vibration, 543: 117363, 2023.
  3. M. A. Biot. "Theory of propagation of elastic waves in a fluid-saturated porous solid. I. Low-frequency range", J. Acoust. Soc. Am., 28(2):168–178, 1956.