Search NASA⌕ Search

SEARCH · Search NASA

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

"Source Term Modeling for Advanced Gas Micro-Reactors"

Maintaining the safety of the public, environment, and operating personnel is the most important factor in designing, operating, maintaining, and decommissioning nuclear reactors. In recent years, there has been a growing interest in the development of micro-reactors employing TRi-structural ISOtropic (TRISO)-coated particle fuel. In gas reactors, TRISO fuel plays an important role in the safety case for high temperature reactors because of the fission product retention properties of the fuel. This ability enables the use of a functional containment strategy for the reactor where multiple barriers are used to prevent fission product release to the environment. Part of the safety analysis of these advanced reactors is the assessment of radionuclide releases under normal and accident conditions through the multiple credited safety barriers. Using conservative assumptions, a mechanistic analysis can be performed to quantify these releases that combines the probabilistic assessment of failure with analytic solutions to radionuclide transport equations. Source term modeling for TRISO fuel has been performed for previous reactor designs; however, these models are outdated, in many cases proprietary, and need updates to be applied to the current state of TRISO fuel technology and alternative gas reactor core configurations [1]. Currently, the only publicly available source term assessment for gas reactors is an expert-based Monte Carlo simulation based on the effectiveness of the fuel kernel, coating layers, and graphite block in a modular high temperature gas reactor [2]. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that could be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. The release is calculated by the diffusion of the key safety important fission products through the kernel, silicon carbide (SiC), graphite for both intact and defective TRISO particles based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system to remove fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. The model then can calculate the fission product release for any transient temperature profile and fission product releases can then be used to assess radiological dose to the workers and the public using conventional dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy (DOE) Advanced Gas Reactor (AGR) TRISO fuel development program. The model is coded in python with inputs and outputs in excel spreadsheets, as well as python plotting utilities to aid in the interpretation of the results. References: [1] INL, NGNP Mechanistic Source Term White Paper, INL-10-17997, July 2010. [2] David A. Petti, Richard R. Hobbins, Peter Lowry, Hans Gougar, “Representative Source Terms and The Influence of Reactor Attributes on Functional Containment in Modular High Temperature Gas-cooled Reactors,” Nuclear Technology, Vol. 184, p. 181-197, Nov. 2013.

07 ISOTOPE AND RADIATION SOURCES↗
NASA NTRS1979

A model for correlating damage and ultrasonic attenuation in composites

Ultrasonic attenuation data measured by the pulse-echo technique have been obtained versus load for two graphite-epoxy laminates and are presented. The attenuation-load data display distinct changes at load values associated with the initiation of transverse cracks in the weakest plies and with the knee in the bilinear stress-strain curve. It was intended to determine if a simplified diffraction model could be used to account for the attenuation changes. The basis of the diffraction model is that transverse cracks or delaminations, or both, initiating and progressing in the specimen under the area covered by the ultrasonic transducer would cause diffraction of the ultrasonic beam and changes in the measured attenuation.

Hayford, D. T.↗
DOE OSTI2026

The stellar origins of 96 Zr excesses in presolar graphites from the Murchison meteorite

Context. Zirconium-96 is a stable isotope that can be synthesized under different neutron-rich nucleosynthetic conditions. Astrophysical models predict its production to occur in various stellar environments: from low-to-intermediate-mass asymptotic giant branch (AGB) stars to massive stars and core-collapse supernovae. Aims. Detections of 96Zr excesses, in combination with other isotopic measurements from presolar grains can provide unique constraints on its stellar origin. Presolar grains are microscopic particles found in primitive Solar System materials, which formed in stellar winds and supernova ejecta. The isotopic composition of each grain can provide us a snapshot of the nucleosynthetic processes that took place during the parent star’s lifetime. Methods. In this study, we measured the stable isotopes of C, N, O, Mo, Zr, and Ru in high-density presolar graphite grains from the Murchison meteorite and found four grains that contain positive isotopic anomalies in 96 Zr carried by their internal subgrains. We analyzed multi-element isotopic datasets from each grain to explore the source of the observed 96 Zr excesses. Results. Comparisons with stellar models indicate that two grains likely condensed in an intermediate-mass AGB star with initial metallicity of Z ≤ 0.014. Their 96 Zr/ 94 Zr ratios also match those predicted for born-again AGB stars undergoing a very late thermal pulse and rapidly accreting white dwarfs. After comparing the relative populations of the aforementioned dust-producing stars, we propose rapidly accreting white dwarfs as a new, and more likely, stellar source for one of the presolar grains. The remaining two grains could have originated in the supernova ejecta of massive stars, due to correlated excesses in the p-nuclides, 92, 94 Mo. Thus, grains with 96 Zr anomalies can have a variety of stellar origins, in agreement with theoretical studies. Conclusions. Our study highlights the importance of multi-element analysis in constraining the types of stars where presolar grains have condensed. These data will help improve our understanding of various nucleosynthesis processes in different stellar phases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Toward model-free temperature diagnostics of warm dense matter from multiple scattering angles

Warm dense matter plays an important role in astrophysical objects and technological applications, but the rigorous diagnostics of corresponding experiments is notoriously difficult. In this work, we present a model-free analysis of x-ray Thomson scattering (XRTS) measurements on isochorically heated graphite obtained at the Linac Coherent Light Source at multiple scattering angles. We demonstrate that the recent imaginary-time thermometry technique works for scattering data that have been measured in both forward and backward scattering geometry. This opens up the way toward a rigorous quantification of nonequilibrium effects in future experiments, where XRTS measurements are being obtained from multiple scattering angles from the same sample.

Equations of state↗
NASA NTRS1986

An analysis on the shapes of ultraviolet extinction curves. I - The 2175 A bump

IUE data are used to investigate the properties of the 2175 A interstellar extinction bump toward 45 reddened Milky Way stars. Using analytic fitting functions, parameters are derived which describe the central position, FWHM, and strength of the bump. The bump position is very stable, with an extreme deviation which is still significantly larger than the measurement uncertainty, indicating true variability. The bump FWHM has a large range of intrinsic values, 0.77/micron to 1.25/micron. The normalized bump strength varies by more than a factor of three. There are no convincing correlations among the three bump parameters, but the width of the bump is strongly correlated with the dust grain environments. Dense quiescent regions yield broad bumps, while diffuse interstellar medium and regions of recent star formation yield narrower bumps. The absence of any correlation between bump width and central position argues strongly against graphite grains in the size range prescribed by the Mathis, Rumpl, and Nordsieck (1977) model as carriers of the bump.

Fitzpatrick, E. L.↗
NASA NTRS2015

An Internal Thermal Environment Model of an Aluminized Solid Rocket Motor with Experimental Validation

Due to the severity of the internal solid rocket motor (SRM) environment, very few direct measurements of that environment exist; therefore, the appearance of such data provides a unique opportunity to assess current thermal/fluid modeling capabilities. As part of a previous study of SRM internal insulation performance, the internal thermal environment of a laboratory-scale SRM featuring aluminized propellant was characterized with two types of custom heat-flux calorimeters: one that measured the total heat flux to a graphite slab within the SRM chamber and another that measured the thermal radiation flux. Therefore, in the current study, a thermal/fluid model of this lab-scale SRM was constructed using ANSYS Fluent to predict not only the flow field structure within the SRM and the convective heat transfer to the interior walls, but also the resulting dispersion of alumina droplets and the radiative heat transfer to the interior walls. The dispersion of alumina droplets within the SRM chamber was determined by employing the Lagrangian discrete phase model that was fully coupled to the Eulerian gas-phase flow. The P1-approximation was engaged to model the radiative heat transfer through the SRM chamber where the radiative contributions of the gas phase were ignored and the aggregate radiative properties of the alumina dispersion were computed from the radiative properties of its individual constituent droplets, which were sourced from literature. The convective and radiative heat fluxes computed from the thermal/fluid model were then compared with those measured in the lab-scale SRM test firings and the modeling approach evaluated.

Martin, Heath T.↗
NASA NTRS1958

Experimental Ablation Cooling

The film shows ablation tests on Teflon, nylon, a 27 percent phenolic resin, Haveg Rocketon, and graphite. Teflon hemisphere-shaped and flat face noses were tested with laboratory-scale ceramic-heated, pilot-model ceramic-heated, and electric-arc-powered air jets. Nylon hemisphere-shaped noses were tested with laboratory-scale ceramic-heated and electric-arc-powered air jets. Phenolic resin hemisphere-shaped noses were tested with laboratory-scaled ceramic-heated air jets. Haveg Rocketon and graphite hemisphere-shaped noses were tested with electric-arc-powered air jets.

Cooling↗
DOE OSTI2024

Interplay of intercalation dynamics and lithium plating in monolithic and architectured graphite anodes during fast charging

Fast charging of high-capacity anodes is challenging due to lithium plating reactions, which lead to poor cycling performance and safety concerns. Thus, accurate predictions of plating onset and an understanding of this electrochemical process are crucial for robust battery design. However, the most commonly used models, based on porous electrode theory (e.g., the pseudo-2D model), are notoriously difficult to calibrate due to their complexity, limiting their predictive power. This work studies the process of lithium plating during fast charging of (small-particle) graphite half-cells by measuring local reaction progression and plating behavior using optical operando techniques. These experiments employ a realistic 1D graphite electrode geometry with commercially-relevant mass loading charged at fast charge rates. It is demonstrated that the local reaction progression and plating onset can not only be predicted accurately with a p2D numerical model, but that these processes follow a simple scaling law. Remarkably, the entire reaction histories of different electrodes charged at different rates (e.g., 160 μm thickness at 0.5C, 111 μm at 1C or 66 μm at 4C) were observed to have self-similar intercalation profiles. It is demonstrated that plating onset is in turn governed by the reaction profile which explains why both processes exhibit the same scaling behavior. Finally, operando measurements of local reaction dynamics are conducted for the first time in electrodes with channeled architectures, quantitatively determining how channels affect reaction uniformity and plating onset. Together, these results reveal underlying simplicity in the complex electrochemical environment of fast charging and lithium plating, improving understanding of this process. These fundamental insights are broadly applicable for design processes, modeling and experimental evaluation of lithium ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Validation of the ERO2.0 code using W7-X and JET experiments and predictions for ITER operation

Abstract The paper provides an overview of recent modelling of global material erosion and deposition in the fusion devices Wendelstein 7-X (W7-X), JET and ITER using the Monte-Carlo code ERO2.0. For validating the modelling tool in a three-dimensional environment, W7-X simulations are performed to describe carbon erosion from the graphite test divertor units, which were equipped in operational phase OP 1.2 and analysed post-mortem. Synthetic spectroscopy of carbon line emission is compared with experimental results from the divertor spectrometer measurement system, showing a good agreement in the e-folding lengths in the radial intensity profiles of carbon. In the case of metallic wall materials, earlier modelling of the Be/W environment in JET and ITER is revisited and extended with an updated set of sputtering and reflection data, as well as including the mixing model for describing the Be/W dynamics in the divertor. Motivated by recent H/D/T isotope experiments in JET, limited and diverted configuration pulses are modelled, showing the expected trend of both Be and W erosion increasing with isotope mass. For the JET diverted configuration pulses, it is shown that Be migrates predominantly to the upper part of the inner divertor where it initially leads to strong W erosion. With longer exposure time, the growth of a Be deposited layer leads to a reduction of W erosion in that region. A similar trend is observed in simulations of the ITER baseline Q = 10 scenario, however with a more symmetric Be migration pattern leading to deposition also on the outer divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Understanding the Competition between Alcohol Formation and Dimerization during Electrochemical Reduction of Aromatic Carbonyl Compounds

The electrochemical reductive dimerization of small aromatic carbonyl compounds derived from lignocellulosic biomass is a crucial C−C coupling reaction for upgrading small molecules to longchain hydrocarbons, particularly in the synthesis of drop-in sustainable aviation fuels. Although other electrochemical reduction reactions of these reactants (i.e., hydrogenation and hydrogenolysis) have undergone extensive mechanistic investigation, the understanding of dimerization remains relatively underdeveloped. Most importantly, there is a lack of understanding of the selectivity-determining step between dimerization and monomer reduction and critical factors that can affect this step. In this study, we provide a comprehensive mechanistic model to explain the competition between dimerization and monomer reduction of benzaldehyde under various conditions. Our model proposes that the selectivity between dimerization and monomer reduction depends on the competition between desorption of a ketyl radical from the electrode, necessary for dimerization, and further reduction of the ketyl radical to an alcohol on the electrode by proton-coupled electron transfer (PCET). Computationally comparing the adsorption/desorption energy and PCET activation barrier energy is challenging because conventional DFT calculations substantially underestimate the PCET kinetic barriers. In this study, we employed constrained DFTbased configuration interaction (CDFT-CI) to obtain a reliable comparison of these energies. Our mechanistic model was tested and supported by experimental results obtained with four electrodes (Cu, Pb, Bi, graphite), three pH conditions (acidic, neutral, basic), and three potentials. Our study offers a coherent mechanistic foundation that can explain how each of these conditions impacts the desorption and PCET processes and the selectivities for dimerization and alcohol production.

09 BIOMASS FUELS↗
DOE OSTI2026

Understanding carbon sourcing and transport originating from the helicon antenna surfaces during high-power helicon discharge in DIII-D Tokamak

The high-power helicon wave system in the DIII-D tokamak could potentially introduce new plasma–material interaction (PMI) challenges owing to rectified RF sheath potentials that develop near the antenna and surrounding plasma-facing components. We present the first application of the STRIPE (Simulated Transport of RF Impurity Production and Emission) framework to helicon-induced PMIs, extending previous STRIPE studies of ICRH antennas by incorporating net erosion, local re-deposition, and three-dimensional global impurity transport. The integrated workflow couples SOLPS, COMSOL, RustBCA, GITR, and GITRm to simulate carbon erosion, re-deposition, and impurity transport for two experimentally constrained DIII-D H-mode helicon operating scenarios with different antenna–plasma gaps, coupled RF powers, and edge plasma conditions. COMSOL predicts rectified RF sheath potentials of 1–5 kV localized near the lower portion of the antenna, where the magnetic field intersects the surface at grazing incidence. Carbon self-sputtering dominates the erosion source, whereas RF-accelerated D + ions contribute approximately 1% of the total gross erosion. The smaller-gap operating scenario exhibits substantially stronger gross erosion, enhanced local re-deposition (∼12%), and a larger confined carbon inventory owing to increased plasma accessibility and broader RF sheath coverage. Comparison with available DIII-D measurements shows no distinct change in the global carbon signal that correlates with the helicon RF pulse, consistent with the simulations indicating that the helicon-generated carbon source remains small compared with the existing background carbon inventory under the present graphite-wall operating conditions. These results demonstrate the capability of STRIPE to integrate RF sheath modeling, plasma transport, surface interaction physics, and three-dimensional impurity transport for the interpretation of helicon-induced PMIs. The study further identifies the principal sources of modeling uncertainty, including grazing-angle RF sheath physics, slow-wave resolution, plasma-background extrapolation, and trace-impurity transport assumptions, providing a framework for future validation and model development.

Kumar, Atul [Oak Ridge National Laboratory (ORNL),↗
NASA NTRS1990

Chemical interactions in the aluminum-carbon and aluminum-silicon carbide systems

XPS was used to investigate the influence of O2 and H2O on the formation of aluminum carbide at Al-C and Al-SiC interfaces. It was determined that dosing the interfaces with H2O catalyzed the formation of aluminum carbide in both systems. This result is consistent with the oxidation model of carbide formation (developed to understand the kinetics of aluminum carbide formation at graphite-aluminum interfaces). These results imply that the formation of aluminum carbide in graphite- and SiC-reinforced metal-matrix composites, which severely degrades the composite mechanical properties by degrading the fiber and interface strength is catalyzed.

Maruyama, Benji↗
NASA NTRS1993

Global/local interlaminar stress analysis of a grid-stiffened composite panel

A global/local procedure for the computation of the interlaminar stress components at the skin wrap, skin core, and wrap core interfaces for an advanced concept stiffened panel, is described. The procedure consists of a global model of two dimensional shell elements that is used to design a grid stiffened panel with blade type stiffeners, a local model of three dimensional solid elements that is used to compute interlaminar stress components, and a scheme devised to assign displacement boundary conditions for a local model that are based on displacement and rotation data of a few nodes of the global model. A global panel was designed according to strength, stiffness, and stability criteria associated with the design of traditional aircraft wing panels. Interlaminar normal and shearing stress components, computed via the local model, were found to be well below typical tensile normal and shearing strengths of a graphite epoxy material.

Wiggenraad, J. F. M.↗
NASA NTRS1999

Evaluation of a Scaling Approach for Stiffened Composite Flat Panels Loaded in Compression

Results of an experimental study of the effects of geometric scaling on a graphite-epoxy flat stiffened-skin panel concept loaded in compression are presented. The scaled models were fabricated using an approach to laminate thickness scaling, referred to as the "ply-thickness method." The structural response and failure characteristics of full-, half-, and quarter-scale specimens fabricated with this scaling concept are discussed and compared. The scaled up failure loads of the half- and quarter-scale models agree well with the failure load of the full-scale prototypes tested. The strength scale effects observed in the results were substantially smaller than those in earlier studies that employed other laminate scaling techniques. The experimental results indicate that failure initiated by crippling and delamination in the cap of the longitudinal stiffeners and skin-stiffener separation at the bond-lines, and was followed by failure of the specimen skin. Nonlinear finite element analysis results correlate well with experimental results up to failure.

Rouse, Marshall↗
NASA NTRS1998

Modeling of the Multiparameter Inverse Task of Transient Thermography

Transient thermography employs preheated surface temperature variations caused by delaminations, cracks, voids, corroded regions, etc. Often, it is enough to detect these changes to declare a defect in a workpiece. It is also desirable to obtain additional information about the defect from the thermal response. The planar size, depth, and thermal resistance of the detected defects are the parameters of interest. In this paper a digital image processing technique is applied to simulated thermal responses in order to obtain the geometry of the inclusion-type defects in a flat panel. A three-dimensional finite difference model in Cartesian coordinates is used for the numerical simulations. Typical physical properties of polymer graphite composites are assumed. Using different informative parameters of the thermal response for depth estimation is discussed.

Plotnikov, Y. A.↗
NASA NTRS1991

Evaluation of some scale effects in the response and failure of composite beams

The feasibility of using scale model testing for predicting full-scale behavior of composite beams loaded in tension and flexure was investigated. Classical laws of similitude were applied to fabricate and test replica model beams to identify scaling effects in the load response, strength, and mode of failure. Experiments were conducted using graphite-epoxy composite beams having different laminate stacking sequences and a range of scaled sizes. Results indicated that the elastic response of scaled composite beams was independent of specimen size. However, a significant scale effect in strength was observed. In addition, a transition in failure mode was observed among scaled beams of certain laminate stacking sequences. Weibull statistical and fracture mechanics based models were applied to predict the strength scale effect since standard failure criteria cannot account for the influence of absolute specimen size in failure.

Jackson, Karen E.↗
NASA NTRS1992

Scale effects in the response and failure of fiber reinforced composite laminates loaded in tension and in flexure

The feasibility of using scale model testing for predicting the full-scale behavior of flat composite coupons loaded in tension and beam-columns loaded in flexure is examined. Classical laws of similitude are applied to fabricate and test replica model specimens to identify scaling effects in the load response, strength, and mode of failure. Experiments were performed on graphite-epoxy composite specimens having different laminate stacking sequences and a range of scaled sizes. From the experiments it was deduced that the elastic response of scaled composite specimens was independent of size. However, a significant scale effect in strength was observed. In addition, a transition in failure mode was observed among scaled specimens of certain laminate stacking sequences. A Weibull statistical model and a fracture mechanics based model were applied to predict the strength scale effect since standard failure criteria cannot account for the influence of absolute specimen size on strength.

Jackson, Karen E.↗