Abstract
Keywords
Introduction
Thermal comfort of vehicle cabin has always been one of the main concerns of car owners, researchers, and the automotive industry, which has a significant influence of passengers’ comfort and driver’s driving. In addition, a large amount of energy is required for automotive heating ventilation and air-conditioning (HVAC) systems to maintain thermal comfort in the cabin. However, the vehicle is characterized by poor thermal insulation and complex space and are susceptible to solar radiation. The convention heat transfer between the vehicle and the environment is greatly increased because of high-speed driving. So the air heat transfer in vehicle cabin is highly non-uniform and transient, which is significantly different from the thermal environment in the building. Therefore, how to scientifically and effectively control the thermal comfort of the cabin and reduce energy waste is an important task of the automotive HVAC system.
Many scholars have studied the thermal environment or thermal comfort in vehicle cabin by numerical simulation or experiment. Chakroun and Al-Fahed 1 studied thermal comfort in cars parked in the sun during the Kuwaiti summer. In addition, the differences of air temperature inside the car between covering the front windshield and the four side windows were studied comparatively. Shojaee et al. 2 analyzed the influence of air supply parameters on the distribution of flow field and temperature field in passenger compartment by computational fluid dynamics (CFD). Zhang et al.3,4 proposed the computing model of the three-dimensional (3D) temperature distribution and flow field in the cabin by the commercial software FLUENT, which was verified with experimental data. The influence of different factors such as the number of passengers, cooling load, and ambient temperature on thermal comfort and energy consumption was also studied. Guan et al.5,6 studied human comfort experimentally by simulating 16 typical climatic conditions in the environment chamber, which was also compared with the highly transient conditions in the cabin. The thermal sensation model, which took into account both physical and psychological factors, was established in the paper. Zhang et al. 7 conducted a study on thermal comfort during the heating period under idling state in the climate of South China, which showed that the mean skin temperature of human body could be used to reflect the human thermal comfort. Chien et al. 8 conducted a study to evaluate the thermal comfort of passengers during cooling. In this research, the temperature field and flow field in the compartment were analyzed by transient numerical analysis and verified by experiments, and the local thermal comfort of passengers was evaluated by the predictive mean vote (PMV) method proposed by Fanger. 9 The results showed that the average temperature decreased rapidly and the flow field became very complicated during the initial cooling period, and began to stabilize after 10 min. The local PMV value indicated that only the feet remained uncomfortable after 20 min. Han and Huang 10 proposed a numerical simulation on an MPV vehicle and analyzed the effects of different solar incident angles, solar radiation intensity, air supply temperature, speed, and the breathing condition of the passenger on thermal comfort. The equivalent temperature was used to evaluate the thermal comfort of human body.
Most of studies on thermal comfort of cabin were conducted under laboratory or outdoor parking conditions, and only a few were conducted under driving conditions.11,12 It is necessary to objectively evaluate the thermal environment and thermal comfort of the passenger compartment under real driving conditions. Moreover, the usage of PMV model of uniform environment and equivalent temperature
Introduction of experiment
Experiment of thermal environment
Due to the uneven and unsteady thermal environment in the cabin, temperature, humidity, wind speed, and mean radiation temperature were measured at multiple positions of the cabin, as shown in Figures 1 and 2. The air temperature and wind speed of the feet, abdomen, and head of the driver were measured at 1st–3rd. The seat temperature of the driver’s hip, waist, and back were measured at 4th–6th. The air temperature and wind speed of the rear passenger’s feet, abdomen, and head were measured at 7th–9th. The seat temperature of the rear passengers’ hip, waist, and back were measured at 10th–12th. Relative humidity in front and back rows were measured at 13th–14th. The mean radiation temperature in the cabin was measured at 15th. The wind speed of outlets of the air conditioner was measured at 16th–19th. There were sensors outside the car, measuring the temperature and humidity of outside air.

Vertical view of measuring point in cabin.

Front view of measuring point in cabin.
Experiment of thermal comfort
Twenty graduate students (10 male students and 10 female students) volunteered to participate in this study as subjects. They all understood the purpose, methods, and procedures of this study. They were asked to stay away from alcohol and strenuous exercise for at least 18 h before the experiment and to wear uniform to avoid clothing differences that might affect test results.
In order to calculate the weighted PMV value and the weighted

Human skin temperature measuring points.
The comfortable questionnaire which is put forward by ASHRAE standard 13 is used to evaluate the human thermal comfort. The seven values are “+3,”“+2,”“+1,”“0,”“−1,”“−2” and “−3,” respectively, reflect the subjective thermal sensation of “hot,”“warm,”“slightly warm,”“neutral,”“slightly cool,”“cool,” and “cold.” This article only conducts a questionnaire on the overall thermal sensation.
Experimental instrument
The main experimental instruments are shown in Table 1.
Main experimental instruments.
Experimental procedure
Before the experiment, each sensor was arranged inside and outside of the vehicle. The temperature sensors need to be shielded by the belt with low solar radiation absorption rate, so as to avoid solar radiation to affect the test accuracy. Volunteers who took part in the experiment rested, sited, read, or walked quietly in office with temperature close to the human thermal comfort zone.
After the arrangement of measuring sensors, all doors of the cabin should be kept open, and shade curtains should be arranged in front of the windshield to avoid interference from solar radiation, so as to keep the thermal environment of the cabin consistent with the external thermal environment. After 1 h, the shade was removed, all doors were closed, and the cabin was exposed in the sun for 1 h. After the exposure test, two testers including a driver and a passenger sitting in the back row got into the car at the same time. And then they immediately turned on the air conditioner under automatic mode. The air outlet direction was adjusted to the level. The temperature, humidity, and air speed data were collected every 5 s, the globe temperature data were collected every 1 min, and the subjective surveys of human comfort were conducted every 5 min.
During the experiment, the air-conditioning was set to run automatically, internal circulation, cooling, and 26°C, no matter the car was under idle or driving state.
Data processing
Mean skin temperature
The mean skin temperature of the upper body and the lower body need to be calculated, respectively, because of the calculation of the thermal evaluation index. The research of Butera 14 on human skin temperature is referred for the calculation of average skin temperature, as shown in formulas (1) and (2)
Weighted PMV model
According to the PMV–PPD model proposed by Fanger, 9 the calculation method of PMV is as follows
where
According to the area-weighted PMV model proposed by Ingersoll et al., 15 PMV of each part of the human body is calculated separately, and then weighted based on the formula (11)
Weighted equivalent temperature teq model
According to SAE J2234-2012, 16 if there is an imaginary closed space with zero wind speed, where the convective and radiant heat transfer between human body and the space are consistent with the real environment, the temperature of the imaginary space is the equivalent temperature of the real environment. The calculation method is as follows
where
According to the calculation of weighted equivalent temperature in SAE J2234-2012,
16
Since

Schematic of relationship between
Analysis and discussion
Analysis of thermal environment
Influence of position on air temperature
The experiment lasts for 120 min. The first 60 min are the preheating stage, and the last 60 min are the cooling and cooling stage. It can be seen from Figures 5 and 6 that the air temperature in the cabin under idle condition is significantly stratified in the vertical direction.

Air temperature change of driver’s position (idle).

Air temperature change of passenger’s position (idle).
During the preheating stage, the air temperature of the driver and rear passenger from top to bottom are distributed from high to low. The maximum temperature difference between the head and the feet occurs on 20 min, which of the driver and rear passenger are 8.6°C and 7.6°C, respectively.
The main reason of the significantly stratified temperature in the vertical direction is that the upper air is strongly influenced by solar radiation. In addition, the air speed of the top is higher than that of the bottom, as is the coefficient of heat transfer with the environment, so the heat flux density of the top is greater than that of the bottom. At the same time, the mean rise rate of the upper air temperature is much higher than that of the lower air temperature due to the heat storage of interior decorations such as car seats and carpets. The passenger position in the back row is slightly less affected by solar radiation than that in the front row because of vehicle structure. Therefore, the heating rate and maximum temperature difference in the back row are slightly less than that in the driver’s position.
During the preheating phase, the vertical temperature difference reaches the peak value first and drops slightly later until the downward trend flattens out, while the lower air temperature rises after 20 min of preheating, the upward trend weakened and finally flattened out. It is the reason that after 20 min, the influence of the difference of temperature rising rate between the upper and lower layers is gradually smaller than that of heat exchange between the upper and lower layers. The heat of the upper air transferred to the lower air is more than that gained from outside. As the temperature difference between the upper and lower air decreases, the heat gradually becomes balanced, and finally, the temperature difference maintains constant.
In the cooling stage, the air temperature of the cabin is obviously stratified, which is different from the preheating stage and is also different between the front and rear. The stratification of air temperature at the driver’s position is the same as that at the preheating stage, but the temperature difference is smaller because of the cooling effect of the air-conditioning. Due to the horizontal setting of the air outlet, the cold air diffuses up and down after blowing to the driver’s body surface, so the upper and the lower air at the driver’s position are effectively cooled, and the cooling rate is relatively close. The cold air can only blow to the upper level of the rear passenger, which is blocked by the front seat. As a result, the air temperature of the rear passenger’s feet and abdomen can hardly drop down, which are higher than the air temperature of his head.
Influence of driving state on air temperature
Figures 7 and 8, respectively, show the air temperature distribution of different parts of the driver and the rear passenger under driving conditions. It can be seen that the air temperature stratification law in the cabin is basically consistent with that under idle speed. However, through comparative analysis with idle speed condition, when the temperature change in cabin tends to be flat, the maximum air temperature difference between layers under driving condition is smaller than that under idle condition. As shown in Table 2, the air temperature is more uniform than that under idle condition. It is the reason that in the process of driving, the convective heat transfer between air and cabin is enhanced. The heat transfer from inside to outside under driving condition is more, when the inside temperature is higher than the outside one, so the cooling rate of cabin is higher than that under idle condition. The heat transfer from outside to inside under driving condition is more, when the outside temperature is lower than the inside one, so the cooling rate of cabin is less than that under idle condition. Combined with the effects of the two processes, the maximum air temperature difference under driving condition is smaller than that under idle condition.

Vertical air temperature change of driver’s position (driving).

Vertical air temperature change of passenger’s position (driving).
Comparison of maximum temperature stratification.
Relative humidity distribution
Figures 9 and 10, respectively, show relative humidity changes of different positions in the cabin under idle and driving condition. In the preheating stage, the relative humidity in the cabin decreases gradually and tends to be stable in the later stage, which decreases from about 50% at the beginning to about 35% at the end. The relative humidity in the front of compartment is lower than that in the rear compartment. The main reason for the great change of humidity in the cabin is that the air temperature rises greatly, while the moisture content stays the same, which causes the relative humidity to decrease greatly.

Relative humidity change of different positions (idle).

Relative humidity change of different positions (driving).
In the cooling stage, the relative humidity in the cabin increases gradually and tends to be stable in the later stage, which increases from about 50% at the beginning to about 35% at the end. It is the reason for the rapid rise of relative humidity in the cabin that not only the air temperature in the cabin drops sharply but also when just enter the cabin, the volunteers sweat a lot due to the high temperature, making the air moisture content greatly increase.
Seat temperature distribution
The change diagram of seat temperature with time, which is mean value of idle and driving conditions, is shown in Figure 11. As can be seen from the figure, the temperature of the back, hips, and waist of the seat is basically the same. In the heating stage, the temperature of the seat is higher than that of the inside air. In the cooling stage, the air temperature drops rapidly, while the seat temperature drops slowly, falling to a temperature close to that of human skin after 80 min. In the heating stage, the seat temperature keeps rising because of the heat storage. When it reaches about 30 min, due to the heat convection transfer between the seat and the lower air, the heat dissipation loss of seat gradually equates with heat storage, and the temperature basically remains unchanged. In the cooling stage, since the seat back, waist, and hip are in contact with the clothing, the heat is dissipated slowly through conduction, resulting in a very slow cooling rate.

Seat temperature changes over time.
Analysis of thermal comfort
Analysis on TSV change
Figure 12 show that the change of subjective TSV of drivers and rear passengers with time under idle and driving conditions. When the participants entered the cabin from a comfortable office environment, although the air conditioner was turned on, the skin temperature rose rapidly in a short time because the air and seat temperature in the cabin was higher after the cabin was exposed to the sun. Therefore, the participants got hotter and hotter, and the TSV kept rising. The temperature in the cabin decreased continuously and then tends to stabilize, so did the TSV.

TSV changes with time.
Through comparative analysis, the TSV of the driver and the rear passenger are basically the same in the TSV rising stage. It is because, the air temperature in the cabin is much higher at this stage and the human thermal sensation is less affected by the cold air of air-conditioning. In the TSV descending stage, the TSV of the driver and the rear passenger decrease at different rates and the TSV of the rear passenger are always higher than that of the driver. It is because, the cold air directly blows to the driver, while the rear passenger is less affected by the cold air due to the obstruction of the seat. Therefore, the cooling effect of the front of the cabin is better than that of the rear.
The duration of TSV rising under driving condition is longer than that under idle condition. It is the reason that the occupants, especially the driver, are more nervous and concentrated in driving, and the metabolic rates of both are higher than that in idle. In the descending stage of TSV, the driver’s TSV are always lower than the rear passenger’s, which was similar to the idling condition. The difference is that the TSV of the driver and rear passenger under driving conditions are always higher than that under the idle condition. The reason is that the metabolic rates of both are higher than that under idle condition and the difference of metabolic rates in the descending stage has more significant influence on TSV than that in the rising stage.
Analysis on absolute deviation of Teq model
As shown in Figures 13 and 14, human thermal sensation of experiment value is compared with that of

Comparison between experimental value and

Comparison between experimental value and
Analysis on absolute deviation of PMV model
As shown in Figures 15 and 16, human thermal sensation of experiment value is compared with that of PMV model under idle and driving conditions. The time trend of the PMV model is basically the same as that of experiment value, but there are deviations between the data of both. Under idle condition, the maximum deviation can reach 0.6, while 92.3% of deviations are within the range of −0.5 to +0.2. It indicates that the prediction deviation of the PMV model under driving condition is larger while the prediction accuracy is lower, compared with that under idle condition.

Comparison between experimental value and PMV model (idle).

Comparison between experiment value and PMV model (driving).
It can be seen from the analysis above that the prediction accuracies of the two models under idle condition are higher than that under driving condition. In the TSV rising stage, the prediction accuracy of PMV is higher than that of the
Analysis on relative deviation of models
In order to facilitate the analysis of relative deviation, the thermal sensation index range of −3 to +3 is adjusted to 0 to +7. Figure 17 shows the relationship on TSV between experimental value and the predicted value of the

Relative deviation on TSV between experimental value and predicated value of

Relative deviation on TSV between experimental value and PMV model.
Previous studies have mostly used PMV model or
Conclusion
In this article, the thermal environment and the human thermal comfort under idling and driving conditions were experimentally studied. The cooling characteristics of the thermal environment of the driver’s and the rear passenger’s positions in the cabin under the cooling condition were analyzed, and then the driver’s thermal comfort was analyzed, based on the PMV model and
The air temperature in the cabin is mainly affected by the space position and driving state. Under idle condition, in the preheating stage, the air temperature of the driver and rear passenger from top to bottom are distributed from high to low, and the heating rate and maximum temperature difference of the rear passenger are slightly lower than that of the drive. Under idle condition, in the cooling stage, the air temperature distribution rule at the driver’s position is the same as in the preheating stage, but the temperature difference is smaller, while the air temperature at the feet and abdomen of the rear passenger are higher than that at the head. Under driving condition, the air temperature distribution rule in the cabin is basically the same as that under idle condition, but the temperature drops faster than that under idle condition and the temperature uniformity is higher than that under idle condition.
The humidity distribution of the cabin is also affected by the space position and driving state. Under idle condition, the relative humidity in the front compartment is lower than that in the rear compartment. In the cooling stage, the relative humidity increases with the progress of cooling and tends to be stable in the later stage. Under driving condition, the relative humidity difference between the front and rear is smaller than that under idle condition.
The temperature of the car seat is affected by solar radiation, which is slightly higher than that of the air in the cabin during the heating stage. During the cooling stage, the temperature of the seats drops slowly, falling to a temperature closed to that of human skin after 80 min.
The time variation trend of TSV of PMV model and
