Experimental measurement of the porosity and the viscous tortuosity of rigid porous material in low frequency

Abstract

In this paper, an acoustic method is presented for measuring the porosity and the viscous tortuosity of air-saturated-porous materials at low frequencies. The proposed method is based on a temporal model of the direct and inverse problem for the reflection of low-frequency waves by homogeneous isotropic slab of porous material having a rigid frame. Reflected coefficient for a slab of porous material is derived from the responses of the medium to an incident acoustic pulse where a simple relation between flow resistivity, porosity, viscous tortuosity and the reflected waves is obtained. A numerical method and efficient tool for the estimation of the porosity and the viscous tortuosity are presented and discussed.

Keywords

Porous materials rigid frame equivalent fluid model low frequency

Introduction

Good acoustic performance is a desirable attribute in almost all types of buildings and is particularly important for residential buildings, places of worship, schools and hospitals. In this context, several studies have been carried out^1–6 to predict the acoustic behavior of the materials used in these buildings. Porous materials such as plastic foams, fibrous, granular materials, rocks and ceramics are among those materials frequently used in our daily lives and in the buildings’ construction to reduce noise and vibration pollution. The effectiveness of the porous acoustic materials in sound insulation is mainly based on their intrinsic properties. The porosity ϕ is one of the most known parameters which appear in the theories of sound propagation in porous media. This parameter intervenes in the description of the viscous coupling between the fluid and the structure. The porosity ϕ refers to the ratio of the fluid volume occupied by the fluid phase to the total volume of porous material. It is also important as one of the several parameters required by acoustical theory to characterize the porous materials. For acoustical materials, its range of values is generally between 0.70 and 0.99. The porosity can be measured directly and there are several methods to do so.^7–11 For example, porosity can be determined by measuring the change in air pressure that occurs with a known volume change in the chamber containing the sample. The practical implementation of the direct methods for measuring the porosity could be complex and expensive. There are other methods that are faster and less expensive to measure the porosity, the so-called indirect methods^12–16 that use sound waves transmitted or reflected by the porous material. The second known parameter which appears in the theories of sound propagation in porous media in low frequency is the viscous tortuosity α₀, namely also inertial factor, introduced by Norris ¹⁷ and Lafarge¹⁸; the viscous tortuosity α₀ has been used recently in the literature¹⁹ to determine the acoustic absorption of multi-periodic composites; it describes the fluid structure interaction in the low frequency, and this parameter has not been yet measured experimentally. Therefore, it is of some importance to devise experimental methods to estimate this parameter.

In this paper, an acoustical method is proposed for measuring the porosity and the viscous tortuosity using reflected waves by air-saturated porous materials at low frequencies (100–400) Hz. The advantage of this method is the use of a simplified expression of the reflection coefficient which depends on the flow resistivity, the porosity and the viscous tortuosity. The results of the inversion method proposed are reasonable and consistent with those found in the literatures.

Acoustical model

The porous material is a bi-phase medium consisting of a solid part saturated with a fluid. Due to the coupling between the solid skeleton and the fluid, the propagation of acoustic waves within the material is well described by the Biot theory.^20,21 Result of this theory – three waves propagate within the material: two compression waves, one related to the fluid phase and the other to the solid phase, and one shear wave in the solid phase. The porous medium is considered rigid^22–24 when the coupling between the solid and fluid phase is so low that little energy is transferred or when the density or stiffness of the solid phase is high, which requires a significant amount of energy to cause a displacement under acoustic excitation; in this case, the wave can be considered to propagate only in the fluid, and the porous medium is considered as a fluid equivalent^22–23 with effective density $\tilde{ρ}$ and bulk modulus $\tilde{K}$ . These effective accounts for viscous and thermal losses reduce the compression wave. In the frequency domain, these loses are related by the two response factors: The dynamic tortuosity of the medium α(ω) given by Johnson et al.²⁵ and the dynamic compressibility of the fluid included in the porous material β(ω) given by Allard and Atalla²⁶ by the relations $\tilde{ρ} (ω) = α (ω) ρ$ and $\tilde{K} (ω) = K / β (ω)$ where ρ and $K$ are the density and bulk modulus of the fluid in free space and ω is the pulsation frequency. The response factors $α (ω)$ and $β (ω)$ represent the deviation from the behavior of the fluid in free space as the frequency changes.

In low frequency range, the expressions of the dynamic tortuosity and the dynamic compressibility of the medium are related to different geometrical parameters of the porous medium which are the static viscous permeability $k_{0}$ , the porosity ϕ, the viscous tortuosity given by Norris¹⁷ $α_{0}$ , the static thermal permeability $k_{0}^{'}$ ^18,27 of the porous material ( $k_{0}^{'} \geq k_{0}$ ) and new parameter $α_{1}$ by the relations^26,27

α (ω) = - \frac{η}{k_{0}} \frac{ϕ}{ρ} \frac{1}{j ω} + α_{0} + α_{1} j ω + O ({(j ω)}^{2})

(1)

β (ω) = γ + \frac{k_{0}^{′}}{η} \frac{(γ - 1) P_{r} ρ}{ϕ} j ω + O ({(j ω)}^{2})

(2)

where

j = \sqrt{- 1}

and

α_{1} = \frac{2 α_{\infty}^{4} k_{0}^{3} ρ}{η Λ^{4} ϕ^{3} p^{3}}

α_{\infty}

is the tortuosity,²⁶ Λ is the viscous characteristic length,²⁵ p is a geometrical parameter introduced by Pride et al.²⁸ and revised by Lafarge¹⁸ and Lafarge et al.²⁷

η is the fluid viscosity, ρ is the saturating fluid density, γ is the adiabatic constant. $P_{r}$ is the Prandtl number. The static viscous permeability $k_{0}$ is related to the specific flow resistivity σ by the relation $k_{0} = η / σ$ . It should be noted that the third term $α_{1}$ written in equation (1) does not result from the exact expansion of $α (ω)$ and $β (ω)$ , this would involve a new geometrical parameter that has not been calculated and measured. This term results from using Pride et al.’s approximate model.²⁸

The objective of this work is proposed to estimate the porosity and the viscous tortuosity of rigid porous medium using the reflected signal at low frequencies. A schematic of the problem under consideration is depicted in Figure 1. It consists of a homogeneous porous material that occupies the region 0 ≤ x ≤ L. An incidence acoustical wave impinges normally on the medium. It generates an acoustic pressure field p(x) and an acoustic velocity field v(x) within the material. To derive the reflection coefficient, it is assumed that the pressure field and flow velocity are continuous at the material boundary

p (0^{-}) = p (0^{+}), p (L^{-}) = p (L^{+})

(3)

v (0^{-}) = ϕ v (0^{+}), v (L^{+}) = ϕ v (L^{-})

(4)

where the ± superscript denotes the limit from right and left, respectively. In these conditions, the reflection coefficient is given by (Appendix 1).

R (ω) = \frac{(1 - Y {(ω)}^{2}) sin h (j k (ω) L)}{2 Y (ω) \cosh (j k (ω) L) + (1 + Y^{2} (ω)) sin h (j k (ω) L)}

(5)

where

Y (ω) = ϕ \sqrt{\frac{β (ω)}{α (ω)}} and k (ω) = ω \sqrt{\frac{ρ α (ω) β (ω)}{K_{a}}}

(6)

Figure 1.

Geometry of the problem.

The Darcy’s regime^15,23 corresponds to the range of frequencies such that viscous skin thickness $δ = \sqrt{2 η / ω ρ}$ is much larger than the radius of the pores r, $\frac{δ}{r} ≫ 1$ . This is also called the low frequency range. By doing the Taylor series expansion of the reflection coefficient (equation (5)) at low frequency, we obtain

R (ω) = \frac{1}{1 + \frac{2 \sqrt{ρ K_{a}}}{σ L}} (1 - j \frac{ω}{ω_{0}} + θ (ω^{2}))

(7)

where

ω_{0} = \frac{σ}{2 γ ρ ϕ} \frac{(1 + \frac{2 \sqrt{ρ K_{a}}}{σ L})}{(1 + \frac{1}{γ ϕ^{2}} \frac{\sqrt{ρ K_{a}}}{σ L} (- α_{0} + γ ϕ^{2} (1 + \frac{1}{3} \frac{σ^{2} L^{2}}{ρ K_{a}})))}

(8)

As a first approximation, at very low frequencies (ω → 0), the reflection coefficient (7) is given by the first term

R (ω) = (\frac{1}{1 + \frac{2 C_{1}}{L C_{2}}}) = \frac{1}{1 + \frac{2 \sqrt{ρ K_{a}}}{L σ}}

(9)

This simplified expression of the reflection coefficient is independent on the frequency and porosity of the medium and depends only on its resistivity σ. As a second approximation, at low frequencies, the reflection coefficient (7) is given by the first and second term

R (ω) = \frac{1}{1 + \frac{2 \sqrt{ρ K_{a}}}{σ L}} (1 - j \frac{ω}{ω_{0}})

(10)

In addition to the resistivity, the expression (10) of the reflection coefficient also depends on the porosity ϕ and the viscous tortuosity α₀.

The incident and scattered fields are related by the reflection operator for the material. These are integral operator represented by

\begin{matrix} p^{r} (x, t) = \int_{0}^{t} \tilde{R} (τ) p^{i} (t - τ + \frac{x}{c_{0}}) d τ \\ = \tilde{R} (t) * p^{i} (t) * δ (t + \frac{x}{c_{0}}) \end{matrix}

(11)

In equation (11), * denotes the convolution operation, $p^{i} (t)$ is the incident field, δ is the Delta function, the functions $\tilde{R} (t)$ are the reflection kernels, and its temporal expression is obtained by taking the inverse Fourier transform of equation (10)

\tilde{R} (t) = (\frac{1}{1 + \frac{2 C_{1}}{L C_{2}}}) (δ (t) - \frac{δ (1, t)}{ω_{0}})

(12)

where

δ (1, t)

represents the first derivative of Dirac function.

Inverted values on the porosity and the viscous tortuosity

The reflected acoustic waves by porous material established in low frequency domain are characterized, in addition to the flow resistivity, by two parameters namely, porosity ϕ and viscous tortuosity α₀, the values of which are crucial for 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 determination of these two parameters is done by minimizing between the experimental reflected signal and the simulated reflected signal given by equation (11) using the expression (12) of the reflected kernel. A numerical solution of the least-square method can be found, which minimizes

U (ϕ, α_{0}) = \sum_{i = 1}^{i = N} {(p^{r} (x, t_{i}) - p_{e x p}^{r} (x, t_{i}))}^{2}

(13)

where

p^{r} (x, t_{i})

represents the discrete set of values of the simulated reflected signal given by equation (11) and

p_{e x p}^{r} (x, t_{i})

is the discrete set of values of the experimental reflected signal. The inverse problem is to find the parametric vector

V = {ϕ, α_{0}}

, which satisfies the conditions

{\begin{array}{l} U (ϕ, α_{0}) \to 0 \\ L V \leq V \leq U V \end{array}

(14)

where

L V

and

U V

are the lower and upper bounds that limit the research domain on the adjustable parametric vector

V

. By definition, the value of the porosity is between 0.and 1. (0.< ϕ < 1.) and the viscous tortuosity α₀ cannot be lower than 1. According to these conditions on these different parameters, the lower and upper bounds in equation (14) can be built from the following constraints

{\begin{array}{l} 0 < ϕ < 1 \\ α_{0} \geq 1 \end{array}

(15)

The inverse problem is solved by the last-square method. For its iterative solution, we used the simplex search method²⁹ which does not require numerical or analytic gradient. The tube length of the experimental setup is adaptable to avoid reflection, for measurements in the frequency range (100–300) Hz, a length of 50 m is sufficient but it is useful to use an anechoic device placed at the end of the pipe. The tube diameter is 5 cm (the cut-off of the tube f_c ∼ 4 kHz). A sound source Driver unit “Brand” constituted by loudspeaker Realistic 40–9000 is used. Tone-bursts are provided by Standard Research Systems Model DS345–30 MHz synthesized by the function generator. The signals are amplified and filtered using model SR 650-Dual channel filter, Standford Research Systems. The incident and reflected signals are measured using the same microphone (Bruel and Kjaer: 4190). The incident signal is measured by putting a total reflector in the same position as the porous sample.

The inverse problem is solved for two cylindrical samples of plastic foam S1 and S2 of a diameter of 5 cm. These materials are frequently used in automotive, aerospace and building applications for thermal and sound absorption. The first sample S1 is less resistive with thickness L₁ = (5.00± 0.10) cm, porosity ϕ₁ = (0.92 ± 0.02) and resistivity σ₁ = (6500 ± 500) Nm⁻⁴s. The second sample S2 is more resistive with thickness L₂ = (4.15 ± 0.10) cm, porosity ϕ₂ = (0.82 ± 0.02) and resistivity σ₂ = (30,000 ± 2000) Nm⁻⁴s. The porosity and the flow resistivity of the samples S1 and S2 are measured using direct methods^{7–10,12–13,30–33} and indirect methods by solving the inverse problem.^14,15 Figure 2(a) and (b) shows the experimental incident signal (dashed line) generated by the loudspeaker in the frequency bandwidth (170–210) Hz, and the experimental reflected signal (solid line) and their spectrum (at right) of the samples S1 and S2, respectively. After solving the inverse problem numerically for the porosity ϕ and the viscous tortuosity α₀, we find the following optimized values presented in Table 1. We present in Figure 3(a) and (b), the variation of the minimization function U with the porosity and viscous tortuosity of the samples S1 and S2. In Figure 4(a) and (b), we show a comparison between an experimental reflected signal and simulated reflected signal for the optimized values of porosity and viscous tortuosity. It can be seen that the agreement between experiment and theory curves is good for both the samples S1 and S2, which let us to conclude that the inversion results are best.

Figure 2.

(a) Experimental incident (dashed line) and reflected (solid line) signals at left and their spectrum at right of the sample (S1). (b) Experimental incident (dashed line) and reflected (solid line) signals at left and their spectrum at right of the sample (S2).

Table 1.

Inverted values on the porosity and viscous tortuosity of the samples (S1) and (S2) in different low frequencies.

Sample	Frequency (Hz)	Viscous tortuosity α₀	Porosity ϕ
S1	180–210	2.7	0.93
S1	90–110	2.5	0.92
S2	180–210	8.0	0.81
S2	90–110	8.0	0.86

Figure 3.

(a) Variation of the cost function U with the porosity ϕ and the viscous tortuosity α₀ for the sample (S1). (b) Variation of the cost function U with the porosity ϕ and the viscous tortuosity α₀ for the sample (S2).

Figure 4.

(a) Comparison between experimental reflected signal (dashed line) and simulated reflected signal (solid line) for the optimized values of the porosity and the viscous tortuosity of the sample (S1). (b) Comparison between experimental reflected signal (dashed line) and simulated reflected signal (solid line) for the optimized values of the porosity and the viscous tortuosity of the sample (S2).

Conclusion

In this work, a simplified expression of the reflection coefficient that depends on the flow resistivity, porosity and viscous tortuosity is established at low frequency, an experimental determination of the porosity and the viscous tortuosity of air-saturated porous sample is given by solving the inverse problem using experimental reflected signals. The simulated signals reconstructed using the optimized values found are in good agreement with the corresponding experimental signals. The most important result in this study is that it is now possible to measure, simultaneously, the porosity and the viscous tortuosity by just using the experimental reflected waves at low frequencies.

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) received no financial support for the research, authorship, and/or publication of this article.

References

Chandra

Tuhkuri

SJ.

Active structural acoustic control of transmitted sound through a double panel partition by weighted sum of spatial gradients. J Low Freq Noise Vib Active Control 2017; 36: 27–42.

Lannace

Acoustic correction of monumental churches with ceramic material: the case of the Cathedral of Benevento (Italy). J Low Freq Noise Vib Active Control 2016; 35: 230–239.

Chandra

Tuhkuri

Reddy

JN.

Active structural acoustic control of a softcore sandwich panel using multiple piezoelectric actuators and Reddy's higher order theory. J. Low Freq Noise Vib Active Control 2015; 34: 385–411.

Bagha

Modak

SV.

Active structural-acoustic control of interior noise in a vibro-acoustic cavity incorporating system identification. J Low Freq Noise Vib Active Control 2017; 36: 261–276.

Poletti

MA.

Active acoustic systems for the control of room acoustics. Noise Vib Worldwide 2013; 44: 10–26.

Asdrubali

New challenges in building acoustics. Build Acoust 2017; 24: 3–4.

Beranek

LL.

Acoustic impedance of porous materials.

J Acoust Soc Am 1942; 13: 248

Leonard

RW.

Simplified porosity measurements. J Acoust Soc Am 1948; 20: 39

Guyon

Oger

Plona

TJ.

Transport properties in sintered porous media composed of two particle sizes. J Phys D 1987; 20: 1637

10.

Champoux

Stinson

Daigle

GA.

Air-based system for the measurement of the porosity. J Acoust Soc Am 1990; 89: 910–916.

11.

Castagnède

Saeid

Moussatov

et al . Reflection and transmission at normal incidence onto air-saturated porous materials and direct measurements based on parametric demodulated ultrasonic waves. Ultrasonics 2006; 44: 221–229.

12.

Delany

Bazley

EN.

Acoustical properties of fibrous absorbent materials. Appl Acoust 1970; 3: 105–116.

13.

Miki

Acoustical properties of porous materials – modifications of Delany-Bazley models. Jacoustsocjpn 1990; 11: 19–24.

14.

Sebaa

Fellah

ZEA

Fellah

et al . Measuring flow resistivity of porous material via acoustic reflected waves. J Appl Phys 2005; 98: 084901.

15.

Fellah

ZEA

Fellah

Sebaa

et al . Measuring flow resistivity of porous materials at low frequencies range via acoustic transmitted waves. J Acoust Soc Am 2006; 119: 1926–1928.

16.

Atallah

Panneton

Inverse acoustical characterization of open cell porous media using impedance tube measurements. Canad Acoust 2015; 33: 11–24.

17.

Norris

AN.

On the viscodynamic operator in Biot’s equations of poroelasticity. J. Wave Mat Interact 1980; 1: 365.

18.

Lafarge

Sound propagation in porous materials having a rigid frame saturated by gas.(in French). PhD Dessertation, Université du Maine, Le Mans, France, 1993.

19.

Lee

Leamy

Nadler

JH.

Frequency band structure and absorption predictions for multi-periodic acoustic composites. J Sound Vib 2010; 329: 1809–1822.

20.

Biot

MA.

The theory of propagation of elastic waves in fluid-saturated porous solid. I. Low frequency range. J Acoust Soc Am 1956; 28: 168–178.

21.

Biot

MA.

The theory of propagation of elastic waves in a fluidsaturated porous solid. II. Higher-frequency range. J Acoust Soc Am 1956; 28: 179–191.

22.

Zwikker

Kosten

CW.

Sound-absorbing materials. Amsterdam: Elsevier, 1949.

23.

Allard

JF.

Propagation of sound in porous media: modeling sound absorbing materials. London: Chapman and Hall, 1993. [Database][Mismatch]

24.

Bardot

Brouard

Allard

JF.

Frame decoupling at low frequencies in thin porous layers saturated by air. J Appl Phys 1996; 79: 8223–8229.

25.

Johnson

Koplik

Dashen

Theory of dynamic permeability and tortuosity in fluid-saturated porous media. J Fluid Mech 1987; 176: 379–402.

26.

Allard

Atalla

Propagation of sound in porous media: modeling sound absorbing materials. Chichester, UK: Wiley, 2009, Chap. 5.

27.

Lafarge

Lemarinier

Allard

et al . Dynamic compressibility of air in porous structures at audible frequencies. J Acoust Soc Am 1995; 102(4): 1995–2006.

28.

Pride

Moran

Gangi

AF.

Drag forces of porous-medium acoustics. Phys Rev 1993; 47: 4964–4978.

29.

Lagarias

Reeds

Wright

et al . Convergence properties of the Nelder-Mead Simplex method in low dimensions. SIAM J Optim 1998; 9: 112–147.

30.

Bies

Hansen

CH.

Flow resistance information for acoustical design. Appl Acoust 1980; 13: 357–391.

31.

del Rey

Alba

Arenas

et al . Evaluation of two alternative procedures for measuring airflow resistance of sound absorbing materials. Arch Acoust 2013; 38: 547–554.

32.

Brown

Bolt

RH.

The measurement of flow resistance of porous acoustic materials. J Acoust Soc Am 1942; 13: 337–344.

33.

Stinson

Daigle

GA.

Electronic system for the measurement of flow resistance. J Acoust Soc Am 1988; 83: 2422–2428.