Coupling Effect of Space-Arrangement and Wall Thermal Resistance on Indoor Thermal Environment of Passive Solar Single-Family Building in Tibet

In areas where solar energy is abundant, such as the Tibetan plateau, passive solar buildings are attracting more and more attention and becoming a popular form of rural building. However, it is often difficult to achieve the satisfactory indoor thermal environment in a local rural passive solar single-family house. In order to improve the indoor thermal environment of passive solar buildings through building design, a systematic study of rural single-family buildings in Tibet was conducted. The basic parameters were investigated on the outdoor thermal environment, space-arrangement, envelope structure, and the indoor thermal environment. The building model considering space-arrangement modes was developed based on the survey data in multi-space passive solar buildings. The general physical and mathematical analysis models of multi-space passive solar buildings were established based on the heat transfer theory. Furthermore, the effects of space-arrangement and exterior wall thermal resistance on indoor air temperature were analyzed by numerical simulation. Results show that the indoor air temperature of the passive solar building is influenced by space-arrangement and wall thermal resistance together. When the space-arrangement of the building model was changed from “north-south through type” (mode a) to “through and separation combination” (mode b) and “north-south separation” (mode c), the indoor air temperature of the living room increased from 8.8 °C to 10.6 °C and 11.6 °C, with increases of 20.5% and 31.8%, respectively. In addition, equally increasing the thermal resistance of exterior walls in different orientations has different effects on the indoor air temperature. In the space-arrangement mode c, comparing with the temperature increment of the living room and bedroom caused by increasing thermal resistance of the south wall and north wall, the temperature increment of the living room caused by increasing thermal resistance of the east/west wall increased by 151.7% and 32.7%, and that of the bedroom increased by 609.1% and 239.1% respectively. This study can provide a reference for the optimal design of passive solar buildings in solar energy abundant areas.


Introduction
There is a growing concern over energy supply and ecological environment conservation in the world. Utilizing a renewable, safe, and eco-friendly energy resource to cater the energy demand of buildings has undoubtedly been regarded as an encouraging solution to achieve sustainable development for buildings.
The Tibet plateau has a peculiar environment, which is characterized by high altitude of over 4000 m, low air temperatures, and transparent and rarified atmosphere with minimal water and dust content. This region has inexhaustible resources of solar energy for annual solar radiation intensity method of no-uniformity insulation envelope configuration was built on a single-space solar building. The heat transfer process of multi-space solar buildings is more complicated than that of single-space buildings. It was indicated that the external boundary thermal conditions varied with the orientation, and an open thermodynamic system which had energy exchanges with the outer surroundings and adjacent spaces in multi-space buildings [29,30]. In other words, in the multi-space solar building, there existed simultaneously diversified indoor and outdoor thermal boundary conditions, and heat transfer between the indoor and outdoor, and between adjacent indoor spaces. Obviously, thermal resistance change of exterior walls in different directions has diverse effects on the indoor thermal environment. The renovation that effectively combines the optimization of the wall thermal performance and the indoor space-arrangement is a useful idea to improve the indoor thermal environment in the passive solar building in winter. However, related research has been scarcely reported.
In this paper, we present results from the field survey and measurement designed to obtain the characteristics of residential buildings and the climatic characteristics in the Lhasa area. Based on the local climatic and building feature, the mathematical models of the multi-space solar buildings were established. The coupling influences of indoor space-arrangement and exterior wall thermal resistance on the indoor thermal environment of multi-space solar building were analyzed. In view of the fact that the indoor air temperature is a key parameter for characterizing the indoor thermal environment, in order to simplify the analysis and evaluation, in this study, the indoor air temperature was adopted as an indoor thermal environment evaluation parameter. The result of this study can provide reference for the optimization of solar building thermal design.

Field Survey and Testing of Solar Buildings
The winter in Lhasa is cold and the local residential buildings are in the heating period from November 15th to March 15th of the following year. To obtain the characteristics and current situation of the local solar houses in the Lhasa area, a field survey was conducted from December 26 to 30,2016. A typical local single-story rural single-family house ( Figure 1) was selected to be tested and surveyed in Lhasa area. The testing building is a traditional north-south oriented passive solar building. The main rooms are located in the south, and the auxiliary rooms in the north. The building envelope parameters were measured, as well as the space layout, indoor and outdoor thermal environment, etc. The survey results show that the 500 mm thick adobe walls delimit the indoor and outdoor space, the living room and the kitchen are separated by a thin plate, and the remaining inner partition walls are 500 mm thick adobe walls. Single glazed windows and wooden doors are used and the exterior walls are painted white. The window-to-wall ratios for the south, east, and north directions are 0.42, 0.20, and 0.10, respectively.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 19 single-space solar building. The heat transfer process of multi-space solar buildings is more complicated than that of single-space buildings. It was indicated that the external boundary thermal conditions varied with the orientation, and an open thermodynamic system which had energy exchanges with the outer surroundings and adjacent spaces in multi-space buildings [29,30]. In other words, in the multi-space solar building, there existed simultaneously diversified indoor and outdoor thermal boundary conditions, and heat transfer between the indoor and outdoor, and between adjacent indoor spaces. Obviously, thermal resistance change of exterior walls in different directions has diverse effects on the indoor thermal environment. The renovation that effectively combines the optimization of the wall thermal performance and the indoor space-arrangement is a useful idea to improve the indoor thermal environment in the passive solar building in winter. However, related research has been scarcely reported.
In this paper, we present results from the field survey and measurement designed to obtain the characteristics of residential buildings and the climatic characteristics in the Lhasa area. Based on the local climatic and building feature, the mathematical models of the multi-space solar buildings were established. The coupling influences of indoor space-arrangement and exterior wall thermal resistance on the indoor thermal environment of multi-space solar building were analyzed. In view of the fact that the indoor air temperature is a key parameter for characterizing the indoor thermal environment, in order to simplify the analysis and evaluation, in this study, the indoor air temperature was adopted as an indoor thermal environment evaluation parameter. The result of this study can provide reference for the optimization of solar building thermal design.

Field Survey and Testing of Solar Buildings
The winter in Lhasa is cold and the local residential buildings are in the heating period from November 15th to March 15th of the following year. To obtain the characteristics and current situation of the local solar houses in the Lhasa area, a field survey was conducted from December 26 to 30, 2016. A typical local single-story rural single-family house ( Figure 1) was selected to be tested and surveyed in Lhasa area. The testing building is a traditional north-south oriented passive solar building. The main rooms are located in the south, and the auxiliary rooms in the north. The building envelope parameters were measured, as well as the space layout, indoor and outdoor thermal environment, etc. The survey results show that the 500 mm thick adobe walls delimit the indoor and outdoor space, the living room and the kitchen are separated by a thin plate, and the remaining inner partition walls are 500 mm thick adobe walls. Single glazed windows and wooden doors are used and the exterior walls are painted white. The window-to-wall ratios for the south, east, and north directions are 0.42, 0.20, and 0.10, respectively. A field measurement was conducted to obtain the indoor and outdoor thermal environment parameters. For the indoor air temperature measurement, the position of the room center 1.5 m from A field measurement was conducted to obtain the indoor and outdoor thermal environment parameters. For the indoor air temperature measurement, the position of the room center 1.5 m from floor was selected as the representative measurement point. The location of the measurement point of the outdoor thermal environment parameter was set in an open space near the testing building. More specifically, the test point of the outdoor air temperature was placed at a position that was not exposed to sunlight and was at a height of 1.5 m from the ground. Measurement of solar radiation intensity was carried out in an open space without any shade. The measuring points of the indoor air temperature were arranged as shown in Figure 2.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 4 of 19 floor was selected as the representative measurement point. The location of the measurement point of the outdoor thermal environment parameter was set in an open space near the testing building. More specifically, the test point of the outdoor air temperature was placed at a position that was not exposed to sunlight and was at a height of 1.5 m from the ground. Measurement of solar radiation intensity was carried out in an open space without any shade. The measuring points of the indoor air temperature were arranged as shown in Figure 2. A 175-H1 type self-recording high precision temperature meter with an accuracy of 0.1 °C was selected for the measurement of the indoor and outdoor air temperature. All the indoor and outdoor thermal environment parameters were measured and recorded every 10 min. The geometric sizes of building and envelope were measured using a laser distance meter. The main information about the testing process is shown in Table 1. A 175-H1 type self-recording high precision temperature meter with an accuracy of 0.1 • C was selected for the measurement of the indoor and outdoor air temperature. All the indoor and outdoor thermal environment parameters were measured and recorded every 10 min. The geometric sizes of building and envelope were measured using a laser distance meter. The main information about the testing process is shown in Table 1. The outdoor thermal environment mainly includes outdoor air temperature and solar radiation intensity. The average diurnal variation of solar radiation intensity and the outdoor air temperature over the testing period is shown in Figures 3 and 4.

Outdoor Thermal Environment
The outdoor thermal environment mainly includes outdoor air temperature and solar radiation intensity. The average diurnal variation of solar radiation intensity and the outdoor air temperature over the testing period is shown in  As shown in Figure 3, the average sunshine duration is approximately 9 h per day in the test period, and the average radiation intensity is 436 W/m 2 . The highest solar radiation intensity is high up to 788 W/m 2 at 14:00 pm. The proportion of scattered radiation to total radiation is approximately 35%. Figure 4 presents that the average, the highest, and the lowest outdoor air temperature are −1.7 °C, 4.9 °C, and −9.8 °C, respectively. Outdoor air temperature fluctuated significantly during the testing day, with a range of 14.7 °C.
The measured results indicate that, in the Lhasa area, the solar radiation intensity is high, the sunshine duration is long, and the outdoor air temperature fluctuates significantly, but the outdoor air temperature is not too low. It is an ideal climatic condition for passive solar buildings to use solar energy to meet the heat need for winter heating [31].

Indoor Thermal Environment
The average diurnal variation of indoor air temperature in each room over the testing period is illustrated in Figure 5. It is obvious that the indoor air temperatures in the living room and kitchen with auxiliary heat sources are higher than those of other rooms without auxiliary heat sources. The indoor air temperature in the living room and the kitchen fluctuates from 9.0 °C to 13.5 °C, and 7.5 °C to 12.0 °C. The indoor air temperature in the hall and bedroom varies in the range of 6.3 °C to 10.7 °C, and 6.4 °C to 11.1 °C. The indoor air temperature in the storage room in the west and middle fluctuates from 1.9 °C to 4.5 °C, and 3.5 °C to 5.3 °C, respectively.

Outdoor Thermal Environment
The outdoor thermal environment mainly includes outdoor air temperature and solar radiation intensity. The average diurnal variation of solar radiation intensity and the outdoor air temperature over the testing period is shown in  As shown in Figure 3, the average sunshine duration is approximately 9 h per day in the test period, and the average radiation intensity is 436 W/m 2 . The highest solar radiation intensity is high up to 788 W/m 2 at 14:00 pm. The proportion of scattered radiation to total radiation is approximately 35%. Figure 4 presents that the average, the highest, and the lowest outdoor air temperature are −1.7 °C, 4.9 °C, and −9.8 °C, respectively. Outdoor air temperature fluctuated significantly during the testing day, with a range of 14.7 °C.
The measured results indicate that, in the Lhasa area, the solar radiation intensity is high, the sunshine duration is long, and the outdoor air temperature fluctuates significantly, but the outdoor air temperature is not too low. It is an ideal climatic condition for passive solar buildings to use solar energy to meet the heat need for winter heating [31].

Indoor Thermal Environment
The average diurnal variation of indoor air temperature in each room over the testing period is illustrated in Figure 5. It is obvious that the indoor air temperatures in the living room and kitchen with auxiliary heat sources are higher than those of other rooms without auxiliary heat sources. The indoor air temperature in the living room and the kitchen fluctuates from 9.0 °C to 13.5 °C, and 7.5 °C to 12.0 °C. The indoor air temperature in the hall and bedroom varies in the range of 6.3 °C to 10.7 °C, and 6.4 °C to 11.1 °C. The indoor air temperature in the storage room in the west and middle fluctuates from 1.9 °C to 4.5 °C, and 3.5 °C to 5.3 °C, respectively. As shown in Figure 3, the average sunshine duration is approximately 9 h per day in the test period, and the average radiation intensity is 436 W/m 2 . The highest solar radiation intensity is high up to 788 W/m 2 at 14:00 pm. The proportion of scattered radiation to total radiation is approximately 35%. Figure 4 presents that the average, the highest, and the lowest outdoor air temperature are −1.7 • C, 4.9 • C, and −9.8 • C, respectively. Outdoor air temperature fluctuated significantly during the testing day, with a range of 14.7 • C.
The measured results indicate that, in the Lhasa area, the solar radiation intensity is high, the sunshine duration is long, and the outdoor air temperature fluctuates significantly, but the outdoor air temperature is not too low. It is an ideal climatic condition for passive solar buildings to use solar energy to meet the heat need for winter heating [31].

Indoor Thermal Environment
The average diurnal variation of indoor air temperature in each room over the testing period is illustrated in Figure 5. It is obvious that the indoor air temperatures in the living room and kitchen with auxiliary heat sources are higher than those of other rooms without auxiliary heat sources. The indoor air temperature in the living room and the kitchen fluctuates from 9.0 • C to 13.5 • C, and 7.5 • C to 12.0 • C. The indoor air temperature in the hall and bedroom varies in the range of 6.3 • C to 10.7 • C, and 6.4 • C to 11.1 • C. The indoor air temperature in the storage room in the west and middle fluctuates from 1.9 • C to 4.5 • C, and 3.5 • C to 5.3 • C, respectively.
The frequency distributions of the indoor air temperature segments in all rooms are illustrated in Figures 6 and 7. For rooms without auxiliary heat supply, it can be seen from Figure 6 that the indoor air temperatures of the south-facing rooms (hall, bedroom) are higher than those in the north-facing rooms (storage rooms). Figure 7 demonstrates that, in the absence of auxiliary heat supply, the indoor air temperatures of the south-facing rooms (hall, bedroom) are above 7 • C with almost 80% frequency and those in the north-facing rooms (two storage rooms) are below 5 • C with the similar frequency.
The above phenomena are ascribed to the north-facing room as "temperature buffer" of the south-facing rooms [32]. The "temperature buffer" effectively reduces the heat transfer losses of the south-facing rooms. The above results show that the difference in ambient conditions of different orientations leads to a significant temperature difference between main and auxiliary rooms in the multi-space solar building. It is worth noting that the indoor air temperatures are low in the testing building. The survey results show that local residents have a good adaptability to such low indoor temperature, but also show an urgent desire to improve the indoor thermal environment. The frequency distributions of the indoor air temperature segments in all rooms are illustrated in Figures 6 and 7. For rooms without auxiliary heat supply, it can be seen from Figure 6 that the indoor air temperatures of the south-facing rooms (hall, bedroom) are higher than those in the north- The frequency distributions of the indoor air temperature segments in all rooms are illustrated in Figures 6 and 7. For rooms without auxiliary heat supply, it can be seen from Figure 6 that the indoor air temperatures of the south-facing rooms (hall, bedroom) are higher than those in the north- The frequency distributions of the indoor air temperature segments in all rooms are illustrated in Figures 6 and 7. For rooms without auxiliary heat supply, it can be seen from Figure 6 that the indoor air temperatures of the south-facing rooms (hall, bedroom) are higher than those in the northfacing rooms (storage rooms). Figure 7 demonstrates that, in the absence of auxiliary heat supply, the indoor air temperatures of the south-facing rooms (hall, bedroom) are above 7 °C with almost 80%

Physical and Mathematical Model
The annual heating period of residential buildings in Lhasa is from 15 November to 15 March of the following year, and January is the coldest month. Based on the data obtained from the survey, the building analysis models were proposed, and the mathematical analysis models of the passive solar building were established. Subsequently, the effects of space-arrangement and wall thermal resistance on the indoor air temperature in the coldest month were analyzed by software.

Building Analysis Model
According to the typical spatial layout and the evolution of rural buildings in the Lhasa area [33], three types of indoor space-arrangement modes were derived and shown in Figure 8 as the building analysis models. In Figure 8, mode a is "north-south through" type, mode b is "transfixion and separation combination" type and mode c is "north-south separation" type. Among them, mode a and mode b represent the transitional modes in the development of indoor space-arrangement of local rural buildings, while mode c represents the current space-arrangement of the modern rural building. The height of the room in all modes is 3.3 m and the parameters of the space-arrangement modes are shown in Table 2.  Figure 8, mode a is "north-south through" type, mode b is "transfixion and separation combination" type and mode c is "north-south separation" type. Among them, mode a and mode b represent the transitional modes in the development of indoor space-arrangement of local rural buildings, while mode c represents the current space-arrangement of the modern rural building. The height of the room in all modes is 3.3 m and the parameters of the space-arrangement modes are shown in Table 2.

Heat Balance on the Outer Surface of the Building Envelope
The external surface of the building envelope directly exchanges heat with the outdoor air and is indirectly affected by solar radiation, including direct solar radiation, sky-scattered radiation, the ground reflected radiation, atmospheric long-wave radiation, and ground long-wave radiation. When these factors simultaneously affect the outer surface of the building envelope, the thermal equilibrium equation [34] is: q S + q R + q B + q g = q 0 + q ra + q ca (1) in which, the q S , q R , q g , q 0 , and q ca are calculated as follows: Formula (6) is the heat transfer equation of the wall and the direction of heat flow is considered in one-dimensional perpendicular to the wall.

Heat Balance on the Inner Surface of the Building Envelope
The inner surface of the building envelope is subjected to conduction heat, direct radiant heat, and radiant heat between surfaces, and convective transfer heat with indoor air. These heat transfers can be incorporated into a heat balance equation. The general formula is stated as following: "Conduction heat" + "convective heat of indoor air" + "mutual radiant heat between surfaces" + "direct radiant heat" = 0.
At the moment n, for the unit area on i-th surface, the heat balance equation is:

Thermal Balance of Transparent Building Envelope
The exterior window and door of the building studied in this paper only are equipped with a double glaze. There exist four glazed surfaces in the outer window and door and consequent four heat balance equations. The heat balance equations on the inner surface of the glass are based on the following assumptions: (1) The heat capacity of the glass is ignored since the glass of the window is very thin; (2) the heat flow direction is perpendicular to the glass surface and one-dimensional; (3) the glass surface is opaque for long-wave radiation; (4) the inner and outer surfaces of the glass are isothermal; (5) the heat of short-wave radiation absorbed by the glass is distributed to the inner and outer surfaces in the same proportions. The heat balance equations of four surfaces are as follows:

Thermal Balance in Indoor Air
The heat loss disparity in every room results in the difference of indoor air temperatures between rooms in natural operating conditions because the rooms in multi-space solar buildings are arranged at diverse locations. The temperature difference between adjacent rooms forces heat transfer through the wall dividing adjoining spaces. In unit time, the rise of indoor air temperature in a room is related to radiation heat exchange, convection heat transfer, heat conduction, cold air infiltration, and other factors. An indoor air thermal equilibrium equation is expressed as following: The above formulas cover the heat balance equation of opaque building envelope inner and outer surfaces, transparent building envelope, and indoor air. These equations together constitute the mathematical model of multi-space passive solar buildings. Temperature indicators in every room will be obtained by solving these equations.

Influence of Space-Arrangement and Exterior Wall Thermal Resistance on Indoor Air Temperature
It can be seen from the previous analysis that, in the design of passive solar buildings, the main rooms such as the bedroom and the living room are usually arranged in the south orientation, and the auxiliary rooms in the north orientation. The indoor space layout combined with the directional characteristics of solar radiation is conducive to improving the thermal environment of the main room. However, it is worth noting that in multi-space passive solar buildings, there is heat transfer between indoor and outdoor, and between adjacent rooms. Therefore, the indoor air temperature of the multi-space solar buildings is affected by space-arrangement and thermal resistance of exterior walls.
In order to study the coupling influence of the space-arrangement mode and the thermal resistance of envelope on the indoor air temperature, an available computation program needs to be selected and used to implement numerical simulations. Design Builder is the first comprehensive interface simulation software developed for EnergyPlus. Owing to these powerful features, it was chosen as the software for this study. The weather data for the Lassa area was downloaded from the official website of the U.S. Department of Energy [35].
According to the parameters of space-arrangement modes in Figure 8 and Table 2, the building numerical models (Figure 9) for the three space-arrangement modes were established by the software of Design Builder and the indoor air temperature change in each room was analyzed. In addition, since the indoor air temperature is the most important parameter affecting the indoor thermal environment, in order to simplify the analysis and evaluation, the indoor air temperature is adopted as the indoor thermal environment evaluation parameter.  Table 3.   Table 3. In the numerical simulation analysis, the non-stationary calculation was performed, and the calculation period was determined to three months (1 December to 1 March of the following year). The influence of space-arrangement and wall thermal resistance on the indoor air thermal environment was investigated by analyzing the change of the average indoor air temperature values over the coldest month (January). In this study, the purpose of determining the three-month period of the numerical calculation is to eliminate the influence of the initial conditions in the non-stationary calculation process on the results.

Influence of the Space-Arrangement Mode on the Indoor Temperature
Space-arrangement modes affect the thermal process of passive solar buildings, and consequently affect the indoor air temperature. According to the parameters in Table 3, the Design Builder software was used to analyze the indoor air temperature in the coldest month (January) of the main room in a multi-space passive solar building. In three space-arrangement modes, the thermal resistances of exterior walls in different orientations were all set to 0.53 m 2 ·K/W, and the average calculated indoor air temperature values of the main rooms in January are illustrated in Figure 10. over the coldest month (January). In this study, the purpose of determining the three-month period of the numerical calculation is to eliminate the influence of the initial conditions in the non-stationary calculation process on the results.

Influence of the Space-Arrangement Mode on the Indoor Temperature
Space-arrangement modes affect the thermal process of passive solar buildings, and consequently affect the indoor air temperature. According to the parameters in Table 3, the Design Builder software was used to analyze the indoor air temperature in the coldest month (January) of the main room in a multi-space passive solar building. In three space-arrangement modes, the thermal resistances of exterior walls in different orientations were all set to 0.53 m 2 ·K/W, and the average calculated indoor air temperature values of the main rooms in January are illustrated in Figure 10.  Figure 10 indicates that the indoor air temperature of the bedroom and living room in mode c is the highest, followed by that of mode b, while that of mode a is the lowest. Apparently, the indoor space-arrangement modes have a significant impact on the thermal environment under the same thermal parameters of the building envelope. Compared to 8.8 °C of the indoor air temperature of the living room in the mode a, that of living room in the model b is 10.6°C, increased by 20.5%, and the indoor air temperature growth rate of the living room is higher than that of the bedroom. When  Figure 10 indicates that the indoor air temperature of the bedroom and living room in mode c is the highest, followed by that of mode b, while that of mode a is the lowest. Apparently, the indoor space-arrangement modes have a significant impact on the thermal environment under the same thermal parameters of the building envelope. Compared to 8.8 • C of the indoor air temperature of the living room in the mode a, that of living room in the model b is 10.6 • C, increased by 20.5%, and the indoor air temperature growth rate of the living room is higher than that of the bedroom. When the space-arrangement mode is changed from mode b to mode c, the increase rate of the air temperature of the bedroom is higher than that in the living room. The air temperature rise from 8.3 • C to 10.8 • C in bedroom 1, and from 8.5 • C to 11.1 • C in bedroom 2, and from 10.6 • C to 11.6 • C in the living room. It is worth noting that the indoor air temperature of bedroom 1, bedroom 2, and living room has increased by 30.1%, 30.6%, and 9.4%, respectively.
Comparing three space-arrangement modes and corresponding indoor air temperatures of the main rooms, it can be found that the north-facing auxiliary rooms in the model c not only enrich the function of the building interior space, but also act as a "temperature buffer" which is beneficial to improve the temperature of the main rooms.

Influence of Exterior Wall Thermal Resistance on Indoor Air Temperature
Due to the temperature difference between adjacent rooms in a multi-space passive solar building, there is not only heat transfer between indoor and outdoor, but also heat transfer between adjacent rooms. Therefore, the change of the thermal resistance of the exterior wall in different orientations has different effects on the indoor air temperature of each room. It is investigated that, in the space-arrangement mode c, the indoor air temperatures of the main rooms varied with the increasing thermal resistance of the exterior walls, and the result is shown in Figure 11.  It can be seen from Figure 11 that increasing the same amount of the thermal resistance of different exterior walls has different effects on the indoor air temperature of the main room. The indoor air temperatures of the mode c were analyzed under the condition of thermal resistance change of different exterior walls. It is deduced that the indoor air temperature of the mode c varies with the thermal resistances of diverse orientation walls. Figure 12 shows the increment of indoor air temperatures of all rooms caused by increasing thermal resistance of exterior walls in different directions. It can be seen from Figure 11 that increasing the same amount of the thermal resistance of different exterior walls has different effects on the indoor air temperature of the main room. The indoor air temperatures of the mode c were analyzed under the condition of thermal resistance change of different exterior walls. It is deduced that the indoor air temperature of the mode c varies with the thermal resistances of diverse orientation walls. Figure 12 shows the increment of indoor air temperatures of all rooms caused by increasing thermal resistance of exterior walls in different directions.
From Figure 12, it can be seen that the effect of the wall thermal resistance on the indoor air temperature is different even if the thermal resistance of walls in different orientations is equally increased. Through the comparative analysis of Figure 12a-c, it can be found that increasing the thermal resistance of east/west walls can significantly increase the indoor air temperature of the main rooms (living room and bedrooms). However, increasing the thermal resistance of the north wall can increase the indoor air temperature of the auxiliary room more obviously. different exterior walls has different effects on the indoor air temperature of the main room. The indoor air temperatures of the mode c were analyzed under the condition of thermal resistance change of different exterior walls. It is deduced that the indoor air temperature of the mode c varies with the thermal resistances of diverse orientation walls. Figure 12 shows the increment of indoor air temperatures of all rooms caused by increasing thermal resistance of exterior walls in different directions. From Figure 12, it can be seen that the effect of the wall thermal resistance on the indoor air temperature is different even if the thermal resistance of walls in different orientations is equally increased. Through the comparative analysis of Figure 12a-c, it can be found that increasing the thermal resistance of east/west walls can significantly increase the indoor air temperature of the main rooms (living room and bedrooms). However, increasing the thermal resistance of the north wall can increase the indoor air temperature of the auxiliary room more obviously. Figure 12a presents that as the thermal resistance of the south wall increases, the indoor air temperature in each room rises slightly, and the temperature rise in the main room is slightly higher than that in the auxiliary room. When the thermal resistance of the south wall increases from 0.53 m 2 ·K/W to 1.86 m 2 ·K/W, the indoor air temperature of the living room and bedroom (the average value of indoor air temperature in bedroom 1 and bedroom 2) increases by 0.29 °C and 0.22 °C, and the temperature increment of auxiliary rooms is less than 0.11 °C. However, under the condition of the same increment of the north wall thermal resistance, the air temperature of the auxiliary rooms increases significantly, while the indoor air temperature of the living room and bedroom only slightly increase by 0.55 °C and 0.46 °C, respectively (Figure 12b). Figure 12c indicates that the same increment of the east/west wall thermal resistance has a greater impact on the indoor air temperature of the main rooms (bedrooms and living room) than those of the south and north walls. When the thermal resistance of the east/west walls is increased from 0.53 m 2 ·K/W to 1.86 m 2 ·K/W, the indoor air temperature of the living room and the bedroom increases by 0.73 °C and 1.56 °C, respectively.  Figure 12a presents that as the thermal resistance of the south wall increases, the indoor air temperature in each room rises slightly, and the temperature rise in the main room is slightly higher than that in the auxiliary room. When the thermal resistance of the south wall increases from 0.53 m 2 ·K/W to 1.86 m 2 ·K/W, the indoor air temperature of the living room and bedroom (the average value of indoor air temperature in bedroom 1 and bedroom 2) increases by 0.29 • C and 0.22 • C, and the temperature increment of auxiliary rooms is less than 0.11 • C. However, under the condition of the same increment of the north wall thermal resistance, the air temperature of the auxiliary rooms increases significantly, while the indoor air temperature of the living room and bedroom only slightly increase by 0.55 • C and 0.46 • C, respectively (Figure 12b). Figure 12c indicates that the same increment of the east/west wall thermal resistance has a greater impact on the indoor air temperature of the main rooms (bedrooms and living room) than those of the south and north walls. When the thermal resistance of the east/west walls is increased from 0.53 m 2 ·K/W to 1.86 m 2 ·K/W, the indoor air temperature of the living room and the bedroom increases by 0.73 • C and 1.56 • C, respectively. Comparing with the temperature rise in the living room and bedroom caused by the thermal resistance increasing of the south wall and north wall, that of the living room caused by the thermal resistance increasing of the east/west wall increases by 151.7% and 32.7%, and that of the bedroom increases by 609.1% and 239.1%, respectively.
The above results show that in the natural working state of multi-space passive solar buildings, there is a significant difference between the heat transfer loss and the passive heat collection in each room, due to the directional characteristics of solar radiation and the complicated heat transfer process between adjacent rooms.

Influence of Wall Thermal Resistance on Indoor Air Temperature Deference between Main and Auxiliary Rooms
During natural operation, the indoor air temperature of the main room in a passive solar building is mainly affected by two aspects: One is the net amount of solar radiant energy collected by the building; the other is the proportion of heat allocated between the main and auxiliary rooms. The previous analysis results show that even if the thermal resistance of different exterior walls changes by the same amount, the temperature difference between the main and auxiliary rooms is different. Therefore, it is necessary to analyze the change trend of the temperature difference between the main and auxiliary rooms under the thermal resistance increment of exterior walls in different orientations. Figure 13 illustrates the variation of indoor air temperature difference between main and auxiliary rooms with the thermal resistance increasing of the exterior wall in the indoor space-arrangement mode c. To simplify the analysis, the indoor air temperature of the main room is defined as the average indoor air temperature of the bedroom 1, bedroom 2, and living room, and the indoor air temperature of the auxiliary room is defined as the average temperature of the storage room 5, storage room 6, and kitchen. Figure 14 presents the indoor air temperature differences between the main and auxiliary rooms when the thermal resistances of the exterior walls in different orientations are increased from 0.53 m 2 ·K/W to 1.86m 2 ·K/W.
From Figure 13 it can found that increasing thermal resistance of the north wall decreases the temperature differences between the main and auxiliary rooms. When the thermal resistance of the north wall is increased from 0.53 m 2 ·K/W to 1.86 m 2 ·K/W, the indoor air temperature differences between the main and auxiliary rooms decline from 6.3 • C to 4.9 • C, and decreases by 22.2%. Due to the temperature increasing of the auxiliary room, the temperature difference between the main and auxiliary rooms decreases with the thermal resistance increasing of the north wall. Compared with the air temperature rise caused by increasing thermal resistance of the south wall and the east/west wall from 0.53 m 2 ·K/W to 1.86 m 2 ·K/W, the thermal resistance increasing of the north wall can increase the indoor air temperature of the auxiliary room by 36.7% and 21.8%. Meanwhile, the temperature values in Figure 14 indicate that the proportion of thermal energy distributed in the main room can be improved by increasing the thermal resistance of the south exterior wall, but the overall indoor thermal environment of the building does not improve significantly. It can be deduced that as the heat transfer loss of the building through the east/west wall decreases with the increase of thermal resistance, the indoor air temperature of the main room increases significantly, while the temperature of the auxiliary room increases slightly. It can be seen from Figure 14 that when the east/west wall adopts the maximum thermal resistance of 1.86m 2 ·K/W, the indoor air temperature of the main room is 12.4 • C, which is 8.8% and 6.9% higher than that of the same thermal resistance of the south wall and the north wall. The above result is attributed to the increase of the east/west wall thermal resistance increasing the net heat gain of the main room, thereby increasing the indoor air temperature of the main room.
as the average indoor air temperature of the bedroom 1, bedroom 2, and living room, and the indoor air temperature of the auxiliary room is defined as the average temperature of the storage room 5, storage room 6, and kitchen. Figure 14 presents the indoor air temperature differences between the main and auxiliary rooms when the thermal resistances of the exterior walls in different orientations are increased from 0.53 m 2 ·K/W to 1.86m 2 ·K/W.  From Figure 13 it can found that increasing thermal resistance of the north wall decreases the temperature differences between the main and auxiliary rooms. When the thermal resistance of the north wall is increased from 0.53 m 2 ·K/W to 1.86 m 2 ·K/W, the indoor air temperature differences between the main and auxiliary rooms decline from 6.3 °C to 4.9 °C, and decreases by 22.2%. Due to the temperature increasing of the auxiliary room, the temperature difference between the main and auxiliary rooms decreases with the thermal resistance increasing of the north wall. Compared with the air temperature rise caused by increasing thermal resistance of the south wall and the east/west wall from 0.53 m 2 ·K/W to 1.86 m 2 ·K/W, the thermal resistance increasing of the north wall can increase the indoor air temperature of the auxiliary room by 36.7% and 21.8%. Meanwhile, the temperature values in Figure 14 indicate that the proportion of thermal energy distributed in the main room can auxiliary rooms with the thermal resistance increasing of the exterior wall in the indoor spacearrangement mode c. To simplify the analysis, the indoor air temperature of the main room is defined as the average indoor air temperature of the bedroom 1, bedroom 2, and living room, and the indoor air temperature of the auxiliary room is defined as the average temperature of the storage room 5, storage room 6, and kitchen. Figure 14 presents the indoor air temperature differences between the main and auxiliary rooms when the thermal resistances of the exterior walls in different orientations are increased from 0.53 m 2 ·K/W to 1.86m 2 ·K/W. 4   From Figure 13 it can found that increasing thermal resistance of the north wall decreases the temperature differences between the main and auxiliary rooms. When the thermal resistance of the north wall is increased from 0.53 m 2 ·K/W to 1.86 m 2 ·K/W, the indoor air temperature differences between the main and auxiliary rooms decline from 6.3 °C to 4.9 °C, and decreases by 22.2%. Due to the temperature increasing of the auxiliary room, the temperature difference between the main and auxiliary rooms decreases with the thermal resistance increasing of the north wall. Compared with the air temperature rise caused by increasing thermal resistance of the south wall and the east/west wall from 0.53 m 2 ·K/W to 1.86 m 2 ·K/W, the thermal resistance increasing of the north wall can increase the indoor air temperature of the auxiliary room by 36.7% and 21.8%. Meanwhile, the temperature values in Figure 14 indicate that the proportion of thermal energy distributed in the main room can be improved by increasing the thermal resistance of the south exterior wall, but the overall indoor thermal environment of the building does not improve significantly. It can be deduced that as the heat transfer loss of the building through the east/west wall decreases with the increase of thermal

Conclusions
(1) The indoor space-arrangement of passive solar buildings has a significant impact on the air temperature. Through the comparison of indoor air temperature in three modes, it is found that the space-arrangement of "north-south separation" (mode c) not only enriches the building space, but also enables the auxiliary room on the north side to function as a "temperature buffer", resulting in a higher temperature in the main room. The rationality of the evolution of the indoor space-arrangements of local passive solar buildings has been confirmed by investigation and analysis. (2) For the space-arrangement of the "north-south separation" (mode c), increasing the thermal resistance of east/west walls has a significant effect on the indoor air temperature of the main room. Compared with the temperature rise in the living room and bedroom induced by increasing thermal resistance of the south walls and north walls, that of the living room induced by increasing thermal resistance of the east/west walls increases by 151.7% and 32.7%, and that of the bedroom increases by 609.1% and 239.1%, respectively. (3) The existence of indoor air temperature difference between the main and auxiliary rooms in the multi-space passive solar building indicates that the optimal distribution of indoor thermal energy in different rooms can be realized through the reasonable design of the building space-arrangement and wall thermal resistance. When the indoor space is arranged as the "north-south separation" (mode c), the effect of equally increasing the thermal resistance of different exterior walls on the indoor air temperature difference between the main and auxiliary rooms is different. Among them, the thermal resistance increase of the east/west wall leads to the indoor air temperature difference increase, while increasing the thermal resistance of the north wall reduces the indoor air temperature difference between the main and auxiliary rooms.

Conflicts of Interest:
The authors declare no conflict of interest.

Nomenclature
The following abbreviations are used in this manuscript: Black-body radiation constant, Equal to 5.67 (W/m 2 ·K 4 ) (cρ) a , (cρ) r Unit heat capacity of outdoor and indoor air (kJ/m 2 ·K) E 0 , E i Long wave radiation intensity of the outer and inner surfaces (W/m 2 ) F k Inner surface area of building envelope (m 2 ) h Heat transfer coefficient of gas between glass layers (W /m 2 ·K) HE s (n) Heat removal capacity of air conditioning system at "n" moment (W) I Dθ , I dθ Direct radiation and sky-scattered radiation absorbed by envelope (W/m 2 ). For vertical plane, θ = 90 0 ; For horizontal plane, θ = 0 0 .

I Rθ
Ground reflected radiation intensity accepted by an inclined plane with a horizontal inclination of θ (W/m 2 ) I SH Total solar radiation intensity to be accepted by the ground (W/m 2 ) K Heat transfer coefficient of envelope (W/m 2 ·K) k 1 , k 2 Heat transfer coefficient of the outer glass and inner glass of a double-glazed window (W/m 2 ·K) L a (n) Amount of air permeability at "n" moment (m 3 /h) N i Total number of inner surfaces of rooms with different building envelopes q S Solar radiation absorbed by the outer surface of the envelope (W/m 2 ) q R Ground reflected radiant heat absorbed by outer surface of the building envelope (W/m 2 ) q B Intensity of long wave radiation absorbed by the outer surface (W/m 2 ) q g Ground radiation heat absorbed by envelope (W/m 2 ) q 0 Heat transfer of the outer surface of the envelope to the inner surface (W/m 2 ) q ca Convective heat transfer of the outer surface of the envelope to the surrounding air (W/m 2 ) q ra Amount of heat radiation from the outer surface of the envelope to the surrounding environment (W/m 2 ) q i (n) The amount of heat transfer obtained by i-th inner surface of the building envelope (W/m 2 ) q τ i (n) Solar radiation and internal disturbance radiation heat in the inner surface of the i-th building envelope (W/m 2 ) q c 1 (n) Heat from illumination, human body heat and device display heat at the n-th moment (W) q c 2 (n) Room heat loss caused by evaporation of water at "n" moment (W) S 1 , S 2 Radiant intensity absorbed by the outer surface and inner surface of the outer glass of a double-glazed window (W/m 2 ) S 3 , S 4 Radiant intensity absorbed by the outer surface and inner surface of the inner glass of a double-glazed window (W/m 2 ) T g Ground temperature (K) t a , t i , t 0 Temperature of outer surface and inner surface of the envelope, and outdoor air temperature (K) t r (n) Indoor air temperature (K) t i (n), t k (n) Inner surface temperature of the i-th and k-th building envelope (K) t k (n), t r (n) Temperature of the inner and outer surfaces of the building envelope at "n" moment (K) V Volume of the room (m 3 ) α Dθ , α dθ Absorptivity of direct solar radiation and the sky-scattering radiation Outer surface and inner surface temperature of the outer glass of a double-glazed window (K) θ 3 , θ 4 Outer surface and inner surface temperature of the inner glass of a double-glazed window (K) ε g Ground emissivity respectively ε ik Degree of system blackness between the i-th and k-th building envelope, as the product of the blackness of the i-th and k-th surface ε i Long-wave emissivity of the i-th surface ϕ g Radiation angle coefficient of envelope outer surface to ground. Value is 0.5 ϕ ik Radiation angle coefficient of the i-th inner surface of the building envelope to the k-th inner surface σ Stephen Boltzmann constant ρ g Average reflectivity of the ground ∆τ Unit time (1 h)