SEARCH · Engineering Papers
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.
Introduction of ART Graphite R&D to include oxidation activities, oxidation resistant graphite, model development, ASME component failure, ASME code development (design rules), ceramic composites, AGC update, molten salt intrusion, split-disk studies, wear testing, and concluding remarks.
Graphite is a critical material used as the negative electrode in lithium-ion batteries. Both natural and synthetic graphites are utilized, with the latter obtained from a range of carbon raw materials. In this paper, efforts to synthesize graphite from coal as a domestic feedstock for synthetic graphite are reported. Domestic coal-derived graphite could address national security and energy issues by standing up domestic supply chains for battery critical materials. The performance in lithium-ion coin cells of this coal derived graphite is compared to a commercial battery-grade graphite. For the first time, a multi-species, multi-reaction (MSMR) modeling technique is applied to synthetic graphite derived from coal. Key thermodynamic, transport, and kinetic parameters are obtained for the coal derived graphite and compared to the same parameters for commercial battery-grade graphite. Modeling of synthetic graphites will allow for virtual evaluation of these materials toward production of domestically sourced graphite.
High-fidelity modeling and simulation were performed for the MIT graphite exponential pile (MGEP) using Monte Carlo neutron transport codes OpenMC and MCNP, and the results were validated by experimental data. The MGEP is being used as the test bed for the design of an autonomous control system for the pile's neutron flux distribution. The main contribution of this work is to generate the training data sets of neutron flux distributions with different locations of control rods that perturb the neutron flux profiles. First, code -to-code cross verification between OpenMC and MCNP was performed to ensure consistency of the numerical modeling within statistical uncertainties. To validate the accuracy of this high-fidelity model, a series of neutron flux measurements were conducted using a Helium-3 (He-3) neutron detector on a mobile platform that is placed inside the pile. Second, the neutron flux profiles were measured in four vertical layers of interest, and compared to the corresponding simulation results. The comparison results shows that the root mean square error is less than 2.5% in the two upper layers, and less than 4.5% in all four measured layers. Here the results validated the accuracy of the modeling and simulation. Finally, the relative change of the neutron flux profiles from moving control rods was analyzed, which identified the layer that has the best sensitivity regarding the control rods movements. Thus, this work identified and provided training data sets of both simulated and experimental neutron flux profiles in the most sensitive layer, paving the path forward to the real-time experimental demonstration of the autonomous control system.
Nuclear-grade graphites are used extensively in the core designs of multiple types of advanced nuclear reactors. In the reactor environment, graphite is exposed over long durations to extreme conditions, including high temperatures, radiation and potentially molten salt and oxygen. Exposure to these conditions can cause several degradation mechanisms in graphite, including nonuniform volumetric strains induced by irradiation and thermal expansion, which lead to stresses that can compromise the performance of graphite components. Evaluating component integrity, predicting component performance over the reactor lifetime, and developing design standards all require robust tools for predicting fracture initiation and propagation in graphite structural components in nuclear reactors. This report documents progress in an ongoing effort to develop modeling and simulation tools in the Grizzly code for predicting the performance of graphite exposed to reactor conditions. Recent developments include a set of thermal and mechanical models that now include the IG-110, NBG-18, and H-451 graphite grades. Improvements have also been made to a nonlinear damaged plasticity model applicable to predicting damage under tension and compression to quasibrittle materials, including graphite. In addition, enhancements have been made to the extended finite element method implementation targeted at simulating graphite fracture. These include new capabilities for crack nucleation in the interior of a solid body, improved treatment of crack nucleation on free surfaces, and more robust modeling approaches for crack growth approaching free surfaces or other cracks.
A better estimation of surface reaction efficiency of semiconductor-grade graphite with atomic nitrogen, as well as the calibration error are calculated using Bayesian updating based on experimental data. Compared with a conventional deterministic model, the stochastic model approach is a powerful tool in the sense that the model is capable of taking into account underlying error correlations among the data quantities. In this paper, we investigate four different stochastic models (called “stochastic system model classes” herein) corresponding to different descriptions of modeling and measurement error structures, given one deterministic physical model. These stochastic system model classes differ in the covariance matrix structure that is used in the uncertainty model to represent uncertainties associated with the physical model and experimental measurements. For each model class, Bayesian inference is used to estimate the posterior probabilities of the physical model parameters as well as of the stochastic model parameters. Model comparison and selection are then applied based on two measures including Bayesian evidence and Bayesian information criterion, as well as the deviance information criterion. Both measures suggest the stochastic model class, which considers that a correlation between errors in two data quantities among different data points is the most plausible. With the stochastic model class, the range of uncertainty in surface reaction efficiency is estimated to be about two orders of magnitude at [Formula: see text].
Lithium plating is commonly observed in anodes charged at fast rates, and can lead to capacity loss and battery safety issues. The increased risk of plating has been attributed to transport limitations, and architectured electrodes may reduce plating risk. However, while theoretical studies have shown that reaction non-uniformity arises due to interplay of transport limitations, anode open circuit voltage behavior and reaction kinetics, its effect on lithium plating has not been studied. We use analytic and numerical simulations to predict onset of plating in graphite anode half-cells at high C-rates and demonstrate how anodes with layered porosities can delay plating. Simplified analytical models identify trends for plating onset and predictions are calibrated against numerical models. A calibrated numerical model of graphite demonstrates qualitative agreement with analytical model predictions. This reaction inhomogeneity mechanism occurs in the absence of lithium ion depletion, indicating that these mechanisms may contribute to capacity loss independently or simultaneously. A bilayer model of graphite exhibits delayed plating onset, and an optimization procedure is presented. This theoretical work presents quantitative and mechanistic insight on how reaction inhomogenity affects lithium metal plating onset and can be used as a guide to engineer anodes resistant to lithium plating.
An experimental method is introduced that uses microprobe Li/Cu reference electrodes to map "safe" charging regimes to avoid Li metal plating on graphite anodes at rates to 8C. With this method, a Li-ion cell closely approaches the Li plating condition but avoids descending into this regime. This "safe" approach is used to demonstrate effects of the microporous separator on Li plating in cells consisting of a layered nickel-manganese -cobalt oxide cathode and graphite anode. We argue that most of the difference in behaviors for different separators arises from time delayed Li + ion percolation through the membrane. In conclusion, our study suggests that explicitly treating kinetics of phase transitions in the lithiated graphite anode is essential for electrochemical models predicting "safe" lines for fast charging cells.
The importance of graphite-hydrogen chemical reactions to fusion, fission, and hydrogen storage applications, combined with the rapidly evolving knowledge on the underlying mechanisms, has led to the development of multiple models to describe hydrogen transport in graphite. Significant differences exist among these models, resulting from discrepancies in the modeling assumptions, intended degree of fidelity, and conditions of applicability. This paper attempts at reconciling these apparent differences by providing a comprehensive description of the constitutive equations governing hydrogen transport in graphite at high-temperature, identifying outstanding gaps in knowledge, illustrating how these different models approach them, and proposing alternative analytical formulations grounded on experimental results from hydrogen-graphite studies. Governing equations, closing relations, and simplifying assumptions are discussed for hydrogen transport at the inter-granular and intra-granular level, accompanied by compiled experimental data and illustrated energy diagrams associated to the proposed transport mechanisms. Analytical formulations are provided to reproduce competing hypotheses on the mechanisms, supporting the development of a range of computational models that can enable resolution of outstanding knowledge gaps through comparative testing against experimental data.
The SAM code is under development and supported by DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. These advanced reactor concepts incorporate novel and improved approaches to achieve safety and economic feasibility. This report summarizes recent efforts and progress in addressing the code capability gaps in SAM for the modeling of the air-ingress phenomenon in High-Temperature Gas-Cooled Reactors (HTGRs). The capability enhancements implemented to the code include: a multi-component flow model to capture the air-helium mixture during air ingress, a 0-D graphite oxidation model to capture the reaction of graphite with oxygen in the air, and an isentropic critical flow model to more accurately predict the de-pressurization of the reactor system due to a small break in the primary loop. Verification cases and demonstrations are provided to showcase these capabilities in SAM. Existing gaps in the code’s capability are also identified.
Explore the source record for details and available documents.
General description of the ART Graphite Program with Status for 2023. This includes list of staff members, and brief description of behavior models, licensing & code, graphite R&D, As-Fab'd Properties, Machanism and Analysis, Irradiation. Also detail FY23 Activites, with AGC Update, Oxidation Activiteis, Baseline STatus, ASME Code Development, VIC Project Update, Molten salt studies. Slide included to highlight each employee contributions for work and papers and NEUPs.
The thermo-mechanical response of an ATJ graphite sample to controlled runaway electron (RE) dissipation, realized in DIII-D, is modelled with a novel work-flow that features the RE orbit code KORC, the Monte Carlo particle transport code Geant4 and the finite element multiphysics software COMSOL. KORC provides the RE striking positions and momenta, Geant4 calculates the volumetric energy deposition and COMSOL simulates the thermoelastic response. Brittle failure is predicted according to the maximum normal stress criterion, which is suitable for ATJ graphite owing to its linear elastic behavior up to fracture and its isotropic mechanical properties. Measurements of the conducted energy, damage topology, explosion timing and blown-off material volume, impose a number of empirical constraints that suffice to distinguish between different RE impact scenarios and to identify RE parameters which provide the best match to the observations.
This study was motivated by the need to incorporate more physics-based mass transfer models in BISON in order to improve fuel particle and compact matrix fission product release prediction. This document details the integration of new modeling capabilities in BISON, including (1) development of sorption mass transfer models; improvement of thermal contact model by adding neon to the gas gap inventory, (3) implementation of PCEA and IG-110 graphite thermal models, (4) validation using Advanced Gas Reactor (AGR)-3/4 compact and capsule rings data. These new capabilities have been shown to enforce the desired physics with satisfactory accuracy. BISON’s predictions of fission product release of AGR-3/4 compacts compare favorably with PARFUME and the concentration profiles across capsule rings show reasonable trend to the experiment measurement.
Flat bands and nontrivial topological physics are two important topics of condensed matter physics. With a unique stacking configuration analogous to the Su–Schrieffer–Heeger model, rhombohedral graphite (RG) is a potential candidate for realizing both flat bands and nontrivial topological physics. Here, in this study, we report experimental evidence of topological flat bands (TFBs) on the surface of bulk RG, which are topologically protected by bulk helical Dirac nodal lines via the bulk-boundary correspondence. Moreover, upon in situ electron doping, the surface TFBs show a splitting with exotic doping evolution, with an order-of-magnitude increase in the bandwidth of the lower split band, and pinning of the upper band near the Fermi level. These experimental observations together with Hartree–Fock calculations suggest that correlation effects are important in this system. Our results demonstrate RG as a platform for investigating the rich interplay between nontrivial band topology, correlation effects, and interaction-driven symmetry-broken states.
Immobilized cobalt phthalocyanine (CoPc) is a highly promising architecture for the six-proton, six-electron reduction of CO 2 to methanol. This electroreduction process relies on proton-coupled electron transfer (PCET) reactions that can occur by sequential or concerted mechanisms. Immobilization on a conductive support such as carbon nanotubes or graphitic flakes can fundamentally alter the PCET mechanisms. We use density functional theory (DFT) calculations of CoPc adsorbed on an explicit graphitic surface model to investigate intermediates in the electroreduction of CO 2 to methanol. Our calculations show that the alignment of the CoPc and graphitic electronic states influences the reductive chemistry. These calculations also distinguish between charging the graphitic surface and reducing the CoPc and adsorbed intermediates as electrons are added to the system. This analysis allows us to identify the chemical transformations that are likely to be concerted PCET, defined for these systems as the mechanism in which protonation of a CO 2 reduction intermediate is accompanied by electron abstraction from the graphitic surface to the adsorbate without thermodynamically stable intermediates. Furthermore, this work establishes a mechanistic pathway for methanol production that is consistent with experimental observations and provides fundamental insight into how immobilization of the CoPc impacts its CO 2 reduction chemistry.
The phase separation dynamics in graphitic anodes significantly affects lithium plating propensity, which is the major degradation mechanism that impairs the safety and fast charge capabilities of automotive lithium-ion batteries. In this study, we present comprehensive investigation employing operando high-resolution optical microscopy combined with non-equilibrium thermodynamics implemented in a multi-dimensional (1D+1D to 3D) phase-field modeling framework to reveal the rate-dependent spatial dynamics of phase separation and plating in graphite electrodes. Here we visualize and provide mechanistic understanding of the multistage phase separation, plating, inter/intra-particle lithium exchange and plated lithium back-intercalation phenomena. A strong dependence of intra-particle lithiation heterogeneity on the particle size, shape, orientation, surface condition and C-rate at the particle level is observed, which leads to early onset of plating spatially resolved by a 3D image-based phase-field model. Moreover, we highlight the distinct relaxation processes at different state-of-charges (SOCs), wherein thermodynamically unstable graphite particles undergo a drastic intra-particle lithium redistribution and inter-particle lithium exchange at intermediate SOCs, whereas the electrode equilibrates much slower at low and high SOCs. These physics-based insights into the distinct SOC-dependent relaxation efficiency provide new perspective towards developing advanced fast charge protocols to suppress plating and shorten the constant voltage regime.
The ASTM International Standard Test Method for Tensile Stress-Strain of Carbon and Graphite (ASTM C749) requires dog-bone shaped specimens. Unfortunately, specimen sizes required in ASTM C749 are much too large for irradiation within a test reactor and are incompatible with irradiation capsule volumes, molten salt degradation facilities, or oxidation apparatus. The ASTM International Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite (ASTM D8289) was developed as an alternative that provides a means for testing tensile properties using smaller specimens when geometric constraints are present. However, the important stress–strain relation throughout the test loading history that could be recorded using ASTM C749 is lost by using ASTM D8289 because of the difficulties involved with using traditional strain gauges or extensometers on such a small specimen. In this work, a digital image correlation (DIC) method, which uses a full-field noncontact surface displacement measurement technique, was applied along with the ASTM D8289 splitting tensile test on eight small samples of graphite grade Mersen 2114 to track their deformation during compressive loading. Initial results showed that the DIC technique can measure the surface displacement/strain on these small (Ø6 × 3 mm) graphite specimens with high accuracy and good repeatability, indicating that the DIC technique, combined with the small-disc splitting tensile test in ASTM D8289, can capture the strain evolution in the test. In conclusion, further comparison of measured strains with analytical and numerical simulation results provides additional understanding of the mechanical behaviour of this nuclear graphite material in complex stress states.
The goal of this document is to provide a set of comprehensive benchmarking neutronics problems for gas-cooled, graphite-moderated reactors. The problems are designed to be representative problems that increase in complexity to test performance of various codes. For the simulation of the benchmarking problems, the Shift Monte Carlo code and MPACT deterministic code are chosen because of their state-of-the-art capabilities, which are optimized for various reactor designs. Because these codes have been developed primarily for light water reactor applications, their application to graphite-moderated, gas-cooled reactors has not been extensively explored. This report attempts to evaluate benchmarks for a Magnox-style reactor, one that is graphite-moderated and gas-cooled. Parameters that are calculated and presented in this benchmark include the eigenvalue, peaking factors, and isotope concentrations. Results show good agreement between Shift and MPACT, and demonstrate the capability of MPACT to model graphite-moderated systems well.