8.11. Poroelastic Acoustics

The mixed displacement-pressure finite element model, based on Biot's theory, solves the poroelastic problem in acoustics (Allard [440]). Sound absorption occurs in the poroelastic materials.

8.11.1. Mixed Displacement and Pressure Formulation for Poroelastic Acoustic Material

In Biot's poroelasticity theory, the frequency-dependent effective densities are defined as:

(8–261)

where:

= porosity
= fluid mass density (air)
= tortuosity
= the bulk density of the solid phase of the poroelastic material
= solid mass density
= fluid resistivity
= fluid dynamic (shear) viscosity
= viscous characteristic length
ω = angular frequency

In the poroelastic acoustic material, the mixed displacement and pressure (u-P) governing equations are as follows:

(8–262)

(8–263)

where:

= solid stress tensor (eliminating the dependency of fluid displacement) defined as:
= displacement of the solid (skeleton)
= strain tensor of the solid
N = complex shear modulus in vacuo (second Lame coefficient) of the solid
= effective density defined as:
p = pressure in the pores
Ks = bulk modulus of the elastic solid (frame)
Kb = bulk modulus of the porous material in vacuo
Kf = bulk modulus of the air defined as:
P0 = static reference pressure
= specific heat ratio of air
Prt = Prandtl number
= thermal characteristic length

For the majority of poroelastic materials used in acoustics, . The Biot's constant is defined by:

(8–264)

In most cases, the Biot's constant equals 1.

8.11.2. The Finite Element Model

Introducing variable transformation,

(8–265)

the finite element formulations are obtained by testing Equation 8–262 and Equation 8–263 using the Galerkin procedure. Equation 8–262 and Equation 8–263 are multiplied by testing function and q, respectively, and integrated over the volume of the domain with some manipulation to yield the following:

(8–266)

(8–267)

where:

= total stress tensor, defined by:
= fluid phase displacement

The fluid phase displacement is cast by:

(8–268)

In poroelastic material, the total displacement is written as:

(8–269)

8.11.3. Boundary Conditions of Poroelastic Acoustics

Boundary conditions are applied to the finite element model to solve poroelastic acoustic equations according to the application.

Natural Boundary Condition

(8–270)

where subscript n = normal direction of the interface.

In a finite element model, the boundary integral terms are not considered with the natural boundary condition.

Fixed Edge

(8–271)

A poroelastic medium is bounded onto a rigid impervious wall. The solid phase displacement and normal relative displacement vanish at the rigid wall.

Free Porous Surface

(8–272)

On the free porous surface, the pressure is set to zero.

Sliding Edge

The sliding edge condition indicates that the normal displacements on an impervious rigid wall are constrained, and the tangential displacements are free.

(8–273)

The boundary integrations vanish because the normal solid displacement is constrained.

Pervious Porous Surface with Permeability K

The boundary condition imposed with the difference in displacements between solid phase, fluid phase, and pressure on the pervious porous surface is written as:

(8–274)

Imposed Pressure

(8–275)

For the case of imposed pressure, the pressure is set to p0.

Imposed Displacement Field

(8–276)

The continuity between the imposed displacement and the solid phase displacement is maintained. The normal displacements between solid phase and fluid phase continues.

8.11.4. Coupling Conditions of Poroelastic Acoustics

The couplings occur at the interfaces between the poroelastic material and other materials.

Poroelastic-Elastic Coupling Conditions

The coupling conditions at the interface between poroelastic and elastic media are given by:

  • Displacement continuity

  • Impervious boundary

  • Stress continuity

where:

= elastic stress tensor

The coupling between poroelastic and elastic media is natural.

Poroelastic-Acoustic Coupling Conditions

The coupling conditions at the interface between poroelastic and acoustic media are given by:

  • Normal stress continuity

  • Total displacement continuity

  • Pressure continuity

where:

pa = pressure in acoustic media
ρ0 = mass density of acoustic media