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 181 records · Page 10

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↗

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↗

Unexpected Diffusion Anisotropy of Carbon Dioxide in the Metal–Organic Framework Zn 2 (dobpdc)

Metal–organic frameworks are promising materials for energy-efficient gas separations, but little is known about the diffusion of adsorbates in materials featuring one-dimensional porosity at the nanoscale. An understanding of the interplay between framework structure and gas diffusion is crucial for the practical application of these materials as adsorbents or in mixed-matrix membranes, since the rate of gas diffusion within the adsorbent pores impacts the required size (and therefore cost) of the adsorbent column or membrane. In this paper, we investigate the diffusion of CO 2 within the pores of Zn 2 (dobpdc) (dobpdc 4– = 4,4'-dioxidobiphenyl-3,3'-dicarboxylate) using pulsed field gradient (PFG) nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations. The residual chemical shift anisotropy for pore-confined CO 2 allows PFG NMR measurements of self-diffusion in different crystallographic directions, and our analysis of the entire NMR line shape as a function of the applied field gradient provides a precise determination of the self-diffusion coefficients. In addition to observing CO 2 diffusion through the channels parallel to the crystallographic c axis (self-diffusion coefficient D ∥ = (5.8 ± 0.1) × 10 –9 m 2 s –1 at a pressure of 625 mbar CO 2 ), we unexpectedly find that CO 2 is also able to diffuse between the hexagonal channels in the crystallographic ab plane (D ⊥ = (1.9 ± 0.2) × 10 –10 m 2 s –1 ), despite the walls of these channels appearing impermeable by single-crystal X-ray crystallography and flexible lattice MD simulations. Finally, observation of such unexpected diffusion in the ab plane suggests the presence of defects that enable effective multidimensional CO 2 transport in a metal–organic framework with nominally one-dimensional porosity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation of fission gas diffusion in UO{sub 2} by vacancy clusters

The low solubility of fission gas in nuclear fuel leads to inter- and intra-granular bubbles. These bubbles affect performance of the fuel rod and increase the fuel temperature through the increase of porosity in the fuel pellets and decreased thermal conductivity across the pellet-clad gap as gas is released into the plenum. The nucleation, growth, and inter-connectivity of the fission gas bubbles ultimately depend on the diffusivity of gas atoms through the fuel matrix. Using mechanisms and data previously established from density functional theory (DFT) and empirical potential calculations, diffusion and reaction models have been derived for xenon atoms, uranium atoms, and vacancies in the substitutional and interstitial locations. Although similar to previous models, the diffusion and reaction rates have been re-derived in a more robust manner. Previous studies have under-predicted the experimentally observed diffusivity of xenon atoms in UO{sub 2}, however it has been previously shown that by including xenon atoms bound to two vacancies, the calculated diffusivity approaches the measured data. The derived model for species diffusion and inter-species reactions has been implemented in the MARMOT phase-field code in a modular structure so as to capture the enhanced diffusion of xenon due to bound gas + vacancy-cluster movement. The MARMOT code allows for the diffusion, reaction, and source rates to be solved simultaneously, resulting in steady state or transient concentrations of each species, and ultimately the diffusion behavior of xenon atoms. By including the diffusion of clusters of different combinations of xenon and vacancies, the calculated diffusivity more closely captures the measured diffusivity of xenon when using preliminary parameters. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Exponential time differencing for the tracer equations appearing in primitive equation ocean models

The tracer equations are part of the primitive equations used in ocean modeling and describe the transport of tracers, such as temperature, salinity or chemicals, in the ocean. Depending on the number of tracers considered, several equations may be added to and coupled to the dynamics system. In many relevant situations, the time-step requirements of explicit methods imposed by the transport and mixing in the vertical direction are more restrictive than those for the horizontal, and this may cause the need to use very small time steps if a fully explicit method is employed. To overcome this issue, we propose an exponential time differencing (ETD) solver where the vertical terms (transport and diffusion) are treated with a matrix exponential, whereas the horizontal terms are dealt with in an explicit way. In this work, we investigate numerically the computational speed-ups that can be obtained over other semi-implicit methods, and we analyze the advantages of the method in the case of multiple tracers.

42 ENGINEERING↗

Enabling fast charging of lithium-ion batteries through secondary- /dual- pore network: Part I - Analytical diffusion model

Battery performance is strongly correlated with electrode microstructural properties. Enabling fast charging of lithium-ion batteries requires improved through-plane ionic diffusion that can be achieved through, among other strategies, structured electrodes with a secondary- or dual-pore network (SPN). In this work, an analytical model investigates the impact of such an SPN on ionic diffusion with a composite electrode, considering various pore-channel geometries and comparing to standard electrodes with identical gravimetric- and volumetric-specific theoretical capacities. Relevant SPN design parameters and tortuosity coefficients are identified according to three optimization objectives that aim to balance the improved overall through-plane diffusion, thanks to the coarse aligned channels, and degraded in-plane diffusion because of the porous matrix densification required to maintain gravimetric- and volumetric-specific theoretical capacities. The model indicates that a relatively low amount of SPN is required and that electrodes with high through-plane tortuosity and low in-plane tortuosity benefit most from such architecture.

25 ENERGY STORAGE↗

Chemical reactions, porosity, and microfracturing in shale during weathering: The effect of erosion rate

The rate of chemical weathering has been observed to increase with the rate of physical erosion in published comparisons of many catchments, but the mechanisms that couple these processes are not well understood. We investigated this question in our research by examining the chemical weathering and porosity profiles from catchments developed on marine shale located in Pennsylvania, USA (Susquehanna Shale Hills Critical Zone Observatory, SSHCZO); California, USA (Eel River Critical Zone Observatory, ERCZO); and Taiwan (Fushan Experimental Forest). The protolith compositions, protolith porosities, and the depths of regolith at these sites are roughly similar while the catchments are characterized by large differences in erosion rate (1–3 mm yr ₋1 in Fushan &laquo 0.2–0.4 mm yr ₋1 in ERCZO &laquo 0.01–0.025 mm yr ₋1 in SSHCZO). The natural experiment did not totally isolate erosion as a variable: mean annual precipitation varied along the erosion gradient (4.2 m yr ₋1 in Fushan > 1.9 m yr ₋1 in ERCZO > 1.1 m yr ₋1 in SSHCZO), so the fastest eroding site experiences nearly twice the mean annual temperature of the other two. Even though erosion rates varied by about 100×, the depth of pyrite and carbonate depletion (defined here as regolith thickness) is roughly the same, consistent with chemical weathering of those minerals keeping up with erosion at the three sites. These minerals were always observed to be the deepest to react, and they reacted until 100% depletion. In two of three of the catchments where borehole observations were available for ridges, these minerals weathered across narrow reaction fronts. On the other hand, for the rock-forming clay mineral chlorite, the depth interval of weathering was wide and the extent of depletion observed at the land surface decreased with increasing erosion/precipitation. Thus, chemical weathering of the clay did not keep pace with erosion rate. But perhaps the biggest difference among the shales is that in the fast-eroding sites, microfractures account for 30–60% of the total porosity while in the slow-eroding shale, dissolution could be directly related to secondary porosity. We argue that the microfractures increase the influx of oxygen at depth and decrease the size of diffusion-limited internal domains of matrix, accelerating weathering of pyrite and carbonate under high erosion-rate conditions. Thus, microfracturing is a process that can couple physical erosion and chemical weathering in shales.

58 GEOSCIENCES↗

Investigation of the ozone formation reaction pathway: Comparisons of full configuration interaction quantum Monte Carlo and fixed-node diffusion Monte Carlo with contracted and uncontracted MRCI

The association/dissociation reaction path for ozone (O 2 + O ↔ O 3 ) is notoriously difficult to describe accurately using ab initio electronic structure theory, due to the importance of both strong and dynamic electron correlations. Experimentally, spectroscopic studies of the highest lying recorded vibrational states combined with the observed negative temperature dependence of the kinetics of oxygen isotope exchange reactions confirm that the reaction is barrierless, consistent with the latest potential energy surfaces. Previously reported potentials based on Davidson-corrected internally contracted multireference configuration interaction (MRCI) suffer from a spurious reef feature in the entrance channel even when extrapolated towards the complete basis set limit. In this work, we report an analysis of comparisons between a variety of electronic structure methods including internally contracted and uncontracted MRCI (with and without Davidson corrections), as well as full configuration interaction quantum Monte Carlo, fixed-node diffusion Monte Carlo, and density matrix renormalization group.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hybrid quantum-classical simulations of magic angle spinning dynamic nuclear polarization in very large spin systems

Solid-state nuclear magnetic resonance can be enhanced using unpaired electron spins with a method known as dynamic nuclear polarization (DNP). Fundamentally, DNP involves ensembles of thousands of spins, a scale that is difficult to match computationally. This scale prevents us from gaining a complete understanding of the spin dynamics and applying simulations to design sample formulations. We recently developed an ab initio model capable of calculating DNP enhancements in systems of up to ~1000 nuclei; however, this scale is insufficient to accurately simulate the dependence of DNP enhancements on radical concentration or magic angle spinning (MAS) frequency. We build on this work by using ab initio simulations to train a hybrid model that makes use of a rate matrix to treat nuclear spin diffusion. We show that this model can reproduce the MAS rate and concentration dependence of DNP enhancements and build-up time constants. Finally, we then apply it to predict the DNP enhancements in core–shell metal-organic-framework nanoparticles and reveal new insights into the composition of the particles’ shells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Superdiffusive transport of energy in one-dimensional metals

Metals in one spatial dimension are described at the lowest energy scales by the Luttinger liquid theory. It is well understood that this free theory, and even interacting integrable models, can support ballistic transport of conserved quantities including energy. In contrast, realistic one-dimensional metals, even without disorder, contain integrability-breaking interactions that are expected to lead to thermalization and conventional diffusive linear response. We argue that the expansion of energy when such a nonintegrable Luttinger liquid is locally heated above its ground state shows superdiffusive behavior (i.e., spreading of energy that is intermediate between diffusion and ballistic propagation), by combining an analytical anomalous diffusion model with numerical matrix-product–state calculations on a specific perturbed spinless fermion chain. Different metals will have different scaling exponents and shapes in their energy spreading, but the superdiffusive behavior is stable and should be visible in time-resolved experiments.

Bulchandani, Vir B.↗

The mechanism of action of pepR, a viral-derived peptide, against Staphylococcus aureus biofilms

Abstract Objectives To investigate the mechanism of action at the molecular level of pepR, a multifunctional peptide derived from the Dengue virus capsid protein, against Staphylococcus aureus biofilms. Methods Biofilm mass, metabolic activity and viability were quantified using conventional microbiology techniques, while fluorescence imaging methods, including a real-time calcein release assay, were employed to investigate the kinetics of pepR activity at different biofilm depths. Results Using flow cytometry-based assays, we showed that pepR is able to prevent staphylococcal biofilm formation due to a fast killing of planktonic bacteria, which in turn resulted from a peptide-induced increase in the permeability of the bacterial membranes. The activity of pepR against pre-formed biofilms was evaluated through the application of a quantitative live/dead confocal laser scanning microscopy (CLSM) assay. The results show that the bactericidal activity of pepR on pre-formed biofilms is dose and depth dependent. A CLSM-based assay of calcein release from biofilm-embedded bacteria was further developed to indirectly assess the diffusion and membrane permeabilization properties of pepR throughout the biofilm. A slower diffusion and delayed activity of the peptide at deeper layers of the biofilm were quantified. Conclusions Overall, our results show that the activity of pepR on pre-formed biofilms is controlled by its diffusion along the biofilm layers, an effect that can be counteracted by an additional administration of peptide. Our study sheds new light on the antibiofilm mechanism of action of antimicrobial peptides, particularly the importance of their diffusion properties through the biofilm matrix on their activity.

Pinto, Sandra N.↗

Secondary Pore Network Design Optimization Analytical Diffusion Model for Lithium Ion Battery

Fast ion transport along the electrode thickness is required for Lithium ion battery (LIB) to sustain continuous high-rate discharge or charge currents, which is a desired feature for the fast charging of electric vehicle batteries. Advance in LIB electrode manufacturing allows to control, to a certain extent, the electrode microstructure architecture in order to improve the macroscopic performances of these devices. For instance, several groups have experimentally tested architectures tailored specifically to increase the through plane ionic diffusion in regards for fast charge application by adding channels or groove-like large pores oriented along the electrode thickness. These voids do not replace the existing pores of the electrode microstructure but consist in a second, coarser, pore network, called the secondary pore network (SPN). The model quantifies the impact of SPN on ionic diffusion with a composite electrode through a simple analytical approach, considering various pore channel geometries and comparing to standard electrodes with identical gravimetric and identical volumetric specific theoretical capacity. The goal is to find optimal SPN design parameters according to three optimization objectives aiming to balance the improved overall through-plane, thanks to the coarse aligned channels, and the degraded in-plane diffusion, due to the porous matrix densification required to maintain gravimetric and volumetric specific theoretical capacity. The analytical model must be considered as a pre-screening tool to quickly investigate large design space, that can be refined later with more complex electrochemical CPU-expensive numerical models.

Usseglio Viretta, Francois↗

Calculation of vaporization rates assuming various rate determining steps: Application to the resistojet

The various steps that could control the vaporization rate of a material are discussed. These steps include the actual vaporization, flow rate of matrix gas, chemical reaction, gas diffusion, and solid state diffusion. The applicable equations have been collected from diverse appropriate sources, and their use is explained. Rate equations are derived for conditions where more than one step is rate controlling. Calculations are made for two model materials: rhenium which vaporizes congruently, and tantalum carbide which vaporizes incongruently. The case of vaporization under thermal gradient conditions is also treated. The existence of a thermal gradient in the resistojet means that the vaporization rate of a material may be only one thousandth of that predicted under isothermal conditions. Calculations show that rhenium might have a 100,000 hr lifetime at temperature in a 2500 C resistojet. Tantalum carbide would have a life of only 660 sec under similar conditions.

May, C. E.↗

Problems and possible solutions involved in hard target calibration of coherent Doppler lidar

The field of radiometry is surveyed with reference to coherent light, since the surfaces to be considered here are irradiated with coherent, polarized light. Definitions of some concepts are modified. In light of these modifications calibration problems are analyzed and solutions suggested. The most important task was to develop hard targets exhibiting minimal specular reflection (mirror-like and retroreflection) and following closely a Lambertian scattering curve. Bistable reflectometer experiments and integrating sphere measurements are used to physically characterize the targets. Optical and electron microscopy are used to physically characterize the targets. Since the BRDF matrix is not measured, this capability must be developed preferably for both bistatic and monostatic reflectometer measurements. The equipment is expensive and not yet developed. Calibration can be simplified if it can be proved that the BRDF matrix is diagonal for a diffuse scatterer.

Anderson, Richard↗

Ignition and combustion characteristics of metallized propellants

Research designed to develop detailed knowledge of the secondary atomization and ignition characteristics of aluminum slurry propellants was started. These processes are studied because they are the controlling factors limiting the combustion efficiency of aluminum slurry propellants in rocket applications. A burner and spray rig system allowing the study of individual slurry droplets having diameters from about 10 to 100 microns was designed and fabricated. The burner generates a near uniform high temperature environment from the merging of 72 small laminar diffusion flames above a honeycomb matrix. This design permits essentially adiabatic operation over a wide range of stoichiometries without danger of flashback. A single particle sizing system and velocimeter also were designed and assembled. Light scattered from a focused laser beam is related to the particle (droplet) size, while the particle velocity is determined by its transit time through the focal volume. Light from the combustion of aluminum is also sensed to determine if ignition was achieved. These size and velocity measurements will allow the determination of disruption and ignition times as functions of drop sizes and ambient conditions.

Turns, S. R.↗

Effect of Specimen Thickness on Mechanical Behavior of SiC/SiC Composites

Potential composite applications in aerospace and transportation application systems have different thickness requirements. For example, space applications such as nozzle ramps or heat exchangers use very thin (less than 1 mm) structures whereas turbine blades need very thick parts greater than or equal to cm). There has been little investigation into the effect of thickness on stress-strain behavior or elevated temperature tensile properties controlled by oxidation. In this study, composites consisting of woven Hi-NicalonTM fibers, a carbon interphase, and CVI Sic matrix were fabricated with different numbers of plies to provide variable thickness. The composites ranged from a single ply (approximately 0.4 mm) to thirty-six plies (approximately 1 cm). Tensile tests were performed at room temperature with acoustic emission used to monitor matrix crack behavior. Elevated temperature tensile stress-rupture tests were performed in air. Considerably different room and elevated temperature tensile behavior was observed that will be discussed with respect to the effect of thickness on matrix crack formation, matrix crack growth and oxidation diffusion kinetics.

Morscher, Gregory N.↗

Microstructural Observations on Advanced Multilayer Protective Coatings for Nuclear Thermal Propulsion

The present research discusses the microstructures of a new multilayered coating concept proposed to solve the “midrange corrosion” problem observed in NbC and ZrC-coated graphite (Gr) nuclear fuels during the NERVA/Rover programs. The concept envisions designing a compliant multilayered coatings architecture to accommodate the thermal expansion mismatch strains between the ZrC outer coating and the Gr/(U,Zr)C fuel matrix along with using a diffusion barrier to minimize carbon and hydrogen diffusion, respectively. Multilayered coatings were deposited on Gr substrates by chemical vapor deposition (CVD). Graphite disks, as well as the inner channels of a 19-hole hexagonal rod 152.4 mm long, were coated by CVD in proof-of concept process demonstration studies. Coated disk specimens were thermally-cycled between ambient and 1900 K for as long as 20 cycles with a 1h hold time at the temperature using both stepped and non-interrupted thermal cycling methods. Detailed cross-sectional microstructural analyses were conducted on the thermally-cycled coated disks. These results showed that the coatings were intact after 20 thermal cycles at 1900 K although the bonding of the overlay Mo coating appeared to be influenced by the oxygen content in the layer. Seven transverse sections were cut along the length from the coated hexagonal rod for microstructural studies. All the nineteen channels were observed to have been coated in all the sections. This paper reports the results of detailed microstructural analyses, compositional line scans and x-ray dot maps conducted on three coated channels from a mid-section of the rod. The measured values of carbon (C) were observed to be exceptionally high in all the specimens, which were attributed to hydrocarbon contamination of the specimen surfaces. It was observed that the oxygen (O) levels were high in all the coating layers even after approximately correcting for the observed high values of C. The source of O contamination of the CVD reactor system was attributed to either the presence of oxychloride impurities in the precursor salts or a pin hole leak in the coating system. It is concluded that high values of O result in poor bonding of the Mo overlay to the ZrC layer. These trials have established that it is possible to deposit multilayered coatings by CVD in 19-hole hexagonal rods 152.4 mm long similar to the NERVA fuel rod design although the quality of the coatings and the uniformity of their thicknesses require further improvement.

Nuclear thermal propulsion, multilayer coatings, m↗

Creating Hydride Rim Structures in Zircaloy 4 Cladding Tubes to Mimic High Burnup

During their long service in commercial light water reactors, zirconium alloy cladding tubes are subject to oxidation and subsequent hydrogen pickup which can degrade the materials mechanical properties [1]. These degradations can lead to premature cladding ruptures during reactor power transients [2] as well as negatively impact the claddings post quench ductility following high temperature steam oxidation in thermal hydraulic transients [3]. Hydrogen has low solid solubility in zirconium and usually precipitates as brittle zirconium hydride platelets in the cladding metal matrix. Due to temperature driven diffusion (Soret effect) and the grain texture of stress relieved zirconium metal tubes, these hydrides precipitate with a circumferential orientation near the outer surface of the cladding tubes forming a hydride rim which can serve as a site for crack initiation during circumferential ‘hoop’ loading of the cladding [4rim formation][5]. In order to better understand the effect of hydride rims on the fracture behavior of stress relieved zirconium alloy cladding tubes a procedure to produce these rim-like structures is being developed at Idaho National Laboratory (INL).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗