key: cord-0883589-scwhwurr authors: Zhu, Shengwei; Demokritou, Philip; Spengler, John title: Experimental and numerical investigation of micro-environmental conditions in public transportation buses date: 2010-03-17 journal: Build Environ DOI: 10.1016/j.buildenv.2010.03.004 sha: 1484ce1aef99cbc619c270a1280a97f1aed00bc9 doc_id: 883589 cord_uid: scwhwurr This study examines both numerically and experimentally the micro-environmental conditions in public transportation buses. A Computational Fluid Dynamics (CFD) model was developed and experimentally validated. The developed CFD model was used to calculate the spatial distributions of the mean age and mean residual lifetime of air in the bus environment and evaluate the efficiency of the bus ventilation system. Additionally, the passengers’ exposures to a variety of environmental conditions were experimentally monitored in “real world” field campaigns using the Harvard University shuttle bus system. Real time continuous monitoring systems were used to assess indoor environmental quality in the buses. It was found that CO levels were very low, while the levels of particulate matter varied and were influenced by the ambient air penetrated into the bus through the operation of the doors and the ventilation system. The CO(2) level was found elevated and greatly affected by occupancy conditions. The elevated CO(2) level indicates that the current bus ventilation is insufficient to dilute air pollutants in the bus especially under heavy occupancy conditions. This lack of sufficient ventilation indicates an elevated risk for airborne transmitted diseases in such a popular public transportation system. Commuting occupies around 7% of the daily time of people who live in urban areas [1] with the Americans spending approximately 5.5% of their typical day inside a vehicle [2] . And the bus transportation system is considered to be one of the most available and popular transportation modes around the world [3] . There are a number of published field studies focusing on air quality and thermal comfort in the bus microenvironment. In a study performed in Munich, Germany showed that the concentrations of PM 10 in buses and trams depended on the ambient sources and the road traffic conditions [4] . In a similar study performed in another urban setting, Guangzhou, China, the parameters such as poor vehicle emission controls, poor vehicle maintenance and high-density traffic conditions were linked to poor indoor air quality in buses. Furthermore, it was found that the concentration levels of carbon monoxide (CO) and PM 10 in buses could be lowered substantially by using a HVAC system [5] . This is due to the "protective" effect of mechanically ventilated and air-conditioned buses which results to a lower indoor penetration factor for a variety of air pollutants (CO, PM 2.5 , PM 10 , etc) [6] . Chan [3] also reported that CO 2 concentration could be up to 10 times higher than outdoor concentrations under overcrowded conditions, which is usually the case during rush hours. However, there are a limited number of studies focusing on the effect of ventilation strategies on diluting or removing air pollutants generated inside the buses. In addition, in terms of health effects linked with the unique microenvironment in the buses, nausea, dizziness and respiratory allergies have been identified in previous studies [3, 7] . It also became apparent in recent years, that airborne transmitted diseases, such as TB, Severe Acute Respiratory Syndrome (SARS), Avian Influenza and Swine Influenza (H1N1), may impose a serious global health and economic burden. Therefore, there is a need to further study and understand the transmission of airborne infectious diseases in the indoor environment, and the bus environment might be one of the most challenging ones because of the very close proximity of the people. In recent years, with the rapid development of computer technology and advanced numerical methods, investigators were also able to successfully use Computational Fluid Dynamics (CFD) models to study the microenvironments in a various of public transportation systems, including aircraft [8e10], trucks [11] and cars [12, 13] . Those advanced numerical methods have been proved an effective tool to obtain the detailed spatial and temporal distributions of temperature, velocity and contaminants in microenvironmental settings, which are usually very difficult and expensive to be obtained in experimental measurements. The thrust of this investigation is two-fold: (a) develop and use a CFD model to understand the bus micro-environmental conditions, especially the ventilation effectiveness and air distribution system, and how they may affect the dilution and removal of any internally generated air pollutant; (b) perform an extensive field campaign using the Harvard University shuttle bus system and monitor continuously the indoor environmental quality as well as the personal exposure levels of passengers for a variety of air pollutants including CO, CO 2 , and particulate matter in its various respirable size fractions such as PM 2.5 , PM 10 and ultrafine particles (UFPs). This bus uses the widely used mixing ventilation method. Six linear air supply diffusers are located near the ceiling. On each side of the bus, there is a two-slot diffuser supplying air downwards and a one-slot diffuser supplying air horizontally (Fig. 1) . The total area of the air supply openings is 0.21 m 2 . The exhaust opening with an area of 0.544 m 2 (0.34 m  1.6 m) is placed in the ceiling close to the rear door. In the CFD simulation, the bus was modeled with all of the doors and windows assumed to be closed. The air infiltration through leakage pathways was also omitted. Fig. 2 shows the geometry of the air diffusers. The total air supply rate was 0.54 m 3 /s, which accounts for an air exchange rate of 57.6 ACH assuming no air recirculation in the bus. 1 In addition, solar and long wave radiation as well as occupancy was not considered in the developed model. A 3D CFD model based on the Reynolds-averaged NaviereStokes (RANS) equations [14] was used to simulate the bus airflow, temperature and scalar contaminant fields. The CFD model is based on the conservation equations of mass, momentum, energy, scalar concentration, turbulent kinetic energy, and dissipation rate of turbulent kinetic energy, that govern the transport phenomena in the bus. The aforementioned governing equations can be closed with appropriate boundary conditions at all of the boundaries such as air supply openings, air exhaust opening, and wall surfaces. The detailed boundary conditions are summarized in Table 1 . For air supply openings, the actual values of air temperature, velocity and turbulence intensity were measured in the field experiments, and used as boundary conditions (see following sections). The turbulence scale was calculated as one-half of the width of the linear opening. Fig. 1(c) shows the location and direction of the supplied airflows in cross section BeB 0 . At the surfaces of floor, windows and lights, the temperatures were fixed based on the data obtained in the field experiment. The rest of wall surfaces were assumed adiabatic. Moreover, the standard wall function [15] was used for the near wall boundary layer. In this study, the spatial distributions of the mean age and mean residual lifetime of air [16] were used as the indices to estimate the ventilation efficiency in the bus. The mean age of air is defined as the mean time that it takes for the air from the supply openings to be transported at a specified location; the mean residual lifetime of air is defined as the mean time that it takes for the air at certain location to be exhausted. Numerically, these distributions can be obtained using the mathematical methods of SVE3 and SVE6 (Scale for Ventilation Efficiency 3 and 6, see Table 2 ) [17e19], under the assumption that the air pollutant ("tracer gas") is generated uniformly and continuously in the bus. Based on the so-called SVE3 method, the mean age of air can be calculated by solving the following transport equation [20] . If q 0 is assumed to be generated uniformly and continuously throughout the bus, the air mass from a supply opening is gradually contaminated as it is regarded proportional to the time elapsed from the time the air mass leaves the supply opening until it reaches the point. The total amount of tracer gas generated in the room divided by the airflow rate of ventilation is termed a "perfectly mixed concentration C 0 (kg/kg)". Under a steady state condition, this is equivalent to the concentration averaged over the exhaust opening. Because the age of air averaged over the exhaust opening is equivalent to the inverse of the air exchange rate, the distribution of trace gas C in the room can be converted to the distribution of the mean age of air. Similarly based on the so-called SVE6 method, the mean residual lifetime of air is computed by applying the same computational procedure as SVE3, except that the time sequence of the flow field is reversed. The CFD numerical grid system used in the model consists of 126,838 triangular surface meshes and 1,062,098 tetrahedral fluid cells, which were created using Gridgen V15.10 at a growth rate of 1.1. A good grid was ensured with over 99.9% of the meshes of less than 0.8 in EquiAngle Skewness, which represented the maximum ratio of the cell's inner angle to the angle of an equilateral element. 1 Note: The total air supply rate of 0.54 m 3 /s was calculated by mutiplying the velocity of the supply airflow by the total area of the air supply openings. It resulted in an air exchange rate of 57.6 ACH based on the assumption that there was 100% outside air supplied in the bus and no air recirculation which may not be the case in real world operations. In addition, the previous studies reported high ACH in automobiles. For example, the air exchange rate was measured between 36.3 and 47.5 ACH in an automobile with windows and doors closed and ventilated on 100% fresh air under stationary conditions [29] , and averagely 92 ACH at a speed of 105 Km/h [30] . Therefore, the air exchange rate of 57.6 ACH agrees with the published studies. The commercial CFD software, Star-CD V4.08, was used to solve the governing equations together with the standard ke3 model and implicit SIMPLE algorithm [21] . The finite volume method with the first-order upwind scheme (UD) was adopted for discretizing the governing equations [22] . In order to validate the developed CFD method, a field experiment was also conducted under stationary conditions using a bus Table 3 ). The averages of the total 140 data points (for 70 min) were recorded. The velocities and temperatures of the supply airflow at the linear diffusers were also measured using a hot-wire anemometer (TSI 964, see Table 3 ) at an interval of 1 s for the duration of 3 min in order to obtain the boundary conditions for the air supply openings. The measurements were performed at four bus locations (in the front and rear compartment at each side of the bus). The average of the 180 data points was used as the value for that location. In addition, turbulence intensity was calculated automatically as the percent of the standard deviation of the turbulent velocity fluctuations to the mean velocity over the measuring period at each location [23] . The averages of those results were also used as the boundary conditions for air supply opening. It should be addressed that while the average values for temperature and velocity were used at the walls and supply inlets, in real life these parameters vary as a function of time based on weather and other operating conditions. The Harvard University shuttle bus system was used for our field investigation study. The Harvard University Shuttle Service is the main transport mode for the students, staff and faculty commuting throughout the Cambridge and Allston campuses of Harvard University. It is estimated that over 25,000 passengers per week use the Harvard University Shuttle Service. Two different routes were selected. Both routes have similar traffic conditions of primary and secondary streets. A typical route consists of four stops and takes approximately 15 min. At each stop, it usually takes approximately 30 s for boarding/unboarding passengers. The age profile of passengers is also unique since 85% of the passengers are under 25 years old on our survey investigation. Four full day field campaigns were performed in the fall of 2008 (November 17th, 19th, 25th and 26th). All buses were Daimler SLF 200, same as that used in the field experiments for CFD model validation. The buses were air-conditioned, and all of their windows were fully closed when in operation. Outdoor air infiltrated through the doors and the ventilation system. The field campaigns started early in the morning before 9:00AM, with a lunch break around noon and continued until 16:30PM in the afternoon. It is worth pointing out that typical traffic patterns in the Cambridge campus area were reflected in our field campaigns. The drivers were also asked to operate as usual without intervening in the temperature setup and the ventilation settings in the bus. Continuous instruments were used to monitor indoor environmental conditions in the bus. Air temperature, relative humidity and the concentration levels of CO 2 , CO and particulate matters in three size fractions, PM 10 and PM 2.5 , and ultrafine particles (UFPs) were continuously measured in the bus. Real time measurements for the above IEQ parameters were recorded at an interval of 10 s. The instruments were placed in two mesh boxes made by coarse wire, which were hanged at the shoulders of two of our field personnel at a height between 0.6 and 1.1 m from the floor. The field personnel remained standing in the middle of the front compartment for the whole trip as shown in Fig. 1(a) . Moreover, since there was only one particle counter used in the Table 2 Scales for ventilation efficiency 3 and 6. Age of air at position X SVE6(X) Residual lifetime at position X C x 0 (X) Concentration at position X where contaminant is uniformly generated throughout a room in total generation rate q Table 3 summarizes the instruments used in the continuous measurements of the aforementioned IEQ parameters. Fig. 3 shows the spatial distribution of air velocity in the bus. As shown in cross section AeA 0 , the air close to the front windshield was heated and moved upward, an indication of a "drafty" condition in the proximity of the driver. Moreover, an eddy was generated in the mid upper region as shown in cross section BeB 0 . In cross section DeD 0 , a large eddy was shown at each side of the bus. These eddies are indicative of the highly turbulent nature of the flow in the bus environment. However, overall, the air velocities were mostly distributed under 0.3 m/s in the bus. This is almost satisfactory based on the ASHRAE Standard 55-2004 [24] , which has a limit of 0.2 m/s for air velocity. Fig. 4 shows the spatial distribution of air temperature in the bus environment. Air temperatures were within AE1 C of the set point (20.2 C) in the passenger compartment with the exception of the area close to the windows and floor. As shown in cross sections BeB 0 and CeC 0 , large air temperature gradients were present close to the floor, primarily because of the high thermal losses and the very poor air circulation under the seats. The same large temperature gradients were also shown close to the front windshield and the floor in the driver's space. there was not a great agreement between measured and simulated results. This is expected since it is a known limitation of the standard wall function used in the model for the wall surfaces covered by a low Reynolds number turbulent flow. This might be also attributed to a cold air infiltration close to the floor, which was not considered in the CFD model. Fig. 6 shows the spatial distribution of the mean age of air in the bus. For the driver's breathing zone, the mean age of air was calculated to be generally in 62e86.8 s. In the front compartment, it was mainly in the range of 49.6e62 s. In the rear compartment, it was mostly less than 62 s (with the exception of the area involved in the developed eddy as shown in Fig. 3(d) ). It is worth to point out that in general, fresh air can reach most of the area covering passengers' breathing zone within 74.4 s. Similarly, Fig. 7 shows the spatial distributions of the mean residual lifetime of air. It was found that the mean residual lifetime of air varied significantly from 0 to over 111.6 s in the bus. As expected, it was minimum in the area directly under the exhaust opening (less than 12.4 s), and increased very quickly with the distance from the air exhaust opening. In conclusion, the polluted air in passengers' breathing zone will be exhausted within 2 min. Table 4 shows the weather conditions occurred during the experimental campaigns. Table 5 summarizes the measured air temperature and relative humidity in the buses for each day of the field campaign. Figs. 8 and 9 show the bus air temperature and relative humidity respectively, as a function of time. Here, only the data for a representative campaign day (Nov. 19th) were plotted in the figures since the results were similar on other days. Moreover, please also note that the number of passengers as a function of time was also presented in the figures. As shown in Tables 4 and 5 , the air temperature in the bus was independent of outdoor weather conditions and the passenger occupancy. Additionally, the door operation showed little influence on the air temperature in the bus as shown in Fig. 8 . This is an indication that the bus HVAC system can maintain really well the temperature set point. On the contrary, as shown in Tables 4 and 5 , relative humidity in the bus was affected greatly by the weather conditions with its daily mean value ranged from 49.8% when it was raining to 12.9% when it was sunny. The high relative humidity on the rainy day was caused by the water drops taken into the bus by the passengers or entering the bus when the doors were opened to load/unload the passengers. In addition, as indicated in Fig. 9 , relative humidity was greatly affected by the opening of bus doors during loading/unloading of the passengers. Also, it was as expected that relative humidity depended on the occupancy load, primary because of the latent heat generation from the human occupancy. The thermal comfort range for winter as it is recommended by ASHRAE standard 55-2004 [24] was presented in Fig. 10 . This was made possible by plotting the experimental data points of the hourly mean air temperature and relative humidity measured in the bus throughout the field campaigns. It is necessary to point out that this is under the assumption that the measured air temperature in the bus is equal to the operative temperature and the clothing conditions of the passengers and drivers is 1 Clo as stated in the ASHRAE standard. Operative temperature was not measured in the field experiments and it can be higher than the air temperature if thermal radiation is considered. Therefore, Fig. 10 is considered just an approximation of the thermal comfort conditions in the bus. As shown in Fig. 10 , most of the results were distributed at the "warmer side" of the thermal comfort zone but still in the comfort zone based on the ASHRAE standard 55-2004. Table 6 summarizes the measured concentration levels of all the air pollutants monitored continuously in the bus (CO, CO 2 , PM 2.5 , PM 10 and UFPs). Similarly, Fig. 11 also illustrates the CO 2 concentration levels and the occupancy load as a function of time. Additionally, Figs. 12 ande15 show the concentration levels of CO, PM 2.5 , PM 10 and UFPs as a function of time, respectively. The door open/ close status was also indicated on those figures. In addition, also only the data for a representative campaign day were plotted in the figures since the results were similar on other days. As it was shown in the aforementioned figures, the CO 2 concentration level in the bus varied widely during the field campaigns. It usually varied from around 1000 ppm, and sometimes reached the levels of over 2000 ppm. It is worth to point out that based on the guidelines of the National Institute for Occupational Safety and Health (NIOSH) for indoor environments, CO 2 concentration levels higher than 1000 ppm is an indication of inadequate ventilation [25] . Furthermore, Fig. 11 reveals that CO 2 concentration levels in the bus were linked with the passenger occupancy. As it was expected, the CO 2 concentration level increased with higher occupancy load. Moreover, the CO 2 concentration level decreased sharply when the doors were open and increased soon after doors were closed in most cases due to the dilution effect from infiltrating air. The CO concentration levels were below 1.9 ppm, much lower than the 8-hour average limit of 9 ppm, which was recommended by the US National Ambient Air Quality Standards (NAAQS). Additionally, as it is shown in Fig. 12 , the CO concentration level generally reached its peak values when the doors were open to load and unload passengers. This is a clear indication of the influence of the traffic conditions in the street on the CO concentration level. The mass concentrations of PM 2.5 and PM 10 were always below 220 mg/m 3 , with the daily mean values to be much lower than the recommendation by the NAAQS (35 mg/m 3 for PM 2.5 and 150 mg/m 3 for PM 10 ). As it is shown in Fig. 13 , PM 2.5 concentration generally did not exceed 20 mg/m 3 except for the "heavy traffic hours" in the morning and afternoon (before 9:30AM and after 14:00PM). It is worth pointing out that PM 2.5 concentrations during peak hours were found to be an order of magnitude higher. This is due to the high traffic conditions during peak hours and PM 2.5 particles are well correlated with mobile sources (traffic conditions). Similar to the CO concentration level, the PM 2.5 concentration level reached its highest values when the doors were open to load and unload passengers, a clear indication of particle penetration from outdoors. Similar trend for the PM 10 concentration level is shown in Fig. 14 . UFPs concentration level also varied widely over time, and reached its peak values when the doors were open, a clear indication of the influence of the traffic conditions. This is expected since UFPs are primarily emitted from combustion engines related to mobile sources such as cars, buses, etc. The daily mean UFPs concentration levels were found to be in the range of 4.0  10 4 e5.7  10 4 pt/cm 3 . This is in line with the findings from the similar studies in other cities in the US (9.5  10 3 e5.3  10 4 pt/ cm 3 in Ann Arbor, MI, 3.0  10 4 e7.5  10 4 pt/cm 3 in Chicago, and 7.0  10 3 e5.0  10 4 pt/cm 3 in Atlanta [26] ). As it was shown in both the simulation results and the experimental results previously presented, there was overall thermal satisfaction based on the ASHRAE 55-2004. However, the low relative humidity levels below 25% on sunny days documented in the field campaigns raise a health concern in terms of the airborne transmission of certain viruses. For example, according to the previous animals studies [27, 28] , the most favorable relative humidity levels for the airborne transmission of influenza virus was in the range of 17e35%. Our field campaign results imply that the low relative humidity found in this microenvironment may elevate the risk for influenza transmission. As outlined in the previous result section, the concentration levels of CO, PM 2.5 and PM 10 were found in compliance with the existing standards for indoor air quality. Their concentration levels were always very small when the doors were closed and generally reached their peak values when the doors were open. Since there are no indoor sources of CO, the only reasonably explanation for this change is that the use of mechanical ventilation in the bus provided a "shielding effect" in terms of particulate matter and CO. This agrees with the results from the previous studies on air-conditioned buses [6] . However, the CO 2 concentration levels were found to be usually high and not in compliance with indoor environmental guidelines and standards. This is an indication of inadequate ventilation that is linked to an increased risk for airborne transmitted diseases. The simulation results also confirmed the lack of ventilation efficiency and found the mean residual lifetime of air in the bus for the breathing zone of passengers to exceed 112 s in some areas of the bus. The developed CFD model can be used to predict the ventilation efficiency in the passenger's breathing zone and evaluate alternate operational scenarios related to the ventilation system. The simulation results of mean age and mean residual lifetime of air indicate very good ventilation efficiency in an empty bus with 100% fresh air supply (no air recirculation), under the conditions different from those in the "real world" field campaigns on IEQ in the buses. Based on both our numerical and field investigation on the bus environment, there is a need to further study and identify ventilation strategies and alternate air distribution methods in order to improve the ventilation efficiency in this important microenvironment. Additionally, due to the very close proximity of people in the bus environment, supplementary air filtration and air disinfection methods such as upper room UV irradiation might be necessary to be explored as a mean to minimize the risk of airborne infectious diseases. We plan to use the developed CFD model to parametrically analyze the impact of the operational conditions of ventilation system, such as air recirculation, air distribution method, etc. The field experiments indicated that the thermal comfort conditions were found on the warmer side of the ASHRAE recommended comfort zone for winter conditions. It was also found that the ventilation system provided a "shielding effect" and protected the passengers from the harmful exposures to ambient air pollutants infiltrating indoors, such as CO and particulate matter. The high CO 2 concentration levels found in the field campaigns is a good indicator of insufficient venitlation in the buses. These results reveal an increased risk for the airborne transmitted diseases. Therefore, alternate ventilation strategies and possible air filtration/purification systems may be needed in order to safeguard the health of people using this public transportation system. Activity patterns of Californians: use of and proximity to indoor pollutant sources The national human activity pattern survey (NHAPS): a resource for assessing exposure to environmental pollutants Commuter exposure and indooreoutdoor relationships of carbon oxides in buses in Hong Kong Dust exposure in Munich public transportation: a comprehensive 4-year survey in buses and trams Exposure level of carbon monoxide and respirable suspended particulate in public transportation modes while PM 10 mass concentration level and door operation condition in the bus as a function of time UFPs particle concentration level and door operation condition in the bus as a function of time (Nov.19th). commuting in urban area of Guangzhou Influence of in-tunnel environment to in-bus air quality and thermal condition in Hong Kong Combined comfort model of thermal comfort and air quality on buses in Hong Kong Numerical simulation of airflow and airborne pathogen transport in aircraft cabins e part I: numerical simulation of the flow field Numerical simulation of airflow and airborne pathogen transport in aircraft cabins e part II: numerical simulation of airborne pathogen transport Novel air distribution systems for commercial aircraft cabins Prediction of indoor thermal comfort in vehicle by a numerical thermoregulation model and CFD Analysis of flow patterns and heat transfer in generic passenger car mini-environment Ventilation characteristics of modeled compact car (part 2) Estimation of local ventilation efficiency and inhaled air quality Progress in the development of Reynolds stress turbulence closure The numerical computation of turbulent flows The use of Moments for ventilation assessing air quality in ventilated room New ventilation efficiency scales based on spatial distribution of contaminant concentration aided by numerical simulation New scales for ventilation efficiency and their application based on numerical simulation of room airflow New scales for evaluating ventilation efficiency as affected by supply and exhaust openings based on spatial distribution of contaminant Study on inhaled air quality in a personal air-conditioning environment using new scales of ventilation efficiency Numerical heat transfer and fluid flow Star-CD Version 4.08-Methodology Turbulence intensity measurements. Application Note AF-141 of TSI Instruments Ltd Standard 55-2004: thermal environment conditions for human occupancy. Atlanta: American Society of Heating, Refrigerating and Air Conditioning Engineers Indoor environmental quality: building ventilation. National Institute for Occupational Safety and Health A multi-city investigation of the effectiveness of retrofit emissions controls in reducing exposutres to particulate matter in school buses. Clean Air Task Force Experimental airborne influenza infection. I. Influence of humidity on survival of virus in air Influenza virus transmission is dependent on relative humidity and temperature Measurement of air exchange rate of stationary vehicles and estimation of in-vehicle exposure Simultaneous measurement of ventilation using tracer gas techniques and VOC concentrations in homes, garages and vehicles We express our deepest gratitude to Harvard University Transport Service for accommodating and assisting during the field campaigns. We would like to acknowledge the contribution of Xuan Du and Anna Algevik, the students of Harvard Extension School, who participated and helped in the study. And we are very grateful for Dr. Kato Shinsuke, from the Institute of Industrial Science, at the Univeristy of Tokyo, who shared his CFD resources with us.