Investigation of Indoor Asymmetric Thermal Radiation in Tibet Plateau: Case Study of a Typical Office Building

<jats:p>The unique climate in cold plateaus leads to long and cold winters, which result in the characteristics and creation of an indoor thermal environment different from that of plain areas. However, there is a lack of detailed research on and evaluation of indoor temperature distributions. This study took an office building in Lhasa as an example to study the indoor non-uniform temperature distributions with radiator and solar radiation. The indoor and outdoor thermal environment parameters were tested. Next, a numerical model was established and verified. On a typical winter weather day, although the average indoor air temperature and radiation temperature in Lhasa are higher than those in Beijing (both are cold areas), the temperature in both is lower than comfortable levels. The indoor vertical air temperature differences are below 3 °C, with a 5% dissatisfaction most of the time. Solar radiation and radiators seriously affect the uniformity of the indoor thermal environment. The radiation asymmetry in Lhasa is significant, and the maximum radiation temperature asymmetry can reach 8.73 °C. In addition, the setting of north-facing windows should be avoided as far as possible in Tibetan areas. Suitable design and evaluation standards should pay attention to the obvious phenomenon of uneven indoor temperature distribution.</jats:p>

section_content: Tibet is called the "Roof of the World" and the "Third Pole of the Earth" [1] due to its average altitude of over 4000 m and unique ecological features [2]. There is a population of approximately 3. 44 million currently residing in Tibet [3]. The climate in the Tibetan Plateau has several peculiar characteristics, such as low atmospheric pressure, low atmospheric humidity, small seasonal temperature fluctuation and wide diurnal temperature range, and strong solar radiation [1]. Compared with plain areas, the annual average temperature in Lhasa and Xigaze is 10 • C-15 • C lower than that in Chongqing, Wuhan, and Shanghai at similar latitudes, but the solar radiation in Lhasa is a third higher than in those cities, and sometimes even twice as high [1]. The indoor thermal environment of the Tibetan Plateau is more susceptible to the outdoor climate, making it more difficult to create a comfortable indoor thermal environment, especially during its long and cold winters. 

The mean indoor air temperature of the Tibet Plateau in winter is very low, only 5-12 • C for rural buildings [4] and 3. 8-7. 8 • C for traditional houses [5] in Lhasa; 2. 8 • C in Kangding [6] ; 4. 16 • C in Nyingchi [7] ; 2. 5 • C in Gannan [8] ; and 6 • C in Qamdo [9]. The indoor air temperature can reach 15 • C with space heating [10] or 11. 5 • C with passive solar design [11]. However, the range of operating temperature during winter in Lhasa varies from 12. 74 • C to 23. 15 • C [10, 12, 13], and the acceptable temperature range is 13. 16 • C-26. 76 • C [9]. The thermal discomfort rates are 66% [11] and 68% [13]. Based on the existing research results, residents feel dry [11, 12], and the resulting discomfort rate ranges from 10% to 57. 4% [4, 7, 9, 12, 14]. However, some researchers found that the indoor thermal environment satisfaction rate is 90% and the low humidity does not have negative impacts [15]. Moreover, the perception of discomfort can also be attributed to the special characteristics of solar radiation, sky background temperature, and outdoor temperature in Tibet, which makes the internal surface temperature vary greatly between walls and leads to variations in radiative heat exchange differences between the human body and each internal surface. Hence, asymmetric radiation is an obstacle to achieving satisfying indoor thermal comfort in Tibet. 

Fanger proposed the concept of "radiant temperature asymmetry (RTA)" based on a series of thermal sensation experiments; the concept is defined as the difference between the radiation temperature of the plane of the two opposite sides for a small plane element or of the environment on opposite sides of a person, and is evaluated by the percentage of dissatisfaction [16, 17]. In winter, RTA in a space may be caused by cold windows, walls, floors, ceilings, and heating panels on the wall, ceiling, and floor. Some studies involved the investigation of the effects of vertical radiant asymmetry. However, their conclusions differ. Chrenko [18] found that the comfort limit for vertical RTA caused by heated ceilings was 2 • C. The work conducted by Fanger et al. [17] has shown that the vertical RTA of 14 • C for cool ceilings and 4 • C for the warm ceilings could be tolerated without adversely affecting comfort, but the results of Loveday and Hodder et al. [19, 20] and Griffiths and McIntyre [21] have shown that there is no effect of vertical radiant asymmetry from cool (warm) ceilings on thermal comfort. Furthermore, the comfort limits for horizontal radiant asymmetries from walls (windows) are different. In the study by McNall and Biddison [22], the occupants felt uncomfortable when the temperature of the warm wall reached 35 • C higher than ambient temperature, and the temperature difference between cool walls and other walls was more than 11. 1 • C. In addition, Fanger [17] recommended that the radiant temperature asymmetries between warm walls and the cool walls should be 23 • C and 10 • C, respectively. Hence, although RTA is an accepted and promoted metric in the thermal comfort evaluation, there are remaining research questions concerning consensus on comfort bounds and systematic understanding in various indoor environmental settings, such as RTA caused by windows and radiators, among others. 

The International Standards Organization (ISO) 7730-2005 [23] and the American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) standard 55-2017 [24] prescribe limits for radiant asymmetry, but the limits are different based on different conditions. The ISO 7730-2005 standard prescribes the limit value of RTA caused by warm ceilings, cool walls, cool ceilings, and warm walls, and the percentage of those dissatisfied as a function of the RTA are 5% (category A and B) and 10% (category C) [23]. However, to better control local discomfort, when a seated occupant wears standard clothing and is exposed to comfortable operative temperature and not in a draught, ASHRAE recommends that the RTA should be less than 5-23 • C by testing the radiation temperatures on opposite planes of a small plane element with a height of 0. 6 m (the height at the 'center' for a seated person) [24]. The comfort limits of the RTA were only applicable to specific geometries and conditions in the experiment. In practice, the angle factors in the cold or warm part of the panel are diverse, and the parameters of the indoor thermal environment are changeable and asymmetrical. Moreover, the RTA of the three-dimensional human body is different from that measured with a planar element, especially in a confined space. As a consequence, defining comfort limits for RTA is not easy work. 

However, there is a lack of relevant studies on the phenomenon of non-uniform indoor temperature distribution and its evaluation in the Tibetan Plateau. The temperature differences of windows, wall, floors, and ceilings can be severe in this area. Especially when radiators are used for space heating, the temperature difference between panels can be further exacerbated. Thus, the influence of asymmetrical temperatures on the thermal environment cannot be ignored. Furthermore, when the radiator is working, the occupiers can potentially suffer from radiant asymmetry due to the cold panels (windows, floor, and walls) and warm panels (windows and radiators) at the same time. In addition, the temperature of windows may vary greatly between the morning and the afternoon, during which cold panels may become warm due to the strong solar radiation. However, these unique situations have not been considered in other studies. Therefore, a field test and numerical simulations were carried out in a typical office building in Lhasa to analyze the non-uniform distribution of the thermal environment. Moreover, the indoor thermal characteristics on a typical-weather day in the plateau (Lhasa) and plain (Beijing) cities in cold regions [25] were compared and analyzed. The degrees of VATD and RTA were quantitatively assessed as well as the corresponding percentage of dissatisfaction. This research will help to better understand, design, build, and evaluate indoor thermal comfort in the Tibetan Plateau.
section_name: Introduction

section_content: In this study, the indoor and outdoor environmental parameters of an office building were collected through on-site measurement, and the CFD simulation results were verified by the measured parameters. Both experimental data and simulation results were used for indoor temperature distribution analysis. In addition, to further study the uniqueness of the indoor thermal environment in the Tibet Plateau in the typical weather day of winter, the same computational model was used and the simulation results of Lhasa and Beijing were analyzed comparatively. At the same time, taking into account human exposure positions when the vertical air temperature difference (VATD) and RTA were evaluated. Figure 1 shows the framework of this study. and walls) and warm panels (windows and radiators) at the same time. In addition, the temperature of windows may vary greatly between the morning and the afternoon, during which cold panels may become warm due to the strong solar radiation. However, these unique situations have not been considered in other studies. Therefore, a field test and numerical simulations were carried out in a typical office building in Lhasa to analyze the non-uniform distribution of the thermal environment. Moreover, the indoor thermal characteristics on a typical-weather day in the plateau (Lhasa) and plain (Beijing) cities in cold regions [25] were compared and analyzed. The degrees of VATD and RTA were quantitatively assessed as well as the corresponding percentage of dissatisfaction. This research will help to better understand, design, build, and evaluate indoor thermal comfort in the Tibetan Plateau.
section_name: Methodology

section_content: In this study, the indoor and outdoor environmental parameters of an office building were collected through on-site measurement, and the CFD simulation results were verified by the measured parameters. Both experimental data and simulation results were used for indoor temperature distribution analysis. In addition, to further study the uniqueness of the indoor thermal environment in the Tibet Plateau in the typical weather day of winter, the same computational model was used and the simulation results of Lhasa and Beijing were analyzed comparatively. At the same time, taking into account human exposure positions when the vertical air temperature difference (VATD) and RTA were evaluated. Figure 1 shows the framework of this study. 

Step one: filed survey
section_name: Methodology

section_content: Step two: on-site measurement
section_name: Measure the indoor and outdoor thermal environment parameters in heating period

section_content: Validate the computational modeling by measured data

Step three: Computational modeling
section_name: Modeling the office room by actual situation and suitable assumptions

section_content: Case selection: the most typical office building in Lhasa Obtain on-site testing permission
section_name: Analysis of simulation results

section_content: Step four: nonuniform temperature analysis

Human exposure positions are considered to analyze the nonuniform temperature distribution using the model in step three

The average indoor temperature, VATD, and RTA are analyzed comparatively between Lhasa and Being at typical weather day Analyzing the applicability of existing regulations to cold plateau areas
section_name: Validation

section_content: Lhasa is located in the southeastern part of the Tibetan Plateau, with an altitude of 3648. 9 m. The geographical coordinates are 91 °06 ′ east longitude and 29 °36 ′ north latitude. According to the "Building Climatic Zone Classification Standard" (GB 50178-93) [25] and the "Design Code for Heating Ventilation and Air Conditioning of Civil Buildings" (GB 50736-2012) [26], Lhasa is a cold area and a heating area. A typical office building was selected in Lhasa (see Figure 2 ) for this research. Due to the abundant solar energy, the window-to-wall ratio of the south wall reached 0. 53; detailed information is
section_name: Field test

section_content: Lhasa is located in the southeastern part of the Tibetan Plateau, with an altitude of 3648. 9 m. The geographical coordinates are 91 • 06 east longitude and 29 • 36 north latitude. According to the "Building Climatic Zone Classification Standard" (GB 50178-93) [25] and the "Design Code for Heating Ventilation and Air Conditioning of Civil Buildings" (GB 50736-2012) [26], Lhasa is a cold area and a heating area. A typical office building was selected in Lhasa (see Figure 2 ) for this research. Due to the abundant solar energy, the window-to-wall ratio of the south wall reached 0. 53; detailed information is shown in Table 1. The height difference between the indoor and outdoor areas is 750 mm, the total building area is 1821. 32 m 2, and the story height is 3 m. Furthermore, this is a mixedstructure building. The wall adopts a concrete solid block 240 mm thick; the thermal properties of the enclosure structures is shown in Table 2. The windows use doubletempered insulating glass with a 9 mm air layer. As shown in Figure 3, all the offices are heated by radiators, which are installed under the south windows. shown in Table 1. The height difference between the indoor and outdoor areas is 750 mm, the total building area is 1821. 32 m 2, and the story height is 3 m. Furthermore, this is a mixed-structure building. The wall adopts a concrete solid block 240 mm thick; the thermal properties of the enclosure structures is shown in Table 2. The windows use doubletempered insulating glass with a 9 mm air layer. As shown in Figure 3, all the offices are heated by radiators, which are installed under the south windows. 

(a) (b) (c) shown in Table 1. The height difference between the indoor and outdoor areas is 750 mm, the total building area is 1821. 32 m 2, and the story height is 3 m. Furthermore, this is a mixed-structure building. The wall adopts a concrete solid block 240 mm thick; the thermal properties of the enclosure structures is shown in Table 2. The windows use doubletempered insulating glass with a 9 mm air layer. As shown in Figure 3, all the offices are heated by radiators, which are installed under the south windows. The field experiments were carried out from 14 to 16 March in 2018. The measured environmental parameters include indoor and outdoor air temperature, humidity, and wall temperature. The following devices were used to measure the variables: a HOBO U12-012 digital temperature and humidity meter and a Testo 830-T4 handheld infrared thermometer. Those can be seen in Figure 4. The specifics of each device are given below:

1. 

HOBO U12-012 digital temperature and humidity meter:

External window 3. 2 5(glass) + 9(air) + 5(glass)

The field experiments were carried out from March 14th to 16th in 2018. The measured environmental parameters include indoor and outdoor air temperature, humidity, and wall temperature. The following devices were used to measure the variables: a HOBO U12-012 digital temperature and humidity meter and a Testo 830-T4 handheld infrared thermometer. Those can be seen in Figure 4. The specifics of each device are given below:

1. HOBO U12-012 digital temperature and humidity meter:

Temperature range: -20-70 °C ; accuracy: ± 0. 21 °C ; humidity range: 1-95%; accuracy: ± 2. 5%. The field measurements were carried out based on the method recommended in ASHRAE 55-2017 [27]. More specifically, for the outdoor environmental parameters, the location of the measuring points was open ground near the office building at a height of 1. 1 m without direct sunlight. For the indoor environmental parameters, a measuring point at a height of 1. 1 m in an occupied office room was selected as the representative point for the measurement of air temperature and relative humidity because the occupants were usually in a sitting position in the office. The measuring points for internal wall temperature were located on the central line of the east and west wall, and the points of the north wall are 30 cm away from the west wall with the height of 50 cm, 75 cm, 100 cm, 125 cm, 150 cm, 175 cm, 200 cm.
section_name: Field Test

section_name: Simulation and Assumptions

section_content: Since there is no ventilation equipment and the door and windows are closed during working hours, the office can be regarded as a closed cavity. Thus, when the indoor natural convection heat transfer was calculated, the Boussinesq approximation for thermal buoyancy was used to deal with the buoyancy term caused by the temperature difference Temperature range: -20 to 70 • C; accuracy: ±0. 21 • C; humidity range: 1-95%; accuracy: ±2. 5%.
section_name: Governing Equations

section_content: Testo 830-T4 handheld infrared thermometer:

Temperature range: -30 to 350 • C; accuracy: ±2%. 

The field measurements were carried out based on the method recommended in ASHRAE 55-2017 [27]. More specifically, for the outdoor environmental parameters, the location of the measuring points was open ground near the office building at a height of 1. 1 m without direct sunlight. For the indoor environmental parameters, a measuring point at a height of 1. 1 m in an occupied office room was selected as the representative point for the measurement of air temperature and relative humidity because the occupants were usually in a sitting position in the office. The measuring points for internal wall temperature were located on the central line of the east and west wall, and the points of the north wall are 30 cm away from the west wall with the height of 50 cm, 75 cm, 100 cm, 125 cm, 150 cm, 175 cm, 200 cm.
section_name: 2.

section_content: Since there is no ventilation equipment and the door and windows are closed during working hours, the office can be regarded as a closed cavity. Thus, when the indoor natural convection heat transfer was calculated, the Boussinesq approximation for thermal buoyancy was used to deal with the buoyancy term caused by the temperature difference [28]. The Boussinesq approximation hypothesis consists of neglected viscous dissipation in fluids, constant physical properties, and fixed density (density variations are considered only for terms related to momentum equation and volume force). The fluid was in compressible flows in a steady-state so that the air density ρ in the gravity terms can be expressed as [28] :

where ρ c is air density at corresponding T c, kg/m 3 ; α is the coefficient of cubical expansion; and T and T c are the air temperature and reference air temperature, • C. 

Reynolds-averaged Navier-Stokes (RANS) was used to model the conservation of momentum energy and mass. Based on the Boussinesq vortex viscosity hypothesis, the basic equations applicable to natural convection heat transfer problems are as follows [29] :

where the u i is the air velocity component in point x i, m/s and x i is the distance of the point x i for origin, m. 

where u j is the air velocity component in point x j, m/s; P is the partial pressure of water vapor in moist air, Pa; µ is dynamic viscosity, (N•s)/m 2 ; µ t is turbulent viscosity, (N•s)/m 2 ; g i is the gravitational acceleration vector, m/s 2 ; and β is the thermal expansion coefficient, K -1. The energy equation is as follows [30] :

where C p is the specific heat of the air, J/(kg•K) and λ e f f is the effective thermal conductivity (W/m• • C), which can be represented as:

where Pr t is turbulent Prandtl number and λ l is the laminar effective thermal conductivity, W/m• • C.
section_name: Governing Equations

section_content: According to the computational model and theoretical analysis, the indoor airflow pattern of the office room is generally of natural convection. The Rayleigh number is usually used as a criterion to distinguish laminar flow from turbulent flow [31]. The numerical simulation method is based on the Rayleigh-averaged equation, and the indoor zero-equation model is chosen for this model. The turbulent viscosity µ t is expressed by Equation ( 6 ) [32]. 

where v is the local mean velocity, m/s and l is the nearest distance from the wall, m. Moreover, the indoor airflow has turbulence characteristics, while the floor, window, and walls are the starting areas of the turbulent flow, which is a laminar or quasi-laminar zone. Hence, the high-Reynolds-number turbulence model is not applicable. Unlike the solution of the high-Reynolds-number turbulence model, the wall function model combines the physical quantity of the wall and the variables of the turbulence development zone in this computational model. The non-dimensional parameters u + and y + can give an accurate description of the issue, which can be calculated by:

where u is the tangential velocity of the wall, m/s; u τ is the friction velocity of the wall, m/s; ρ w is the density of the wall, kg/m 3 ; and y is the normal distance from the wall to the first layer of mesh, m. The results show that the y + was within the range of 30 to 150. 

In this study, air does not participate in the radiation heat transfer in space, and the surfaces are assumed to be diffused. In addition, the radiation model is used to study the radiation heat transfer of the enclosure structure, the radiator, and the human body when the occupier is in different positions. The discrete ordinate radiation model is as follows [33] :

where the r is position vector; s is a directional vector; s is scattered direction; s is the length path, m; a is absorption coefficient; n is convert coefficient; σ s is scattering coefficient; σ is the Stefan-Boltzmann constant, 5. 672 × 10 -8 W/m 2 K 4 ; I is radiation intensity depending on location ( r ) and direction ( s ), W/m 2 ; Φ is phase function; and Ω is the spatial angle. 

The solver of the CFD model used the finite volume method, because it constructs discrete equations from a physical standpoint, while each discrete equation defines the conservation of a physical variable for a finite volume. Hence, the equations are always conservational. In addition, the SIMPLE algorithm that is sequentially close to the approximate solutions of linear algebraic equations is used.
section_name: Solution Methods

section_content: To investigate the indoor thermal environment characteristics in Lhasa, a 3 dimensional computational office model was built based on the real scenario. The geometry and surface mesh are shown in Figure 5, and the model was constructed and computed in Fluent Airpak 3. 0. 16. The computational model consists of walls, windows, seated persons, standing persons, lights, etc., and the detailed parameters of each item are shown in Table 3. The default person modules were used as the standing and sitting person models in this study. The numbers in Table 3 correspond to Figure 5. In addition, the computers were hidden under the desks. In order to determine the suitable meshing scheme, the gird independence tests were carried out by different numbers of elements. Five grid sizes were considered: 0. 1 m, 0. 07 m, 0. 05 m, 0. 04 m, and 0. 03 m. Furthermore, the temperatures of points (1. 95, 1. 1, 3), (1. 95, 1. 1, 3. 45) and (1. 95, 1. 1, 3. 9) are connected in Figure 6. The temperature of the three points changes little when the grid decreases from 0. 05 m to 0. 03 m. This shows that the impact of grid size on the simulation results can be ignored when the grid size is equal to or smaller than 0. 05 m. The results showed that the computational domain was divided into 279,146 small hexahedron cells with finite volumes under optimal grid size. An acceptable In order to determine the suitable meshing scheme, the gird independence tests were carried out by different numbers of elements. Five grid sizes were considered: 0. 1 m, 0. 07 m, 0. 05 m, 0. 04 m, and 0. 03 m. Furthermore, the temperatures of points (1. 95, 1. 1, 3), (1. 95, 1. 1, 3. 45) and (1. 95, 1. 1, 3. 9) are connected in Figure 6. The temperature of the three points changes little when the grid decreases from 0. 05 m to 0. 03 m. This shows that the impact of grid size on the simulation results can be ignored when the grid size is equal to or smaller than 0. 05 m. The results showed that the computational domain was divided into 279,146 small hexahedron cells with finite volumes under optimal grid size. An acceptable grid quality was obtained in which the elements less than 0. 15 were severely distorted (Figure 6 ). The convergence criteria were the solutions of nonlinear equations and relative errors with tolerances from 1 × 10 -7 to 1 × 10 -3. The residuals of momentum were below 1× 10 -3. The simulation converged in 1123 iterations. In order to determine the suitable meshing scheme, the gird independence tests were carried out by different numbers of elements. Five grid sizes were considered: 0. 1 m, 0. 07 m, 0. 05 m, 0. 04 m, and 0. 03 m. Furthermore, the temperatures of points (1. 95, 1. 1, 3), (1. 95, 1. 1, 3. 45) and (1. 95, 1. 1, 3. 9) are connected in Figure 6. The temperature of the three points changes little when the grid decreases from 0. 05 m to 0. 03 m. This shows that the impact of grid size on the simulation results can be ignored when the grid size is equal to or smaller than 0. 05 m. The results showed that the computational domain was divided into 279,146 small hexahedron cells with finite volumes under optimal grid size. An acceptable grid quality was obtained in which the elements less than 0. 15 were severely distorted (Figure 6 ). The convergence criteria were the solutions of nonlinear equations and relative errors with tolerances from 1 × 10 -7 to 1 × 10 -3. The residuals of momentum were below 1× 10 -3. The simulation converged in 1123 iterations.
section_name: Computational Modeling

section_content: Lhasa is 3650 m above sea level. The atmospheric pressure and air density are 64. 504 kpa and 0. 764 kg/m 3, respectively. In addition, when setting the solar radiation, geographic coordinates of 91 • 13 east longitude and 29 • 66 north latitude were filled in. The outdoor environment is set according to the measured parameters. However, the operating environment parameters were changed when the computational model was used to simulate the indoor thermal environment on a typical-weather day in Lhasa or Beijing. The boundary of the enclosure structure is a wall, which reflects the influence of outdoor temperature through outdoor temperature and convective heat transfer coefficient, and adopts a non-slip and standard wall function wall treatment method. The detailed thermal properties of the building's structure are displayed in Table 2 and the parameters of the windows are shown in Table 4. According to the results of the field survey, the influence of curtains is ignored when calculating the heat transfer process of the outer window. The boundary condition of the indoor heat source is block; the heat dissipation and temperature are shown in Table 3.
section_name: Operating Environment and Boundary Conditions

section_name: Results

section_content: According to the measured data of the office room, the outdoor air temperature in Lhasa ranged between 0. 05 • C and 17. 96 • C, and the indoor air temperature between 17. 01 • C and 24. 05 • C. It is indicated in Figure 7 that, similar to the outdoor air temperature, the indoor air temperature rose from morning to afternoon and fell at about 5 p. m., but the daily range of the outdoor air temperature was much greater than that of the indoor temperature. This is mainly due to the thermal inertia of the building, the use of radiators, and the fact that the windows and door were kept almost closed. Moreover, the average indoor air relative humidity stayed stable, around 21%, while the outdoor air relative humidity fluctuated between 6. 11% to 77. 21%. 

ing environment parameters were changed when the computational model was used to simulate the indoor thermal environment on a typical-weather day in Lhasa or Beijing. The boundary of the enclosure structure is a wall, which reflects the influence of outdoor temperature through outdoor temperature and convective heat transfer coefficient, and adopts a non-slip and standard wall function wall treatment method. The detailed thermal properties of the building's structure are displayed in Table 2 and the parameters of the windows are shown in Table 4. According to the results of the field survey, the influence of curtains is ignored when calculating the heat transfer process of the outer window. The boundary condition of the indoor heat source is block; the heat dissipation and temperature are shown in Table 3.
section_name: Measured Data and Analysis

section_content: According to the measured data of the office room, the outdoor air temperature in Lhasa ranged between 0. 05 °C and 17. 96 °C, and the indoor air temperature between 17. 01 °C and 24. 05 °C. It is indicated in Figure 7 that, similar to the outdoor air temperature, the indoor air temperature rose from morning to afternoon and fell at about 5 p. m., but the daily range of the outdoor air temperature was much greater than that of the indoor temperature. This is mainly due to the thermal inertia of the building, the use of radiators, and the fact that the windows and door were kept almost closed. Moreover, the average indoor air relative humidity stayed stable, around 21%, while the outdoor air relative humidity fluctuated between 6. 11% to 77. 21%. height, and the largest wall temperature difference among the locations from 50 cm to 200 cm was 2. 6 ℃. The average internal wall temperatures of the east wall, west wall, and north wall were 15. 53 ℃, 14. 6 ℃, and 14. 5 ℃, respectively. However, the temperature of the window was not measured. The wall temperature difference could be a hint of indoor air stratification.
section_name: Measured Data and Analysis

section_content: To examine the accuracy of the computational model, the results from March 14 were compared. The input ambient environment parameters are the measured data. The solar load model parameters include local time and position. The sunshine fraction and ground reflectance were 1 and 0. 2, respectively. Figure 9
section_name: Model Verification

section_content: To examine the accuracy of the computational model, the results from 14 March were compared. The input ambient environment parameters are the measured data. The solar load model parameters include local time and position. The sunshine fraction and ground reflectance were 1 and 0. 2, respectively. Figure 9 presents the relative errors of air temperature on the test point between the simulation and the experiment at different times. The measured air temperatures at 11:30, 14:30, and 17:30 were 19. 32 • C, 22. 79 • C, and 23. 67 • C, respectively. The simulated results agreed well with the measured data. The calculated relative errors were within 2%. The west wall temperatures are shown in Figure 10 ; the simulated and measured wall temperatures had the same change trend, with the height rising from 50 cm to 200 cm. The simulated wall temperature was higher than that of the measured temperature by 2. 62 • C. This was because the thermal inertia of the envelope was not considered in the simulation. The actual ground temperature was lower, and the measured temperature stratification was more obvious. Nevertheless, the relative error of the wall temperature at different times was less than 15%. These results illustrate that the developed CFD model can reflect the indoor temperature with sufficient accuracy. 

respectively. The simulated results agreed well with the measured data. The calculated relative errors were within 2%. The west wall temperatures are shown in Figure 10 ; the simulated and measured wall temperatures had the same change trend, with the height rising from 50 cm to 200 cm. The simulated wall temperature was higher than that of the measured temperature by 2. 62 ℃. This was because the thermal inertia of the envelope was not considered in the simulation. The actual ground temperature was lower, and the measured temperature stratification was more obvious. Nevertheless, the relative error of the wall temperature at different times was less than 15%. These results illustrate that the developed CFD model can reflect the indoor temperature with sufficient accuracy.
section_name: Model Verification

section_content: The simulated air temperature distribution at 2:30 pm on March 14th is shown in Figure 11. The horizontal temperature field at the height of 0. 6 m is not evenly distributed. On the one hand, the air temperature was affected by heat sources, resulting in higher air temperature near the heat source areas. On the other hand, with the heat fluxes from solar radiation and radiators, the indoor air temperature in the southern area was higher than in the northern area. According to the simulation results, the air temperature had obvious stratification in the vertical direction, and the temperature rises with the increase in height, which had the same trend as the internal wall temperature. In addition, the temperature of the south-facing window was much higher than that of the north-facing window. respectively. The simulated results agreed well with the measured data. The calculated relative errors were within 2%. The west wall temperatures are shown in Figure 10 ; the simulated and measured wall temperatures had the same change trend, with the height rising from 50 cm to 200 cm. The simulated wall temperature was higher than that of the measured temperature by 2. 62 ℃. This was because the thermal inertia of the envelope was not considered in the simulation. The actual ground temperature was lower, and the measured temperature stratification was more obvious. Nevertheless, the relative error of the wall temperature at different times was less than 15%. These results illustrate that the developed CFD model can reflect the indoor temperature with sufficient accuracy.
section_name: Simulation Results

section_content: The simulated air temperature distribution at 2:30 pm on March 14th is shown in Figure 11. The horizontal temperature field at the height of 0. 6 m is not evenly distributed. On the one hand, the air temperature was affected by heat sources, resulting in higher air temperature near the heat source areas. On the other hand, with the heat fluxes from solar radiation and radiators, the indoor air temperature in the southern area was higher than in the northern area. According to the simulation results, the air temperature had obvious stratification in the vertical direction, and the temperature rises with the increase in height, which had the same trend as the internal wall temperature. In addition, the temperature of the south-facing window was much higher than that of the north-facing window.
section_name: Simulation Results

section_content: The simulated air temperature distribution at 2:30 pm on 14 March is shown in Figure 11. The horizontal temperature field at the height of 0. 6 m is not evenly distributed. On the one hand, the air temperature was affected by heat sources, resulting in higher air temperature near the heat source areas. On the other hand, with the heat fluxes from solar radiation and radiators, the indoor air temperature in the southern area was higher than in the northern area. According to the simulation results, the air temperature had obvious stratification in the vertical direction, and the temperature rises with the increase in height, which had the same trend as the internal wall temperature. In addition, the temperature of the south-facing window was much higher than that of the north-facing window.
section_name: Simulation Results

section_content: Beijing and Lhasa are both located in cold areas, with an altitude of 43. 5 m and geographic coordinates of 116 °25 ′ east longitude and 39 °54 ′ north latitude. A comparison between these two cities can easily explain the difference in indoor thermal environment between the plateau and plain cold areas. Since the model was verified by the measured data in Section 3, it can be used to conduct a more in-depth analysis of the indoor thermal environment in Lhasa and Beijing. For example: What is the indoor thermal environment on a typical day in winter? It can be seen from the results of the measured and simulated temperatures that the vertical and horizontal indoor temperature differences are relatively large. Therefore, this section mainly analyzes the non-uniform distribution of indoor temperature comparatively.
section_name: Discussion

section_content: The meteorological parameters of typical weather days in Lhasa and Beijing in January are shown in Figure 12 [36]. As can be seen from Figure 12a, the outdoor air temperature is below 0 ℃ in the forenoon in the two cities. Furthermore, the temperature in Lhasa ranges from -7. 6 ℃ to 5. 8 ℃, with a diurnal range of 13. 4 ℃, but the temperature in Beijing varies from -6. 6 ℃ to 1. 2 ℃ with a diurnal range of 7. 6 ℃. The time zone in Lhasa is two hours later than that in Beijing. Thus, the solar radiation intensity in Lhasa is two hours behind the actual values in Figure 12b. This is almost identical to Beijing, but Lhasa's direct solar radiation intensity is significantly higher than that of Beijing. As a result, the global solar radiation intensity in Lhasa is higher. At midday, the global solar radiation intensity in Lhasa is 696 W/m 2, but in Beijing it is 411 W/m 2. Consequently, the outdoor temperature in Lhasa exceeds that of Beijing in the afternoon, resulting in a higher daily average temperature in Lhasa. Due to the time difference, the highest measured temperatures in Lhasa and Beijing occurred at 17:00 and 15:00, respectively.
section_name: Typical-Weather Day in Lhasa and Beijing

section_content: Beijing and Lhasa are both located in cold areas, with an altitude of 43. 5 m and geographic coordinates of 116 • 25 east longitude and 39 • 54 north latitude. A comparison between these two cities can easily explain the difference in indoor thermal environment between the plateau and plain cold areas. Since the model was verified by the measured data in Section 3, it can be used to conduct a more in-depth analysis of the indoor thermal environment in Lhasa and Beijing. For example: What is the indoor thermal environment on a typical day in winter? It can be seen from the results of the measured and simulated temperatures that the vertical and horizontal indoor temperature differences are relatively large. Therefore, this section mainly analyzes the non-uniform distribution of indoor temperature comparatively.
section_name: Discussion

section_content: The meteorological parameters of typical weather days in Lhasa and Beijing in January are shown in Figure 12 [36]. As can be seen from Figure 12a, the outdoor air temperature is below 0 • C in the forenoon in the two cities. Furthermore, the temperature in Lhasa ranges from -7. 6 • C to 5. 8 • C, with a diurnal range of 13. 4 • C, but the temperature in Beijing varies from -6. 6 • C to 1. 2 • C with a diurnal range of 7. 6 • C. The time zone in Lhasa is two hours later than that in Beijing. Thus, the solar radiation intensity in Lhasa is two hours behind the actual values in Figure 12b. This is almost identical to Beijing, but Lhasa's direct solar radiation intensity is significantly higher than that of Beijing. As a result, the global solar radiation intensity in Lhasa is higher. At midday, the global solar radiation intensity in Lhasa is 696 W/m 2, but in Beijing it is 411 W/m 2. Consequently, the outdoor temperature in Lhasa exceeds that of Beijing in the afternoon, resulting in a higher daily average temperature in Lhasa. Due to the time difference, the highest measured temperatures in Lhasa and Beijing occurred at 17:00 and 15:00, respectively.
section_name: Typical-Weather Day in Lhasa and Beijing

section_content: The RTA value depends on the view factor between the human body and the surrounding surface. Different parts of the human body, the surfaces of the building envelope, and the heat source will produce different net radiative heat transfers, resulting in inconsistent radiant temperatures of various parts of the human body surface. The difference between the radiation temperature of the plane of the two opposite sides for the human body is RTA. Hence, six indoor points were selected as the human body position to study the non-uniform thermal distribution under the influence of a unique climate, as shown in Figure 13. In this model, the seated persons are facing east (X-axis reverse). The numerical model settings of the building envelope, radiator, and lamps are consistent with Section 2. 2.
section_name: Human Exposure Positions

section_content: Figure 14 shows the variation in the average air temperature and radiation temperature on a horizontal surface at a height of 0. 6 m. The indoor air temperature and radiation temperature have the same changing trend, and the difference between the two is small. Corresponding to the change in outdoor air temperature, in the morning, the air temperature in Beijing is higher, and a greater daily air temperature difference is evident in Lhasa. However, the maximum air temperature in the two cities does not exceed 15 • C, which is lower than the requirements of the "Design Standard for Heating and Ventilation of Civil Buildings for Tibet" (DBJ 540002-2016) [37] and the "Design Code for Heating Ventilation and Air Conditioning of Civil Buildings" (GB 50736-2012) [26], which are 18-24 • C. Moreover, the indoor air temperature is also lower than the neutral temperature that is proposed in existing research conclusions [10, 12, 13]. This shows that the current heating plan cannot make rooms reach a comfortable state. Even in March, indoor personnel still need to wear down jackets to keep warm in the morning. The average minimum and maximum indoor temperatures at a height of 0. 6 m in Lhasa (Ta1) are 5. 98 • C and 14. 7 • C, respectively. It is worth noting that the outdoor temperature in Lhasa is up to 5. 7 • C higher than in Beijing, and the solar radiation is also stronger. However the indoor temperature in Lhasa is up to 4. 17 • C higher than that of Beijing. Furthermore, the standard deviation (SD) of the radiation temperature is between 2. 5 • C and 3 • C, and the SD of the air temperature is around 1. 5 • C. In addition, the SD of the indoor air temperature at a height of 0. 6 m in Beijing (Ta2) is similar to that of Ta1, but the SD of the radiation temperature at a height of 0. 6 m (Tr2) is higher than that of Lhasa (Tr1) afternoon. Solar radiation has a greater impact on the radiant temperature of different indoor areas, and further analysis of non-uniform indoor radiant temperature is necessary. Therefore, the design load of central heating public buildings should match the actual situation.
section_name: Average Temperature

section_content: Figure 14 shows the variation in the average air temperature and radiation temperature on a horizontal surface at a height of 0. 6 m. The indoor air temperature and radiation temperature have the same changing trend, and the difference between the two is small. Corresponding to the change in outdoor air temperature, in the morning, the air temperature in Beijing is higher, and a greater daily air temperature difference is evident in Lhasa. However, the maximum air temperature in the two cities does not exceed 15 ℃, which is lower than the requirements of the "Design Standard for Heating and Ventilation of Civil Buildings for Tibet" (DBJ 540002-2016) [37] and the "Design Code for Heating Ventilation and Air Conditioning of Civil Buildings" (GB 50736-2012) [26], which are 18-24 ℃. Moreover, the indoor air temperature is also lower than the neutral temperature that is proposed in existing research conclusions [10, 12, 13]. This shows that the current heating plan cannot make rooms reach a comfortable state. Even in March, indoor personnel still need to wear down jackets to keep warm in the morning. The average minimum and maximum indoor temperatures at a height of 0. 6 m in Lhasa (Ta1) are 5. 98 ℃ and 14. 7 ℃, respectively. It is worth noting that the outdoor temperature in Lhasa is up to 5. 7 ℃ higher than in Beijing, and the solar radiation is also stronger. However the indoor temperature in Lhasa is up to 4. 17 ℃ higher than that of Beijing. Furthermore, the standard deviation (SD) of the radiation temperature is between 2. 5 ℃ and 3 ℃, and the SD of the air temperature is around 1. 5 ℃. In addition, the SD of the indoor air temperature at a height of 0. 6 m in Beijing (Ta2) is similar to that of Ta1, but the SD of the radiation temperature at a height of 0. 6 m (Tr2) is higher than that of Lhasa (Tr1) afternoon. Solar radiation has a greater impact on the radiant temperature of different indoor areas, and further analysis of non-uniform indoor radiant temperature is necessary. Therefore, the design load of central heating public buildings should match the actual situation. Even in the morning during March and wearing leather boots, the feet still feel uncomfortably cold. However, the body does not feel cold. In addition, as shown in Figure 11, the phenomenon of the thermal stratification of indoor air is obvious. Therefore, the VATD is a problem worth analyzing. According to the model configuration provided by the "Evaluation Standard for the Indoor Thermal Environment in Civil Buildings" (GB/T 50785-2012) [38] and ISO 7730-2005 [23], when the human body is in a sitting position, the air temperature around the head and ankles are set to 1. 1 m and 0. 1 m, respectively. In addition, the percentage of dissatisfaction 2 caused by the VATD is calculated according to [23, 38] : 100 Even in the morning during March and wearing leather boots, the feet still feel uncomfortably cold. However, the body does not feel cold. In addition, as shown in Figure 11, the phenomenon of the thermal stratification of indoor air is obvious. Therefore, the VATD is a problem worth analyzing. According to the model configuration provided by the "Evaluation Standard for the Indoor Thermal Environment in Civil Buildings" (GB/T 50785-2012) [38] and ISO 7730-2005 [23], when the human body is in a sitting position, the air temperature around the head and ankles are set to 1. 1 m and 0. 1 m, respectively. In addition, the percentage of dissatisfaction LPD 2 caused by the VATD is calculated according to [23, 38] :

where ∆t a,v is the VATD between head and ankle, • C. ISO 7730-2005 stipulates that the above formula should only be used ∆t a,v < 8 • C [23]. 

According to the simulation results, the air temperature around the ankle and head of the human body at different positions is similar to the changing trend in indoor temperature at a height of 0. 6 m. Lhasa has a wider range of temperature changes, and the air temperature at various locations in the morning is lower than in Beijing. Since the outdoor temperature of Lhasa and Beijing in the morning is below 0 • C, the indoor temperature is very low, even near the radiator (position 1 (P1) and position 6 (P6)). When the solar radiation and outdoor temperature are low, the building exterior envelopes are in a heat-loss state. The temperature in different positions is low, and the temperature difference is small. Because this area is far away from the radiators, the air temperature around the ankle in position 3 (P3) and position 4 (P4) is lower than in other positions. The lowest temperatures near the ankle in Lhasa and Beijing are 3. 71 • C and 4. 33 • C, respectively. When the solar radiation is strong, the indoor temperature is relatively high. The temperature around the head in position 2 (P2) can reach 15. 8 • C. The temperature difference between different positions in the room is obvious. This is because part of the indoor area receives heat from solar energy. In addition, the south outer wall (window) is in a state of heat gain, while the north outer wall (window) is still in a state of heat loss. 

Figure 15a calculates the VATDs between the ankle and the head, and Figure 15b lists the corresponding LPD 2. The vertical temperature difference between Lhasa and Beijing is less than 3 • C most of the time, and the LPD 2 is also within 5%. Generally, P3 has the largest vertical temperature difference due to its proximity to the west and north outer wall. The maximum vertical temperature differences between Lhasa and Beijing are 3. 12 • C and 3. 2 • C, and the corresponding dissatisfaction rates are 4. 87% and 5. 21% (in Figure 15b ). According to the requirements of ISO 773-2005 [23], the VATD should be less than 3 • C (category B) or 4 • C (category C), and the corresponding LPD 2 is 10%. Moreover, the strictest limit for LPD 2 is 10% in GB/T 50785-2012 [38]. Therefore, the indoor vertical temperature difference in Lhasa and Beijing meets the requirements of ISO 773-2005 and GB/T50785-2012. However, ASHRAE 55-2017 [27] stipulates that LPD 2 should be less than 5%. 

℃ (category B) or 4 ℃ (category C), and the corresponding 2 is 10%. Moreover, the strictest limit for 2 is 10% in GB/T 50785-2012 [38]. Therefore, the indoor vertical temperature difference in Lhasa and Beijing meets the requirements of ISO 773-2005 and GB/T50785-2012. However, ASHRAE 55-2017 [27] stipulates that 2 should be less than 5%.
section_name: Average Temperature

section_content: According to the concept of RTA proposed by Fanger [18, 19], human RTA can be evaluated in three orthogonal directions: chest and back (X-axis), head and ankles (Y-axis), and left and right arms (Z-axis). Equation ( 11 ) can calculate the dissatisfaction rate caused by asymmetric radiation due to cold walls (windows). However, this equation is only applicable to the asymmetry of the left and right arms [23]. 

where 4 is the percentage of dissatisfied people due to radiant asymmetry, %. ∆ is RTA, ℃. 

The RTA of the human body in different directions is shown in Figure 16. In general, the radiant asymmetry in Lhasa is more obvious than in Beijing. The RTA in the X direction is the average radiant temperature difference between the chest and the back. A negative value indicates that the back temperature of the human body is higher. As shown in Figure 16a, the RTA of indoor position 5 (P5) in the two cities is the largest. The RTA in Lhasa ranges from 1. 1 ℃ to 2. 53 ℃ and the RTA in Beijing ranges from 1. 6 ℃ to 2. 8 ℃. In addition, the radiation temperature of the back at P6 is higher than that of the chest. The RTA value of Lhasa in this position is greater, with a maximum value of -1. 6 ℃. On the
section_name: Radiant Temperature Asymmetry

section_content: According to the concept of RTA proposed by Fanger [18, 19], human RTA can be evaluated in three orthogonal directions: chest and back (X-axis), head and ankles (Y-axis), and left and right arms (Z-axis). Equation (11) can calculate the dissatisfaction rate caused by asymmetric radiation due to cold walls (windows). However, this equation is only applicable to the asymmetry of the left and right arms [23]. 

where LPD 4 is the percentage of dissatisfied people due to radiant asymmetry, %. ∆t pr is RTA, • C. The RTA of the human body in different directions is shown in Figure 16. In general, the radiant asymmetry in Lhasa is more obvious than in Beijing. The RTA in the X direction is the average radiant temperature difference between the chest and the back. A negative value indicates that the back temperature of the human body is higher. As shown in Figure 16a, the RTA of indoor position 5 (P5) in the two cities is the largest. The RTA in Lhasa ranges from 1. 1 • C to 2. 53 • C and the RTA in Beijing ranges from 1. 6 • C to 2. 8 • C. In addition, the radiation temperature of the back at P6 is higher than that of the chest. The RTA value of Lhasa in this position is greater, with a maximum value of -1. 6 • C. On the one hand, the back receives more heat radiation from the radiators. On the other hand, the back receives solar energy. The RTA value of the third position in Beijing is similar to that of Lhasa, but after 16:00, the radiation temperature of the chest is lower than that of the back. It may be that the outdoor temperature drops, causing the cold radiation of the north outer window to be stronger than that of the west outer wall. The RTA value in other positions is positive because the person is facing the east inner wall and their back is facing the west outer wall. Taking the simulation result at 10:00 in Lhasa as an example, it is possible to further analyze the human body radiation temperature at indoor location 4. Figure 17 The RTA in the Z direction refers to the average radiation difference between the left arm and the right arm. The distribution trend in RTA at different positions in Lhasa and Beijing is the same most of the time, but there is a big difference at 9:00 and 18:00, as shown in Figure 16b. To take P1 as an example, the average radiant temperature of the left arm in Lhasa is higher than that of the right arm at 9:00, while the temperature of the right arm is higher in Beijing. At this juncture, the solar radiation intensity in Lhasa is 0 W/m 2, and the outdoor air temperature is as low as -6. 7 • C. The south outward window is the heat loss surface, which can be regarded as a cold radiant panel. The human body receives more cold radiation from the external window than heat radiation from the radiator. However, the total solar radiation intensity at 9:00 in Beijing is 125 W/m 2 and the outdoor air temperature is -5. 1 • C. The right arm of the human body in P1 receives more heat than the left arm. At 18:00, the RTA of Lhasa and Beijing at location 1 are -3. 5 • C and -0. 57 • C, respectively. In addition, the average radiation temperature of the right arm in Beijing is higher after 16:00. This is because the sunrise and sunset times are different, since they are affected by the time difference between the two cities. Overall, the RTA in the Z direction of Lhasa is higher, reaching -3. 95 • C, with and LPD 4 of 0. 52% at P4. 

The RTA in the Y direction refers to the average radiant temperature difference between the head and ankle, as shown in Figure 16c. The ankle receives more heat than the head when the human body is seated in P1 and P6. Therefore, the radiation temperature of the ankle is higher than that of the head, resulting in a negative RTA in the Y direction. Lhasa has the largest RTA at P1 in the Y direction, which can reach -2. 79 • C, with LPD 4 of 0. 35%. The RTA of Beijing at P1 is relatively small, with a maximum of -1. 86 • C. When individuals are in P2 and P5, they are far away from the outer window and receive less cold radiation. Therefore, the average radiant temperature of the head is relatively high, and the RTA in Lhasa is positive in the Y direction. In addition, the average radiant temperature of the head in P4 in Lhasa is lower than that of the ankle when the solar radiation is strong. This can be explained by the fact that the radiation temperature on the side of the head near the radiators and the south outer window is higher, and the radiation temperature on the other sides is lower, especially the left side near the north outer window. Nevertheless, the ankle is far from the outer window, and its average radiation temperature is relatively high. Due to the relatively weak solar radiation in Beijing, the influence of solar energy on indoor radiant temperature is small. Coupled with the influence of outdoor temperature, the average radiant temperature of P4 in Beijing is lower than that in Lhasa. At noon, the average radiant temperature of the head and ankle of P4 in Lhasa are 12. 6 • C and 13. 56 • C, while the average radiant temperatures of Beijing at that this position are 12. 2 • C and 11. 96 • C, respectively. The cold radiation caused by the north outer window on the head has a weaker effect on the average radiation temperature than the indoor thermal stratification. Therefore, the RTA of Beijing P4 in the Y direction is positive. 

Taking the simulation result at 10:00 in Lhasa as an example, it is possible to further analyze the human body radiation temperature at indoor location 4. Figure 17 is the radiation temperature cloud diagram of the front chest and back at indoor location 4 in Lhasa. At this time, the average radiation temperatures of the human chest and back are 9. 54 • C and 7. 32 • C, respectively. Affected by the cold radiation from the north outer window, the radiation temperature of the central axis of the chest and back increases with the decrease in height. In addition, the radiant temperature on the left side of the human body (near the north outer window) is much lower than on the right side, as is clearly manifested by the head, arms, back, and legs. Figure 14 is also a cloud map of the radiation temperature on the left and right sides of the human body at P4 in Lhasa. The difference in radiation temperature between the left and right arms of the human body is more evident. The average radiation temperature of the left arm and right arm are 5. 94 • C and 9. 8 • C, respectively. The radiant temperature of the head and ankles on the same side is little different. However, the low temperature of the north outer window causes the temperature on the left side of the head to be significantly lower than that on the right side, and the maximum radiant temperature difference reaches 8. 73 • C, with an LPD 4 of 2. 66%. 

However, the low temperature of the north outer window causes the temperature on the left side of the head to be significantly lower than that on the right side, and the maximum radiant temperature difference reaches 8. 73 ℃, with an 4 of 2. 66%. In summary, there are certain differences in the indoor thermal environment between Lhasa and Beijing. Under the same envelope structure and heating conditions, the average indoor temperature and average radiant temperature in Lhasa are higher, and the nonuniform of indoor temperature distribution in Lhasa is also more significant. The indoor temperature in both cities is below comfortable levels. The indoor vertical temperature difference between Lhasa and Beijing meets the requirements of ISO 7730-2005 [23] and GB/T50785-2012 [38]. In general, the indoor RTA in Lhasa is larger than in Beijing and lower than the limits set by the current standard. The RTA in the Z direction is the largest. This is due to the large heat loss of the north outer window on the left side of the human body, while the right side can obtain the heat from the radiator and the solar energy passing through the south outer window. Lhasa's position near the north outer window received strong cold radiation, resulting in higher RTA during certain periods. This shows that direct solar radiation through windows helps to improve thermal comfort. Due to thermal stratification, the air temperature at the height of 1. 1 m in Lhasa and Beijing is higher than at 0. 1 m. The average radiant temperature of the head in indoor P1 and P6 is lower than the average radiant temperature of the ankle. In addition, the average radiation temperature of the head in P4 is lower than that of the ankle during periods of strong solar radiation in Lhasa. This shows that the cold radiation from outer windows in the cold plateau area has a greater impact on the indoor thermal environments to which human bodies are subjected, even when the temperature is high during the day.
section_name: Radiant Temperature Asymmetry

section_content: At present, there are two standards for the design and evaluation of indoor thermal environment in Tibet. The regulations on the lighting, ventilation, heat preservation, and heat insulation of building layouts, orientations, and forms can be found in the " Design Standard for Green Building in the Tibet Autonomous Region" (DBJ540001-2018) [39]. Moreover, the standard also requires that indoor thermal comfort should consider dry bulb temperature, air humidity, airflow speed, average radiant temperature difference, and indoor personnel activities and clothing conditions. However, there is no clear description of the concept, calculation method, and limit of the average radiant temperature difference. The "Design Standard for Heating and Ventilation of Civil Buildings for Tibet" (DBJ540002-2016) stipulates that the main room should be 18-24 ℃ and that the wind speed of the heating indoor activity area should be lower than 0. 3 m/s [37]. The national indoor thermal environment design and evaluation standards include GB 50736-2012 and In summary, there are certain differences in the indoor thermal environment between Lhasa and Beijing. Under the same envelope structure and heating conditions, the average indoor temperature and average radiant temperature in Lhasa are higher, and the nonuniform of indoor temperature distribution in Lhasa is also more significant. The indoor temperature in both cities is below comfortable levels. The indoor vertical temperature difference between Lhasa and Beijing meets the requirements of ISO 7730-2005 [23] and GB/T50785-2012 [38]. In general, the indoor RTA in Lhasa is larger than in Beijing and lower than the limits set by the current standard. The RTA in the Z direction is the largest. This is due to the large heat loss of the north outer window on the left side of the human body, while the right side can obtain the heat from the radiator and the solar energy passing through the south outer window. Lhasa's position near the north outer window received strong cold radiation, resulting in higher RTA during certain periods. This shows that direct solar radiation through windows helps to improve thermal comfort. Due to thermal stratification, the air temperature at the height of 1. 1 m in Lhasa and Beijing is higher than at 0. 1 m. The average radiant temperature of the head in indoor P1 and P6 is lower than the average radiant temperature of the ankle. In addition, the average radiation temperature of the head in P4 is lower than that of the ankle during periods of strong solar radiation in Lhasa. This shows that the cold radiation from outer windows in the cold plateau area has a greater impact on the indoor thermal environments to which human bodies are subjected, even when the temperature is high during the day.
section_name: Analysis of Standards

section_content: At present, there are two standards for the design and evaluation of indoor thermal environment in Tibet. The regulations on the lighting, ventilation, heat preservation, and heat insulation of building layouts, orientations, and forms can be found in the " Design Standard for Green Building in the Tibet Autonomous Region" (DBJ540001-2018) [39]. Moreover, the standard also requires that indoor thermal comfort should consider dry bulb temperature, air humidity, airflow speed, average radiant temperature difference, and indoor personnel activities and clothing conditions. However, there is no clear description of the concept, calculation method, and limit of the average radiant temperature difference. The "Design Standard for Heating and Ventilation of Civil Buildings for Tibet" (DBJ540002-2016) stipulates that the main room should be 18-24 • C and that the wind speed of the heating indoor activity area should be lower than 0. 3 m/s [37]. The national indoor thermal environment design and evaluation standards include GB 50736-2012 and GB/T 50785-2012. In GB 50736-2012, the indoor thermal environment with a temperature of 18-22 • C, relative humidity of more than 30%, and wind speed of less than 0. 2 m/s is classified as level I. When the indoor temperature is 22-24 • C and the wind speed is lower than 0. 2 m/s, which can be considered as meeting level II [26]. The evaluation of the overall indoor thermal environment for GB/T 50785-2012 is consistent with GB 50736-2012. Furthermore, the grades and limits of local thermal comfort are supplemented in GB/T 50785-2012, but the impact of radiant asymmetry on local comfort is not included [38]. ISO 7730-2005 and ASHRAE 55-2017 present comprehensive and detailed regulations on local thermal comfort. The limits of the vertical temperature difference and RTA of category A for ISO 7730-2005 are 2 • C and 10 • C (cool wall), and the corresponding percentages of dissatisfied people are 3% and 5%, respectively [23]. In addition, the upper limits of ASHRAE standard 55-2017 for vertical temperature difference (sitting posture) and RTA (cool wall) are 3 • C and 14 • C, respectively, and the corresponding dissatisfaction rate is 5%. However, the RTA limits of the ASHRAE 55 standard are applicable when the ambient temperature is within thermal comfort levels for a lightly clothed person engaged in near-sedentary physical activity. In addition, they may underestimate local discomfort when the temperature is lower than comfortable levels. However, the results presented by Su et al. show that the RTA limit under radiator heating is 1. 6 • C [40]. 

It can be seen that corresponding indoor thermal environment standards in cold plateau areas help to guide the design of buildings and indoor thermal environment parameters. However, the guidance for the design of indoor thermal parameters in the unique climate of the alpine region is not perfect. The indoor thermal parameters are limited, but their particularity is not considered. The non-uniform distribution of indoor temperatures in the cold plateau regions is caused by low outdoor temperatures, strong solar radiation, and heating. However, the local thermal discomfort caused by non-uniform indoor temperatures is not mentioned in the local standards. Moreover, when the indoor temperature is low, the local discomfort of the human body is more sensitive. Although detailed regulations on non-uniform indoor local thermal comfort can be found in ISO 7730 and ASHRAE 55, the limits are set according to certain indoor conditions and the relative position between the human body and the environment. They do not apply to the evaluation of indoor local thermal comfort in Tibetan Plateau areas.
section_name: Analysis of Standards

section_content: A unique indoor thermal environment is created by the climatic characteristics of the Tibetan Plateau, but a comprehensive understanding of the indoor thermal environment cannot be found in relevant research. In this study, the non-uniform distribution of radiant temperature was evaluated and analyzed comparatively. In short, the following conclusions can be drawn:

(1) According to the field measurement results in March, the selected office room stayed within the comfortable temperature range. Furthermore, the temperature of the inner wall increased with the height of the location due to thermal stratification. 

(2) On a typical-weather day in the winter in Lhasa and Beijing, the maximum indoor average temperature at the height of 0. 6 m is 14. 7 • C. The VATD of Lhasa and Beijing are below 3 • C most of the time. However, the asymmetry of the indoor radiant temperature in Lhasa is more notable. The maximum RTA of the human body and the head, respectively, are -3. 95 • C and 8. 73 • C. 

(3) Due to the strong solar radiation and low temperature in the Tibetan Plateau, indoor thermal comfort is greatly affected by the ambient environment. Architectural designs should try to avoid setting north-facing windows in this area to reduce heat loss. 

(4) A comprehensive standard should be compiled based on the characteristics of non-uniform indoor thermal environments in cold plateau areas for building thermal environment definition, calculation, design, and evaluation. 

Moreover, future developments of this research may include:

(1) A study of the indoor thermal environment characteristics of residential buildings. It is severely cold, particularly in winter, in the Tibetan Plateau. 

(2) A subjective questionnaire survey of overall thermal sensation and local thermal sensation at the same time as environmental measurements including air temperature, air relative humidity, air speed, and radiant temperature.
section_name: Conclusions

section_name: Institutional Review Board Statement: Not applicable.

section_content: Funding: This research was funded by " National Natural Science Foundation of China (NSFC), grant number 52008132 ", " Basic and Applied Basic Research Foundation of Guangdong Province, China, grant number. 2019A1515111194", and " Shenzhen Science and Technology Program, grant number. RCBS20200714114921062".

section_content: Informed Consent Statement: Not applicable.
section_name: Institutional Review Board Statement: Not applicable.

section_name: Data Availability Statement:

section_content: The study did not report any data.
section_name: Data Availability Statement:

section_name: Conflicts of

section_content: Author Contributions: Conceptualization, M. W. and P. S. ; methodology, M. W. ; software, M. W. ; validation, M. W. ; formal analysis, P. S. ; investigation, M. W. ; resources, P. S. ; data curation, M. W. ; writing-original draft preparation, M. W. ; writing-review and editing, P. S. ; visualization, M. W. ; supervision, P. S. ; project administration, P. S. ; funding acquisition, P. S. All authors have read and agreed to the published version of the manuscript.