Advances in Materials Science and Engineering

Volume 2015, Article ID 896035, 10 pages

http://dx.doi.org/10.1155/2015/896035

## Mechanical Parameters Effects on Acoustic Absorption at Polymer Foam

^{1}Hyperfrequencies and Semiconductor Laboratory, Faculty of Sciences of Technology, Mentouri Brothers University of Constantine, P.B. 325, 25017 Constantine, Algeria^{2}Faculty of Medical Science, Constantine 3 University, P.B. 125, 25000 Constantine, Algeria

Received 20 June 2015; Revised 9 August 2015; Accepted 10 August 2015

Academic Editor: Belal F. Yousif

Copyright © 2015 Lyes Dib et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### 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 [2–7]. 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 [8–10]. 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 [11–13].

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.