Search NASASearch

SEARCH · Search NASA

Results for “gaussian method”

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 127 records · Page 7

Error mitigation in variational quantum eigensolvers using tailored probabilistic machine learning

Quantum computing technology has the potential to revolutionize the simulation of materials and molecules in the near future. A primary challenge in achieving near-term quantum advantage is effectively mitigating the noise effects inherent in current quantum processing units (QPUs). This challenge is also decisive in the context of quantum-classical hybrid schemes employing variational quantum eigensolvers (VQEs) that have attracted significant interest in recent years. In this paper, we present a method that employs parametric Gaussian process regression (GPR) within an active learning framework to mitigate noise in quantum computations, focusing on VQEs. Our approach, grounded in probabilistic machine learning, exploits a custom prior based on the VQE ansatz to capture the underlying correlations between VQE outputs for different variational parameters, thereby enhancing both accuracy and efficiency. We demonstrate the effectiveness of our method on a two-site Anderson impurity model and a eight-site Heisenberg model, using the IBM open-source quantum computing framework, Qiskit, showcasing substantial improvements in the accuracy of VQE outputs while reducing the number of direct QPU energy evaluations. This paper contributes to the ongoing efforts in quantum-error mitigation and optimization, bringing us a step closer to realizing the potential of quantum computing in quantum matter simulations. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING

User's manual for the Gaussian windows program

'Gaussian Windows' is a method for exploring a set of multivariate data, in order to estimate the shape of the underlying density function. The method can be used to find and describe structural features in the data. The method is described in two earlier papers. I assume that the reader has access to both of these papers, so I will not repeat material from them. The program described herein is written in BASIC and it runs on an IBM PC or PS/2 with the DOS 3.3 operating system. Although the program is slow and has limited memory space, it is adequate for experimenting with the method. Since it is written in BASIC, it is relatively easy to modify. The program and some related files are available on a 3-inch diskette. A listing of the program is also available. This user's manual explains the use of the program. First, it gives a brief tutorial, illustrating some of the program's features with a set of artificial data. Then, it describes the results displayed after the program does a Gaussian window, and it explains each of the items on the various menus.

Jaeckel, Louis A.

Uncertainty quantification of a physics-informed model based on sparse identification of a Thermal Energy Distribution System

Integrated energy systems (IES)s are crucial for enhancing the economy and efficiency of power generation sources (e.g., nuclear energy) necessary to unleash American energy dominance. These systems can be integrated with thermal energy storage (TES) and intermittent renewable energies to optimize overall energy use, peak-load regulation, and demand-side responses. However, the stabilization of energy generation, transport, and utilization introduces operational complexities that exceed the challenges of managing each sub-component individually. Currently, though IESs rely on human operators for efficiency and stability, reducing human error risk and enhancing performance through automation is highly desirable. Recent advances at Idaho National Laboratory have demonstrated successful control of the Thermal Energy Distributed System (TEDS). However, the automatic control system depends on a deterministic Sparse Identification of Nonlinear Dynamics with Control (SINDyC) model, which are trained based on simulation data from physics-based simulations. Because of uncertainties in physics-based simulation, SINDyC model results in large discrepancies against experimental data and cannot be reliably used in automatic control. In this paper, we present an innovative approach to address these discrepancies by quantifying uncertainties and developing a more robust model. We first generated trajectories by using first-principles physics codes to encapsulate the experiment. Next, we trained thousands of models by randomly sampling these trajectories. We then collapsed all those models into one probabilistic SINDyC by fitting a multivariate Gaussian distribution onto the resulting coefficient’s distribution. Despite its simplicity, our approach successfully produced 95% confidence intervals that captured the experimental trajectories. It even did so with a higher probability and better U-pooling score across six of the seven relevant quantities of interest (QoIs), as compared to other classical approaches. In conclusion, ongoing research is focusing on generating new experimental trajectories to validate this approach, and on employing Bayesian calibration to refine parametric uncertainties and guide future model development efforts.

22 - GENERAL STUDIES OF NUCLEAR REACTORS

Attitude and vibration control of a large flexible space-based antenna

The problem of control systems synthesis is considered for controlling the rigid body attitude and elastic motion of a large deployable space based antenna. Two methods for control systems synthesis are considered. The first method utilizes the stability and robustness properties of the controller consisting of torque actuators and collocated attitude and rate sensors. The second method is based on the linear quadratic Gaussian (LQG) control theory. A combination of the two methods, which results in a two level hierarchical control system, is also briefly discussed. The performance of the controllers is analyzed by computing the variances of pointing errors, feed misalignment errors and surface contour errors in the presence of sensor and actuator noise.

Joshi, S. M.

Orphaned oil and gas well methane emission rates quantified using Gaussian plume inversions of ambient observations

Abstract. Annually, ∼ 3.6 million abandoned oil and gas wells in the US emit a combined ∼ 2.6 Tg methane (CH4), adversely affecting climate and regional air quality. However, these estimates depend on emission factors derived from measuring subpopulations of wells that vary by orders of magnitude due to very limited field sampling and poorly characterized distributions. Currently, US protocols to remediate orphaned wells lacks standardized quantification methods needed to both prioritize plugging and account for emission reductions. Therefore, sensitive, reliable, affordable, and scalable CH4 flux quantification methods are needed. We report the use of a simple Gaussian plume method where the dispersion parameters are constrained by in situ ground measurements of CH4 concentration at four locations 7.5–49 m downwind of the orphan well as well as local winds to estimate the leak rate from an orphan well. We derive a flux of 10.53 ± 1.16 kg CH4 h−1 during a venting procedure in April 2023 that agrees with the directly measured volumetric flow rate of 9.00 ± 0.25 kg CH4 h−1. This is 71 % greater than the 5.3 kg CH4 h−1 flux measured 7 months prior. Additionally, we discovered a secondary leak through the surface casing inferred as 0.43–0.67 kg CH4 h−1 both by our ground Gaussian analysis and by transecting the plume with an uncrewed aerial system (UAS). We show that in situ determination of the dispersion parameters used in our Gaussian inversions allows us to measure methane emissions to 15 % accuracy, significantly reducing errors when compared to the standard practice of assuming stability class. Our results help develop simpler methods and protocols for robust orphan well emission quantification that can be used for reporting.

Follansbee, Emily

Uncertainty guided online ensemble for non-stationary data streams in fusion science

Machine Learning (ML) is poised to play a pivotal role in the development and operation of next-generation fusion devices. Fusion data shows non-stationary behavior with distribution drifts, resulted by both experimental evolution and machine wear-and-tear. ML models assume stationary distribution and fail to maintain performance when encountered with such non-stationary data streams. Online learning techniques have been leveraged in other domains, however it has been largely unexplored for fusion applications. In this paper, we investigate online learning for continuous adaptation to drifting data streams in the prediction of Toroidal Field (TF) coils deflection at the DIII-D fusion facility. We further address the short-term performance degradation inherent to standard online learning, which arises because ground truth is unavailable at prediction time. To mitigate this issue, we propose an uncertainty-guided online ensemble framework. The method leverages the Deep Gaussian Process Approximation (DGPA) for calibrated uncertainty estimation and uses these uncertainty measures to guide a meta-algorithm that aggregates predictions from learners trained over different historical horizons. Our results show that online learning reduces prediction error by 80% compared to a static model. The online ensemble and the proposed uncertainty-guided ensemble further reduce error by approximately 6%, and 10% respectively, relative to standard single-model online learning, while also providing calibrated uncertainty estimates to support operational decision-making.

AI

Simulations of Attosecond Metallization in Quartz and Diamond Probed with Inner-Shell Transient Absorption Spectroscopy

When dielectrics are hit with intense infrared (IR) laser pulses, transient metalization can occur. The initial attosecond dynamics behind this metallization are not entirely understood. Therefore, simulations are needed to understand this process and to help interpret experimental observations of it, such as with attosecond transient absorption (ATA). In this paper, we present first-principles simulations of ATA based on bulk-mimicking clusters and real-time time-dependent density functional theory (RT-TDDFT), with Koopmans-tuned range-separated hybrid functionals and Gaussian basis sets. Our method gives good agreement with the experiment for the breakdown threshold in silica and diamond. This breakdown voltage corresponds to a Keldysh parameter of approximately one and thus involves a transition to a regime where the dynamics are driven by tunneling. Pumping at an amplitude just below this value causes a mixture of multiphoton and tunneling excitations across the band gap to occur. The computed extreme ultraviolet and X-ray attosecond transient spectra also agree well with the experiment and show a decrease in optical density due to the transient population of the conduction band from the IR field. First-principles approaches such as this are valuable for interpreting the complicated modulations in a spectrum and for guiding future attosecond experiments on solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Inference of response functions with the help of machine-learning algorithms

Response functions are a key quantity to describe the near-equilibrium dynamics of strongly interacting many-body systems. Recent techniques that attempt to overcome the challenges of calculating these ab initio have employed expansions in terms of orthogonal polynomials. We employ a neural network prediction algorithm to reconstruct a response function 𝑆⁡(𝜔) defined over a range in frequencies 𝜔. Here, we represent the calculated response function as a truncated Chebyshev series whose coefficients can be optimized to reduce the representation error. We compare the quality of response functions obtained using coefficients calculated using a neural network (NN) algorithm with those computed using the Gaussian integral transform (GIT) method. In the regime where only a small number of terms in the Chebyshev series are retained, we find that the NN scheme outperforms the GIT method.

Kurkcuoglu, Doga Murat [Fermi National Accelerator

A matrix equation solution by an optimization technique

The computer solution of matrix equations is often difficult to accomplish due to an ill-conditioned matrix or high noise levels. Two methods of solution are compared for matrices of various degrees of ill-conditioning and for various noise levels in the right hand side vector. One method employs the usual Gaussian elimination. The other solves the equation by an optimization technique and employs a function minimization subroutine.

Johnson, M. J.

Attitude and vibration control of a large flexible space-based antenna

Control systems synthesis is considered for controlling the rigid body attitude and elastic motion of a large deployable space-based antenna. Two methods for control systems synthesis are considered. The first method utilizes the stability and robustness properties of the controller consisting of torque actuators and collocated attitude and rate sensors. The second method is based on the linear-quadratic-Gaussian control theory. A combination of the two methods, which results in a two level hierarchical control system, is also briefly discussed. The performance of the controllers is analyzed by computing the variances of pointing errors, feed misalignment errors and surface contour errors in the presence of sensor and actuator noise.

Joshi, S. M.

N(+)-N long-range interaction energies and resonance charge exchange

The aerothermodynamic studies of proposed space missions require atmospheric charge-transfer data. N2(+) eigenstate energies are calculated with use of the complete-active-space self-consistent-field method with an extended Gaussian basis set. The N(+)-N charge-exchange cross section, determined from these energies, agrees with merged-beam measurements. This contradicts the previous theoretical conclusion. A simple physical description of the long-range interaction is presented and should expedite future charge-transfer studies.

Stallcop, J. R.

Aeroelastic control of oblique-wing aircraft

The U.S. Navy and NASA are currently involved in the design and development of an unsymmetric-skew-wing aircraft capable of 65 deg wing sweep and flight at Mach 1.6. A generic skew-wing aircraft model was developed for 45 deg wing skew at a flight condition of Mach 0.70 and 3048 m altitude. At this flight condition the aircraft has a wing flutter mode. An active implementable control law was developed using the linear quadratic Gaussian design technique. A method of modal residualization was used to reduce the order of the controller used for flutter suppression.

Burken, J. J.

Aeroelastic control of oblique-wing aircraft

The U.S. Navy and NASA are currently involved in the design and development of an unsymmetric-skew-wing aircraft capable of 65 deg wing sweep and flight at Mach 1.6. A generic skew-wing aircraft model was developed for 45 deg wing skew at a flight condition of Mach 0.70 and 3048 m altitude. At this flight condition the aircraft has a wing flutter mode. An active implementable control law was developed using the linear quadratic Gaussian design technique. A method of modal residualization was used to reduce the order of the controller used for flutter suppression.

Burken, J. J.

Modern CACSD using the Robust-Control Toolbox

The Robust-Control Toolbox is a collection of 40 M-files which extend the capability of PC/PRO-MATLAB to do modern multivariable robust control system design. Included are robust analysis tools like singular values and structured singular values, robust synthesis tools like continuous/discrete H(exp 2)/H infinity synthesis and Linear Quadratic Gaussian Loop Transfer Recovery methods and a variety of robust model reduction tools such as Hankel approximation, balanced truncation and balanced stochastic truncation, etc. The capabilities of the toolbox are described and illustated with examples to show how easily they can be used in practice. Examples include structured singular value analysis, H infinity loop-shaping and large space structure model reduction.

Chiang, Richard Y.

On three dimensional object recognition and pose-determination: An abstraction based approach

A method of computing reliable Gaussian and mean curvature sign-map descriptors from the polynomial approximation of surfaces was demonstrated. Such descriptors which are invariant under perspective variation are suitable for hypothesis generation. A means for determining the pose of constructed geometric forms whose algebraic surface descriptors are nonlinear in terms of their orienting parameters was developed. This was done by means of linear functions which are capable of approximating nonlinear forms and determining their parameters. It was shown that biquadratic surfaces are suitable companion linear forms for cylindrical approximation and parameter estimation. The estimates provided the initial parametric approximations necessary for a nonlinear regression stage to fine tune the estimates by fitting the actual nonlinear form to the data. A hypothesis-based split-merge algorithm for extraction and pose determination of cylinders and planes which merge smoothly into other surfaces was developed. It was shown that all split-merge algorithms are hypothesis-based. A finite-state algorithm for the extraction of the boundaries of run-length regions was developed. The computation takes advantage of the run list topology and boundary direction constraints implicit in the run-length encoding.

Quek, Kok How Francis

Analysis of EM penetration into and scattering by electrically large open waveguide cavities using Gaussian beam shooting

Gaussian beam (GB) representation methods are used to analyze the electromagnetic coupling into and the scattering by a large nonuniform cavity. The aperture field in the cavity is decomposed into beams using the Gabor expansion, and shooting techniques are then employed. The method is illustrated only for the two-dimensional (2-D) case. The GBs are tracked axially using the rules of beam optics which ignore any beam distortion upon reflection at the walls. The effects of beam distortion are not significant for relatively slowly varying waveguide cavities. The field scattered into the exterior by the termination within the cavity is found using a reciprocity integral formulation which requires a knowledge of the beam fields near the termination. Numerical results based on this GB approach are presented and compared with results based on an independent reference solution.

Burkholder, Robert J.

Implementation of a Collision Probability Prediction Technique for Constellation Maneuver Planning

On March 22, 2006, the Space Technology 5 (ST5) constellation spacecraft were successfully delivered to orbit by a Pegasus XI, launch vehicle. An unexpected relative motion experienced by the constellation after orbit insertion brought about a problem. Soon after launch the observed relative position of the inert rocket body was between the leading and the middle spacecraft within the constellation. The successful planning and execution of an orbit maneuver that would create a fly-by of the rocket body was required to establish the.formation. This maneuver would create a close approach that needed to conform to predefined collision probability requirements. On April 21, 2006, the ST5 "155" spacecraft performed a large orbit maneuver and successfully passed the inert Pegasus 3rd Stage Rocket Body on April 30, 2006 15:20 UTC at a distance of 2.55 km with a Probability of Collision of less than 1.0E-06. This paper will outline the technique that was implemented to establish the safe planning and execution of the fly-by maneuver. The method makes use of Gaussian distribution models of state covariance to determine underlying probabilities of collision that arise under low velocity encounters. Specific numerical examples used for this analysis are discussed in detail. The mechanics of this technique are explained to foster deeper understanding of the concepts presented and to improve existing processes for use in future constellation maneuver planning.

Concha, Marco a.