Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Robust optimization”

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Convergence Analysis of the Alternating Anderson–Picard Method for Nonlinear Fixed-Point Problems

Anderson acceleration (AA) has been widely used to solve nonlinear fixed-point problems due to its rapid convergence. This work focuses on a variant of AA in which multiple Picard iterations are performed between each AA step, referred to as the Alternating Anderson–Picard (AAP) method. Furthermore, despite introducing more “slow” Picard iterations, this method has been shown to be efficient and even more robust in both linear and nonlinear cases. However, there is a lack of theoretical analysis for AAP in the nonlinear case. In this paper, we address this gap by establishing the equivalence between AAP and a multisecant-GMRES method that uses GMRES to solve a multisecant linear system at each iteration. From this perspective, we show that AAP “converges” to the Newton-GMRES method. Specifically, as the residual approaches zero, the multisecant matrix, the approximate Jacobian inverse, the search direction, and the optimization gain of AAP converge to their counterparts in the Newton-GMRES method. These connections provide insights for analyzing the asymptotic convergence properties of AAP. Consequently, we show that AAP is locally 𝑞-linear convergent and provide an upper bound for the convergence factor of AAP. To validate the theoretical results, numerical examples are provided.

Anderson acceleration↗
DOE OSTI2020

UNRAVEL: A Robust Modular Velocity Dealiasing Technique for Doppler Radar

Unfold Radar Velocity (UNRAVEL) is an open-source modular Doppler velocity dealiasing algorithm for weather radars. UNRAVEL is an algorithm that does not need external reference velocity data, making it easily applicable. The proposed algorithm includes 11 core modules and 2 dealiasing strategies. UNRAVEL is an iterative algorithm. The goal is to build the dealiasing results starting with the strictest possible continuity tests in azimuth and range and, after each step, relaxing the parameters to include more results from a progressively growing number of reference points. UNRAVEL also has modules that perform 3D continuity checks. Thanks to this modular design, the number of dealiasing strategies can be expanded in order to optimize the dealiasing results. While the first driver dealiases Doppler velocity from each tilt independently from one another, the second driver also performs a three-dimensional continuity check of the velocity using successive elevations. The proposed dealiasing algorithm is tested using severe weather data from an S-band Doppler radar that have been aliased to mimic aliased radial velocity patterns that would be observed by a C-band Doppler radar. Artificially aliasing S-band data permits creation of a reference to which the performance of various dealiasing techniques can be compared. Comparisons show that UNRAVEL consistently outperforms other established dealiasing algorithms for the test period selected in this work.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Compact fiber-coupled narrowband two-mode squeezed light source

Quantum correlated states of light, such as squeezed states, are a fundamental resource for the development of quantum technologies, as they are needed for applications in quantum metrology, quantum computation, and quantum communications. It is thus critical to develop compact, efficient, and robust sources to generate such states. Here, we report on a compact, narrowband, fiber-coupled source of two-mode squeezed states of light at 795 nm based on four-wave mixing (FWM) in an 85 Rb atomic vapor. The source is designed in a small modular form factor, with two input fiber-coupled beams, the seed and pump beams required for the FWM, and two output fibers, one for each of the modes of the squeezed state. The system is optimized for low pump power (135 mW) to achieve a maximum intensity-difference squeezing of 4.4 dB after the output of fibers at an analysis frequency of 1 MHz. Furthermore, the narrowband nature of the source makes it ideal for atomic-based quantum sensing and quantum networking configurations that rely on atomic quantum memories. Such a source paves the way for a versatile and portable platform for applications in quantum information science.

Jain, Umang [University of Oklahoma, Norman, OK (U↗
DOE OSTI2026

Machine learning of factors for improving oyster hatchery production

Oyster aquaculture and restoration in the Chesapeake Bay are vital, yet hatcheries frequently struggle with inconsistent larval growth and sudden mass mortality events. Unpredictable disruptions in larval production cause large economic losses, represent a perceived risk to growers, and impede industry expansion. To better understand associations between production yield and its potential predictors, we applied machine learning (random forest, and neural network) and statistical (generalized additive model) models to a comprehensive dataset of environmental, water quality, and operational parameters from a Maryland oyster hatchery, aiming to identify key yield predictors and develop a robust forecasting tool. We used recursive Boruta algorithm for variable selection, pinpointing critical predictors, and employed cross-validation to fine-tune model settings. Shapley value analysis offered crucial insights into model interpretations, highlighting week number, Normalized Difference Vegetation Index, salinity, turbidity, and fecundity as primary drivers of yield variability. For low-yield cases, salinity-related variables were particularly important. Our findings provide an early warning system for potential production downturns, empowering hatchery operators to make data-driven decisions for optimizing water conditions, feeding schedules, and broodstock management. By boosting predictability and efficiency, this research directly supports economic stability of the oyster industry and ecological health of the Chesapeake Bay.

Vishwakarma, Srishti [Oak Ridge National Laborator↗
DOE OSTI2025

Efficient on-chip terahertz generation and detection with GaN photoconductive emitters

Photoconductive emitters for terahertz generation hold promise for highly efficient down-conversion of optical photons because it is not constrained by the Manley-Rowe relation. Existing terahertz photoconductive devices, however, faces limits in efficiency due to the semiconductor properties of commonly used GaAs materials. Here, we demonstrate that large bandgap semiconductor GaN, characterized by its high breakdown electric field, facilitates the highly efficient generation of terahertz waves in a coplanar stripline waveguide. Towards this goal, we investigated the excitonic contribution to the electro-optic response of GaN under static electric field both through experiments and first-principles calculations, revealing a robust excitonic Stark shift. Using this electro-optic effect, we developed a novel ultraviolet pump-probe spectroscopy for in-situ characterization of the terahertz electric field strength generated by the GaN photoconductive emitter. Our findings show that terahertz power scales quadratically with optical excitation power and applied electric field over a broad parameter range. We achieved an optical-to-terahertz conversion efficiency approaching 100% within the 0.03–1 THz bandwidth at the highest bias field (116 kV/cm) in our experiment. Further optimization of GaN-based terahertz generation devices could achieve even greater optical-to-terahertz conversion efficiencies.

Uzundal, Can B↗
DOE OSTI2021

Improving reproducibility in synchrotron tomography using implementation-adapted filters

For reconstructing large tomographic datasets fast, filtered backprojection-type or Fourier-based algorithms are still the method of choice, as they have been for decades. These robust and computationally efficient algorithms have been integrated in a broad range of software packages. The continuous mathematical formulas used for image reconstruction in such algorithms are unambiguous. However, variations in discretization and interpolation result in quantitative differences between reconstructed images, and corresponding segmentations, obtained from different software. This hinders reproducibility of experimental results, making it difficult to ensure that results and conclusions from experiments can be reproduced at different facilities or using different software. In this paper, a way to reduce such differences by optimizing the filter used in analytical algorithms is proposed. These filters can be computed using a wrapper routine around a black-box implementation of a reconstruction algorithm, and lead to quantitatively similar reconstructions. Use cases for this approach are demonstrated by computing implementation-adapted filters for several open-source implementations and applying them to simulated phantoms and real-world data acquired at the synchrotron. Our contribution to a reproducible reconstruction step forms a building block towards a fully reproducible synchrotron tomography data processing pipeline.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Voltage positioning using co-optimization of controllable grid assets in radial networks

With increasing penetration of solar PV, some distribution feeders are experiencing highly variable net-load flows and even reverse flows. To optimize distribution systems under such conditions, the scheduling of mechanical devices, such as OLTCs and capacitor banks, needs to take into account forecasted solar PV and actual grid conditions. However, these legacy switching assets are operated on a daily or hourly timescale, due to the wear and tear associated with mechanical switching, which makes them unsuitable for real-time control. Therefore, there is a natural timescale- separation between these slower mechanical assets and the responsive nature of inverter-based resources. In this paper, we present a network admissible convex formulation for holistically scheduling controllable grid assets to position voltage optimally against solar PV. An optimal hourly schedule is presented that utilizes mechanical resources to position the predicted voltages close to nominal values, while minimizing the use of inverter-based resources (i.e., DERs), making them available for control at a faster time-scale (after the uncertainty reveals itself). A convex, inner approximation of the OPF problem is adapted to a mixed-integer linear program that minimizes voltage deviations from nominal (i.e., maximizes voltage margins). Here, the resulting OPF solution respects all the network constraints and is, hence, robust against modeling simplifications. Simulation based analysis on IEEE distribution feeders validates the approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Uncovering the basis of protein-protein interaction specificity with a combinatorially complete library

Protein-protein interaction specificity is often encoded at the primary sequence level. However, the contributions of individual residues to specificity are usually poorly understood and often obscured by mutational robustness, sequence degeneracy, and epistasis. Using bacterial toxin-antitoxin systems as a model, we screened a combinatorially complete library of antitoxin variants at three key positions against two toxins. This library enabled us to measure the effect of individual substitutions on specificity in hundreds of genetic backgrounds. These distributions allow inferences about the general nature of interface residues in promoting specificity. We find that positive and negative contributions to specificity are neither inherently coupled nor mutually exclusive. Further, a wild-type antitoxin appears optimized for specificity as no substitutions improve discrimination between cognate and non-cognate partners. By comparing crystal structures of paralogous complexes, we provide a rationale for our observations. Collectively, this work provides a generalizable approach to understanding the logic of molecular recognition.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Quantification of trace iodine using laser-induced breakdown spectroscopy for real-time monitoring of nuclear off-gas streams

This study evaluated the potential of laser-induced breakdown spectroscopy (LIBS) for real-time monitoring of trace gas-phase iodine, which is an element of high significance in nuclear applications due to its long radioactive half-life (as iodine-129), volatility, and biological impact. In anticipation of iodine evolving into off-gas systems in molten salt reactor and nuclear fuel recycling applications, this research aimed to assess LIBS performance in flowing argon and helium matrices; optimize measurement parameters using a multichannel spectrometer; and perform calibrations to assess predictive capabilities and limits of detection (LODs). Experimental results successfully measured gas-phase iodine in flowing argon and helium; however, trace iodine was not detected in air. Optimal delay times were determined to be 10 µs for argon and 1 µs for helium, which are consistent with the expected shorter plasma lifetime in helium relative to argon. An emission line survey was provided with the 206.16, 804.37, 902.24, and 905.83 nm peaks, which were identified as the strongest emission peaks. Calibration models were successfully built in both helium and argon, achieving LODs down to 3 ppm in helium and 5 ppm in argon. The iodine emission at 905.83 nm emerged as the most robust for calibration and was subsequently applied to a time series dataset in argon. The predictive trace confirmed the feasibility of employing LIBS for continuous, online quantification of trace iodine in flowing gas systems.

Andrews, Hunter B. [Oak Ridge National Laboratory ↗
DOE OSTI2020

The role of iron in magnetic damping of Mg(Al,Fe) 2 O 4 spinel ferrite thin films

© 2020 Author(s). We have investigated magnesium aluminum ferrite thin films with a range of iron concentrations and identified the optimal iron content to obtain high crystalline quality thin films with the low magnetic damping required for spin current-based applications. Epitaxial MgAl 2-x FexO4 films with 0.8 < x < 2.0 were grown by pulsed laser deposition on single crystal MgAl2O4 substrates and were characterized structurally and magnetically. We find that the x = 1.5 composition minimizes the room-temperature magnetic damping with a typical Gilbert damping parameter of α eff = 1.8 × 10-3. This minimized damping is governed by a competition between the more robust magnetic ordering with increased iron content, x, and the more defective structure due to larger film-substrate lattice mismatch with increased iron content. The temperature-dependent magnetization curves indicate that Tc is suppressed below room temperature for iron content x ≤ 1.2 and eventually suppressed entirely for x = 0.8. X-ray magnetic circular dichroism results indicate that for all x the magnetic moment is dominated by Fe 3 + cations distributed in a 60:40 ratio on the octahedral and tetrahedral sites, with minimal contribution from Fe 2 + cations. Films with x = 1.4-1.6 exhibit very strong ferromagnetic resonance and low Gilbert damping with α eff = (1.8-6) × 10-3, making them ideal candidates for microwave and spintronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Chemical substitution induced half-metallicity in CrMnSb (1– x ) P x

Half-metallic Heusler alloys have been intensively studied in recent years due to their potential applications in spin-based devices, e.g., in magnetic tunnel junctions. Yet, their properties may be very sensitive to the choice of the substrates, i.e., to the epitaxial strain and interface properties. Here, we report the results of our computational work on the half-Heusler compound CrMnSb (1–x) P x . In particular, we demonstrate that the parent compound CrMnSb is close to a half-metallic material at the optimized lattice parameter, with the onset of the half-metallic bandgap a few meV above the Fermi energy. Moreover, although it undergoes a half-metallic transition under a uniform compression of ~1.5%, such a transition is absent under epitaxial strain. At the same time, we show that a half-metallic transition could be induced by a chemical substitution of Sb with P, which results in a volume reduction of the unit cell. In particular, 50% substitution of Sb with P leads to a robust half-metallicity in CrMnSb (1–x) P x , with 100% spin polarization being retained at a large range of epitaxial strain. Furthermore, our results indicate that CrMnSb0.5P0.5 could be grown on different types of substrates, e.g., GaAs, without its electronic properties being detrimentally affected by biaxial strain. In addition, CrMnSb 0.5 P 0.5 exhibits a fully compensated ferrimagnetic alignment, which could be potentially useful in applications where stray magnetic fields are undesirable.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid–liquid mixtures

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

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

Minimizing the impacts of the ammonia economy on the nitrogen cycle and climate

Ammonia (NH 3 ) is an attractive low-carbon fuel and hydrogen carrier. However, losses and inefficiencies across the value chain could result in reactive nitrogen emissions (NH 3 , NO x , and N 2 O), negatively impacting air quality, the environment, human health, and climate. A relatively robust ammonia economy (30 EJ/y) could perturb the global nitrogen cycle by up to 65 Mt/y with a 5% nitrogen loss rate, equivalent to 50% of the current global perturbation caused by fertilizers. Moreover, the emission rate of nitrous oxide (N 2 O), a potent greenhouse gas and ozone-depleting molecule, determines whether ammonia combustion has a greenhouse footprint comparable to renewable energy sources or higher than coal (100 to 1,400 gCO 2 e/kWh). The success of the ammonia economy hence hinges on adopting optimal practices and technologies that minimize reactive nitrogen emissions. We discuss how this constraint should be included in the ongoing broad engineering research to reduce environmental concerns and prevent the lock-in of high-leakage practices.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Covariance Shaping Over Riemannian Manifolds for Massive MIMO Communication

Acquiring accurate instantaneous channel state information (CSI) is a challenging aspect of massive multi-input multi-output (MIMO) communication. Utilizing statistical information, such as channel covariance matrix, to design statistical beamforming vectors is robust when compared to instantaneous CSI. In this paper, we propose a novel MIMO covariance shaping scheme over Riemannian manifolds. It serves as an effective statistical beamforming solution to a number of close proximity user equipment (UE) that are undergoing substantial channel correlation. Proposed algorithm exploits the Hermitian positive definite nature of covariance matrices lying over Riemannian manifold. We introduce Wasserstein distance function as a Riemannian metric to measure distances between channel covariance matrices. Furthermore, K-means clustering technique is utilized to effectively identify the optimal shape of effective optimal covariance matrices. Our findings suggest that maximizing the geodesic distance between covariance matrices ultimately leads to a corresponding increase in the network throughput, as determined by the beamforming vector used to shape the covariance matrices. Simulation results validate that the proposed solution converges faster than Euclidean-based state-of-the-art, while maintaining the same computational complexity. Finally, the sum rate performance asymptotically achieves full capacity for two-UE case and more than 96% of the upper bound exhaustive search benchmark for multi-UE scenario.

42 ENGINEERING↗
DOE OSTI2020

Distributed Intrusion Detection System using Semantic-based Rules for SCADA in Smart Grid

Cyber-physical system (CPS) security for the smart grid enables secure communication for the SCADA and wide-area measurement system data. Power utilities world-wide use various SCADA protocols, namely DNP3, Modbus, and IEC 61850, for the data exchanges across substation field devices, remote terminal units (RTUs), and control center applications. Adversaries may exploit compromised SCADA protocols for the reconnaissance, data exfiltration, vulnerability assessment, and injection of stealthy cyberattacks to affect power system operation. In this paper, we propose an efficient algorithm to generate robust rule sets. We integrate the rule sets into an intrusion detection system (IDS), which continuously monitors the DNP3 data traffic at a substation network and detects intrusions and anomalies in real-time. To enable CPS-aware wide-area situational awareness, we integrated the methodology into an open-source distributed-IDS (D-IDS) framework. The D-IDS facilitates central monitoring of the detected anomalies from the geographically distributed substations and to the control center. The proposed algorithm provides an optimal solution to detect network intrusions and abnormal behavior. Different types of IDS rules based on packet payload, packet flow, and time threshold are generated. Further, IDS testing and evaluation is performed with a set of rules in different sequences. The detection time is measured for different IDS rules, and the results are plotted. All the experiments are conducted at Power Cyber Lab, Iowa State University, for multiple power grid models. After successful testing and evaluation, knowledge and implementation are transferred to field deployment.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Mitigating the noise of DESI mocks using analytic control variates

In order to address fundamental questions related to the expansion history of the Universe and its primordial nature with the next generation of galaxy experiments, we need to model reliably large-scale structure observables such as the correlation function and the power spectrum. Cosmological N-body simulations provide a reference through which we can test our models, but their output suffers from sample variance on large scales. Fortunately, this is the regime where accurate analytic approximations exist. To reduce the variance, which is key to making optimal use of these simulations, we can leverage the accuracy and precision of such analytic descriptions using Control Variates (CV). The power of control variates stems from utilizing inexpensive but highly correlated surrogates of the statistics one wishes to measure. The stronger the correlation between the surrogate and the statistic of interest, the larger the variance reduction delivered by the method. We apply two control variate formulations to mock catalogs generated in anticipation of upcoming data from the Dark Energy Spectroscopic Instrument (DESI) to test the robustness of its analysis pipeline. Our CV-reduced measurements offer a factor of 5-10 improvement in the measurement error compared with the raw measurements. We explore the relevant properties of the galaxy samples that dictate this reduction and comment on the improvements we find on some of the derived quantities relevant to Baryon Acoustic Oscillation (BAO) analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A simple method for obtaining heat capacity coefficients of minerals

Abstract Heat capacity data are unavailable or incomplete for many minerals at geologically relevant temperatures. Despite the availability of entropy and enthalpy values in numerous thermodynamic tables (even sometimes at elevated temperatures), there remains need for extrapolation beyond, or interpolation between, temperatures. This approach inevitably results in estimates for entropy and enthalpy values because the heat capacity coefficients required for optimal thermodynamic treatment are less frequently available. Here we propose a simple method for obtaining heat capacity coefficients of minerals. This method requires only the empirically measured temperature-specific heat capacity for calculation via a matrix algorithm. The system of equations solver is written in the Python computing language and has been made accessible in an online repository. Thermodynamically, the solution to a system of equations represents the heat capacity coefficients that satisfy the mineral-specific polynomial. Direct coefficient calculation will result in more robust thermodynamic data, which are not subject to fitting uncertainties. Using hematite as an example, this method provides results that are comparable to conventional means and is applicable to any solid material. Coefficients vary within the traditional large 950 K temperature interval, indicating that best results should instead utilize a smaller 400 K temperature interval. Examples of large-scale implications include the refinement of geothermal gradient estimation in rapidly subsiding sedimentary basins or metamorphic and hydrothermal evolution.

Geochemistry & Geophysics↗
DOE OSTI2021

Material Identification Using Dual Energy X-ray Absorptiometry

Two implementations of dual energy X-ray absorptiometry were studied to identify materials using X-ray attenuation data taken with the Digital Radiography and Computed Tomography (DRCT) systems that were developed for the Recovered Chemical Materiel Directorate (RCMD). Maitrejean et al.’s approach utilizes eigen effects through Principal Component Analysis, while Osipov et al.’s approach proposed a physics-based method. Both approaches approximate mass attenuation coefficients of materials as a linear combination of basis functions (eigen effects) or physics-based equations. A set of coefficients {a 1 , a 2 , a 3 } or {B, D} were found by parameter optimization in EXCEL Solver. The identification parameters, {$\frac{a_{2}}{a_{1}}$, $\frac{a_{3}}{a_{1}}$} or estimated effective atomic number $\hat{Z}$ from {B, D}, were calculated to identify material of an aluminum 8 step wedge and a steel 8 step wedge in X-ray radiography images taken by a DRCT system. Maitrejean et al.’s approach was unable to provide reliable $\frac{a_{3}}{a_{1}}$ ratio values for identification of materials. Osipov et al.’s approach was found to be more robust in identify materials with a semi-empirical formula derived from test results in this study.

36 MATERIALS SCIENCE↗