Research on the Sound Absorption Performance of Porous Asphalt Concrete with Different Air Voids Based on the Finite Element Models

: This paper aims to investigate the effect of the void structure of porous asphalt concrete (PAC) on the sound absorption performance. The sound absorption coefﬁcient (SAC) spectra of PAC with various voids were measured using the transfer function method and the air void structure was identiﬁed from Computed Tomography (CT) scanning images. The ﬁnite element model (FEM) of the void structure was established via CT cross-sectional images to research the noise reduction mechanism of PAC. At different frequencies, the sound pressure level distribution and SAC were calculated. The following are the main conclusions: (1) As the void content of PAC increases, the SAC spectrum curve becomes higher, more rightward and wider, while the peak values, their corresponding frequencies and mean values increase. (2) The SAC decreases with the average number of voids within each image and increases with the void equivalent diameter and fractal dimension; there exist linear relationships between these three parameters and the SAC. (3) The sound absorption capacity of the void with larger size and fractal dimension is higher than the smaller one and the voids connected to the upper and lower surfaces provide more noise reduction. (4) The FEMs of PAC can predict the SAC of PAC.


Introduction
Road traffic has been rapidly developed in the past decades.However, traffic noise pollution has endangered the physical and emotional health of drivers and inhabitants near roads.Long-term exposure to noise may even cause health issues including tinnitus and hypertension [1].
Common methods for reducing road noise include installing sound barriers, planting green belts, and constructing porous asphalt pavements.However, there are some shortcomings with the first and second methods, such as more time and money are needed.At the annual meetings of the Transportation Research Board (TRB), many opinions were shared to deal with the traffic noise problems [2][3][4].It is widely recognized that aiming at the PAC design is better than the research on tire structure because the former costs less and still has massive noise reduction.
Therefore, many road engineers have focused on porous asphalt pavements.Porous asphalt pavement is an environmentally friendly pavement.Due to the rich air voids, the drainage ability, frictional property and sound absorption ability of porous asphalt pavement are higher than those of asphalt concrete (AC) pavement.The void structure of the road surface can effectively reduce noise, which spreads from the interface of the road and the tire.Thus, it is significant to study the air void structure of porous asphalt concrete (PAC) [5].
During engineering, it is inconvenient to obtain the air void structure in PAC; therefore, most existing studies have focused on the effects of macroscopic indicators on the noise reduction capability of PAC.Dai and Lou [6] conducted near-field noise measurements on asphalt concrete pavement with a nominal maximum particle size of 13 mm (AC-13) and porous asphalt concrete pavements with a nominal maximum particle size of 13 mm (PAC-13), ascertaining the sound-decreasing level.Luong et al. [7] analyzed the influences of different contents of rubber additives on the sound absorption coefficient (SAC) of AC-13, whose results showed a positive relationship.By analyzing the mix grading, component percentage, aggregate shape, and other factors, Losa and Leandri [8] predicted the SAC of PAC-13 using the Neithalath model.Tang [9] verified that a double-layer porous asphalt pavement structure still has a high noise reduction performance under semirigid subgrade conditions.Although the test results of noise value could be used to compare the effect of reducing sound between PAC pavement and AC pavement, the noise absorption mechanism of PAC is difficult to be found by the macroscopic evaluation of noise reduction value.
The sound absorption capability of PAC is highly affected by the size and morphology of voids.Accordingly, it is important to accurately obtain the void structure, which supports further research on the noise reduction effect of PAC.Industrial computed tomographic (CT) scanning is the primary technology used to acquire a microscopic view of PAC structures with voids.Masad et al. [10] used a mathematical method of multiple regression to predict the air void distribution in PAC and estimate the water permeability by void morphology, which shows adequate accuracy.Pei et al. [11] investigated the fractional dimension of voids and void content, revealing a positive correlation between them.Jiang et al. [12] introduced image processing techniques to investigate the relationship between the macroscopic performance of PAC and microscopic features, which revealed that they were closely related.Gao et al. [13] imposed CT scanning images to determine the true depth of structural voids to adjust the sand spreading method for PAC.Wang et al. [14] adopted a 3D reconstruction technique to obtain void parameters, including void surface area, void volume and equivalent void diameter, by which a prediction model for SAC, which possesses high simulation accuracy, was created.Sun et al. [15][16][17][18] employed finite element software to examine the sound transmission loss (STL) of cross-sectional images; the result was very similar to the test values, demonstrating the correctness of the simulation.Some research established the simplified void models using acoustic-electric analogy theory or Helmholtz resonance theory [19,20], but it is essential to acquire the material's properties and void characteristics, which were difficult to obtain.The mentioned researches focus mainly on the void characteristic, but the effect of the void on sound absorption has been less concerned, which is the key to revealing the noise reduction mechanism of PAC.
In conclusion, previous research studies have provided most of the macroscopic indicators for the noise reduction design of PAC, and the relationship between void parameters and noise reduction impact has been investigated from a microscopic perspective.However, how to reveal the noise reduction mechanism of PAC with voids is still a challenge to be addressed.This research is aimed at solving the above problems.Two highlight points were proposed: one is the statistical models between void parameters of PAC and SAC and the other is the FEM established by CT images, which shows high accuracy in SAC prediction.Both two highlights are helpful to the further study of the sound absorption capability of PAC and useful to the practical design of low-noise pavement.
In this paper, the relationship between the air voids and the sound absorption capability of PAC was examined.First, the transfer function method was used to obtain the SAC spectrum of PAC with different void contents.Second, the void structure, identified from CT scanning images, was used to analyze the microscopic parameters.Third, statistical models between microscopic parameters and SAC were established.Last, a finite element model based on CT section images was built to predict the SAC spectrum of PAC and explore the noise reduction mechanism of porous asphalt pavement at the microscopic scale.

Materials
The technical indexes of the aggregates are shown in Table 1.To model high-quality PAC samples, styrene butadiene styrene (SBS)-modified asphalt combined with a highviscosity asphalt (HVA) modifier at a ratio of 92:8 (asphalt:modifier) was used in this study.Tables 2 and 3 show the results of the HVA, SBS-modified asphalt and composite-modified asphalt (both SBS and HVA were used to modify the asphalt) technical indicator tests.The technical indexes of the materials meet the requirement of Chinese standards [21,22].

Aggregate Gradation
The passage rates of the 2.36 mm sieve (P 2.36 ) and 4.75 mm sieve (P 4.75 ) have a significant impact on the void content of PAC-13 [23].The P 2.36 and P 4.75 are closely related to the maximum theoretical relative density by which the void content is affected.Therefore, four PAC-13 gradations with target void contents of 16%, 18%, 22% and 24% were designed by altering P 2.36 and P 4.75 , which were numbered A1, A2, A3 and A4, respectively.The gradations and the void content of samples (A1~A4) are shown in Table 4; the asphalt-stone ratio was calculated by using the leak test and the Cantabro test.

Test Scheme
The main raw materials of PAC models include aggregate and asphalt.In this test, basalt, basalt machine-made sand and limestone powder were used as the coarse aggregate, fine aggregate and mineral powder filler, respectively.The aggregate was sieved using a sieve shaker.After the cleaning and drying, the aggregate was set in a heating cabinet at 160 • C for 90 min.The SBS-modified asphalt was heated to a liquid state before the models were made.At the beginning of the test, the aggregate was poured into the concrete mixing machine singly.After 60 s, the weighed HVA modifier was put into the machine and stirred together with the aggregate.After 30 s, the quantitative SBS-modified asphalt was finally poured.After 90 s, the PAC was completed.At last, the concrete was taken from the machine to fill the rutting plate mold which was heated to 165 • C. The heating and mixing schemes were carried out above 160 • C.
To verify the basic performance of the PAC, the Marshall test pieces were made by the aggregate gradation of the A4 sample.The main tested results of Marshall stability, dynamic stability, fly-away loss and Marshall flow are 6.97 kN, 9719 times/mm, 14.1% and 27.8 mm.All the technical indexes of the Marshall sample satisfy the requirements of Chinese specifications [22].
The dimensions of all the rutting plate models were 300 mm × 300 mm × 40 mm.Each model was cut into 9 cuboids with the same size of 100 mm × 100 mm × 40 mm.The sample at the center of the rutting plate models after cutting was chosen as the study sample for further investigation in this study.The samples were also named A1~A4, the same as the corresponding gradation in Table 4.

Sound Absorption Coefficient Test
The SAC is frequently regarded as the sound absorption performance evaluation parameter.By using the transfer function method [24], the SAC spectra of samples A1~A4 were measured at different frequencies.Figure 1 shows the instrument for the SAC measurements.

Sound Absorption Coefficient Test
The SAC is frequently regarded as the sound absorption performance evaluation pa rameter.By using the transfer function method [24], the SAC spectra of samples A1~A were measured at different frequencies.Figure 1 shows the instrument for the SAC mea urements.

Result Analysis
Although the traffic noise spectrum has an extremely wide frequency distributio range (20 Hz to 20,000 Hz), the frequency band from 500 Hz to 1500 Hz is the noise fre quency range of traffic noise energy [25].As a result, the sound absorption performanc of samples in the frequency range of 500-1500 Hz was investigated in this study.Figure and Table 4 show the analysis results of the SAC test.The SAC results of samples A1~A4 are shown in Figure 2 and Table 5.The peak ab sorption coefficients of A1~A4 are 0.45, 0.59, 0.88 and 0.93, which are found at 725 Hz, 91 Hz, 1050 Hz and 1140 Hz, respectively.In most frequency ranges, the SAC of samples A and A4 exceeded 0.5, resulting in good noise reduction effects.It can be concluded tha the SAC of PAC gradually increases with increasing frequency, reaches a peak, and the decreases as the frequency increases.The peak SAC corresponds to the frequency ban mostly located in the range of 700 to 1200 Hz.According to research on the dominan traffic noise frequency, PAC can effectively reduce noise caused by small cars.As the voi

Result Analysis
Although the traffic noise spectrum has an extremely wide frequency distribution range (20 Hz to 20,000 Hz), the frequency band from 500 Hz to 1500 Hz is the noise frequency range of traffic noise energy [25].As a result, the sound absorption performance of samples in the frequency range of 500-1500 Hz was investigated in this study.Figure 2 and Table 4 show the analysis results of the SAC test.
corresponding gradation in Table 4.

Sound Absorption Coefficient Test
The SAC is frequently regarded as the sound absorption performance evaluation parameter.By using the transfer function method [24], the SAC spectra of samples A1~A4 were measured at different frequencies.Figure 1 shows the instrument for the SAC measurements.

Result Analysis
Although the traffic noise spectrum has an extremely wide frequency distribution range (20 Hz to 20,000 Hz), the frequency band from 500 Hz to 1500 Hz is the noise frequency range of traffic noise energy [25].As a result, the sound absorption performance of samples in the frequency range of 500-1500 Hz was investigated in this study.Figure 2 and Table 4 show the analysis results of the SAC test.The SAC results of samples A1~A4 are shown in Figure 2 and Table 5.The peak absorption coefficients of A1~A4 are 0.45, 0.59, 0.88 and 0.93, which are found at 725 Hz, 910 Hz, 1050 Hz and 1140 Hz, respectively.In most frequency ranges, the SAC of samples A3 and A4 exceeded 0.5, resulting in good noise reduction effects.It can be concluded that the SAC of PAC gradually increases with increasing frequency, reaches a peak, and then decreases as the frequency increases.The peak SAC corresponds to the frequency band mostly located in the range of 700 to 1200 Hz.According to research on the dominant traffic noise frequency, PAC can effectively reduce noise caused by small cars.As the void The SAC results of samples A1~A4 are shown in Figure 2 and Table 5.The peak absorption coefficients of A1~A4 are 0.45, 0.59, 0.88 and 0.93, which are found at 725 Hz, 910 Hz, 1050 Hz and 1140 Hz, respectively.In most frequency ranges, the SAC of samples A3 and A4 exceeded 0.5, resulting in good noise reduction effects.It can be concluded that the SAC of PAC gradually increases with increasing frequency, reaches a peak, and then decreases as the frequency increases.The peak SAC corresponds to the frequency band mostly located in the range of 700 to 1200 Hz.According to research on the dominant traffic noise frequency, PAC can effectively reduce noise caused by small cars.As the void content increases, the band of the SAC spectrum widens, the peak and mean values of SAC increase and the frequency corresponding to the peak increases, which causes the overall spectrum curves to show a trend of moving up and to the right.This means that the overall noise reduction effect improves.The characteristics of voids are critical for understanding the PAC noise reduction mechanism.In this study, the CT images of four samples were obtained using industrial CT scanning technology.The CT scanning equipment was set parallel to the height of the samples, as shown in Figure 3.The scan spacing was set to 0.80 mm in the vertical X and vertical Y directions and the images were obtained after rejecting the ends.Figure 4a shows a CT image of sample A2.
content increases, the band of the SAC spectrum widens, the peak and mean values of SAC increase and the frequency corresponding to the peak increases, which causes the overall spectrum curves to show a trend of moving up and to the right.This means that the overall noise reduction effect improves.

Computed Tomography Image Segmentation
The characteristics of voids are critical for understanding the PAC noise reduction mechanism.In this study, the CT images of four samples were obtained using industrial CT scanning technology.The CT scanning equipment was set parallel to the height of the samples, as shown in Figure 3.The scan spacing was set to 0.80 mm in the vertical X and vertical Y directions and the images were obtained after rejecting the ends.Figure 4a shows a CT image of sample A2.The Otsu method is a simple and efficient image binarization algorithm without losing precision [26].This method was utilized to segment the CT images and identify void structures in the images.Using the CT image in Figure 4a as an example, the processing result is displayed in Figure 4b, demonstrating that the concrete and voids can be correctly and fully separated.Among all the analysis results of the CT images, the aggregates were not misdiagnosed as voids.Therefore, the quality of this void identification can satisfy the requirements of this investigation.content increases, the band of the SAC spectrum widens, the peak and mean values of SAC increase and the frequency corresponding to the peak increases, which causes the overall spectrum curves to show a trend of moving up and to the right.This means that the overall noise reduction effect improves.

Computed Tomography Image Segmentation
The characteristics of voids are critical for understanding the PAC noise reduction mechanism.In this study, the CT images of four samples were obtained using industrial CT scanning technology.The CT scanning equipment was set parallel to the height of the samples, as shown in Figure 3.The scan spacing was set to 0.80 mm in the vertical X and vertical Y directions and the images were obtained after rejecting the ends.Figure 4a shows a CT image of sample A2.The Otsu method is a simple and efficient image binarization algorithm without losing precision [26].This method was utilized to segment the CT images and identify void structures in the images.Using the CT image in Figure 4a as an example, the processing result is displayed in Figure 4b, demonstrating that the concrete and voids can be correctly and fully separated.Among all the analysis results of the CT images, the aggregates were not misdiagnosed as voids.Therefore, the quality of this void identification can satisfy the requirements of this investigation.The Otsu method is a simple and efficient image binarization algorithm without losing precision [26].This method was utilized to segment the CT images and identify void structures in the images.Using the CT image in Figure 4a as an example, the processing result is displayed in Figure 4b, demonstrating that the concrete and voids can be correctly and fully separated.Among all the analysis results of the CT images, the aggregates were not misdiagnosed as voids.Therefore, the quality of this void identification can satisfy the requirements of this investigation.

Void Characteristics Analysis
The CT images were processed to identify the void structures in the samples [27][28][29][30].In this paper, the analysis of void characteristics was conducted using CT images that are vertical to the X-axis and the Y-axis.The characteristic parameters of the void were identified by writing programs in MATLAB, including N as the average number of voids within each image, a as the average area of individual voids and d as the equivalent diameter of voids.The calculations are as follows: where N i is the total number of voids in the sample's i-th layer cross-sectional image; n is the total number of layers in the CT scan; A i is the total area of voids in the sample's ith layer cross-sectional image; A is the average area of voids in all cross-sectional pictures of each sample and l is the average length of each void.
To compute the void content of CT scan identification, the total area of voids within each cross-sectional picture was transformed into the void content of each layer, and the average value was determined.
The aforementioned formulas only describe the size and number of voids, not their morphological properties.Because the internal voids of PAC have obvious fractal properties [31], in this paper, the fractal dimension was used to evaluate the internal complexity of voids, with the box-counting method in fractal theory as the quantification criterion.The CT images were divided by a square grid with a side length δ; the number of boxes N(δ) occupied by the image boundary with the corresponding side length δ was counted separately.The side lengths of the square grid were then changed, which generates different side lengths with corresponding grid numbers.The side lengths δ and the grid numbers N(δ) were linearly regressed in logarithmic coordinates, with the slope D of the regression line representing the fractal dimension measured using the box-counting method.The formulations are as follows: taking the logarithm of both sides of Equation ( 14) simultaneously, that is, Equation (15).
where D i is the fractal dimension within the i-th layer of the section images of the sample and D B is the average fractal dimension within all section images of the sample.Table 6 shows the microscopic void characteristics of samples A1~A4.The void content obtained by image recognition technology is very similar to the void content measured by the volume method.The recognition void content is slightly smaller than that obtained by actual measurement, which is due to the limitations of CT scanning accuracy and algorithm recognition accuracy.There are micro-voids that are difficult to identify, which results in a smaller identification void content.For samples A1~A4, Figure 5 shows the relationship between the void content determined and the void structural parameters, including the equivalent diameter and the average number of voids within each image.The measured void content is represented by the horizontal coordinate x, while the double vertical coordinates y 1 and y 2 represent the void equivalent diameter and the average number of voids within each image, respectively.For y 1 and y 2 , linear fits are performed and the resulting relational equation and correlation coefficient R 2 are shown in Figure 5.The number of voids inside the PAC decreases as the void content increases but the void equivalent diameter continuously increases, which reveals that the void content of PAC increases not because the number of voids increases but because the void size increases.For samples A1~A4, Figure 5 shows the relationship between the void content determined and the void structural parameters, including the equivalent diameter and the average number of voids within each image.The measured void content is represented by the horizontal coordinate , while the double vertical coordinates  1 and  2 represent the void equivalent diameter and the average number of voids within each image, respectively.For  1 and  2 , linear fits are performed and the resulting relational equation and correlation coefficient  2 are shown in Figure 5.The number of voids inside the PAC decreases as the void content increases but the void equivalent diameter continuously increases, which reveals that the void content of PAC increases not because the number of voids increases but because the void size increases.The relationships between the above parameters and the peak and mean values of the SAC (named  in the figures) of the PAC are shown in Figure 6a,b.The results follow the same pattern as the improvement in noise reduction performance caused by increasing the macroscopic void content.Although an increase in void content reduces the number of voids, the fitted results show that the sensitivity of the SAC with the equivalent diameter is greater than that of the number of voids.This indicates that while the number of voids decreases, the equivalent diameter demonstrates a greater increase in the noise reduction effect, which proves that an increase in void diameter can increase the incident acoustic waves and provide a larger internal surface area for sound energy reflection and frictional attenuation.When the two void parameters change, the peak value of the SAC changes more than the mean value.It can be seen that as long as the PAC's design void content requirements (16-25%) are met, PAC may provide good noise reduction in most frequency bands.The relationships between the above parameters and the peak and mean values of the SAC (named α in the figures) of the PAC are shown in Figure 6a,b.The results follow the same pattern as the improvement in noise reduction performance caused by increasing the macroscopic void content.Although an increase in void content reduces the number of voids, the fitted results show that the sensitivity of the SAC with the equivalent diameter is greater than that of the number of voids.This indicates that while the number of voids decreases, the equivalent diameter demonstrates a greater increase in the noise reduction effect, which proves that an increase in void diameter can increase the incident acoustic waves and provide a larger internal surface area for sound energy reflection and frictional attenuation.When the two void parameters change, the peak value of the SAC changes more than the mean value.It can be seen that as long as the PAC's design void content requirements (16-25%) are met, PAC may provide good noise reduction in most frequency bands.

Sample
The geometric complexity of the void structure is described by the fractal dimension.Figure 7 shows the relationship between the void content and the average fractal dimension of all images of each sample.As the void content increases, D B increases, showing that as the complexity of the void morphology improves, the structure of the internal voids of the PAC becomes more complicated and the internal roughness rises.
As shown in Figure 8, the noise reduction impact improves as the fractal dimension grows, which results from a more complicated void structure, leading to a greater chance of acoustic wave reflections when the waves enter the interior void structure.Due to the increase in the fractal dimension, the interface between the acoustic waves and the void surface is increased, providing a larger friction area.This causes greater viscous loss generated between the air and the void wall, resulting in a higher SAC.
The fractal dimension of the voids can explain the variation in SAC with d and N.As the void content increases, N decreases, which seems to reduce the space where sound energy is lost.However, the increase in d increases the sound wave incidence and D B increases with increasing void content, which enhances the amount of sound energy loss inside the structure; the overall sound absorption and noise reduction performance of the structure is improved.
The three characteristics, d, N and D B , are strongly correlated with the SAC.These characteristics may be used to design PAC pavement sound absorption performance in the future.The geometric complexity of the void structure is described by the fractal dimension.Figure 7 shows the relationship between the void content and the average fractal dimension of all images of each sample.As the void content increases,   increases, showing that as the complexity of the void morphology improves, the structure of the internal voids of the PAC becomes more complicated and the internal roughness rises.As shown in Figure 8, the noise reduction impact improves as the fractal dimension grows, which results from a more complicated void structure, leading to a greater chance of acoustic wave reflections when the waves enter the interior void structure.Due to the increase in the fractal dimension, the interface between the acoustic waves and the void surface is increased, providing a larger friction area.This causes greater viscous loss generated between the air and the void wall, resulting in a higher SAC.The geometric complexity of the void structure is described by the fractal dimension.Figure 7 shows the relationship between the void content and the average fractal dimension of all images of each sample.As the void content increases,   increases, showing that as the complexity of the void morphology improves, the structure of the internal voids of the PAC becomes more complicated and the internal roughness rises.As shown in Figure 8, the noise reduction impact improves as the fractal dimension grows, which results from a more complicated void structure, leading to a greater chance of acoustic wave reflections when the waves enter the interior void structure.Due to the increase in the fractal dimension, the interface between the acoustic waves and the void surface is increased, providing a larger friction area.This causes greater viscous loss generated between the air and the void wall, resulting in a higher SAC.The fractal dimension of the voids can explain the variation in SAC with  and .As the void content increases,  decreases, which seems to reduce the space where sound energy is lost.However, the increase in  increases the sound wave incidence and   increases with increasing void content, which enhances the amount of sound energy loss

Theory of Finite Element Models
The internal interconnected voids of PAC and the damping of the material itself perform the functions of sound absorption and noise reduction.As a result, not only must the continuity equation for fluid momentum conservation (named the Navier-Stokes equation, N-S equation) be considered in this problem but also the structural dynamics equations.The following equations are available for a compressible Newtonian fluid [32]: where v is the speed of sound; p is the sound pressure; ρ is the air density and µ is the aerodynamic viscosity.From left to right, each part of the formula represents the inertial force, acoustic pressure, viscous force and external force.
Equations ( 17) and ( 18) should be solved simultaneously; they represent momentum conservation and mass conservation, respectively.By solving these two equations under specified boundary conditions, the fluid flow rate and pressure inside a given structure can be obtained.The analytical solution is often difficult to obtain accurately due to the complexity of the equations and the structure itself.As a result, solving the equations necessitates the simplification of the assumed conditions and the structure with the application of numerical methods [33].
When air is assumed to be a non-viscous compressible and homogeneous fluid, the N-S equation can be simplified to Equation (19), which means the acoustic fluctuation equation [34].
In the finite element method, the discretized sound field fluctuation equations and structural dynamics equations are Equations ( 11) and (12).
The acoustic pressure term acting at the air-pavement interface is added to the righthand side of the structural dynamic Equation ( 12) to couple the fluid to the solid and obtain Equation ( 13): Then, Equation ( 13) can be simplified into equation (15).
[M e ] Combining Equation (11) with Equation ( 15), the acoustic-solid coupling Equation ( 16) can be obtained as follows: is the fluid stiffness matrix; [M e ] is the mass matrix of the solid structure; [C e ] is the damping matrix of the solid structure; [K e ] is the stiffness matrix of the solid structure; {U e } is the solid displacement matrix; {p e } is the fluid displacement matrix and {F e } is the exogenic force.

Construction of the Finite Element Model
Sound absorption in porous asphalt pavements is a common acoustic-solid coupling problem that is influenced by a variety of factors, such as air temperature, humidity and vehicle loading, in real-world situations.The majority of existing research has focused on the fields of pressure acoustics or fluctuating acoustics using air as the analysis medium [35][36][37].The damping of pavement structures also has an impact on sound absorption.The sound pressure acts on the structure, causing structural vibrations.At the same time, the structural vibrations cause the surrounding air to vibrate, causing sound waves to radiate outward.The mass coupling matrix and the stiffness coupling matrix are used in Equation ( 16) to represent this pair of interactions.
Numerical solutions were obtained by using the finite element method in this paper.COMSOL Multiphysics is a finite element simulation program that can handle a coupled multi-physics field.This software was used in this study to couple pressure acoustics and solid mechanics to obtain an acoustic-structural boundary.In the simulated results, including the calculated sound pressure and sound pressure level distribution, the noise reduction effect of PAC was analyzed under different void contents, which explains the noise reduction mechanism of PAC from a microscopic perspective.
The void structure of PAC is complicated and the distribution of void locations is extremely random.To obtain an accurate void structure, the processed CT images of the PAC samples were used as the import model.For each sample, as the image processing and void identification operations described in Chapter 3.1, the bitmap images were obtained.The bitmap images in this paper were programmed to be converted into vector images in AutoCAD format, which were used as the import models by the COMSOL CAD interface.The result is shown in Figure 9.The media involved in the simulation are air and PAC.Based on existing studies [38][39][40], the required material parameters are summarized, as shown in Table 7.The media involved in the simulation are air and PAC.Based on existing studies [38][39][40], the required material parameters are summarized, as shown in Table 7.
Because material parameters change with external temperature, humidity and other factors, the following assumptions were made to clarify the study's scope and accurately describe the noise absorption of PAC.

1.
The acoustic calculations in this paper are solved under linear acoustic conditions.The acoustic pressure is very small compared to the atmospheric pressure because it is the residual pressure after the atmospheric pressure perturbation.Therefore, in this paper, a compressional wave under a small acoustic perturbation is considered.

2.
All the air particles conduct simple harmonic vibration at the acoustic equilibrium position.In this paper, it is assumed that no air flow occurs in the air flow field, i.e., no wind condition.Therefore, the N-S equation can be simplified to a linear N-S equation.

3.
The void wall surface is completely filled with air.In the model, the air is in continuous and complete contact with the void structure.4.
The adiabatic condition is applied to the air model boundary.The thermal energy generated by friction between the air and the void wall raises the temperature of the PAC but it does not exchange heat with the air.Only the sound absorption effect at 293.15 K is taken into account.5.
The incident direction of the sound source is selected as the Z direction.The type and size of the sound source have no effect on the final calculation results because the material's SAC is its own property.In this paper, the incident source is a plane wave perpendicular to the surface of the sample model with a pressure amplitude of 1 Pa.The model of the measuring instrument was established and the transfer function method was used to obtain the SAC, in which the sizes of the instrument, the position and the spacing of the microphone were designed according to the actual instrument to visualize the difference between the simulation model and the measured sound absorption effect.The model is shown in Figure 10.The working principle is that the two microphones measure the incident and reflected sound pressures, which are used to calculate the transfer function.The SAC is calculated from the incident sound reflection factor determined by the transfer function.The calculation process is shown as follows.
If the side of the sample in contact with the acoustic wave is the reference surface ( = 0), then the distance of microphone 1 from the reference surface is  1 and the distance of microphone 2 from the reference surface is  2 , as the Figure .10 shown.The incident sound pressure   and the reflected sound pressure   can be expressed as: =   ̂− 0  (18) where   ̂ is the amplitude of   on the reference plane;   ̂ is the amplitude of   on the reference plane;  0 =  0 ′ −  0 ′′ is the complex wavenumber.The sound pressures  1 and  2 at the two microphone positions are: The working principle is that the two microphones measure the incident and reflected sound pressures, which are used to calculate the transfer function.The SAC is calculated from the incident sound reflection factor determined by the transfer function.The calculation process is shown as follows.
If the side of the sample in contact with the acoustic wave is the reference surface (x = 0), then the distance of microphone 1 from the reference surface is x 1 and the distance of microphone 2 from the reference surface is x 2 , as the Figure .10 shown.The incident sound pressure p i and the reflected sound pressure p R can be expressed as: p R = pR e −jk 0 x (18 where pi is the amplitude of P i on the reference plane; pR is the amplitude of P R on the reference plane; k 0 = k 0 − jk 0 is the complex wavenumber.
The sound pressures p 1 and p 2 at the two microphone positions are: p 2 = pi e jk 0 x 2 + pR e −jk 0 x 2 (20) If the distance between microphones 1 and 2 is x 1 − x 2 = s, then the incident and reflected wave transfer functions H 1 and H R are: If the incident acoustic reflection factor is r, then pR = r pi , and the total sound field transfer function H 12 can be derived by combining Equations ( 19) and (20) as follows: Combining Equation (21) with Equation ( 22), we can obtain r as Equation ( 24): Therefore, the SAC α can be calculated by r, as shown in Equation ( 25): where r r is the real part of the incident acoustic reflection factor and r i is the imaginary part of the incident acoustic reflection factor.The pressure acoustic module and the structural mechanics module were added to the model after it was established; the acoustic-solid coupling interface for both is generated by the COMSOL automatically, after which the mesh can be divided.Notably, if the size of the divided cells is too large, then the final model calculation results will be distorted and will have difficulty converging [24].An acoustic medium cell was used in this study, whose size is smaller than 1/6 of the wavelength corresponding to the calculated maximum frequency (1500 Hz). Figure 11 shows the grid created by dividing the imported CT image in Figure 9, in which the blue part represents the air and the gray part represents the PAC model.generated by the COMSOL automatically, after which the mesh can be divided.Notably, if the size of the divided cells is too large, then the final model calculation results will be distorted and will have difficulty converging [24].An acoustic medium cell was used in this study, whose size is smaller than 1/6 of the wavelength corresponding to the calculated maximum frequency (1500 Hz). Figure 11 shows the grid created by dividing the imported CT image in Figure 9, in which the blue part represents the air and the gray part represents the PAC model.

Analysis of the Simulation Results
In this analysis, the model was swept in steps of 25 Hz in the frequency range of 500-1500 Hz after meshing was completed.As calculated from parts of CT images of A1~A4, Figure 12 shows the distribution of the sound pressure levels for four different void content PAC structures at different frequencies.

Analysis of the Simulation Results
In this analysis, the model was swept in steps of 25 Hz in the frequency range of 500-1500 Hz after meshing was completed.As calculated from parts of CT images of A1~A4, Figure 12 shows the distribution of the sound pressure levels for four different void content PAC structures at different frequencies.
Figure 12 shows that PAC-13 has a good sound absorption effect and the sound pressure level has a good grading level.The sound absorption performance of PAC first increases and then decreases with increasing frequency, which is consistent with the results of the macroscopic tests.In the light of the perforated plate principle [31] and the Helmholtz resonator principle [32], the void structures can be regarded as sound absorbers.When the frequency is close to the resonant frequency, the SAC of PAC is larger than the frequency deviating from the resonant frequency.Because the noise-reducing capability of small voids is limited, noise reduction is primarily provided by the large voids.When the void content of PAC increases, the average internal void size and void depth of PAC increase, which improves the sound absorption effect of the whole structure by increasing the space of acoustic vibration.Due to the larger distribution range of void size, there is an obvious noise reduction effect in the middle and high frequency bands.
Figure 12 shows that the voids connecting the surface and the surrounding concrete have superior sound absorption levels and sound waves are more likely to be reflected when transmission occurs in concrete near the interior closed voids.The surface-connected voids are the key structures that provide the PAC noise reduction effect, as seen by these figures.The capacity of PAC to reduce noise is also influenced by the internal complexity of the voids.The void morphology grows increasingly complex as the void content increases, enhancing the overall sound absorption capability.However, the complexity of voids is constrained by their pore size and length.With a small void content, the void structure has a large fractal dimension; however, the small size is more likely to lead to acoustic energy concentration, which is not conducive to energy dissipation.With a large void content, the internal structure of voids being very rough, sound energy can be significantly reduced.
In summary, as the void content rises, the internal voids become more easily joined to create a whole, resulting in noise reduction.The complexity of the void has a favorable effect on the sound absorption performance of PAC.The surface-connected voids are an important structure with respect to the sound absorption performance.However, the void size and the complexity of the voids often affect each other and work together to influence the sound absorption ability of PAC.In a word, as the void size increases, the sound absorption performance of PAC improves.
After importing all models, the COMSOL built-in commands, real (acpr_p) and imag (acpr_p), were used to obtain the real part and the imaginary part sound pressure at the positions of microphones 1 and 2; the SAC was calculated using Equations ( 1)-(9) and the SAC spectrum of each image was obtained.Each sample can export spectra of 50 images, which were averaged to obtain the simulated SAC spectrum, as shown in Figure 13.The simulation results show the same trend as the measured results, which increase to peak, then consequently decrease.The peak SACs of the simulated A1~A4 samples are 0.42, 0.51, 0.57 and 0.69 at 760 Hz, 975 Hz, 1115 Hz and 1225 Hz, respectively.The structural characteristic frequencies correspond not only to the peak but also to the trough.The overall resonance of the structure at certain frequencies causes a decrease in the absorption coefficient, which is the trough in the spectrum, as shown in Figure 13b,d for the SAC trough at approximately 700 Hz.The simulated results show the same trend as the test values in the 500-1500 Hz range.
In order to analyze the errors between the test and simulation results, maximum errors and average errors were calculated using Equations ( 26)-( 28): where δ i is the absolute value of the error between test result and simulation result in a single frequency; α 1 is the test SAC of PAC; α 2 is the simulation SAC of PAC; δ is the average error of all frequencies; n is the total number of each frequency, the taken values space is 5 Hz in this paper and δ max is the maximum error at all frequencies.Furthermore, the measured results were compared to the mean values of the SAC obtained from the simulations.The measured and simulated SACs of A1 to A4 are shown in Figure 14.Both the tested and simulated results have a strong linear relationship with the void content, with correlation coefficients of 0.9368 and 0.9074, respectively.The simulated values are smaller than the measured values, with the simulated mean values of the absorption coefficients of samples A1 to A4 being 74.36%, 74.47%, 76.77% and 75.01% of the test values, respectively.Combined with Figures 12 and 13, the simulation results are more similar to the experimental results at low frequencies but less so at high frequencies.One reason is that the effect of inner roughness on the loss of acoustic waves grows as the frequency increases.The acquired void structure in this study after the CT scanning, image recognition and vectorization conversion causes the recognition accuracy to be limited, resulting in smaller simulation results than test results at high frequencies.Another reason is that the scanning interval between each cross-sectional image will truncate the surfaceconnected pores, which are cut into some surface-connected voids and inner voids in each image.The damping qualities of the material and the structure of the voids are primarily responsible for the noise reduction effect of PAC.The material properties can be simulated fully during the simulation process.However, the simulation of voids is limited, which reduces the simulation accuracy.The goal of follow-up research is to determine more about 3D void modeling.frequencies.One reason is that the effect of inner roughness on the loss of acoustic waves grows as the frequency increases.The acquired void structure in this study after the CT scanning, image recognition and vectorization conversion causes the recognition accuracy to be limited, resulting in smaller simulation results than test results at high frequencies.
Another reason is that the scanning interval between each cross-sectional image will truncate the surface-connected pores, which are cut into some surface-connected voids and inner voids in each image.The damping qualities of the material and the structure of the voids are primarily responsible for the noise reduction effect of PAC.The material properties can be simulated fully during the simulation process.However, the simulation of voids is limited, which reduces the simulation accuracy.The goal of follow-up research is to determine more about 3D void modeling.The SAC and the structural characteristic frequencies of the test results are similar to those of the simulated results.All the simulated average values of the SAC account for approximately 75% of the measured values.By the comparison between test and simulated results, the reasonableness of the model is proved, which shows that the PAC acoustic model that was created using the CT scanning methodology and the finite element approach is practical.

Conclusions
In this paper, CT scanning technology was used to obtain slice images of the PAC samples.Programming was used to identify the voids of the images.The relationship between the characteristic parameters of the voids and the SAC was investigated.The processed CT images were imported into the finite element software to build the acoustic models.Based on the model, the noise reduction mechanism of the void structure was analyzed.The main conclusions are as follows: 1.The SAC spectrum of PAC tends to increase to the peak value and then decrease.As the void content increases, the SAC spectrum curve shows a trend that becomes higher, more rightward and wider, while the peak values and their corresponding frequencies increase.The overall noise reduction effect of PAC is improved; 2. The SAC decreases with an increasing average number of voids and increases with an increasing equivalent diameter of voids and average fractal dimension; 3. The finite element models are built using CT images.At different frequencies, the simulated distribution of the sound pressure level is calculated, which can explain the noise reduction mechanism of PAC; The SAC and the structural characteristic frequencies of the test results are similar to those of the simulated results.All the simulated average values of the SAC account for approximately 75% of the measured values.By the comparison between test and simulated results, the reasonableness of the model is proved, which shows that the PAC acoustic model that was created using the CT scanning methodology and the finite element approach is practical.

Conclusions
In this paper, CT scanning technology was used to obtain slice images of the PAC samples.Programming was used to identify the voids of the images.The relationship between the characteristic parameters of the voids and the SAC was investigated.The processed CT images were imported into the finite element software to build the acoustic models.Based on the model, the noise reduction mechanism of the void structure was analyzed.The main conclusions are as follows: 1.
The SAC spectrum of PAC tends to increase to the peak value and then decrease.
As the void content increases, the SAC spectrum curve shows a trend that becomes higher, more rightward and wider, while the peak values and their corresponding frequencies increase.The overall noise reduction effect of PAC is improved; 2.
The SAC decreases with an increasing average number of voids and increases with an increasing equivalent diameter of voids and average fractal dimension; 3.
The finite element models are built using CT images.At different frequencies, the simulated distribution of the sound pressure level is calculated, which can explain the noise reduction mechanism of PAC; 4.
The finite element models can predict the SAC of PAC.At low frequencies, the results are accurate; at high frequencies, the results are smaller than the test results but still acceptable.

Figure 3 .Figure 4 .
Figure 3.This is a figure.Schemes follow the same formatting.

Figure 3 .
Figure 3.This is a figure.Schemes follow the same formatting.

Figure 3 .Figure 4 .
Figure 3.This is a figure.Schemes follow the same formatting.

Figure 5 .
Figure 5. Relationship between ,  and the void content.

Figure 5 .
Figure 5. Relationship between d, N and the void content.

Figure 6 .
Figure 6.The relationship between the SAC and the parameters of size and number: (a) Relationship between SAC and ; (b) Relationship between SAC and .

Figure 7 .
Figure 7. Relationship between   and the void content.

Figure 6 .Figure 6 .
Figure 6.The relationship between the SAC and the parameters of size and number: (a) Relationship between SAC and N; (b) Relationship between SAC and d.

Figure 7 .
Figure 7. Relationship between   and the void content.

Figure 7 . 21 Figure 8 .
Figure 7. Relationship between D B and the void content.Appl.Sci.2022, 12, x FOR PEER REVIEW 10 of 21

Figure 8 .
Figure 8. Relationship between the SAC and D B .
U e + [C e ] .U e + [K e ]{U e } = {F e } + s = ρ 0 [R e ] T is the mass coupling matrix; K f s = −[R e ] is the stiffness coupling matrix; M p e is the fluid mass matrix; K p e

21 Figure 10 .
Figure 10.Model of the impedance tube. 1 : the distance of microphone 1 from the reference surface,  2 : the distance of microphone 2 from the reference surface.

Figure 10 .
Figure 10.Model of the impedance tube.x 1 : the distance of microphone 1 from the reference surface, x 2 : the distance of microphone 2 from the reference surface.

Figure 12 .
Figure 12.Sound pressure level distribution of PAC with different void contents at different frequencies: (a) sample A1; (b) sample A2; (c) sample A3 and (d) sample A4.

Figure 12 .
Figure 12.Sound pressure level distribution of PAC with different void contents at different frequencies: (a) sample A1; (b) sample A2; (c) sample A3 and (d) sample A4.

Table 1 .
Technical indexes of the aggregate.

Table 2 .
Technical indexes of the HVA high-viscosity modifier.

Table 3 .
Technical indexes of the modified asphalt.

Table 5 .
Statistical results of the sound absorption coefficient α.
4. Study of the Microscopic Void Characteristics of PAC 4.1.Computed Tomography Image Segmentation

Table 5 .
Statistical results of the sound absorption coefficient α.

Table 5 .
Statistical results of the sound absorption coefficient α.

Table 6 .
Void characteristics of PAC.

Table 7 .
Material parameters involved in the simulation.
Because material parameters change with external temperature, humidity and other

Table 7 .
Material parameters involved in the simulation.