Numerical Study on the Suitability of Passive Solar Heating Technology Based on Differentiated Thermal Comfort Demand

2022-07-30 01:31XiaonaFanQinZhaoGuochenSangandYiyunZhu

Xiaona Fan,Qin Zhao,Guochen Sang and Yiyun Zhu

School of Civil Engineering and Architecture,Xi’an University of Technology,Xi’an,710048,China

ABSTRACT Indoor thermal comfort and passive solar heating technologies have been extensively studied.However,few studies have explored the suitability of passive solar heating technologies based on differentiated thermal comfort demands.This work took the rural dwellings in Northwest China as the research object.First,the current indoor and outdoor thermal environment in winter and the mechanism of residents’differentiated demand for indoor thermal comfort were obtained through tests,questionnaires,and statistical analysis.Second,a comprehensive passive optimized design of existing buildings was conducted,and the validity of the optimized combination scheme was explored using DesignBuilder software.Finally,the suitability of passive solar heating technology for each region in Northwest China was analyzed based on residents’differentiated demand for indoor thermal comfort.The regions were then classified according to the suitability of the technology for these.The results showed that the indoor heating energy consumption was high and the indoor thermal environment was not ideal,yet the solar energy resources were abundant.Indoor comfort temperature indexes that match the functional rooms and usage periods were proposed.For the buildings with the optimized combination scheme,the average indoor temperature was increased significantly and the temperature fluctuation was decreased dramatically.Most regions in Northwest China were suitable for the development of passive solar heating technology.Based on the obtained suitability of the technology for the regions of Northwest China,these were classified into most suitable,more suitable,less suitable,and unsuitable regions.

KEYWORDS Rural dwellings;differentiated thermal comfort demand;optimization and simulation;statistical analysis;passive solar heating technology;suitability

1 Introduction

Humans spend approximately 70–90%of their lifetime indoors,particularly the sick,weak,disabled,and elderly[1,2].Thus,with the progress of science and technology and the improvement in living standards,higher demands have been put forward for the thermal comfort of the indoor environment[3,4].However,the improvement in indoor thermal comfort is generally accompanied by a rapid increase in building energy consumption[5,6].It has been reported that buildings account for 40%of the total global energy consumption,and this proportion continues to increase[7,8].This study aims to effectively address this contradictory problem to possess a comfortable and energy-efficient indoor thermal environment.

Indoor temperature plays a vital role in the achievement of a comfortable indoor thermal environment.Additionally,building standards[9,10]emphasize that the energy demand of temperature control systems can be evaluated only after the indoor comfort temperature is determined.Evidently,it is crucial to determine the indoor comfort temperature index for developing comfortable and energy-efficient buildings.In thermal comfort research,however,the existing theories regarding indoor comfort temperature are poorly developed.Most of the existing studies are based on the assumption that residents have similar demands for indoor thermal comfort during day and night[11–13].However,several studies have demonstrated that the human body has different demands for indoor comfort temperature in the active and sleeping states[14–16].Additionally,previous studies generally regarded the indoor space of a building in its entirety[13,17,18].Although this assumption is conducive to simplifying the analysis process in building thermal engineering,it does not consider the differentiated thermal comfort demands of residents for different functional rooms[19,20].To summarize,the indoor comfort temperature varies with functional room and usage period.However,few discussions on this difference were undertaken in previous studies.This is highly unfavorable to the realization of energy-saving and comfortable buildings.Therefore,this study explored the differentiated indoor comfort temperature indexes by adopting a research method that differentiates functional rooms and usage periods.

The indoor temperature directly affects indoor thermal comfort and building energy consumption[21].Meanwhile,the space form,space dimensions,thermal parameters,construction practices,and passive thermal response of the building directly influence its thermal performance and solar thermal utilization[22–25].Thereby,these affect indoor temperature.However,previous studies have focused on exploring the role of individual influencing factors on the indoor thermal environment rather than the synergistic improvement effect of multi-factors[26–30].To summarize,it is crucial to thoroughly consider each influencing factor in the design of schemes for improving the indoor thermal environment.

As is widely known,the application of passive solar heating technology can enable buildings to achieve certain heating functions without the aid of mechanical power.Given this,many scholars have conducted studies on passive technologies.It has been demonstrated that the application of passive solar heating technology can significantly improve the indoor environment and reduce building energy consumption[31,32].Nevertheless,the application effect of this technology depends significantly on the regional climatic conditions[33,34].In addition,it is noteworthy that few documents are available to guide the promotion of passive solar heating technology in rural areas of China.In particular,there are rarely specifications or studies that can guide designers in the selection of suitable region and the corresponding optimal type[35–37].Undoubtedly,the suitability zoning of passive solar heating technology in rural areas would contribute significantly to the realization of comfortable and energy-efficient rural dwellings.

As a rapidly-developing country,China’s construction industry is undergoing rapid development and innovation.Nevertheless,the existing schemes for improving the thermal environment of rural dwellings are imperfect.The problems of low indoor thermal comfort and high heating energy consumption continue to be prominent[38–40].The main reason is the deficiency of scientific theoretical guidance for rural buildings[41,42].Most regions in Northwest China are located in cold or severe cold climate zones.Here,the climate conditions are harsh,and the ecological environment is vulnerable[43].As a result,the above-mentioned problems of the indoor environment are particularly prominent in the rural areas of Northwest China.However,most of Northwestern China has abundant solar energy resources.Considering the above analysis,this study takes the rural dwellings in Northwest China as the research object to explore the suitability of passive solar heating technology.

The investigation of the suitability of passive solar heating technology involves an analysis of the building thermal process under geographical zoning.An issue that needs to be focused upon while discussing building thermal processes is the heat transfer in the building envelope.To clarify the heat transfer condition of the envelope,it is necessary to solve the differential equation of heat conduction.The solving of the equation involves its discretization.At present,the commonly used discretization methods include the finite-volume method(FVM),finite-element method(FEM),and finite-difference method(FDM).However,a comprehensive comparison revealed that the FVM has a low computational accuracy[44].The FEM has limitations such as complex procedures and substantial calculations[45].It is noteworthy that the FDM exhibits strong applicability,computational efficiency,optional accuracy,and convenience to implement on a computer[46].In summary,the FDM is a highly effective numerical calculation method.

In addition,the two main specialized computer-aided engineering software programs for building thermal simulation are EnergyPlus and DesignBuilder(DB).EnergyPlus is a building energy simulation engine developed by the U.S.Department of Energy and Lawrence Berkeley National Laboratory.It is based on the building heat balance equation and the FDM to analyze building thermal performance,solar energy utilization scheme,and indoor temperature fields under unsteady heat transfer[47].It is noteworthy that with the development of science and technology,a variety of user-graphic interface software has emerged for EnergyPlus.Of these,and DB is one of the most powerful[47].Compared with EnergyPlus,DB has the characteristics of rapid modeling,convenient use,and powerful functions.Additionally,DB passed the ASHRAE test in 2006.This indicates that the simulation results of DB are consistent with those of EnergyPlus.Hence,DB is adopted in this study to simulate and analyze the suitability of passive solar heating technology.

To develop a comfortable and energy-efficient indoor thermal environment,this study undertook the following tasks with the rural dwellings in Northwest China as the research object:

·The current indoor and outdoor thermal environment in winter and the mechanism of residents’differentiated demand for indoor thermal comfort were obtained through tests,questionnaires,and statistical analysis.

·The indoor comfort temperature indexes that match the functional rooms and usage periods were explored.

·A comprehensive passive optimized design of existing buildings is conducted,and the validity of the optimized combination scheme was explored using DB.

·The suitability of passive solar heating technology for each region in Northwest China was analyzed using the obtained differentiated comfort temperature indexes.Furthermore,the regions were classified according to the suitability of the technology.

2 Methods

Ningxia is located in inland Northwest China(Fig.1a).The architectural forms and climate characteristics here are representative of Northwest China[7].Therefore,this study takes the rural dwellings in Ningxia,Northwest China,as the research object.We conducted a field survey on rural dwellings in Ningxia,Northwest China,during January 16–20,2015,and January 15–17,2017.The site photographs of the test and questionnaire are shown in Fig.2.

Figure 1:Location of the study regions.(a)Location of Ningxia in China map;(b)Location of 28 representative regions in Northwest China

Figure 2:Site photographs of the test and questionnaire.(a)Recording of test data;(b)Questionnaire for face-to-face;(c)Testing of solar radiation intensity;(d)Testing of indoor air temperature and relative humidity;(e)Testing of indoor black-globe temperature

2.1 Tests

2.1.1 Test Building

The field survey revealed two types of local rural dwellings:earth buildings(7%)and brickconcrete buildings(93%).A majority of the local dwellings were approximately square.The internal and external walls were made of solid clay bricks.The thicknesses of the internal and external walls were 370 mm and 240 mm,respectively.The external walls were decorated with 10 mm-thick white tiles or fair-faced concrete.The internal walls were painted with lime mortar.The roof was generally single sloped or double sloped.Only the bedroom had a suspended ceiling(gypsum board).The windows were generally made of single-glass aluminum alloy.The door and indoor door were made of aluminum alloy and wood,respectively.A representative brickconcrete dwelling(Fig.3a)was selected for the thermal environment test.The layout of indoor and outdoor measuring points is shown in Fig.3b.

Figure 3:Measured dwelling and layout of measuring points.(a)Image of measured dwelling;(b)Layout of indoor and outdoor measuring points

2.1.2 Test Contents and Instruments

The test site was located in Heilin Village,Ningxia,Northwest China.The data were collected from January 16 to 20,2015.The solar radiation intensity,relative humidity and temperature of air,air velocity,and black-globe temperature were measured.Considering the complexity involved in the measurement,the mean radiation temperature was assumed to be approximately equal to the black-globe temperature[48].The instruments used are presented in Table 1.

Table 1:Instruments for measuring parameters of the indoor and outdoor thermal environment

The measurement duration of solar radiation intensity was 10 h,and the acquisition time interval was 1 h.The measuring point was located outdoors without shelter.The monitoring period for the relative humidity and temperature of the air was 24 h,and the acquisition time interval was 10 min.The measuring points are approximately 1.1 m above the ground.

The mean radiant temperature and air velocity were not recorded continuously owing to the inconvenience of placing the black-globe thermometer and hot-bulb anemometer in the respondent’s home for a long time.The mean radiation temperature and air velocity were measured every 2 h.The globe thermometer and hot-bulb anemometer were placed approximately 1.1 m above the ground.

2.2 Questi onnaire

Verbal and written informed consent was obtained from the participants prior to the questionnaire survey.The participants were villagers who had lived in the region for a long time.To ensure that the participants completed the written questionnaire conveniently and effectively,each question in it was explained before they executed it.In the experiment,the participants were asked to maintain daily habits and emotional stability.They were paid after they had participated in the entire experiment.In addition,to induce the participants to support this experiment,we assured them that positive performers would be provided additional rewards.

The data were collected during January 16–20,2015 and January 15–17,2017.A total of 610 participants were received face-to-face questionnaires.Finally,583 valid questionnaires were recovered,which corresponded to an effective recovery rate of 95.6%.The questionnaires included basic life information and subjective perceptions of the participants.The thermal sensation votes(TSVs)were determined using the ASHRAE 7-point scale[49]:cold(–3),cool(–2),slightly cool(–1),neutral(0),slightly warm(+1),warm(+2),and hot(+3).The thermal sensation votes corresponding to acceptability were–1 to+1.The thermal preference votes used a three-level metric:cooler(–1),no change(0),and warmer(+1).

2.3 Optimization

2.3.1 Optimized Design

The thermal parameters and construction methods of the existing rural houses in Ningxia were optimized by combining the building standards[50,51]and field survey.The optimized thermal parameters and construction methods are presented in Table 2.Furthermore,the forms and dimensions of the indoor spaces in existing buildings were optimized based on the field survey and building standards[50,51].The results are displayed in Fig.4.

Table 2:Optimized thermal parameters and construction methods

Note:According to reference[20],the heat transfer resistance of the internal and external surface of the wall was 0.11(m2·K)/W and 0.04(m2·K)/W,respectively.The thermal conductivity correction coefficient of the insulation layer was 1.1.

To investigate the impact of different types of passive solar houses on the indoor thermal environment,three types of passive solar houses were introduced in the design.The optimized floor plans of these three types of houses are shown in Fig.4.The parameters of the representative components of passive solar houses are presented in Table 3.In addition,to facilitate comparative analysis,the rooms with the least favorable heat gain on the north and south sides were defined as room1 and room2,respectively.The direct-gain window passive solar house,Trombe wall+direct-gain window passive solar house,and attached sunspace passive solar house were denoted as mode1,mode2,and mode3,respectively.

2.3.2 Simulated Regions

Northwest China includes the three provinces Shaanxi,Gansu and Qinghai and the two autonomous regions of Ningxia and Xinjiang(Fig.1a).Here,the outdoor climatic conditions vary substantially across regions[52].28 representative regions in Northwest China were selected in this study by combining solar radiation classes,building thermal zone,and administrativegeographical divisions[51,53](see Fig.1b and Table 4).To further facilitate the analysis of the suitability law of passive solar heating technology in Northwest China,8 typical regions were selected from the 28 regions[54].The 8 regions and the corresponding meteorological data are shown in Table 5[10,55].

Table 3:Parameters of representative components of passive solar buildings

Table 4:Representative regions in Northwest China

Table 5:8 typical regions and the corresponding meteorological data

Figure 4:Optimized floor plans of three types of passive solar houses.(a)Direct-gain window passive solar house;(b)Trombe wall+direct-gain window passive solar house;(c)Attached sunspace passive solar house

Table 3:(continued)Passive solar house style Component Size/m Trombe wall+direct-gain window passive solar house Glass curtain wall 4.0×3.0 4.3×3.0 Air layer 0.08 The top vent 0.5×0.4×2 0.5×0.5×2 The bottom vent 1.0×0.4 1.2×0.4 The south window 1.5×2.2 1.8×2.2 Attached sunspace passive solar house Sunspace 4.5×1.5×3.5 4.8×1.5×3.5 The south window 4.0×2.8 4.3×2.8

2.3.3 Classi fication of Thermal Comfort Demands

The residents’demand for indoor thermal comfort was classified based on the responses to the questionnaire(see Table 6).

Table 6:Classification of residents’thermal comfort demands

2.4 Si mulation

2.4.1 Mathematical Model

Under natural operation,the indoor thermal environment of a building is determined by the building thermal process.This process involves a series of heat balance equations and thermal transfer equations.The heat balance equations include those of the internal(external)surface of the non-transparent envelope,internal(external)surface of the transparent envelope,and indoor air.

The heat balance equation of the external surface of the non-transparent envelope is expressed as[56]

whereqSis the solar radiation heat absorbed by the external surface of the envelope(W/m2),qRis the ground-reflected radiation heat absorbed by the external surface of the envelope(W/m2),qBis the atmospheric long-wave radiant heat absorbed by the external surface of the envelope(W/m2),qgis the ground radiation heat absorbed by the envelope(W/m2),q0is the heat transfer from the external surface to the internal surface of the envelope(W/m2),qcais the convective heat transfer between the external surface of the envelope and surrounding air(W/m2),andqrais the radiant heat between the external surface of the envelope and surrounding environment(W/m2).

The heat balance equation for the internal surface of the non-transparent envelope is expressed as[56]

whereqi(n)is the heat transfer obtained by the internal surface of the i-thenvelope at time n(W/m2),αciis the convective heat transfer coefficient of the i-thinternal surface of the envelope[W/(m2·K)],tr(n)is the indoor air temperature at time n(°C),Niis the total number of internal surfaces of different envelopes in the room,Cbis the blackbody radiation constant,with a value of 5.67[W/(m2·K)],φikis the radiation angle coefficient between the i-thand k-thinternal surfaces of the envelope,ti(n)andtk(n)are the temperature of the i-thand k-thinternal surface of the envelope,respectively(°C),εikis the system emissivity between the i-thand k-thinternal surfaces of the envelope,which is equal to the product of the emissivity of the i-thand k-thsurfaces itself.

The heat balance equation of the internal and external surfaces of the transparent envelope can be expressed as(for double-glazed windows as an example)[56]

whereEoandEiare the long-wave radiation intensity of the external and internal surfaces,respectively(W/m2),αoandαiare the air convection heat transfer coefficients of the external and internal surfaces,respectively[W/(m2·K)],toandtiare the outdoor and indoor air temperatures,respectively(°C),θiis the temperature of the air around the i-thsurface(°C),εiis the long-wave emissivity of the i-thsurface,Siis the intensity of short-wave or long-wave radiation absorbed by the i-thsurface(W/m2),kiis the heat transfer coefficient of the i-thlayer of glass[W/(m2·K)],his the heat transfer coefficient of the glass interlayer[W/(m2·K)],andσis the Stefan–Boltzmann constant[W/(m2·K4)].

The heat balance equation of indoor air is[56]

whereFkis the area of the k-thinternal surface of the envelope(m2),is the unit capacity of the kth internal surface of the envelope and the indoor air[kJ/(m3·K)],is the convective heat dissipation from lighting,human sensible heat and equipment sensible heat at time n(W),is the sensible heat of the room that is consumed by water evaporation at time n(W),HEs(n)is the sensible heat removal of the air conditioning system at time n(W),La(n)is the air infiltration capacity at time n(m3/h),Δτis the time increment,(cp)aand(cp)rare the unit heat capacity of indoor and outdoor air,respectively[kJ/(m3·K)],andVis the volume of the room(m3).

An issue that should be focused upon in the solution of the heat balance equation is the heat transfer by the envelope.

The indoor temperature varies continuously owing to the variation in outdoor weather conditions with time.Therefore,the heat transfer phenomenon of the envelope structure is considered to be an unsteady heat transfer process.The differential equation of unsteady heat conduction in the three-dimensional(3D)Cartesian coordinate system is expressed as[21]

wherex,yandzare the 3D Cartesian system coordinates,λiis the thermal conductivity along thei-axis direction[W/(m·K)],ρis the density(kg/m3),cis the specific heat capacity[W/(kg·K)],τis the time(s),qvis the heat flux(W/m2),andTis the temperature(°C).

However,considering that the height(length)and width of the external envelope are 8–10 times the thickness,the unsteady heat transfer of the external envelope is considered as a onedimensional unsteady heat transfer process.The one-dimensional differential equation of unsteady heat conduction is expressed as[56]

whereais the thermal diffusivity of the material,which is equal to(m2/h),λis the thermal conductivity of the material[W/(m2·K)],cis the specific heat capacity of the material[kJ/(kg·K)],andρis the density of the material(kg/m3).

Engineering calculation problems can be solved effectively and rapidly using the FDM.Using this method,the temperature gradient from the center of the k–1-th(k-th)floor to the center of the k-th(k+1-th)floor at time m is expressed as

whereΔxkis the thickness of the k-ththin layer of the envelope(k=0,1,2,3...n),is the central temperature of the k-thlayer at time m.

Based on Eqs.(9)and(10),we obtain

whereΔτis the time increments.

Eq.(11)is substituted into Eq.(9)to get

whereλkis the thermal conductivity of the k-thmaterial layer of the envelope[W/(m2·K)],ckis the specific heat capacity of the k-thmaterial layer of the envelope[kJ/(kg·K)],ρkis the density of the i-thmaterial layer of the envelope(kg/m3).

Because the envelope is a homogeneous material layer,Eq.(12)can be simplified as

whereak,k-1,ak,k,andak,k+1are dependent on the thermal properties and dimensions of the envelope material.

To determine the specific heat transfer conditions,the initial and boundary conditions need to be specified in addition to the heat conduction equation.

The boundary condition is[56]

whereLis the total thickness of the envelope.

The initial condition is[56]

The finite difference equation system(matrix expression form)of one-dimensional unsteady thermal conductivity is obtained according to Eqs.(15)–(17).

Further,the sequential Gaussian elimination method was used to solve the equation system to obtain the temperature of the surface of the envelope at each moment.

2.4.2 Development of DB Model

The process of the simulation using DB is as follows:(1)establishment of 3D model,(2)meshing,(3)setting of parameters,(4)setting of boundary conditions,and(5)verification of the model.It was worth noting that the meshing was done automatically by the DB in the running background.Therefore,this step was not shown during the simulation.

(1)3D models

The 3D models for the simulation were established according to Fig.4 in Section 2.3.1.The results are shown in Fig.5.

Figure 5:Developed 3D models.(a)Direct-gain window passive solar building;(b)Trombe wall+direct-gain window passive solar building;(c)Attached sunspace passive solar building

(2)Setting of parameters

The construction parameters of the building are shown in Table 2.The internal disturbance parameters of the building are listed in Table 7.During the simulation,the timestep was set to 12.The grid size was automatically selected by DB.

Table 7:Internal disturbance parameters of the building

(3)Setting of boundary conditions

China is located in the northern hemisphere of the earth.The day with the shortest sunshine time and the lowest solar altitude angle is the winter solstice.From the perspective of passive heat collection,the solar radiation conditions on this day are the worst for the year.Therefore,the winter solstice day was considered as the simulated weather date to explore the suitability of passive solar heating technology for each region in Northwest China.

(4)Validation of the model

The correctness of the model was verified to compare the tested data with the simulated data.The tested data were the indoor temperature of the tested dwelling on January 19,2016.The comparison results between the simulated data and the tested data are presented in Fig.6.

As shown in Fig.6,the measured temperature was highly consistent with the simulation temperature.The R-square values of the main and secondary functional rooms were 0.97 and 0.95,respectively.The deviation between the simulated and measured values are likely to have resulted from(1)the objective errors of the data test instruments and(2)the effect of infiltration of cold air during the measurement because of closing or opening of windows and doors.Given these reasons,the error is considered acceptable.Therefore,the numerical model established by DB is considered valid.

Figure 6:Comparison between measured and simulated results

3 Results and Discussion

3.1 Thermal Envi ronment

During the test,only the bedroom had a simple heating device(coal stove).It was operated from 6:00 am to 10:00 pm.The meteorological data of a sunny day(January 19,2015)during the test period were selected for analysis in this study.

The variations in indoor and outdoor meteorological parameters during the test are shown in Fig.7.

As shown in Fig.7,the outdoor temperature showed an increasing trend from 10:00 am to 6:00 pm,and a decreasing trend from 6:00 pm to 10:00 am the next day.The outdoor temperature varied from–7.6°C to 2.5°C,with a temperature fluctuation of 10.1°C and an average temperature of–2.7°C.Additionally,the total solar radiation intensity was 290.8 W/m2on average.The peak value was 544.8 W/m2,which appeared at approximately 1:00 pm.The average direct radiation intensity was 216.7 W/m2,which accounted for 74.5%of the total radiation intensity.The indoor relative humidity of the main bedroom(main functional room)and storage room(secondary functional room)varied within the ranges 45%–52%and 30%–43%,respectively.This implied that the indoor relative humidity satisfied the thermal comfort demands of the human body[50].Fig.7 also shows that the average indoor temperature of main and secondary functional rooms from 6:00 am to 10:00 pm was 10.5°C and 5.8°C,respectively.From 10:00 pm to 6:00 am,the average indoor temperature of the main functional room was 7.3°C.The indoor temperature obtained in this field test is similar to that obtained in the previous study for Northwest China[58–60].That is,the indoor temperature is low,and the temperature varies substantially across functional rooms and usage periods.

In addition,the intermittent test results demonstrated that there was a marginal difference between the indoor mean radiation temperature and indoor air temperature(±0.5°C).The indoor airflow was significantly weak(within 0.1 m/s,which is lower than the human body’s air-flowperception threshold of 0.2 m/s[7]).

Figure 7:Variations in solar radiation intensity,air temperature,and relative humidity

During the data measurement,we also surveyed the indoor heating in winter.The survey showed that most families heated only the bedrooms.Each dwelling consumed approximately 3t of coal through the coal stove during the heating period.When all the rooms were heated,the heating energy consumption would be 2–3 times the present energy consumption.

To sum up,the outdoor temperature was low,and the temperature difference in winter in this area was large.Notwithstanding the high energy consumption for heating,the indoor thermal environment was still not ideal.However,the area was rich in solar energy resources.Given this,passive solar buildings can be developed according to local conditions.

3.2 Residents’Di fferentiated Demands for Indoor Thermal Comfort

The determination of indoor comfort temperature index was influenced by the regional climate,residents’behavior habits,psychological activities,etc.[61].

3.2.1 Thermal Resistance

Using a research method that differentiates functional rooms and usage periods,a detailed investigation was carried out on the clothing conditions of rural residents in Northwest China in winter.Further,the thermal resistance of residents’clothing during the activity period was obtained according to the ASHRAE Standard 55-2013.The results are presented in Fig.8.

Figure 8:Distribution ratio of the thermal resistance of residents’clothing during the activity period

Fig.8 illustrates that the clothing thermal resistance of residents in the main functional rooms was in the range of 1.5 to 1.7 clo(58.4%),with an average value of 1.6 clo.The clothing thermal resistance in the secondary functional rooms was in the range of 1.7 to 1.9 clo(57.8%),with an average of 1.8 clo.The main reason for this difference is that the region belongs to the cold region,where indoor heating is required in winter.However,most rural dwellings heat only the main functional rooms owing to economic constraints.Consequently,residents need to wear a thick coat as additional clothing whenever they enter the secondary functional rooms.Furthermore,the thermal resistance of clothing obtained in this study is significantly higher than the 0.9 clo specified by ASHRAE Standard 55-2013.The main reason is that the indoor thermal environment in winter is not ideal,and the economic capability of local residents to improve the thermal environment is weak.Therefore,local residents need to wear more clothes to ensure their thermal comfort.

In addition,the survey found that during the nighttime sleep period,the local rural residents were used to wearing pajamas and covering thick quilts.The thermal resistance of pajamas and thick quilts were 0.38 clo and 2.86 clo,respectively[49,62].

In summary,the thermal resistance of winter clothing of rural residents in Northwest China was large.Besides,the thermal resistance of residents’clothing varies significantly across functional rooms.The residents’demand for thermal resistance varies substantially across usage periods.

3.2.2 Activity Rate

The activity rate of residents in a room is defined as the number of residents in the room as a percentage of the total number of individuals in all the rooms.The activity rate of residents in different functional rooms was investigated at intervals of 1 h.The activity rate of residents in different functional rooms during the activity period is presented in Fig.9.

Figure 9:Activity rate of residents in the rooms during the activity period

Fig.9 indicates that the activity rate for the main functional rooms during the activity period(from 6:00 am to 10:00 pm) showed a decreasing–increasing–decreasing–increasing trend with an average value of 78%.Correspondingly,the activity rate for the secondary functional rooms showed an increasing–decreasing–increasing–decreasing trend,with an average value of 22%.

Moreover,the survey found that during the sleep period(from 10:00 pm to 6:00 am the next day),only the main functional rooms constituted the activity space of the local rural residents.Accordingly,the residents’ demands for indoor thermal comfort of secondary functional rooms at night are not discussed in the subsequent analysis.Referring to the ASHRAE Standard 55-2013,the metabolic rate of residents during sleep at night was 45 W/m2.The survey also found that the local rural residents were in the period of agricultural leisure during winter.During the activity period,the activities of the residents in the main functional rooms were watching TV,chatting,and sleeping for brief periods,whereas those of the residents in the secondary activity rooms were cooking,sewing,knitting,and maintenance.The metabolic rate of residents in main and secondary functional rooms was 85 W/m2and 95 W/m2,respectively,according to the ASHRAE Standard 55-2013.

To summarize,the activity rate of the same functional room varied across usage periods.Furthermore,the activity rates of different functional rooms varied across the same usage periods.

3.2.3 Thermal Acceptance Rate

The temperature was grouped with 0.5°C as the group distance.For each temperature group,the thermal acceptance rate (the ratio of the number of residents with voting values of –1 to+1 as a percentage of the total number of voters) was calculated.Finally,a polynomial fit was performed on the thermal acceptance rate and indoor air temperature.The fitting results areshown in Figs.10 and 11.

Figure 10:Acceptance temperature of main functional rooms.(a)The activity period;(b)The sleep period

Figure 11:Acceptance temperature of the secondary functional rooms

The acceptance rate of 80%of residents is regarded as the acceptable range of thermal environment[49,50,63].Fig.10 shows that the lower limits of the acceptance temperature for the main functional rooms in winter during the activity period and sleep period were 14.9°C and 12.9°C,respectively.Fig.11 shows that the lower limit of the acceptable temperature for the secondary functional room during the activity period was 11.8°C.Evidently,the acceptance temperature varied across functional rooms and usage periods.

3.2.4 Preferred Temperature

First,the temperature was grouped within 0.5°C as the group difference.Second,the percentage of people who expect the e nvironment to be hotter or colder than the current environment was counted,and linear regression was performed with the air temperature,respectively.The results are shown in Figs.12 and 13.

Figure 12:Preferred temperature of the main functional rooms.(a)The activity period;(b)The sleep period

Figure 13:Preferred temperature of the secondary functional rooms

Fig.12 shows that the preferred temperature for the main functional rooms during the activity and sleep periods were 15.5°C and 13.3°C,respectively.The difference in preferred temperature between these two periods was 2.2°C.Fig.13 displays that the preferred temperature for the secondary functional rooms during the activity period was 12.3°C.Based on the above analysis,the preferred temperature of the same functional room varied across usage periods.Furthermore,the preferred temperature differs across functional rooms for the same usage period.

3.3 Determination of Thermal Comfort Indexes

The thermal comfort demands of the local rural residents varied across functional rooms during the same period.Besides,the thermal comfort demands of the same functional room varied across usage periods.Thus,the indoor thermal comfort evaluation indexes based on residents’differentiated thermal demand should be explored.

The widely used indoor thermal comfort evaluation indexes are Operative Temperature(Top),Standard Effective Temperature(SET*),Predicted Mean Vote(PMV),and Subjective Temperature(Tsub).Topwas proposed as an indoor thermal comfort evaluation index based on the apparent heat transfer between the surrounding environment and human body by convection and radiation[64].The greater the indoor air velocity,the closer is the operating temperature to the air temperature.However,the indoor air velocity of rural dwellings in winter is generally significantly low.This implies that there is a limitation to the use of Topas an indoor thermal comfort evaluation index.The SET*was put forward as an indoor thermal comfort evaluation index based on a dynamic two-node model of the human temperature regulation[65].It is perfect as an evaluation index.However,the difficulty of calculating skin temperature and wettability hinders the acquisition of SET*[35].Fanger[66]recommended the use of PMV-PPD to evaluate indoor thermal comfort.This index integrates a variety of factors that affect human thermal comfort,wherefore it is considered to be a fairly comprehensive evaluation index.However,many scholars consider that it is difficult for most rural dwellings in China to achieve the standard of PMV=–0.5 to 0.5[7].McIntyre[67]recommended the use of Tsubas an evaluation index of indoor thermal comfort.Tsubwas defined as the temperature of the standard environment when the subjective thermal sensation of a person in a certain actual environment was the same as in the standard environment.This implies that the Tsubcan be used to guide the evaluation and design of indoor thermal environments.The standard environment was a uniform enclosed space withTa=Tr,v=0.1 m/s,relative humidity of 50%,and air temperature equal to the mean radiation temperature[49].

Based on the above analysis,it is considered that Tsubwas effective and feasible as an indoor thermal environment evaluation index.It can be calculated using Eq.(20)[49].

whereTs ubis the subjective temperature(°C),Rclois the total clothing thermal resistance(clo),andMis the metabolic rate(%).

The Tsubunder different conditions was obtained based on Section 3.2 and Eq.(20).The results are listed in Table 8.

Table 8 demonstrates that the comfort temperature of the main functional rooms during the activity period(6:00 am to 10:00 pm)was 15.1°C.The comfort temperature during the sleep period(10:00 pm to 6:00 am the next day)was 12.9°C.The comfortable temperature of the secondary functional rooms during the activity period was 12.0°C.

Table 8:Subjective temperature under different conditions

Meanwhile,the comfort temperature was estimated using Griffiths’method,as shown in Eq.(21)[68,69].The results are shown in Table 9.

Table 9:Comfort temperature by the Griffiths method

whereTcis comfort temperature(°C),Tiis indoor air temperature(°C),0 is‘neutral’thermal sensation vote,andαis Griffiths constant,with a value of 0.5[32,68].

A comparison of Tables 8 and 9 shows that the thermal comfort temperature indexes calculated by subjective temperature are basically consistent with those estimated by Griffiths’method.Therefore,it is considered that the comfort temperature obtained by subjective temperature is reliable and valid.

This study aims to develop a comfortable and energy-efficient indoor thermal environment.Given this,this study proposes that during the activity period(from 6:00 am to 10:00 pm),the average indoor temperature of the main and secondary functional rooms should be at least 15°C and 12°C,respectively.During the sleep period(from 10:00 pm to 6:00 am the next day),the average indoor temperature of the main functional rooms should be at least 13°C.

A comparison of the thermal comfort temperature indexes proposed in this study with the results of previous studies in Northwest China,it is found that the temperature range of the thermal comfort indexes is basically consistent[70–72].However,many studies have shown that the indoor comfort temperature varies with the functional rooms and usage periods,but this difference has not been taken into account in previous studies.The studies would have been more interesting if these had involved the exploration of indoor comfort temperature indicators that matched the functional rooms and usage periods.

3.4 Numerical Simul ation

3.4.1 Comparison of Indoor Temperature before and after Adopti on of Optimized Combination Scheme

There was no heating source in the room during the numerical simulation.Comparison of indoor temperature before and after optimization.The results are presented in Table 10.

Table 10:Comparison of indoor temperature between optimized and nonoptimized buildings

Table 10 shows that the average indoor temperature of the main and secondary functional rooms in the optimized building was higher than the measured indoor temperature by at least 3.6°C and 4.0°C,respectively.Table 10 shows that compared with the measured indoor temperature,the average indoor temperature of the main and secondary functional rooms in the optimized building was increased,at least by 3.6°C and 4.0°C,respectively.The temperature fluctuation was reduced by at least 32.6%and 75.0%,respectively.Evidently,the average indoor temperature of the building with the passive optimized combination scheme was increased significantly,and temperature fluctuation was decreased dramatically without an increase in energy consumption.

3.4.2 Suitability of Passive Solar Heating Technology

For brevity,the simulation results of 8 typical regions are expounded.The simulation results of the other 20 representative regions are listed in Table 11.

Table 11:Simulation results of 20 representative regions

Note:AP:6:00 am to 10:00 pm,and PA:10:00 pm to 6:00 am.

Most suitable region:Dunhuang and Jiuquan are located in the I-Cold and I-Severe Cold regions,respectively,which are extremely rich in solar energy resources.Longnan is located in the II-HSCW region,where the average outdoor temperature on the winter solstice day is 6.4°C.Although the solar radiation intensity in Longnan is marginally lower than that in Dunhuang and Jiuquan,its average outdoor temperature is significantly higher than that of Dunhuang and Jiuquan.

As shown in Fig.14,the indoor temperature of the three types of passive solar houses could satisfy the thermal comfort standard proposed by this research.During the activity period(from 6:00 am to 10:00 pm),the average indoor temperature of the main and secondary functional rooms was at least 15°C and 12°C,respectively.During the sleeping period(10:00 pm to 6:00 am the next day),the average indoor temperature of the main functional rooms was at least 13°C.Owing to the buffering effect provided by the attached sunspace,the temperature fluctuation in the room with attached sunspace is significantly smaller than that in the room with Trombe wall+direct-gain window or direct-gain window.Prior studies have also shown that the indoor thermal comfort of passive solar houses is affected not only by the average temperature,but also by temperature fluctuation[19,20,73].Therefore,it is strongly recommended to adopt passive solar heating technology in these regions.The passive solar house with attached sunspace is preferred by designers.The other two types of passive houses are available as alternatives.

More suitable region:Yinchuan is located in the II-Cold region,which is rich in solar energy resources.The average outdoor temperature on the winter solstice day is–4.0°C.

Figure 14:Indoor temperature curves of different passive solar houses.(a)Dunhuang;(b)Jiuquan;(c)Longnan

Fig.15 demonstrates that only the passive solar house with attached sunspace could satisfy the requirement that the average indoor temperature of the main and secondary functional rooms from 6:00 am to 10:00 pm was at least 15°C and 12°C,respectively.The average indoor temperature of the main functional rooms from 10:00 pm to 6:00 am the next day was at least 13°C.If the passive solar house with direct-gain window or Trombe wall+direct-gain window is adopted,auxiliary heating should be provided to the room,especially for the passive solar house with the direct-gain window.Thus,it is recommended to adopt passive solar heating technology in this region,and the design scheme of the passive solar house design with attached sunspace should be selected.

Less suitable region:Xining is located in the II-Severe Cold region,where the average outdoor temperature on the winter solstice day is–6.9°C.The region has abundant solar energy resources.Xi’an is located in the III-Cold region.The average outdoor temperature on the winter solstice day is 1.6°C,which is marginally higher than that in Xining.However,the level of solar energy resources here is lower than that in Xining.Hanzhong is located in the III-HSCW region.The average outdoor temperature on the winter solstice day is higher,with a value of 5.3°C.

Figure 15:Indoor temperature curves of different passive solar houses in Yinchuan

Fig.16 illustrates that the indoor temperature of the passive solar house with attached sunspace was higher than that of the passive solar house with Trombe wall+direct-gain window or direct-gain window.However,no matter the type of passive solar house was considered,an individual passive heating technology could not satisfy the indoor thermal comfort demands of rural residents in Northwest China.That is,it could not satisfy the demand that the average indoor temperature of the main and secondary functional rooms from 6:00 am to 10:00 pm was at least 15°C and 12°C,respectively.The average indoor temperature of the main functional rooms from 10:00 pm to 6:00 am the next day was at least 13°C.However,the solar energy resources in these regions have a relatively high potential for application.Thus,it is recommended that these regions develop passive solar heating technology to the extent feasible.To satisfy indoor thermal comfort standards,active solar heating technology may be adopted on the basis of the development of passive solar heating technology.

Figure 16:(Continued)

Figure 16:Indoor temperature of different passive solar houses.(a)Xining;(b)Xian;(c)Hanzhong

Unsuitable region:Karamay is located in the III-Severe Cold region.The average outdoor temperature on the winter solstice day is exceptionally low with–12.1°C.

Fig.17 indicates that regardless of the type of passive solar room design adopted in this study,the indoor temperature cannot satisfy the requirement of at least 12°C.In fact,all the types of passive solar house in this study are far from satisfying this requirement,including the passive solar house with attached sunspace.The main reason is the low outdoor temperature and solar radiation intensity in this region.To summarize,active solar heating technology can be considered for this region to achieve indoor thermal comfort.

Figure 17:Indoor temperature of different passive solar houses in Karamay

Based on the above analysis,the suitability of passive solar heating technology in different regions of Northwest China was obtained.It can be learned that most regions in Northwest China are suitable for the development of passive solar heating technology.Moreover,over 70%of the regions in Northwest China comply with the following law:I-C is the most suitable region,II-Cold is the more suitable region,II-Severe Cold and III-Cold are the less suitable regions,and III-Severe Cold is the unsuitable region.However,this regularity is not apparent in certain regions of Xinjiang and Qinghai(italics in Table 11).This may be owing to the differences in topography[54].However,this difference is not discussed in this study and will be considered in the follow-up study.

4 Conclusion

This work took the rural dwellings in Northwest China as the research object.The current indoor and outdoor thermal environment in winter and the mechanism of residents’differentiated demand for indoor thermal comfort were analyzed.On this basis,suitable differentiated indoor comfort temperature indexes and the passive optimized combination scheme for improving the indoor thermal environment were explored.Finally,the suitability of passive solar heating technology for each region of Northwest China was analyzed based on residents’differentiated demand for indoor thermal comfort.The conclusions are summarized below:

(1)In winter,the outdoor temperature was low and the temperature difference between day and night was large in the rural areas of Northwest China.Although the heating energy consumption was high,the indoor thermal environment was still not ideal.However,the area was rich in solar energy resources.

(2)It was proposed that during the activity period(from 6:00 am to 10:00 pm),the average indoor temperature of the main and secondary functional rooms should be at least 15°C and 12°C,respectively.During the sleep period(from 10:00 pm to 6:00 am the next day),the average indoor temperature of the main functional rooms should be at least 13°C.

(3)The indoor thermal environment of the buildings with optimized combination schemes was improved significantly.The average indoor temperature of the main and secondary functional rooms increased by at least 3.6°C and 4.0°C,respectively,and the temperature fluctuation decreased by at least 32.6%and 75.0%,respectively.

(4)The suitability of passive solar heating technology for each region of Northwest China was obtained.According to the suitability,the regions were classified into most suitable,more suitable,less suitable,and unsuitable regions.

Funding Statement:This research was supported by the National Natural Science Foundation of China(Grant Nos.52078419 and 51678483).This research was supported by the Doctoral Dissertation Innovation Fund of Xi’an University of Technology(310–252072116).

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.