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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Phase change material-to-refrigerant heat exchangers: Experimental validation and uneven melting analysis

Thermal energy storage (TES) using phase change materials (PCMs) enables load shifting and reduces mismatches between the building's thermal demand and the heat pump (HP) system's thermal capacity. While most TES studies have focused on single-phase heat transfer fluids, PCM-to-refrigerant heat exchangers (PRHX), in which the refrigerant undergoes liquid-vapor phase change, remain largely underexplored. Here, to address this gap, this study systematically investigates a shell-and-tube PRHX operating as the condenser in an HP through experiments and simulations. A two-dimensional enthalpy-based finite-volume PRHX model was developed and validated against experimental data, including refrigerant and PCM temperature profiles, PRHX capacity, and system-level performance. Results show three characteristic stages: Stage I, the condenser temperature rapidly increases with PCM in the solid state absorbing sensible heat; Stage II, stable operation during PCM phase change; and Stage III, performance decline once part of the PCM becomes fully melted. Heat transfer analysis revealed that uneven PCM melting along the condenser length was driven primarily by variations in the refrigerant-side heat transfer coefficient, rather than local approach temperature differences. Further study demonstrated that minimizing refrigerant outlet subcooling led to more uniform PCM melting and extended Stages I and II by up to 89%. This paper provides a validated PRHX model, clarifies the mechanisms of uneven PCM melting, and highlights subcooling control as an effective strategy to improve TES performance in HP systems without requiring secondary loops.

Graphite matrix↗
DOE OSTI2025

Experiment-modeling studies comparing energy dissipation in the DIII-D SAS and SAS-VW divertors

Recent DIII-D experiments on Small Angle Slot (SAS) divertors have confirmed that a combination of divertor closure and target shaping can enhance cooling across the divertor target and increase energy dissipation, but with significant dependence on B T (toroidal magnetic field) direction. In these novel divertors, the roles of closure, target shaping, drifts, and scale lengths are all interconnected in optimizing dissipation, with the separatrix electron density n eSEP being the key parameter associated with the level of dissipation/detachment. After modifying the original flat-targeted graphite SAS to include a V shape with a tungsten coating on the outer side of the divertor (SAS-VW), matched series of discharges were run to compare to detailed SOLPS-ITER modeling. Experimentally, when run as designed with the outer strike point at the slot vertex, SAS-VW requires nearly identical n eSEP for detachment as the original SAS, with little difference in dissipation for the new geometry. This is in contrast to (1) earlier modeling predictions that a small change of the SAS geometry to a V shape should enhance dissipation at the same n eSEP for magnetic configurations having better H-mode access (ion B × ∇B drift directed into the divertor), and (2) despite the achievement of significantly higher (2-7x) neutral pressures and compression in the SAS-VW slot. Comparisons of experimental density scans to the most recent SOLPS-ITER modeling with ExB drifts show reasonable agreement for dissipation/detachment onset when using separatrix density as the independent parameter. In order to help understand the discrepancy in modeled vs actual performance for the new configuration, additional measurements varying gas injection location and impurity injection were undertaken. In-slot D 2 gas fueling is more effective (5–22 %) in promoting detachment, in accord with modeling. In-slot impurity injection (N 2 or Ne) can yield 30 % lower core Z eff and 15 % less confinement degradation after detachment compared to main chamber puffing, as well as relatively lower tungsten leakage from the divertor. Modeling can also reproduce the improved detachment seen as the strike point moves inboard of the slot vertex. While we can explain the effects of the most important parameters causing energy dissipation in these slot divertors, it remains that many aspects of their behavior cannot be accurately modeled using state-of-art codes such as SOLPS-ITER. This is of concern for future model-driven designs utilizing similar V-shaped geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Design and performance evaluation of a dual-circuit thermal energy storage module for air conditioners

We present experimental results and a validated numerical model of a dual-circuit phase-change thermal energy storage module for air conditioners. The module incorporates a phase-change material encapsulated in compressed expanded natural graphite foam. We used n-tetradecane as the PCM with a transition temperature (~4.5 °C) suitable for air-conditioning applications. Heat exchange to and from the module is accomplished through two fluid loops operating as a heat source and sink embedded inside multiple slabs of the composite material. This dual-circuit design enables easier integration with air-conditioning equipment and provides enhanced flexibility in system operation as compared to the state-of-the-art thermal storage systems. When integrated with an air-conditioner, this design will enable peak-load shaving and enhances operational efficiency. The thermal storage device was designed for a nominal storage capacity of ~ 3.5 kWh. We evaluated the heat transfer and energy storage performance of this device using standalone heat transfer experiments to estimate key thermal resistances and identify design improvements before integration with an air conditioner. The numerical model of the heat exchanger uses a combination of discretized and lumped parameter approaches to maintain a balance between accuracy and computational expense. Our analyses show that the geometric features and integration of fluid tubes are key contributors to the thermal contact resistance between the fluid and the thermal storage material, and consequently, to the overall performance of the thermal storage module. Our standalone experiments also identified important operating scenarios in which this thermal storage module can be used for air-conditioning in buildings.

25 ENERGY STORAGE↗
DOE OSTI2021

Unraveling compacted graphite evolution during solidification of cast iron using in-situ synchrotron X-ray tomography

In spite of many years of research, the physical phenomena leading to the evolution of compacted graphite (CG) during solidification is still not fully understood. In particular, it is unknown how highly branched CG aggregates form and evolve in the semi-solid, and how local microstructural variations at micrometer length scale affect this growth process. We present here the first time-resolved synchrotron tomography combined with a bespoke high-temperature environmental cell that allows direct observation of the evolution of CG and relates this dynamic process to the local surrounding microstructures in a cast iron sample during repeated melting and solidification. Distinct processes are identified for the formation of CG involving the nucleation, growth, development of branches and interconnection of graphite particles, ultimately evolving into highly branched graphite aggregates with large sizes and low sphericities. CG is found to nucleate with a spheroidal or a plate-like shape, developing branches induced by high carbon concentration, e.g. thin melt channels. Additionally, CG grows much faster than spheroidal graphite during subsequent cooling in solid state. The direct visualization of the dynamic solidification process provides unprecedented new insights into formation mechanisms of CG and correlating factors such as local microstructural variations, and guides the development of CG iron solidification models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A combined SECM and electrochemical AFM approach to probe interfacial processes affecting molecular reactivity at redox flow battery electrodes

Redox flow batteries are attractive technologies for grid energy storage since they use solutions of redox-active molecules that enable a superior scalability and the decoupling of power and energy density. However, the reaction mechanisms of the redox active components at RFB electrodes are complex, and there is currently a pressing need to understand how interfacial processes impact the kinetics and operational reversibility of RFB systems. Here, we developed a combined electrochemical imaging methodology rooted in scanning electrochemical microscopy (SECM) and atomic force microscopy (AFM) for exploring the impact of electrode structure and conditioning on the electron transfer properties of model redox-active dialkoxybenzene derivatives, 2,5-di-tert-butyl-1,4-bis(2-methoxyethoxy)benzene (C1) and 2,3-dimethyl-1,4-dialkoxybenzene (C7). Using AFM and secondary-ion mass spectrometry (SIMS), we observed the formation of interfacial films with distinct mechanical properties compared to those of cleaved graphitic surfaces, and exclusively during reduction of electrogenerated radical cations. These films had an impact on the median rate and distribution of the electron transfer rate constant at the basal plane of multilayer and single layer graphene electrodes, displaying kinetically-limited values that did not yield the activation expected per the Butler–Volmer model with a transfer coefficient ~0.5. These changes were dependent on redoxmer structure: furthermore, SECM showed strong attenuation of C7 kinetics by a surface layer on MLG and SLG, while C1 kinetics were only affected by SLG. SECM and AFM results together show that these limiting films operate exclusively on the basal plane of graphite, with the edge plane showing a relative insensitivity to cycling and operation potential. Ultimately, this integrated electrochemical imaging methodology creates new opportunities to understand the unique role of interfacial processes on the heterogeneous reactivity of redoxmers at electrodes for RFBs, with a future role in elucidating phenomena at high active concentrations and spatiotemporal variations in electrode dynamics.

25 ENERGY STORAGE↗
NASA NTRS1990

Interlaminar shear stress effects on the postbuckling response of graphite-epoxy panels

The influence of shear flexibility on overall postbuckling response was assessed, and transverse shear stress distributions in relation to panel failure were examined. Nonlinear postbuckling results are obtained for finite element models based on classical laminated plate theory and first-order shear deformation theory. Good correlation between test and analysis is obtained. The results presented analytically substantiate the experimentally observed failure mode.

Engelstad, S. P.↗
NASA NTRSDate not supplied

Toughening of Boron Carbide Composites by Hierarchical Microstructuring

Due to a unique combination of properties including high hardness, low density, chemical and thermal stability, semi-conductivity, and high neutron absorption, boron carbide (B4C) is a potential candidate for various applications involving extreme environment. However, B4C’s current application is limited because of its low fracture toughness. In this study, a hierarchical microstructure design with features including TiB2 grains and graphite platelets was used to toughen B4C by simultaneously utilizing multiple toughening mechanisms including crack deflection, bridging, and micro-crack toughening. Using field-assisted sintering technology (FAST), B4C composites with dense and hierarchical microstructure were fabricated. Previously, the fracture toughness of fabricated B4C composites was measured at micro-scale using micro-indentation to have up to 56% improvement. In this work, the B4C composites’ fracture toughness was characterized at macro-scale using four-point bending methods and compared with previous results obtained at micro-scale. Micromechanics modeling of fracture behaviors for B4C-TiB2 composites was also performed to evaluate the contributions from experimentally observed toughening mechanisms. From four-point bending tests, B4C composites reinforced with both TiB2 grains (~15 vol%) and graphite platelets (~8.7 vol%) exhibited the highest fracture toughness enhancement from 2.38 to 3.65 MPa∙m^(1/2). The measured values were lower than those obtained using micro-indentation but maintained the general trends. The discrepancy between the indentation and four-point bending test results originated from the complex deformation behaviors triggered by the high contact load during indentation tests. Through micromechanics modeling, introduced thermal residual stress due to thermal expansion mismatch between B4C and TiB2, and weak interphases at B4C-TiB2 boundaries were identified as the main causes for experimentally observed toughness enhancement. These results proved the effectiveness of hierarchical microstructure designs for B4C toughening and can provide reference for the future design of B4C composites with optimized microstructures for further fracture toughness enhancement.

composites↗
DOE OSTI2021

Time-Resolved X-ray Operando Observations of Lithiation Gradients across the Cathode Matrix and Individual Oxide Particles during Fast Cycling of a Li-Ion Cell

Lithiated transition metal oxides serve as active materials in the positive electrode (cathode) of lithium-ion cells. During electrochemical cycling, lithium ions intercalate and deintercalate into these oxide particles. This behavior causes two types of lithiation gradients to emerge: (i) a bulk gradient across the depth of the cathode matrix (averaged over individual oxide particles) and (ii) a microscopic gradient across the particles themselves, which also depends on their location in the electrode. Here we show how both gradients can be studied using operando X-ray diffraction during 4C charge and 4C discharge. The oxide (de)lithiation is estimated from the unit cell parameters by indexing the X-ray diffraction spectra. By fitting the lithiation profiles with orthogonal polynomials, the bulk gradients across the electrode thickness are quantified. These gradients develop as the current flows through the cell and dissipate during open-circuit and potentiostatic-hold periods. Further details of lithiation dynamics can be obtained through shape analysis of the Bragg peaks. In particular, from electrochemical model simulations, we show that the width and skewness of the (003) peak track (de)lithiation fronts moving across the individual oxide particles.

(003) peak↗
DOE OSTI2025

Summary of the 1st Reactor Graphite Workshop, 8-9 July, 2025

Graphite is widely used as moderator in many fission reactors and other criticality systems. The understanding of its short- and long-term behavior in fission environments is crucial not only for the operation and optimization of existing reactors but also for the efficient design and safe deployment of future advanced reactors. It is thus of utmost importance to incorporate in nuclear data libraries accurate evaluations of the interaction of thermal neutrons with graphite used in applications. This is however a complex task as there are challenges in defining the detailed characterizations of the studied material, appropriate modeling, and optimal validation suite. This will only be achieved by the concerted effort of experts of all relevant areas.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1981

Progress report on the high resolution spectrograph for the Space Telescope

The instrument design is complete and many components and subsystems have been manufactured and tested. Solar-blind, multichannel, pulse counting detectors for two ultraviolet spectral bands have been developed and satisfactory flight units have been chosen. A large, lightweight graphite-epoxy optical bench, utilizing fingerplate joints with no metallic parts in order to minimize thermal expansion, has been built. Tests of the engineering model torque-motor-driven grating carrousel indicate that this critical system exceeds tight specifications for accurate and repeatable angular positioning. A complement of large, high-frequency gratings and an echelle is being completed in several laboratories. Software resident in a computer on the spacecraft will control the spectrograph with the capability to alter procedures in response to real-time evaluation of the data.

Brandt, J. C.↗
NASA NTRS1995

A 75 Watt, 59 to 64 Ghz Space TWT

The design and performance of a high efficiency 70 Watt, 59 to 64 GHz traveling-wave tube (TWT) is described. This intersatellite communications TWT, designated the Hughes Aircraft Company Model 961HA, employs a coupled-cavity slow-wave structure with a two-step velocity taper and a four-stage isotropic graphite multistage depressed collector (MDC). This TWT achieved outstanding beam transmission, 96% with RF. This is the highest power and widest bandwidth device of this type. The tube has 45 dB of gain and 30% efficiency. This TWT is conduction cooled and weighs 13 lbs.

Limburg, Helen↗
NASA NTRS1989

Infrared emission from ultracompact H II regions

Models of circumstellar dust shells around ultracompact (UC) H II regions were constructed that accurately fit the observed IR flux distributions. The models assume spherically symmetric dust shells illuminated by stars whose bolometric luminosity is inferred from the integrated FIR flux densities. Assuming ionization by a single zero age main sequence (ZAMS) star, the relations of Panagia were used to infer the stellar radius and effective temperature for a given luminosity. The grain mixture in the dust shell consists of bare graphite and silicate grains with the optical properties of Draine and Lee and the size distribution of Mathis et al. The computer code of Wolfire et al was used to solve the radiative transfer equations through a spherical dust shell. The model provides monochromatic luminosities, dust temperatures, and opacities through the shell. Aside from the stellar and dust properties, the only other input parameters to the model are the distance to the shell, the form of its density distribution, and its outer radius. Predictions of the model are compared with observations of a typical UC H II region and the run of dust temperature with radius and the optical depth with frequency are discussed.

Churchwell, ED↗
NASA NTRS1986

Advanced composites for large Navy spacecraft

An overview is given of work conducted on contract for the Naval Sea Systems Command. The objective of this contract was to provide direction for the development of high modulus graphite reinforced metal matrix composites. These advanced materials can have a significant effect on the performance of a spacecraft before, during and after an evasive maneuver. The work conducted on this program was organized into seven technical tasks. Task 1 was development of a generic Navy spacecraft model. Finite element models of candidate structural designs were developed. In Task 2, the finite-element model(s) of the structure were used to conduct analytical assessments involving conventional materials, resin matrix composites and metal matrix composites (MMC). In Task 3 and 4, MMC material design, fabrication and evaluation was conducted. This consisted of generating material designs and developing a data base for a broad range of graphite reinforced MMC materials. All material was procured according to specifications which set material quality and material property standards. In Task 5, a set of evasive maneuvering requirements were derived and used in Task 6 to conduct analytical simulations. These analytical simulations used current SOA material properties and projected material properties to provide an indication of key payoffs for material development. In Task 7, a set of material development recommendations was generated.

Davis, William E.↗
NASA NTRS2009

Status of the NEXT Long-Duration Test After 23,300 Hours of Operation

The NASA s Evolutionary Xenon Thruster (NEXT) program is developing the next-generation ion propulsion system with significant enhancements beyond the state-of-the-art in ion propulsion to provide future NASA science missions with enhanced mission capabilities at a low total development cost. As part of a comprehensive thruster service life assessment utilizing both testing and analyses, a Long-Duration Test (LDT) was initiated in June 2005, to verify the NEXT propellant throughput capability to a qualification-level of 450 kg, 1.5 times the anticipated throughput requirement of 300 kg per thruster from mission analyses. The LDT is being conducted with a modified, flight-representative NEXT engineering model ion thruster, designated EM3. As of July 2009, the thruster has accumulated 23,300 h of operation with extensive durations at the following input powers: 6.9, 4.7, 1.1, and 0.5 kW. The thruster has processed 427 kg of xenon surpassing the NSTAR propellant throughput demonstrated during the extended life testing of the Deep Space 1 flight spare ion thruster and approaching the NEXT development qualification throughput goal. The NEXT LDT has demonstrated a total impulse of 16.0 10(exp 6) N/s; the highest total impulse ever demonstrated by an ion thruster. Thruster performance tests are conducted periodically over the entire NEXT throttle table with input power ranging 0.5 to 6.9 kW. Thruster performance parameters including thrust, input power, specific impulse, and thruster efficiency have been nominal with little variation to date. The NSTAR first-failure mode, accelerator aperture erosion leading to electron backstreaming, has been mitigated in the NEXT design. The severe NSTAR discharge cathode assembly erosion has been mitigated by a graphite keeper in the NEXT thruster. Tracking of the NEXT first failure mode, charge-exchange ion impingement on the accelerator grid causing hexagonal groove erosion, is consistent with model predictions and indicates thruster life greater than or equal to 750 kg throughput. This paper presents the status, performance data, and wear characteristics of the NEXT LDT to date.

Herman, Daniel A.↗
DOE OSTI2025

Layers Can Be Deceiving: A Hopping Model for Small Molecule Diffusion in TATB Crystal

Sorption of small molecule gases in materials can play a significant role in their long‐term stability and compatibility within multi‐material assemblies. While many material properties of the insensitive high explosive TATB (1,3,5‐triamino‐2,4,6‐trinitrobenzene) are well understood, very little is known regarding its permeability to gases. TATB crystal exhibits a graphitic‐like layered packing structure that evokes a mental schema in which the layers form nanoscopic channels, but it is unclear whether this structure promotes gas transport. Here, we use molecular dynamics (MD) simulations to predict transport of small molecules through TATB single crystal. An approach to fit classical force fields is developed to model TATB interactions with H 2 O, He, Ne, and Ar, which is then combined with steered MD to probe gas transport along selected directions in the crystal. We find that small molecule transport occurs via a hopping mechanism that exhibits distinct jumps between interstitial sites and is substantially faster normal to the layers as compared to through them. This result stems from the finding that intralayer junctions between adjacent TATB molecules are the most stable interstitial sites and that energetic barriers are lower for hopping between adjacent layers. An empirical model for diffusion rate based on the MD data shows that the rate decays exponentially with increasing molecular radius and is negligibly small for all molecules larger than He, including common atmospheric gases. These findings have implications for the interpretation of experiments that measure surface area, material response to extreme conditions, and are expected to help constrain models for material aging.

36 MATERIALS SCIENCE↗
NASA NTRS1986

An in-situ measurement of particulates from solid rocket motors fired in space

Current models exist that predict the damage caused by the impact of aluminum oxide exhaust particles as well as their lifetime in useable space. In these models, two necessary inputs are the size and flux of the particles. An experiment, referred to as the Plume Witness Plate, was designed for the Remote Manipulator System of the space shuttle orbiter to measure in-situ the flux and material effects of a solid rocket motor (SRM) firing in space. Five different types of samples were used to provide a broad range of substances: (1) fused quartz glass (representative of orbiter windows); (2) germanium micrometeroid capture cells; (3) orbiter HRTS tiles from the thermal protection system; (4) Kapton foil; and (5) metallic disks of aluminum, copper, titanium, graphite epoxy, and gold. The analyses of the data show excellent agreement with ground-based SRM firings in terms of particle size distribution and mass distribution. The Particle Impact Damage Integrator computer model used to calculate potential damage of orbiter surfaces by SRM exhaust plumes agrees favorable with the results in terms of particle size and velocity distributions though it may be conservative by as much as 20%.

Alred, J. W.↗
DOE OSTI2024

Embedding Neural Thermal Scattering (NeTS) Modules in SERPENT for Higher Fidelity Advanced Reactor Analysis

When a neutron born in fission thermalizes to the order of $k$ $B$ $T$, it’s de-Broglie wavelength and energy approach the order of inter-atomic spacing and elementary lattice oscillations, respectively. $S$($a,β,t$) or the scattering law, uuantify these temperature-dependent crystallographic contributions to total cross section (or reaction rate). In a Monte Carlo analysis, cumulative distribution functions (CDFs) of $S$($a,β,t$) are loaded to memory from “A Compact ENDF” (ACE) files for stochastically selecting thermal scattered neutron trajectories. In this work, novel neural thermal scattering (NeTS) modules for $S$($a,β,t$) CDFs are designed, trained, serialized and embedded within SERPENT using Python’s limited C-API for on-the-fly deployment of crystalline graphite $S$($a,β,t$) sampling. Torchscript tracing and Numba just-in-time (JIT) compilation streamline neural inference on NVIDIA GPUs with CUDA libraries. Demonstrations of bare sphere thermalization of fast and thermal sources show excellent agreement between embedded NeTS in SERPENT and MCNP. With an explicit model of the reactor, NeTS can predict on-the-fly changes in TREAT neutron spectra as a function of local temperature, which can serve to improve transient and accident predictions in a multiphysics analysis framework. This framework can be further extended to account on-the-fly for changes in local graphitic microstructure to scattering cross sections, and outlines a novel coupling of modern machine learning with state-of-the-art reactor physics methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Assessing Electrolyte Fluorination Impact on Calendar Aging of Blended Silicon-Graphite Lithium-Ion Cells Using Potentiostatic Holds

Silicon-based lithium-ion batteries have started to meet cycle life metrics, but they exhibit poor calendar life. Here, electrolyte fluorination impact on calendar fade of blended silicon-graphite anodes is explored using a LiPF 6 in EC:EMC:FEC electrolyte vs LiBOB in EC:EMC electrolyte. We utilize a combined experimental-modeling approach applying potentiostatic voltage holds (V-hold) to evaluate electrolyte suitability for calendar life in a shortened testing timeframe (~2 months). Our theoretical framework deconvolutes the irreversible parasitic capacity losses (lithium lost to the solid electrolyte interphase) from the V-hold electrochemical data. Unfluorinated electrolyte (dominant LiBOB reduction) exhibits higher cell resistance as compared to fluorinated electrolyte (dominant FEC reduction). Both systems have similar irreversible capacities during the voltage hold duration with slower rate of parasitic capacity loss for the LiBOB system. Extrapolation of the parasitic losses to end of life capacity fade of 20% shows LiBOB electrolyte outperforming LiPF 6 electrolyte in calendar life. The results demonstrate the applicability of the V-hold protocol as a rapid material screening tool providing semi-quantitative calendar lifetime estimates.

25 ENERGY STORAGE↗