SEARCH · Engineering Papers

Results for “reactive transport”

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.

SourceResearch connectionFrom yearThrough year

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

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Reaction Kernel Structure of a Slot Jet Diffusion Flame in Microgravity

Diffusion flame stabilization in normal earth gravity (1 g) has long been a fundamental research subject in combustion. Local flame-flow phenomena, including heat and species transport and chemical reactions, around the flame base in the vicinity of condensed surfaces control flame stabilization and fire spreading processes. Therefore, gravity plays an important role in the subject topic because buoyancy induces flow in the flame zone, thus increasing the convective (and diffusive) oxygen transport into the flame zone and, in turn, reaction rates. Recent computations show that a peak reactivity (heat-release or oxygen-consumption rate) spot, or reaction kernel, is formed in the flame base by back-diffusion and reactions of radical species in the incoming oxygen-abundant flow at relatively low temperatures (about 1550 K). Quasi-linear correlations were found between the peak heat-release or oxygen-consumption rate and the velocity at the reaction kernel for cases including both jet and flat-plate diffusion flames in airflow. The reaction kernel provides a stationary ignition source to incoming reactants, sustains combustion, and thus stabilizes the trailing diffusion flame. In a quiescent microgravity environment, no buoyancy-induced flow exits and thus purely diffusive transport controls the reaction rates. Flame stabilization mechanisms in such purely diffusion-controlled regime remain largely unstudied. Therefore, it will be a rigorous test for the reaction kernel correlation if it can be extended toward zero velocity conditions in the purely diffusion-controlled regime. The objectives of this study are to reveal the structure of the flame-stabilizing region of a two-dimensional (2D) laminar jet diffusion flame in microgravity and develop a unified diffusion flame stabilization mechanism. This paper reports the recent progress in the computation and experiment performed in microgravity.

Takahashi, F.↗

NASA NTRS2003

Quadrature Moments Method for the Simulation of Turbulent Reactive Flows

A sub-filter model for reactive flows, namely the DQMOM model, was formulated for Large Eddy Simulation (LES) using the filtered mass density function. Transport equations required to determine the location and size of the delta-peaks were then formulated for a 2-peak decomposition of the FDF. The DQMOM scheme was implemented in an existing structured-grid LES solver. Simulations of scalar shear layer using an experimental configuration showed that the first and second moments of both reactive and inert scalars are in good agreement with a conventional Lagrangian scheme that evolves the same FDF. Comparisons with LES simulations performed using laminar chemistry assumption for the reactive scalar show that the new method provides vast improvements at minimal computational cost. Currently, the DQMOM model is being implemented for use with the progress variable/mixture fraction model of Pierce. Comparisons with experimental results and LES simulations using a single-environment for the progress-variable are planned. Future studies will aim at understanding the effect of increase in environments on predictions.

Raman, Venkatramanan↗

DOE OSTI2021

Characterizing the Influence of Organic Polymers on the Specific Reactivity of Particulate Remedial Amendments

Commercially available particulate amendments demonstrate high reactivity for effective treatment of water soluble organic and inorganic contaminants in laboratory studies; however, transport of these particles is constrained in the subsurface. In many field applications, particulate amendments are mixed with organic polymers to enhance mobility for direct push applications or stabilize suspensions for high mass loadings. As such, the interactions between particulate amendments, organic polymers and contaminant species need to be systematically investigated to properly understand mechanistic processes that facilitate predictive performance metrics for specific applications in situ . In this study, batch experiments were conducted to quantify the effects of organic polymers (xanthan gum, guar gum, and sodium alginate), polymer concentration (800 and 4,000 mg/L), and aging (up to 28 days) on chromate treatment rate and capacity by two classes of amendments: reductants [granular zero-valent iron (gZVI), micron-ZVI (mZVI), sulfur modified iron (SMI)], and an adsorbent (bismuth sub-nitrate). When particulate amendments were suspended in polymer solutions, reductants retained between 84–100% of the amendment treatment capacity. Conversely, the adsorbent maintained 63–97% relative treatment capacity of the no-polymer control. Polymer solutions had a more pronounced impact on the rate of chromate removal; first order rates of chemical reduction decreased by as much as 70% and adsorption by up to 81% relative to the no-polymer controls. Polymer–amendment aging experiments also showed decreased Cr(VI) treatment capacity; reductants decreased by as much as 24% and adsorption decreased by as much as 44% after 28 days of incubation. While polymer suspensions are needed to aid the injection of particulate amendments into the subsurface, the results from this study indicate potential losses of treatment capacity and a decrease in the rate of remedial performance due to the physical and chemical interactions between polymer suspensions and reactive particulate amendments. Simple batch systems provide baseline characterization of tripartite interactions for the removal of Cr(VI). Additional work is needed to quantify the full impact of polymers on remedial outcomes under site relevant conditions at field scale.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1970

Polymer degradation rate control of hybrid rocket combustion

Polymer degradation to small fragments is treated as a rate controlling step in hybrid rocket combustion. Both numerical and approximate analytical solutions of the complete energy and polymer chain bond conservation equations for the condensed phase are obtained. Comparison with inert atmosphere data is very good. It is found that the intersect of curves of pyrolysis rate versus interface temperature for hybrid combustors, with the thermal degradation theory, falls at a pyrolysis rate very close to that for which a pressure dependence begins to be observable. Since simple thermal degradation cannot give sufficient depolymerization at higher pyrolysis rates, it is suggested that oxidative catalysis of the process occurs at the surface, giving a first order dependence on reactive species concentration at the wall. Estimates of the ratio of this activation energy and interface temperature are in agreement with best fit procedures for hybrid combustion data. Requisite active species concentrations and flux are shown to be compatible with turbulent transport. Pressure dependence of hybrid rocket fuel regression rate is thus shown to be describable in a consistent manner in terms of reactive species catalysis of polymer degradation.

Stickler, D. B.↗

DOE OSTI2021

Role of surface diffusion in formation of unique reactivity for graphite oxidation: Time-resolved measurements in a pulsed diffusion reactor

Quantification of oxidation kinetics is essential to develop graphitic materials for diverse applications: from refractories found in gas-cooled nuclear reactors to catalysts needed for chemical manufacturing. In this work, using well-defined highly oriented pyrolytic graphite, low-pressure isotopic transient experiments combined with controlled annealing periods, we resolve the role of surface diffusion and quantify oxidation kinetics with nanomole-precision. We observe an unexpected increase in reactivity following annealing which is explained by the role of surface diffusion increasing the probability for trapping mobile oxygen at more reactive edge sites. Here, the locus of adsorption and spillover to the basal plane is distinct from the trapping location creating a more active oxygen species. Isotopic products reflect the population dynamics of oxygen added at the edge and surface diffusion that relocates basal plane oxygen to more reactive edge sites. Since this process proceeds in parallel with direct oxidation reactions, it is not likely to be observed using steady-state or conventional ‘bulk’ characterization techniques. Our unique time-resolved non-equilibrium measurement in a well-defined transport regime, enables observation of three distinct behaviors: short-term deactivation due to the balance of rates in oxygen supply/product formation, reactivity increases due to surface diffusion and longer-term reactivity increase with oxygen accumulation.

36 MATERIALS SCIENCE↗

NASA NTRS1995

The fate of atmospheric phosgene and the stratospheric chlorine loadings of its parent compounds: CCl4, C2Cl4, C2HCL3, CH3CCl3, and CHCl3

A study of the tropospheric and stratospheric cycles of phosgene is carried out to determine its fate and ultimate role in controlling the ozone depletion potentials of its parent compounds. Tropospheric phosgene is produced from the OH-initiated oxidation of C2Cl4, CH3CCl3, CHCl3, and C2HCl3. Simulations using a two-dimensional model indicate that these processes produce about 90 pptv/yr of tropospheric phosgene with an average concentration of about 18 pptv, in reasonable agreement with observations. We estimate a residence time of about 70 days for tropospheric phosgene, with the vast majority being removed by hydrolysis in cloudwater. Only about 0.4% of the phosgene produced in the troposphere avoids wet removal and is transported to the stratosphere, where its chlorine can be released to participate in the catalytic destruction of ozone. Stratospheric phosgene is produced from the photochemical degradation of CCl4, C2Cl4, CHCl3, and CH3CCl3 and is removed by photolysis and downward transport to the troposphere. Model calculations, in good agreement with observations, indicate that these processes produce a peak stratospheric concentration of about 25-30 pptv at an altitude of about 25 km. In contrast to tropospheric phosgene, stratospheric phosgene is found to have a lifetime against photochemical removal of the order of years. As a result, a significant portion of the phosgene that is produced in the stratosphere is ultimately returned to the troposphere, where it is rapidly removed by clouds. This phenomenon effectively decreases the amount of reactive chlorine injected into the stratosphere and available for ozone depletion from phosgene's parent compounds. A similar phenomenon due to the downward transport of stratospheric COFCl produced from CFC-11 is estimated to cause a 7% decrease in the amount of reactive chlorine injected into the stratosphere from this compound. Our results are potentially sensitive to a variety of parameters, most notably the rate of reaction of phosgene with sulfate aerosols. However, on the basis of the observed vertical distribution of COCl2, we estimate that the reaction of COCl2 with sulfate aerosol most likely has a gamma less than 5 x 10(exp -5) and, as a result, has a negligible impact on the stratospheric chlorine loadings of the phosgene parent compounds.

Kindler, T. P.↗

DOE OSTI2023

Asymptotic hydrogen redistribution analysis in yttrium-hydride-moderated heat-pipe-cooled microreactors using DireWolf

Yttrium hydride (YH x ) is one of the materials being considered for moderating thermal and epithermal nuclear microreactors. One potential issue with YH x use is that the hydrogen redistributes in the hydride when thermal and concentration gradients are present. This hydrogen redistribution leads to spatial gradients in the hydrogen concentration, thus affecting neutron transport in the reactor. Here, by building upon observations in prior works, this paper aims to gain a better understanding of the reactivity feedback associated with such hydrogen redistributions. In particular, we wish to understand the sign (+/–) of the hydrogen redistribution neutronic feedback, its order of magnitude, and its underlying physical causes. To achieve this goal, the DireWolf multiphysics software driver was used to solve the coupled radiation transport, heat transfer, heat pipe two-phase flow, and hydrogen redistribution equations for the Simplified Microreactor Benchmark Assessment (SiMBA) problem, a full-core microreactor numerical benchmark developed at Idaho National Laboratory.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Quantifying the effects of artificial aging on the ignition and self-propagating reactions of Ni(V)/Al multilayers

Bimetallic, reactive multilayers are uniformly structured materials composed of alternating nanoscale layers that may be ignited to produce self-propagating intermetallic-formation reactions. When reactive multilayers age, there is a change in local composition and loss of stored chemical energy due to mass transport and rearrangement at the interfaces. Here, to quantify the long-term reliability of commercial Ni(V)/Al multilayers, the effects of accelerated aging on both ignition sensitivity and self-propagating reactions have been examined. Thermally aged samples were characterized using transmission electron microscopy, differential scanning calorimetry, and laser ignition combined with high-speed videography. The analytically quantifiable nature of both continuous wave laser ignition and reactive wave propagation enabled calculations of Arrhenius rate constants for as-received and various heat-treated multilayers. With heat treatment, there is a change in the intermixed thickness and interfacial chemistry that decreases the activation energy for point ignition but increases it for self-propagating reactions. This finding implies that with increased thermal aging, the reaction becomes easier to facilitate in the solid state but harder in the liquid phase.

Energetic material↗

DOE OSTI2020

Influence of acid leaching surface treatment on indentation cracking of soda lime silicate glass

Past work has shown that water or acid soaking treatments can increase the mechanical strength of soda lime silicate (SLS) glasses. In this work, we show that acid leaching treatments result in an increase in the apparent crack resistance of the acid-leached surface of SLS glass during indentation. Vickers indentation tests in controlled environments show a humidity dependence of radial cracking, suggesting that the transport of water through the leached layer plays a critical role in the propagation of cracks to the glass surface. Here, molecular dynamics simulations with reactive force fields indicate that the leached surface layer can undergo pressure-induced mechanochemical reactions during indentation, which increases the bridging oxygen connectivity in the silica network of the leached layer. Such structural changes can hinder transport of water molecules from the environment to the subsurface crack tip. Based on experimental observations and simulation results, a new hypothesis is proposed that mechanochemical restructuring in the leached layer in response to the applied load may lower the transport kinetics of molecular water to critical flaws, resulting in an apparent enhancement in the crack resistance of the acid-leached surface of SLS glass.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Electrical Interrogation of Thickness‐Dependent Multiferroic Phase Transitions in the 2D Antiferromagnetic Semiconductor NiI 2

Abstract 2D magnetic materials hold promise for quantum and spintronic applications. 2D antiferromagnetic materials are of particular interest due to their relative insensitivity to external magnetic fields and higher switching speeds compared to 2D ferromagnets. However, their lack of macroscopic magnetization impedes detection and control of antiferromagnetic order, thus motivating magneto‐electrical measurements for these purposes. Additionally, many 2D magnetic materials are ambient‐reactive and electrically insulating or highly resistive below their magnetic ordering temperatures, which imposes severe constraints on electronic device fabrication and characterization. Herein, these issues are overcome via a fabrication protocol that achieves electrically conductive devices from the ambient‐reactive 2D antiferromagnetic semiconductor NiI 2 . The resulting gate‐tunable transistors show band‐like electronic transport below the antiferromagnetic and multiferroic transition temperatures of NiI 2 , revealing a Hall mobility of 15 cm 2 V −1 s −1 at 1.7 K. These devices also allow direct electrical probing of the thickness‐dependent multiferroic phase transition temperature of NiI 2 from 59 K (bulk) to 28 K (monolayer).

2D magnets↗

DOE OSTI2023

Source Characterization of Volatile Organic Compounds at Carlsbad Caverns National Park

Carlsbad Caverns National Park (CAVE), located in southeastern New Mexico, experiences elevated ground-level ozone (O 3 ) exceeding the National Ambient Air Quality Standard (NAAQS) of 70 ppbv. It is situated adjacent to the Permian Basin, one of the largest oil and gas (O&G) producing regions in the US. In 2019, the Carlsbad Caverns Air Quality Study (CarCavAQS) was conducted to examine impacts of different sources on ozone precursors, including nitrogen oxides (NO x ) and volatile organic compounds (VOCs). Here, we use positive matrix factorization (PMF) analysis of speciated VOCs to characterize VOC sources at CAVE during the study. Seven factors were identified. Three factors composed largely of alkanes and aromatics with different lifetimes were attributed to O&G development and production activities. VOCs in these factors were typical of those emitted by O&G operations. Associated residence time analyses (RTA) indicated their contributions increased in the park during periods of transport from the Permian Basin. These O&G factors were the largest contributor to VOC reactivity with hydroxyl radicals (62%). Two PMF factors were rich in photochemically generated secondary VOCs; one factor contained species with shorter atmospheric lifetimes and one with species with longer lifetimes. RTA of the secondary factors suggested impacts of O&G emissions from regions farther upwind, such as Eagle Ford Shale and Barnett Shale formations. The last two factors were attributed to alkenes likely emitted from vehicles or other combustion sources in the Permian Basin and regional background VOCs, respectively.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Light regulation of light-harvesting antenna size substantially enhances photosynthetic efficiency and biomass yield in green algae

One of the major factors limiting biomass productivity in algae is the low thermodynamic efficiency of photosynthesis. The greatest thermodynamic inefficiencies in photosynthesis occur during the conversion of light into chemical energy. At full sunlight the light-harvesting antenna captures photons at a rate nearly 10 times faster than the rate-limiting step in photosynthetic electron transport. Excess captured energy is dissipated by non-productive pathways including the production of reactive oxygen species. Substantial improvements in photosynthetic efficiency have been achieved by reducing the optical cross-section of the light-harvesting antenna by selectively reducing chlorophyll b levels and peripheral light-harvesting complex subunits. Smaller light-harvesting antenna, however, may not exhibit optimal photosynthetic performance in low or fluctuating light environments. Additionally, we describe a translational control system to dynamically adjust light-harvesting antenna sizes for enhanced photosynthetic performance. Furthermore, by expressing a chlorophyllide a oxygenase (CAO) gene having a 5' mRNA extension encoding a Nab1 translational repressor binding site in a CAO knockout line it was possible to continuously alter chlorophyll b levels and correspondingly light-harvesting antenna sizes by light-activated Nab1 repression of CAO expression as a function of growth light intensity. Significantly, algae having light-regulated antenna sizes had substantially higher photosynthetic rates and two-fold greater biomass productivity than the parental wild-type strains as well as near wild-type ability to carry out state transitions and non-photochemical quenching. These results have broad implications for enhanced algae and plant biomass productivity.

09 BIOMASS FUELS↗

DOE OSTI2020

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Assessment of Current MACCS Capabilities for Modeling Atmospheric Physical and Chemical Transformations

The physical and chemical transformation during atmospheric transport of radionuclides released into the environment has the possibility of impacting consequence modeling results. Accordingly, this report identifies physical and chemical transformations that may occur following release of chemically reactive radioactive species, how those transformations may affect modeling of consequences of release to the atmosphere and identifies current capabilities – in both MACCS and other state-of-practice atmospheric transport and dispersions models– to model those transformations. It was found that the inclusion of physical and chemical transformations is currently very limited in current state-of-practice codes for atmospheric dispersion of radionuclides. State-of-practice atmospheric dispersion codes appear to be typically limited to simulating either physical-chemical transformations or radioactive transformations, but not both. A state-of-practice atmospheric dispersion code capable of performing parallel physical, chemical, and radioactive transformation was not identified. A few atmospheric dispersion codes capable of modeling physical and chemical transport of specific species such as tritium or uranium hexafluoride were identified. Consequently, there is currently no information available that clearly suggests updates to the MACCS code are needed to bring it up to state-of-practice. However, investigations concluded that the MACCS computational framework can currently accommodate multiple physical/chemical forms in one simulation. Additionally, with some major assumptions, the computational framework in MACCS can accommodate parallel physical-chemical and radioactive transformations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

NASA NTRS1999

Measurements of Acidic Gases and Aerosol Species Aboard the NASA DC-8 Aircraft During the Pacific Exploratory Mission in the Tropics (PEM-Tropics A)

We received funding to provide measurements of nitric acid (HNO3), formic acid (HCOOH), acetic acid (CH3COOH), and the chemical composition of aerosols aboard the NASA Ames DC-8 research aircraft during the PEM-Tropics A mission. These measurements were successfully completed and the final data resides in the electronic archive (ftp-gte.larc.nasa.gov) at NASA Langley Research Center. For the PEM-Tropics A mission the University of New Hampshire group was first author of four different manuscripts. Three of these have now appeared in the Journal of Geophysical Research-Atmospheres, included in the two section sections on PEM-Tropics A. The fourth manuscript has just recently been submitted to this same journal as a stand alone paper. All four of these papers are included in this report. The first paper (Influence of biomass combustion emissions on the distribution of acidic trace gases over the Southern Pacific basin during austral springtime) describes the large-scale distributions of HNO3, HCOOH, and CH3COOH. Arguments were presented to show, particularly in the middle tropospheric region, that biomass burning emissions from South America and Africa were a major source of acidic gases over the South Pacific basin. The second paper (Aerosol chemical composition and distribution during the Pacific Exploratory Mission (PEM) Tropics) covers the aerosol aspects of our measurement package. Compared to acidic gases, O3, and selected hydrocarbons, the aerosol chemistry showed little influence from biomass burning emissions. The data collected in the marine boundary layer showed a possible marine source of NH3 to the troposphere in equatorial areas. This source had been speculated on previously, but our data was the first collected from an airborne platform to show its large-scale features. The third paper (Constraints on the age and dilution of Pacific Exploratory Mission-Tropics biomass burning plumes from the natural radionuclide tracer Pb-210) utilized the unexpectedly high concentrations of Pb-210 in the combustion plumes to estimate their ages and mixing along the transport route to the South Pacific basin. The final paper in the group (Tropospheric reactive-odd nitrogen over the South Pacific in austral springtime). This paper provides a summary of reactive nitrogen during PEM-Tropics A, with HNO3 and PAM showing the most impact from combustion emission.

Talbot, Robert W.↗

NASA NTRS2018

GMI-IPS: Python Processing Software for Aircraft Campaigns

NASA's Atmospheric Tomography Mission (ATom) seeks to understand the impact of anthropogenic air pollution on gases in the Earth's atmosphere. Four flight campaigns are being deployed on a seasonal basis to establish a continuous global-scale data set intended to improve the representation of chemically reactive gases in global atmospheric chemistry models. The Global Modeling Initiative (GMI), is creating chemical transport simulations on a global scale for each of the ATom flight campaigns. To meet the computational demands required to translate the GMI simulation data to grids associated with the flights from the ATom campaigns, the GMI ICARTT Processing Software (GMI-IPS) has been developed and is providing key functionality for data processing and analysis in this ongoing effort. The GMI-IPS is written in Python and provides computational kernels for data interpolation and visualization tasks on GMI simulation data. A key feature of the GMI-IPS, is its ability to read ICARTT files, a text-based file format for airborne instrument data, and extract the required flight information that defines regional and temporal grid parameters associated with an ATom flight. Perhaps most importantly, the GMI-IPS creates ICARTT files containing GMI simulated data, which are used in collaboration with ATom instrument teams and other modeling groups. The initial main task of the GMI-IPS is to interpolate GMI model data to the finer temporal resolution (1-10 seconds) of a given flight. The model data includes basic fields such as temperature and pressure, but the main focus of this effort is to provide species concentrations of chemical gases for ATom flights. The software, which uses parallel computation techniques for data intensive tasks, linearly interpolates each of the model fields to the time resolution of the flight. The temporally interpolated data is then saved to disk, and is used to create additional derived quantities. In order to translate the GMI model data to the spatial grid of the flight path as defined by the pressure, latitude, and longitude points at each flight time record, a weighted average is then calculated from the nearest neighbors in two dimensions (latitude, longitude). Using SciPya's Regular Grid Interpolator, interpolation functions are generated for the GMI model grid and the calculated weighted averages. The flight path points are then extracted from the ATom ICARTT instrument file, and are sent to the multi-dimensional interpolating functions to generate GMI field quantities along the spatial path of the flight. The interpolated field quantities are then written to a ICARTT data file, which is stored for further manipulation. The GMI-IPS is aware of a generic ATom ICARTT header format, containing basic information for all flight campaigns. The GMI-IPS includes logic to edit metadata for the derived field quantities, as well as modify the generic header data such as processing dates and associated instrument files. The ICARTT interpolated data is then appended to the modified header data, and the ICARTT processing is complete for the given flight and ready for collaboration. The output ICARTT data adheres to the ICARTT file format standards V1.1. The visualization component of the GMI-IPS uses Matplotlib extensively and has several functions ranging in complexity. First, it creates a model background curtain for the flight (time versus model eta levels) with the interpolated flight data superimposed on the curtain. Secondly, it creates a time-series plot of the interpolated flight data. Lastly, the visualization component creates averaged 2D model slices (longitude versus latitude) with overlaid flight track circles at key pressure levels. The GMI-IPS consists of a handful of classes and supporting functionality that have been generalized to be compatible with any ICARTT file that adheres to the base class definition. The base class represents a generic ICARTT entry, only defining a single time entry and 3D spatial positioning parameters. Other classes inherit from this base class; several classes for input ICARTT instrument files, which contain the necessary flight positioning information as a basis for data processing, as well as other classes for output ICARTT files, which contain the interpolated model data. Utility classes provide functionality for routine procedures such as: comparing field names among ICARTT files, reading ICARTT entries from a data file and storing them in data structures, and returning a reduced spatial grid based on a collection of ICARTT entries. Although the GMI-IPS is compatible with GMI model data, it can be adapted with reasonable effort for any simulation that creates Hierarchical Data Format (HDF) files. The same can be said of its adaptability to ICARTT files outside of the context of the ATom mission. The GMI-IPS contains just under 30,000 lines of code, eight classes, and a dozen drivers and utility programs. It is maintained with GIT source code management and has been used to deliver processed GMI model data for the ATom campaigns that have taken place to date.

Damon, M. R.↗

DOE OSTI2022

Identification of efficient and stable solid-solid electrochemical interfaces for Li-S batteries. Final Report

The objective of this work is to develop a better understanding of the interactions between battery components that drive the observed battery performance. Using computational tools at the atomistic level, we elucidate molecular, crystal, and amorphous structures, interfacial structures and reactions, and electrochemical behavior. We focus on Li metal anodes, liquid electrolytes, solid electrolytes, S cathodes, and metal-oxide cathodes, and their corresponding interfaces. Within liquid electrolytes, we study their bulk structure and dynamics, as a function of salt concentration. The analysis allows to understand the formation of localized high concentrated regions as the solvent is substituted by a weakly interacting diluent species. Analysis of the electron distribution yields the reactive behavior at interfaces. Moreover, we employ tools that allow the characterization of ion and electron transport at interfaces, during the processes of solvation and desolvation occurring while cycling. For Sulfur cathodes we elucidate synthesis processes and the role of carbon-sulfur and nitrogen-sulfur interactions not only to provide electronic conductivity but also to define alternative reaction pathways during discharge of Li-S batteries. Full-cell simulations offer new insights regarding the integration of all the battery components. Our studies also include solid electrolytes and their interfaces with Li metal, S-C, and metal-oxide electrodes.

25 ENERGY STORAGE↗