Search NASA⌕ Search

SEARCH · Search NASA

Results for “matrix diffusion”

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 OSTI2018

Developing a Novel Hydrogen Sponge with Ideal Binding Energy and High Surface Area for Practical Hydrogen Storage

This Phase I (5 quarters) research project was to examine the validity of a new class of boron-containing polymer (B-polymer) frameworks, serving as the adsorbents for the practical onboard H2 storage applications. Three B-polymer frameworks were synthesized and investigated, which include B-poly(butyenylstyrene) (B-PBS) framework (A), B-poly(phenyldiacetyene) (B-PPDA) framework (B), and B-poly(phenyltriacetylene) (B-PPTA) framework (C). They are 2-D polymer structures with the repeating cyclic units that spontaneously form open morphology and the B-doped (p-type) π-electrons delocalized surfaces. The ideal B-polymer framework shall exhibit open micropores (pore size in the range of 1-1.5nm) with high surface area (>3000 m 2 /g), and the B-dopants in the conjugated framework shall provide high surface energy for interacting with H 2 molecules (an ideal H 2 binding energy in the range of 15-25 kJ/mol). The pore size distribution and H2 binding energy were investigated at both Penn State and NREL laboratories. So far, the experimental results show the successful synthesis of B-polymer frameworks with the relatively well-defined planar (2-D) structures. The intrinsically formed porous morphology exhibits a broad pore size distribution (in the range of 0.5-10 nm) with specific surface area (~1000 m 2 /g). The miss-alignment between 2-D layers may block some micropore channels and limit gas diffusion throughout the entire matrix. In addition, the 2-D planar conjugated structure may also allow free π-electrons delocalization throughout the framework, which significantly reduces the acidity of B-moieties (electron-deficiency).The resulting 2-D B-polymer frameworks only exhibit a small increase of H 2 binding energy in the range of 8-9 KJ/mole (quite constant over the whole sorption range).

08 HYDROGEN↗
DOE OSTI2026

Impact of Annealing on Superconducting Properties of Tantalum Thin Films

The pursuit for higher-coherence superconducting qubits has led to the adoption of new materials and innovative fabrication techniques. Among the commonly used superconducting materials, tantalum-based or tantalum-capped superconducting qubits have demonstrated the longest lifetimes, outperforming niobium devices even in the presence of their native oxides. Given the losses introduced by the surface oxide, one method to further improve performance is to reduce or eliminate the native oxide on tantalum through UHV annealing. In this work, we examine the effect of dissolved oxygen as a function of annealing temperature on the superconducting properties of tantalum thin films. Dissolved oxygen is quantified using atom probe tomography (APT) and secondary ion mass spectrometry (TOF-SIMS) and correlated with electromagnetic transport measurements performed using a Physical Property Measurement System (PPMS). This study observes that annealing decomposes the native oxide and drives oxygen diffusion into the film matrix. Increasing oxygen content decreases the superconducting transition temperature (Tc) and the residual resistivity ratio (RRR). The critical field decreases with the incorporation of oxygen, while the film’s coherence length decreases relatively by 5%. These results, as well as findings from ongoing studies linking fabrication processes to the structural and chemical properties of the superconducting film and Josephson junction, provide insight into effective strategies for improving the performance of superconducting devices.

Shinde, Maithilee [IIT, Chicago; Fermilab]↗
DOE OSTI2023

Impact of phase state and non-ideal mixing on equilibration timescales of secondary organic aerosol partitioning

Abstract. Evidence has accumulated that secondary organic aerosols (SOAs) exhibit complex morphologies with multiple phases that can adopt amorphous semisolid or glassy phase states. However, experimental analysis and numerical modeling on the formation and evolution of SOA still often employ equilibrium partitioning with an ideal mixing assumption in the particle phase. Here we apply the kinetic multilayer model of gas–particle partitioning (KM-GAP) to simulate condensation of semi-volatile species into a core–shell phase-separated particle to evaluate equilibration timescales of SOA partitioning. By varying bulk diffusivity and the activity coefficient of the condensing species in the shell, we probe the complex interplay of mass transfer kinetics and the thermodynamics of partitioning. We found that the interplay of non-ideality and phase state can impact SOA partitioning kinetics significantly. The effect of non-ideality on SOA partitioning is slight for liquid particles but becomes prominent in semisolid or solid particles. If the condensing species is miscible with a low activity coefficient in the viscous shell phase, the particle can reach equilibrium with the gas phase long before the dissolution of concentration gradients in the particle bulk. For the condensation of immiscible species with a high activity coefficient in the semisolid shell, the mass concentration in the shell may become higher or overshoot its equilibrium concentration due to slow bulk diffusion through the viscous shell for excess mass to be transferred to the core phase. Equilibration timescales are shorter for the condensation of lower-volatility species into semisolid shell; as the volatility increases, re-evaporation becomes significant as desorption is faster for volatile species than bulk diffusion in a semisolid matrix, leading to an increase in equilibration timescale. We also show that the equilibration timescale is longer in an open system relative to a closed system especially for partitioning of miscible species; hence, caution should be exercised when interpreting and extrapolating closed-system chamber experimental results to atmosphere conditions. Our results provide a possible explanation for discrepancies between experimental observations of fast particle–particle mixing and predictions of long mixing timescales in viscous particles and provide useful insights into description and treatment of SOA in aerosol models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Effective Lanthanide Diffusivity through U-Zr Metallic Fuel

Fuel-cladding chemical interaction (FCCI) is one of the main concerns for the performance of U-Zr metallic fuels. Although FCCI results in complex phase transformations and morphological changes in both the fuel and cladding sides of the interface, the formation of a brittle layer in the cladding side of the interface (sometimes referred to as wastage) has one of the most significant impacts on the cladding mechanical integrity. This brittle layer is associated with the formation of intermetallic compounds between the cladding constituents Fe, Cr and the lanthanide fission product species. The formation of this cladding wastage layer has been the focus of recent efforts to develop a mechanistic modeling framework, including production of lanthanides in the fuel, transport through the fuel to the fuel-cladding interface, and phase transformation to form intermetallic phases. To calculate the rate of transport of lanthanides through the fuel, a multi-scale computational approach has been used. Atomistic calculations have been performed to determine the diffusivities of Nd, the most prevalent lanthanide species, through the solid fuel matrix and along pore surfaces. These diffusivities have been used as input parameters for a mesoscale model to determine an effective diffusion coefficient that accounts for the formation of porosity, its interconnection, and infiltration with bond sodium from the fuel-cladding gap. The effective diffusivity has been used in engineering-scale simulations using the BISON fuel performance code, and BISON simulations have been compared to representative validation cases from the EBR-II reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Multiscale Modeling of CO 2 Migration and Trapping in Fractured Reservoirs with Validation by Model Comparison and Real-Site Applications (Final Technical Report)

This report documents the accomplishments achieved during the project titled “Multiscale Modeling of CO 2 Migration and Trapping in Fractured Reservoirs with Validation by Model Comparison and Real-Site Applications” funded by the US Department of Energy, Office of Fossil Energy. The objectives of the project were to develop modeling capabilities for predicting CO 2 and brine migration in fractured reservoirs during geologic carbon storage and to investigate the feasibility of carbon storage in fractured reservoirs using the newly developed modeling capabilities. To achieve these objectives new mass transfer functions were developed to adapt existing dual-continuum approaches – a commonly used approach in hydrocarbon reservoir modeling – to the CO 2 -brine system. In the dual-continuum approach the undisturbed rock matrix and the fractures are modeled as two separate continua, which are coupled through mass exchange of CO 2 and brine between the two continua. Mass transfer functions represent the mass exchange between fractures and the rock matrix, without requiring highly resolved numerical grids of individual fractures, thus allowing large modeling domains relevant to geologic carbon storage questions (tens of meters to hundreds of kilometers). Mass transfer functions for the exchange of both free-phase CO 2 and brine between fractures and rock matrix were developed for both the injection phase and the post-injection phase. During the injection phase the fractures are filled with CO 2 , while the rock matrix is initially filled with brine. Therefore, a mass transfer function was developed that is based on vertical displacement of brine by CO 2 in the rock matrix due to gravitational forces. Once CO 2 injection ceases and CO 2 migrates away from the direct vicinity of the injection site, the fractures will begin to re-fill with brine, while the rock matrix will have high CO 2 saturations. For these conditions a mass transfer functions was developed that is based on spontaneous imbibition, as the pore spaces in the rock matrix are much smaller than in the fractures. For formations with low permeability or high capillary entry pressure of the rock matrix, free-phase CO 2 may not be able to enter the rock matrix. However, free-phase CO 2 at the interface between fractures and rock matrix can dissolve into brine resident in the rock matrix. Therefore, additional mass transfer functions were developed based on aqueous-phase diffusion of CO 2 through the rock matrix. These mass transfer functions are based on existing solutions for diffusion and were expanded to be applicable to a variety of block shapes relevant to dual-continuum modeling. Implementations of the newly developed mass transfer functions showed that they give accurate results when compared to highly resolved models. In addition to developing new transfer functions, a new vertically-integrated dual-continuum modeling framework was developed. This development takes advantage of fast vertical segregation of CO 2 and brine in the fractures due to the density difference between CO 2 and brine and the high permeability of the fractures. Therefore, a vertically-integrated vertical equilibrium model for the fractures is coupled with a vertically-distributed model for the rock matrix. The newly developed approach significantly reduces the computational effort, but is only accurate for fracture permeabilities of 100 mD or higher. Finally, the newly developed models were applied to real and hypothetical geologic carbon storage sites to investigate the feasibility of using fractured reservoirs for carbon storage. The general results were, that storage of CO 2 – both as free-phase and dissolved – is feasible, but that storage capacities are likely lower than in unfractured systems. While the high permeability of the fractures reduces the injection pressure (or permits higher injection rates), fast migration of CO 2 may limit the overall storage volume, because CO 2 may escape laterally with not enough of the rock matrix being swept for effective storage.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

Diffusive hydrodynamics from integrability breaking

We describe here the crossover from generalized to conventional hydrodynamics in nearly integrable systems. Integrable systems have infinitely many conserved quantities, which spread ballistically, in general. When integrability is broken, only a few of these conserved quantities survive. The remaining conserved quantities are generically transported diffusively; we further derive a compact and general diffusion equation for these. The diffusion constant depends on the matrix elements of the integrability-breaking perturbation; for a certain class of integrability-breaking perturbations, including long-range interactions, the diffusion constant can be expressed entirely in terms of generalized hydrodynamic data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

One-Dimensional Anomalous Diffusion of Gold Nanoparticles in a Polymer Melt

We investigated the dynamics of polymer-grafted gold nanoparticles loaded into polymer melts using x-ray photon correlation spectroscopy. For low molecular weight host matrix polymer chains, normal isotropic diffusion of the gold nanoparticles is observed. For larger molecular weights, anomalous diffusion of the nanoparticles is observed that can be described by ballistic motion and generalized Levy walks, similar to those often used to discuss the dynamics of jammed systems. Under certain annealing conditions, the diffusion is one-dimensional and related to the direction of heat flow during annealing and is associated with an dynamic alignment of the host polymer chains. Molecular dynamics simulations of a single gold nanoparticle diffusing in a partially aligned polymer network semiquantitatively reproduce the experimental results to a remarkable degree. Here, the results help to showcase how nanoparticles can under certain circumstances move rapidly in polymer networks.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Transport properties of C and O in UN fuels

Uranium nitride fuel is considered for fast reactors (GEN-IV generation and space reactors) and for light water reactors as a high-density fuel option. Despite this large interest, there is a lack of information about its behavior for in-pile and out-of-pile conditions. From the present literature, it is known that C and O impurities have significant influence on the fuel performance. Here we perform a systematic study of these impurities in the UN matrix using electronic-structure calculations of solute-defect interactions and microscopic jump frequencies. These quantities were calculated in the DFT + U approximation combined with the occupation matrix control scheme, to avoid convergence to metastable states for the 5f levels. The transport coefficients of the system were evaluated with the self-consistent mean-field theory. It is demonstrated that carbon and oxygen impurities have different diffusion properties in the UN matrix, with O atoms having a higher mobility, and C atoms showing a strong flux coupling anisotropy. The kinetic interplay between solutes and vacancies is expected to be the main cause for surface segregation, as incorporation energies show no strong thermodynamic segregation preference for (001) surfaces compared with the bulk.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Modeling the influence of preferential flow on the spatial variability and time-dependence of mineral weathering rates

Inferences of weathering rates from laboratory and field observations suggest significant scale and time-dependence. Preferential flow induced by heterogeneity (manifest as permeability variations or discrete fractures) has been suggested as one potential mechanism causing scale/time-dependence. In this paper, we present a quantitative evaluation of the influence of preferential flow on weathering rates using reactive transport modeling. Simulations were performed in discrete fracture networks (DFNs) and correlated random permeability fields (CRPFs), and compared to simulations in homogeneous permeability fields. The simulations reveal spatial variability in the weathering rate, multidimensional distribution of reactions zones, and the formation of rough weathering interfaces and corestones due to preferential flow. In the homogeneous fields and CRPFs, the domain-averaged weathering rate is initially constant as long as the weathering front is contained within the domain, reflecting equilibrium-controlled behavior. The behavior in the CRPFs was influenced by macrodispersion, with more spread-out weathering profiles, an earlier departure from the initial constant rate and longer persistence of weathering. DFN simulations exhibited a sustained time-dependence resulting from the formation of diffusion-controlled weathering fronts in matrix blocks, which is consistent with the shrinking core mechanism. A significant decrease in the domain-averaged weathering rate is evident despite high remaining mineral volume fractions, but the decline does not follow a math formula dependence, characteristic of diffusion, due to network scale effects and advection-controlled behavior near the inflow boundary. Finally, the DFN simulations also reveal relatively constant horizontally averaged weathering rates over a significant depth range, challenging the very notion of a weathering front.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Self-learning kinetic Monte Carlo simulations of diffusion in ferromagnetic α -Fe–Si alloys

Diffusion in α-Fe-Si alloys is studied using AKSOME, an on-lattice self-learning KMC code, in the ferromagnetic state. Si diffusivity in the α-Fe matrix were obtained with and without the magnetic disorder in various temperature ranges. In addition we studied vacancy diffusivity in ferromagnetic α-Fe at various Si concentrations up to 12.5at.% in the temperature range of 350–550 K. The results were compared with available experimental and theoretical values in the literature. Local Si-atom dependent activation energies for vacancy hops were calculated using a broken-model and were stored in a database. The migration barrier and prefactors for Si-diffusivity were found to be in reasonable agreement with available modeling results in the literature. Magnetic disorder has a larger effect on the prefactor than on the migration barrier. Prefactor was approximately an order of magnitude and the migration barrier a tenth of an electron-volt higher with magnetic disorder when compared to a fully ferromagnetic ordered state. In addition, the correlation between various have a larger effect on the Si-diffusivity extracted in various temperature range than the magnetic disorder. In the case of vacancy diffusivity, the migration barrier more or less remained constant while the prefactor decreased with increasing Si concentration in the disordered or A2-phase of Fe-Si alloy. Important vacancy-Si/Fe atom exchange processes and their activation barriers were also identified and discuss the effect of energetics on the formation of ordered phases in Fe-Si alloys.

AKSOME↗
DOE OSTI2019

Relationship between Surface Self-Diffusion and Bubble Mobility in Solids: Theory and Atomistic Simulation

Diffusion processes in solids are characterized by complex mechanisms occurring at the atomistic level. The relevant theoretical concepts are sufficiently developed in the case of the diffusion of point defects. At the same time, the diffusion of objects such as clusters of vacancies or nanometer-size cavities in the crystal lattice is described only in terms of continuous-medium approximation. This paper attempts to link the existing theory of diffusion of bubbles in a crystalline matrix with the atomistic description of the process on the basis of the molecular dynamics method. For simplicity, the diffusion of cavities in a two-dimensional Lennard-Jones lattice is considered. A controlled molecular dynamics method is proposed to speed up the direct computation of the diffusion rate of bubbles. Based on the results of simulation, a practical interpretation is given to a parameter of the continuum theory such as the rate of surface self-diffusion. The applicability range of the continuum theory is demonstrated, and principles for the refinement of the theory in the case of minimum-size bubbles are formulated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Diffusion Behavior of Mn and Si Between Liquid Oxide Inclusions and Solid Iron-Based Alloy at 1473 K

In order to clarify the changes in the composition of oxide inclusions in steel, the effect of the metal and oxide composition on the reaction between solid Fe-based alloys and liquid multi-component oxides was investigated using the diffusion couple method at 1473 K. The measured concentration gradients of Mn and Si in the metal indicated that Mn diffused into the metal from the oxide, while the diffusion of Si occurred in the opposite direction. In addition, the MnO content in the oxide decreased with heat treatment time, while the SiO{sub 2} content increased. The compositional changes in both phases indicated that the Mn content in the metal near the interface increased with heat treatment with decreasing MnO content in the oxide. Assuming local equilibrium at the interface, the calculated [Mn]{sup 2}/[Si] ratio at the interface in equilibrium with the oxide increased with increases in the MnO/SiO{sub 2} ratio in the oxide. The difference in the [Mn]{sup 2}/[Si] ratios between the interface and the metal matrix increased, which caused the diffusion of Mn and Si between the multi-component oxide and metal. By measuring the diffusion lengths of Mn and Si in the metal, the chemical diffusion coefficients of Mn and Si were obtained to calculate the composition changes in Mn and Si in the metal. The calculated changes in Mn and Si in the metal agreed with the experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Methane Diffusion and Adsorption in Shale Rocks: A Numerical Study Using the Dusty Gas Model in TOUGH2/EOS7C-ECBM

Gas production from shale gas reservoirs plays a significant role in satisfying increasing energy demands. Compared with conventional sandstone and carbonate reservoirs, shale gas reservoirs are characterized by extremely low porosity, ultra-low permeability and high clay content. Slip flow, diffusion, adsorption and desorption are the primary gas transport processes in shale matrix, while Darcy flow is restricted to fractures. Understanding methane diffusion and adsorption, and gas flow and equilibrium in the low-permeability matrix of shale is crucial for shale formation evaluation and for predicting gas production. Modeling of diffusion in low-permeability shale rocks requires use of the Dusty gas model (DGM) rather than Fick’s law. The DGM is incorporated in the TOUGH2 module EOS7C-ECBM, a modified version of EOS7C that simulates multicomponent gas mixture transport in porous media. Also included in EOS7C-ECBM is the extended Langmuir model for adsorption and desorption of gases. In this study, a column shale model was constructed to simulate methane diffusion and adsorption through shale rocks. The process of binary CH 4 - N 2 diffusion and adsorption was analyzed. A sensitivity study was performed to investigate the effects of pressure, temperature and permeability on diffusion and adsorption in shale rocks. The results show that methane gas diffusion and adsorption in shale is a slow process of dynamic equilibrium, which can be illustrated by the slope of a curve in CH4mass variation. The amount of adsorption increases with the pressure increase at the low pressure, and the mass change by gas diffusion will decrease due to the decrease in the compressibility factor of the gas. With the elevated temperature, the gas molecules move faster and then the greater gas diffusion rates make the process duration shorter. The gas diffusion rate decreases with the permeability decrease, and there is a limit of gas diffusion if the permeability is less than 1.0×10 -15 m 2 . In conclusion, the results can provide insights for a better understanding of methane diffusion and adsorption in the shale rocks so as to optimize gas production performance of shale gas reservoirs.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2019

Thermodynamics of Dual Doping in Quantum Dots

Dual doping is a powerful way to tailor the properties of semiconductor quantum dots (QDs) arising out of host–dopant and dopant–dopant interactions. Nevertheless, it has seldom been explored due to a variety of thermodynamic challenges, such as the differential bonding strength and diffusion constant within the host matrix that integrates with the host in dissimilar ways. This work discusses the challenges involved in administering them within the constraints of one host under similar conditions of temperature, time, and chemical parameters such as solubility and reactivity using CoPt-doped CdS QDs as a model system. In addition, the various forces in play, such as Kirkendall diffusion, solid- and liquid-state diffusion, hard acid soft base interaction with the host, and the effect of lattice strain due to lattice mismatch, are studied to understand the feasibility of the core to doped transformation. Furthermore, these findings suggest a potential approach for manipulating the properties of semiconductors by dual doping engineering.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2015

Formation of Ge nanoparticles in SiO{sub x}N{sub y} by ion implantation and thermal annealing

Germanium nanoparticles embedded within dielectric matrices hold much promise for applications in optoelectronic and electronic devices. Here we investigate the formation of Ge nanoparticles in amorphous SiO{sub 1.67}N{sub 0.14} as a function of implanted atom concentration and thermal annealing temperature. Using x-ray absorption spectroscopy and other complementary techniques, we show Ge nanoparticles exhibit significant finite-size effects such that the coordination number decreases and structural disorder increases as the nanoparticle size decreases. While the composition of SiO{sub 1.67}N{sub 0.14} is close to that of SiO{sub 2}, we demonstrate that the addition of this small fraction of N yields a much reduced nanoparticle size relative to those formed in SiO{sub 2} under comparable implantation and annealing conditions. We attribute this difference to an increase in an atomic density and a much reduced diffusivity of Ge in the oxynitride matrix. These results demonstrate the potential for tailoring Ge nanoparticle sizes and structural properties in the SiO{sub x}N{sub y} matrices by controlling the oxynitride stoichiometry.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Chloride diffusivity, chloride threshold, and corrosion initiation in reinforced alkali-activated mortars: Role of calcium, alkali, and silicate content

The aim of this study is to investigate systematically the chloride diffusivity and chloride threshold of a wide range of calcium-rich and fly ash-dominated alkali-activated samples in light of their compositional differences. To this end, the effects of various fly ash (FA)-to-slag ratios, of alkali concentrations and of silicate content in the activator were investigated. The electrochemical aspects of the passive samples were also assessed. Results show the prominent role of calcium in the matrix to reduce the chloride diffusivity. While higher alkali concentration increased the porosity and chloride diffusivities in general, lower modulus ratios provided considerably better performance in the FA-dominated samples. Chloride threshold values range between 0.19 (wt% binder mass) for calcium-rich mortars fabricated at low levels of alkalinities and 0.69 for FA-dominated mortars fabricated with highly alkaline activators. Half-cell potential and polarization resistance of alkali-activated samples were in general lower than their Portland cement counterparts.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Constraining the Reaction Rate of Criegee Intermediates with Carboxylic Acids during the Multiphase Ozonolysis of Aerosolized Alkenes

Criegee intermediates (CIs) are central to the ozonolysis of unsaturated organic compounds, controlling the formation of small fragmentation reaction products and higher molecular weight oligomeric compounds through competing unimolecular and bimolecular pathways. In particular, the reaction of thermalized CI with carboxylic acids (RCOOH) to produce oligomeric α-acyloxyalkyl hydroperoxides (AAHPs) is of interest to the chemical transformation of organic compounds in the atmosphere. In the gas phase, the CI + RCOOH reaction proceeds at supercollisional rates, consistent with a barrierless reaction mechanism. However, the rate of this reaction in condensed phases remains uncertain. Here, we report the results of aerosol flow tube studies of the ozonolysis of submicron organic aerosol containing mixtures of an internal n-alkene (Z-9-tricosene, tri) and a saturated acid (2-hexyldecanoic acid, HDA). The decay of the reactants and formation of reaction products are monitored using atmospheric pressure chemical ionization mass spectrometry (APCI–MS). A 6-fold decrease in the reactive uptake of ozone is observed with increasing acid concentration, consistent with a previous study using alcohol additives. The decay of HDA with ozone exposure shows that HDA scavenges CIs at rates comparable to competing unimolecular pathways involving CIs. The decay kinetics of HDA are analyzed using a simple kinetic model in order to constrain the ratio of the CI + RCOOH and unimolecular CI rate constants. Empirical fitting to this model, combined with theoretical estimates of the unimolecular loss rate, yields a CI + RCOOH rate constant of 1.85 ± 0.27 × 10 –19 cm 3 molec –1 s –1 , which is six orders of magnitude slower than the expected diffusion limit in a liquid organic matrix. Stochastic kinetic simulations are conducted with different values of the CI + RCOOH rate constant to validate this result and show that the diffusion-limited rate is indeed much too rapid to reproduce the experimental results.

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

Determination of the effective diffusivity of water in a poly (methyl methacrylate) membrane containing carbon nanotubes using kinetic Monte Carlo simulations

A kinetic Monte Carlo (kMC) simulation algorithm is developed for computing the effective diffusivity of water molecules in a poly(methyl methacrylate) (PMMA) matrix containing carbon nanotubes (CNTs) at several loadings. The simulations are conducted on a cubic lattice to the bonds of which rate constants are assigned governing the elementary jump events of water molecules from one lattice site to another. Lattice sites belonging to PMMA domains of the membrane are assigned different rates than lattice sites belonging to CNT domains. Values of these two rate constants are extracted from available numerical data for water diffusivity within a PMMA matrix and a CNT pre-computed on the basis of independent atomistic molecular dynamics simulations, which show that water diffusivity in CNTs is 3 orders of magnitude faster than in PMMA. Our discrete-space, continuum-time kMC simulation results for several PMMA-CNT nanocomposite membranes (characterized by different values of CNT length L and diameter D and by different loadings of the matrix in CNTs) demonstrate that the overall or effective diffusivity, D{sub eff}, of water in the entire polymeric membrane is of the same order of magnitude as its diffusivity in PMMA domains and increases only linearly with the concentration C (vol. %) in nanotubes. For a constant value of the concentration C, D{sub eff} is found to vary practically linearly also with the CNT aspect ratio L/D. The kMC data allow us to propose a simple bilinear expression for D{sub eff} as a function of C and L/D that can describe the numerical data for water mobility in the membrane extremely accurately. Additional simulations with two different CNT configurations (completely random versus aligned) show that CNT orientation in the polymeric matrix has only a minor effect on D{sub eff} (as long as CNTs do not fully penetrate the membrane). We have also extensively analyzed and quantified sublinear (anomalous) diffusive phenomena over small to moderate times and correlated them with the time needed for penetrant water molecules to explore the available large, fast-diffusing CNT pores before Fickian diffusion is reached.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗