**Transient Acoustic Wave Propagation in Porous Media**

Zine El Abiddine Fellah, Mohamed Fellah and Claude Depollier

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/55048

## **1. Introduction**

126 Modeling and Measurement Methods for Acoustic Waves and for Acoustic Microdevices

(3):482-488

New York.

8345-8350

118:1157-1159.

Acoustics, vol. 22: 123-131.

of surface waves, J. Appl. Phys., vol. 49, No. 1.

[49] Pustelny T (1997) Surface Acoustic Wave Techniques for the Real-Surfaces Investigations of InP (111), InAs (111) Semiconductors, Acustica-acta acoustics, vol. 83

[50] Matthews H (1997) Surface wave filters Design, construction and use, J Wiley and Sons,

[51] D'Amico A, Palma A, Verona E (1982), Hydrogen sensor using a palladium coated

[52] Jakubik W, Urbańczyk M, Kochowski S, Nadolski M (2001) Surface acoustic wave sensor for hydrogen detection with a two-layered structure, Molecular and Quantum

[53] Jakubik W (2007), Investigation of thin film structure of WO3 and WO3 with Pd for hydrogen detection in a surface acoustic wave sensor system, Thin Solid Films, 515:

[54] Tiersten H, Sinha B.K (1978) A perturbation analysis of the attenuation and dispersion

[57] Hejczyk T, Urbańczyk M, Jakubik W (2010), Analytical Model of Semiconductor Sensor

[58] Hejczyk T, Urbańczyk M, Jakubik W (2010), Numerical Results of Modeling Semiconductor Sensor Layers in SAW Gas Sensors, Acta Physica Polonica A, vol

surface acoustic wave delay-line, IEEE Ultrasonics Symposium: 308-311.

[55] White R (1970), Surface elastic waves, Proc. IEEE, vol.58, No8:1238-1277.

[56] Williams D E, Hilger A, (1987), Solid state gas sensors, Bristol and Philadelphia.

Layers in SAW Gas Sensors, Acta Physica Polonica A, vol 118:1148-1152.

A porous material is a two-phase medium consisting of a solid part (skeleton) and a fluid part (pore space). During the propagation of a sound wave in such a medium, interactions between these two phases of different nature take place, giving various physical properties that are unusual in classical media. The large contact area between solid and fluid, which is the main characteristic of porous media induces new phenomena of diffusion and transport in the fluid, in relation to micro-geometry of the pore space. Many applications are concerned with understanding the behavior of acoustic waves in such media. In geophysics, we are interested in the propagation of acoustic waves in porous rocks, for information on soil composition and their fluid content. Oil companies have greatly contributed to the study of acoustic properties of natural porous media. In medicine, the characterization of porous media such as trabecular bone, is useful for diagnosing osteoporosis, bone disease that is manifested by the deterioration of bone microarchitecture. Acoustic characterization of materials is often achieved by measuring the attenuation coefficient and phase velocity in the frequency domain [1] or by solving the direct and inverse problems directly in the time domain [2–4]. The attractive feature of a time domain based approach is that the analysis is naturally limited by the finite duration of ultrasonic pressures and is consequently the most appropriate approach for the transient signal. The objective of this chapter is to show the most recent theoretical and experimental methods developed by the authors for the acoustic characterization of porous materials. The direct and inverse scattering problems are solved in time domain using experimental reflected and transmitted signals. The physical parameters of the porous medium are recovered by solving the inverse problem at the asymptotic domain corresponding to the high frequency range (ultrasound), and at the viscous domain (low frequency range). Figures 1 and 2 give an example of porous materials commonly used in the characterization. The figure 1 shows a sample of air-saturated plastic foam used for sound absorption, and figure 2 shows a sample of cancellous bone used in the diagnosis of the disease of osteoporosis.

©2013 Fellah et al., licensee InTech. This is an open access chapter distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0),which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. ©2013 Fellah et al., licensee InTech. This is a paper distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. © 2013 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

2 Will-be-set-by-IN-TECH 128 Modeling and Measurement Methods for Acoustic Waves and for Acoustic Microdevices Transient Acoustic Wave Propagation in Porous media <sup>3</sup>

finite set of parameters low and high frequencies may be sufficient to describe the dynamic susceptibilities. For a given frequency, one can define a boundary layer thickness, respectively viscous and thermal *δ* = (2*η*/*ωρ <sup>f</sup>*)1/2, *δ* = (2*κ*/*ρ <sup>f</sup> Cpω*)1/2, where *η* is the viscosity, *ρ <sup>f</sup>* the fluid density, *κ* the coefficient of thermal conduction, *Cp* the specific heat at constant pressure and *ω* the angular frequency. The quantities *δ* and *δ* are respectively the penetration depths of isothermal viscous and rotational movement. The boundary layer is small if the penetration depth is small compared to a "radius" characteristic pore. It is then at high frequencies. The opposite situation corresponds to low frequencies. In the high frequencies, inertial effects are described by the tortuosity *α*∞, the visco-thermal losses by viscous and thermal lengths Λ et

, respectively. At low frequencies, inertial effects are represented by the viscous tortuosity

There are several geometrical parameters that can be associated with a porous structure. Some depend only on the shape of the porous solid and do not vary when it expands or changes the scale of the solid. This is the case of the porosity and tortuosity. Others depend on the dimensions, in the case of viscous and thermal permeabilities that vary in shape as an equal

• **Porosity** *φ*: The porosity of a porous material saturated with a fluid is defined by the ratio of the volume of fluid saturating the porous medium to the total volume of the sample. This quantity is dimensionless and can vary from 0 (non-porous solid material) to 1 (free

• **Tortuosity** *α*∞: The tortuosity is a geometric parameter describing the appearance of tortuous pores, and inertial coupling between the fluid and the structure of the material in high-frequency regime. Tortuosity takes its lowest value (*α*<sup>∞</sup> = 1) in the case of porous materials having a straight pores, and large values (*α*∞ = 1.5, 2) for the most resistive

• **Inertial factor** *α*<sup>0</sup> **(low frequency tortuosity)**: The inertial factor *α*<sup>0</sup> corresponds to the low frequency approximation of the dynamic tortuosity given by Norris [10]; *<sup>α</sup>*<sup>0</sup> <sup>=</sup> *<sup>&</sup>lt;***v**(**r**)**<sup>2</sup>***<sup>&</sup>gt;*

• **The specific flow resistivity** *σ*: When a porous material placed in a tube sealed to the walls, is traversed by a fluid, it is a pressure difference between the two free faces. This means that a proportional relationship exists between flow rate *Qv* and pressure drop Δ*P* given by the Δ*P* = *RQv*, the coefficient R is the resistance to fluid flow and can therefore be defined as the ratio Δ*P*/*Qv*. For a homogeneous porous material, the resistance is proportional to

where *<* **v**(**r**) *>* is the average velocity of the viscous fluid for direct current flow within an elementary volume, small compared to the relevant wavelength, but large compared to the individual grains/pores of the solid. Note that the tortuosity *α*∞ (high frequency parameter) has a similar definition than *α*0, the difference is that for *α*∞, **v**(**r**) corresponds

<sup>0</sup>: This low frequency parameter is the thermal counterpart of *α*0. Its

*<sup>&</sup>lt;θ*0*>*<sup>2</sup> , where *<sup>θ</sup>*<sup>0</sup> is the scaled response excess temperature field

<sup>0</sup>. Porosity *φ* is a

Transient Acoustic Wave Propagation in Porous Media 129

*<sup>&</sup>lt;***v**(**r**)*>***<sup>2</sup>** ,

*α*0, and visco-thermal exchanges by the visco-thermal permeabilities *k*<sup>0</sup> and *k*

surface, and viscous and thermal characteristic lengths that vary as lengths.

parameter playing an important role both in high in low frequency.

**2.2. Parameters describing the porous materials**

to the velocity of the perfect incompressible fluid.

<sup>0</sup> <sup>=</sup> *<sup>&</sup>lt;θ*<sup>2</sup> <sup>0</sup>*>*

Λ

fluid).

material.

• **Thermal tortuosity** *α*

definition is given by: *α*

in the low frequency range.

**Figure 1.** Air saturated Plastic foam

**Figure 2.** Human cancellous bone sample

## **2. Porous materials in rigid structure: Equivalent fluid model**

## **2.1. Introduction**

The air saturated porous materials are often used as a sound absorber in the fight against noise pollution. The polyurethane foam, felt or glass wools are three examples of common materials used in aerospace and automotive industries and in building construction. Generally, when a sound wave propagates through a porous material saturated with air, the structure remains stationary and non-deformable with respect to the acoustic excitation, this is due to the heaviness and stiffness of the skeleton of the structure relative in air. This is called porous rigid structure and, in this case, we use the equivalent fluid model where viscous and thermal effects are disconnected [5]. The fluid-structure interactions responsible for the sound attenuation (particularly important in porous media), are described by two dynamic susceptibilities, the dynamic tortuosity and compressibility. The dynamic tortuosity [6] decribes the viscous and inertial effects, the dynamic compressibility [7–9] represents the thermal effects. A prediction of the acoustic behavior of porous material requires the determination of both dynamic susceptibilities. The latter can depend only on the physical characteristics of the fluid and the geometry of the fluid domain. Sound propagation in the saturating fluid involves both transport and diffusion characteristics. The phenomenon of transport is the fact that a global movement fluid/solid is induced during the propagation. The phenomenon of diffusion is related to the expansion/compression of the fluid. They result in a change in temperature fluid/solid (the latter remaining at room temperature) driving a heat diffusion, not negligible if the fluid is a gas. Moreover, in the case of sound propagation, the frequency of movement is important. It is clear that the character of the movement is changed to different frequencies, the fluid probe different aspects of micro geometry. In principle one should expect that there are an infinite number of geometric parameters. But on the other hand, a practical viewpoint, the information provided by a finite set of parameters low and high frequencies may be sufficient to describe the dynamic susceptibilities. For a given frequency, one can define a boundary layer thickness, respectively viscous and thermal *δ* = (2*η*/*ωρ <sup>f</sup>*)1/2, *δ* = (2*κ*/*ρ <sup>f</sup> Cpω*)1/2, where *η* is the viscosity, *ρ <sup>f</sup>* the fluid density, *κ* the coefficient of thermal conduction, *Cp* the specific heat at constant pressure and *ω* the angular frequency. The quantities *δ* and *δ* are respectively the penetration depths of isothermal viscous and rotational movement. The boundary layer is small if the penetration depth is small compared to a "radius" characteristic pore. It is then at high frequencies. The opposite situation corresponds to low frequencies. In the high frequencies, inertial effects are described by the tortuosity *α*∞, the visco-thermal losses by viscous and thermal lengths Λ et Λ , respectively. At low frequencies, inertial effects are represented by the viscous tortuosity *α*0, and visco-thermal exchanges by the visco-thermal permeabilities *k*<sup>0</sup> and *k* <sup>0</sup>. Porosity *φ* is a parameter playing an important role both in high in low frequency.

## **2.2. Parameters describing the porous materials**

2 Will-be-set-by-IN-TECH

**Figure 1.** Air saturated Plastic foam

**Figure 2.** Human cancellous bone sample

**2.1. Introduction**

**2. Porous materials in rigid structure: Equivalent fluid model**

The air saturated porous materials are often used as a sound absorber in the fight against noise pollution. The polyurethane foam, felt or glass wools are three examples of common materials used in aerospace and automotive industries and in building construction. Generally, when a sound wave propagates through a porous material saturated with air, the structure remains stationary and non-deformable with respect to the acoustic excitation, this is due to the heaviness and stiffness of the skeleton of the structure relative in air. This is called porous rigid structure and, in this case, we use the equivalent fluid model where viscous and thermal effects are disconnected [5]. The fluid-structure interactions responsible for the sound attenuation (particularly important in porous media), are described by two dynamic susceptibilities, the dynamic tortuosity and compressibility. The dynamic tortuosity [6] decribes the viscous and inertial effects, the dynamic compressibility [7–9] represents the thermal effects. A prediction of the acoustic behavior of porous material requires the determination of both dynamic susceptibilities. The latter can depend only on the physical characteristics of the fluid and the geometry of the fluid domain. Sound propagation in the saturating fluid involves both transport and diffusion characteristics. The phenomenon of transport is the fact that a global movement fluid/solid is induced during the propagation. The phenomenon of diffusion is related to the expansion/compression of the fluid. They result in a change in temperature fluid/solid (the latter remaining at room temperature) driving a heat diffusion, not negligible if the fluid is a gas. Moreover, in the case of sound propagation, the frequency of movement is important. It is clear that the character of the movement is changed to different frequencies, the fluid probe different aspects of micro geometry. In principle one should expect that there are an infinite number of geometric parameters. But on the other hand, a practical viewpoint, the information provided by a

There are several geometrical parameters that can be associated with a porous structure. Some depend only on the shape of the porous solid and do not vary when it expands or changes the scale of the solid. This is the case of the porosity and tortuosity. Others depend on the dimensions, in the case of viscous and thermal permeabilities that vary in shape as an equal surface, and viscous and thermal characteristic lengths that vary as lengths.


#### 4 Will-be-set-by-IN-TECH 130 Modeling and Measurement Methods for Acoustic Waves and for Acoustic Microdevices Transient Acoustic Wave Propagation in Porous media <sup>5</sup>

the length *d* of the sample and inversely to his section *S*. Then can be defined, the resistivity as the specific resistance to the passage of a fluid. Denoted *σ*, the resistivity is defined by *σ* = (*S*/*d*)*R*.

where *j*

in the frequency domain is given by:

 <sup>1</sup> <sup>−</sup> <sup>1</sup> *jx*

*α*(*ω*) = *α*<sup>∞</sup>

to the tortuosity *α*∞.

*β*(*ω*) = *γ* − (*γ* − 1)/

the parameter *k*�

permeability.

where *k*�

(5):

<sup>2</sup> <sup>=</sup> <sup>−</sup>1, *<sup>η</sup>* is the dynamic viscosity of the fluid, *<sup>ω</sup>* the angular frequency, *<sup>ρ</sup> <sup>f</sup>* the

*<sup>ρ</sup> <sup>f</sup> <sup>k</sup>*<sup>0</sup> *<sup>j</sup><sup>ω</sup>* , *<sup>ω</sup>* −→ 0, (2)

*<sup>σ</sup>* (with the dimension of a surface), *σ*

Transient Acoustic Wave Propagation in Porous Media 131

and *<sup>M</sup>* <sup>=</sup> <sup>8</sup>*k*0*α*<sup>∞</sup>

<sup>0</sup>/*φ*, the quantities Λ� and *k*�

, when *ω* −→ 0. (5)

<sup>0</sup>*Pr ηφ*

and *<sup>M</sup>*� <sup>=</sup> <sup>8</sup>*k*�

0 *<sup>φ</sup>*Λ�<sup>2</sup> .

, *ω* −→ ∞, (4)

*<sup>φ</sup>*Λ<sup>2</sup> . (3)

<sup>0</sup>/*φ* for

density of the fluid. In the low frequency limit, the expression of the dynamic tortuosity

is the specific resistance passage of the fluid. The meaning of the quantities *α*∞, Λ, *k*<sup>0</sup> have been given previously. Johnson et al [6] imply that these three geometric parameters provide sufficient information about the microgeometry. The proposed general expression for *α*(*ω*) is

A model similar to the previous one has been set for the frequency dependence of the thermal exchanges between the fluid and the structure, because of the heat capacity of the solid part is

thermal effects playing respectively the role of Λ and *k*<sup>0</sup> for viscous effects. In a high-frequency

where *Pr* is the number of Prandt and *γ* the adiabatic constant. As before, this behavior implies that the thermal boundary layer thickness *δ*� becomes small compared to any characteristic radius of curvature of the interface. One can see that there is no equivalent of the tortuosity *α*<sup>∞</sup> for thermal effects, at very high frequency *β*(*ω*) tends to 1 as *α*(*ω*) tends

<sup>0</sup>*Pr*

analogous to the parameter *k*0, but adapted to the thermal problem. The proposed general expression for *β*(*ω*) is then that of the simplest analytic function satisfying the limits ((4) and

<sup>0</sup>, which has the same size (area) that the Darcy permeability of *k*0, is a parameter

where *<sup>x</sup>*� <sup>=</sup> *ωρ <sup>f</sup> <sup>k</sup>*�

<sup>0</sup> introduced by Lafarge called thermal permeability by analogy to the viscous

*ηφ jω*

1/2 <sup>1</sup>

*jω*

1/2

where *<sup>x</sup>* <sup>=</sup> *ωα*∞*<sup>ρ</sup> <sup>f</sup>*

*σφ*

*<sup>α</sup>*(*ω*) = <sup>−</sup> *ηφ*

where, *<sup>k</sup>*<sup>0</sup> is the viscous permeability given by: *<sup>k</sup>*<sup>0</sup> <sup>=</sup> *<sup>η</sup>*

*2.3.2. The equivalent model for the dynamic compressibility*

*<sup>β</sup>*(*ω*) = <sup>1</sup> <sup>−</sup> <sup>2</sup>(*<sup>γ</sup>* <sup>−</sup> <sup>1</sup>)

In a low frequency limit, Lafarge [9] showed that:

 <sup>1</sup> <sup>−</sup> <sup>1</sup> *j x*�

significant. The result is a two-parameter model, Λ� and *k*�

limit, Allard and Champoux [7, 8] showed the following behavior of *β*(*ω*):

 *η Prρ <sup>f</sup>*

Λ�

*<sup>β</sup>*(*ω*) = *<sup>γ</sup>* <sup>+</sup> (*<sup>γ</sup>* <sup>−</sup> <sup>1</sup>)*<sup>ρ</sup> <sup>f</sup> <sup>k</sup>*�

<sup>1</sup> <sup>−</sup> *<sup>M</sup>*� <sup>2</sup> *j x*�

then that of the simplest analytic function verifying these limits:

<sup>1</sup> <sup>−</sup> *<sup>M</sup>* 2 *j x*


## **2.3. Models of dynamic susceptibilities**

There are several models describing the inertial viscous and thermal exchanges between the fluid and the structure. We present here the most elaborate.

### *2.3.1. The model of Johnson et al. for the dynamic tortuosity*

Johnson et al [6] have proposed a simple model describing the function *α*(*ω*), which is the dynamic tortuosity, when the porous structure is saturated with a Newtonian viscous fluid. This model can be applied to the case of our porous structures saturated with air. Two model parameters, the tortuosity *α*∞ and the viscous characteristic length Λ are involved in a high-frequency limit:

$$a(\omega) = a\_{\infty} \left( 1 - \frac{2}{\Lambda} \left( \frac{\eta}{j\omega \rho\_f} \right)^{1/2} \right), \qquad \omega \longrightarrow \infty. \tag{1}$$

where *j* <sup>2</sup> <sup>=</sup> <sup>−</sup>1, *<sup>η</sup>* is the dynamic viscosity of the fluid, *<sup>ω</sup>* the angular frequency, *<sup>ρ</sup> <sup>f</sup>* the density of the fluid. In the low frequency limit, the expression of the dynamic tortuosity in the frequency domain is given by:

$$\mathfrak{a}(\omega) = -\frac{\eta \phi}{\rho\_f k\_0 j \omega'} , \quad \omega \longrightarrow 0,\tag{2}$$

where, *<sup>k</sup>*<sup>0</sup> is the viscous permeability given by: *<sup>k</sup>*<sup>0</sup> <sup>=</sup> *<sup>η</sup> <sup>σ</sup>* (with the dimension of a surface), *σ* is the specific resistance passage of the fluid. The meaning of the quantities *α*∞, Λ, *k*<sup>0</sup> have been given previously. Johnson et al [6] imply that these three geometric parameters provide sufficient information about the microgeometry. The proposed general expression for *α*(*ω*) is then that of the simplest analytic function verifying these limits:

$$a(\omega) = a\_{\infty} \left( 1 - \frac{1}{j\mathbf{x}} \sqrt{1 - \frac{M}{2}j\mathbf{x}} \right) \quad \text{where} \quad \mathbf{x} = \frac{\omega a\_{\infty} \rho\_f}{\sigma \phi} \quad \text{and} \quad M = \frac{8k\_0 a\_{\infty}}{\phi \Lambda^2}. \tag{3}$$

#### *2.3.2. The equivalent model for the dynamic compressibility*

4 Will-be-set-by-IN-TECH

• **Viscous permeability** *k*0: Static viscous permeability is related to the flow resistivity by the relation *k*<sup>0</sup> = *η*/*σ*, where *η* is the dynamic viscosity of the fluid. Viscous permeability is a geometrical parameter, homogeneous to a surface (it is expressed in *m*2), independent of the nature of the saturating fluid and depending only on the internal geometry of the porous material. This parameter represents the effective pore section for the flow of fluid, it is connected to the low frequency behavior of the fluid exchange between viscous and

• **Viscous characteristic length** Λ: The work of Johnson et al [6] in 1987 on the characterization of viscous effects at high frequencies has introduced the viscous

> � *<sup>S</sup> <sup>u</sup>*<sup>2</sup> *dS* �

incompressible perfect fluid, *S* the area of the interface between solid and fluid phases and *V* the volume of fluid. The length Λ is a geometrical parameter. This definition applies to a smooth interface fluid/solid, and for a low boundary layer thickness to the radius of curvature characteristic of the interface. When the pore surface has singularities (peaks), this definition of the characteristic length is no longer valid. The parameter Λ is an indicator of the size of the narrow neck of the pore, i.e. the privileged place of viscous

analogy with Johnson *et al* [6], a geometric parameter called thermal characteristic length

contact surface between the solid and fluid. The factor 2 is introduced so that in the case of cylindrical pores, Λ� is the pore radius. The length Λ� is an indicator of the size of large

similarities between these two effects. Starting from the equation of heat diffusion, Lafarge [9] introduced the equivalent of Darcy's law for the temperature excess, by introducing a

There are several models describing the inertial viscous and thermal exchanges between the

Johnson et al [6] have proposed a simple model describing the function *α*(*ω*), which is the dynamic tortuosity, when the porous structure is saturated with a Newtonian viscous fluid. This model can be applied to the case of our porous structures saturated with air. Two model parameters, the tortuosity *α*∞ and the viscous characteristic length Λ are involved in a

�1/2⎞

� *η jωρ <sup>f</sup>*

*<sup>V</sup> dV* . In other words <sup>Λ</sup>� is twice the ratio of fluid volume on the total

<sup>0</sup>: The study of viscous and thermal effects shows that there are

*<sup>V</sup> <sup>u</sup>*<sup>2</sup> *dV* where *<sup>u</sup>* is the speed of a microscopic

: In 1991, Allard and Champoux [8], introduced by

<sup>0</sup>, which is the inverse of the constant trapping.

⎠ , *ω* −→ ∞. (1)

<sup>Λ</sup> =

*σ* = (*S*/*d*)*R*.

exchanges.

given by <sup>2</sup>

saturating the porous structure.

• **Thermal characteristic length** Λ�

� *<sup>S</sup> dS* �

pores, privileged place of heat exchange.

new parameter: the thermal permeability *k*�

fluid and the structure. We present here the most elaborate.

*2.3.1. The model of Johnson et al. for the dynamic tortuosity*

*α*(*ω*) = *α*<sup>∞</sup>

⎛ <sup>⎝</sup><sup>1</sup> <sup>−</sup> <sup>2</sup> Λ

**2.3. Models of dynamic susceptibilities**

<sup>Λ</sup>� =

• **Thermal permeability** *k*�

high-frequency limit:

characteristic length Λ given by <sup>2</sup>

the length *d* of the sample and inversely to his section *S*. Then can be defined, the resistivity as the specific resistance to the passage of a fluid. Denoted *σ*, the resistivity is defined by

> A model similar to the previous one has been set for the frequency dependence of the thermal exchanges between the fluid and the structure, because of the heat capacity of the solid part is significant. The result is a two-parameter model, Λ� and *k*� <sup>0</sup>/*φ*, the quantities Λ� and *k*� <sup>0</sup>/*φ* for thermal effects playing respectively the role of Λ and *k*<sup>0</sup> for viscous effects. In a high-frequency limit, Allard and Champoux [7, 8] showed the following behavior of *β*(*ω*):

$$\beta(\omega) = 1 - \frac{2(\gamma - 1)}{\Lambda'} \left(\frac{\eta}{P\_t \rho\_f}\right)^{1/2} \left(\frac{1}{j\omega}\right)^{1/2}, \qquad \omega \longrightarrow \infty,\tag{4}$$

where *Pr* is the number of Prandt and *γ* the adiabatic constant. As before, this behavior implies that the thermal boundary layer thickness *δ*� becomes small compared to any characteristic radius of curvature of the interface. One can see that there is no equivalent of the tortuosity *α*<sup>∞</sup> for thermal effects, at very high frequency *β*(*ω*) tends to 1 as *α*(*ω*) tends to the tortuosity *α*∞.

In a low frequency limit, Lafarge [9] showed that:

$$\mathfrak{G}(\omega) = \gamma + \frac{(\gamma - 1)\rho\_f k\_0' P\_r}{\eta \phi j \omega}, \qquad \text{when} \qquad \omega \longrightarrow 0. \tag{5}$$

where *k*� <sup>0</sup>, which has the same size (area) that the Darcy permeability of *k*0, is a parameter analogous to the parameter *k*0, but adapted to the thermal problem. The proposed general expression for *β*(*ω*) is then that of the simplest analytic function satisfying the limits ((4) and (5):

$$\beta(\omega) = \gamma - (\gamma - 1) / \left[ 1 - \frac{1}{j \cdot \mathbf{x'}} \sqrt{1 - \frac{M'}{2} j \mathbf{x'}} \right] \quad \text{where} \quad \mathbf{x'} = \frac{\omega \rho\_f k\_0' P\_r}{\eta \phi} \quad \text{and} \quad M' = \frac{8 k\_0'}{\phi \Lambda'^2}.$$

the parameter *k*� <sup>0</sup> introduced by Lafarge called thermal permeability by analogy to the viscous permeability.

#### 6 Will-be-set-by-IN-TECH 132 Modeling and Measurement Methods for Acoustic Waves and for Acoustic Microdevices Transient Acoustic Wave Propagation in Porous media <sup>7</sup>

#### *2.3.3. Extensions of models by pride and Lafarge*

The development of Johnson et al in high frequency for the tortuosity dynamics stops after the term √ 1 *<sup>j</sup><sup>ω</sup>* , Pride *et al* [11] proposed a second correction term which is <sup>1</sup> *<sup>j</sup><sup>ω</sup>* . The asymptotic expansion of the dynamic tortuosity now is written as follows:

$$a(\omega) = a\_{\infty} \left( 1 - \frac{2}{\Lambda} \left( \frac{\eta}{j\omega \rho\_f} \right)^{1/2} - \frac{\sigma \phi (1 - \wp)}{j\omega \rho\_f \alpha\_{\infty}} \right), \qquad \omega \longrightarrow \infty. \tag{6}$$

The Pride parameter *℘* [11] is a geometrical parameter, dimensionless, connected to a correction term *α*<sup>0</sup> playing the role of inertia of fluid at low frequencies:

$$\wp = \frac{M}{4\left(\frac{a\_0}{a\_\infty} - 1\right)}, \qquad \text{where} \qquad M = \frac{8k\_0 a\_\infty}{\phi \Lambda^2}.$$

The expansion of low frequency dynamic tortuosity is given by:

$$\mathfrak{a}(\omega) = -\frac{\eta\phi}{\rho\_f k\_0 j\omega} + \mathfrak{a}\_{0\prime} \qquad \omega \longrightarrow 0 \tag{7}$$

One can notice that the models of Pride and Lafarge are reduced to those of Johnson and Allard when *℘* = 1 and *℘*� = 1. We will see later the influence of these parameters on the

The use of transient signals in the acoustic characterization of porous media is widely used experimentally. These signals have a wide frequency content making some complicated frequency approaches. The frequency methods are very effective for monochromatic signals, however, for transient signals, the temporal approach [2, 12–15] is best suited for several

• Time domain analysis is naturally bounded (limited) by the finite duration of impulses. • Time-modeling is often easier to develop because it is closer to the reality of experiments. • In many situations, the introduction of the "time" parameter make easier the analysis of

• For certain applications, it is fast because it avoids back and forth operations between time

• It provides an elegant solution of the direct problem necessary for the resolution of the

In return for these benefits, it is necessary to use a new mathematical formalism. In the frequency domain approach, the relevant quantities for the characterization of porous media are functions of frequency which can be interpreted as susceptibilities. While frequency methods lead to the frequency responses of the medium, in the temporal approach, one is

In the asymptotic regime corresponding to high frequencies, the fluid-structure interactions are described by expressions (1) et (4) according to the Johnson-Allard model. Writing the equations in the time domain is equivalent to taking the inverse Fourier transform of (1) and (4). The temporal equivalent of *jω* is *∂*/*∂t*, while the temporal equivalent of *jω* is a fraction0al derivative of order 1/2. The definition of fractional derivative of order *ν* is given

> *t* 0

where Γ is the Eulerian function of the second kind [17]. From the definition (9), the expressions of response factors *α*(*ω*) and *β*(*ω*) are then given in the time domain by [2]:

−→ *<sup>α</sup>*∞*δ*(*t*) + <sup>1</sup>

−→ *δ*(*t*) + *Kaχth*(*t*),

*ρ f*

*χv*(*t*),

(*<sup>t</sup>* <sup>−</sup> *<sup>u</sup>*)−*ν*−1*x*(*u*)*du*, (9)

Transient Acoustic Wave Propagation in Porous Media 133

**2.4. Temporal modeling of propagation: asymptotic regime**

• It is well suited for comparison between theory and experiment.

**2.5. Concept of fractional derivative and propagation equation**

*<sup>D</sup>ν*[*x*(*t*)] = <sup>1</sup>

*<sup>α</sup>*(*ω*) *<sup>t</sup>*

*<sup>β</sup>*(*ω*) *<sup>t</sup>*

Γ(−*ν*)

propagation.

reasons:

by [16]:

experimental results.

inverse problem.

and frequency domains by FFT.

interested in these impulse responses.

The new expression of the analytic function taking into account developments of (6) and (7) is given by:

$$\mathfrak{a}(\omega) = \mathfrak{a}\_{\infty} \left( 1 - \frac{1}{j\chi} \left( 1 - \wp + \wp \sqrt{1 - \frac{M}{2\wp^2} j\,\mathrm{x}} \right) \right) \dots$$

The value of *℘* depends on the geometry of the pores, for example in the case of cylindrical pores of circular section, it is shown that *℘* = 3/4. For rectangular sections, one obtains values close to 3/4.

For thermal effects, Lafarge [9] suggested a similar expression for the dynamic thermal compressibility, also introducing a dimensionless parameter *℘*� and gives a correction to the Allard and Champoux model at high frequencies. The expression of the dynamic compressibility is then given by:

$$\beta(\omega) = 1 - (\gamma - 1) \left[ \frac{2}{\Lambda'} \left( \frac{\eta}{Pr\rho\_f} \right)^{1/2} \left( \frac{1}{j\omega} \right)^{1/2} + \left( \frac{\eta}{\rho\_f Pr j\omega} \right) \cdot \left( \frac{(1 - \wp')\phi}{k\_0'} - \frac{4}{\Lambda'^2} \right) \right], \omega \longrightarrow \infty. \tag{8}$$

The parameter *℘*� is connected to *α*� <sup>0</sup>, thermal equivalent of *α*0, which describes the thermal inertia of fluid at low frequency:

$$
\wp' = \frac{M'}{4(\alpha'\_0 - 1)}.
$$

The general expression for the dynamic compressibility is then given by:

$$\beta(\omega) = \gamma - (\gamma - 1) / \left[ 1 - \frac{1}{j \cdot \mathbf{x'}} \left( 1 - \wp' + \wp' \sqrt{1 - \frac{M'}{2 \wp'^2} j \mathbf{x'}} \right) \right].$$

One can notice that the models of Pride and Lafarge are reduced to those of Johnson and Allard when *℘* = 1 and *℘*� = 1. We will see later the influence of these parameters on the propagation.

### **2.4. Temporal modeling of propagation: asymptotic regime**

6 Will-be-set-by-IN-TECH

The development of Johnson et al in high frequency for the tortuosity dynamics stops after

The Pride parameter *℘* [11] is a geometrical parameter, dimensionless, connected to a

The new expression of the analytic function taking into account developments of (6) and (7) is

1 − *℘* + *℘*

The value of *℘* depends on the geometry of the pores, for example in the case of cylindrical pores of circular section, it is shown that *℘* = 3/4. For rectangular sections, one obtains values

For thermal effects, Lafarge [9] suggested a similar expression for the dynamic thermal compressibility, also introducing a dimensionless parameter *℘*� and gives a correction to the Allard and Champoux model at high frequencies. The expression of the dynamic

> *η ρ <sup>f</sup> Prjω*

� .

� (<sup>1</sup> <sup>−</sup> *<sup>℘</sup>*�

*k*� 0

<sup>0</sup>, thermal equivalent of *α*0, which describes the thermal

<sup>1</sup> <sup>−</sup> *<sup>M</sup>*� <sup>2</sup>*℘*�<sup>2</sup> *j x*�

)*φ*

− 4 Λ�<sup>2</sup>

> �� .

�⎤

⎦ , *ω* −→ ∞.

(8)

�

<sup>1</sup> <sup>−</sup> *<sup>M</sup>* <sup>2</sup>*℘*<sup>2</sup> *j x*

<sup>−</sup> *σφ*(<sup>1</sup> <sup>−</sup> *<sup>℘</sup>*) *jωρ <sup>f</sup> α*<sup>∞</sup>

� , where *<sup>M</sup>* <sup>=</sup> <sup>8</sup>*k*0*α*<sup>∞</sup>

⎞

*<sup>φ</sup>*Λ<sup>2</sup> .

+ *α*0, *ω* −→ 0 (7)

�� . *<sup>j</sup><sup>ω</sup>* . The asymptotic

⎠ , *ω* −→ ∞. (6)

*<sup>j</sup><sup>ω</sup>* , Pride *et al* [11] proposed a second correction term which is <sup>1</sup>

�1/2

*2.3.3. Extensions of models by pride and Lafarge*

*α*(*ω*) = *α*<sup>∞</sup>

expansion of the dynamic tortuosity now is written as follows:

� *η jωρ <sup>f</sup>*

correction term *α*<sup>0</sup> playing the role of inertia of fluid at low frequencies:

*<sup>α</sup>*(*ω*) = <sup>−</sup> *ηφ*

� <sup>1</sup> <sup>−</sup> <sup>1</sup> *jx*

�1/2 � <sup>1</sup>

*jω*

The general expression for the dynamic compressibility is then given by:

� <sup>1</sup> <sup>−</sup> <sup>1</sup> *j x*�

�1/2 + �

*<sup>℘</sup>*� <sup>=</sup> *<sup>M</sup>*� 4(*α*� <sup>0</sup> − 1) .

�

1 − *℘*� + *℘*�

�

*ρ <sup>f</sup> k*<sup>0</sup> *jω*

�

⎛ <sup>⎝</sup><sup>1</sup> <sup>−</sup> <sup>2</sup> Λ

*<sup>℘</sup>* <sup>=</sup> *<sup>M</sup>* 4 � *<sup>α</sup>*<sup>0</sup> *<sup>α</sup>*<sup>∞</sup> − 1

*α*(*ω*) = *α*<sup>∞</sup>

The expansion of low frequency dynamic tortuosity is given by:

the term √

given by:

close to 3/4.

*β*(*ω*)=1−(*γ*−1)

compressibility is then given by:

⎡ ⎣ 2 Λ�

The parameter *℘*� is connected to *α*�

inertia of fluid at low frequency:

� *η Prρ <sup>f</sup>*

*β*(*ω*) = *γ* − (*γ* − 1)/

1

The use of transient signals in the acoustic characterization of porous media is widely used experimentally. These signals have a wide frequency content making some complicated frequency approaches. The frequency methods are very effective for monochromatic signals, however, for transient signals, the temporal approach [2, 12–15] is best suited for several reasons:


In return for these benefits, it is necessary to use a new mathematical formalism. In the frequency domain approach, the relevant quantities for the characterization of porous media are functions of frequency which can be interpreted as susceptibilities. While frequency methods lead to the frequency responses of the medium, in the temporal approach, one is interested in these impulse responses.

#### **2.5. Concept of fractional derivative and propagation equation**

In the asymptotic regime corresponding to high frequencies, the fluid-structure interactions are described by expressions (1) et (4) according to the Johnson-Allard model. Writing the equations in the time domain is equivalent to taking the inverse Fourier transform of (1) and (4). The temporal equivalent of *jω* is *∂*/*∂t*, while the temporal equivalent of *jω* is a fraction0al derivative of order 1/2. The definition of fractional derivative of order *ν* is given by [16]:

$$[D^{\nu}]\mathbf{x}(t)] = \frac{1}{\Gamma(-\nu)} \int\_{0}^{t} (t-u)^{-\nu-1} \mathbf{x}(u) du,\tag{9}$$

where Γ is the Eulerian function of the second kind [17]. From the definition (9), the expressions of response factors *α*(*ω*) and *β*(*ω*) are then given in the time domain by [2]:

> *<sup>α</sup>*(*ω*) *<sup>t</sup>* −→ *<sup>α</sup>*∞*δ*(*t*) + <sup>1</sup> *ρ f χv*(*t*), *<sup>β</sup>*(*ω*) *<sup>t</sup>* −→ *δ*(*t*) + *Kaχth*(*t*),

#### 8 Will-be-set-by-IN-TECH 134 Modeling and Measurement Methods for Acoustic Waves and for Acoustic Microdevices Transient Acoustic Wave Propagation in Porous media <sup>9</sup>

where *δ*(*t*) is the Dirac distribution and the operators, *Ka* the modulus of the fluid, *χv*(*t*) and *χth*(*t*) are given by:

$$\chi\_{\upsilon}(t) = \frac{2\rho\_f \alpha\_{\infty}}{\Lambda} \sqrt{\frac{\eta}{\pi \rho\_f}} \, t^{-1/2} \Big/$$

$$\chi\_{th}(t) = \frac{2(\gamma - 1)}{K\_d \Lambda'} \sqrt{\frac{\eta}{\pi P\_r \rho\_f}} \, t^{-1/2} \Big/$$

In this model, the time convolution of *t* <sup>−</sup>1/2 with a function, is interpreted as an operator of fractional derivative. The acoustic fields satisfy the Euler equation and the constitutive equation [2] given by:

$$
\rho\_f \alpha\_{\infty} \,\partial\_t \mathbf{v}(\mathbf{r}, t) + \int\_0^t \chi\_{\mathcal{V}}(t - t') \,\partial\_t \mathbf{v}(\mathbf{r}, t') \,dt' = -\nabla p(\mathbf{r}, t), \tag{10}
$$

respectively, The coefficient *<sup>A</sup>* gives the wave velocity *<sup>c</sup>* <sup>=</sup> 1/

*<sup>p</sup>*(*x*, *<sup>t</sup>*) = *<sup>t</sup>*

0

*G*(*x*, *t* − *t*

To solve the propagation equation and obtain the Green function of the medium, we solve the equation (13) using the method of Laplace transform [18], taking into account the following

> ∞ 0

*<sup>p</sup>*(*x*, *<sup>t</sup>*) = <sup>L</sup>−<sup>1</sup> [*P*(*x*, *<sup>z</sup>*)] .

*z*

where *H*(*t*) is the Heaviside jump function, the Laplace transform of the wave equation (13)

 1 *<sup>c</sup>*<sup>2</sup> <sup>+</sup>*<sup>B</sup>*

*P*(*x*, *z*)=−

where *p*(*x*, *t*) is the sound pressure in the porous medium and *p*(0, *t*) the incident signal.

**2.6. Solution of the propagation equation: Green function (Johnson-Allard**

*<sup>p</sup>*(*x*, *<sup>t</sup>*)|*t*=<sup>0</sup> <sup>=</sup> 0 and *<sup>∂</sup><sup>p</sup>*

Denote by *P*(*x*, *z*) the Laplace transform of *p*(*x*, *t*) defined by:

<sup>L</sup> [*δ*(*t*)] <sup>=</sup> 1, <sup>L</sup> [*H*(*t*)] <sup>=</sup> <sup>1</sup>

*P*(*x*, *z*) = L [*p*(*x*, *t*)] =

�

*∂t*

exp(−*zt*)*p*(*x*, *t*)*dt*,

and L

*π z*

 1 √*t* = *π z* ,

*zp*(*x*, 0)+ *<sup>∂</sup>p*(*x*, 0)

*∂t*

−*Cp*(*x*, 0).

(15)

)*p*(0, *t*)*dt*�

,

changes the speed *c*<sup>0</sup> =

**model)**

is:

*∂*2*P*(*x*, *z*) *<sup>∂</sup>x*<sup>2</sup> <sup>−</sup>*z*<sup>2</sup>

the porous medium [3, 18], defined by:

initial conditions of causality:

and the inverse Laplace transform by:

Using the following relationship:

 1 *<sup>c</sup>*<sup>2</sup> <sup>+</sup> *C z* +*B π z* 

the porous medium. *α*∞ plays the role of the refractive index of the porous medium that

*Ka*/*ρ <sup>f</sup>* in free space to *c*<sup>0</sup> = *c*/

The other coefficients depend essentially on characteristic lengths Λ and Λ� and express the fluid-structure visco-thermal interactions. The constant *B* describes the dispersion of the signal and *C* the attenuation of the amplitude of the wave (without dispersion). The propagation equation (13) describes the evolution of the acoustic wave inside the porous material. The boundary conditions have not been introduced, the porous medium is assumed to be infinite. The coefficients of the equation of propagation do not depend on the porosity, in fact, porosity appears at the interfaces of the porous material via the relationship of continuity of acoustic flow (the flow rate is equal to the porosity multiplied by the speed). This parameter appears naturally at the reflection and transmission of an acoustic wave through a slice of porous material. The solution of the wave equation (13) is given by the Green function *G* of

*ρ <sup>f</sup> α*∞/*Ka* in the air saturating

Transient Acoustic Wave Propagation in Porous Media 135

<sup>√</sup>*α*<sup>∞</sup> in the porous medium.


$$\frac{1}{K\_a}\partial\_t p(\mathbf{r}, t) + \int\_0^t \chi\_{th}(t - t')\,\partial\_t p(\mathbf{r}, t')\,dt' = -\nabla.\mathbf{v}(\mathbf{r}, t). \tag{11}$$

These constitutive relations in the time domain satisfy the principle of causality. In these equations, *p* is the acoustic pressure, **v** the particle velocity. The term *α*∞*δ*(*t*) reflects the instantaneous response of the porous medium and describes the inertial coupling between fluid and structure. For instantaneous responses, we mean that the time response is smaller than the characteristic time scale of the change in the acoustic field. The susceptibilities *χv* and *χth* are operators of memory that determine the dispersion of the medium.

We assume that the medium varies with the thickness *x* only, and that the incident wave is plane and normal to the surface of the material. Sound pressure is denoted by *p*(*x*, *t*). We assume that the pressure field is zero for moments earlier to 0. The wave equation for acoustic pressure field of a porous medium having a dispersive rigid structure is obtained from the constitutive equations (10, 11), and is of the form:

$$\left[\partial\_{\mathbf{x}}^{2}p(\mathbf{x},t) - \frac{1}{c\_{0}^{2}}\left[a\_{\infty}\partial\_{t}^{2}p(\mathbf{x},t) + \left(a\_{\infty}K\_{a}\chi\_{th} + \frac{\chi\_{v}}{\rho\_{f}} + c\_{0}^{2}\chi\_{th}\*\chi\_{v}\right)\*\partial\_{t}^{2}p(\mathbf{x},t)\right] = 0,\tag{12}$$

where *c*<sup>0</sup> = (*Ka*/*ρ <sup>f</sup>*)1/2 is the sound speed in the free fluid. The following notation is used for convolution integrals of two causal functions *f*(*x*,*t*) and *g*(*x*, *t*);

$$[f \ast g](\mathbf{x}, t) = \int\_0^t f(\mathbf{x}, t - t') g(\mathbf{x}, t') dt'.$$

The propagation equation (12) can be written as:

$$\frac{\partial^2 p(\mathbf{x}, t)}{\partial \mathbf{x}^2} - A \frac{\partial^2 p(\mathbf{x}, t)}{\partial t^2} - B \int\_0^t \frac{\partial^2 p(\mathbf{x}, t) / \partial t'^2}{\sqrt{t - t'}} dt' - C \frac{\partial p(\mathbf{x}, t)}{\partial t} = 0,\tag{13}$$

where the coefficients *A*, *B* and *C* are constant and given by:

$$A = \frac{\rho\_f a\_{\infty}}{K\_a}, \quad B = \frac{2a\_{\infty}}{K\_a} \sqrt{\frac{\rho\_f \eta}{\pi}} \left(\frac{1}{\Lambda} + \frac{\gamma - 1}{\sqrt{Pr \Lambda'}}\right) \quad \mathbb{C} = \frac{4a\_{\infty}(\gamma - 1)\eta}{K\_a \Lambda \Lambda' \sqrt{Pr}}, \quad \mathbb{C} = \frac{2a\_{\infty}(\gamma - 1)\eta}{\sqrt{Pr \Lambda'}}$$

respectively, The coefficient *<sup>A</sup>* gives the wave velocity *<sup>c</sup>* <sup>=</sup> 1/ *ρ <sup>f</sup> α*∞/*Ka* in the air saturating the porous medium. *α*∞ plays the role of the refractive index of the porous medium that changes the speed *c*<sup>0</sup> = *Ka*/*ρ <sup>f</sup>* in free space to *c*<sup>0</sup> = *c*/ <sup>√</sup>*α*<sup>∞</sup> in the porous medium. The other coefficients depend essentially on characteristic lengths Λ and Λ� and express the fluid-structure visco-thermal interactions. The constant *B* describes the dispersion of the signal and *C* the attenuation of the amplitude of the wave (without dispersion). The propagation equation (13) describes the evolution of the acoustic wave inside the porous material. The boundary conditions have not been introduced, the porous medium is assumed to be infinite. The coefficients of the equation of propagation do not depend on the porosity, in fact, porosity appears at the interfaces of the porous material via the relationship of continuity of acoustic flow (the flow rate is equal to the porosity multiplied by the speed). This parameter appears naturally at the reflection and transmission of an acoustic wave through a slice of porous material. The solution of the wave equation (13) is given by the Green function *G* of the porous medium [3, 18], defined by:

$$p(\mathbf{x},t) = \int\_0^t \mathbf{G}(\mathbf{x}, t - t') p(\mathbf{0}, t) dt',$$

where *p*(*x*, *t*) is the sound pressure in the porous medium and *p*(0, *t*) the incident signal.

### **2.6. Solution of the propagation equation: Green function (Johnson-Allard model)**

To solve the propagation equation and obtain the Green function of the medium, we solve the equation (13) using the method of Laplace transform [18], taking into account the following initial conditions of causality:

$$p(\mathbf{x},t)|\_{t=0} = 0 \quad \text{and} \quad \frac{\partial p}{\partial t}|\_{t=0} = 0. \tag{14}$$

Denote by *P*(*x*, *z*) the Laplace transform of *p*(*x*, *t*) defined by:

$$P(\mathbf{x}, z) = \mathcal{L}\left[p(\mathbf{x}, t)\right] = \int\_0^\infty \exp(-zt)p(\mathbf{x}, t)dt,$$

and the inverse Laplace transform by:

$$p(\mathbf{x}, t) = \mathcal{L}^{-1} \left[ P(\mathbf{x}, z) \right].$$

Using the following relationship:

8 Will-be-set-by-IN-TECH

where *δ*(*t*) is the Dirac distribution and the operators, *Ka* the modulus of the fluid, *χv*(*t*) and

 *η πρ <sup>f</sup> t* <sup>−</sup>1/2,

> *η πPrρ <sup>f</sup> t* <sup>−</sup>1/2,

<sup>−</sup>1/2 with a function, is interpreted as an operator

) *dt*� = −∇*p*(**r**, *t*), (10)

) *dt*� = −∇.**v**(**r**, *t*). (11)

2*ρ <sup>f</sup> α*<sup>∞</sup> Λ

*Ka*Λ�

*χv*(*t* − *t* � ) *∂t***v**(**r**, *t* �

*χth*(*t* − *t* � ) *∂<sup>t</sup> p*(**r**, *t* �

and *χth* are operators of memory that determine the dispersion of the medium.

of fractional derivative. The acoustic fields satisfy the Euler equation and the constitutive

These constitutive relations in the time domain satisfy the principle of causality. In these equations, *p* is the acoustic pressure, **v** the particle velocity. The term *α*∞*δ*(*t*) reflects the instantaneous response of the porous medium and describes the inertial coupling between fluid and structure. For instantaneous responses, we mean that the time response is smaller than the characteristic time scale of the change in the acoustic field. The susceptibilities *χv*

We assume that the medium varies with the thickness *x* only, and that the incident wave is plane and normal to the surface of the material. Sound pressure is denoted by *p*(*x*, *t*). We assume that the pressure field is zero for moments earlier to 0. The wave equation for acoustic pressure field of a porous medium having a dispersive rigid structure is obtained from the

*<sup>α</sup>*∞*Kaχth* <sup>+</sup> *<sup>χ</sup><sup>v</sup>*

where *c*<sup>0</sup> = (*Ka*/*ρ <sup>f</sup>*)1/2 is the sound speed in the free fluid. The following notation is used for

*f*(*x*, *t* − *t* � )*g*(*x*, *t* � )*dt*� .

*∂*<sup>2</sup> *p*(*x*, *t*)/*∂t*

<sup>Λ</sup> <sup>+</sup> *<sup>γ</sup>* <sup>−</sup> <sup>1</sup> <sup>√</sup>*Pr*Λ�

�2 <sup>√</sup>*<sup>t</sup>* <sup>−</sup> *<sup>t</sup>*� *dt*� <sup>−</sup> *<sup>C</sup>∂p*(*x*, *<sup>t</sup>*)

*<sup>C</sup>* <sup>=</sup> <sup>4</sup>*α*∞(*<sup>γ</sup>* <sup>−</sup> <sup>1</sup>)*<sup>η</sup> Ka*ΛΛ�

 *t* 0

> *t* 0

> > 1

*ρ f* + *c*<sup>2</sup>

<sup>0</sup>*χth* ∗ *χ<sup>v</sup>*

 <sup>∗</sup> *<sup>∂</sup>*<sup>2</sup> *<sup>t</sup> p*(*x*, *t*)

*<sup>∂</sup><sup>t</sup>* <sup>=</sup> 0, (13)

<sup>√</sup>*Pr* ,

= 0, (12)

*<sup>χ</sup>th*(*t*) = <sup>2</sup>(*<sup>γ</sup>* <sup>−</sup> <sup>1</sup>)

 *t* 0

 *t* 0

*χv*(*t*) =

*χth*(*t*) are given by:

equation [2] given by:

*∂*2

*<sup>x</sup> <sup>p</sup>*(*x*, *<sup>t</sup>*) <sup>−</sup> <sup>1</sup>

*c*2 0

*∂*<sup>2</sup> *p*(*x*, *t*)

*<sup>A</sup>* <sup>=</sup> *<sup>ρ</sup> <sup>f</sup> <sup>α</sup>*<sup>∞</sup> *Ka*

 *α*∞*∂*<sup>2</sup>

In this model, the time convolution of *t*

1 *Ka*

constitutive equations (10, 11), and is of the form:

The propagation equation (12) can be written as:

*<sup>∂</sup>x*<sup>2</sup> <sup>−</sup> *<sup>A</sup> <sup>∂</sup>*<sup>2</sup> *<sup>p</sup>*(*x*, *<sup>t</sup>*)

where the coefficients *A*, *B* and *C* are constant and given by:

, *<sup>B</sup>* <sup>=</sup> <sup>2</sup>*α*<sup>∞</sup>

*Ka*

*<sup>t</sup> p*(*x*, *t*) +

convolution integrals of two causal functions *f*(*x*,*t*) and *g*(*x*, *t*);

[ *f* ∗ *g*](*x*, *t*) =

*<sup>∂</sup>t*<sup>2</sup> <sup>−</sup> *<sup>B</sup>*

*ρ <sup>f</sup> η π*

*ρ <sup>f</sup> α*<sup>∞</sup> *∂t***v**(**r**, *t*) +

*∂<sup>t</sup> p*(**r**, *t*) +

$$
\mathcal{L}\left[\delta(t)\right] = 1, \qquad \qquad \mathcal{L}\left[H(t)\right] = \frac{1}{z} \qquad \text{and} \qquad \mathcal{L}\left[\frac{1}{\sqrt{t}}\right] = \sqrt{\frac{\pi}{z}}.
$$

where *H*(*t*) is the Heaviside jump function, the Laplace transform of the wave equation (13) is:

$$\frac{\partial^2 P(\mathbf{x}, z)}{\partial \mathbf{x}^2} - z^2 \left( \frac{1}{c^2} + \frac{\mathbb{C}}{z} + B\sqrt{\frac{\pi}{z}} \right) P(\mathbf{x}, z) = -\left( \frac{1}{c^2} + B\sqrt{\frac{\pi}{z}} \right) \left( z p(\mathbf{x}, 0) + \frac{\partial p(\mathbf{x}, 0)}{\partial t} \right) - \mathbb{C} p(\mathbf{x}, 0). \tag{15}$$

#### 10 Will-be-set-by-IN-TECH 136 Modeling and Measurement Methods for Acoustic Waves and for Acoustic Microdevices Transient Acoustic Wave Propagation in Porous media <sup>11</sup>

Taking into account the initial conditions (14), equation (15) simplifies to

$$\frac{\partial^2 P(\mathbf{x}, z)}{\partial \mathbf{x}^2} - z^2 \left( \frac{1}{c^2} + \frac{\mathbf{C}}{z} + B\sqrt{\frac{\pi}{z}} \right) P(\mathbf{x}, z) = 0\_\mathbf{v}$$

which is a second order differential equation with real coefficients constant. The general solution of equation is:

$$P(x,z) = e^{-\frac{x}{\varepsilon}}\sqrt{f(z)}\,\varphi(z) + e^{\frac{x}{\varepsilon}}\sqrt{f(z)}\,\psi(z),$$

where *ϕ*(*z*) and *ψ*(*z*) are two independent functions of *x* and:

$$f(z) = z^2 \left( 1 + \frac{\mathbb{C}c^2}{z} + \mathrm{B}c^2 \sqrt{\frac{\pi}{z}} \right) = z(z + b' \sqrt{z} + c'). \tag{16}$$

In equation (16), the constants *b*� and *c*� are positive and are given by:

$$b' = \mathcal{B}c^2 \sqrt{\pi} \quad \text{and} \quad c' = \mathcal{C}c^2. \tag{17}$$

When Δ<sup>2</sup> *>* 0, the general solution of the propagation equation is given by:

<sup>−</sup> *<sup>b</sup>*�2*x*<sup>2</sup> <sup>16</sup>*c*<sup>2</sup>(*τ*−*x*/*c*)

1 *<sup>ξ</sup>*3/2 � <sup>1</sup> −1 exp

⎛

� *μ*Δ

<sup>−</sup> *<sup>b</sup>*�2*x*<sup>2</sup> <sup>16</sup>*c*<sup>2</sup>(*t*−*x*/*c*)

�1/2 *<sup>∂</sup>*−1/2

*∂*−1*x*(*t*) *<sup>∂</sup>t*−<sup>1</sup> <sup>=</sup>

> 0 *χ*� *<sup>v</sup>*(*t* − *t* � ) *∂t***v**(**r**, *t* �

0 *χ*� *th*(*t* − *t* � ) *∂<sup>t</sup> p*(**r**, *t* �

⎜⎝

×

1 (*τ*−*x*/*c*)3/2 exp �

1 �(*<sup>τ</sup>* <sup>−</sup> *<sup>ξ</sup>*)<sup>2</sup> <sup>−</sup> *<sup>x</sup>*2/*c*<sup>2</sup>

In this case, the Green function [18] is given by:

1 (*t*−*x*/*c*)3/2 exp �

0, if 0 ≤ *t* ≤ *x*/*c*,

�

⎛

⎜⎝<sup>−</sup>

0, if 0 ≤ *t* ≤ *x*/*c*,

�

**2.7. Solution of the propagation equation: Green function (Pride-Lafarge model)** Extensions of Pride-Lafarge (6,8) are expressed in the time domain [3] using the definition of

*<sup>∂</sup>t*−1/2 <sup>+</sup> *σφ*(<sup>1</sup> <sup>−</sup> *<sup>℘</sup>*)

*<sup>∂</sup>t*−1/2 <sup>+</sup>

*α*˜(*t*) and *β*˜(*t*) are the operators of tortuosity and compressibility. In these equations, the

� *t* 0 *x*(*t* � )*dt*� ,

�1/2 *<sup>∂</sup>*−1/2

*ρ <sup>f</sup> α*<sup>∞</sup>

� *η ρ <sup>f</sup> Pr*

+Δ � *<sup>τ</sup>*−*x*/*<sup>c</sup>* <sup>0</sup> *h*�

> � *μ*Δ

�(*<sup>τ</sup>* <sup>−</sup> *<sup>ξ</sup>*)<sup>2</sup> <sup>−</sup> *<sup>x</sup>*2/*c*<sup>2</sup> <sup>+</sup> *<sup>b</sup>*�

+ Δ � *<sup>t</sup>*−*x*/*<sup>c</sup>* <sup>0</sup> *h*�

> *∂*−<sup>1</sup> *∂t*−<sup>1</sup>

⎞ ⎠ ,

� �

<sup>−</sup>1/2 represents the fractional derivative. In this case, the basic

(1 − *℘*�

*k*� 0 )*φ*

<sup>−</sup> <sup>4</sup> Λ�<sup>2</sup>

) *dt*� = −∇*p*(**r**, *t*), (24)

) *dt*� = −∇.**v**(**r**, *t*), (25)

� *∂*−<sup>1</sup> *∂t*−<sup>1</sup> ⎞ ⎠ ,

(*ξ*)*dξ* �

�(*<sup>τ</sup>* <sup>−</sup> *<sup>ξ</sup>*)<sup>2</sup> <sup>−</sup> *<sup>x</sup>*2/*c*<sup>2</sup> <sup>+</sup> *<sup>b</sup>*�

16*ξ*

Transient Acoustic Wave Propagation in Porous Media 137

(*τ* − *ξ*)

<sup>8</sup>*<sup>ξ</sup>* <sup>−</sup> <sup>1</sup>

(*ξ*)*dξ* � �2

*p*(0, *t*−*τ*)*dτ*, *t > x*/*c*.

(*τ* − *ξ*)

⎞

⎟⎠

, if *<sup>t</sup> <sup>&</sup>gt; <sup>x</sup>*/*c*. (23)

(22)

�2

*μdμ* �<sup>1</sup> <sup>−</sup> *<sup>μ</sup>*<sup>2</sup>

⎞

⎟⎠

.

*p*(*x*, *t*)=

where

*h*�

⎧ ⎨ ⎩

(*ξ*) = <sup>−</sup> <sup>1</sup>

*G*(*x*, *t*) =

*x c* � *t x*/*c* � *<sup>b</sup>*� 4 <sup>√</sup>*<sup>π</sup>*

4*π*3/2

⎧ ⎨ ⎩

fractional derivative:

*α*˜(*t*) = *α*<sup>∞</sup>

operator *∂*<sup>−</sup>1/*∂t*

*x c* � *<sup>b</sup>*� 4 <sup>√</sup>*<sup>π</sup>*

⎛

*<sup>β</sup>*˜(*t*) = *<sup>δ</sup>*(*t*)+(*<sup>γ</sup>* <sup>−</sup> <sup>1</sup>)

and semi-operator *∂*<sup>−</sup>1/2/*∂t*

equations are given by:

<sup>⎝</sup>*δ*(*t*) + <sup>2</sup>

Λ

� *η ρ f*

⎛ ⎝ 2 Λ�

� *η Prρ <sup>f</sup>*

<sup>−</sup><sup>1</sup> represents the time integral:

*<sup>ρ</sup> <sup>f</sup> <sup>α</sup>*<sup>∞</sup> *<sup>∂</sup>t***v**(**r**, *<sup>t</sup>*) + � *<sup>t</sup>*

*<sup>∂</sup><sup>t</sup> <sup>p</sup>*(**r**, *<sup>t</sup>*) + � *<sup>t</sup>*

1 *Ka*

Retaining the finite solution at infinity, which corresponds to the physical solution of our problem, we have:

*P*(*x*, *z*) = *e* − *x c* √*f*(*z*) *ϕ*(*z*).

The solution of equation (13) is the inverse Laplace transform of *P*(*x*, *z*). We obtain:

$$p(\mathbf{x},t) = \mathcal{L}^{-1}\left(e^{-\frac{\mathbf{z}}{\varepsilon}\sqrt{f(z)}}\rho(z)\right) = \mathcal{L}^{-1}\left(e^{-\frac{\mathbf{z}}{\varepsilon}\sqrt{f(z)}}\right) \* \mathcal{L}^{-1}\left(\rho(z)\right).$$

The core problem is the calculation of the inverse Laplace transform of the term <sup>L</sup>−<sup>1</sup> *e*<sup>−</sup> *<sup>x</sup> c* √*f*(*z*) . This has been well studied in reference [18]. By letting:

$$
\Delta^2 = b'^2 - 4c' \,\text{.}\tag{18}
$$

it is easy to check that Δ<sup>2</sup> is always positive in the Johnson-Allard model. The expression of *f*(*z*) (Eq. 16) can be written as

$$f(z) = \left(z + \frac{b'}{2}\sqrt{z}\right)^2 - \left(\frac{\Delta\sqrt{z}}{2}\right)^2. \tag{19}$$

When Δ = 0, the solution of the propagation equation(13) is given [18] by:

$$p(\mathbf{x},t) = \begin{cases} 0, & \text{if } \quad 0 \le t \le \mathbf{x}/c, \\ \frac{1}{4\sqrt{\pi}} \frac{b^t \mathbf{x}}{c} \int\_{\mathbf{x}/c}^t \frac{1}{(\tau - \mathbf{x}/c)^{3/2}} \exp\left(-\frac{b^{t2}\mathbf{x}^2}{16c^2(\tau - \mathbf{x}/c)}\right) p(0, t-\tau) d\tau, & \text{if } \quad t > \mathbf{x}/c, \end{cases} \tag{20}$$

where *<sup>p</sup>*(0, *<sup>t</sup>*) = <sup>L</sup>−<sup>1</sup> (*ϕ*(*z*)). In this case the Green function is given [18, 19] by:

$$G(\mathbf{x},t) = \begin{cases} 0, & \text{if } \quad 0 \le t \le \mathbf{x}/c, \\ \frac{1}{4\sqrt{\pi}} \frac{b^\circ \mathbf{x}}{c} \frac{1}{(t-\mathbf{x}/c)^{3/2}} \exp\left(-\frac{b^2 \mathbf{x}^2}{16c^2(t-\mathbf{x}/c)}\right), & \text{if } \qquad t > \mathbf{x}/c, \end{cases} \tag{21}$$

When Δ<sup>2</sup> *>* 0, the general solution of the propagation equation is given by:

$$p(\mathbf{x},t) = \begin{cases} 0, & \text{if} \quad 0 \le t \le \mathbf{x}/c, \\ \frac{\mathbf{x}}{c} \int\_{\mathbf{x}/c}^{t} \left( \frac{b'}{4\sqrt{\pi}} \frac{1}{(\tau - \mathbf{x}/c)^{3/2}} \exp\left(-\frac{b'^2 \mathbf{x}^2}{16c^2(\tau - \mathbf{x}/c)}\right) + \Delta \int\_{0}^{\tau - \mathbf{x}/c} h'(\xi) d\xi \right) p(0, t - \tau) d\tau, t > \mathbf{x}/c. \end{cases} (22)$$

where

10 Will-be-set-by-IN-TECH

which is a second order differential equation with real coefficients constant. The general

<sup>√</sup>*f*(*z*)*ϕ*(*z*) + *<sup>e</sup>*

*π z* 

Retaining the finite solution at infinity, which corresponds to the physical solution of our

− *x c* √*f*(*z*)

The core problem is the calculation of the inverse Laplace transform of the term

<sup>=</sup> <sup>L</sup>−<sup>1</sup> *e* − *x c* √*f*(*z*) 

*x c* √*f*(*z*)

= *z*(*z* + *b*�

*ϕ*(*z*).

*P*(*x*, *z*) = 0,

*ψ*(*z*),

*<sup>b</sup>*� <sup>=</sup> *Bc*2√*<sup>π</sup>* and *<sup>c</sup>*� <sup>=</sup> *Cc*2. (17)

2

<sup>√</sup>*<sup>z</sup>* <sup>+</sup> *<sup>c</sup>*�

∗ L−<sup>1</sup> (*ϕ*(*z*)).

, (18)

. (19)

*<sup>p</sup>*(0, *<sup>t</sup>* <sup>−</sup> *<sup>τ</sup>*)*dτ*, if *<sup>t</sup> <sup>&</sup>gt; <sup>x</sup>*/*c*, (20)

, if *<sup>t</sup> <sup>&</sup>gt; <sup>x</sup>*/*c*, (21)

). (16)

Taking into account the initial conditions (14), equation (15) simplifies to

*<sup>P</sup>*(*x*, *<sup>z</sup>*) = *<sup>e</sup>*<sup>−</sup> *<sup>x</sup>*

In equation (16), the constants *b*� and *c*� are positive and are given by:

where *ϕ*(*z*) and *ψ*(*z*) are two independent functions of *x* and:

 1 + *Cc*<sup>2</sup> *z*

*f*(*z*) = *z*<sup>2</sup>

*<sup>p</sup>*(*x*, *<sup>t</sup>*) = <sup>L</sup>−<sup>1</sup>

*f*(*z*) (Eq. 16) can be written as

1 4 <sup>√</sup>*<sup>π</sup> b*� *x c t <sup>x</sup>*/*<sup>c</sup>* <sup>1</sup>

*G*(*x*,*t*) =

 *e* − *x c*

*f*(*z*) =

(*τ*−*x*/*c*)3/2 exp

0, if 0 <sup>≤</sup> *<sup>t</sup>* <sup>≤</sup> *<sup>x</sup>*/*c*,

0, if 0 <sup>≤</sup> *<sup>t</sup>* <sup>≤</sup> *<sup>x</sup>*/*c*,

1 4 <sup>√</sup>*<sup>π</sup> b*� *x <sup>c</sup>* <sup>1</sup>

 *z* + *b*� 2 √*z* 2 − Δ √*z* 2

When Δ = 0, the solution of the propagation equation(13) is given [18] by:

where *<sup>p</sup>*(0, *<sup>t</sup>*) = <sup>L</sup>−<sup>1</sup> (*ϕ*(*z*)). In this case the Green function is given [18, 19] by:

(*t*−*x*/*c*)3/2 exp

 1 *<sup>c</sup>*<sup>2</sup> <sup>+</sup> *C z* + *B π z* 

*c*

*P*(*x*, *z*) = *e*

<sup>√</sup>*f*(*z*)*ϕ*(*z*)

The solution of equation (13) is the inverse Laplace transform of *P*(*x*, *z*). We obtain:

. This has been well studied in reference [18]. By letting:

<sup>Δ</sup><sup>2</sup> <sup>=</sup> *<sup>b</sup>*�<sup>2</sup> <sup>−</sup> <sup>4</sup>*c*�

it is easy to check that Δ<sup>2</sup> is always positive in the Johnson-Allard model. The expression of

<sup>−</sup> *<sup>b</sup>*�2*x*<sup>2</sup> <sup>16</sup>*c*<sup>2</sup>(*τ*−*x*/*c*)

<sup>−</sup> *<sup>b</sup>*�2*x*<sup>2</sup> <sup>16</sup>*c*<sup>2</sup>(*t*−*x*/*c*)

+ *Bc*<sup>2</sup>

*∂*2*P*(*x*, *z*) *<sup>∂</sup>x*<sup>2</sup> <sup>−</sup> *<sup>z</sup>*<sup>2</sup>

solution of equation is:

problem, we have:

<sup>L</sup>−<sup>1</sup> *e*<sup>−</sup> *<sup>x</sup> c* √*f*(*z*) 

*p*(*x*, *t*) =

$$h'(\xi) = -\frac{1}{4\pi^{3/2}} \frac{1}{\sqrt{(\tau-\xi)^2 - \mathbf{x}^2/c^2}} \frac{1}{\xi^{3/2}} \int\_{-1}^1 \exp\left(-\frac{\left(\mu\Delta\sqrt{(\tau-\xi)^2 - \mathbf{x}^2/c^2} + b'(\tau-\xi)\right)^2}{16\xi}\right) d\xi$$

$$\times \left(\frac{\left(\mu\Delta\sqrt{(\tau-\xi)^2 - \mathbf{x}^2/c^2} + b'(\tau-\xi)\right)^2}{8\xi} - 1\right) \frac{\mu d\mu}{\sqrt{1-\mu^2}}.$$

In this case, the Green function [18] is given by:

$$G(\mathbf{x},t) = \begin{cases} 0, & \text{if } \quad 0 \le t \le \mathbf{x}/\mathbf{c} \\ \frac{\mathbf{x}}{\mathbf{c}} \left( \frac{\mathbf{b}^{\prime}}{4\sqrt{\pi}} \frac{1}{(t-\mathbf{x}/\mathbf{c})^{3/2}} \exp\left(-\frac{\mathbf{b}^{2}\mathbf{x}^{2}}{16\mathbf{c}^{2}(t-\mathbf{x}/\mathbf{c})}\right) + \Delta \int\_{0}^{t-\mathbf{x}/\mathbf{c}} h^{\prime}(\boldsymbol{\xi}) d\boldsymbol{\xi} \right), & \text{if } \quad t > \mathbf{x}/\mathbf{c}. \end{cases} \tag{23}$$

#### **2.7. Solution of the propagation equation: Green function (Pride-Lafarge model)**

Extensions of Pride-Lafarge (6,8) are expressed in the time domain [3] using the definition of fractional derivative:

$$\begin{split} \tilde{\alpha}(t) &= a\_{\infty} \left( \delta(t) + \frac{2}{\Lambda} \left( \frac{\eta}{\rho\_f} \right)^{1/2} \frac{\partial^{-1/2}}{\partial t^{-1/2}} + \frac{\sigma \phi (1-\wp)}{\rho\_f \mathfrak{a}\_{\infty}} \frac{\partial^{-1}}{\partial t^{-1}} \right), \\ \tilde{\beta}(t) &= \delta(t) + (\gamma - 1) \left( \frac{2}{\Lambda'} \left( \frac{\eta}{P\_r \rho\_f} \right)^{1/2} \frac{\partial^{-1/2}}{\partial t^{-1/2}} + \left( \frac{\eta}{\rho\_f P\_r} \right) \left[ \frac{(1-\wp')\mathfrak{\phi}}{k\_0'} - \frac{4}{\Lambda'^2} \right] \frac{\partial^{-1}}{\partial t^{-1}} \right), \end{split}$$

*α*˜(*t*) and *β*˜(*t*) are the operators of tortuosity and compressibility. In these equations, the operator *∂*<sup>−</sup>1/*∂t* <sup>−</sup><sup>1</sup> represents the time integral:

$$\frac{\partial^{-1}\mathfrak{x}(t)}{\partial t^{-1}} = \int\_0^t \mathfrak{x}(t')dt',$$

and semi-operator *∂*<sup>−</sup>1/2/*∂t* <sup>−</sup>1/2 represents the fractional derivative. In this case, the basic equations are given by:

$$
\rho\_f \mathbf{a}\_{\infty} \,\partial\_t \mathbf{v}(\mathbf{r}, t) + \int\_0^t \chi\_v'(t - t') \,\partial\_t \mathbf{v}(\mathbf{r}, t') \,dt' = -\nabla p(\mathbf{r}, t), \tag{24}
$$

$$\frac{1}{K\_a} \partial\_t p(\mathbf{r}, t) + \int\_0^t \chi\_{th}'(t - t') \,\partial\_t p(\mathbf{r}, t') \, dt' = -\nabla \cdot \mathbf{v}(\mathbf{r}, t), \tag{25}$$

#### 12 Will-be-set-by-IN-TECH 138 Modeling and Measurement Methods for Acoustic Waves and for Acoustic Microdevices Transient Acoustic Wave Propagation in Porous media <sup>13</sup>

where the operators *χ*� *<sup>v</sup>*(*t*) and *χ*� *th*(*t*) are given by:

$$\begin{split} \chi'\_{\upsilon}(t) &= \frac{2\rho\_f \mathfrak{a}\_{\infty}}{\Lambda} \sqrt{\frac{\eta}{\pi \rho\_f}} t^{-1/2} + \sigma \phi (1 - \wp) \frac{\mathfrak{d}^{-1}}{\partial t^{-1}}, \\\chi'\_{th}(t) &= \frac{2(\gamma - 1)}{\mathcal{K}\_a \Lambda'} \sqrt{\frac{\eta}{\pi \mathcal{P}\_r \rho\_f}} t^{-1/2} + \left( \frac{\eta(\gamma - 1)}{\mathcal{K}\_a \rho\_f \mathcal{P}\_r} \right) \left[ \frac{(1 - \wp') \mathfrak{d}}{k'\_0} - \frac{4}{\Lambda'^2} \right] \frac{\partial^{-1}}{\partial t^{-1}}. \end{split}$$

In these expressions, the convolutions express the dispersive nature of the porous material. They take into account the memory effect, where the medium response to excitation of the wave is not instantaneous but takes some time.

The propagation equation can be easily obtained from (24, 25). The propagation equation obtained has exactly the same form as the equation (13), the only difference appears at the coefficient *C* = <sup>4</sup>*α*∞(*γ*−1)*<sup>η</sup> Ka*ΛΛ� <sup>√</sup>*Pr* + *σφ*(1−*℘*) *Ka* <sup>+</sup> *<sup>α</sup>*∞(*γ*−1)*<sup>η</sup> Kak*� <sup>0</sup>*Pr* (1−*℘*� )*φ k*� 0 <sup>−</sup> <sup>4</sup> Λ�<sup>2</sup> , However, we will see that its solution is different. Since the coefficient *C* exchange rate, the coefficients *c*� and Δ given by equation ((17, 18) also change. With extensions of Pride-Lafarge, the coefficient Δ<sup>2</sup> can be negative [3], it gives:

$$
\Delta^{'2} = -\Delta^2 \Longrightarrow \Delta^{'2} = 4c' - b'^2. \tag{26}
$$

where

A = Δ� *μ* 

*<sup>F</sup>*1(*t*) + *<sup>t</sup>*−*x*/*<sup>c</sup>*

**2.8. The reflection and transmission operators**

*pr*

*pt*

The expressions of *R*˜ and *T*˜ are given by (see Ref. [4])

 ∑ *n*≥0

<sup>2</sup> ∑ *n*≥0

> <sup>√</sup>*α*<sup>∞</sup> <sup>−</sup> *<sup>φ</sup>* <sup>√</sup>*α*<sup>∞</sup> <sup>+</sup> *<sup>φ</sup>*

<sup>√</sup>*α*<sup>∞</sup> (*<sup>φ</sup>* <sup>+</sup> <sup>√</sup>*α*∞)<sup>2</sup> *<sup>G</sup>*

the interfaces of the material:

*<sup>R</sup>*˜(*t*) = −*<sup>φ</sup>* <sup>+</sup> <sup>√</sup>*α*<sup>∞</sup>

*<sup>T</sup>*˜(*t*) = <sup>4</sup>*<sup>φ</sup>*

(

are simplified (Ref. [4] ) as follows:

*<sup>φ</sup>* <sup>+</sup> <sup>√</sup>*α*<sup>∞</sup>

<sup>√</sup>*α*<sup>∞</sup>

*R*˜(*t*) =

*<sup>T</sup>*˜(*t*) = <sup>4</sup>*<sup>φ</sup>*

<sup>√</sup>*α*<sup>∞</sup> <sup>+</sup> *<sup>φ</sup>*)

The Green function is given by [3]:

(*t* − *ξ*)

<sup>0</sup> *H*(*ξ*, *t*)*dξ* + *j*

the material. These are operators [4] represented in integral form:

0

0

*R*˜(*τ*)*p<sup>i</sup>*

*T*˜(*τ*)*p<sup>i</sup>*

*<sup>φ</sup>* <sup>−</sup> <sup>√</sup>*α*<sup>∞</sup> *<sup>φ</sup>* <sup>+</sup> <sup>√</sup>*α*<sup>∞</sup>

 *<sup>φ</sup>* <sup>−</sup> <sup>√</sup>*α*<sup>∞</sup> *<sup>φ</sup>* <sup>+</sup> <sup>√</sup>*α*<sup>∞</sup>

(*x*, *<sup>t</sup>*) = *<sup>t</sup>*

(*x*, *<sup>t</sup>*) = *<sup>t</sup>*

<sup>2</sup> <sup>−</sup> *<sup>x</sup>*2/*c*<sup>2</sup> and <sup>B</sup> <sup>=</sup> *<sup>b</sup>*�

*<sup>G</sup>*(*t*, *<sup>x</sup>*) = 0, if 0 <sup>≤</sup> *<sup>t</sup>* <sup>≤</sup> *<sup>x</sup>*/*c*,

*<sup>t</sup>*−*x*/*<sup>c</sup>*

For a slice of porous material occupying the domain 0 ≤ *x* ≤ *L*, the incident and scattered fields are connected by scattering operators (i.e. operators of reflection and transmission) of

In equations (30) and (31) functions *R*˜ and *T*˜ represent the kernel of reflection and transmission operators, respectively. These are independently-owned operators of the incident field and depend only on the material properties. To express the operators of reflection and transmission, we assume [4] that the field of acoustic pressure and flow are continuous at

*p*(0+, *t*) = *p*(0−, *t*), *p*(*L*−, *t*) = *p*(*L*+, *t*)

2*<sup>n</sup> G t* + *L c*0

where *G* is the Green function of the medium. These expressions take into account the multiple *n* reflections within the porous material. Given the high attenuation of acoustic waves in air-saturated porous media, multiple reflections are negligible. For reflections at the interfaces *x* = 0 and *x* = *L*, the expressions of the operators of reflection and transmission

<sup>√</sup>*α*∞(

(

<sup>√</sup>*α*<sup>∞</sup> <sup>−</sup> *<sup>φ</sup>*)

<sup>√</sup>*α*<sup>∞</sup> <sup>+</sup> *<sup>φ</sup>*)<sup>3</sup> *<sup>G</sup>*

*<sup>δ</sup>*(*t*) <sup>−</sup> <sup>4</sup>*<sup>φ</sup>*

 *t* + *L c* , *L c* 

2*<sup>n</sup> G t*, 2*n L c* − *G* 

*t* − *τ* +

*<sup>t</sup>* <sup>−</sup> *<sup>τ</sup>* <sup>−</sup> *<sup>L</sup>*

*x c*0 

*<sup>c</sup>* <sup>−</sup> (*<sup>x</sup>* <sup>−</sup> *<sup>L</sup>*) *c*0

*v*(0−, *t*) = *φv*(0+, *t*), *v*(*L*+, *t*) = *φv*(*L*−, *t*) (32)

,(2*n* + 1)

*L c* .

 *t*, 2*L c* 

(*t* − *ξ*).

Transient Acoustic Wave Propagation in Porous Media 139

<sup>0</sup> <sup>Ξ</sup>(*ξ*, *<sup>t</sup>*)*d<sup>ξ</sup>* if *<sup>t</sup> <sup>&</sup>gt; <sup>x</sup>*/*c*. (29)

*dτ*, (30)

*t*,(2*n* + 2)

. (34)

*L c* ,

, (33)

*dτ*. (31)

In this case, *f*(*z*) given by equation (16) can be written in the form:

$$f(z) = \left(z + \frac{b'}{2}\sqrt{z}\right)^2 + \left(\frac{\Delta'\sqrt{z}}{2}\right)^2. \tag{27}$$

Using expressions (26 and 27) and analytical calculations given in [3], we obtain the solution of the wave equation when Δ<sup>2</sup> is negative:

$$p(\mathbf{x},t) = \begin{cases} 0, & \text{if } 0 \le t \le \mathbf{x}/c\\ \int\_{\mathbf{x}/c}^{t} \left[ F\_{\mathbf{f}}(\mathbf{r}) + \int\_{0}^{\mathbf{r}-\mathbf{x}/c} H(\boldsymbol{\xi},\mathbf{r}) d\boldsymbol{\xi} + j \int\_{0}^{\mathbf{r}-\mathbf{x}/c} \Xi(\boldsymbol{\xi},\mathbf{r}) d\boldsymbol{\xi} \right] p(0,t-\mathbf{r}) d\mathbf{r} & \text{if } \quad t > \mathbf{x}/c \end{cases} (28)$$

with

$$F\_1(\tau) = \frac{1}{4\sqrt{\pi}} b^{\prime} \frac{\chi}{c} \frac{1}{(\tau - \mathfrak{x}/c)^{3/2}} \exp\left(-\frac{b^{\prime 2} \mathfrak{x}^2}{16c^2(\tau - \mathfrak{x}/c)}\right), \quad \dot{f}^2 = -1.1$$

The functions *H*(*ξ*, *τ*) and Ξ(*ξ*, *τ*) are given by the following equations where *τ* is replaced by *t*:

$$\begin{split} H(\xi,t) &= \frac{\Delta'}{4\pi\sqrt{\pi}} \frac{x}{c} \frac{1}{\sqrt{(t-\xi)^2 - x^2/c^2}} \frac{1}{\xi^{3/2}} \times \\ &\int\_{-1}^{1} \left[\frac{\mathcal{A}\mathcal{B}}{4\xi} \cos\left(\frac{\mathcal{A}\mathcal{B}}{8\xi}\right) - \left(\frac{\mathcal{B}^2 - \mathcal{A}^2}{8\xi} - 1\right) \sin\left(\frac{\mathcal{A}\mathcal{B}}{8\xi}\right)\right] \exp\left(-\frac{\mathcal{B}^2 - \mathcal{A}^2}{16\xi}\right) \frac{\mu d\mu}{\sqrt{1-\mu^2}} \end{split}$$

and

$$\begin{split} \Xi(\xi,t) &= -\frac{\Delta'}{4\pi\sqrt{\pi}} \frac{\mathbf{x}}{c} \frac{1}{\sqrt{(t-\xi)^2 - \mathbf{x}^2/c^2}} \frac{1}{\xi^{3/2}} \\ &\int\_{-1}^{1} \left[ \left(\frac{\mathcal{B}^2 - \mathcal{A}^2}{8\xi} - 1\right) \cos\left(\frac{\mathcal{A}\mathcal{B}}{8\xi}\right) + \frac{\mathcal{A}\mathcal{B}}{4\xi} \sin\left(\frac{\mathcal{A}\mathcal{B}}{8\xi}\right) \right] \exp\left(-\frac{\mathcal{B}^2 - \mathcal{A}^2}{16\xi}\right) \frac{\mu d\mu}{\sqrt{1-\mu^2}}, \end{split}$$

where

12 Will-be-set-by-IN-TECH

<sup>−</sup>1/2 <sup>+</sup> *σφ*(<sup>1</sup> <sup>−</sup> *<sup>℘</sup>*) *<sup>∂</sup>*−<sup>1</sup>

In these expressions, the convolutions express the dispersive nature of the porous material. They take into account the memory effect, where the medium response to excitation of the

The propagation equation can be easily obtained from (24, 25). The propagation equation obtained has exactly the same form as the equation (13), the only difference appears at the

its solution is different. Since the coefficient *C* exchange rate, the coefficients *c*� and Δ given by equation ((17, 18) also change. With extensions of Pride-Lafarge, the coefficient Δ<sup>2</sup> can be

Using expressions (26 and 27) and analytical calculations given in [3], we obtain the solution

0, if 0 <sup>≤</sup> *<sup>t</sup>* <sup>≤</sup> *<sup>x</sup>*/*c*, *<sup>t</sup>*

*<sup>τ</sup>*−*x*/*<sup>c</sup>*

The functions *H*(*ξ*, *τ*) and Ξ(*ξ*, *τ*) are given by the following equations where *τ* is replaced by

1 *<sup>ξ</sup>*3/2 <sup>×</sup>

> sin

1 *ξ*3/2

<sup>0</sup> Ξ(*ξ*, *τ*)*dξ*

<sup>−</sup> *<sup>b</sup>*�2*x*<sup>2</sup> <sup>16</sup>*c*2(*<sup>τ</sup>* − *<sup>x</sup>*/*c*)

> A B 8*ξ*

A B 8*ξ*

 exp 

 exp 

*Ka* <sup>+</sup> *<sup>α</sup>*∞(*γ*−1)*<sup>η</sup> Kak*� <sup>0</sup>*Pr*

<sup>2</sup> <sup>=</sup> <sup>−</sup>Δ<sup>2</sup> <sup>=</sup><sup>⇒</sup> <sup>Δ</sup>�

*<sup>∂</sup>t*−<sup>1</sup> ,

<sup>−</sup> <sup>4</sup> Λ�<sup>2</sup> 

2

 , *j*

(1 − *℘*�

*k*� 0 )*φ*

<sup>2</sup> <sup>=</sup> <sup>4</sup>*c*� <sup>−</sup> *<sup>b</sup>*�2. (26)

<sup>−</sup> <sup>4</sup> Λ�<sup>2</sup>  *∂*−<sup>1</sup> *<sup>∂</sup>t*−<sup>1</sup> ,

, However, we will see that

. (27)

*p*(0, *t*−*τ*)*dτ* if *t > x*/*c*.

<sup>2</sup> <sup>=</sup> <sup>−</sup>1.

 *μdμ* <sup>1</sup> <sup>−</sup> *<sup>μ</sup>*<sup>2</sup>

 *μdμ* <sup>1</sup> <sup>−</sup> *<sup>μ</sup>*<sup>2</sup>

,

<sup>−</sup> <sup>B</sup><sup>2</sup> − A<sup>2</sup> 16*ξ*

<sup>−</sup> <sup>B</sup><sup>2</sup> − A<sup>2</sup> 16*ξ*

(28)

*η*(*γ* − 1) *Kaρ <sup>f</sup> Pr*

> (1−*℘*� )*φ k*� 0

*th*(*t*) are given by:

where the operators *χ*�

*χ*� *<sup>v</sup>*(*t*) =

*χ*�

coefficient *C* =

*p*(*x*, *t*)=

with

*t*:

and

negative [3], it gives:

*<sup>v</sup>*(*t*) and *χ*�

 *η πρ <sup>f</sup> t*

> *η πPrρ <sup>f</sup> t* <sup>−</sup>1/2 +

+ *σφ*(1−*℘*)

Δ�

*f*(*z*) =

In this case, *f*(*z*) given by equation (16) can be written in the form:

 *z* + *b*� 2 √*z* 2 + Δ� √*z* 2

<sup>0</sup> *H*(*ξ*, *τ*)*dξ*+*j*

1 (*<sup>τ</sup>* <sup>−</sup> *<sup>x</sup>*/*c*)3/2 exp

1 (*<sup>t</sup>* <sup>−</sup> *<sup>ξ</sup>*)<sup>2</sup> <sup>−</sup> *<sup>x</sup>*2/*c*<sup>2</sup>

<sup>B</sup><sup>2</sup> − A<sup>2</sup>

1 (*<sup>t</sup>* <sup>−</sup> *<sup>ξ</sup>*)<sup>2</sup> <sup>−</sup> *<sup>x</sup>*2/*c*<sup>2</sup>

> + A B <sup>4</sup>*<sup>ξ</sup>* sin

A B 8*ξ*

<sup>8</sup>*<sup>ξ</sup>* <sup>−</sup> <sup>1</sup>

2*ρ <sup>f</sup> α*<sup>∞</sup> Λ

wave is not instantaneous but takes some time.

 <sup>4</sup>*α*∞(*γ*−1)*<sup>η</sup> Ka*ΛΛ� <sup>√</sup>*Pr*

of the wave equation when Δ<sup>2</sup> is negative:

*<sup>F</sup>*1(*τ*)+ *<sup>τ</sup>*−*x*/*<sup>c</sup>*

*x*/*c* 

*<sup>H</sup>*(*ξ*, *<sup>t</sup>*) = <sup>Δ</sup>�

<sup>Ξ</sup>(*ξ*, *<sup>t</sup>*) = <sup>−</sup> <sup>Δ</sup>�

<sup>B</sup><sup>2</sup> − A<sup>2</sup>

 1 −1

 1 −1  A B <sup>4</sup>*<sup>ξ</sup>* cos

*<sup>F</sup>*1(*τ*) = <sup>1</sup> 4 <sup>√</sup>*<sup>π</sup> b*� *x c*

> 4*π* <sup>√</sup>*<sup>π</sup> x c*

> > 4*π* <sup>√</sup>*<sup>π</sup> x c*

<sup>8</sup>*<sup>ξ</sup>* <sup>−</sup> <sup>1</sup>

A B 8*ξ*

 −

 cos

*th*(*t*) = <sup>2</sup>(*<sup>γ</sup>* <sup>−</sup> <sup>1</sup>) *Ka*Λ�

$$\mathcal{A} = \Delta' \mu \sqrt{(t - \xi)^2 - \mathfrak{x}^2/c^2} \qquad \text{and} \qquad \mathcal{B} = b'(t - \xi).$$

The Green function is given by [3]:

$$G(t, \mathbf{x}) = \begin{cases} 0, & \text{if } \quad 0 \le t \le \mathbf{x}/c \\ F\_1(t) + \int\_0^{t-\mathbf{x}/c} H(\xi, t) d\xi + j \int\_0^{t-\mathbf{x}/c} \Xi(\xi, t) d\xi & \text{if } \qquad t > \mathbf{x}/c. \end{cases} \tag{29}$$

#### **2.8. The reflection and transmission operators**

For a slice of porous material occupying the domain 0 ≤ *x* ≤ *L*, the incident and scattered fields are connected by scattering operators (i.e. operators of reflection and transmission) of the material. These are operators [4] represented in integral form:

$$p^r(\mathbf{x}, t) = \int\_0^t \tilde{\mathcal{R}}(\mathbf{r}) p^i \left( t - \mathbf{r} + \frac{\mathbf{x}}{c\_0} \right) d\mathbf{r},\tag{30}$$

$$p^t(\mathbf{x}, t) = \int\_0^t \Upsilon(\tau) p^i \left( t - \tau - \frac{L}{c} - \frac{(\mathbf{x} - L)}{c\_0} \right) d\tau. \tag{31}$$

In equations (30) and (31) functions *R*˜ and *T*˜ represent the kernel of reflection and transmission operators, respectively. These are independently-owned operators of the incident field and depend only on the material properties. To express the operators of reflection and transmission, we assume [4] that the field of acoustic pressure and flow are continuous at the interfaces of the material:

$$\begin{aligned} p(0^+,t) &= p(0^-,t), & p(L^-,t) &= p(L^+,t) \\ v(0^-,t) &= \phi v(0^+,t), & v(L^+,t) &= \phi v(L^-,t) \end{aligned} \tag{32}$$

The expressions of *R*˜ and *T*˜ are given by (see Ref. [4])

$$
\begin{split}
\tilde{\mathcal{R}}(t) &= \left(\frac{-\phi + \sqrt{\kappa\_{\infty}}}{\phi + \sqrt{\kappa\_{\infty}}}\right) \sum\_{n\geq 0} \left(\frac{\phi - \sqrt{\kappa\_{\infty}}}{\phi + \sqrt{\kappa\_{\infty}}}\right)^{2n} \left[G\left(t, 2n\frac{L}{c}\right) - G\left(t, (2n+2)\frac{L}{c}\right)\right], \\
\tilde{\mathcal{T}}(t) &= \frac{4\phi\sqrt{\kappa\_{\infty}}}{\left(\sqrt{\kappa\_{\infty}} + \phi\right)^{2}} \sum\_{n\geq 0} \left(\frac{\phi - \sqrt{\kappa\_{\infty}}}{\phi + \sqrt{\kappa\_{\infty}}}\right)^{2n} G\left(t + \frac{L}{c\_{0}}, (2n+1)\frac{L}{c}\right).
\end{split}
$$

where *G* is the Green function of the medium. These expressions take into account the multiple *n* reflections within the porous material. Given the high attenuation of acoustic waves in air-saturated porous media, multiple reflections are negligible. For reflections at the interfaces *x* = 0 and *x* = *L*, the expressions of the operators of reflection and transmission are simplified (Ref. [4] ) as follows:

$$\tilde{\mathcal{R}}(t) = \frac{\sqrt{\alpha\_{\infty}} - \phi}{\sqrt{\alpha\_{\infty}} + \phi} \delta(t) - \frac{4\phi\sqrt{\alpha\_{\infty}}(\sqrt{\alpha\_{\infty}} - \phi)}{(\sqrt{\alpha\_{\infty}} + \phi)^3} G\left(t, \frac{2L}{c}\right),\tag{33}$$

$$\hat{T}(t) = \frac{4\phi\sqrt{\kappa\_{\infty}}}{(\phi + \sqrt{\kappa\_{\infty}})^2} G\left(t + \frac{L}{c}, \frac{L}{c}\right). \tag{34}$$

#### 14 Will-be-set-by-IN-TECH 140 Modeling and Measurement Methods for Acoustic Waves and for Acoustic Microdevices Transient Acoustic Wave Propagation in Porous media <sup>15</sup>

The first term on the second member of the equation (33):

(( <sup>√</sup>*α*<sup>∞</sup> <sup>−</sup> *<sup>φ</sup>*)/( <sup>√</sup>*α*<sup>∞</sup> <sup>+</sup> *<sup>φ</sup>*)) *<sup>δ</sup>*(*t*) is equivalent to the instantaneous reflection response of the porous material. This term corresponds to the wave reflected by the first interface x = 0 . It depends only on the porosity and tortuosity of the material. The reflected wave to the first interface has the advantage not to be dispersive, but simply attenuated. This shows that it is possible to measure the porosity and tortuosity of the porous material by measuring just the first reflected wave.

The second term of equation (33): <sup>−</sup>4*<sup>φ</sup>* <sup>√</sup>*α*∞( <sup>√</sup>*α*∞−*φ*) ( <sup>√</sup>*α*∞+*φ*)<sup>3</sup> *<sup>G</sup> t*, <sup>2</sup>*<sup>L</sup> c* corresponds to the reflection by the second interface *x* = *L*. This term depends on the Green function of the medium that describes the propagation and scattering of the acoustic wave having made a round trip in the slab of porous material. Green's function depends on the tortuosity, and viscous and thermal characteristic lengths Λ and Λ� material, but does not depend on the porosity. Experimentally, this second contribution to the debate can not be measured for low-resistive materials, because the acoustic signal is very attenuated.

1 2 3 4

Time (s)

0 1 2 3 4

Frequency (Hz)

1 2 3 4

Time (s)

sample M2 less resistive, having the following parameters: thickness 0.8 cm, tortuosity *α*<sup>∞</sup> = 1.05, viscous characteristic length Λ = 300 *μ* m, thermal characteristic length Λ 900 = *μ* m,

**Figure 5.** Comparison between simulated transmitted signal corresponding to the real part of the solution of Pride-Lafarge (Eq. 28) for *℘* = *℘* = 0.7 (solid line) and simulated signal corresponding to the

solution of Johnson-Allard (Eq. 22) for *℘* = *℘* = 1 (dotted line) for the sample M1.

−0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8

−15

−0.02 −0.015 −0.01 −0.005 0 0.005 0.01 0.015 0.02

Amplitude (a.u.)

**Figure 4.** Spectrum of the incident signal

−10

Amplitude (dB)

−5

0

**Figure 3.** Incident signal

Amplitude (a.u.)

x 10 −5

Transient Acoustic Wave Propagation in Porous Media 141

x 10 5

x 10−5

Let us study the sensitivity [20] of porosity on the transmission operator *T*˜(*t*). Taking the derivative of *T*˜(*t*) versus the porosity *φ*, we obtain:

$$\frac{\partial \tilde{T}}{\partial \phi} = \frac{4\sqrt{\alpha\_{\infty}}(\sqrt{\alpha\_{\infty}} - \phi)}{\left(\sqrt{\alpha\_{\infty}} + \phi\right)^3} G\left(t + \frac{L}{c}, \frac{L}{c}\right),$$

when *<sup>φ</sup>* −→ <sup>√</sup>*α*∞, the derivative *<sup>∂</sup>T*˜ *∂φ* −→ 0. This shows that when the porous medium tends to a free fluid, or for weakly resistive porous material with low values of porosity and tortuosity (near 1), the sensitivity of the porosity of the transmitted wave tends to zero. More generally, for a wide range of air saturated porous materials, the term *<sup>∂</sup>T*˜ *∂φ* is very small. Finally, we conclude that the operator of transmission depends on all parameters, but the low sensitivity of the porosity, makes impossible the determination of this parameter from transmitted data.

#### **2.9. Numerical simulations**

Let us consider a sample of air saturated porous M1 having the following physical properties: thickness 0.8 cm, tortuosity *α*<sup>∞</sup> = 1.5, viscous characteristic length Λ = 25 *μ* m, thermal characteristic length Λ = 75*μ* m, specific resistance to fluid flow *σ* = 200,000 Nm <sup>−</sup><sup>4</sup> s, porosity *φ* = 0.82 and thermal permeability *k*� <sup>0</sup> <sup>=</sup> 2.77 10−<sup>10</sup> <sup>m</sup> 2. A simulated transmitted signal was calculated from the equation (31). The input signal is given in figure 3 and its spectrum in Figure 4. Figure 5 shows a comparison between two signals, the first (solid line) corresponds to the real part of the solution (28) of Pride-Lafarge model when *℘* = *℘*� = 0.7 and the second (dotted) corresponds to the solution (22) model of Johnson-Allard (*℘* = *℘*� = 1). We choose the same values for *℘* and *℘*� to simplify the study. This specific value of 0.7 is only valid for porous materials with circular pores. In the general case, the values of *℘* and *℘*� may be different from 0.7. Note that it is possible to have positive values of Δ (Eq. 18) for other values of *℘* and *℘*� . From Figure 5, we notice a significant change in the wave amplitude. By increasing *℘* and *℘*� from 0.7 to 1, the wave amplitude increases by 65% of its original value. This result can be predicted by the fact that when the values of *℘* and *℘*� increases, the coefficient *C* decreases, and hence the wave amplitude decreases due to interactions; inertial, viscous and thermal between fluid and structure. This phenomenon is much more important for resistive porous materials. Figure 6 shows the same comparison as Figure 5 for another

**Figure 3.** Incident signal

14 Will-be-set-by-IN-TECH

porous material. This term corresponds to the wave reflected by the first interface x = 0 . It depends only on the porosity and tortuosity of the material. The reflected wave to the first interface has the advantage not to be dispersive, but simply attenuated. This shows that it is possible to measure the porosity and tortuosity of the porous material by measuring just the

<sup>√</sup>*α*∞−*φ*)

<sup>√</sup>*α*∞+*φ*)<sup>3</sup> *<sup>G</sup>*

the second interface *x* = *L*. This term depends on the Green function of the medium that describes the propagation and scattering of the acoustic wave having made a round trip in the slab of porous material. Green's function depends on the tortuosity, and viscous and thermal characteristic lengths Λ and Λ� material, but does not depend on the porosity. Experimentally, this second contribution to the debate can not be measured for low-resistive materials, because

Let us study the sensitivity [20] of porosity on the transmission operator *T*˜(*t*). Taking the

<sup>√</sup>*α*<sup>∞</sup> <sup>−</sup> *<sup>φ</sup>*)

a free fluid, or for weakly resistive porous material with low values of porosity and tortuosity (near 1), the sensitivity of the porosity of the transmitted wave tends to zero. More generally,

conclude that the operator of transmission depends on all parameters, but the low sensitivity of the porosity, makes impossible the determination of this parameter from transmitted data.

Let us consider a sample of air saturated porous M1 having the following physical properties: thickness 0.8 cm, tortuosity *α*<sup>∞</sup> = 1.5, viscous characteristic length Λ = 25 *μ* m, thermal characteristic length Λ = 75*μ* m, specific resistance to fluid flow *σ* = 200,000 Nm <sup>−</sup><sup>4</sup> s, porosity

calculated from the equation (31). The input signal is given in figure 3 and its spectrum in Figure 4. Figure 5 shows a comparison between two signals, the first (solid line) corresponds to the real part of the solution (28) of Pride-Lafarge model when *℘* = *℘*� = 0.7 and the second (dotted) corresponds to the solution (22) model of Johnson-Allard (*℘* = *℘*� = 1). We choose the same values for *℘* and *℘*� to simplify the study. This specific value of 0.7 is only valid for porous materials with circular pores. In the general case, the values of *℘* and *℘*� may be different from 0.7. Note that it is possible to have positive values of Δ (Eq. 18) for other

By increasing *℘* and *℘*� from 0.7 to 1, the wave amplitude increases by 65% of its original value. This result can be predicted by the fact that when the values of *℘* and *℘*� increases, the coefficient *C* decreases, and hence the wave amplitude decreases due to interactions; inertial, viscous and thermal between fluid and structure. This phenomenon is much more important for resistive porous materials. Figure 6 shows the same comparison as Figure 5 for another

. From Figure 5, we notice a significant change in the wave amplitude.

<sup>3</sup> *G*

 *t* + *L c* , *L c* ,

*∂φ* −→ 0. This shows that when the porous medium tends to

<sup>0</sup> <sup>=</sup> 2.77 10−<sup>10</sup> <sup>m</sup> 2. A simulated transmitted signal was

<sup>√</sup>*α*<sup>∞</sup> <sup>+</sup> *<sup>φ</sup>*)

<sup>√</sup>*α*∞(

(

<sup>√</sup>*α*∞(

(

for a wide range of air saturated porous materials, the term *<sup>∂</sup>T*˜

<sup>√</sup>*α*<sup>∞</sup> <sup>+</sup> *<sup>φ</sup>*)) *<sup>δ</sup>*(*t*) is equivalent to the instantaneous reflection response of the

 *t*, <sup>2</sup>*<sup>L</sup> c* 

corresponds to the reflection by

*∂φ* is very small. Finally, we

The first term on the second member of the equation (33):

((

<sup>√</sup>*α*<sup>∞</sup> <sup>−</sup> *<sup>φ</sup>*)/(

first reflected wave.

The second term of equation (33): <sup>−</sup>4*<sup>φ</sup>*

the acoustic signal is very attenuated.

when *<sup>φ</sup>* −→ <sup>√</sup>*α*∞, the derivative *<sup>∂</sup>T*˜

**2.9. Numerical simulations**

values of *℘* and *℘*�

*φ* = 0.82 and thermal permeability *k*�

derivative of *T*˜(*t*) versus the porosity *φ*, we obtain:

*∂T*˜ *∂φ* <sup>=</sup> <sup>4</sup>

**Figure 4.** Spectrum of the incident signal

**Figure 5.** Comparison between simulated transmitted signal corresponding to the real part of the solution of Pride-Lafarge (Eq. 28) for *℘* = *℘* = 0.7 (solid line) and simulated signal corresponding to the solution of Johnson-Allard (Eq. 22) for *℘* = *℘* = 1 (dotted line) for the sample M1.

sample M2 less resistive, having the following parameters: thickness 0.8 cm, tortuosity *α*<sup>∞</sup> = 1.05, viscous characteristic length Λ = 300 *μ* m, thermal characteristic length Λ 900 = *μ* m,

**Figure 6.** Comparison between simulated transmitted signal corresponding to the real part of the solution of Pride-Lafarge (Eq. 28) for *℘* = *℘*� = 0.7 (solid line), and the transmitted signal corresponding to the simulated Allard-Johnson solution (Eq. 22) for *℘* = *℘*� = 1 (dotted line) for the sample M2.

1 2 3 4 5

Time (s)

**Figure 8.** Simulated transmitted signal corresponding to the imaginary part of the solution (28) model of

tan *θ* = *G*2/*G*<sup>1</sup> � 1. This leads to the conclusion that all components of the signal have the

The expression (34) shows that the transmitted wave depends on the porosity *φ*, and parameters describing the high frequency acoustic propagation in the material (tortuosity *α*∞

porosity (in the previous paragraph) showed that the effect of propagation described by the Green's function altogether mask the effect of interfaces: recall that the Green function *G* is independent of porosity and the effect of the interface appears through the term <sup>4</sup>*<sup>φ</sup>*

in the expression of the transmission operator (34). The low sensitivity of the porosity in transmission makes its determination impossible. Solving the inverse problem, using the transmitted waves allowed the determination of the tortuosity and viscous and thermal characteristic lengths. Note that the ratio of two lengths was set at 3 [20–23]. Currently the only method for the separate determination of the two lengths is based on the saturation of the porous by two different fluids [24]. The two characteristic lengths are acting the same way on the dispersion and attenuation of the transmitted signal. So it is mathematically impossible to trace two unknowns giving the same physical effect on the waveform. We limit ourselves in our work to set an arbitrary ratio between the two lengths. The tortuosity acts on the speed of the wavefront (delay of the transmitted signal with respect to the incident signal) and the attenuation of the amplitude of the waveform. While the characteristic lengths affect the dispersion and attenuation of the signal. A basic inverse problem associated with a slab of porous material may be stated as follows: from the measurements of the transmitted signals outside the slab, find the values of the parameters of the medium. As shown in the previous section, the solution of the direct problem is the system of two operators expressed as

of the slab is based on a fitting procedure: find the values of the parameters *α*∞, Λ and Λ� such that the transmitted and reflected signal describes the scattering problem in the best possible way (e.g., in the least-squares sense). Experiments were performed in air using a pairs of broadband Ultran transducers NCG200-D13 with a central frequency of 195 kHz and

. The inversion algorithm for finding the values of the parameters

same phase (very small *θ*, which is therefore a physical factor not essential.

−4 −3 −2 −1 0 1 2 3 4 5 x 10−10

and viscous and thermal characteristic lengths Λ and Λ�

Amplitude(a.u.)

Pride-Lafarge.

**2.10. Transmitted wave**

functions on *α*∞, Λ and Λ�

x 10−5

Transient Acoustic Wave Propagation in Porous Media 143

). The study of the sensitivity of the

<sup>√</sup>*α*<sup>∞</sup> (*φ*+√*α*∞)<sup>2</sup>

**Figure 7.** Spectrum of the transmitted signal simulated for *℘* = *℘*� = 0.7 Eq. (28) (solid line) and simulated transmitted signal spectrum for *℘* = *℘*� = 1 (Eq. 22) (poitillés).

specific resistance to the passage of fluid *<sup>σ</sup>* <sup>=</sup>20000*Nm*−4 s, porosity *<sup>φ</sup>* <sup>=</sup> 0.96 and thermal permeability *k*� <sup>0</sup> <sup>=</sup> 2.77 10−<sup>9</sup> <sup>m</sup> 2. In Figure 6, the influence of the parameters *<sup>℘</sup>* and *<sup>℘</sup>*� on the attenuation is smaller than that in Figure 5. We can conclude that the parameters *℘* and *℘*� play an important role in the acoustic attenuation, especially for resistive media.

It is possible to see from Figures 5 and 6, when the values of *℘* and *℘*� change from 0.7 to 1, the waveform changes only at the amplitude of the wave, but it is not dispersed. Figure 7 shows the spectra of two simulated signals given in Figure 5. From the spectra of two simulated signals, we can see they have the same bandwidth, which means that there is no dispersion. This last result shows that *℘* and *℘*� play an important role in attenuation the acoustic wave but not on its dispersion. Figure 8 shows the imaginary part of the solution (28) for *℘* = *℘*� = 0.7 (sample M1). The amplitude of the imaginary part of the solution is very small compared to the real part (figure 5), this is why only the real part of the solution corresponding to the physical solution is taken into account when a comparison with the experiment is performed. It is possible to write the Green function given by equation (29) as *<sup>G</sup>* <sup>=</sup> *<sup>G</sup>*<sup>1</sup> <sup>+</sup> *iG*<sup>2</sup> <sup>=</sup><sup>|</sup> *<sup>G</sup>* <sup>|</sup>e*iθ*,

**Figure 8.** Simulated transmitted signal corresponding to the imaginary part of the solution (28) model of Pride-Lafarge.

tan *θ* = *G*2/*G*<sup>1</sup> � 1. This leads to the conclusion that all components of the signal have the same phase (very small *θ*, which is therefore a physical factor not essential.

#### **2.10. Transmitted wave**

16 Will-be-set-by-IN-TECH

1 2 3 4

Time (s)

0 0.5 1 1.5 2 2.5

specific resistance to the passage of fluid *<sup>σ</sup>* <sup>=</sup>20000*Nm*−4 s, porosity *<sup>φ</sup>* <sup>=</sup> 0.96 and thermal

the attenuation is smaller than that in Figure 5. We can conclude that the parameters *℘* and *℘*�

It is possible to see from Figures 5 and 6, when the values of *℘* and *℘*� change from 0.7 to 1, the waveform changes only at the amplitude of the wave, but it is not dispersed. Figure 7 shows the spectra of two simulated signals given in Figure 5. From the spectra of two simulated signals, we can see they have the same bandwidth, which means that there is no dispersion. This last result shows that *℘* and *℘*� play an important role in attenuation the acoustic wave but not on its dispersion. Figure 8 shows the imaginary part of the solution (28) for *℘* = *℘*� = 0.7 (sample M1). The amplitude of the imaginary part of the solution is very small compared to the real part (figure 5), this is why only the real part of the solution corresponding to the physical solution is taken into account when a comparison with the experiment is performed. It is possible to write the Green function given by equation (29) as *<sup>G</sup>* <sup>=</sup> *<sup>G</sup>*<sup>1</sup> <sup>+</sup> *iG*<sup>2</sup> <sup>=</sup><sup>|</sup> *<sup>G</sup>* <sup>|</sup>e*iθ*,

**Figure 7.** Spectrum of the transmitted signal simulated for *℘* = *℘*� = 0.7 Eq. (28) (solid line) and

play an important role in the acoustic attenuation, especially for resistive media.

Frequency (Hz)

<sup>0</sup> <sup>=</sup> 2.77 10−<sup>9</sup> <sup>m</sup> 2. In Figure 6, the influence of the parameters *<sup>℘</sup>* and *<sup>℘</sup>*� on

**Figure 6.** Comparison between simulated transmitted signal corresponding to the real part of the solution of Pride-Lafarge (Eq. 28) for *℘* = *℘*� = 0.7 (solid line), and the transmitted signal corresponding to the simulated Allard-Johnson solution (Eq. 22) for *℘* = *℘*� = 1 (dotted line) for the sample M2.

−0.4

−20

simulated transmitted signal spectrum for *℘* = *℘*� = 1 (Eq. 22) (poitillés).

−15

−10

Amplitude (dB)

permeability *k*�

−5

0

−0.2

0

0.2

Amplitude (a.u.)

0.4

0.6

x 10−5

x 105

The expression (34) shows that the transmitted wave depends on the porosity *φ*, and parameters describing the high frequency acoustic propagation in the material (tortuosity *α*∞ and viscous and thermal characteristic lengths Λ and Λ� ). The study of the sensitivity of the porosity (in the previous paragraph) showed that the effect of propagation described by the Green's function altogether mask the effect of interfaces: recall that the Green function *G* is independent of porosity and the effect of the interface appears through the term <sup>4</sup>*<sup>φ</sup>* <sup>√</sup>*α*<sup>∞</sup> (*φ*+√*α*∞)<sup>2</sup> in the expression of the transmission operator (34). The low sensitivity of the porosity in transmission makes its determination impossible. Solving the inverse problem, using the transmitted waves allowed the determination of the tortuosity and viscous and thermal characteristic lengths. Note that the ratio of two lengths was set at 3 [20–23]. Currently the only method for the separate determination of the two lengths is based on the saturation of the porous by two different fluids [24]. The two characteristic lengths are acting the same way on the dispersion and attenuation of the transmitted signal. So it is mathematically impossible to trace two unknowns giving the same physical effect on the waveform. We limit ourselves in our work to set an arbitrary ratio between the two lengths. The tortuosity acts on the speed of the wavefront (delay of the transmitted signal with respect to the incident signal) and the attenuation of the amplitude of the waveform. While the characteristic lengths affect the dispersion and attenuation of the signal. A basic inverse problem associated with a slab of porous material may be stated as follows: from the measurements of the transmitted signals outside the slab, find the values of the parameters of the medium. As shown in the previous section, the solution of the direct problem is the system of two operators expressed as functions on *α*∞, Λ and Λ� . The inversion algorithm for finding the values of the parameters of the slab is based on a fitting procedure: find the values of the parameters *α*∞, Λ and Λ� such that the transmitted and reflected signal describes the scattering problem in the best possible way (e.g., in the least-squares sense). Experiments were performed in air using a pairs of broadband Ultran transducers NCG200-D13 with a central frequency of 195 kHz and a bandwith of 6 dB extending from 140 kHz to 250 kHz. Pulses of 400 V are provided by a 5058PR Panametrics pulser/receiver. The received signals are filtered above 1 MHz to avoid high frequency noise. Electronic interference is eliminated by 1000 acquisition averages. The experimental setup is shown in Fig. 9. Consider a sample of plastic foam, of thicknesses

6.6 6.8 7 7.2 7.4

Time (s)

**Figure 11.** Comparison between experimental transmitted signal (solid line) and simulated transmitted

[25]. This method is based on a temporal model of direct and inverse scattering problems affecting the propagation of transient ultrasonic waves in a homogeneous isotropic slab of porous material with a rigid frame. Generally, porosity and tortuosity can be evaluated simultaneously by the wave reflected at the first interface at oblique incidence [26–29] but this is not possible at normal incidence [4, 30]. Viscous and thermal characteristic lengths can be estimated only by transmitted waves. Porosity cannot be determined from transmitted waves because of its weak sensitivity in this mode [4]. The advantage of the proposed method is that all the parameters can be determined at normal incidence (the ratio between viscous and thermal lengths is fixed as in classical methods [21–24] based on transmission measurement). The main principle of this method is the experimental detection of reflected contributions from the first and the second interface of the medium. The properties of these two contributions are used to estimate the four acoustical parameters needed for ultrasonic propagation in porous material with a rigid frame, by solving the inverse problem. Studying the sensitivity of each reflected wave parameter demonstrates the importance of each contribution) first and second interface! for the inversion. Numerical and experimental validation for weak resistive

This section concerns the ultrasonic characterization of air-saturated porous materials by solving the inverse problem using experimental data. It is generally easy to solve the inverse problem via transmitted waves, obtaining optimized values of tortuosity, viscous and thermal characteristic lengths, but this is not the case for the porosity because of its weak sensitivity in the transmitted mode. The reflection mode is an alternative to the transmission mode, in that it gives a good estimation of porosity and tortuosity by processing the data relative to measurements of the wave reflected by the first interface. The viscous and thermal characteristic lengths cannot be determined via the first interface reflection. The wave reflected by the second interface can be experimentally detected only for the weakly resistive porous materials. In this case, the characteristic lengths can be estimated. But for common air-saturated porous materials, the second reflection is very damped and its experimental detection is difficult. We solve the inverse problem [31] numerically by the least-squares method, using both reflected and transmitted experimental data. We determine

−0.04 −0.03 −0.02 −0.01 0 0.01 0.02 0.03 0.04

air-saturated industrial plastic foams is given to validate this method.

**2.12. Reflected and transmitted waves**

Amplitude (V)

signal (dashed line).

x 10 −4

Transient Acoustic Wave Propagation in Porous Media 145

**Figure 9.** Experimental set-up of the ultrasonic measurements in transmitted mode.

0.7 ± 0.01*cm*. The measured signal generated by the transducer (measured without the sample) is given in Fig. 10 as a dashed line. The measured transmitted signal (obtained with the sample inserted) is given in the same figure (Fig. 11) as a solid line. By solving the inverse problem using the experimental transmitted data, we find the following optimized values: *α*<sup>∞</sup> = 1.26 and Λ = 66*μm*. In Figs. 11, we compare the experimental transmitted signal and simulated transmitted signal for the optimized values of tortuosity and viscous characteristic length. The difference between experiment and theory is slight, which leads us to conclude that the physical parameters are well-identified.

**Figure 10.** Experimental incident signal (dashed line) and transmitted signal (solid line) using the pair of transducers Ultran NCG200-D13 (140 kHz - 250 kHz).

## **2.11. Determination of transport parameters in air-saturated porous materials via reflected ultrasonic waves**

A method for measuring transport parameters in porous materials simultaneously, using measurements of waves reflected at the first and second interface, has been proposed

**Figure 11.** Comparison between experimental transmitted signal (solid line) and simulated transmitted signal (dashed line).

[25]. This method is based on a temporal model of direct and inverse scattering problems affecting the propagation of transient ultrasonic waves in a homogeneous isotropic slab of porous material with a rigid frame. Generally, porosity and tortuosity can be evaluated simultaneously by the wave reflected at the first interface at oblique incidence [26–29] but this is not possible at normal incidence [4, 30]. Viscous and thermal characteristic lengths can be estimated only by transmitted waves. Porosity cannot be determined from transmitted waves because of its weak sensitivity in this mode [4]. The advantage of the proposed method is that all the parameters can be determined at normal incidence (the ratio between viscous and thermal lengths is fixed as in classical methods [21–24] based on transmission measurement). The main principle of this method is the experimental detection of reflected contributions from the first and the second interface of the medium. The properties of these two contributions are used to estimate the four acoustical parameters needed for ultrasonic propagation in porous material with a rigid frame, by solving the inverse problem. Studying the sensitivity of each reflected wave parameter demonstrates the importance of each contribution) first and second interface! for the inversion. Numerical and experimental validation for weak resistive air-saturated industrial plastic foams is given to validate this method.

#### **2.12. Reflected and transmitted waves**

18 Will-be-set-by-IN-TECH

a bandwith of 6 dB extending from 140 kHz to 250 kHz. Pulses of 400 V are provided by a 5058PR Panametrics pulser/receiver. The received signals are filtered above 1 MHz to avoid high frequency noise. Electronic interference is eliminated by 1000 acquisition averages. The experimental setup is shown in Fig. 9. Consider a sample of plastic foam, of thicknesses

> Triggering Sample

**Figure 9.** Experimental set-up of the ultrasonic measurements in transmitted mode.

0.7 ± 0.01*cm*. The measured signal generated by the transducer (measured without the sample) is given in Fig. 10 as a dashed line. The measured transmitted signal (obtained with the sample inserted) is given in the same figure (Fig. 11) as a solid line. By solving the inverse problem using the experimental transmitted data, we find the following optimized values: *α*<sup>∞</sup> = 1.26 and Λ = 66*μm*. In Figs. 11, we compare the experimental transmitted signal and simulated transmitted signal for the optimized values of tortuosity and viscous characteristic length. The difference between experiment and theory is slight, which leads us to conclude

6.7 6.8 6.9 7 7.1 7.2

Time (s)

**Figure 10.** Experimental incident signal (dashed line) and transmitted signal (solid line) using the pair

**2.11. Determination of transport parameters in air-saturated porous materials via**

A method for measuring transport parameters in porous materials simultaneously, using measurements of waves reflected at the first and second interface, has been proposed

x 10 −4

Transducers

Computer

Digital oscilloscope

High frequency filtering Pre-amplifier

pulse generator

that the physical parameters are well-identified.

−0.15 −0.1 −0.05 0 0.05 0.1 0.15

Amplitude (V)

of transducers Ultran NCG200-D13 (140 kHz - 250 kHz).

**reflected ultrasonic waves**

This section concerns the ultrasonic characterization of air-saturated porous materials by solving the inverse problem using experimental data. It is generally easy to solve the inverse problem via transmitted waves, obtaining optimized values of tortuosity, viscous and thermal characteristic lengths, but this is not the case for the porosity because of its weak sensitivity in the transmitted mode. The reflection mode is an alternative to the transmission mode, in that it gives a good estimation of porosity and tortuosity by processing the data relative to measurements of the wave reflected by the first interface. The viscous and thermal characteristic lengths cannot be determined via the first interface reflection. The wave reflected by the second interface can be experimentally detected only for the weakly resistive porous materials. In this case, the characteristic lengths can be estimated. But for common air-saturated porous materials, the second reflection is very damped and its experimental detection is difficult. We solve the inverse problem [31] numerically by the least-squares method, using both reflected and transmitted experimental data. We determine simultaneously all the physical parameters intervening in the propagation. The minimization between experiment and theory is made in the time domain. The inverse problem is well posed, and its solution is unique. As with the classic ultrasonic approach for characterizing porous material saturated with one gas, the characteristic lengths are estimated by assuming a given ratio between them. Tests are performed using industrial plastic foams. Experimental and numerical results, and prospects are discussed [31].

solution of the inverse problem. The propagation equation is given (along the *x*−axis) by:

*<sup>∂</sup><sup>t</sup>* <sup>+</sup> *<sup>B</sup> <sup>∂</sup>*<sup>3</sup> *<sup>p</sup>*(*x*, *<sup>t</sup>*)

*ρ <sup>f</sup> k*<sup>0</sup> ,

> � <sup>−</sup> <sup>2</sup> *π* <sup>√</sup>*<sup>π</sup>*

−

*γα*∞*ρ*<sup>2</sup> *<sup>f</sup>* <sup>Λ</sup><sup>2</sup> � *<sup>α</sup>*<sup>0</sup> *<sup>α</sup>*<sup>∞</sup> − 1 �3

> � ∞ *x c* √*b*�

⎤ ⎦

(*t*) fields are related in time domain by the transmission

*du*, *�* =

*<sup>∂</sup><sup>t</sup>* <sup>=</sup> 0, (37)

, *<sup>A</sup>* <sup>=</sup> *ηφ*

0 2 *Pr* 2*ρ f* 2

*ηφk*0*Ka*

<sup>−</sup> *<sup>x</sup>*2<sup>Δ</sup> 16*c*<sup>2</sup> *b*� *t*

)ℵ(*ζ*,*s*)*ds*,

) �

*<sup>ζ</sup>*<sup>2</sup> <sup>−</sup> *<sup>x</sup>*<sup>2</sup> *<sup>c</sup>*<sup>2</sup> *<sup>b</sup>*�

*D*(*τ*)*G*(*t* − *τ* + *L*/*c*, *L*/*c*)*dτ*, − ∗ (36)

<sup>√</sup>*s*(*<sup>s</sup>* <sup>+</sup> 1/2*b*�

*<sup>∂</sup>t*<sup>3</sup> <sup>=</sup> 0, (35)

Transient Acoustic Wave Propagation in Porous Media 147

4*ηKa*

�

*yd y* �<sup>1</sup> <sup>−</sup> *<sup>y</sup>*<sup>2</sup> ,

*F*(*ζ*,*t*)

*<sup>ζ</sup>*<sup>2</sup> <sup>−</sup> *<sup>x</sup>*<sup>2</sup> *<sup>c</sup>*<sup>2</sup> *b*� *dζ* ⎤ ⎦ ,

�

. (38)

*η* �

*φ*3*γσ*3*ρ <sup>f</sup>* ,

, *<sup>B</sup>* <sup>=</sup> <sup>1</sup>

*ρ <sup>f</sup> k*0*γπ <sup>η</sup>* .

*<sup>∂</sup>t*<sup>2</sup> <sup>−</sup> *<sup>A</sup> <sup>∂</sup>p*(*x*, *<sup>t</sup>*)

0

*α*�

The solution of this equation gives the Green function [34] of the porous material:

�

<sup>0</sup>(*γ* − 1)*k*�

*∂*<sup>2</sup> *p*(*x*, *t*) *<sup>∂</sup>x*<sup>2</sup> <sup>−</sup> <sup>1</sup>

1 *<sup>c</sup>*<sup>2</sup> <sup>=</sup> *<sup>ρ</sup> <sup>f</sup> Ka* �

*B* = −

*<sup>G</sup>*(*x*, *<sup>t</sup>*) = <sup>1</sup>

*<sup>F</sup>*(*ζ*,*t*) = � <sup>∞</sup>

<sup>ℵ</sup>(*ζ*,*s*) = � <sup>1</sup>

*<sup>T</sup>*˜(*t*) = � *<sup>t</sup>*

*<sup>D</sup>*(*t*) = <sup>2</sup>

0

*�*√*<sup>π</sup>*

function of the porous material [33, 35, 36]:

convolution operation. The operator *D*(*t*) is given by

exp � *<sup>t</sup> B*2 �

The incident *p<sup>i</sup>*

equation

where *d* = *σφγ*

*<sup>D</sup>*(*t*) = <sup>−</sup> <sup>8</sup>

*B*2 � 1 + *t B*2 �

<sup>√</sup>*<sup>π</sup> x c* √ *b*� ⎡ ⎣ √Δ 4 *b*�

0

−1 cos ⎡ ⎣ √Δ <sup>2</sup>*b*� *<sup>ζ</sup>*

*b*� = *Bc*2, *c*� = *Ac*2, Δ = 1 + 4*b*�

*c*2

*α*0(*γ* − 1)*k*�

exp(−*s t*)

(*t*) and transmitted *p<sup>t</sup>*

1 *<sup>t</sup>*3/2 � <sup>∞</sup> 0

*∂*<sup>2</sup> *p*(*x*, *t*)

*<sup>α</sup>*0*<sup>γ</sup>* <sup>−</sup> (*<sup>γ</sup>* <sup>−</sup> <sup>1</sup>)*Prk*�

0*Prρ*<sup>2</sup> *f Kaηφ* <sup>+</sup>

*k*0

1 *<sup>t</sup>*3/2 exp �

<sup>√</sup>*s*(*<sup>s</sup>* <sup>+</sup> 1/2*b*�

*c*� .

scattering operator 31. Its expression is given by the relation given in Ref. [34]:

�

√*d*

1 *<sup>t</sup>*3/2 exp �

The transmission scattering operator is *<sup>T</sup>*˜(*t*) = *<sup>D</sup>*(*t*) <sup>∗</sup> *<sup>G</sup>*(*L*,*t*), where \* denotes the time

� +

2 <sup>√</sup>*<sup>π</sup>*

Erfc � <sup>√</sup>*<sup>t</sup> B*

*<sup>u</sup>* exp �

In the very low frequency range, the propagation equation 35 is reduced to a diffusive

− *u �* <sup>−</sup> *<sup>u</sup>*<sup>2</sup> 4*t* �

*∂p*(*x*, *t*)

*Ka* the diffusion constant. The solution of this diffusive equation gives the Green

4 *B* <sup>√</sup>*π<sup>t</sup>*

<sup>−</sup> *<sup>x</sup>*<sup>2</sup> *<sup>d</sup>* 4*t*

> � 2 *t <sup>B</sup>*<sup>2</sup> <sup>+</sup> <sup>1</sup>

�

�

<sup>√</sup>*<sup>s</sup>* <sup>−</sup> *<sup>y</sup>*

� *u*<sup>2</sup> <sup>2</sup>*<sup>t</sup>* <sup>−</sup> <sup>1</sup>

> *∂*<sup>2</sup> *p*(*x*, *t*) *<sup>∂</sup>x*<sup>2</sup> <sup>−</sup> *<sup>d</sup>*

*<sup>G</sup>*(*x*, *<sup>t</sup>*) = *<sup>x</sup>*

## **2.13. Ultrasonic characterization of air-saturated double-layered porous media**

This section concerns a time-domain model of transient wave propagation in double-layered porous materials [32]. An analytical derivation of reflection and transmission scattering operators is given in the time domain. These scattering kernels are the medium's responses to an incident acoustic pulse. The expressions obtained take into account the multiple reflections occurring at the interfaces of the double-layered material. The double-layered porous media consist of two slabs of homogeneous isotropic porous materials with a rigid frame. Each porous slab is described by a temporal equivalent fluid model, in which the acoustic wave propagates only in the fluid saturating the material. In this model, the inertial effects are described by the tortuosity; the viscous and thermal losses of the medium are described by two susceptibility kernels which depend on the viscous and thermal characteristic lengths. the inverse problem is solved using experimental reflected signals at normal incidence. The sensitivity of porosity, tortuosity, and viscous characteristic length of each layer is studied showing their effect on the reflected interface waveforms. The inverse problem is solved numerically by the least-squares method. Five parameters are inverted: porosity and tortuosity of the two layers and the viscous characteristic length of the first layer. The minimization of the discrepancy between experimental and theoretical data is made in the time domain. The inverse problem is shown to be well posed and its solution to be unique. Experimental results for waves reflected by the interfaces of the double-layered porous material are given and compared with theoretical predictions [32].

### **2.14. Acoustic characterization of porous materials at low frequency range**

In this frequency range, the viscous forces are important throughout the fluid, the cycle of compression/expansion in the porous medium is slow enough to allow heat exchange between fluid and structure. At the same time, the temperature of the structure is practically unchanged by the passage of the sound wave due to the large value of its specific heat conducting relation to that of the fluid: the structure acts as a thermostat. In this case the isothermal compression is directly applicable. This range corresponds to frequencies where the viscous skin depth *δ* = (2*η*/*ωρ <sup>f</sup>*)1/2 is larger than the pore radius.

An acoustic method [33–37] using transmitted and reflected waves is proposed for measuring static viscous permeability *k*0, flow resistivity *σ*, static thermal permeability *k* <sup>0</sup> and the inertial factor *α*<sup>0</sup> (low frequency tortuosity), of porous materials having a rigid frame at low frequencies. Flow resistivity of porous material is defined as the ratio between the pressure difference across a sample and the velocity of flow of air through that sample per unit cube. The static thermal permeability of porous material is a geometrical parameter equal to the inverse trapping constant of the solid frame. The inertial factor describes the fluid structure interactions in the low frequency range (1-3 kHz). These parameters are determined from the solution of the inverse problem. The propagation equation is given (along the *x*−axis) by:

20 Will-be-set-by-IN-TECH

simultaneously all the physical parameters intervening in the propagation. The minimization between experiment and theory is made in the time domain. The inverse problem is well posed, and its solution is unique. As with the classic ultrasonic approach for characterizing porous material saturated with one gas, the characteristic lengths are estimated by assuming a given ratio between them. Tests are performed using industrial plastic foams. Experimental

**2.13. Ultrasonic characterization of air-saturated double-layered porous media**

porous material are given and compared with theoretical predictions [32].

the viscous skin depth *δ* = (2*η*/*ωρ <sup>f</sup>*)1/2 is larger than the pore radius.

**2.14. Acoustic characterization of porous materials at low frequency range**

In this frequency range, the viscous forces are important throughout the fluid, the cycle of compression/expansion in the porous medium is slow enough to allow heat exchange between fluid and structure. At the same time, the temperature of the structure is practically unchanged by the passage of the sound wave due to the large value of its specific heat conducting relation to that of the fluid: the structure acts as a thermostat. In this case the isothermal compression is directly applicable. This range corresponds to frequencies where

An acoustic method [33–37] using transmitted and reflected waves is proposed for measuring

inertial factor *α*<sup>0</sup> (low frequency tortuosity), of porous materials having a rigid frame at low frequencies. Flow resistivity of porous material is defined as the ratio between the pressure difference across a sample and the velocity of flow of air through that sample per unit cube. The static thermal permeability of porous material is a geometrical parameter equal to the inverse trapping constant of the solid frame. The inertial factor describes the fluid structure interactions in the low frequency range (1-3 kHz). These parameters are determined from the

<sup>0</sup> and the

static viscous permeability *k*0, flow resistivity *σ*, static thermal permeability *k*

This section concerns a time-domain model of transient wave propagation in double-layered porous materials [32]. An analytical derivation of reflection and transmission scattering operators is given in the time domain. These scattering kernels are the medium's responses to an incident acoustic pulse. The expressions obtained take into account the multiple reflections occurring at the interfaces of the double-layered material. The double-layered porous media consist of two slabs of homogeneous isotropic porous materials with a rigid frame. Each porous slab is described by a temporal equivalent fluid model, in which the acoustic wave propagates only in the fluid saturating the material. In this model, the inertial effects are described by the tortuosity; the viscous and thermal losses of the medium are described by two susceptibility kernels which depend on the viscous and thermal characteristic lengths. the inverse problem is solved using experimental reflected signals at normal incidence. The sensitivity of porosity, tortuosity, and viscous characteristic length of each layer is studied showing their effect on the reflected interface waveforms. The inverse problem is solved numerically by the least-squares method. Five parameters are inverted: porosity and tortuosity of the two layers and the viscous characteristic length of the first layer. The minimization of the discrepancy between experimental and theoretical data is made in the time domain. The inverse problem is shown to be well posed and its solution to be unique. Experimental results for waves reflected by the interfaces of the double-layered

and numerical results, and prospects are discussed [31].

$$\begin{split} \frac{\partial^2 p(\mathbf{x}, t)}{\partial \mathbf{x}^2} - \frac{1}{c^2} \frac{\partial^2 p(\mathbf{x}, t)}{\partial t^2} - A \frac{\partial p(\mathbf{x}, t)}{\partial t} + B \frac{\partial^3 p(\mathbf{x}, t)}{\partial t^3} &= 0, \\ \frac{1}{c^2} = \frac{\rho\_f}{K\_d} \left( a\_0 \gamma - \frac{(\gamma - 1) P\_r k\_0'}{k\_0} \right), \quad A = \frac{\eta \phi}{\rho\_f k\_0'}, \\ B = -\frac{a\_0 (\gamma - 1) k\_0' P\_r \rho\_f^2}{K\_d \eta \phi} + \frac{a\_0' (\gamma - 1) k\_0'^2 P\_r^2 \rho\_f^2}{\eta \phi k\_0 K\_d} - \frac{\gamma a\_0 \rho\_f^2 \Lambda^2 \left( \frac{a\_0}{a\_\infty} - 1 \right)^3}{4 \eta K\_d} \end{split}$$

The solution of this equation gives the Green function [34] of the porous material:

$$\begin{split} G(x,t) &= \frac{1}{\sqrt{\pi}} \frac{x}{c} \sqrt{b'} \left[ \frac{\sqrt{\Delta}}{4b'} \frac{1}{b^{3/2}} \exp\left(-\frac{x^2 \Delta}{16c^2 b' t}\right) - \frac{2}{\pi \sqrt{\pi}} \int\_{\frac{x}{c}\sqrt{b'}}^{\infty} \frac{F(\zeta, t)}{\sqrt{\zeta^2 - \frac{x^2}{c^2} b'}} d\zeta \right], \\ F(\zeta, t) &= \int\_0^{\infty} \exp(-st) \sqrt{s} (s + 1/2b') \aleph(\zeta, s) ds, \\ \aleph(\zeta, s) &= \int\_{-1}^1 \exp\left[\frac{\sqrt{\Delta}}{2b'} \zeta \sqrt{s} - y\sqrt{s} (s + 1/2b') \sqrt{\zeta^2 - \frac{x^2}{c^2} b'}\right] \frac{y d\, y}{\sqrt{1 - y^2}}, \\ b' &= Bc^2, c' = Ac^2, \Delta = 1 + 4b'c'. \end{split}$$

The incident *p<sup>i</sup>* (*t*) and transmitted *p<sup>t</sup>* (*t*) fields are related in time domain by the transmission scattering operator 31. Its expression is given by the relation given in Ref. [34]:

$$T(t) = \int\_0^t D(\tau)G(t-\tau+L/c, L/c)d\tau, -\ast \tag{36}$$

$$D(t) = \frac{2}{\mathcal{O}\sqrt{\pi}}\frac{1}{t^{3/2}}\int\_0^\infty \left(\frac{u^2}{2t}-1\right)u\exp\left(-\frac{u}{\mathcal{O}}-\frac{u^2}{4t}\right)du, \qquad \mathcal{O} = \sqrt{\frac{\rho\_f k\_0 \gamma \pi}{\eta}}.$$

In the very low frequency range, the propagation equation 35 is reduced to a diffusive equation

$$\frac{\partial^2 p(\mathbf{x}, t)}{\partial \mathbf{x}^2} - d \frac{\partial p(\mathbf{x}, t)}{\partial t} = 0,\tag{37}$$

where *d* = *σφγ Ka* the diffusion constant. The solution of this diffusive equation gives the Green function of the porous material [33, 35, 36]:

$$G(\mathbf{x}, t) = \frac{\mathbf{x}\sqrt{d}}{2\sqrt{\pi}} \frac{1}{t^{3/2}} \exp\left(-\frac{\mathbf{x}^2 d}{4t}\right). \tag{38}$$

The transmission scattering operator is *<sup>T</sup>*˜(*t*) = *<sup>D</sup>*(*t*) <sup>∗</sup> *<sup>G</sup>*(*L*,*t*), where \* denotes the time convolution operation. The operator *D*(*t*) is given by

$$D(t) = -\frac{8}{B^2} \left( 1 + \frac{t}{B^2} \right) \exp\left(\frac{t}{B^2}\right) \text{Erfc}\left(\frac{\sqrt{t}}{B}\right) + \frac{4}{B\sqrt{\pi t}} \left( 2\frac{t}{B^2} + 1 \right), \qquad B = \frac{1}{\eta} \sqrt{\phi^3 \eta \sigma^3 \rho\_{f'}}$$

where Erfc is the complementary error function and *ρ <sup>f</sup>* is the fluid density.

The inverse problem is to find the parameters *k*� and *α*<sup>0</sup> [37] which minimize numerically the discrepancy function *U*(*k*� , *α*0) = ∑*i*=*<sup>N</sup> <sup>i</sup>*=<sup>1</sup> (*p<sup>t</sup> exp*(*x*, *ti*) <sup>−</sup> *<sup>p</sup><sup>t</sup>* (*x*, *ti*))2, wherein *p<sup>t</sup> exp*(*x*, *ti*)*i*=1,2,...*<sup>n</sup>* is the discrete set of values of the experimental transmitted signal and *p<sup>t</sup>* (*x*, *ti*)*i*=1,2,...*<sup>n</sup>* the discrete set of values of the simulated transmitted signal predicted from Eq. (31). The inverse problem is solved numerically by the least-square method. For its iterative solution, we used the simplex search method (Nedler Mead) [38] which does not require numerical or analytic gradients. Experiments are performed in a guide [33–37] (pipe), having a diameter of 5 cm. The experimental set up is given in Fig. 12. For measuring the static thermal permeability and the inertial factor, a pipe of 3 m long, since the frequencies used in the experiment are between 1kHz and 4 kHz. However, for measuring the viscous permeability (and flow resistivity) [33, 35, 36], a pipe of 50 m long must be used (50Hz-1kHz) to avoid the reflections at its end. In this case, it is not important to keep the pipe straight, it can be rolled in order to save space without perturbations on experimental signals (the cut-off frequency of the tube *fc* ∼ 4kHz). A sound source Driver unit "Brand" constituted by loudspeaker Realistic 40-9000 is used. Bursts are provided by synthesized function generator Standford Research Systems model DS345-30MHz. The signals are amplified and filtered using model SR 650-Dual channel filter, Standford Research Systems. The signals (incident and transmitted) are measured using the same microphone (Bruel&Kjaer, 4190) in the same position in the tube. The incident signal is measured without porous sample, however, the transmitted signal is measured with the porous sample. Consider a cylindrical sample of plastic foam M of diameter 5 cm, porosity *φ* = 0.85 and thickness 4.15 cm. The viscous permeability *k*<sup>0</sup> of the porous sample is measured by solving the inverse problem in time domain at viscous domain (very low frequency range)[33, 35, 36]. The obtained value of *<sup>k</sup>*<sup>0</sup> is (0.44 <sup>±</sup> 0.02) <sup>×</sup> <sup>10</sup>−<sup>9</sup> *<sup>m</sup>*2. The static thermal permeability is measured using classic Kundt tube (continuous frequency) [39] (continuous frequency) obtaining the value of *k*� <sup>0</sup> = (1.54 <sup>±</sup> 0.05) <sup>×</sup> <sup>10</sup>−9, with a ratio of 3.5 between *<sup>k</sup>*� 0 and *k*0. The tortuosity *α*<sup>∞</sup> and the viscous characteristic length Λ are measured in the high frequency range [4] obtaining the values: *α*<sup>∞</sup> = 1.2 and Λ = 65*μm*. The value of the thermal tortuosity *α*� <sup>0</sup> is approximated by the relation *α*� <sup>0</sup> � *α*0/*α*∞. Fig. 13 shows the experimental incident signal (solid line) generated by the loudspeaker in the frequency bandwidth (2.5 - 3.5) kHz, and the experimental transmitted signal (dashed line). After solving the inverse problem numerically for the thermal permeability *k*� <sup>0</sup> and the inertial factor *α*<sup>0</sup> , we find the following optimized values: *k*� <sup>0</sup> = (1.06 <sup>±</sup> 0.2) <sup>×</sup> <sup>10</sup>−9*m*<sup>2</sup> and *<sup>α</sup>*<sup>0</sup> <sup>=</sup> 2.15 <sup>±</sup> 0.15. We present in Fig. 14 the variation of the minimization function *U* with the thermal permeability *k*� <sup>0</sup> and tortuosity *α*0. The obtained ratio of 2.4 between the inverted *k*� <sup>0</sup> and *k*<sup>0</sup> is in adequacy with the ratio given in literature [9, 39], which is found generally to be between 2 and 4 for plastic foams. In Fig. 15, we show a comparison between an experimental transmitted signal and simulated transmitted signal for the optimized values of thermal permeability and tortuosity. The difference between the two curves is slight, which leads us to conclude that the optimized values of the thermal permeability and inertial factor are correct. This study has been carried out, in the frequency bandwidths (1.5 - 2.5) kHz and has also given good results.

**Figure 12.** Experimental setup of acoustic measurements.

Amplitude (a.u.)

−0.3 −0.2 −0.1 0 0.1 0.2 0.3

0.3 0.4 0.5

0.25

U

**Figure 14.** Minima of the inverse problem.

0.3

U

0.088 0.089 0.09 0.091

Transient Acoustic Wave Propagation in Porous Media 149

Time (s)

1.5 2 2.5

α0

k'0 /k0

skeleton and the saturating fluid. The fluid-structure interactions then play an essential role in the propagation. The study of theses effects has been largely developed by Biot [40, 41] for applications in the domain of oil exploration since 1950. This semi-phenomenological theory provides a rigorous description of the propagation of acoustical waves in porous

1.5 2 2.5 3

**Figure 13.** Experimental incident signal (solid line) and experimental transmitted signal (dashed line).

### **3. Modified Biot theory**

When the structure of a porous material is not rigid, the equivalent fluid model, on which the previous sections are based, no longer applies because the waves propagate in both the

**Figure 12.** Experimental setup of acoustic measurements.

22 Will-be-set-by-IN-TECH

The inverse problem is to find the parameters *k*� and *α*<sup>0</sup> [37] which minimize numerically

discrete set of values of the simulated transmitted signal predicted from Eq. (31). The inverse problem is solved numerically by the least-square method. For its iterative solution, we used the simplex search method (Nedler Mead) [38] which does not require numerical or analytic gradients. Experiments are performed in a guide [33–37] (pipe), having a diameter of 5 cm. The experimental set up is given in Fig. 12. For measuring the static thermal permeability and the inertial factor, a pipe of 3 m long, since the frequencies used in the experiment are between 1kHz and 4 kHz. However, for measuring the viscous permeability (and flow resistivity) [33, 35, 36], a pipe of 50 m long must be used (50Hz-1kHz) to avoid the reflections at its end. In this case, it is not important to keep the pipe straight, it can be rolled in order to save space without perturbations on experimental signals (the cut-off frequency of the tube *fc* ∼ 4kHz). A sound source Driver unit "Brand" constituted by loudspeaker Realistic 40-9000 is used. Bursts are provided by synthesized function generator Standford Research Systems model DS345-30MHz. The signals are amplified and filtered using model SR 650-Dual channel filter, Standford Research Systems. The signals (incident and transmitted) are measured using the same microphone (Bruel&Kjaer, 4190) in the same position in the tube. The incident signal is measured without porous sample, however, the transmitted signal is measured with the porous sample. Consider a cylindrical sample of plastic foam M of diameter 5 cm, porosity *φ* = 0.85 and thickness 4.15 cm. The viscous permeability *k*<sup>0</sup> of the porous sample is measured by solving the inverse problem in time domain at viscous domain (very low frequency range)[33, 35, 36]. The obtained value of *<sup>k</sup>*<sup>0</sup> is (0.44 <sup>±</sup> 0.02) <sup>×</sup> <sup>10</sup>−<sup>9</sup> *<sup>m</sup>*2. The static thermal permeability is measured using classic Kundt tube (continuous frequency) [39] (continuous

and *k*0. The tortuosity *α*<sup>∞</sup> and the viscous characteristic length Λ are measured in the high frequency range [4] obtaining the values: *α*<sup>∞</sup> = 1.2 and Λ = 65*μm*. The value of the thermal

incident signal (solid line) generated by the loudspeaker in the frequency bandwidth (2.5 - 3.5) kHz, and the experimental transmitted signal (dashed line). After solving the inverse

the ratio given in literature [9, 39], which is found generally to be between 2 and 4 for plastic foams. In Fig. 15, we show a comparison between an experimental transmitted signal and simulated transmitted signal for the optimized values of thermal permeability and tortuosity. The difference between the two curves is slight, which leads us to conclude that the optimized values of the thermal permeability and inertial factor are correct. This study has been carried

When the structure of a porous material is not rigid, the equivalent fluid model, on which the previous sections are based, no longer applies because the waves propagate in both the

in Fig. 14 the variation of the minimization function *U* with the thermal permeability *k*�

out, in the frequency bandwidths (1.5 - 2.5) kHz and has also given good results.

*exp*(*x*, *ti*) <sup>−</sup> *<sup>p</sup><sup>t</sup>*

(*x*, *ti*))2, wherein *p<sup>t</sup>*

<sup>0</sup> = (1.54 <sup>±</sup> 0.05) <sup>×</sup> <sup>10</sup>−9, with a ratio of 3.5 between *<sup>k</sup>*�

<sup>0</sup> = (1.06 <sup>±</sup> 0.2) <sup>×</sup> <sup>10</sup>−9*m*<sup>2</sup> and *<sup>α</sup>*<sup>0</sup> <sup>=</sup> 2.15 <sup>±</sup> 0.15. We present

<sup>0</sup> � *α*0/*α*∞. Fig. 13 shows the experimental

<sup>0</sup> and the inertial factor *α*<sup>0</sup> , we find the

<sup>0</sup> and *k*<sup>0</sup> is in adequacy with

*exp*(*x*, *ti*)*i*=1,2,...*<sup>n</sup>*

(*x*, *ti*)*i*=1,2,...*<sup>n</sup>* the

0

<sup>0</sup> and

*<sup>i</sup>*=<sup>1</sup> (*p<sup>t</sup>*

where Erfc is the complementary error function and *ρ <sup>f</sup>* is the fluid density.

, *α*0) = ∑*i*=*<sup>N</sup>*

is the discrete set of values of the experimental transmitted signal and *p<sup>t</sup>*

the discrepancy function *U*(*k*�

frequency) obtaining the value of *k*�

following optimized values: *k*�

**3. Modified Biot theory**

<sup>0</sup> is approximated by the relation *α*�

tortuosity *α*0. The obtained ratio of 2.4 between the inverted *k*�

problem numerically for the thermal permeability *k*�

tortuosity *α*�

**Figure 13.** Experimental incident signal (solid line) and experimental transmitted signal (dashed line).

**Figure 14.** Minima of the inverse problem.

skeleton and the saturating fluid. The fluid-structure interactions then play an essential role in the propagation. The study of theses effects has been largely developed by Biot [40, 41] for applications in the domain of oil exploration since 1950. This semi-phenomenological theory provides a rigorous description of the propagation of acoustical waves in porous

wherein *P*, *Q* and *R* are generalized elastic constants which are related, via gedanken experiments, to other, measurable quantities, namely *φ* (porosity), *Kf* (bulk modulus of the pore fluid), *Ks* (bulk modulus of the elastic solid) and *Kb* (bulk modulus of the porous skeletal frame). *N* is the shear modulus of the composite as well as that of the skeletal frame. The

<sup>3</sup> *<sup>N</sup>*, *<sup>Q</sup>*<sup>=</sup> (1−*φ*<sup>−</sup> *Kb*

3(1 − 2*νb*)

The Young modulus and the Poisson ratio of the solid *Es*, *ν<sup>s</sup>* and of the skeletal frame *Eb*, *ν<sup>b</sup>*

The temporal operators *<sup>ρ</sup>*<sup>11</sup>(*t*), *<sup>ρ</sup>*<sup>12</sup>(*t*) et *<sup>ρ</sup>*<sup>22</sup>(*t*) are the operators of mass coupling between

Biot coefficients *ρmn*, are the "mass factors" and are connected to the densities of solid (*ρs*) and

As in the case of an elastic solid, the wave equations of dilatational and rotational waves can be obtained using scalar and vector displacement potentials, respectively. Two scalar potentials

−→*<sup>u</sup>* <sup>=</sup> −→∇<sup>Φ</sup>*s*, −→*<sup>U</sup>* <sup>=</sup> −→∇<sup>Φ</sup>*<sup>f</sup>* .

<sup>1</sup> <sup>−</sup><sup>1</sup>

The result of solving this system of equations is the existence of two distinct longitudinal modes called slow mode and fast mode. The system (41) can be expressed on the basis of slow

> *λ*1(*t*) 0 0 *λ*2(*t*)

<sup>−</sup>1 1 *<sup>∂</sup>*3/2

*<sup>∂</sup>t*3/2 <sup>Φ</sup>*<sup>s</sup>* Φ*f* =

> <sup>Φ</sup><sup>1</sup> Φ<sup>2</sup> ,

<sup>1</sup>−*φ*<sup>−</sup> *Kb Ks* +*φ Ks Kf*

*Ks* )*φKs*

, *<sup>ρ</sup>*<sup>12</sup>(*t*) = *<sup>ρ</sup>*<sup>12</sup> <sup>−</sup> <sup>2</sup>*φα*<sup>∞</sup>

, *<sup>N</sup>* <sup>=</sup> *Eb*

, *<sup>R</sup>*<sup>=</sup> *<sup>φ</sup>*2*Ks* <sup>1</sup>−*φ*<sup>−</sup> *Kb*

Transient Acoustic Wave Propagation in Porous Media 151

2(1 + *νb*)

 *ρ <sup>f</sup> η πt*

 *P Q Q R <sup>∂</sup>*<sup>2</sup> *∂x*<sup>2</sup>  Φ*<sup>s</sup>* Φ*f* ,

(41)

1/2 ,

Λ

*Ks* <sup>+</sup> *<sup>φ</sup> Ks Kf* .

. (40)

equations which explicitly relate *P*, *Q*, and *R* to *φ*, *Kf* , *Ks*, *Kb*, and *N* are given by

, *Kb* <sup>=</sup> *Eb*

1/2

1/2 .

+ 4

depend on the generalized elastic constant *P*, *Q* and *R* via the relations:

fluid and structure and are given in the high frequency range by [47]:

 *ρ <sup>f</sup> η πt*

 *ρ <sup>f</sup> η πt*

fluid (*ρ <sup>f</sup>*) by *ρ*<sup>11</sup> + *ρ*<sup>12</sup> = (1 − *φ*)*ρ<sup>s</sup>* and *ρ*<sup>12</sup> + *ρ*<sup>22</sup> = *φρ <sup>f</sup>* , *ρ*<sup>12</sup> = −*φρ <sup>f</sup>*(*α*<sup>∞</sup> − 1).

In this case the equations 39 become for a propagation along the *x* axis:

 *η ρ f*

2*φρ <sup>f</sup> α*<sup>∞</sup> Λ

*∂*2 *∂x*<sup>2</sup>  <sup>Φ</sup><sup>1</sup> Φ<sup>2</sup> = 

2*φα*∞ Λ

2*φα*∞ Λ

*P*=

(1−*φ*)(1−*φ*<sup>−</sup> *Kb*

<sup>1</sup>−*φ*<sup>−</sup> *Kb Ks* +*φ Ks Kf*

*Ks*

*Ks* <sup>=</sup> *Es*

*<sup>ρ</sup>*<sup>11</sup>(*t*) = *<sup>ρ</sup>*<sup>11</sup> <sup>+</sup>

*<sup>ρ</sup>*<sup>22</sup>(*t*) = *<sup>ρ</sup>*<sup>22</sup> <sup>+</sup>

for the frame and the fluid Φ*<sup>s</sup>* and Φ*<sup>f</sup>*

 Φ*<sup>s</sup>* Φ*f* +

and fast waves Φ<sup>1</sup> and Φ<sup>2</sup> respectively by:

*∂t*<sup>2</sup>

*3.1.1. Longitudinal waves*

 *<sup>ρ</sup>*<sup>11</sup> *<sup>ρ</sup>*<sup>12</sup> *<sup>ρ</sup>*<sup>12</sup> *<sup>ρ</sup>*<sup>22</sup> *<sup>∂</sup>*<sup>2</sup> )*Ks* + *φ Ks Kf Kb*

3(1 − 2*νs*)

**Figure 15.** Comparison between the experimental transmitted signal (solid line) and the simulated transmitted signals (dashed line) using the reconstructed values of *k*� <sup>0</sup> and *α*0.

media saturated by a compressible viscous fluid. Such diphasic materials are supposed to be elastic and homogeneous. Biot's theory was initially introduced for petroleum prospecting and research. Due to its very general and rather fundamental character, it has been applied in various fields of acoustics such as geophysics, underwater acoustics, seismology, ultrasonic characterization of bones, etc. Cancelous bone is a porous material consisting of a matrix of solid trabeculae filled with soft bone marrow. The interaction between ultrasound and bone is highly complex. Modeling ultrasonic propagation through trabecular tissue has been considered using porous media theories, such as Biot's theory. This theory derives the equations of motion for each phase (i.e. the solid frame and the fluid) based on energy considerations which include the inertial, potential and viscous coupling between the two phases. For an isotropic porous medium, three different bulk modes are predicted, i.e. two compression waves and one shear wave. One compressional wave, the so-called wave of the first type or fast longitudinal wave, and the transverse wave are similar to the two bulk waves observed in an anisotropic linear elastic solid. The other longitudinal wave, called a wave of the second kind, or slow wave, is a highly damped and very dispersive mode. It is diffusive at low frequencies and propagative at high frequencies. In this section, the direct and inverse scattering problems for the ultrasonic propagation in porous material having an elastic frame are solved [42–46]. An experimental application on human cancelous bone is discussed [42, 43].

#### **3.1. Model**

The equations of motion of the frame and fluid are given by the Euler equations applied to the Lagrangian density. Here −→*u* and −→*U* are the displacements of the solid and fluid phases. The equations of motion are:

$$
\begin{split}
\widetilde{\rho}\_{11}(t) \* \frac{\partial^2 \overrightarrow{\boldsymbol{u}}^{\flat}}{\partial t^2} + \widetilde{\rho}\_{12}(t) \* \frac{\partial^2 \overrightarrow{\boldsymbol{U}}}{\partial t^2} &= P \overrightarrow{\nabla} (\overrightarrow{\nabla} \cdot \overrightarrow{\boldsymbol{u}}^{\flat}) + Q \overrightarrow{\nabla} (\overrightarrow{\nabla} \cdot \overrightarrow{\boldsymbol{U}}) - N \overrightarrow{\nabla} \wedge (\overrightarrow{\nabla} \wedge \overrightarrow{\boldsymbol{u}}^{\flat}), \\
\widetilde{\rho}\_{12}(t) \* \frac{\partial^2 \overrightarrow{\boldsymbol{u}}^{\flat}}{\partial t^2} + \widetilde{\rho}\_{22}(t) \* \frac{\partial^2 \overrightarrow{\boldsymbol{U}}}{\partial t^2} &= Q \overrightarrow{\nabla} (\overrightarrow{\nabla} \cdot \overrightarrow{\boldsymbol{u}}^{\flat}) + R \overrightarrow{\nabla} (\overrightarrow{\nabla} \cdot \overrightarrow{\boldsymbol{U}}).
\end{split}
\tag{39}
$$

wherein *P*, *Q* and *R* are generalized elastic constants which are related, via gedanken experiments, to other, measurable quantities, namely *φ* (porosity), *Kf* (bulk modulus of the pore fluid), *Ks* (bulk modulus of the elastic solid) and *Kb* (bulk modulus of the porous skeletal frame). *N* is the shear modulus of the composite as well as that of the skeletal frame. The equations which explicitly relate *P*, *Q*, and *R* to *φ*, *Kf* , *Ks*, *Kb*, and *N* are given by

$$P = \frac{(1 - \phi)(1 - \phi - \frac{K\_b}{K\_s})K\_s + \phi \frac{K\_s}{K\_f}K\_b}{1 - \phi - \frac{K\_b}{K\_s} + \phi \frac{K\_s}{K\_f}} + \frac{4}{3}N\_s \quad Q = \frac{(1 - \phi - \frac{K\_b}{K\_s})\phi K\_s}{1 - \phi - \frac{K\_b}{K\_s} + \phi \frac{K\_s}{K\_f}} \prime \quad R = \frac{\phi^2 K\_s}{1 - \phi - \frac{K\_b}{K\_s} + \phi \frac{K\_s}{K\_f}} \prime$$

The Young modulus and the Poisson ratio of the solid *Es*, *ν<sup>s</sup>* and of the skeletal frame *Eb*, *ν<sup>b</sup>* depend on the generalized elastic constant *P*, *Q* and *R* via the relations:

$$K\_s = \frac{E\_s}{3(1 - 2\nu\_s)}, \qquad K\_b = \frac{E\_b}{3(1 - 2\nu\_b)}, \qquad N = \frac{E\_b}{2(1 + \nu\_b)}.\tag{40}$$

The temporal operators *<sup>ρ</sup>*<sup>11</sup>(*t*), *<sup>ρ</sup>*<sup>12</sup>(*t*) et *<sup>ρ</sup>*<sup>22</sup>(*t*) are the operators of mass coupling between fluid and structure and are given in the high frequency range by [47]:

$$\begin{split} \widetilde{\rho}\_{11}(t) &= \rho\_{11} + \frac{2\phi\alpha\_{\infty}}{\Lambda} \left( \frac{\rho\_f \eta}{\pi t} \right)^{1/2}, \quad \widetilde{\rho}\_{12}(t) = \rho\_{12} - \frac{2\phi\alpha\_{\infty}}{\Lambda} \left( \frac{\rho\_f \eta}{\pi t} \right)^{1/2}, \\\widetilde{\rho}\_{22}(t) &= \rho\_{22} + \frac{2\phi\alpha\_{\infty}}{\Lambda} \left( \frac{\rho\_f \eta}{\pi t} \right)^{1/2}. \end{split}$$

Biot coefficients *ρmn*, are the "mass factors" and are connected to the densities of solid (*ρs*) and fluid (*ρ <sup>f</sup>*) by *ρ*<sup>11</sup> + *ρ*<sup>12</sup> = (1 − *φ*)*ρ<sup>s</sup>* and *ρ*<sup>12</sup> + *ρ*<sup>22</sup> = *φρ <sup>f</sup>* , *ρ*<sup>12</sup> = −*φρ <sup>f</sup>*(*α*<sup>∞</sup> − 1).

#### *3.1.1. Longitudinal waves*

24 Will-be-set-by-IN-TECH

0.088 0.089 0.09 0.091

<sup>0</sup> and *α*0.

Time (s)

media saturated by a compressible viscous fluid. Such diphasic materials are supposed to be elastic and homogeneous. Biot's theory was initially introduced for petroleum prospecting and research. Due to its very general and rather fundamental character, it has been applied in various fields of acoustics such as geophysics, underwater acoustics, seismology, ultrasonic characterization of bones, etc. Cancelous bone is a porous material consisting of a matrix of solid trabeculae filled with soft bone marrow. The interaction between ultrasound and bone is highly complex. Modeling ultrasonic propagation through trabecular tissue has been considered using porous media theories, such as Biot's theory. This theory derives the equations of motion for each phase (i.e. the solid frame and the fluid) based on energy considerations which include the inertial, potential and viscous coupling between the two phases. For an isotropic porous medium, three different bulk modes are predicted, i.e. two compression waves and one shear wave. One compressional wave, the so-called wave of the first type or fast longitudinal wave, and the transverse wave are similar to the two bulk waves observed in an anisotropic linear elastic solid. The other longitudinal wave, called a wave of the second kind, or slow wave, is a highly damped and very dispersive mode. It is diffusive at low frequencies and propagative at high frequencies. In this section, the direct and inverse scattering problems for the ultrasonic propagation in porous material having an elastic frame are solved [42–46]. An experimental application on human cancelous bone is

The equations of motion of the frame and fluid are given by the Euler equations applied to the Lagrangian density. Here −→*u* and −→*U* are the displacements of the solid and fluid phases.

−→∇.

−→*<sup>u</sup>* ) + *<sup>Q</sup>*−→∇(

−→∇( −→∇.

−→*u* ) + *R*

−→∇.

−→*<sup>U</sup>* ) <sup>−</sup> *<sup>N</sup>*−→∇ ∧ (

−→∇ ∧ −→*<sup>u</sup>* ),

−→*<sup>U</sup>* ). (39)

−→∇.( −→∇.

**Figure 15.** Comparison between the experimental transmitted signal (solid line) and the simulated

−0.04 −0.03 −0.02 −0.01 0 0.01 0.02 0.03 0.04 0.05

transmitted signals (dashed line) using the reconstructed values of *k*�

discussed [42, 43].

The equations of motion are:

*<sup>ρ</sup>*<sup>11</sup>(*t*) <sup>∗</sup> *<sup>∂</sup>*2−→*<sup>u</sup>*

*<sup>ρ</sup>*<sup>12</sup>(*t*) <sup>∗</sup> *<sup>∂</sup>*2−→*<sup>u</sup>*

*<sup>∂</sup>t*<sup>2</sup> <sup>+</sup> *<sup>ρ</sup>*<sup>12</sup>(*t*) <sup>∗</sup> *<sup>∂</sup>*2−→*<sup>U</sup>*

*<sup>∂</sup>t*<sup>2</sup> <sup>+</sup> *<sup>ρ</sup>*<sup>22</sup>(*t*) <sup>∗</sup> *<sup>∂</sup>*2−→*<sup>U</sup>*

*<sup>∂</sup>t*<sup>2</sup> <sup>=</sup> *<sup>P</sup>*

*<sup>∂</sup>t*<sup>2</sup> <sup>=</sup> *<sup>Q</sup>*−→∇(

**3.1. Model**

Amplitude (a.u.)

As in the case of an elastic solid, the wave equations of dilatational and rotational waves can be obtained using scalar and vector displacement potentials, respectively. Two scalar potentials for the frame and the fluid Φ*<sup>s</sup>* and Φ*<sup>f</sup>*

$$
\overrightarrow{\boldsymbol{\mu}} = \overrightarrow{\nabla} \boldsymbol{\Phi}\_{\sf s\nu} \qquad \overrightarrow{\boldsymbol{\mathcal{U}}} = \overrightarrow{\nabla} \boldsymbol{\Phi}\_{\sf f}.
$$

In this case the equations 39 become for a propagation along the *x* axis:

$$
\begin{pmatrix}
\rho\_{11}\,\rho\_{12} \\
\rho\_{12}\,\rho\_{22}
\end{pmatrix}
\frac{
\partial^2}{\partial t^2}
\begin{pmatrix}
\Phi\_s \\
\Phi\_f
\end{pmatrix}
+
\frac{2\phi\rho\_f\alpha\_{\infty}}{\Lambda}
\sqrt{\frac{\eta}{\rho\_f}}
\begin{pmatrix}
1 & -1 \\
\end{pmatrix}
\frac{
\partial^{3/2}}{\partial t^{3/2}}
\begin{pmatrix}
\Phi\_s \\
\Phi\_f
\end{pmatrix}
=
\begin{pmatrix}
P \ Q \\
Q \ R
\end{pmatrix}
\frac{
\partial^2}{\partial \mathbf{x}^2}
\begin{pmatrix}
\Phi\_s \\
\Phi\_f
\end{pmatrix}
\tag{41}
$$

The result of solving this system of equations is the existence of two distinct longitudinal modes called slow mode and fast mode. The system (41) can be expressed on the basis of slow and fast waves Φ<sup>1</sup> and Φ<sup>2</sup> respectively by:

$$
\frac{\partial^2}{\partial x^2} \begin{pmatrix} \Phi\_1\\ \Phi\_2 \end{pmatrix} = \begin{pmatrix} \widetilde{\lambda}\_1(t) & 0\\ 0 & \widetilde{\lambda}\_2(t) \end{pmatrix} \begin{pmatrix} \Phi\_1\\ \Phi\_2 \end{pmatrix} .
$$

#### 26 Will-be-set-by-IN-TECH 152 Modeling and Measurement Methods for Acoustic Waves and for Acoustic Microdevices Transient Acoustic Wave Propagation in Porous media <sup>27</sup>

where � *λ*1(*t*) and � *λ*2(*t*) are operators of time corresponding to the eigenvalues of the system (41). Their expressions are given by:

where

where

*F*3(*ω*) = *ρ <sup>f</sup> c*<sup>0</sup>

**Inverse problem**

transmitted *p<sup>t</sup>*

*<sup>τ</sup>*<sup>5</sup> = *R*�

*Fi*(*ω*) = (<sup>1</sup> <sup>+</sup> *<sup>φ</sup>* (�*i*(*ω*) <sup>−</sup> <sup>1</sup>))

*<sup>F</sup>*1(*ω*) cosh

*ρ*<sup>11</sup> − *Q*�

<sup>T</sup> (*ω*), which is the Fourier transform of *<sup>T</sup>*˜ is given by [42]

*l λ*1(*ω*) 

Ψ(*ω*) = 2(*Z*1(*ω*)*Z*4(*ω*) − *Z*2(*ω*)*Z*3(*ω*)),

and the coefficients *Z*1(*ω*), *Z*2(*ω*), *Z*3(*ω*) and *Z*4(*ω*) by

*<sup>ρ</sup>*<sup>12</sup> *<sup>τ</sup>*<sup>6</sup> = *<sup>A</sup>*

For a slab of cancellous bone occupying the region 0 <sup>≤</sup> *<sup>x</sup>* <sup>≤</sup> *<sup>L</sup>*, the incident *<sup>p</sup><sup>i</sup>*

<sup>T</sup> (*ω*) = *<sup>j</sup>ω*2*<sup>ρ</sup> <sup>f</sup> <sup>c</sup>*0*F*4(*ω*)

*<sup>λ</sup>i*(*ω*) <sup>Ψ</sup>*i*(*ω*) sinh *l λi*(*ω*)

The functions *λ*1(*ω*), *λ*2(*ω*), �1(*ω*) and �2(*ω*) are the Fourier transform of

and �2(*t*), respectively. The coefficients Ψ1(*ω*), Ψ2(*ω*) and Ψ(*ω*) are given by

*R*� + *Q*� 

*T*˜ given by Eq. 31. In the frequency domain, the expression of the transmission coefficient

(*jωρ <sup>f</sup> <sup>c</sup>*0*F*4(*ω*))<sup>2</sup> <sup>−</sup> (*jωF*3(*ω*) <sup>−</sup> <sup>1</sup>)

<sup>+</sup>*F*2(*ω*) cosh

Ψ1(*ω*) = *φZ*2(*ω*) − (1 − *φ*)*Z*4(*ω*), Ψ2(*ω*)=(1 − *φ*)*Z*3(*ω*) − *φZ*1(*ω*),

*Z*1(*ω*)=(*P* + *Q*�1(*ω*))*λ*1(*ω*), *Z*2(*ω*)=(*P* + *Q*�2(*ω*))*λ*2(*ω*), *Z*3(*ω*)=(*Q* + *R*�1(*ω*))*λ*1(*ω*), *Z*4(*ω*)=(*Q* + *R*�2(*ω*))*λ*2(*ω*).

**3.2. Ultrasonic characterization of human cancellous bone using the Biot theory:**

Osteoporosis is a degenerative bone disease associated with biochemical and hormonal changes in the ageing body. These changes perturb the equilibrium between bone apposition and bone removal, resulting in a net decrease in bone mass. This leads to a modification of the structure (porosity, trabecular thickness, connectivity etc.) and, to a lesser extent, the composition (mineral density) of the bone. These changes result in a decrease of the mechanical strength of bone and in an increase of the risk of fracture. Osteoporosis mainly affects the trabecular bone (located at the hip, vertebrae or heels for instance). Early clinical detection of this pathological condition is very important to insure proper treatment. The primary method currently used for clinical bone assessment is based on x-ray absorptiometry, and measures total bone mass at a particular anatomic site. Because other factors, such as architecture, also appear to have a role in determining an individual's risk of fracture, ultrasound is an alternative to X-rays that has generated much attention. In addition to their potential for conveying the architectural aspects of bone, ultrasonic techniques also may have advantages in view of their use as a nonionising radiation and inherently lower costs,

, *<sup>τ</sup>*<sup>7</sup> =

(*t*) fields are related in the time domain by the transmission scattering operator

*l λ*2(*ω*)

2 Ψ(*ω*) *R*�

*ρ*<sup>12</sup> − *Q*�

, *i* = 1, 2.

*<sup>ρ</sup>*22 .

Transient Acoustic Wave Propagation in Porous Media 153

<sup>2</sup> , (43)

, *<sup>F</sup>*4(*ω*)= *<sup>F</sup>*1(*ω*)+*F*2(*ω*).

*λ*1(*t*),

*λ*2(*t*), �1(*t*)

(*t*) and

$$\begin{split} \widetilde{\lambda}\_{1}(t) &= \frac{1}{2} \left( \tau\_{1} - \sqrt{\tau\_{1}^{2} - 4\tau\_{3}} \right) \frac{\partial^{2}}{\partial t^{2}} + \frac{1}{2} \left( \tau\_{2} - \frac{\tau\_{1}\tau\_{2} - 2\tau\_{4}}{\sqrt{\tau\_{1}^{2} - 4\tau\_{3}}} \right) \frac{\partial^{3/2}}{\partial t^{3/2}} \\ &- \frac{1}{4} \left( \frac{\tau\_{2}^{2}}{\sqrt{\tau\_{1}^{2} - 4\tau\_{3}}} - \frac{(\tau\_{1}\tau\_{2} - 2\tau\_{4})^{2}}{2(\tau\_{1}^{2} - 4\tau\_{3})^{3/2}} \right) \frac{\partial}{\partial t}, \\\\ &\qquad \cdot \frac{1}{4} \left( \tau\_{1} - \boxed{\tau\_{1}^{2} - \tau\_{2}} - \tau\_{1} \left( \tau\_{2} - \tau\_{3} - 2\tau\_{4} \right) \right) \partial^{3/2}. \end{split}$$

$$\begin{split} \widetilde{\lambda}\_{2}(t) &= \frac{1}{2} \left( \tau\_{1} + \sqrt{\tau\_{1}^{2} - 4\tau\_{3}} \right) \frac{\partial^{2}}{\partial t^{2}} + \frac{1}{2} \left( \tau\_{2} + \frac{\tau\_{1}\tau\_{2} - 2\tau\_{4}}{\sqrt{\tau\_{1}^{2} - 4\tau\_{3}}} \right) \frac{\partial^{3/2}}{\partial t^{3/2}} \\ &+ \frac{1}{4} \left( \frac{\tau\_{2}^{2}}{\sqrt{\tau\_{1}^{2} - 4\tau\_{3}}} - \frac{(\tau\_{1}\tau\_{2} - 2\tau\_{4})^{2}}{2(\tau\_{1}^{2} - 4\tau\_{3})^{3/2}} \right) \frac{\partial}{\partial t} \end{split}$$

with:

$$\begin{aligned} \tau\_1 &= \mathbf{R}' \rho\_{11} + \mathbf{P}' \rho\_{22} - 2\mathbf{Q}' \rho\_{12}, \quad \tau\_2 = A \left( \mathbf{P}' + \mathbf{R}' + 2\mathbf{Q}' \right), \quad \tau\_3 = \left( \mathbf{P}' \mathbf{R}' - \mathbf{Q}'^2 \right) \left( \rho\_{11} \rho\_{22} - \rho\_{12}^2 \right), \\ \text{and} \quad \tau\_4 &= A \left( \mathbf{P}' \mathbf{R}' - \mathbf{Q}'^2 \right) \left( \rho\_{11} + \rho\_{22} - 2\rho\_{12} \right). \end{aligned}$$

Coefficients *R*� , *P*� and *Q*� are given by: *R*� = *<sup>R</sup> PR*−*Q*<sup>2</sup> , *<sup>Q</sup>*� <sup>=</sup> *<sup>Q</sup> PR*−*Q*<sup>2</sup> , and *<sup>P</sup>*� <sup>=</sup> *<sup>P</sup> PR*−*Q*<sup>2</sup> . The system of equations (41) shows that the slow and fast waves obey the same propagation equation developed in the framework of the equivalent fluid model (equation 13). The eigenvectors (1, ��1(*t*)) and (1, ��2(*t*)) associated with the eigenvalues � *λ*1(*t*) et � *λ*2(*t*) link the potential solid and fluid Φ*<sup>s</sup>* and Φ*<sup>f</sup>* , respectively to slow and fast waves Φ<sup>1</sup> and Φ<sup>2</sup> by the following relations:

$$
\begin{pmatrix} \Phi\_{\mathbf{s}} \\ \Phi\_{f} \end{pmatrix} = \begin{pmatrix} 1 & 1 \\ \widetilde{\otimes}\_{1}(t) \ \widetilde{\otimes}\_{2}(t) \end{pmatrix} \begin{pmatrix} \Phi\_{1} \\ \Phi\_{2} \end{pmatrix} . \tag{42}
$$

the analytical expressions of temporal operators ��1(*t*) and ��2(*t*) are given by:

$$
\begin{split}
\tilde{\mathbb{G}}\_{1}(t) &= \frac{\tau\_{1} - 2\tau\_{5} - \sqrt{\tau\_{1}^{2} - 4\tau\_{3}}}{2\tau\_{7}} + \frac{t^{-1/2}}{4\tau\_{7}^{2}\sqrt{\pi}} \times \\
&\left[ \left( \tau\_{2} - 2\tau\_{6} - \frac{\tau\_{1}\tau\_{2} - 2\tau\_{4}}{\sqrt{\tau\_{1}^{2} - 4\tau\_{3}}} \right) 2\tau\_{7} + \left( \tau\_{1} - 2\tau\_{5} - \sqrt{\tau\_{1}^{2} - 4\tau\_{3}} \right) 2\tau\_{6} \right].
\end{split}
$$

$$\begin{aligned} \widetilde{s}\_2(t) &= \frac{\tau\_1 - 2\tau\_5 + \sqrt{\tau\_1^2 - 4\tau\_3}}{2\tau\_7} + \frac{t^{-1/2}}{4\tau\_7^2\sqrt{\pi}} \times \\ &\left[ \left( \tau\_2 - 2\tau\_6 + \frac{\tau\_1\tau\_2 - 2\tau\_4}{\sqrt{\tau\_1^2 - 4\tau\_3}} \right) 2\tau\_7 + \left( \tau\_1 - 2\tau\_5 + \sqrt{\tau\_1^2 - 4\tau\_3} \right) 2\tau\_6 \right]. \end{aligned}$$

where

26 Will-be-set-by-IN-TECH

� *∂*<sup>2</sup> *<sup>∂</sup>t*<sup>2</sup> <sup>+</sup>

<sup>−</sup> (*τ*1*τ*<sup>2</sup> <sup>−</sup> <sup>2</sup>*τ*4)<sup>2</sup> 2(*τ*<sup>2</sup>

> � *∂*<sup>2</sup> *<sup>∂</sup>t*<sup>2</sup> <sup>+</sup>

<sup>−</sup> (*τ*1*τ*<sup>2</sup> <sup>−</sup> <sup>2</sup>*τ*4)<sup>2</sup> 2(*τ*<sup>2</sup>

*P*� + *R*� + 2*Q*�

The system of equations (41) shows that the slow and fast waves obey the same propagation equation developed in the framework of the equivalent fluid model (equation 13). The

potential solid and fluid Φ*<sup>s</sup>* and Φ*<sup>f</sup>* , respectively to slow and fast waves Φ<sup>1</sup> and Φ<sup>2</sup> by the

� 1 1 ��1(*t*) ��2(*t*)

<sup>1</sup> − <sup>4</sup>*τ*3)3/2

<sup>1</sup> − <sup>4</sup>*τ*3)3/2

1 2 ⎛

� , *τ*<sup>3</sup> =

*λ*2(*t*) are operators of time corresponding to the eigenvalues of the system

<sup>⎝</sup>*τ*<sup>2</sup> − �

⎞ ⎠ *∂ ∂t* ,

<sup>⎝</sup>*τ*<sup>2</sup> <sup>+</sup> �

⎞ ⎠ *∂ ∂t*

*PR*−*Q*<sup>2</sup> , *<sup>Q</sup>*� <sup>=</sup> *<sup>Q</sup>*

� � <sup>Φ</sup><sup>1</sup> Φ<sup>2</sup> �

*τ*1*τ*<sup>2</sup> − 2*τ*<sup>4</sup> *τ*2 <sup>1</sup> − 4*τ*<sup>3</sup>

*τ*1*τ*<sup>2</sup> − 2*τ*<sup>4</sup> *τ*2 <sup>1</sup> − 4*τ*<sup>3</sup>

> � *P*�

⎞ ⎠

⎞ ⎠

*<sup>R</sup>*� <sup>−</sup> *<sup>Q</sup>*�<sup>2</sup>

*∂*3/2 *∂t*3/2

� �

*PR*−*Q*<sup>2</sup> , and *<sup>P</sup>*� <sup>=</sup> *<sup>P</sup>*

*λ*1(*t*) et �

, (42)

� 2*τ*<sup>6</sup> ⎤ ⎦ ,

� 2*τ*<sup>6</sup> ⎤ ⎦ ,

*<sup>ρ</sup>*11*ρ*<sup>22</sup> <sup>−</sup> *<sup>ρ</sup>*<sup>2</sup>

12 � ,

*PR*−*Q*<sup>2</sup> .

*λ*2(*t*) link the

*∂*3/2 *∂t*3/2

1 2 ⎛

where �

with:

*τ*<sup>1</sup> =*R*�

and *τ*<sup>4</sup> = *A*

following relations:

��1(*t*) =

��2(*t*) =

Coefficients *R*�

*λ*1(*t*) and �

� *<sup>λ</sup>*1(*t*) = <sup>1</sup>

� *<sup>λ</sup>*2(*t*) = <sup>1</sup>

*ρ*11+*P*�

� *P*�

(41). Their expressions are given by:

2 � *τ*<sup>1</sup> − � *τ*2 <sup>1</sup> − 4*τ*<sup>3</sup>

⎛ ⎝

�

*τ*2

*τ*2

*<sup>ρ</sup>*12, *<sup>τ</sup>*<sup>2</sup> = *<sup>A</sup>* �

, *P*� and *Q*� are given by: *R*� = *<sup>R</sup>*

� Φ*<sup>s</sup>* Φ*f* � =

> � *τ*2 <sup>1</sup> − 4*τ*<sup>3</sup>

� *τ*2 <sup>1</sup> − 4*τ*<sup>3</sup>

2*τ*<sup>7</sup>

<sup>⎝</sup>*τ*<sup>2</sup> <sup>−</sup> <sup>2</sup>*τ*<sup>6</sup> <sup>+</sup> �

2*τ*<sup>7</sup>

<sup>⎝</sup>*τ*<sup>2</sup> <sup>−</sup> <sup>2</sup>*τ*<sup>6</sup> <sup>−</sup> �

(*ρ*11+*ρ*22−2*ρ*12).

eigenvectors (1, ��1(*t*)) and (1, ��2(*t*)) associated with the eigenvalues �

the analytical expressions of temporal operators ��1(*t*) and ��2(*t*) are given by:

*τ*1*τ*<sup>2</sup> − 2*τ*<sup>4</sup> *τ*2 <sup>1</sup> − 4*τ*<sup>3</sup>

*τ*1*τ*<sup>2</sup> − 2*τ*<sup>4</sup> *τ*2 <sup>1</sup> − 4*τ*<sup>3</sup>

+ *t* −1/2 4*τ*<sup>2</sup> 7 <sup>√</sup>*<sup>π</sup>* ×

+

⎞

⎠ 2*τ*<sup>7</sup> +

*t*−1/2 4*τ*<sup>2</sup> 7 <sup>√</sup>*<sup>π</sup>* ×

⎠ 2*τ*<sup>7</sup> +

⎞

�

�

*τ*<sup>1</sup> − 2*τ*<sup>5</sup> −

*τ*<sup>1</sup> − 2*τ*<sup>5</sup> +

� *τ*2 <sup>1</sup> − 4*τ*<sup>3</sup>

� *τ*2 <sup>1</sup> − 4*τ*<sup>3</sup>

2 *τ*2 <sup>1</sup> − 4*τ*<sup>3</sup>

2 *τ*2 <sup>1</sup> − 4*τ*<sup>3</sup>

− 1 4

> 2 � *τ*<sup>1</sup> + � *τ*2 <sup>1</sup> − 4*τ*<sup>3</sup>

⎛ ⎝

�

�

+ 1 4

*ρ*22−2*Q*�

*<sup>R</sup>*� <sup>−</sup> *<sup>Q</sup>*�<sup>2</sup>

*τ*<sup>1</sup> − 2*τ*<sup>5</sup> −

*τ*<sup>1</sup> − 2*τ*<sup>5</sup> +

⎡ ⎣ ⎛

⎡ ⎣ ⎛

$$\pi\_5 = \left(\mathbb{R}'\rho\_{11} - \mathbb{Q}'\rho\_{12}\right) \quad \pi\_6 = A\left(\mathbb{R}' + \mathbb{Q}'\right), \quad \pi\_7 = \left(\mathbb{R}'\rho\_{12} - \mathbb{Q}'\rho\_{22}\right) \dots$$

For a slab of cancellous bone occupying the region 0 <sup>≤</sup> *<sup>x</sup>* <sup>≤</sup> *<sup>L</sup>*, the incident *<sup>p</sup><sup>i</sup>* (*t*) and transmitted *p<sup>t</sup>* (*t*) fields are related in the time domain by the transmission scattering operator *T*˜ given by Eq. 31. In the frequency domain, the expression of the transmission coefficient <sup>T</sup> (*ω*), which is the Fourier transform of *<sup>T</sup>*˜ is given by [42]

$$\mathcal{T}(\omega) = \frac{j\omega 2\rho\_f c\_0 \mathcal{F}\_4(\omega)}{\left(j\omega \rho\_f c\_0 \mathcal{F}\_4(\omega)\right)^2 - \left(j\omega \mathcal{F}\_5(\omega) - 1\right)^2},\tag{43}$$

where

$$F\_i(\omega) = \left(1 + \phi\left(\Im\_i(\omega) - 1\right)\right) \sqrt{\lambda\_i(\omega)} \frac{\Psi\_i(\omega)}{\sinh\left(l\sqrt{\lambda\_i(\omega)}\right)} \frac{2}{\Psi(\omega)}, \qquad i = 1, 2.$$

$$F\_3(\omega) = \rho\_f c\_0 \left(F\_1(\omega) \cosh\left(l\sqrt{\lambda\_1(\omega)}\right) + F\_2(\omega) \cosh\left(l\sqrt{\lambda\_2(\omega)}\right)\right), \ F\_4(\omega) = F\_1(\omega) + F\_2(\omega).$$

The functions *λ*1(*ω*), *λ*2(*ω*), �1(*ω*) and �2(*ω*) are the Fourier transform of *λ*1(*t*), *λ*2(*t*), �1(*t*) and �2(*t*), respectively. The coefficients Ψ1(*ω*), Ψ2(*ω*) and Ψ(*ω*) are given by

$$\begin{aligned} \Psi\_1(\omega) &= \phi Z\_2(\omega) - (1 - \phi) Z\_4(\omega), & \Psi\_2(\omega) &= (1 - \phi) Z\_3(\omega) - \phi Z\_1(\omega), \\ \Psi(\omega) &= 2(Z\_1(\omega) Z\_4(\omega) - Z\_2(\omega) Z\_3(\omega)), \end{aligned}$$

and the coefficients *Z*1(*ω*), *Z*2(*ω*), *Z*3(*ω*) and *Z*4(*ω*) by

$$\begin{aligned} Z\_1(\omega) &= (P + Q \circledS\_1(\omega))\lambda\_1(\omega), & Z\_2(\omega) &= (P + Q \circledS\_2(\omega))\lambda\_2(\omega), \\ Z\_3(\omega) &= (Q + R \circledS\_1(\omega))\lambda\_1(\omega), & Z\_4(\omega) &= (Q + R \circledS\_2(\omega))\lambda\_2(\omega). \end{aligned}$$

### **3.2. Ultrasonic characterization of human cancellous bone using the Biot theory: Inverse problem**

Osteoporosis is a degenerative bone disease associated with biochemical and hormonal changes in the ageing body. These changes perturb the equilibrium between bone apposition and bone removal, resulting in a net decrease in bone mass. This leads to a modification of the structure (porosity, trabecular thickness, connectivity etc.) and, to a lesser extent, the composition (mineral density) of the bone. These changes result in a decrease of the mechanical strength of bone and in an increase of the risk of fracture. Osteoporosis mainly affects the trabecular bone (located at the hip, vertebrae or heels for instance). Early clinical detection of this pathological condition is very important to insure proper treatment. The primary method currently used for clinical bone assessment is based on x-ray absorptiometry, and measures total bone mass at a particular anatomic site. Because other factors, such as architecture, also appear to have a role in determining an individual's risk of fracture, ultrasound is an alternative to X-rays that has generated much attention. In addition to their potential for conveying the architectural aspects of bone, ultrasonic techniques also may have advantages in view of their use as a nonionising radiation and inherently lower costs,

#### 28 Will-be-set-by-IN-TECH 154 Modeling and Measurement Methods for Acoustic Waves and for Acoustic Microdevices Transient Acoustic Wave Propagation in Porous media <sup>29</sup>

compared with x-ray densitometric methods. Although ultrasonic methods appear promising for noninvasive bone assessment, they have not yet fulfilled their potential. Unfortunately, a poor understanding of the ultrasound interaction with bone has become one of the obstacles preventing it from being a fully developed diagnostic technique. Despite extensive research on the empirical relationship between ultrasound and the bulk properties of bone, the mechanism of how ultrasound physically interacts with bone is still unclear. Since trabecular bone is an inhomogeneous porous medium, the interaction between ultrasound and bone is a highly complex phenomenon. Modelling ultrasonic propagation through trabecular tissue has been considered using porous media theories, such as Biot's theory. As seen in the previous section, within the framework of the modified Biot theory, the propagation of ultrasonic waves in a slab of cancellous bone is conditioned by many parameters: porosity *φ*, tortuosity *α*∞, viscous characteristic length Λ, fluid viscosity *η*, Young's modulus of the elastic solid *Es*, Young's modulus of porous skeletal frame *Eb*, Poisson's ratio of the elastic solid *νs*, Poisson's ratio of the porous skeletal frame *νb*, the solid density *ρs*, the bulk modulus of the saturating fluid *Kf* and the fluid density *ρ <sup>f</sup>* . It is therefore important to develop new experimental methods and efficient tools for their estimation. The basic inverse problem associated with the slab of cancellous bone may be stated as follows: from measurements of the signal transmitted outside the slab, find the values of the medium's parameters. Solving the inverse problem for all the Biot parameters using only the transmitted experimental data is difficult, if not impossible. To achieve this task, requires more experimental data for obtaining a unique solution. For this reason, in this contribution we limit the inversion to the five parameters: *Eb*, *νb*, *φ*, *α*<sup>∞</sup> and Λ. In our previous paper [42], we studied the sensitivity of transmitted waveforms to variations of *φ*, *α*<sup>∞</sup> and Λ. The sensitivity of *Eb* and *ν<sup>b</sup>* has been examined in Ref. [43]

wherein *p<sup>t</sup>*

and *p<sup>t</sup>*

*exp*(*x*, *ti*)*i*=1,2,...*<sup>n</sup>* is the discrete set of values of the experimental transmitted signal

(*x*, *ti*)*i*=1,2,...*<sup>n</sup>* the discrete set of values of the simulated transmitted signal. A large variation range is applied of each estimating parameter value in solving the inverse problem. The variation range of the parameters is ; *α*<sup>∞</sup> ∈ [1, 2], Λ ∈ [1, 200]*μm*, *φ* ∈ [0.5, 0.99], *ν* ∈ [0.1, 0.5] and *Eb* ∈ [0.5, 5]GPa. The variations of the cost function with the physical parameters present one clear minimum corresponding to the mathematical solution of the inverse problem. This shows that the inverse problem is well posed mathematically, and that the solution is unique. The minima, corresponding to the solution of the inverse problem, are clearly observed for each parameter. After solving the inverse problem, we find the following optimized values: *φ* = 0.64, *α*<sup>∞</sup> = 1.018, Λ = 10, 44*μm*, *ν<sup>b</sup>* = 0.28 and *Eb* = 4.49 GPa. Using these values, we present in Figs. 16-19 the variations in the discrepancy function U with respect to two values of the inverted parameters. To show clearly the solution of the

8

**Figure 16.** Variation of the minimization function *U* with the viscous characteristic length Λ and the

4

**Figure 17.** Variation of the minimization function *U* with the the Young modulus of the skeletal frame

inverse problem, the variation of U in Figs. 16-19 is given only around the minima values of the inverted parameters. In Fig. 20, a comparison is made between the experimental transmitted signal and the simulated transmitted signals using the reconstructed values of

4.5

5

Eb (GPa)

−100.4

−100.5

Log (U)

5

(GPa) <sup>Λ</sup> (μm)

4

Transient Acoustic Wave Propagation in Porous Media 155

4.5

<sup>1</sup> 1.02 1.04 1.06

α∞

Eb

10

12

−100.4

−100.5

Log(U)

Young modulus of the skeletal frame *Eb*.

*Eb* and the tortuosity *α*∞.

As an application of this model, some numerical simulations are compared with experimental results. Experiments are performed in water using two broadband Panametrics A 303S plane piezoelectric transducers with a central frequency of 1 MHz in water, and diameter of 1cm. 400 V Pulses are provided by a 5058PR Panametrics pulser/receiver. Electronic interference is removed by averaging 1000 acquisitions. The experimental setup is shown in Fig. 9. The parallel-faced samples were machined from femoral heads and femoral necks of human cancellous bone. The liquid in the pore space (blood and marrow) is removed from the bone sample and substituted by water. The size of the ultrasound beam is very small compared to the size of the specimens. The emitting transducer insonifies the sample at normal incidence with a short (in time domain) pulse. When the pulse hits the front surface of the sample, a part is reflected, a part is transmitted as a fast wave, and a part is transmitted as a slow wave. When any of these components, travelling at different speeds, hit the second surface, a similar effect takes place: a part is transmitted into the fluid, and a part is reflected as a fast or slow wave. The experimental transmitted waveforms are travelling through the cancellous bone in the same direction as the trabecular alignment (x direction). The fluid characteristics are: bulk modulus *Kf* = 2.28 GPa, density *ρ <sup>f</sup>* = 1000Kgm−3, viscosity *η* = 10−3*Kg*.*m*.*s*−1. Consider a sample of human cancellous bone M1 (femoral neck) of thickness 11.2 mm and solid density *ρ<sup>s</sup>* = 1990Kgm−3. The Young's modulus *Es* = 13*GPa* and Poisson ratio *ν<sup>s</sup>* = 0.3 of the solid bone are taken from the literature [48]. Figure 9 shows the experimental incident signal. The inverse problem is solved by minimizing the function *U*(*φ*, *α*∞, Λ, *Eb*, *νs*) given by:

$$\mathcal{U}(\phi\_{\prime}\alpha\_{\infty\prime}\Lambda\_{\prime}E\_{b\prime}\nu\_{b}) = \sum\_{i=1}^{i=n} (p^{t}\_{\exp}(\chi\_{\prime}t\_{i}) - p^{t}(\chi\_{\prime}t\_{i}))^{2}\Lambda\_{\prime}$$

wherein *p<sup>t</sup> exp*(*x*, *ti*)*i*=1,2,...*<sup>n</sup>* is the discrete set of values of the experimental transmitted signal and *p<sup>t</sup>* (*x*, *ti*)*i*=1,2,...*<sup>n</sup>* the discrete set of values of the simulated transmitted signal. A large variation range is applied of each estimating parameter value in solving the inverse problem. The variation range of the parameters is ; *α*<sup>∞</sup> ∈ [1, 2], Λ ∈ [1, 200]*μm*, *φ* ∈ [0.5, 0.99], *ν* ∈ [0.1, 0.5] and *Eb* ∈ [0.5, 5]GPa. The variations of the cost function with the physical parameters present one clear minimum corresponding to the mathematical solution of the inverse problem. This shows that the inverse problem is well posed mathematically, and that the solution is unique. The minima, corresponding to the solution of the inverse problem, are clearly observed for each parameter. After solving the inverse problem, we find the following optimized values: *φ* = 0.64, *α*<sup>∞</sup> = 1.018, Λ = 10, 44*μm*, *ν<sup>b</sup>* = 0.28 and *Eb* = 4.49 GPa. Using these values, we present in Figs. 16-19 the variations in the discrepancy function U with respect to two values of the inverted parameters. To show clearly the solution of the

28 Will-be-set-by-IN-TECH

compared with x-ray densitometric methods. Although ultrasonic methods appear promising for noninvasive bone assessment, they have not yet fulfilled their potential. Unfortunately, a poor understanding of the ultrasound interaction with bone has become one of the obstacles preventing it from being a fully developed diagnostic technique. Despite extensive research on the empirical relationship between ultrasound and the bulk properties of bone, the mechanism of how ultrasound physically interacts with bone is still unclear. Since trabecular bone is an inhomogeneous porous medium, the interaction between ultrasound and bone is a highly complex phenomenon. Modelling ultrasonic propagation through trabecular tissue has been considered using porous media theories, such as Biot's theory. As seen in the previous section, within the framework of the modified Biot theory, the propagation of ultrasonic waves in a slab of cancellous bone is conditioned by many parameters: porosity *φ*, tortuosity *α*∞, viscous characteristic length Λ, fluid viscosity *η*, Young's modulus of the elastic solid *Es*, Young's modulus of porous skeletal frame *Eb*, Poisson's ratio of the elastic solid *νs*, Poisson's ratio of the porous skeletal frame *νb*, the solid density *ρs*, the bulk modulus of the saturating fluid *Kf* and the fluid density *ρ <sup>f</sup>* . It is therefore important to develop new experimental methods and efficient tools for their estimation. The basic inverse problem associated with the slab of cancellous bone may be stated as follows: from measurements of the signal transmitted outside the slab, find the values of the medium's parameters. Solving the inverse problem for all the Biot parameters using only the transmitted experimental data is difficult, if not impossible. To achieve this task, requires more experimental data for obtaining a unique solution. For this reason, in this contribution we limit the inversion to the five parameters: *Eb*, *νb*, *φ*, *α*<sup>∞</sup> and Λ. In our previous paper [42], we studied the sensitivity of transmitted waveforms to variations of *φ*, *α*<sup>∞</sup> and Λ. The sensitivity of *Eb* and *ν<sup>b</sup>* has been examined in

As an application of this model, some numerical simulations are compared with experimental results. Experiments are performed in water using two broadband Panametrics A 303S plane piezoelectric transducers with a central frequency of 1 MHz in water, and diameter of 1cm. 400 V Pulses are provided by a 5058PR Panametrics pulser/receiver. Electronic interference is removed by averaging 1000 acquisitions. The experimental setup is shown in Fig. 9. The parallel-faced samples were machined from femoral heads and femoral necks of human cancellous bone. The liquid in the pore space (blood and marrow) is removed from the bone sample and substituted by water. The size of the ultrasound beam is very small compared to the size of the specimens. The emitting transducer insonifies the sample at normal incidence with a short (in time domain) pulse. When the pulse hits the front surface of the sample, a part is reflected, a part is transmitted as a fast wave, and a part is transmitted as a slow wave. When any of these components, travelling at different speeds, hit the second surface, a similar effect takes place: a part is transmitted into the fluid, and a part is reflected as a fast or slow wave. The experimental transmitted waveforms are travelling through the cancellous bone in the same direction as the trabecular alignment (x direction). The fluid characteristics are: bulk modulus *Kf* = 2.28 GPa, density *ρ <sup>f</sup>* = 1000Kgm−3, viscosity *η* = 10−3*Kg*.*m*.*s*−1. Consider a sample of human cancellous bone M1 (femoral neck) of thickness 11.2 mm and solid density *ρ<sup>s</sup>* = 1990Kgm−3. The Young's modulus *Es* = 13*GPa* and Poisson ratio *ν<sup>s</sup>* = 0.3 of the solid bone are taken from the literature [48]. Figure 9 shows the experimental incident signal. The

inverse problem is solved by minimizing the function *U*(*φ*, *α*∞, Λ, *Eb*, *νs*) given by:

*i*=*n* ∑ *i*=1 (*pt*

*exp*(*x*, *ti*) <sup>−</sup> *<sup>p</sup><sup>t</sup>*

(*x*, *ti*))2,

*U*(*φ*, *α*∞, Λ, *Eb*, *νb*) =

Ref. [43]

**Figure 16.** Variation of the minimization function *U* with the viscous characteristic length Λ and the Young modulus of the skeletal frame *Eb*.

**Figure 17.** Variation of the minimization function *U* with the the Young modulus of the skeletal frame *Eb* and the tortuosity *α*∞.

inverse problem, the variation of U in Figs. 16-19 is given only around the minima values of the inverted parameters. In Fig. 20, a comparison is made between the experimental transmitted signal and the simulated transmitted signals using the reconstructed values of

30 Will-be-set-by-IN-TECH 156 Modeling and Measurement Methods for Acoustic Waves and for Acoustic Microdevices Transient Acoustic Wave Propagation in Porous media <sup>31</sup>

**Figure 18.** Variation of the minimization function *U* with the viscous characteristic length Λ and the tortuosity *α*∞.

9 9.5 10

Time (s)

88 90 92 94 96 98 100

21) reproduce correctly the experimental transmitted signals. This leads us to conclude that

*LMA, CNRS, UPR 7051, Aix-Marseille Univ, Centrale Marseille, F-13402 Marseille Cedex 20, France*

*Laboratoire de Physique Théorique, Faculté de Physique, USTHB, BP 32 El Alia, Bab-Ezzouar, 16211,*

*Universite du Maine UMR CNRS 6613 Laboratoire d'Acoustique de l'Universite du Maine UFR STS*

**Figure 21.** Comparison between the experimental transmitted signal (solid line) and the simulated transmitted signals (dashed line) using the reconstructed values of *α*∞, *φ*, Λ, *ν<sup>b</sup>* and *Eb* (sample M2).

this method is well adapted for the characterization of cancellous bone.

Time (μs)

**Figure 20.** Comparison between the experimental transmitted signal (solid line) and the simulated transmitted signals (dashed line) using the reconstructed values of *α*∞, *φ*, Λ, *ν<sup>b</sup>* and *Eb* (sample M1).

−0.2 −0.15 −0.1 −0.05 0 0.05 0.1 0.15 0.2

−0.3

*Avenue O. Messiaen 72085 Le Mans CEDEX 09 France*

**Author details**

Mohamed Fellah

Claude Depollier

*Alger- Algeria*

Zine El Abiddine Fellah

−0.2

−0.1

Amplitude (V)

0

0.1

0.2

0.3

Amplitude (V)

x 10 −5

Transient Acoustic Wave Propagation in Porous Media 157

**Figure 19.** Variation of the minimization function *U* with the porosity *φ* and the viscous characteristic length Λ.

*α*∞, *φ*, Λ, *ν<sup>b</sup>* and *Eb*. The difference between the two curves is small, which leads us to conclude that the optimized values of the physical parameters are correct. The fast and slow waves predicted by the Biot theory are easily detected in the transmitted signal. The slow wave seems to be less attenuated than the fast wave. In the other applications, the slow wave is generally more-attenuated and dispersive than the fast wave. We usually observe the opposite phenomena for cancellous bone samples; this can be explained by the different orders of magnitude of the physical parameters (high porosity, low tortuosity...etc). Let us now solve the inverse problem for sample M2 (femoral neck) of thickness 12 mm. Using another sample of cancellous bone (femoral head) M2 of thickness 10.2 mm. The results after solving the inverse problem are: *φ* = 0.72, *α*<sup>∞</sup> = 1.1, Λ = 14.97*μm*, *ν<sup>b</sup>* = 0.22 and *Eb* = 3.1*GPa*. In Fig. 21, we compare the experimental transmitted signal to the transmitted simulated signal using reconstruct values of the physical parameters. The correlation between the two curves is excellent. The simulated transmitted signals obtained using optimized values (Figs. 20 and

**Figure 20.** Comparison between the experimental transmitted signal (solid line) and the simulated transmitted signals (dashed line) using the reconstructed values of *α*∞, *φ*, Λ, *ν<sup>b</sup>* and *Eb* (sample M1).

**Figure 21.** Comparison between the experimental transmitted signal (solid line) and the simulated transmitted signals (dashed line) using the reconstructed values of *α*∞, *φ*, Λ, *ν<sup>b</sup>* and *Eb* (sample M2).

21) reproduce correctly the experimental transmitted signals. This leads us to conclude that this method is well adapted for the characterization of cancellous bone.

## **Author details**

30 Will-be-set-by-IN-TECH

1

α∞ <sup>Λ</sup>(μm)

8

**Figure 18.** Variation of the minimization function *U* with the viscous characteristic length Λ and the

0.55 0.6 0.65 0.7 0.75 <sup>8</sup>

**Figure 19.** Variation of the minimization function *U* with the porosity *φ* and the viscous characteristic

*α*∞, *φ*, Λ, *ν<sup>b</sup>* and *Eb*. The difference between the two curves is small, which leads us to conclude that the optimized values of the physical parameters are correct. The fast and slow waves predicted by the Biot theory are easily detected in the transmitted signal. The slow wave seems to be less attenuated than the fast wave. In the other applications, the slow wave is generally more-attenuated and dispersive than the fast wave. We usually observe the opposite phenomena for cancellous bone samples; this can be explained by the different orders of magnitude of the physical parameters (high porosity, low tortuosity...etc). Let us now solve the inverse problem for sample M2 (femoral neck) of thickness 12 mm. Using another sample of cancellous bone (femoral head) M2 of thickness 10.2 mm. The results after solving the inverse problem are: *φ* = 0.72, *α*<sup>∞</sup> = 1.1, Λ = 14.97*μm*, *ν<sup>b</sup>* = 0.22 and *Eb* = 3.1*GPa*. In Fig. 21, we compare the experimental transmitted signal to the transmitted simulated signal using reconstruct values of the physical parameters. The correlation between the two curves is excellent. The simulated transmitted signals obtained using optimized values (Figs. 20 and

9

10

−100.4

−100.5

Log (U)

tortuosity *α*∞.

length Λ.

−100.49

−100.52

Log (U)

−100.55

1.02

9

φ Λ (μm)

10

1.04

Zine El Abiddine Fellah *LMA, CNRS, UPR 7051, Aix-Marseille Univ, Centrale Marseille, F-13402 Marseille Cedex 20, France*

Mohamed Fellah *Laboratoire de Physique Théorique, Faculté de Physique, USTHB, BP 32 El Alia, Bab-Ezzouar, 16211, Alger- Algeria*

Claude Depollier *Universite du Maine UMR CNRS 6613 Laboratoire d'Acoustique de l'Universite du Maine UFR STS Avenue O. Messiaen 72085 Le Mans CEDEX 09 France*

*LUNAM Universite du Maine. UMR CNRS 6613 Laboratoire d'Acoustique de l'Universite du MaineUFR STS Avenue O. Messiaen 72085 Le Mans CEDEX 09 France*

[19] Fellah ZEA, Depollier C, Fellah M (2001), Direct and inverse scattering problem in porous material having a rigid frame by fractional calculus based method, J. Sound.

Transient Acoustic Wave Propagation in Porous Media 159

[20] Fellah ZEA, Ogam E, Wirgin A, Fellah M, Depollier C et Lauriks W, Ultrasonic characterization of porous materials: Inverse problem, J. Sound and Vib. 203:746-759. [21] Leclaire P, Kelders L, Lauriks W, Glorieux C, Thoen J (1996), Determination of the viscous characteristic length in air-filled porous materials by ultrasonic attenuation

[22] Ayrault C, Moussatov A, Castagnède B, Lafarge D, (1999) Ultrasonic characterization of plastic foams via measurements with static pressure variations, Applied Physics Letters

[23] Moussatov A, Ayrault C, Castagnède B (2001), Porous material characterization ultrasonic method for estimation of tortuosity and characteristic length using a

[24] Leclaire P, Kelders L, Lauriks W, Brown NR, Melon M, Castagnède B (1996), Determination of viscous and thermal characteristic lengths of plastic foams by

[25] Fellah ZEA, Depollier C, Berger S, Lauriks W, Trompette P, Chapelon JY (2003), Detemination of transport parameters in air saturated porous materials via ultrasonic

[26] Fellah ZEA, Berger S, Lauriks W, Depollier C, Chapelon JY (2003), Inverse problem in air-saturated porous media via reflected waves, Rev. Sci. Inst. 74: (5), p. 2871. [27] Fellah ZEA, Berger S, Lauriks W, Depollier C, Aristégui C, Chapelon JY (2003), Measuring the porosity and the tortuosity of porous materials via reflected waves at

[28] Fellah ZEA, Berger S, Lauriks W, Depollier C, Trompette P, Chapelon JY (2003), Ultrasonic measuring of the porosity and tortuosity of air- saturated random packings

[29] Fellah ZEA, Mitri FG, Depollier D, Berger S, Lauriks W Chapelon JY (2003), Characterization of porous materials having a rigid frame via reflected waves, J. Appl.

[30] Fellah ZEA, Berger S, Lauriks W, Depollier C, Fellah M (2003), Measurement of the porosity of porous materials having a rigid frame via reflected waves : A time domain

[31] Fellah ZEA, Mitri FG, Fellah M, Ogam E, Depollier C (2007), Ultrasonic characterization of porous absorbing materials: Inverse problem, J. Sound. Vib. 302: 746-759. [32] Fellah ZEA, Wirgin A, Fellah M, Sebaa N, Depollier C, Lauriks W (2005), A time-domain model of transient acoustic wave propagation in double-layered porous media, J.

[33] Fellah ZEA, Fellah M, Sebaa N, Lauriks W, Depollier C (2006) , Measuring flow resistivity of porous materials at low frequencies range via acoustic transmitted waves,

[34] Fellah ZEA, Fellah M, Mitri FG, Sebaa N, Lauriks W, Depollier C (2007), Transient acoustic wave propagation in air-saturated porous media at low frequencies, J. Appl.

[35] Sebaa N, Fellah ZEA, Fellah M, Lauriks W, Depollier C (2005), Measuring flow resistivity

of porous material via acoustic reflected waves, J. Appl. Phys. 98: 084901.

ultrasonic measurements in helium and air, J. Appl. Phys. 80:, p. 2009.

reflected waves, J. Acoust. Soc. Am. 113: (5), 2561-2569.

oblique incidence, J. Acoust. Soc. Am. 113: (5), p. 2424.

analysis with fractional derivatives, J. App. Phys, 93: p. 296.

Vib. 244: (2), p. 359-366.

74: 3224.

measurements, J. Acoust. Soc. Am. 99: p. 1944.

barometric chamber, Ultrasonic 39: 195.

of beads, J. Appl. Phys. 93: p. 9352.

Acoust. Soc. Am. 118: 661-670.

J. Acoust. Soc. Am. 119: 1926-1928.

Phys. 94: p. 7914.

Phys. 102: 084906.

## **4. References**


[19] Fellah ZEA, Depollier C, Fellah M (2001), Direct and inverse scattering problem in porous material having a rigid frame by fractional calculus based method, J. Sound. Vib. 244: (2), p. 359-366.

32 Will-be-set-by-IN-TECH

*LUNAM Universite du Maine. UMR CNRS 6613 Laboratoire d'Acoustique de l'Universite du*

[1] Allard JF (1993), Propagation of Sound in Porous Media : Modeling Sound Absorbing

[2] Fellah ZEA, Depollier C (2000) Transient acoustic wave propagation in rigid porous

[3] Fellah ZEA, Depollier C, Fellah M, Lauriks W et Chpaleon JY (2005), Influence of dynamic tortuosity and compressibility on the propagation of transient waves in porous

[4] Fellah ZEA, Fellah M, Lauriks W, Depollier C (2003), Direct inverse scattering of transient acoustic waves by a slab of rigid porous material., J. Acoust. Soc. Am. 113:

[7] Allard JF, Champoux Y (1992), New empirical equations for sound propagation in rigid

[8] Champoux Y, Allard JF (1191), Dynamic tortuosity and bulk modulus in air-saturated

[9] Lafarge D, Lemarnier P, Allard JF, Tarnow V (1996), Dynamic compressibility of air in

[10] Norris AN (1986), On the viscodynamic operator in Biot's equations of poroelasticity, J.

[11] Pride SR, Morgan FD, Gangi AF (1993), Drag forces of porous media acoustics. Phys.

[12] Szabo TL (1994), Time domain wave equations for lossy media obeying a frequency

[13] Szabo TL (1995), Causal theories and data for acoustic attenuation obeying a frequency

[14] Norton V, Novarini JC (2003), Including dispersion and attenuation directly in time domain for wave propagation in isotropic media, J. Acoust. Soc. Am. 113: p. 3024. [15] Chen W, Holm S (2003), Modified Szabo's wave equation models for lossy media

[16] Samko SG, Kilbas AA, Marichev OI (1993), Fractional Integrals and Derivatives Theory

[17] Abramowitz M, Stegun TA (1968), Handbook of mathematical functions, National bureau of standards, U.S. Governement Printing Office, Washington, DC. Chap. 6, p.

[18] Fellah ZEA, Fellah M, Lauriks W, Depollier C, Angel Y, Chapelon JY (2003), Solution in time domain of ultrasonic propagation equation in porous material, Wave motion, 38:

porous structures at audible frequencies Journ. Acoust. Soc. Am, 102: 4.

power law, Journal of the Acoustical Society of America 97: p. 14.

obeying frequency power law, J. Acoust. Soc. Am. 113: p. 3024.

and Applications, Gordon and Breach Publishers, Amsterdam.

[5] Zwikker C, Kosten CW (1949) Sound Absorbing Materials, Elsevier , New York . [6] Johnson D.L, Koplik J, Dashen R (1987) Theory of dynamic permeability and tortuosity

in fluid-saturated porous media. J. Fluid. Mech. 176: 379-402.

frame fibrous materials . Acoust. Soc. Am. 91: p. 3346.

media: A time domain approach, J. Acoust. Soc. Am. 107: (2), 683-688.

*MaineUFR STS Avenue O. Messiaen 72085 Le Mans CEDEX 09 France*

Materials, Chapman and Hall, London.

media, Wave Motion 41 (2): 145-161.

porous media J. Appl. Phys. 70: p. 1975.

power law, J. Acoust. Soc. Am. 96: p. 491.

Wave Mat. Interact. 1: 365-380.

Rev. B, 47: p. 4964.

358.

p. 151-163.

**4. References**

(1) 61-72.

	- [36] Fellah ZEA, Fellah M, Mitri FG, Sebaa N, Depollier C, Lauriks W, Measuring permeability of porous materials at low frequency range via acoustic transmitted waves, Rev. Sci. Instrum. 78: 114902-10.
	- [37] Sadouki M, Fellah M, Fellah ZEA, Ogam E, Sebaa N, Mitri FG, Depollier C, Measuring static thermal permeability and inertial factor of rigid porous materials, J. Acoust. Soc. Am. 130: 2627-2630.
	- [38] Lagarias JC, Reeds JA, Wright MH, Wright PE (1998), Convergence properties of the Nedler-Mead Simplex Method in low dimensions SIAM J. Optim. 9: 112-147.
	- [39] Henry M, Allard JF (1997), Acoustical measurement of the trapping constant of foams with open cells, C. R. Acad. Sci. Paris 325: 331-338.
	- [40] Biot MA (1956), The theory of propagation of elastic waves in fluid-saturated porous solid. I. Low frequency range, J. Acoust. Soc. Am. 28: 168.
	- [41] Biot MA (1956), The theory of propagation of elastic waves in fluid-saturated porous solid. I. Higher frequency range, J. Acoust. Soc. Am. 28: 179.
	- [42] Fellah ZEA, Chapelon JY, Berger S, Lauriks W and Depollier C (2004), Ultrasonic wave propagation in human cancellous bone: Application of Biot theory, J. Acoust. Soc. Am. 116: (1), 61-73.
	- [43] Sebaa N, Fellah ZEA, Fellah M, Ogam E, Wirgin A, Mitri FG, Depollier C, Lauriks W (2006), Ultrasonic characterization of human cancellous bone using the Biot theory: Inverse problem, J. Acoust. Soc. Amer. 120: 1816-1824.
	- [44] Ogam E, Fellah ZEA, Sebaa N, Groby JP (2011), Non-ambiguous recovery of Biot poroelastic parameters of cellular panels using ultrasonic waves. J. Sound and Vibr. 330(6):1074-1090.
	- [45] Ogam E, Fellah ZEA (2011), The direct and inverse problems of an airsaturated poroelastic cylinder submitted to acoustic radiation. AipAdvances, 1(3):032174.
	- [46] Ogam E, Depollier C, Fellah ZEA (2010), The direct and inverse problems of an air-saturated porous cylinder submitted to acoustic radiation. Rev. Sci. Inst. 81(094902):9.
	- [47] Sebaa N, Fellah ZEA, Lauriks W, Depollier C (2006), Application of fractional calculus to ultrasonic wave propagation in human cancellous bone, Signal Processing 10: 2668-2677, DOI : 10.1016/j.sigpro.2006.02.015.
	- [48] Ashman RB, Rho JY (1988), Elastic modulus of trabecular bone material, J. Biomech. 21: 177-181.

© 2013 Hosseini Fouladi et al., licensee InTech. This is an open access chapter distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

© 2013 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution,

© 2013 Hosseini Fouladi et al., licensee InTech. This is a paper distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use,

distribution, and reproduction in any medium, provided the original work is properly cited.

and reproduction in any medium, provided the original work is properly cited.

**Utilizing Malaysian Natural Fibers as** 

Previous researchers [1, 2, 3, 4] successfully produced sound absorption panels using agricultural wastes. These panels play important role both in noise absorption and heat insulation preserving the comfort of indoor living spaces. Yang et al. [4] produced rice straw-wood particle composite boards. They found that the sound absorption coefficient of these boards is higher than other wood-based materials in the 500-8000 Hz frequency range. Reason was the low specific gravity of composite boards having high amount of porosity compared to wood-based materials [4]. Another study by Davern [5] aimed at producing airspace layers and examined the influence of porosity on the acoustic properties of materials. He found that the porosity of the perforated plate and the density of the porous material would significantly affect the acoustic impedance and sound absorption coefficient of the panel, in which case, the frequency band near the resonance frequency achieved high acoustic absorption. In addition to Davern's study, Lee and Chen [6] reported that the acoustic absorption of multi-layer materials is better with a perforated plate backed with airspaces. Other usages of natural fibres are in reduction of sound propagation in automotive interior spaces, or to improve the control of outdoor noise propagation [7, 8]. Recent studies show that researchers are focusing on coir fibre and oil palm fibre in replacing synthetic-based fibres for sound absorption applications simply due to their

Acoustics, as science, deals with sound creation and transmission through materials. Sound and light waves share same vibrational system as these waves creating pressure during propagation due to the nature of fluctuation in the material. However, unlike electromagnetic waves, sound cannot travel through a vacuum. Consequently, sound wave has become important in many applications in materials science, medicine, dentistry,

Mohammad Hosseini Fouladi, Mohamed H. Nassir, Masomeh Ghassem, Marwan Shamel, Sim Yeng Peng, Sin Yi Wen, Pang Zong Xin and Mohd Jailani Mohd Nor

Additional information is available at the end of the chapter

abundance in tropical countries such as Malaysia.

**Sound Absorber** 

http://dx.doi.org/10.5772/53197

**1. Introduction** 
