TY - JOUR AU1 - Li, Weijian AU2 - Chen, Jiqing AU3 - Lan, Fengchong AB - Introduction Intelligent and automatic systems in vehicles have been receiving attention for recent years. In such an autonomous vehicle, the HVAC (Heating, Ventilation, and Air Conditioning) system is required to enable self-serving passengers’ thermal comfort with energy-saving strategies. The thermal state of the cabin is the result of its interaction with the external environment, such as solar radiation, air as well as dust [1–3], and fails to meet human comfort objectives [4, 5]. Without the cabin pre-conditioning, passengers entering the cabin have to endure long periods of thermal discomfort and intense self-regulation, as well as toxic volatiles [6, 7]. The HVAC is also bound to take time to satisfy comfort objectives with additional energy. How to achieve efficient and accurate HVAC pre-conditioning controlling according to the different initial thermal states in the cabin is a problem to be solved. Actually, city layout on landscapes planning makes local climates strongly different between road and outdoor semi-open spaces, where parked cars experience various thermal states in cabins, although it is the same entire climate at one area district. As known in statistics, the ratio of the parking time is about 10–20 to the driving during the automobile serving period [8], while more than one-third are open-air parked [9]. Although HVAC systems of the building already apply pre-conditioning, the vehicle cabin has a more changeable thermal environment and complex interactions than buildings [10]. In such a narrow car space, the different material properties and sizes of structural parts have various thermal responses to environmental conditions. An inappropriate value might be set without a correct perception of the cabin environment [11], resulting in dissatisfaction and energy waste [12, 13]. Therefore, cabin pre-conditioning strategies must consider the dynamically changing thermal states in cabins. Many studies have been carried out to improve the thermal environment in vehicles [14, 15]. Healthy ventilation and temperature control strategies are applied to reduce the threat of harmful gases and viruses in a vehicle cabin [15]. Supply control of variable air volumes helps to meet different thermal demand objectives in cabin [16]. Thermal condition strategies based on thermal loads and thermal comfort can effectively reduce the energy consumption of electric vehicles [17–19]. Related strategies focus on thermal comfort and energy management when the vehicle environment stabilizes, while effective management strategies for the initial thermal environment is lacking. For cabin pre-conditioning, serval studies applied additional structural to reduce the heat accumulated by parked vehicles [20–24]. Pan et al. [20] and Lee et al. [22] by adding a pre-ventilation system, were to reduce cabin thermal load with fan and solar sunroof. Maan et al. [24] proposed a device to cool the vehicle cabin with compressed air by vaporizing endothermic, but the safety of the air tank was another problem to be worth considering. Related studies about cabin pre-conditioning with fixed cooling parameters cannot adapt to the dynamically changing urban climate, which would provide an unsatisfactory outcome once the environmental conditions deviate from the scope of their setting conditions. As reviewed, a systematic understanding of urban climate on vehicle cabins and perfect cooling strategies are lacking. It is necessary to for pre-conditioning identify the cabin thermal environment characteristics formed by climate conditions. Evaluating passengers’ satisfaction with thermal comfort is another issue for the cabin pre-conditioning. In the narrow cabin, the thermal states and seat surface contact with the passenger both directly affect the passenger’s thermal sensation. The impact of tactility factors cannot be ignored to human comfort. An effective cabin pre-conditioning should meet the human thermal needs for these factors, especially the surface temperature, air quality, and energy consumption. The PMV- PPD (Predicted Mean Vote- Predicted Percentage Dissatisfied) is most commonly used for assessing passenger satisfaction with air thermal state [25, 26], but the evaluation methods for other factors are still needed. How to evaluate the passenger’s satisfaction with the cabin pre-conditioning comprehensively troubles the planning of the strategy. A suitable pre-conditioning decision method is one of the subjects for the intelligent vehicles adapting to changing external environment conditions. Reliable data on cabin thermal characteristics is essential to pre-conditioning decisions. The data in the study come from experiments conducted in a tropical desert climate which collect the thermal properties outside and inside the cabin. The relationship of thermal responses between the inside and outside cabin was analyzed by thermal transfer. According to cabin air temperature varied from hot to neutral, the cabin thermal environmental variable would be divided into multiple levels corresponding to the cooling amplitude in pre-conditioning. A data-driven decision model with Cubic Support Vector Machines is proposed to identify the cabin thermal environment level and determine pre-conditioning schemes. Considering the multiple satisfaction objectives, a comprehensive evaluation index CEI is suggested to evaluate the passengers’ satisfaction on both thermal comfort, seat temperature, air quality, and energy consumption. CEI effectively provides assessments for pre-conditioning parameter schemes. Materials and methods Cabin pre-conditioning strategy based on thermal environment levels Driverless technology and intelligent vehicles place new demands on the management of the thermal environment in the vehicle cabin. When passengers need to use parked vehicles, pre-conditioning of the HVAC is asked to eliminate the worse environment in the cabin formed during the parking process. Among these, the cabin pre-conditioning with HVAC means that HVAC cools or heats the cabin environment before the passenger enters with energy-saving, thermal comfort and efficiency, activating in advance via network reservation commands. Fig 1 shows an outline of the HVAC pre-conditioning implementation with environment data driven. The pre-conditioning system obtains environmental parameters and predicts component temperatures firstly. Then, the level of the thermal environment state is identified based on measured parameters, and the corresponding level of pre-conditioning solution in the database is called up to set the cooling parameters of the HVAC. The implementation of multi-satisfaction settings parameters for pre-conditioning in different scenarios is the main focus of this research. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. Implementation process of the vehicle HVAC pre-conditioning system. The system was designed for internet-enabled vehicles; therefore, passengers send network commands to activate parked vehicle systems for ventilation. The pre-conditioning system runs 10 minutes before passenger demand. The vehicle system specifies the optimal air conditioning parameters based on the current environmental conditions to the satisfy the passengers. https://doi.org/10.1371/journal.pone.0266672.g001 All differences in geographic, spatial, and temporal dimensions, as well as time and seasons variations, lead to differences in the impact of weather on the cabin environment. It is laborious and inefficient to develop an optimal solution of pre-conditioning for each condition needing multiple iterative operations. Therefore, the study constructs a database of thermal environmental variables and classifies them with the research methods in Fig 2. In this way, it is only necessary to develop control strategies for each category level and call them from the database, without repeating the analysis. This data-driven method provides effective solutions and achieves promising results through data training and classification with a simple model structure and fast calculation. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. Research method of automobile cabin pre-conditioning. The system levels and classifies thermal environment. The database develops solutions for each thermal environment level. Cabin pre-conditioning invokes solutions based on the identified thermal environment level. https://doi.org/10.1371/journal.pone.0266672.g002 The research ideas of cabin pre-conditioning parameters decision have two main parts: candidate schemes formulation and data-driven decision model construction. On the one hand, typical climate scenarios from the experimental environment data set are selected to level multiple climate scenes. Parameter combinations of pre-conditioning are designed for corresponded levels with a proposed comprehensive evaluation index. On the other hand, the study suggests a classification algorithm model for identifying the climate level for the inputted environmental conditions and matching a suitable pre-conditioning parameters scheme. The system obtains the thermal environment states in an automotive cabin coupling experimental climate data and machine learning algorithms. The data-driven decision model matches satisfactory cooling solutions according to hot severity levels of the thermal environments with variables relationship in Fig 3. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. The data flow within the data-driven decision model. A classification model identifies and matches schemes based on the cabin temperature and humidity, seat and wall temperature and climate conditions. https://doi.org/10.1371/journal.pone.0266672.g003 Comprehensive evaluation index proposed for human multi-satisfaction The effect factors on the passenger satisfaction of cabin cooling is varied, especially the high-temperature inner walls and poor air quality inside. A suitable pre-conditioning scheme should meet passenger’s multiple satisfaction needs. However, the PMV model focuses on the thermal state of the air. Therefore, a comprehensive evaluation index (CEI), new subjective evaluation method, is proposed to characterize the passengers’ satisfaction on considering multi-satisfaction objectives. The proposed CEI combines the thermal comfort, surface-cooling state, air quality, and energy consumption as represented by Eq (1), (1) where CTC is the thermal comfort scores; Cs is the surface acceptance scores; Ca is the air quality scores; Cp is the energy consumption scores; kTC, ks, ka, kp are the weights coefficient of different items respectively. In this study, kTC = 0.35, ks = 0.2, ka = 0.15, kp = 0.3 that thermal comfort improving is prioritize and energy-saving is second. Compared with air quality, touching the seat and steering wheel in hot is more sensitive to human skin. (1) Thermal comfort The thermal comfort scores depend on the passenger’s thermal sensation and thermal comfort. These can be described by PMV-PPD indexes in EN ISO 14505 criterion, which are well-known variables. The global thermal sensation of people is divided into 7 levels, from -3 (cold) to 0 (neutral) to 3 (hot). These parameters are calculated as shown by Eqs (2) and (3) [26], (2) (3) where M is the metabolism (W/m2); W is the human body work (W/m2), considering that the passenger walks into the car with 116 W/m2; Pca is the cabin air pressure (Pa); Pa is the partial pressure of water vapor around the body, (kPa); fcl is the clothing area factor; Tcl is the clothes temperature (°C); Tr is the radiation temperature (°C); hc is the convection coefficient (Wm−2K−1). Relevant researches concluded that the human is comfortable and satisfied with the thermal environment in which PPD < 10 [27]. With the increase of PPD, thermal conditions can lead to increased discomfort for occupants. Therefore, the inverse ratio function is used to determine the relationship between PPD and CTC, and the value of CTC is normalized from 0 to 1. The thermal comfort scores CTC is proposed as follows [27], (4) (2) Surface-cooling state The surface temperature of the seat in narrow cabins is a local thermal condition that profoundly affects passenger. The hot surface would cause discomfort to the tactile sensation of passengers, and even skin-damaging during prolonged exposure. Medical clinical studies found that 42°C is the critical temperature to cause moderate temperature burns and over 70°C will cause skin high-temperature burns [28–30]. Therefore, during the cabin pre-conditioning, the temperature of the human main contact surface should decrease lower than 42°C to avoid skin moderate temperature burns. As the seating surface is the main contact area for passengers which absorbs a lot of heat, the satisfaction on seat state is necessary to be assessed. The model set 40°C as the critical temperature of satisfactory surface cooling down considering medical research [31] and the satisfaction increases as the surface cooling, so the Cs is proposed as, (5) (3) Air quality The releasing of harmful substances from car interiors is accelerated by high temperature conditions, harming the occupants’ physical health [32]. The external circulation is the main way of cabin ventilation during cooling, impacting the air quality significantly [33]. The air circulation rate β is a parameter to control the mixing ratio of ambient air and cabin circulating air. A low circulation rate β can obtain a high percentage of outside fresh air. Therefore, the air quality scores Ca in CEI is proposed as, (6) (4) Energy consumption Actuality, the energy consumption of HVAC is calculated by the performance compression curve of the compressor and the change in enthalpy actuality. To simplify the calculation, the study approximates that energy consumption is inversely proportional to the HVAC outlet airflow temperature Tac [18]. High energy consumption would reduce users’ satisfaction. With normalization, the energy consumption satisfaction item is proposed as, (7) where compressor power Pc is 14 kW; normalized coefficient np is 2.4. Model for cabin thermal environment Generally, the airflow and other local surfaces within the cabin have different thermal responses to various levels of environmental loads. Therefore, it is a premise of cabin pre-conditioning strategy formulation to analyze the impact of environmental factors and supply air characteristics on the cabin thermal characteristics. A model for cabin thermal environment is established to determine the thermal environment characteristic and cabin pre-conditioning effects of the considered schemes. Combining with the proposed CEI, suitable solutions for various climate levels are explored. (1) Thermal transfer between vehicle cabin and its surroundings The external environment affects the cabin’s thermal environment through convection and solar radiation. And the air supplied from the HVAC system makes the dynamic non-uniform change of the cabin thermal environment more intense. To study the cooling effects of different air conditions, the model of the thermal transfer in the cabin is built in Fig 4. According to the heat transfer relationship of the cabin, a thermal model is constructed to analyze the cabin thermal environment. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. The physical relationship of cabin thermal transfer. Solar radiation (ΔQsolar −w, ΔQsolar −s) is an important factor to the cabin temperature, the convection of external airflow (ΔQatm−s) as well. As the uneven airflow from HVAC, the cabin is divided into two thermal zones: the front row’s region (zone 1) and the rear row’s region (zone 2). During cooling, thermal zones are affected by energy and mass exchanges: 1) the convection (ΔQa−w, ΔQa−s) and cooling mass flow (Δmac) from air-conditioning; 2) the solar radiation load (ΔQsolar −w, ΔQsolar −s) received; 3) air circulation between zones (Δma1−a2) and in HVAC (Δmcir); 4) thermal exchanges with outside (ΔQatm−s). https://doi.org/10.1371/journal.pone.0266672.g004 (2) Energy balance in the cabin The cabin thermal model is performed with Lumped Capacitance Method. The temperature distribution in each thermal zone is assumed to be uniform. Dynamic equations are constructed to calculate the air temperature and humidity in the cabin, including the mass balance equation and energy balance equation. The mathematical model is derived from the energy conservation and mass balance in thermal zone 1 and 2. The heat exchange of air, wall, and seat are both considered. The equation can be given as follow, while the specific heat capacity ci of the air is set as a constant. The distribution ratio δ is the proportion of air from the vents in the dashboard (zone 1) that directly enters the rear area (zone 2). (8) (9) (10) (11) Qw1−w2 is thermal conduction of the body wall between zone 1 and zone 2. The thermal convective load Qca affects the air temperature of the cabin. The seats and walls receive the gain Qsolar −w, Qsolar −s of solar radiation while radiating their own heat Qw−s, Qs−w to the surroundings and cooling with convection Qca−w, Qca−s. The calculation method of the above thermal load can be obtained in the literature [17, 21]. (3) Moisture mass balance in the cabin The moisture mass in the cabin is affected by the water vapor content from HVAC. Regarding the cabin as the control volume, the moisture balance can be expressed as, (12) (13) The fresh air and circulating air are mixed to enter the evaporator with an enthalpy imix. The energy conservation and mass balance occur in the evaporator, and supplied air enthalpy iac can be expressed [17], (14) (15) where mwa is the mass of water condensed on the cooling surface of the evaporator. Classification model for environment levels identification Identifying the current thermal environment level of the cabin is a prerequisite for targeted and satisfactory pre-conditioning work. The research classifies the real-time weather conditions and the cabin temperature with classification algorithms to clarify which level of the cabin environment is in. The classification algorithm determines a reasonable parameter scheme based on the cabin climate level. The algorithms obtain objective classification rules through the analysis of the training set of known categories, mapping each attribute set x to a pre-defined class label y, to predict the category of the new data. The study surveys Support Vector Machines (SVM) [34], K-nearest neighbor (KNN) [35], and Decision Tree algorithms [34] for thermal environment levels identification when cabin pre-conditioning on. The selected methods are well established and their basic principles can be obtained from the relevant literature. By comparing the effectiveness of the methods, the method with the highest accuracy in identifying the thermal environment class is selected. Cabin pre-conditioning strategy based on thermal environment levels Driverless technology and intelligent vehicles place new demands on the management of the thermal environment in the vehicle cabin. When passengers need to use parked vehicles, pre-conditioning of the HVAC is asked to eliminate the worse environment in the cabin formed during the parking process. Among these, the cabin pre-conditioning with HVAC means that HVAC cools or heats the cabin environment before the passenger enters with energy-saving, thermal comfort and efficiency, activating in advance via network reservation commands. Fig 1 shows an outline of the HVAC pre-conditioning implementation with environment data driven. The pre-conditioning system obtains environmental parameters and predicts component temperatures firstly. Then, the level of the thermal environment state is identified based on measured parameters, and the corresponding level of pre-conditioning solution in the database is called up to set the cooling parameters of the HVAC. The implementation of multi-satisfaction settings parameters for pre-conditioning in different scenarios is the main focus of this research. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. Implementation process of the vehicle HVAC pre-conditioning system. The system was designed for internet-enabled vehicles; therefore, passengers send network commands to activate parked vehicle systems for ventilation. The pre-conditioning system runs 10 minutes before passenger demand. The vehicle system specifies the optimal air conditioning parameters based on the current environmental conditions to the satisfy the passengers. https://doi.org/10.1371/journal.pone.0266672.g001 All differences in geographic, spatial, and temporal dimensions, as well as time and seasons variations, lead to differences in the impact of weather on the cabin environment. It is laborious and inefficient to develop an optimal solution of pre-conditioning for each condition needing multiple iterative operations. Therefore, the study constructs a database of thermal environmental variables and classifies them with the research methods in Fig 2. In this way, it is only necessary to develop control strategies for each category level and call them from the database, without repeating the analysis. This data-driven method provides effective solutions and achieves promising results through data training and classification with a simple model structure and fast calculation. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. Research method of automobile cabin pre-conditioning. The system levels and classifies thermal environment. The database develops solutions for each thermal environment level. Cabin pre-conditioning invokes solutions based on the identified thermal environment level. https://doi.org/10.1371/journal.pone.0266672.g002 The research ideas of cabin pre-conditioning parameters decision have two main parts: candidate schemes formulation and data-driven decision model construction. On the one hand, typical climate scenarios from the experimental environment data set are selected to level multiple climate scenes. Parameter combinations of pre-conditioning are designed for corresponded levels with a proposed comprehensive evaluation index. On the other hand, the study suggests a classification algorithm model for identifying the climate level for the inputted environmental conditions and matching a suitable pre-conditioning parameters scheme. The system obtains the thermal environment states in an automotive cabin coupling experimental climate data and machine learning algorithms. The data-driven decision model matches satisfactory cooling solutions according to hot severity levels of the thermal environments with variables relationship in Fig 3. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. The data flow within the data-driven decision model. A classification model identifies and matches schemes based on the cabin temperature and humidity, seat and wall temperature and climate conditions. https://doi.org/10.1371/journal.pone.0266672.g003 Comprehensive evaluation index proposed for human multi-satisfaction The effect factors on the passenger satisfaction of cabin cooling is varied, especially the high-temperature inner walls and poor air quality inside. A suitable pre-conditioning scheme should meet passenger’s multiple satisfaction needs. However, the PMV model focuses on the thermal state of the air. Therefore, a comprehensive evaluation index (CEI), new subjective evaluation method, is proposed to characterize the passengers’ satisfaction on considering multi-satisfaction objectives. The proposed CEI combines the thermal comfort, surface-cooling state, air quality, and energy consumption as represented by Eq (1), (1) where CTC is the thermal comfort scores; Cs is the surface acceptance scores; Ca is the air quality scores; Cp is the energy consumption scores; kTC, ks, ka, kp are the weights coefficient of different items respectively. In this study, kTC = 0.35, ks = 0.2, ka = 0.15, kp = 0.3 that thermal comfort improving is prioritize and energy-saving is second. Compared with air quality, touching the seat and steering wheel in hot is more sensitive to human skin. (1) Thermal comfort The thermal comfort scores depend on the passenger’s thermal sensation and thermal comfort. These can be described by PMV-PPD indexes in EN ISO 14505 criterion, which are well-known variables. The global thermal sensation of people is divided into 7 levels, from -3 (cold) to 0 (neutral) to 3 (hot). These parameters are calculated as shown by Eqs (2) and (3) [26], (2) (3) where M is the metabolism (W/m2); W is the human body work (W/m2), considering that the passenger walks into the car with 116 W/m2; Pca is the cabin air pressure (Pa); Pa is the partial pressure of water vapor around the body, (kPa); fcl is the clothing area factor; Tcl is the clothes temperature (°C); Tr is the radiation temperature (°C); hc is the convection coefficient (Wm−2K−1). Relevant researches concluded that the human is comfortable and satisfied with the thermal environment in which PPD < 10 [27]. With the increase of PPD, thermal conditions can lead to increased discomfort for occupants. Therefore, the inverse ratio function is used to determine the relationship between PPD and CTC, and the value of CTC is normalized from 0 to 1. The thermal comfort scores CTC is proposed as follows [27], (4) (2) Surface-cooling state The surface temperature of the seat in narrow cabins is a local thermal condition that profoundly affects passenger. The hot surface would cause discomfort to the tactile sensation of passengers, and even skin-damaging during prolonged exposure. Medical clinical studies found that 42°C is the critical temperature to cause moderate temperature burns and over 70°C will cause skin high-temperature burns [28–30]. Therefore, during the cabin pre-conditioning, the temperature of the human main contact surface should decrease lower than 42°C to avoid skin moderate temperature burns. As the seating surface is the main contact area for passengers which absorbs a lot of heat, the satisfaction on seat state is necessary to be assessed. The model set 40°C as the critical temperature of satisfactory surface cooling down considering medical research [31] and the satisfaction increases as the surface cooling, so the Cs is proposed as, (5) (3) Air quality The releasing of harmful substances from car interiors is accelerated by high temperature conditions, harming the occupants’ physical health [32]. The external circulation is the main way of cabin ventilation during cooling, impacting the air quality significantly [33]. The air circulation rate β is a parameter to control the mixing ratio of ambient air and cabin circulating air. A low circulation rate β can obtain a high percentage of outside fresh air. Therefore, the air quality scores Ca in CEI is proposed as, (6) (4) Energy consumption Actuality, the energy consumption of HVAC is calculated by the performance compression curve of the compressor and the change in enthalpy actuality. To simplify the calculation, the study approximates that energy consumption is inversely proportional to the HVAC outlet airflow temperature Tac [18]. High energy consumption would reduce users’ satisfaction. With normalization, the energy consumption satisfaction item is proposed as, (7) where compressor power Pc is 14 kW; normalized coefficient np is 2.4. Model for cabin thermal environment Generally, the airflow and other local surfaces within the cabin have different thermal responses to various levels of environmental loads. Therefore, it is a premise of cabin pre-conditioning strategy formulation to analyze the impact of environmental factors and supply air characteristics on the cabin thermal characteristics. A model for cabin thermal environment is established to determine the thermal environment characteristic and cabin pre-conditioning effects of the considered schemes. Combining with the proposed CEI, suitable solutions for various climate levels are explored. (1) Thermal transfer between vehicle cabin and its surroundings The external environment affects the cabin’s thermal environment through convection and solar radiation. And the air supplied from the HVAC system makes the dynamic non-uniform change of the cabin thermal environment more intense. To study the cooling effects of different air conditions, the model of the thermal transfer in the cabin is built in Fig 4. According to the heat transfer relationship of the cabin, a thermal model is constructed to analyze the cabin thermal environment. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. The physical relationship of cabin thermal transfer. Solar radiation (ΔQsolar −w, ΔQsolar −s) is an important factor to the cabin temperature, the convection of external airflow (ΔQatm−s) as well. As the uneven airflow from HVAC, the cabin is divided into two thermal zones: the front row’s region (zone 1) and the rear row’s region (zone 2). During cooling, thermal zones are affected by energy and mass exchanges: 1) the convection (ΔQa−w, ΔQa−s) and cooling mass flow (Δmac) from air-conditioning; 2) the solar radiation load (ΔQsolar −w, ΔQsolar −s) received; 3) air circulation between zones (Δma1−a2) and in HVAC (Δmcir); 4) thermal exchanges with outside (ΔQatm−s). https://doi.org/10.1371/journal.pone.0266672.g004 (2) Energy balance in the cabin The cabin thermal model is performed with Lumped Capacitance Method. The temperature distribution in each thermal zone is assumed to be uniform. Dynamic equations are constructed to calculate the air temperature and humidity in the cabin, including the mass balance equation and energy balance equation. The mathematical model is derived from the energy conservation and mass balance in thermal zone 1 and 2. The heat exchange of air, wall, and seat are both considered. The equation can be given as follow, while the specific heat capacity ci of the air is set as a constant. The distribution ratio δ is the proportion of air from the vents in the dashboard (zone 1) that directly enters the rear area (zone 2). (8) (9) (10) (11) Qw1−w2 is thermal conduction of the body wall between zone 1 and zone 2. The thermal convective load Qca affects the air temperature of the cabin. The seats and walls receive the gain Qsolar −w, Qsolar −s of solar radiation while radiating their own heat Qw−s, Qs−w to the surroundings and cooling with convection Qca−w, Qca−s. The calculation method of the above thermal load can be obtained in the literature [17, 21]. (3) Moisture mass balance in the cabin The moisture mass in the cabin is affected by the water vapor content from HVAC. Regarding the cabin as the control volume, the moisture balance can be expressed as, (12) (13) The fresh air and circulating air are mixed to enter the evaporator with an enthalpy imix. The energy conservation and mass balance occur in the evaporator, and supplied air enthalpy iac can be expressed [17], (14) (15) where mwa is the mass of water condensed on the cooling surface of the evaporator. Classification model for environment levels identification Identifying the current thermal environment level of the cabin is a prerequisite for targeted and satisfactory pre-conditioning work. The research classifies the real-time weather conditions and the cabin temperature with classification algorithms to clarify which level of the cabin environment is in. The classification algorithm determines a reasonable parameter scheme based on the cabin climate level. The algorithms obtain objective classification rules through the analysis of the training set of known categories, mapping each attribute set x to a pre-defined class label y, to predict the category of the new data. The study surveys Support Vector Machines (SVM) [34], K-nearest neighbor (KNN) [35], and Decision Tree algorithms [34] for thermal environment levels identification when cabin pre-conditioning on. The selected methods are well established and their basic principles can be obtained from the relevant literature. By comparing the effectiveness of the methods, the method with the highest accuracy in identifying the thermal environment class is selected. Thermal environment characteristic experiments for parked vehicles The pre-conditioning strategy uses a hierarchical approach to rate the thermal states in the cabin, so it is required to clarify how the thermal states are caused by climate conditions. The climate conditions, such as solar radiation, ambient air temperature and humidity, vary when vehicles are parked in open-space, shade and underground. Cars in open parking have the most comprehensive range of environmental parameters experienced by vehicles. Therefore, the study experimentally investigated the cabin thermal environment characteristics formed by climate conditions for an open parked vehicle. The experiment provides support for rating the thermal environment levels and model verification. Experimental design and data collection The tested vehicle is a sport utility vehicle (SUV) with five doors, black. The front end of the vehicle faces the south without any shelter to ensure the most severe thermal environment inside the car. The natural climate exposure experiment was carried out in Jeddah, Saudi Arabia (21°3’ N, 39°1’ E), in a tropical desert climate. The duration of the experiment was for a whole year, from March 2016 to February 2017. The experiment examined air temperature, humidity, precipitation, total solar radiation of ambient environment and cabin environment, the thermal characteristic of vehicle structural as well. Ambient environmental parameters were collected via a meteorological test platform. The cabin temperature data of wall surfaces were logged using T-type thermocouples (± 0.1°C) with the Agilent 34970A data collector. The temperature measurement points on the car body were set as shown in Fig 5. The case study measured the variables of solar radiation, atmosphere temperature, and relative humidity. The sensors were installed as shown in Fig 6. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. Arrangement of thermocouples in the cabin. The temperature measurement points on the car body were set to analyze the thermal load of the vehicle. (a): Measurement points on the car body. (b): Measurement points in the cabin. https://doi.org/10.1371/journal.pone.0266672.g005 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. Experimental scenario and related layout. A solar radiation meter was set to obtain the amount of solar radiation entering the cabin through the windows. Air temperature and humidity were measured by a sensor suspended from the roof. (a): Dashboard’s thermocouples and radiometer. (b): Air temperature and humidity sensor. https://doi.org/10.1371/journal.pone.0266672.g006 Experimental results The measured parameters were recorded every 5 minutes, and a total of 104,955 sets of valid experimental data were obtained from the experiment. Experiments were conducted to obtain the characteristics of the thermal load on the cabin during parking in a tropical desert climate. It has a clear picture of the maximum and minimum temperature states that exist in the vehicle, helping to define and differentiate levels of thermal environments. Under the tropical desert climate, the vehicle cabin was at a high temperature most of the time. The statistics of the air temperature in the cabin are shown in Fig 7. From Fig 7(a), the range from 27°C to 28°C has the most samples, accounting for 8.3%. The proportion of air temperature over 26°C accounts for 73%, which means that the thermal environment of the cabin is unsatisfactory most time. From June to September, cabin air temperature both exceeded 25°C, and even in winter exceeds 30°C under sun exposure in Fig 7(b), confirming the necessity for cabin pre-conditioning. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. Cabin air temperature Tca obtained by experiment. The frequency in (a) is described as ; while average cabin temperature f(tm) = 48.08 sin(0.172tm + 0.257)+ 11.06 sin(0.86tm + 1.65). tm is the month. (a): Frequency of occurrence of cabin air temperature. (b): Change in cabin air temperature by month. https://doi.org/10.1371/journal.pone.0266672.g007 For a vehicle parked in the open air, its thermal environment can be divided into three phases according to the changes in solar radiation: increasing radiation, decreasing radiation and none radiation. Fig 8(a) presents thermal characteristics of the cabin environment as solar radiation increases to decrease. It is worth noting that the highest cabin air temperature Tca occurs at the beginning of the solar radiation waning, while the lowest is at the appearance of solar radiation, due to the thermal storage effect of the material. During decreasing radiation, the greenhouse effect maintains the occupant compartment at a higher temperature than the same level of solar radiation during radiation increase. And during the none radiation phase, the occupant compartment temperature continues to decrease to a minimum, close to the ambient temperature with thermal equilibrium, shown in Fig 8(b). The RH max occurs at the lowest temperature and the RH mini occurs at the highest temperature. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 8. Thermal load of main components during parking. During the radiation increasing, the relationship between the solar radiation and the components temperature follows as: atmosphere, Tatm = 0.006086I1.117 + 25.95; cabin air, Tca = 0.006086I1.117 + 25.95; the steering wheel, TSw = 0.1161I0.9124 + 20.69; front seat, Ts1 = 2.475I4 + 0.3333I3−9.094I2 + 10.64I + 50.25; rear seat, Ts2 = −0.1648I2 + 10.35I + 40.8; During the radiation decreasing, the correlation follows as: atmosphere ; steering wheel TSw = 1.739I0.4881 + 33.89; front seat Ts1 = −1.886I4 + 0.8018I3 + 4.752I2 + 6.391I + 53.98; rear seat Ts2 = −0.8943I4 + 1.891I3 + 1.865I2 + 1.704I + 51.24. In (b), the atmospheric temperature is described as and ; cabin air temperature are and Tacmin = ; atmospheric RH are described as RHatm max = and ; cabin RH are described as and 6.219tm + 28.13. (a): The influence of changes in solar radiation. (b): Maximum and minimum of air temperature and relative humidity. https://doi.org/10.1371/journal.pone.0266672.g008 The highest temperature measured in the crew cabin during the experiment was on 15 July, as shown in Fig 9. The experimental results in Fig 9 show that the air temperature in the cabin increases with the atmosphere due to the solar radiation. Due to the thermal insulation effect of the car interior, the cabin air temperature remains at a relatively high level in the early part of the evening. The state of time is therefore also an influencing factor for the cabin environment to be considered. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 9. Temperature and humidity distribution on July 15. The car suffered from the worst heat load at noon, with a maximum Tca of 68.66°C. https://doi.org/10.1371/journal.pone.0266672.g009 Experimental design and data collection The tested vehicle is a sport utility vehicle (SUV) with five doors, black. The front end of the vehicle faces the south without any shelter to ensure the most severe thermal environment inside the car. The natural climate exposure experiment was carried out in Jeddah, Saudi Arabia (21°3’ N, 39°1’ E), in a tropical desert climate. The duration of the experiment was for a whole year, from March 2016 to February 2017. The experiment examined air temperature, humidity, precipitation, total solar radiation of ambient environment and cabin environment, the thermal characteristic of vehicle structural as well. Ambient environmental parameters were collected via a meteorological test platform. The cabin temperature data of wall surfaces were logged using T-type thermocouples (± 0.1°C) with the Agilent 34970A data collector. The temperature measurement points on the car body were set as shown in Fig 5. The case study measured the variables of solar radiation, atmosphere temperature, and relative humidity. The sensors were installed as shown in Fig 6. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. Arrangement of thermocouples in the cabin. The temperature measurement points on the car body were set to analyze the thermal load of the vehicle. (a): Measurement points on the car body. (b): Measurement points in the cabin. https://doi.org/10.1371/journal.pone.0266672.g005 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. Experimental scenario and related layout. A solar radiation meter was set to obtain the amount of solar radiation entering the cabin through the windows. Air temperature and humidity were measured by a sensor suspended from the roof. (a): Dashboard’s thermocouples and radiometer. (b): Air temperature and humidity sensor. https://doi.org/10.1371/journal.pone.0266672.g006 Experimental results The measured parameters were recorded every 5 minutes, and a total of 104,955 sets of valid experimental data were obtained from the experiment. Experiments were conducted to obtain the characteristics of the thermal load on the cabin during parking in a tropical desert climate. It has a clear picture of the maximum and minimum temperature states that exist in the vehicle, helping to define and differentiate levels of thermal environments. Under the tropical desert climate, the vehicle cabin was at a high temperature most of the time. The statistics of the air temperature in the cabin are shown in Fig 7. From Fig 7(a), the range from 27°C to 28°C has the most samples, accounting for 8.3%. The proportion of air temperature over 26°C accounts for 73%, which means that the thermal environment of the cabin is unsatisfactory most time. From June to September, cabin air temperature both exceeded 25°C, and even in winter exceeds 30°C under sun exposure in Fig 7(b), confirming the necessity for cabin pre-conditioning. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. Cabin air temperature Tca obtained by experiment. The frequency in (a) is described as ; while average cabin temperature f(tm) = 48.08 sin(0.172tm + 0.257)+ 11.06 sin(0.86tm + 1.65). tm is the month. (a): Frequency of occurrence of cabin air temperature. (b): Change in cabin air temperature by month. https://doi.org/10.1371/journal.pone.0266672.g007 For a vehicle parked in the open air, its thermal environment can be divided into three phases according to the changes in solar radiation: increasing radiation, decreasing radiation and none radiation. Fig 8(a) presents thermal characteristics of the cabin environment as solar radiation increases to decrease. It is worth noting that the highest cabin air temperature Tca occurs at the beginning of the solar radiation waning, while the lowest is at the appearance of solar radiation, due to the thermal storage effect of the material. During decreasing radiation, the greenhouse effect maintains the occupant compartment at a higher temperature than the same level of solar radiation during radiation increase. And during the none radiation phase, the occupant compartment temperature continues to decrease to a minimum, close to the ambient temperature with thermal equilibrium, shown in Fig 8(b). The RH max occurs at the lowest temperature and the RH mini occurs at the highest temperature. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 8. Thermal load of main components during parking. During the radiation increasing, the relationship between the solar radiation and the components temperature follows as: atmosphere, Tatm = 0.006086I1.117 + 25.95; cabin air, Tca = 0.006086I1.117 + 25.95; the steering wheel, TSw = 0.1161I0.9124 + 20.69; front seat, Ts1 = 2.475I4 + 0.3333I3−9.094I2 + 10.64I + 50.25; rear seat, Ts2 = −0.1648I2 + 10.35I + 40.8; During the radiation decreasing, the correlation follows as: atmosphere ; steering wheel TSw = 1.739I0.4881 + 33.89; front seat Ts1 = −1.886I4 + 0.8018I3 + 4.752I2 + 6.391I + 53.98; rear seat Ts2 = −0.8943I4 + 1.891I3 + 1.865I2 + 1.704I + 51.24. In (b), the atmospheric temperature is described as and ; cabin air temperature are and Tacmin = ; atmospheric RH are described as RHatm max = and ; cabin RH are described as and 6.219tm + 28.13. (a): The influence of changes in solar radiation. (b): Maximum and minimum of air temperature and relative humidity. https://doi.org/10.1371/journal.pone.0266672.g008 The highest temperature measured in the crew cabin during the experiment was on 15 July, as shown in Fig 9. The experimental results in Fig 9 show that the air temperature in the cabin increases with the atmosphere due to the solar radiation. Due to the thermal insulation effect of the car interior, the cabin air temperature remains at a relatively high level in the early part of the evening. The state of time is therefore also an influencing factor for the cabin environment to be considered. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 9. Temperature and humidity distribution on July 15. The car suffered from the worst heat load at noon, with a maximum Tca of 68.66°C. https://doi.org/10.1371/journal.pone.0266672.g009 Results and discussion Environmental levels division for pre-conditioning scenario Natural environmental conditions are variable and uncertain. This system levels multiple climate scenes corresponding to variables in the cabin air temperature and provides them with several tailored pre-conditioning solutions. Based on experimental results, the thermal environment levels in the vehicle cabin are graded by the magnitude of the cabin air temperature Tca to satisfy most of the real scenarios. The experimental data of Tca are ranked in descending order and classified into 6 levels of thermal environment varied from hot, warm to neutral levels according to the average intervals. Representative specify thermal environment scenarios are selected to character environment levels in Table 1 and Fig 10 considering the differences in solar radiation intensity I and ambient temperature Tatm over time. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 10. Characteristics of scenarios for different thermal environment levels. Tca and RHca are respectively the air temperature and relative humidity of the cabin, measured at the front and rear rows. Ts1 and Ts2 are the temperatures of the front and rear seat cushion, while Tw is derived from the average wall temperature in the thermal zone, including car door interior trim panel, pillars, the car roof, and floor. https://doi.org/10.1371/journal.pone.0266672.g010 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. The typical scenarios for different thermal environment levels. https://doi.org/10.1371/journal.pone.0266672.t001 Six levels indicate various pre-conditioning requirements of the cabin environment. The first two scenarios mainly occur in the afternoon with high solar radiation, while the latter two are mainly during the morning. The study sets up a simulation analysis to explore the suitable pre-conditioning schemes corresponding to these 6 thermal environmental levels. Pre-conditioning schemes for various thermal environment levels Requirements and objectives for cabin pre-conditioning. The developed cabin pre-conditioning schemes should achieve a satisfactory thermal environment within the specified time (10 minutes) using as low energy consumption as possible. External fresh air also should be considered to minimize the toxic gas content accumulated. This allows passengers to be as comfortable as possible in the cabin when they enter. Generally, it takes about 10 minutes for citizens to travel from their workplaces to their parking places [36]. Therefore, the pre-conditioning schemes should meet the two objectives. Firstly, the cabin pre-conditioning should be completed within 10 minutes; secondly, the CEI of the thermal environment after cabin pre-conditioning should be high as possible that satisfy passengers. Boundary conditions. The dynamic thermal responses during cooling are simulated with the established model of cabin thermal environment under various conditions and control strategies. In simulations, the settings of vehicle dimensions were taken from manufacturer data, 3100mm × 1700mm × 1300mm cabin sizes, 180ml HVAC compressor. The major material parameters including the thermodynamic and optical properties were consistent with previous related research [5, 37, 38]. The environmental parameters in Table 1 were used as boundary conditions for the different environmental levels. The external environmental factors were simplified as a constant value during the pre-conditioning. The cooling performance of the HVAC parameter combinations was analyzed sequentially for each thermal environment level. The simulations conditions respectively selected a value from the candidate ranges of TR, vac, and β shown in Table 2 to form a parameter combination [TR, vac, β] for cabin pre-conditioning. The refrigerant temperature TR, velocity of supplied air vac, and air circulation rate β are used as the control variables of HVAC to control the Tca and Δmca. The TR depends on the performance compression curve of the HVAC compressor, so TR needs an electronic control of the compressor to be managed. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. The typical scenarios for different thermal environment levels. https://doi.org/10.1371/journal.pone.0266672.t002 The velocity of supplied air selected in Table 2 corresponds to the blower power Pb being 100 W, 60 W, and 20 W. It is worth noting that when setting β=100%, β=0% is used first for 60 seconds, after which β=100% takes effect. This setup is to improve the air quality and reduce the toxic gas content in the cabin as possible. Parameter schemes for different thermal environment levels. To meet multi satisfaction objectives, the CEI proposed is applied to evaluate the satisfaction on various cabin pre-conditioning schemes. Through multiple simulations and comparisons with CEI, the research has formulated 6 pre-conditioning schemes for cabin thermal environment, as shown in Table 3. For Level 6, the airflow is supplied from the blower without cooling. Each combination scheme obtains the highest CEI index to match each scenario level. All the selected parameter schemes have excellent comprehensive performance and consider multi-satisfaction objectives that more realistically meet passengers’ needs for vehicle using. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. The typical scenarios for different thermal environment levels. https://doi.org/10.1371/journal.pone.0266672.t003 During the scheme analysis of Level 1, the cooling effect of different parameter combinations is shown in Fig 11. As the β decreases, the cooling effect would decrease leading to a higher air temperature. The main reason is that external circulation improves air quality by introducing ambient air but the ambient heat is also introduced into the system. The mass flow rate of the supply air also affects the cooling effect. With a small supplied air volume, the thermal environment in the cabin is still hot with Tca1 being 44.6°C and Ts1 being 59.7°C at 10 minutes. For Level 1 scenario, Tca1 is still over 24°C for any combinations. Also, Ts is still higher than 40°C, which will cause a slightly hot touch to the passengers. Limited by the performance of the HVAC system of a vehicle, a longer pre-conditioning time would be a more effective solution for more satisfactory cooling results for these extreme scenarios. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 11. Influence of different parameter combinations on scenario 1. Compared with others, the parameter scheme [6, 7.5, 100] wins the highest score, with a CEI being -0.099. https://doi.org/10.1371/journal.pone.0266672.g011 The parameters in Table 3 are conducted for satisfaction analysis. Level 1 has the most severe thermal environment experienced in the cabin. The CEI of all schemes are less than 0 for the Level 1 scenario, which means unsatisfied items still exist. Compared with other schemes, the parameter scheme [6, 7.5, 100] wins the highest score, with a CEI being -0.099 in Fig 11. Therefore, the scheme [TR=6°C, vac=7.5m/s, β=100%] is the better choice to solve those thermal environments in Level 1. The cooling effect in 10 minutes of different combinations for Level 2 is presented in Fig 12. With constant Tac, the reduction of β would fresh the cabin air but increase Tca. With multiple groups of simulations, the parameter combinations [6, 7.5, 100] and [8, 7.5, 100] enable to make cabin environment satisfactory. As the air circulation rate increases, the PPD index of the thermal environment decreases. In combinations [6, 7.5, 100] and [8, 7.5, 100], the PPD of the front thermal zone is less than 10%, and the back row is also close to 10%. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 12. Influence of different parameter combinations in 7.5 m/s on scenario 2. The parameter combination [8, 7.5, 100] wins the higher score. https://doi.org/10.1371/journal.pone.0266672.g012 According to the CEI model, the combination [8, 7.5, 100] wins the higher score, chosen as the pre-conditioning schemes to solve the thermal environment scenario in Level 2. Due to limited space, the simulation results of other scenarios are shown in Table A1 in S2 Appendix. Training and verification of decision model To accurately and effectively classify thermal environment levels and pre-conditioning decisions, data-driven methods are applied to identify inputted environmental conditions. The classification algorithms: support vector machines (SVM), Decision Tree and Random Forest (RF), and weighted K-Nearest Neighbor (KNN) are compared in environment level identification. The predicted environmental level would be used to match the corresponding parameter scheme to achieve pre-conditioning. To train the classification model, a training set is first generated from the experimental data in Section 3. To avoid the sample points selection aimlessly, the optimal Latin hypercube experimental design method was applied to obtain 300 samples as the training set, as shown in Fig 13. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 13. Tca and schemes distribution of the sample point. The 300 samples selected consider the uniformity and representativeness of the sample points in the thermal environment state. https://doi.org/10.1371/journal.pone.0266672.g013 The models were trained in toolboxes of MATLAB. To avoid over-fitting, k-fold Cross Validation was set to 10 [30]. The t, Tca, Tatm, Ts, Tw, RHa, RHatm, and PMV were used as the classification input labels, and the output was the serial number of the level. The accuracy is used to evaluate the performance of the classification models in Table 4. The confusion matrix for classification in Fig 14 shows that classification mistakes mainly happen between level 1 and 2 and between level 4 and 5. To test the classification models, 62 groups of test data are used and the test results are shown in Table 4. Compared with the analysis of the above classification models, the SVM using cubic kernel function has the highest accuracy in training, while it also has well accuracy in testing. Therefore, the study applies the trained Cubic SVM model to identify the cabin thermal environment level. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 14. Confusion matrix for classification models training. Classification mistakes mainly happen between scheme 1 and scheme 2 and between scheme 4 and scheme 5. (a): Cubic SVM. (b): RF. (c): Weighted KNN. https://doi.org/10.1371/journal.pone.0266672.g014 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. Accuracy of different classification models. https://doi.org/10.1371/journal.pone.0266672.t004 Synthesis of the above analysis process, a complete cabin pre-conditioning strategy decision model is established. Firstly, measure and predict the cabin environment variables with the weather data inputted; then input the predicted variables and weather data to Cubic SVM to predict its climate level category and match the required scheme. Case studies for cabin pre-conditioning Case studies are conducted to verify the effect of the pre-conditioning strategies matching by the Cubic SVM model. The case studies are chosen and analyzed in Table 5. The analysis results show that when the environmental characteristics are between the boundaries of the environment levels, matching scenarios and solutions are probably inconsistent with the actual optimal solution; otherwise, the model can effectively give the ideal solution. All the case studies keep thermal comfort index PMV within a suitable range, which is the most concerning element to passengers. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 5. Case study for pre-conditioning analysis. https://doi.org/10.1371/journal.pone.0266672.t005 Cases 1 to 6 apply the corresponding ideal parameters schemes for cabin pre-conditioning, with the cabin air temperature changes shown in Fig 15. With the specified scheme predicted, the pre-conditioning system adjusts the cabin air temperature to a suitable temperature range within 10 minutes. The pre-conditioning parameters developed for different thermal environment levels meet the thermal comfort needs of passengers entering the cabin. For Case 7 to Case 9, although the decision scheme of these cases is inconsistent with the ideal, the pre-conditioning effect is also comfortable and satisfactory under the solution provided by the model. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 15. The pre-conditioning results for case study. All case studies achieve a CEI rating of 0 or higher, bringing cabin heat down to an acceptable value. https://doi.org/10.1371/journal.pone.0266672.g015 According to the above analysis results, when the parked vehicle applies the cabin pre-conditioning scheme according to its scenario level, it can be expected to maximize the overall benefits to effectively improve the thermal comfort. Environmental levels division for pre-conditioning scenario Natural environmental conditions are variable and uncertain. This system levels multiple climate scenes corresponding to variables in the cabin air temperature and provides them with several tailored pre-conditioning solutions. Based on experimental results, the thermal environment levels in the vehicle cabin are graded by the magnitude of the cabin air temperature Tca to satisfy most of the real scenarios. The experimental data of Tca are ranked in descending order and classified into 6 levels of thermal environment varied from hot, warm to neutral levels according to the average intervals. Representative specify thermal environment scenarios are selected to character environment levels in Table 1 and Fig 10 considering the differences in solar radiation intensity I and ambient temperature Tatm over time. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 10. Characteristics of scenarios for different thermal environment levels. Tca and RHca are respectively the air temperature and relative humidity of the cabin, measured at the front and rear rows. Ts1 and Ts2 are the temperatures of the front and rear seat cushion, while Tw is derived from the average wall temperature in the thermal zone, including car door interior trim panel, pillars, the car roof, and floor. https://doi.org/10.1371/journal.pone.0266672.g010 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. The typical scenarios for different thermal environment levels. https://doi.org/10.1371/journal.pone.0266672.t001 Six levels indicate various pre-conditioning requirements of the cabin environment. The first two scenarios mainly occur in the afternoon with high solar radiation, while the latter two are mainly during the morning. The study sets up a simulation analysis to explore the suitable pre-conditioning schemes corresponding to these 6 thermal environmental levels. Pre-conditioning schemes for various thermal environment levels Requirements and objectives for cabin pre-conditioning. The developed cabin pre-conditioning schemes should achieve a satisfactory thermal environment within the specified time (10 minutes) using as low energy consumption as possible. External fresh air also should be considered to minimize the toxic gas content accumulated. This allows passengers to be as comfortable as possible in the cabin when they enter. Generally, it takes about 10 minutes for citizens to travel from their workplaces to their parking places [36]. Therefore, the pre-conditioning schemes should meet the two objectives. Firstly, the cabin pre-conditioning should be completed within 10 minutes; secondly, the CEI of the thermal environment after cabin pre-conditioning should be high as possible that satisfy passengers. Boundary conditions. The dynamic thermal responses during cooling are simulated with the established model of cabin thermal environment under various conditions and control strategies. In simulations, the settings of vehicle dimensions were taken from manufacturer data, 3100mm × 1700mm × 1300mm cabin sizes, 180ml HVAC compressor. The major material parameters including the thermodynamic and optical properties were consistent with previous related research [5, 37, 38]. The environmental parameters in Table 1 were used as boundary conditions for the different environmental levels. The external environmental factors were simplified as a constant value during the pre-conditioning. The cooling performance of the HVAC parameter combinations was analyzed sequentially for each thermal environment level. The simulations conditions respectively selected a value from the candidate ranges of TR, vac, and β shown in Table 2 to form a parameter combination [TR, vac, β] for cabin pre-conditioning. The refrigerant temperature TR, velocity of supplied air vac, and air circulation rate β are used as the control variables of HVAC to control the Tca and Δmca. The TR depends on the performance compression curve of the HVAC compressor, so TR needs an electronic control of the compressor to be managed. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. The typical scenarios for different thermal environment levels. https://doi.org/10.1371/journal.pone.0266672.t002 The velocity of supplied air selected in Table 2 corresponds to the blower power Pb being 100 W, 60 W, and 20 W. It is worth noting that when setting β=100%, β=0% is used first for 60 seconds, after which β=100% takes effect. This setup is to improve the air quality and reduce the toxic gas content in the cabin as possible. Parameter schemes for different thermal environment levels. To meet multi satisfaction objectives, the CEI proposed is applied to evaluate the satisfaction on various cabin pre-conditioning schemes. Through multiple simulations and comparisons with CEI, the research has formulated 6 pre-conditioning schemes for cabin thermal environment, as shown in Table 3. For Level 6, the airflow is supplied from the blower without cooling. Each combination scheme obtains the highest CEI index to match each scenario level. All the selected parameter schemes have excellent comprehensive performance and consider multi-satisfaction objectives that more realistically meet passengers’ needs for vehicle using. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. The typical scenarios for different thermal environment levels. https://doi.org/10.1371/journal.pone.0266672.t003 During the scheme analysis of Level 1, the cooling effect of different parameter combinations is shown in Fig 11. As the β decreases, the cooling effect would decrease leading to a higher air temperature. The main reason is that external circulation improves air quality by introducing ambient air but the ambient heat is also introduced into the system. The mass flow rate of the supply air also affects the cooling effect. With a small supplied air volume, the thermal environment in the cabin is still hot with Tca1 being 44.6°C and Ts1 being 59.7°C at 10 minutes. For Level 1 scenario, Tca1 is still over 24°C for any combinations. Also, Ts is still higher than 40°C, which will cause a slightly hot touch to the passengers. Limited by the performance of the HVAC system of a vehicle, a longer pre-conditioning time would be a more effective solution for more satisfactory cooling results for these extreme scenarios. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 11. Influence of different parameter combinations on scenario 1. Compared with others, the parameter scheme [6, 7.5, 100] wins the highest score, with a CEI being -0.099. https://doi.org/10.1371/journal.pone.0266672.g011 The parameters in Table 3 are conducted for satisfaction analysis. Level 1 has the most severe thermal environment experienced in the cabin. The CEI of all schemes are less than 0 for the Level 1 scenario, which means unsatisfied items still exist. Compared with other schemes, the parameter scheme [6, 7.5, 100] wins the highest score, with a CEI being -0.099 in Fig 11. Therefore, the scheme [TR=6°C, vac=7.5m/s, β=100%] is the better choice to solve those thermal environments in Level 1. The cooling effect in 10 minutes of different combinations for Level 2 is presented in Fig 12. With constant Tac, the reduction of β would fresh the cabin air but increase Tca. With multiple groups of simulations, the parameter combinations [6, 7.5, 100] and [8, 7.5, 100] enable to make cabin environment satisfactory. As the air circulation rate increases, the PPD index of the thermal environment decreases. In combinations [6, 7.5, 100] and [8, 7.5, 100], the PPD of the front thermal zone is less than 10%, and the back row is also close to 10%. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 12. Influence of different parameter combinations in 7.5 m/s on scenario 2. The parameter combination [8, 7.5, 100] wins the higher score. https://doi.org/10.1371/journal.pone.0266672.g012 According to the CEI model, the combination [8, 7.5, 100] wins the higher score, chosen as the pre-conditioning schemes to solve the thermal environment scenario in Level 2. Due to limited space, the simulation results of other scenarios are shown in Table A1 in S2 Appendix. Requirements and objectives for cabin pre-conditioning. The developed cabin pre-conditioning schemes should achieve a satisfactory thermal environment within the specified time (10 minutes) using as low energy consumption as possible. External fresh air also should be considered to minimize the toxic gas content accumulated. This allows passengers to be as comfortable as possible in the cabin when they enter. Generally, it takes about 10 minutes for citizens to travel from their workplaces to their parking places [36]. Therefore, the pre-conditioning schemes should meet the two objectives. Firstly, the cabin pre-conditioning should be completed within 10 minutes; secondly, the CEI of the thermal environment after cabin pre-conditioning should be high as possible that satisfy passengers. Boundary conditions. The dynamic thermal responses during cooling are simulated with the established model of cabin thermal environment under various conditions and control strategies. In simulations, the settings of vehicle dimensions were taken from manufacturer data, 3100mm × 1700mm × 1300mm cabin sizes, 180ml HVAC compressor. The major material parameters including the thermodynamic and optical properties were consistent with previous related research [5, 37, 38]. The environmental parameters in Table 1 were used as boundary conditions for the different environmental levels. The external environmental factors were simplified as a constant value during the pre-conditioning. The cooling performance of the HVAC parameter combinations was analyzed sequentially for each thermal environment level. The simulations conditions respectively selected a value from the candidate ranges of TR, vac, and β shown in Table 2 to form a parameter combination [TR, vac, β] for cabin pre-conditioning. The refrigerant temperature TR, velocity of supplied air vac, and air circulation rate β are used as the control variables of HVAC to control the Tca and Δmca. The TR depends on the performance compression curve of the HVAC compressor, so TR needs an electronic control of the compressor to be managed. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. The typical scenarios for different thermal environment levels. https://doi.org/10.1371/journal.pone.0266672.t002 The velocity of supplied air selected in Table 2 corresponds to the blower power Pb being 100 W, 60 W, and 20 W. It is worth noting that when setting β=100%, β=0% is used first for 60 seconds, after which β=100% takes effect. This setup is to improve the air quality and reduce the toxic gas content in the cabin as possible. Parameter schemes for different thermal environment levels. To meet multi satisfaction objectives, the CEI proposed is applied to evaluate the satisfaction on various cabin pre-conditioning schemes. Through multiple simulations and comparisons with CEI, the research has formulated 6 pre-conditioning schemes for cabin thermal environment, as shown in Table 3. For Level 6, the airflow is supplied from the blower without cooling. Each combination scheme obtains the highest CEI index to match each scenario level. All the selected parameter schemes have excellent comprehensive performance and consider multi-satisfaction objectives that more realistically meet passengers’ needs for vehicle using. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. The typical scenarios for different thermal environment levels. https://doi.org/10.1371/journal.pone.0266672.t003 During the scheme analysis of Level 1, the cooling effect of different parameter combinations is shown in Fig 11. As the β decreases, the cooling effect would decrease leading to a higher air temperature. The main reason is that external circulation improves air quality by introducing ambient air but the ambient heat is also introduced into the system. The mass flow rate of the supply air also affects the cooling effect. With a small supplied air volume, the thermal environment in the cabin is still hot with Tca1 being 44.6°C and Ts1 being 59.7°C at 10 minutes. For Level 1 scenario, Tca1 is still over 24°C for any combinations. Also, Ts is still higher than 40°C, which will cause a slightly hot touch to the passengers. Limited by the performance of the HVAC system of a vehicle, a longer pre-conditioning time would be a more effective solution for more satisfactory cooling results for these extreme scenarios. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 11. Influence of different parameter combinations on scenario 1. Compared with others, the parameter scheme [6, 7.5, 100] wins the highest score, with a CEI being -0.099. https://doi.org/10.1371/journal.pone.0266672.g011 The parameters in Table 3 are conducted for satisfaction analysis. Level 1 has the most severe thermal environment experienced in the cabin. The CEI of all schemes are less than 0 for the Level 1 scenario, which means unsatisfied items still exist. Compared with other schemes, the parameter scheme [6, 7.5, 100] wins the highest score, with a CEI being -0.099 in Fig 11. Therefore, the scheme [TR=6°C, vac=7.5m/s, β=100%] is the better choice to solve those thermal environments in Level 1. The cooling effect in 10 minutes of different combinations for Level 2 is presented in Fig 12. With constant Tac, the reduction of β would fresh the cabin air but increase Tca. With multiple groups of simulations, the parameter combinations [6, 7.5, 100] and [8, 7.5, 100] enable to make cabin environment satisfactory. As the air circulation rate increases, the PPD index of the thermal environment decreases. In combinations [6, 7.5, 100] and [8, 7.5, 100], the PPD of the front thermal zone is less than 10%, and the back row is also close to 10%. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 12. Influence of different parameter combinations in 7.5 m/s on scenario 2. The parameter combination [8, 7.5, 100] wins the higher score. https://doi.org/10.1371/journal.pone.0266672.g012 According to the CEI model, the combination [8, 7.5, 100] wins the higher score, chosen as the pre-conditioning schemes to solve the thermal environment scenario in Level 2. Due to limited space, the simulation results of other scenarios are shown in Table A1 in S2 Appendix. Training and verification of decision model To accurately and effectively classify thermal environment levels and pre-conditioning decisions, data-driven methods are applied to identify inputted environmental conditions. The classification algorithms: support vector machines (SVM), Decision Tree and Random Forest (RF), and weighted K-Nearest Neighbor (KNN) are compared in environment level identification. The predicted environmental level would be used to match the corresponding parameter scheme to achieve pre-conditioning. To train the classification model, a training set is first generated from the experimental data in Section 3. To avoid the sample points selection aimlessly, the optimal Latin hypercube experimental design method was applied to obtain 300 samples as the training set, as shown in Fig 13. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 13. Tca and schemes distribution of the sample point. The 300 samples selected consider the uniformity and representativeness of the sample points in the thermal environment state. https://doi.org/10.1371/journal.pone.0266672.g013 The models were trained in toolboxes of MATLAB. To avoid over-fitting, k-fold Cross Validation was set to 10 [30]. The t, Tca, Tatm, Ts, Tw, RHa, RHatm, and PMV were used as the classification input labels, and the output was the serial number of the level. The accuracy is used to evaluate the performance of the classification models in Table 4. The confusion matrix for classification in Fig 14 shows that classification mistakes mainly happen between level 1 and 2 and between level 4 and 5. To test the classification models, 62 groups of test data are used and the test results are shown in Table 4. Compared with the analysis of the above classification models, the SVM using cubic kernel function has the highest accuracy in training, while it also has well accuracy in testing. Therefore, the study applies the trained Cubic SVM model to identify the cabin thermal environment level. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 14. Confusion matrix for classification models training. Classification mistakes mainly happen between scheme 1 and scheme 2 and between scheme 4 and scheme 5. (a): Cubic SVM. (b): RF. (c): Weighted KNN. https://doi.org/10.1371/journal.pone.0266672.g014 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. Accuracy of different classification models. https://doi.org/10.1371/journal.pone.0266672.t004 Synthesis of the above analysis process, a complete cabin pre-conditioning strategy decision model is established. Firstly, measure and predict the cabin environment variables with the weather data inputted; then input the predicted variables and weather data to Cubic SVM to predict its climate level category and match the required scheme. Case studies for cabin pre-conditioning Case studies are conducted to verify the effect of the pre-conditioning strategies matching by the Cubic SVM model. The case studies are chosen and analyzed in Table 5. The analysis results show that when the environmental characteristics are between the boundaries of the environment levels, matching scenarios and solutions are probably inconsistent with the actual optimal solution; otherwise, the model can effectively give the ideal solution. All the case studies keep thermal comfort index PMV within a suitable range, which is the most concerning element to passengers. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 5. Case study for pre-conditioning analysis. https://doi.org/10.1371/journal.pone.0266672.t005 Cases 1 to 6 apply the corresponding ideal parameters schemes for cabin pre-conditioning, with the cabin air temperature changes shown in Fig 15. With the specified scheme predicted, the pre-conditioning system adjusts the cabin air temperature to a suitable temperature range within 10 minutes. The pre-conditioning parameters developed for different thermal environment levels meet the thermal comfort needs of passengers entering the cabin. For Case 7 to Case 9, although the decision scheme of these cases is inconsistent with the ideal, the pre-conditioning effect is also comfortable and satisfactory under the solution provided by the model. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 15. The pre-conditioning results for case study. All case studies achieve a CEI rating of 0 or higher, bringing cabin heat down to an acceptable value. https://doi.org/10.1371/journal.pone.0266672.g015 According to the above analysis results, when the parked vehicle applies the cabin pre-conditioning scheme according to its scenario level, it can be expected to maximize the overall benefits to effectively improve the thermal comfort. Conclusion The thermal states of the vehicle cabin fail to meet the human comfort objectives under the environmental conditions like solar radiation during parking. Based on passenger multi-target satisfaction needs, the study has developed a data-driven decision model to determine cabin pre-conditioning solutions efficiently and satisfactorily at different environmental conditions. In addition, the model quantifies passenger satisfaction considering comprehensive factors during cabin pre-conditioning. Notably, the following work is conducted: The research rates the thermal states levels for cabin pre-conditioning considering diversity of thermal environmental conditions. A data-driven decision model with the Cubic SVM algorithm was established to identify the thermal environment levels and match satisfactory pre-conditioning solutions. The characteristics of cabin thermal environment are analyzed from experiments in tropical desert climate. With the model for cabin thermal environment and reliable data on cabin thermal characteristics, the thermal responses on cabin air and structural under the effects of climate conditions and HVAC airflow is analyzed. Considering the multiple satisfaction objectives, a comprehensive evaluation index CEI model is proposed to evaluate passengers’ satisfaction in pre-conditioning on both thermal comfort, seat temperature, air quality, and energy consumption. With the assessments for passengers’ satisfaction in pre-conditioning, cabin pre-conditioning strategies are designed to different environment levels based on CEI. The performance of different temperatures, air volume and internal circulation ratios were compared for cabin pre-conditioning. The high cabin temperature increases the proportional demand for internal circulation, and the high ventilation speed shows high efficiency and energy efficiency. The decision algorithm of SVM, Decision Tree, and KNN are compared in environment levels identification. The Cubic SVM algorithm shows high effectiveness in identify thermal environment levels with 92.3% accuracy. Several scene cases are carried out and both thermal environments are improved within multi-satisfaction objectives. Supporting information S1 Appendix. Equations derivation process. https://doi.org/10.1371/journal.pone.0266672.s001 (PDF) S2 Appendix. Simulation results analysis. https://doi.org/10.1371/journal.pone.0266672.s002 (PDF) Acknowledgments We would like to thank collaborators from the State Key Laboratory of Environment Apparatus Research Institute for Conducting tests and collecting experimental data. TI - Environmental conditions driven method for automobile cabin pre-conditioning with multi-satisfaction objectives JF - PLoS ONE DO - 10.1371/journal.pone.0266672 DA - 2022-05-23 UR - https://www.deepdyve.com/lp/public-library-of-science-plos-journal/environmental-conditions-driven-method-for-automobile-cabin-pre-6bb5J4mNO2 SP - e0266672 VL - 17 IS - 5 DP - DeepDyve ER -