1,350
Views
2
CrossRef citations to date
0
Altmetric
Original Articles

Effects of volatility, viscosity, and non-ideality on particle–particle mixing timescales of secondary organic aerosols

ORCID Icon, ORCID Icon & ORCID Icon
Pages 411-426 | Received 30 Jul 2023, Accepted 30 Aug 2023, Published online: 26 Sep 2023

Figures & data

Figure 1. Schematic of the two populations implementation in KM-GAP. Both populations consist of a near-surface bulk layer (nsb) and the number of bulk layers (b1-L; L = 10 in this study). Additionally, each particle has a sorption layer containing adsorbed and desorbing species and a near-surface gas phase to account for potential gas phase diffusion, which is not pictured. There is one shared gas phase between the two populations. Each population consists of different nonvolatile species, while population 1 initially contains semi-volatile species in equilibrium with the gas phase.

Figure 1. Schematic of the two populations implementation in KM-GAP. Both populations consist of a near-surface bulk layer (nsb) and the number of bulk layers (b1-L; L = 10 in this study). Additionally, each particle has a sorption layer containing adsorbed and desorbing species and a near-surface gas phase to account for potential gas phase diffusion, which is not pictured. There is one shared gas phase between the two populations. Each population consists of different nonvolatile species, while population 1 initially contains semi-volatile species in equilibrium with the gas phase.

Figure 2. Particle–particle mixing timescale (τmix) of two populations where population 1 is liquid (Db1 = 10−8 cm2 s−1) and population 2 is (a) liquid (Db2 = 10−8 cm2 s−1), (b) semi-solid (Db2 = 10−15 cm2 s−1), and (c) ultraviscous semi-solid (Db2 = 10−18 cm2 s−1).

Figure 2. Particle–particle mixing timescale (τmix) of two populations where population 1 is liquid (Db1 = 10−8 cm2 s−1) and population 2 is (a) liquid (Db2 = 10−8 cm2 s−1), (b) semi-solid (Db2 = 10−15 cm2 s−1), and (c) ultraviscous semi-solid (Db2 = 10−18 cm2 s−1).

Figure 3. The final equilibrium number concentration of the semi-volatile species in population 2 with different volatilities (C0) and activity coefficients of semi-volatile species in population 2 (γ2).

Figure 3. The final equilibrium number concentration of the semi-volatile species in population 2 with different volatilities (C0) and activity coefficients of semi-volatile species in population 2 (γ2).

Figure 4. (a) Particle–particle mixing timescale (τmix) with various diffusivities and activity coefficients of semi-volatile species (C0 = 10 µg m−3) in population 2 (γ2). (b) τmix when the bulk diffusivity (Db) of a semi-volatile species with C0 = 10 µg m−3 is varied in both populations. In these simulations, the semi-volatile species is ideally miscible in both populations.

Figure 4. (a) Particle–particle mixing timescale (τmix) with various diffusivities and activity coefficients of semi-volatile species (C0 = 10 µg m−3) in population 2 (γ2). (b) τmix when the bulk diffusivity (Db) of a semi-volatile species with C0 = 10 µg m−3 is varied in both populations. In these simulations, the semi-volatile species is ideally miscible in both populations.

Figure 5. Particle–particle mixing timescale (τmix) when mass loading of population 2 (COA2) is varied in the range of 0.1–100 µg m−3, while mass loading of population 2 (COA1) is held constant at 40 µg m−3.

Figure 5. Particle–particle mixing timescale (τmix) when mass loading of population 2 (COA2) is varied in the range of 0.1–100 µg m−3, while mass loading of population 2 (COA1) is held constant at 40 µg m−3.

Figure 6. Simulations of toluene SOA mixing experiments and the data from Ye et al. (Citation2016, Citation2018). Plotted is the mass fraction of SVOCs evaporated from population 1 in population 2. (a) Mass fraction of toluene SOA (population 1) in D-toluene SOA (population 2). RH legend indicates the RH at the end of the experiment. For each experiment, the RH begins at 10%. For RH 20%, the RH is increased at 1 h; for RH 40%, it is increased at 2 h; and for RH 90%, the experimental RH was increased continually, reaching 90% at 1 h, while in the model this is achieved by a stepwise increase in RH. (b) Mass fraction of toluene SOA (population 1) in limonene SOA (population 2) initially at 10% RH with an increase to 30% at 2 h. (c) Mass fraction of toluene SOA (population 1) in β-caryophyllene SOA (population 2) initially at 10% RH with an increase to 30% at 2 h.

Figure 6. Simulations of toluene SOA mixing experiments and the data from Ye et al. (Citation2016, Citation2018). Plotted is the mass fraction of SVOCs evaporated from population 1 in population 2. (a) Mass fraction of toluene SOA (population 1) in D-toluene SOA (population 2). RH legend indicates the RH at the end of the experiment. For each experiment, the RH begins at 10%. For RH 20%, the RH is increased at 1 h; for RH 40%, it is increased at 2 h; and for RH 90%, the experimental RH was increased continually, reaching 90% at 1 h, while in the model this is achieved by a stepwise increase in RH. (b) Mass fraction of toluene SOA (population 1) in limonene SOA (population 2) initially at 10% RH with an increase to 30% at 2 h. (c) Mass fraction of toluene SOA (population 1) in β-caryophyllene SOA (population 2) initially at 10% RH with an increase to 30% at 2 h.