Skip to content

Commit

Permalink
Add EF_MPP model to docs/theory
Browse files Browse the repository at this point in the history
  • Loading branch information
jakep72 committed Feb 22, 2024
1 parent f80fc75 commit 93dc9ff
Show file tree
Hide file tree
Showing 4 changed files with 105 additions and 0 deletions.
99 changes: 99 additions & 0 deletions docs/Theory/EF_MPP_Model.md
Original file line number Diff line number Diff line change
@@ -0,0 +1,99 @@
# Equivalent Fluid MPP Model

The equivalent fluid microperforate model estimates the input parameters to the [JCA](https://jakep72.github.io/acoustipy/Theory/JCA_Model/) model based on the perforate diameter $(d)$, center-to-center distance $(b)$, and thickness $(t)$.

The internal [_calc_dynamics](https://jakep72.github.io/acoustipy/AcousticTMM/#src.acoustipy.TMM.AcousticTMM._calc_dynamics) method is then used to determine the dynamic mass density $(\tilde{\rho})$ and dynamic bulk modulus $(\widetilde{K})$ using the [JCA](https://jakep72.github.io/acoustipy/Theory/JCA_Model/) model.

# Estimation of JCA parameters

\[
\phi = \frac{\pi}{4}\Bigg(\frac{d}{b}\Bigg)^2
\]

\[
\sigma = \frac{32\eta}{\phi d^2}
\]

\[
\Lambda = \frac{d}{2}
\]

\[
\Lambda^{\prime} = \frac{d}{2}
\]

\[
\tau = 1+\frac{2*fok}{t}
\]

where $fok$ is:

\[
fok = \frac{4d}{3\pi}(1-1.13eps-0.09eps^2+0.27eps^3)
\]

and $eps$ is:

\[
eps = 2\sqrt{\frac{\phi}{\pi}}
\]

The [Add_MPP_EF_Layer](https://jakep72.github.io/acoustipy/AcousticTMM/#src.acoustipy.TMM.AcousticTMM.Add_MPP_EF_Layer) method then converts The modified dynamic mass density and bulk modulus to the characteristic impedence $(Z_{c})$ and wavenumber $(k_{c})$ for use in the layer transfer matrix.

\[
Z_{c} = \sqrt{\tilde{\rho}\widetilde{K}}
\]

\[
k_{c} = {\omega}\sqrt{\frac{\tilde{\rho}}{\widetilde{K}}}
\]

## Model Parameters:

##### Using the following nomenclature --- Symbol = [Units] (name)

###### Model Specific
\[
d = \Bigg[m\Bigg]\tag{perforate diameter}
\]

\[
b = \Bigg[m\Bigg]\tag{center-to-center distance}
\]

\[
t = \Bigg[m\Bigg]\tag{layer thickness}
\]

###### JCA Parameters

\[
\sigma = \Bigg[\frac{Pa*s}{m^2}\Bigg]\tag{static airflow resistivity}
\]

\[
\phi = \Bigg[unitless\Bigg]\tag{porosity}
\]

\[
\tau = \Bigg[unitless\Bigg]\tag{tortuosity}
\]

\[
\Lambda = \Bigg[{\mu}m\Bigg]\tag{viscous characteristic length}
\]

\[
\Lambda^{\prime} = \Bigg[{\mu}m\Bigg]\tag{thermal characteristic length}
\]



## Defining Other Symbols:

##### Using the following nomenclature --- Symbol = [Units] (name)


\[
\eta = \Bigg[Pa*s\Bigg]\tag{viscosity of air}
\]
2 changes: 2 additions & 0 deletions docs/Theory/JCAL_Model.md
Original file line number Diff line number Diff line change
@@ -1,3 +1,5 @@
# JCAL Model

The Johnson-Champoux-Allard-Lafarge model is a six parameter model consisting of the static airflow resistivity $(\sigma)$, porosity $(\phi)$, tortuosity $(\tau)$, viscous characteristic length $(\Lambda)$, thermal characteristic length $(\Lambda^\prime)$, and thermal permeability $(k_{0}^\prime)$.

The acoustipy implementation for the JCA, JCAL, and JCAPL models are all based on the implementation from [APMR](https://apmr.matelys.com/PropagationModels/MotionlessSkeleton/JohnsonChampouxAllardPrideLafargeModel.html). The equations described below can be found in the [_calc_dynamics](https://jakep72.github.io/acoustipy/AcousticTMM/#src.acoustipy.TMM.AcousticTMM._calc_dynamics) method.
Expand Down
2 changes: 2 additions & 0 deletions docs/Theory/JCAPL_Model.md
Original file line number Diff line number Diff line change
@@ -1,3 +1,5 @@
# JCAPL Model

The Johnson-Champoux-Allard-Pride-Lafarge model is an eight parameter model consisting of the static airflow resistivity $(\sigma)$, porosity $(\phi)$, tortuosity $(\tau)$, viscous characteristic length $(\Lambda)$, thermal characteristic length $(\Lambda^\prime)$, thermal permeability $(k_{0}^\prime)$, thermal tortuosity $(\alpha_{0}^\prime)$, and viscous tortuosity $(\alpha_{0})$.

The acoustipy implementation for the JCA, JCAL, and JCAPL models are all based on the implementation from [APMR](https://apmr.matelys.com/PropagationModels/MotionlessSkeleton/JohnsonChampouxAllardPrideLafargeModel.html). The equations described below can be found in the [_calc_dynamics](https://jakep72.github.io/acoustipy/AcousticTMM/#src.acoustipy.TMM.AcousticTMM._calc_dynamics) method.
Expand Down
2 changes: 2 additions & 0 deletions docs/Theory/JCA_Model.md
Original file line number Diff line number Diff line change
@@ -1,3 +1,5 @@
# JCA Model

The Johnson-Champoux-Allard model is a five parameter model consisting of the static airflow resistivity $(\sigma)$, porosity $(\phi)$, tortuosity $(\tau)$, viscous characteristic length $(\Lambda)$ and thermal characteristic length $(\Lambda^\prime)$.

The acoustipy implementation for the JCA, JCAL, and JCAPL models are all based on the implementation from [APMR](https://apmr.matelys.com/PropagationModels/MotionlessSkeleton/JohnsonChampouxAllardPrideLafargeModel.html).
Expand Down

0 comments on commit 93dc9ff

Please sign in to comment.