**Journal of Applied Mathematics and Physics**

Vol.07 No.07(2019), Article ID:93575,23 pages

10.4236/jamp.2019.77097

Model for Frequency Dependence of Thermal Permeability in Order to Quantify the Effects of Thermal Exchange on Wave Propagation in Multi Layered Porous Medium

Allen Teagle-Hernandez^{1,2}, Ortwin Ohtmer^{1}, Duyen Nguyen^{1}^{ }

^{1}Mecahnical and Aerospace Engineering Department, California State University, Long Beach, USA

^{2}Vibro-Acoustic Department, HSEA Tech. Inc., Nagoya, Japan

Received: April 12, 2019; Accepted: July 3, 2019; Published: July 10, 2019

ABSTRACT

The ability to quantify and predict the energy absorption/transmission characteristics of multi-layered porous medium is imperative if one is involved in the automotive, launch vehicle, commercial aircraft, architectural acoustics, petroleum exploration, or even in modeling human tissue. A case in point, the first four aforementioned fields rely on effective Noise and Vibration (NV) development for their commercial success. NV development requires the setting of NV targets at different system levels. The targets are then translated to Transmission Loss (TL), Insertion Loss (IL), and absorption (Alpha) performance for the multi-layered porous materials being utilized. Thus, it behooves to have a thorough understanding of the physics behind the energy dissipating mechanism of the material that entails the effects of the fluid meandering through the pores of the material and its interaction with the structural skeleton. In this section of the project the focus is on the thermal interchange that occurs within the porous medium. Via the acoustic modeling at the micro/macro level it is shown how this thermal exchange affects the acoustic compressibility within the porous material. In order to obtain a comprehensive approach the ensuing acoustic modeling includes the effects due to relaxation process, thus bulk viscosity and instantaneous entropy functions (effects due to vibration of diatomic molecules of air) are incorporated into the equation. The instantaneous entropy functions are explained by means of the Boltzmann’s distribution, partition function, and quantum states. The concept of thermal length and its connection to thermal permeability is clarified. Lastly, the results for TL calculations employing the aforementioned thermal exchange into the Transfer Matrix Method with finite size correction, (FTMM), pertaining to a simple multi-layered material is compared with experimentally obtained data.

**Keywords:**

Porous Medium, Fluid Structure Interaction, Thermal Permeability, Thermal Length, Partition Function, Quantum Mechanics, Asymptotic Expansion, ket/bra Vectors

1. Introduction

Teagle et al. [1] derived a coupled set of fluid/structure equations for a porous medium applying asymptotic and homogenization techniques. Via this modeling it is established in [1] that there are mainly 3 modes of energy transformation: 1) The first mode is through the connection between the micro and macro structural framework of the porous skeleton, 2) the second is via the viscous boundary layer, and 3) the third interaction is through thermal (entropy) boundary layer. The combination of the viscous boundary layer and how tortuous the porous material is, results in the encapsulation of the fluid medium. This, in turn, changes the apparent mass of the structural medium. Details pertaining to these encapsulating phenomena can be found in the work by Johnson, et al. [2]. In [2] the concepts of tortuosity, viscous length, and viscous permeability are explained.

This current work concentrates in the third interaction mode. There is a boundary layer where the flow experiences changes in entropy, therefore a loss/gain of energy due to thermal exchange will occur. This thermal exchange also changes the acoustic bulk modulus at certain frequencies and in turn changing the speed of sound. This is depicted through Equation (31) found in [1]

$-i{\omega}^{*}{p}_{o}^{*}(\gamma +(\gamma -1)Pr{{k}^{\prime}}^{\ast}i{\omega}^{\ast})=-{\nabla}_{y}^{*}\cdot {v}_{1}^{*}+i{\omega}^{*}{\nabla}_{x}^{*}\cdot {u}_{o}^{*}+i{\omega}^{*}{\nabla}_{x}^{*}\cdot {w}^{*}(x,y)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}{D}_{f}$ (1)

From this equation it is understood that the compressibility factor, ${\beta}_{a}$ , is defined as

${\beta}_{a}=\gamma +\left(\gamma -1\right)i\omega \frac{Pr}{\mu \varphi}\rho {k}^{\prime}$ (2)

where ${\kappa}_{a}=\frac{1}{{\beta}_{a}}$ is the bulk modulus of the air inside the pore. ${k}^{\prime}$ is the

thermal permeability and is defined as

$\varphi {T}_{o}=-\frac{{k}^{\prime}(\omega )}{\kappa}\frac{\partial {P}_{o}}{\partial t}$ (3)

${T}_{o},{P}_{o}$ are the first terms of the asymptotic expansion (macro scale) pertaining to the temperature and pressure respectively, their equations are derived in [1]. The bulk of this work is to derive a frequency dependent general expression for the thermal permeability, ${k}^{\prime}(\omega )$ . In order to attain this goal a good understanding of the physics behind the thermal exchange process has to be accomplished. This is achieved by: 1) analyzing the intensity escaping the main acoustic mode or 2) deriving a partial differential equation incorporating the porous wall impedance due to the thermal/viscous boundary layer; both are accomplished here.

This work also models an additional energy absorbing phenomena and that is the relaxation process. This is the hysteresis effect due to compression/decompress- sion cycle of the molecules. The application of quantum mechanics, explained via the construct of Hilbert Spaces along with its dual space, facilitates the modeling of this hysteresis effect.

The main interest of the current authors is to calculate transmission loss (TL), Insertion loss (IL), and absorption coefficient (Alpha) for multi-layered porous medium. Therefore, calculation results for TL pertaining to a layered material, consisting of 2 different foams, are shown. This calculation is achieved by incorporating all of the aforementioned energy phenomena into the Finite-Cor- rection Transfer Matrix Method (FTMM) [3] [4] [5].

In addition to the goals mentioned in the last paragraph a parallel effort is underway to lay down the mathematical foundation to understand how the propagating energy transfigures as it travels through the porous medium. This creates a tool that will allow design teams to make efficient absorption comparison of porous mediums consisting of different rheological substances.

2. Formulation

2.1. Basic Statistical Mechanics: Partition Function/Boltzmann Distribution/Quantum States

One of the goals is to obtain a comprehensive set of fluid equations that incorporates any phenomena that leads to energy depletion. This means that in addition to viscous and thermal effects that transpire within their respective boundary layers, a quantification of the hysteresis effects due to the delay in thermal response should be obtained and incorporated. The most efficient tool to describe this process and for that matter any process is to apply statistical mechanic techniques. Once the Partition Function, Z, is obtained any macro level characteristic pertaining to the system can be derived. In this study Boltzmann’s

version of the combinatorial arrangement, $\stackrel{\u02dc}{C}=\frac{N!}{{{\displaystyle \prod}}_{i-1}^{m}\text{\hspace{0.05em}}{n}_{i}!}$ , is used. N represents

the number of states, and n_{i} are the number of states at each energy level. Applying Sterling’s formula, this arrangement number becomes,

$C=\mathrm{ln}(\stackrel{\u02dc}{C})=-{\displaystyle {\sum}_{i}{p}_{i}\mathrm{ln}({p}_{i})}$ , where
${p}_{i}=\frac{{n}_{i}}{N}$ , is the probability of being in the i^{th}

energy state that has energy
${E}_{i}$ . The statistical mechanic definition of Entropy, s, is the measure of information plus its distribution and is defined as finding the p_{i}s that will maximize the arrangement number, C, under the following constraints:
${\sum}_{i}{p}_{i}}=1$ , and
${\sum}_{i}{p}_{i}{E}_{i}}=E$ (E = Average Energy). From this mathematical framework, the temperature,
$\stackrel{\xaf}{T}$ , plays the role of a Lagrange Multiplier. The main results are that

${p}_{i}=\frac{{e}^{-\beta {E}_{i}}}{Z}$ (4)

where Z is the partition function defined as

$Z={\displaystyle {\sum}_{i}{e}^{-\beta {E}_{i}}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\beta =\frac{1}{\stackrel{\xaf}{\stackrel{\xaf}{T}}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{\xaf}{T}=kT\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}k=Boltzmann\text{\hspace{0.17em}}Constant$ (5)

and note that the units for $\stackrel{\xaf}{T}$ (“theorist” temperature) is in Joules. It is easily shown that the average energy and the Helmholtz Free Energy, A, are represented respectively as,

$E=-\frac{\partial \mathrm{ln}(Z)}{\partial \beta}$ and $A=E-TS=-\stackrel{\xaf}{T}\mathrm{ln}(Z)$ (6)

In quantum mechanics (see Section 2.2 below) the probability density is

defined as $p\left(r,t\right)=\psi \stackrel{\xaf}{\psi}={\left|\psi \right|}^{2}$ , where $\psi $ is a state vector of the system. Applying Schrodinger’s equation, $i\hslash \frac{\partial \Psi}{\partial t}=\stackrel{^}{H}\Psi (r,t)$ , and Hamilton’s operator,

defined in Equation (19a), the following differential equation is derived which expresses the propagation of probability

$\frac{\partial p}{\partial t}+\nabla J=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}J=\frac{i\hslash}{2m}\left[\psi \nabla \stackrel{\xaf}{\psi}-\stackrel{\xaf}{\psi}\nabla \psi \right]$

$J$ is the probability flux. For adiabatic processes the divergence of the probability flux is shown to be zero,
$\nabla J=0$ ; therefore the probability density is stationary with respect to time. This means that the p_{i}s are adiabatic invariant (which in turn means that entropy is constant) and thus the following expression for pressure holds

$\begin{array}{l}P\stackrel{adiabatic}{\to}{\frac{\partial E}{\partial V}|}_{S}=-{\frac{\partial E}{\partial V}|}_{T}+{\frac{\partial S}{\partial V}|}_{T}{\frac{\partial E}{\partial S}|}_{V}=-{\frac{\partial E}{\partial V}|}_{T}+{\frac{\partial S}{\partial V}|}_{T}T\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}={-\frac{\partial (E-TS)}{\partial V}|}_{T}={\frac{T\partial \mathrm{ln}(Z)}{\partial V}|}_{T}\end{array}$ (7)

The average potential between the particles is defined as ${U}_{0}$ and under classical mechanic conditions the partition function is

$Z={Z}_{0}(\beta )\left[1-\frac{\beta {N}^{2}{U}_{0}}{2V}\right]$ (8)

${Z}_{0}\left(\beta \right)={\left(\frac{e}{\rho}\right)}^{N}{\left(\frac{2m\pi}{\beta}\right)}^{\frac{3N}{2}}$ is the partition function pertaining to the Ideas Gas.

In order to obtain a general expression for the average energy of the gas, Equation (6) is applied

$E=-\frac{\partial \mathrm{ln}\left(Z\right)}{\partial \beta}=N\frac{3}{2}\stackrel{\xaf}{T}+\frac{{N}^{2}}{2V}{U}_{0}=\left[\frac{3}{2}\stackrel{\xaf}{T}+\frac{\stackrel{\u02dc}{\rho}}{2}{U}_{0}\right]N$ (9)

Diatomic molecules like O_{2} and N_{2} can be modeled as stiff oscillators. As a simple representative example, a model of a simple oscillator with stiffness k will be studied. The total energy for this system is

${E}_{T}=\frac{{p}^{2}}{2m}+\frac{k{x}^{2}}{2}$ (10)

Via the Boltzmann distribution, the Partition Function pertaining to the oscillator, ${Z}_{os}$ , is

${Z}_{os}={\displaystyle \iint {e}^{-\beta \frac{{p}^{2}}{2m}}{e}^{\frac{k{x}^{2}}{2}}dpdx}={\displaystyle \int {e}^{-\beta \frac{{p}^{2}}{2m}}dp}{\displaystyle \int {e}^{\frac{k{x}^{2}}{2}}dx}$ (11)

By applying Gaussian integrals the Partition Function for an undamped oscillator is

${Z}_{os}=\sqrt{\frac{2m\pi}{\beta}}\sqrt{\frac{2\pi}{\beta k}}=\frac{2\pi}{\beta}\sqrt{\frac{m}{k}}=\frac{2\pi}{{\omega}_{n}}\frac{1}{\beta}$ (12)

As before, Equation (6) is used in order to obtain the average energy of the oscillating system

$-\frac{\partial \mathrm{ln}\left({Z}_{os}\right)}{\partial \beta}=-\frac{\partial \left[\mathrm{ln}\left(\frac{2\pi}{{\omega}_{n}}\right)-\mathrm{ln}\left(\beta \right)\right]}{\partial \beta}=\frac{1}{\beta}=\stackrel{\xaf}{T}$ (13)

The analysis above, especially Equation (11), assumes classical mechanic conditions. Equation (13) states that the average energy of the oscillating system is $\stackrel{\xaf}{T}$ which, based on prior analysis, is connected to kinetic energy. This result is somewhat suspect since it states that regardless of how stiff the system is, the average energy is always $\stackrel{\xaf}{T}$ . In classical mechanics, the path integrals measured with respect to the Max Planck constant, $\hslash $ , are relatively large. An observable in classical mechanics is simply a real function in Euclidean space. In quantum mechanics every possible state of a given system corresponds to a separable Hilbert Space over the complex number field. Additionally, to every physical observable there corresponds, in the Hilbert space, a linear Hermitian operator that has a complete set of orthogonal eigenvectors. The following limit applies to systems oscillating at relatively high energy levels,

${\mathrm{lim}}_{n\to \infty}\frac{\Delta E}{{E}_{n}}={\mathrm{lim}}_{n\to \infty}\hslash \omega /\left[\left(n+\frac{1}{2}\right)\hslash \omega \right]\to 0$ thus, in classical mechanics, it can

be assumed that the energy forms a continuous spectrum; as compared to quantum mechanics where there is now a quantization (discretization) of the energy states [6]. In quantum mechanics a unit of energy is proportional to the natural frequency of the system by the constant,
$\hslash $ . The n^{th} energy state is defined as

${E}_{n}=n\hslash \omega $ (14)

Thus, the Partition Function is

${Z}_{q}={\displaystyle {\sum}_{n}{e}^{-\beta n\hslash \omega}}={\displaystyle {\sum}_{n}{\left({e}^{-\beta \hslash \omega}\right)}^{n}}=\frac{1}{1-{e}^{-\beta \hslash \omega}}$

Equation (6) is utilized once more in order to obtain the average energy of the stiff oscillator

$E=-\frac{\partial \mathrm{ln}\left(Z\right)}{\partial \beta}=-\frac{1}{Z}\frac{\partial Z}{\partial \beta}=\frac{\hslash \omega {e}^{-\beta \hslash \omega}}{1-{e}^{-\beta \hslash \omega}}$ (15)

At high temperatures $\beta \ll 1$ and therefore Equation (15) becomes

$E=\frac{\hslash \omega {e}^{-\beta \hslash \omega}}{1-{e}^{-\beta \hslash \omega}}~\frac{\hslash \omega}{\beta \hslash \omega}=\frac{1}{\beta}=\stackrel{\xaf}{T}$ (16)

Therefore at high temperature levels there is enough energy such that the quantization of energy is not required. The criteria to follow are the following

$\beta \hslash \omega >1\to Quantum\text{\hspace{0.17em}}Mechanics\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\beta \hslash \omega <1\to Classical\text{\hspace{0.17em}}Mechanics$ (17)

The cross-over point is $\hslash \omega =\stackrel{\xaf}{T}$ . Lastly, it is possible to extract the Energy’s standard deviation from the Partition Function. It is easy to show that

$\langle {E}^{2}\rangle =\frac{1}{Z}\frac{{\partial}^{2}Z}{\partial {\beta}^{2}}$ . What follows is the following definition

$\langle {\left(\Delta E\right)}^{2}\rangle =\langle {E}^{2}\rangle -{\langle E\rangle}^{2}=\frac{1}{Z}\frac{{\partial}^{2}Z}{\partial {\beta}^{2}}-\frac{1}{{Z}^{2}}{\left(\frac{\partial Z}{\partial \beta}\right)}^{2}=\frac{{\partial}^{2}\mathrm{ln}\left(Z\right)}{\partial {\beta}^{2}}=-\frac{\partial E}{\partial \beta}=-\frac{\partial E}{\partial T}\frac{\partial \stackrel{\xaf}{T}}{\partial \beta}$ (18a)

Recognize that $\frac{\partial E}{\partial \stackrel{\xaf}{T}}={C}_{V}$ , the specific heat constant at constant volume and $\frac{\partial \stackrel{\xaf}{T}}{\partial \beta}=-{\stackrel{\xaf}{T}}^{2}$

$\langle {\left(\Delta E\right)}^{2}\rangle ={\stackrel{\xaf}{T}}^{2}{C}_{V}$ or $\langle \Delta E\rangle =\stackrel{\xaf}{T}\sqrt{{C}_{V}}$ fluctuation-susceptibility Theorem (18b)

The specific heat constant at constant volume is proportional to standard deviation of the energy.

2.2. Relaxation Process: Fluid Hysteresis/Bra-Ket Vectors, Schrodinger’s Equation

The vibration due to the collision of molecules extract energy from the passing wave but releases the energy after some delay, thus the relaxation process. This delay, akin to mechanical damping systems leads to hysteretic energy losses. The internal energy can be partitioned into translated (defined relative to the average flow velocity), rotational kinetic energy and energy due to molecular natural frequency,

$u={u}_{tr}+{u}_{rot}+{\displaystyle {\sum}_{\nu}{u}_{\nu}}$

${u}_{\nu}$ quantum energy level/mass

The temperatures, T_{tr} and T_{rot}, related to translational and rotational energies are larger than the temperatures related to molecular oscillation. They usually satisfy the second criteria in Equation (17) and therefore can be considered using Classical Theory. For the molecular vibration case,
$\hslash \omega >k{T}_{\nu}$ , and, based on Equation (17), quantum mechanics should be applied.

The first postulate in quantum theory states that everything that can be known about the state of a system can be extracted from its state vector (wave function) and this is represented as a vector in Hilbert Space; these are denoted as ket-vectors $|\Psi (x)\rangle $ . The bra-vectors are denoted as $\langle \Psi (x)|$ , the vector space of bras is the dual space(one-to-one correspondence with the space of functionals) to the kets vector space. The Riesz representation theorem allows the following notation $\langle \Psi |\Psi \rangle $ .

The 1-D Hamilton operator is given as

$\stackrel{^}{H}=\frac{{\stackrel{^}{p}}^{2}}{2m}+V\left(\stackrel{^}{x}\right)=-\frac{{h}^{2}}{2m}\frac{{\partial}^{2}}{\partial {x}^{2}}+V\left(\stackrel{^}{x}\right)$ (19)

where $\stackrel{^}{p}$ is the Linear Momentum operator, $\stackrel{^}{p}=-i\hslash \frac{\partial}{\partial x}$ , and $\stackrel{^}{x}$ is the position

operator.
$\stackrel{^}{p}$ and
$\stackrel{^}{x}$ are both observables and thus are hermitian operators on the Hilbert Space, therefore
$\stackrel{^}{H}$ has a purely discrete spectrum of eigenvalues,
${E}_{n}$ , that has a complete set of eigenstates,
$|{E}_{n}(x)\rangle $ , i.e. (
$\stackrel{^}{H}|{E}_{n}(x)\rangle ={E}_{n}|{E}_{n}(x)\rangle $ ). Hence, any state or ket-vector,
$|\Psi (x)\rangle $ , can be expanded by the set of eigenstates,
$|\Psi (x)\rangle ={\displaystyle {\sum}_{n}{a}_{n}|{E}_{n}(x)\rangle}$ . Due to Schrodinger’s equation the eigenstates pertaining to the energy operator are stationary with respect to time. For a diatomic molecule the oscillating like characteristics between the 2 atoms can be simply represented by a spring constant K. Given a particular natural frequency, f_{e}, K can be represented as
$K=4{\pi}^{2}{f}_{e}^{2}m$ . The potential energy pertaining to the simple oscillator is represented as

$\u0424\left(x\right)=-{\displaystyle {\int}_{0}^{\infty}Fdx}={\displaystyle {\int}_{0}^{\infty}Kxdx}={\displaystyle {\int}_{0}^{\infty}4{\pi}^{2}{f}_{e}^{2}mxdx}=2{\pi}^{2}{f}_{e}^{2}m{x}^{2}$

Therefore the potential operator becomes $\stackrel{^}{V}\left(\stackrel{^}{x}\right)=\frac{1}{2}m{\omega}^{2}{\stackrel{^}{x}}^{2}$ . Hamilton’s operator can be pseudo separated as follows

$\begin{array}{c}\stackrel{^}{a}=\sqrt{\frac{m}{2}}\omega \stackrel{^}{x}+\frac{i}{\sqrt{2m}}\stackrel{^}{p}=\frac{\sqrt{\hslash \omega}}{2}\left[{\left(\hslash /m\omega \right)}^{1/2}\frac{\partial}{\partial x}+{\left(m\omega /\hslash \right)}^{1/2}\stackrel{^}{x}\right]\\ =\frac{\sqrt{\hslash \omega}}{2}\left(\frac{\partial}{\partial \eta}+\stackrel{^}{\eta}\right)\end{array}$ (20a)

$\begin{array}{c}{\stackrel{^}{a}}^{*}=\sqrt{\frac{m}{2}}\omega \stackrel{^}{x}-\frac{i}{\sqrt{2m}}\stackrel{^}{p}=\frac{\sqrt{\hslash \omega}}{2}\left[-{\left(\hslash /m\omega \right)}^{1/2}\frac{\partial}{\partial x}+{\left(m\omega /\hslash \right)}^{1/2}\stackrel{^}{x}\right]\\ =\frac{\sqrt{\hslash \omega}}{2}\left(-\frac{\partial}{\partial \eta}+\stackrel{^}{\eta}\right)\end{array}$ (20b)

where $\stackrel{^}{\eta}={\left(m\omega /\hslash \right)}^{1/2}\stackrel{^}{x}$ . It can be shown that the commutator $[\stackrel{^}{a},{\stackrel{^}{a}}^{*}]=\hslash \omega $ and therefore

$\stackrel{^}{H}={\stackrel{^}{a}}^{*}\stackrel{^}{a}+\frac{1}{2}\hslash \omega =\stackrel{^}{a}{\stackrel{^}{a}}^{*}-\frac{1}{2}\hslash \omega $ (20c)

Based on Equation (20c) the eigenvalue equation can be written as

${\stackrel{^}{a}}^{*}\stackrel{^}{a}|{E}_{n}(x)\rangle =\left({E}_{n}-\frac{1}{2}\hslash \omega \right)|{E}_{n}(x)\rangle ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{^}{a}{\stackrel{^}{a}}^{*}|{E}_{n}(x)\rangle =\left({E}_{n}+\frac{1}{2}\hslash \omega \right)|{E}_{n}(x)\rangle $ (21)

Built on the ideas in Equation (21) the concept of annihilators, $\stackrel{^}{a}$ , and creators, ${\stackrel{^}{a}}^{*}$ , are established; the following relationships exists

$\stackrel{^}{a}|{E}_{n}(x)\rangle =|{E}_{n-1}(x)\rangle $ corresponding eigenvalue ${E}_{n-1}={E}_{n}-\hslash \omega $ } annihilator

${\stackrel{^}{a}}^{*}|{E}_{n}(x)\rangle =|{E}_{n+1}(x)\rangle $ corresponding eigenvalue ${E}_{n+1}={E}_{n}+\hslash \omega $ } creator

Another possible solution, tied to the annihilator equation, is that $\stackrel{^}{a}|{E}_{0}(x)\rangle =0$ . Here, $|{E}_{0}(x)\rangle $ is known as the zero-point wave-function and ${E}_{0}$ is the zero-point energy level. Equation (20a) pertaining to the zero-point state results in

$\frac{\sqrt{\hslash \omega}}{2}\left(\frac{\partial}{\partial \eta}+\stackrel{^}{\eta}\right){\psi}_{0}=0$ , $\therefore \frac{\partial {\psi}_{0}}{\partial \eta}+\stackrel{^}{\eta}{\psi}_{0}=0$ normalized solution ${\psi}_{0}={\left(\frac{m\omega}{\pi \hslash}\right)}^{1/4}{e}^{-{\eta}^{2}/2}$

The eigenstates pertaining to the higher energy levels can be generated by applying the Riesz representation notation.

$\begin{array}{c}\langle {\psi}_{n+1}|{\psi}_{n+1}\rangle =\langle {\stackrel{^}{a}}^{*}{\psi}_{n}|{\stackrel{^}{a}}^{*}{\psi}_{n}\rangle =\langle {\psi}_{n}|\stackrel{^}{a}{\stackrel{^}{a}}^{*}{\psi}_{n}\rangle =\langle {\psi}_{n}|\left(\stackrel{^}{H}+\frac{\hslash \omega}{2}\right){\psi}_{n}\rangle \\ =\hslash \omega (n+1)\langle {\psi}_{n}|{\psi}_{n}\rangle \end{array}$

Thus, if ${\psi}_{n}$ is normalized so is

${\psi}_{n+1}={\left(\hslash \omega (n+1)\right)}^{-1/2}{\stackrel{^}{a}}^{*}{\psi}_{n}={\left(n+1\right)}^{-1/2}\frac{1}{2}\left(-\frac{\partial}{\partial \eta}+\stackrel{^}{\eta}\right){\psi}_{n}$

As a consequence, the recurring formula for the eigenfunctions/states becomes

${\psi}_{n}={\left(2n!\right)}^{-1/2}{\left(-\frac{d}{d\eta}+\eta \right)}^{n}{e}^{-\frac{{\eta}^{2}}{2}}={\left(2n!\right)}^{-\frac{1}{2}}{\left(\frac{m\omega}{\pi \hslash}\right)}^{\frac{1}{4}}{H}_{n}\left(\eta \right){e}^{-\frac{{\eta}^{2}}{2}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=0,1,2,\dots $ (22)

where ${H}_{n}\left(\eta \right)$ are the Hermite polynomials. Additionally due to Equation

(19a) the pertaining eigenvalues of $\stackrel{^}{H}$ are ${E}_{n}=\hslash \omega \left(n+\frac{1}{2}\right)$

An alternative approach is to apply Schrodinger’s Equation $i\hslash \frac{\partial \Psi}{\partial t}=\stackrel{^}{H}\Psi (r,t)$ .

Plugging in the definition for the Hamilton Operator defined in Equation (19a) results in the following simple oscillator equation

$\frac{{d}^{2}{\Psi}_{v}}{d{x}^{2}}+\frac{2m}{{\hslash}^{2}}\left({E}_{j}-\Phi \left(x\right)\right){\Psi}_{v}=\frac{{d}^{2}{\Psi}_{v}}{d{x}^{2}}+\frac{2m}{{\hslash}^{2}}\left({E}_{j}-\frac{1}{2}m{\left(\omega x\right)}^{2}\right){\Psi}_{v}=0$ (23)

Setting

$z=\sqrt{\frac{m\omega}{\hslash}}x$ , $\alpha =\frac{2m{E}_{j}}{{\hslash}^{2}}$ , and ${\Psi}_{v}={e}^{\frac{-{z}^{2}}{2}}w$ (24)

Schrodinger’s Equation simplifies to

$\frac{{d}^{2}w}{d{z}^{2}}-2z\frac{dw}{dz}+\left(\frac{\alpha}{{\left(\frac{dz}{dx}\right)}^{2}}-1\right)w=\frac{{d}^{2}w}{d{z}^{2}}-2z\frac{dw}{dz}+\left(\frac{2{E}_{j}}{\hslash \omega}-1\right)w=0$ (25)

When $\frac{2{E}_{j}}{\hslash \omega}-1=2j=0,1,2,3,\dots $ the solutions to this equation are the Hermite

Polynomials. Note that the solution is equivalent to Equation (22). Applying this definition and solving for ${E}_{j}$ , the energy levels are defined as

${E}_{j}=\hslash \omega \left(j+\frac{1}{2}\right)=j=0,1,2,\dots $ (26)

Due to the assumption of independence of the different energies the partition function pertaining to internal energies can be represented as

${Z}_{int}={Z}_{rot}{Z}_{v}{Z}_{chem}$ (27)

The partition function for the vibrational portion pertaining to a particular type of molecule v (e.g. O_{2}, N_{2}) is represented as

${Z}_{v}={\displaystyle {\sum}_{{v}_{j}}{e}^{-\left[\frac{\left({\u03f5}_{{v}_{j}}-{\u03f5}_{{v}_{0}}\right)}{k{T}_{v}}\right]}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\beta =\frac{1}{\stackrel{\xaf}{\stackrel{\xaf}{T}}}=\frac{1}{k{T}_{v}}$ (28)

${\u03f5}_{{v}_{0}}$ is the ground vibrational state energy (obtained by setting j = 0 in Equation (26)). Using Equation (26) in (28), the vibrational partition function can be rewritten as

${Z}_{v}={\displaystyle {\sum}_{j=0}^{\infty}{e}^{-\left[\frac{\hslash \omega}{k{T}_{v}}\right]j}}=\frac{1}{1-{e}^{-\left[\frac{\hslash \omega}{k{T}_{v}}\right]}}=\frac{1}{1-{e}^{-\frac{{\theta}_{e}^{*}}{{T}_{v}}}}$ (29)

Such that ${\theta}_{e}^{*}=\frac{\hslash \omega}{k}$ .

Applying Equation (6), an expression for the average energy is obtained

${u}_{int}=R{T}^{2}{\left(\frac{\partial \mathrm{ln}{Z}_{int}}{\partial T}\right)}_{V}$ (30)

Therefore the vibrational internal energy becomes

$\left({u}_{v}-{u}_{{v}_{0}}\right)=-R{T}_{v}{}^{2}\left(\frac{{e}^{-\frac{{\theta}_{e}^{*}}{{T}_{v}}}}{1-{e}^{-\frac{{\theta}_{e}^{*}}{{T}_{v}}}}\left(-\frac{{\theta}_{e}^{*}}{{T}_{v}{}^{2}}\right)\right)=\frac{R{\theta}_{e}^{*}}{{e}^{\frac{{\theta}_{e}^{*}}{{T}_{v}}}-1}$ (31)

Differentiating this equation with respect to temperature one obtains [7] [8]

$\frac{d{u}_{v}}{d{T}_{v}}={C}_{{v}_{v}}=\frac{{\theta}_{e}^{*}{}^{2}}{{T}_{v}^{2}}\frac{R{e}^{\frac{{\theta}_{e}^{*}}{{T}_{v}}}}{{\left({e}^{\frac{{\theta}_{e}^{*}}{{T}_{v}}}-1\right)}^{2}}$ (32)

Applying the Boltzmann distribution definition for entropy

$s=-{\displaystyle {\sum}_{i}{p}_{i}\mathrm{ln}\left({p}_{i}\right)}=\beta E+\mathrm{ln}(Z(\; \beta \; )\; )$

results in

${s}_{v}=R\left(-\frac{\mathrm{ln}\left(1-{e}^{-\frac{{\theta}_{e}^{*}}{{T}_{v}}}\right)}{R}+\frac{\frac{{\theta}_{e}^{*}}{k{T}_{v}}}{{e}^{\frac{{\theta}_{e}^{*}}{{T}_{v}}}-1}\right)$ (33)

Recall that ${T}_{v}\ll {T}_{tr},{T}_{rot}$ and the temperatures pertaining to Translational and rotational energy levels are high enough to be considered by classical

mechanics. The combination of kinetic energies $\left(\frac{{p}^{2}}{2m}+\frac{{L}^{2}}{2I}\right)$ and positions are

considered as possible states, the application of the Boltzmann distribution results in the following equation for average energies.

$-{\sigma}_{n}=p-{\mu}_{B}\nabla \cdot v$ such that ${\mu}_{B}=\frac{p}{{\beta}_{rot}{N}_{c}}\frac{{u}_{rot}^{2}}{{u}^{2}}=bulk\text{\hspace{0.17em}}viscosity$ (34)

The instantaneous change of entropy can be represented in the following form

$Tds=du+pd{\rho}^{-1}+{\displaystyle {\sum}_{\upsilon}{A}_{\upsilon}d{T}_{\upsilon}}$ (35)

The $\sum}_{\upsilon}{A}_{\upsilon}d{T}_{\upsilon$ terms represent the small deviation from equilibrium. The affinities ${A}_{\upsilon}$ are to be determined. The portion devoid of vibration energy is denoted from now on with subscript fr (frozen), and the portion pertaining to vibrational energy is indicated with subscript v [9].

$s={s}_{fr}\left({u}_{tr}+{u}_{rot},{\rho}^{-1}\right)+{\displaystyle {\sum}_{\upsilon}{s}_{\upsilon}({T}_{\upsilon})}$ (36)

$Td{s}_{fr}=d\left({u}_{tr}+{u}_{rot}\right)+pd{\rho}^{-1}$ (37)

Set

${u}_{tr}+{u}_{rot}={C}_{V,fr}T$ assuming ideal gas $P=\rho RT$ (38)

${C}_{V,fr}=\frac{R}{\gamma -1}$ , ${C}_{V,fr}$ = coefficient of specific heat at constant volume. (39)

Substituting Equation (38) into (37) and solve for ${s}_{fr}$ the following expression is obtained

${s}_{fr}={C}_{V,fr}\mathrm{ln}\left(u-{\displaystyle {\sum}_{\upsilon}{u}_{\upsilon}}\right)+R\mathrm{ln}{\rho}^{-1}+const\text{\hspace{0.17em}}recall\text{\hspace{0.17em}}that\text{\hspace{0.17em}}u-{\displaystyle {\sum}_{\upsilon}{u}_{\upsilon}}={u}_{tr}+{u}_{rot}$ (40)

Utilizing Equation (37) - (40) and applying them into Equation (35) yields

${A}_{\upsilon}=T\left(\frac{-{C}_{V,fr}}{{u}_{tr}+{u}_{rot}}\frac{d{u}_{\upsilon}}{d{T}_{\upsilon}}+\frac{1}{{T}_{\upsilon}}\frac{d{u}_{\upsilon}}{d{T}_{\upsilon}}\right)=\left(\frac{T}{{T}_{\upsilon}}-1\right)\frac{d{u}_{\upsilon}}{d{T}_{\upsilon}}=\left(\frac{T}{{T}_{\upsilon}}-1\right){C}_{V,v}$ (41)

The expression in Equation (32) can be substituted into Equation (41). Since ${T}_{\upsilon}\ll {\theta}_{e}^{*}$ , ${C}_{V,v}$ can deduce to

${C}_{V,v}=\frac{d{u}_{\upsilon}}{d{T}_{\upsilon}}=\frac{{n}_{\upsilon}}{n}R{\left(\frac{{\theta}_{e}^{*}}{{T}_{\upsilon}}\right)}^{2}{e}^{-\frac{{\theta}_{e}^{*}}{{T}_{\upsilon}}}$ (42)

${C}_{V,v}$ act as the specific heat pertaining to each molecular type v.

The inclusion of these energy terms into the standard energy equations results in

$\frac{Du}{Dt}-{\sigma}_{n}\frac{D{\rho}^{-1}}{Dt}=T\frac{D{s}_{v}}{Dt}-{\displaystyle {\sum}_{\upsilon}{s}_{\upsilon}\frac{D{T}_{\upsilon}}{Dt}}-{\mu}_{B}\nabla \cdot v\frac{D{\rho}^{-1}}{Dt}$ (43)

Applying the conservation of mass and the heat flux the entropy balance equations become

$\rho \frac{Ds}{Dt}+\nabla \cdot \frac{q}{T}=\frac{1}{T}{\mu}_{B}{\left(\nabla \cdot v\right)}^{2}+\frac{1}{2T}\mu \stackrel{\leftrightarrow}{D\nabla v}:\stackrel{\leftrightarrow}{D\nabla v}+\frac{\kappa}{{T}^{2}}{\left(\nabla T\right)}^{2}+\frac{\rho}{T}{\displaystyle {\sum}_{\upsilon}{A}_{\upsilon}\frac{D{T}_{\upsilon}}{Dt}}$ (44)

The Navier Stokes Equation with the inclusion of the bulk viscosity becomes

$\rho \frac{Dv}{Dt}=-\nabla p+\nabla \left({\mu}_{B}\nabla \cdot v\right)+\mu \nabla \cdot \stackrel{\leftrightarrow}{D\nabla v}$ (45)

Applying conservation of energy results in the following expression for the vibration energy

$\frac{D{u}_{\upsilon}}{Dt}={C}_{V,v}\frac{D{T}_{\upsilon}}{Dt}={n}_{\upsilon}{N}_{c\upsilon}\Delta {\u03f5}_{\upsilon}$ (46)

${N}_{c\upsilon}$ represents the number of collisions a molecule of type $\upsilon $ has per unit time, ${n}_{\upsilon}$ is the number of molecules of type $\upsilon $ , and $\Delta {\u03f5}_{\upsilon}$ is the average energy gained per collision. It is known that most of the high temperature resides in the translational and rotational portions, therefore $\Delta {\u03f5}_{\upsilon}$ can best be expressed as

$\Delta {\u03f5}_{\upsilon}={\beta}_{\upsilon}k(T-{T}_{\upsilon})$ (47)

and yields

$\frac{D{T}_{\upsilon}}{Dt}=\left(T-{T}_{\upsilon}\right)\frac{1}{{\tau}_{\upsilon}}$ such that ${\tau}_{\upsilon}=\frac{{C}_{V,v}}{{n}_{\upsilon}{N}_{c\upsilon}k{\beta}_{\upsilon}}$ (48)

More elaborate modeling for $\Delta {\u03f5}_{\upsilon}$ can be achieved by applying Discrete Boltzmann Equations but these are outside the scope of this project [10]. Solving the differential equation for ${T}_{\upsilon}$ results in

${T}_{\upsilon}=T-(T-{T}_{\upsilon}(0)){e}^{-t/{\tau}_{\upsilon}}$ (49)

If T is changed by an amount $\Delta T$ it will take a time of ${\tau}_{\upsilon}$ for the incremental change in ${T}_{\upsilon}$ is $(1-{e}^{-1})\Delta t$ . Therefore ${\tau}_{\upsilon}$ is called the relaxation time. And it is this relaxation time that introduces the hysteretic energy losses [11].

The energy conservation dissipation equation is in the form of

$\frac{\partial w}{\partial t}+\nabla \cdot I=-D$

where

$w=\frac{1}{2}{\rho}_{o}{v}^{2}+\frac{1}{2}\frac{{p}^{2}}{{\rho}_{o}{c}^{2}}+\frac{1}{2}{\left(\frac{\rho T}{{C}_{p}}\right)}_{o}{s}_{fr}^{2}+{\displaystyle {\sum}_{\upsilon}\frac{1}{2}{\left(\frac{\rho {C}_{V,v}}{T}\right)}_{o}{T}_{\upsilon}^{2}}$ (50a)

$I=pv-{\mu}_{B}v\left(\nabla \cdot v\right)-\mu v\cdot \left(\stackrel{\leftrightarrow}{D\nabla v}\right)-\kappa {T}_{0}^{-1}T\nabla T$ (50b)

$\begin{array}{l}D={\mu}_{B}{\left(\nabla \cdot v\right)}^{2}+\frac{1}{2}\mu \left(\stackrel{\leftrightarrow}{D\nabla v}\right):\left(\stackrel{\leftrightarrow}{D\nabla v}\right)+\kappa {T}_{0}^{-1}{\left(\nabla T\right)}^{2}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+{\displaystyle {\sum}_{\upsilon}{\left(\frac{{\rho}_{o}{C}_{V,v}}{T{\tau}_{\upsilon}}\right)}_{o}{\left(T-{T}_{\upsilon}\right)}^{2}}\end{array}$ (50c)

From the coupled system of acoustic equation, asymptotic techniques result in different dispersion relationships for 3 different regimes; for the viscous and entropy boundary layers, and the equations pertaining to the acoustic mode. These translate to different set of partial differential equations pertaining to the 3 different regimes. The vibrational relaxation terms will mostly affect the acoustic mode regime. In order to obtain the new dispersion relationship a solution in the form of ${e}^{ikx}{e}^{-i\omega t}$ is plugged into the equation of motion (see Appendix A.1 - A.5), the ensuing is obtained

${k}^{2}\stackrel{^}{p}=\left({\omega}^{2}+i\omega {k}^{2}\frac{{\mu}_{B}+\frac{4}{3}\mu}{{\rho}_{0}}\right)\stackrel{^}{\rho}$ (51a)

$i\omega {\rho}_{0}{\stackrel{^}{s}}_{fr}=\left[{\left(\frac{\kappa}{T}\right)}_{o}{k}^{2}-{\left(\frac{i\omega \rho}{T}\right)}_{0}{\displaystyle {\sum}_{v}\frac{{C}_{V,v}}{1-i\omega {\tau}_{v}}}\right]\stackrel{^}{T}$ (51b)

$\stackrel{^}{\rho}=\frac{\stackrel{^}{p}}{{c}^{2}}-\frac{{\left(\beta T/{C}_{p}\right)}_{0}^{2}\left({\stackrel{^}{s}}_{fr}/\stackrel{^}{T}\right)\stackrel{^}{p}}{1-{\left(T/{C}_{P}\right)}_{0}\left({\stackrel{^}{s}}_{fr}/\stackrel{^}{T}\right)}$ (51c)

Plugging this last equation into Equation (51a) the dispersion relation is obtained [12]

$\frac{{k}^{2}}{{\omega}^{2}+i\omega {k}^{2}\frac{{\mu}_{B}+\frac{4}{3}\mu}{{\rho}_{0}}}=\frac{1}{{c}^{2}}-\frac{\left(\gamma -1\right){T}_{0}\left({\stackrel{^}{s}}_{fr}/T\right)/{C}_{p}{c}^{2}}{1-{\left(T/{C}_{p}\right)}_{0}\left({\stackrel{^}{s}}_{fr}/\stackrel{^}{T}\right)}$ (51d)

For the acoustic mode and the frequency range of interest the following

approximations exist: ${k}^{2}={\left(\frac{\omega}{c}\right)}^{2},\frac{{C}_{V,v}}{{C}_{p}}\ll 1,\frac{\omega \kappa}{{\rho}_{0}{C}_{p}{c}^{2}}\ll 1,\frac{\omega \mu}{{\rho}_{0}{c}^{2}}\ll 1$ the last 2

expressions can be found in Kirchhoff’s dispersion relation and were considered small in [14]. Given these assumptions, the following asymptotic expansion is obtained

$k=\frac{\omega}{c}+i\frac{{\omega}^{2}\mu}{2{\rho}_{0}{c}^{3}}\left[\frac{4}{3}+\frac{{\mu}_{B}}{\mu}+\frac{\left(\gamma -1\right)\kappa}{{C}_{p}\mu}\right]+\frac{1}{2}\left(\gamma -1\right)\frac{\omega}{c}{\displaystyle {\sum}_{v}\frac{{C}_{V,v}/{C}_{p}}{1-i\omega {\tau}_{v}}}$ (51e)

The absorption coefficients are the imaginary portion of Equation (51e)

$\frac{{\omega}^{2}\mu}{2{\rho}_{0}{c}^{3}}\left[\frac{4}{3}+\frac{{\mu}_{B}}{\mu}+\frac{\left(\gamma -1\right)\kappa}{{C}_{p}\mu}\right]+\frac{1}{2}\left(\gamma -1\right)\frac{{\omega}^{2}}{c}{\displaystyle {\sum}_{v}\frac{{\tau}_{v}{C}_{V,v}/{C}_{p}}{1+{\left(\omega {\tau}_{v}\right)}^{2}}}$ (51f)

where

${\alpha}_{v}=\frac{\pi}{2\lambda}\left(\gamma -1\right)\frac{{\omega}^{2}}{c}\frac{{\tau}_{v}{C}_{V,v}/{C}_{p}}{1+{\left(\omega {\tau}_{v}\right)}^{2}}$ (51g)

Recall that $\lambda =\frac{2\pi c}{\omega}$ . At relatively low frequencies, $\omega \ll \frac{1}{{\tau}_{\upsilon}}$ Equation (51g) indicates that

${\alpha}_{v}\to \frac{{\tau}_{\upsilon}}{c}\frac{{C}_{V,v}}{{C}_{p}}\left(\gamma -1\right){\omega}^{2}$ quadratic increase w.r.t. $\omega $ (51h)

and at relatively high frequencies, $\omega \gg \frac{1}{{\tau}_{\upsilon}}$ , results in the following

${\alpha}_{v}\to \frac{{\tau}_{\upsilon}}{c}\frac{{C}_{V,v}}{{C}_{p}}\left(\gamma -1\right)$ approaches a constant value (51i)

2.3. Energy Loss: Thermal Boundary Layer-Bessel Functions

In [1], multi-scales technique is used on the coupled fluid/structure set of equations. Equation (16), in that publication, represents the partial differential equation that connects fluctuations of temperature and pressure in the fluid medium. The dimensional form of this equation is

${\nabla}_{y}{}^{2}{T}_{0}+i\omega {\rho}_{0}\frac{Pr}{\upsilon}{T}_{0}=\frac{i\omega {P}_{0}}{\kappa}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}{D}_{f}$ (52)

Before deriving the thermal exchange for a general porous medium (with pores of arbitrary geometry), the thermal exchange for a typical cylindrical pore with cross-sectional radius R will be analyzed. The form this solution takes will be used as a pattern to emulate when considering solutions for the general case. Equation (52) in cylindrical form is

$\frac{1}{r}\frac{d}{dr}\left(r\frac{d{T}_{0}}{dr}\right)+i\omega {\rho}_{0}\frac{Pr}{\upsilon}{T}_{0}=\frac{i\omega {P}_{0}}{\kappa}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{in}\text{\hspace{0.17em}}{D}_{f}$ (53)

The Bessel function generator is defined as

${e}^{ix\mathrm{sin}\left(\theta \right)}={J}_{0}\left(x\right)+i2{J}_{1}\left(x\right)\mathrm{sin}\theta +2{J}_{2}\left(x\right)\mathrm{cos}2x+\dots $ (54a)

The ${J}_{0}$ Bessel function is efficiently represented via Hankel’s integral

${J}_{0}\left(x\right)=\frac{1}{\pi}{\displaystyle {\int}_{-1}^{1}\frac{{e}^{ixt}}{\sqrt{\left(1-{t}^{2}\right)}}dt}=\frac{2}{\pi}{\displaystyle {\int}_{0}^{1}\frac{\mathrm{cos}\left(xt\right)}{\sqrt{\left(1-{t}^{2}\right)}}dt}$ (54b)

It can be deduced from Equation (54a) that

$\frac{2n}{x}{J}_{n}\left(x\right)={J}_{n-1}\left(x\right)+{J}_{n+1}\left(x\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}2\frac{\partial {J}_{n}\left(x\right)}{\partial x}={J}_{n-1}\left(x\right)-{J}_{n+1}\left(x\right)$ (55a)

The ensuing 2 equations are obtained by combining the equations in (55a)

${J}_{n-1}\left(x\right)=\frac{n}{x}{J}_{n}\left(x\right)+\frac{\partial {J}_{n}\left(x\right)}{\partial x}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{J}_{n+1}\left(x\right)=\frac{n}{x}{J}_{n}\left(x\right)-\frac{\partial {J}_{n}\left(x\right)}{\partial x}$ (55b)

Applying Equation (55b) results in

$\frac{d}{dx}\left\{x{J}_{1}(x)\right\}=x{J}_{0}(x)$ (55c)

Equation (55c) is integrated from 0 to R, this generates an equation that will be utilized shortly

${\int}_{0}^{R}x{J}_{0}\left(\gamma x\right)dx}=\frac{R}{\gamma}{J}_{1}\left(\gamma \right)$ (55d)

Teagle [1] showed that ${P}_{0}$ is independent of the micro (fast) scale y. Therefore

the particular solution for Equation (53) is $\frac{\upsilon}{\kappa Pr}{P}_{0}=\frac{{\rho}_{0}}{{C}_{p}}{P}_{0}$ . The homogeneous solution is ${A}_{1}{J}_{0}\left(\sqrt{\frac{i\omega {\rho}_{0}Pr}{\mu}}r\right)$ , where ${A}_{1}$ is a derived so that boundary conditions

at the wall of the cylinder pores are met, ${T}_{0}\left(R\right)=0$ . Combining the particular and homogenous solutions and satisfying boundary conditions, the following expression for temperature is obtained

${T}_{0}\left(r\right)=\frac{{\rho}_{0}}{{C}_{p}}{P}_{0}\left(1-\frac{{J}_{0}\left(\chi r\right)}{{J}_{0}\left(\chi R\right)}\right)\ni \chi =\sqrt{\frac{i\omega {\rho}_{0}Pr}{\mu}}$ (56)

For statistical energy analysis calculations, average temperature of the pore cross-section suffices, therefore

$\begin{array}{c}\langle {T}_{0}\rangle =\frac{{\displaystyle {\int}_{0}^{R}2\pi {T}_{0}\left(r\right)rdr}}{\pi {R}^{2}}=\frac{{\rho}_{0}}{{C}_{p}}{P}_{0}\left(1-\frac{2{J}_{1}\left(\chi R\right)}{R\chi {J}_{0}\left(\chi R\right)}\right)\\ =\frac{{\rho}_{0}}{{C}_{p}}{P}_{0}\left(1-\left(1-i\right)\frac{2{J}_{1}\left(\chi R\right)}{R{J}_{0}\left(\chi R\right)}\sqrt{\frac{\mu}{\omega {\rho}_{0}Pr}}\right)\end{array}$ (57)

For the more general case, the focus is on analyzing the energy flux escaping the acoustic mode into the entropy boundary layer. Johnson et al. [2] accomplished similar work for the viscous boundary layer.

Trilling [13] gave a detailed explanation of the separation technique utilized to obtain 3 dispersion relations from the Navier-Stokes Equation. Applying asymptotic expansion for each of these relationships, it is possible to derive 3 sets of simple partial differential equations along with their pertinent conditions. The process polarized the linear Navier-Stokes equation into 3 different flow regimes. The superposition principle can be applied to obtain the total velocity. The flow regimes or modes can be characterized by the following:

Vorticity Mode field

${\nabla}^{2}{v}_{vor}=\frac{\rho}{\mu}\frac{\partial {v}_{vor}}{\partial t}$ , $\nabla \cdot {v}_{vor}=0$ , ${s}_{vor}={p}_{vor}={T}_{vor}=0$ (58)

It is important to point out that Equation (58) is a classical diffusion equation

The entropy Mode field

${\nabla}^{2}{s}_{ent}=\frac{\rho {c}_{P}}{\kappa}\frac{\partial {s}_{ent}}{\partial t}$ (59a)

$\begin{array}{l}{p}_{ent}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{v}_{ent}={\left(\frac{\beta T\kappa}{\rho {c}_{P}}\right)}_{o}\nabla {s}_{ent},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\nabla \times {v}_{ent}=0,\\ {{T}^{\prime}}_{ent}=\left(\frac{T}{{c}_{P}}\right){s}_{ent},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{{\rho}^{\prime}}_{ent}={\left(\frac{\beta T\rho}{{c}_{P}}\right)}_{o}{s}_{ent}\end{array}$ (59b)

Acoustic Mode field

$\rho \frac{\partial {v}_{ac}}{\partial t}=-\nabla {p}_{ac},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{s}_{ac}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\nabla \times {v}_{ac}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{{T}^{\prime}}_{ac}={\left(\frac{\beta T}{\rho {c}_{P}}\right)}_{o}{p}_{ac},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{{\rho}^{\prime}}_{ac}=\frac{{p}_{ac}}{{c}^{2}}$ (60)

Due to the superposition principle the microscopic velocity flow field can be written as

$v={v}_{vor}+{v}_{ent}+{v}_{ac}$ (61)

The dispersion for the vorticity and entropy modes are

${k}_{vor}^{2}=\frac{i\omega \rho}{\mu}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{k}_{ent}^{2}=\frac{i\omega \rho {c}_{P}}{\kappa}$ (62)

Therefore the vorticity and entropy mode field extend to $\frac{1}{\left|{k}_{vor}\right|}$ and $\frac{1}{\left|{k}_{ent}\right|}$ so

${\delta}_{vor}={\left(\frac{2\mu}{\omega \rho}\right)}^{1/2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\delta}_{ent}={\left(\frac{2\kappa}{\omega \rho {c}_{P}}\right)}^{1/2}=\frac{{\delta}_{vor}}{{\left(Pr\right)}^{1/2}}$ (63)

Near the solid surface it can be assumed that the solution has the form

$\Psi \left(x,y,z\right)=\Psi \left(x,y,0\right){e}^{-\left(1-i\right)z/{\delta}_{ent}}$ . (64)

where this relation satisfies the following equation

$\frac{{\partial}^{2}\Psi (x,y,z)}{\partial z}=\frac{-i2}{{\delta}^{2}}\Psi (x,y,z)$ . Additionally, as in the aforementioned

cylindrical case, the temperature fluctuation at the wall is zero. The temperature conditions in Equations (58) (29b) (60) produces the following

${{T}^{\prime}}_{ac}+{T}_{ent}=0$ (65)

Applying Equation (64) into Equations (58), (59a), (59b) and splitting the “del” operator into tangential and normal components (where the tangential portion will have a subscript T) yields the following

${\nabla}_{T}\cdot {v}_{vor,T}\left(x,y\right)-\left(1-i\right){v}_{vor}\left(x,y\right)\cdot \frac{n}{{\delta}_{vor}}=0$ (66)

${v}_{ent,T}\left(x,y\right)=\frac{\beta \kappa}{\rho {c}_{P}}{\nabla}_{T}{T}_{ent}=0$ (67)

${v}_{ent}\cdot n=-\frac{\beta \kappa}{\rho {c}_{P}}\left(1-i\right)\frac{{T}_{ent}}{{\delta}_{ent}}$ (68)

At z = 0 plane Equation (61) will give the following boundary condition

${v}_{wall}={v}_{vor}+{v}_{ent}+{v}_{ac}$ (69)

and applying the ${\nabla}_{T}\cdot $ (horizontal divergence) along with Equation (67), Equation (69) can be rewritten the following way

${\nabla}_{T}\cdot {v}_{ac,T}\left(x,y\right)+\left(1-i\right){v}_{vor}\left(x,y\right)\cdot \frac{n}{{\delta}_{vor}}=0$ (70)

Remember that

$\gamma -1=\frac{T{\stackrel{\u02dc}{\beta}}^{2}\gamma {K}_{T}}{\rho {C}_{p}}=\frac{T{\stackrel{\u02dc}{\beta}}^{2}{K}_{S}}{\rho {C}_{p}}=\frac{T{\stackrel{\u02dc}{\beta}}^{2}{c}^{2}}{{C}_{p}}$ ${K}_{T}/{K}_{S}$

= isothermal/adiabatic bulk modul (71)

Additionally, using Equation (68) and (70) the subsequent conditions are obtained in the direction normal to the surface

${v}_{wall}\cdot n=-(1+i)\frac{{\delta}_{vor}}{2}{\nabla}_{T}\cdot {v}_{ac,T}\left(x,y\right)+(1-i)(\gamma -1)\frac{\omega}{{c}^{2}}\frac{{\delta}_{ent}}{2}\frac{{p}_{ac}}{\rho}+{v}_{ac}\cdot n$ (72)

Recall that the main goal is to obtain an expression for the energy that is being dissipated and this is dictated by the net power flowing out of the acoustic (ac) mode close to the boundary. Multiplying Equation (72) by the complex conjugate of the acoustic mode pressure, ${p}_{ac}^{*}$ , and solving for $-{p}_{ac}^{*}{v}_{ac}\cdot n$ , which is the intensity leaving the acoustic mode, provides the following equation

$\begin{array}{l}-{\left({I}_{ac}\cdot n\right)}_{av}\\ =\frac{1}{2}Re\left({p}_{ac}^{*}{v}_{wall}\cdot n\right)-\frac{{\delta}_{vor}}{4}\nabla \cdot \left\{Re\left[(1+i){p}_{ac}^{*}{v}_{ac,T}\left(x,y\right)\right]\right\}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{{\delta}_{vor}}{4}Re\left[(1+i){\nabla}_{T}{p}_{ac}^{*}\cdot {v}_{ac,T}\left(x,y\right)\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+Re\left[\frac{1}{2}(1-i)(\gamma -1)\frac{\omega}{{c}^{2}}\frac{{\delta}_{ent}}{2}\frac{{\left|{p}_{ac}\right|}^{2}}{\rho}\right]\end{array}$ (73)

Recall that the relative velocity at the wall is equal to zero. From the first relation in Equation (60) it can be deduced that

${\nabla}_{T}{p}_{ac}^{*}=i\omega \rho {v}_{ac,T}^{*}\left(x,y\right)$ (74)

The second term in the right hand side of Equation (73) will average out to zero given enough surface area. Equation (74) will be applied in the third term of Equation (73) to finally obtain the following relation

${\left(\frac{{d}^{2}E}{dAdt}\right)}_{diss}={\left(\frac{\omega \rho \mu}{2}\right)}^{1/2}{\left({\text{v}}_{ac,T}^{2}\left(x,y\right)\right)}_{av}+(\gamma -1){\left(\frac{\omega \rho \kappa}{2{c}_{P}}\right)}^{1/2}\frac{{\left|{p}_{ac}\right|}_{av}^{2}}{\rho {c}^{2}}$ (75)

For the mechanical dissipation the first term in the right hand side is used by [2]. The second term will be utilized in this study to treat the thermal dissipation and derive the acoustic bulk modulus.

The thermal counterpart of dynamic permeability and tortuosity are

$\varphi {T}_{o}=-\frac{{k}^{\prime}(\omega )}{\kappa}\frac{\partial {P}_{o}}{\partial t}$ (76)

${\alpha}^{\prime}(\omega ){c}_{P}{\rho}_{f}\frac{\partial {T}_{o}}{\partial t}=\frac{\partial {P}_{o}}{\partial t}$ (77)

where ${T}_{o}$ and ${P}_{o}$ satisfy the following

(78)

The real constant acts like a thermal form of a static thermal permeability, i.e. the thermal counterpart to Darcy’s static equation of.

A volume of porous sample enclosed by a thin membrane is considered. In this particular case, it is assumed that the thin membrane acts like an adiabatic wall. The interest is to know where the singularities to and lie. The existence of a singularity means that there exists a nontrivial solution to the above microscopic equation while the macroscopic input is zero, i.e.. This is equivalent to saying that there is no mechanical power entering the element. Therefore Equation (78) becomes

(79)

To guarantee that no thermal power is entering the element, the following has to hold

(80)

Since the membrane is adiabatic and since conditions listed in Equation (78) are being satisfied the integral in Equation (80) is consequentially fulfilled. Pre multiplying the Fourier equation of heat diffusion equation, Equation (79), by the complex conjugate, assuming a solution of the form, and then perform a volume integral in terms of the microscopic variable yields

(81)

Making use of the usual identification

(82)

(83)

So

(84)

The first term in the right hand side of Equation (84) goes to zero because of boundary condition found in Equation (80). Inserting this result in Equation (81) and solving for the following is obtained

(85)

Due to the shifting rule, this indicates that the singularity for lies in the negative imaginary axis. In the meantime consider Equation (77) where the boundary condition is that of. Additionally, no mechanical power is entering through the adiabatic wall. The aim is to find where the singularities of lie. Equation (78) is multiplied by. A volume integral results in

(86)

In [1] it is shown that is independent of the micro scale variable and therefore it is considered as an external source of power for this thermal element. Applying Equation (82), (83), and the Green’s theorem to the first term in Equation (86), the same result of Equation (85) is obtained; again, the singularities of also lie in the negative imaginary axis.

Following the same steps performed in [2] in deriving the relation for viscous length, the last term in the right hand side of Equation (73)

is studied. It is understood that this term repre-

sents the thermal portion of the energy flux exiting the main acoustic mode. The third term in the RHS of Equation (50a) is the thermal energy. Multiplying this term by the speed of sound, c, the energy flux entering the pore is attained. The ensuing equation represents the ratio of thermal energy that is leaving the acoustic mode to that of thermal energy that is entering the thermal element.

(87)

Recall temperature condition and from Equation (59b),

. Additionally from Equation (60),. Lastly,

applying Equation (71), Equation (87) simplifies to

(88)

In Equation (88),. The integration takes place at the micro level and it has been established that is independent of the micro scale variable, y, and thus the first equality sign in Equation (88). is the Thermal Length and is

defined as. Using Equation (57) as a pattern solution, the general solution

for average temperature is defined as (assuming ideal gas conditions)

(89)

In Equation (89), as increases the thinner the entropy boundary layer, , is and therefore the limit as the acoustic mode condition (ideal gas)

applies for the whole cross-section. This equation for is

plugged into Equation (3). The following expression for thermal permeability

is obtained

(90)

The second term in Equation (90) is obtained by applying the definition of. This term can asymptotically (as)be represented as

(91)

where, is the thermal permeability as. A combined asymptotic representation that interpolates thermal permeability for the whole frequency range is represented as

(92)

Finally, this is the term that is used in Equation (2) to obtain the compressibility factor of the fluid of the porous medium.

3. Conclusions and Closing Remarks: Absorption

Consider the acoustic medium near the wall of a tube (the x axis is taken to be along the tube). Set as the angle of incident and recall that the solution to the wave equation tangent to the tube wall will have the form of

therefore. Also recognize that

. With these definitions, the term in Equation (73) can be represented as

(93)

An expression for the wall impedance can be given as

(94)

The objective is to obtain a simple partial differential equation that incorporates the impedance effects of the viscous and thermal boundary layer. Following techniques introduced by Crandall et al. [15], The variational form of the Helmholtz equation is

(95)

The is a segment of the perimeter of the tube and

Applying the definition of impedance, the following is obtained

(96)

and applying this to the variational Helmholtz Equation yields

(97)

As a first order approximation, it is assumed that the parameters do not vary with respect to the cross-sectional area

(a) (b)

Figure 1. Left: TL comparison (dB); Right: TL difference(dB) measured data is the reference.

Figure 2. Multi layer configuration.

perimeter length (98)

Performing the variation generates

(99)

For the case pertaining to constant cross-sectional area, the dispersion relationship is

(100)

is the absorption

coefficient of the wall. It is interesting to note the similarity in form between the imaginary part of Equations (100) and (88).

Equation (51g) pertaining to the absorption due to the relaxation process, , combined with the effect of the thermal permeability, , on the bulk modulus, , of the porous medium are incorporated into Transmission Matrix Method with finite size correction via Green’s functions techniques (FTMM) [16]. The parameters to this study are listed in Table 1 in [1]. Figure 1 shows the comparison between experimental data and the FTMM calculations for the multi-layer configuration shown in Figure 2. Further analysis, especially a parameter study on the different relaxation and thermal components is currently being done. This study mainly assumes air within the porous medium but there is current interest in applying these to pores with different rheological make up.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Teagle-Hernandez, A., Ohtmer, O. and Nguyen, D. (2019) Model for Frequency Dependence of Thermal Permeability in Order to Quantify the Effects of Thermal Exchange on Wave Propagation in Multi Layered Porous Medium. Journal of Applied Mathematics and Physics, 7, 1440-1462. https://doi.org/10.4236/jamp.2019.77097

References

- 1. Teagle-Hernandez, A., Ohtmer, O. and Nguyen, D. (2018) Mathematical Modeling of Porous Medium for Sound Absorption Simulations: Application of Multi-Scales and Homogenization. Journal of Applied Mathematics and Physics, 6, 2705-2717. https://doi.org/10.4236/jamp.2018.612224
- 2. Johnson, D. and Koplik, J. (1987) Theory of Dynamic Permeability and Tortuosity in Fluid-Saturated Porous Media. Journal of Fluid Mechanics, 176, 379-402. https://doi.org/10.1017/S0022112087000727
- 3. Thomasson, S.I. (1980) On the Absorption Coefficient. Acustica, 44, 265-273.
- 4. Brouard, B., Lafarge, D. and Allard, J.F. (1995) A General Method of Modeling Sound Propagation in Layered Media. J. Sound Vib., 183, 129-142. https://doi.org/10.1006/jsvi.1995.0243
- 5. Folds, D. and Loggins, C.D. (1977) Transmission and Reflection of Ultrasonic Waves in Layered Media. J. Acoust. Soc. Amer., 62, 1102-1109. https://doi.org/10.1121/1.381643
- 6. Feynman, R. and Hibbs, A. (1965) Quantum Mechanics and Integrals, McGraw-Hill.
- 7. ter Haar, D. (1954) Elements of Statistical Mechanics. Rinehart, New York. https://doi.org/10.1119/1.1933869
- 8. Schiff, L.I. (1955) Quantum Mechanics. McGraw-Hill, New York.
- 9. Landau, L. and Lifshitz, E. (1959) Fluid Mechanics. Pergamon.
- 10. Bellomo, N. and Gatignol, R. Lecture Notes on the Discretization of the Boltzmann Equation. World Scientific.
- 11. Herzfield, K.F. and Rice, F.O. (1928) Dispersion and Absorption of High Frequency Sound Waves. Phys. Rev., 31, 691-695. https://doi.org/10.1103/PhysRev.31.691
- 12. Lighthill, M.J. (1956) Viscosity Effects n Sound Waves of Finite Amplitude, Batchelor and Davies (Eds.) Surveys in Mechanics. Cambridge University Press, London.
- 13. Trilling, L. (1955) On Thermally Induced Sound Fields. J. Acoust. Soc. Am, 27, 425-431. https://doi.org/10.1121/1.1907920
- 14. Truesdell, C. (1953) Precise Theory of the Absorption and Dispersion of Forced Plane Infinitesimal Waves According to the Navier-Stokes Equations. J. Ration. Mech. Anal., 2, 643-730. https://doi.org/10.1512/iumj.1953.2.52032
- 15. Crandall, S.H., Karnopp, D.C., Kurtz Jr., E.F. and Pridmore-Brown, D.C. (1968) Dynamics of Mechanical and Electromechanical Systems. McGraw-Hill, New York.
- 16. Nguyen, D. and Teagle-Hernadez, A. (2019) Application of Transfer Matrix Method, Green’s Functions, and Variational Technique to Predict Diffuse Absorption Coefficient. CalTech-SOCAMS, Pasadena.

Appendix A: Linear Acoustic Equation

(A.1)

(A.2)

(A3)

(A4)

(A.5)

Nomenclature

= Fluid Pressure, Macro Scale

= Porosity

= Viscosity

= Domain Occupied by Fluid

C_{p} = specific heat at constant pressure C_{v} = Specific heat at constant volume

= Coefficient of thermal Conductivity

= Prandlt Number =

, = density of the Fluid

T = Temperature deviation

= Vibrational temperature pertaining to molecular type v

= theorist temperature, kT, k = Boltzmann constant

= Dynamic Viscous Permeability

= Dynamic Thermal Permeability

= asterisk superscript means conjugate

= special tilde 2^{nd} rank tensor

= Entropy, = Entropy without molecular vibration,

= molecular vibration entropy

= Thermal expansion coefficient =