This chapter provides a temporal method for measuring the porosity and the tortuosity of air-saturated porous materials using experimental reflected waves. The direct problem of reflection and transmission of acoustic waves by a slab of porous material is studied. The equivalent fluid model has considered in which the acoustic wave propagates only in the pore-space. Since the acoustic damping in air-saturated porous materials is important, only the reflected waves by the first interface are taken into account, and the multiple reflections are neglected. The study of the sensitivity analysis shows that porosity is much more sensitive than tortuosity to reflection, especially when the incident angle is less than its critical value, at which the reflection coefficient vanishes. The inverse problem is solved using experimental data at a different incidence angle in reflection. Some advantages and perspectives of this method are discussed.
- porous material
Porosity [1, 2, 3] is one of the most important parameters for describing the acoustic propagation in porous materials. This parameter intervenes in the propagation phenomena at all frequencies. Porosity is the relative fraction, by volume, of the air contained within the material. Air-saturated porous materials [1, 2, 4] as plastic foams, fibrous or granular materials are of great interest for a wide range of industrial applications. Figure 1 gives an example of air-saturated porous material commonly used in the sound absorption (passive control). These materials are frequently used in the automotive and aeronautics industries and in the building trade. Beranek  has developed an apparatus for measuring the porosity of air saturated porous materials. This device was based on the equation of state for ideal gases at constant temperature i.e., Boyle’s law. Porosity can be determined by measuring the change in air pressure occurring with a known change in volume of the chamber containing the sample. In this apparatus, both pressure change and volume change are monitored using a U-shaped fluid-filled manometer. Leonard  has given an alternate dynamic method for measuring porosity. Other techniques using water as the pore-filling fluid, rather than air, are common in geophysical studies [7, 8]. Mercury has been used as the pore-filling fluid in other applications . However, for many materials, the introduction of liquids into the material is not appropriate. Recently, a similar device to that of Beranek , involving the use of an electronic pressure transducer, was introduced by Champoux et al. . This device can be used to measure very slight changes in pressure accurately, and the output can be recorded by a computer.
Generally, the most methods used for measuring the porosity cited previously do not use acoustic waves. Here, we present an ultrasonic method for measuring porosity using ultrasonic reflected waves by the porous material. The direct and inverse problem is solved in the time domain using experimental reflected data. The inverse problem is solved directly in time domain using the waveforms. The attractive feature of a time domain-based approach [11, 12, 13, 14, 15, 16] 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.
In the acoustics of porous materials, a distinction can be made between two situations depending on whether the frame is moving or not. In the general case (when the frame is moving), as for cancellous bone and rock and all porous media saturated by a liquid, the dynamics of the waves due to the coupling between the solid frame and the fluid are clearly described by the Biot theory [17, 18]. In air-saturated porous media, as plastic foams or fibrous materials that are used in sound absorption, the structure is generally motionless and the waves propagate only in the fluid. This case is described by the model of equivalent fluid, which is a particular case in the Biot model, in which the interactions between the fluid and the structure are taken into account in two frequency response factors: the dynamic tortuosity of the medium given by Johnson et al.  and the dynamic compressibility of the air included in the porous material given by Allard and Champoux . In the frequency domain, these factors multiply the density of the fluid and its compressibility, respectively, and represent the deviation from the behavior of the fluid in free space as the frequency increases. In the time domain, they act as operators, and in the high-frequency approximation, their expressions are given by 
In these equations, is the Dirac function, is the Prandtl number, and are fluid viscosity and fluid density, respectively, and is the adiabatic constant. The most important physical parameters of the model are the medium’s tortuosity initially introduced by Zwikker and Kosten  and viscous and the thermal characteristic lengths and introduced by Johnson et al.  and Allard and Champoux . 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 in Samko et al. 
where is the Gamma function.
In this framework, the basic equations of our model can be expressed as follows
where denotes the time convolution operation, is the acoustic pressure, is the particle velocity and is the bulk modulus of the air. The first equation is the Euler equation, the second is a constitutive equation obtained from the equation of mass conservation associated with the behavior (or adiabatic) equation.
In the plane (), the constitutive equation (4) can be written as
where and are the components of particle velocity along axes and . In these equations, the integrals express the fractional derivatives, defined mathematically by convolutions and describe the dispersive nature of the porous material. They take into account the memory effects due to the fluid-structure interactions and to the viscous and thermal losses in the medium.
The problem geometry is shown in Figure 2. A homogeneous porous material occupies the region . This medium is assumed to be isotropic and to have a rigid frame. A short sound pulse impinges at oblique incidence on the medium from the left, giving rise to an acoustic pressure field and an acoustic velocity field within the material, which satisfies the system of Eq. (5) that can be written as:
In the region , the incident pressure wave is given by
where is the velocity of the free fluid (); .
In the region , the pressure wave is given by
To simplify the system of Eq. (5), we can then use the following property
where has the following form:
and , and .
If the incident sound wave is launched in region , then the expression of the pressure field in the region on the left of the material is the sum of the incident and reflected fields
Here, is the field in the region , and denote the incident and reflected fields, respectively.
The incident and reflected fields are related by the scattering operator (i.e., the reflection operator) for the material. This is an integral operator represented by
In Eq. (22), the function is the reflection kernel for incidence from the left. Note that the lower limit of integration in Eq. (22) is set to 0, which is equivalent to assuming that the incident wavefront impinges on the material at
Expression of the reflection-scattering operator taking into account the n-multiple reflections in the material is given by
Generally , in air-saturated porous materials, acoustic damping is very important, and the multiple reflections are thus negligible inside the material. So, by taking into account only the first reflections at interfaces and , the expression of the reflection operator will be given by
where is the instantaneous response of the porous material corresponding to the reflection contribution at the first interface (). is equivalent to reflection at the interface , which is the bulk contribution to reflection. The part of the wave corresponding to is not subjected to dispersion but simply multiplied by the factor .
The reflection coefficient at the first interface vanishes for a critical angle
Figure 3 shows the variation of the reflection coefficient at the first interface with the incident angle , for a fixed porosity value , and for different values of tortuosity (solid line), (star), (dashed line), (dash-dot line) and (circle).
Figure 4 shows the variation of with the incident angle, for a fixed tortuosity of , and for different values of porosity (solid line), (star), (dashed line), (dash-dot line) and (circle).
When the incident angle is , the reflection coefficient decreases slowly with the incident angle, and when it is , the reflection coefficient increases quickly with the angle. It can be seen also from Figures 2 and 3 that the sensitivity of porosity variation is more important than the sensitivity of tortuosity on the reflection coefficient at the first interface.
3. Inverse problem
The inverse problem is solved using transmitted waves and an estimation of the tortuosity, viscous and thermal characteristic length is given. However, the porosity cannot be inverted in transmission since its sensitivity is low. In this work, we determine porosity and tortuosity by solving the inverse problem for waves reflected by the first interface and by taking into account experimental data concerning all measured incident angles (Figure 4). The propagation of acoustic waves in a slab of porous material in the high-frequency asymptotic domain is characterized by four parameters: porosity , tortuosity , viscous characteristic length , and thermal characteristic length , the values of which are crucial to the behavior of sound waves in such materials. It is of some importance to work out new experimental methods and efficient tools for their estimation. The basic inverse problem associated with the slab may be stated as follows: from the measurement of the signals transmitted and/or reflected outside the slab, find the values of the medium’s parameters. The inverse problem is to find values of parameters , , which minimize the function:
where represents the discrete set of values of the experimental reflected signal for different incident angles , is the reflection coefficient at the first interface, and is the experimental incident signal. The term represents the predicted reflected signal. The inverse problem is solved numerically by the least-square method.
Let us consider three samples of plastic foam M1, M2, and M3. Their tortuosity and porosity were measured using classical methods [22, 23, 24, 25, 26, 27, 28, 29, 30, 31] (Table 1). We solved the inverse problem for these samples via waves reflected at the first interface and for different incident angles.
Experiments were performed in air with two broadband Ultran NCT202 transducers with a 190 kHz central frequency in air and a 6 dB bandwidth extending from 150 to 230 kHz. Pulses of 400 V were provided by a 5052PR panametrics pulser/receiver. An optical goniometer was used to position the transducers. The received signals were amplified to 90 dB and filtered above 1 MHz to avoid high-frequency noise. Electronic interference was removed by 1000 acquisition averages. The experimental setup is shown in Figure 5.
Figures 6 and 7 show a variation of the cost function, , with tortuosity and porosity, for samples M1, M2, and M3, respectively. The reconstructed values of porosity and tortuosity corresponding to the positions of the minima of these cost functions are given in Table 2.
A comparison is given in Figures 8–10, between simulated reflection coefficient at the first interface using reconstructed values of porosity and tortuosity (solid line) and experimental data of the reflection coefficient at the first interface (circle) for plastic foams M1, M2, and M3, respectively.
The correspondence between experiment and theory is good, which leads us to conclude that this method based on the solution of the inverse problem is appropriate for estimating the porosity and tortuosity of porous materials with a rigid frame. This method is very interesting comparing the classical one using non-acoustic waves.
In this chapter, an inverse determination of porosity and tortuosity is given using experimental reflected acoustic waves at different incidence angles. The inverse problem is solved numerically by the least-square method. The obtained values of porosity and tortuosity are close to those given using classical (non-acoustical) methods. Generally, the tortuosity is inverted easily using transmitted waves, and this is not the case for porosity because of its weak sensitivity in transmitted mode. This method is very interesting from its simplicity and constitutes an alternative to the usual method involving the use of a porosimeter introduced by Beranek  and improved by Champoux et al.  or the other ultrasonic methods based on transmitted mode.
This method is valuable for porous materials having a rigid frame as air-saturated media, in which the equivalent fluid model is used. However, for liquid-saturated materials as rock or cancellous bone, this method cannot be used and should be adapted using the general Biot model [17, 18]. We hope, in the future, to extend this method to porous media with an elastic frame saturated with viscous fluid, in order to estimate other parameters that play an important role in acoustic propagation.
Allard JF. Propagation of Sound in Porous Media: Modeling Sound Absorbing Materials. London: Chapman and Hall; 1993
Zwikker C, Kosten CW. Sound Absorbing Materials. New York: Elsevier; 1949
Beranek LL. Acoustic impedance of porous materials. The Journal of the Acoustical Society of America. 1942; 13:248-260
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
Allard JF, Champoux Y. New empirical equations for sound propagation in rigid frame fibrous materials. Acoustical Society of America. 1992; 91:3346
Leonard RW. Simplified porosity measurements. Journal of the Acoustical Society of America. 1948; 20:39-34
Guyon E, Oger L, Plona TJ. Transport properties in sintered porous media composed of two particles sizes. Journal of Physics D: Applied Physics. 1987; 20:1637-1644
Johnson DL, Plona TJ, Scala C, Psierb F, Kojima H. Tortuosity and acoustic slow waves. Physical Review Letters. 1982; 49:1840-1844
Van Brakel J, Modry S, Svata M. Mercury porosimetry: State of the art. Powder Technology. 1981; 29:1-12
Champoux Y, Stinson MR, Daigle GA. Air-based system for the measurement of porosity. The Journal of the Acoustical Society of America. 1991; 89:910-916
Fellah ZEA, Depollier C. Transient acoustic wave propagation in rigid porous media: A time domain approach. Journal of the Acoustical Society of America. 2000; 107(2):683-688
Fellah ZEA, Depollier C, Fellah M, Lauriks W, Chpaleon JY. Influence of dynamic tortuosity and compressibility on the propagation of transient waves in porous media. Wave Motion. 2005; 41(2):145-161
Fellah ZEA, Fellah M, Lauriks W, Depollier C. Direct inverse scattering of transient acoustic waves by a slab of rigid porous material. Journal of the Acoustical Society of America. 2003; 113(1):61-72
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
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
Samko SG, Kilbas AA, Marichev OI. Fractional Integrals and Derivatives Theory and Applications. Amsterdam: Gordon and Breach Publishers; 1993
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, Mitri FG, Depollier C, Berger S, Lauriks W, Chapelon JY. Characterization of porous materials with a rigid frame via reflected waves. Journal of Applied Physics. 2003; 94:7914-7922
Fellah ZEA, Ogam E, Wirgin A, Fellah M, Depollier C, Lauriks W. Ultrasonic characterization of porous materials: Inverse problem. Journal of Sound and Vibration. 2007; 302:746-759
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
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. Detemination 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