Abstract

Polymer foams have acoustic absorption properties that play an important role in reducing noise level. When the skeleton is set to motion, it is necessary to use generalized Biot-Allard model which takes into account the deformation of the skeleton and the fluid and the interactions between them. The aim of this work is to study the quality of acoustic absorption in polyurethane foam and to show the importance of the structural vibration of this foam on the absorption by varying mechanical parameters (Young’s modulus , Poisson’s coefficient , structural damping factor , and the density ). We calculated the absorption coefficient analytically using classical Biot formulation (, ) and numerically using Biot mixed formulation (, ) in 3D COMSOL Multiphysics. The obtained results are compared together and show an excellent agreement. Afterwards, we studied the effect of varying each mechanical parameter independently on the absorption in interval of ±20%. The simulations show that these parameters have an influence on the sound absorption around the resonance frequency .

1. Introduction

Porous materials are materials well known for their promising applications in many areas, for example, in automotive and aeronautics; they are mainly used to reduce noise level. According to their frame state, they can be classified to three types: elastic, rigid, or limp. Metallic foams and fiber layers are common examples of materials having limp or rigid frame, respectively. Because of the huge rigidity in their frames, only longitudinal waves can propagate inside the fluid phase. The “equivalent fluid” model is often used to model these types of materials [1]. This model is characterized by the effective density and the compression modulus. Many works have been done to evaluate these effective properties to predict the behavior of these types of materials [27]. Polymer foams (polyurethane) are well known porous materials with elastic frame. In a poroelastic medium, acoustic wave propagation is described by a generalized Biot-Allard model [810]. Unlike rigid or limp materials, waves can propagate in both phases of the poroelastic medium, that is, a longitudinal acoustic wave in the fluid phase and both longitudinal and transversal waves in the solid phase.

In this work, we are interested in studying the acoustic behavior of poroelastic materials (polymer foams), specifically polyurethane. To know the quality of absorption of this material, it is necessary to have all the properties that define every phase, that is, fluid and solid, and the interconnection between them and to know other parameters that can influence this quality. A lot of works have been done to study these properties, such as the porosity, resistivity, tortuosity, and viscous and thermal characteristic lengths [1113].

In this paper, we calculate the surface impedance and the absorption coefficient versus frequency, and we study the influence of the mechanical parameters on the quality of the absorption. We use the classical (, ) formulation of Biot-Allard model for the analytical calculations and a combination of (, ) formulation and finite element approach for the numerical calculation using COMSOL environment [14, 15].

2. Modeling of Sound Absorbing Materials

In this section, we will present two important formulations, which with we can study and predict the acoustic behavior of poroelastic medium. These formulations are the classic formulation of Biot known as displacement-displacement formulation (, ) which uses 6 variables for 3D space, 3 for solid phase and 3 for fluid phase, and a mixed formulation that uses 4 space variables; it had been developed by Atalla et al. [16]. This formulation is used to describe the fluid phase and the acoustic pressure in the pores. It is called displacement-pressure formulation (, ).

2.1. The (, ) Formulation for Poroelastic Material

For a monochromatic acoustic wave, of pulsation , incident on porous medium with elastic structure, the wave equation in the solid skeleton, and saturating fluid are obtained from the energetic considerations [1]. With the conventional , the equations can be written as follows:(i)In the solid phase, (ii)In the fluid phase,

and are, respectively, the displacement vector in the structure and the macroscopic displacement vector in the saturating fluid. The coefficients, , , and in (3), are defined from mass coupling factors, , , and in (4), and from viscous coupling parameter in (5). ConsiderMass coupling factors [17] are as follows: Viscous coupling parameter [18] is as follows: is the porosity of the considered medium. The elasticity coefficients , , , and introduced by Biot model are evaluated by three Gedenken experiments [1]. In the case where the material which composes the solid skeleton is less compressible, the four elasticity coefficients can be approximated by the following relations: is the compressibility modulus of the solid frame (in vacuum) which can be evaluated by formulation (7) [1]. is the compressibility modulus of the fluid contained in the material pores and corresponds to the dynamic compressibility established for the equivalent fluid (porous material with rigid frame) in the Johnson et al. model [3] (8). Consider and are, respectively, Young’s modulus and the Poisson coefficient of the deformable solid matrix. , and are the characteristic parameters of the studied porous medium. represents the density of the solid frame and is the effective density defined in the case of the equivalent fluid (10). Consider

2.2. Acoustic Wave Propagation

Two categories of waves propagate at a time across the solid skeleton of the material and in the saturating fluid inside the pores: compression waves and shear waves. In order to determine the characteristics associated to these two types of waves, the displacement vectors and are replaced in (1) and (2) with a scalar potential (), in the case of compression waves, and with a vector potential in the case of shear waves. The calculations are detailed in [1].

The results show that only two compression waves propagate simultaneously in the fluid phase and the solid phase of the porous material with elastic structure. The medium is then characterized by two wavenumbers and in (12) and four characteristic impedances and in (14) in the saturating fluid as well as and in (15) in the deformable matrix. Consider where(i)In the fluid phase, (ii)In the solid phase, The reports and (16) between the speed in the solid frame and the speed in air, for the two compression waves, indicate in which medium, solid or fluid, the waves propagate preferentially: In contrast, only one shear wave propagates in the two mediums composing the porous material. These characteristics are given by In air at 18°C, atmospheric pressure  Pa, with density of fluid , sound wave speed , ration of specific heats , the Prandtl number , and air viscosity [1].

2.3. The () Mixed Formulation for Poroelastic Material

From Biot equations, Atalla et al. [16] have implanted an equivalent mixed formulation (). This formulation is valid only for harmonic oscillations. It is derived from the classic formulation which is equivalent mathematically. The equilibrium modified equations (for small harmonic oscillations) are as follows: where the tilde symbol () indicates that the associated physical property is complex and frequency dependent. In (18), is the angular frequency; and , respectively, denote the solid macroscopic displacement vector and the fluid sound pressure. denotes the modified partial stress tensor associated with the skeleton particle and only depends on the displacement of the solid phase. stands for the porosity defined as the ratio between the volume of the fluid phase and the total volume of the porous material, and , , and are given in the previous section. is the effective density given by . The coefficient is given by .

3. Modeling Poroelastic Materials in COMSOL Multiphysics

In this part, we will focus on the implementation of two equations from the Biot mixed formulation () proposed by Atalla in COMSOL Multiphysics. We choose this formulation because it allows reducing the number of liberty degrees to four by a node instead of six liberty degrees when we use the classic Biot formulation. This gives a considerable reduction in calculation time. COMSOL is a tool for finite element analysis designed specifically to treat the multiphysics problems. The user combines a couple of predefined physics modules in COMSOL and introduces additional coupling terms to the constitutive equations. Of course, each physics mode can be used individually in the case of resolution of nonmultiphysics classic problems. Concerning porous materials, COMSOL does not provide a specific module. But, with the fact that these materials have a rigid structure and they can be modeled as equivalent fluid, the fluid dynamic module and acoustic module in COMSOL can be used. Moreover, porous materials have an elastic structure that can not be modeled multiphysically, that is, using solid mechanics and fluid mechanics modules to model solid and fluid phase, respectively. This is not only for not knowing the coupling terms but also for the high coupling in the equilibrium equations. As a result, the need of using EDP module is necessary to implement, for example, either the classic formulation or Biot mixed formulation, which are made in the form of differential equations. In COMSOL, the general form of PDE (for temporal analysis) must be expressed in the following matrix form: where is the flux vector matrix and is the right part of the vector (the two can be functions of spatial coordinates, the unknown variables , and/or their derivatives in space), and is the gradient/divergence operator. The dimensions of these quantities are as follows: Here is the number of equations equal to the unknowns number (), whereas is the space dimension and then depends of the problem; it can be 1, 2, or 3. In Cartesian coordinates, the gradient/divergence operator vector , for , is defined as follows: The limit conditions in the case of PDE in the general form are as follows: where the vector and can be functions of space coordinates, the unknown variable , and/or their spatial derivatives, whereas is the normal unit vector outgoing from the limit surface. These are, respectively, the limit conditions of Dirichlet and Neumann. The term in the Neumann limit condition is a synonym of Lagrange multiplier. To analyze the harmonic behavior of a porous medium, we use the mixed formulation of Biot (18). The latter depends on 4 variables (): the displacements of the solid phase, , and the fluid pressure, . From these two equations, the matrices and , constituting the form (19), are identified as [14, 15, 19] According to the definition in [16], the expression of can be written as where is the Lamé coefficient for the elastic solid.

Define and can be rewritten in detail as

3.1. Poroelastic/Air Coupling

In the case of poroelastic medium bound to an acoustic medium, (27) describes the continuity conditions of the total normal stress, acoustic pressure, and fluid flow. ConsiderHere, is the pressure in the acoustic medium, is the total stress tensor in the poroelastic material, is the displacement of the fluid phase in the () formulation, and is the outward normal unit vector. The detailed expressions for and were given by Atalla et al. [16]. After some substitution, the vectors and can be expressed as

4. Acoustic Absorption Coefficient

In this section, we consider a structure composed of a poroelastic medium glued on a rigid wall (the rigid wall is a condition to the rear boundary); the surface impedance of this structure was introduced by Allard [1] and given by and are the characteristic impedances of the poroelastic medium and the fluid, respectively [1]. are the ratios between the speeds in poroelastic material and fluid, respectively [1]. However, in COMSOL environment, is defined as the ratio of the acoustic pressure and the total velocity at the impinged face [20], and it can be written as The sound absorption of a poroelastic layer glued to a rigid wall and submitted to plane acoustic wave propagating in the air at the surface of the layer at normal incidence is calculated from the surface impedance and the impedance of air as follows [1, 7, 21]:

5. Numerical Results and Discussion

The acoustic properties of polyurethane foams are given in Table 1 [22], the thickness of this foam is 16 mm, and one of its extremities is glued to a rigid wall, while the other is excited with a normal incidence by a monochromatic pressure wave of pulsation . We have calculated analytically the real and imaginary parts of surface impedance, shown in Figures 1 and 2, and the absorption coefficient versus frequency (Figure 3). In order to validate the obtained results we made a comparison with results calculated with COMSOL in 3D. Indeed, we found an excellent agreement between the results calculated by the two methods. Since the foam skeleton is set to motion, which depends on the frequency, the study of the absorption coefficient is made in three characteristic frequency bands centered around the resonance frequency (32) of the skeleton, which is the quart-wave frequency for the materials glued to a rigid wall. In the vicinity of this frequency, the rigidity of the frame can have a huge influence on the absorption coefficient. The foam fluid phase (air) is very light and has compressibility modulus much lesser than the skeleton’s, which allows the estimation of this frequency by simply considering the properties of the frame under vacuum as follows [23]: where is the thickness of the porous layer (here  mm). is the skeleton density (i.e., the density of the porous material in vacuum). The estimation of quart wave resonance frequency for the layer made of this foam is  kHz. From this frequency, we can find three characteristic zones, the low frequencies range [LF] [1, ], the medium frequencies range [MF] [, 2 ], and the high frequencies range [HF] [2 , ]. In Figure 3, the low frequencies zone (between 1 and 615 Hz) and medium frequencies (between 615 and 2462 Hz), small absorption has been observed, while, in the high frequencies zone (between 2462 and  Hz), the absorption coefficient gets a maximal value for certain frequencies and drops slightly each time reaching the maximum.

5.1. Mechanical Behavior of Porous Materials

The skeleton of porous materials generally has a mechanical behavior of a viscoelastic type at room temperature and in the audible frequency range 20 Hz to 20 kHz. Therefore, its response to mechanical stress depends on the time or the stress frequency (or pulse ) and also the temperature; the two variables are closely linked: the apparent stiffness of the polymer decreases with increasing temperature or as the frequency decreases. As part of this paper, the polymeric foam (polyurethane) will be subject to normal conditions of temperature and pressure ( = 18°C and  Pa).

In the case of small deformations, the behavior can be considered linear and described by Hooke’s law using complex variables [24]. Using a vector representation of strain and stress field, we have where , and , respectively, are the stress, strain, and the complex matrix of elasticity of the solid phase (the tilde symbol for a complex and frequency-dependent quantity). In the case of an isotropic model, the elastic matrix can be characterized from the Young modulus and the complex Poisson coefficient. Consider where is the loss factor defined as the ratio of the instantaneous response (real part of the matrix) on the quadrature phase response (imaginary part). If we consider that the skeleton is isotropic-transverse, in this case, five elastic coefficient are necessary [25, 26]. and are, respectively, the moduli of elasticity in the longitudinal and transverse directions. is shear modulus in the plane . is the Poisson coefficient proportional to the deformation in the transverse plane generated by a deformation in the longitudinal plane, and is the Poisson coefficient in the transverse plane.

Because of the difficulty of determining the above five parameters and the real directions of symmetry, we will focus in this paper on an isotropic model. This assumption allows us the use of an analytical model to describe the vibroacoustic behavior of the porous material. In addition, Poisson’s coefficient is considered real and independent of frequency [27]: this is fixed to 0.3 for polymers foams such as polyurethane. So, from this assumption, (34) are written as follows:

5.2. Influence of the Mechanical Parameters on the Acoustic Absorption

In this part, we discuss the importance of the structure vibration of the polyurethane foam on its absorption properties and study the mechanical parameters effects on these properties (i.e., the Young modulus , structural loss coefficient , the Poisson coefficient , and the density ). Figure 4 shows the variations of the absorption coefficient by varying the Young modulus by ; we notice that this variation is found around the frequency of resonance given by (32); when the value of the Young modulus increases or decreases by , the maximum amplitude of resonance increases and decreases with a value of 0.05, we also notice that there is a shift in frequency of this maximum of a value of 100 Hz towards the higher frequencies in the case where the Young modulus increases and towards the lower frequencies when it decreases, without any modification in the quality factor (Figure 5).

In the same way concerning the influence of the Poisson coefficient on the absorption, the same variation has been observed (Figure 6) except there is a small increase in the resonance amplitude as the value of the Poisson coefficient increases by 20% (Figure 7).

In Figure 8, a variation of 0.02 of the resonance amplitude has been observed; a decrease or an increase occurs depending on the variation in the value of the material density, with a resonance frequency shift of 50 Hz (Figure 9), while, in Figure 10, we distinguish a small variation observed at level of the resonance amplitude without any shift in the resonance frequency (Figure 11).

6. Conclusion

In this paper, we presented a study of the efficiency of acoustic absorption in poroelastic medium of type polymer foam (a case of polyurethane foam). Using the generalized Biot-Allard model and combined formulation of Atalla et al., the real and imaginary parts were calculated. The obtained results show an excellent agreement between analytical calculations and those numericals performed using COMSOL environment. These results confirm that the foam has a very good performance in absorbing noise at medium and high frequencies and less performance in low frequencies. Afterwards, simulations have been done in COMSOL Multiphysics to study the effect of mechanical parameters such as the Young modulus , the Poisson coefficient , structural damping factor , and density of material on the absorption coefficient. The simulations showed that the Young modulus, the Poisson coefficient, and the density of material play an important role in the absorption quality of the polyurethane foam particularly around the resonant frequency, while the structural damping factor plays a minor role in the acoustic absorption. Determining experimentally the mechanical parameters is needed to better understand and predict the acoustic behavior of this foam.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.