Thermal Insulation Performance of Radiant Barrier Roofs for Rural Buildings in the Qinba Mountains

Strengthening the roof insulation is the main method for improving the indoor thermal environment of rural buildings during summer. )is study proposed a radiant barrier roof with a roof insulation structure based on climatic conditions in the Qinba Mountains and the characteristics of rural buildings, in which radiant barriers were installed on the underside of the roof in the attic space to enhance roof insulation. To study the thermal insulation performance of radiant barrier roofs, aluminum foil bubble composites were installed on the underside of the roof in the attic space in the Qinba Mountains, and comparative experiments were carried out. Based on the results, the average temperature of the second-floor bedrooms dropped by 1.1°C, thus proving the effectiveness of thermal insulation with radiant barrier roofs.)e thermal insulation performance of radiant barrier roofs made of aluminum foil bubble composite was studied through numerical simulation analysis. Based on the findings of this study, radiant heat transfer contributed 65% to the thermal insulation effect when aluminum foil bubble composites were applied. Based on the analysis of thermal performance and cost, it is recommended to use radiant barrier materials with an emissivity of less than 0.3 and a thermal resistance of no more than 2mK/W. Furthermore, the application of radiant barrier roofs made of aluminum foil bubble composites can have a positive effect on thermal insulation during winter. )is facilitates the widespread use of radiant barrier roofs for rural buildings throughout hot-summer/cold-winter zones.


Introduction
At present, energy consumption in buildings accounts for 36% of total energy consumption in China, while the area of rural buildings and residences is approximately three times that of urban residences [1,2]. Rural buildings have indoor thermal environments of poor quality. When the per capita income of villagers increases to a certain level as the economy grows, energy consumption in rural buildings will rise in a nonlinear manner, and the total energy consumption of rural buildings will exceed that of urban areas. erefore, faced with the rapid development of urbanization-driven construction, improving indoor thermal environments in rural buildings based on local conditions is the main way to realize green and livable rural areas and has important strategic significance for building a sustainable and harmonious society.
Enclosure structures are an important factor affecting indoor thermal environments during summer. e roof is the main source of heat loss in the enclosure structure of buildings, constituting 70% of total heat loss, because roofs are exposed to solar radiation all day long [3][4][5]. erefore, it is especially important to strengthen roof insulation structures. At present, studies on roof insulation design primarily focus on superinsulating materials [6,7], green roofs [8,9], natural ventilation techniques [10,11], reflection techniques [12,13], and other techniques [14][15][16]. A type of smooth and thin aluminum foil layer with low emissivity and high reflectivity is used in reflection techniques to reduce radiant heat transfer in air [17][18][19]. Reflection techniques include cold roofs and radiant barriers, both of which reduce radiant heat transfer in air to lower the quantity of heat flowing into the living space through the roof [20,21]. e energy-saving design of buildings in hot-summer/ cold-winter zones in China is mainly based on heat protection during summer and also takes into account thermal insulation during winter. e Qinba Mountains are located in a northern high-latitude region within a hot-summer/ cold-winter zone in China. In such zones, the temperature difference between day and night during summer is small, and solar radiation is stronger in the same season, whereas winter temperature is low, and solar radiation is still acceptable in the same season; therefore, these zones have a certain potential to utilize such a design. Moreover, the air in these zones contains large amounts of moisture, resulting in typical regional climatic conditions, including hot and humid summer, as well as wet and cold winter. According to surveys, most local houses use rafter tile-type roofs in the attic space for roof insulation. Heat transfer takes place mainly in the air layer within the attic space via convection and radiation. e roof tiles conduct heat to the inner surface of the roof tiles owing to exposure to sunlight during the day. en, heat in the roof tiles is transferred to the cooler roof panels via radiation and convection. Due to the large temperature difference between the roof tiles and the roof panels during summer, radiant heat transfer constitutes a higher proportion than convective heat transfer [22].
Regulating radiant heat transfer in the attic space reduces the quantity of heat entering indoors through the roof. Despite being significantly effective, cold roofs will be covered with dust and moss after being used for a few years, which, in turn, reduces the reflectivity of the roof tiles and enhances living space heating ( Figure 1); therefore, this technique is not suitable for these zones. Installing radiant barriers on the bottom surface of the roof tiles in the attic space can effectively prevent dust accumulation and reduce radiant heat transfer in the air layer within the attic space [23][24][25][26].
Applying radiant barriers to the roof in the attic space is a feasible thermal insulation measure. To assess the actual effects of radiant barriers in the attic space, Ferreira and Corvacho [27] conducted field tests on two buildings in northern Portugal, and the results showed that radiant barriers can effectively reduce maximum indoor temperature. Medina [28] found that a small change in the emissivity of radiant barriers will result in a significant change in heat flux during summer, while reducing heat flowing into indoor space during winter. In recent years, to improve indoor thermal comfort, certain new aluminum foil composites have appeared in radiation barriers, and a layer of insulating materials has been added in the middle of radiant barriers.
is kind of insulating material further reduces the quantity of heat flowing into the attic space through the roof tiles. Kosny et al. [29] tested the thermal insulation effect of the new aerogel-based radiant barriers installed in the attic space and found that the thermal load in the attic space installed with radiant barriers decreased significantly. Reduced thermal load is the result of the combined effect of low emissivity and thermal resistance; however, their study failed to point out the contribution of aerogel-based radiant barriers to radiant heat transfer.
Based on the abovementioned research background, this study proposes a type of roof insulation structure, namely, a radiant barrier roof, which consists of a pitched roof, radiant barriers, and an attic space. Figure 2 shows the structure of the radiant barrier roof.
is roof insulation structure is expected to be extended to rural buildings in hot-summer/ cold-winter zones. is is because residents will also install a ceiling in their homes, which can also form an attic-like space, for decoration purposes, even though most rural buildings currently do not have an attic space. To promote the use of radiant barrier roofs in the Qinba Mountains and assist in designing them, it is necessary to study the thermal insulation performance of radiant barrier roofs and conduct a parametric analysis of the emissivity and thermal resistance of radiant barriers.
To study the thermal insulation performance of a radiant barrier roof, the research team conducted a field experiment in the indoor thermal environment of rural buildings in the Qinba Mountains in the summer of 2019 as well as studied the thermal insulation performance of radiant barrier roofs made of aluminum foil bubble composite through numerical simulation analysis and drew corresponding conclusions.
is study provides new ideas and methods of thermal insulation for rural buildings situated in the Qinba Mountains during summer and is of guiding value to rural constructions in the context of rural revitalization.

Field Experiment.
A comparative experiment was conducted on the thermal performance of radiant barriers under the climatic conditions of Jinxing Village, Hanzhong City, from July 27 to 28, 2019. e experiment was performed on two parallel buildings with similar characteristics (Figure 3). Both buildings were very similar in terms of structure and could be used in the comparative experiment. Nevertheless, they were not identical: Residence 1 had a smaller attic area than Residence 2 (the plan dimensions of Residence 2 were 12.6 m × 8.9 m, while the plan dimensions of Residence 1 were 14.4 m × 9.5 m). To seek conditions that were as similar as possible, the second floors of both buildings were uninhabited, and the windows were closed. e type of roof involved was a pitched roof with a slope of 30°, erected by directly installing roof tiles on the rafters without an insulation layer. Meanwhile, the outer wall was a 240 mm thick solid clay brick wall, with plaster applied to the inner and outer walls but without a thermal insulation layer. Considering that aluminum foil bubble composites are easy to install, low cost, and high efficiency, radiant barriers were made of aluminum foil bubble composites. In Residence 1, aluminum foil bubble composites were applied to the inner surface of the roof tiles in the attic space and directly covered the rafters. By contrast, the structure of the attic space in Residence 2 remained unchanged ( Figure 4). e content of the experiment comprised indoor and outdoor air temperatures, wall temperature, and solar radiation intensity. e test equipment was placed at similar positions in both residences. Air temperature was measured using the TR-72U hygrothermograph, with a test accuracy of ±0.2°C. ree hygrothermographs were used in Residence 1, placed in the attic space, second-floor bedroom, and outdoors, where the outdoor hygrothermograph was placed in a location away from direct sunlight. Two hygrothermographs were used in Residence 2, placed in the attic space and second-floor bedroom. Wall temperature was measured using the LR5021 thermocouple, with a test accuracy of ±0.5°C. e experiment was performed continuously for 24 hours, and data were collected at one-hour intervals. In Residence 1, thermocouples were used to test temperatures at the top of the south-facing roof tiles, the bottom of the roof tiles, the top of the radiant barriers, the bottom of the radiant barriers, the attic floor, the bedroom ceiling, and the outer walls of the attic (in all four directions). Since both Residence 1 and Residence 2 have very similar structures, only the outer surface of the enclosure structure in Residence 1 was measured. In Residence 2, thermocouples were used to measure the wall temperatures of the bottom surface of the south-facing roof tiles, the attic floor, and the bedroom ceiling. Figure 5 shows the layout of the measuring points, while Figure 6 displays the location of some thermocouples. Meanwhile, solar radiation intensity was measured using the TBD-1 radiometer and TBQ-2 pyranometer, where the measurement range of the pyranometer is 0 to 2000 W/m2, and the sensitivity coefficient of the pyranometer is 8.789 μV/(W·m −2 ). Data were collected at one-hour intervals, while the measuring points were located in an outdoor open space with an unobstructed surrounding area.

Numerical
Analysis. ANSYS CFX (SIMULIA Inc., USA) is currently one of the most comprehensive, versatile, and widely used types of CFD software in China. ANSYS CFX is the world's only large-scale commercial software that uses a fully implicit coupling algorithm. Pressure-velocity coupling is achieved using the semi-implicit method for pressurelinked equations (SIMPLE) algorithm. e software adopts the element-based finite volume method and absorbs the numerical accuracy of the finite element method by ensuring the conservative features of the finite volume method. ANSYS CFX has a common physical model that includes flow, heat transfer, and radiation. Advanced algorithms, a wide range of physical models, and powerful pre-and postprocessing capabilities enable ANSYS CFX to demonstrate excellent performance in terms of result accuracy, computational stability, and flexibility. erefore, ANSYS CFX was employed in this study. Experimental data were used for validating the CFX model, while the validated model was used for quantifying the contribution of radiant heat transfer and heat conduction to the reduction of peak Mathematical Problems in Engineering ceiling temperature when aluminum foil bubble composites were applied. In addition, this study also simulated the effect of emissivity and thermal resistance on the peak temperature of the bedroom ceiling. Last, this study investigated the effects of the application of aluminum foil bubble composites to radiant barrier roofs during winter.

Physical Model and Material Composition.
Considering that the subject of this study is the attic space, the model was properly simplified to reduce the number of model calculations, while a three-dimensional model was constructed to test the attic space and the second floor (onemeter tall) of the buildings. Since both test residences had very similar structures, only the three-dimensional model for Residence 1 was constructed to represent two similar residential buildings. e model performed grid generation using hexahedral and tetrahedral grids, where the number of grid points was 685,200. Figure 7 shows the grid of the model. Fluid density is related to pressure and temperature.
When a fluid is in a natural convection state, the density of the fluid mainly depends on temperature, but the effect of pressure is negligible. To enable the governing equations to converge more easily, the Boussinesq model was employed to consider the effect of air lift on the calculation results, namely, ρ � ρ 0 (1 − βΔT)g [30]. e Monte Carlo model was adopted in the radiation model, while the most commonly used k − ε standard model was employed in the turbulence model [31,32]. e material dimensions and thermal parameters for the building structures were set according to Tables 1 and 2, where d is the thickness, ρ is the density, λ is the thermal conductivity, c is the specific heat capacity, r is the thermal resistance, and ε is the surface emissivity.

Control Equations.
Indoor heat transfer follows three basic conservation laws, namely, the laws of conservation of energy, conservation of mass, and conservation of momentum. ese are described numerically through the general control equation as follows: where (zρ φ/zt) is the unsteady term; div(pvφ) is the convection term; div(Γ.grad φ) is the diffusion term; Γ is the generalized diffusion coefficient; and S φ is the source term.
e turbulence dissipation equation (ε-equation) is expressed as follows: where G k and G b are the kinetic energy generation terms arising from average velocity gradient and buoyancy. For compressible fluids, where t is the temperature, Pr t is the Prandtl number, and Pr t � 0.85 in this study. Y M represents the contribution of arterial dilatation in compressible turbulence. C 1ε , C 2ε , and C 3ε are empirical constants, while σ k and σ ε are the Prandtl numbers corresponding to turbulence kinetic energy, k, and dissipation rate, ε, respectively. e values of these parameters are shown as follows: In radiation calculations, the effective radiant energy of a surface generally includes emitted energy and reflected energy on the surface. e quantity of emitted energy on the surface is determined by the emissivity and temperature of the surface, while the quantity of reflected energy on the surface is determined by the incident energy projected from the surrounding surface and the reflectivity of the surface itself.
e effective radiant energy of the surface can be expressed by the following equation: where q ok is the effective radiant energy of surface k; ε k is the emissivity of surface k; i.e., 5.67 × 10 −8 (W/m 2 K 2 ); σ is the Boltzmann constant; ρ k is the reflectivity of surface k; and q ik is the radiant energy projected to surface k.

Boundary Conditions.
In the study of the thermal insulation performance of radiant barriers during summer, the boundary conditions for the enclosure structure were derived from measured data. e outer surface temperatures of the four walls and roof tiles tested in Residence 1 were fitted to fourth-order Fourier series using MATLAB software (MathWorks, Inc., Natick, MA, USA) and then used as the boundary conditions for the wall temperature of the enclosure structure in the numerical model ( Figure 8). Meanwhile, in the study of the thermal insulation performance of radiant barriers during winter, the boundary conditions for the enclosure structure were derived from the Hanzhong typical meteorological year database [33]. Meteorological data on the winter solstice were selected for their research significance owing to the fact that the winter solstice is the day with the shortest sunshine duration in the year. As the Qinba Mountains are relatively rich in solar resources, the constructed enclosure structure not only experiences heat exchange with outdoor air but is also affected by solar radiation. To reflect the local climatic conditions more accurately, outdoor sol-air temperatures in different orientations were used as boundary conditions for the enclosure structure. Outdoor sol-air temperatures in different orientations on the winter solstice were calculated based on meteorological data. Outdoor sol-air temperature includes not only the heating effect of direct radiation, scattered radiation, ground-reflected radiation, and outdoor air temperature on the enclosure structure but also involves the scattering effect of effective long-wave radiation.  According to [34], the outdoor sol-air temperature of the enclosure structure in different orientations during winter in Hanzhong is calculated using the following equation: where t Z and t a are outdoor sol-air temperature and outdoor air temperature (°C), respectively, while q S and q R are the quantity of heat from ground-reflected radiation and solar radiation absorbed by the outer surface of the enclosure structure (W/m 2 ), respectively. e absorption rate of the outer surface of the enclosure structure is 0.7; q e is the effective radiation (W/m 2 ); and α a is the overall heat transfer coefficient of the surface. Table 3shows the calculation results for outdoor air temperature as well as hourly outdoor sol-air temperatures in different directions on the winter solstice. Hourly sol-air temperatures were fitted to fourth-order Fourier series using MATLAB (Figure 8). Based on Figure 8, the fourthorder Fourier series-based outdoor sol-air temperature curves were basically consistent with the discrete values of hourly outdoor sol-air temperature. erefore, the fitted external disturbance functions of the enclosure structure during winter and summer were used as the boundary conditions of the simulation. e initial indoor air temperature during summer was taken from the corresponding measured data, while the initial indoor air temperature during winter, which was derived from thermal comfort indicators, was 12°C.

Results and Discussion
Test data were used to validate the accuracy of the ANSYS CFX model and the effectiveness of the radiant barrier roofs. Numerical model analysis was performed to quantify the contribution of radiant heat transfer and heat conduction to the reduction of peak ceiling temperature. is study also simulated and analyzed the effect of the reflectivity and thermal resistance of radiant barriers on the temperature of the bedroom ceiling. In addition, this study investigated the effect of the application of aluminum foil bubble composites during winter.

Experimental Results.
With respect to the field experiment conducted on two similar buildings in Jinxing Village, Hanzhong City, this section exhibits the temperature of the roof tiles, the air temperature in the attic space, the temperature of the attic floor, the ceiling temperature, and the temperature of the second-floor bedroom in the test residences. Figures 9 and 10 show the outdoor conditions during the field test period (from July 27, 2019, to July 28, 2019). During the test period, the peak total radiation intensity was 974 W/m 2 , and the average radiation intensity was 676 W/m 2 , while the direct solar radiation intensity accounted for 76% of the total radiation intensity. Based on the experimental results, the solar radiation intensity was higher during summer. e average outdoor air temperature was 29.5%, while the daily temperature varied by 5°C. Figure 11shows the temperature of the roof, attic floor, and ceiling of both residences. It can be seen from this figure that there was not much difference in the bottom surface temperature of the roof tiles in both buildings; however, with the use of radiant barriers, there was a significant difference in the inner surface temperature of the roof. e maximum bottom surface temperature of the radiant barriers in Residence 2 was approximately 15.1°C lower than that in Residence 1. is goes to show that the thermal insulation effect of the radiant barriers was significant. At the same time, it can also be seen that, in

Mathematical Problems in Engineering
Residence 1, where radiant barriers were installed, the maximum bottom surface temperature of the radiant barriers was approximately 17.2°C lower than the bottom surface temperature of the roof tiles. ere are two reasons to explain the large temperature differences. One of the main reasons is that as aluminum foil bubble composites are adhered to the rafters, small air cavities remain between the aluminum foil and the clay tiles. e presence of air cavities enables both sides of the aluminum foil to be fully utilized, thereby effectively reducing radiant heat transfer from the roof tiles to the attic floor. On the other hand, since the radiant barriers block heat transfer to the attic floor, the heat trapped in the attic space will heat up the inner surface of the roof tiles, thus leading to an increase in the temperature of the bottom of the clay tiles. e bottom surface temperature of the attic floor was always higher than the temperature of the floor without radiant barriers. Table 4 summarizes the characteristics of the inner surface temperature of the roof and bottom surface temperature of the floor (ceiling). Table 5 shows the maximum, minimum, and average temperatures of the attic space and bedroom. Figure 12 displays the temperature variation in the attic space and bedroom. Based on Figure 12, the use of radiant barriers significantly improved the thermal stability of the attic space, resulting in a 7°C drop in the peak temperature of the attic space during daytime, a 2.1°C decrease in the average temperature of the attic space, and a slight increase in the temperature of the attic space during nighttime. e bedroom temperature in Residence 1 was always lower than the bedroom temperature in Residence 2, while the average bedroom temperature dropped by 1.1°C.

Validation of the CFX Model with Measured Data.
is section compares the estimations made using the CFX model and the experimental data recorded in both test residences. is study constructed two numerical models using the parameters in Tables 1 and 2 and performed numerical simulations, in which all temperatures in the CFX simulation were highly consistent with the experimental data. To illustrate the CFX simulation results, Figure 13 compares the air temperatures in the attic space of Residence 1 and Residence 2 obtained from the simulation and the experiment. Table 6 shows the average difference between the simulation results and the experiment results. According to this table, the largest average temperature difference was 0.3°C.  To understand the specific contribution of each heat transfer mode, this study simulated the variation of peak temperature of the second-floor ceiling under three operating conditions (Table 7). Operating Conditions A, B, and C represented the conditions that aluminum foil bubble composites were applied, that there was no thermal resistance in aluminum foil bubble composites, and that aluminum foil bubble composites were not applied, respectively. By comparing both Operating Conditions A and C, it can be found that aluminum foil bubble composites can effectively reduce the peak temperature of the indoor ceiling, which decreased by 1.7°C. is is the result of the combined effect of both radiant heat transfer and heat conduction. According to the comparison between Operating Conditions B and C, aluminum foil bubble composite, with thermal resistance remaining constant, not only changed the emissivity of the ceiling but also reduced the peak ceiling temperature by 1.1°C. is is the result of radiant heat transfer. Based on these results, it can be concluded that radiant heat transfer and heat conduction account for 65% and 35%, respectively, in the thermal insulation effect of aluminum foil   bubble composites. erefore, the low emissivity of aluminum foil bubble composite has a greater impact on its thermal insulation performance.

Effect of Emissivity and ermal Resistance on Ceiling
Temperature. In this model, the initial thermal resistance of aluminum foil bubble composites used was 0.075 m 2 K/W, and the emissivity values of aluminum foil bubble composites used were 0.1, 0.3, 0.5, 0.7, and 0.9. Figure 14 shows the variation of the maximum temperature of the secondfloor bedroom ceiling with emissivity. As expected, emissivity has a greater effect on ceiling temperature. An emissivity of less than 0.3 can effectively reduce peak ceiling temperature, which also demonstrates the advantage of applying low-emissivity materials to the inner surface of the roof. e initial emissivity of aluminum foil bubble composites used was 0.25, while the thermal resistance values of aluminum foil bubble composites used were R � 0.5, 2, 3.5, and 5. Figure 15 exhibits the effect of thermal resistance of the radiant barriers on peak ceiling temperature. Based on this figure, it can be observed that as thermal resistance increases, the inner surface temperature of the roof decreases significantly. When thermal resistance continues to increase from 2 m 2 K/W, the temperature of the indoor ceiling decreases slowly. erefore, it is not necessary to have high thermal resistance. When thermal resistance is less than 2 m 2 K/W, thermal resistance has a greater effect on peak temperature. Based on the analysis of thermal performance and cost, low thermal resistance results in better performance than high thermal resistance. Hence, it is recommended to use radiant barrier materials with an emissivity of less than 0.3 and a thermal resistance of no more than 2 m 2 K/W.

ermal Insulation Performance of Radiant Barrier Roof during Winter.
Heat entering indoors through the roof has to be reduced during summer. However, a roof requires more heat to enter indoors during winter. erefore, the use of radiant barriers is beneficial to heat proofing during summer but is unfavorable to heating during winter. Hence, this study investigated the thermal insulation performance of radiant barrier roofs during winter. Table 8shows the simulated air temperature in the attic space during winter. As thermal resistance is present in aluminum foil bubble composites, air temperature in the attic space where aluminum foil bubble composites were used was always higher than that in the attic space without radiant barriers. e maximum and minimum air temperatures in the attic space rose by 3.3°C and 5.4°C, respectively. e use of aluminum foil bubble composites enhanced the stability of air temperature in the attic space. Hence, aluminum foil bubble composites are favorable to roof insulation during winter.

Conclusion
is study investigated the thermal insulation performance of radiant barrier roofs. e following conclusions were drawn: (1) e results proved the superiority of the thermal insulation performance of radiant barrier roofs. e average bedroom temperature in the residence to which aluminum foil bubble composites were applied dropped by 1.1°C compared to the traditionally constructed residence. (2) In aluminum foil bubble composites, radiant heat transfer contributed 65% to the thermal insulation effect of radiant barriers. It also explained the   advantage of applying low-emissivity materials to the inner surface of the roof. Based on the analysis of thermal performance and cost, low thermal resistance resulted in better performance than high thermal resistance. It is recommended to use radiant barrier materials with an emissivity of less than 0.3 and a thermal resistance of no more than 2 m 2 K/W. (3) Aluminum foil bubble composites have a positive insulating effect during winter due to the presence of thermal resistance, thereby facilitating the widespread use of radiant barrier roofs in buildings located in hot-summer/cold-winter zones.

Data Availability
e data used to support the findings of this study are included within the supplementary information file.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.