Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion coefficient”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗
DOE OSTI2020

Residually Stressed Bimaterial Beam Specimen for Measuring Environmentally Assisted Crack Growth

Background: Subcritical crack growth can occur in a brittle material when the stress intensity factor is smaller than the fracture toughness if an oxidizing agent (such as water) is present at the crack tip. Objective: Here we present a novel bi-material beam specimen which can measure environmentally assisted crack growth rates. The specimen is “self-loaded” by residual stress and requires no external loading. Methods: Two materials with different coefficient of thermal expansion are diffusion bonded at high temperature. After cooling to room temperature a subcritical crack is driven by thermal residual stresses. A finite element model is used to design the specimen geometry in terms of material properties in order to achieve the desired crack tip driving force. Results: The specimen is designed so that the crack driving force decreases as the crack extends, thus enabling the measurement of the crack velocity versus driving force relationship with a single test. The method is demonstrated by measuring slow crack growth data in soda lime silicate glass and validated by comparison to previously published data. Conclusions: The self-loaded nature of the specimen makes it ideal for measuring the very low crack velocities needed to predict brittle failure at long lifetimes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Machine learning-based ethylene and carbon monoxide estimation, real-time optimization, and multivariable feedback control of an experimental electrochemical reactor

Electrochemical reduction of CO 2 gas is a novel CO 2 utilization technique that has the potential to mitigate the global climate crisis caused by anthropogenic CO 2 emissions, and enable the large-scale storage of energy generated from renewable sources in the form of carbon-based chemicals and fuels. However, due to the complexity of the electrochemical reactions, the explicit first-principles models for CO2 reduction are not available yet, and there has been a limited effort to develop process modeling, optimization and control of CO 2 electrochemical reactors. To this end, a rotating cylinder electrode (RCE) reactor has been constructed at UCLA to understand the mass transfer and reaction kinetics effects separately on the productivity. In the RCE reactor, the applied potential strongly influences the reaction energetics and the electrode rotation speed affects the hydrodynamic boundary layer and modifies the film mass transfer coefficient, which involves convective and diffusive transport. Further, the present work aims to develop a multi-input multi-output (MIMO) control scheme for the RCE reactor that integrates techniques from artificial and recurrent neural network modeling, nonlinear optimization, and process controller design. Specifically, production rates of two products from the experimental reactor, ethylene and carbon monoxide, are controlled by manipulating two inputs, applied potential and catalyst rotation speed. Process dynamics and controllability are analyzed, a feedback control strategy is designed and the controllers are tuned accordingly. The experimental electrochemical cell is employed to gather data for process modeling and implement the multivariable control system. Finally, the experimental results are presented which demonstrate excellent closed-loop performance by the control system and regulation of the outputs at three different set-points including an economically-optimal set-point.

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

Quantifying heavy quark transport coefficients with an improved transport model

We report that the heavy-flavor transport coefficients contain important information on the strong interaction at finite temperatures. The extraction of these numbers from experimental data requires dynamical modeling of heavy-flavor transport that is coupled to realistic medium evolution. Furthermore, meaningful extractions necessitate both a faithful implementation of the physical inputs to be tested and the quantification of model uncertainty. For these purposes, we have developed a partonic transport model LIDO. It has an improved treatment of in-medium parton bremsstrahlung, which has been calibrated to theoretical calculations in a simple medium to reduce modeling uncertainty. Regarding the interaction between heavy quark and the medium, few-body perturbative scatterings are applied to large-momentum transfer (q) processes, while a diffusion equation models the dynamics of small-q processes. Such a separation restricts the explicit use of medium quasi-particles to large-q processes only. Another advantage is that deviations from the leading-order probe-medium coupling can be parametrized as an additional contribution to the diffusion constant. The heavy quark transport coefficients are then extracted with uncertainty estimation from a Bayesian analysis including both the RHIC and the LHC data. The results are found to be consistent with earlier extraction of the light-quark transport coefficients at high momentum and be comparable with lattice calculations of the heavy-flavor diffusion constant in the static limit at low momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Leveraging Sequential Doping of Semiconducting Polymers to Enable Functionally Graded Materials for Organic Thermoelectrics

With the ability to modulate electronic properties through molecular doping coupled with ease in processability, semiconducting polymers are at the forefront in enabling organic thermoelectric devices for thermal energy management. In contrast to uniform thermoelectric material properties, an alternative route focuses on functionally graded materials (FGMs) where one spatially controls and optimizes transport properties across the length of a thermoelectric material. While primarily studied in the context of inorganic materials, the concept of FGMs for organic thermoelectrics has not been explored. Herein, we introduce how molecular doping of semiconducting polymers enables spatial compositional control of thin-film FGMs. Specifically, we utilize sequential vapor doping of poly[2,5-bis(3-tetradecylthiophen-2-yl) thieno [3,2-b]thiophene] (PBTTT) with the small molecule acceptor 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ) to fabricate the simplest form of FGMs – double- segmented thin films. The two thin-film segments are of equal length (7.5 mm), but each set to different doping levels. Our study focuses on understanding the thermoelectric properties (Seebeck coefficient, α, and electronic conductivity, σ) and structural properties (through X-ray scattering, UV-vis-NIR spectroscopy, and Raman spectroscopy) within and across the two segments. We observe the presence of a small diffuse interfacial region of 0.5 – 1 mm between the two segments where the doping level and transport properties vary continuously. Despite the diffuse interface, the measured effective Seebeck coefficient (αeff) across the two segments is simply the average of α within each segment. Importantly, this experimental result is consistent with reported mathematical models describing αeff in graded thermoelectric materials. Finally, our results demonstrate the facile fabrication and characterization of functionally graded organic thermoelectric materials, providing guidelines for further development on more complex FGMs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles investigation of lanthanides diffusion in HCP zirconium via vacancy-mediated transport

The diffusion of lanthanide fission products plays an important role in the growth of the fuel-cladding chemical interaction (FCCI) region in metallic fuels. The use of a Zr interdiffusion barrier may mitigate the transport of lanthanides from the fuel to the cladding, but the efficacy of such a liner is not yet known. In this paper, the stability and vacancy-mediated diffusion of La, Ce, Pr, and Nd in hexagonal close-packed (HCP) Zr is investigated via density functional theory (DFT) calculations and self-consistent mean field (SCMF) analysis. DFT is used to calculate the formation, binding, and migration energies of vacancies and vacancy-solute pairs. The DFT energetics are used in the KineCluE code to calculate the Onsager transport coefficients. La is found to be the fastest diffusing species in HCP Zr and experiences an almost isotropic diffusion behavior. The other three species (Ce, Pr, and Nd) demonstrate anisotropic diffusion where the diffusion in the basal planes is significantly faster than that along the c-axis. Further, the calculated lanthanide diffusivities in HCP Zr are fitted to an Arrhenius relation and the activation energies and prefactors are reported for the first time. Furthermore, the vacancy drag and the segregation tendencies were analyzed using the calculated off-diagonal transport coefficients. According to our vacancy-mediated diffusion model, lanthanides are expected to be enriched at vacancy sinks at low temperatures, while at high temperatures, lanthanides are depleted at sinks and will preferably diffuse into the bulk. The enrichment/depletion transition temperature depends on the diffusion direction (basal or axial) and hence will be controlled by the grain texture and orientation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Mechanical and Thermophysical Properties of 3D-Printed SiC-FY20

In the Transformational Challenge Reactor (TCR), the fuel blocks consist of an additively-manufactured silicon carbide (SiC) matrix and uranium nitride tristructural isotropic (UN TRISO) fuel particles, which are stacked to form fuel columns. The SiC matrix is manufactured using binderjet 3D printing followed by loading the TRISO fuel particles and the chemical vapor infiltration (CVI) process. Because the fuel matrix is a primary component of the TCR core and its response to mechanical and thermal loads during operation is one of the most influential factors on the integrity of TCR core, testing and evaluation have focused on producing mechanical and thermophysical properties data for the binderjet/CVI SiC. Mechanical and thermophysical properties were measured from various types of specimens printed for two or three orientations, which included equibiaxial flexural failure strength, elastic constants, thermal diffusivity and conductivity, density, and the coefficient of thermal expansion. Flexural failure strength datasets showed similar Weibull distributions regardless of sample variants including different orientations. The mean failure strengths of the 3D-printed SiC variants were in the range of 286–306 MPa, which are 22–27% lower than that of the CVD SiC. Thermophysical test results showed that specific heat and thermal expansion are not sensitive to the build directions of SiC samples, while thermal diffusivity is highly dependent on the build direction and can be correlated to the anisotropic character of the 3D-printed SiC. This report also includes discussions on the uniaxial tensile properties of the as-printed SiC before CVI and on ongoing efforts for irradiation effects studies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Conditional Karhunen-Loève expansion for uncertainty quantification and active learning in partial differential equation models

We use a conditional Karhunen-Lo` eve (KL) model to quantify and reduce uncertainty in a stochastic partial differential equation (SPDE) problem with partially-known space-dependent coefficient, Y (x). We assume that a small number of Y (x) measurements are available and model Y (x) with a KL expansion. We achieve reduction in uncertainty by conditioning the KL expansion coefficients on measurements. We consider two approaches for conditioning the KL expansion: In Approach 1, we condition the KL model first and then truncate it. In Approach 2, we first truncate the KL expansion and then condition it. We employ the conditional KL expansion together with Monte Carlo and sparse grid collocation methods to compute the moments of the solution of the SPDE problem. Uncertainty of the problem is further reduced by adaptively selecting additional observation locations using two active learning methods. Method 1 minimizes the variance of the PDE coefficient, while Method 2 minimizes the variance of the solution of the PDE. We demonstrate that conditioning leads to dimension reduction of the KL representation of Y (x). For a linear diffusion SPDE with uncertain log-normal coefficient, we show that Approach 1 provides a more accurate approximation of the conditional log-normal coefficient and solution of the SPDE than Approach 2 for the same number of random dimensions in a conditional KL expansion. Furthermore, Approach 2 provides a good estimate for the number of terms of the truncated KL expansion of the conditional field of Approach 1. Finally, we demonstrate that active learning based on Method 2 is more efficient for uncertainty reduction in the SPDE’s states (i.e., it leads to a larger reduction of the variance) than active learning using Method 2.

Conditioned Karhunen-Lo` eve expanion, machine lea↗
DOE OSTI2023

High‐Throughput Evaluation of Hardening Coefficients of Eight Alloying Elements in Magnesium

Liquid–solid diffusion couples (LSDCs) are employed to generate a composition gradient in the single‐phase hexagonal closed‐packed (hcp) solid solution with compositions up to the solubility limit of various solutes in Mg. Nanoindentation scanning across the composition gradient in LSDCs allows effective evaluation of composition‐dependent hardness of eight alloying elements (Al, Ca, Ce, Gd, Li, Sn, Y, and Zn) in the hcp Mg phase. The hardening coefficients, an indicator of the potency of solid‐solution hardening, are evaluated from the measured composition‐hardness data and correlated with various materials properties such as atomic radius, shear modulus, and elastic modulus of the solutes. The rank of hardening potency of Al, Gd, Sn, Y, and Zn measured by nanoindentation is in good agreement with that measured by microindentation reported in the literature. The hardening coefficient (potency) from the strongest to the weakest is Ce > Ca > Y ≈ Gd > Zn > Al ≈ Sn > Li in Mg‐based hcp binary solid solutions. The hardening coefficient is found to be closely correlated with the strengthening potency.

36 MATERIALS SCIENCE↗
DOE OSTI2022

New Particle Formation and Growth Dynamics for α-Pinene Ozonolysis in a Smog Chamber and Implications for Ambient Environments

Secondary organic aerosol (SOA) plays an important role in the Earth’s energy balance and air quality. Understanding the formation of SOA accompanying the oxidation of volatile organic compounds (VOCs) remains incomplete because of the complex reactions and partitioning mechanisms. In this study, we investigated new particle formation (NPF) of SOA derived from α-pinene ozonolysis in the Harvard Environmental Chamber using the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) with the condensing species as two products, one of lower volatility (LVOC, 0.01 μg m –3 ) and one of higher volatility (SVOC, 1.0 μg m –3 ). A module for the classical nucleation theory was included. The surface tension (σ), bulk diffusivity (D b ), and mass accommodation coefficient (α) were varied until convergence between the simulated and observed particle number size distributions. α of LVOC was constrained as 0.03–0.1 with the reported σ of α-pinene-ozonolysis-derived SOA and Tolman surface correction (~23.0–27.5 mN m –1 ) and Db from 10 –15.5 to >10 –11 cm 2 s –1 . The observed particle number size distribution of α-pinene ozonolysis was well-simulated to derive the physical parameters, including σ, Db, and α, suggesting the possible application to group the reported complex oxidation products and further application for other atmospheric-related oxidations for quantifying NPF.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

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

Improved analytic modeling of the linear Rayleigh–Taylor instability with plasma transport

The Rayleigh–Taylor (R–T) instability with transport effects has long been studied analytically in the neutral fluid regime. How this picture may differ in plasmas has been given less consideration, however. Leveraging previous numerical simulations with plasma viscosity and diffusion of R–T at a binary plasma interface, we show how plasma transport and kinetic effects alter the linear R–T dispersion relation. Our insight informs a much improved analytical dispersion relation, which better matches the simulation data than previously established models. Additionally, we demonstrate the importance of employing asymptotically correct viscosity coefficients in fluid codes for accurately capturing plasma R–T growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Stable Rate of Slip Along the Karakax Section of the Altyn Tagh Fault from Observation of Interglacial and Postglacial Offset Morphology and Surface Dating

Digital elevation maps obtained using TanDEM-X and Pleiades data combined with newly obtained surface age estimates using cosmogenic radionuclide (CRN) and optically simulated luminescence (OSL) methods are used to quantify the slip rate along the western section of the Altyn Tagh Fault in southern Xinjiang. The reconstruction of the conical shape of massive alluvial fans inferred to be from the Eemian (115±7 ka) from CRN dating shows consistent left-lateral offsets of 300±20 m, yielding a slip rate of 2.6±0.3 mm/year. Successive episodes of incision have left cut terraces inset in wide canyons, 10–25 m below the fans' surface. The incision was followed by the deposition of a broad terrace of early Holocene age, which is reincised by modern stream channels. Near the village of Shanxili, a 200 m wide valley is partially dammed by a shutter ridge displaced by the fault. A fill terrace deposited upstream from the ridge has an OSL age of 8.8±0.6 ka. The 23±2 m offset of the riser incising the terrace indicates a minimum postdepositional movement on the fault, yielding a Holocene rate of 2.6±0.5 mm/year, consistent with the 115-ka average slip rate. Furthermore, scarp degradation analysis using mass diffusion reveals a nonlinear relationship between fault displacement and degradation coefficient along the progressively exposed fault scarp, a pattern suggesting either seismic clustering or variable diffusion rate since the Eemian. Together with the Gozha Co-Longmu Co fault to the south, the Karakax section of the Altyn Tagh Fault contributes to the eastward movement of the western corner of Tibet.

58 GEOSCIENCES↗
DOE OSTI2021

Small to Full-Height Scale Comparisons of Cesium Ion Exchange Performance with Crystalline Silicotitanate

Crystalline silicotitanate (CST) ion exchanger is planned to be used to remove cesium (137Cs) from the aqueous phases of Hanford tank wastes in preparation for vitrification at the Waste Treatment and Immobilization Plant (WTP) Pretreatment Facility. Column scale up testing was conducted to evaluate performance of Cs exchange onto the CST. Batch contact testing was conducted to assess exchange kinetics of Cs exchange at four different CST PSDs. Column testing was conducted at three column sizes, small (2.5% full height), medium (12% full height), and full height, to assess Cs load performance behavior. Testing at the small scale was compared to actual Hanford tank waste testing to verify the validity of simulant tests to accurately represent full height column performance. A change in CST particle size was essential in scaling the small column dynamics up to full scale. Based on these results, a determination of intraparticle and film diffusion impacts on overall mass transfer coefficients will allow future modeling of the breakthrough performance at a range of process conditions.

Westesen, Amy M.↗
DOE OSTI2023

Measurements of He-Ar, He-N 2 , and He-air thermal creep slip and accommodation coefficients at high temperatures

Understanding gas flows in capillaries has many applications in modeling the transport of gases in nano-structured, porous, or fractured media. A network of capillaries can often approximate these media, and also information on gas-surface interactions obtained from capillary experiments can be used for modeling flows in these media. Experimental data on flows of different mixtures of He-Ar, He-N 2 , and He-air in the slip regime and in stainless steel capillaries at high temperatures were obtained by using a two-bulb apparatus. An accurate expression for the thermal creep slip coefficients with the Lennard-Jones potential parameters and diffuse-specular reflection gas-surface interaction conditions were then used to obtain the accommodation coefficients for the different gases. The experimental data are best described with values of accommodation coefficients in the range of 0.1–0.3 for He and 0.5–1.0 for Ar, N 2 , and air. Further, the use of values in this range is suggested for modeling gaseous flows in capillaries and nano-structured, porous, or fractured media if other direct measured values for a particular medium are unavailable.

42 ENGINEERING↗
DOE OSTI2021

Mechanical and Thermophysical Properties of 3D-Printed SiC before and after Neutron Irradiation – FY21

This report presents the mechanical and thermophysical properties of 3D-printed SiC before and after neutron irradiation that have been evaluated to assess the fuel matrix material for the Transformational Challenge Reactor (TCR). The TCR fuel form consists of an additively-manufactured silicon carbide (SiC) matrix and uranium nitride tristructural isotropic (UN TRISO) fuel particles, which is manufactured through a newly developed processing route combining binderjet 3D printing, TRISO fuel particle loading, and chemical vapor infiltration (CVI). Because the fuel matrix is a primary component of the TCR core and its response to mechanical and thermal loads during operation is one of the most influential factors on the integrity of TCR core, testing and evaluation have focused on producing mechanical and thermophysical properties data for the binderjet/CVI SiC. Baseline mechanical and thermophysical properties were measured from the disk specimens printed for different and sizes orientations, which included equibiaxial flexural failure strength, elastic constants, thermal diffusivity and conductivity, density, and the coefficient of thermal expansion. Flexural failure strength datasets showed similar Weibull distributions regardless of sample variants including different orientations. The mean failure strengths of the 3D-printed SiC variants were in the range of 280–310 MPa, which are slightly lower than that of the chemical vapor deposition (CVD) SiC. Thermophysical test results showed that specific heat and thermal expansion are not sensitive to the build direction of SiC samples, while thermal conductivity is highly dependent on the build direction and can be correlated to the anisotropic character of the 3D-printed SiC. Neutron irradiation tests were carried out on the 3D-printed 6-mm diameter SiC disk specimens. Irradiation was carried to 2.3 dpa over a temperature range of 360–880°C. No significant degradation in strength was observed in SiC after irradiations in various conditions and with different orientations. Anisotropy that had been observed in the thermal conductivity of 3D-printed SiC prior to irradiation vanished after irradiation as the irradiation defect thermal resistivity accumulated in the material. Electron microscopy of the microstructure after neutron irradiation showed distinct defect morphologies in the heterogenous material, but no evidence for irradiation-induced cracking or degradation in the microstructure was observed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Creating Novel Gas Separation Constructs by Manipulating Free-Volume Distributions in Polymer-Grafted Nanoparticles

Polymer membranes are critical to sustainability applications, and in this work we focused on one key application, the efficient separation of gas mixtures. Despite their widespread use, important challenges remain in the use of polymer membranes in this context – i.e., finding membrane materials with selectively improved transport of the desired mixture component(s) while possessing enhanced strength and improved aging behavior relative to the best current benchmarks. The important separation figures of merit are the gas flux, which is proportional to the permeability, P i = D i × S i (Di is the gas diffusivity and S i its solubility coefficient) and selectivity (i.e., gas purity,α ij = P i /P j ). The goal is to simultaneously maximize P i and α ij . Most research to date has empirically targeted the development of new glassy polymers with improved separation performance. These include thermally rearranged (TR) polymers and polymers of intrinsic microporosity (PIMs).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Estimation of secondary organic aerosol viscosity from explicit modeling of gas-phase oxidation of isoprene and α -pinene

Abstract. Secondary organic aerosols (SOA) are major components of atmospheric fine particulate matter, affecting climate and air quality. Mounting evidence exists that SOA can adopt glassy and viscous semisolid states, impacting formation and partitioning of SOA. In this study, we apply the GECKO-A (Generator of Explicit Chemistry and Kinetics of Organics in the Atmosphere) model to conduct explicit chemical modeling of isoprene photooxidation and α-pinene ozonolysis and their subsequent SOA formation. The detailed gas-phase chemical schemes from GECKO-A are implemented into a box model and coupled to our recently developed glass transition temperature parameterizations, allowing us to predict SOA viscosity. The effects of chemical composition, relative humidity, mass loadings and mass accommodation on particle viscosity are investigated in comparison with measurements of SOA viscosity. The simulated viscosity of isoprene SOA agrees well with viscosity measurements as a function of relative humidity, while the model underestimates viscosity of α-pinene SOA by a few orders of magnitude. This difference may be due to missing processes in the model, including autoxidation and particle-phase reactions, leading to the formation of high-molar-mass compounds that would increase particle viscosity. Additional simulations imply that kinetic limitations of bulk diffusion and reduction in mass accommodation coefficient may play a role in enhancing particle viscosity by suppressing condensation of semi-volatile compounds. The developed model is a useful tool for analysis and investigation of the interplay among gas-phase reactions, particle chemical composition and SOA phase state.

54 ENVIRONMENTAL SCIENCES↗