Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “successive linear approximation”

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 OSTI2026

Stochastic Microgrid Scheduling With Chance‐Constrained Resilience Consideration

Traditionally, it is assumed that microgrids transition seamlessly from grid‐connected operation to islanded mode in the event of sudden main grid outages. In reality, the islanding process, especially unintentional islanding, is rarely seamless. Instead, it is subject to voltage and frequency fluctuations caused by the instantaneous disconnection of the point of common coupling (PCC) switch, variations in loads and renewable generation output and even the protection tripping of distributed energy resources (DERs). To mitigate these fluctuations and facilitate a smooth islanding process, we propose a stochastic microgrid scheduling model that incorporates chance‐constrained resilience measures. Specifically, the resilience measure is defined as the probability of successful islanding (PSI), that is, the probability that a microgrid can mitigate the generation‐demand imbalance caused by the disconnection of the PCC switch, variations in load and renewable generation and DER tripping. This measure is modelled using chance constraints. Unlike existing reliability and resilience indices, which typically neglect the possibility of microgrid/DER failure under extreme events and assume their survival while primarily focussing on reducing impact duration or magnitude, the proposed PSI‐based framework explicitly addresses microgrid and DER survival during the islanding transition. The formulated nonlinear chance constraints are approximated using a multiinterval approach and equivalently represented as a mixed‐integer linear programming (MILP) formulation. Case study results validate the proposed method, showing that the PSI estimation error is reduced to less than 8%, compared to approximately 28% with existing methods. Various sensitivity analyses on the DER tripping rate and PSI settings were performed to validate the robustness of the proposed method. In particular, the necessity of accounting for DER tripping in the PSI calculation was demonstrated.

chance constrained optimization↗
DOE OSTI2024

Constrained Local Approximate Ideal Restriction for Advection-Diffusion Problems

Herein this paper focuses on developing a reduction-based algebraic multigrid (AMG) method that is suitable for solving general (non)symmetric linear systems and is naturally robust from pure advection to pure diffusion. Initial motivation comes from a new reduction-based AMG approach, $\ell \text{AIR}$ (local approximate ideal restriction), that was developed for solving advection-dominated problems. Though this new solver is very effective in the advection-dominated regime, its performance degrades in cases where diffusion becomes dominant. This is consistent with the fact that in general, reduction-based AMG methods tend to suffer from growth in complexity and/or convergence rates as the problem size is increased, especially for diffusion-dominated problems in two or three dimensions. Motivated by the success of $\ell \text{AIR}$ in the advective regime, our aim in this paper is to generalize the AIR framework with the goal of improving the performance of the solver in diffusion-dominated regimes. To do so, we propose a novel way to combine mode constraints as used commonly in energy-minimization AMG methods with the local approximation of ideal operators used in $\ell \text{AIR}$. The resulting constrained $\ell \text{AIR}$ algorithm is able to achieve fast scalable convergence on advective and diffusive problems. In addition, it is able to achieve standard low complexity hierarchies in the diffusive regime through aggressive coarsening, something that was previously difficult for reduction-based methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A parametric study of slow dynamic nonlinear elasticity with comparisons to models

Several phenomenological models that aspire to quantitative description of anomalous nonlinear mesoscopic elasticity are reviewed and compared with laboratory measurements. This class of nonlinearity, best known perhaps for slow dynamics and aging, is seen widely in imperfectly consolidated granular solids but is not well understood. Typical slow dynamic tests show that a modest conditioning oscillatory "pump" strain depresses material stiffness, which then recovers like the logarithm of time after conditioning ceases. Several phenomenological models based on physical arguments have been proposed that predict the material stiffness response to arbitrary pump strain histories during conditioning and recovery. Approximate closed form and numerical solutions to the models are presented that predict the quantitative influence of three key pump parameters: the pump's strain amplitude, the pump's strain rate, and the pump’s duration. Laboratory measurements on Berea sandstone, concrete and a confined single aluminum bead find that slow dynamic responses are linear in pump strain and independent of pump frequency. Measurements also show that, after pump-off, stiffness recovers over times far longer than the pump duration. These observations and others are compared to model predictions. One of the considered models, based on a picture of fast brittle damage and slow healing, successfully matches all these behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Temperature Evolution of the Activation Barriers Leads to Meyer−Neldel Rules for Structural Relaxation and Transport in Polymers

Understanding activation barriers controlling structural relaxation in glass-forming liquids, molecular transport, and ionic conductivity in amorphous polymers is a grand challenge of fundamental scientific and materials engineering interest across disciplines. Over decades, intriguing but puzzling empirical correlations between the elementary time scale of activated barrier crossing and the apparent Arrhenius activation energy, the so-called Meyer−Neldel (MN) rules, have been discovered in diverse liquids and glasses. Here, in this study, we formulate and successfully apply a new experimental analysis and an explicitly dynamical theoretical framework which provides an understanding of the origin, validity, and failure of such correlations, that bridge and unify the three fields of structural relaxation, molecular transport, and ionic conductivity in liquids and quenched glasses. Distinct quasi-universal laws are predicted in equilibrated liquids and nonequilibrium glasses, consistent with experiments. Our analysis reveals that even if the relaxation appears Arrhenius over a limited temperature range, the physical activation barrier is generally temperature-dependent in polymeric systems even below glass transition temperature. In addition, we show that the approximate validity of classical MN rules hinges on a linear temperature dependence of this barrier and the temperature range probed in experiments. Our findings are relevant for controlling the activation barrier in functional soft polymeric materials relevant to molecular separations, barrier coatings, and charge transport, and also provide new constraints on the theoretical understanding of the mechanism underlying slow activated dynamics in glass-forming condensed matter.

Meyer-Neldel rules↗
DOE OSTI2025

Computational Algorithms for Unit Commitment with AC Power Flows (Final Report)

Security-constrained unit commitment (SCUC) is a key component in power system operations. When AC power flow constraints are considered in the SCUC model (AC-SCUC), the problem becomes extremely difficult due to its discrete and non-convex nature, as described in “Grid Optimization Competition Challenge 3 Problem Formulation (GOCC)”. There are four main challenges: (i) Discrete decisions regarding unit online/offline status and start-up/shut-down procedures for every single unit. The number of discrete decision variables increases considerably when a system integrates multiple generators; (ii) Configuration-based combined-cycle formulations, and multi-commodity models that include ramping products, spin/non-spin products, and regulation up/down products. The combined-cycle units introduce additional discrete decision variables and auxiliary service products further complicate the model by connecting multi-commodity products’ continuous and discrete variables; (iii) SCUC models with AC power flow constraints are far more complex due to massive bilinear terms in the large-scale nonlinear power balance equations. The nonlinear power balance equations are further complicated by the discrete step control variables of shunts; (iv) N − 1 contingency analysis. The size of the model increases linearly with the number of contingencies considered, greatly increasing the size of the optimization model. Accordingly, there is an emergent need to develop a robust algorithm capable of deriving a high-quality solution in a short time and passing through contingency tests simultaneously. In this project, we explore innovative techniques to address this challenging problem by integrating advanced polyhedral theory, approximation methods, relaxation strategies, decomposition techniques, and parallel computing. Each technique approaches the problem from a different perspective, leveraging its specific strengths to tackle distinct challenges. Each individual method has demonstrated its effectiveness in the PI’s previous research. Their integration is expected to significantly reduce the computational time required to solve the proposed complex problem. Successful completion of this project has the potential to transform the industry by enhancing optimization solvers capable of handling large-scale day-ahead energy market clearing models within strict time constraints, while incorporating AC power flow constraints. This advancement will lead to reduced overall generation costs and, consequently, increased social welfare.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Impedance Response Influenced by Variability in the Random Distribution of Physical Properties of Coated Materials in Two-Dimensional Space

Heterogeneous physical characteristics of a system featuring a single-layer film on a metallic surface have been explored via its impedance response. The Nyquist plot showed a distorted semicircle, indicative of the system's unique distribution characteristics. Utilizing a copula-based probability method, a two-dimensional deterministic impedance model was successfully integrated, accounting for spatial physical properties such as permittivity and electrical conductivity. This strategy enabled in-depth exploration and mechanistic quantification of a broad spectrum of properties. A quantitative understanding of impedance signal alterations, characterized by normally or log-normally correlated variables, was achieved through the variation in aspect ratio and characteristic frequency of the impedance spectra. Log-normally distributed electrical properties provided a superior representation of the distorted impedance spectra. As coefficient of variation (CV) values fluctuated, the aspect ratio and characteristic frequency showed heightened sensitivity to log-normal permittivity compared to log-normal electrical conductivity. Notably, a marked positive linear correlation between electrical properties resulted in an impedance response that approximated perfect semicircular spectra. The variability in the electrical properties' distribution was demonstrated by considering the correlation coefficient between electrical conductivity and the z-direction position. Furthermore, the highest aspect ratio of the impedance spectra was observed when the electrical conductivity was randomly distributed across the z-direction space.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Gas Stopper Developments for Improved Purity and Intensity of Low-Energy, Rare Isotope Ion Beams (Final Technical Report)

This final technical report summarizes the work of the Michigan State University (MSU) team supported by grant # DE-SC0021423 awarded by the Office of Nuclear Physics, Department of Energy. Objectives: The successful fulfillment of the FRIB science mission hinges on ensuring the availability of fast, stopped, and reaccelerated beams consisting of rare isotopes. This project's research and development focus was dedicated to supporting the advancement and creation of a cutting-edge linear gas stopper. The primary aim is to efficiently convert the high-intensity fast beams of rare isotopes provided by FRIB into high-quality, low-energy beams. These beams are essential for conducting stopped beam experiments or for subsequent reacceleration. The overarching goal is to advance technology, aiming to increase the beam rate capability of the linear gas stopper for medium-to-heavy-mass rare isotopes by more than tenfold compared to the currently most effective gas stopper in operation, and to improve the purity of the delivered beams. Project Description: The existing technology employed in gas stopping devices designed for low-energy, rare-isotope beams presents limitations in both the purity of extracted beams and the intensities of injected beams. These limitations are incompatible with the requirements of the recently commissioned rare isotope beam facility, FRIB. Our research and development efforts, aligned with the previously outlined objectives, focused on addressing the most critical aspects for enhancing beam-rate capability and purity. Specifically, advanced particle-in-cell simulations were developed and integrated into a simulation pipeline to explore the efficacy of multi-layer RF carpets on increasing ion transport efficiency with high incoming beam rates that generate space charge fields which can limit it. We also explored the possibility of using a collision-induced-dissociation (CID) gas cell to break up molecular contaminant ions that are generated during the stopping process. A prototype CID gas cell was constructed and tested with beams from an offline ion source, validating the concept with the successful demonstration of breaking of molecular ions. The outcome of this research enabled the formulation a conceptual design for a next-generation linear gas stopping device specifically tailored for FRIB. This device is envisioned to deliver rare-isotope-ion beams at a rate of 10 8 particles per second or higher, accompanied by advancements in purity. Methods employed: This project leverages advancements in technologies initially designed for the Advanced Cryogenic Gas Stopper (ACGS), the current state-of-the-art linear gas stopper, through the use of new simulations and beam purification via collision-induced-dissociation. The methods include: 1. Development of a prototype low-energy, low-pressure CID gas stopper. This prototype features a thin, approximately 20 nm, Si 3 N 4 entrance window designed for dissociating stable and rare-isotope molecular ions. The goal is to enhance beam purification and overall efficiency. 2. Creation of Particle-in-cell (PIC) simulations to assess the advantages of multi-layer RF carpets and multi-point extraction for ion transport efficiency. These simulations rely on the 3DCylPIC package, specifically designed for studying devices of this nature. The goal is to quantify and mitigate ion transport losses due to space charge generated in the stopping process of large numbers of ions. 3. Perform ion transport simulations across an RF carpet using an 8-phase travelling wave and evaluate its performance. Compared to the 4-phase RF carpets used in ACGS, the 8-phase carpets will double the wavelength of the generated traveling wave allowing for larger maximum RF amplitudes that could result in improved ion transport efficiency for high-intensity incoming beams when large space charge fields are present. Impact: Tackling the primary challenges associated with transforming high-energy projectile fragment beams into low-energy beams—specifically, addressing efficiency and purity—holds significant promise for advancing FRIB science. This advancement will particularly impact precision mass measurements, laser spectroscopy of short-lived nuclei, and studies in astrophysics and nuclear reactions using reaccelerated beams. These domains play a crucial role in addressing key questions outlined in the 2023 NSAC long-range plan, spanning nuclear structure, nuclear astrophysics, and fundamental symmetries. Additionally, they contribute to addressing 10 out of the 17 benchmarks identified by the NRC RISAC. The development of a next-generation gas stopping device capable of delivering low-energy, rare-isotope beams at a rate of 10 8 particles per second, or more, with high purity holds the potential to unlock experiments that would otherwise be unfeasible. Furthermore, it is expected to reduce the time required for experiments at FRIB, thereby maximizing scientific output. The research and development activities performed as part of this project bolstered essential competencies at FRIB in beam physics and ion source technologies, provided valuable training opportunities for junior scientists.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Deep-freeze graph training for latent learning

Scientific and engineering advances are primarily driven by multi-tier conceptual constructs and conditional theoretical frameworks. The theories allow predictions of hypothetical system responses, given a set of approximate conditions (ranges of applicability) imposed on latent parameters that cannot be measured directly. Learning to estimate the latent variables (Latent Learning) helps to pinpoint the anticipated range-edge anomalies and improves the confidence in interpretation, interpolation and extrapolation of limited experimental data. Due to high dimensionality and extreme non-linearity of the materials science problems, very large datasets are typically required for conventional data-driven model development. The vital experimental data collection, particularly on microstructural phases, is very challenging, which makes it difficult to compile a high-quality database. Incorporation of the domain knowledge into the computational graph structure, initialization and optimization processes presents a viable mechanism for developing accurate models, with limited datasets. Furthermore, this study successfully utilized the approach to build the Deep Freeze Graph (DeepFreG) by mapping known causality relationships and by digitizing empirical domain knowledge for Latent Learning (LL), with specific applications in materials science.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Internal Standard Triggered-Parallel Reaction Monitoring Mass Spectrometry Enables Multiplexed Quantification of Candidate Biomarkers in Plasma

Despite advances in proteomic technologies, clinical translation of plasma biomarkers remains low, partly due to a major bottleneck between the discovery of candidate biomarkers and costly clinical validation studies. Due to a dearth of multiplexable assays, generally only a few candidate biomarkers are tested, and the validation success rate is accordingly low. Previously, mass spectrometry-based approaches have been used to fill this gap but feature poor quantitative performance and were generally limited to hundreds of proteins. Here, we demonstrate the capability of an internal standard triggered-parallel reaction monitoring (IS-PRM) assay to greatly expand the numbers of candidates that can be tested with improved quantitative performance. The assay couples immunodepletion and fractionation with IS-PRM and was developed and implemented in human plasma to quantify 5176 peptides representing 1314 breast cancer biomarker candidates. Characterization of the IS-PRM assay demonstrated the precision (median % CV of 7.7%), linearity (median R 2 > 0.999 over 4 orders of magnitude), and sensitivity (median LLOQ < 1 fmol, approximately) to enable rank-ordering of candidate biomarkers for validation studies. Using three plasma pools from breast cancer patients and three control pools, 893 proteins were quantified, of which 162 candidate biomarkers were verified in at least one of the cancer pools and 22 were verified in all three cancer pools. The assay greatly expands capabilities for quantification of large numbers of proteins and is well suited for prioritization of viable candidate biomarkers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Simulating large one-dimensional neutral-atom quantum systems

While abstract models of quantum computation assume a closed system of two-level states, practical quantum devices inevitably couple to the environment in some way, creating sources of noise. Understanding the tolerance to noise of specific quantum algorithms run on specific devices is important for determining the feasibility of quantum computing in the current noisy intermediate-scale quantum era. Of particular interest is understanding the noise sensitivity of these devices as more qubits are added to the system. Classical simulations are a useful tool to understand the effects of this noise, but direct classical simulations of open quantum systems are burdened by an exponentially growing cost in the number of qubits and a large local Hilbert space dimension. For onedimensional, shallow circuits, using tensor networks can replace this exponential cost with a linear one and simulate far wider systems than what would normally be available. In this paper, we describe a tensor network simulation of a neutral atom quantum system under the presence of noise, while introducing a purity-preserving truncation technique that compromises between the simplicity of the matrix product state and the positivity of the matrix product density operator. We apply this simulation to a near-optimized iteration of the quantum approximate optimization algorithm on a transverse field Ising model in order to investigate the influence of large system sizes on the performance of the algorithm. We find that while circuits with a large number of qubits fail more often under noise that depletes the qubit population, their outputs on a successful measurement are just as robust under Rydberg atom dissipation or qubit dephasing as smaller systems. However, such circuits might not perform as well under coherent multiqubit errors such as Rydberg atom crosstalk. We also find that the optimized parameters are especially robust to noise, suggesting that a noisier quantum system can be used to find the optimal parameters before switching to a cleaner system for measurements of observables.

Allen, James↗
DOE OSTI2023

A Unification Framework for Euclidean and Hyperbolic Graph Neural Networks

Hyperbolic neural networks have recently gained significant attention due to their promising results on several graph problems including node classification and link prediction. The primary reason for this success is the effectiveness of hyperbolic space in capturing the inherent hierarchy of graph datasets. However, they are limited in terms of generalization, scalability, and have inferior performance when applied to non-hierarchical datasets. In this paper, we take a completely different perspective for modeling hyperbolic networks and answer the following question: is an Euclidean model able to approximate a function or behavior in the hyperbolic space? Extending the universal approximation theory developed for Euclidean models, We draw an analogy from the hyperbolic components to the Euclidean counterparts and conclude that, in order to capture hierarchical features, it is possible to generalize hyperbolic models to be a special case of Euclidean models with the proposed Pseudo-Poincaré technique. We applied our non-linear hyperbolic normalization to the current state-of-the-art homogeneous and multi-relational graph networks and demonstrate significant improvements in performance compared to both Euclidean and hyperbolic counterparts. The primary impact of this work lies in its ability to capture hierarchical features in the Euclidean space, and thus, can replace hyperbolic networks without any loss in performance metrics while simultaneously leveraging the power of Euclidean networks such as interpretability and efficient execution of various model components.

Khatir, Mehrdad↗
DOE OSTI2025

PID-Regulated Heating System for PIP-II Reference Line

The Proton Improvement Project-2 centers on building a new superconducting linear particle accelerator (Linac) at Fermilab. At the heart of the accelerator is the reference line, a critical system that defines the ideal path for the particle beam as it passes through magnets, RF cavities, and other beamline elements. Temperature stability is crucial for the reliable operation of RF components, such as mixers and filters. Fluctuations affect key performance parameters like conversion loss, isolation, and linearity. To mitigate any drift caused by ambient temperature changes, a heating plate assembly is utilized to maintain key components at a controlled temperature of 40°C. The system utilizes an aluminum 36”x36”x0.5” heat plate powered by a MOSFET-based control circuit, delivering approximately 460 W of thermal energy through a resistor array. Real-time temperature feedback is provided by a PT100 Resistance Temperature Detector (RTD), which interfaces with a Proportional–Integral–Derivative (PID) control algorithm to maintain closed-loop temperature regulation. The control signal actively modulates the gate voltage of an N channel MOSFET, dynamically adjusting power delivery in response to deviations from the temperature setpoint. Simulations and LTspice models validate the functionality and responsiveness of the circuit under varying conditions. The prototype has successfully demonstrated stable thermal control, paving the way for integration into the PIP-II infrastructure. The final design will feature an expanded resistor array, as well as communication with a PLC for continuous data acquisition and diagnostics. This work directly supports Fermilab’s broader mission by contributing to the stability and reliability of core accelerator systems, enhancing the precision of particle beam delivery for future physics experiments.

Mosher, Alexander [Fermilab]↗
DOE OSTI2023

High-brightness electron injectors for high-duty cycle X-ray free electron lasers

The successful development in the last two decades of X-ray free electron lasers (FELs) with their revolutionary brightness performance has been tightly dependent on the parallel development of electron guns and injectors capable of providing the high-brightness electron beams required by FELs lasing at these short wavelengths. The ultimate brightness delivered by a linear accelerator (linac) is already set at its injector and the remaining part of the accelerator can be only designed to preserve the injector performance. The technology to be used for the accelerator part of an X-Ray FEL strongly depends on the duty-cycle at which the FEL operates. Normal-conducting, room-temperature, copper-based radio frequency (RF) technology is typically used for low duty-cycles of up to approximately 10 −3 . For higher duty-cycles and up to continuous wave (CW) operation, the linac must rely on superconductive RF technology because, with the higher duty-cycle, the increasingly higher power dissipated in normal conducting RF structures becomes excessive for the warm technology. The situation changes in the lower energy part of the accelerator, where injector schemes, based on direct current, normal-conducting, and superconducting RF electron guns, are demonstrating the beam quality performance required by high-duty-cycle X-ray FELs. In this paper we start with a description of the requirements for such injectors, followed by an overview of the pursued technologies and schemes, and by a discussion on the main differences in terms of beam dynamics between low and high duty-cycle injectors.

Sannibale, Fernando↗
DOE OSTI2023

An electrified boom actuation system with energy regeneration capability driven by a novel electro-hydraulic unit

This paper presents a novel electrified solution for linear actuation systems, which has undergone a successful in-vehicle demonstration, marking a significant advancement in this field. Aligned with the growing trend towards electrification, this solution holds promise for a wide range of off-road applications. The showcased electrohydraulic actuation architecture enables efficient four-quadrant operation and incorporates energy recovery capabilities. At the core of the system lies a distinctive electro-hydraulic unit, integrating a fixed hydraulic pump and a variable-speed electric motor into a compact and powerful entity. A specialized hydraulic system is designed to efficiently transform electrical energy from the electro-hydraulic unit into boom cylinder actuation. The sizing of the system is determined by the boom operation of a commercial compact loader. The maximum power output of the electro-hydraulic system is approximately 15 kW. The components of the proposed system are integrated into the reference vehicle, and comprehensive performance tests are conducted. The results of the tests performed on the vehicle with the boom lifting and lowering cycle indicate an overall system efficiency of 60%. The electrified solution consumes only about one-third of the energy used by the conventional boom actuation system on the reference machine, while exhibiting energy savings exceeding 75% during specific work phases. The energy savings come mainly from energy regeneration and minimal loss from throttling through the specialized hydraulic circuit design. Finally, the successful demonstration on the reference vehicle underscores the practicality and effectiveness of the proposed electrified solution in current mobile hydraulic applications

42 ENGINEERING↗
DOE OSTI2026

Variability of MHD instabilities in benign termination of high-current runaway electron beams in the JET and DIII-D tokamaks

Benign termination, in which magnetohydrodynamic (MHD) instabilities deconfine runaway electrons (REs) following hydrogenic injections, is a promising strategy for mitigating dangerous RE loads after disruptions. Recent experiments on the Joint European Torus (JET) have explored this scenario at higher pre-disruptive plasma currents than are achievable on other devices, revealing challenges in obtaining benign terminations at I p ≥2.5 MA. This work analyzes the evolution of these high-current RE beams and their terminating MHD events using fast magnetic sensor measurements and EFIT equilibrium reconstructions for approximately 40 JET and 20 DIII-D tokamak discharges. On JET, unsuccessful non-benign terminations occur at low edge safety factor (q edge ≈ 2), and are preceded by intermittent, non-terminating MHD events at higher rational qedge. Trends in the internal inductance I i indicate more peaked RE current profiles in the high-I p non-benign population, which may hinder successful recombination through re-ionization of the companion plasma. In contrast, benign terminations on JET typically occur at higher q edge ≥3 and exhibit less peaked RE current profiles. DIII-D displays a broader range of terminating edge safety factors, again correlated with the measured values. Across both tokamaks, the RE current peaking is therefore found to determine which MHD instability boundary is encountered, a result confirmed by linear resistive MHD modeling with the CASTOR3D code. Measured growth rates are similar for benign and non-benign cases, indicating that ideal MHD timescales at low density after hydrogenic injection do not alone explain efficient RE deconfinement. Instead, non-benign cases are most readily characterized by their comparably lower overall MHD perturbation amplitudes δB. These observations suggest that the interplay between ideal and resistive dynamics governs the termination process, with implications for extrapolating benign RE termination to high-I p reactor scenarios.

MHD instabilities↗
DOE OSTI2024

High-Fidelity and High-Performance Computational Simulations for Rapid Design Optimization of Sulfur Thermal Energy Storage (CRADA Report)

NREL and Element 16 collaborated on sulfur thermal energy storage modeling using NREL’s high performance computing (HPC) resources to assist its application in industrial processes. Industrial process heat (IPH) accounts for ~70% of US manufacturing energy use and is primarily produced by fossil fuel combustion. Approximately, 1500 TWht (~60% Terawatt hour thermal) of IPH demand is in the temperature range of 100-300°C. Industrial applications in this temperature range include drying, hydrothermal processing, thermal enhanced oil recovery, food and beverage, bioethanol production, etc. Cost-effective thermal energy storage (TES) that increases the utilization of waste and renewable heat (solar, geothermal, etc.) could provide significant energy savings and reliable heat sources, decrease emissions, and increase US manufacturing competitiveness through reductions in fuel consumption. This HPC4EI project facilitated Element 16’s development of low-cost and high-impact molten sulfur TES for dispatchable IPH. The development of a high-fidelity model validated by experimental data and HPC simulations enabled the successful resolution of the complex interplay between fluid dynamics and heat transfer processes during transient operation of sulfur TES, overcoming the numerical challenges posed by the non-linear temperature-dependent physical properties of sulfur. The project helped accelerate Element 16’s molten sulfur TES product design and support its broad applications.

25 ENERGY STORAGE↗
DOE OSTI2025

Randomized Algorithms for Symmetric Nonnegative Matrix Factorization

Symmetric Nonnegative Matrix Factorization (SymNMF) is a technique in data analysis and machine learning that approximates a matrix with a product of a nonnegative, low-rank matrix and it transpose. To design faster and more scalable algorithms for SymNMF we develop two randomized algorithms for its computation. The first method uses randomized matrix sketching to compute an initial low-rank approximation to the input matrix and proceeds to uses this as a low-rank input to rapidly compute a SymNMF. The second methods uses randomized leverage score sampling to approximately solve constrained least squares problems. Many successful methods for SymNMF rely on (approximately) solving sequences of constrained least squares problems. Here, we prove theoretically that leverage score sampling can approximately solve constrained least squares problems to e-accuracy. Finally we demonstrate both methods work in practice by applying them to graph clustering tasks on large real world data sets. These experiments show that our methods approximately maintain solution quality and achieve significant speed ups for both large dense and large sparse problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Multi-Parameter Optical Fiber Sensing of Humidity, CH4, CO2, and Corrosion

In this work, the previously demonstrated capability of the optical fiber sensor (OFS) for successful monitoring of humidity has been extended to monitor the humidity, CH4, and CO2 with different gas composition based on the strain produced along the single mode fiber (SMF) sensor. This is enabled by absorption of H2O/gases on to the commercially available polyacrylate coated jacketed portion of the fiber resulting a change in strain. Under equilibrium, a differential microstrain was observed along the jacketed portion of the SMF with N2, CH4, and CO2 at different humidity conditions, while the unjacketed portion of the fiber was used only for sensing pressure/temperature induced strain. In case of N2 at 800 psig pressure, the observed microstrain was approximately 80, 65, 50 and 35 µε at 100, 75.0, 46.8, and 23.4 RH% respectively. Comparatively, a microstrain of approximately 95 µε was observed with 100% RH CH4 which demonstrates that SMF produces a measurable CH4 response alongside water. Similarly, the observed microstrain with CO2 was approximately 105, 90, 85, 80, and 70 µε at 100, 75.0, 46.8, 23.4, and 0 RH% respectively. Linear regression and principal component analysis of these dataset provided deconvolution of the impact of strain from H2O, CH4, and CO2 enabling good cross sensitivity. Additionally, modified OFS comprising of Fe coated fiber section was employed to monitor corrosion, using Fe as corrosion proxy, under harsh corrosive environment.

Mainali, Badri↗