Engineering Papers⌕ Search

SEARCH · Engineering Papers

Search Engineering Papers

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

A framework and tool for designing cost-effective, resilient, and circular net-zero supply chains under uncertainty with an application to multilayer plastic films

While 55% of Fortune 500 companies have committed to achieving net-zero emissions and/or zero-waste operations by 2035, only 2% are currently on track, revealing a critical gap between ambition and action. Designing supply chains that reduce both emissions and waste is a complex non-intuitive, multi-objective challenge, compounded by the high costs of new technologies and the need for resilient, profitable solutions. This paper aims to address this challenge by presenting a generic framework and multi-objective optimization formulation for designing cost-effective, circular, and resilient supply chains under uncertainty, implemented through a user-friendly decision-support tool with intuitive data visualization capabilities, enabling communication of results to both technical and non-technical stakeholders. We demonstrate the application of this framework in the context of multilayer plastic films (barrier films), which are widely used in food packaging and composite materials. The model quantifies trade-offs across three objectives: minimizing global warming potential, maximizing circularity, and minimizing cost. A key contribution of this work is the explicit modeling of technological resilience, the ability of supply chains to maintain function under disruption. In the cost-minimization case, the resilience constraint makes the design approximately three times more expensive in the short-term metric, but shifts the system from relying on a single recovery pathway to a portfolio of four recovery pathways, improving the robustness of the optimization solution under uncertainty. Lastly, we introduce TranZero, a decision-support tool that integrates material flow analysis, hotspot identification, and optimization-based scenario planning to support net-zero and circularity decisions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

PET waste- and bio-derived imine vitrimers for shape-memory, intrinsic flame-retardant, and recyclable carbon fiber composites

Developing circular multifunctional vitrimers and carbon fiber–reinforced polymers (CFRPs) that are simultaneously recyclable, mechanically robust, and intrinsically flame retardant remains a major challenge. Here, in this study, we report multifunctional vitrimers and their carbon fiber–reinforced vitrimer (CFRV) composites, where the vitrimer design integrates closed-loop recyclability, enhanced interfacial adhesion, and intrinsic flame retardancy within a single materials platform. The vitrimer matrix is synthesized from post-consumer polyethylene terephthalate (PET) waste and a vanillin-derived phosphorus-containing crosslinker, forming an imine-based network. The resulting vitrimer resin exhibits high tensile strength, thermal healability, repeated reprocessability, programmable shape memory, and rapid chemical depolymerization under mild conditions. Amine-functionalized carbon fibers significantly improve fiber–matrix interfacial bonding, yielding CFRVs with tensile strengths up to 789 MPa and complete recovery of structurally intact fibers after chemical recycling. The phosphorus-rich aromatic network further imparts intrinsic flame retardancy, enabling self-extinguishing behavior without external additives. This work advances a materials design paradigm for next-generation multifunctional, sustainable vitrimers and CFRVs, while simultaneously addressing the recycling challenges associated with both plastic and CFRP waste.

Bio-derived crosslinker↗
DOE OSTI2026

An immersed interface method for microstructure-scale electrochemical battery models: numerical formulation and performance portable implementation

We present the numerical formulation, verification, and performance portable implementation of an immersed interface method for microstructure scale electrochemical modeling of batteries. The innovation in this approach is the resolution of chemical species and electrostatic potential discontinuities at active interfaces without the use of interface conforming unstructured grids. A unified formulation on Cartesian grids for all domains (electrodes and electrolyte) is used with interfacial flux conditions applied using volume fraction or “color” function gradients. We have developed one dimensional and two dimensional test cases with analytic solutions for electrochemical modeling using which we verified the consistency and accuracy of our scheme. Our solver is also validated against solutions from a macroscale model and an unstructured multi-subdomain solver for a full lithium ion cell. We then demonstrated the utility of our solver on an image-based complex battery electrode microstructure at high charging rate. Our technique also exhibits good scalability on distributed memory architectures using central processing units (CPU), with problem sizes up to 1.8 billion degrees of freedom and with 5400 ranks. Initial performance studies of our open-source performance portable solver showed about 70 times speed up using a graphics processing unit (GPU) compared to single compute core for a problem with 4 million cells.

25 ENERGY STORAGE↗
DOE OSTI2026

Fast solvers for tokamak fluid models with PETSc

Multigrid (MG) is widely recognized as a highly effective solver for the model problem, the Laplacian, but textbook MG fails on most problems of interest. MG methods have been applied to complex, real-world applications with careful consideration of the physical model and discretization. In this work we develop the first step in applying MG methods to science and engineering relevant magnetohydrodynamics (MHD) tokamak models in the M3D-C1 (https://m3dc1.pppl.gov) fusion energy science code. The semi-implicit time integrator in M3D-C1 is composed of many linear solves. The implicit advance of the momentum equation is the most challenging and is the focus of this work. The current production solver in M3D-C1 is a block Jacobi (BJ) preconditioner within a Krylov solver, where blocks group degrees of freedom on planes of constant toroidal coordinate. BJ convergence degrades as the number of planes increases due to the spectral properties of the matrix preconditioned with BJ. The partially magnetic field-aligned, regular toroidal grid structure in M3D-C1 is amenable to semi-coarsening geometric MG in the toroidal direction. This paper develops such a solver and demonstrates competitive performance on a runaway electron model of a SPARC (https://cfs.energy/technology/sparc) disruption, and superior robustness on a stellarator model on which the BJ solver fails to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Magnus method for electronic structure calculations at extreme conditions

We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

general physics↗
DOE OSTI2026

Hydropower potential derived from streamflow extremes for Alaska, USA

Alaska is an expansive region known for its abundant natural resources, including thousands of miles of streams and rivers. These rivers represent potential opportunities for future hydropower development that could provide reliable energy supply for local communities. There is limited long-term high temporal resolution streamflow data available for the region, making data-driven estimates of potential hydropower and its variability across the state challenging. This study provides a novel data-driven approach for hydropower capacity estimation across Alaska. We use supervised machine learning to develop a relationship between the daily and peak flow duration curves in order to augment the size of our dataset from 44 sites to 67 sites. We perform a stochastic hydropower estimation across the 67 sites and identify approximately 1000 MW of total potential hydropower capacity distributed across these sites. Our study provides the first step towards more comprehensive hydropower estimation for this critical region, highlighting the need for future work integrating high-resolution spatial data, community needs, and economic constraints in estimates of potential hydropower development in Alaska.

Hydropower↗
DOE OSTI2026

(CrMnCoNiTi) 3 O 4 high-entropy spinel oxide as a high-performance electrode for supercapacitors

In this study, spinel-structured (CrMnCoNiTi)₃O₄ high-entropy oxides (HEOs) have been successfully synthesized using the solution combustion method, and their performance as supercapacitor electrode materials was investigated. These HEOs exhibit excellent performance, stemming from the unique high-entropy effect and lattice distortion. This structure effectively buffers volumetric strain during cycling, conferring the material with extremely high structural stability and specific capacity. The material demonstrates a high specific capacity of 491.5 C∙g⁻¹ at a current density of 0.5 A∙g⁻¹. Furthermore, it shows extraordinary cycling stability: after 10,000 charge/discharge cycles at 5 A∙g⁻¹, the capacity retention rate is 99%, and the Coulombic efficiency consistently remains close to 100%. Ex-situ structural characterizations further reveal that this excellent cycling durability originates from a cooperative multination valence reconstruction mechanism within an entropy-stabilized spinel framework. Additionally, a symmetric supercapacitor assembled using this material and 1 M KOH electrolyte successfully extended the working voltage to 1.2 V. This research highlights the great potential of high-entropy oxides in synergizing high energy density and ultra-long cycling life, laying a solid foundation for the development of next-generation high-performance supercapacitor electrode materials.

High entropy oxides↗
DOE OSTI2026

Direct regeneration of degraded LiFePO4 cathodes via a separator-enabled prelithiation strategy

A persistent challenge in lithium-ion batteries is the loss of active lithium due to the solid electrolyte interphase (SEI) formation and associated side reactions. While prelithiation employing lithium replenishment separator (LRS) has been proven effective in compensating for lithium loss, previous studies have largely been accompanied by gas evolution or solid residue formation during the prelithiation process. To surmount this challenge, we present a LRS based on 4-fluoro-1,2-dihydroxybenzene lithium salt (LiDF), capable of mitigating lithium loss while producing decomposition products that integrate directly into the electrolyte as functional additives which can assist with the stability of the SEI, free from gas or solid formation, thus establishing a sustainable and environmentally benign strategy for lithium compensation. Incorporation of the LRS enables the pristine LiFePO4||graphite (Gr) full cell to achieve 10.8% higher capacity than the cell with a polypropylene separator (PPS) after 200 cycles at 0.5C. Remarkably, the degraded LiFePO4 (D-LFP)||Gr full cell with the LRS exhibits a 135.8% capacity improvement over the PPS-based cell after 500 cycles. These findings establish the LRS as a powerful approach for both boosting high-performance lithium-ion batteries and recovering the capacity of degraded batteries.

Tao, Fujun↗
DOE OSTI2026

Insights from a coupled thermo-hydro-mechanical analysis of a layered high-temperature thermal energy storage reservoir

Coupled thermal-hydraulic-mechanical (THM) modeling is applied to investigate the performance of a seasonal high-temperature aquifer thermal energy storage operation based on data and conditions from current site investigations at the Geostorage Forsthaus pilot project in Bern (Switzerland). The model includes subhorizontal sand lenses of various lengths and dips that are embedded in a low permeability clay matrix. Thermal energy storage is simulated by seasonal injection and withdrawal of hot (up to 90 °C) water from a main well, with reservoir pressure regulated by two auxiliary wells at a distance of about 70 m from the main well. The results show how targeted injection into deeper permeable storage formations, along with active deep well pressure control, can effectively minimize geomechanical impact and the potential risk of damaging subsurface storage and sealing formations, or even surface facilities. With such pressure control, the subsurface mechanical responses are dominated by thermal strain and stress, which can be monitored with subsurface fiber optics. The study demonstrates how coupled THM modeling can be applied for the design of a safe and efficient thermal energy storage operation, and how subsurface fiber optic monitoring can be applied for performance confirmation, allowing for more confident operational forecasting.

Rutqvist, Jonny↗
DOE OSTI2026

Arm and shoulder muscle segmentation in axial MRI with UNet deep learning model

Quantifying individual upper-limb muscle volumes from MRI provides key insight into muscle-specific strength, deficits, and adaptations. Manual delineation is the gold standard but time‑intensive, and the performance of current deep learning approaches, particularly for small or anatomically complex muscles, remains incompletely characterized. We evaluated a state‑of‑the‑art deep learning framework across the entire upper limb and analyzed factors governing segmentation performance, with attention to the forearm. Three previously published MRI datasets (1.5 T, 3D GRE T1‑weighted; total n = 39) spanning young, middle‑aged, and older adults were curated and quality‑checked, including expert manual segmentations for 31 muscles. Following multiclass mask reconstruction, we trained three 3D nnU‑Net multiclass models matched to the muscle subsets present across datasets, using five‑fold cross‑validation and a composite Dice Similarity Coefficient (DSC) + cross entropy loss. Segmentation accuracy was assessed with DSC. Performance varied across muscles (mean DSC = 0.806 ± 0.098), ranging from 0.920 (Deltoid) to 0.461 (Extensor pollicis brevis). In uncertainty‑weighted regressions, muscle volume was positively associated with DSC (R2 = 0.36, p < 0.001), whereas training segmentation count and muscle orientation showed negligible associations (R2 ≤ 0.06). A weighted mixed‑effects model identified volume as the strongest evaluated predictor, explaining 23.9% of variance in DSC; orientation and training count each contributed <1%, leaving 61.5% unexplained. These results indicate that deep learning–based segmentation can accurately quantify muscle volume for many upper‑limb muscles but remains constrained for small, low‑contrast forearm muscles.

Gillespie, Samuel↗
DOE OSTI2026

TANTE: Time-adaptive operator learning via neural Taylor expansion

Operator learning for time-dependent partial differential equations (PDEs) has seen rapid progress in recent years, enabling efficient approximation of complex spatiotemporal dynamics. However, most existing methods rely on fixed time step sizes during rollout, which limits their ability to adapt to varying temporal complexity and often leads to error accumulation. In this work, we propose the Time-Adaptive Transformer with Neural Taylor Expansion (TANTE), a novel operator-learning framework that produces continuous-time predictions with adaptive step sizes. TANTE predicts future states by performing a Taylor expansion at the current state, where neural networks learn both the higher-order temporal derivatives and the local radius of convergence. This allows the model to dynamically adjust its rollout based on the local behavior of the solution, thereby reducing cumulative error and improving computational efficiency. We demonstrate the effectiveness of TANTE across a wide range of PDE benchmarks, achieving superior accuracy and adaptability compared to fixed-step baselines, delivering accuracy gains of 60-80 % and speed-ups of 30-40 % at inference time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Re-evaluating probable maximum precipitation estimates: sensitivity to transposition domains and storm rotation using modern datasets

This study examines the sensitivity of Probable Maximum Precipitation (PMP) estimates to key methodological decisions embedded in the legacy approach adopted in the U.S. National Weather Service Hydrometeorological Reports No. 51 and No. 52. Although widely used for infrastructure design and risk regulation, fundamental aspects of PMP estimation—such as storm sample size, transposition domain, maximization procedures, and storm rotation—remain poorly constrained and lack formal guidance. Using the Red Rock watershed in Iowa as a case study, and leveraging the 2002–2023 NOAA Analysis of Record for Calibration (AORC) precipitation dataset, we systematically evaluate how each methodological choice, individually and in combination, influences PMP estimates. Our findings demonstrate that PMP is not a fixed physical upper bound but rather a modeling construct shaped heavily by user-defined assumptions. Notably, PMP values derived from modern gridded rainfall datasets can be substantially higher than the legacy estimate used in the original spillway design for Red Rock Dam. Decisions regarding storm sample size, domain extent, climatological window, and particularly storm rotation all contributed to higher PMP estimates. Storm rotation alone—a loosely constrained element in the current PMP practice—can amplify PMP by more than 25%. These results reveal the lack of standardized bounds in current PMP workflows and the need for systematic sensitivity and uncertainty analysis. As PMP estimation shifts toward probabilistic approaches, incorporating physically meaningful storm attributes will be key to developing more transparent, defensible methods for dam safety and climate-resilient infrastructure.

Probable maximum precipitation↗
DOE OSTI2026

Magnetic structure of topological nodal-line semimetal EuSb 2

Here, we investigated the magnetic structure of the magnetic topological nodal-line semimetal EuSb 2 using x-ray resonant magnetic scattering at the Eu 𝐿 3 absorption edge. Forbidden magnetic reflections confirm antiferromagnetic ordering with propagation vector 𝐪 = (0.5,0,0) below 𝑇 N ≈ 27 K. Azimuthal-angle-dependent measurements reveal that the Eu moments are aligned along the monoclinic 𝑏 axis. Comparison between the measured 𝐐-dependent magnetic Bragg peak intensities and calculations based on representation analysis identifies the magnetic structure as the A-type antiferromagnetic 𝛤 4 representation with a (++−−) arrangement along the 𝑎 direction. These results establish the microscopic magnetic structure of EuSb 2 and provide a basis for understanding the interplay between antiferromagnetism and topological electronic states in this system.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Impact of Crystalline Phases on Low-Activity Waste Glass Durability: Insights from PCT and VHT

During vitrification of nuclear wastes, slow cooling along the container centerline promotes crystalline phase formation, which can alter residual glass composition and reduce chemical durability. This study investigates the effects of crystalline phases on the chemical durability of low-activity waste (LAW) borosilicate glasses using the product consistency test (PCT) and vapor hydration test (VHT) on container centerline cooled (CCC) samples. A preliminary model (R2 = 0.88) was developed to predict CCC PCT responses based on glass composition, PCT data from quenched glasses, and measured crystal fractions. Using the latest LAW glass dataset, the feasibility of predictive modeling is evaluated, limitations in current data and methods are identified, and challenges for improving model accuracy are discussed to guide future data collection and model development.

borosilicate glass↗
DOE OSTI2026

Stoichiometric effects on grain growth in zirconium carbide coatings for high-temperature nuclear fuel

Interest in coated particle fuel for space nuclear propulsion (SNP) has expanded in recent years due to successful demonstrations of the resiliency of the coatings to extreme environments. For SNP applications, the coating layer for the particle design needs to be able to withstand exposure to high temperature hydrogen during operating conditions. ZrC has been proposed as a protective layer, however, it is important to understand the high temperature behavior to ensure adequate protection to this fuel. In this study, surrogate ZrC coated particles were heat treated at 1900 °C up to 300 min, to examine how the microstructure evolves when exposed to high temperature. Scanning electron microscopy and electron backscatter diffraction (EBSD) were conducted to determine grain size and grain boundary character and orientation to determine the degree of change in the ZrC layer post heat treatment. Raman spectroscopy provided insight to understand how the as-fabricated stoichiometry of each sample contributed to the differences in grain growth behavior. Despite the as-fabricated samples showing a similar initial grain size and grain boundary character, the samples with a higher amount of excess carbon exhibited smaller grain areas and slower growth rates when exposed to 1900 °C. In conclusion, this investigation details the as-fabricated microstructure of the ZrC layer, specifically grain size, evolved under high temperature as this can impact the performance of the fuel under operating conditions.

EBSD↗
DOE OSTI2026

Effects of high-dose neutron irradiation at light-water reactor relevant temperature on the mechanical properties of SiC/SiC composites

For this study, the neutron dose-dependent evolutions and the underlying mechanisms of properties of SiC fiber-reinforced SiC matrix (SiC/SiC) composites at a temperature relevant to light-water reactors (∼600 K) were investigated and analyzed. Chemical vapor infiltrated (CVI) SiC/SiC composites reinforced with Hi-Nicalon Type S or Tyranno SA3 fiber were neutron-irradiated to doses up to 30 dpa. The irradiated composites retained their flexural strengths. The thermal diffusivity and dimensional changes were mostly retained from 2.0 to 30.2 dpa. Discrepancies in irradiation responses among CVI SiC/SiC composites from different sources were found. Additional microstructural analysis using Raman spectroscopy and numerical analysis on irradiation effect on residual stress were used to explain how the microstructural variables, especially of carbon interphases, affect the mechanical properties in the 30–40 dpa dose range.

Continuous fiber-reinforced ceramic matrix composi↗
DOE OSTI2026

Comparative neutron-irradiation effects on thermal conductivity degradation and dimensional stability of TiC, TiB 2 , and ZrB 2 at 200–1000 °C

Ultra-high-temperature ceramics (UHTCs), including TiC, Ti 11 B₂, and Zr 11 B₂, show great potential for plasma-facing components due to their excellent high-temperature properties prior to irradiation. However, their response to neutron irradiation remains insufficiently understood, limiting robust assessment of their viability for fusion energy applications. Here, this study examines the thermal conductivity, dimensional stability and microstructure of TiC, TiB₂, and ZrB₂ following neutron irradiation at temperatures of 200–1000 °C and fast neutron fluences of 2.0 × 10 25 to 1.1 × 10 26 n/m 2 (E > 0.1 MeV). Lattice swelling measured by synchrotron X-ray diffraction in all three UHTCs was maximized at 200 °C and decreased with increasing irradiation temperature, with no evidence of amorphization observed at 200 °C. Above 600 °C, significant macroscopic volume swelling was observed in irradiated Ti 11 B₂ and Zr 11 B₂, but not in TiC, likely due to cavity formation in the diborides. The post-irradiation thermal conductivity, measured at the irradiation temperature, ranged from 28 to 45 W/m·K, representing a 34–45% reduction relative to the unirradiated material. Notably, neutron-irradiated UHTCs exhibit recoverable thermal conductivity at elevated temperatures, comparable to ferritic–martensitic steels and potentially superior to W when transmutation effects are considered, highlighting promise for shielding or armor plasma-facing components. At 600 °C, both thermal conductivity degradation and lattice swelling saturated at doses exceeding 2–4 dpa.

fusion materials↗