SEARCH · Engineering Papers

Results for “iterative methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Predicting Li-Ion Battery Capacity Fade Using Early-Life Data and a Hybrid Data-Driven Gaussian Process-Bayesian Regression Approach

Accurately predicting Li-ion battery capacity trajectories using early-life data can dramatically improve battery-life understandings and be used to rapidly evaluate design/cost/performance trade-offs when developing new battery materials. Accurate early-life predictions enable researchers to quickly iterate over cell designs and material precursor properties without consistently cycling cells to failure. To this end, we present a toolbox that uses a combined Gaussian Process and Bayesian regression approach that capitalizes on signals other than just capacity (e.g., dQ/dV, voltage drops) to rapidly predict capacity-fade trajectories. The prediction tool uses Bayesian regression to fit functional forms, e.g., power law, sigmoids, etc., to predict capacity-fade dynamics. By fitting functional forms, the capacity fade can be interrogated at any point in the future, allowing for early cell-failure prediction. Additionally, Bayesian regression allows for accurate uncertainty estimates that account for cell-to-cell variability (aleatoric uncertainty) and the lack of observation data (epistemic uncertainty). By only using early cycle data to predict the capacity fade trajectory, uncertainty bounds at end-of-life can be extremely large. The large uncertainty bounds are further exacerbated because there is no systematic way to define the prior distribution of the functional forms' parameters. We improve our the predicted trajectory confidence interval of our predicted trajectory using two methods. First, we shows that a small amount of held-out cycling data is sufficientuse some train cells, that have been cycled to failure to derive information regarding the appropriate prior distributions for the functional forms' parameters of the functional form, effectively leading to data-driven priors.. We propose constructing the data-driven priors by first running a Bayesian regression starting with uninformed priors to generate intermediate cell-specific posterior parameter distributions. These posterior distributions are combined using a Ggaussian mixture model for each parameter to create the data-driven priors. These mixture models serve as the data-driven prior distributions for the parameters for. Second, we derive multiple features, e.g., C_dchg 0.5 DoD 0.5, log (|mean(dQ/dV_(w_3-w_0 ) (V)|), etc., from the train cellsheld-out cycling data, identify which the features are that best predicting capacity at early/mid-life cycles, and then create Ggaussian process regression models that are used for predicting capacity at early/mid-life cycles for the test cells (see blue dots with error bars in Fig 1b). Finally, these predicted data-points are used in addition to the actual early cycle data capacity fade to construct the Bayesian regression trajectory for the test cell s. Notably. We note that these two methods are complementary and can be combined with each other. We evaluate the performance of our proposed method on an testing open-source dataset from Iowa State University and Iowa Lakes Community College (ISU-ILCC). This dataset comprises of 251 nickel-manganese-cobalt/graphite Lithium-ion cells that are cycled under 63 different conditions. We compute the mean average percentage error (MAPE) and negative log predictive density (NLPD) to quantify the efficacy of our method. Our initial findings suggest that, when only few observations are available, for test cells, when using only Bayesian regression with uninformed priors, a power law functional provides the most accurate predictions. with very few data points. However, asHowever, a the number of data points increases, a twin sigmoidal function becomes more accurate as the number of observations further increases. We also find that using as little as 10% of the data set towards generating data-driven priors can lead to significant improvement in prediction accuracy when using early cycle data. Lastly, we found that augmenting early-cycle data with Gaussian process-predicted capacity data for Bayesian regression greatly improves the prediction accuracy. We will present a comprehensive comparison of our methods to other methods available in the literature and apply this method to additional battery datasets.

42 ENGINEERING↗

DOE OSTI2025

$\mathrm{SageNet}$: Fast Neural Network Emulation of the Stiff-amplified Gravitational Waves from Inflation

Accurate modeling of the inflationary gravitational waves (GWs) requires time-consuming, iterative numerical integrations of differential equations to take into account their backreaction on the expansion history. To improve computational efficiency while preserving accuracy, we present the Stiff-amplified Gravitational-wave Emulator Network (SageNet), a deep learning framework designed to replace conventional numerical solvers (code available at https://github.com/YifangLuo/SageNet). SageNet employs a long short-term memory architecture to emulate the present-day energy density spectrum of the inflationary GWs with possible stiff amplification, Ω GW (f). Trained on a data set of 25,689 numerically generated solutions, SageNet allows accurate reconstructions of Ω GW (f) and generalizes well to a wide range of cosmological parameters; 90.9% of the test emulations with randomly distributed parameters exhibit errors of under 4%. In addition, SageNet demonstrates its ability to learn and reproduce the artificial, adaptive sampling patterns in numerical calculations, which implement denser sampling of frequencies around changes in spectral indices in Ω GW (f). The dual capability of learning both physical and artificial features of the numerical GW spectra establishes SageNet as a robust alternative to exact numerical methods. Finally, our benchmark tests show that SageNet reduces the computation time from tens of seconds to milliseconds, achieving a speedup of ∼10 4 times over standard CPU-based numerical solvers with the potential for further acceleration on GPU hardware. These capabilities make SageNet a powerful tool for accelerating Bayesian inference procedures for extended cosmological models. In a broad sense, the SageNet framework offers a fast, accurate, and generalizable solution to modeling cosmological observables whose theoretical predictions demand costly differential equation solvers.

Astronomy data modeling↗

DOE OSTI2021

Flattening the Radial Temperature Profile across the Transformational Challenge Reactor Core

The Transformational Challenge Reactor (TCR) program is demonstrating an agile development approach to advanced nuclear reactor design, which has traditionally utilized a linear design process. In leveraging artificial intelligence, additive manufacturing, advanced materials, and cutting-edge modeling and simulation, the TCR program aims to minimize the high cost and lengthy deployment timelines now standard in the nuclear industry. Within a relatively short period of time, a robust and mature advanced gas-cooled reactor was iteratively designed under the TCR program, using these cutting-edge technologies. The TCR is a 3 MWt gas-cooled microreactor fueled with uranium nitride (UN) tristructural isotropic (TRISO) fuel particles. Though manufactured via traditional means, these UN TRISO particles are loaded into additively manufactured silicon carbide (SiC) cans [4]. Once loaded with TRISO particles, the SiC cans are densified using a chemical vapor infiltration process. The additively manufactured SiC enables significantly more freedom in the design of the fuel form than could ever be achieved using traditionally manufactured SiC. The helium coolant, pressurized to 5 MPa, enters the core at 300°C and nominally exits it at 500°C. Typically, the most thermally limiting components in any reactor design are the fuel assemblies in the core center. To provide a wide thermal margin in these central fuel assemblies, the flow may be biased toward the center of the core to more effectively cool these fuel assemblies with more power deposition and flatten the core’s radial temperature distribution. An analytical fluid model of the TCR core was developed to explore methods for biasing the flow away from the cooler outer fuel assemblies and towards the hotter inner ones. Higher-fidelity models developed in STAR-CCM+ 2020.3.1, a computational fluid dynamics code, were then utilized to verify the analytical model’s findings.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Optimization of the number and locations of the calibration stations needed to monitor soil moisture using distributed temperature sensing systems: A proof-of-concept study

The single-probe heat-pulse (SPHP) technique combined with the Fiber-optic Distributed Temperature Sensing (DTS) technology can offer novel high-resolution measurements of soil moisture (θ) over spatial scales ranging from several centimeters to several kilometers. However, the key limitation of this method is in obtaining the calibration relationship between θ and soil thermal conductivity (λ) across a specific field. In a previous study, a new methodology using a Gaussian processes model was presented to account for the spatial variability in the λ-θ relationship. The model aggregated θ measurements from soil moisture sensors scattered over the SPHP transect with the corresponding DTS λ measurements at their locations. In this study, a novel methodology is tested to optimize the number and locations of soil moisture sensors required to account for the spatial variability of the λ - θ relationship to achieve higher accuracy from the SPHP technique. The proposed methodology utilizes hierarchical clustering to analyze the information contained in the spatial structure of the SPHP measurements as the soil dries from a nearly-saturated condition. The proposed methodology was tested using data from a field in Oklahoma. Monte-Carlo simulation was performed to validate the performance of the proposed methodology. The predictions obtained from the proposed methodology resulted in θ measurements accuracy comparable to those obtained from the 10% best Monte-Carlo iterations of randomly assigned soil moisture locations. Further, this study demonstrates that the proposed methodology is more efficient than the traditional practice of randomly spreading calibration soil moisture sensors along the SPHP transect.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Iterated projected position algorithm for constructing exponentially localized generalized Wannier functions for periodic and nonperiodic insulators in two dimensions and higher

Localized bases play an important role in understanding electronic structure. In periodic insulators, a natural choice of localized basis is given by the Wannier functions which depend on a choice of unitary transform known as a gauge transformation. Over the past few decades, there have been many works that have focused on optimizing the choice of the gauge so that the corresponding Wannier functions are maximally localized or reflect some symmetry of the underlying system. In this work, we consider fully nonperiodic materials where the usual Wannier functions are not well defined and gauge optimization is impractical. To tackle the problem of calculating exponentially localized generalized Wannier functions in both periodic and nonperiodic systems, we discuss the ‘iterated projected position (IPP)” algorithm. Here, the IPP algorithm is based on matrix diagonalization and therefore unlike optimization-based approaches, it does not require initialization and cannot get stuck at a local minimum. Furthermore, the IPP algorithm is guaranteed by a rigorous analysis to produce exponentially localized functions under certain mild assumptions. We numerically demonstrate that the IPP algorithm can be used to calculate exponentially localized bases for the Haldane model, the Kane-Mele model (in both Z 2 invariant even and Z 2 invariant odd phases), and the p x + ip y model on a quasicrystal lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements

Here, we present a computationally efficient approach to perform systematically convergent real-space all-electron Kohn-Sham density functional theory calculations for solids using an enriched finite element (FE) basis. The enriched FE basis is constructed by augmenting the classical FE basis with atom-centered numerical basis functions, comprising of atomic solutions to the Kohn-Sham problem. Notably, to improve the conditioning, we orthogonalize the enrichment functions with respect to the classical FE basis, without sacrificing the locality of the resultant basis. In addition to improved conditioning, this orthogonalization procedure also renders the overlap matrix block diagonal, greatly simplifying its inversion. Subsequently, we use a Chebyshev polynomial based filtering technique to efficiently compute the occupied eigenspace in each self-consistent field iteration. We demonstrate the accuracy and efficiency of the proposed approach on periodic unit cells and supercells. The benchmark studies show a staggering 130× speedup of the orthogonalized enriched FE basis over the classical FE basis. We also present a comparison of the orthogonalized enriched FE basis with the linearized augmented plane-wave + local orbitals basis, both in terms of accuracy and efficiency. Notably, we demonstrate that the orthogonalized enriched FE basis can handle large system sizes of ~10 000 electrons. Finally, we observe good parallel scalability of our implementation with 92% efficiency at 22× speedup for a system with 620 electrons

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Active Learning Framework

Machine learning (ML) of interatomic potentials show great promise to accelerate scientific simulation, e.g., by emulating expensive computations at a high accuracy but much reduced computational cost. Training datasets are calculated from computationally expensive ab initio quantum mechanics methods, density functional theory (DFT). Trained on this data, an ML model can be very successful in predicting energy and forces for new atomic configurations. A critical factor is the quality and diversity of the training dataset. Thus, a highly automated approach to dataset construction based on active learning framework is designed suitable for material physics. The active learning scheme begins with fully randomized atomic configurations. Then, many Molecular Dynamics (MD) trajectories are simulated using current ML potentials, where each MD trajectory is initialized to a random disordered configuration. The temperature is varied in order to diversify the sampled configuration during these simulations. The variance of predictions for eight neural networks within an ensemble is analyzed to determine whether the model is operating as expected. This helps in determining whether collecting more data would be helpful to the model by checking the ensemble variance is greater than the threshold. In this case, the MD trajectory is terminated and the final atomic configuration is placed on a queue (SQL database) for DFT calculations and added to training dataset. Periodically, ML model is retrained to the updated training model. This Active Learning loop is iterated until the cost of MD simulations becomes prohibitively expensive. The MD simulations will hopefully be sufficiently robust to support nucleation after many active learning iterations. In this sense, active learning scheme must automatically discover the important low energy and nonequilibrium physics.

Nebgen, Benjamin↗

DOE OSTI2023

Sustainable bioleaching of lithium-ion batteries for critical metal recovery: Process optimization through design of experiments and thermodynamic modeling

Recycling spent lithium-ion batteries (LIBs) could alleviate supply risks for critical metals and be less harmful to the environment compared to new production of metals from mining. Developing a cost-effective LIB bioleaching process could be a promising alternative to traditional energy-intensive recycling technologies. Here, this study aimed to optimize bioleaching conditions for maximum economic competitiveness through design of experiments using iterative response surface methodology (RSM), assisted by thermodynamic modeling. The optimal condition was identified as 2.5% pulp density in 75 mM gluconic acid biolixiviant at 55°C for 30 h which could recover 57%–84% of nickel, 71%–86% of cobalt, and 100% of lithium and manganese, yielding a 17%–26% net profit margin. The recommended pulp density and acid concentrations, together with the observed metal solubilization, were supported by thermodynamic modeling predictions. Our study demonstrated that combining RSM with thermodynamic simulations could be a powerful tool for optimizing bioleaching conditions.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Neural network-based classification and regression of magnetohydrodynamic modes in tokamaks

We present a machine learning-based magnetohydrodynamic (MHD) classifier and regressor that utilizes real or complex-valued 3D magnetic sensor array data to determine neoclassical tearing mode (NTM) onset times in tokamaks with millisecond accuracy. The input dataset consists of poloidal profiles of complex Fourier amplitudes with an n = 1 toroidal mode number from 144 human-labeled ITER Baseline Scenario discharges in the DIII-D tokamak, spanning both tearing-dominated and sawtooth-dominated regimes. Since m, n = 2,1 NTMs frequently emerge alongside sawteeth at the same frequency in this scenario, the focus is on isolating the m = 1 and m = 2 components of the n = 1 MHD mode near the tearing onset. To improve model regularization and prediction stability, singular value decomposition was applied to balance the sawtooth and tearing datasets. The enriched datasets facilitated training neural networks that learn the key distinguishing features of sawtooth and tearing modes in the poloidal profiles of their magnetic amplitude and phase. When the modes occur independently, the networks achieve perfect classification due to the modes’ distinct characteristics and low measurement noise. In the more experimentally relevant case where both modes coexist, the networks maintain exceptional performance across key metrics. Tests on synthetic data with known ground truth demonstrate the superior accuracy of the neural network trained on complex-valued input compared to models using real amplitude, phase, or pseudo-complex data, achieving both a mean time delay and standard deviation below 1 ms. Notably, standard linear regression methods fitting the dominant singular modes to the data closely match the neural network’s performance. Applying these methods across a broad range of H-mode scenarios will enable future studies to systematically identify dominant NTM triggers as scenario-specific variables, paving the way for more effective tearing mode avoidance strategies in future fusion reactor designs.

machine learning↗

DOE OSTI2025

Terra-Populus v0.1: A Python Library for LandScan High-Definition Population Analysis and Modeling

The terra-populus library is designed for use by the LandScan HD technical team, offering a streamlined set of tools for generating and updating LandScan HD datasets from foundational building-level data, referred to as 'molecules,' provided by the building-level attribution team. This document serves as the primary technical documentation for terra-populus. Version 0.1 of the library includes the core modeling components necessary for LandScan HD production. It enables the generation of the LandScan HD Baseline dataset as well as corresponding confidence measures for the occupancy rates used. Parameters have been included for incorporating damaged building indicators and changes in population, to faciliate the creation of rapid updates for LandScan HD. Future iterations of terra-populus will introduce tools for creating a confidence index, and quantifying and propagating uncertainty, facilitating the creation of probabilistic LandScan HD outputs. This report provides an overview of the tools available in the library and the corresponding code implementations. One of the key advancements implemented in terra-populus is a redefinition of the atomic modeling unit for LandScan HD. Traditionally, the LandScan HD vector analytical framework has generated population estimates at the building sub-component (molecule) level. However, terra-populus adopts a building-level modeling approach. This shift is an operational decision aimed at aligning LandScan HD outputs with confidence measures, which are computed and validated at the building level (confidence measures are not included in this version of terra-populus, aside from those associated with the occupancy rates). Additional advancements to the LandScan HD modeling, as implemented by terra-populus, include a minimum population value parameter and an auto assignment of building floor counts. The population minimum value was implemented to prevent buildings and subsequent LandScan HD pixels that contained small values that may not rasterize in production. An 'auto' value has been included as a method for dealing with buildings lacking floor count information, where it is the average floor count of all other buildings with a residential building use type tag. The logic behind this is to remain consistent with the current logic employed for dealing with building use type null instances, where a null use type is defaulted to residential since it is the most common building type. The auto logic is intended to apply the most common building floor count of the most common type of buildings. The tools provided in terra-populus represent a significant step forward in improving the efficiency, reproducibility, and transparency of the LandScan HD modeling process. As the library evolves, it will continue to serve as a foundational resource for high-resolution population modeling.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

DOE OSTI2025

Autonomous organic synthesis for redox flow batteries via flexible batch Bayesian optimization

Traditional trial-and-error methods for materials discovery are inefficient to meet the urgent demands posed by the rapid progression of climate change. This urgency has driven the increasing interest in integrating robotics and machine learning into materials research to accelerate experimental learning. However, idealized decision-making frameworks to achieve maximum sampling efficiency are not always compatible with high-throughput experimental workflows inside a laboratory. For multi-step chemical processes, differences in hardware capacities can complicate the digital framework by introducing constraints on the maximum number of samples in each step of the experiment, hence causing varying batch sizes in variable selection within the same batch. Therefore, designing flexible sampling algorithms is necessary to accommodate the multi-step synthesis with practical constraints unique to each high-throughput workflow. In this work, we designed and employed three strategies on a high-throughput robotic platform to optimize the sulfonation reaction of redox-active molecules used in flow batteries. Our strategies adapt to the multi-step experimental workflow, where their formulation and heating steps are separate, causing varying batch size requirements. By strategically sampling using clustering and mixed-variable batch Bayesian optimization, we were able to iteratively identify optimal conditions that maximize the yields. Our work presents a flexible approach that allows tailoring the machine learning decision-making to suit the practical constraints in individual high-throughput experimental platforms, followed by performing resource-efficient yield optimization using available open-source Python libraries.

Tamura, Clara [Univ. of Washington, Seattle, WA (U↗

DOE OSTI2025

Resolving the Coverage Dependence of Surface Reaction Kinetics with Machine Learning and Automated Quantum Chemistry Workflows

Microkinetic models for catalytic systems require estimation of many thermodynamic and kinetic parameters that can be calculated for isolated species and transition states using ab initio methods. However, the presence of nearby coadsorbates on the surface can dramatically alter these thermodynamic and kinetic parameters causing them to be dependent on species coverage fractions. As there are combinatorially many coadsorbed configurations on the surface, computing the coverage dependence of these parameters is far less straightforward. We present a framework for generating and applying machine learning models to predict coverage-dependent parameters for microkinetic models. Our toolkit enables automatic calculation and evaluation of coadsorbed configurations allowing us to sample 2,000 coadsorbed adsorbates and transition states (TSs) for a diverse set of 9 reactions on Cu(111), a challenging surface, with four possible coadsorbates. This dataset was then used to train subgraph isomorphic decision trees (SIDTs) to predict the stability and association energy of configurations. We were able to achieve mean absolute errors (MAEs) of 0.106 eV on adsorbates, 0.172 eV on TSs, and due to natural error cancellation in SIDTs for relative properties, 0.130 eV on reaction energies and 0.180 eV on activation barriers. In conclusion, we describe how to use these models to predict coverage-dependent corrections for adsorbates and TSs and demonstrate on H*, HO*, and O* comparing the generated SIDT model with an iteratively refined version.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

The Grassmann interpolation method for spin-unrestricted open-shell systems

The recently reported Grassmann interpolation (G-Int) method [J. A. Tan and K. U. Lao, J. Chem. Phys. 158, 051101 (2023)] has been extended to spin-unrestricted open-shell systems. In contrast to closed-shell systems, where G-Int has to be performed only once since the α and β density matrices are the same, spin-unrestricted open-shell systems require G-Int to be performed twice—one for the α spin and another for the β spin density matrix. In this work, we tested the performance of G-Int to the carbon monoxide radical cation CO●+ and nickelocene complex, which have the doublet and triple ground states, respectively. We found that the Frobenius norm errors associated with the interpolations for the α and β spin density matrices are comparable for a given molecular geometry. These G-Int density matrices, when used as an initial guess for a self-consistent field (SCF) calculation, outperform the conventional SCF guess schemes, such as the superposition of atomic densities, purified superposition of atomic densities, core Hamiltonian, and generalized Wolfsberg–Helmholtz approximation. Depending on the desired accuracy, these G-Int density matrices can be used to directly evaluate the SCF energy without performing SCF iterations. In addition, the spin-unrestricted G-Int density matrices have been used for the first time to directly calculate the atomic charges using the Mulliken and ChElPG population analysis.

Chemistry↗

DOE OSTI2020

Extremum seeking for optimal control problems with unknown time-varying systems and unknown objective functions

We consider the problem of optimal feedback control of an unknown, noisy, time-varying, dynamic system that is initialized repeatedly. Examples include a robotic manipulator which must perform the same motion, such as assisting a human, repeatedly and accelerating cavities in particle accelerators which are turned on for a fraction of a second with given initial conditions and vary slowly due to temperature fluctuations. In this paper, we present an approach that applies to systems of practical interest. The method presented here is model independent; does not require knowledge of the objective function; is robust to measurement noise; is applicable for any set of initial conditions; is applicable to simultaneously controlling an arbitrary number of parameters; and may be implemented with a broad range of continuous or discontinuous functions such as sine or square waves. For systems with convex cost functions we prove that our algorithm will produce controllers that approach the minimal cost. For linear systems we reproduce the cost minimizing linear quadratic regulator optimal controller that could have been designed analytically had the system and cost function been known. We demonstrate the effectiveness of the algorithm with simulation studies of noisy and time-varying systems.

42 ENGINEERING↗

DOE OSTI2022

Fully Distributed Multi-Material Magnetic Sensing Structures for Multiparameter DAS Applications

This dissertation demonstrates the first of its kind distributed magnetic field sensor based on a fiber optic distributed acoustic sensing (DAS) scheme. Ferromagnetic nickel and Metglas® were dispersed internally within a fiber optic preform and then drawn on an in-house fiber optic draw tower to lengths in the kilometers. Due to the close proximity of the ferromagnetic metals and fiber optic core, the magnetostrictive strain response of the ferromagnetic materials when exposed to a magnetic field would perturbate within the fiber cladding and transfer that strain, internally, to the fiber optic core. Strain resulting from the magnetostrictive effect allows the DAS based sensor to accurately translate strain into readable magnetic field data. Due to the high sensitivity seen in this sensor design, multiparameter sources, acoustic and magnetic fields, were tested and validated and a three dimensional magnetic-field vector sensor was proposed. Numerical analysis of the novel sensor design was first implemented using COMSOL Multiphysics, where inputs such as magnetostrictive element shape, size, distance, and number were first investigated. Upon optimizing system constraints, the sensor design was further modified such that single mode operation was consistent across multiple fiber draws while retaining high strain transfer from the ferromagnetic elements to the fiber optic core. Ferromagnetic material selection was evaluated as a function of the saturation magnetostriction constants and a total of 4 modules were used to fully characterize the complex physics involved in this sensor design. All fabrication and testing were performed in-house using a full scale 3-story fiber draw tower and custom environmental testing stations to imitate naturally occurring events such as magnetic or acoustic point sources. A unique stacking method was used to embed ferromagnetic nickel and Metglas® into a fiber optic preform which when combined with a custom fiber draw process resulted in consistent multi-material fibers drawn to lengths of 1-km. In-house testing facilities included different types of electromagnetic generators, in addition to a soil test bed, and an outdoor test bed which allowed 100 meters of fiber to be tested simultaneously. All tested sensors demonstrated high strain transfer capabilities on the order of 0.01-10 μϵ depending on the materials used, ferromagnetic rod number, and core to metal spacing. Due to the sensitivity of the system the difference between AC and DC was distinct, and directional magnetostriction was studied. Transverse and longitudinal magnetic wave propagation was controlled through a solenoid and rectangular Helmholtz coil, both built in-house. A three-dimensional magnetic field vector sensor was proposed due to the success of the magnetic field sensor, and a design was proposed and initially tested to validate direction as a function of field strength and distance. To summarize, this dissertation explores the first fully distributed magnetic field sensor using DAS based techniques and one of the first multi-material fiber draw processes which can produce consistent single mode fiber up to 1-km. Due to extensive FEA modeling, multiple iterations of the magnetic sensor were fully characterized and an equation describing the relationship between sensor design and strain transfer has been created and validated experimentally. Multi-parameter tests including acoustic and magnetic fields were implemented and an algorithm was developed to separate the mixed signals. Finally, a test was performed to demonstrate the feasibility of sensing magnetic fields directionally. Cumulative results demonstrate a high-quality sensor alternative to current designs which may surpass other magnetic sensors due to innate multi-parameter capabilities, in addition to the inexpensive production cost and extremely long operating lengths.

02 PETROLEUM↗

DOE OSTI2021

Follow the power—pathways to steady-state tokamak reactors

Abstract Diagramming pathways of dimensionless power is a potent method for extrapolating between operating points on present-day tokamaks and future burning plasma devices. The heat transport power, current drive power, H-mode threshold power and other ‘plasma physics’ powers can be expressed in dimensionally correct (or normalized) form as Pa 3/4 , where P is the power and a is the plasma minor radius, with the relative gyroradius ( ρ * ) dependence ranging from gyro-Bohm-like for transport ( P a 3 / 4 ∝ ρ * − 3 / 2 ), Bohm-like for current drive ( P a 3 / 4 ∝ ρ * − 5 / 2 ) and worse than Bohm-like for H-mode threshold ( P a 3 / 4 ∝ ρ * − 3 ). The D–T fusion power cannot be normalized in the same fashion since it is governed by nuclear physics, but at fixed B T it scales like P a 3 / 4 ∝ ρ * − 9 / 2 . Other ‘mixed physics’ powers can be incorporated into the dimensionless power framework by holding B T fixed in the same manner. Diagramming these dimensionless powers vs ρ * shows how the pathway to a steady-state reactor can be optimized relative to various operational boundaries. Using a steady-state hybrid discharge with β N = 3.2 from DIII-D as the starting point, a multi-parameter optimization finds an attractive pathway to steady-state operation on ITER using 76 MW of current drive power (fusion gain of Q fus = 8), along with a pathway to Q fus = 20 in a JET-sized steady-state reactor with B T = 10 T.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

RegularizedOptimization.jl: A Julia framework for regularized and nonsmooth optimization

RegularizedOptimization.jl is a Julia package that implements families of quadratic regularization and trust-region methods for solving the nonsmooth optimization problem $^{\textrm{minimize}}_{𝑥∈ℝ^𝑛}$ 𝑓(𝑥) + ℎ(𝑥) subject to 𝑐(𝑥) = 0, (1) where 𝑓 ∶ ℝ 𝑛 → ℝ and 𝑐 ∶ ℝ 𝑛 → ℝ 𝑚 are continuously differentiable, and ℎ ∶ ℝ 𝑛 → ℝ∪{+∞} is lower semi-continuous. The nonsmooth objective ℎ can be a regularizer, such as a sparsity inducing penalty, model simple constraints, such as 𝑥 belonging to a simple convex set, or can be a combination of both. All 𝑓, ℎ, and 𝑐 can be nonconvex. RegularizedOptimization.jl provides a modular and extensible framework for solving (1), and developing novel solvers. Currently, the following solvers are implemented: • Trust-region solvers TR and TRDH (Aravkin et al., 2022; Leconte & Orban, 2025) • Quadratic regularization solvers R2, R2DH and R2N (Aravkin et al., 2022; Diouane, Habiboullah, et al., 2024) • Levenberg-Marquardt solvers LM and LMTR (Aravkin et al., 2024) used when 𝑓 is a least-squares residual. • Augmented Lagrangian solver AL (De Marchi et al., 2023). All solvers rely on first derivatives of 𝑓 and 𝑐, and optionally on their second derivatives in the form of Hessian-vector products. If second derivatives are not available, quasi-Newton approximations can be used. In addition, the proximal mapping of the nonsmooth part ℎ, or adequate models thereof, must be evaluated. At each iteration, a step is computed by solving a subproblem of the form (1) inexactly, in which 𝑓, ℎ, and 𝑐 are replaced with appropriate models around the current iterate. The solvers R2, R2DH, and TRDH are particularly well suited to solve the subproblems, though they are general enough to solve (1). All solvers are allocation-free, so re-solves incur no additional allocations. To illustrate our claim of extensibility, a first version of the AL solver was implemented by an external contributor. Furthermore, a nonsmooth penalty approach, described in Diouane, Gollier, et al. (2024), is currently being developed, that relies on the library to efficiently solve the subproblems.

Gollier, Maxence [Polytechnique Montréal, QC (Cana↗

DOE OSTI2026

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER↗