Study on the Influence of Clogging on the Cooling Performance of Permeable Pavement

Permeable pavement is often known as “cool road”. However, the cooling performance will be weakened due to clogging. In this paper, the temperature field distribution model of asphalt pavement was obtained by Green’s function. Gradations of porous asphalt mixture were designed to obtain different porosities, and the thermal properties of specimens with different porosities were tested and calculated. The simulation test was carried out to obtain the heating curves, which were used to verify the accuracy of the temperature model by comparing the results of the theoretical calculation to results of the test. The daily solar radiation intensity and air temperature changing functions were plugged into the model to calculate the temperature at the bottom of the middle surface. In this way, the simplified model of void fraction and cooling performance of the porous asphalt pavement was obtained. The results showed that the temperature at the bottom of the middle surface for permeable pavement was lower than that for traditional asphalt pavements. The gap was between 0.29 to 2.75 ◦C and it increased as the porosity of permeable pavement increased.


Introduction
As an important part of a "sponge city", permeable pavement plays a significant role in regulating local climate, accelerating groundwater recycling and reducing runoff on road surface.Permeable pavement is also known as "cool road".Due to its porous structure, the maximum temperature in the permeable pavement is lower than that in traditional asphalt pavement in summer [1].The maximum temperature of asphalt pavement affects the service and safety performance.Problems such as rutting and deflection are more likely to occur as the maximum temperature rises, leading to instability failures [2][3][4][5].The difference between the maximum temperatures of the permeable asphalt pavement and traditional asphalt pavement is defined as the cooling performance of permeable asphalt pavement.In the long-term use of permeable pavement, the study found that it was easily affected by the environment and the compaction conditions, leading to void clogging [6][7][8], weakening its permeability characteristics, and also affecting the cooling performance.
The numerical analysis method is often used to predict the temperature distribution internal to the pavement [9][10][11][12].This method is based on finite element or finite difference method.In these numerical models, pavement structures are discretized into elements for solving governing equations.However, the accuracy of the calculation results is highly dependent on the density of the mesh and total number of elements in the model.
Different from the numerical analysis method, the analytical method can be used to obtain the rigorous solutions of pavement temperature fields with appropriate assumptions.Additionally, various transfer methods were used in the analytical method.With this approach, each layer of the pavement was treated as isotropic [13].Based on the Fourier-Biot heat conduction equation and transform methods such as Laplace transform method and Green Function method, an analytical solution to heat conduction can be obtained very fast by solving the initial boundary value problems.The solution is usually an equation of instant temperature, thermal properties of materials in each layer and the heat flux on the surface of pavement [14][15][16].Comparing with the numerical analysis method, the analytical method is more convenient and accurate.So the analytical method was adopted in the paper.
According to the analytical method, the changing rule between thermal properties and porosity of the asphalt mixture should be studied first.Then the maximum internal temperatures of the pavement can be obtained, based on which the changing rule of cooling performance with void clogging can also be obtained.

Objective
This paper aims to establish a prediction model of the cooling performance decline for layered permeable pavements due to clogging problems.In order to achieve the objective, both theoretical analysis and laboratory test were performed.The Green's function method was used in the prediction of pavement temperature fields.The thermal properties of porous asphalt mixture with different porosities were tested in laboratory.Finally, the cooling performance was represented by the difference between the maximum temperatures in the permeable asphalt pavement and traditional pavement.

Heat Conduction Model of Multilayer Pavement Structure
Because the horizontal dimension of each layer of the pavement is much larger than the vertical direction, the heat conduction problem of the pavement can be regarded as the heat conduction problem of the large flat plate.Therefore, in the prediction of the temperature field, the method of one-dimensional heat conduction was adopted [17][18][19].The pavement structure is shown in Figure 1.At the surface of the pavement z = z 0 = 0.The surface is considered as a mixed boundary where the pavement receives the solar radiation, radiates part of the heat into the atmosphere and exchanges heat with the air.When the depth is infinitely deep, the temperature can be considered as a constant.
Water 2018, 10, x FOR PEER REVIEW 2 of 16 as Laplace transform method and Green Function method, an analytical solution to heat conduction can be obtained very fast by solving the initial boundary value problems.The solution is usually an equation of instant temperature, thermal properties of materials in each layer and the heat flux on the surface of pavement [14][15][16].Comparing with the numerical analysis method, the analytical method is more convenient and accurate.So the analytical method was adopted in the paper.
According to the analytical method, the changing rule between thermal properties and porosity of the asphalt mixture should be studied first.Then the maximum internal temperatures of the pavement can be obtained, based on which the changing rule of cooling performance with void clogging can also be obtained.

Objective
This paper aims to establish a prediction model of the cooling performance decline for layered permeable pavements due to clogging problems.In order to achieve the objective, both theoretical analysis and laboratory test were performed.The Green's function method was used in the prediction of pavement temperature fields.The thermal properties of porous asphalt mixture with different porosities were tested in laboratory.Finally, the cooling performance was represented by the difference between the maximum temperatures in the permeable asphalt pavement and traditional pavement.

Heat Conduction Model of Multilayer Pavement Structure
Because the horizontal dimension of each layer of the pavement is much larger than the vertical direction, the heat conduction problem of the pavement can be regarded as the heat conduction problem of the large flat plate.Therefore, in the prediction of the temperature field, the method of one-dimensional heat conduction was adopted [17][18][19].The pavement structure is shown in Figure 1.At the surface of the pavement z = z0 = 0.The surface is considered as a mixed boundary where the pavement receives the solar radiation, radiates part of the heat into the atmosphere and exchanges heat with the air.When the depth is infinitely deep, the temperature can be considered as a constant.According to the theory of heat transfer, heat conduction and heat convection that occurred due to direct contact between fluid and solid are called convention heat transfer, which can be calculated according to Newton cooling formula: where qr is the heat flux when convective heat transfer occurs, W•m −2 •s −1 ; hr is convective heat transfer coefficient, W/(m 2 •K).Tsurf is the temperature of the surface, K; and Tair is the temperature of the air, K.
The energy entering the road consists of short wave radiation and long wave radiation.Meanwhile, the road will radiate outward in a long wave way.The total amount of long wave radiation energy can be expressed as the following formula: According to the theory of heat transfer, heat conduction and heat convection that occurred due to direct contact between fluid and solid are called convention heat transfer, which can be calculated according to Newton cooling formula: where q r is the heat flux when convective heat transfer occurs, W•m −2 •s −1 ; h r is convective heat transfer coefficient, W/(m 2 •K).T surf is the temperature of the surface, K; and T air is the temperature of the air, K.

of 16
The energy entering the road consists of short wave radiation and long wave radiation.Meanwhile, the road will radiate outward in a long wave way.The total amount of long wave radiation energy can be expressed as the following formula: where E L is the amount of long wave radiation energy, W•m −2 •s −1 ; σ sb is Stefan-Boltzmann constant, 5.67 × 10 −8 W•m −2 •K −4 ; and ε a is atmospheric long wave emissivity.
The long wave radiation [9] of asphalt pavement emission can be calculated according to the following formula: where q L is the outward amount of pavement radiation, W•m −2 •s −1 ; and ε is pavement emissivity.So the heat flux into the pavement system can be expressed as: where q is the heat flux into the pavement; E g is the total solar shortwave radiation, W•m −2 •s −1 ; ε e is the absorptivity of solar shortwave radiation; and ε f is the absorptivity of atmospheric long wave radiation.

Model of Temperature Field in Asphalt Pavement Based on Green's Function
Green's function method is widely used to solve heat conduction problems, especially when the boundary conditions are nonhomogeneous.It can be used to solve single dielectric problems, as well as inhomogeneous problems for composite dielectrics.The heat conduction problem can be expressed as follows: where T i (z,t) is the distribution of temperature in i-th layer, K; z is the depth, m; and α i is the thermal diffusivity of i-th layer, m 2 •s −1 .
The boundary conditions can be expressed as follows: where k i is the thermal conductivity of i-th layer, W•m −1 •K −1 ; f 1 (t) is the heat flux into the pavement; f 2 (t) is the distribution of temperature as z = z m .Construct the solution of T i (z,t) as the form as follows: where φ i (z) and ψ i (z) are functions of z; f 1 (t) and f 2 (t) are functions of t; and θ i (z,t) is function of z and t.
In addition, ψ i (z) should satisfy the steady state heat conduction problem given as Equation ( 17), subjected to the boundary conditions Equations ( 18)- (21).
Because there is no heat source inside the pavement, the solution of θ i (z,t) can be constructed as θ i (z,t) = Z(z) Γ(t); Equations (31) and (32) can be obtained as θ i (z,t) substituted into Equations ( 22) and (23).
So θ i (z,t) can be expressed as follows: Equation ( 22) can be deduced as: where According to the heat transfer theory, once the Green's function is determined, the temperature field of the pavement can be determined directly as:

High Viscosity Modified Asphalt
Different from other studies, the high viscosity additive (HVA) and styrene-butadiene-styrene (SBS) block copolymer modified asphalt were used in this test.The indexes of asphalt and HVA were measured by the methods in Chinese standards.The results of SBS modified asphalt test are shown in Table 1.The test results of HVA are shown in Table 2 and the indexes of high viscosity modified asphalt modified by high viscosity additives are shown in Table 3.

Mix and Structure Design
Asphalt concrete (AC) is widely used in the traditional pavement.Different from the traditional pavement, the concrete used in the permeable pavement is often called porous asphalt concrete (PAC).Besides, the maximum particle size of upper-surface is usually 13 mm in China, and for the mid-surface, it is usually 20 mm.In addition, for better simulation effect, the test specimens were all made up of two layers.Traditional pavement specimen was made up of 4 cm thick AC13 mixture and 6 cm thick AC20 mixture.The permeable pavement specimens were made up of 4 cm thick PAC13 mixture and 6 cm thick AC20 mixture.Additionally, the different clogging situations were simulated by porosities ranged from 15% to 24%.

Mix Design
The present study shows that the porosity of PAC13 is greatly influenced by the aggregate passing proportion of 4.75 mm (P 4.75 ) and 2.36 mm (P 2.36 ) [20].So the porosity of PAC13 mixture was mainly adjusted by controlling P 4.75 and P 2.36 .The gradation is shown in Table 4 and the gradation curves are shown in Figure 2.

Mix Design
The present study shows that the porosity of PAC13 is greatly influenced by the aggregate passing proportion of 4.75 mm (P4.75) and 2.36 mm (P2.36) [20].So the porosity of PAC13 mixture was mainly adjusted by controlling P4.75 and P2.36.The gradation is shown in Table 4 and the gradation curves are shown in Figure 2.  In addition, the gradations of aggregate blends for asphalt mixture AC13 and AC20 are listed in Tables 5 and 6.

Structure Design
A total of 6 experimental groups were set up in the test.The size of the specimens was 30 cm × 30 cm × 10 cm.They were divided into 2 types: traditional group (AC) and permeable group (S1-S5).The specific combination form is shown in Table 7.Additionally, in order to facilitate later experiment, every gradation of asphalt mixture was also made into rutting plate specimens in the size of 30 cm × 30 cm × 5 cm (Labeling AC13, AC20 and PAC13-1-PAC13-5).In addition, the gradations of aggregate blends for asphalt mixture AC13 and AC20 are listed in Tables 5 and 6.

Structure Design
A total of 6 experimental groups were set up in the test.The size of the specimens was 30 cm × 30 cm × 10 cm.They were divided into 2 types: traditional group (AC) and permeable group (S1-S5).The specific combination form is shown in Table 7.Additionally, in order to facilitate later experiment, every gradation of asphalt mixture was also made into rutting plate specimens in the size of 30 cm × 30 cm × 5 cm (Labeling AC13, AC20 and PAC13-1-PAC13-5).The temperature field in the pavement needed to be monitored in the later tests.Thus, the thermocouple sensors should be embedded in the specimen preparation process.The length of the thermocouples inside the specimen was 15 cm, and the depth of 2 cm, 4 cm, 7 cm and 10 cm respectively, which was shown in Figure 3. Two thermocouples were embedded in each depth, and the temperature of each layer is expressed as the average value of the measured values of the two sensors.In order to facilitate the embedment of the sensor, the customized model in Figure 4a was made.The specimens were prepared by the method of layered compaction.After the lower layer was cooled, the upper layer was added and compacted.The temperature field in the pavement needed to be monitored in the later tests.Thus, the thermocouple sensors should be embedded in the specimen preparation process.The length of the thermocouples inside the specimen was 15 cm, and the depth of 2 cm, 4 cm, 7 cm and 10 cm respectively, which was shown in Figure 3. Two thermocouples were embedded in each depth, and the temperature of each layer is expressed as the average value of the measured values of the two sensors.In order to facilitate the embedment of the sensor, the customized model in Figure 4a was made.The specimens were prepared by the method of layered compaction.After the lower layer was cooled, the upper layer was added and compacted.

Determination of Thermal Properties
Parameters in the heat conduction model included mass density, specific heat and thermal conductivity of the mixture.The mass density of the mixture was expressed by the bulk density.In addition, in the study of the thermally physical properties of asphalt mixtures, Zou [21] proposed that the specific heat can be expressed in a parallel model, which is shown as follows: where Ac is the specific heat of asphalt mixtures; Ai is the specific heat of the i-th component; Bi is mass fraction of the i-th component.
The mass densities and specific heat of each gradation are shown in Table 8.The temperature field in the pavement needed to be monitored in the later tests.Thus, the thermocouple sensors should be embedded in the specimen preparation process.The length of the thermocouples inside the specimen was 15 cm, and the depth of 2 cm, 4 cm, 7 cm and 10 cm respectively, which was shown in Figure 3. Two thermocouples were embedded in each depth, and the temperature of each layer is expressed as the average value of the measured values of the two sensors.In order to facilitate the embedment of the sensor, the customized model in Figure 4a was made.The specimens were prepared by the method of layered compaction.After the lower layer was cooled, the upper layer was added and compacted.

Determination of Thermal Properties
Parameters in the heat conduction model included mass density, specific heat and thermal conductivity of the mixture.The mass density of the mixture was expressed by the bulk density.In addition, in the study of the thermally physical properties of asphalt mixtures, Zou [21] proposed that the specific heat can be expressed in a parallel model, which is shown as follows: where Ac is the specific heat of asphalt mixtures; Ai is the specific heat of the i-th component; Bi is mass fraction of the i-th component.
The mass densities and specific heat of each gradation are shown in Table 8.

Determination of Thermal Properties
Parameters in the heat conduction model included mass density, specific heat and thermal conductivity of the mixture.The mass density of the mixture was expressed by the bulk density.In addition, in the study of the thermally physical properties of asphalt mixtures, Zou [21] proposed that the specific heat can be expressed in a parallel model, which is shown as follows: where A c is the specific heat of asphalt mixtures; A i is the specific heat of the i-th component; B i is mass fraction of the i-th component.
The mass densities and specific heat of each gradation are shown in Table 8.The thermal conductivity of the mixture was related to the conductivity of aggregate, asphalt, mineral powder and porosity.Williamson [22] once put forward a formula: where k a , k b , k v and k w are the thermal conductivity of the aggregate, asphalt, mineral powder and air respectively.G, h, i and j are volume fraction of each component.The thermal conductivities [23] of the components were listed in Table 9.In addition, the thermal conductivities of the specimens were also tested by heat flow meter method, the results of the test and Williamson's formula were shown in Table 10.From Table 10, for PAC13 mixtures, the calculation results of Williamson formula and the measured results decreased with the increase of porosity.The increasing trend is consistent, so it was considered whether there was linear relationship between them.To facilitate the analysis and comparison, the data of Table 10 is shown in Figure 5, and the relationship between the tested and calculated values was obtained by linear fitting.
Can be seen from Figure 5, it can be linear fitted well between test value and the theoretical value of the thermal conductivity of PAC mixture.So in this paper, a new model was proposed based on the Williamson formula to explain the relationship between thermal conductivity of the mixture and that of the components, which is shown as Equation (43):

Test Method
The infrared lamp was used to simulate the solar radiation in summer.Considering the size ratio of the vertical and horizontal directions of the road surface, the sides and the bottom of the rutting plate specimens can be set as an adiabatic surface while the upper surface was set as a mixed surface to exchange heat with the outside.
As the temperature of the road surface in summer was about 65-70 °C, the target of the equilibrium temperature of the heating test was set in the gap.The AC specimen was used to adjust the height of the infrared light.When the specimen temperature reaches the equilibrium range, the height of the infrared lamp would be fixed.The result showed that when the height of the infrared lamp was controlled at 28 cm above the surface and the test lasting time was 60 min, the temperature of the specimen surface would reach the target range and kept stable.
The heating tests were conducted after the height was determined.The infrared lamp should be cooled after each heating test to avoid the influence of the remaining warmth on the later test.The experimental devices were shown in Figure 6a.Additionally, in order to prevent the influence of wind, the enclosure was installed around the specimen, shown in Figure 6b.

Data Comparison
Several hypotheses were adopted in the theoretical calculation.
(1) The heating power on the infrared light remained constant during the test; (2) The sides and the bottom of the specimen were adiabatic boundaries, and the upper surface of the specimen is a mixed boundary of fluid and solid; (3) Each layer of asphalt mixture was isotropic materials; (4) The thermal properties of asphalt mixture remained constant during the test; (5) There was no thermal resistance between two linked layers.

Test Method
The infrared lamp was used to simulate the solar radiation in summer.Considering the size ratio of the vertical and horizontal directions of the road surface, the sides and the bottom of the rutting plate specimens can be set as an adiabatic surface while the upper surface was set as a mixed surface to exchange heat with the outside.
As the temperature of the road surface in summer was about 65-70 • C, the target of the equilibrium temperature of the heating test was set in the gap.The AC specimen was used to adjust the height of the infrared light.When the specimen temperature reaches the equilibrium range, the height of the infrared lamp would be fixed.The result showed that when the height of the infrared lamp was controlled at 28 cm above the surface and the test lasting time was 60 min, the temperature of the specimen surface would reach the target range and kept stable.
The heating tests were conducted after the height was determined.The infrared lamp should be cooled after each heating test to avoid the influence of the remaining warmth on the later test.The experimental devices were shown in Figure 6a.Additionally, in order to prevent the influence of wind, the enclosure was installed around the specimen, shown in Figure 6b.

Test Method
The infrared lamp was used to simulate the solar radiation in summer.Considering the size ratio of the vertical and horizontal directions of the road surface, the sides and the bottom of the rutting plate specimens can be set as an adiabatic surface while the upper surface was set as a mixed surface to exchange heat with the outside.
As the temperature of the road surface in summer was about 65-70 °C, the target of the equilibrium temperature of the heating test was set in the gap.The AC specimen was used to adjust the height of the infrared light.When the specimen temperature reaches the equilibrium range, the height of the infrared lamp would be fixed.The result showed that when the height of the infrared lamp was controlled at 28 cm above the surface and the test lasting time was 60 min, the temperature of the specimen surface would reach the target range and kept stable.
The heating tests were conducted after the height was determined.The infrared lamp should be cooled after each heating test to avoid the influence of the remaining warmth on the later test.The experimental devices were shown in Figure 6a.Additionally, in order to prevent the influence of wind, the enclosure was installed around the specimen, shown in Figure 6b.

Data Comparison
Several hypotheses were adopted in the theoretical calculation.
(1) The heating power on the infrared light remained constant during the test; (2) The sides and the bottom of the specimen were adiabatic boundaries, and the upper surface of the specimen is a mixed boundary of fluid and solid; (3) Each layer of asphalt mixture was isotropic materials; (4) The thermal properties of asphalt mixture remained constant during the test; (5) There was no thermal resistance between two linked layers.

Data Comparison
Several hypotheses were adopted in the theoretical calculation.
(1) The heating power on the infrared light remained constant during the test; (2) The sides and the bottom of the specimen were adiabatic boundaries, and the upper surface of the specimen is a mixed boundary of fluid and solid; (3) Each layer of asphalt mixture was isotropic materials; (4) The thermal properties of asphalt mixture remained constant during the test; (5) There was no thermal resistance between two linked layers.
Water 2018, 10, 299 11 of 16 The AC specimen was set as the matched group.With the thermally physical properties of the AC specimen plugged into Equation (40) and surface temperature curve as the reference, the power of the infrared radiation was determined to be 0.157 W/mm 2 .Then the thermal properties of the remaining specimens were also plugged into Equation (40).The results of theoretical calculation and tests are shown in the following figures.
Firstly, it can be found that the theoretical calculation results are similar with the actual results.So it can be determined that the model was suitable with the heat transfer of porous asphalt mixture, and it can be used to predict the temperature of the whole pavement.
In addition, it can be seen from Figures 7 and 8 that the temperature on the surface of the specimen rose fastest.As the test went on, the temperature rising curve of the specimen gradually became gentle and tended to be stable.In addition, by comparing the temperature among different depths, it can be concluded that the temperature conduction in the asphalt pavement had the hysteresis in the vertical direction.With the increase in depth, the hysteresis phenomenon became more obvious.Taking AC specimen as an example, the surface temperature changes rapidly in the first 5 min.But at the bottom, the temperature had almost no change in 30 min.
Water 2018, 10, x FOR PEER REVIEW 11 of 16 The AC specimen was set as the matched group.With the thermally physical properties of the AC specimen plugged into Equation (40) and surface temperature curve as the reference, the power of the infrared radiation was determined to be 0.157 W/mm 2 .Then the thermal properties of the remaining specimens were also plugged into Equation (40).The results of theoretical calculation and tests are shown in the following figures.
Firstly, it can be found that the theoretical calculation results are similar with the actual results.So it can be determined that the model was suitable with the heat transfer of porous asphalt mixture, and it can be used to predict the temperature of the whole pavement.
In addition, it can be seen from Figures 7 and 8 that the temperature on the surface of the specimen rose fastest.As the test went on, the temperature rising curve of the specimen gradually became gentle and tended to be stable.In addition, by comparing the temperature among different depths, it can be concluded that the temperature conduction in the asphalt pavement had the hysteresis in the vertical direction.With the increase in depth, the hysteresis phenomenon became more obvious.Taking AC specimen as an example, the surface temperature changes rapidly in the first 5 min.But at the bottom, the temperature had almost no change in 30 min.The AC specimen was set as the matched group.With the thermally physical properties of the AC specimen plugged into Equation (40) and surface temperature curve as the reference, the power of the infrared radiation was determined to be 0.157 W/mm 2 .Then the thermal properties of the remaining specimens were also plugged into Equation (40).The results of theoretical calculation and tests are shown in the following figures.
Firstly, it can be found that the theoretical calculation results are similar with the actual results.So it can be determined that the model was suitable with the heat transfer of porous asphalt mixture, and it can be used to predict the temperature of the whole pavement.
In addition, it can be seen from Figures 7 and 8 that the temperature on the surface of the specimen rose fastest.As the test went on, the temperature rising curve of the specimen gradually became gentle and tended to be stable.In addition, by comparing the temperature among different depths, it can be concluded that the temperature conduction in the asphalt pavement had the hysteresis in the vertical direction.With the increase in depth, the hysteresis phenomenon became more obvious.Taking AC specimen as an example, the surface temperature changes rapidly in the first 5 min.But at the bottom, the temperature had almost no change in 30 min.The temperature of each depth of the specimens can also be obtained by the curves shown in the figures above.The temperature in the depth of 4 cm and 10 cm of each specimen was extracted and is listed in Table 11.It can be visually seen from Table 11 that in both depth of 4 cm and 10 cm, the temperatures in permeable pavements were lower than the traditional pavement.It showed that the cooling effect still existed in the permeable pavement under the extreme conditions of no wind on the surface or water in the void.Through the comparison of single layered PAC asphalt specimens, it can be found that with the increase of porosity, both of the temperature values in the two depths were showing a downward trend.The results showed that the cooling effect of drainage asphalt increased with the increase of void fraction.As the porosity of PAC13 layer changed from 16.68% to 23.05%, the maximum temperature of permeable asphalt pavement at the depth of 4 cm and 10 cm were 0.4-2.4°C and 0.2-2.5 °C lower than the traditional pavement respectively.

The Cooling Performance of the Temperature Field with Different Porosities
The pavement was usually divided into six layers.The upper surface was AC13/PAC13, the middle surface was AC20, and the lower surface was AC25.The order and thickness of each layer was shown in Figure 9, where CTB stands for cement treated base, LS stands for lime-stabilized soil and SG stands for subgrade.The temperature of each depth of the specimens can also be obtained by the curves shown in the figures above.The temperature in the depth of 4 cm and 10 cm of each specimen was extracted and is listed in Table 11.It can be visually seen from Table 11 that in both depth of 4 cm and 10 cm, the temperatures in permeable pavements were lower than the traditional pavement.It showed that the cooling effect still existed in the permeable pavement under the extreme conditions of no wind on the surface or water in the void.Through the comparison of single layered PAC asphalt specimens, it can be found that with the increase of porosity, both of the temperature values in the two depths were showing a downward trend.The results showed that the cooling effect of drainage asphalt increased with the increase of void fraction.As the porosity of PAC13 layer changed from 16.68% to 23.05%, the maximum temperature of permeable asphalt pavement at the depth of 4 cm and 10 cm were 0.4-2.4• C and 0.2-2.5 • C lower than the traditional pavement respectively.

The Cooling Performance of the Temperature Field with Different Porosities
The pavement was usually divided into six layers.The upper surface was AC13/PAC13, the middle surface was AC20, and the lower surface was AC25.The order and thickness of each layer was shown in Figure 9, where CTB stands for cement treated base, LS stands for lime-stabilized soil and SG stands for subgrade.The temperature of each depth of the specimens can also be obtained by the curves shown in the figures above.The temperature in the depth of 4 cm and 10 cm of each specimen was extracted and is listed in Table 11.It can be visually seen from Table 11 that in both depth of 4 cm and 10 cm, the temperatures in permeable pavements were lower than the traditional pavement.It showed that the cooling effect still existed in the permeable pavement under the extreme conditions of no wind on the surface or water in the void.Through the comparison of single layered PAC asphalt specimens, it can be found that with the increase of porosity, both of the temperature values in the two depths were showing a downward trend.The results showed that the cooling effect of drainage asphalt increased with the increase of void fraction.As the porosity of PAC13 layer changed from 16.68% to 23.05%, the maximum temperature of permeable asphalt pavement at the depth of 4 cm and 10 cm were 0.4-2.4°C and 0.2-2.5 °C lower than the traditional pavement respectively.

The Cooling Performance of the Temperature Field with Different Porosities
The pavement was usually divided into six layers.The upper surface was AC13/PAC13, the middle surface was AC20, and the lower surface was AC25.The order and thickness of each layer was shown in Figure 9, where CTB stands for cement treated base, LS stands for lime-stabilized soil and SG stands for subgrade.The thermally physical properties of each layer were shown in Table 12.

Layer
Density (kg/m 3 ) Specific Heat (J After the structural form was determined, it was necessary to ascertain the thermal physical parameters of each layer, the solar radiation intensity and the variation law of the daily air temperature.The amount of solar radiation [12] in the model was expressed as the form of q(t) by Equation (44).
where q 0 is the maximum radiation intensity at midday, q 0 = 0.131 mQ, m = 12/c; Q is the total amount of solar radiation in the day, J/m 2 ; c is the actual effective sunshine time, h; and ω is angular frequency, ω = 2π/24, rad.The changing rule of the daily air temperature [12] was shown in Equation ( 27) where T a is the air temperature, • C; T a is the average value of the temperature in a day, T a = (T max + T min )/2, • C; T m is the daily temperature variation, T m = (T max − T min )/2, • C; t 0 is the initial phase; and T max = 40 • C, T min = 26 • C, t 0 = 3.
Because the temperature at the bottom of the middle surface has the greatest influence on rutting, the change of temperature in the depth of 10 cm within 24 h was studied.The calculation results of AC and S5 were shown in Figure 10.
Water 2018, 10, x FOR PEER REVIEW 13 of 16 The thermally physical properties of each layer were shown in Table 12.After the structural form was determined, it was necessary to ascertain the thermal physical parameters of each layer, the solar radiation intensity and the variation law of the daily air temperature.The amount of solar radiation [12] in the model was expressed as the form of q(t) by Equation (44).
where q0 is the maximum radiation intensity at midday, q0 = 0.131 mQ, m = 12/c; Q is the total amount of solar radiation in the day, J/m 2 ; c is the actual effective sunshine time, h; and ω is angular frequency, ω = 2π/24, rad.The changing rule of the daily air temperature [12] was shown in Equation ( 27) where Ta is the air temperature, °C; a is the average value of the temperature in a day, a = (Tmax + Tmin)/2, °C; Tm is the daily temperature variation, Tm = (Tmax − Tmin)/2, °C; t0 is the initial phase; and Tmax = 40 °C, Tmin = 26 °C, t0 = 3.
Because the temperature at the bottom of the middle surface has the greatest influence on rutting, the change of temperature in the depth of 10 cm within 24 h was studied.The calculation results of AC and S5 were shown in Figure 10.It can be seen from Figure 10 that the temperature would be more stable and the maximum temperature would lower in the permeable pavement.It indicated that the permeable pavement could reduce the internal maximum temperature and be more environmental-friendly.The maximum temperatures at the bottom of middle surface of each pavement type were calculated and It can be seen from Figure 10 that the temperature would be more stable and the maximum temperature would lower in the permeable pavement.It indicated that the permeable pavement could reduce the internal maximum temperature and be more environmental-friendly.The maximum temperatures at the bottom of middle surface of each pavement type were calculated and shown in Table 13.Additionally, the cooling performance was represented by the reduced temperature between AC and permeable pavements.From the theoretical calculation results, it can be seen that the cooling performance would be better with the increase of porosity.The maximum value of cooling performance could reach 1.18 • C.So it was an effective method to contrast the urban heat island effect by adopting permeable pavement in the building of sponge cities.The results in Table 13 were shown in Figure 11.
Water 2018, 10, x FOR PEER REVIEW 14 of 16 shown in Table 13.Additionally, the cooling performance was represented by the reduced temperature between AC and permeable pavements.From the theoretical calculation results, it can be seen that the cooling performance would be better with the increase of porosity.The maximum value of cooling performance could reach 1.18 °C.So it was an effective method to contrast the urban heat island effect by adopting permeable pavement in the building of sponge cities.The results in Table 13 were shown in Figure 11.In addition, for the permeable pavement, as the porosity was 16.68%, the maximum temperature of the bottom of the middle surface was only 0.29 °C lower than that of traditional pavement.There was almost no cooling performance, and the porosity of permeable pavement was usually 20%, whose cooling performance was about 0.63 °C.So it can be judged that the cooling performance of the permeable pavement would reduce by over 50% as the porosity changed from 20% to 16.68%.Additionally, there was a linear relationship between the cooling performance and porosity in the single layered permeable pavements, which was shown in Figure 12.In addition, for the permeable pavement, as the porosity was 16.68%, the maximum temperature of the bottom of the middle surface was only 0.29 • C lower than that of traditional pavement.There was almost no cooling performance, and the porosity of permeable pavement was usually 20%, whose cooling performance was about 0.63 • C.So it can be judged that the cooling performance of the permeable pavement would reduce by over 50% as the porosity changed from 20% to 16.68%.Additionally, there was a linear relationship between the cooling performance and porosity in the single layered permeable pavements, which was shown in Figure 12.
Water 2018, 10, x FOR PEER REVIEW 14 of 16 shown in Table 13.Additionally, the cooling performance was represented by the reduced temperature between AC and permeable pavements.From the theoretical calculation results, it can be seen that the cooling performance would be better with the increase of porosity.The maximum value of cooling performance could reach 1.18 °C.So it was an effective method to contrast the urban heat island effect by adopting permeable pavement in the building of sponge cities.The results in Table 13 were shown in Figure 11.In addition, for the permeable pavement, as the porosity was 16.68%, the maximum temperature of the bottom of the middle surface was only 0.29 °C lower than that of traditional pavement.There was almost no cooling performance, and the porosity of permeable pavement was usually 20%, whose cooling performance was about 0.63 °C.So it can be judged that the cooling performance of the permeable pavement would reduce by over 50% as the porosity changed from 20% to 16.68%.Additionally, there was a linear relationship between the cooling performance and porosity in the single layered permeable pavements, which was shown in Figure 12.

Figure 1 .
Figure 1.Heat transfer diagram of multilayered pavement system.

Figure 1 .
Figure 1.Heat transfer diagram of multilayered pavement system.

Figure 6 .
Figure 6.(a) The experimental devices; (b) The enclosure around the specimen.

Figure 6 .
Figure 6.(a) The experimental devices; (b) The enclosure around the specimen.

Figure 6 .
Figure 6.(a) The experimental devices; (b) The enclosure around the specimen.

Figure 7 .Figure 7 .
Figure 7. Results of theoretical calculation and tests of AC.

Figure 7 .Figure 8 .
Figure 7. Results of theoretical calculation and tests of AC.

Figure 8 .
Figure 8.(a) Results of theoretical calculation and tests of S1; (b) Results of theoretical calculation and tests of S2; (c) Results of theoretical calculation and tests of S3; (d) Results of theoretical calculation and tests of S4; (e) Results of theoretical calculation and tests of S5.

Figure 8 .
Figure 8.(a) Results of theoretical calculation and tests of S1; (b) Results of theoretical calculation and tests of S2; (c) Results of theoretical calculation and tests of S3; (d) Results of theoretical calculation and tests of S4; (e) Results of theoretical calculation and tests of S5.

Figure 10 .
Figure 10.Temperature at the bottom of the middle surface.

Figure 10 .
Figure 10.Temperature at the bottom of the middle surface.

Figure 11 .
Figure 11.Cooling performance of the structures.

Figure 12 .
Figure 12.Relationship between the cooling performance and porosity.

Figure 11 .
Figure 11.Cooling performance of the structures.

Figure 11 .
Figure 11.Cooling performance of the structures.

Figure 12 .
Figure 12.Relationship between the cooling performance and porosity.

Figure 12 .
Figure 12.Relationship between the cooling performance and porosity.

Table 1 .
Properties of SBS modified asphalt.

Table 2 .
Properties of high viscosity additive.

Table 3 .
Properties of high viscosity modified asphalt.

Table 4 .
Gradation of aggregate blends for asphalt mixture PAC13.

Table 4 .
Gradation of aggregate blends for asphalt mixture PAC13.

Table 7 .
Structure design of specimens.

Table 7 .
Structure design of specimens.

Table 7 .
Structure design of specimens.

Table 8 .
Specific heat and density of 5 cm thick specimens.

Table 9 .
Thermal conductivities of the components.

Table 11 .
Temperature in the specimens.

Table 11 .
Temperature in the specimens.

Table 11 .
Temperature in the specimens.

Table 12 .
Thermally physical properties of each layer.

Table 12 .
Thermally physical properties of each layer.

Table 13 .
Cooling performance of all structures.

Table 13 .
Cooling performance of all structures.

Table 13 .
Cooling performance of all structures.