Mathematical Modeling of Porous Medium for Sound Absorption Simulations: Application of Multi-Scales and Homogenization ()
1. Introduction
NV development begins by assigning NV targets for different systems. Using an automobile development cycle as an example, this translates to assigning NV targets to systems like dashboard, floor, roof, trunk, and other systems of the automobile. These targets dictate TL, IL and Alpha performances that are projected to different parts that the system comprises of. These performance levels are usually derived via a Hybrid Statistical Energy Analysis (HSEA) technique [1] [2] [3]. The aggregate sound attenuation performance of each component is affected by the multi-layer porous materials that are utilized to manufacture it. NV specialists run simulations in order to decide the optimum combination.
Material suppliers furnish layer parameters such as porosity, tortuosity, resistivity, foam bulk modulus, skeleton Young’s modulus, viscous and thermal lengths that are then applied to calculate flat sample TL, IL, and Alpha. Quite often, OEM engineers plug in values pertaining to these parameters, but have little understanding of the energy absorbing or sound blocking mechanism. Decisions are based on past experiences. The difficulty is to find a work that encompasses the gamut of equations required and clearly explains the physics behind transmission and absorption of energy. This work builds a bridge between the gaps in order to obtain a more fundamental understanding. In this part of the series the goal is to derive the coupled fluid/structural equations. The macro-scale equations are obtained by applying the multi scales and homogenization techniques. The scales of the pores are small compared to the macroscales. The process will show how energy is being dissipated due to fluid/structure interaction; the dilatational/compressional interplay between the acoustic/fluid and structural medium.
The encapsulation of the acoustic medium due to how tortuous the foam material is and the mass entrapped in the viscous boundary layer are incorporated in the above equations. There is also a boundary layer where the flow is not adiabatic, 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 paper ties the aforementioned physical phenomena to the parameters of tortuosity, dynamic viscous and thermal permeability, viscous as well as thermal lengths. The derivation of the equation bands together ideas and techniques presented in [4] [5] [6] [7]. Part 2 of this series details how to incorporate the foam equations with other type of layers, e.g. panels, air gap, other foams, via the Transfer Matrix Method (TMM) and Finite Transfer Matrix Method (FTMM) in order to calculate TL, IL, and Alpha.
2. Formulation
2.1. Basic Tensor Calculus and Notation
This body of work takes advantage of tensor notation and tensor calculus. This notational generalization makes it easy to describe the physics regardless of the abstract methodology used to define the dynamics of the system. In order to clearly define ideas, the Einstein summation convention is used in which covariant index followed by the identical contravariant index is implicitly summed over. Variables written in bold are tensors of rank 1, i.e. vectors. In the meantime, a bold letter with a special tilde as shown here,
, is a tensor of rank 2. The covariant derivative of a contravariant vector,
is given as
(1)
Additionally, the use of the divergence of a second rank tensor,
, is defined as
(2)
In Equation (1) and (2),
, is the Chritoffel symbol of the second kind.
2.2. Fluid-Structure Interaction: Two-Scale Expansion
The pore composition of the porous medium has amicroscopic length l, while the macroscopic length is designated L. The material is statistically homogenous at the macroscopic scale. The ratio of the two length scale is set as
where
. The different acoustic and vibration variables will depend on two scales: a slow scale x and a fast scale
. The two-scale expansion is applied to the coupled fluid-structure equations. For relatively small viscosity values, conditions found in acoustics, the intrinsic viscous relaxation time is
, here
is
the kinematic viscosity
. In order to correctly compare orders of asymptotic
expansion the dimensionless scaled terms of the equations of motion are derived and hence the right scaling factors are required. Due to the dynamic viscous permeability conditions, acoustic pressure and structural stresses are scaled with
respect to
. Concurrently, due to the dynamic thermal permeability conditions the temperature fluctuation is scaled with respect to
.
The dimensionless form of the coupled Navier-Stokes and structural equations along with the fluid/structure boundary conditions are given as
(3)
(4)
(5)
(6)
(7)
(8)
(9)
The dimensionless/scaled version of each variable is denoted with an asterisk. For each of the variables in Equations (3)-(9) are represented by a two scale expansion
. The covariant derivative becomes
. Applying the expansion and collecting all terms of O(1/) the following are obtained
(10)
(11)
(12)
(13)
Collecting all O(1) terms the following key equations are obtained
(14)
(15)
(16)
(17)
(18)
To better explain the interaction conditions the relative displacement
is introduced,
(19)
and is incorporated into Equations (14) and (15)
(20)
(21)
(22)
Thus, the new boundary condition becomes
at
In the meantime, to solve the coupled partial differential Equations (20)-(22),it is assumed that there exist a rank 1 operator,
, a second rank tensor
and a scalar
that act as transfer functions or operators that satisfy the following
(23)
(24)
(25)
In [8] [9], a more detailed analysis at the microscopic level is constructed. There, it is shown that the microscopic dynamics produces exponential decaying effects on the transfer functions
and
. Fluid equations pertaining to the viscous boundary layer relates this decaying effects on viscous length. For the thermal counterpart, the decaying effects due to the entropy boundary layer are related to thermal length. Homogeneity is obtained by averaging at the microscopic level, i.e. with respect to y, thus the following relation is obtained
(26)
Comparing this equation and the definition of dynamic permeability will determine that
Therefore
(27)
Taking the inverse of
, Equation (27) can be rearranged as
(28)
is known as Dynamic tortuosity. Variables with an overbar, e.g.
, are ones that have been averaged at the microscopic level.
The dimensional form of Equation (25) is
(29)
By taking the microscopic average of the above equation one will obtain
such that
(30)
Replacing the temperature with this relationship to pressure in Equation (21), the following is obtained
(31)
The term
is the compressibility factor in dimensionless form. Set
(bulk modulus for the acoustic medium) and apply this in the last equation above
(31')
After taking the microscopic average of (31), the following expression is obtained
(32)
2.3. Structural: Finite Element ? Representative Volume Element
The structural skeleton of the porous medium will be split up into identical unit cells where the macroscopic variables will represent conditions at the boundary, while microscopic variables will be tied to conditions interior of the unit cell. The boundary conditions and structural characteristics of the unit cell will be assumed statistically periodic. Special boundary conditions and interior forces are applied in order to obtain a set of Basis functions that will expand the interior solutions. A set of equations that will relate boundary stresses to interior strains will be easily obtained by representing
(33)
In order to obtain a relationship between macroscopic strain and macroscopic stresses and considering the impedance from the micro-scale portion the following finite element form is used
(34)
are the elements pertaining to interior nodes,
and
are the coupling elements, while
are the elements pertaining to the boundary nodes. The assumption that there exist a potential energy for the fluid-structure problem forces
and
to be transpose operators from each other. From the top part of the above matrix formulation an expression for
as a function of
can be obtained
(35)
thus
(36)
The interior basis functions can be obtained by solving the following matrix formulation
(37)
where
is a special tensor, of rank 2, that transforms the scalar,
, pressure into its tensor form and applies it normal to the structural surface. It should be clear that
,
,
, and
are 4th order contravariant tensors. The stress component that is transferred to the boundary nodes will be equal to
(38)
By combining the two concepts the equations for microscopic strains and macroscopic stress is
(39)
(40)
From the acoustic relations, Equation (20) can be placed into the following form
(41)
Taking the microscopic average of this equation and utilizing Green’s theorem Equation (41) becomes
(42)
Similarly, the same procedure is done to the structural Equation, (17)
(43)
Due to the boundary conditions at the interface, Equation (7), the foregoing equations becomes
(44)
where
.
2.4. Coupled Equations
Equation (28) will be rewritten as
(45)
Reintroducing dimensions to the compressibility factor
one will get
(46)
The final dimensional form pertaining to the acoustic is
(47)
Plugging in Equation (39) into (47), the dimensional and averaged structural/acoustical expression is
(48)
(49)
(50)
where
(51)
The grouping of Equations (40), (45), (44) and (50) develops one form of the coupled equation
(52)
(53)
(54)
(55)
The term
and the expression for po will be subtracted from Equation (54)
(56)
Equation (56) in its current form can model anisotropic material. A simpler expression can be obtained by assuming that the material is statistically isotropic. Therefore the first term in the right hand side can be expressed in the following form after taking advantage of basic concepts of symmetric tensors
(57)
Additionally, in [6] it is shown that
(58)
Set
and because of (58).
. Analyzing how Equation (58) was derived,
represents the proportion of fluid pressure that produces the same strains as the total stress.
Plugging Equation (57) and the definition for
into Equation (56) the subsequent equation is obtained
(59)
and
are new Lame’ constants of the elastic portion when the porous material is drained where as
is the contained Lame’ constant sometimes referred as
in the literature. Equation (59) is in the exact form as Equation (30) found in Biot and Willis’ paper [10]. Equation (55) can also be expressed, using the definition for
, as
(60)
It is desired to define the equations above as a function of uo and Uo. Hence, the relative displacement can be expressed as
. Equation (53) becomes
(61)
(62)
Now, subtracting (61) from (52) one obtains
(63)
Applying the same definition of
into Equation (53) the acoustic portion of the coupled equation of motion turns into
(64)
Set
(65)
(66)
(67)
so that (63) and (64) can be expressed in a simpler form
(68)
(69)
Plug into (68) and (69) the definitions for
and
found in Equations (54) and (55), along with the replacements for
and
one arrives to the following equation
(70a)
(70b)
and
(71a)
(71b)
Equation 70(a) and Equation 70(b) and Equation 71(a) and Equation 71(b) are the Equation which will be utilized to study wave propagation through porous medium. Note, that the form of the equations are in the same format as Biot’s Equation [11]
3. Results and Conclusion
A Simple application of Equations (70) and (71) is applied to the simple layer configuration shown in Figure 1. The figure shows a multi layer system where there is a plate that is glued to foam 1. The parameters pertaining to plate1 and foam1 are listed in Table 1. This simulation also applied the following parameters that are not listed in the table: Viscous length = 40 μm, Thermal length = 80 m. The goal of this simulation is to calculate Random Incidence Transmission Loss (TL). This is achieved by applying Finite Size correction Transfer Matrix Method (FTMM). The Finite Size Correction is achieved by incorporating a Green’s Function integration to the radiation efficiency. Equations (70) and (71) are used to simulate how the structural and fluid stresses will change as the acoustic wave travels through the panel. The results are shown in Figure 2. The results are compared to measured results and also to the empirical formulas of Delaney and Bazley [12] [13]. The graph shows that the calculated results come within 0.4 dB.
Table 1. Parameters pertaining to plate and foam.
Figure 2. Top: TL comparison (dB); Bottom: TL difference c(dB) measured data is the reference.
The current method used in order to calculate TL requires the use of various measured variables. Parameters like tortuosity, viscous/thermal lengths, tortuosity, flow resistivity are all measured. Additionally, these measurements are being done at different times at different facilities that have disparate techniques and equipment. This makes running simulations challenging. One of the goals of going through this mathematical derivation is to find systematic techniques to reduce the disparity in these parameters or at least get a better understanding of the physics.
Nomenclature
= trace of a tensor
= averaged at the microscopic level
= Fluid Velocity
= Structural displacement
= Fluid displacement
= Fluid Pressure
= Porosity
= Viscosity
= Second Viscosity
= Domain Occupied by Fluid
= Domain Occupied by Structure
= Boundary of Domain Occupied by Fluid
= Boundary of Domain Occupied by Structure
= Cp/Cv Cp = specific heat at constant pressure Cv = Specific heat at constant volume
= Coefficient of thermal Conductivity
= Prandlt Number =
= unit normal pointing into the solid
= Stress Tensor in the Fluid
= Stress Tensor in the structure
= Elastic fourth ranked Contravariant tensor, operates on
= density of the Fluid
= density of the Structure
= Lame Parameters
= Coefficient of thermal Conductivity
= Temperature deviation
= Dynamic Viscous Permeability
= Dynamic Thermal Permeability
= Acoustic bulk Modulus
= asterisk superscript means scaled variable (dimensionless)