Next Article in Journal
Numerical Investigations on Seismic Behavior of Segmental Assembly of Concrete Filled Steel Tube Piers with External Replaceable Energy-Dissipating Links
Next Article in Special Issue
Fretting Fatigue Life Prediction of Dovetail Structure Based on Plastic Effect and Sensitivity Analysis of Influencing Factors
Previous Article in Journal
An Innovative Method to Analyze the Hydraulic Fracture Reopening Pressure of Hot Dry Rock
Previous Article in Special Issue
Estimation of the Acoustic Properties of the Random Packing Structures of Granular Materials: Estimation of the Sound Absorption Coefficient Based on Micro-CT Scan Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Model for Estimating the Sound Absorption Coefficient in Grid Network Structures

1
FUKOKU Co., Ltd., 6 Showa Chiyoda-machi, Oura-gun, Gunma 370-0723, Japan
2
Department of Engineering, Niigata University, Ikarashi 2-nocho 8050, Nishi-ku, Niigata 950-2181, Japan
3
Graduate School of Science and Technology, Niigata University, Ikarashi 2-nocho 8050, Nishi-ku, Niigata 950-2181, Japan
*
Author to whom correspondence should be addressed.
Materials 2023, 16(3), 1124; https://doi.org/10.3390/ma16031124
Submission received: 2 December 2022 / Revised: 21 January 2023 / Accepted: 24 January 2023 / Published: 28 January 2023
(This article belongs to the Special Issue Advanced Materials Structures for Sound and Vibration Damping)

Abstract

:
Although grid network structures are often not necessarily intended to absorb sound, the gaps between the rods that make up the grid network are expected to have a sound absorption effect. In this study, the one-dimensional transfer matrix method was used to develop a simple mathematical model for accurately estimating the sound absorption coefficient of a grid network structure. The gaps in the grid network structure were approximated as the clearance between two parallel planes, and analysis units were derived to consider the exact geometry of the layers. The characteristic impedance and propagation constant were determined for the approximated gaps and treated as a one-dimensional transfer matrix. The transfer matrix obtained for each layer was used to calculate the sound absorption coefficient. The samples were fabricated from light-curing resin by using a Form2 3D printer from Formlabs. The measurement results showed that a sound absorption coefficient of 0.81 was obtained at the peak when seven layers were stacked. A sensitivity analysis was carried out to investigate the influence of the rod diameter and pitch. The simulated values tended to be close to the experimental values. The above results indicate that the mathematical model used to calculate the sound absorption coefficient is sufficiently accurate to predict the sound absorption coefficient for practical application.

1. Introduction

Grid network structures formed by a large number of round or square rods are often used as screens in machine openings and in natural or forced ventilation inlets and outlets. Various experimental studies have investigated the fluid flow characteristics in grid network structures, such as the oscillatory pressure drop when airflow is applied to a mesh screen [1], the use of a wire mesh as a catalyst for soot mitigation [2] and the pressure loss of wire mesh filters [3].
Yamamoto et al. [4] numerically analyzed the acoustic properties of a frame structure, which is similar to a grid network structure. Satoh et al. [5] compared the predicted and experimental results for sound waves incident perpendicular to a group of cylinder axes, including the sound absorption coefficient. In the grid network structure, such sound waves traverse a direction where the size of the gap between the cylinders varies continuously, which is considered to facilitate a high sound absorption coefficient.
In general, thin grid network materials are often treated as acoustically transparent objects. Therefore, they are used in practical applications or experiments to hold acoustic materials. However, if the grid network has a low aperture ratio, an air layer behind it or several stacked layers, its sound absorption may be significant. In a previous study, Iizuka et al. [6] demonstrated the sound absorption effect of laminated wire mesh structures. Dias and Monaragala [7] performed a similar study on the sound absorption of knitted structures and on metamaterials consisting of laminated thin sheets with apertures [8,9,10]. Dias and Monaragala [7] carried out a similar study on the sound absorption of knitted structures.
Although grid network structures are often not necessarily intended to absorb sound, the gaps between the rods that make up the grid network are expected to have a sound absorption effect. Therefore, predicting acoustic properties such as the sound absorption coefficient of a grid network structure is useful in noise engineering because it is a common engineering shape. Being able to predict the sound absorption effect based on the geometric dimensions of the grid network structure and the physical properties of the gas may help in the development of a compact and simple sound absorption mechanism for various applications. It is expected to provide basic knowledge for many applications, such as filters, sifters, and screens. It can also provide insights into the noise control properties of heaters, heat exchangers, catalysts, and other devices. Furthermore, the grid network structure can function as a reflection barrier for electromagnetic waves and as a vibration damper; thus, it can be used as a device that combines these functions with sound reduction.
Although the finite element method [11,12] is useful in the analysis of acoustic metamaterials, the simpler transfer matrix method was used in this work. The one-dimensional transfer matrix method was used to develop a simple mathematical model for accurately estimating the sound absorption coefficient of a grid network structure. The propagation constants and characteristic impedances were derived by using the Navier–Stokes equations and other equations to consider the viscosity of the air. In addition, a simulated analysis was performed to account for the continuous change in the cross-sectional area when the sound waves are incident perpendicular to the rods in a group of cylinder axes [5]. Moreover, a simulated analysis corresponding to the continuous change in the cross-sectional area is attempted by geometrically estimating the surface area of the wall surface and the volume of the void portion in each segmented element correctly and applying it to a two-plane approximation. Finally, the mathematical model was validated by 3D-printing samples and comparing the measured sound absorption coefficients with the predicted values. The measurement results showed that a sound absorption coefficient of 0.81 was obtained at the peak when seven layers were stacked. A sensitivity analysis was performed to evaluate the effects of parameters such as the diameter and pitch on the predicted values.

2. Experimental Validation

2.1. Measurement Samples

Figure 1 shows a schematic of the test samples used for the experimental measurements. Figure 2 shows a photograph of the samples. Table 1 presents the sample specifications. The samples were fabricated from acrylic-based light-curing resin by using a Form2 3D printer from Formlabs Inc., (Somerville, MA, USA). The rods comprising the grid network had a length of 25.7 mm, and the diameter and pitch were varied to evaluate their effects on the sound absorption coefficient. Four samples were produced with different dimensions. The 3D printer used a layer pitch of 0.025 mm, which allows a fabrication error of less than 0.05 mm in rod diameter. However, during the initial fabrication runs, the fabrication error of the rod diameters was biased to be several tens of micrometers smaller. To compensate for this, the diameter of the rod on the computer-aided design (CAD) drawing was increased by several tens of micrometers. This resulted in a finished product with more accurate dimensions, as shown in Table 1. The room temperature at the time of fabrication was 20 °C, and the time required to fabricate one sample was approximately 7 h.
The samples were inserted into a sample holder made of aluminum alloy with internal dimensions of 25.7 mm per side for measurement of the sound absorption coefficient. The end faces of the rods and the gap between the bottoms of the sample holder and sample were filled with Vaseline to eliminate factors other than the sample that could contribute to sound absorption.

2.2. Measurement Equipment

A Brüel & Kjær Type 4206 (Brüel & Kjær Sound & Vibration Measurement, Nærum, Denmark)two-microphone acoustic impedance tube system was used to measure the sound absorption coefficient. Figure 3a shows the measurement system. Sound waves with a sinusoidal signal were generated by a signal generator built into the fast Fourier transform (FFT) analyzer and were radiated into the impedance tube by a loudspeaker. The transfer function between the sound pressure signals of two microphones attached to the impedance tube was measured by a FFT analyzer. An Onosokki DS-3000 FFT (Ono Sokki Co., Ltd., Yokohama, Japan) analyzer was used for the measurement. The measured transfer function was then used to calculate the normal incident sound absorption coefficient in accordance with ISO 10534-2 [13]. The derivation of the normal incident sound absorption coefficient is described in detail in ISO 10534-2.
The critical frequency at which a plane wave forms depends on the inner diameter of the impedance tube. In this study, impedance tubes with an inner diameter of 29 mm were used because the sound absorption coefficient is not high at low frequencies. However, a rectangular tube with side lengths of 25.7 mm was used as the sample holder to ensure uniform rod lengths. Thus, a conversion tube was used to smoothly change the cross-sectional shape of the impedance tube from circular to square. Figure 3b shows a perspective view of the conversion tube. The diagonal length of the square cross-section was approximately 36.3 mm, so the limiting frequency of the plane wave was determined as approximately 5400 Hz.

3. Simulated Analysis

3.1. Transfer Matrix of an Acoustic Element

The transfer matrix method was used to calculate the simulated value of the sound absorption coefficient of a grid network structure [14,15]. Figure 4 shows the Cartesian coordinate system for a gap of thickness b between two parallel planes 1 and 2. The transfer matrix method can be used to calculate the sound pressure and volume velocity for such a gap based on the one-dimensional wave equation. The relationship between the sound pressure p1 and volume velocity Su1 at plane 1 and the sound pressure p2 and volume velocity Su2 at plane 2 can be expressed by using the four-terminal constant t:
[ p 1 S u 1 ] = [ cos h   γ l Z c S sin h   γ l S Z c sin h   γ l cos h   γ l ] [ p 2 S u 2 ] = [ t 11 t 12 t 21 t 22 ] [ p 2 S u 2 ]
where Zc is the characteristic impedance, γ is the propagation constant, S is the area of the gap and l is the length between planes 1 and 2.

3.2. Propagation Constant and Characteristic Impedance Considering Attenuation

The propagation constant and characteristic acoustic impedance of small tubes considering attenuation due to air viscosity have been studied for tubes with circular [16,17,18] and equilateral-triangle cross-sections [19].
In this study, Stinson’s method [19] was applied to consider the attenuation of sound waves by the loss of friction due to the viscosity of the boundary layer near the wall. The propagation constant γ and characteristic acoustic impedance Zc of a gap between two surfaces, as shown in Figure 4, can be derived by solving the Navier–Stokes equations, continuity equation, gas state equation, energy equation and dissipative function representing heat transfer. Here, air was assumed to be a compressible fluid, and the air viscosity was assumed to be constant.
The propagation constant γ can be expressed as follows:
γ = κ ( κ 1 ) B ( s σ ) B ( s ) ,   B ( x ) = 1 tan h x x ,   s = b 2 ρ 0 ω μ
where κ is the specific heat ratio, σ is the square root of the Prandtl number (0.8677), μ is the viscosity (1.869 × 10−5 Pa·s), ω is the angular frequency and ρ0 is the density of the gas between two planes. Then, the characteristic impedance Zc can be expressed as follows:
Z c = p + u ¯ +
u+ and p+ are the particle velocity and sound pressure, respectively, of the travelling wave:
u ¯ + = c 0 B ( s ) j κ ( γ k ) β e γ x
p + = P s β e γ x
where k is the wavenumber and Ps is the atmospheric pressure (1.013 × 105 Pa).

3.3. Mathematical Model of the Grid Network Structure

The transfer matrix method was used to calculate the sound absorption coefficient of the grid network structure when the sound wave is incident perpendicular to the axial direction of the rods of the grid network structure. Three analysis units were considered to represent different positions in the grid network structure in relation to the tube wall.

3.3.1. Analysis Unit Surrounded by Four Rods

Figure 5 shows a model of the grid network structure. The first analysis unit (i.e., Unit I) represents an area surrounded by four rods (red rectangle). To obtain the transfer matrix, Unit I is divided into n equal layers in the direction of sound incidence (i.e., x-axis direction), as shown in Figure 6. The rod diameter is Dp, so the thickness of one layer is Dp/n. For the actual analysis, 100 divisions were used because this is when the simulated values converged sufficiently. Then, the clearance between two parallel planes can be approximated by using Allard’s method [20] so that the volume of the void space and the surface area of the rod portion are each equal. Figure 7 approximates a single layer as two parallel planes. The surface area of a wall constituted by one rod is S1. As shown in Figure 8, the central angle θn of the arc constituting S1 can be expressed as
θ n = cos 1 { 1 ( x r + D p n r ) } cos 1 ( 1 x r )
where r is the radius of the rod.
Then, the surface area S1 is given by
S 1 = ( P 2 r sin θ ) r d θ
where P is the pitch of the rod. The volume V of the void space is given by
V = ( P 2 r 2 x 2 ) 2 d x
The thickness in the direction of the sound wave incidence is equal to Dp/n, which is the thickness of one layer. The distance b between the two planes can be expressed as
b = V 2 S 1
where the volume between two parallel planes is V. An open-ended correction was added to the grid network structure because the voids can be regarded as orifices in the x-axis direction. For the cross-section perpendicular to the direction of sound wave incidence, the orifice was assumed to be in the x-axis direction. Then, an open-ended correction length of 0.4 times the radius of a circle equivalent to the area of the smallest gap is applied [21,22]:
Δ l = P D p π × 0.4
The open-ended correction length is added to the length l of the relevant layer.

3.3.2. Analysis Unit Surrounded by Three Rods and Tube Wall

Figure 9 shows an analysis unit surrounded by three rods and a tube wall (Unit II). Figure 10 shows the approximation of Unit II as two parallel planes. The lateral surface area S1 is derived in the same manner as in Unit I. The lateral surface area S2 is derived in the same manner as in Unit I. The side area S2 of the rod in contact with the tube wall is derived as follows:
S 2 = { L t ( N p 1 ) P 2 r sin θ } r d θ
where Lt is the length of one side of the sample holder (i.e., Lt = 25.7 mm).
Figure 11 shows the area S3 of the inner wall of the sample holder, which is the rectangular area comprising the pitch between the rods and layer thickness minus the area of overlap between the rectangle and end faces of the rods. In other words, Aar (area enclosed by the red border) can be calculated as follows:
A a r = x 1 x 2 r 2 x 2 d x
Then, the surface area S3 can be expressed as follows:
S 3 = D p n × P 2 A a r
The volume V′ in Figure 12 can be expressed as follows:
V = { L t ( N p 1 ) P 2 r 2 x 2 } ( P 2 r 2 x 2 ) d x
Then, the distance b between two planes is derived as follows:
b = V S 1 + 2 S 2 + S 3 2

3.3.3. Analysis Unit Surrounded by Two Rods and Tube Wall

Figure 12 shows an analysis unit surrounded by two rods and a tube wall (Unit III). Figure 13 shows the approximation of Unit III as two parallel planes. The lateral area S2 of the rod tangential to the tube wall is derived in the same way as in Unit II. The area of the tube wall S3′ is derived as follows:
S 3 = D p n × L t ( N p 1 ) P 2 A a r
The volume V″ in Figure 13 can be expressed as follows:
V = { L t ( N p 1 ) P 2 r 2 x 2 } 2 d x
Then, the distance b between the two planes can be derived as follows:
b = V S 2 + S 3

3.4. Transfer Matrix of the Grid Network Structure

3.4.1. Transfer Matrix of Analysis Units

The two parallel planes of each layer in an analysis unit correspond to the planes in Figure 4. Thus, Equations (2) and (3) can be used to determine the propagation constant and characteristic impedance for each layer. The propagation constant and characteristic impedance can then be substituted into Equation (1) to obtain the transfer matrices T1, T2, T3, …, Tn for the layers, where the subscript 1 represents the layer at the plane incident to the sound wave, and the subscript n represents the number of layers. These transfer matrices can then be cascaded, as shown in Figure 14a, to obtain the transfer matrix for an analysis unit.
Here, Tu is the transfer matrix of Unit I, Tw is the transfer matrix of Unit II and Tf is the transfer matrix of Unit III.

3.4.2. Transmission Matrix of the Whole Sample

Figure 15 shows that the transfer matrices of each analysis unit are connected in parallel along the y-axis and then along the z-axis. This allows the transfer matrix Tl per grid network layer to be derived.
As shown in Figure 15 and Figure 16, the transfer matrices TM and TN are calculated for each analysis unit in parallel in the y-axis direction and for one row in the y-axis direction, respectively. The number of parallel connections of Tw in TM and Tu in TN is one less than the number of rods in the y-axis direction. As shown in Figure 17, the calculated TM and TN are connected in parallel in the z-axis direction. Then, the transfer matrix Tl of the entire sample can be calculated. In this case, the number of TN in parallel is one less than the number of rods in the z-axis direction.
Next, the transfer matrix Tall for the entire sample is derived. As shown in Figure 18, the transfer matrix Tl for one grid network layer and the transfer matrix Ta corresponding to the air layer are connected in an alternating cascade. This allows the entire transfer matrix to be derived.

3.5. Derivation of the Sound Absorption Coefficient

The transfer matrix can be used to derive the sound absorption coefficient. The end of the gap was assumed to be a rigid wall, so the particle velocity u2 = 0. Then, Equation (1) can be transformed as follows:
[ p 1 S u 1 ] = [ t 11 t 12 t 21 t 22 ] [ p 2 0 ] = [ t 11 p 2 t 21 p 2 ]
As shown in Figure 4, p0 and u0 are the sound pressure and particle velocity, respectively, just outside plane 1. If p0 = p1 and S0u0 = Su1, then Equation (19) can be used to obtain the specific acoustic impedance Z0 looking inside the sample from the plane of incidence:
Z 0 = p 0 u 0 = p 0 u 0 S 0 S 0 = p 1 u 1 S S 0 = t 11 t 21 S 0
The relationship between the specific acoustic impedance Z0 and reflectance R is expressed as follows:
R = Z 0     ρ 0 c 0 Z 0 + ρ 0 c 0
Equation (21) can then be used to obtain the sound absorption coefficient α:
α = 1     | R | 2

4. Results

4.1. Experimental and Simulated Values of Sound Absorption Coefficient

The experimental and simulated sound absorption coefficients were compared for each sample. Figure 19 shows the experimental and simulated values for the rod diameters Dp = 2.3 and 2.5 mm and pitch P = 3.5 mm. Figure 20 shows the results for the rod diameters Dp = 2.8 and 3.0 mm and pitch P = 4.2 mm.
For all samples, the experimental and simulated values generally showed similar trends. The simulated values showed high prediction accuracy, but tended to be lower than the experimental values at peak frequencies. A similar tendency was previously observed by Sakamoto et al. [23] when the one-dimensional transfer matrix method was used to estimate the sound absorption coefficient for other cross-sectional shapes. The difference between the simulated and experimental values may be attributed to the simulated calculation underestimating the sound wave attenuation. The reason for this can be that the particle velocity distribution of the sound wave between the two surfaces was averaged in this simulated analysis.
When sound waves are incident to a gap, the distribution of the x-axial component of the particle velocity in the z-axis direction within the gap generally has a distribution with respect to the distance from the wall. In contrast, the present simulated analysis assumed the x-axial component of the particle velocity in the gap between two planes to have a uniform distribution in the z-axis direction. This assumption may have led to the difference between the experimental and simulated values.

4.2. Sensitivity Analysis

In this section, simulation results are presented for different wire diameters and pitches.
For the sensitivity analysis, the sample shown in Figure 2b (i.e., Dp = 2.5 mm, pitch P = 3.5 mm, number of grid network layers N1 = 7) was used as the baseline. Figure 21 shows the calculated results when the rod diameter Dp was varied. The corresponding experimental result is shown for reference. When Dp was incremented 0.1 mm larger than the actual value (i.e., Dp = 2.6 mm), the simulated and experimental values were almost identical. Therefore, the difference between the simulated and experimental values can be attributed to a difference in the rod diameter of approximately 0.1 mm. One way to reduce the difference between the experimental and simulated values to an acceptable level is to adjust the rod diameter used in the simulated calculation.
Figure 21 also shows that the peak value of the sound absorption coefficient increased with an increasing rod diameter, and the peak frequency tended to shift lower. This may be attributed to the increasing influence of the boundary layer as the gap becomes smaller. In general, the peak becomes sharper as the aperture ratio decreases. Here, the sound absorption peak also became sharper as the rod diameter increased (i.e., as the aperture ratio decreased).
Figure 22 shows the calculated results when the pitch P was varied. The size of the sample tube was changed to match the change in pitch. A smaller pitch tended to increase the peak value of the sound absorption coefficient. This may be because a smaller pitch decreases the gap, which again increases the influence of the boundary layer.
Figure 23 shows the calculated results when the number of grid network layers Nl was varied. The peak value increased with the number of layers, and the peak frequency tended to shift lower. This may be because an increase in the number of grid network layers corresponds to an increase in the thickness of the sound-absorbing material. For a single grid network layer, the sound absorption coefficient was difficult to confirm.
Figure 24 shows the calculation results for one grid network layer when the thickness of the back air layer was varied. Increasing the thickness of the back air layer decreased the peak value of the sound absorption coefficient, and the peak frequency tended to shift lower. This confirmed that a sound absorption effect could be achieved by the back air layer even when there was only one grid network layer.
These results show that the parameters used in the simulated analysis can be varied to explore their effects on the sound absorption coefficient. The comparison between the experimental and simulated results show that the observed trends for the changes in the sound absorption coefficient are reasonable.

5. Conclusions

A mathematical model of a grid network structure consisting of rods was developed to estimate the sound absorption coefficient, and the simulated results were compared with experimental measurements. The gaps in the grid network structure were approximated as the clearance between two parallel planes, and analysis units were derived to consider the exact geometry of the layers. The characteristic impedance and propagation constant were determined for the approximated gaps and treated as a one-dimensional transfer matrix. The transfer matrix obtained for each layer was used to calculate the sound absorption coefficient.
The simulated values tended to be close to the experimental values, but generally underestimated them. This may be because the simulated analysis underestimated the sound wave attenuation. For the peak frequency, the simulated analysis showed high prediction accuracy for the peak frequency. The measurement results show that the rod diameter of 2.5 mm, pitch of 3.5 mm, and the presence of seven layers resulted in a sound absorption coefficient of 0.81 at the first peak.
A sensitivity analysis was carried out to investigate the influence of the rod diameter and pitch. The peak sound absorption coefficient tended to increase with an increasing rod diameter and decreasing pitch. This was attributed to the increasing influence of the boundary layer as the gaps become smaller. The above results indicate that the mathematical model used to calculate the sound absorption coefficient is sufficiently accurate to predict the sound absorption coefficient for practical applications. Thus, it was found that the difference between the simulated and experimental values was approximately equivalent to 0.1 mm in diameter. Therefore, it is possible to reduce the difference between the experimental and calculated values to a practically acceptable level by adjusting the diameter of the rod.

Author Contributions

Conceptualization, S.S.; data curation, T.I., T.S. and K.T.; formal analysis, T.S., T.I., K.I. and K.T.; project administration, S.S.; software, K.I.; supervision, S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by FUKOKU Co., Ltd.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Acknowledgments

This study was conducted in collaboration with FUKOKU Co., Ltd.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, T.S.; Cheng, P. Oscillatory pressure drops through a woven-screen packed column subjected to a cyclic flow. Cryogenics 1996, 36, 333–341. [Google Scholar] [CrossRef]
  2. Banús, E.D.; Sanz, O.; Milt, V.G.; Miró, E.E.; Montes, M. Development of a stacked wire-mesh structure for diesel soot combustion. Chem. Eng. J. 2014, 246, 353–365. [Google Scholar] [CrossRef]
  3. Sun, H.; Bu, S.; Luan, Y. A high-precision method for calculating the pressure drop across wire mesh filters. Chem. Eng. Sci. 2015, 127, 143–150. [Google Scholar] [CrossRef]
  4. Yamamoto, T.; Maruyama, S.; Terada, K.; Izui, K.; Nishiwaki, S. A generalized macroscopic model for sound-absorbing poroelastic media using the homogenization method. Comput. Methods Appl. Mech. Eng. 2011, 200, 251–264. [Google Scholar] [CrossRef]
  5. Satoh, T.; Sakamoto, S.; Unai, S.; Isobe, T.; Iizuka, K.; Tasaki, K.; Nitta, I.; Shintani, T. Sound-absorption coefficient of a pin-holder structure for sound waves incident in the direction perpendicular to the cylinder’s axis. Noise Control Eng. J. 2022, 70, 136–149. [Google Scholar] [CrossRef]
  6. Iizuka, K.; Sakamoto, S.; Sato, T.; Tasaki, K.; Nitta, I.; Mizuno, C. Theoretical Estimation and Experiment of Sound Absorption Coefficient for Mesh Structure. In Proceedings of the JSME, No.227-1, Paper No. C032, Kanazawa, Japan, 5 March 2022. [Google Scholar]
  7. Dias, T.; Monaragala, R. Sound absorption in knitted structures for interior noise reduction in automobiles. Meas. Sci. Technol. 2006, 17, 2499. [Google Scholar] [CrossRef]
  8. Gao, N.; Tang, L.; Deng, J.; Lu, K.; Hou, H.; Chen, K. Design, fabrication and sound absorption test of composite porous metamaterial with embedding I-plates into porous polyurethane sponge. Appl. Acoust. 2021, 175, 107845. [Google Scholar] [CrossRef]
  9. Gao, N.; Wu, J.; Lu, K.; Zhong, H. Hybrid composite meta-porous structure for improving and broadening sound absorption. Mech. Syst. Signal Process. 2021, 154, 107504. [Google Scholar] [CrossRef]
  10. Gao, N.; Wang, B.; Lu, K.; Hou, H. Teaching-learning-based optimization of an ultra-broadband parallel sound absorber. Appl. Acoust. 2021, 178, 107969. [Google Scholar] [CrossRef]
  11. Yang, F.; Wang, E.; Shen, X.; Zhang, X.; Yin, Q.; Wang, X.; Yang, X.; Shen, C.; Peng, W. Optimal Design of Acoustic Metamaterial of Multiple Parallel Hexagonal Helmholtz Resonators by Combination of Finite Element Simulation and Cuckoo Search Algorithm. Materials 2022, 15, 6450. [Google Scholar] [CrossRef] [PubMed]
  12. Yang, X.; Yang, F.; Shen, X.; Wang, E.; Zhang, X.; Shen, C.; Peng, W. Development of Adjustable Parallel Helmholtz Acoustic Metamaterial for Broad Low-Frequency Sound Absorption Band. Materials 2022, 15, 5938. [Google Scholar] [CrossRef] [PubMed]
  13. ISO 10534-2; Acoustics−Determination of Sound Absorption Coefficient and Impedance in Impedance Tubes−Part 2: Transfer-Function Method. International Organization for Standardization: Geneva, Switzerland, 1998.
  14. Sasao, H. A guide to acoustic analysis by Excel-Analysis of an acoustic structural characteristic-(4) Analysis of the duct system silencer by Excel. J. Soc. Heat. Air-Cond. Sanit. Eng. Jpn. 2007, 81, 51–58. (In Japanese) [Google Scholar]
  15. Suyama, E.; Hirata, M. Attenuation Constant of Plane Wave in a Tube: Acoustic Characteristic Analysis of Silencing Systems Based on Assuming of Plane Wave Propagation with Frictional Dissipation Part 1. J. Acoust. Soc. Jpn. 1979, 35, 165–170. (In Japanese) [Google Scholar]
  16. Sakamoto, S.; Hoshino, A.; Sutou, K.; Sato, T. Estimating Sound-Absorption Coefficient and Transmission Loss by the Di-mentions of Bundle of Narrow Holes (Comparison between Theoretical Analysis and Experiments). Trans. Jpn. Soc. Mech. Eng. Ser. C 2013, 79, 4164–4176. (In Japanese) [Google Scholar] [CrossRef] [Green Version]
  17. Tijdeman, H. On the propagation of sound waves in cylindrical tubes. J. Sound Vib. 1975, 39, 1–33. [Google Scholar] [CrossRef]
  18. Stinson, M.R. The propagation of plane sound waves in narrow and wide circular tubes, and generalization to uniform tubes of arbitrary cross—Sectional shape. J. Acoust. Soc. Am. 1991, 89, 550–558. [Google Scholar] [CrossRef]
  19. Stinson, M.R.; Champou, Y. Propagation of sound and the assignment of shape factors in model porous materials having simple pore geometries. J. Acoust. Soc. Am. 1992, 91, 685–695. [Google Scholar] [CrossRef]
  20. Allard, J.F.; Atalla, N. Propagation of Sound in Porous Media: Modeling Sound Absorbing Materials, 2nd ed.; Wiley: Hoboken, NJ, USA, 2009; pp. 45–54. [Google Scholar]
  21. Bolt, H.R.; Labate, S.; Ingård, U. The acoustic reactance of small circular orifices. J. Acoust. Soc. Am. 1949, 21, 94–97. [Google Scholar] [CrossRef]
  22. Benade, A.H. Measured End Corrections for Woodwind Tone—Holes. J. Acoust. Soc. Am. 1967, 41, 1609. [Google Scholar] [CrossRef]
  23. Sakamoto, S.; Higuchi, K.; Saito, K.; Koseki, S. Theoretical analysis for sound-absorbing materials using layered narrow clearances between two planes. J. Adv. Mech. Des. Syst. Manuf. 2014, 8, JAMDSM0036. [Google Scholar] [CrossRef]
Figure 1. Schematic of a test sample. Dp: rod diameter; P: rod pitch; and Nm: number of grids.
Figure 1. Schematic of a test sample. Dp: rod diameter; P: rod pitch; and Nm: number of grids.
Materials 16 01124 g001
Figure 2. Test samples in sample holder: (a) Dp = 2.3 mm, P = 3.5 mm, Nm = 64; (b) Dp = 2.5 mm, P = 3.5 mm, Nm = 64; (c) Dp = 2.8 mm, P = 4.2 mm, Nm = 49; (d) Dp = 3.0 mm, P = 4.2 mm, and Nm = 49.
Figure 2. Test samples in sample holder: (a) Dp = 2.3 mm, P = 3.5 mm, Nm = 64; (b) Dp = 2.5 mm, P = 3.5 mm, Nm = 64; (c) Dp = 2.8 mm, P = 4.2 mm, Nm = 49; (d) Dp = 3.0 mm, P = 4.2 mm, and Nm = 49.
Materials 16 01124 g002
Figure 3. Schema of absorption coefficient measurement: (a) two-microphone impedance tube and (b) converter tube to convert the cross-sectional shape from circular to square.
Figure 3. Schema of absorption coefficient measurement: (a) two-microphone impedance tube and (b) converter tube to convert the cross-sectional shape from circular to square.
Materials 16 01124 g003
Figure 4. Cartesian coordinate system for the parallel clearance between a pair of planes.
Figure 4. Cartesian coordinate system for the parallel clearance between a pair of planes.
Materials 16 01124 g004
Figure 5. Analysis unit (Unit I) of the grid network structure.
Figure 5. Analysis unit (Unit I) of the grid network structure.
Materials 16 01124 g005
Figure 6. Analysis unit (Unit I) divided into layers of equal thickness in the x-axis direction.
Figure 6. Analysis unit (Unit I) divided into layers of equal thickness in the x-axis direction.
Materials 16 01124 g006
Figure 7. Approximation of analysis unit (Unit I) as the clearance between two planes.
Figure 7. Approximation of analysis unit (Unit I) as the clearance between two planes.
Materials 16 01124 g007
Figure 8. Derivation of surface area S1.
Figure 8. Derivation of surface area S1.
Materials 16 01124 g008
Figure 9. Analysis unit (Unit II) of the grid network structure.
Figure 9. Analysis unit (Unit II) of the grid network structure.
Materials 16 01124 g009
Figure 10. Approximation of Unit II as the clearance between two planes.
Figure 10. Approximation of Unit II as the clearance between two planes.
Materials 16 01124 g010
Figure 11. Derivation of surface area S3.
Figure 11. Derivation of surface area S3.
Materials 16 01124 g011
Figure 12. Analysis unit (Unit III) of the grid network structure.
Figure 12. Analysis unit (Unit III) of the grid network structure.
Materials 16 01124 g012
Figure 13. Approximation of Unit III as the clearance between two planes.
Figure 13. Approximation of Unit III as the clearance between two planes.
Materials 16 01124 g013
Figure 14. Equivalent circuit analysis: (a) cascade connecting Tn, (b) summation of all Tn for Unit I, (c) summation of all Tn for Unit II, and (d) summation of all Tn for Unit III.
Figure 14. Equivalent circuit analysis: (a) cascade connecting Tn, (b) summation of all Tn for Unit I, (c) summation of all Tn for Unit II, and (d) summation of all Tn for Unit III.
Materials 16 01124 g014
Figure 15. Parallel connection of transfer matrices in a grid network structure.
Figure 15. Parallel connection of transfer matrices in a grid network structure.
Materials 16 01124 g015
Figure 16. Equivalent circuit of parallel connections in the y-axis direction: (a) Tw and Tf and (b) Tu and Tw.
Figure 16. Equivalent circuit of parallel connections in the y-axis direction: (a) Tw and Tf and (b) Tu and Tw.
Materials 16 01124 g016
Figure 17. Equivalent circuit of the whole sample as parallel connections of TM and TN in the z-axis direction.
Figure 17. Equivalent circuit of the whole sample as parallel connections of TM and TN in the z-axis direction.
Materials 16 01124 g017
Figure 18. Equivalent circuit of analysis unit: (a) cascade connecting Tl, Ta and Tall; and (b) cross-section of the grid network structure.
Figure 18. Equivalent circuit of analysis unit: (a) cascade connecting Tl, Ta and Tall; and (b) cross-section of the grid network structure.
Materials 16 01124 g018
Figure 19. Comparison between experimental and simulated values (Dp = 2.3, 2.5 mm, P = 3.5 mm). The rod diameters used in the calculations are the measured values given in Table 1.
Figure 19. Comparison between experimental and simulated values (Dp = 2.3, 2.5 mm, P = 3.5 mm). The rod diameters used in the calculations are the measured values given in Table 1.
Materials 16 01124 g019
Figure 20. Comparison between experimental and simulated values (Dp = 2.8, 3.0 mm, P = 4.2 mm). The rod diameters used in the calculations are the measured values given in Table 1.
Figure 20. Comparison between experimental and simulated values (Dp = 2.8, 3.0 mm, P = 4.2 mm). The rod diameters used in the calculations are the measured values given in Table 1.
Materials 16 01124 g020
Figure 21. Comparison between experimental and simulated values with changing rod diameter (Dp ≂ 1.5–3.0 mm, P = 3.5 mm).
Figure 21. Comparison between experimental and simulated values with changing rod diameter (Dp ≂ 1.5–3.0 mm, P = 3.5 mm).
Materials 16 01124 g021
Figure 22. Comparison between experimental and simulated values with changing pitch (P ≂ 3.3–3.7 mm, Dp = 2.5 mm).
Figure 22. Comparison between experimental and simulated values with changing pitch (P ≂ 3.3–3.7 mm, Dp = 2.5 mm).
Materials 16 01124 g022
Figure 23. Comparison between experimental and simulated values with changing number of layers (Nl ≂ 1–9).
Figure 23. Comparison between experimental and simulated values with changing number of layers (Nl ≂ 1–9).
Materials 16 01124 g023
Figure 24. Simulated values with a changing length of back air space (1–9 mm) (Nl = 1, Dp = 2.5 mm, P = 3.5 mm).
Figure 24. Simulated values with a changing length of back air space (1–9 mm) (Nl = 1, Dp = 2.5 mm, P = 3.5 mm).
Materials 16 01124 g024
Table 1. Specifications of test samples.
Table 1. Specifications of test samples.
Diameter of Rods [mm]
(Measured)
Pitch of Rods [mm]Height of Test Sample [mm]Number of Grids NmNumber of Layers NlCorrespondence to Figure
2.3 (2.31)3.525.7647Figure 2a
2.5 (2.51)Figure 2b
2.8 (2.80)4.2496Figure 2c
3.0 (3.02)Figure 2d
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Satoh, T.; Sakamoto, S.; Isobe, T.; Iizuka, K.; Tasaki, K. Mathematical Model for Estimating the Sound Absorption Coefficient in Grid Network Structures. Materials 2023, 16, 1124. https://doi.org/10.3390/ma16031124

AMA Style

Satoh T, Sakamoto S, Isobe T, Iizuka K, Tasaki K. Mathematical Model for Estimating the Sound Absorption Coefficient in Grid Network Structures. Materials. 2023; 16(3):1124. https://doi.org/10.3390/ma16031124

Chicago/Turabian Style

Satoh, Takamasa, Shuichi Sakamoto, Takunari Isobe, Kenta Iizuka, and Kastsuhiko Tasaki. 2023. "Mathematical Model for Estimating the Sound Absorption Coefficient in Grid Network Structures" Materials 16, no. 3: 1124. https://doi.org/10.3390/ma16031124

APA Style

Satoh, T., Sakamoto, S., Isobe, T., Iizuka, K., & Tasaki, K. (2023). Mathematical Model for Estimating the Sound Absorption Coefficient in Grid Network Structures. Materials, 16(3), 1124. https://doi.org/10.3390/ma16031124

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop