Research on Indoor Thermal Environment Analysis and Optimization Strategy of Rural Dwellings around Xi’an Based on PET Evaluation

: Rural dwellings are an important group of residential buildings in China. With the continuous development of rural construction in China, the contradiction between the pursuit of a simple material space and the villagers’ demand for living quality, especially the indoor thermal comfort of rural dwellings, has become increasingly prominent. Therefore, it is particularly important to study the optimization strategies of the indoor thermal environment in rural dwellings. Current research on optimizing the indoor thermal environment of rural dwellings mainly focuses on analyzing the impact of individual factors, such as the envelope structures, building constructions, and building technology applications, but there is a lack of strategy development based on the comprehensive evaluation. This study aims to analyze the combined effects of multiple design elements on the indoor thermal environment and propose a comprehensive optimization strategy for rural dwellings. This study selects the rural dwellings around Xi’an as an example and establishes a basic model of the rural dwellings around Xi’an through ﬁeld investigation and software simulation. Then, through univariate and compound-variable simulations, we analyze the inﬂuence of changes in passive architectural design indicators on the indoor physiological equivalent temperature (PET) of rural dwellings and obtain a comprehensive design indicator optimization strategy. This strategy can improve the indoor thermal comfort in winter and summer, especially in winter, achieving an average increase of 4.17 ◦ C in the winter PET value and an average decrease of 0.66 ◦ C in summer. This provides a reference for the design and renovation of rural dwellings in Xi’an and other rural areas in the cold regions of China.


Introduction 1.Overview of Rural Dwelling Development
Rural dwellings are the main body of rural settlements and an important group of residential buildings [1].In 2020, the national rural residential area accounted for approximately 61.23% of the total rural construction area in China [2].With China's urbanization and new rural construction, the living standards in rural areas have been continuously improved [3], and rural residential construction has achieved good results.In 2020, the area of new and rebuilt rural residential buildings totaled 756 million square meters, the per capita residential building area reached 34.3 square meters [2], and the per capita disposable income reached 17,131.5 yuan [4].However, the current rural construction process still focuses on the construction of the physical space [5], and pays less attention to the quality of the living space [6].The contradiction between this situation and the villagers' increasing pursuit of quality of life is becoming increasingly serious [7].Most of the current rural dwellings are houses self-built by villagers, accounting for approximately 89% [8], and a few are new houses planned and constructed by local governments and design units in the process of new rural construction [9].Due to the lack of design resources in the self-built Sustainability 2023, 15, 7889 2 of 25 houses of the villagers [10], the self-built houses pay more attention to the appearance and shape, while ignoring the comfort of the indoor environment [11].Designers often give priority to economy in the process of planning and design, and do not pay sufficient attention to the regional environment and indoor comfort of the planned houses [12].This has led to the neglect of thermal comfort factors in rural dwellings, especially the inability to scientifically grasp the changing relationship of the indoor thermal environment of rural dwellings under different design conditions, so that it is impossible to effectively formulate corresponding optimal design strategies for different conditions [13].This is an urgent problem that needs to be solved in the current rural housing construction context, and it is also particularly important for the comprehensive realization of a beautiful and livable rural revitalization strategy [14].

Research on the Indoor Thermal Environment of Rural Dwellings
Aiming at the neglect of the indoor thermal environment in the development of rural dwellings, related research is also gradually being carried out [15].China's research on residential thermal environment-related issues is mainly focused on urban dwellings.The research on the thermal environment of urban dwellings has made a great deal of progress, and a large number of optimization strategies have been proposed [16].Existing research mainly includes the following: 1. Optimization of the design of the spatial layout, material structure, and facilities of urban residential buildings [17,18].2. Analysis of the influence of factors such as the scale, floor area ratio, and greening conditions of urban residential areas on the indoor thermal environment of residential buildings [19,20].3. Research on passive low-energy houses and near-zero-energy houses in response to the national green building development plan [21,22].
With the development of new rural construction and the emergence of indoor thermal environment problems in rural dwellings, many scholars have also turned their research attention to rural areas, conceiving ideas and experimenting to apply the research and optimization methods of urban residences to rural residences [23].However, rural dwellings are quite different from urban dwellings in terms of building site selection, orientation, spatial layout, building materials, and structure selection due to the strong regional characteristics and the fact that most of them are built independently by villagers.This leads to the inapplicability of many existing research methods, and the relatively shallow development of related research [24].
Research on the indoor thermal environment of rural dwellings can be roughly divided into two aspects: analyzing the influence of housing design elements on the indoor thermal environment and analyzing the design characteristics suitable for the climate environment of the rural area [25].In terms of analyzing the influence of housing design elements on the indoor thermal environment, Shao et al. [26] conducted a study on the impact of parameterized changes in the thickness of insulation layers under multiple building envelope structures on the indoor thermal environment, energy-saving effects, and economic costs of rural dwellings throughout their lifecycle.They determined the optimal thickness of insulation layers for three different life spans of 10, 20, and 30 years.Chi et al. [27] compared the effects of different orientations and window-to-wall ratios on the indoor temperature and air velocity of rural dwellings, and determined the optimal window-towall ratio range for each orientation, with a span of 20 • .Zhang et al. [28] studied the effects of various types of attached sunrooms on the indoor temperature, energy efficiency, and cost of rural dwellings, and selected the most suitable sunroom structure for the local climate environment.Compared with dwellings without sunrooms, the optimized indoor temperature increased by 2.38 • C, and the energy-saving rate was 44.8%.Zhang et al. [29] conducted an experimental study on the indoor air temperature and CO 2 concentration of rural dwellings using six different mechanical ventilation methods, and determined the best air supply mode for improving the indoor thermal environment and air quality.
Overall, current research mainly focuses on the analysis of single elements, such as the envelope structures, building constructions, and building technology applications, but there is a lack of strategy development based on the comprehensive evaluation.
In terms of analyzing the design characteristics of regional rural dwellings, the research regions are mainly divided into the cold and severely cold regions of Northern China and the hot summer and cold winter regions of Southern China [30].Among them, the cold and severely cold regions not only occupy a large geographical area in China, but the main research focus is also on improving the indoor thermal environment in winter and optimizing heating energy consumption.This is closely related to China's development of energy-efficient buildings and green architecture, and therefore receives great attention [31].In the cold and severely cold regions, the northwest region is particularly affected due to its location deep inland and its high terrain, leading to arid climate conditions.Additionally, the low economic level limits its rural development, resulting in the poor thermal performance of residential buildings, and particularly severe indoor thermal environment problems in rural dwellings [32].Xi'an, as the most economically developed city in the northwest region, has seen rapid growth in rural housing construction in recent years, resulting in a widespread phenomenon of new construction and renovation in surrounding rural areas.In 2020, the construction area of rural housing in Xi'an was approximately 4.0733 million square meters, an increase of approximately 27.87% compared to 2019 [33].
In summary, current research on the analysis and optimization of the indoor and outdoor environments of rural dwellings has focused more on the analysis of single design elements, but lacks a comprehensive analysis of multiple design elements, such as the floor space, building area, building form, and roof form.The optimization strategies derived from these studies are often one-sided and cannot form a complete set of design optimization strategies [34].Therefore, further research is needed on how to identify and extract multiple design elements and analyze their degrees of influence on the indoor thermal environment of housing, in order to derive comprehensive optimization strategies for rural dwellings.In the research on the design characteristics of regional rural dwellings, the characteristics of rural dwellings in the northwest region and the surrounding rural areas of Xi'an have great research value.However, due to the constraints of economic development and other practical factors, relevant research is difficult to carry out, and research on optimization strategies for the indoor thermal environment is even scarcer [35].
Therefore, this study selected rural dwellings around Xi'an as the research object.By combining a field investigation and data analysis, the study analyzed the influence of design elements on the indoor thermal environment of houses, and proposed comprehensive design strategies that could optimize the indoor thermal environment, providing design references for future rural construction projects.

Current Situation of Rural Dwellings around Xi'an
The rural dwellings around Xi'an are mostly in the form of narrow courtyards of traditional Guanzhong houses.In terms of the spatial layout, the traditional homestead is a long and narrow rectangle with a length-to-width ratio of about 2:1 [36,37].The short side of the homestead is located along the village road.The courtyard is composed of a gatehouse near the courtyard gate, flanking houses on both sides, and a main house on the inner side [37].The main house is the primary living space and includes a living room and bedrooms.The interior of the main house has a three-bay layout, with the living room in the center and the bedrooms on both sides.The gatehouse and flanking houses are auxiliary rooms that include storerooms, kitchens, and washrooms [38].With the development of rural areas in Xi'an, the dwellings have roughly retained the original spatial layout, while the gatehouse and flanking houses have gradually simplified and merged [39].The common spatial form of existing dwellings is shown in Figure 1.In terms of spatial dimensions, the width of the existing residential homestead is mostly 9-10 m, and the length is about 20 m.The width of each bay of the main house is mostly 3 m, and some of the central bays are widened.The length of the main house is mostly 5-7 m.The width of the auxiliary rooms is mostly 3 m [37].
In terms of the building structure, traditional dwellings are made of raw soil and mostly have one floor.The roof of traditional dwellings is the pitched roof or semi-pitched roof unique to Guanzhong dwellings [40].With the development of rural areas, the newly built dwellings are all brick-concrete structures.The traditional raw soil buildings have gradually disappeared.Based on the brick-concrete structures, many existing dwellings have added cladding materials such as cement and ceramic tiles to the exterior walls [39].The number of floors in existing dwellings has increased to two, and the pitched roof and semi-pitched roof have been replaced by a flat roof [40].This has resulted in some differences in the types of rural dwellings, as shown in Table 1.The common building structure of existing dwellings is shown in Figure 1.There are also a large number of thermal environmental problems in rural dwellings in Xi'an.They are mainly reflected in the following: 1.Many traditional building practices that adapt to the local climate and can improve the indoor thermal environment are abandoned.2. Due to insufficient construction technology, the heat insulation effect of rural dwellings is poor.3.There is a lack of sustainable housing design, and there are a great deal of over-design phenomena in rural areas [41].In terms of the building structure, traditional dwellings are made of raw soil and mostly have one floor.The roof of traditional dwellings is the pitched roof or semi-pitched roof unique to Guanzhong dwellings [40].With the development of rural areas, the newly built dwellings are all brick-concrete structures.The traditional raw soil buildings have gradually disappeared.Based on the brick-concrete structures, many existing dwellings have added cladding materials such as cement and ceramic tiles to the exterior walls [39].The number of floors in existing dwellings has increased to two, and the pitched roof and semi-pitched roof have been replaced by a flat roof [40].This has resulted in some differences in the types of rural dwellings, as shown in Table 1.The common building structure of existing dwellings is shown in Figure 1.
Table 1.Types of rural dwellings around Xi'an (self-made by the author based on [37,39]).There are also a large number of thermal environmental problems in rural dwellings in Xi'an.They are mainly reflected in the following:

Structure
1.Many traditional building practices that adapt to the local climate and can improve the indoor thermal environment are abandoned.2. Due to insufficient construction technology, the heat insulation effect of rural dwellings is poor.3.There is a lack of sustainable housing design, and there are a great deal of over-design phenomena in rural areas [41].
Due to the economic level and building conditions in rural areas around Xi'an, the heating and cooling equipment level of rural houses is relatively low [11].For example, during winter, there is no centralized heating system for rural houses as there is in urban areas, and many villagers still use stoves as heating equipment.Due to the limitations of personal economic levels and mindsets, villagers usually try to minimize the use of heating and cooling equipment to save costs [13].In this case, it is difficult to use heating and cooling equipment to solve the thermal environment problems of rural dwellings in Xi'an.Therefore, it is very urgent to improve the indoor thermal environment by optimizing the existing design features of rural dwellings [42].

Evaluation Methods for the Indoor Thermal Environment
The commonly used methods for evaluating thermal comfort include predicted mean vote (PMV), physiological equivalent temperature (PET) and standard effective temperature (SET) [43].Among them, PET is a human biometeorological parameter that describes individual thermal perception.Its definition is the air temperature at which, in a typical indoor setting (without wind and solar radiation), the heat budget of the human body is balanced with the same core and skin temperature as under the complex outdoor conditions to be assessed [44].In this study, the indoor thermal environment of rural residences is greatly influenced by outdoor conditions due to the lower level and less frequent use of heating and cooling equipment.In this case, PMV can only calculate the integrated index of human thermal perception of the indoor environment, rather than provide a real description of the human thermal status.SET can only express the influence of the indoor thermal environment on human heat dissipation and cannot directly represent the thermal sensation of individuals.PET, on the other hand, can directly represent the core and skin temperatures of the human body, providing a more realistic description of human thermal perception and a more intuitive reflection of the indoor thermal environment [45,46].Therefore, this study is based on the PET evaluation method to analyze the indoor thermal environment of rural residences.

Field Investigation Method
A number of rural dwellings around Xi'an were selected for the field investigation, and their architectural status and spatial design indicators were measured, so as to obtain the spatial layout, design mode, and design parameters of the rural dwellings around Xi'an, which were then used for simulation modeling.

Selection of Investigation Objects
Several villages around Xi'an were selected for investigation based on factors such as location and development situation.All the dwellings in these villages were examined.According to the classification of existing rural dwellings around Xi'an in Table 1, the specific types of rural dwellings in terms of building structure, number of floors, roof type, and external wall cladding materials were identified through the examination.Several dwellings covering all types were selected as the investigation objects.

Investigation Method 1. Current Status Record of Architectural Design
We recorded the spatial layout, orientation, structure, door and window material, wall material, floor material, floor material, roof type, roof cladding material, and other building status indicators of the selected dwellings.

Mapping of Architectural Space Dimensions
We measured the dimensions of courtyards and indoor rooms, the height of buildings, the thickness of walls and floors, the dimensions and thickness of doors and windows, etc., and drew plane sketches.

Software Simulation Method
Modeling and software simulation were carried out based on the measured results to determine the impact of various design indicators' changes on the indoor thermal environment, and a comprehensive optimization strategy was established based on the optimal results.

Selection of Simulation Software
Common simulation software can be divided into two parts: the simulation engine that calculates indoor thermal environment indicators and the software platform for building modeling based on the simulation engine.In terms of simulation engines, commonly used engines include EnergyPlus, IESVE, and DeST [47].Among them, IESVE and DeST can only handle linear and steady-state systems, that is, simulating the thermal performance of buildings under indoor thermal environment conditions that do not vary with time.EnergyPlus can handle various transient systems where multiple thermal parameters change over time, and simulate indoor thermal environments under the interference of equipment, personnel, and other factors [48].DeST only supports the CAD modeling platform, IESVE supports SketchUp and Revit platforms, while EnergyPlus can support multiple platforms, such as CAD, SketchUp, Revit, and Rhino [49].A summary of the differences between the simulation engines is shown in Table 2. Therefore, EnergyPlus was selected as the simulation engine in this study.In terms of software platforms, compared to platforms such as SketchUp and Revit, Rhino can accurately create building design and thermal models through parametric modeling, and achieve real-time coordination between modeling and EnergyPlus engine calculations [47].Therefore, this study selected Rhino as the modeling platform and used Rhino software based on the EnergyPlus engine for the simulation.

Simulation Method 1. Boundary Condition Settings
The simulation time was set to typical winter and summer days.The outdoor climate conditions were set according to the Xi'an meteorological data provided by EnergyPlus's climate database.The building thermal conditions, personnel object conditions, and PET calculation conditions were set based on the design status of rural housing obtained through the field investigation, in combination with relevant standards, such as "GB/T 50176-2016 Technical Standard for Thermal Design of Civil Buildings" and "DB/J 61/T91-2014 Energy-Saving Technical Standard for Rural Residential Buildings in Xi'an Area" [50,51].

Establishment of the Basic Model
The spatial layout and dimensions of the basic model were determined based on the results of the field investigation and the spatial layout of existing rural dwellings around Xi'an.The design values of each indicator of the basic model were determined based on the design status and dimension values of the investigation objects.

Univariate Simulation and Evaluation
Based on the extraction of the design elements of the measured dwellings, some design elements that have an influence on the indoor thermal environment were selected as simulation indicators for the univariate simulation.The indoor PET values corresponding to the basic model were simulated under the variation of a single indicator.The influence of each variable indicator on the indoor thermal environment was evaluated based on the change in PET compared to the basic value.The optimal solution for each variable indicator was obtained, and the design indicators were evaluated based on their influence, thus determining the indicators that need further research on the compound impact on PET under mutual interference.

Composite-Variable Simulation and Optimization Strategy
Based on the evaluation results, the indicators for the composite-variable simulation were obtained.The indicators were combined into multiple working conditions, and the indoor PET values corresponding to each working condition were simulated.Based on the change in PET compared to the basic value for each working condition, the compound working condition with the best PET optimization result was obtained, and thus the indoor thermal environment optimization strategy for the basic model was determined.
The flowchart of the software simulation method is shown in Figure 2.
Sustainability 2023, 15, x FOR PEER REVIEW 7 of 27 simulation indicators for the univariate simulation.The indoor PET values corresponding to the basic model were simulated under the variation of a single indicator.The influence of each variable indicator on the indoor thermal environment was evaluated based on the change in PET compared to the basic value.The optimal solution for each variable indicator was obtained, and the design indicators were evaluated based on their influence, thus determining the indicators that need further research on the compound impact on PET under mutual interference.

Composite-Variable Simulation and Optimization Strategy
Based on the evaluation results, the indicators for the composite-variable simulation were obtained.The indicators were combined into multiple working conditions, and the indoor PET values corresponding to each working condition were simulated.Based on the change in PET compared to the basic value for each working condition, the compound working condition with the best PET optimization result was obtained, and thus the indoor thermal environment optimization strategy for the basic model was determined.
The flowchart of the software simulation method is shown in Figure 2.

Investigation Objects
The rural areas in Xi'an are mainly under the jurisdiction of various towns in the city.In order to select the investigation objects, it is necessary to first select the towns around Xi'an [52].After screening the towns around Xi'an, Taiyigong Street was selected.Taiyigong Street is one of the nine provincial-level demonstration towns for rural revitalization in Xi'an, which has strong representativeness for the current development status of rural areas in Xi'an.Taiyigong Street is located along the south ring road of Xi'an, about 27.2 km away from the city center.There are 11 villages under Taiyigong Street.Among them, Shachang Village and Guanjia Village were selected as the rural investigation objects.
After the investigation, it was found that the construction conditions of dwellings in the two villages were roughly the same, and the spatial layout was consistent with the overall characteristics of rural dwellings around Xi'an.In terms of dwelling type, the dwellings in both villages were brick-concrete structures.The main differences included the number of floors, roof type, and exterior wall cladding materials.In addition, due to the current situation of village road and homestead planning, there were also differences in the orientation of dwellings in the two villages.The specific types are shown in Table 3.Typical rural dwellings covering each type were selected as the investigation objects.Five households were selected from Shachang Village and six households were selected from Guanjia Village, totaling eleven households as investigation objects.The distribution and numbering of each household are shown in Figures 3 and 4.

Investigation Objects
The rural areas in Xi'an are mainly under the jurisdiction of various towns in the city.In order to select the investigation objects, it is necessary to first select the towns around Xi'an [52].After screening the towns around Xi'an, Taiyigong Street was selected.Taiyigong Street is one of the nine provincial-level demonstration towns for rural revitalization in Xi'an, which has strong representativeness for the current development status of rural areas in Xi'an.Taiyigong Street is located along the south ring road of Xi'an, about 27.2 km away from the city center.There are 11 villages under Taiyigong Street.Among them, Shachang Village and Guanjia Village were selected as the rural investigation objects.
After the investigation, it was found that the construction conditions of dwellings in the two villages were roughly the same, and the spatial layout was consistent with the overall characteristics of rural dwellings around Xi'an.In terms of dwelling type, the dwellings in both villages were brick-concrete structures.The main differences included the number of floors, roof type, and exterior wall cladding materials.In addition, due to the current situation of village road and homestead planning, there were also differences in the orientation of dwellings in the two villages.The specific types are shown in Table 3.Typical rural dwellings covering each type were selected as the investigation objects.Five households were selected from Shachang Village and six households were selected from Guanjia Village, totaling eleven households as investigation objects.The distribution and numbering of each household are shown in Figures 3 and 4.

Architectural Design and Construction Element Measurement
According to the measured results, we sorted out the architectural form features and material structures of each household, as shown in Table A1 in Appendix A. All households were brick-concrete structures with no enclosure insulation layer.The main differences included orientations, exterior wall cladding materials, roof types, roof cladding materials, and window materials.Except for houses 7 and 8, the exterior wall cladding material of all households was ceramic tile.House 7 was cement, while house 8 had no exterior wall cladding material.The roof types and cladding materials can be simplified into three forms: traditional tile pitched roof, cement-covered flat roof, and color steel roofs added for shading on cement flat roofs.Houses 5, 8, and 11 had tile pitched roofs, while houses 3 and 7 had color steel roofs installed.Window materials were divided into window frame materials and window glass types.Window frame materials included

Architectural Design and Construction Element Measurement
According to the measured results, we sorted out the architectural form features and material structures of each household, as shown in Table A1 in Appendix A. All households were brick-concrete structures with no enclosure insulation layer.The main differences included orientations, exterior wall cladding materials, roof types, roof cladding materials, and window materials.Except for houses 7 and 8, the exterior wall cladding material of all households was ceramic tile.House 7 was cement, while house 8 had no exterior wall cladding material.The roof types and cladding materials can be simplified into three forms: traditional tile pitched roof, cement-covered flat roof, and color steel roofs added for shading on cement flat roofs.Houses 5, 8, and 11 had tile pitched roofs, while houses 3 and 7 had color steel roofs installed.Window materials were divided into window frame materials and window glass types.Window frame materials included

Architectural Design and Construction Element Measurement
According to the measured results, we sorted out the architectural form features and material structures of each household, as shown in Table A1 in Appendix A. All households were brick-concrete structures with no enclosure insulation layer.The main differences included orientations, exterior wall cladding materials, roof types, roof cladding materials, and window materials.Except for houses 7 and 8, the exterior wall cladding material of all households was ceramic tile.House 7 was cement, while house 8 had no exterior wall cladding material.The roof types and cladding materials can be simplified into three forms: traditional tile pitched roof, cement-covered flat roof, and color steel roofs added for shading on cement flat roofs.Houses 5, 8, and 11 had tile pitched roofs, while houses 3 and 7 had color steel roofs installed.Window materials were divided into window frame materials and window glass types.Window frame materials included wood and aluminum alloy, and window glass types were all single layers.Types of functional room were divided into five categories: living room, bedroom, kitchen, washroom, and storeroom.The numbers of the first four types were similar, while the number of storerooms greatly varied.

Building Space Dimensions
We calculated the dimensional data of each household's building space and the window-to-wall ratio (WWR) of the main living room based on the measurement results, as shown in Table A2 in Appendix A. Except for houses 3, 5, 8, and 11, the length of the main living rooms was the same as the width of the homestead.The width and floor height of the main living rooms were generally similar.Most of the houses had a two-story main living room, while houses 4, 6, 8, and 10 had a one-story main living room.The windows of the main living rooms in each household faced the front and back of the building, and there were no windows on the sides.The window-to-wall ratio of each household significantly varied.

Building Plans
Based on the measurement results, we drew the floor plans of each household, as shown in Table 4.According to the floor plan, it can be seen that the spatial layout of each household generally conformed to the spatial layout of existing rural dwellings around Xi'an.Except for house 3, the main living room of each household was located near the short side of the homestead.The auxiliary rooms of each household consisted of one or more rooms and were all located near the entrance of the courtyard.

Extraction of Rural Dwellings' Design Elements
Based on the results of the investigation into the current state of rural housing design, we identified the design elements that were decisive for the design of the main living rooms and extracted the numerical or typological ranges of each element, as shown in Table 5.This will be used for the establishment of the basic model and the setting of simulation indicators in the subsequent simulation process.The simulation time was set to the winter solstice (21 December) and the summer solstice (21 June).Based on the meteorological data of Xi'an in the EnergyPlus climate database, we obtained the parameters for outdoor climate conditions used in the simulation.The main parameters' average values during the simulation period are shown in Table 6.The materials involved in the actual construction process and subsequent optimization process can be divided into structural materials and window materials.Their corresponding thermal parameters are shown in Tables 7 and 8. To simulate the thermal environment without equipment use, the effects of heating, cooling, and other equipment on the indoor thermal environment were excluded in the modeling, as well as the impacts of other dwellings and facilities outdoors.According to the actual situation, the model assumes no shading devices for windows in winter and summer, and the average air change rate was 0.5 times/h.The investigation found that the elderly were the main occupants of the investigation objects.Therefore, the simulation process took a 60-year-old male as the personnel object.The personnel object's height was set to 175 cm, weight was 75 kg, and activity level was sitting.The clothing level was set according to the clothing worn by elderly residents during the investigation.The specific parameters are shown in Table 9.The calculation method for PET simulation was to set a calculation grid that covers the entire area of each floor at a certain height from the floor level in the main living room.The length and width of each cell on the grid were both 0.5 m.The height from the floor level was set to 0.6 m, which is the measurement point height for sitting activity.The calculation period was 24 h, and the PET value of each cell on the grid was calculated every hour.The average PET value of all cells on the grid was taken as the indoor PET simulation result for every hour.

Basic Model Establishment
Based on the field investigation results, we established the basic model of the main living room of rural dwellings.Firstly, based on the measured plans, we simplified the spatial layout of the investigation objects.Combined with the spatial layout of existing rural dwellings around Xi'an, we developed the outdoor courtyard and indoor spatial layout of the basic model, as shown in Figure 5. Based on the spatial dimensions of existing rural dwellings around Xi'an and the measurement results of the investigation objects, we set the dimension values of the basic model that were approximate to the actual situation.The width of the homestead was 10 m, and the length was 18 m.The length of the main living space was 10 m, and the width was 6 m.The width of the central bay of the main living room was 4 m, and the width of bays on both sides was 3 m.The width of the auxiliary room was 3 m, and the length was 6 m.The detailed dimension values are shown in the dimension labeling of Figure 5.The other indicator design values of the basic model were also determined by values and types that were approximate to the actual situation, and the detailed values are shown in Table 10.The digital model of the basic model is shown in Figure 6.  10.The digital model of the basic model is shown in Figure 6.

Univariate Simulation Process
We simulated the change in the PET value under the change in a single indicator, and then analyzed the influence of each indicator on the indoor thermal environment.The range of indicators that need to be simulated and compared is shown in Table 11.We simulated the indoor PET values corresponding to the changes in various indicators of the basic model during the simulation period, as shown in Table A3 in Appendix A. According to the simulation results, the indicators with better optimization results in winter were as follows: south orientation, 3.0 m floor height, inner and outer insulation   11.We simulated the indoor PET values corresponding to the changes in various indicators of the basic model during the simulation period, as shown in Table A3 in Appendix A. According to the simulation results, the indicators with better optimization results in winter were as follows: south orientation, 3.0 m floor height, inner and outer insulation layers, and a 40% south-facing window-to-wall ratio.The indicators with better optimization results in summer were as follows: south orientation, color steel roof, and a 30% south-facing window-to-wall ratio.

Evaluation of Univariate Simulation Results
For the simulation results, we took the average value of the PET data of a single indicator after a change in winter as T wi ( • C), the average value of the basic PET data of the corresponding indicator as T wb ( • C), and the optimal value, ∆T w ( • C), of the indicator change relative to the basic indicator can be obtained, where: In the same way, taking the average value of the PET data of a single indicator after a change in summer as T si ( • C), and the average value of the basic PET data of the corresponding indicator as T sb ( • C), the optimal value, ∆T s ( • C), of the indicator change relative to the basic indicator can be obtained, where: It is better to increase the PET value in winter, and to decrease it in summer.Therefore, for the optimal value of a certain indicator, its ∆T w should be a positive number, and ∆T s should be a negative number.From this, it can be seen that the value of its comprehensive optimal value, ∆T ( • C), should be: From this, the PET optimization results can be obtained, as shown in Table 12.From analyzing and evaluating the single-variable simulation results, it can be concluded that the optimal solution variables for each indicator in the PET optimization results were as follows: south direction, 3.0 m floor height, color steel roof, cement exterior wall cladding material, external insulation, 40%S10%N window-to-wall ratio, and aluminum alloy double-glazed window material.Through the evaluation of the univariate simulation results, it can be seen that among the various indicators, the south direction, the addition of a color steel roof, and the addition of an external insulation layer had a clear improvement effect, the effect of changing the cladding material of the exterior wall and the window material was not obvious, and the effect of changing the floor height and the window-to-wall ratio requires further study.Therefore, in the original basic model, the orientation was determined to be south, the colored steel roof and external insulation layer were added, and the simulation indicators of the composite simulation were obtained as shown in Table 13.According to the simulation results, we obtained the design indicators for the optimization strategy.Compared to the original basic model, the improved indicators of the optimization strategy included: building length of 12.6 m, width of 7.8 m, floor height of 2.7 m, a south window-to-wall ratio of 10%, a north window-to-wall ratio of 10%, a color steel roof, 60 mm-thick external insulation, a cement exterior wall cladding material, and an aluminum alloy double glass window material.As a result of the changes in length and width, the building area and room area of the optimized building also changed.A comparison between the detailed indicator values of the optimized building and the original values is shown in Table 17.Compared with the basic model before simulation, this strategy increased the average PET value by 4.17 • C in winter and decreased it by 0.66 • C in summer.Xi'an belongs to the cold regions, and the optimization results of the PET value in winter were better, which proved that the optimization strategy is more feasible.

Comparison with Previous Studies
Through our research, we know that the typical rural residential space around Xi'an is three longitudinal rooms symmetrically arranged side-by-side, and the building form is a two-story, flat-roofed, brick-concrete building.The indicators that had a greater impact on the indoor thermal environment of rural dwellings around Xi'an were the building shape, enclosure structure, and window-to-wall ratio.These are the same conclusions as obtained in previous studies [27,39,53,54].After this study, it can be further found that the indoor thermal environment performance is better when the length and width of the house are larger, the floor height is smaller, the insulation layer of the enclosure structure is thicker, and the windows facing the south and north are smaller.

Limitations of the Study
This study has some limitations in the field investigation and simulation process.In the field investigation part, due to the limitation of the number of people and time, the research objects were limited, and only 11 dwellings that could cover the rural dwelling types around Xi'an City were selected.It is impossible to obtain generalized measurement data through small-sample surveys.In terms of the simulation, there were some factors that affected the accuracy of the simulation results.For example, the external environmental parameters referred to in the simulation were the average environmental parameters of Xi'an City, which are somewhat different from the microclimate environment in the actual village.Hence, the simulation process was unable to express the influence of other users of the house on the simulated experience objects, etc.The limitations of the existing process will be improved in the subsequent research.Due to the current status of rural dwellings around Xi'an, this study mainly analyzed the influence of design elements on the indoor thermal environment, without calculating the building energy consumption.In the subsequent research, energy consumption can also be included as an evaluation indicator to derive multi-objective optimization strategies.

Economic Issues of Optimization Strategies
The optimization strategy proposed in this study has undergone a large number of modifications compared with the initial conditions of the basic model to ensure the maximum PET optimization results.In the actual construction process, due to the limited economic level in rural areas and the fact that most rural houses are self-built by villagers, there are many economic problems in the optimization and renovation of houses.The optimization strategy should consider the construction costs required for the transformation of various indicators, compare this with the PET optimization results, and choose the most cost-effective and economical optimization strategy according to the actual situation of the user.Subsequent research will also calculate the required material costs for all working conditions in the composite simulation and obtain the cost and the variation value of the basic model.At the same time, a coupling analysis will be carried out with the simulated PET value to obtain the coupling results and cost-effectiveness of the PET optimization value and the cost of each working condition and establish a dataset for obtaining a multiobjective optimization strategy.

Possible Future Research Directions
This study analyzed the influence of various design indicators on the indoor thermal environment of rural dwellings around Xi'an through field investigation and simulation.Based on the results, we developed a comprehensive optimization design strategy.In future studies, it is possible to actually construct the optimized rural dwellings based on the optimization strategy through practical methods.Through the post-occupancy evaluation of constructed dwellings, the actual influence and effectiveness of the optimization strategy on the indoor thermal environment of rural dwellings can be verified.Furthermore, the indicator values of the optimization strategy can serve as a reference for setting or updating local energy-saving design standards in the future.For example, as mentioned in the "Energy-Saving Technical Standard for Rural Residential Buildings in Xi'an Area", the length of a new rural residence should not exceed 18 m, and the floor height should not exceed 3 m [51].The length of the optimization strategy is 12.6 m and the floor height is 2.7 m. Additionally, the thickness of the insulation layer for the external insulation structure in the standard is 20 mm, while the thickness of the optimization strategy is 60 mm [51].These values can provide supplementary optimized values for the implementation of relevant standards in the future.
In the current composition of residents in rural areas of Xi'an, the majority of those who permanently reside in the countryside are elderly.Children and adults have higher mobility, as they work or study in urban areas during weekdays and only stay in the rural areas during weekends or holidays [41].For this population of residents, it is relatively difficult to develop thermal adaptability to the indoor thermal environment of rural dwellings.In the future promotion of the rural revitalization work, a large number of people will return to rural areas, and many personnel engaged in revitalization work will reside in rural areas.These people also lack thermal adaptability to rural dwellings [55].Therefore, after achieving further progress in the rural revitalization work and stabilizing the number of permanent residents in rural areas, the impact of thermal adaptability on the optimization design strategy for rural dwellings can be further considered.

Conclusions
This study took rural dwellings surrounding Xi'an City as an example, conducted field investigation on rural dwellings, simulated the influence of different design indicators on indoor PET during winter and summer, and obtained the composite effect of multiple design indicators on indoor PET.The study drew the following conclusions.
1. Through field investigation, it was found that the 11 rural dwellings in the study had similar characteristics in terms of homestead shape, courtyard spatial layout, building structure, building form, and functional room types, which was consistent with the current characteristics of rural dwellings around Xi'an.There were differences in the building orientation, number of floors, building materials, and building spatial dimensions.The value or type range of each design element was extracted through measurements.Based on the measurement results, the basic model of rural dwellings around Xi'an was established.2. This study simulated the indoor PET of the basic model on the winter solstice and summer solstice days under the variation of single design indicators, including the building orientation, floor height, roof type, exterior wall cladding material, insulation layer setting method, window-to-wall ratio, and window material.The optimal values for each indicator that resulted in the best overall indoor thermal environment in winter and summer were as follows: south-facing orientation, 3.0 m floor height, color steel roof, cement exterior wall cladding material, external insulation, 40% south window-to-wall ratio, 10% north window-to-wall ratio, and aluminum alloy double glass window material.3. Through composite-variable simulation, this study obtained the indoor PET simulation, resulting in 8400 working conditions corresponding to various building lengths, widths, heights, south window-to-wall ratios, north window-to-wall ratios, and insulation layer thicknesses.By comparing the results, a comprehensive design optimization strategy for rural dwellings was obtained.Compared to the basic model, this strategy optimized the building size, roof type, insulation layer setting method, window-towall ratio, and material.The optimized building had a length of 12.6 m, width of 7.8 m, floor height of 2.7 m, with a south window-to-wall ratio of 10%, and a north window-to-wall ratio of 10%.It used a color steel roof, a 60 mm-thick external insulation layer, a cement exterior wall cladding material, and an aluminum alloy double glass window material.By optimizing the design indicators, the indoor thermal envi-

Figure 1 .
Figure 1.Design characteristics of existing rural dwellings around Xi'an (self-made by the author based on [36,37]).(a) Spatial layout.(b) Building construction.

Figure 1 .
Figure 1.Design characteristics of existing rural dwellings around Xi'an (self-made by the author based on [36,37]).(a) Spatial layout.(b) Building construction.

Figure 2 .
Figure 2. Flowchart of the software simulation method.

Figure 2 .
Figure 2. Flowchart of the software simulation method.

Figure 3 .
Figure 3. General plan and household selection of Shachang Village (Numbers are the serial number of each household).

Figure 4 .
Figure 4. General plan and household selection of Guanjia Village (Numbers are the serial number of each household).

Figure 3 .Figure 3 .
Figure 3. General plan and household selection of Shachang Village (Numbers are the serial number of each household).

Figure 4 .
Figure 4. General plan and household selection of Guanjia Village (Numbers are the serial number of each household).

Figure 4 .
Figure 4. General plan and household selection of Guanjia Village (Numbers are the serial number of each household).
living room was 4 m, and the width of bays on both sides was 3 m.The width of the auxiliary room was 3 m, and the length was 6 m.The detailed dimension values are shown in the dimension labeling of Figure 5.The other indicator design values of the basic model were also determined by values and types that were approximate to the actual situation, and the detailed values are shown in Table

Figure 5 .Table 10 .Figure 5 .
Figure 5. Spatial layout of the basic model.Table 10.Design values of the basic model.Indicator Value Number of rooms 9 Number of main rooms Living room 1 Bedroom 2

Figure 6 .
Figure 6.Digital model of the basic model.

4. 3 .
Univariate Simulation Results and Evaluation 4.3.1.Univariate Simulation Process We simulated the change in the PET value under the change in a single indicator, and then analyzed the influence of each indicator on the indoor thermal environment.The range of indicators that need to be simulated and compared is shown in Table

Table 2 .
Comparison of different simulation engines.

Table 3 .
Main classification of rural dwellings' design types.

Table 3 .
Main classification of rural dwellings' design types.
Sustainability 2023, 15, x FOR PEER REVIEW 8 of 27 Exterior wall cladding material Ceramic tile, cement, brick
Sustainability 2023, 15, x FOR PEER REVIEW 12 of 30 Plan Number

Table 5 .
Main living room design elements and ranges.

Table 6 .
The main parameters' average values of outdoor climate conditions.

Table 7 .
Thermal parameters of structural materials.

Table 8 .
Thermal parameters of window materials.

Table 9 .
Personnel object condition settings.

Table 10 .
Design values of the basic model.
Figure 6.Digital model of the basic model.

Table 16 .
Compound-variable PET optimization results (top ten items after sorting).

Table 17 .
Optimization strategy design indicators.