Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffused interface model”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Understanding LIB Battery Electrodes Through Classical Electrochemical Interface Theory

Classical models of the electrochemical interface can be applied to complex and dynamic electrodes to understand important mechanisms that are relevant to device-level phenomena like calendar aging or cycle life. Here, we show that negative alloy electrodes with wide electrochemical windows cannot be assumed to have static interfaces throughout all states of charge. Under high states of charge, the interface -including the solid electrolyte interphase - has characteristics that resemble a highly polarized electrode far from the point of zero charge. Under these conditions, the interface can be understood through the Helmholtz model. At lower states of charge, the interface capacitance decreases and the space charge layer length increases. This transformation resembles a metal electrode approaching the point of zero charge which is understood through the Gouy-Chapman-Stern-Grahame model of the electrochemical interface. This transformation happens at both silicon and carbon-coated silicon interfaces. From voltage-limited cycling experiments, in the limit of the diffuse double layer, the silicon electrode impedance rises significantly in the first few cycles. This rise is related to the continuous electrochemical reduction of electrolyte components from an interface that is not electronically screened from the electrolyte. In other words, an electrode interface for batteries should resemble a 'Helmholtz-like' structure.

anode↗
DOE OSTI2023

Modeling and Characterization of Surface Discharges in Insulating Material for Spacers: Electrode Shape, Discharge Mode, and Revision of the Creepage Concept

In the design of MV AC and DC spacers, the predominant factors are surface and interface conditions. Design is generally carried out on specifications and standards which are based on long-term experience and lab testing. However, the diffusion of power electronics with a trend to increase electric field, switching frequency, and rise time to achieve higher power density calls for an innovative, global approach to optimized insulation system design. A new methodology, based on field simulation, discharge modeling, and partial discharge inception measurements, called the three-leg approach, can form the basis to optimize insulation design for any type of supply voltage waveform. This paper focuses on the influence of the type of electrode on the inception and phenomenology of surface discharges and, as a consequence, on the interpretation of the results used for application of the three-leg approach. It is demonstrated that a typical electrode system used for insulating material testing can generate both gas and surface discharges at the triple point, when the electrodes have a smooth profile that is used to avoid corona or flashover. Hence, testing partial discharge may not provide a straightforward indication of the surface discharge inception and, thus, be partially misleading for insulation design. Another takeover is that such analysis must benefit from PD testing tools endowed with analytics able to provide automatic identification of the type of defect generating PD, i.e., internal, surface, and corona, since design and remedy actions can be taken, and adequate insulating materials developed, only knowing the type of source generating PD. Hence, testing partial discharge may not provide a straightforward indication of surface discharge inception and, thus, be partially misleading for insulation design. In addition to the importance of the three-leg approach to favor reliable and optimized design of insulation systems, there is a clear need to have a PD testing tool endowed with analytics. It should preferably be able to provide automatic identification of the type of defect generating PD, i.e., internal, surface, and corona.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Development of Corrosion Resistant Coatings for Structural Materials for Liquid Fueled Molten Salts Reactors Applications (Final Report)

Current structural alloys, code certified for the MSRs temperature ranges, contain high levels of chromium, making them highly susceptible to molten halide salt corrosion. One potential solution to circumvent the need for code certification of novel alloys, which is expensive and time consuming, is to design claddings that will protect the underlying code certified materials from corrosion damage during operation. In this work, we examined the corrosion of Ni and Cu electroplated, Ni and CuNi weld overlay, Mo-laser clad, and carburized claddings on SS316H for use in molten salt reactor environments. After characterization of the cladded materials (Task 1), static corrosion tests were used to assess corrosion resistance of the claddings in typical molten fluoride salts (Task 2). The corrosion tests were performed in molten FLiNaK at 700°C up to 1000 hours. Pre- and post-corrosion Scanning Electron Microscopy (SEM), energy-dispersive-spectroscopy (EDS), Transmission Electron Microscopy (TEM), X-ray Diffraction (XRD) and glow-discharge-optical-emission-spectroscopy (GDOES) were performed on cladding cross-sections and surfaces to evaluate degradation. The Cu and Ni electroplated samples, as well as the carburized samples, showed excellent corrosion resistance relative to the bare SS316H. To assess high temperature cladding stability, ageing experiments were performed at temperatures up to 900 C in inert atmosphere for the electroplated and carburized samples and a diffusion model was developed to predict long term cladding behavior (Task 3). It was found that the Cu cladding was basically insensitive to the high temperature ageing, except for small secondary phases forming at the interface. On the other hand, the Ni electroplated cladding experienced significant interdiffusion. Nevertheless, the gain in corrosion resistance for a 100μm Ni electroplated cladding is phenomenal, with more than 50% reduction in chromium dissolution from the substrate material for the first 15 years of salt exposure at 700°C. To assess radiation resistance and phase stability, high-temperature 4MeV Ni heavy ion irradiation was performed across the cladding/substrate interface up to 50 displacement-per-atom (DPA) at 500°C and 700°C to assess the phase stability and irradiation behaviors of the electroplated systems (Task 4). The Ni and Cu electroplated systems did not experience void swelling at the contrary to the SS316H substrate due to their nanocrystalline nature. Some level of recrystallization was observed in the cladding, as well as radiation induced segregation and enhanced diffusion. The interface acted as a potent sink for point defects with the presence of a void denuded zone. The mechanical properties of the claddings were assessed using thermal shock, micro-indentation, nano-indentation, and four-point bend testing experiments (Task 5). These experiments were performed on the claddings in as-received, corroded, irradiated, and thermally aged states to determine the effects of typical molten salt reactor environments on cladding integrity. No significant mass loss was observed after repeated thermal shocks. While the electroplated samples softened after high-temperature ageing, irradiation hardening compensate this effect, such that there is little different with the as-received materials. The results of these experiments suggest that the electroplated (copper/nickel) show the greatest promise for application in molten salt reactor development. These claddings prevented any chromium dissolution from occurring during the static corrosion experiments. Additionally, they demonstrated favorable properties for high temperature diffusion, phase stability, interfacial mechanical properties, and irradiation resistance. Weld-overlay cladding are also of great interest and could reach properties similar to the electroplated samples upon optimization (especially with multiple weld passes). The carburized SS316H showed excellent behavior as well, but more work needs to be done to assess their long-term stability. Finally, the Mo-laser clad system was not pursued further due to the manufacturing challenges. It is also worth noting that the Ni cladding systems should behave relatively well in terms of weldability since the Ni-weld overlay microstructure and associated corrosion rate are sound. While the NiCu weld-overlay has even lower corrosion rates than the Ni-weld overlay, more studies on welding of Cu-electroplated systems are necessary since Cu clusters are known to embrittle steels.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unravelling Microstructure Selection in an Additively Manufactured Eutectic High‐Entropy Alloy

High-entropy alloys (HEAs) are promising candidates for advanced structural applications due to their excellent mechanical properties. Additive manufacturing (AM), with its rapid solidification conditions, enables the creation of unique nonequilibrium microstructures. To fully leverage the synergy between AM and HEAs, understanding how processing affects structure and properties is essential. Here, how solidification rate influences microstructure evolution and phase transformation pathway in laser additively manufactured AlCrFe2Ni2 eutectic HEAs is investigated. By increasing the laser scan speed and hence the solidification rate, distinct solidification modes evolving from coupled eutectic to anomalous eutectic and eventually to single-phase solidification are revealed. These transitions result in distinct microstructures and a wide range of mechanical properties. Thermodynamic modeling and molecular dynamics simulations reveal that low cooling rates allow for sufficient atomic diffusion and phase separation, facilitating coupled eutectic growth. In contrast, rapid cooling suppresses diffusion and destabilizes the solid–liquid interface, promoting anomalous or single-phase solidification. This integrated experimental and computational approach provides a multiscale understanding of solidification mechanisms in HEAs and underscores how kinetic effects can over-ride thermodynamic predictions under nonequilibrium conditions. Furthermore, these results demonstrate that AM can serve as a powerful tool to design HEAs with tailored microstructures and properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Distinguishing Surface and Bulk Reactivity: Concentration-Dependent Kinetics of Iodide Oxidation by Ozone in Microdroplets

Iodine oxidation reactions play an important role in environmental, biological, and industrial contexts. The multiphase reaction between aqueous iodide and ozone is of particular interest due to its prevalence in the marine atmosphere and unique reactivity at the air–water interface. Here, we explore the concentration dependence of the I – + O 3 reaction in levitated microdroplets under both acidic and basic conditions. To interpret the experimental kinetics, molecular simulations are used to benchmark a kinetic model, which enables insight into the reactivity of the interface, the nanometer-scale subsurface region, and the bulk interior of the droplet. For all experiments, a kinetic description of gas- and liquid-phase diffusion is critical to interpreting the results. We find that the surface dominates the iodide oxidation kinetics under concentrated and acidic conditions, with the reactive uptake coefficient approaching an upper limit of 10 –2 at pH 3. In contrast, reactions in the subsurface dominate under more dilute and alkaline conditions, with inhibition of the surface reaction at pH 12 and an uptake coefficient that is 10× smaller. In conclusion, the origin of a changing surface mechanism with pH is explored and compared to previous ozone-dependent measurements.

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

Challenges and approaches to interpretive modeling of boundary plasma and neutral transport in a closed, pumped divertor

An experimental discharge from the DIII-D tokamak is modeled using the SOLPS-ITER code suite and compared against measurements in the pumped and relatively closed upper divertor. Uncertainties of boundary plasma simulations are identified by attempting to match code inputs to experimental conditions, including iteratively solving transport coefficients to match upstream experimental profiles using varying quantities of core particle flux, different pumping models, and various assumptions of ion thermal transport. Simulated boundary conditions for particle injection at the core interface are shown to be relevant to the plasma solution at the divertor targets, even if upstream transport is modified so that plasma profiles are comparatively similar, although seperatrix density is not held constant. When upstream plasma profiles are matched to experimental measurements by varying diffusive transport coefficients, using either poloidally symmetric or ballooning structure, the model finds a majority of injected energy being transported radially off the computational domain, in conflict with experimental radiated power measurements and heat flux measurements at the divertor target. Imposing a maximum thermal diffusivity or radially shifting the experimental separatrix location of the fitted profiles to increase power conducted to the targets by increasing the upstream electron temperature does not significantly modify this result. Including a thermalizing plenum volume in the simulation domain is shown to maintain the experimental volumetric pumping rate without knowing the neutral energy distribution incident on the pump duct a priori. By modifying transport parameters to match different assumptions for ion temperature, downstream neutral pressure changes by more than a factor of two, suggesting that attention to ion thermal transport may be a critical parameter for simulations to accurately resolve recycling and neutral transport, particularly in a closed divertor geometry. In addition to quantifying various modeling uncertainties, this work motivates both further experimental study and modeling improvements to improve predictive capabilities.

divertor↗
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 OSTI2020

Scaling Laws of Exciton Recombination Kinetics in Low Dimensional Halide Perovskite Nanostructures

Carrier recombination is a crucial process governing the optical properties of a semiconductor. Although various theoretical approaches have been utilized to describe carrier behaviors, a quantitative understanding of the impact of defects and interfaces in low dimensional semiconductor systems is still elusive. Here, we develop a model system consisting of chemically tunable, highly luminescent halide perovskite nanocrystals to illustrate the role of carrier diffusion and material dimensionality on the carrier recombination kinetics and luminescence efficiency. Our advanced synthetic methods provide a well-controlled colloidal system consisting of nanocrystals with different aspect ratios, halide compositions, and surface conditions. Using this system, we reveal the scaling laws of photoluminescence quantum yield and radiative lifetime with respect to the aspect ratio of nanocrystals. The scaling laws derived herein are not only a phenomenological observation but proved a powerful tool disentangling the carrier dynamics of microscopic systems in a quantitative and interpretable manner. The investigation of our model system and theoretical formulation bring to light the dimensionality, as a hidden constraint on carrier dynamics, and identify the diffusion length as an important parameter that distinguishes nanoscale and macroscale carrier behaviors. Furthermore, the conceptual distinction in carrier dynamics in different dimensionality regimes informs new design rules for optical devices where complex microstructures are involved.

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

Precipitation-site competition in duplex stainless steels: Cu clusters vs spinodal decomposition interfaces as nucleation sites during thermal aging

Competing microstructural evolution mechanisms can exist simultaneously when duplex stainless steels are operating for several decades in a high temperature service environment. Such competition between different microstructural evolution pathways can be difficult to ascertain using simple model alloy systems necessitating detailed structural and compositional analysis of phase transformation mechanisms in complex alloys Thus, in this study duplex stainless steels with complex but well understood chemistries were used to investigate the relative importance of different heterogeneous nucleation sites – specifically spinodal decomposition and Cu clustering – on precipitation of Ni-Si-Mn particles during long term thermal aging. Precipitation of Ni-Si-Mn particles in duplex stainless steels and ferritic steels during thermal aging and irradiation is known to have profound effect on mechanical properties. Using duplex stainless steels with custom modified compositions along with atom probe tomography characterization and first-passage kinetic Monte Carlo (FPKMC) simulations it is revealed that while the interface between Cr and Fe formed during spinodal decomposition can be a pathway for diffusion of Ni, Si, Mn, and Cu, it is not a site that encourages precipitation of Ni-Si-Mn rich particles. Instead, the presence of a higher concentration of Cu in these steels leads to formation of small Cu clusters with high energy interfaces that act as nucleation sites for Ni-Si-Mn particles because of its strong interaction particularly with Ni increasing the chemical driving force for precipitation. These results will help inform predictive models for the use of duplex stainless steels and other precipitation-hardened alloys for extended operation at high temperatures.

Duplex stainless steel, Thermal aging, Precipitati↗
DOE OSTI2024

Measuring thermal diffusivity and gap conductance in uranium nitride and Zircaloy relevant for microreactor applications

Heat transfer across nuclear fuels and structural interfaces is an important factor for evaluating the performance of nuclear power systems. Specifically, heat generated as nuclear fuel fissions must be transported through the cladding material and through the reactor to reach the steam turbine for power generation. As new microreactor designs emerge, maximizing the efficiency of this heat transfer process becomes crucial to make them commercially viable. This article examines thermal diffusivity and gap conductance in uranium nitride (UN) fuel and Zircaloy-4 (Zry4) cladding using light flash analysis (LFA). Thermal diffusivity measurements were made on monolithic UN pellets and Zry4 exposed to carbon at peak operating temperatures of microreactors and show that carbon ingress has a minimal effect on thermal diffusivity when compared with identical materials not exposed to carbon. Evaluation of gap conductance at the UN-Zry4 interface was done using one-dimensional two-layer thermal transport models as a function of applied pressure. Here the results show that increasing pressure on the UN-Zry4 interface leads to gains in gap conductance per unit area in fuel-cladding assemblies at microreactor operating temperatures. While many other variables are expected to influence UN-Zry4 interfacial gap conductance (e.g. contact surface roughness, porosity, localized heating, environmental gas pressure), the work offers a demonstration of using a conventional LFA apparatus to determine this parameter at elevated temperatures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Strain-Driven Surface Reactivity in Magnesium-Ion Battery Cathodes

Successful deployment of a Mg-ion battery requires cathodes that can achieve reversible Mg intercalation and high energy density. Recent theoretical and experimental studies indicated that the overall transport is likely limited by sluggish Mg transport at the cathode-electrolyte interface and not Mg diffusion through bulk. Here, in this work, we investigated the surface electrochemical activity of Mg ions by using a spinel-structured manganese oxide thin-film model system and in situ X-ray scattering. In combination with post-mortem microscopy analysis, we found that magnesium insertion was more favorable than subsequent extraction near the surface of the Mg x Mn 2 O 4 film, resulting in overmagnesiation, and eventually amorphization of the surface. This structural irreversibility and high overpotential required for Mg extraction could explain significant voltage hysteresis and Mg surface enrichment previously observed in bulk cathodes. Density functional theory calculations suggested that the tendency for the Mg surface enrichment could be associated with Mg diffusion kinetics, which varies with the strain state evolved due to constrained film volume change during Mg insertion and extraction. Particularly, out-of-plane Mg migration was predicted to be favorable in the tensile strain rather than in the compressive case.

25 ENERGY STORAGE↗
DOE OSTI2022

Uptake of N 2 O 5 by aqueous aerosol unveiled using chemically accurate many-body potentials

The reactive uptake of N 2 O 5 to aqueous aerosol is a major loss channel for nitrogen oxides in the troposphere. Despite its importance, a quantitative picture of the uptake mechanism is missing. Here we use molecular dynamics simulations with a data-driven many-body model of coupled-cluster accuracy to quantify thermodynamics and kinetics of solvation and adsorption of N 2 O 5 in water. The free energy profile highlights that N 2 O 5 is selectively adsorbed to the liquid–vapor interface and weakly solvated. Accommodation into bulk water occurs slowly, competing with evaporation upon adsorption from gas phase. Leveraging the quantitative accuracy of the model, we parameterize and solve a reaction–diffusion equation to determine hydrolysis rates consistent with experimental observations. We find a short reaction–diffusion length, indicating that the uptake is dominated by interfacial features. The parameters deduced here, including solubility, accommodation coefficient, and hydrolysis rate, afford a foundation for which to consider the reactive loss of N 2 O 5 in more complex solutions.

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

Micro-continuum approach for mineral precipitation

Abstract Rates and extents of mineral precipitation in porous media are difficult to predict, in part because laboratory experiments are problematic. It is similarly challenging to implement numerical methods that model this process due to the need to dynamically evolve the interface of solid material. We developed a multiphase solver that implements a micro-continuum simulation approach based on the Darcy–Brinkman–Stokes equation to study mineral precipitation. We used the volume-of-fluid technique in sharp interface implementation to capture the propagation of the solid mineral surface. Additionally, we utilize an adaptive mesh refinement method to improve the resolution of near interface simulation domain dynamically. The developed solver was validated against both analytical solution and Arbitrary Lagrangian–Eulerian approach to ensure its accuracy on simulating the propagation of the solid interface. The precipitation of barite (BaSO 4 ) was chosen as a model system to test the solver using variety of simulation parameters: different geometrical constraints, flow conditions, reaction rate and ion diffusion. The growth of a single barite crystal was simulated to demonstrate the solver’s capability to capture the crystal face specific directional growth.

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

The nanoscale distribution of copper and its influence on charge collection in CdTe solar cells

For decades, Cu has been the primary dopant in CdTe photovoltaic absorbers. Typically, Cu acceptor concentrations in these devices are on the order of 1 1014 cm-3, which has made it notoriously difficult to directly correlate nanoscale Cu distributions to the local charge transport properties of these devices. To measure and correlate these properties, measurement techniques require high sensitivity to elemental concentration, large penetration depth, and operando compatibility. Techniques such as secondary-ion mass spectroscopy and X-ray energy dispersive spectroscopy are widely adopted to measure Cu concentrations, but they are limited by penetration depth, sensitivity, or spatial resolution. Additionally, they lack the operando capabilities required to correlate one-to-one Cu concentrations to electrical performance. In this work, correlative X-ray microscopy is used to investigate the spatial distribution of Cu and its impact on charge collection through the depth and breadth of CdTe photovoltaic devices. Plan-view, nanoscale X-ray fluorescence maps clearly demonstrate the spatial segregation of copper around regions thought to be CdTe grain boundaries. Complementary cross-section imaging unveils the transition of the maximum charge-collection efficiency from the ZnTe-CdTe interface to the CdS-CdTe interface as a function of Cu incorporation. The copper concentration through the depth of the CdTe layer is characterized by slow and fast diffusion components, and cross-section charge-transport modeling shows that the experimentally obtained charge collection can be explained by the modeled acceptor distribution through the depth of the CdTe layer.

14 SOLAR ENERGY↗
DOE OSTI2022

The Critical Stack Pressure to Alter Void Generation at Li/Solid-Electrolyte Interfaces during Stripping

The lithium stripping process generates vacancies, which may accumulate as voids and lead to uneven current distribution and dendrite growth in the following plating cycles. A stack pressure is typically required during stripping, but how to optimize the stack pressure is not clear. In this work, extremely lithiophilic Li/Li 2 O and lithiophobic Li/LiF interfaces were used to reveal the combining effect of interface interaction and stack pressure induced lithium creep on the stripping critical current density (CCD). A multiscale simulation scheme with Density Functional Theory (DFT), kinetic Monte Carlo (KMC) simulations, and an analytical model was developed. The analytical model predicted lithiophobic interfaces require a higher stack pressure than lithiophilic interfaces to reach the same CCD. The KMC simulations also showed higher stack pressure is needed at lithiophobic interfaces to accelerate Li vacancy diffusion into the bulk and maintain a flat surface. This stack pressure needs to be high enough to alter the Li forward-and-backward hopping barriers at the interface. This multiscale simulation scheme illustrates the importance to include the chemical-mechanical effects during Li stripping morphology evolution. It can be used to design ideal interlayer coating materials to maintain a flat Li surface during cycling.

25 ENERGY STORAGE↗
DOE OSTI2025

Selective Recovery of Critical Minerals from Simulated Electronic Wastes Via Reaction‐Diffusion Coupling

Abstract Atom‐ and energy‐efficient chemical separations are urgently needed to meet the surging demand for critical materials that has strained supply chains and threatened environmental damage. In this study, we used reaction‐diffusion coupling to separate iron, neodymium, and dysprosium ions from model feedstocks of permanent magnets, which are typically found in electronic wastes. Feedstock solutions were placed in contact with a hydrogel loaded with potassium hydroxide and/or dibutyl phosphate, resulting in complex precipitation patterns as the various metal ions diffused into the reaction medium. Specifically, we observed the precipitation of up to 40 mM of iron from the feedstock, followed by the enrichment of 73 % dysprosium, and the extraction of >95 % neodymium product at a further distance from the solution‐gel interface. We designed a series of experiments and simulations to determine the relevant ion diffusivities, D Nd =5.4×10 −10 and D Dy =5.1×10 −10 m 2 /s, and precipitation rates, k Nd =1.0×10 −5 and k Dy =5.0×10 −3 m 9 mol −3 s −1 , which enabled a numerical model to be established for predicting the distribution of products in the reaction medium. Our proof‐of‐concept study validates reaction‐diffusion coupling as an effective and versatile approach for critical materials separations, without relying on ligands, membranes, resins, or other specialty chemicals.

Wang, Qingpu [Physical and Computational Sciences ↗
DOE OSTI2021

Mass accommodation and gas–particle partitioning in secondary organic aerosols: dependence on diffusivity, volatility, particle-phase reactions, and penetration depth

Abstract. Mass accommodation is an essential process for gas–particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as being incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface–bulk exchange, it can be important for viscous semi-solid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs–Sutugin approximation of mass-transport kinetics at the gas–particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Model for Simulating Aerosol Interactions and Chemistry, MOSAIC; Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semi-solid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in the case of low-volatility compounds, whereas it can decrease by several orders of magnitude in the case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs–Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient, regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models. While kinetic limitations may not be critical for partitioning into liquid SOA particles in the planetary boundary layer (PBL), the effects are likely important for amorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT) as well as in the PBL at low relative humidity and low temperature.

54 ENVIRONMENTAL SCIENCES↗