Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KM-GAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KM-GAP is based on the PRA model framework (Pöschl-Rudich-Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modelled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmo- spheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270 K is close to unity. Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for eðcient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone

Equilibration timescale of atmospheric secondary organic aerosol partitioning

[1] Secondary organic aerosol (SOA) formed from partitioning of oxidation products of anthropogenic and biogenic volatile organic compounds (VOCs) accounts for a substantial portion of atmospheric particulate matter. In describing SOA formation, it is generally assumed that VOC oxidation products rapidly adopt gas-aerosol equilibrium. Here we estimate the equilibration timescale, τ_(eq), of SOA gas-particle partitioning using a state-of-the-art kinetic flux model. Τ_(eq) is found to be of order seconds to minutes for partitioning of relatively high volatility organic compounds into liquid particles, thereby adhering to equilibrium gas-particle partitioning. However, τ_(eq) increases to hours or days for organic aerosol associated with semi-solid particles, low volatility, large particle size, and low mass loadings. Instantaneous equilibrium partitioning may lead to substantial overestimation of particle mass concentration and underestimation of gas-phase concentration

Fatigue Life Prediction of Welded Joint by Microstructure-based Simulation

This paper proposes a numerical framework to predict fatigue life on welded joints by integrating several computational techniques. The framework consists of five steps: i) materials properties estimation; ii) welding simulation using thermo-mechanical finite element method; iii) macroscopic stress field analysis under cyclic loading; iv) mesoscopic stress field analysis using crystal plasticity finite element method (CPFEM); v) analysis of fatigue crack growth. The total number of cycles to failure is eventually obtained by the sum of initiation life calculated by CPFEM and propagation life calculated by X-FEM. A fatigue life of butt joint is evaluated by the proposed method. The results demonstrated the possibility of evaluating the fatigue life and its scattering by the proposed framework

Mass Accommodation of Water: Bridging the Gap Between Molecular Dynamics Simulations and Kinetic Condensation Models

The condensational growth of submicrometer aerosol particles to climate relevant sizes is sensitive to their ability to accommodate vapor molecules, which is described by the mass accommodation coefficient. However, the underlying processes are not yet fully understood. We have simulated the mass accommodation and evaporation processes of water using molecular dynamics, and the results are compared to the condensation equations derived from the kinetic gas theory to shed light on the compatibility of the two. Molecular dynamics simulations were performed for a planar TIP4P-Ew water surface at four temperatures in the range 268–300 K as well as two droplets, with radii of 1.92 and 4.14 nm at T = 273.15 K. The evaporation flux from molecular dynamics was found to be in good qualitative agreement with that predicted by the simple kinetic condensation equations. Water droplet growth was also modeled with the kinetic multilayer model KM-GAP of Shiraiwa et al. [ Atmos. Chem. Phys. 2012, 12, 2777]. It was found that, due to the fast transport across the interface, the growth of a pure water droplet is controlled by gas phase diffusion. These facts indicate that the simple kinetic treatment is sufficient in describing pure water condensation and evaporation. The droplet size was found to have minimal effect on the value of the mass accommodation coefficient. The mass accommodation coefficient was found to be unity (within 0.004) for all studied surfaces, which is in agreement with previous simulation work. Additionally, the simulated evaporation fluxes imply that the evaporation coefficient is also unity. Comparing the evaporation rates of the mass accommodation and evaporation simulations indicated that the high collision flux, corresponding to high supersaturation, present in typical molecular dynamics mass accommodation simulations can under certain conditions lead to an increase in the evaporation rate. Consequently, in such situations the mass accommodation coefficient can be overestimated, but in the present cases the corrected values were still close to unity with the lowest value at ≈0.99

Discontinuities in hygroscopic growth below and above water saturation for laboratory surrogates of oligomers in organic atmospheric aerosols

Discontinuities in apparent hygroscopicity below and above water saturation have been observed for organic and mixed organic–inorganic aerosol particles in both laboratory studies and in the ambient atmosphere. However, uncertainty remains regarding the factors that contribute to observations of low hygroscopic growth below water saturation but enhanced cloud condensation nuclei (CCN) activity for a given aerosol population. Utilizing laboratory surrogates for oligomers in atmospheric aerosols, we explore the extent to which such discontinuities are influenced by organic component molecular mass and viscosity, non-ideal thermodynamic interactions between aerosol components, and the combination of these factors. Measurements of hygroscopic growth under subsaturated conditions and the CCN activity of aerosols comprised of polyethylene glycol (PEG) with average molecular masses ranging from 200 to 10 000 g mol⁻¹ and mixtures of PEG with ammonium sulfate (AS) were conducted. Experimental results are compared to calculations of hygroscopic growth at thermodynamic equilibrium conducted with the Aerosol Inorganic Organic Mixtures Functional groups Activity Coefficients (AIOMFAC) model, and the potential influence of kinetic limitations on observed water uptake was further explored through estimations of water diffusivity in the PEG oligomers. Particle-phase behavior, including the prevalence of liquid–liquid phase separation (LLPS), was also modeled with AIOMFAC. Under subsaturated relative humidity (RH) conditions, we observed little variability in hygroscopic growth across PEG systems with different molecular masses; however, an increase in CCN activity with increasing PEG molecular mass was observed. This effect is most pronounced for PEG–AS mixtures, and, in fact, an enhancement in CCN activity was observed for the PEG10000–AS mixture as compared to pure AS, as evidenced by a 15 % reduction in critical activation diameter at a supersaturation of 0.8 %. We also observed a marked increase in apparent hygroscopicity for mixtures of higher molecular mass PEG and AS under supersaturated conditions as compared to subsaturated hygroscopic growth. AIOMFAC-based predictions and estimations of water diffusivity in PEG suggest that such discontinuities in apparent hygroscopicity above and below water saturation can be attributed, at least in part, to differences in the sensitivity of water uptake behavior to surface tension effects. There is no evidence that kinetic limitations to water uptake due to the presence of viscous aerosol components influenced hygroscopic growth. For the systems that display an enhancement in apparent hygroscopicity above water saturation, LLPS is predicted to persist to high RH. This indicates a miscibility gap and is likely to influence bulk-to-surface partitioning of PEG at high RH, impacting droplet surface tension and CCN activity. This work provides insight into the factors likely to be contributing to discontinuities in aerosol water-uptake behavior below and above water saturation that have been observed previously in the ambient atmosphere

Gas-particle partitioning of atmospheric aerosols: interplay of physical state, non-ideal mixing and morphology

Atmospheric aerosols, comprising organic compounds and inorganic salts, play a key role in air quality and climate. Mounting evidence exists that these particles frequently exhibit phase separation into predominantly organic and aqueous electrolyte-rich phases. As well, the presence of amorphous semi-solid or glassy particle phases has been established. Using the canonical system of ammonium sulfate mixed with organics from the ozone oxidation of α-pinene, we illustrate theoretically the interplay of physical state, non-ideality, and particle morphology affecting aerosol mass concentration and the characteristic timescale of gas–particle mass transfer. Phase separation can significantly affect overall particle mass and chemical composition. Semi-solid or glassy phases can kinetically inhibit the partitioning of semivolatile components and hygroscopic growth, in contrast to the traditional assumption that organic compounds exist in quasi-instantaneous gas–particle equilibrium. These effects have significant implications for the interpretation of laboratory data and the development of improved atmospheric air quality and climate models

Optimization of process models for determining volatility distribution and viscosity of organic aerosols from isothermal particle evaporation data

The composition of organic aerosol under different ambient conditions as well as their phase state have been a subject of intense study in recent years. One way to study particle properties is to measure the particle size shrinkage in a diluted environment at isothermal conditions. From these measurements it is possible to separate the fraction of low-volatility compounds from high-volatility compounds. In this work, we analyse and evaluate a method for obtaining particle composition and viscosity from measurements using process models coupled with input optimization algorithms. Two optimization methods, the Monte Carlo genetic algorithm and Bayesian inference, are used together with process models describing the dynamics of particle evaporation. The process model optimization scheme in inferring particle composition in a volatility-basis-set sense and composition-dependent particle viscosity is tested with artificially generated data sets and real experimental data. Optimizing model input so that the output matches these data yields a good match for the estimated quantities. Both optimization methods give equally good results when they are used to estimate particle composition to artificially test data. The timescale of the experiments and the initial particle size are found to be important in defining the range of values that can be identified for the properties from the optimization

Kinetic limitations in gas-particle reactions arising from slow diffusion in secondary organic aerosol

The potential for aerosol physical properties, such as phase, morphology and viscosity/diffusivity, to affect particle reactivity remains highly uncertain. We report here a study of the effect of bulk diffusivity of polycyclic aromatic hydrocarbons (PAHs) in secondary organic aerosol (SOA) on the kinetics of the heterogeneous reaction of particle-borne benzo[a]pyrene (BaP) with ozone. The experiments were performed by coating BaP-ammonium sulfate particles with multilayers of SOA formed from ozonolysis of α-pinene, and by subsequently investigating the kinetics of BaP loss via reaction with excess ozone using an aerosol flow tube coupled to an Aerodyne Aerosol Mass Spectrometer (AMS). All reactions exhibit pseudo-first order kinetics and are empirically well described by a Langmuir–Hinshelwood (L-H) mechanism. The results show that under dry conditions (RH 1 × 10^(−12) for dry (RH < 5%), 50% RH and 70% RH conditions, respectively. These results clearly indicate that slow diffusion of reactants through SOA coats under specific conditions can provide shielding from gas-phase oxidants, enabling the long-range atmospheric transport of toxic trace species, such as PAHs and persistent organic pollutants

Size distribution dynamics reveal particle-phase chemistry in organic aerosol formation

Organic aerosols are ubiquitous in the atmosphere and play a central role in climate, air quality, and public health. The aerosol size distribution is key in determining its optical properties and cloud condensation nucleus activity. The dominant portion of organic aerosol is formed through gas-phase oxidation of volatile organic compounds, so-called secondary organic aerosols (SOAs). Typical experimental measurements of SOA formation include total SOA mass and atomic oxygen-to-carbon ratio. These measurements, alone, are generally insufficient to reveal the extent to which condensed-phase reactions occur in conjunction with the multigeneration gas-phase photooxidation. Combining laboratory chamber experiments and kinetic gas-particle modeling for the dodecane SOA system, here we show that the presence of particle-phase chemistry is reflected in the evolution of the SOA size distribution as well as its mass concentration. Particle-phase reactions are predicted to occur mainly at the particle surface, and the reaction products contribute more than half of the SOA mass. Chamber photooxidation with a midexperiment aldehyde injection confirms that heterogeneous reaction of aldehydes with organic hydroperoxides forming peroxyhemiacetals can lead to a large increase in SOA mass. Although experiments need to be conducted with other SOA precursor hydrocarbons, current results demonstrate coupling between particle-phase chemistry and size distribution dynamics in the formation of SOAs, thereby opening up an avenue for analysis of the SOA formation process