Modeling of Humidity in Passenger Cars Equipped with Mechanical Ventilation

: This paper presents research on humidity in a passenger car cabin with the use of supply ventilation without cooling the air. Based on the tests carried out and the humidity balance in the car, a model was developed for changing the humidity in the passenger car cabin as a function of time. The study of thermohumid conditions was carried out in two passenger cars. During the tests, the heating and cooling functions were turned o ﬀ . The relative humidity and temperature were measured outside the car before and after driving the car and in the supply air duct and inside the passenger car cabin while driving the car. The tests were carried out for an average temperature range from 20 to 42.9 ◦ C. In order to develop a model of humidity changes as a function of time, a humidity balance was prepared. Human-generated humidity in the car cabin depends mainly on the temperature inside the car and the age of the person and can range from 20 to 180 g / (h × person) for an adult in the temperature range of 20–43 ◦ C, while for a child under six years old the humidity ranges from 8 to 19.5 g / (h × person) in the temperature range 22–34 ◦ C. A formula of humidity generated by an adult and a child aged six years old was obtained as a function of temperature inside a passenger car. Based on the experimental research and the model developed, the humidity generated by a single adult and a six-year-old child in the car was determined. The developed model can be used in the automatic airﬂow adjustment systems in passenger cars.


Introduction
Air humidity is one of the parameters of air quality that significantly affects the comfort of driving a car and human health [1][2][3]. The main sources of moisture in the passenger car cabin are the driver and passengers inside the passenger car. The generation of moisture by humans depends primarily on human physical intensity and can range from about 30 g/h during sleep to over 300 g/h during intense physical effort [4]. Air quality in car cabins has been studied by many researchers [5,6], and a few of them have analyzed moisture. In Aquino et al. [5], an analysis of fungal contamination of air-conditioning filters collected from 21 cars was performed, and the study showed a population of fungi in all samples taken (100%), including toxin-forming fungi such as Penicillium, Fusarim, and Aspergillus, which indicates that air quality indoors may be a health hazard. These fungi cause severe allergic and pulmonary respiratory diseases. For health reasons and to ensure thermal comfort, the relative humidity of air at a temperature of 20-22 • C can be 30-65%. At a temperature of about 26 • C, the relative humidity should not exceed 55%. When the relative humidity exceeds 70%, there are conditions that cause the formation of fungi and mold that are harmful to people and can destroy the structure of the car. According to the results of a survey [7], 47.7% of the people from the surveyed population spent one to two hours in means of transport; in larger cities (over 500,000 inhabitants), this time was longer, and about 20% of the respondents from large cities spent 2-5 h a day in cars. Due to

Materials and Methods
The study of thermohumid conditions was carried out in two passenger cars-one a type of station wagon marked A and the other a type of hatchback marked B. The maximum number of people allowed in the tested cars together with the driver is five people. The car speed in all measurement series was maintained at a constant level of 70 km/h. Both cars are equipped with air heating and cooling functions.
A simplified diagram of the ventilation system is shown in Figure 1. The blown air is first taken in from the outside through the inlet (1) through the backward-curved centrifugal fan (2), then the air flows through the air filter (3) and heat exchangers (4) in which the air can be heated or cooled, and the air is then supplied through the channels (5-7) inside the car with the option of switching off individual channels. The cooling function is carried out by a thermodynamic cycle, while heating is most often done before the heat is collected from the engine. During the tests, the heating and cooling functions Energies 2020, 13, 2987 3 of 14 were turned off. The relative humidity of the supply air in such a system is equal to the outside air humidity. Despite the heating function being turned off, the supply air had a higher temperature than the outside air, which was caused by heat gains from the running engine and solar radiation.
Energies 2020, 13, x FOR PEER REVIEW 3 of 15 most often done before the heat is collected from the engine. During the tests, the heating and cooling functions were turned off. The relative humidity of the supply air in such a system is equal to the outside air humidity. Despite the heating function being turned off, the supply air had a higher temperature than the outside air, which was caused by heat gains from the running engine and solar radiation.
The main elements included in the passenger car ventilation system: 1-air intake, 2-supply fan, 3-air filter, 4-heat exchangers (cooler, heater), 5-air supply duct under the glass, 6-duct air supply towards the driver and passenger, 7-air supply channel towards the driver's and passengers' legs, 8-car windscreen.
The relative humidity and temperature in the tested passenger cars were recorded while driving the car with the windows closed. The research was conducted in Poland in the spring and the summer when mechanical ventilation in passenger cars is intensively exploited. Measurements were made outside of an urban area on a national road, where the car's velocity varied from 70 to 90 km/h. The relative humidity and temperature were measured outside the car before and after driving the car as well as in the supply air duct and inside the passenger car cabin while driving the car using a Testo-435 recorder, whose parameters are shown in [39]. During the measurement, the airflow was set only before the windscreen in channel 5; the other channels 6 and 7 were closed ( Figure 1). Additionally, in order to determine the air change rate inside the car cabin, the flow velocity in the supply air duct (5) was measured using a Testo-435 recorder. The air change rate (ACH) n was determined according to the following relationship: where Q (m 3 /s) is the volumetric air flow, V (m 3 ) is the volume of the passenger car cabin, t (s) is the time, and um (m/s) is the average air velocity in the supply air duct. Table 1 presents the characteristics of the persons in the car, and the body surface area is determined from the DuBois formula: The temperature and relative humidity were measured in the center of the passenger car between the front seats. Table 2 presents the measurement series number, car designation, volume of the passenger car cabin, the number of people in the car, the number of adults and children aged six, the driving time, air change rate, and the date of measurement.   Figure 1. The main elements included in the passenger car ventilation system: 1-air intake, 2-supply fan, 3-air filter, 4-heat exchangers (cooler, heater), 5-air supply duct under the glass, 6-duct air supply towards the driver and passenger, 7-air supply channel towards the driver's and passengers' legs, 8-car windscreen.
The relative humidity and temperature in the tested passenger cars were recorded while driving the car with the windows closed. The research was conducted in Poland in the spring and the summer when mechanical ventilation in passenger cars is intensively exploited. Measurements were made outside of an urban area on a national road, where the car's velocity varied from 70 to 90 km/h. The relative humidity and temperature were measured outside the car before and after driving the car as well as in the supply air duct and inside the passenger car cabin while driving the car using a Testo-435 recorder, whose parameters are shown in [39]. During the measurement, the airflow was set only before the windscreen in channel 5; the other channels 6 and 7 were closed ( Figure 1). Additionally, in order to determine the air change rate inside the car cabin, the flow velocity in the supply air duct (5) was measured using a Testo-435 recorder. The air change rate (ACH) n was determined according to the following relationship: where . Q (m 3 /s) is the volumetric air flow, V (m 3 ) is the volume of the passenger car cabin, t (s) is the time, and u m (m/s) is the average air velocity in the supply air duct. Table 1 presents the characteristics of the persons in the car, and the body surface area is determined from the DuBois formula: The temperature and relative humidity were measured in the center of the passenger car between the front seats. Table 2 presents the measurement series number, car designation, volume of the passenger car cabin, the number of people in the car, the number of adults and children aged six, the driving time, air change rate, and the date of measurement.  Table 2. Basic data of the tested passenger cabins.

Number of Measurement Series
Car Room Volume  Table 3 presents the assessment of thermohumidity conditions in the passenger car cabin based on the average temperature and relative humidity values. The temperature of the air blown into the passenger car cabin is much higher than the outside air temperature, which is caused by heat gains from solar radiation and heat gains from the internal combustion engine. In the case of car A, the largest increase in temperature by 64% was obtained for the air change rate equal to 4.66 1/h (measuring series No. 4), while the smallest increase in the temperature of the supply air in relation to the outside temperature was 16% for the air change rate equal to 83.01 1/h (series measuring No. 8). Air flowing at a lower velocity through heated ventilation ducts achieves a higher temperature. A larger flow through the ventilation system means a lower temperature of the air blown into the car. According to the Mollier chart, air heating in this way runs along the constant absolute humidity line [10,40]. The average absolute humidity inside the car cabin increased compared to the average absolute humidity of the supply air, which is caused by the generation of moisture by the driver and passengers. The increase in absolute humidity inside the car cabin depends on the air change rate, the number of people inside the car, and the temperature inside the car. The highest increase in the average absolute humidity inside the car in relation to the average absolute humidity of the supply air was recorded for measurement series No. 10 (n = 6.21 1/h, T in = 41.6 • C), while the smallest increase in the absolute humidity inside the car was recorded for measurement series No. 8 (n = 83.01 1/h, T in = 31.8 • C). High temperatures inside the car are primarily attributed to solar radiation [41]. According to measurements made by [41] in Athens, the maximum temperatures on clear days in a parked car reached up to 76 • C.  1 Air quality assessed on the basis of developed thermohumidity categories in [13,34], 2 possibility of hyperthermia in one-year-old children [40][41][42].

Results and Discussion
Energies 2020, 13, 2987 6 of 14 The assessment of thermohumidity conditions (Table 3) in the passenger car cabin was made in accordance with the guidelines in [13,34] based on the average values of the temperature and the relative humidity in the passenger car cabin. In the presented research, the "intermediate" and "bad" grades prevailed. Bad thermohumidity conditions were caused by too low ventilation efficiency (3-5 series) and by high temperature inside the car (9 series). The "very bad" category was also added to the assessment, which means the possibility of hyperthermia in one-year-old children when the temperature is above 40 • C, based on studies in [40][41][42]. The "very bad" category was received by the series 9, 10, and 11, in which the average temperature was around 40 • C. These take into account the guidelines from [31][32][33][34] where comfort can be obtained at a temperature of 22-24 • C and a humidity of 40-60%; using these guidelines and numbers for our case, comfort is only ensured in series 1 and 2. Figure 2 presents graphic results of absolute humidity measurements for 15 measurement series as a function of time. Universal trends can be observed in the charts. In the case of cars ventilated before driving, an increase in absolute humidity was observed when driving a passenger car. The increase in humidity then depends on the number of passengers and the air change rate as well as the relative humidity of the air supplied to the car, where the greater the number of passengers and the smaller the air change rate and the higher relative humidity of the supply air, then the gradient of relative humidity inside the car is higher. The air change rate in a car has the greatest impact on the absolute humidity increase in a car. For example, for the same two adults in car A for series 3-5 with an average air change rate equal to 4.7 1/h, the increase in absolute humidity is about 23% from the initial value, while for the air change rate equal to 29.7 1/h (series No. 6), the absolute humidity value is approximately constant. The temperature inside the car has a significant impact on the increase in the absolute humidity in the car. The largest increase in absolute humidity by about 30% in relation to the initial value of all measurement series is for series No. 10, namely for the temperature inside the car that is equal to 41.6 • C. It should be emphasized here that the tests were carried out without supply cooling. In the case of cooling, the humidity in a car usually has a decreasing trend as a result of dehumidifying the air [27].
In the case of cars that were not ventilated before driving (series 7, 8, 11) in which people had previously stayed, after switching on the ventilation for a greater flow while driving the car, a decrease in absolute humidity can be observed, where the humidity is removed by mechanical ventilation outside. In measurement series 1-6, 9-10, 12-15, the car was ventilated for two hours before driving. After the initial increase or decrease in absolute humidity, the absolute humidity content stabilizes at a more or less constant level. Figure 3 shows the relationship between the time of absolute humidity stabilization t s and the sum of the weight w of people in the car as well as the air change rate n in the car. The time to stabilize the absolute humidity in the car depends primarily on the air change rate in the car. The longest absolute humidity stabilization time is on average t s = 47 min for the air change rate from 4.17 1/h to 6.21 1/h, while the shortest stabilization time values (2 < t s < 8 min) of absolute humidity are reached above n = 27 1/h. In the case of the sum of the weight of people in the car, the minimum values of absolute humidity stabilization time are achieved for weights over 164 kg and reached the average value of t s = 5.3 min. Figure 4 shows the absolute humidity as a function of time and temperature for the case of opening windows while driving the car, and the absolute humidity graph of the outside air. After opening the windows, a significant decrease in the absolute humidity inside the car can be observed compared to the absolute humidity of the outside air. In the case of the temperature inside the car cabin, a decrease in temperature was also observed. Opening the windows of the car while driving it quickly removes excess moisture from the car cabin; however, it causes great discomfort in the form of a draft, which can be particularly bothersome in the event of low outside air temperatures. When driving a car with open windows (Figure 5), the absolute humidity in the passenger car cabin is about equal to the absolute humidity of the outside air.  (7)).   (7)).    (7)).   All measurements were made on a rainless day, and therefore the impact of precipitation on car humidity is not the subject of the study in this publication.
In order to develop a model of humidity changes as a function of time, a humidity balance was prepared: where V (m 3 ) is the volume of the car, t (h) is the time, ω (g/m 3 ) is the absolute humidity, Qn (g/h) is the humidity supplied from the outside or discharged from the car by supply ventilation, and Qg (g/h) is the humidity generated by the driver and passengers. It should be noted that the moisture buffering by materials, which in many models has been included [35][36][37][38], is omitted here. The humidity that is supplied from the outside or removed from the car by the supply ventilation is described by the following formula: where n is the air change rate, ωa (g/m 3 ) is the absolute external humidity, and ω (g/m 3 ) is the absolute humidity inside the car. The humidity generated by the driver and passengers is described by the following relationship (Equation (5)): where qg (g/(h × person)) is the humidity emitted by one passenger and m is the number of persons (person). All measurements were made on a rainless day, and therefore the impact of precipitation on car humidity is not the subject of the study in this publication.
In order to develop a model of humidity changes as a function of time, a humidity balance was prepared: where V (m 3 ) is the volume of the car, t (h) is the time, ω (g/m 3 ) is the absolute humidity, Q n (g/h) is the humidity supplied from the outside or discharged from the car by supply ventilation, and Q g (g/h) is the humidity generated by the driver and passengers. It should be noted that the moisture buffering by materials, which in many models has been included [35][36][37][38], is omitted here. The humidity that is supplied from the outside or removed from the car by the supply ventilation is described by the following formula: where n is the air change rate, ω a (g/m 3 ) is the absolute external humidity, and ω (g/m 3 ) is the absolute humidity inside the car. The humidity generated by the driver and passengers is described by the following relationship (Equation (5)): Energies 2020, 13, 2987 10 of 14 where q g (g/(h × person)) is the humidity emitted by one passenger and m is the number of persons (person). After substituting the relationship in Equation (4) to the formula in Equation (3), the absolute humidity balance takes this form: After the integration of Equation (5), the absolute humidity model is obtained as a function of time t: Figure 2 compares the absolute humidity obtained from the model in Equation (7) with the experiment. The relative error was determined according to the following relationship (Equation (8)): where ω exp is absolute humidity measured in in a passenger car while ω num is absolute humidity determined using Equation (7). The relative error is 9.7%, while the average relative error is 1.5%. Based on the measurement data and the simulation according to the formula in Equation (7), the humidity emissions by an adult and a six-year-old child were determined. The relationship between the humidity generated by an adult and a six-year-old child as a function of the interior temperature of a passenger car cabin is shown in Figure 6. The humidity generated by an adult in the temperature range 21.75-32.2 • C increases linearly and increases from 28 to 48 g/h; when above 32.2 • C, there is a significant increase in humidity generated by an adult from 48 g/h to 180 g/h. It should be noted here that the maximum value of humidity generated by an adult occurred at 42.85 • C, at which there is a risk of hyperthermia [40][41][42]. For five-year-old children, only the temperature range from 21.75 to 34.01 • C was examined, in which the humidity generated by the six-year-old child increased from 8 to 19.5 g/h. In the temperature range from 22 to 29 • C, the humidity generation by adults in a car is on average 4 times greater than by children aged six years old, while the body area of adults is 2.7 larger than the body areas for children aged six years old ( Table 1). As a result of approximation, an approximate formula of the humidity generated by an adult q ga and a child aged six years old q gc was obtained as a function of the temperature inside a passenger car: q ga , q gc (T) = a + bT 0.5 + cT 1.0 + dT 1.5 1 + eT 0.5 + f T 1.0 + gT 1.5 , The shape coefficients of the function of the formula in Equation (8) are presented in Table 4. Based on the formula in Equation (9), it is possible to estimate the humidity generated by humans with a maximum relative error of 1% for an adult and 2% for a six-year-old child.  The regulation of the supply ventilation flow rate in both the A and B cars is done by changing the rotational speed of the impeller and by setting the dampers in the supply channels. During all measurements, the ventilation flow rate was regulated by changing the rotational speed of the impeller, while the dampers were set in one unchanged position. Figure 7 shows the measured power consumption of the fan at the set impeller speed as a function of the air change rate for the supply ventilation for cars A and B. Based on the measured characteristics of the fan power as a function of the air change rate, it is possible to estimate the energy E consumed by the fan over a given time t: During one hour, the fan's energy consumption for the ranges of air change rates at 27-32 1/h, 56-64 1/h, and 83 1/h were on average 19 Wh, 63 Wh, and 150 Wh, respectively. The presented model can be used in modern air quality controllers in the car cabin. The automatic adjustment system based on the number of passengers and the outside air humidity on the basis of the model can adjust the airflow to the car cabin to maintain adequate air humidity in the The regulation of the supply ventilation flow rate in both the A and B cars is done by changing the rotational speed of the impeller and by setting the dampers in the supply channels. During all measurements, the ventilation flow rate was regulated by changing the rotational speed of the impeller, while the dampers were set in one unchanged position. Figure 7 shows the measured power consumption of the fan at the set impeller speed as a function of the air change rate for the supply ventilation for cars A and B. Based on the measured characteristics of the fan power as a function of the air change rate, it is possible to estimate the energy E consumed by the fan over a given time t: During one hour, the fan's energy consumption for the ranges of air change rates at 27-32 1/h, 56-64 1/h, and 83 1/h were on average 19 Wh, 63 Wh, and 150 Wh, respectively. The regulation of the supply ventilation flow rate in both the A and B cars is done by changing the rotational speed of the impeller and by setting the dampers in the supply channels. During all measurements, the ventilation flow rate was regulated by changing the rotational speed of the impeller, while the dampers were set in one unchanged position. Figure 7 shows the measured power consumption of the fan at the set impeller speed as a function of the air change rate for the supply ventilation for cars A and B. Based on the measured characteristics of the fan power as a function of the air change rate, it is possible to estimate the energy E consumed by the fan over a given time t: During one hour, the fan's energy consumption for the ranges of air change rates at 27-32 1/h, 56-64 1/h, and 83 1/h were on average 19 Wh, 63 Wh, and 150 Wh, respectively. The presented model can be used in modern air quality controllers in the car cabin. The automatic adjustment system based on the number of passengers and the outside air humidity on the basis of the model can adjust the airflow to the car cabin to maintain adequate air humidity in the The presented model can be used in modern air quality controllers in the car cabin. The automatic adjustment system based on the number of passengers and the outside air humidity on the basis of the model can adjust the airflow to the car cabin to maintain adequate air humidity in the car.