Abstract

The hydronic snow melting pavement (HSMP) system is an environmentally friendly, clean, and sustainable alternative to traditional approaches for bridge deck snowmelt. The objective of this paper is to investigate the temperature field and thermal responses of HSMP by a three-dimensional finite element (3D FE) model based on the thermal-fluid coupling method. Considering the full fluid domain, the model simulates the dynamic temperature field and obtains the dynamic heat load of HSMP. Results show that the model factually simulates the decrease of fluid temperature along the pipe. The insulation of bridge deck bottom reduces heat loss, and heating demand can be lowered. Due to the ambient temperature changes, preheating is an effective approach of energy conservation and the start time of heating proposed is at 10 : 00 to 16 : 00. The flow velocity has a slight influence, and the recommended magnitude is 0.6 m/s. The shallower pipe embedded depth and the narrower pipe spacing can improve the surface temperature of HSMP and the uniformity of snowmelt process. Under ambient temperature loads, the maximum principle tensile stress of HSMP is induced at the contact interface between pipes and surrounding concrete, and the magnitude is greater than that of the conventional pavement. Under ambient temperature loads and fluid circulation, heated pipes can effectively prevent thermal shrinkage cracking and extend service life of HSMP. With increase of the pipe embedded depth and decrease of the pipe spacing, the chance of thermal shrinkage cracks decreases. For both optimum snowmelt efficiency and thermal cracking reduction, the pipe embedded depth of 7 cm and pipe spacing of 10 cm are recommended for HSMP with the inlet fluid temperature of 15°C.

1. Introduction

To prevent snow accumulation and ice formation on pavement surfaces and improve traffic safety during winter, the applications of chemical salt are the most common means. As reported, annually nearly 35 million tons of salt was used to melt snow and deice ice worldwide, and this dosage has apparently risen in recent years [1]. However, the usage of salt is ineffective when the ambient temperature falls below −3.9°C [2]. Moreover, chemical salt accelerates the deterioration of pavement structure [3], corrosion of steel [4], and pollutes soil and water [5]. As an environmentally friendly, clean, and sustainable alternative, a thermal snowmelt technology can provide one of the solutions [6]. The thermal method includes hydronic pipes [2, 7, 8], electric heating cables [9, 10], and electrically conductive concrete [11, 12]. Among these three technologies, the hydronic heating system in which the heat is extracted from geothermal, solar energy, and industrial waste heat has received the most attention due to its certain inherent advantages such as the improvement of energy efficiency and reduction in emissions [6, 13].

The hydronic snow melting pavement (HSMP) system is that heated fluid circulates through a series of serpentine or parallel embedded pipes inside the pavement to transfer heat to the pavement and consequently melt snow on the top surface. Therefore, the cleaner and more sustainable hydronic snow melting system has received increasing attention in recent years. Research on the hydronic snowmelt system has been performed since the introduction of the technology in 1948 in Klamath Falls, Oregon, USA [14, 15]. Balbay and Esen [6, 7] investigated the feasibility of a ground source heat pump system consisting of the vertical type single U-borehole heat exchangers with three different depths (30 m, 60 m, and 90 m) for snow melting on pavement and bridge decks through experiments and the finite element models and obtained congruence between the experimental data and simulations. Liu et al. [16] built a hydronic snowmelt system which could successfully remain the surface snow-free even during severe weather conditions. Rees et al. [17] developed a two-dimensional transient model, and Liu et al. [18] modified the model to simulate a bridge deck hydronic snowmelt system. The modified model could predict the surface temperature and outlet temperature, but the results were slightly higher than the experimental results.

In the open literature, a number of models could predict the surface temperature of hydronic snowmelt pavement, while only a few [1, 13] focused on the fluid domain, which was considered as the mean heat output [19, 20]. The surface temperature was even along the pipe. However, it is untrue since the fluid temperature decreases along the pipe. For the low fluid temperature, the end of the pipe circuit should be controlled according to the lowest surface temperature. It is necessary to develop a model that takes the full 3D domain into account [21]. The dynamic temperature field of HSMP during all day is influenced by the external climatic conditions. Li and Hong [22] reported the heat load of HSMP obtained from the dynamic temperature field had an effect on the efficiency of snow melt. Therefore, further studies are needed to develop an appropriate control system for the low temperature HSMP system.

The effects of embedded pipes on the mechanical properties of HSMP have been the subject of many studies. Numerical models have been developed and validated by experimental results, which are applicable to perform design with respect to the structural aspect [15]. Van Vliet et al. [23] performed different numerical analyses to quantify the stress and strain in the structural layer, and showed that the shallow depth hollow system could be applied to airfield pavements. Tan et al. [24] compared the responses of traditional versus snowmelt airfield rigid pavement under aircraft loads, temperature loads, and their coupling effects by the means of field tests and a three-dimensional finite element (3D FE) model. The results showed that the embedded pipe did not change the mechanical-induced responses of the pavement. However, there are few reports about the influence of the pipe layouts on the thermal responses of HSMP.

The present paper is aimed to evaluate the effects of heat pipes on the temperature field and thermal responses of conventional pavement and HSMP. Based on the thermal-fluid coupling method, a 3D FE model of HSMP was developed by applying ANSYS. The model simulated the dynamic temperature field under ambient temperature and heat pipes and obtained the temperature load of HSMP. In order to obtain the optimal performance of the HSMP system, influence factors such as boundary conditions of the bottom, start time of heating, flow velocity, pipe embedded depth, and pipe spacing on the thermal analyses of HSMP were studied. Simultaneously, the thermal responses of conventional pavement and HSMP were simulated to identify the effects of pipe layout. Finally, the recommended pipe embedded depth and pipe spacing were proposed for the HSMP with the inlet fluid temperature of 15°C. The study presented would provide theoretical support for the design and application of the bridge deck snow melting system.

2. Methods

Aiming to evaluate the effects of the heat pipe on the HSMP temperature distribution and thermal responses, the transient conduction and dynamic heat load are taken into account in the thermal analyses and responses for the low temperature (15°C) HSMP system. In this study, the thermal-fluid coupling method models the heat conduction between pipes and concrete by applying the thermal-fluid pipe element, which can fully reflect temperature change law of the place near the pipe or far from the pipe and consider temperature changes of the fluid along pipes [25].

2.1. Assumptions of the Model

The model of the HSMP system is performed in pipes and concrete materials for bridge decks. The important assumptions applied in this model are as follows [10, 26]:(1)The heat transfer process of the hydronic snow melting system is transient(2)Neglect the thermal contact resistance between different materials(3)The velocity profile is fully developed within the entire pipe section(4)Materials are uniform, continuous, and isotropic elastic(5)The interface between materials is fixed

2.2. Mathematical Model of HSMP

The heat conduction in the surface course is a three-dimensional and transient conduction process. The energy equation in the surface course can be written by using Fourier’s law as follows:where T is temperature (°C); τ is time (s); x, y, and z are the coordinates; α = λ/(ρc); λ is the thermal conductivity (W/m·°C); ρ is density (kg/m3); and c is the specific heat capacity (J/kg·°C).

Figure 1 shows the heat exchange of the coupled thermal-fluid pipe. The concrete and pipes are simulated by 3D solid elements and the thermal-fluid pipe element, respectively. The thermal-fluid pipe element includes two primary nodes (I, J) and two supplemental nodes (K, L). The heat convection of the pipe and concrete is simulated by the coupling between the supplemental nodes of the thermal-fluid pipe element and the nodes of the concrete element.

The governing equation for heat transfer and conservation of energy is described aswhere is the specific heat capacity matrix; is the node temperature rate vector; is the conduction matrix; is the node temperature vector; is the node heat flux vector; A is the cross-sectional area of the pipe; is the density of water; is the flow rate of water; is the specific heat capacity of water; is the thermal conductivity of water; L is the element length; ; D is the hydraulic diameter; k4 and k5 represent the energy changes attributable to the mass transport that is lumped at the outlet nodes, by which the fall in water temperature along the pipe is calculated; k2 and k3 represent the heat convection of water; and is the convective heat transfer coefficient of water and is calculated by [27].where Nu is the Nusselt number; Re is the Reynolds number, ; Pr is the Prandtl number, ; is the flow velocity of water; is the viscosity of water; and d is the diameter of pipe. According to McAdams [28], a = 0.023; b = 0.8; when the fluid is heated, c = 0.4; when the fluid is cooled, c = 0.3. There are the service conditions of equation (3), i.e., Re = 10000–120000; Pr = 0.7–120; and the temperature difference is less than 30°C. In this paper, the outer diameter and inner diameter of the pipe is 0.022 m and 0.02 m, respectively. The length of the pipe is between 20-30 m, and the flow velocity of heating water is 0.6–2 m/s. Thus, we get L/d = 1000–1500, Pr = 7.7, and Re = 10909–36364, which meet the use conditions of equation (3).

2.3. Fluid Temperature Calculations along Pipes

The geometrical model of fluid temperature calculations along the pipe is shown in Figure 2.

In the numerical simulation of fluid-structure interactions, the fluid in heating pipes is considered to be one-dimensional steady flow fluid. According to Fourier’s law, the heat flux of concrete in contact with pipes is . During , heat exchanges between fluid and concrete (from section W1 to section W2) are expressed as follows:(1)The heat flux of the concrete interface transfers from the fluid and is calculated using the following equation:(2)The fluid heat flux at the inlet section W1 is calculated as follows:(3)The fluid heat flux at the outlet section W2 is calculated by the following equation:where and are the inlet fluid temperature and outlet fluid temperature, respectively.(4)The fluid heat flux difference between section W1 and section W2 is calculated as follows:where, is the fluid temperature; l is the coordinate along the pipe; AP is the flow area.(5)Based on the energy balance, the governing equation is described as

The fluid temperature difference is calculated as

Due to the smaller volume of flows and slight change of fluid temperature, equation (9) is simplified as

The inlet temperature of segments is , by which i segment fluid temperature is calculated as

2.4. Establishment of the 3D FE Model

The object of this study is shown in Figure 3. In Figure 3(a), serpentine pipes are embedded in the pavement. In Figure 3(b), from top to bottom, the structure of the bridge component is an asphalt layer (AL, 5 cm in thickness), concreter layer (CL, 13 cm in thickness), and reinforced concrete body. A point is above the pipe, and B point is above the centre line between two adjacent pipes. Thermal parameters and elastic parameters of HSMP are shown in Table 1.

A 3D FE model of HSMP was developed to solve the thermal analyses and responses by applying ANSYS. The meshing for the whole model and pipes is presented in Figure 4. The concrete, pipe, and pavement surface (bottom) are modelled by Solid70, Fluid116, and Surf152, respectively. The concrete closest to pipes, where the temperature gradient is steep, are more finely meshed to ensure the temperature and responses to be accurately predicted, as shown in Figure 4(c). The HSMP model domain is the cover area of one of the serpentine pipes (width of 6 m and length depending on the pipe spacing). Design parameters of pipe layouts are shown in Table 2.

2.5. Weather Data

The temperature field of HSMP is dynamic during the 1st day, which produces the dynamic heating loads. To ensure the efficiency of HSMP, the control design needs to know the conditions of the heating pipes to start. In this paper, the periodic ambient condition was adopted. The air temperature is expanded into periodic functions in the linear combination form of two sinusoidal functions, involving a 24 h cycle, as shown in the following equation [29]:where is the maximum value of air temperature (°C); is the minimum value of air temperature (°C); is the angular frequency, ; is the time (s). In the study, , , and wind speed are 0°C, −5°C, and 5 m/s, respectively. Solar radiation, snowfall, and humidity are neglected.

2.6. Boundary Conditions

As reported in Ref. [30], the energy balance at the HSMP surface is related to the atmospheric factors, and the comprehensive heat transfer coefficient is adopted in the study. The sides are adiabatic boundaries [13]. The bottom of CL is adiabatic and convective for HSMP with and without the insulated bottom, respectively.

3. Results and Discussion

3.1. Transient Temperature Field of HSMP

This section presents influence of factors (e.g., insulation, start time of heating, flow velocity, pipe embedded depth, and pipe spacing) on HSMP surface temperature distribution. The relationship between surface temperature and factors were analysed.

3.1.1. Effect of Insulation

(1) Surface Temperature without Fluid Circulation. The time when the temperature of the model domain becomes periodic is called the transition time. In this study, the calculated transition time of the model domain is 72 h, and the periodic temperature is used as the initial temperature of HSMP. The minimum and maximum temperatures of HSMP without heating are at 6 : 00 and 16 : 00, respectively. The heat flux of HSMP is positive from 10 : 00 to 20 : 00, and the other time is negative. The results are in accordance with Li and Hong [22]. It is necessary to control the heating system to prevent deck surface freezing when HSMP is in a state of heat loss or deck surface temperature is below 0°C.

Figures 5 and 6 show the temperature field of HSMP with the insulating bottom and of HSMP without the insulating bottom, respectively. The simulated temperature field is similar to the measured data in Ref. [31]. With the insulation at the bottom of HSMP, the HSMP domain temperature gradient in thickness direction is largest near the top surface due to the heat transfer. The maximum temperature difference is about 0.6°C regardless of at 6 : 00 and 16 : 00. Without the heat transfer on the bottom of HSMP, the temperature is higher than that of the noninsulated bottom at 6 : 00. Without insulating the bottom of HSMP, the temperature gradient is steepest near the surface and bottom, respectively. Due to the heat transfer on the surface and the bottom, the maximum temperature difference (approximately 0.3°C) is smaller than that of the insulated bottom, and the temperature is higher than that of the insulated bottom at 16 : 00.

(2) Surface and Pipe Temperature with Fluid Circulation. HSMP (the pipe spacing of 30 mm and the pipe embedded depth of 11 mm) is heated by the fluid circulation (fluid temperature of 15°C and the flow velocity of 0.6 m/s). When the pipes are heated at 6 : 00, it takes 4.25 h and 4.75 h to rise to above 0°C for the surface of HSMP with the insulated bottom and of HSMP without the insulated bottom, respectively.

Figure 7 indicates the surface temperature (above 0°C) distribution of HSMP. The surface thermal field is approximate to sinusoidal distribution along the x direction and decreases linearly along the pipe. The surface temperature of HSMP without the insulated bottom is lower than that of with the insulated bottom, and it easily generates the stripe phenomena in the snowmelt process. Figure 8 shows the pipe temperature difference between the inlet and the outlet is about 2°C. Due to the shorter length of the serpentine pipe, the bottom boundary condition has slight effect on the pipe temperature. The optimization of pipe layout is obtained through the 3D FE model, which can reduce the cost of experiments.

As shown in Figure 9, the heat flux of the insulated bottom gradually increases and remains constant (about 20 kJ/(h·m3)) with the fluid circulation, whereas the heat flux of the uninsulated bottom rapidly rises. As Liu et al. [32] reported, conductive and convective heat losses exist at the surface. Temperature difference and wind speed have obvious influences on the heat losses of HSMP. Therefore, insulation can be very important to significantly reduce heat loss, and heating demand can be lowered.

3.1.2. Effect of Start Time of Heating

Due to ambient temperature changes, it is necessary to preheat the HSMP at an approriate time for satisfactory snowmelt effect, and it is an effective approach of energy conservation. As Hu [33] reported, HSMP obtained more energy and made great impact on snowmelt through proper preheating. The hours required for the surface temperature above 0°C when HSMP is heated at different points in time are presented in Figure 10. The curve is approximate to the cosine function. The hours at 22 : 00 to 4 : 00 is more than that of other points in time. The heating hours at 4 : 00 to 10 : 00 constantly decrease with the rise in ambient temperature. The heating hours at 10 : 00 to 16 : 00 is least (3.5 h). The heating hours at 16 : 00 to 22 : 00 gradually increase with the drop in ambient temperature. Therefore, the start time of heating proposed is at 10 : 00 to 16 : 00.

3.1.3. Effect of Flow Velocity

The fluid is turbulent and meets equation (3) when the flow velocity is greater than the critical flow velocity (0.55 m/s). Figure 11 shows the temperature time-history curve of a certain point at the surface when the flow velocity increases from 0.6 to 2 m/s. The HSMP surface temperature increases slightly as flow velocity increases, which means that the flow velocity increase does not have a great impact on the heat transfer efficiency, identical to Ref. [34]. Therefore, the effect of flow velocity on the HSMP surface temperature is negligible, and the recommended magnitude is 0.6 m/s.

3.1.4. Effect of Pipe Embedded Depth

The pipe embedded depth represents the heat transfer distance from the pipe to the HSMP surface. Figure 12 shows the surface temperature distribution with the different (7 cm, 9 cm, 11 cm) pipe embedded depth and the pipe spacing (15 cm and 30 cm) after heating for 4.25 h. As the pipe embedded depth increases, the heat transfer distance increase results in the surface temperature decrease, and the temperature curve is more smooth. It needs more time to transfer the energy of heated pipes to the surface of HSMP as the pipe embedded depth increases. Therefore, the deeper the pipe embedded depth is, the less energy the surface of HSMP obtains. The reason is that the distance between A point and the pipe has almost no difference with the distance between B point and the pipe when the pipe embedded depth increases. Thus, the A and B point temperatures tend to be consistent and the temperature curve becomes flattening. In Figure 12(a), as the pipe embedded depth increases from 7 cm to 11 cm, the A and B point temperatures significantly descend, and the surface maximum temperature difference between two adjacent pipes decreases from 0.61°C to 0.11°C. In Figure 12(b), as the pipe embedded depth increases from 7 cm to 11 cm, the A point temperature obviously decreases, while the B point temperature has slight change, and the surface maximum temperature difference between two adjacent pipes decreases from 3.36°C to 1.45°C. Results are in good agreement with the literature data [23]. It shows that appropriate decrease of pipe embedded depth and spacing can improve the surface temperature of HSMP and uniformity of the snowmelt process for the low temperature HSMP system.

3.1.5. Effect of Pipe Spacing

The pipe spacing is the horizontal distance between two adjacent pipes that affects the heat-transfer areas of the heat pipes. Figure 13 shows the surface temperature distribution with different (10 cm, 15 cm, and 20 cm) pipe spacing and pipe embedded depth (7 cm and 11 cm) after heating for 4.25 h. As the pipe spacing increases, the heat-transfer area increase results in the surface temperature decrease, and the temperature curve becomes more steep. As the pipe spacing increases from 10 cm to 20 cm, the surface maximum temperature difference between two adjacent pipes decreases from 1°C to 0.03°C at the pipe embedded depth of 7 cm and from 1.7°C to 0.14°C at the pipe embedded depth of 11 cm, respectively. For a given size of HSMP, the narrower pipe spacing represents more number of pipes and smaller heat-transfer areas. As a result, the surface temperature increase and the temperature curve is more flattening. In other words, the larger pipe spacing presumes less number of pipes and larger heat-transfer areas. As mentioned in Refs. [35, 36], the narrower the pipe spacing is, the more uniform the surface temperature is. As shown in Figures 12 and 13, conclusions can be drawn that the the surface maximum temperature difference between two adjacent pipes become greater when the pipe embedded depth is shallower and the pipe spacing is wider.

It can be inferred that the pipe layouts have dominated influence on the uniformity of HSMP surface temperature. The shallower the pipe embedded depth and the narrower the pipe spacing, the higher the surface temperature and the lower the stripe distribution in the snowmelt process. Table 3 lists the numbers of hours for HSMP surface temperature below 0°C in different pipe layouts. The pipe embedded depth of 7 cm and pipe spacing of 10 cm were recommended for both convenience of construction and optimum snowmelt efficiency.

3.2. Thermal Response Analyses

In this section, the influence of pipe layouts on HSMP thermal responses was analysed and compared with the conventional pavement (without embedded pipes) under temperature load based on the 3D FE model. As the layer interface and contact interface between pipes and surrounding concrete almost have no effect on thermal responses [24, 37], the continuity condition is used at all interfaces. The boundary conditions for the sides and bottom of HSMP are simulated as fixed constraints.

3.2.1. Thermal Response of HSMP under Ambient Temperature Loads

For the sake of investigating the difference of thermal responses between HSMP and conventional pavement under ambient temperature loads, a 3D FE model of conventional pavement with the same size was modelled. The temperature distribution of pavement at 6 : 00 is considered as temperature loads. The maximum principle tensile stress is generated at the bottom of the conventional pavement under a positive temperature gradient load, and the magnitude of AL and CL are 0.4 MPa and 2.2 MPa, respectively. However, the HSMP stress distribution features (e.g., maximum stress position, distribution shape, and magnitude) are evidently different with the conventional pavement. As shown in Figure 14, the stress distribution shape of HSMP (the pipe spacing of 15 mm and the pipe embedded depth of 7 mm) is approximate to the diffraction fringe. The maximum principle tensile stress is occurred at the contact interface between pipes and surrounding concrete, and the magnitude is greater than that of the conventional pavement. The magnitude of AL and CL are 0.43 MPa and 2.78 MPa, respectively. The results are in good agreement with Refs. [24, 38]. Due to the greater line expansion coefficient, the shrinkage deformation of pipes is greater than that of concrete. The concrete around pipes is pulled by the pipes, due to which the greater tensile stress is occurred. As shown in Figure 15, the maximum principle stress and vertical displacement of the pipe are the greatest at the bottom and the top, respectively, due to the free boundary of the HSMP surface.

Table 4 and 5 show the maximum principle stress of HSMP with different pipe layouts under ambient temperature loads. The pipe shrinkage results in the special thermal responses of HSMP. The thermal response of CL is more significantly influenced than that of AL due to the pipe embedded. With the increase of pipe embedded depth and pipe spacing, the maximum principle stress of HSMP has little influence and is in accordance with Ref. [39].

3.2.2. Thermal Response of HSMP with Fluid Circulation

In this section, the thermal response of HSMP is simulated when HSMP suffers the ambient temperature (at 6 : 00) and fluid circulation (15°C). In Figure 16, the maximum principle compressive stress and principle tensile stress of the concrete are generated on the contact interface of the top and bottom of the pipe and of the left and right of the pipe, respectively. Due to the pipe thermal expansion, the concrete above the top or below the bottom of the pipe and the concrete between adjacent pipes are subjected to pressure and tension, respectively. It can be found from Figure 16(b) that a smaller maximum principle tensile stress (1.523 MPa) of HSMP with fluid circulation is induced and the magnitude reduces to 54.8% comparing to conventional pavement. As reported in Ref. [24], the thermal stress of HSMP is attributed to the heat transfer and exchange interaction between fluid circulation and concrete.

Figure 17 illustrates the maximum principle stress and vertical displacement of the pipe. Due to the steep temperature gradient of the pipe and the weaker constraint of concrete above the pipe, the extreme values are at the pipe bottom and the pipe top, respectively. Compared with Figure 15, the stress and displacement distribution of the pipe are obviously different and HSMP with fluid circulation is in a state of “internal dilatation and external shrinkage.”

In Table 6, with the increase of pipe embedded depth, the maximum principle tensile stress and the maximum principle compressive stress decreases and increases, respectively. The deeper the pipe embedded is, the greater the constrain is and the greater the compressive stress of the pipe is. It implies that the chance of thermal shrinkage cracks decreases with increasing the pipe embedded depth. In Table 7, with the increase of the pipe spacing, the maximum principle stress increases. For a given size of HSMP, the larger pipe spacing represents less number of pipes and greater temperature gradients, which increases the thermal shrinkage cracking risk. Therefore, the deeper and narrower the pipe is embedded, the smaller the principle tensile stress of HSMP is. Overall, heated pipes can effectively prevent thermal shrinkage cracks and prolong service life of HSMP. Considering the snowmelt efficiency, the pipe embedded depth of 7 cm and pipe spacing of 10 cm were proposed.

4. Conclusions

For the bridge deck hydronic snow melting system utilizing low temperature water, a 3D FE model of HSMP was built. Influence factors on the temperature distribution and thermal responses of HSMP were studied, respectively. The following conclusions are drawn from this study:(1)The 3D FE model of HSMP based on the thermal-fluid coupling method was developed to solve the heat transmission and thermal responses of HSMP.(2)Particularly, the insulation, start time of heating, pipe embedded depth, and pipe spacing have significant impacts on the snowmelt performance, while the flow velocity has a slight influence.(3)The pipe layout has little impact on the thermal response of HSMP under ambient temperature loads. However, the principle tensile stress and the chance of thermal shrinkage cracks decreases for HSMP with fluid circulation as the pipe embedded depth is deeper and the pipe spacing is narrower.(4)The pipe embedded depth of 7 cm and pipe spacing of 10 cm for HSMP with the inlet fluid temperature of 15°C are recommended for both optimum snowmelt efficiency and the reduction of thermal shrinkage cracking risk.

Data Availability

The data used to support the findings of this study are included within the supplementary information files.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the Doctoral Fund of the Ministry of Education of China (grant no. 20130205110014), the Natural Science Basic Research Plan in Shaanxi Province of China (grant no. 2014JM1005), the Science and Technology Program of the Department of Transportation in Shaanxi Province of China (grant no. 13-16k), and the Science and Technology Planning Project in Henan Province of China (grant no. 182102311091).

Supplementary Materials

In the “insulated bottom” folder, two text files are, respectively, the command-flow of the 3D FE model of HSMP (pipe spacing of 30 mm and pipe embedded depth of 11 cm) and conventional pavement (without embedded pipes) with the insulation at the bottom. In the “noninsulated bottom” folder, two text files are, respectively, the command-flow of the 3D FE model of HSMP (pipe spacing of 30 mm and pipe embedded depth of 11 cm) and conventional pavement (without embedded pipes) without the insulation at the bottom. In the “spacing of 10 mm” folder, three text files are, respectively, the command-flow of the 3D FE model of HSMP with pipe spacing of 10 mm and different pipe embedded depths (7 mm, 9 mm, and 11 cm). In the “spacing of 15 mm” folder, three text files are, respectively, the command-flow of the 3D FE model of HSMP with pipe spacing of 15 mm and different pipe embedded depths (7 mm, 9 mm, and 11 cm). In the “spacing of 20 mm” folder, three text files are, respectively, the command-flow of the 3D FE model of HSMP with pipe spacing of 20 mm and different pipe embedded depths (7 mm, 9 mm, and 11 cm). In the “spacing of 30 mm” folder, three text files are, respectively, the command-flow of the 3D FE model of HSMP with pipe spacing of 30 mm and different pipe embedded depths (7 mm, 9 mm, and 11 cm). (Supplementary Materials)