SEARCH · Search NASA

Results for “Stability Analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Grid Integration of Renewable Energy and Energy Storage

Grid integration of renewable energy and energy storage requires forward-looking planning process, and increased emphasizes on reliability, resilience, and equity. Power-electronics based energy generation including solar, wind, distributed energy resources (DERs), and various types of grid-tied energy storage and emerging loads, are reshaping grid operator's understanding on interconnection level performance and responses. This paper will present the ongoing work at PNNL related to power electronics R&D, energy modeling and analysis, and a wide spectrum of grid stability studies and technologies in support of grid integration of renewable energy and energy storage.

Power Electronics, grid integration of renewable e↗

DOE OSTI2025

Risk Assessment of EIC Central Detector (ePIC) Solenoid Magnet (MARCO)

As part of the BNL-JLab-CEA Electron Ion Collider (EIC) collaboration, the design of a 2 T, 2.8 m bore diameter, 3.8 m long conduction cooled superconducting detector magnet design is completed. Such magnet will be employed at the interaction region of ePIC for physics experiment. The magnet is a passive shielded solenoidal magnet system consisting of 3 coils wound with specially designed conductor using NbTi Rutherford type cable in copper stabilized channel. This paper describes the risk analysis as the part of Failure Modes and Effects Analysis (FMEA) that was carried out as a team to identify their various failure modes and risks associated with the magnets system. In conclusion, this FMEA is intended to become the content of the designed document as an integral part of the engineering assessment and the statement of work for the potential vendors towards design and built.

Ghoshal, Probir K. [Thomas Jefferson National Acce↗

DOE OSTI2025

Experimental test of the ratio method for nuclear-reaction analysis

Nuclear halos are exotic quantal structures observed far from stability. They are mostly studied through reactions. The ratio of angular cross sections for breakup and scattering is predicted to be independent of the reaction process and to be very sensitive to the halo structure. We test this new observable experimentally for the first time on the collision of 11 Be on C at 22.8 MeV/nucleon and using existing data on Pb at 19.1 MeV/nucleon. The theoretical predictions are verified, which offers the possibility to develop a new spectroscopic tool to study nuclear structure far from stability.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

EIC ESR beam position monitors thermal and thermomechanical study preliminary design

The Electron Ion Collider (EIC) Electron Storage Ring (ESR) Beam Position Monitors (BPM) are designed to monitor the electron beam orbit position on a turn-by-turn and bunch-by-bunch basis in the vacuum chamber. Large variations in ESR fill pattern and bunch intensity imposes a large dynamic range on the BPM system design. And given the large ESR beam current - up to 2.5 A with 1160 x 20 nC bunches – the heating by synchrotron radiation (SR) and beam impedance is a significant concern. The high energy of the electrons – up to 18 GeV – leads to scat tering of the SR that produces high radiation doses and as sociated heating of structures away from the SR impact area. Because of this intense heating, and stringent position stability requirements, an in-depth engineering analysis was carried out. This paper will review the analysis made of the preliminary design and its main outcomes.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment (Final Report)

This report documents the FY25 Theory and Simulation Performance Target (TSPT) of developing an integrated modeling framework for fusion reactor design and assessment (FREDA). Over Q1-Q4, new capabilities were developed across both plasma and engineering domains and demonstrated on an example representation of a Compact Advanced Tokamak with a Dual Cooled Lead Lithium blanket. This represents a first-of-a-kind demonstration of coupled core-to-wall-to-engineering for a reactor. Self-consistent CESOL workflows were applied to provide core, pedestal, and SOL prediction; new modules were developed for energetic particle stability (FAR3D) and transport (TGLF-EP) analysis; and boundary plasma modeling (SOLPS-ITER, BOUT++/Hermes-3) was expanded to evaluate wall and divertor heat fluxes and interface with engineering thermal analysis. A parameterized CAD tool, TRACER, was expanded to generate medium-fidelity divertor, blanket, and coil geometries; OpenFOAM and Diablo workflows were applied for first-wall and divertor thermal analyses with helium cooling; and reduced-order models were created for high-mass-flux divertor cooling. Magnet multiphysics capabilities were verified between Elmer, Diablo, and a new MFEM-based solver, and workflows enable stress, thermal, and neutron-fluence analysis of TF coils with neutronics-driven heating. Nuclear and blanket analysis workflows were demonstrated, including tritium breeding, transport, and CFD-informed thermo-mechanical assessment. Preliminary multi-fidelity uncertainty quantification workflows were applied to boundary modeling codes and shown to achieve variance reductions with fewer high-fidelity boundary simulations. Key findings highlight the challenges of resolving the ITEP gap to find suitable balance between wall and divertor loads, neutron heating, and practical limits of PFC cooling. Next step priorities are to develop automated workflows to check boundary code convergence and detachment, implement tighter physics-engineering CAD provenance tracking, and inclusion of plasma-material interface models for SLAG and tungsten cracking behavior. Collectively, these developments establish sophisticated capabilities for predictive, multi-fidelity, whole-device modeling that integrates plasma physics, materials, magnets, and nuclear engineering to guide pathways to viable Fusion Pilot Plant design points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Understanding Defect Activation and Kinetics in Next Generation CdTe Absorbers

This project investigated the role of defects and dopants in cadmium telluride (CdTe) solar cells to improve efficiency and long-term stability. Advanced X-ray microscopy and spectroscopy techniques were used to map composition, identify defect structures, and measure charge collection at the nanoscale. These experimental results were combined with computational modeling to understand how dopants such as copper and arsenic interact with other elements during processing. The study found that many dopants become electrically inactive due to the formation of complexes with chlorine and oxygen, limiting carrier concentration and device performance. It also showed that selenium, used to enhance efficiency, can migrate within the material during both manufacturing and operation, altering local structure and potentially affecting stability. The project developed new diagnostic and analysis tools for studying photovoltaic materials and provided a mechanistic understanding of key performance limitations. These insights support the development of more efficient, stable, and cost-effective CdTe solar technologies for large-scale energy deployment.

14 SOLAR ENERGY↗

DOE OSTI2026

Experimental Studies of Nonlinear Integrable Optics with Elliptic Potentials

The stable transport of charged particle beams is the core challenge in designing and operating accelerators. The current paradigm for transverse focusing is based on quadrupole and dipole elements, described as a linear integrable Hamiltonian by Courant and Snyder. The operational limits of high intensity accelerators are often determined by collective instabilities within the beam. Nonlinear elements may be added to suppress certain forms of these collective effects, but restrict the range of stable trajectories. These shortcomings motivate extending to a nonlinear integrable system which can offer the benefits of suppressing collective instabilities without limiting the stable trajectories. Danilov and Nagaitsev proposed a novel nonlinear integrable system with elliptical potentials, which could be implemented as magnetic elements in an accelerator. Such a system has been implemented for practical verification in the integrable optics test accelerator (IOTA), a small storage ring constructed for beam dynamics studies. Electron beam studies in IOTA treat the low-emittance beam as a macroparticle for detailed probing of the expected single particle dynamics. This dissertation details new measurements of the nonlinear integrable system relevant for practical implementation. Turn-by-turn measurements of the kicked beam responses are used for phase space reconstruction and analysis of the predicted nonlinear dynamics. The stability and aperture for various nonlinear configurations are measured with beam losses. Amplitude dependent detuning, a core figure of merit for suppressing instabilities, is measured and compared with high fidelity particle tracking simulations. The synchrotron radiation images of circulating beam allow direct measurements of the topology and lifetime in configurations where nonlinear focusing terms dominate.

Wieland, John [Michigan U.] (ORCID:000000028971852↗

DOE OSTI2025

Characterization of the degradation of gamma-irradiated elastomers using Raman spectroscopy

This report presents key findings from Raman spectroscopic analysis of gamma-irradiated rubber samples extracted from a laminated lead-damped rubber (LDR) seismic isolation device. The samples were exposed to gamma radiation from a 60Co source in a Foss Therapy Services gamma irradiator, reaching absorbed doses up to 1600 kGy. A distinct threshold near 400 kGy was identified, beyond which significant spectral changes were observed. Two Raman peaks - at approximately 425 cm-1 and 2440 cm-1 - were tracked as a function of dose using Gaussian fitting. The 425 cm-1 peak, attributed to sulfur–sulfur (S–S) bond stretching (resulting from vulcanization of the rubber), exhibited a dose-dependent upshift, indicating radiation-induced crosslinking within the sulfur-based polymer network. Conversely, the 2440 cm-1 peak, likely associated with vibrational modes of additives or impurities, showed a downward shift with increasing dose, suggesting chain scission and degradation of non-rubber constituents. These results provide first-of-a-kind insights into the microstructural evolution of elastomers under high-dose gamma irradiation and establish a preliminary dose threshold for significant degradation. Future work will incorporate multi-modal characterization—including Fourier Transform Infrared (FTIR) spectroscopy, scanning electrom microscopy (SEM) of the rubber surface morphology, thermogravimetric analysis (TGA) to determine changes in thermal stability, and mechanical testing—to correlate molecular-level changes with macroscopic performance of these elastomers as damping media in seismic isolation devices. These findings are expected to provide regulatory guidance and design criteria for qualifying low-damping rubber seismic isolators in advanced nuclear reactor applications.

36 - MATERIALS SCIENCE↗

DOE OSTI2025

Impact of the isotope effect on divertor detachment and pedestal structure in DIII-D

Dedicated H-mode hydrogen experiments have been conducted at DIII-D and are compared to similar deuterium plasmas under a range of divertor conditions, revealing significant effects of the isotope mass on divertor and pedestal conditions. H-mode confinement is found to be systematically lower in hydrogen and experiences a worse degradation with increased fueling. This coincides with a higher ne,sep in hydrogen for similar fueling rates. When ne,ped is matched in hydrogen and deuterium as are the ne and Te pedestal structures, the target conditions deviate strongly from each other. An analysis of two major factors in the pedestal structure, pedestal stability and neutral fueling, shows they are not sufficient to explain either the pedestal differences between isotopes or the hydrogen pedestal structure respectively. An analysis of hydrogen and deuterium density ramps demonstrates a 20% higher ion current and a 30% higher upstream density at detachment onset in hydrogen. The characteristics of target heat flux and divertor radiation are shown to change significantly with the isotope mass in each divertor operational regime. The results shown here indicate divertor conditions trend with ion mass in a manner beneficial for divertor performance in future devices.

Mattes, Ray [Univ. of Tennessee, Knoxville, TN (Un↗

DOE OSTI2026

Effects of stabilizer strength on the delamination strength of ReBCO coated conductors

Delamination strengths (DSs) of ReBCO tapes were evaluated using a transverse tensile anvil test. This study first demonstrated that the electromechanical DS is nearly equivalent to the purely mechanical DS for ReBCO tapes exhibiting abrupt delamination behavior. Mechanical DS was then examined for tapes with and without stabilizer edges to assess the influence of stabilizer strength. Scanning electron microscopy/energy-dispersive x-ray spectroscopy analysis revealed the delamination sites and crack propagation in ReBCO tapes. Further investigation revealed that the DS of ReBCO tapes is primarily determined by two factors: (1) interfacial bonding strength between the layers and (2) mechanical strength of the stabilizer edges. These findings point to a potential approach to enhancing DS by increasing the strength of the stabilizer (e.g. by optimizing the Cu electroplating process).

Cu stabilizer↗

DOE OSTI2026

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗

DOE OSTI2025

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES↗

DOE OSTI2025

Quench Analysis of EIC Central Detector (ePIC) Solenoid Magnet (MARCO)

The 2.0 T superconducting solenoid magnet (MARCO) being designed to be built for the EIC central detector (ePIC). Here, the magnet consists of three set of coil in a vacuum impregnated coil pack. Each coil pack consists of a total 554 turns winding wrapped with 6 layers of winding per coil. The coil case assembly is cooled by supercritical helium at 4.6 K. This manuscript presents quench analysis carried out on the magnet system along with the conductor stability. For the normal operation the magnet at 4.6 K and at normal operating current, the quench evaluation is carried out at 4.7 K for design assessment. The quench analysis carried out using commercial multi-physics FEA tool along an equivalent single coil quench model and the requirements for the magnet protection is also presented.

Ghoshal, Probir K. [Thomas Jefferson National Acce↗

DOE OSTI2026

Enhanced electrochemical performance and extended cycling of resorcinol-formaldehyde derived N-doped carbon xerogel for alkali metal-ion (Li/Na/K) batteries

Resorcinol formaldehyde-derived carbon xerogel (RFC) is a versatile material with tuneable properties, synthesized through a simple sol-gel method. This study presents nitrogen-doped RF carbon xerogel (N-RFC) with 11.8 at% nitrogen doping, offering a microporous architecture ideal for alkali metal-ion (Li, Na, K) batteries. The porous N-doped framework enhances electrochemical performance by improving ion transport, increasing active storage sites, and significantly boosting metal-ion adsorption, particularly through pyrrolic nitrogen, as revealed by first-principles calculations supported by XPS analysis. N-RFC anodes showed excellent cycling stability, high-capacity retention, and fast charge/discharge capabilities, rendering them suitable for commercial applications. Notably, the N-RFC anode demonstrates high-rate long-term cycling stability, retaining its capacity of 83.5 % (188 mAh/g at 2 C-rate) and 50 % (133 mAh/g at 1250 mA/g) over 1000 cycles for Li and Na-ion batteries, respectively, favorable for commercial battery applications. Additionally, N-RFC demonstrates a reversible capacity of 120 mAh/g after 394 cycles with a retention of 82 % for K-ion batteries. The ability of the material to accommodate larger ions like Na + and K + further emphasizes its versatility and potential application in diverse alkali metal-ion battery systems.

25 ENERGY STORAGE↗

DOE OSTI2024

Identification of Suitable Vacuum Gas Oils as Plasticizers Using HT-GC × GC-HRMS

Vacuum gas oils (VGOs) have long been heralded as effective plasticizers due to their high boiling-points and lubricating properties. One of these VGOs, HyVac Oil 93050, is a paraffinic plasticizer that is currently utilized in polymer–plasticizer explosive formulations, while this oil effectively increases the elasticity and decreases the sensitivity for current formulations. We seek to identify alternative VGOs with similar density, viscosity, molecular composition, and impurities for future formulations. In this study, 18 VGOs, including HyVac Oil 93050, were initially evaluated for their density and viscosity. The oils were then ranked based on physical characteristics and analyzed for molecular composition using high-temperature comprehensive two-dimensional gas chromatography with high-resolution time-of-flight mass spectrometry (HT-GC × GC-HRMS). Further, the HT-GC × GC-HRMS chromatograms of the top 10 most similar VGOs to HyVac Oil 93050, according to density, were compared utilizing the Pearson correlation coefficient. Three of the oils with densities similar to those of HyVac Oil 93050 had a Pearson correlation within the lot-to-lot variation of HyVac Oil 93050. For the top 3 candidates, Pearson correlation was then utilized as a feature selection technique to discover significant chemical differences. Positive chemical ionization (PCI) and negative chemical ionization (NCI) HT-GC × GC-HRMS chromatograms were also evaluated and aided in the discovery and identification of several key compounds including additives that acted as stabilizers for the VGOs. After this further chemical analysis, two of the original 17 potential VGOs were determined to be physically and chemically similar to HyVac Oil 93050.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗

DOE OSTI2025

Lattice Structure and Dynamics of Sparse Molecular Crystals: OsO 4 and RuO 4

OsO 4 and RuO 4 are molecular oxides with unique tetrameric structures and rare +8 oxidation states. Accurately modeling their properties remains challenging for density functional theory (DFT) due to weak intertetramer interactions, which standard functionals fail to capture. Here, in this work, we show that the van der Waals (vdW)-corrected density functional (vdW-DF-optB86b) provides structural parameters that are much closer to experimental values than the standard generalized gradient approximation, with volume predictions that fall within the experimentally observed range. Phonon band structure analysis shows that the inclusion of vdW interactions stabilizes soft phonon modes, highlighting the importance of dispersion corrections for accurate predictions of lattice dynamics. Experimental measurements of the phonon density of states for OsO 4 , obtained via inelastic neutron scattering, demonstrate good agreement with our vdW-DF-optB86b calculations. These results validate OsO 4 and RuO 4 as valuable benchmarks for structural and vibrational calculations via vdW-corrected DFT methods and offer insights for studying the broader class of sparse molecular materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Single-Phase Spinel NiCo 2 O 4 as Highly Active and Stable Electrocatalysts for Urea Oxidation Reaction in Urea Electrolysis

Exploring and designing a stable and active catalyst for the urea electro-oxidation reaction (UOR, CO(NH 2 ) 2 + 6OH – → CO 2 + N 2 + 5H 2 O + 6e – ) is crucial for the long-term sustainability of ecological systems and clean energy production. We found that spinel NiCo 2 O 4 is a stable and active electrocatalyst for UOR at a relatively low anodic potential without triggering the competing oxygen evolution reaction (OER). A urea electrolysis cell (CO(NH 2 )2 + H 2 O → CO 2 + N 2 + 2H 2 ) utilizing a spinel NiCo 2 O 4 anode and a commercial Pt cathode was further characterized through galvanostatic polarization tests, demonstrating excellent structural stability at various current densities. Post-mortem analysis of long-term urea electrolysis measurements suggested that NiCo 2 O 4 electrocatalysts maintained a stable spinel structure. However, redistribution of Ni 3+ to Ni 2+ valence on the catalyst surface was observed, in contrast to the intact Co valence, indicating that (i) Ni sites are active toward urea adsorption and sequential electro-oxidation; (ii) while urea oxidation proceeds primarily through the direct electro-oxidation mechanism, chemical reactions between the Ni 3+ site and urea occur during long-term electrochemical UOR operation. Density functional theory (DFT) simulations were used to calculate the adsorption energies of urea molecules on NiO, Co 3 O 4 , and NiCo 2 O 4 , revealing the importance of regulating the configuration of adsorbed urea molecules on the NiCo 2 O 4 surface.

36 MATERIALS SCIENCE↗