Next Article in Journal
Experimental Study on the Influence of Different Factors on the Mechanical Properties of a Soil–Rock Mixture Solidified by Micro-Organisms
Next 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
Previous Article in Journal
Experimental Investigations and Effect of Nano-Powder-Mixed EDM Variables on Performance Measures of Nitinol SMA
Previous Article in Special Issue
Optimal Design of Acoustic Metamaterial of Multiple Parallel Hexagonal Helmholtz Resonators by Combination of Finite Element Simulation and Cuckoo Search Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Models and Experiments on the Acoustic Properties of Granular Packing Structures (Measurement of Tortuosity in Hexagonal Close-Packed and Face-Centered Cubic Lattices)

1
Department of Engineering, Niigata University, Ikarashi 2-no-cho 8050, Nishi-ku, Niigata 950-2181, Japan
2
Graduate School of Science and Technology, Niigata University, Ikarashi 2-no-cho 8050, Nishi-ku, Niigata 950-2181, Japan
*
Author to whom correspondence should be addressed.
Materials 2022, 15(20), 7393; https://doi.org/10.3390/ma15207393
Submission received: 6 September 2022 / Revised: 18 October 2022 / Accepted: 19 October 2022 / Published: 21 October 2022
(This article belongs to the Special Issue Advanced Materials Structures for Sound and Vibration Damping)

Abstract

:
In this study, the sound absorption characteristics of hexagonal close-packed and face-centered cubic lattices were estimated by theoretical analysis. Propagation constants and characteristic impedances were obtained by dividing each structure into elements perpendicular to the incident direction of sound waves and by approximating each element to a clearance between two parallel planes. Consequently, the propagation constant and the characteristic impedance were treated as a one-dimensional transfer matrix in the propagation of sound waves, and the normal incident sound absorption coefficient was calculated by the transfer matrix method. The theoretical value of the sound absorption coefficient was derived by using the effective density applied to the measured tortuosity. As a result, the theoretical value was becoming closer to the measured value. Therefore, the measured tortuosity is reasonable.

1. Introduction

The structures filled with granular materials are used to reduce noise, such as in low-noise pavements [1] and ballast tracks [2], owing to their acoustic characteristics. These continuous clearance structures exhibit acoustic properties similar to those of porous materials, and their acoustic properties change depending on the layer thickness, particle size, and packing structure. Therefore, the prediction of the acoustic properties of the structures filled with granular materials based on the particle size, packing structure, and physical properties of gases is useful for their engineering applications.
Various studies have been conducted on the sound absorption characteristics of granular packing structures, including experimental studies on the acoustic properties of loosely packed granular materials [3], the prediction of the acoustic properties of face-centered cubic structures [4], the numerical analysis of multiple regularly packed structures using commercial software [5], and the study of the sound absorption characteristics of random close-packed materials [6]. In addition, the acoustic properties of granular materials [7], the sound absorption coefficient in a powder bed [8,9], and sound absorption due to gaps in the packing structure of simple cubic and hexagonal lattices [10] have been studied.
In this study, the sound absorption characteristics of hexagonal close-packed and face-centered cubic lattices were estimated by theoretical analysis. Propagation constants and characteristic impedances were obtained by dividing each structure into elements perpendicular to the incident direction of sound waves and approximating each element to a clearance between two parallel planes. Consequently, the propagation constant and the characteristic impedance were treated as a one-dimensional transfer matrix in the propagation of sound waves, and the normal incident sound absorption coefficient was calculated by the transfer matrix method. This method measures the sound absorption coefficient using a one-dimensional (1D) wave equation from simplified geometrical information of the granular material. The program and the number of calculations are so simple that they cannot be compared to the finite element method (FEM), for example, and the calculation can be completed in a fraction of a second on a standard PC.
We also measured the tortuosity, which is used to represent the complexity of the path through which sound waves propagate in the structure. The theoretical value of the sound absorption coefficient was derived by using the effective density applied to the measured tortuosity which is one of the Biot parameters. The drawback of this method may be that the theoretical value includes a measurement of tortuosity, but instead, a simplified method for measuring tortuosity is proposed.
In the measurement of the sound absorption coefficient, the normal incident sound absorption coefficient of each sample was measured using a two-microphone impedance tube. The comparison between experimental and theoretical values was reported.

2. Samples and Measuring Device Used to Measure the Sound Absorption Coefficient

2.1. Transmission Loss Measurement

In this study, two types of packing structures, hexagonal close-packed and face-centered cubic lattices, were investigated. The samples used for the measurement of sound absorption are shown in Figure 1a–d, and the sample specifications are shown in Table 1. Stainless steel spheres with diameters of d = 4 mm and 8 mm were used as particles.
In the preparation of the measurement sample, the granular material was packed regularly with the sample holder. The sample holder was regularly arranged with a hemispherical convex part on the wall, and a regular packing structure was formed by placing a granular body in a predetermined position. The sample holder was fabricated in photocurable resin using a photocurable 3D printer Form 2 manufactured by Formlabs Inc. (Somerville, MA, USA).

2.2. Measurement Equipment for Sound Absorption Coefficient

A 4206-type two-microphone impedance tube made by Brüel and Kjær was used to measure the sound absorption coefficient. The configuration of the device is shown in Figure 2. The sample was enclosed in an impedance tube, and a sine wave signal was output by an internal signal generator in the DS-3000 fast Fourier transform (FFT) analyzer manufactured by Ono Sokki, and the transfer function between the sound pressure signals of two microphones attached to the impedance tube was measured by an FFT analyzer. Using the measured transfer function, the normal incident sound absorption coefficient was calculated in accordance with ISO 10534-2. The critical frequency of the plane wave depended on the inner diameter of the acoustic tube. In this study, a small tube with an inner diameter of 29 mm was used because the sound absorption coefficient in the low-frequency range was not high. Therefore, the measurement was performed in the range of 500–6400 Hz.

3. Measurement Method and Results of Tortuosity

3.1. Overview of Tortuosity

Tortuosity, which is a Biot parameter, is the ratio of the average length of voids in the poroelastic material to the thickness of the sound-absorbing material. When a sound wave passes through a sound-absorbing material with a complicated internal clearance structure, the tortuosity represents the complexity of the path of the sound wave. In this study, the tortuosity for each packed structure was measured using an ultrasonic sensor. In general, the tortuosity α is expressed as Equation (1) using the sound velocity c0 in the air and the apparent sound velocity c in the packed structure [11]:
α = c 0 c 2

3.2. Tortuosity Measurement

The tortuosity was calculated by Equation (1) from the square of the ratio of two sound velocities: the sound velocity in the air without a sample and the sound velocity transmitted through the sample. Ultrasonic waves were output from an ultrasonic sensor on the transmitter side using a sinusoidal signal generated by a signal generator. Ultrasonic waves propagated in the sample were measured using an ultrasonic sensor on the receiving side, and the waveform was observed using an oscilloscope. The propagation time of ultrasonic waves from the transmitter to the receiver was calculated from the comparison between the original waveform obtained by the signal generator and that of ultrasonic waves propagated in the sample, and the apparent sound velocity c in the sample in Equation (1) was calculated.
Figure 3 shows the configuration of the tortuosity measurement device. Ultrasonic sensors with center frequencies of 32.7 kHz, 40 kHz, 58 kHz, 110 kHz, 150 kHz, 200 kHz, and 300 kHz were used. First, the tortuosity α at each frequency was measured for each packed structure. Then, the reciprocal of the square root of the frequency used for the measurement was determined as the value of the horizontal axis, and the tortuosity α at each frequency obtained by the measurement was plotted as the value of the vertical axis. The linear approximation of these point clouds leads to a soaring straight line using the least-squares method. When the frequency of the approximate line is set to infinity, the limiting value of the tortuosity, i.e., the y-intercept of the graph, becomes the tortuosity α of the packing structure.
At higher frequencies, ultrasonic sensors have a lower signal-to-noise ratio (S/N ratio or SNR) due to a lower conversion efficiency and a higher attenuation of sound waves in the air. To improve the S/N ratio and measurement accuracy, 150 measurements were added synchronously. The signal was measured in the amplitude direction with a resolution of 16 bits.
Table 2 and Figure 4 show the results of tortuosity measurements. As mentioned above, the y-intercept of the approximate straight line is the tortuosity α of each packing structure. At 40 kHz, the deviation was excluded for the hexagonal structure. The tortuosity of each packed structure was 1.44 for the hexagonal close-packed structure and 1.43 for the face-centered cubic lattice.
Generally, measuring tortuosity requires expensive equipment, and outsourcing the measurement is not cost-effective. In addition, experimental methods with very high ultrasonic frequencies make it difficult to measure tortuosity on rigid granular objects such as those used in this study (where the rigid surfaces have a significant size relative to the wavelength) because ultrasonic waves are reflected.
To avoid the abovementioned problems, this apparatus uses several ultrasonic transducers with low frequencies. Other devices are very common. This simple method of measuring tortuosity may be useful for rigid granular materials such as those used in this study.

4. Theoretical Analysis

4.1. Analysis Units and Element Division

The gap in the packed structure of granular materials was analyzed by the transfer matrix method based on the one-dimensional wave equation. In this section, the outline of the theoretical analysis is explained. Figure 5 shows the analysis units for each packing structure. In these structures, the cross-sectional shape changes periodically in the x-, y-, and z-axis directions. The range, whereby the change in the cross-sectional shape was completed in one cycle, was used as the analysis unit. Namely, the analysis unit was defined as the area surrounded by the broken line in Figure 5. In addition, the analysis unit located at the upper edge of the sample in Figure 5a,b (i.e., the incident surface of the sound wave) in the x direction is shaped, as shown in Figure 6a,b, respectively.
In the packing structure of the granular materials, the shape of the cross-section changes continuously in the plane perpendicular to the direction of sound wave propagation. Therefore, the analysis unit was divided into 100 parts in the direction perpendicular to the x-axis. Figure 7a shows the analysis unit division method for the face-centered cubic lattice shown in Figure 5b. The analysis units shown in Figure 5a and Figure 6a,b were also divided, as in the case of Figure 7a. The number of divisions of the analysis unit, n = 100, is the value at which the theoretical value of the normal incident sound absorption coefficient converges.
For the clearance between particles in the divided element shown in Figure 7a, the real surface area Sn and real volume Vn of the clearance were geometrically calculated and approximated to the thickness of the clearance between two planes bn, as shown in Figure 7b and Equation (2).
For the divided elements in contact with the sample holder, as shown in Figure 8a,b, the area Sh of the wall of the sample holder was also taken into consideration, and it was approximated to the clearance between two planes. As shown in Figure 8a,b, when the clearance thickness is set to bn’, Sn, Sh, and Vn are expressed as Equation (3):
b n = 2 V n S n
b n = 2 V n S n + S h

4.2. Derivation of the Surface Area of a Sphere in a Divided Element

In this section, we describe the derivation method of the surface area Sn of spheres in the divided element for the hexagonal close-packed and face-centered cubic lattices used in Equations (2) and (3).
Using the radius r of the sphere and the length Zunit of the dividing element in the x-axis direction, the surface area Sn of the sphere at the dividing element can be obtained by Equation (4):
S n = 2 π r × Z u n i t
As shown in Figure 8, the area Sh of the blue portion of the sample holder wall was considered for the divided element in contact with the sample holder.

4.3. Derivation of the Volume of the Clearance in a Divided Element

In this section, we describe the derivation method of the volume Vn of clearances in the divided element of the hexagonal close-packed and face-centered cubic lattice structures.
Vn, for the hexagonal close-packed structure, shows the clearance in the range surrounded by the broken line in Figure 5a and is derived geometrically using the dimensions of the radius r of the sphere and the length Zunit of the divided element in the x-axis direction, as follows:
V n = 2 3 r 2 × Z u n i t 2 π r 2 x 2 d x
Vn, for the face-centered cubic structure, shows the clearance in the range surrounded by the broken line in Figure 5b, which is geometrically derived as follows:
V n = 8 r 2 × Z u n i t 4 π r 2 x 2 d x

4.4. Propagation Constant and Characteristic Impedance Considering Tortuosity

Propagation constants and characteristic impedances in the clearance between two planes approximated in Section 4.1 were obtained by considering the attenuation of sound waves. Propagation constants and characteristic impedances considering the viscosity of the air in tubes were studied by Tijdeman [12] and Stinson [13] for circular tubes, Stinson [14] for equilateral triangles, and Beltman [15] for rectangular tubes. In addition, Allard [16] considered the effect of tortuosity. In this study, the methods reported by Stinson [14] and Allard [16] were applied.
The Cartesian coordinate system was used, as shown in Figure 9, and the effective density ρs and compressibility Cs were obtained from a three-dimensional analysis using Equations (7) and (8) [14], respectively, using the Navier–Stokes equations, the equation of state of gas, the continuity equation, the energy equation, and the dissipation function representing heat transfer, where ρ0 is the density of the air, λs is an intermediary variable, bn is the clearance thickness between two planes, ω is the angular frequency, η is the viscosity of the air, κ is the specific heat ratio of the air, P0 is the atmospheric pressure, and Npr is the Prandtl number.
ρ s = ρ 0 1 tanh j λ s j λ s 1 , λ s = b n 2 ω ρ 0 η
C s = 1 κ P 0 1 + κ 1 tanh j N p r λ s j N p r λ s
By using the effective density ρs multiplied by the tortuosity α, the propagation constant and characteristic impedance considering the tortuosity can be obtained [16]. Therefore, the propagation constant γ and the characteristic impedance Zc can be expressed by the following equations [16] in terms of the effective density ρs and compressibility Cs when the tortuosity α is considered:
γ = j ω α ρ s C s
Z c = α ρ s C s

4.5. Transfer Matrix

The clearance between two planes was analyzed by the transfer matrix method with respect to sound pressure and volumetric velocity based on the one-dimensional wave equation. Figure 10 shows a schematic diagram showing one element in the x-direction for the clearance between two planes shown in Figure 9. Using the characteristic impedance, the propagation constant; the cross-sectional area of the clearance, S; and the length of the divided element, l, obtained in Section 4.1,Section 4.2,Section 4.3 and Section 4.4 the transfer matrix, Tn; and the four-terminal constants of the acoustic tube element, A to D, can be calculated using Equation (11):
T n = cosh γ l Z c S s i n h γ l S Z c s i n h γ l cosh γ l = A B C D
Plane 1 is the incident surface of sound waves, and Plane 2 is the transmission surface of sound waves. The sound pressure and the particle velocity can be expressed as p1 and u1 at Plane 1 and p2 and u2 at Plane 2, respectively, and the transfer matrix can be expressed as Equation (12):
p 1 S u 1 = A B C D p 2 S u 2
When Equation (12) was applied to the clearance between the two planes obtained in Section 4.2 and Section 4.3, the transfer matrix was obtained at each divided element. Since each divided element is continuous in the x-axis direction, the transfer matrix Tunit and Ttop of the analysis unit and the analysis unit at the upper edge of the sample were obtained by cascading the transfer matrix of each divided element based on the equivalent circuit shown in Figure 11.
Next, as shown in Figure 12, the transfer matrix corresponding to the analysis unit was connected based on the equivalent circuit of the whole sample, and the transfer matrix of the whole sample was derived.
First, the analysis unit that was arranged in the x-axis direction, i.e., the incident direction of sound waves, was cascaded. Next, the transfer matrix of the whole sample was obtained by connecting the cascade-connected transfer matrix that aligned on the yz plane, the plane perpendicular to the incident direction of the sound wave, in parallel.

4.6. Normal Incident Sound Absorption Coefficient

The sound absorption coefficient was calculated from the transfer matrix Tall obtained in Section 4.5. Since the end of the specimen used in this study was a rigid wall, the particle velocity u2 = 0, Equation (12), can be deformed as shown in Equation (13), yielding Equation (14):
p 1 S u 1 = A B C D p 2 0
p 1 S u 1 = A p 2 C p 2
When the sound pressure and the particle velocity just outside of Plane 1 are p0 and u0, the specific acoustic impedance Z0 from the incident surface of the sample is expressed as follows:
Z 0 = p 0 u 0
Therefore, according to p0 = p1, S0u0 = Su1, and Equation (15), the specific acoustic impedance Z0 of the sample is expressed as
Z 0 = p 0 u 0 = p 0 u 0 S 0 S 0 = p 1 u 1 S S 0 = A C S 0
where S/S0 is the aperture ratio of the sample shown in Table 1.
The relation between the specific acoustic impedance Z0 and the reflectance R is expressed by the following equation:
R = Z 0 ρ 0 c 0 Z 0 + ρ 0 c 0
According to Equation (17), the theoretical value of the normal incident sound absorption coefficient α of the sample is expressed as follows:
α = 1 R 2

5. Comparison between Measured and Theoretical Values

In this section, the theoretical value of the normal incident sound absorption coefficient is compared with the measured value for each sample. Three theoretical values are presented for Figure 13a–d: first, for the case where the tortuosity is not considered (unity); for the tortuosity measured in Section 3.2 (hexagonal close-packed lattice: 1.44, face-centered cubic lattice: 1.43); and for the sound absorption coefficient which derived from the tortuosity obtained by the numerical analysis in [5] (HCP: 1.65, FCC: 1.66). Here, the tortuosity measured in Section 3.2 shows smaller values than those obtained from the numerical analysis in [5].
Figure 13a,b show the hexagonal close-packed structure, and Figure 13c,d show the comparison of the face-centered cubic lattice with grain sizes of d = 4 mm and d = 8 mm, respectively.
First, we compared the measured values with theoretical values without considering the tortuosity. As shown in Figure 13a–d, the peak frequency of the theoretical value, which did not consider the tortuosity, compared with the measured value that appeared on the high-frequency side. Except for Figure 13a, the theoretical sound absorption peak values, without considering the tortuosity, were lower than the measured values.
Next, we focused on the theoretical values that considered the tortuosity. The peak frequency shifted to a lower frequency in Figure 13a–d, and the peak sound absorption value increased in all cases (except for Figure 13a).
As a result, the difference between the theoretical and measured values decreased when the tortuosity was considered. A larger tortuosity than unity means a longer path of sound waves that propagate in the sample, similar to the increase in sample thickness. Therefore, it is considered that the peak frequency of the theoretical value moved to the low-frequency side.
In general, the sound absorption coefficient of porous materials is greatly influenced by thickness [17,18]. In other words, the sound absorption peak occurs at a frequency where the thickness of the porous sound-absorbing material corresponds to one-quarter wavelength of the sound wave. Conversely, if the apparent sound velocity in the material decreases because of boundary layer friction or tortuosity, the peak sound absorption frequency will decrease, as this is equivalent to an increase in the apparent thickness of the material. That is, an increase in tortuosity decreases the peak sound absorption frequency and increases the sound absorption curve in most cases. For the same reason, in Figure 13b–d, it is considered that the theoretical sound absorption peak value increased, owing to the consideration of the tortuosity.
Next, we discuss the tortuosity in each sample. As shown in Section 3, the measured tortuosities in the hexagonal close-packed and face-centered cubic lattices were 1.44 and 1.43, respectively.
By considering one of the two types of tortuosities, the theoretical values of the peak frequency and the sound absorption peak were in good agreement with the measured values for the hexagonal close-packed lattice with a grain size d = 4 mm, as shown in Figure 13a,b. Here, the theoretical values, considering the measured tortuosity in Section 3.2, are closer to the experimental values. However, when the grain size of the hexagonal close-packed structure was d = 8 mm, the peak frequency shifted to a lower frequency with the consideration of either of the tortuosities, but it was lower than the peak frequency side of the measured value because the division and assembly of the sample holder led to errors in the measured values. For the measured hexagonal close-packed structure with a grain size of d = 8 mm, the number of spheres used was small. To avoid interference between the particles during the sample assembly, the holder was divided into three parts perpendicular to the incident direction of the sound wave. Dimensional errors between the divided holders may have caused gaps in the spheres, thereby affecting the measured values. On the other hand, for the grain size d = 4 mm, the accuracy of the theoretical value was improved by considering the measured tortuosity in Section 3.2.
In the case of the face-centered cubic lattice with a grain size of d = 4 mm, as shown in Figure 13c, the theoretical values were fairly close to the measured values when either tortuosity was considered. When the two types of tortuosities were compared, the theoretical values from the tortuosity measured in Section 3.2 agreed very well with the measured values. When the grain size of the face-centered cubic lattice was d = 8 mm, as shown in Figure 13d, by considering either of the tortuosities, the theoretical values of both the peak frequency and the sound absorption peak approached the measured values, but the degree of agreement was smaller than that in Figure 13c. In terms of the agreement between theoretical (applying the tortuosity in Section 3.2) and experimental values, the case with a grain size of 4 mm (Figure 13a,c) is better than the case with a grain size of 8 mm (Figure 13b,d). This is because the method used to estimate the attenuation of sound waves in gaps [14] is more accurate for smaller gaps [19].

6. Conclusions

The normal incidence sound absorption coefficients of hexagonal close-packed and face-centered cubic lattices, which are typical packing structures of granular materials, were theoretically calculated. In the derivation of the propagation constant and the characteristic impedance, the measured tortuosity was taken into account. The results were compared with the experimental data. The following conclusions were obtained:
  • In both packing structures, the real area of the granular surface and the real volume of the clearance were obtained geometrically and analyzed theoretically.
  • In both packing structures, the peak frequency tended to appear at a higher frequency than the measured value when the tortuosity was not considered.
  • In the theoretical sound absorption, the peak value was higher when the tortuosity was considered compared to that without the consideration of the tortuosity (the peak frequency moved to a lower frequency). As a result, the theoretical value was becoming closer to the measured value.

Author Contributions

Conceptualization, S.S. (Shuichi Sakamoto); Data curation, K.T. and S.S. (Shotaro Seino); Formal analysis, K.S.; Project administration, S.S. (Shuichi Sakamoto); Software, K.S.; Supervision, S.S. (Shuichi Sakamoto). All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Japan Society for the Promotion of Science (JSPS) KAKENHI Grant No. 20K04359.

Institutional Review Board Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sandberg, U. Low noise road surfaces—A state-of-the-art review. J. Acoust. Soc. Jpn. 1999, 20, 1–17. [Google Scholar] [CrossRef] [Green Version]
  2. Zhang, X.; Thompson, D.; Jeong, H.; Squicciarini, G. The effects of ballast on the sound radiation from railway track. J. Sound Vib. 2017, 339, 137–150. [Google Scholar] [CrossRef]
  3. Voronina, N.N.; Horoshenkov, K.V. A new empirical model for the acoustic properties of loose granular media. Appl. Acoust. 2003, 64, 415–432. [Google Scholar] [CrossRef]
  4. Gasser, S.; Paun, F.; Bréchet, Y. Absorptive properties of rigid porous media: Application to face centered cubic sphere packing. J. Acoust. Soc. Am. 2005, 117, 2090–2099. [Google Scholar] [CrossRef] [PubMed]
  5. Lee, C.-Y.; Leamy, M.J.; Nadler, J.H. Acoustic absorption calculation in irreducible porous media: A unified computational approach. J. Acoust. Soc. Am. 2009, 126, 1862–1870. [Google Scholar] [CrossRef] [PubMed]
  6. Dung, V.V.; Panneton, R.; Gagné, R. Prediction of effective properties and sound absorption of random close packings of monodisperse spherical particles: Multiscale approach. J. Acoust. Soc. Am. 2019, 145, 3606–3624. [Google Scholar] [CrossRef] [Green Version]
  7. Sakamoto, S.; Sakuma, Y.; Yanagimoto, K.; Watanabe, S. Basic Study for the Acoustic Characteristics of Granular Material (Normal Incidence Absorption Coefficient for Multilayer with Different Grain Diameters). J. Environ. Eng. 2012, 7, 12–22. [Google Scholar] [CrossRef] [Green Version]
  8. Sakamoto, S.; Yamaguchi, K.; Ii, K.; Takakura, R.; Nakamura, Y.; Suzuki, R. Theoretical and experiment analysis on the sound absorption characteristics of a layer of fine lightweight powder. J. Acoust. Soc. Am. 2019, 146, 2253–2262. [Google Scholar] [CrossRef]
  9. Sakamoto, S.; Takakura, R.; Suzuki, R.; Katayama, I.; Saito, R.; Suzuki, K. Theoretical and Experimental Analyses of Acoustic Characteristics of Fine-grain Powder Considering Longitudinal Vibration and Boundary Layer Viscosity. J. Acoust. Soc. Am. 2021, 149, 1030–1040. [Google Scholar] [CrossRef]
  10. Sakamoto, S.; Ii, K.; Katayama, I.; Suzuki, K. Measurement and Theoretical Analysis of Sound Absorption of Simple Cubic and Hexagonal Lattice Granules. Noise Control Eng. J. 2021, 69, 401–410. [Google Scholar] [CrossRef]
  11. Allard, J.F.; Castagnede, B.; Henry, M. Evaluation of tortuosity in acoustic porous materials saturated by air. Rev. Sci. Instrum. 1994, 65, 754. [Google Scholar] [CrossRef] [Green Version]
  12. Tijdeman, H. On the propagation of sound waves in cylindrical tubes. J. Sound Vib. 1975, 39, 1–33. [Google Scholar] [CrossRef]
  13. 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]
  14. 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]
  15. Beltman, W.M.; van der Hoogt, P.J.M.; Spiering, R.M.E.J.; Tijdeman, H. Implementation and experimental validation of a new viscothermal acoustic finite element for acousto-elastic problems. J. Sound Vib. 1998, 216, 159–185. [Google Scholar] [CrossRef]
  16. Allard, J.F. Propagation of sound in Porous Media Modeling Sound Absorbing Materials. J. Acoust. Soc. Am. 1994, 95, 2785. [Google Scholar] [CrossRef] [Green Version]
  17. Tabana, E.; Khavanina, A.; Jafarib, A.J.; Faridanc, M.; Tabrizid, A.K. Experimental and mathematical survey of sound absorption performance of date palm fibers. Heliyon 2019, 5, e01977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Duan, C.; Cui, G.; Xu, X.; Liu, P. Sound absorption characteristics of a high-temperature sintering porous ceramic material. Appl. Acoust. 2012, 73, 865–871. [Google Scholar] [CrossRef]
  19. 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, 16. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Test samples: (a) Hexagonal close-packed lattice (d = 4 mm); (b) Hexagonal close-packed lattice (d = 8 mm); (c) Face-centered cubic lattice (d = 4 mm); (d) Face-centered cubic lattice (d = 8 mm).
Figure 1. Test samples: (a) Hexagonal close-packed lattice (d = 4 mm); (b) Hexagonal close-packed lattice (d = 8 mm); (c) Face-centered cubic lattice (d = 4 mm); (d) Face-centered cubic lattice (d = 8 mm).
Materials 15 07393 g001
Figure 2. Configuration diagram of a two-microphone impedance tube for the absorption coefficient measurement.
Figure 2. Configuration diagram of a two-microphone impedance tube for the absorption coefficient measurement.
Materials 15 07393 g002
Figure 3. Configuration diagram of the tortuosity measurement.
Figure 3. Configuration diagram of the tortuosity measurement.
Materials 15 07393 g003
Figure 4. Measurement results of the tortuosity.
Figure 4. Measurement results of the tortuosity.
Materials 15 07393 g004
Figure 5. Analysis unit in each structure: (a) Hexagonal close-packed lattice; (b) Face-centered cubic lattice.
Figure 5. Analysis unit in each structure: (a) Hexagonal close-packed lattice; (b) Face-centered cubic lattice.
Materials 15 07393 g005
Figure 6. Analysis unit at the top of the sample: (a) Hexagonal close-packed lattice; (b) Face-centered cubic lattice.
Figure 6. Analysis unit at the top of the sample: (a) Hexagonal close-packed lattice; (b) Face-centered cubic lattice.
Materials 15 07393 g006
Figure 7. Divided element approximated to the clearance between two planes: (a) Divided element (face-centered cubic); (b) Approximated clearance between two planes.
Figure 7. Divided element approximated to the clearance between two planes: (a) Divided element (face-centered cubic); (b) Approximated clearance between two planes.
Materials 15 07393 g007
Figure 8. Approximated clearance between two planes at the element in contact with the sample holder: (a) Divided element in contact with the sample holder on two sides; (b) Divided element in contact with the sample holder on one side.
Figure 8. Approximated clearance between two planes at the element in contact with the sample holder: (a) Divided element in contact with the sample holder on two sides; (b) Divided element in contact with the sample holder on one side.
Materials 15 07393 g008
Figure 9. Cartesian coordinate system for the parallel clearance between the two planes.
Figure 9. Cartesian coordinate system for the parallel clearance between the two planes.
Materials 15 07393 g009
Figure 10. Sound incident area, incident plane, and transmission plane of approximated clearance between two planes in Figure 7b and Figure 9.
Figure 10. Sound incident area, incident plane, and transmission plane of approximated clearance between two planes in Figure 7b and Figure 9.
Materials 15 07393 g010
Figure 11. Equivalent circuit in the analysis unit (cascade connecting the transfer matrix of each element).
Figure 11. Equivalent circuit in the analysis unit (cascade connecting the transfer matrix of each element).
Materials 15 07393 g011
Figure 12. Equivalent circuit of the whole sample (parallel connection of the cascaded Tunit).
Figure 12. Equivalent circuit of the whole sample (parallel connection of the cascaded Tunit).
Materials 15 07393 g012
Figure 13. Comparison between the experimental and calculated values (considering tortuosity in Section 3.2 and tortuosity in Reference (Lee et al. 2009) [5]) of the peak frequency and the sound absorption coefficient: (a) Hexagonal close-packed lattice (d = 4 mm); (b) Hexagonal close-packed lattice (d = 8 mm); (c) Face-centered cubic lattice (d = 4 mm); and (d) Face-centered cubic lattice (d = 8 mm).
Figure 13. Comparison between the experimental and calculated values (considering tortuosity in Section 3.2 and tortuosity in Reference (Lee et al. 2009) [5]) of the peak frequency and the sound absorption coefficient: (a) Hexagonal close-packed lattice (d = 4 mm); (b) Hexagonal close-packed lattice (d = 8 mm); (c) Face-centered cubic lattice (d = 4 mm); and (d) Face-centered cubic lattice (d = 8 mm).
Materials 15 07393 g013aMaterials 15 07393 g013b
Table 1. Properties of test samples.
Table 1. Properties of test samples.
Packing StructureDiameter [mm]Length [mm]Aperture Ratio of Sample HolderFilling RateMeasured TortuosityCorrespondence to Figure
Hexagonal close-packed4270.570.741.441a
8270.670.741.441b
Face-centered cubic4220.580.741.431c
8210.850.741.431d
Table 2. Measurement results of the tortuosity.
Table 2. Measurement results of the tortuosity.
Frequency [kHz]Inverse of the Square Root of Frequency
[ 1 / Hz ]
Distance between Sensors [mm]Transmission Time
[ms]
Tortuosity
Hexagonal
Close-Packed
Face-Centered Cubic
32.70.005593951.2294.784.48
400.005395--3.35
580.0041523951.1994.313.98
1100.0030153451.0443.593.36
1500.0025822040.6263.112.84
2000.0022362040.6182.562.67
3000.0018262040.6152.561.93
0--1.441.43
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sakamoto, S.; Suzuki, K.; Toda, K.; Seino, S. Mathematical Models and Experiments on the Acoustic Properties of Granular Packing Structures (Measurement of Tortuosity in Hexagonal Close-Packed and Face-Centered Cubic Lattices). Materials 2022, 15, 7393. https://doi.org/10.3390/ma15207393

AMA Style

Sakamoto S, Suzuki K, Toda K, Seino S. Mathematical Models and Experiments on the Acoustic Properties of Granular Packing Structures (Measurement of Tortuosity in Hexagonal Close-Packed and Face-Centered Cubic Lattices). Materials. 2022; 15(20):7393. https://doi.org/10.3390/ma15207393

Chicago/Turabian Style

Sakamoto, Shuichi, Kyosuke Suzuki, Kentaro Toda, and Shotaro Seino. 2022. "Mathematical Models and Experiments on the Acoustic Properties of Granular Packing Structures (Measurement of Tortuosity in Hexagonal Close-Packed and Face-Centered Cubic Lattices)" Materials 15, no. 20: 7393. https://doi.org/10.3390/ma15207393

APA Style

Sakamoto, S., Suzuki, K., Toda, K., & Seino, S. (2022). Mathematical Models and Experiments on the Acoustic Properties of Granular Packing Structures (Measurement of Tortuosity in Hexagonal Close-Packed and Face-Centered Cubic Lattices). Materials, 15(20), 7393. https://doi.org/10.3390/ma15207393

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