Research Article  Open Access
Qin Zhao, Xiaona Fan, Qing Wang, Guochen Sang, Yiyun Zhu, "Research on EnergySaving Design of Rural Building Wall in Qinba Mountains Based on Uniform Radiation Field", Mathematical Problems in Engineering, vol. 2020, Article ID 9786895, 16 pages, 2020. https://doi.org/10.1155/2020/9786895
Research on EnergySaving Design of Rural Building Wall in Qinba Mountains Based on Uniform Radiation Field
Abstract
How to create a healthy and comfortable indoor environment without causing a substantial increase in energy consumption has become a strategic problem that the development of all countries must face and solve. According to the climatic conditions of Qinba Mountains in China, combined with the characteristics of local rural residential buildings and residents’ living habits, the field survey and theoretical analysis were used to study the thermal environment status and the heating energy consumption condition of local rural residential buildings. The thermal design method of walls for the local rural energysaving buildings based on the indoor uniform radiation field was explored by using the outdoor comprehensive temperature function expressed by the fourthorder harmonic Fourier series as the boundary condition of the wall thermal analysis. ANSYS CFX was adopted to study the suitability of the energysaving wall structure designed by the above method. The results show that the indoor thermal environment of local rural residential buildings in winter is not ideal and the heating energy consumption is high, but this area has the geographical advantage to develop solar energy buildings. It is proposed that the indoor thermal comfort temperature of local rural residential buildings in winter should not be lower than 14°C. When the internal surface temperature of the external walls in different orientations are equally based on the design principle of uniform radiation field, the heat transfer coefficient of the east wall, the west wall, and the north wall of the local rural residential buildings is 1.13 times, 1.06 times, and 1.14 times of the south wall heat transfer coefficient, respectively. The energysaving structural wall with KPI porous brick as the main material and the south wall heat transfer coefficient of 0.9 W/(m^{2}·K) is the most suitable energysaving wall for local rural residential buildings.
1. Introduction
Many research studies have shown that good environmental conditions can enhance people’s productivity and happiness [1–4]. People spend more than 90% of their time indoors, especially for the old, the weak, the sick, and the disabled [5, 6]. Therefore, a comfortable indoor thermal environment is necessary. However, with the improvement of indoor thermal comfort, the proportion of building energy consumption is also increasing [7–9]. At present, energy consumption in buildings already accounts for nearly 40% of total energy consumption in China, a figure which is still rising [10–12]. If the rapid growth of building energy consumption cannot be controlled, the present development model will bring negative effects to the sustainable development of China’s national economy, which is still the top priority of developing countries. Therefore, how to create a healthy, comfortable, and safe building indoor environment without causing a substantial increase in energy consumption has become a strategic problem that countries must face and solve.
Currently, more than 560 million people are living in China’s rural areas [13], with a total building area of 24 billion m^{2} [14], accounting for 50% of the total building area in China [15]. According to the projections of the World Bank, even if China maintains a stable urbanization rate, its rural population will still be more than 500 million by 2030 [16]. If the thermal design and construction of buildings are not rationalized, the energy consumption in rural residential buildings will continue to increase in the future [17–19]. Therefore, promoting and accelerating the sustainable development of rural buildings have a strategic meaning for balancing the contradiction between building energy efficiency and indoor thermal environment improvement.
The energysaving of the envelope is considered to be one of the most effective ways of building energysaving [20–23]. The external wall is the main part of the building envelope [24, 25]. Among the energy consumption of the building envelope, the heat consumption of the external wall accounts for the largest proportion [26, 27]. Thus, it is an important measure to reduce building energy consumption and improve the indoor thermal environment by optimizing the thermal design and structural design of external walls. So far, many scholars have conducted research on the energysaving structure design of the wall [28–31] and indicated that indoor thermal environment quality can be effectively improved, and building energy consumption can be significantly reduced if external walls can be rationalized in design and construction.
Generally, the heat transfer coefficient of the wall with different orientations is considered to be the constant in building thermal design [31, 32]. Nevertheless, for some areas with high solar radiation intensity, the solar radiation intensity varies greatly in different orientations [33, 34]. In the thermal design of the wall, the same heat transfer coefficient is adopted for the external walls with different orientations, which greatly simplifies the thermal design process of the wall. However, the above design method neglects the influence of solar radiation difference in different orientations on wall heat transfer and also fails to consider the economics of the wall structure and the influence of the wall thermal performance on the indoor thermal comfort [31].
The effect of the wall on the indoor thermal environment and energy consumption of the building is influenced by the thermal performance, the construction mode, and the construction materials of the wall. The walls with the same heat transfer coefficient but different construction methods or construction materials have different influences on the indoor thermal environment and building energy consumption. However, in the research of building external walls, there is a lack of research on structural optimization of energysaving walls based on the uniform radiation field, especially the research methods of steadystate theoretical calculation and dynamic numerical simulation [31, 35–37]. Many scholars have conducted researches on rural residential buildings [38–40]. These researches show that it is important to fully consider the climate characteristics and residents’ living habits in building design, which can enable buildings to be environmentally friendly, lowenergy consuming, and sustainable on a longer term. Therefore, it is an important measure to reduce the building energy consumption and improve the indoor thermal environment of the rural residential buildings by reasonably controlling the thermal parameters, construction methods, and construction materials of the external walls according to the specific climatic conditions, building characteristics, and residents’ living habits.
The Qinba mountains are located in the northern high latitude area of the hot summer and cold winter area, between the Qinling mountains and the Baba mountains, including three cities, Hanzhong, Ankang, and Shangluo [41, 42]. The area is hot and humid in summer and cold and wet in winter [43]. The temperature in the lowest month is similar to that in the cold area, and the minimum outdoor temperature can reach about −10°C [44, 45]. The annual sunshine hours are from 1400∼2200 h, and the annual solar radiation amount is from 4190∼5016 MJ·m^{−2}·a^{−1} [46]. Consequently, the thermal effect of indoor and outdoor environment on the building envelope in this area is quite different from that of other areas with hot summer and cold winter. The amount of rural buildings in Qinba mountains area is large and widely distributed. However, the building type is simple and the structure is extensive, which contributes to the poor thermal environment indoor in winter and high heating energy consumption [43, 47]. Therefore, problems such as poor indoor thermal environment and high heating energy consumption are particularly prominent in this area.
Ankang area is located in the intersection zone of Sichuan, Chongqing, Hubei, and Shaanxi provinces, and its climate characteristics and residential forms are representative in Qinba mountains [48, 49]. According to the climate conditions of Qinba mountains, combined with the characteristics of rural residential buildings and the living habits of residents in Ankang area, the field survey and theoretical analysis were used to study the indoor and outdoor thermal environment status and the heating energy consumption condition of local rural residential buildings and explore the thermal design method of walls for the local rural energysaving buildings based on the indoor uniform radiation field. In addition, ANSYS CFX was adopted to study the suitability of the energysaving wall structure designed by the above method. The research results can provide theoretical and technical reference for the construction of rural residential buildings in Qinba mountains and other similar areas and play an important role in promoting the sustainable development of rural buildings.
2. Current Situation of Indoor and Outdoor Thermal Environment and Heating înergy Consumption
2.1. Object Selection
We conducted a field survey of rural residential buildings in Sanyuan Village, Ankang City, and Shaanxi Province of China from January 1 to 3, 2016.
According to the survey, it is found that most of the rural residential buildings in Qinba Mountains are twostory independent buildings with households as units. Most of the local residential buildings face south, with a long eastwest and short northsouth layout. The external walls are mostly 240mmthick solid clay brick wall. The roof is flat or in double slope. Most of the windows are in single frame and made of aluminum alloy or plastics. Based on the above analysis, a representative twostory independent residential building (Figure 1) is selected as the following measurement and analysis.
2.2. Data Acquisition
A field measurement was conducted to obtain the indoor and outdoor thermal environment parameters. The outdoor temperature measuring point is located in the outdoor shade. The measuring point of the solar radiation intensity is arranged outdoors with no shelter. The indoor air temperature measuring point is located 1.5 m above the floor in the center of the room. The monitoring period for the solar radiation intensity measuring points is 10 h, and the acquisition time interval is 1 h. The monitoring period for the air temperature measuring points is 24 h, and the acquisition time interval is 10 min.
We are concerned that the hot wire anemometer and globe thermometer will cause inconvenience to the residents’ life if they are placed indoors for a long time. Therefore, the measurements of the mean radiant temperature and air velocity are conducted by intermittent tests. The measuring points of the indoor air velocity and mean radiation temperature are placed at the height of 1.0 m above the floor in the center of the room. The main information about the test instruments is shown in Table 1.

2.3. Results and Analysis
During the test, the outdoor weather conditions were similar every day. Therefore, the test data on January 2, 2016 are selected for the following analysis.
The changes of outdoor air temperature and solar radiation intensity are shown in Figures 2 and 3, respectively.
Figure 2 shows that from 5 am to 4:00 pm, the outdoor temperature is on the rise, while from 4:00 pm to 12 pm, the outdoor temperature shows a decreasing trend. The outdoor temperature varies from 1.3°C to 6.6°C with the temperature fluctuation of 5.3°C, and the mean temperature of 3.6°C. It can be seen that the outdoor air temperature in this area is lower and the temperature fluctuation is slightly larger in winter.
Figure 3 shows that the local sunshine duration is about 10 h. The peak total radiation intensity is 500.7 W/m^{2}, and the average radiation intensity is 234.6 W/m^{2}. The direct radiation intensity accounts for 64.2% of the total radiation intensity. To conclude, the region has a long sunshine time and relatively high solar radiation intensity, which has a geographical advantage to develop solar buildings.
The changes of indoor air temperature, as shown in Table 2.

Table 2 shows the maximum, minimum, and average temperatures of the master bedrooms and the secondary bedrooms. It is obvious that the master bedrooms are higher than that of the secondary bedrooms, and the indoor temperature of the bedrooms on the second floor is lower than that on the first floor. The fluctuation range of the indoor temperature of the master and secondary bedrooms on the first floor is from 8.7°C to 12.1°C and from 8.0°C to 10.8°C, with the average temperature of 10.5°C and 9.8°C, respectively. The fluctuation range of the indoor temperature of the master and secondary bedrooms on the second floor is from 8.1°C to 12.8°C and 7.5°C to 11.6°C, with the average temperature of 10.0°C and 9.4°C, respectively. To sum up, the indoor temperature in this area is low and fluctuates greatly in winter.
Through the intermittent tests, it is found that the mean radiation temperature is approximately equal to the air temperature. The indoor air velocity varies from 0.03 m/s to 0.12 m/s, and the average value is 0.1 m/s.
In addition, we also surveyed the heating status of local rural residential buildings. It is found that the heating energy for local rural residential buildings in winter is mainly coal, and each household consumes about 3.2 t of coal during the heating period. Even so, most of the residents in this area are still concentrated in the same room in winter. If the residents adopt the auxiliary heating method of the compartment, the heating energy consumption will be 1.5 to 2.0 times of the existing energy consumption. The survey also shows that the percentage of the residents who are satisfied and general with the current heating method and heating effect is 14.1% and 43.7%, respectively. Besides, 42.2% of residents express dissatisfaction with the current heating method and heating effect. They think that the indoor heating energy consumption is high, but the indoor thermal comfort is still poor and the indoor air quality is low.
In summary, the indoor thermal environment of the existing local rural residential buildings is not ideal. The residents’ satisfaction is low, and the heating energy consumption is high. However, this area is rich in solar energy resources and has the potential to develop solar energy buildings.
3. Design Method of EnergySaving Wall Based on Uniform Radiation Field
3.1. Determination of Indoor and Outdoor Calculated Temperature
3.1.1. Indoor Design Calculation Temperature
The common international comfort indexes are operative temperature (To), effective temperature (ET), new effective temperature (SET), predicted mean vote (PMV), and subjective temperature (Tsub).
Yaglon et al. [50] proposed the ET as the indoor thermal comfort evaluation index, but this index had the disadvantage of estimating the impact of humidity on thermal sensation. This evaluation index is more perfect as an evaluation index when the effective temperature ET is further improved to the standard effective temperature SET. It is inconvenient to apply to the design because it is intricate to calculate the skin temperature and wettability [51]. PMV is put forward as the evaluation index of indoor thermal comfort. Although the indicator has been considered perfect, it is difficult for most areas of China to reach the standard of PMV = −0.5∼0.5, and it is also considered unnecessary in rural areas of China [52]. T_{o} is put forward as the evaluation index of the indoor thermal environment. However, this index can only qualitatively evaluate the indoor thermal comfort and does not guide the thermal design of buildings well [53]. Mcintyre [54] brought forward the T_{sub} as an evaluation index of the indoor thermal comfort. The T_{sub} refers to the temperature of a uniform closed environment that can produce the same thermal sensation as the actual environment. In this indoor thermal environment, the air temperature and the mean radiation temperature are equal to each other. The air velocity is 0.1 m/s, and the air relative humidity is 50%.
Based on the above analysis, it is considered that T_{sub} as the evaluation index of the indoor thermal environment is effective and feasible. Within the allowable error range, the T_{sub} can be calculated as follows [54]:where T_{sub} is the subjective temperature; M is the metabolic rate (W/m^{2}); and R_{clo} is the thermal resistance of clothing, and its calculation formula is as follows [54]:where R_{clo} is the total clothing thermal resistance and R_{i} is the thermal resistance for each piece of clothing in units of clo (1 clo = 0.155 (°C·m^{2})/W).
According to the survey on the living habits and clothing of rural residents in Qinba mountains, the following conclusions can be drawn:(1)In winter, the local rural residents are in the period of agricultural leisure. The metabolism rate of the local rural residents is calculated according to the average labor intensity. According to ASHRAE standard 552013 [55], the metabolism rate of the local rural residents is 85 W/m^{2}.(2)According to the literature [56], combined with the dressing thickness of different genders and different ages of the local rural residents, the thermal resistance of local residents’ winter clothing is 1.7 clo.
According to the above formulas (1)∼(3), the T_{sub} is 14.4°C. In winter, it is suggested that the indoor thermal comfort temperature of rural buildings in Qinba mountains should not be lower than 14°C based on the comprehensive analysis of comfort, energy conservation, and economy.
3.1.2. Outdoor Calculated Temperature
The comprehensive effect of outdoor air temperature and solar radiation should be considered in the thermal calculation of buildings in areas with abundant solar energy resources. Therefore, the outdoor comprehensive temperature is more accurate than the outdoor air temperature to describe the outdoor thermal effect of different orientation envelope. The calculation formula of outdoor comprehensive temperature is as follows [57]:where t_{sa} is the outdoor comprehensive temperature; t_{a} is the outdoor air temperature, °C; a_{a} is the total heat transfer coefficient of building wall surface with the value of 23.0 W/(m^{2}·°C) [57]; is the effective radiation or night radiation, W/m^{2}; q_{S} is the ground reflected radiant heat absorbed by the outer surface of building envelope, W/m^{2}; q_{R} is the solar radiation heat absorbed by the outer surface of building envelope, W/m^{2}. The absorption rate of sunlight on the outer surface of the building envelope is 0.7 [57].
The winter solstice has the shortest sunshine time and the lowest solar altitude angle. So, it is typical to select the weather data of this day to analyze the outdoor thermal effect of buildings. On the winter solstice, the outdoor air temperature in Ankang area and the comprehensive outdoor temperature with different orientations are shown in Table 3.

3.2. Heat Transfer Coefficient of External Wall of EnergySaving Residential Buildings
3.2.1. Heat Balance Equation of EnergySaving Residential Buildings
In areas with abundant solar energy resources, the indoor thermal condition of buildings is not only affected by factors such as building construction mode and thermal parameters but also affected by solar radiation intensity and outdoor air temperature.
In the heating period, the indoor heat gain is composed of solar radiation heat obtained from heat collection components as well as indoor residual heat and auxiliary heat consumption. The heat loss includes the heat transfer loss of envelope and the heat consumption of air infiltration. The heat balance equation is as follows [58]:where NLC is the static load factor of the room, W/°C; A_{j} is the area of the envelope j, m^{2}; n is the number of air changes in the room with the value of 1.0 times/h [59]; V_{f} is the internal volume of the room, m^{3}; is the average density of outdoor air with the value of 1.27 kg/m^{3} [59]; c_{p} is the constant pressure specific heat capacity of air with the value of 1.008 kJ/(kg·°C) [59]; is the calculated temperature of indoor air, °C; is the calculated temperature of outdoor air, °C; is the monthly average daily radiation of the total daily solar radiation transmitted to the room by the rest of the envelope except the heat collecting component of the room, W; is the average value of heat generated by indoor people, lighting and nondedicated heating equipment in the room, W; Q_{aux} is the auxiliary heat consumption, W; is the window area of the heat collecting component, m^{2}; is the total average daily radiation amount projected on the window surface, W/m^{2}; η is the efficiency of the heat collecting component with the value of 0.13 [59]; is the equivalent temperature of the outdoor solar radiation and the longwave effective radiation of the envelope i, °C.
3.2.2. Analysis of External Wall Heat Transfer Coefficient
According to formula (5), the calculation formula of net average heat transfer coefficient of external wall is obtained, as shown in the following formula:where , is the outdoor comprehensive temperature of the envelope j, °C; mb is the number of external walls and roof; and nbmb is the total number of external doors and windows of nonheat collecting components.
(1) Characteristics of Representative Rural Residential Buildings in Ankang Area. According to the literature [60] and our field survey, the plans and basic parameters of the representative rural residential buildings in Ankang area are determined, as shown in Figure 4 and Table 4, respectively.
(a)
(b)
 
(2) Correlation Analysis of External Wall Heat Transfer Coefficient Based on Uniform Radiation. Many scholars have studied the influence of indoor nonuniform radiation on human thermal comfort. The results show that asymmetric radiation temperature has a significant adverse effect on human thermal comfort. For the rural buildings with strong solar radiation in Qinba mountains, the discomfort of the human body caused by the indoor nonuniform radiation field will be more serious if the same thermal parameters are adopted for the enclosure structures with different orientations. Therefore, it is an effective way to reasonably determine the thermal resistance of envelopes with different orientations to solve the problems of indoor comfort and energy saving of rural buildings.
Taking the indoor calculated temperature and the comprehensive temperature of the air in different orientations as the boundary conditions, the correlation of the heat transfer resistance for each orientation of the external wall is analyzed based on the principle of uniform radiation, as shown in the following equations:where θ_{i.S}, θ_{i.E}, , and θ_{i.N} are the inner surface temperature of the south, east, west, and north walls, respectively, °C:where is the indoor design calculation temperature, °C; R_{oj} is the heat transfer resistance to external walls and roofs with different orientations, (m^{2}·K)/W; and is the comprehensive outdoor temperature corresponding to the envelope, °C.
From equations (7)∼(8) and Table 3, the correlation of the heat transfer coefficients of external walls with different orientations can be obtained:where U_{j} is the heat transfer coefficient of external walls with different orientations, W/(m^{2}·K).
It can be seen from equation (9) that the heat transfer coefficient of the external walls with different orientations has a quantitative proportional relationship. In order to simplify the calculation process, the heat transfer coefficient U_{S} of the south wall is determined as the basic heat transfer coefficient.
(3) Analysis of Relationship between Auxiliary Heat Consumption and Heat Transfer Coefficient of External Wall. The relationship between the basic heat transfer coefficient and the auxiliary heat consumption can be obtained by taking the indoor and outdoor calculated temperature obtained from the above analysis and the basic parameters of Table 4 into equations (6) and (9), as shown in equation (10) and Figure 5:
Figure 5 shows that the basic heat transfer coefficient has a linear relationship with indoor auxiliary heat consumption. Therefore, for the rural residential buildings which have determined the basic parameters, window wall ratio, and window materials, it is an effective method to reduce the building energy consumption by reasonably controlling the basic heat transfer coefficient of the external wall of the building.
The literature [59] shows that the heat consumption of rural residential buildings in Qinba Mountains during the winter solstice is 61.54 W/m^{2} under the state of natural operation state. It can be seen from Figure 2 that, under the premise of building energy saving of 50%, the basic heat transfer coefficient of the external wall of rural residential buildings in Qinba mountains is 0.9 W/(m^{2}·K). The heat transfer coefficients of external walls with different orientations are as shown in Table 5.

4. Structure Suitability Analysis of EnergySaving Wall
4.1. Selection of Numerical Analysis
Many engineering problems can be summarized as the mathematical problems of solving the control equations under given boundary conditions. However, only a few problems in which the properties of the equations are relatively simple and the geometric boundary conditions are quite regular can be solved analytically. Thus, many researchers have developed effective numerical methods, such as finite element method (FEM) [62, 63], meshless method [64–66], and boundary element method (BEM) [67, 68], when it is impossible to solve engineering problems by the analytical method.
The BEM is to divide the unit only on the boundary of the definition domain and use the function that meets the control equation to approximate the boundary conditions. The boundary element method reduces the dimensions of the problem by one degree, which leads to significant advantages in terms of ease of discretization of the domain and of reducing the overall time to perform the analysis of objects with complex geometries. Although BEM can save the preprocessing time, it requires a basic solution when establishing boundary integral equations, and the singular integral of boundary integral equation is difficult to solve complex structural problems [62]. As an alternative to meshbased methods, meshless methods have attracted much attention in recent years due to (1) the lack of dependence on mesh and avoiding meshing; (2) remeshing and mesh distortion (in large deformation problems); and (3) easy preparation of the initial data. Meshless methods are used to establish a system of algebraic equations for the whole problem domain without the use of a predefined mesh. Meshless methods use a set of nodes scattered within the problem domain as well as sets of nodes scattered on the boundaries of the domain to represent the problem domain and its boundaries. These sets of scattered nodes do not form a mesh, which means that no information on the relationship between the nodes is required, at least for field variable interpolation. Because there is no need to create a mesh and the nodes can be created by a computer in a fully automated manner, the time an engineer would spend on conventional mesh generation can be saved. This can translate to substantial cost and time savings in modeling and simulation projects. But meshless methods developed so far are not really ideal and fail in some categories. For example, some meshless methods require background cells for the integration of system matrices derived from the weak form over the problem domain, and some methods that do not require a mesh at all are less stable and less accurate [68].
The FEM as a general discretization procedure of continuum problems posed by mathematically defined statements has been one of the major numerical solution techniques. The finite element method requires division of the problem domain where the partial differential governing equations are defined into meshes, and each mesh is called a finite element. The FEM not only has high calculation accuracy but also can be applied to various complex geometric models, complex material characteristics, and complex boundary conditions. Another major advantage of the finite element method is that a general purpose computer program can be developed easily to analyze various kinds of problems. In the field of architecture, ANSYS CFX is one of the most comprehensive and universal finite element analysis software. The software has general physical models including flow, heat transfer, and radiation. The software adopts the elementbased finite volume method and absorbs the numerical accuracy of the finite element method. Therefore, ANSYS CFX was adopted for the followup study in this paper.
4.2. Model Establishment
4.2.1. Mathematical Model
The model is controlled by the k∼ε equation of turbulent natural convection [69] and the IMMERSOL (IMMER sed SOLids) equation of Stefan–Boltzmann law [70]:(1)Standard k∼ε model: where Gk is the turbulent kinetic energy generated by the laminar velocity gradient, J; Gb is the turbulent kinetic energy generated by buoyancy, J; Ym is the fluctuation caused by the transitional diffusion; ∂_{j} is the turbulent Prandtl number of the j equation; C1, C2, and C3 all are constants; S is a custom physical quantity; k is turbulent pulsating kinetic energy, J; ε is the dissipation rate of turbulent pulsating kinetic energy, %.(2)IMMERSOL (IMMER sed SOLids) model of Stefan–Boltzmann law. When T_{3} is the solid temperature, the equation is as follows: where c is the specific heat capacity of the solid, J/kg·W; λ is the thermal conductivity of the solid, W/m·K; is the heat flux per unit volume of the solid, W/kg; and t is the time, s.
When T_{3} is the radiation temperature between solids, the equation is as follows:where E_{3} is the net value of the radiant energy gained or lost per unit volume, W/m^{2}; E_{12} is the average radiant force of the first and second fluid surface, W/m^{2}; a is the absorptivity of the medium; S is the scattering coefficient of the medium; and L is the distance between two solid radiation surfaces, s.
4.2.2. Physical Model
The physical model of representative rural residential buildings in Ankang area is shown in Figure 6.
In order to facilitate the numerical simulation, the physical model is simplified as follows: (1) the indoor air body of all the rooms is considered as a whole; (2) the external envelope and internal partition lines are assumed to be isotropic onedimensional thermal conductors whose thermal parameters are calculated by the weighted average of the thermal parameters of each material layer; (3) the influence of external door, window type, and cold air infiltration on the indoor air temperature is ignored; (4) the indoor door is set default to be open, the external door and window to be closed, and the bright window part of the external door to be vent; (5) the influence of indoor human activities on the indoor air temperature is ignored; (6) the sloping roof is simplified to a flat roof with an air layer.
4.3. Mesh Generation
In order to simplify the calculation process of the model, only key parts such as external walls, roofs, doors, and windows are encrypted in the study. The size of the calculating model is given in Figure 4 and Table 4. The total number of the grid formed is 306036, and the minimum grid size is 0.02 m. The physical model of the encrypted building is shown in Figure 7.
4.4. Parameter Setting
Through field survey, it is found that the external walls of rural residential buildings in Ankang area are mostly 240mmthick solid clay brick walls without thermal insulation layer. A few buildings also adopt sandwich insulation walls with solid clay brick, limes and brick, and KP1 porous brick. According to the thermal design method of the wall proposed above, combined with the materials and their combination modes of the local rural building walls, three typical energysaving structural walls are selected to conduct analysis in this paper, which are category I, category II, and category III, respectively. Category I is an energysaving structural wall with KPI porous brick as the main material and a basic heat transfer coefficient of 0.9 W/(m^{2}·K). Category II is an energysaving structural wall with KPI porous brick as the main material and a basic heat transfer coefficient of 0.9 W/(m^{2}·K). Category III is an energysaving structural wall with KPI porous brick as the main material and a basic heat transfer coefficient of 0.9 W/(m^{2}·K). The construction methods and thermal parameters of different categories of energysaving walls are shown in Table 6.
 
Note. According to the literature [71], the heat transfer resistance for the inner surface of the wall is 0.11 (m^{2}·k)/W. The heat transfer resistance of the inner surface of the wall is 0.04 (m^{2}·k)/W, and the thermal conductivity correction coefficient of the insulation layer is 1.1. 
4.5. Setting of Boundary Conditions
Under the premise that indoor temperature is a certain value, the unstable heat transfer of building envelope is determined by the change of the outdoor disturbance and the response of the envelope to the external disturbance. The outdoor disturbances, such as outdoor air temperature and solar radiation intensity acting on the building envelope, take 24 hours or 365 days as a cycle. Although the meteorological conditions generally show periodic changes, the changes are random and irregular. Because the complex disturbance curve can be expressed as a group of simple periodic functions by using the Fourier series expansion (formula (15)), the outdoor comprehensive temperature described by the Fourier function is used in this research:where A_{o} is the zeroorder outdoor disturbance, which is considered to be the average value of the outdoor comprehensive temperature in the calculation period; A_{n} is the amplitude of the outdoor disturbance of the nth sine wave; is the frequency of the external disturbance of the nth sine wave, = 2πn/T (rad); φ_{n} is the initial phase of the external disturbance of the nth sine wave; and is the period of the function, h.
In the calculation of the building heat transfer process, the sampling interval of outdoor disturbance is 1 h, j is the sequence number of sampling value, and n is the order of harmonic with the value of 4. The expression of harmonic function is as follows:
According to Table 3 and formula (16), the fourthorder Fourier series expression of outdoor comprehensive temperature in different orientations on the winter solstice in Ankang area is obtained.
The comprehensive outdoor temperature function in south orientation t_{sa·S} is as follows:
The comprehensive outdoor temperature function in south orientation t_{sa·E} is
The comprehensive outdoor temperature function in south orientation is
The comprehensive outdoor temperature function in south orientation t_{sa·N} is
4.6. Model Verification
The accuracy of numerical simulation is affected by the setting of boundary conditions, the density of mesh division, and the convergence of iteration. Therefore, it is necessary to verify the accuracy of the model. The comparison between the measured data and the simulated data is shown in Figure 8.
Figure 8 shows that the variation trend of the simulation results is generally consistent with the measured data, and the relative error between them is within 8%. The main reasons for the error are the simplification of the physical model before the numerical simulation and the accuracy errors of the testing instruments. Therefore, it is effective and feasible for the physical model to be applied in the following numerical simulation analysis.
4.7. Results and Analysis
The energysaving walls with the same heat transfer coefficient but different construction methods or construction materials have different storage and release capacity for heat energy. Therefore, this paper studies the structural suitability of energysaving wall through the analysis of indoor temperature field.
Taking the above three energysaving walls as examples, the indoor temperature of rooms on different floors is simulated. The indoor temperature of the numerical simulation is compared with the indoor temperature of the test, and the results are shown in Figure 9.
(a)
(b)
Figure 9 shows that the energysaving walls with the same heat transfer coefficient but different construction methods or construction materials have significantly different effects on indoor temperature. Compared with the original buildings, the indoor temperature of the buildings with the above energysaving external wall is significantly improved, and the fluctuation of the indoor temperature is slightly reduced.
When the construction method of building wall adopts category I, the average indoor temperature of the master bedroom and the secondary bedroom on the first floor increased by 13.3% and 15.3%, respectively, reaching 11.9°C and 11.3°C. The average indoor temperature of the master bedroom and the secondary bedroom on the second floor rose by 11.0% and 13.3%, respectively, reaching 11.1°C and 10.65°C. When the construction method of building wall adopts category II, the average indoor temperature of the master bedroom and the secondary bedroom on the first floor increased by 12.4% and 13.3%, respectively, reaching 11.8°C and 11.1°C. The average indoor temperature of the master bedroom and the secondary bedroom on the second floor rose by 8.7% and 11.7%, respectively, reaching 11.1°C and 10.65°C.
Compared with category I and category II, when the construction method of building wall adopts category III, the indoor temperature increase of each room reaches the maximum. The average indoor temperature of the master bedroom and secondary bedroom on the first floor increases by 14.3% and 15.3%, respectively, reaching 12.0°C and 11.3°C. The average indoor temperature of the master bedroom and secondary bedroom on the second floor increases by 14.0% and 14.9%, respectively, reaching 11.4°C and 10.8°C. These temperature indexes are close to the lower limit of indoor comfort temperature of the rural residential buildings in winter [17]. It can be seen that, under the climatic conditions of the Qinba mountains, the indoor thermal comfort can be significantly improved and the building energy consumption can be effectively reduced when the energysaving structural wall with KPI porous brick as the main material and the heat transfer coefficient of the south wall of 0.9 W/(m^{2}·K) is selected.
5. Conclusion
(1)In winter, the indoor thermal environment of rural residential buildings in Qinba mountains is not ideal and the heating energy consumption is large. However, this area is relatively rich in solar energy resources and has the geographical advantage to develop solar energy buildings.(2)According to the climatic conditions of Qinba mountains and the living habits of local rural residents, it is suggested that the indoor thermal comfort temperature of rural residential buildings in winter should not be lower than 14°C.(3)When the internal surface temperature of the external walls in different orientations are equally based on the design principle of uniform radiation field, the heat transfer coefficient of the east wall, the west wall, and the north wall of the local rural residential buildings are 1.13 times, 1.06 times, and 1.14 times of the south wall heat transfer coefficient, respectively.(4)The indoor temperature of the buildings with the energysaving wall is significantly improved. When the main material of the energysaving wall is KP1 porous brick and the heat transfer coefficient of the south wall is 0.9 W/(m^{2}·K), the average indoor temperature of the master bedroom and secondary bedroom on the first floor increases by 14.3% and 15.3%, respectively, reaching 12.0°C and 11.3°C. Besides, the average indoor temperature of the master bedroom and secondary bedroom on the second floor increases by 14.0% and 14.9% respectively, reaching 11.4°C and 10.8°C.
Data Availability
All the data and models used to support the findings of this study are included within the paper.
Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was supported by the National Natural Science Foundation of China (Grant no. 51678483).
References
 S. P. Corgnati, M. Filippi, and S. Viazzo, “Perception of the thermal environment in high school and university classrooms: subjective preferences and thermal comfort,” Building and Environment, vol. 42, no. 2, pp. 951–959, 2007. View at: Publisher Site  Google Scholar
 D. G. Valeria, D. P. Osvaldo, and D. C. Michele, “Indoor environmental quality and pupil perception in Italian primary schools,” Building and Environment, vol. 56, no. 335, pp. 335–345, 2012. View at: Publisher Site  Google Scholar
 N. E. Klepeis, W. C. Nelson, W. R. Ott et al., “The nationalhuman activity pattern survey (NHAPS): a resource for assessing exposure to environment pollutants,” Journal of Exposure Analysis and Environment Epidemiology, vol. 11, no. 335, pp. 231–252, 2001. View at: Publisher Site  Google Scholar
 C. P. Le, N. Bonvallot, P. Glorennecet et al., “Indoor environment and children’s health: recent developments in chemical, biological, physical and social aspects,” International Journal of Hygiene and Environmental Health, vol. 215, no. 1, pp. 1–18, 2011. View at: Publisher Site  Google Scholar
 F. Wu, D. Jacobs, C. Mitchell, D. Miller, and M. H. Karol, “Improving indoor environmental quality for public health: impediments and policy recommendations,” Environmental Health Perspectives, vol. 115, no. 6, pp. 953–957, 2007. View at: Publisher Site  Google Scholar
 O. Hänninen, R. Sorjamaa, P. Lipponen, J. Cyrys, T. Lanki, and J. Pekkanen, “Aerosolbased modelling of infiltration of ambient PM2.5 and evaluation against populationbased measurements in homes in Helsinki, Finland,” Journal of Aerosol Science, vol. 66, no. 1, pp. 111–122, 2013. View at: Publisher Site  Google Scholar
 R. Pacheco, J. Ordóñez, and G. Martínez, “Energy efficient design of building: a review,” Renewable and Sustainable Energy Reviews, vol. 16, no. 6, pp. 3559–3573, 2012. View at: Publisher Site  Google Scholar
 W. He, Q. Y. Jiang, J. Ji, and W. Wei, “A study on thermal performance, thermal comfort in sleeping environment and solar energy contribution of solar Chinese Kang,” Energy and Buildings, vol. 58, no. 3, pp. 66–75, 2013. View at: Publisher Site  Google Scholar
 D. Song and Y.J. Choi, “Effect of building regulation on energy consumption in residential buildings in Korea,” Renewable and Sustainable Energy Reviews, vol. 16, no. 1, pp. 1074–1081, 2012. View at: Publisher Site  Google Scholar
 G. Liu, S. Zheng, P. Xu, and T. Zhuang, “An ANPSWOT approach for ESCOs industry strategies in Chinese building sectors,” Renewable and Sustainable Energy Reviews, vol. 93, no. 1, pp. 90–99, 2018. View at: Publisher Site  Google Scholar
 T. Huo, H. Ren, X. Zhang et al., “China's energy consumption in the building sector: a Statistical YearbookEnergy Balance Sheet based splitting method,” Journal of Cleaner Production, vol. 185, no. 6, pp. 665–679, 2018. View at: Publisher Site  Google Scholar
 Y. Zhang, C.Q. He, B.J. Tang, and Y.M. Wei, “China's energy consumption in the building sector: a life cycle approach,” Energy and Buildings, vol. 94, no. 5, pp. 240–251, 2015. View at: Publisher Site  Google Scholar
 M. Evans, S. Yu, B. Song, Q. Deng, J. Liu, and A. Delgado, “Building energy efficiency in rural China,” Energy Policy, vol. 64, no. 1, pp. 243–251, 2014. View at: Publisher Site  Google Scholar
 M. Shan, P. Wang, J. Li, G. Yue, and X. Yang, “Energy and environment in Chinese rural buildings: situations, challenges, and intervention strategies,” Building and Environment, vol. 91, no. 9, pp. 271–282, 2015. View at: Publisher Site  Google Scholar
 B.j. He, L. Yang, M. Ye, B. Mou, and Y. Zhou, “Overview of rural building energy efficiency in China,” Energy Policy, vol. 69, no. 4, pp. 385–396, 2014. View at: Publisher Site  Google Scholar
 World Bank, China 2030: Building a Modern, Harmonious, and Creative High Income Society, World Bank Publications, Washington, DC, USA, 2013.
 M. Ma, W. Cai, and Y. Wu, “China act on the energy efficiency of civil buildings (2008): a decade review,” Science of The Total Environment, vol. 651, no. 2, pp. 42–60, 2018. View at: Publisher Site  Google Scholar
 L. Jiang, X. P. Chen, and B. Xue, “Features, driving forces and transition of the household energy consumption in China: a review,” Sustainability, vol. 11, no. 4, pp. 1–20, 2019. View at: Publisher Site  Google Scholar
 Y. Y. Zhu, X. N. Fan, G. C. Sang, Q. Zhao, and X. L. Cui, “Study on indoor thermal environment of rural dwellings in Northwest China in winter,” Journal of Physics: Conference Series, vol. 1176, no. 4, pp. 1–9, 2019. View at: Publisher Site  Google Scholar
 S. Mirrahimi, M. F. Mohamed, L. C. Haw, N. L. N. Ibrahim, W. F. M. Yusoff, and A. Aflaki, “The effect of building envelope on the thermal comfort and energy saving for highrise buildings in hothumid climate,” Renewable and Sustainable Energy Reviews, vol. 53, no. 1, pp. 1508–1519, 2016. View at: Publisher Site  Google Scholar
 N. Prescott, “Energy saving through building insulation and airtightness,” Proceedings of the Institution of Civil Engineers  Energy, vol. 161, no. 2, pp. 51–55, 2008. View at: Publisher Site  Google Scholar
 L. Zhang and Q. F. Shi, “The comparison of energy saving characteristic of traditional residential envelope structure in southeastern shanxi,” Applied Mechanics and Materials, vol. 409410, no. 1, pp. 624–629, 2013. View at: Publisher Site  Google Scholar
 S. Lu, X. Tang, L. Ji, and D. Tu, “Research on energysaving optimization for the performance parameters of ruralbuilding shape and envelope by TRNSYSGenOpt in hot summer and cold winter zone of China,” Sustainability, vol. 9, no. 2, p. 294, 2017. View at: Publisher Site  Google Scholar
 Y. Z. Tian and Y. Yu, “Analysis of anshan existing residential building exterior wall energy saving reconstruction,” Advanced Materials Research, vol. 10041005, no. 1, pp. 1565–1569, 2014. View at: Publisher Site  Google Scholar
 X. N. Fan, Y. Y. Zhu, and G. C. Sang, “Optimization design of enclosure structure of solar energy building in Northwest China based on indoor zoning,” Advances in Engineering Research, vol. 143, no. 1, pp. 30–40, 2017. View at: Google Scholar
 J. Zhang and S. F. Dong, “Comprehensive measures should be taken for exterior window thermal insulating energysaving,” Applied Mechanics and Materials, vol. 361363, no. 1, pp. 327–330, 2013. View at: Publisher Site  Google Scholar
 R. Broun, H. Babaizadeh, A. Zakersalehi, and G. Menzies, “Integrated life cycle energy and greenhouse gas analysis of exterior wall systems for residential buildings,” Sustainability, vol. 6, no. 12, pp. 8592–8603, 2014. View at: Publisher Site  Google Scholar
 X. Wang, N. Wang, X. Liu, and R. Shi, “Energysaving analysis for the modern wing of the art institute of chicago and green city strategies,” Renewable and Sustainable Energy Reviews, vol. 73, no. 6, pp. 714–729, 2017. View at: Publisher Site  Google Scholar
 T. Nakaya, M. Yamasaki, and Y. Sasaki, “Influence of wall composition on thermal environment of wooden houses,” Journal of Wood Science, vol. 60, no. 2, pp. 117–126, 2014. View at: Publisher Site  Google Scholar
 B. Yu, J. Hou, W. He et al., “Study on a highperformance photocatalyticTrombe wall system for space heating and air purification,” Applied Energy, vol. 226, no. 9, pp. 365–380, 2018. View at: Publisher Site  Google Scholar
 X. L. Cui, Y. K. Zhang, G. C. Sang et al., “Coupling effect of spacearrangement and wall thermal resistance on indoor thermal environment of passive solar singlefamily building in tibet,” Applied Sciences, vol. 9, no. 17, pp. 1–19, 2019. View at: Publisher Site  Google Scholar
 JGJ2695, Energy Conservation Design Standard for Heating Residential Buildings, China Architecture Building Press, Beijing, China, 1996.
 J. I. Jiménez and Y. Castro, “Solar radiation on sloping surfaces with different orientations in granada, Spain,” Solar Energy, vol. 28, no. 3, pp. 257–262, 1982. View at: Publisher Site  Google Scholar
 M. Alabiso, F. Parrini, and R. Sidri, “Hourly solar radiation on Sloping planes with different orientations from measured daily radiation on the horizontal surfaces,” Environmental Software, vol. 4, no. 1, pp. 35–40, 1989. View at: Publisher Site  Google Scholar
 Y. Y. Zhu, Y. F. Liu, and J. P. Liu, “Thermal design of energy saving buildings in rural areas of northwest under the indoor uniform radiation filed,” Acta Energiae Solaris Sinica, vol. 32, no. 7, pp. 1034–1039, 2011. View at: Google Scholar
 L. Zhang, X.H. Liu, and Y. Jiang, “Experimental evaluation of a suspended metal ceiling radiant panel with inclined fins,” Energy and Buildings, vol. 62, no. 7, pp. 522–529, 2013. View at: Publisher Site  Google Scholar
 Y. Y. Zhu, Q. Zhang, J. P. Liu, and Q. Zhao, “Research on thermal stability of configuration of residential building energy saving exterior wall in Northwestern areas,” Journal of Xi’an University of Technology, vol. 27, no. 1, pp. 46–50, 2011. View at: Google Scholar
 Y. Q. Cui, N. H. Sun, H. B. Cai, and S. M. Li, “Indoor temperature improvement and energysaving renovations in rural houses of China’s cold region—a case study of Shandong Province,” Energies, vol. 13, no. 4, pp. 1–26, 2020. View at: Publisher Site  Google Scholar
 S. Gou, Z. Li, Q. Zhao, V. M. Nik, and J.L. Scartezzini, “Climate responsive strategies of traditional dwellings located in an ancient village in hot summer and cold winter region of China,” Building and Environment, vol. 86, no. 4, pp. 151–165, 2015. View at: Publisher Site  Google Scholar
 C. Mitterer, H. M. Künzel, S. Herkel, and A. Holm, “Optimizing energy efficiency and occupant comfort with climate specific design of the building,” Frontiers of Architectural Research, vol. 1, no. 3, pp. 229–235, 2012. View at: Publisher Site  Google Scholar
 C. Deng, H. Bai, X. Ma, T. Zhao, S. Gao, and X. Huang, “Spatiotemporal differences in the climatic growing season in the Qinling Mountains of China under the influence of global warming from 1964 to 2015,” Theoretical and Applied Climatology, vol. 138, no. 34, pp. 1899–1911, 2019. View at: Publisher Site  Google Scholar
 J. Li, X. Feng, J. Yin, and F. Chen, “Change analysis of spring vegetation greenup date in Qinba mountains under the support of spatiotemporal data cube,” Journal of Sensors, vol. 2020, no. 1, Article ID 6413654, pp. 1–12, 2020. View at: Publisher Site  Google Scholar
 Y. Y. Zhu, J. Z. Zhao, G. C. Sang, and Q. Zhao, “Thermal performance and energysaving construction for rural architecture in Qinba Mountain Areas,” Journal of Ningxia University (Natural Science Edition), vol. 38, no. 4, pp. 426–12438, 2017. View at: Google Scholar
 C. H. Deng, H. Y. Bai, S. Cao, T. Zhao, and X. P. Ma, “Differences and variations in the elevationdependent climatic growing season of the northern and southern slopes of the Qinling Mountains of China from 19852015,” Theoretical and Applied Climatology, vol. 137, no. 12, pp. 1159–1169, 2018. View at: Publisher Site  Google Scholar
 Q. Y. Zhang and H. Joe, Chinese Standard Weather Databank for Buildings, Machine Industry Press, Beijing, China, 2014.
 JGJ/T 267–2012, Technical Code for Passive Solar Buildings, China Architecture Building Press, Beijing, China, 2012.
 L. Zhao and J. Xu, “Architectural design of residential buildings in Qinba mountains based on knowledge mapping,” Case Studies in Thermal Engineering, vol. 14, no. 9, pp. 1–7, 2019. View at: Publisher Site  Google Scholar
 M. M. Yu, Y. Y. Ren, S. H. Liu, and D. M. Ma, “Southern Shanxi Ankang region residence geographic features and its inheritance and development,” Architecture and Culture, vol. 132, no. 3, pp. 109–111, 2015. View at: Google Scholar
 L. W. Dong, “Climate suitability analysis of villages and towns in Southern Shanxi based on climate consultant—taking the Ankang area as an example,” Architecture and Culture, vol. 180, no. 3, pp. 105106, 2013. View at: Google Scholar
 F. C. Houghtona and C. P. Yaglou, “Determining equal comfort lines,” Journal of the American Society of Heating and Ventilating Engineers, vol. 29, no. 1, pp. 165–176, 1923. View at: Google Scholar
 J. Liu, L. Wang, Y. Yoshino, and Y. Liu, “The thermal mechanism of warm in winter and cool in summer in China traditional vernacular dwellings,” Building and Environment, vol. 46, no. 8, pp. 1709–1715, 2011. View at: Publisher Site  Google Scholar
 J. P. Liu and Y. Liu, “Thermal design of a zero energy cavedwelling solar house,” Acta Energiae Solaris Sinica, vol. 20, no. 3, pp. 302–310, 1999. View at: Google Scholar
 T. Galloway, Solar House: A Guide for the Solar Designer, Architectural Press, Paris, England, 2004.
 D. A. Mcintyre, Indoor Climate, Applied Science Publisher, London, UK, 1980.
 ANSI/ASHRAE Standard 552013, Thermal Environmental Conditions for Human Occupancy, American Society of Heating, Refrigerating and Air conditioning Engineers, Inc., Atlanta, GA, USA, 2013.
 P. O. Fanger, Thermal Comfort, Danish Technical, Copenhagen, Denmark, 1970.
 Q. S. Yan and Q. Z. Zhao, Building Thermal Process, China Architecture & Building Press, Beijing, China, 1986.
 Y. Z. Li, Handbook of Thermal Design for Passive Solar House, Shanghai Scientific & Technical Publishers, Shanghai, China, 1989.
 Q. Wang, Study on Low Energy Consumption Structure System of Qinba Mountainous Building, Xi’an University of Technology, Xi’an, China, 2018.
 GB/T508242013, Design Standard for Energy Efficiency of Rural Residential Buildings, China Architecture Building Press, Beijing, China, 1993.
 06J6071, National Building Standard Design Atlas, China Planning Press, Beijing, China, 2006.
 H. Gao and G. Wei, “Numerical solution of potential problems using radial basis reproducing kernel particle method,” Results in Physics, vol. 13, no. 2, Article ID 102122, 2019. View at: Publisher Site  Google Scholar
 H. Gao and G. Wei, “Complex variable meshless manifold method for transient heat conduction problems,” International Journal of Applied Mechanics, vol. 9, no. 5, Article ID 1750067, pp. 1–20, 2017. View at: Publisher Site  Google Scholar
 Z. J. Meng, H. Cheng, L. D. Ma et al., “The dimension splitting elementfree Galerkin method for 3D transient heat conduction problems,” Physics, Mechanics & Astronomy, vol. 62, no. 4, Article ID 040711, pp. 1–12, 2019. View at: Publisher Site  Google Scholar
 J.F. Wang and Y.M. Cheng, “A new complex variable meshless method for transient heat conduction problems,” Chinese Physics B, vol. 21, no. 12, Article ID 120206, pp. 1–9, 2012. View at: Publisher Site  Google Scholar
 G. Zhang, F. Zhang, and H. Ren, “Accuracy analysis of interpolating elementfree galerkin (IEFG) method in solving transient heat conduction problems,” International Journal of Applied Mechanics, vol. 8, no. 6, Article ID 1650078, pp. 1–18, 2016. View at: Publisher Site  Google Scholar
 M. J. Peng and Y. M. Cheng, “A boundary elementfree Galerkin (IEFG) for twodimensional potential problems,” Engineering Analysis with Boundary Elements, vol. 36, no. 5, pp. 77–82, 2009. View at: Publisher Site  Google Scholar
 S. Mukherjee and Y. Liu, “The Boundary element method,” International Journal of Computational Methods, vol. 10, no. 6, Article ID 1350037, pp. 1–91, 2013. View at: Publisher Site  Google Scholar
 Z. Q. Zhang, Z. J. Wang, and L. M. Lian, “Numercal SsimuIation for indoor thermal environment in residential buildings,” Building Energy and Environment, vol. 23, no. 6, pp. 88–92, 2004. View at: Google Scholar
 K. Matsunawa, H. Iizuka, and S. Tanabe, “Development and application of an underfloor airconditioning system with improved outlets for a ‘smart' building in Tokyo,” Ashrae Transactions, vol. 101, no. 1, pp. 887–901, 1995. View at: Google Scholar
 GB 501762016, Code for Thermal Design of Civil Building, China Architecture Building Press, Beijing, China, 2016.
Copyright
Copyright © 2020 Qin Zhao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.