Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational 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.

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

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Counterdiabatic Driving with Performance Guarantees

Counterdiabatic (CD) driving has the potential to speed up adiabatic quantum state preparation by suppressing unwanted excitations. However, existing approaches either require intractable classical computations or are based on approximations that do not have performance guarantees. We propose and analyze a nonvariational, system-agnostic CD expansion method and analytically show that it converges exponentially quickly in the expansion order. In finite systems, the required resources scale inversely with the spectral gap, which we argue is asymptotically optimal. To extend our method to the thermodynamic limit and suppress errors stemming from high-frequency transitions, we leverage finite-time adiabatic protocols. In particular, we show that a time determined by the quantum speed limit is sufficient to prepare the desired ground state, without the need to optimize the adiabatic trajectory. Numerical tests of our method on the quantum Ising chain show that our method can outperform state-of-the-art variational CD approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A simple method for mitigating error in the fixed-offset of a double-crystal monochromator

In order to maintain a constant monochromatic synchrotron X-ray beam height for all energies, the separation between the crystals of a double-crystal monochromator is typically adjusted, via translation of the second crystal, while X-ray energy is varied, via rotation of the crystal pair. The ability to accurately translate the second crystal requires precise knowledge of the separation between the two crystals and, when present, crystal miscuts. Here, a simple method for calibrating the crystal gap from measured variation in the X-ray beam height that eliminates error in the fixed beam offset is provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Using distributed power electronics-based devices to improve the voltage and frequency stability of distribution systems

Methods include, in response to a line frequency variation of a power grid, adjusting a voltage setpoint of a voltage regulator coupled to the power grid at a grid edge to maintain a voltage at the grid edge, wherein the adjusting the regulated voltage setpoint is configured to reduce the line frequency variation to stabilize the line frequency of the power grid. Apparatus include a voltage regulator configured to couple to a power grid at a grid edge and to maintain a voltage at the grid edge, wherein the voltage regulator is further configured to adjust a voltage setpoint of a voltage regulator in response to a line frequency variation of the power grid to reduce the line frequency variation and stabilize the line frequency of the power grid.

Du, Wei↗
DOE OSTI2023

Economies of Numbers Formulations for Optimal Process Family Design of Carbon Capture Systems

Enabling rapid, widespread deployment of many process variants across a range of design requirements is key to ensuring the widespread availability of critical decentralized technologies such as carbon capture systems. Unfortunately, traditional process design methods are not efficient for designing many similar variations of a process rapidly across a wide range of system requirements. Rather than considering each process variant as a separate design task, we have proposed viewing the entire set of variants as a “family.” In our work, we present rigorous optimization formulations that simultaneously design each variant within the family from a platform of available unit modules, largely inspired by product family and platform design.

Stinchfield, Georgia↗
DOE OSTI2024

Accelerating Application Bulk Synchronous Writes in HPC Environments

High-bandwidth storage tiers are becoming more common for their capability to absorb high-rate, bursty I/Os. Notably, the designs of these fast storage tiers differ from system to system. The variation of these layers and non-uniform methods of access can pose chal- lenges for applications seeking to run at multiple HPC facilities. Therefore, in this work, we present Spectral, a rapid-output ab- straction library to accelerate application, bulk-synchronous writes on HPC systems. We design Spectral to enable applications to use high-bandwidth storage, such as node-local storage and dis- tributed, write-caches (e.g., burst buffers) transparently without requiring modifications to the application or file system source code. The key idea is to allow applications to spend most of the time performing productive work and to not require any source code changes for maximum portability on different HPC archi- tectures. Spectral internally re-routes write-only files through available, high-performance I/O resources before ultimately mi- grating them to the shared global parallel file system. For instance, on Summit, Spectral transparently places application outputs on node-local storage and then utilizes asynchronous migration to the center-wide GPFS file system. We evaluate Spectral on the Summit HPC system (1024 nodes) using the IOR benchmark and real scientific applications. Spectral shows linear performance scaling, improving application write performance by over an order of magnitude when compared to GPFS.

Khan, Awais↗
DOE OSTI2026

Renormalon subtracted nonrelativistic QCD for heavy hadron systems

We present a renormalon-subtracted formulation of potential nonrelativistic QCD (pNRQCD) for precision spectroscopy of heavy hadron systems, combining variational and Green's function Monte Carlo (VMC/GFMC) methods with NNLO static two- and three-body potentials. Minimal renormalon subtraction (MRS) systematically sums leading factorially growing terms, thereby stabilizing perturbative convergence and reducing renormalization-scale dependence. We tune charm and bottom quark masses to spin-averaged $1S$ quarkonium states and predict $Ω_{ccc}$, $Ω_{ccb}$, $Ω_{cbb}$, and $Ω_{bbb}$ baryon masses, as well as QCD-stable baryons containing top quarks. NNLO MRS results undershoot lattice QCD by 125--175~MeV, with fractional differences decreasing as $\sim 1/m_Q$, consistent with neglected $\mathcal{O}(1/m_Q)$ corrections. Applying these methods to unequal-mass fully-heavy tetraquarks, we determine the critical heavy-to-light mass ratio for binding and compute binding energies across the mass-ratio landscape.

Assi, Benoît [Cincinnati U.; Fermilab]↗
DOE OSTI2021

Adaptive variational quantum eigensolvers for highly excited states

Highly excited states of quantum many-body systems are central objects in the study of quantum dynamics and thermalization that challenge classical computational methods due to their volume-law entanglement content. In this work, we explore the potential of variational quantum algorithms to approximate such states. We propose an adaptive variational quantum eigensolver (VQE) for excited states (X) that self-generates a variational ansatz for arbitrary eigenstates of a many-body Hamiltonian H by attempting to minimize the energy variance with respect to H. We benchmark the method by applying it to an Ising spin chain with integrable and nonintegrable regimes, where we calculate various quantities of interest, including the total energy, magnetization density, and entanglement entropy. We also compare the performance of adaptive VQE-X to an adaptive variant of the folded-spectrum method. For both methods, we find a strong dependence of the algorithm's performance on the choice of operator pool used for the adaptive construction of the ansatz. In particular, an operator pool including long-range two-body gates accelerates the convergence of both algorithms in the nonintegrable regime. Here, we also study the scaling of the number of variational parameters with system size, finding that an exponentially large number of parameters may be necessary to approximate individual highly excited states. Nevertheless, we argue that these methods lay a foundation for the use of quantum algorithms to study finite-energy-density properties of many-body systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Variational AutoEncoders Reveal Intensifying GPP Extremes in Continental United States based on CESM2 Simulations

Climate extremes significantly impact terrestrial carbon cycle dynamics, necessitating robust methods for detecting and analyzing anomalous behavior in plant productivity. This study presents a novel application of variational autoencoders (VAE) for identifying extreme events in gross primary productivity (GPP) from Community Earth System Model version 2 simulations across four AR6 regions in the Continental United States. We compare VAE-based anomaly detection with traditional singular spectral analysis (SSA) methods across three time periods: 1850-80, 1950-80, and 2050-80 under SSP5-8.5 scenario. The VAE architecture employs three dense layers and a latent space with input sequence length of 12 months, training on normalized GPP time series to reconstruct the GPP and identify anomalies based on reconstruction errors. Extreme events are defined using 5th percentile thresholds applied to both VAE and SSA anomalies. Results demonstrate strong regional agreement between VAE and SSA methods in spatial patterns of extreme event frequencies, despite VAE consistently producing higher threshold values (179-756 GgC for VAE vs. 100-784 GgC for SSA across regions and periods). Both methods reveal increasing magnitudes and frequencies of negative carbon cycle extremes toward 2050-80, particularly in Western and Central North America. The VAE approach shows comparable performance to established SSA techniques while offering computational advantages and enhanced capability for capturing non-linear temporal dependencies in carbon cycle variability. This research demonstrates the potential of deep learning approaches for extremes detection and provides a foundation for improved understanding of future carbon cycle risks under future conditions.

Sharma, Bharat [ORNL] (ORCID:0000000266982487)↗
DOE OSTI2021

Unpacking the drivers of diurnal dynamics of sun-induced chlorophyll fluorescence (SIF): Canopy structure, plant physiology, instrument configuration and retrieval methods

Sun-induced chlorophyll fluorescence (SIF) from spaceborne sensors is a promising tool for global carbon cycle monitoring, but its application is constrained by insufficient understanding of the drivers underlying diurnal SIF dynamics. SIF measurements from ground-based towers can reveal diurnal SIF dynamics across biomes and environmental conditions; however, meaningful interpretation of diurnal variations requires disentangling impacts from canopy structure, plant physiology, instrument configuration and retrieval methods, which often interact with and confound each other. This study aims to unpack these drivers using 1) concurrent ground and airborne canopy-scale and leaf-scale measurements at a corn field, 2) a mechanistic SIF model that explicitly considers the dynamics of photochemistry (via the fraction of open photosystem II reaction centers, qL) and photoprotection (via nonphotochemical quenching, NPQ) as well as their interactive dependence on the sub-canopy light environment, and 3) cross-comparison of SIF instrument configurations and retrieval methods. We found that crop row orientations and sun angles can introduce a distinctive midday dip in SIF in absence of stress, due to a midday drop of absorbed photosynthetically active radiation (APAR) when crop rows are north-south oriented. Canopy structure caused distinctive responses in both qL and NPQ at different positions within the vertical canopy that collectively influenced fluorescence quantum yield (Φ F ) at the leaf scale. Once integrated at the canopy scale, diurnal dynamics of both APAR and canopy escape probability (ε) are critical for accurately shaping diurnal SIF variations. While leaf-level qL and NPQ exhibited strong diurnal dynamics, their influence was attenuated at the canopy scale due to opposing effects on SIF at different canopy layers. Furthermore, different system configurations (i.e., bi-hemispherical vs. hemispherical-conical) and retrieval methods can bias the SIF magnitude and distort its diurnal shapes, therefore confounding the interpretation of inherent strength and dynamics of SIF emission. Our findings demonstrate the importance of crop row structures, interactive variations in canopy structure and plant physiology, instrument configuration, and retrieval method in shaping the measured dynamics of diurnal SIF. This study highlights the necessity to account for these factors to accurately interpret satellite SIF, and informs future synthesis work with different SIF instrumentation and retrieval methods across sites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A unified framework of stabilized finite element method for solving the Boltzmann transport equation

This paper presents a unified framework of stabilized finite element method for solving the Boltzmann transport equation. Unlike the traditional Petrov-Galerkin finite element method which modifies the test function to construct the stabilization term, we derive the stabilization methods from the standard Galerkin weak form with Sub-grid scale model. The basic idea of this method is to decompose the unknowns into its numerical solution and residual, with an approximation for the residual and embeds it in the Galerkin weak form to yield a stabilized variational formula. Different approximations of the residual lead to different stabilization methods, all the frequently used stabilized methods, including the Streamline Upwinding Petrov-Galerkin (SUPG) method, Galerkin/Least-Square (GLS) method, and Algebraic Sub-Grid Scale (ASGS) method can be obtained from this framework. The similarities and differences of the different approximations are compared in this paper. The numerical results show that the behaviors of the different methods area similar with the same stabilization parameter, and all these stabilized techniques can obtain a correct and stable solution. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Estimators and Fusers for Fiber Delay Estimation Using Environmental Measurements

The properties of deployed network fiber are affected by environmental factors due to their exposure to the elements. Particularly for quantum networks, the resultant delay variations may have significant impacts due to the extreme sensitivity of synchronization, coincidence counting, and other critical operations. In this paper, the delays of 15 km aerial-inground fiber connections are measured, and effects due to temperature, humidity and wind speed are analyzed over multiple periods spanning four seasons of a year. Machine learning methods are first utilized to reveal surprisingly pronounced effects of humidity on the delay, in addition to the expected temperature and its seasonal variations. Estimator and fusion methods are developed to estimate the delay using temperature, humidity and wind speed measurements, by utilizing smooth Gaussian Process Regression (GPR) and nonsmooth Ensemble of Trees (EOT) methods. Measurements from winter and summer periods are temporally fused using twelve different methods, and eight methods provide estimates for the delay throughout the year with median test errors under 1.28%. The results reveal distinct temperature-humidity trends across the seasons, and the ability of estimator and temporal fusion methods to exploit them for estimating the delay. These results constitute a case study of machine learning analytical results, wherein generalization equations explain the performance of various estimator and fuser methods.

Rao, Nageswara [ORNL] (ORCID:0000000234085941)↗
DOE OSTI2022

The Effect of Cell Geometry and Trigger Method on the Risks Associated with Thermal Runaway of Lithium-Ion Batteries

Consideration of thermal runaway heat output variability is paramount for the development of safe lithium-ion battery assemblies. This study utilizes data gathered from fractional thermal runaway calorimetry (FTRC) experiments to conduct a comparative analysis of thermal runaway heat output for three cell formats (18650, 21700, and 33600) as a function of trigger method (heaters, internal short-circuiting device, and nail penetration). The analysis is based on comparisons for the calculated total energy yield, fractional energy yield, heat rate, and heat flux. This study reveals that nail penetration tends to result in higher thermal runaway heat output for larger cells (21700 & 33600); these experiments also tended to result in higher fractions of the total energy being released through the cell body. The smaller cells (18650) did not appear to have significant variation in heat output as a function of trigger method. This finding suggests that, for this cell type, worst-case scenario heat output could be achievable in assembly level testing regardless of the utilized trigger method. This study also demonstrates successful translation of FTRC results, as recorded in the Battery Failure Databank, into meaningful analysis that breaks down the influence of specific conditions on thermal runaway heat output.

calorimetry↗
DOE OSTI2022

General-Purpose Bayesian Tensor Learning With Automatic Rank Determination and Uncertainty Quantification

A major challenge in many machine learning tasks is that the model expressive power depends on model size. Low-rank tensor methods are an efficient tool for handling the curse of dimensionality in many large-scale machine learning models. The major challenges in training a tensor learning model include how to process the high-volume data, how to determine the tensor rank automatically, and how to estimate the uncertainty of the results. While existing tensor learning focuses on a specific task, this paper proposes a generic Bayesian framework that can be employed to solve a broad class of tensor learning problems such as tensor completion, tensor regression, and tensorized neural networks. We develop a low-rank tensor prior for automatic rank determination in nonlinear problems. Our method is implemented with both stochastic gradient Hamiltonian Monte Carlo (SGHMC) and Stein Variational Gradient Descent (SVGD). We compare the automatic rank determination and uncertainty quantification of these two solvers. We demonstrate that our proposed method can determine the tensor rank automatically and can quantify the uncertainty of the obtained results. We validate our framework on tensor completion tasks and tensorized neural network training tasks.

Bayesian inference↗
DOE OSTI2025

The electronic structure, crystal fields, and magnetic anisotropy in RECo 5 magnets

The current progress in describing rare-earth-based magnets' electronic structure and magnetic properties is discussed. We use several currently popular electronic structure methods to show the typical values of critical parameters that define the physics of RECo 5 (RE = rare earth atom) materials. The magnetic moments and magnetic anisotropy of 4f atoms are obtained using several approaches, including anisotropic 4f-charge density-constrained DFT and DFT+HI methods. We also suggest the introduction of "penalty" functional for obtaining correct variational total energy in the traditional Hund's rule-constrained DFT-based techniques. The applicability and future extensions are discussed. The proposed combination of methods is potentially suitable for high-throughput computational searches of new rare-earth-containing magnetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Time-lapse electrical resistivity and induced polarization monitoring of a simulated ecosystem-scale coastal flooding experiment

Accurate estimation of the spatio-temporal variation in soil saturation and salt concentration is valuable for calibrating and evaluating Earth system models as well as detecting important eco-hydrological interactions. Such representative models play a vital role in investigating the stability of ecosystems facing changing hydrologic disturbance regimes. Ecosystem-scale field experiments are a useful way to capture hydrological disturbance and gain a deeper understanding of the complex mechanisms driving environmental changes. However, previous studies have not explored the use of quantitative imaging at a large spatial scale to reproduce breakthrough curves and inform system response and recovery following hydrologic disturbance events. Therefore, in this study, non-invasive geophysical methods were used to capture the spatio-temporal variation in subsurface saturation and salt concentrations during an ecosystem-scale coastal forest flooding experiment.

coastal↗
DOE OSTI2025

A critical review of the history of fabricating monolithic U-Mo fuel plates

The fabrication of monolithic uranium-molybdenum alloy fuels, specifically those developed for high-performance research and test reactors, began in the early 2000s. The primary fuel form consists of uranium alloyed with 10-wt% molybdenum in a thin foil coated with zirconium and encapsulated in aluminum-6061 cladding. Over the years, the process has evolved with different types of casting, heat treatments, rolling schedules, and cladding applications. This review examines the history of these fabrication processes and their impact on microstructure and fuel-swelling performance. Even though various fabrication methods were used, we found little correlation between fabrication variation and fuel swelling. This insensitive relationship between fabrication variation and fuel swelling is primarily due to inhomogeneous microstructures that formed during casting and grain refinement that occurred during rolling. We conclude that the fabrication processes we examined produced similar microstructures, indicating that the fuel microstructure is somewhat insensitive to the fabrication parameters evaluated. However, the relatively small amount of historic data, such as those for grain sizes, limited this analysis. More recently fabricated materials, such as those from ongoing irradiation experiment, Mini-Plate-1 and Mini-Plate-2, were also excluded from this analysis and are intended to be reviewed separately. The findings, that fuel microstructure is somewhat insensitive to the fabrication parameters, do not imply that any fabrication method is acceptable, given the uncertainties in data and fuel-swelling observations. For example, only arc melting and vacuum induction melting casting processes were previously explored in the historic fabrication efforts. In conclusion, the findings should not be extrapolated to other casting processes.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

An Electrical Resistance Diagnostic for Conductivity Monitoring in Laser Powder Bed Fusion

With the growing interest in metal additive manufacturing using laser powder bed fusion (LPBF), there is a need for advanced in-situ nondestructive evaluation (NDE) methods that can dynamically monitor manufacturing process-related variations, that can be used as a feedback mechanism to further improve the manufacturing process, leading to parts with improved microstructural properties and mechanical properties. Current NDE techniques either lack sensitivity beyond build layer, are costly or time-consuming, or are not compatible for in-situ integration. In this research, we develop an electrical resistance diagnostic for in-situ monitoring of powder fused regions during laser powder bed fusion printing. The technique relies on injecting current into the build plate and detecting voltage differences from conductive variations during printing using a simple, cheap four-point electrode array directly connected to the build plate. A computational model will be utilized to determine sensitivities of the approach, and preliminary experiments will be performed during the printing process to test the overall approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Prime factorization using quantum variational imaginary time evolution

The road to computing on quantum devices has been accelerated by the promises that come from using Shor’s algorithm to reduce the complexity of prime factorization. However, this promise hast not yet been realized due to noisy qubits and lack of robust error correction schemes. Here we explore a promising, alternative method for prime factorization that uses well-established techniques from variational imaginary time evolution. We create a Hamiltonian whose ground state encodes the solution to the problem and use variational techniques to evolve a state iteratively towards these prime factors. We show that the number of circuits evaluated in each iteration scales as \(O(n^{5}d)\) , where n is the bit-length of the number to be factorized and d is the depth of the circuit. We use a single layer of entangling gates to factorize 36 numbers represented using 7, 8, and 9-qubit Hamiltonians. We also verify the method’s performance by implementing it on the IBMQ Lima hardware to factorize 55, 65, 77 and 91 which are greater than the largest number (21) to have been factorized on IBMQ hardware.

97 MATHEMATICS AND COMPUTING↗