Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Graphite modeling”

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 343 records · Page 19

On the Efficacy of Repeat Voltage Holds for Conditioning and Calendar Life Testing of Graphite and Silicon Cells

Voltage-hold (V-hold) protocols have shown promise toward calendar lifetime analysis of cells with graphite (Gr) and silicon (Si) anodes. In this work, repeat V-holds are performed on Gr and Si cells paired with lithium iron phosphate cathodes to delineate their beneficial role in formation and conditioning. We find that V-hold at the top of charge supplements constant current cycling in conditioning the cell to higher capacities for both Gr and Si cells after the first V-hold. A reduced order model provides the irreversible capacity proportions of each V-hold. With each repeat V-hold, parasitic loss of lithium to the solid electrolyte interphase (SEI) decreases on both Gr and Si cells. Gr cells show the square-root-of-time capacity loss behavior within 200 h of V-hold, indicative of its fast relaxation and low impact of reference performance test cycles on the SEI growth. Lifetime estimates from repeat V-holds on Gr can reach years. Si exhibits longer transition times from kinetic to diffusion-limited SEI growth, evidenced by the 400 h and 200 h holds showing square-root-of-time and linear behavior, respectively. Lifetime predictions from repeat V-holds on Si only reach 1–2 months, highlighting its limitations. Recommended duration of V-holds for Si cells should be ≥400 h.

25 ENERGY STORAGE↗

In situ synchrotron investigation of degenerate graphite nodule evolution in ductile cast iron

Ductile cast irons (DCIs) are of increasing importance in the renewable energy and transportation sectors. The distribution and morphology of the graphite nodules, in particular the formation of degenerate features during solidification, dictate the mechanical performance of DCIs. In situ high-speed synchrotron X-ray tomography was used to capture the evolution of graphite nodules during solidification of DCI, including degenerate features and the effect of the carbon concentration field. The degeneration of nodules is observed to increase with re-melting cycles, which is attributed to Mg-loss. The dendritic primary austenite and carbon concentration gradients in the surrounding liquid phase were found to control nodule morphology by locally restricting and promoting growth. A coupled diffusion-mechanical model was developed, confirming the experimentally informed hypothesis that protrusions form through liquation cracking of the austenite shell and subsequent localised growth. In conclusion, these results provide valuable insights into the solidification kinetics of cast irons, supporting the design of advanced alloys.

36 MATERIALS SCIENCE↗

The presence of Ru metal modifies the behavior of deuterium bound to pyridinic nitrogen in doped graphene-like materials

We used a combination of x-ray photoelectron spectroscopy (XPS) and density functional theory (DFT) to investigate how the presence of a ruthenium (Ru) metal substrate modifies the N-D bond strength of the pyridinic nitrogen species in nitrogen (N)-doped graphene (Gr) on Ru(0001) compared to metal-free highly-oriented pyrolytic graphite (HOPG) substrate. The N-dopants were introduced through low-energy N+/N2+ irradiation. N is embedded in the carbon materials in two predominant configurations: as graphitic N (GN: N substituted in the hexagonal C lattice) and pyridinic N (PN: substitutional N adjacent to a C vacancy). XPS data showed that upon atomic deuterium (D) exposure at 220 K, the PN peak shifted to a higher binding energy by +1.2 eV for HOPG and +1.0 eV for N-doped graphene on Ru(0001), while the GN peak remained unchanged, indicating that the D atoms bound solely to pyridinic N. 85% of PN sites on HOPG can be saturated with D atoms, whereas only ~30% of pyridinic N sites are able to bind D atoms in N-doped Gr/Ru(0001). Our DFT analysis shows that this difference is due to the coordination of PN to Ru atoms, which necessitates bond cleavage of the N-Ru interaction prior to D atom adsorption. D begins desorbing from N-HOPG at ~573 K and is fully desorbed by ~973 K, whereas desorption from N-doped graphene on Ru(0001) begins at ~290 K, with complete desorption observed at ~700 K, indicating that the Ru metal weakens the N-D bond strength. We also studied a high-surface-area, layered, and porous N-doped carbon material, internally labeled NC900, which was synthesized by pyrolysis of glucose and graphitic carbon nitride (g-C3N4) at 1173 K. D exposure caused a +1.1 eV PN peak shift in NC900, with no change in the GN peak. Moreover, the NC900 exhibited the same desorption behavior as HOPG, demonstrating that well-defined model systems can effectively capture the behavior of more complex N-doped carbon materials.

Alupothe Gedara, Buddhika S.↗

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗

Deconvoluting the benefits of porosity distribution in layered electrodes on the electrochemical performance of Li-ion batteries

Thick Lithium Ion Battery (LIB) electrodes suffer from poor rate capability and high ionic impedance due to their thickness and mesostructure. Therefore, optimizing thick electrode architectures becomes crucial. In this work, we report a systematic assessment of the ionic resistance in heterogeneous porous electrodes through the combination of computational simulations using a 4D-resolved model and experimental measurements. The first part of the study is devoted to a general assessment of Electrochemical Impedance Spectroscopy (EIS) spectra, mapping the impact of ionic and electronic resistances on the overall impedances of uncalendered and calendered LiNi 1/3 Mn 1/3 Co 1/3 O 2 , LiFePO 4 and graphite electrodes. In the second part, in silico-generated electrodes with different porosities are used in computational EIS simulations to analyze the impact of the electrode porosity on the ionic impedance. As expected, the results show that a lower porosity leads to a higher ionic impedance because of a higher electrode tortuosity factor. Furthermore, in silico-generated electrodes with different porosities were stacked and assembled to create heterogeneities of porosity along the thickness, and used in computational EIS and galvanostatic discharge simulations. The computational results show that the porosity heterogeneity along the electrode thickness has a significant effect on the ionic impedance and capacity of the electrode. The electrode architecture with progressively decreasing porosity from separator to current collector shows the highest performance, a trend validated by our in house experimental EIS and galvanostatic discharge also reported in this manuscript. Overall, we conclude that the ionic resistance in a thick electrode can be effectively reduced through proper tuning of the porosity heterogeneity. The proposed heterogeneous electrode architectures presented here could enormously help building efficient thick electrodes for LIBs.

25 ENERGY STORAGE↗

BMINN: Learning chemical potentials and parameters from voltage data for multi-phase battery modeling

Free-energy landscapes and chemical potentials govern the dynamics of phase transitions, transport, and stability in functional materials, yet they remain experimentally inaccessible under realistic operating conditions. Here we introduce a Bayesian model-integrated neural network (BMINN) that embeds physics-based formulations of non-autonomous partial differential-algebraic equations into probabilistic learning. This approach reconstructs hidden thermodynamics directly from macroscopic current-voltage data, providing quantitative access to metastable states, staging transitions, and energy barriers without synchrotron probes. Demonstrated on lithium-graphite electrodes, BMINN recovers full Gibbs free-energy landscapes with fidelity validated against operando X-ray diffraction. The framework generalizes across dynamical regimes, enabling accurate voltage prediction, internal state estimation, and inference of governing parameters. Beyond batteries, BMINN exemplifies a broadly applicable strategy for learning missing physics in multiphase, non-equilibrium systems, offering a new pathway to uncover hidden thermodynamic functions across condensed matter and materials physics.

25 ENERGY STORAGE↗

Data-driven design of electrolyte additives supporting high-performance 5 V LiNi 0.5 Mn 1.5 O 4 positive electrodes

LiNi 0.5 Mn 1.5 O 4 (LNMO) is a high-capacity spinel-structured material with an average lithiation/de-lithiation potential at ca. 4.6–4.7 V vs Li + /Li, far exceeding the stability limits of electrolytes. An efficient way to enable LNMO in lithium-ion batteries is to reformulate an electrolyte composition that stabilizes both graphitic (Gr) negative electrode with solid-electrolyte-interphase and LNMO with cathode-electrolyte-interphase. In this study, we select and test a diverse collection of 28 single and dual additives for the Gr||LNMO battery system. Subsequently, we train machine learning models on this dataset and employ the trained models to suggest 6 binary compositions out of 125, based on predicted final area-specific-impedance, impedance rise, and final specific-capacity. Such machine learning-generated new additives outperform the initial dataset. This finding not only underscores the efficacy of machine learning in identifying materials in a highly complicated application space but also showcases an accelerated material discovery workflow that directly integrates data-driven methods with battery testing experiments.

batteries↗

Insights into the thermal history of carbon fibers using Raman spectroscopy and a novel kinetic model

Carbon fibers and carbon fiber composites are applied in high-performance applications, but a key consideration for application is their relative sensitivity to oxidative environments. To enable in-situ characterization of carbon fibers exposed to oxidative conditions, the Raman spectral response of T700 carbon fibers that have been exposed to a variety of dwell temperatures is reported herein with dwell times reaching up to 1 month. Here, we evaluate the spectra holistically by using integrated absolute difference analysis. By combining this analysis with straightforward kinetic models, we connect the total Raman spectral response to the temperature-time curve that could yield such a shift in spectral parameters. Our work connects the Raman spectral response of carbon fibers to their thermal history and can easily be extended to other graphitic materials, such as nuclear graphite.

36 MATERIALS SCIENCE↗

SCALE inventory and reactivity analysis as part of the Hermes 2021 PSAR review

The readiness of SCALE for comprehensive studies of pebble-bed reactors has been demonstrated through detailed analysis of a fluoride salt–cooled, high-temperature pebble-bed reactor (PB-FHR). The methods developed for pebble-bed reactor modeling in SCALE, particularly for inventory generation, have proven effective in gaining insights into the reactor physics of this advanced reactor. Excellent agreement with another code package has been observed, further highlighting SCALE’s strong performance. The SCALE results supported the US Nuclear Regulatory Commission’s construction permit application review of the Hermes low-power PB-FHR demonstration reactor. A SCALE model of the Hermes reactor was developed at Oak Ridge National Laboratory using information from the Preliminary Safety Analysis Report (PSAR) and supplemented with publicly available data. SCALE reactivity coefficient simulations reproduced PSAR results within 1σ statistical uncertainties. Sensitivity studies emphasized the importance of graphite specifications for accurate keff predictions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Tritium Transport Phenomena in Molten-Salt Reactors

In this work, we review phenomena relevant to tritium transport in molten-salt reactors, which produce tritium from lithium and beryllium salts at significantly higher levels than other reactor types. Modeling of such phenomena began following MSRE operations, and these early models attempted to predict measured tritium distributions in the MSRE, accounting for turbulent mass-transport processes (using established heat transfer correlations), permeation through a variety of metal structures such as heat exchanger tubes, and transport to and from bubbles introduced into the salt by gas sparging. The models reasonably reproduced the MSRE data, but did so best when the permeability of structures was reduced by about a factor of 1,000. This issue does not appear to have been conclusively resolved, and all of the more recent attempts to model tritium transport in molten salts appear to make use of the same methodology. MSRE remains, however, essentially our only source of integral tritium transport data relevant to MSRs. Here, we generalize the MSRE approach to permeation in order to include potential rate-limiting effects at interfaces, as well the effects of added hydrogen. Appropriately non-dimensionalized, this system of equations identifies two dimensionless numbers whose relative values clearly delineate the conditions under which mass transport, surface effects, permeation, and hydrogen swamping are expected to become rate-limiting. We also describe the preliminary conceptual design of a forced convection FLiBe loop, into which tritium would be introduced for the purpose of providing validation data for such a model. The primary purpose of this is to investigate the coupled transport phenomena described above and identify those that are rate-limiting in MSRs. Additional test section configurations are described that would address other transport phenomena relevant to MSRs, including bubbly flows and graphite interactions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

3D TRISO particle-explicit compact meshing

The TRI-structural ISOtropic (TRISO) layered fuel particle is a robust nuclear fuel form offering enhanced safety and performance for advanced reactor concepts, including high-temperature gas-cooled reactors and other Generation IV designs. These poppy-seed-sized particles are embedded in a graphite matrix to form fuel elements that must withstand elevated temperatures and high burn-up levels. The heterogeneous nature of these fuel elements — comprising thousands of randomly distributed TRISO particles — produces complex stress fields and thermal gradients that one- and two-dimensional models cannot accurately capture. While three-dimensional modeling has improved predictions of dimensional changes, internal pressure buildup, and fission product transport under irradiation, current approaches rely on homogenized material properties that are known to have considerable divergence from experimental observations. This work presents a methodology for optimized random packing of TRISO fuel compacts and full three-dimensional mesh generation within the BISON fuel performance code, with each particle coating layer individually discretized. The resulting mesh was demonstrated through heat conduction simulations under representative in-reactor operating conditions, showing strong agreement with expected behavior. This capability enables detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating layer failures — all of which directly govern fuel performance and safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Radial Deconsolidation and Leach-Burn-Leach of AGR-3/4 Compact 1-4 and 10-4

The Advanced Gas Reactor (AGR) Fuel Development and Qualification Program third and fourth irradiation experiments (AGR-3/4), originally planned as separate tests, were combined in one test train for irradiation in the Advanced Test Reactor at Idaho National Laboratory (INL). The irradiation test began on December 14, 2011 and ended on April 12, 2014 (Collin 2016). The originally planned AGR-3 and AGR-4 irradiation experiments were both focused on obtaining data on fission product transport to support the improvement of modeling. The AGR-3 experimental plan was focused on gaseous and metallic fission product release from the kernels and diffusion in the coatings during irradiation and postirradiation safety testing. The AGR-4 experimental plan was focused on diffusivities and sorptivities in the compact matrix and reactor graphite (Petti et al. 2005). These two goals were combined in the AGR-3/4 irradiation, which consisted of twelve independently monitored capsules that each contained four AGR-3/4 compacts in a single stack surrounded by an inner ring of matrix or graphite and an outer ring of graphite. There were two capsule types: a standard capsule and a “fuel body”, in which the outer graphite ring included a floor and cap that fully enclosed the fuel (Stempien et al. 2018a). The fuel body design supported post-irradiation safety testing of the intact fuel and ring assembly to provide data on fission product transport and release from matrix and graphite at accident temperatures (Demkowicz 2017).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Radial Deconsolidation and Leach-Burn-Leach of AGR-3/4 Compacts 8-4 and 7-4

The Advanced Gas Reactor (AGR) Fuel Development and Qualification Program’s third and fourth irradiation experiments (AGR-3/4), originally planned as separate tests, were combined into one test train for irradiation in the Advanced Test Reactor at Idaho National Laboratory (INL). The irradiation test began on December 14, 2011, and ended on April 12, 2014 (Collin 2016). The originally planned AGR-3 and AGR-4 irradiation experiments were both to be focused on obtaining fission product transport data to support improvements in modeling. The AGR-3 experimental plan was focused on gaseous and metallic fission product release from the kernels and diffusion in the coatings during irradiation and postirradiation safety testing. The AGR-4 experimental plan was focused on diffusivities and sorptivities in the compact matrix and reactor graphite (Petti et al. 2005). These goals were combined in the AGR-3/4 irradiation experiment, which consisted of twelve independently monitored capsules that each contained four AGR-3/4 compacts in a single stack surrounded by an inner ring of matrix or graphite and an outer ring of graphite. There were two capsule types: a standard capsule and a fuel body, in which the outer graphite ring included a floor and cap to fully enclose the fuel (Stempien et al. 2018a). The fuel body design supported post-irradiation safety testing of the intact fuel and ring assembly to obtain data on fission product transport and release from matrix and graphite at accident temperatures (Demkowicz 2017).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Hydrogen Burning of 29 Si and Its Impact on Presolar Stardust Grains from Classical Novae

Abstract Presolar stardust grains found in primitive meteorites are believed to retain the isotopic composition of stellar outflows at the time of grain condensation. Therefore, laboratory measurements of their isotopic ratios represent sensitive probes for investigating open questions related to stellar evolution, stellar explosions, nucleosynthesis, mixing mechanisms, dust formation, and galactic chemical evolution. For a few selected presolar grains, classical novae have been discussed as a potential source. For SiC, silicate, and graphite presolar grains, the association is based on the observation of small N ( 12 C)/ N ( 13 C) and N ( 14 N)/ N ( 15 N) number abundance ratios compared to solar values, and abundance excesses in 30 Si relative to 29 Si, as previously predicted by models of classical novae. We report on a direct measurement of the 29 Si(p, γ ) 30 P reaction, which strongly impacts simulated δ 29 Si values from classical novae. Our new experimental 29 Si(p, γ ) 30 P thermonuclear reaction rate differs from previous results by up to 50% in the classical nova temperature range ( T = 100–400 MK), while the rate uncertainty is reduced by up to a factor of 3. Using our new reaction rate in Monte Carlo reaction network and hydrodynamic simulations of classical novae, we estimate δ 29 Si values with much reduced uncertainties. Our results establish δ 29 Si values measured in presolar grains as a sensitive probe for assessing their classical nova paternity. We also demonstrate that δ 30 Si values from nova simulations are currently not a useful diagnostic tool unless the large uncertainty of the 30 P(p, γ ) 31 S reaction rate can be significantly reduced.

79 ASTRONOMY AND ASTROPHYSICS↗

Computational Analysis of Anode and Cathode Structuring Effects on Charge and Discharge in Graphite|LiNi 0.6 Mn 0.2 Co 0.2 O 2 Batteries

Structured electrodes (SEs) improve the rate capability of Lithium-ion batteries by engineering micrometer-scale electrolyte regions into the electrode, promoting rapid ionic transport. Prior research has focused on structuring one electrode (anode or cathode) with an analysis on either the charge or discharge performance. We present a holistic study using three-dimensional models to investigate the isolated effects of structuring either electrode and the combined effects of structuring both electrodes on the charge and discharge capacity of single-layer cells at 4 C and 6 C. Volumetric and gravimetric discharge energy density (Wh/L stack and Wh/kg stack ) and charge capacity (Ah/kg stack and Ah/L stack ) are evaluated for multi-layer pouch cell stacks. Pairing SE anodes with SE cathodes demonstrated improvements up to 15% in discharge Wh/kg stack and up to 33% in charge Ah/kg stack over a conventional cell; Energy required to charge per Ah/kg stack was improved by 13%–14%. SE cathodes paired with a conventional anode exhibited improvements of 0.3%–22% across all performance metrics evaluated. Conversely, pairing a SE anode with a conventional cathode demonstrated improved charge capacity up to 13% but showed a 2%–23% lower discharge energy density. The importance of aligning SEs in a cell from a performance and manufacturing perspective is also analyzed.

25 ENERGY STORAGE↗

Development of Nanocrystalline Graphite from Lignin Sources

Carbon composites are attractive to a variety of high-impact applications, such as carbon fibers, batteries, and vehicle parts, due to their multifunctional properties. The properties of carbon are highly dependent on the allotrope the carbon takes and the functionality, impurities, and defects contained within the structure. The increase in demand for sustainable carbon sources in energy storage devices motivates interest in understanding synthesis parameters of lignin value-added products. Also, as the dependence on oil for fuel decreases, alternative sources for carbon in many applications will be needed. In this work, the thermochemical conversion of lignin powders from different feedstocks was evaluated via small and wide-angle X-ray scattering techniques to resolve the amorphous, disordered, and crystalline domains present in the lignin carbons. Scattering analyses indicated an evolution of hierarchical structures along with an increase in ordered domains as a function of carbonization temperature. Qualitative and quantitative methods were used to describe isotropic scattering intensity profiles at multiple length scales. The use of power law models in the mesoscopic region served as the basis to describe morphological changes related to structural features, for example, graphene stacking, degree of roughness, and surface fractals. Kraft softwood and switchgrass produced carbon powder with the most crystalline domains and the least surface roughness. Softwoods reached the highest degree of crystallinity followed by switchgrass samples and had less variability in particle sizes. These results suggest lignin carbons extracted from softwoods and switchgrass are viable substitutes for graphite. Interpretation of X-ray scattering data from lignin carbon powders elucidates feedstock- and processing-dependent morphological features across multiple length scales providing a straightforward framework to evaluate the feasibility of leveraging lignin carbons for producing tunable application-specific materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Leveraging Curvature on N–Doped Carbon Materials for Hydrogen Storage

Carbon sorbent materials have shown great promise for solid-state hydrogen (H 2 ) storage. Modification of these materials with nitrogen (N) dopants has been undertaken to develop materials that can store H 2 at ambient temperatures. In this work density functional theory (DFT) calculations are used to systematically probe the influence of curvature on the stability and activity of undoped and N-doped carbon materials toward H binding. Specifically, four models of carbon materials are used: graphene, [5,5] carbon nanotube, [5,5] D 5d -C 120 , and C 60 , to extract and correlate the thermodynamic properties of active sites with varying degrees of sp 2 hybridization (curvature). From the calculations and analysis, it is found that graphitic N-doping is thermodynamically favored on more pyramidal sites with increased curvature. In contrast, it is found that the hydrogen binding energy is weakly affected by curvature and is dominated by electronic effects induced by N-doping. These findings highlight the importance of modulating the heteroatom doping configuration and the lattice topology when developing materials for H 2 storage.

08 HYDROGEN↗

Towards a Benchmark Experiment with the Compact Nuclear Power Source

The Compact Nuclear Power Source (CNPS) was a high-assay low enriched uranium (HALEU) tristructural isotropic (TRISO)-fueled, graphite-moderated microreactor constructed in 1987 at Los Alamos National Laboratory. The reactor, conceived as a power source for short-range radar stations, was designed to be "walk-away safe," and was cooled by heat pipes and ambient air. Though the project was formally cancelled after the fuel and moderator material had been received, its potential to meaningfully advance the body of critical and integral data was evident, and critical experiments with a mock-up of the reactor proceeded at TA-18 until its disassembly in 1991. The reflector and some components of the core would later go on to see service as part of the New Production Reactor Modular High-Temperature Gas-Cooled Reactor (NP-MHTGR) critical experiments. In light of the progress made by groups like Westinghouse and X-Energy towards contemporary graphite-moderated microreactors, the system remains an attractive candidate for the basis of a benchmark experiment even today. Uncertainties remain, however, in some features of the system–inconsistencies in the dimensions and composition of the reactor’s components as described in literature and implemented in current computational models. The identification, quantification, and to the extent possible, the minimization of these uncertainties is a crucial task on the path towards creating a benchmark based on the CNPS system. The present work seeks to initiate this process with a survey of available documentation and improvements to the accuracy of CNPS neutronics models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗