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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Understanding particle size effect on fast-charging behavior of graphite anode using ultra-thin-layer electrodes

Extreme fast charging (≤ 15 min) of lithium-ion batteries is highly desirable to accelerate mass-market adoption of electric vehicles. However, significant capacity fading, as well as safety issues due to the lithium plating caused by the fast charging rate, limit its implementation. In this study, we investigated the fast-charging capability of graphite materials with various particle sizes. To eliminate the Li + ion concentration gradient effect across the thickness of the electrode, ultra-thin-layer graphite electrodes were developed to investigate the "real" fast-charging capability of graphite at the particle level. Electrochemical assessments as well as microscopic characterizations revealed that smaller particles exhibited superior fast-charging performance, featuring enhanced capacity reversibility, faster charging rate, and less lithium plating under the same fast-charging conditions. It is shown that small-particle graphite (mean radius of 3.3 μm) could withstand a 4C charge (to 80 % state-of-charge) without plating, with minimal plating occurring at 6C. Thicker particles exhibited plating at lower C-rates. Since the experimental data could not directly explain whether intra-particle diffusion limitations or interfacial reaction limitations dominated the plating mechanism, the pseudo-2-dimensional model was used to evaluate the most likely plating mechanism. The model suggested that particle-level diffusion is the dominant mechanism contributing to plating at high rates. Finally, this work provides comprehensive insights into the particle size effects on fast-charging capability, offering a better understanding of fast-charging behavior and valuable guidance for designing optimal electrode architecture for high-rate lithium-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2023

On the Impact of Mechanics on Electrochemistry of Lithium-Ion Battery Anodes

Abstract Models exploring electrochemistry-mechanics coupling in liquid electrolyte lithium-ion battery anodes have traditionally incorporated stress impact on thermodynamics, bulk diffusive transport, and fracture, while stress-kinetics coupling is more explored in the context of all solid-state batteries. Here, we showcase the existence of strong link between active particle surface pressure and reaction kinetics affecting performance even in liquid electrolyte systems. Traction-free and immobile particle surface mechanical boundary conditions are used to delineate the varying pressure magnitudes in graphite host during cycling. Both tensile and compressive stresses are generated in traction-free case, while a fixed surface subjects the entire particle to a compression state. Pressure magnitudes are nearly two to three orders of magnitude higher for the latter resulting in significant depression of open circuit potential and improvement of exchange current densities compared to stress-free state. The results demonstrate the need for incorporating stress-kinetics linkage in models and provide a rationale for putting battery electrodes under compression to improve kinetics.

25 ENERGY STORAGE↗
NASA NTRSDate not supplied

In-situ Nanoscale Ablation

An understanding of the ablation of carbon-based materials is crucial to modeling the behavior of atmospheric entry spacecrafts equipped with thermal protection systems (TPS).Carbon is the backbone of TPS systems such as PICA (phenolic-impregnated carbon ablator).Therefore in the present work we study ablation of highly oriented pyrolytic graphite (HOPG)in oxygen at temperatures up to 1000°C done within a gas reaction cell housed in a scanning transmission electron microscope (STEM). Observation of the HOPG oxidation on specimens sectioned parallel and normal to the carbon basal planes in the presence of oxygen are reported.Pitting caused by residual oxygen and/or platinum particles from the sectioning process was observed before oxygen gas flow was established at 750°C in the specimen sectioned with the carbon basal planes parallel to the beam direction. Introduction of oxygen flow caused rapid oxidation moving in a uniform front that completely consumed the HOPG in approximately 1.5minutes. The specimen sectioned with the carbon basal planes normal to the beam direction did not show the same pitting phenomenon but exhibited rapid oxidation at 1000°C that proceeded in a uniform front and completed in approximately 1 minute. All specimens tested had a husk resembling the original specimen shape left over after oxidation. It is concluded that this husk is most likely ash impurity from the HOPG or impurities from the sectioning process or E-chip. A residual gas analyszer (RGA) was successfully used to monitor gas flows during the experiments but was yet unsuccessful in monitoring oxidation gas products produced during the experiment.Future studies will optimize the RGA setup in order to maximize the potential for detecting these species. These in-situ studies successfully show how this highly ordered carbon ablates on a nano- and micro-scopic length scale and can be used to provide fundamental understanding of carbon ablation that can be used to design the next-generation of TPS systems.

Nanoscale Ablation↗
DOE OSTI2020

Effect of Anode Porosity and Temperature on the Performance and Lithium Plating During Fast-Charging of Lithium-Ion Cells

Twenty-four single-layer approximate to 32 mAh pouch cells are tested to determine the effect of electrode porosity on lithium plating. Twelve cells contain a graphite electrode that is 26% porous, and 47% for the other twelve. The cells are cycled using a 6-C charge and a C/2 discharge protocol at temperatures in the range of 20-50 degrees C. A macro-homogeneous electrochemical model and microstructure analysis tool set are used to help interpret experimental observations for the effect of anode porosity and ambient temperature on fast-charging performance. Additionally, comparison between the two also highlights gaps in current theoretical understanding that need to be addressed. In post-test examination, lithium plating is seen in all cells, regardless of porosity. Elevated temperature is shown to reduce the amount of lithium plating and improve initial fast-charge capacity, but also changes the rate of other, less well-understood degradation mechanisms. Apparent kinetic rate laws, At + Bt 1/2 , where A and B are constants, can be fit to most of the capacity loss and resistance increase data. The relative magnitudes of A and B change with temperature and porosity. The capacity loss data at 50 degrees C from the high-porosity cells are fit by a logistics rate law.

25 ENERGY STORAGE↗
NASA NTRS1983

Shock processing of interstellar grains

Theoretical and observational evidence is presented that shock processing of interstellar dust grains by supernova blast waves affects both heavy element depletions and ultraviolet extinction curves. By coupling a realistic model of grain sizes and populations with a radiative shock code, significant grain destruction at velocities as low as 40 km/s is demonstrated. Nonthermal sputtering and grain-grain collisions destroy relatively more large grains than small, and more silicates than graphite. Consequently, both the 2175 A extinction 'bump' and the far-ultraviolet normalized extinction are increased in strength. Ultraviolet extinction studies with the International Ultraviolet Explorer of nine stars near three supernova remnants (the Monoceros Loop, Shajn 147, and Vela) exhibit strong 2175 A bumps and normal or high far-ultraviolet extinction. Diffuse bands, if they are created by small grains, should show little correlation with such activity.

Seab, C. G.↗
DOE OSTI2025

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↗
DOE OSTI2021

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↗
DOE OSTI2026

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.↗
NASA NTRS1983

The relative importance of aerosol scattering and absorption in remote sensing

The relative importance of aerosol optical thickness and absorption is illustrated through computing radiances for radiative transfer models. The radiance of sunlight reflected from models of the earth-atmosphere system is computed as a function of the aerosol optical thickness and its albedo of single scattering; it is noted that the albedo varies from 0.6 in urban environment to nearly 1 in areas with low graphitic carbon content. The calculations are applied to the example of satellite measurements of biomass. It is found that when surface classifications are made by means of clustering techniques the presence of gradients in the aerosol optical properties results in the dispersion of points in the plot correlating radiances viewed in two different directions. Finally, though such a remote sensing parameter as contrast is weakly affected by aerosol absorption, it is highly dependent on its optical thickness.

Fraser, R. S.↗
NASA NTRS1999

Analytical Prediction of Damage Growth in Notched Composite Panels Loaded in Axial Compression

A progressive failure analysis method based on shell elements is developed for the computation of damage initiation and growth in stiffened thick-skin stitched graphite-epoxy panels loaded in axial compression. The analysis method involves a step-by-step simulation of material degradation based on ply-level failure mechanisms. High computational efficiency is derived from the use of superposed layers of shell elements to model each ply orientation in the laminate. Multiple integration points through the thickness are used to obtain the correct bending effects through the thickness without the need for ply-by-ply evaluations of the state of the material. The analysis results are compared with experimental results for three stiffened panels with notches oriented at 0, 15 and 30 degrees to the panel width dimension. A parametric study is performed to investigate the damage growth retardation characteristics of the Kevlar stitch lines in the pan

Ambur, Damodar R.↗
DOE OSTI2025

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↗
DOE OSTI2022

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↗
DOE OSTI2026

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↗
DOE OSTI2025

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↗
NASA NTRS1979

An experimental and analytical investigation of the rail shear-test method as applied to composite materials

This report presents the results from an experimental and analytical investigation of the stress distributions occurring in a rail shear test. The effects of non-uniform stresses induced by differential thermal expansion, rail flexibility and specimen aspect ratio on measured shear modulus and ultimate strength of composite laminates are shown. A two-dimensional linearly elastic finite element model was used to analytically determine how various geometric parameters influenced the magnitude and distribution of inplane normal and shear stresses in a tensile rail shear specimen. Rail shear tests were conducted at room temperature and 589 K (600 F) on selected graphite-polyimide composite laminates using two titanium rail configurations. The analysis and test methods are discussed, and the results of the effects of the various parameters on shear modulus and ultimate strength are presented.

Garcia, R.↗
DOE OSTI2023

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↗
DOE OSTI2024

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↗
NASA NTRS1982

Predictions for partial and monolayer coverages of O2 on graphite

Monolayer properties of O2 on graphite are calculated using a pattern recognition, optimization scheme. Equilibrium monolayers are predicted at two different densities with properties in agreement with recent X-ray diffraction, specific heat, and neutron scattering data. Properties of the extremely low density regime are calculated using a model based upon a distribution of two-dimensional O2 clusters. The results are consistent with experimental evidence.

Pan, R. P.↗