This chapter provides different models for the acoustic wave propagation in porous materials having a rigid and an elastic frames. The direct problem of reflection and transmission of acoustic waves by a slab of porous material is studied. The inverse problem is solved using experimental reflected and transmitted signals. Both high- and low-frequency domains are studied. Different acoustic methods are proposed for measuring physical parameters describing the acoustic propagation as porosity, tortuosity, viscous and thermal characteristic length, and flow resistivity. Some advantages and perspectives of this method are discussed.
- acoustic porous materials
- viscous and thermal charactertistic lengths
- fractional derivatives
More than 50 years ago, Biot [1, 2] proposed a semi-phenomenological theory which provides a rigorous description of the propagation of acoustic waves in porous media saturated by a compressible viscous fluid. 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. Biot’s theory describes the motion of the solid and the fluid, as well as the coupling between the two phases. The loss of acoustic energy is due mainly to the viscosity of the fluid and the relative fluid-structure movement. The model predicts that the acoustic attenuation, as well as the speed of sound, depends on the frequency and elastic constants of the porous material, as well as porosity, tortuosity, permeability, etc. The theory predicts two compressional waves: a fast wave, where the fluid and solid move in phase, and a slow wave where fluid and solid move out of phase. Johnson et al.  introduced the concept of tortuosity or dynamic permeability which has better described the viscous losses between fluid and structure in both high and low frequencies.
Air-saturated porous materials such as plastic foams or fibrous materials are widely used in passive control and noise reduction. These materials have interesting acoustic properties for sound absorption, and their use is quite common in the building trade and automotive and aeronautical fields. The determination of the physical parameters of the medium from reflected and transmitted experimental data is a classical inverse scattering problem.
Pulse propagation in porous media is usually modeled by synthesizing the signal via a Fourier transform of the continuous wave results. On the other hand, experimental measurements are usually carried out using pulses of finite bandwidth. Therefore, direct modeling in the time domain is highly desirable [4–10]. The temporal and frequency approaches are complementary for studying the propagation of acoustic signals. For transient signals, the temporal approach is the most appropriate because it is closer to the experimental reality and the finite duration of the signal. However, for monochromatic harmonic signals, the frequency approach is the most suitable .
Fractional calculus has been used in the past by many authors as an empirical method to describe the viscoelastic properties of materials (e.g., see Caputo  and Bagley and Torvik ). The fact that acoustic attenuation, stiffness, and damping in porous materials are proportional to the fractional powers of frequency [4, 5, 7, 9, 10] suggests that fractional-order time derivatives could describe the propagation of acoustic waves in these materials.
In this chapter, acoustic wave propagation in porous media is studied in the high- and the low-frequency range. The direct and inverse scattering problems are solved for the mechanical characterization of the medium. The general Biot model applied to porous materials having elastic structure is treated, and also the equivalent fluid model, used for air-saturated porous materials ( Figures 1 and 2 ).
2. Porous materials with elastic frame
In porous media, the equations of motion of the frame and fluid are given by the Euler equations applied to the Lagrangian density. Here, and are the displacements of the solid and fluid phases. The equations of motion are given by [1, 2]
where , , and are the generalized elastic constants, is the porosity, is the bulk modulus of the pore fluid, is the bulk modulus of the elastic solid, and is the bulk modulus of the porous skeletal frame. is the shear modulus of the composite as well as that of the skeletal frame. The equations which explicitly relate , , and to , , , , and are given by
is the “mass coefficients” which are related to the densities of solid ( ) and fluid ( ) phases by
The Young modulus and the Poisson ratio of the solid and and of the skeletal frame and depend on the generalized elastic constant , , and via the relations:
The mass coupling parameter between the fluid and solid phases is always negative:
where is the tortuosity of the medium. The damping of the acoustic wave in porous material is essentially due to the viscous exchanges between the fluid and the structure. To express the viscous losses, the dynamic tortuosity is introduced  given by
where , is the angular frequency, is the fluid resistivity, is the viscous permeability, and is the viscous characteristic length given by Johnson et al. . The ratio of the sizes of the pores to the viscous skin depth thickness gives an estimation of the parts of the fluid affected by the viscous exchanges. In this domain of the fluid, the velocity distribution is perturbed by the frictional forces at the interface between the viscous fluid and the motionless structure. At high frequencies, the viscous skin thickness is very thin near the radius of the pore . The viscous exchanges are concentrated in a small volume near the surface of the frame . The expression of the dynamic tortuosity is given by 
The range of frequencies such that viscous skin thickness is much larger than the radius of the pores
is called the low-frequency range. For these frequencies, the viscous forces are important everywhere in the fluid. When , the expression of the dynamic tortuosity becomes
where is the microscopic velocity. The angle brackets represent the average of the random variable over the sample of material. In the time domain, and in the high-frequency domain, the dynamic tortuosity (Eq. 7) acts as the operator, and its expression is given by 
is the Dirac function. In this model the time convolution of with a function is interpreted as a semi-derivative operator according to the definition of the fractional derivative of order given by Samko et al. :
where and is the gamma function. A fractional derivative acts as a convolution integral operator and no longer represents the local variations of the function. The properties of fractional derivatives and fractional calculus are given by Samko et al. .
The introduction of the tortuosity operator (Eq. 11) in Biot’s Eqs. (1) and (2) to describe the inertial and viscous interactions between fluid and structure will express the propagation equations in the time domain. When is used instead of in Eqs. (1) and (2), the equations of motion (1) and (2) will be written as 
In these equations, the temporal operators , , and represent the mass coupling operators between the fluid and solid phases and are given by
where is given by Eq. (11).
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, and , are defined for compressional waves giving
where , is the Laplacian, and represents the fractional derivative following the definition given by Eq. (12).
Two distinct longitudinal modes called fast and slow waves are obtained by the resolution of the eigenvalue problem of the matrix of Biot (Eq. (14)). On a basis of fast and slow waves and , one can have
where and are the “eigenvalue operators” of the Biot matrix (Eq. (14)). Their expressions are given by
Their corresponding eigenvectors are
Coefficients , , and are given by
The fast and slow waves and are obeying to the following propagation equations along the axis:
where the coefficients , , and are constants, respectively, given by
where Eq. (18) is a fractional propagation equations  in time domain of the fast and slow waves, respectively. These equations describe the attenuation and the spreading of the temporal signal propagating inside the porous material. These fractional propagation equations have been solved and well-studied in the case of rigid porous materials using the equivalent fluid model.
3. Porous materials with rigid frame
In the acoustics of porous media, two situations can be distinguished: elastic and rigid frame materials. In the first case, the Biot [1, 2] theory is best suited. In the second case, the acoustic wave cannot vibrate the structure. The equivalent fluid model is then used, in which the acoustic wave propagates inside the saturating fluid [8, 11]. The equations for the acoustics in the equivalent fluid model are given by
In these relations, is the acoustic pressure. The first equation is the Euler equation, and the second one is a constitutive equation obtained from the equation of mass conservation associated with the behavior (or adiabatic) equation. These equations can be obtained from the Biot Eqs. (1, 2) by canceling the solid displacement. Assuming that the porous medium studied is homogeneous and has a linear elasticity, we obtain easily the following wave equation (propagation along the axis) for the acoustic pressure in a lossless porous material:
In Eq. (20), the viscous and thermal losses that contribute to the sound damping in acoustic materials are not described. The thermal exchanges are generally negligible near viscous effects in the porous materials obeying to the Biot theory, this is not the case for air-saturated porous materials using the equivalent fluid model. To take into account the fluid-structure exchanges, the density and compressibility of the fluid are “renormalized” by the dynamic tortuosity and the dynamic compressibility , via the relations and , giving the following wave equation in frequency domain (Helmholtz equation) for a lossy porous material:
The thermal exchanges to the fluid compressions-dilatations are produced by the wave motion. The parts of the fluid affected by the thermal exchanges can be estimated by the ratio of a microscopic characteristic length of thermal skin depth thickness ( is the fluid viscosity; is the Prandtl number).
The expression of the dynamic compressibility is given by
where is the adiabatic constant, the magnitude introduced by Lafarge  called thermal permeability by analogy to the viscous permeability, and is the thermal characteristic length. The low-frequency approximation of  is given by
where , which has the same size (area) that of Darcy’s permeability of , is a parameter analogous to the parameter but is adapted to the thermal problem.
In a high-frequency limit, Allard and Champoux  showed the following behavior of :
In the time domain (using the convention ), we obtain the following fractional propagation equation:
In this equation, the term is interpreted as a semi-derivative operator following the definition of the fractional derivative of order , given by Samko and coll. . The solution of the wave Eq. (26) with suitable initial and boundary conditions is by using the Laplace transform. is the medium’s Green function  given by
where has the following form:
Let us consider a homogeneous porous material which occupies the region ; the expressions of the reflection and transmission coefficients in the frequency domain are given by
These expressions are simplified by taking into account the reflections at the interfaces and ; the expressions of the reflection and transmission operators are given in time domain by
where is the Dirac function and is the Green function of the medium given by Eq. (27). In the next sections, we will use the reflected and transmitted waves for solving the inverse problem in order to characterize the porous materials.
3.1. Ultrasonic measurement of porosity, tortuosity, and viscous and thermal characteristic lengths via transmitted waves
The experimental setup consists of two transducers broadband Ultran NCT202 with a central frequency of 190 kHz in air and a bandwidth of 6 dB extending from 150 to 230 kHz . A pulser/receiver 5058PR Panametrics sends pulses of 400 V. The high-frequency noise is avoided by filtering the received signals above 1 MHz. Electronic interference is eliminated by 1000 acquisition averages. The experimental setup is shown in Figure 3 . The inverse problem is to find the parameters , , , and which minimize numerically the discrepancy function wherein is the discrete set of values of the experimental transmitted signal and is the discrete set of values of the simulated transmitted signal predicted from Eq. (33). The least squares method is used for solving the inverse problem using the simplex search method (Nelder-Mead)  which does not require numerical or analytic gradients.
Consider a sample of plastic foam M1, of thicknesses . Sample M1 was characterized using classic methods [21–31] and gave the following physical parameters , , , and . Figure 4 shows the experimental incident signal (dashed line) generated by the transducer and the experimental transmitted signal (solid line). After solving the inverse problem simultaneously for the porosity , tortuosity , and viscous and thermal characteristic lengths and , we find the following optimized values: , , , and . The values of the inverted parameters are close to those obtained by conventional methods [21–31]. We present in Figures 5 and 6 the variation of the minimization function with the porosity, tortuosity, viscous characteristic length, and the ratio between and . In Figure 7 , we show a comparison between an experimental transmitted signal and simulated transmitted signal for the optimized values of , , , and . The difference between the two curves is small, which leads us to conclude that the optimized values of the physical parameters are correct.
3.2. Measuring flow resistivity of porous material via acoustic reflected waves at low-frequency domain
In the low-frequency domain, the viscous forces are important everywhere in all the fluid saturating the porous material. The thermal exchanges between fluid and structure are favored by the slowness of the cycle of expansion and compression in the material. The temperature of the frame is practically unchanged by the passage of the sound wave because of the high value of its specific heat: the frame acts as a thermostat; the isothermal compressibility is directly applicable. In this domain, the viscous skin thickness is much larger than the radius of the pores
For a wave traveling along the direction , the generalized forms of the basic Eqs. (19) in the time domain are now
where the Euler equation is reduced to Darcy’s law which defines the static flow resistivity . The wave equation in time domain is given by
The fields which are varying in time, the pressure, the acoustic velocity, etc. follow a diffusion equation with the diffusion constant:
The diffusion constant is connected to Darcy’s constant (called also the viscous permeability) by the relation
where is the fluid viscosity.
The expression of the reflection coefficient in Laplace domain (put for obtaining the frequency domain of ), is given by 
The development of these expressions in exponential series leads to the reflection coefficient:
The multiple reflections in the material are taken into account in these expressions. As the attenuation is high in the porous materials, the multiple reflection effects are negligible. Let us consider the reflections at the interfaces and :
The reflection scattering operator is calculated by taking the inverse Laplace transform of the reflection coefficient.
We infer  that
where erf is the error function. By putting
Using the relation
which with the variable change , yields
The reflection scattering operator is then given by
3.2.1. Acoustic parameter sensitivity
Consider a sample of porous material having a physical parameters that correspond to quite common acoustic materials, as follows: thickness , porosity , flow resistivity , and radius of the pore . Let us study the sensitivity of the main parameters using numerical simulations of waves reflected by a porous material. Fifty percent variation is applied to the physical parameters (flow resistivity and porosity ).
To obtain the simulated reflected waves, we use the incident signal given in Figure 8 (dashed line). The result (reflected wave) is the wave given in the same figure ( Figure 8 ) in solid line. The spectra of the two waves (incident and reflected) are given in Figure 9 . From Figure 8 , we can see that there is just an attenuation of the reflected wave without dispersion, since the two waves have the same spectral bandwidth ( Figure 9 ). Figure 8 shows the results obtained after reducing flow resistivity by of its initial value. The wave in dashed line corresponds to the simulated reflected signal for and the second one (solid line) to . The values of the porosity and thickness have been kept constant. When the flow resistivity is reduced, the amplitude of reflected wave decreases by 30% of its initial value. Physically, by reducing the flow resistivity, the medium is less resistive, since the viscous effects become less important in the porous material, and thus the amplitude of the reflected wave decreases. No change is observed in the reflected wave when reducing the porosity by of its initial value. We can conclude that the porosity has no significant sensitivity in reflected mode.
For the propagation of transient signals at low frequency, a guide (pipe) , having a diameter of 5 cm and of length 50 m, is chosen. The pipe can be rolled without perturbations on experimental signals (the cutoff frequency of the tube ). The same microphone (Brüel & Kjær, 4190) is used for measuring the incident and reflected signals. Burst is provided by synthesized function generator Stanford Research Systems model DS345-30 MHz. A sound source driver unit “Brand” constituted by loudspeaker Realistic 40-9000 is used. The incident signal is measured by putting a total reflector in the same position than the porous sample. The experimental setup is shown in Figure 10 . Consider a cylindrical sample of plastic foam M1 of flow resistivity value . This value is obtained using the method of Bies and Hansen . The sample M1 has a diameter of 5 cm and a thickness of 3 cm. Figure 11 shows the experimental incident wave (solid line) generated by the loudspeaker in the frequency bandwidth (35–75) Hz, and the experimental reflected signal (dashed line), with their spectra. There is no dispersion, since the two signals have practically the same bandwidth. The minimization of the function gives the solution if the inverse problem:
where and represent the discrete set of values of the experimental reflected signal and of the simulated reflected signal, respectively. The optimized value of is obtained by solving the inverse problem. The variation of the minimization function with the flow resistivity is given in Figure 12 . A comparison between experiment and theory is given in Figure 13 . The difference between theory and experiment is slight, which leads us to conclude that the optimized value of the flow resistivity is good.
This alternative acoustic method has the advantage of being simple and effective since it requires the use of only one microphone and therefore no calibration problem. In addition, this approach is different from conventional methods (Bies and Hansen ) that involve the use of fluid flow measurement techniques and pressure differences. The mathematical analysis of the reflected wave at low frequency is quite simple, because this wave is not propagative in the medium but simply diffusive (having the same frequency band with the incident signal). The wave reflected by the resistive materials has the advantage of being easily detectable experimentally compared to the transmitted wave.
Acoustic propagation in porous media involves a large number of physical parameters when the structure is elastic. This number is reduced when the structure is rigid, because the mechanical part does not intervene and thus remains only the acoustic part. The study of high and low frequencies separately solves the inverse problem and characterizes the porous materials in the domain of influence of the physical parameters. The proposed methods are simple and effective and allow an acoustic characterization of porous materials using transmitted or reflected experimental waves.
Biot MA. The theory of propagation of elastic waves in fluid-saturated porous solid. I. Low frequency range. The Journal of the Acoustical Society of America. 1956; 28:168
Biot MA. The theory of propagation of elastic waves in fluid-saturated porous solid. I. Higher frequency range. The Journal of the Acoustical Society of America. 1956; 28:179
Johnson DL, Koplik J, Dashen R. Theory of dynamic permeability and tortuosity in fluid-saturated porous media. Journal of Fluid Mechanics. 1987; 176:379-402
Szabo TL. Time domain wave equations for lossy media obeying a frequency power law. The Journal of the Acoustical Society of America. 1994; 96:491
Szabo TL. Causal theories and data for acoustic attenuation obeying a frequency power law. Journal of the Acoustical Society of America. 1995; 97:14
Norton V, Novarini JC. Including dispersion and attenuation directly in time domain for wave propagation in isotropic media. The Journal of the Acoustical Society of America. 2003; 113:3024
Chen W, Holm S. Modified Szabo’s wave equation models for lossy media obeying frequency power law. The Journal of the Acoustical Society of America. 2003; 113:3024
Fellah ZEA, Depollier C. Transient acoustic wave propagation in rigid porous media: A time domain approach. The Journal of the Acoustical Society of America. 2000; 107(2):683-688
Fellah ZEA, Fellah M, Lauriks W, Depollier C, Angel Y, Chapelon JY. Solution in time domain of ultrasonic propagation equation in porous material. Wave Motion. 2003; 38:151-163
Fellah ZEA, Depollier C, Fellah M. Direct and inverse scattering problem in porous material having a rigid frame by fractional calculus based method. Journal of Sound and Vibration. 2001; 244(2):359-366
Allard JF. Propagation of Sound in Porous Media : Modeling Sound Absorbing Materials. London: Chapman and Hall; 1993
Caputo M. Vibration of an infinite plate with a frequency dependent Q. The Journal of the Acoustical Society of America. 1976; 60:634-639
Bagley RL, Torvik PJ. On the fractional calculus model of viscoelastic behavior. Journal of Rheology. 1983; 30:133-155
Lafarge D, Lemarnier P, Allard JF, Tarnow V. Dynamic compressibility of air in porous structures at audible frequencies. Journal of the Acoustical Society of America. 1996; 102:4
Norris AN. On the viscodynamic operator in Biot’s equations of poroelasticity. Journal of Wave-Material Interaction. 1986; 1:365-380
Samko SG, Kilbas AA, Marichev OI. Fractional Integrals and Derivatives Theory and Applications. Amsterdam: Gordon and Breach Publishers; 1993
Fellah M, Fellah ZEA, Mitr FG, Ogam E, Depollier C. Transient ultrasound propagation in porous media using biot theory and fractional calculus: Application to human cancellous bone. The Journal of the Acoustical Society of America. 2013; 133(4):683-688
Allard JF, Champoux Y. New empirical equations for sound propagation in rigid frame fibrous materials. Journal of the Acoustical Society of America. 1992; 91:3346-3353
Fellah ZEA, Sadouki M, Fellah M, Mitri FG, Ogam E, Depollier C. Simultaneous determination of porosity, tortuosity, viscous and thermal characteristic lengths of rigid porous materials. Journal of Applied Physics. 2013; 114:204902-204905
Lagarias JC, Reeds JA, Wright MH, Wright PE. Convergence properties of the Nelder–mead simplex method in low dimensions. SIAM Journal on Optimization. 1998; 9:112Ű147
Leclaire P, Kelders L, Lauriks W, Glorieux C, Thoen J. Determination of the viscous characteristic length in air-filled porous materials by ultrasonic attenuation measurements. The Journal of the Acoustical Society of America. 1996; 99:1944
Ayrault C, Moussatov A, Castagnède B, Lafarge D. Ultrasonic characterization of plastic foams via measurements with static pressure variations. Applied Physics Letters. 1999; 74:3224
Moussatov A, Ayrault C, Castagnède B. Porous material characterization ultrasonic method for estimation of tortuosity and characteristic length using a barometric chamber. Ultrasonics. 2001; 39:195
Leclaire P, Kelders L, Lauriks W, Brown NR, Melon M, Castagnède B. Determination of viscous and thermal characteristic lengths of plastic foams by ultrasonic measurements in helium and air. Journal of Applied Physics. 1996; 80:2009
Fellah ZEA, Depollier C, Berger S, Lauriks W, Trompette P, Chapelon JY. Determination of transport parameters in air saturated porous materials via ultrasonic reflected waves. The Journal of the Acoustical Society of America. 2003; 113(5):2561-2569
Fellah ZEA, Berger S, Lauriks W, Depollier C, Chapelon JY. Inverse problem in air-saturated porous media via reflected waves. Review of Scientific Instruments. 2003; 74(5):2871
Fellah ZEA, Berger S, Lauriks W, Depollier C, Aristégui C, Chapelon JY. Measuring the porosity and the tortuosity of porous materials via reflected waves at oblique incidence. The Journal of the Acoustical Society of America. 2003; 113(5):2424
Fellah ZEA, Berger S, Lauriks W, Depollier C, Trompette P, Chapelon JY. Ultrasonic measuring of the porosity and tortuosity of air- saturated random packings of beads. Journal of Applied Physics. 2003; 93:9352
Fellah ZEA, Mitri FG, Depollier D, Berger S, Lauriks W, Chapelon JY. Characterization of porous materials having a rigid frame via reflected waves. Journal of Applied Physics. 2003; 94:7914
Fellah ZEA, Berger S, Lauriks W, Depollier C, Fellah M. Measurement of the porosity of porous materials having a rigid frame via reflected waves : A time domain analysis with fractional derivatives. Journal of Applied Physics. 2003; 93:296
Fellah ZEA, Mitri FG, Fellah M, Ogam E, Depollier C. Ultrasonic characterization of porous absorbing materials: Inverse problem. Journal of Sound and Vibration. 2007; 302:746-759
Sebaa N, Fellah ZEA, Fellah M, Lauriks W, Depollier C. Measuring flow resistivity of porous material via acoustic reflected waves. Journal of Applied Physics. 2005; 98:084901
Bies DA, Hansen CH. Flow resistance information for acoustical design. Applied Acoustics. 1980; 13:357-391