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.

SourceResearch connectionFrom yearThrough year

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

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Gain balanced nonlinear optical interferometer

A nonlinear fiber interferometer is disclosed suitable for fiber sensor and other applications. A first nonlinear fiber section amplifies probe and conjugate sidebands of a pump through four-wave mixing. A second section introduces a phase shift to be measured, for example from a sensor. A third nonlinear fiber section amplifies with phase-sensitive gain to increase signal-to-noise ratio. Based on phase-sensitive output power of probe and/or conjugate components, the phase shift can be measured. Superior performance can be obtained by balancing gain between the (first and third) nonlinear sections. Non-fiber, for example photonic integrated circuit, embodiments are disclosed. Differential sensing, alternative detection schemes, sensing applications, associated methods, and other variations are disclosed.

Lukens, Joseph M.↗

DOE OSTI2023

Gain balanced nonlinear optical interferometer

Lukens, Joseph M.↗

DOE OSTI2025

Nonunitary Variational Quantum Eigensolver with the Localized Active Space Method and Cost Mitigation

Accurately describing strongly correlated systems with affordable quantum resources remains a central challenge for quantum chemistry applications on near and intermediate term quantum computers. The localized active space self-consistent field (LASSCF) approximates the complete active space self-consistent field (CASSCF) by generating active space-based wave functions within specific fragments while treating interfragment correlation with mean-field approach, hence is computationally less expensive. Hardware-efficient ansatzes (HEA) offer affordable and shallower circuits, yet they often fail to capture the necessary correlation. Previously, Jastrow-factor-inspired nonunitary qubit operators were proposed to use with HEA for variational quantum eigensolver (VQE) calculations (so-called nuVQE), as they do not increase circuit depths and recover correlation beyond the mean-field level for Hartree–Fock initial states. Here, in this study, we explore running nuVQE with LASSCF as the initial state. The method, named LAS-nuVQE, is shown to recover interfragment correlations, reach chemical accuracy with a small number of gates (<70) in both H 4 and square cyclobutadiene (C 4 H 4 ), and produces more accurate energetics than its HEA counterparts at all circuit depths. To further address the inherent symmetry-breaking in HEA, we implemented spin-constrained LAS-nuVQE to extend the capabilities of HEA further and show spin-pure results for square cyclobutadiene. We also mitigate the increased measurement overhead of nuVQE via Pauli grouping and shot-frugal sampling, reducing measurement costs by up to 2 orders of magnitude compared to ungrouped operator, and show that one can achieve better accuracy with a small number of shots (10 3–4 ) per one expectation value calculation compared to noiseless simulations with one or two orders of magnitude more shots. Finally, wall clock time estimates show that, with our measurement mitigation protocols, nuVQE becomes a cheaper and more accurate alternative than vanilla VQE with HEA. Taken together, these developments illustrate a practical pathway toward performing multireference chemical simulations with accuracy and affordable resources on today’s quantum hardware, achieving both accuracy and affordability in challenging correlated systems.

Wang, Qiaohong [Univ. of Chicago, IL (United State↗

DOE OSTI2025

Phenome‐to‐genome insights for evaluating root system architecture in field studies of maize

Abstract Understanding the genetic basis of root system architecture (RSA) in crops requires innovative approaches that enable both high‐throughput and precise phenotyping in field conditions. In this study, we evaluated multiple phenotyping and analytical frameworks for quantifying RSA in mature, field‐grown maize in three field experiments. We used forward and reverse genetic approaches to evaluate >1700 maize root crowns, including a diversity panel, a biparental mapping population, and maize mutant and wild‐type alleles at two known RSA genes,DEEPER ROOTING 1(DRO1) andRootless1(Rt1). We show the utility of increasing the dimensionality of traditional two‐dimensional (2D) techniques, referred to as the “2D multi‐view” method, to improve the capture of whole root system information for mapping genetic variation influencing RSA. Comparison of univariate and multivariate genome‐wide association study (GWAS) approaches revealed that multivariate traits were effective at dissecting complex RSA phenotypes and identifying pleiotropic quantitative trait loci (QTLs). Overall, three‐dimensional (3D) root models generated from X‐ray computed tomography and digital phenotyping captured a larger proportion of RSA trait variations compared to other methods of root phenotyping, as evidenced by both genome‐wide and single‐gene analyses. Among the individual root traits, root pulling force emerged as a highly heritable estimate of RSA that identified the largest number of shared QTLs with 3D phenotypes. Our study shows that integrating complementary phenotyping technologies helps to provide a more comprehensive understanding of the genetic architecture of RSA in field‐grown maize.

Genetics & Heredity↗

DOE OSTI2021

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Parametric Comparative Analysis between Virtual Synchronous Generator and Droop-based Inertia for Inverter-Based Microgrids

This paper presents a parametric comparative analysis between the virtual synchronous generator (VSG) method and the droop control method to emulate inertia in the voltage-source inverter (VSI). Droop controllers are commonly used to regulate sharing power in microgrids and distribute power generation proportionally among VSI’s depending on their rated power. Additionally, VSG has been used to regulate the Rate-of-Change-of-Frequency (RoCoF) of the microgrid using virtual inertia. Although both methods can be used to regulate the frequency variation, the influence of each method on the closed loop eigenvalues is not the same. In this work, the transient response of the frequency is analyzed for each method to determine their advantages and disadvantages regarding frequency regulation in microgrid applications. The results were verified by conducting experimental trials using VSI’s. These experiments demonstrated that VSG is more suitable for regulating RoCoF and frequency nadir than droop controllers since it provides inertial support and improves frequency response.

Campo-Ossa, Daniel D.↗

DOE OSTI2020

Integrating Predictions for Improving Defect Classification Accuracy in NDT-based Assessment of Concrete - 20229

There is an increasing need to create predictive models for defect classification in concrete using the output of non-destructive testing (NDT) techniques. Recent advancement of machine-learning algorithms has offered several techniques for developing classification models for different types of data sets. However, the performance of these algorithms is very uncertain, mainly when applied to small and noisy datasets. For example, when human access is limited (e.g., nuclear facility), robot-based NDT is preferred. But compared to manual tests with humans present on site, the data sets are small, and more noise can exist. Therefore, it is imperative to develop new approaches to ensure a consistently high classification accuracy for inadequate data sets. This study explores the classification performance on NDT dataset using classifiers from different machine-learning algorithms, namely k-Nearest Neighbor (kNN), Decision Tree, Naive Bayes, Logistic Regression, and Support Vector Machine (SVM). The authors further integrated the predictions from these classifiers using proposed methods. The integration strategy combines the output of the classifiers based on two different measures, accuracy, and performance (ACC and PERF), using equations such as sum, average, and square-root-of-sums-of-squares (SRSS). Our results reveal varying classification accuracies across individual classifiers with different misclassifications across the test data set. The integration strategy provided significant improvement in the classification accuracy compared to the individual classifiers. Furthermore, the results indicate minimal variation across the integration methods as compared to the variation across the individual classifiers. To conclude, prediction integration offers a unique approach for combining the output of multiple classifiers to create redundancies with the potential of achieving high classification performance and improved reliability in predictive models for defect detection in concrete. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Capturing the Relative‐Humidity‐Sensitive Gas–Particle Partitioning of Organic Aerosols in a 2D Volatility Basis Set

Abstract Aerosol water affects the physicochemical properties and mass concentration of organic aerosols (OA), but it is typically omitted by air quality, weather, and climate models. We compare two classes of simplified models to estimate the OA water uptake and gas–particle partitioning of organic compounds. One class uses a single‐hygroscopicity‐parameter (κ) approach while the other is based on the reduced‐complexity Binary Activity Thermodynamics (BAT) model. We show that a BAT‐based two‐dimensional volatility basis set (VBS) model always predicts a higher OA mass concentration at elevated relative humidity (RH), for example, ∼16% at 80% RH, than any variation of theκ‐based method considered—even when BAT‐VBS predicts a lower water uptake. The main reason being that the BAT‐VBS model captures variations in effective saturation mass concentration of organics (C*) with RH, a feature that other VBS methods lack. The BAT‐VBS framework offers an efficient, RH‐sensitive treatment for reduced‐complexity OA modeling.

Geology↗

DOE OSTI2020

Variational Localized Active Space Self-Consistent Field Method

Not provided.

Chemistry↗

DOE OSTI2024

Residual Stress Measurements of Cr-coated Cladding

Accident-tolerant fuel concepts for light-water reactor applications have been developed and tested in diverse research programs around the world. Industry teams have developed coated cladding concepts using different coating application methods, such as physical vapor deposition and coldspraying. Recently, variations in thickness and application method were investigated, and their effects on the microstructure and mechanical properties were reported. The results suggest an increase in performance in an unirradiated state, however, it remains unclear what phenomenon drove this performance. One explanation to this performance is associated with the coating process imparting a residual stress. Residual stress analysis was performed on coatings produced with various process parameters to understand how a thin layer of Cr impacts the mechanical properties after coating application and after simulated operating conditions. The gathered information is intended to address safety concerns or support margin identification related to accident-tolerant fuels and to support subsequent modeling and simulation efforts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Geometric electrostatic particle-in-cell algorithm on unstructured meshes

We present a geometric particle-in-cell (PIC) algorithm on unstructured meshes for studying electrostatic perturbations with frequency lower than electron gyrofrequency in magnetized plasmas. In this method, ions are treated as fully kinetic particles and electrons are described by the adiabatic response. The PIC method is derived from a discrete variational principle on unstructured meshes. To preserve the geometric structure of the system, the discrete variational principle requires that the electric field is interpolated using Whitney 1-forms, the charge is deposited using Whitney 0-forms and the electric field is computed by discrete exterior calculus. The algorithm has been applied to study the ion Bernstein wave (IBW) in two-dimensional magnetized plasmas. The simulated dispersion relations of the IBW in a rectangular region agree well with theoretical results. In a two-dimensional circular region with fixed boundary condition, the spectrum and eigenmode structures of the IBW are obtained from simulations. We compare the energy conservation property of the geometric PIC algorithm derived from the discrete variational principle with that of previous PIC methods on unstructured meshes. The comparison shows that the new PIC algorithm significantly improves the energy conservation property.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Reweighting Underlying Event and Colour Reconnection parameter variations in Sherpa

We propose and validate a new method to trace the impact of parameter variations in the simulation of multi-parton interactions and colour reconnections in the Sherpa event generator. They are reflected, at an event-by-event basis, through relative weights with respect to the central production parameters that give rise to the generated events and distributions. Our method facilitates the tuning of the Monte Carlo event generator at a dramatically reduced computational cost, alleviates parameter sensitivity studies, and enables robust quantification of parametric uncertainties on-the-fly, one of the missing ingredients for future simulations of high-energy particle collisions. The method can easily be adapted to and implemented in other event generators. To illustrate its potential, we here consider combined tunes of the multi-parton-interaction and colour-reconnection models in Sherpa using LHC proton-proton collision data at $\sqrt{s}=7\,\text{TeV}$. We furthermore calibrate the energy-scaling behaviour of dimensionful model parameters based on $\sqrt{s}=13\,\text{TeV}$ LHC data and Tevatron data taken at $\sqrt{s}=1.96\,\text{TeV}$.

Pabst, Moritz [Gottingen U.]↗

DOE OSTI2023

Simulating scalar field theories on quantum computers with limited resources

Here we present a quantum algorithm for implementing φ 4 lattice scalar field theory on qubit computers. The field is represented in the discretized field amplitude basis. The number of qubits and elementary gates required by the implementation of the evolution operator is proportional to the lattice size. The algorithm allows efficient φ 4 state preparation for a large range of input parameters in both the normal and broken-symmetry phases. The states are prepared using a combination of variational and adiabatic evolution methods. First, the ground state of a local Hamiltonian, which includes the φ 4 self-interaction, is prepared using short variational circuits. Next, this state is evolved by switching on the coupling between the lattice sites adiabatically. The parameters defining the local Hamiltonian are adjustable and constitute the input of our algorithm. We present a method to optimize these parameters in order to reduce the adiabatic time required for state preparation. For preparing broken-symmetry states, the adiabatic evolution problems caused by crossing the phase transition critical line and by the degeneracy of the broken-symmetry ground state can be addressed using an auxiliary external field which gradually turns off during the adiabatic process. We show that the time dependence of the external field during the adiabatic evolution is important for addressing the broken-symmetry ground state degeneracy. The adiabatic time dependence on the inverse error tolerance can be reduced from quadratic to linear by using a field strength that decreases exponentially in time relative to one that decreases linearly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Synchronous and Concurrent Multidomain Computing Method for Cloud Computing Platforms

We present a numerical method for synchronous and concurrent solution of transient elastodynamics problem where the computational domain is divided into subdomains that may reside on separate computational platforms. Here, this work employs the variational multiscale discontinuous Galerkin (VMDG) method to develop interdomain transmission conditions for transient problems. The fine-scale modeling concept leads to variationally consistent coupling terms at the common interfaces. The method admits a large class of time discretization schemes, and decoupling of the solution for each subdomain is achieved by selecting any explicit algorithm. Numerical tests with a manufactured solution problem show optimal convergence rates. The energy history in a free vibration problem is in agreement with that of the solution from a monolithic computational domain.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Deeply uncertain: comparing methods of uncertainty quantification in deep learning algorithms

We present a comparison of methods for uncertainty quantification (UQ) in deep learning algorithms in the context of a simple physical system. Three of the most common uncertainty quantification methods - Bayesian Neural Networks (BNN), Concrete Dropout (CD), and Deep Ensembles (DE) - are compared to the standard analytic error propagation. We discuss this comparison in terms endemic to both machine learning ("epistemic" and "aleatoric") and the physical sciences ("statistical" and "systematic"). The comparisons are presented in terms of simulated experimental measurements of a single pendulum - a prototypical physical system for studying measurement and analysis techniques. Our results highlight some pitfalls that may occur when using these UQ methods. For example, when the variation of noise in the training set is small, all methods predicted the same relative uncertainty independently of the inputs. This issue is particularly hard to avoid in BNN. On the other hand, when the test set contains samples far from the training distribution, we found that no methods sufficiently increased the uncertainties associated to their predictions. This problem was particularly clear for CD. In light of these results, we make some recommendations for usage and interpretation of UQ methods.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Analysis of Long-Term Quality Control Data for a 137 Cs Dosimetry Calibration Source

Strict quality assurance programs are required for many radiological applications, but these seldom exist for verifying dosimetry calibration sources. After initial characterization of a dosimetry calibration facility, quality control procedures are recommended to ensure the early detection of any changes or malfunctions. These also result in refined knowledge about average dose rate and experimental variations in dose delivery. This paper describes the implementation of a phase I quality control protocol for a 137 Cs dosimetry calibration source and includes an analysis of the resulting data collected over a 24-mo period. During this time, substantial data was collected to establish trial control limits. Air kerma rate measurements were obtained using an ion chamber and were adjusted for decay, corrected for ambient temperature, pressure and humidity, and then analyzed using quality control charts. Three variations of rational subgrouping methods were used in order to find assignable causes of error, and Nelson's Rules were followed to detect any non-random statistical variations. Measurements were subgrouped according to same-day measurements in order to detect positional errors as well as atmospheric correction errors. Additionally, measurements were subgrouped according to analogous experimental setups in order to detect failure in equipment or incorrect settings. Both were analyzed using the X-bar and R chart method. Similarly, individuals and moving ranges charts were used to carefully examine each position in order to observe any situational errors that may occur which include timing, positional, or interference errors. Each method was successful in identifying unique out-of-control data points that occurred during the phase I application of forming control limits. Furthermore, over the 24-mo period, enough data points were deemed in-control to establish reliable trial limits. Future experiments will include the phase II application of gaining more reliable measurements in order to fine-tune the limits, as well as performing a designed experiment, where variables are purposefully changed in order to test the variation of the data.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2022

Non‐Linear Dimensionality Reduction With a Variational Encoder Decoder to Understand Convective Processes in Climate Models

Abstract Deep learning can accurately represent sub‐grid‐scale convective processes in climate models, learning from high resolution simulations. However, deep learning methods usually lack interpretability due to large internal dimensionality, resulting in reduced trustworthiness in these methods. Here, we use Variational Encoder Decoder structures (VED), a non‐linear dimensionality reduction technique, to learn and understand convective processes in an aquaplanet superparameterized climate model simulation, where deep convective processes are simulated explicitly. We show that similar to previous deep learning studies based on feed‐forward neural nets, the VED is capable of learning and accurately reproducing convective processes. In contrast to past work, we show this can be achieved by compressing the original information into only five latent nodes. As a result, the VED can be used to understand convective processes and delineate modes of convection through the exploration of its latent dimensions. A close investigation of the latent space enables the identification of different convective regimes: (a) stable conditions are clearly distinguished from deep convection with low outgoing longwave radiation and strong precipitation; (b) high optically thin cirrus‐like clouds are separated from low optically thick cumulus clouds; and (c) shallow convective processes are associated with large‐scale moisture content and surface diabatic heating. Our results demonstrate that VEDs can accurately represent convective processes in climate models, while enabling interpretability and better understanding of sub‐grid‐scale physical processes, paving the way to increasingly interpretable machine learning parameterizations with promising generative properties.

54 ENVIRONMENTAL SCIENCES↗