Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS CALCULATION”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Efficient calculation of self magnetic field, self-force, and self-inductance for electromagnetic coils with rectangular cross-section

Abstract For designing high-field electromagnets, the Lorentz force on coils needs to be computed in order to design suitable support structures, and the inductance should be computed to evaluate the stored energy and dynamics. Also, the magnetic field and its variation inside the conductor is of interest for computing stress and strain, and due to superconducting quench limits. For these force, inductance, energy, and internal field calculations, the coils cannot be naively approximated as infinitesimally thin filaments due to divergences when the source and evaluation points coincide, so more computationally demanding calculations are usually required, resolving the finite cross-section of the conductors. Here, we present a new alternative method that enables the internal magnetic field vector, self-force, and self-inductance to be computed rapidly and accurately within a 1D filament model. The method is applicable to coils for which the curve center-line can have general noncircular shape, as long as the conductor width is small compared to the radius of curvature. This paper extends a previous calculation for circular-cross-section conductors (Hurwitz et al 2024 IEEE Trans. Magn. ) to consider the case of rectangular cross-section. The reduced model is derived by rigorous analysis of the singularity, regularizing the filament integrals such that they match the true high-dimensional integrals at high coil aspect ratio. The new filament model exactly recovers analytic results for a circular coil, and is shown to accurately reproduce full finite-cross-section calculations for a non-planar coil of a stellarator magnetic fusion device. Due to the efficiency of the model here, it is well suited for use inside design optimization.

Landreman, Matt (ORCID:000000027233577X)↗
DOE OSTI2025

Evaluating methods to reduce duration of near-threshold fatigue crack growth rate measurements for low-alloy steels in hydrogen gas

Measurement of the near-threshold fatigue crack growth rate (da/dN) vs. stress-intensity factor range (∆K) relationship in hydrogen gas is essential for maximizing the calculated design fatigue life of high-pressure hydrogen storage vessels. However, such measurements are rarely performed, since the low cyclic loading frequencies applied in standard practice lead to prohibitively protracted test durations. The objective of this study was to demonstrate two means for reducing test durations when measuring near-threshold da/dN vs. ΔK relationships under decreasing ΔK for low-alloy pressure vessel steels in hydrogen gas: 1) imposing steeper K-gradients relative to the recommended limits in standards such as ASTM E647, and 2) increasing cyclic loading frequency relative to typical values applied during fatigue crack growth testing of low-alloy steels in hydrogen gas. Recognizing that steeper K-gradients could amplify loading-history effects, test methods employing this approach were designed to mitigate such effects by either maintaining constant K max or gradually increasing the K-gradient as the threshold was approached. Although the varying K-gradient method was vulnerable to loading-history effects in the form of plasticity-induced crack closure, particularly at lower stress ratio (R) and higher starting K max values, these effects could be compensated by applying the adjusted compliance ratio (ACR) method. Here, it was demonstrated that steeper K-gradients in concert with increased cyclic loading frequency reduced the duration of near-threshold fatigue crack growth tests in hydrogen gas by more than 99% relative to standard practices.

Fatigue threshold↗
DOE OSTI2025

Drought adaptation index (DAI) based on BLUP as a selection approach for drought-resilient switchgrass germplasm

This study introduces a Drought Adaptation Index (DAI), derived from Best Linear Unbiased Prediction (BLUP), as a method to assess drought resilience in switchgrass (Panicum virgatum L.). A panel of 404 genotypes was evaluated under drought-stressed (CV) and well-watered (UC) conditions over four consecutive years (2019–2022). BLUP-estimated biomass yields were used to calculate the DAI, which enabled classification of genotypes into four adaptation groups: very well-adapted, well-adapted, adapted, and unadapted. The DAI was compared with conventional drought tolerance indices, including the Stress Susceptibility Index (SSI), Stress Tolerance Index (STI), Geometric Mean Productivity (GMP), and Yield Stability Index (YSI). Correlation analyses demonstrated strong agreement between DAI and these indices, supporting its validity and consistency. Biplot analyses using the Genotype plus Genotype-by-Environment Interaction (GGE) and Additive Main Effects and Multiplicative Interaction (AMMI) models revealed significant genotype-by-environment interactions (GEI) and identified J222.A, J463.A, and J295.A. A as high-performing genotypes, with J222.A exhibiting greater yield stability across treatments and years. Additionally, DAI isoline curves provided a graphical representation of differential genotype performance under drought and control conditions. These visualizations aided in distinguishing genotypes with stable and superior biomass yield across contrasting environments. Overall, the BLUP-based DAI is a robust and practical selection tool that improves the accuracy of identifying drought-resilient, high-yielding switchgrass genotypes. Its integration into breeding programs offers a comprehensive framework for improving biomass productivity and stress adaptation under variable climatic conditions. The application of DAI supports the development of climate-resilient cultivars and contributes to sustainable bioenergy and forage production systems.

BLUP↗
DOE OSTI2025

Chromium enhances the mechanical performance of 3d transition metal high entropy alloys

A comprehensive examination of the compositional effects on the deformation behavior of high entropy alloys (HEAs) was conducted through nanoindentation, indentation creep, and stress relaxation experiments at ambient temperature. The alloys investigated included NiCoFe, NiCoCr, NiCoCrFe, and NiCoCrFeMn, all of which underwent identical thermomechanical processing. From the experimental results, NiCoCr exhibited the highest maximum shear stress for dislocation nucleation (13 GPa) and nanoindentation hardness (3.8 GPa) among the alloys tested. It also showed the lowest stress sensitivity (0.02) and activation volume (3–4 b 3 ) under creep and stress relaxation conditions, indicative of its superior plastic flow properties. Density functional theory (DFT) calculations further revealed an uneven charge density and larger bonding-length variation in NiCoCr due to the presence of Cr, which increased lattice distortion, impeding dislocation movement. Furthermore, in this study, post-deformation microscopy using transmission electron microscopy (TEM) revealed a high density of stacking faults and twins in NiCoCr that effectively enhanced its creep strength. The results highlight the significance of specific elemental effects, particularly from Cr in this study, over configurational entropy effect (i.e., compositional complexity) in governing the deformation microstructure and mechanical properties of HEAs. Finally, these insights will be instrumental for the design and optimization of advanced alloys for load-bearing applications.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Probing the elastic modulus and hardness of superionic boron cluster solid electrolytes

Polyhedral borane salts are highly tunable compounds with the potential to be used as solid electrolytes due to their high Li and Na superionic conductivity and relatively wide electrochemical stability window. In considering their application to all solid-state batteries, their mechanical properties play a critical role in their operation due to stresses the solid electrolyte is subjected to during battery operation. Density functional theory (DFT) calculations have provided an initial assessment of bulk and elastic moduli for some selected boron cluster solid electrolytes, but direct measurements are still scarce. Here, in this paper, we report the elastic moduli and hardnesses of LiCB 11 H 12 , LiCB 9 H 10 , NaCB 11 H 12 , and NaCB 9 H 10 for the first time using nanoindentation continuous stiffness method (CSM) measurements under inert atmosphere. Experimental modulus values ranging 8.8–12.3 GPa and hardness values ranging 0.15–0.38 GPa are lower than that of other inorganic solid-state electrolytes such as oxide- and sulfide-based solid electrolytes. DFT calculations on the expected moduli of these compounds are also presented and discussed. Analysis of the indentation plasticity index reveals that these compounds have higher plasticity index compared to other common oxide and sulfide solid electrolytes. According to the calculated Pugh's ratio, all compounds in this study except LiCB 11 H 12 are considered ductile.

25 ENERGY STORAGE↗
DOE OSTI2025

Tuning the structure and superconductivity of SrNi 2⁢ P 2 by Rh substitution

The compound SrNi 2⁢ P 2 is unique among the ThCr 2 ⁢Si 2 class since it exhibits a temperature-induced transition upon cooling from an uncollapsed tetragonal (ucT) state to a one-third-collapsed orthorhombic (tcO) state where one out of every three P-rows bond across the Sr layers. This compound is also known for exhibiting bulk superconductivity below 1.4 K at ambient pressure. Here, in this paper, we report on the effects of Rh substitution in Sr⁢(Ni 1-x ⁢Rh x ) 2 ⁢P 2 on the structural and superconducting properties. We studied the variation of the nearest P-P distances as a function of the Rh fraction at room temperature, as well as its temperature dependence for selected compositions. We find that increasing the Rh fraction leads to a decrease in the transition temperature between the ucT and tcO states, until a full suppression of the tcO state for x ≥ 0.166. The superconducting transition first remains nearly insensitive to the Rh fraction, and then it increases to 2.3 K after the tcO state is fully suppressed. These results are summarized in a phase diagram, built upon the characterization by energy dispersive x-ray spectroscopy, x-ray diffraction, resistance, magnetization, and specific heat measurements done on crystalline samples with varying Rh content. The relationship between band structure, crystal structure, and superconductivity is discussed based on previously reported band structure calculations on SrRh 2 ⁢P 2 . Moreover, the effect of Rh fraction on the stress-induced structural transitions is also addressed by means of strain-stress studies done by uniaxial compression of single-crystalline micropillars of Sr⁢(Ni 1-x ⁢Rh x ) 2 ⁢P 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Deconvoluting thermomechanical effects in X-ray diffraction data using machine learning

X-ray diffraction is ideal for probing the sub-surface state during complex or rapid thermomechanical loading of crystalline materials. However, challenges arise as the size of diffraction volumes increases due to spatial broadening and because of the inability to deconvolute the effects of different lattice deformation mechanisms. Here, we present a novel approach that uses combinations of physics-based modeling and machine learning to deconvolve thermal and mechanical elastic strains for diffraction data analysis. The method builds on a previous effort to extract thermal strain distribution information from diffraction data. The new approach is applied to extract the evolution of the thermomechanical state during laser melting of an Inconel 625 wall specimen which produces significant residual stress upon cooling. A combination of heat transfer and fluid flow, elasto-plasticity and X-ray diffraction simulations is used to generate training data for machine-learning (Gaussian process regression, GPR) models that map diffracted intensity distributions to underlying thermomechanical strain fields. First-principles density functional theory is used to determine accurate temperature-dependent thermal expansion and elastic stiffness used for elasto-plasticity modeling. The trained GPR models are found to be capable of deconvoluting the effects of thermal and mechanical strains, in addition to providing information about underlying strain distributions, even from complex diffraction patterns with irregularly shaped peaks.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evaluation of Maximum Allowable Working Pressure and Svensson Burst Pressure Recommended in API 579-1 2021 Edition

Abstract API 579-1/ASME FFS-1 2021 Edition provides the minimum wall thickness, the maximum allowable working pressure (MAWP), and the membrane stress equations for thin and thick-walled cylindrical shells subject to internal pressure in Section 2C.3.3.1 of Appendix 2C – Thickness, MAWP, and Stress Equations for an FFS Assessment. The minimum wall thickness and MAWP are determined using the hoop stress and the Tresca yield criterion. Section 2C.7 – Estimation of Burst Pressure newly added the Svensson method for calculating burst pressure of cylindrical shells under internal pressure, where the plastic yielding is characterized by the von Mises yield criterion. For thin-walled cylinders, the von Mises flow solution of burst pressure in Equation (2C.179) was recommended. For thick-walled cylinders, an implicit burst pressure solution in an integral equation (2C.176) was recommended. But this integral equation is inconvenient to use in practice. It is well known that the classic plasticity theory includes the Tresca and von Mises yield criteria, with the Tresca criterion predicting a lower bound solution and the von Mises criteria predicting an upper bound solution. In addition, the present author developed an average shear stress yield criterion that can determine more accurate limit and burst pressures for thin and thick-walled cylinders. This work uses these three yield criteria to evaluate the minimum required wall thickness, MAWP and Svensson burst pressure recommended in the API 579 code.

burst pressure↗
DOE OSTI2025

Evaluation of Maximum Allowable Working Pressure and Svensson Burst Pressure Recommended in API 579-1 2021 Edition

ABSTRACT API 579-1/ASME FFS-1 2021 Edition provides the minimum wall thickness, the maximum allowable working pressure (MAWP), and the membrane stress equations for thin and thick-walled cylindrical shells subject to internal pressure in Section 2C.3.3.1 of Appendix 2C – Thickness, MAWP, and Stress Equations for an FFS Assessment. The minimum wall thickness and MAWP are determined using the hoop stress and the Tresca yield criterion. Section 2C.7 – Estimation of Burst Pressure newly added the Svensson method for calculating burst pressure of cylindrical shells under internal pressure, where the plastic yielding is characterized by the von Mises yield criterion. For thin-walled cylinders, the von Mises flow solution of burst pressure in Equation (2C.179) was recommended. For thick-walled cylinders, an implicit burst pressure solution in an integral equation (2C.176) was recommended. But this integral equation is inconvenient to use in practice. It is well known that the classic plasticity theory includes the Tresca and von Mises yield criteria, with the Tresca criterion predicting a lower bound solution and the von Mises criteria predicting an upper bound solution. In addition, the present author developed an average shear stress yield criterion that can determine more accurate limit and burst pressures for thin and thick-walled cylinders. This work uses these three yield criteria to evaluate the minimum required wall thickness, MAWP and Svensson burst pressure recommended in the API 579 code.

burst pressure↗
DOE OSTI2025

Planar Defect Layers Template a High-Pressure InBi Polymorph

The short- and long-range order of III–V materials under high pressure has long been the subject of debate, with advancements in structural characterization leading to significant revisions to the accepted structural models. Despite these revisions, previous high-pressure structural assignments in the In–Bi system include the site-disordered β-Sn structure type, a structure type demonstrated to be nonexistent in analogous III–V systems. While X-ray diffraction is consistent with site disordering in InBi at high pressure, cluster expansion calculations indicate that disordering requires temperatures above 3000 K. Here, we propose InBi as a model material for studying unique high-pressure planar defects due to its highly anisotropic stress-dependent properties and structure. Specifically, we identify two sets of planar defects that mimic the diffraction pattern of a site disordered β-Sn structure type and are compatible with the calculated disorder barrier. We derive these defects by symmetry relations over crystallographic transitions. Density functional theory calculations of the proposed defects suggest that these defects are stabilized by diminishing interlayer separations with pressure. Further, we find that one of the proposed defects closely resembles a bulk high-pressure phase of InBi, InBi-ϵ, and we assert that the proposed defects order upon heating, acting as a template for InBi-ϵ growth. The proposed defects and their electronic structure provide a basis for the trend of superconducting critical temperature with increasing pressure. These methods for identifying defects are generalizable to other materials with reports of site disorder at high pressure, prompting a broader search for related high-pressure defects.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Neural network potentials with effective charge separation for non-equilibrium dynamics of ionic solids: a ZnO case study

Developing neural network potentials (NNPs) accurate under non-equilibrium dynamics is challenging, as such systems require extensive sampling beyond equilibrium phases. Here we construct high-fidelity NNPs for zinc oxide (ZnO), a polymorphic ionic solid, using density functional theory (DFT) reference data. To efficiently capture transitional configurations, we combine enhanced-sampling molecular dynamics with empirical potentials, data distillation, and pretraining on short-range atomic energies (A-Train), followed by transfer learning with DFT-relabeled datasets. This hierarchical approach improves transferability across polymorphs and stress states. We further introduce effective charge separation, treating long-range Coulombic terms analytically while short-range residual interactions are learned by the NNP. The optimal effective charges fall in the range 0.5–1.0 q e , consistent with dielectric-screened values derived from formal charges but distinct from Bader estimates. Motivated by this observation, we propose a simple data-driven protocol in which effective charges are optimized by comparing DFT reference energies with explicit Coulomb calculations, without additional NNP training. This strategy improves accuracy and transferability in DFT-level predictions of energies, forces, and stress. Together, these results provide a practical charge-selection framework for robust NNP development in ionic solids, enabling reliable simulation of polymorphic phase transformations and non-equilibrium dynamics.

Chemistry↗
DOE OSTI2025

Anisotropic structure and optical engineering of strontium titanate and zirconia responding to sequential hydrogen and helium irradiation

Unraveling the correlation between strain engineering with anisotropic optical properties via nanoscale defects gradually evolved into a strategy for fundamental studies and technological applications, yet it remains understudied in functional complex oxides compared to metals and semiconductors. Here, the methodology of strain engineering for the determination of the lattice parameters, parallel/normal to the sample surface, in the individual layers of single-crystalline superlattices is derived, which is based on analysis of high-angle X-ray diffraction measurements in combination with diffraction reciprocal space mapping. With modeling that takes into account the effect of elastic properties, two elastically anisotropic materials, SrTiO 3 and ZrO 2 , have been compared in terms of defect-induced elastic strain caused by individual and sequential H + and He 2+ irradiation. The anisotropic lattice swelling with corresponding refractive index, and obstructive behavior of elastic strain recovery are demonstrated in SrTiO 3 , while approaching strain behaviors accompanied by isotropic refractive index distribution in both orientations are confirmed in ZrO 2 . Under sequential He 2+ /H + irradiation, pre-existing He-vacancy complexes acted as vacancy traps, preferentially capturing H + clusters to induce lattice distortion and enhance absorption. Nanohardness increments (ΔH) calculated by the DBH model matched nanoindentation results, confirming that sequential irradiation generated higher indentation yield stress than individual irradiations.

Anisotropic expansion↗
DOE OSTI2025

A Combined Computational and Mathematical Analysis of Interconnect Fatigue Potential in Photovoltaic Modules

A finite element model of a 60-cell monocrystalline silicon glass-polymer photovoltaic module was simulated with ±1.0 kPa and ±2.4 kPa loads applied to the glass to calculate the deformation under load. Cell-to-cell displacements were used to approximate interconnect strain and stress. A mathematical fatigue cycle life relation was fitted to data for the interconnect material (copper), to generate a life prediction at each interconnect location based on the local stress means, reversal extents, and amplitudes. Interconnect stress was found to be significantly asymmetric about zero despite symmetric positive and negative module loads due to laminate thickness offsets about the neutral plane and the effects of module framing. Cycle life results indicated that interconnect fatigue failure was unlikely to occur over a 30-year lifetime of conservative wind and snow load cycles since the typical cell design feature of leaving some unconstrained length between the cell edge and first solder pad increases the effective gauge length and decreases the stress levels below the material endurance limit. Follow-up analyses found that 3.6 mm and 6.4 mm were the minimum unconstrained lengths required to survive the assumed lifetime of wind and snow cycles, respectively, confirming that typical industrial module constructions with 8–15 mm unconstrained lengths should survive conservatively. Notably, large magnitude, low-cycle snow loading was consistently the limiting factor requiring a longer unconstrained interconnect length. Finally, insights and workflows from this study inform module interconnection design limits for survival against mechanical fatigue in deployment environments.

14 SOLAR ENERGY↗
DOE OSTI2025

Understanding the deformation behavior of the γ rich transformative Fe 38.5 Mn 20 Co 20 Cr 15 Si 5 Cu 1.5 complex concentrated alloy using in situ synchrotron diffraction

In situ tensile testing coupled with synchrotron x-ray diffraction was used to study the deformation behavior of metastability-engineered Fe 38.5 Mn 20 Co 20 Cr 15 Si 5 Cu 1.5 complex concentrated alloy. Monitoring the evolution of phase fraction and strain hardening response allowed the determination of true critical stress for the onset of the transformation-induced plasticity (TRIP) to be ∼375 MPa, preceded by slip starting at ∼255 MPa. In situ EBSD was used to validate the critical stress for transformation at the microstructure level and observe slip traces to confirm prior slip activity before the transformation. Further, a modeling framework based on stacking fault energy (SFE) was developed to predict the critical stress for transformation. Modeling suggested the SFE of the alloy to fall nearly 15 mJ/m 2 , which agrees well with SFE values calculated using synchrotron peak shifting (12 mJ/m 2 ) and thermodynamic calculation (11 mJ/m 2 ). As a result of γ-fcc to ε-hcp phase transformation, new {0002} ε planes emerged parallel to unaligned {111} γ planes with the tensile loading following S-N orientation relationship. Such selective emergence of new diffraction rings corresponding to ε phase is understood based on the reorientation of γ crystals with reference to tensile axis. In conclusion, this approach can be extended to effectively design alloys based on critical stress required for activating different deformation mechanisms to further push the limits of the strength-ductility envelope.

Complex concentrated alloy↗
DOE OSTI2026

Quaternary i-MAX Phases (Mo 2/3 RE 1/3 ) 2 AlC (RE: Dy, Tb, Er): Experimental Characterization and First-Principles Insights into their Fundamental Properties

Rare earth (RE)-based materials have unique electronic, magnetic, and optical properties, leading to the recent discovery of atomically layered solids with the chemical formula (M' 2/3 RE 1/3 ) 2 AlC, which have since garnered significant attention in the scientific community. This study aims to synthesize, characterize, and investigate the structural and thermal stability of the RE i-MAX phases. We prepared i-MAX phases using molybdenum (Mo) as M′ and RE elements as Dy, Tb, and Er, namely (Mo 2/3 Dy 1/3 ) 2 AlC, (Mo 2/3 Tb 1/3 ) 2 AlC, and (Mo 2/3 Er 1/3 ) 2 AlC. Structural characterization through x-ray diffraction (XRD) and Raman spectroscopy confirms the formation of the RE-based i-MAX phase, along with the presence of minor impurity phases in the alloys. Thermogravimetric analysis (TGA) conducted up to 1000°C under ambient conditions reveals that the i-MAX phases remain thermally stable up to approximately 450°C, beyond which oxidation leads to a noticeable weight gain in all samples. Differential scanning calorimetry (DSC) measurements during heating and cooling cycles show endothermic and exothermic peaks for (Mo 2/3 Dy 1/3 ) 2 AlC i-MAX in the 410–420°C range, indicating a temperature-induced minor atomic arrangement. In contrast, these peaks are absent in the Tb- and Er-based i-MAX phases. These findings offer valuable insights into the thermal behavior and stability of these i-MAX phases under thermal stress, contributing to a deeper understanding of their unique properties. Furthermore, first-principles density functional theory (DFT) calculations were performed to investigate the electronic and optical properties of the i-MAX phases. The results reveal their metallic nature, with pronounced contributions from Mo and RE elements near the Fermi level and within the conduction band.

Rare earth↗
DOE OSTI2024

Multi-stage load partitioning in additively manufactured Al6061+TiC nanocomposite characterized by in-situ neutron diffraction

Metal matrix composites (MMCs), combining metal matrix and ceramic particles, exhibit high mechanical strength and stiffness compared to conventional metallic materials. During fusion-based additive manufacturing (AM), MMCs undergo a rapid heating and cooling thermal history, which affects the interfacial bonding, thermal misfit stress and mechanical load transfer behavior between the metal matrix and particles, thus impacting the bulk mechanical performance. Here, we investigate the deformation dynamics and phase-specific load transfer in fusion-based AMed Al6061 +TiC MMCs through in-situ neutron diffraction. By calculating the phase-specific lattice strains, a multi-stage load transfer and deformation behavior during uniaxial compression is revealed: (1) elastic deformation in Stage I for both phases; (2) sudden stress rebalance between two phases in Stage II, evidenced by a sudden decrease in Al lattice strain and increase in TiC; (3) active load carrying by both phases in Stage III, with both phases experiencing an increase in lattice strains. Further, the deformation mechanism of each stage is deducted by correlating the evolutions of the interphase stresses and peaks’ broadening and intensity of the Al matrix. The local plastic deformation of the Al matrix near the phase interface is triggered and leads to stress rebalancing in Stage II. The global plastic deformation subsequently propagates throughout the Al matrix in Stage III, and the composite is further hardened through both dislocation multiplication and load transfer. The findings offer valuable insights into deformation and load-sharing behavior in AMed MMCs.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Plasma-Assisted Pre-Chamber Ignition System for Highly Dilute Stoichiometric Heavy-Duty Natural Gas Engines (Final Technical Report)

This project explored advanced ignition technologies to significantly enhance efficiency and reduce operating costs for heavy-duty natural gas engines operating at stoichiometric conditions, while meeting ultra-low NOx emission standards. The main goal was to develop and validate a plasma-assisted pre-chamber ignition system that could deliver at least a 2% increase in brake thermal efficiency (BTE) and a 4% decrease in total cost of ownership (TCO) compared to a typical multi-cylinder engine with three-way catalyst aftertreatment, ensuring compatibility with the expected 2027 EPA/CARB regulations. In the first half of the project, the research team concentrated on developing and testing plasma-assisted pre-chamber ignition using nanosecond pulsed discharges. Extensive experiments were conducted in an optically accessible rapid-compression and expansion machine, a constant-volume chamber, and an optical single-cylinder engine. Experiments were coupled with CFD simulations. The work produced unique insights into pre-chamber flame formation, jet ignition, dilution effects, and flame quenching at pressures, temperatures, and dilution levels relevant to engines. Although plasma-assisted ignition showed promise in controlled lab settings, the research also identified fundamental and practical challenges when applying this technology to real engine conditions. Midway through the project, a crucial pivot was made, guided by three key findings. First, the power electronics required for nanosecond plasma discharges were found to be too costly for commercial use, undermining the project’s cost-of-ownership goals. Second, nanosecond plasma ignition was highly sensitive to turbulent flow in the pre-chamber, resulting in lower ignition reliability than traditional spark under engine-like conditions. Third, achieving a truly diffuse low-temperature plasma at high pressures near top dead center was not possible, reducing the anticipated chemical enhancement benefits. These results collectively suggested that continuing with plasma-assisted ignition was unlikely to meet both efficiency and cost objectives. In response, the project shifted focus to a more realistic approach: enhancing traditional spark-based pre-chamber ignition with significantly less spark energy. Using insights gained earlier in the project, the team redesigned the pre-chamber to maintain high dilution tolerance and quick combustion, even with lower ignition energy. Testing confirmed that with optimized pre-chamber design and combustion timing, a lower-energy spark could reliably ignite highly diluted stoichiometric mixtures, reduce burn time, and boost thermal efficiency. Final engine testing and techno-economic analysis verified that this revised approach successfully achieved the project goals. The optimized pre-chamber ignition system provided over a 2% increase in calculated brake thermal efficiency compared to the baseline engine. Notably, the lower ignition energy and simplified hardware reduced component stress, extended maintenance intervals, and lowered the total cost of ownership. When used with stoichiometric operation and traditional three-way aftertreatment, the system remained compatible with near-zero NOx emissions targets without increasing cost or complexity in the emissions control system. In summary, although the project deviated from its initial plasma-assisted ignition idea, the work produced a more practical and commercially viable solution. The results show that precisely optimized, low-energy pre-chamber spark ignition can significantly improve efficiency and reduce overall ownership costs for heavy-duty natural gas engines. This directly aligns with DOE goals for cleaner, more efficient, and cost-effective transportation technologies.

03 NATURAL GAS↗
DOE OSTI2025

Evaluating design safety margins in the American Society of Mechanical Engineers graphite core components design-by-analysis assessments

Graphite is an important material being used for core components in next-generation high-temperature gas-cooled nuclear reactors. The selection of graphite grade for a specific Designer is a complex task, dependent on reactor conditions, component functionality, and required reliability. The American Society of Mechanical Engineers (ASME) provides two semi-probabilistic design-by-analysis assessments to evaluate graphite core components against design reliability targets. The simplified assessment uses a 2-parameter Weibull distribution to describe the graphite grade’s tensile-strength distribution to establish component stress limits. The full assessment uses the 3-parameter Weibull distribution and a modified Weakest-Link Theory approach to calculate a component design probability of failure. The paper defines recommended assessment rules, which are the as-written simplified assessment and the full assessment with parameter lower bounds, the modulus update with threshold reduction, and the 2027 grouping rules. Code rules are applied to three grades: 2114, IG-110, and NBG-18. The baseline margin calculation is developed using the experimental tensile dogbone specimen. Percent margin is defined as the percent reduction in the median experimental load to obtain the allowable load per ASME assessments. Under the recommended rules, the SRC–1 margin in the simplified assessment ranged from 40.2 % to 52.7 % among the grades in this study and from 36.1 % to 49.8 % in the full assessment. The full assessment only decreases the margins by 2.5–4.5 % for the SRC-1 components and 0–1.5 % for the SRC-2 components for this baseline case. Margin is inversely related to material median strength (i.e., the strongest grade, 2114, has the lowest margin).

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗