Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

ReaxFF-based molecular dynamics study of bio-derived polycyclic alkanes as potential alternative jet fuels

This work investigates the initial stages of the pyrolysis of HtH-1 (C 18 H 32 ; 2,2,7,7,8a,8b-hexamethyl-dodecahydrobiphenylene) and HtH-2 (C 18 H 34 ; 1,1',3,3,3',3'-hexamethyl-1,1'-bi(cyclohexane)), which are bio-derived polycyclic alkanes and potential jet fuels, using ReaxFF force field based molecular dynamics (MD) simulations. Global Arrhenius parameters, such as activation energies and pre-exponential factors, are calculated and used to analyze the overall decomposition kinetics of the fuels. HtH-1 decomposes faster than HtH-2 at the same temperature and density conditions, and they have a faster decomposition rate compared to some existing jet-fuels, such as JP-10. A systematic reaction analysis framework developed in this work is applied to determine a temperature-dependent decomposition mechanism. At lower temperature, the central C-C bond connecting the two cyclohexane rings is dominantly broken in both HtH-1 and HtH-2. However, C-CH 3 bond breaking becomes dominant with increasing temperature due to the large increase in entropy during this reaction. Major products from HtH-1 are C 5 H 8 and C 4 H 8 , and those from HtH-2 are C 4 H 8 and C 2 H 4 . The major products predict that HtH-1 has a higher sooting tendency than HtH-2, which is consistent with measurements. The impact of HtH-2 on the pyrolysis of HtH-1 is also investigated in their binary mixtures. HtH-1 and HtH-2 decompose by unimolecular reactions, and they rarely interact with each other during the pyrolysis of the mixtures. Furthermore, this work demonstrates that ReaxFF can be used to investigate pyrolysis and combustion chemistry of existing or future fuels and to contribute to the development of their chemical kinetic models without any a priori input and chemical intuition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Multi-stage heat release of multi-component fuels: Insights and implications for advanced engine operation

Multi-stage heat release (MSHR) is a unique phenomenon typically seen in lean/diluted and low- to intermediate-temperature combustion. Despite its relevance to advanced engine operation, the MSHR of multi-component fuels has barely been quantified. This study aims to characterize the MSHR of multi-component fuels in a rapid compression machine (RCM) at conditions representative of advanced combustion engines. New experimental data are first reported in the RCM at an equivalence ratio of 0.4, pressure of 60 bar and temperatures from 702 to 795 K for a research grade, multi-component gasoline surrogate, termed PACE-20. Here, experiments confirm the existence of MSHR for PACE-20 at all temperatures, and reveal the strong inhibiting effect of temperature on MSHR, where increasing temperature inhibits MSHR by suppressing first stage heat release and promoting second and third stage heat release. A detailed chemical kinetic model is also adopted to model the experiments, with good agreement observed. The response of MSHR characteristics to changes in different engine operating parameters (i.e., temperature, pressure, equivalence ratio and CO 2 dilution level) and fuel compositions (i.e., mole fraction of n-pentane, n-heptane, isooctane, cyclopentane, 1-hexene, ethanol and aromatics in PACE-20) is further evaluated via statistical analysis coupling quasi-random sampling, extensive computer experiments and high-dimensional model representation. The change in MSHR is characterized through 8 quantities of interest (QOIs), i.e., duration and extent of each heat release stage, and their standard deviations. The analysis highlights the dependence of the QOIs on the individual parameters as well as their interplays, with temperature exhibiting the strongest impact among all parameters. Furthermore, it is demonstrated that if MSHR is to be facilitated and combustion phasing is to be extended to enable engine operation at higher compression ratios, it is recommended to operate the engine at low intake temperatures, boosted intake pressures and high CO 2 dilutions (i.e., high EGR levels) with gasoline fuels containing more n-pentane and less aromatics.

42 ENGINEERING↗
DOE OSTI2022

Water film-driven Mn (oxy)(hydr)oxide nanocoating growth on rhodochrosite

Minerals exposed to moist air stabilize thin water films that drive a score of chemical reactions of great importance to water-unsaturated terrestrial environments. In this study, we identified Mn (oxy)(hydr)oxide nanocoatings formed by the dissolution, oxidation and precipitation of Mn in oxygenated water films grown on rhodochrosite (MnCO 3 ) microparticles. Nanocoatings that could be identified by vibrational spectroscopy, X-ray diffraction, X-ray photoelectron spectroscopy, and (scanning and transmission) electron microscopy formed in water films containing the equivalent of at least 7 monolayers (~84 H 2 O/nm 2 ). These films were formed by exposing microparticles to moist air with at least 50% relative humidity (RH). Films of neutral pH reacted up to 14% of the Mn II located in the topmost ~5 nm region of the microparticles in atmospheres of up to 90% RH for 7 d. These reactions produced MnOOH, birnessite (MnO 2 ) and hausmannite (Mn 3 O 4 ) nanoparticles of low crystallinity, while exposure to atmospheric air for 1 yr. converted only 2% of Mn II in this region to MnOOH. In contrast, reactions in alkaline water films converted up to ~75% of the Mn II but only after 16 d of reaction. These films produced MnOOH and MnO 2 of low crystallinity, as well as crystalline hausmannite. Kinetic modeling of the time-resolved growth of the Mnsingle bondO stretching vibrational bands of these nanocoatings revealed two concurrent reaction processes. A 1 rst -order process was assigned to nucleation events terminating only after a few hours, and a 0-order process was assigned to the sustained growth of nanocoatings from these nuclei over longer reaction time. By identifying nanocoatings formed by water film-driven reactions on rhodochrosite, our study adds new insight into mineralogical transformations relevant to anoxic–oxic boundaries in water-unsaturated environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Microstructure-based prediction of thermal aging strength reduction factors for grade 91 ferritic-martensitic steel

This study aimed to develop a microstructure-based mechanistic approach to address the long-term thermal aging effect on yield stress and ultimate tensile strength and to provide a physical basis for developing thermal aging factors for G91 for a design life of 60 years. Several heats of G91 steel were examined. Controlled aging experiments were conducted on two heats at temperatures of 550, 600, and 650 °C for times up to ~64,000 h. Specimens were also taken from archived creep-tested specimens of G91 and from a tube removed from the Kingston coal-fired power plant to obtain data that cover a wide range of temperatures and for aging times up to 155,000 h. Thermal aging caused significant subgrain recovery, coarsening of M 23 C 6 carbides at subgrain boundaries and MX carbonitrides within subgrains. The intermetallic Laves phase forms during aging and grows rapidly. The growth rate of the subgrain width and MX mean size during thermal aging were described by kinetic models. Thermal aging results in the reductions in the yield stress and the ultimate tensile strength of G91. The effects of thermal aging on the reductions in yield stress and ultimate tensile strength were well described by the microstructure-strength model that considers three superimposed strengthening mechanisms, namely, sub-boundary strengthening, MX precipitation hardening, and Mo solid solution strengthening. Furthermore, the model was independently validated by the ASME Code values and is being used by the ASME to develop the aging-induced strength reduction factors for G91 steel for a design life of 60 years.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Methanol oxidation up to 100 atm in a supercritical pressure jet-stirred reactor

Methanol (CH 3 OH) has attracted considerable attention as a renewable fuel or fuel additive with low greenhouse gas emissions. Methanol oxidation was studied using a recently developed supercritical pressure jet-stirred reactor (SP-JSR) at pressures of 10 and 100 atm, at temperatures from 550 to 950 K, and at equivalence ratios of 0.1, 1.0, and 9.0 in experiments and simulations. Here, the experimental results show that the onset temperature of CH 3 OH oxidation at 100 atm is around 700 K, which is more than 100 K lower than the onset at 10 atm and this trend cannot be predicted by the existing kinetics models. Furthermore, a negative temperature coefficient (NTC) behavior was clearly observed at 100 atm at fuel rich conditions for methanol for the first time. To understand the observed temperature shift in the reactivity and the NTC effect, we updated some key elementary reaction rates of relevance to high pressure CH 3 OH oxidation from the literature and added some new low-temperature reaction pathways such as CH 2 O + HO 2 = HOCH 2 O 2 (RO 2 ), RO 2 + RO 2 = HOCH 2 O (RO) + HOCH 2 O (RO) + O 2 , and CH 3 OH + RO 2 = CH 2 OH + HOCH 2 O 2 H (ROOH). Although the model with these updates improves the prediction somewhat for the experimental data at 100 atm and reproduces well high-temperature ignition delay times and laminar flame speed data in the literature, discrepancies still exist for some aspects of the 100 atm low-temperature oxidation data. In addition, it was found that the pressure-dependent HO 2 chemistry shifts to lower temperature as the pressure increases such that the NTC effect at fuel-lean conditions is suppressed. Therefore, as shown in the experiments, the NTC phenomenon was only observed at the fuel-rich condition where fuel radicals are abundant and the HO 2 chemistry at high pressure is weakened by the lack of oxygen resulting in comparatively little HO 2 formation.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Shock tube/laser absorption measurement of the rate constant of the reaction: H 2 O 2 + CO 2 2OH + CO 2

We address the role of the linear mixing rule in the kinetics of the H 2 O 2 decomposition system by reporting the rate constant for H 2 O 2 + M = 2OH + M (M = Ar and CO 2 ) in the temperature range of 1087-1234 K at low pressures in a mixture of 20% CO 2 in Argon. The reaction rate constant was inferred from H 2 O concentrations monitored by using a laser-absorption spectroscopy-based water diagnostic. To the best of our knowledge, this is the first measurement of the rate constant of this reaction in a mixture to be reported in literature. A significant discrepancy was found between the rate constants derived using the traditional linear mixing rule and the reduced pressure linear mixing rule. Furthermore, this discrepancy can have serious implications on the predictive accuracy of these kinetic models, especially under conditions relevant to the operation of supercritical CO 2 (sCO 2 ) power cycles that rely on oxy-fuel combustion in a working fluid comprised almost entirely of CO 2 .

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Enhancing $\phi$-sensitivity of ignition delay times through dilution of fuel-air mixture

The high $\phi$-sensitivity (η) of ignition delay time (τ IDT ) is one of the desirable fuel properties for the high-load extension of advanced compression ignition engine. Recent studies revealed the effectiveness of a high dilution rate (x D ) for enhancing η at the engine-relevant conditions. This study aims to quantify the effect of dilution on the η of isooctane using a combined experiment-simulation approach. The ignition delay of the isooctane/air mixture was measured with an Advanced Fuel Ignition Delay Analyzer (AFIDA) over the temperature range of 623 - 923 K at 10 bar pressure with global $\phi$ = 0.3 - 0.6, with and without 28.6% of additional N 2 dilution. For precise evaluation of experimental η, the facility effect of the AFIDA experiment was characterized with three-dimensional computational fluid dynamics (3-D CFD) simulation. The temperature in the combustion chamber from 3-D CFD indicated a substantial temporal dependency, varying up to ~52 K by charge-cooling of fuel injection and heat transfer from the wall. We introduced the dimensionless number θ(t) for characterizing the temporal profile of chamber temperature. Consideration of the facility effect using θ(t) resulted in better agreement between the experimental ..eta.. and zero-dimensional (0-D) kinetics simulation. The refined η were then further utilized to quantify the effectiveness of dilution to η The extent of η enhancement with dilution strategy was maximized at the low-temperature chemistry regime, increasing η by 77% with a 28.6% dilution rate. Further analysis on the dilution effect was carried out using 0-D kinetics simulation, revealing the critical dimensionless numbers relevant to the effectiveness of dilution to η enhancement. Here, this study is the first experiment-simulation combined research to quantify the effect of dilution on η, facilitating the kinetics model refinement for better reproduction of $\phi$-sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Radon kinetics in a natural indoor radon chamber

An unusual 180 m 3 storage room in the basement of a two-story laboratory building is unventilated, and separated from occupiable rooms by double steel doors. The space completely borders on soil through the concrete floor and two of its concrete walls. The room also contains a separate inner chamber with 1 m thick concrete walls designed to damp vibrations in the room above it. The space boasts a relatively high radon level, 1083 Bq m –3 , which varies with local outdoor environmental conditions. Measurements were made of radon concentrations at various locations and heights within the facility. More than a year of continuous radon concentration data corresponding to a single location are also available, along with measurements of indoor and outdoor pressure, temperature, and humidity. Data were also collected with as many as five fans placed in different locations and cycled on for variable time periods. First order linear kinetic models were created to explain the observed approaches to steady state due to changing conditions and wash-out resulting from intentional ventilation. Results demonstrate a good fit between changes in the radon concentration level and the developed compartmental models. Furthermore, no significant differences were observed between radon concentration at different locations or heights in the chamber.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Influence of Polymer Type and Carbon Nanotube Properties on Carbon Nanotube/Polymer Nanocomposite Biodegradation

The interaction of anaerobic microorganisms with carbon nanotube/polymer nanocomposites (CNT/PNC) will play a major role in determining their persistence and environmental fate at the end of consumer use when these nano-enabled materials enter landfills and encounter wastewater. Motivated by the need to understand how different parameters (i.e., polymer type, microbial phenotype, CNT characteristics) influence CNT/PNC biodegradation rates, we have used volumetric biogas measurements and kinetic modeling to study biodegradation as a function of polymer type and CNT properties. In one set of experiments, oxidized multiwall carbon nanotubes (O-MWCNTs) with a range of CNT loadings 0 – 5% w/w were incorporated into poly-e-caprolactone (PCL) and polyhydroxyalkanoates (PHA) matrices and subjected to biodegradation by an anaerobic microbial community. For each CNT/PNC, complete polymer biodegradation was ultimately observed, although the rate of biodegradation was inhibited above certain critical CNT loadings dependent upon the polymer type. Higher loadings of pristine MWCNTs were needed to decrease the rate of polymer biodegradation compared to O-MWCNTs, an effect ascribed principally to differences in CNT dispersion within the polymer matrices. Above ertain CNT loadings, a CNT mat of similar shape to the initial PNC was formed after polymer biodegradation, while below this threshold, CNT aggregates fragmented in the media. In situations where biodegradation was rapid, methanogen growth was disproportionately inhibited compared to the overall microbial community. Analysis of the results obtained from this study indicates that the inhibitory effect of CNTs on polymer biodegradation rate is greatest under conditions (i.e., polymer type, microbial phenotype, CNT dispersion) where biodegradation of the neat polymer is slowest. This new insight provides a means to predict the environmental fate, persistence and transformations of CNT-enabled polymer materials.

Frank, Benjamin↗
DOE OSTI2021

Ethylene oxidation on unpromoted silver catalysts: Reaction pathway and selectivity analysis using DFT calculations

Chemical conversions in catalytic partial oxidation processes of light hydrocarbons are responsible for the production of numerous industrial chemicals, plastics, and intermediates. These processes are relatively expensive to perform, and are typically operated at high thermodynamic inefficiency, so the development of novel, highly efficient catalysts would prove to be very cost effective. Herein, our study focused on surface catalytic mechanisms of the ethylene oxide (EO) formation process. Periodic plane-wave Density Functional Theory (DFT) methods were used to analyze related reaction mechanisms on the Ag(111) surface facet with low coverage. Energetic changes of related species and pathways were calculated. Key surface species are identified to suggest the factors for the observed selectivity during EO formation. Our results are consistent with previous kinetic modeling efforts in the literature which did not employ DFT analysis. Lastly, our study demonstrates how fundamental theoretical investigations and multi-scale modeling techniques are currently impacting the advancement of rational catalyst design and microkinetic modeling techniques in the light hydrocarbon processing industry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

A novel sub-pilot-scale direct-contact ultrasonic dehydration technology for sustainable production of distillers dried grains (DDG)

DDG is a major source of protein, calcium, phosphorus, and sulfur is arguably the most important byproduct of the bioethanol industry with increasing demand over the past few years. Reducing energy consumption in the DDG production process and energy recovery from DDG is vital for sustainable bioethanol productions. In this paper, a novel direct-contact multi-frequency, multimode, and modulated (MMM) ultrasonic dryer (US) was developed for the first time and has been applied in dehydration of wet distillers' grain (WDG). Ultrasonic drying (US) was combined with a convective airflow (HA) at different temperatures of 25 (room temperature), 50 and 70 °C to evaluate the impact of US, HA, and US + HA on drying kinetics, activation energy, chemical compositions, microstructure, and color of DDG. Semi-empirical kinetic models were developed and evaluating drying performances showed that the application of ultrasound significantly enhanced the drying rate and decreased the drying time (by 46%), especially at low drying temperatures. The activation energy for moisture removal in the presence of ultrasound was about 50% of that without ultrasound. The final dried distillers' grains product processed by ultrasonic drying had a brighter color, a higher available protein, a higher digestible protein (the lowest acid detergent insoluble crude protein), and a better surface profile with no compromise on minerals and fiber contents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Anhydrous volatile fatty acid extraction through omniphobic membranes by hydrophobic deep eutectic solvents: Mechanistic understanding and future perspective

Volatile fatty acids (VFAs) derived from arrested anaerobic digestion (AD) can be recovered as a valuable commodity for value-added synthesis. However, separating VFAs from digestate with complex constituents and a high-water content is an energy-prohibitive process. This study developed an innovative technology to overcome this barrier by integrating deep eutectic solvents (DESs) with an omniphobic membrane into a membrane contactor for efficient extraction of anhydrous VFAs with low energy consumption. Here, a kinetic model was developed to elucidate the mechanistic differences between this novel omniphobic membrane-enabled DES extraction and the previous hydrophobic membrane-enabled NaOH extraction. Experimental results and mechanistic modeling suggested that VFA extraction by the DES is a reversible adsorption process facilitating subsequent VFA separation via anhydrous distillation. High vapor pressure of shorter-chain VFAs and low Nernst distribution coefficients of longer-chain VFAs contributed to DES-driven extraction, which could enable continuous and in-situ recovery and conversion of VFAs from AD streams.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Reconfigurable structural color by reversible switching of colloidal discoid liquid crystal alignment

Combining the effects of alternating-current (AC) electric fields and sedimentation reconfigures colloidal discoids between planar and homeotropic orientational alignments, thereby controlling the material’s structural color. Here, we self-assemble micrometer-size polystyrene discoids in an isopropanol-water mixture. After sedimentation, the discoids adopt a homeotropic alignment, with the minor axis perpendicular to the substrate. Adding an AC electric field (1 kHz) of ≥0.50 V switches the discoids to planar alignment—with the minor axis parallel to the substrate—within ∼100 s. Removing the field switches them back to homeotropic alignment within ∼300 s. Kinetic modeling of the field-induced torques yields good agreement with these measurements. The peak wavelength and intensity of the crystal’s diffraction response shift significantly upon reconfiguration; scattering simulation predicts these shifts. The reconfiguration is maintained for at least 10 cycles. This method is a simple, scalable avenue to reconfigure the optical properties of colloidal crystals produced from simple dielectric spheroids.

colloidal crystals↗
DOE OSTI2026

Spectrally accurate, reverse-mode differentiable bounce-averaging algorithm and its applications

We present a fast, spectrally (exponentially) accurate, automatically differentiable bounce-averaging algorithm that is used to simplify kinetic models. Using this algorithm, implemented in the DESC stellarator optimisation suite, we can perform efficient optimisation of many objectives to improve stellarator performance, such as the effective ripple 𝜖 eff metric for the neoclassical transport coefficient in the low collisionality regime and proxies for energetic particle confinement. For the first time, we optimise a finite-beta stellarator to directly reduce neoclassical ripple transport using reverse-mode differentiation. This ensures the computational cost of differentiation is independent of the number of controllable parameters.

fusion plasma↗
DOE OSTI2022

Prospects and Limitations of Predicting Fuel Ignition Properties from Low-Temperature Speciation Data

Using chemical kinetic modeling and statistical analysis, we investigate the possibility of correlating key chemical “markers”–typically small molecules–formed during very lean (φ ~ 0.001) oxidation experiments with near-stoichiometric (φ ~ 1) fuel ignition properties. One goal of this work is to evaluate the feasibility of designing a fuel-screening platform, based on small laboratory reactors that operate at low temperatures and use minimal fuel volume. Buras et al. [Combust. Flame2020,216, 472–484] have shown that convolutional neural net (CNN) fitting can be used to correlate first-stage ignition delay times (IDTs) with OH/HO2 measurements during very lean oxidation in low-T flow reactors with better than factor-of-2 accuracy. In this work, we test the limits of applying this correlation-based approach to predict the low-temperature heat release (LTHR) and total IDT, including the sensitivity of total IDT to the equivalence ratio, φ. We demonstrate that first-stage IDT can be reliably correlated with very lean oxidation measurements using compressed sensing (CS), which is simpler to implement than CNN fitting. LTHR can also be predicted via CS analysis, although the correlation quality is somewhat lower than for first-stage IDT. In contrast, the accuracy of total IDT prediction at φ = 1 is significantly lower (within a factor of 4 or worse). Furthermore, these results can be rationalized by the fact that the first-stage IDT and LTHR are primarily determined by low-temperature chemistry, whereas total IDT depends on low-, intermediate-, and high-temperature chemistry. Oxidation reactions are most important at low temperatures, and therefore, measurements of universal molecular markers of oxidation do not capture the full chemical complexity required to accurately predict the total IDT even at a single equivalence ratio. As a result, we find that φ-sensitivity of ignition delay cannot be predicted at all using solely correlation with lean low-T chemical speciation measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Coupling of Lipid Peroxidation and Criegee Intermediate Mediated Autoxidation in the Heterogeneous Oxidation of Linoleic Acid Aerosols

Autoxidation is an established mechanism for the degradation of organic molecules, which is relevant in the atmosphere, combustion processes, the environment, and the rancidification of lipids (commonly known as lipid peroxidation). Autoxidation proceeds via radical chain reactions involving hydroxyl (•OH), peroxy (RO2•), and alkoxy radicals, which are also prominent oxidants in the atmosphere. Recent reports have provided evidence for an alternative autoxidation mechanism driven instead by Criegee intermediates (CIs), which are produced from the reaction of β-hydroxy peroxy radicals (β-OH-RO2•). This work evaluates the contributions of these two mechanisms in the •OH initiated heterogeneous oxidation of linoleic acid (LA) aerosols. Reaction kinetics and product distributions are monitored using a vacuum ultraviolet photoionization aerosol time-of-flight mass spectrometer. To explain the observed kinetics, a kinetic model is developed that incorporates both the conventional peroxidation and alternative CI-mediated autoxidation mechanisms. We observe that the CI-mediated autoxidation pathways enhance the heterogeneous autoxidation rate, while the peroxidation reactions, although present, contributes less to the overall oxidation rate. α-Acyloxyalkyl hydroperoxides (AAHPs) are identified as key indicators for bimolecular reactions of CI with LA, highlighting the role of LA as a CI scavenger. Moreover, the measured functionalized LA products with hydroxyl or carbonyl group(s), serve as markers for the peroxidation reactions. In summary, this work presents a quantitative framework to understand the coupled reaction network of •OH, RO2•, β-OH-RO2• radicals, and CI in driving heterogeneous autoxidation, which is crucial for understanding degradation mechanisms of organic molecules in the environment and atmosphere.

Criegee intermediate↗
DOE OSTI2026

Dramatic Slowdown of Bimolecular Reactions of Criegee Intermediates in Organic Aerosols

Although Criegee intermediates (CIs) play a central role in driving the transformation of organic molecules, their condensed-phase reaction kinetics remains poorly understood. While CIs undergo unimolecular decomposition to produce OH radicals, stabilized CIs participate in a host of other bimolecular reactions forming larger molecular weight products, including secondary ozonides (kSOZ), alkoxyalkyl hydroperoxides (kAlAHP), and acyloxyalkyl hydroperoxides (kAcAHP). Here, we demonstrate, using prior literature and new measurements of heterogeneous organic aerosol reactions using O3 or OH radicals, that these condensed-phase bimolecular reaction rate coefficients kSOZ (10-20 cm3·molecule-1·s-1), kAlAHP, and kAcAHP (10-19 cm3·molecule-1·s-1) are 7-10 orders of magnitude slower than analogous gas-phase reactions (∼ 10-10-10-13 cm3·molecule-1·s-1). This contrasts with current evidence that the rate coefficients for unimolecular steps in the condensed phase are similar in magnitude to those in the gas phase. Kinetic models using gas-phase-like bimolecular coefficients consistently overestimate the importance of bimolecular relative to unimolecular pathways while underpredicting the aerosol reactivity by up to 30-fold. The slowing of bimolecular reactivity in aerosols is likely due to solvation effects and the increased steric hindrance of larger molecular weight CIs and their scavengers in the condensed phase. These results challenge the conventional view that the high reactant density of an organic aerosol should favor bimolecular over unimolecular reactions, underscoring the need for more accurate CI condensed-phase rate constants for reliable modeling of oxidative processes across a diverse set of atmospheric and environmental systems.

Zeng, Meirong↗