A CFD study of the transport and fate of airborne droplets in a ventilated office: The role of droplet−droplet interactions

• Coulomb and Lennard−Jones forces were considered for droplet interactions. • The net droplet interactions were repulsive. • Repulsive droplet interactions increased the transport of droplets. • Repulsive droplet interactions significantly modified the fate of droplets. Previous studies reported that specially designed ventilation systems provide good air quality and safe environment by removing airborne droplets that contain viruses expelled by infected people. These water droplets can be stable in the environment and remain suspended in air for prolonged periods. Encounters between droplets may occur and droplet interactions should be considered. However, the previous studies focused on other physical phenomena (air flow, drag force, evaporation) for droplet transport and neglected droplet interactions. In this work, we used computational fluid dynamics (CFD) to simulate the transport and fate of airborne droplets expelled by an asymptomatic person and considered droplet interactions. Droplet drag with turbulence for prediction of transport and fate of droplets indicated that the turbulence increased the transport of 1 μm droplets, whereas it decreased the transport of 50 μm droplets. In contrast to only considering drag and turbulence, consideration of droplet interactions tended to increase both the transport and fate. Although the length scale of the office is much larger than the droplet sizes, the droplet interactions, which occurred at the initial stages of release when droplet separation distances were shorter, had a significant effect in droplet fate by considerably manipulating the final locations on surfaces where droplets adhered. Therefore, it is proposed that when an exact prediction of transport and fate is required, especially for high droplet concentrations, the effects of droplet interactions should not be ignored.     

Extension of the EUROS Integrated Air Quality Model to Fine Particulate Matter by Coupling to CACM/MADRID 2

The European Operational Smog (EUROS) integrated air quality modelling system has been extended to model fine particulate matter (PM). From an extended literature study, the Caltech Atmospheric Chemistry Mechanism and the Model of Aerosol Dynamics, Reaction, Ionisation and Dissolution were selected and recently coupled to EUROS. Currently, modelling of mass and chemical composition of aerosols in two size fractions (PM_2.5 and PM_10–2.5) is possible. The chemical composition is expressed in terms of seven components: ammonium, nitrate, sulphate, elementary carbon, primary inorganic compounds, primary organic compounds and secondary organic compounds. Calculated PM₁₀ concentrations and chemical composition are presented for two summer months of the year 2003 (1 July to 31 August).     

Size distributions of fine and ultrafine particles in the city of Strasbourg: correlation between number of particles and concentrations of NO(x) and SO(2) gases and some soluble ions concentration determination

An Electrical Low Pressure Impactor (ELPI) was used during spring and autumn 2003 in the centre of Strasbourg for the measurement of atmospheric aerosols size distribution. The concentration of NO(x) and SO(2) in air was simultaneously measured with specific analysers. Samples were collected in the range 0.007-10 microm in equivalent aerodynamic diameter size. Number distributions are representative of a pollution originating from urban traffic with a particle size distribution exhibiting a nucleation mode below 29 nm and an accumulation mode around 80 nm in size. A mean particle density equal to 39000+/-35000 total particles per cm(3) with a size ranging from 7 to 10 microm was obtained after a sampling period of 2 weeks in spring. About 86.9% of the number of particles have an aerodynamic diameter below 0.1 microm and 13.1% between 0.1 and 1 microm. Correlation coefficients between the number of particles impacted on each ELPI plate and gas concentrations (SO(2) and NO(x)) showed that the numbers of particles with diameter between 0.10 and 0.62 microm are highly related to the NO(x) concentration. This result indicates that particles are traffic induced since NO(x) is mainly emitted by cars as shown by measurements on various sites. Particles are less clearly correlated to the SO(2) concentration. Particle analysis on different ELPI plates for a sampling period of 2 weeks in autumn showed high level of soluble NO(3)(-), SO(4)(2-) and NH(4)(+) ions. Indeed, up to 90% b.w. of these three species were found in the particle range 0.1-1 microm. The formation of particulate NH(4)NO(3) is favoured by high NO(x) concentration, which induces the formation of gaseous HNO(3).     

The transformation of outdoor ammonium nitrate aerosols in the indoor environment

Recent studies associate particulate air pollution with adverse health effects; however, the exposure to indoor particles of outdoor origin is not well characterized, particularly for individual chemical species. We conducted a field study in an unoccupied, single-story residence in Clovis, California to provide data and analyses to address issues important for assessing exposure. We used real-time particle monitors both outdoors and indoors to quantify nitrate, sulfate, and carbon particulate matter of particle size 2.5 μm or less in diameter (PM-2.5). The results show that measured indoor ammonium nitrate concentrations were significantly lower than would be expected based solely on penetration and deposition losses. The additional reduction can be attributed to the transformation indoors of ammonium nitrate into ammonia and nitric acid gases, which are subsequently lost by deposition and sorption to indoor surfaces. A mass balance model that accounts for the kinetics of ammonium nitrate evaporation was able to reproduce measured indoor ammonium nitrate and nitric acid concentrations, resulting in a fitted value of the deposition velocity for nitric acid of 0.56 cm s⁻¹. The results indicate that indoor exposure to outdoor ammonium nitrate in Central Valley of California are small, and suggest that exposure assessments based on total particle mass measured outdoors may obscure the actual causal relationships for indoor exposure to particles of outdoor origin.     

Prediction of minimum ignition energy of aerosols using flame kernel modeling combined with flame front propagation theory

Flammable aerosols have created many fire and explosion hazards in the process industry, but the flammability of aerosols has not been fully understood. The minimum ignition energy has been widely used as an indicator for flammability of combustible mixtures, but the amount of experimental data on the minimum ignition energy of aerosols is very limited. In this work, the minimum ignition energy of tetralin aerosols is predicted using an integrated model. The model applies the flame front propagation theory in aerosol systems to the growth of the flame kernel, which was created during the spark discharge in the ignition process. The aerosol minimum ignition energy was defined as the minimum level of energy in the initial flame kernel to maintain the kernel temperature above the minimum ignition temperature of 1073 K specific for tetralin aerosols during the kernel growth. The minimum ignition energy obtained in the model is influenced by the fuel-air equivalence ratio and the size of the aerosol droplets. For tetralin aerosols of 40 μm diameter, E_min decreases significantly from 0.32 mJ to 4.3 × 10 e⁻³ mJ when the equivalence ratio rises from 0.57 to 1.0. For tetralin aerosols of 0.57 equivalence ratio, E_min increases from as 0.09 mJ to 0.32 mJ when the droplet diameter rises from 10 μm to 60 μm. The trends are in agreement with previous experimental observations. The method used in current work has the potential to prediction of the minimum ignition energy of aerosol.     

Particle-size distribution and gas/particle partitioning of atmospheric polybrominated diphenyl ethers in urban areas of Greece

Ambient concentrations, gas/particle partitioning and particle-size distribution of polybrominated diphenyl ethers (PBDEs) were investigated in two urban areas (Athens and Heraklion) of Greece. Atmospheric (gas + particle) concentrations of ∑PBDE varied from 21 to 30 pg m⁻³ in the center of Athens and from 4 to 44 pg m⁻³ in the suburbs of Heraklion. A predominance of particulate PBDEs was observed in Athens (71-76% in particles), whereas the opposite was evident in Heraklion (69-92% in gas phase). In both urban areas, PBDE particle-size distribution featured a distinct enrichment in smaller particles. A similar trend was also observed in aerosols of a background marine site. For all sampling sites, more than 46% of ∑PBDE was associated with particles of <0.57 μm in diameter. Our results imply that particulate PBDEs may have long atmospheric residence time and they may be capable of reaching the deeper parts of the human respiratory system.     

Estimating fugitive dust emission rates using an environmental boundary layer wind tunnel

Emissions from fugitive dust due to erosion of “natural” wind-blown surfaces are an increasingly important part of PM₁₀ (particulate matter with sizes of 10 μm aerodynamic diameter) emission inventories. These inventories are particularly important to State Implementation Plans (SIP), the plan required for each state to file with the Federal government indicating how they will comply with the Federal Clean Air Act (FCAA). However, techniques for determining the fugitive dust contribution to over all PM₁₀ emissions are still in their developmental stages. In the past, the methods have included field monitoring stations, specialized field studies and field wind-tunnel studies. The measurements made in this paper allow for systematic determination of PM₁₀ emission rates through the use of an environmental boundary layer wind tunnel in the laboratory. Near surface steady-state concentration profiles and velocity profiles are obtained in order to use a control volume approach to estimate emission rates. This methodology is applied to soils retrieved from the nation's single largest PM₁₀ source, Owens (dry) Lake in California, to estimate emission rates during active storm periods. The estimated emission rates are comparable to those obtained from field studies and lend to the validity of this method for determining fugitive dust emission rates.     

Performance evaluation of NOAA-EPA developmental aerosol forecasts

To aid air quality model development and assess air quality forecasts, the Meteorological Development Laboratory (MDL) provided categorical verification metrics for developmental aerosol predictions. The National Air Quality Forecasting Capability (NAQFC) generated 48 h (of) gridded hourly developmental predictions for the lower 48 states (CONUS) domain in 12 km horizontal spacing. The NAQFC uses the North American Mesoscale (NAM) model with EPA’s Community Multiscale Air Quality (CMAQ) model to produce predictions of ground level aerosol concentrations. We used bilinear interpolation to calculate predicted daily maximum values at the location of the observation sites. We compared these interpolated predicted values to the observed daily maximum to produce 2 × 2 contingency tables, with a threshold of 40 μg/m³ during the months of March–August, 2007. The model showed some degree of skill in predicting aerosol exceedances. These results are preliminary as the NAQFC model for aerosol prediction is in the developmental stage. A more comprehensive performance evaluation will be accomplished in 2008, when more data become available. Our verification metrics included categorical analyses for Fraction Correct (FC) or percent correct (FC × 100), Threat Score (TS) or Critical Success Index (CSI), Probability of Detection (POD), and the False Alarm Rate (FAR), Mean Absolute Error (MAE) and mean algebraic error or bias, where bias is forecast minus observation. Graphic products included weekly statistics for the CONUS displayed in the form of bar charts, scatterplots, and graphs. In addition, we split the CONUS into six geographic regions and provided regional statistics on a monthly basis. MDL produced spatial maps of daily 1-h maximum predicted aerosol values overlaid with the corresponding point observations. MDL also provided spatial maps of the daily maximum of the 24-h running average. We derived the 24-h running average from the 1-h average predicted aerosol values and observations.     

Simulation of aerosol direct radiative forcing with RAMS-CMAQ in East Asia

The air quality modeling system RAMS-CMAQ is developed to assess aerosol direct radiative forcing by linking simulated meteorological parameters and aerosol mass concentration with the aerosol optical properties/radiative transfer module in this study. The module is capable of accounting for important factors that affect aerosol optical properties and radiative effect, such as incident wave length, aerosol size distribution, water uptake, and internal mixture. Subsequently, the modeling system is applied to simulate the temporal and spatial variations in mass burden, optical properties, and direct radiative forcing of diverse aerosols, including sulfate, nitrate, ammonium, black carbon, organic carbon, dust, and sea salt over East Asia throughout 2005. Model performance is fully evaluated using various observational data, including satellite monitoring of MODIS and surface measurements of EANET (Acid Deposition Monitoring Network), AERONET (Aerosol Robotic Network), and CSHNET (Chinese Sun Hazemeter Network). The correlation coefficients of the comparisons of daily average mass concentrations of sulfate, PM2.5, and PM10 between simulations and EANET measurements are 0.70, 0.61, and 0.64, respectively. It is also determined that the modeled aerosol optical depth (AOD) is in congruence with the observed results from the AERONET, the CSHNET, and the MODIS. The model results suggest that the high AOD values ranging from 0.8 to 1.2 are mainly distributed over the Sichuan Basin as well as over central and southeastern China, in East Asia. The aerosol direct radiative forcing patterns generally followed the AOD patterns. The strongest forcing effect ranging from −12 to −8 W m⁻² was mainly distributed over the Sichuan Basin and the eastern China’s coastal regions in the all-sky case at TOA, and the forcing effect ranging from −8 to −4 W m⁻² could be found over entire eastern China, Korea, Japan, East China Sea, and the sea areas of Japan