Search NASA⌕ Search

SEARCH · Search NASA

Results for “truncation uncertainty”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Design of robust failure detection filters

An essential aspect of the design of control systems for large, flexible spacecraft is fault tolerance. Because it is anticipated that a large number of sensors and actuators will be required to realize good control over these assemblies, the detection and isolation of component failures cannot be based on direct comparisons among replicated components. Instead, the notion of 'analytic redundancy' must be employed for the FDI function. Unfortunately this makes the FDI function sensitive to modeling errors which are certain to exist in the large space structure problem due to model truncation and parameter uncertainty. This paper addresses the robustness to model error of one method of FDI residual generation - the failure detection filter. Initial designs were found to be extremely sensitive to modeling error. The sources of this sensitivity are analyzed and modifications to the design are suggested. The improved filter is shown to have much better visibility of the failure signatures relative to the background due to modeling error.

San Martin, A. M.↗
DOE OSTI2026

Connecting ground-state properties of 6 Li to each other and to scattering data

We examine the relationship between the asymptotic normalization coefficient (ANC) of 6 Li and other low-energy observables in the α–deuteron system. Our analysis uses a set of calculations carried out within the ab initio no core shell model with continuum (NCSMC) using a variety of inter-nucleon interactions and basis sizes, and yielding 6 Li deuteron separation energies between 1.3 and 1.8 MeV (Hebborn et al 2022 Phys. Rev. Lett. 129 042503). These NCSMC calculations show that the square of the ANC is strongly correlated with the separation energy over this range. In this work, we investigate the origin of this correlation using the phenomenological R-matrix, a single-channel potential and a perturbative approach. We show that this correlation occurs because the depth of the α–deuteron central potential changes by only a small relative amount as the separation energy varies. We then investigate if the ANC can be accurately extracted from α–deuteron phase shifts in an ideal case in which low-energy data are available and there are no experimental errors. We find that both R-matrix and Coulomb-modified effective-range theory (CM-ERE) yield extracted ANCs close to, although not exactly equal to, the NCSMC value, provided the extrapolation is constrained by the known position of the bound-state pole and at least three terms are included in the fit function. The R-matrix approach converges faster than the CM-ERE as the number of parameters increases and is also more robust against the inclusion of low-energy and high-energy phase shift data. Finally, our study also shows that a naive quantification of uncertainties by comparing different truncations used in both theories is not accurate, and suggests the accuracy of ANCs extracted from phase shift data needs further investigation.

R-matrix↗
DOE OSTI2023

Auxiliary field diffusion Monte Carlo calculations of magnetic moments of light nuclei with chiral effective field theory interactions

We calculate the magnetic moments of light nuclei $( A < 20 )$ using the auxiliary field diffusion Monte Carlo method and local two- and three-nucleon forces with electromagnetic currents from chiral effective field theory. For all nuclei under consideration, we also calculate the ground-state energies and charge radii. We generally find a good agreement with experimental values for all of these observables. For the electromagnetic currents, we explore the impact of employing two different power counting schemes, and study theoretical uncertainties stemming from the truncation of the chiral expansion order by order for select nuclei within these two approaches. In conclusion, we find that it is crucial to employ consistent power counting schemes for interactions and currents to achieve a systematic order-by-order convergence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1985

Numerical and experimental evaluation for single-axis control of an LSS laboratory model

One of the major problems of LSS attitude control design stems from the modeling uncertainties due to modal truncation errors and modal parameter estimation errors. This modeling problem is investigated by evaluating poles and zeros of the frequency response functions, described in terms of constrained or unconstrained modes. The frequency domain approach is utilized to evaluated the results of modal analyses and direct output feedback controller experiments using a simple flexible spacecraft model consisting of a rigid primary body and flexible aluminum beam(s) with a control moment gyro as an actuator. The modal data are compared under the conditions that the model is constrained as the canti-lever and it is rotationally free on a single-axis air bearing table. In addition, the open-loop responses are experimentally examined under the impulsive disturbance to the appendage. The results of the control experiments show that the first vibration mode of the appendage and the interacted rotational motion of the primary body are damped in a short time as predicted by the numerical analysis, to produce damping ratio of 10-20%.

Ohkami, Y.↗
DOE OSTI2023

Revisiting proton–proton fusion in chiral effective field theory

We calculate the S-factor for proton–proton fusion using chiral effective field theory interactions and currents. By performing order-by-order calculations with a variety of chiral interactions that are regularized and calibrated in different ways, we assess the uncertainty in the S-factor from the truncation of the effective field theory expansion and from the sensitivity of the S-factor to the short-distance axial current determined from three- and four-nucleon observables. We find that S(0) = (4.100 ± 0.024(syst) ± 0.013(stat) ± 0.008(g$_{A}$)) × 10$^{−23}$ MeV fm$^{2}$, where the three uncertainties arise, respectively, from the truncation of the effective field theory expansion, use of the two-nucleon axial current fit to few-nucleon observables and variation of the axial coupling constant within the recommended range. The increased value of S(0) compared to previous calculations is mainly driven by an increase in the recommended value for the axial coupling constant and is in agreement with a recent analysis based on pionless effective field theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Analysis of multipacting threshold sensitivity to the random distributions of the secondary electron yield parameters

The way multipacting develops, depends strongly on the secondary emission property of the surface material. The knowledge of secondary electron yield is crucial for accurate prediction of the multipacting threshold. Variations in secondary electron yield parameters from experimental measurements create uncertainty, stemming from handling and surface preparation, and these uncertainties significantly affect multipacting threshold predictions. Despite their significance, the previous studies on the multipacting phenomenon did not adequately address the effect of an assumed random distribution of the secondary emission parameters on the multipacting threshold. Therefore, this paper aims to provide a comprehensive statistical study on how the different random distributions of the secondary emission parameters and, as a result, the uncertainty in the secondary electron yield affect multipacting thresholds. We focus on three commonly used distributions, namely uniform, normal, and truncated normal distributions, to define the uncertainty of random inputs. We use the chaos polynomial expansion method to determine how much each of the random parameters contributes to the multipacting threshold uncertainty. Additionally, we calculate Sobol sensitivity indices to evaluate the impact of the individual parameters or groups of parameters on the model outputs and study how different random distributions of these parameters affected the Sobol index results.

physics↗
DOE OSTI2021

Fiducial $q_T$ resummation of color-singlet processes at N$^3$LL+NNLO

We present a framework for $q_T$ resummation at N 3 LL+NNLO accuracy for arbitrary color-singlet processes based on a factorization theorem in SCET. Our implementation CuTe-MCFM is fully differential in the Born kinematics and matches to large-$q_T$ fixed-order predictions at relative order $α^2_s$. It provides an efficient way to estimate uncertainties from fixed-order truncation, resummation, and parton distribution functions. In addition to W ± , Z and H production, also the diboson processes $\gamma\gamma$, Z$\gamma$, ZH and W ± H are available, including decays. We discuss and exemplify the framework with several direct comparisons to experimental measurements as well as inclusive benchmark results. In particular, we present novel results for $\gamma\gamma$ and Z$\gamma\gamma$ at N 3 LL+NNLO and discuss in detail the power corrections induced by photon isolation requirements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Bayesian analysis of nucleon-nucleon scattering data in pionless effective field theory

We perform Bayesian model calibration of two-nucleon (NN) low-energy constants (LECs) appearing in an NN interaction based on pionless effective field theory (πEFT). The calibration is carried out for potentials constructed using naive dimensional analysis in NN relative momenta (p) up to next-to-leading order [NLO, O(p 2 )] and next-to-next-to-next-to-leading order [N3LO, O(p 4 )]. We consider two classes of pionless πEFT potential: one that acts in all partial waves and another that is dominated by s-wave physics. The two classes produce broadly similar results for calibrations to NN data up to E lab = 5 MeV. Our analysis accounts for the correlated uncertainties that arise from the truncation of the pionless πEFT. We simultaneously estimate both the πEFT LECs and the parameters that quantify the truncation error. This permits the first quantitative estimates of the pionless πEFT breakdown scale, Λ b : the 95% intervals are Λ b ∈[50.11,63.03] MeV at NLO and Λ b ∈[72.27,88.54] MeV at N3LO. Furthermore, invoking naive dimensional analysis for the NN potential, therefore, does not lead to consistent results across orders in pionless πEFT. This exemplifies the possible use of Bayesian tools to identify inconsistencies in a proposed EFT power counting.

Bayesian methods↗
DOE OSTI2021

A Brief Account of Steven Weinberg’s Legacy in ab initio Many-Body Theory

In this contribution to the special issue “Celebrating 30 years of Steven Weinberg’s papers on Nuclear Forces from Chiral Lagrangians,” we emphasize the important role chiral effective field theory has played in leading nuclear physics into a precision era. To this end, we share our perspective on a few of the recent advances made in ab initio calculations of nuclear structure and nuclear matter observables, as well as Bayesian uncertainty quantification of effective field theory truncation errors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Structure and dynamics of open-shell nuclei from spherical coupled-cluster theory

We extend the spherical coupled-cluster ab initio method for open-shell nuclei where two nucleons are removed from a shell subclosure. Following the recent implementation of the two-particle-attached approach [Phys. Rev. C 110, 044306 (2024)], we focus on the two-particle-removed method. Using the equations-of-motion framework, we address both nuclear structure and dipole response functions by coupling coupled-cluster theory with the Lorentz integral transform technique. We perform calculations using chiral interactions, including three-nucleon forces, and estimate many-body uncertainties by comparing different coupled-cluster truncation schemes. Here, we validate our approach by studying ground-state energies, excited states, and electric dipole polarizabilities in the oxygen and calcium isotopic chains. For binding energies and selected low-lying excited states, we achieve an accuracy comparable to that of the established closed-shell coupled-cluster theory and generally agree with experiment. Finally, we underestimate experimental data for electric dipole polarizabilities, particularly in calcium isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

How Well Do We Know the Neutron-Matter Equation of State at the Densities Inside Neutron Stars? A Bayesian Approach with Correlated Uncertainties

Here, we introduce a new framework for quantifying correlated uncertainties of the infinite-matter equation of state derived from chiral effective field theory (𝜒⁢EFT ). Bayesian machine learning via Gaussian processes with physics-based hyperparameters allows us to efficiently quantify and propagate theoretical uncertainties of the equation of state, such as 𝜒⁢EFT truncation errors, to derived quantities. We apply this framework to state-of-the-art many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the 𝜒⁢EFT expansion. This produces the first statistically robust uncertainty estimates for key quantities of neutron stars. We give results up to twice nuclear saturation density for the energy per particle, pressure, and speed of sound of neutron matter, as well as for the nuclear symmetry energy and its derivative. At nuclear saturation density, the predicted symmetry energy and its slope are consistent with experimental constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Reliability of stiffened structural panels: Two examples

The reliability of two graphite-epoxy stiffened panels that contain uncertainties is examined. For one panel, the effect of an overall bow-type initial imperfection is studied. The size of the bow is assumed to be a random variable. The failure mode is buckling. The benefits of quality control are explored by using truncated distributions. For the other panel, the effect of uncertainties in a strain-based failure criterion is studied. The allowable strains are assumed to be random variables. A geometrically nonlinear analysis is used to calculate a detailed strain distribution near an elliptical access hole in a wing panel that was tested to failure. Calculated strains are used to predict failure. Results are compared with the experimental failure load of the panel.

Stroud, W. Jefferson↗
DOE OSTI2023

Muon capture on the deuteron in chiral effective field theory

Here, we consider the capture of a muon on a deuteron. An uncertainty analysis of the dominant channels is important for a careful analysis of forthcoming experimental data. We quantify the theoretical uncertainties of chiral effective-field-theory predictions of the muon-deuteron capture rate from the relevant neutron-neutron partial wave channels in the final state. We study the dependence on the cutoff used to regularize the interactions, low-energy constants calibrated using different fitting data and strategies, and truncation of the effective-field-theory expansion of the currents. Combining these approaches gives as an estimate of $Γ^{1/2}_{μd}$ = 399.1 ± 7.6 ± 4.4 s –1 for capture from the atomic doublet state, and $Γ^{3/2}_{μd}$= 12.31 ± 0.47 ± 0.04 s –1 for capture from the quartet state and the first and second uncertainties given here are due to the effective field theory truncation error and the uncertainty in the axial radius, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Real-time evolution of Anderson impurity models via tensor network influence functionals

In this work, we present and analyze two tensor network-based influence functional approaches for simulating the real-time dynamics of quantum impurity models such as the Anderson model. Via comparison with recent numerically exact simulations, we show that such methods accurately capture the long-time nonequilibrium quench dynamics. The two parameters that must be controlled in these tensor network influence functional approaches are a time discretization (Trotter) error and a bond dimension (tensor network truncation) error. We show that the actual numerical uncertainties are controlled by an intricate interplay of these two approximations, which we demonstrate in different regimes. Our work opens the door to using these tensor network influence functional methods as general impurity solvers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0

Subglacial drainage models represent water flow at the ice–bed interface through coupled distributed and channelized systems to determine water pressure, discharge, and drainage system geometry. While they are used to understand processes such as the relationship between surface melt and ice flow, the number of uncertain model parameters and the computational cost of running models makes it difficult to adequately explore the high-dimensional parameter space and evaluate uncertainty in model predictions. Here, we develop Gaussian process (GP) emulators that make fast predictions with associated uncertainty of subglacial drainage model outputs. Using a truncated principal component (PC) basis representation, we construct a GP emulator for diurnally averaged subglacial water pressure. We also explore emulation of scalar variables describing drainage efficiency and configuration. We train the emulators using ensembles of up to 512 simulations varying eight parameters of the Glacier Drainage System (GlaDS) model on a synthetic domain intended to represent an ice-sheet margin. The emulators make predictions ∼ 1000 times faster than GlaDS simulations, with errors <3 % for the water pressure field and ∼ 5 %–9 % for drainage efficiency and configuration. We apply the emulators to explore the eight-dimensional parameter space by computing variance-based parameter sensitivity indices, finding that three parameters (ice flow coefficient, bed bump aspect ratio, and the subglacial cavity system conductivity) explain 90 % of the variance in modelled water pressure in response to parameter changes. The GP emulator approach described here is well suited to integrating observational data with models to make calibrated, credible predictions of subglacial drainage.

58 GEOSCIENCES↗
DOE OSTI2021

B-PINNs: Bayesian Physics-informal Neural Networks for Forward and Inverse PDE Problems with Noisy Data

We propose a Bayesian physics-informed neural network (B-PINN) to solve both forward and inverse nonlinear problems described by partial differential equations (PDEs) and noisy data. In this Bayesian framework, the Bayesian neural network (BNN) combined with a PINN for PDEs serves as the prior while the Hamiltonian Monte Carlo (HMC) or the variational inference (VI) could serve as an estimator of the posterior. B-PINNs make use of both physical laws and scattered noisy measurements to provide predictions and quantify the aleatoric uncertainty arising from the noisy data in the Bayesian framework. Compared with PINNs, in addition to uncertainty quantification, B-PINNs obtain more accurate predictions in scenarios with large noise due to their capability of avoiding overfitting. We conduct a systematic comparison between the two different approaches for the B-PINNs posterior estimation (i.e., HMC or VI), along with dropout used for quantifying uncertainty in deep neural networks. Our experiments show that HMC is more suitable than VI with mean field Gaussian approximation for the B-PINNs posterior estimation, while dropout employed in PINNs can hardly provide accurate predictions with reasonable uncertainty. Finally, we replace the BNN in the prior with a truncated Karhunen-Loève (KL) expansion combined with HMC or a deep normalizing flow (DNF) model as posterior estimators. The KL is as accurate as BNN and much faster but this framework cannot be easily extended to high-dimensional problems unlike the BNN based framework.

Non linear PDEs, Noisy data, Bysian physics inform↗
DOE OSTI2024

Physics-coupled data-driven design of high-temperature alloys

We present a materials design loop, which streamlines physics-coupled machine learning (ML) surrogate models to discover new alloy chemistries with improved properties. The efficacy is demonstrated by discovering a high-temperature alumina-forming austenitic (AFA) stainless steel with enhanced creep, followed by experimental validation. The ML models have been trained using a well-curated, highly consistent experimental dataset augmented with synthetic microstructural features from a computational thermodynamic approach. We have populated a large number of hypothetical AFA alloys to explore the high-dimensional composition space and have predicted their creep properties by providing the same synthetic input features obtained from the trained ML models. Uncertainties from the ML training were taken as thresholds for truncating predicted results to identify alloys with improved or deteriorated creep. Individual elemental compositions have been determined via probability density distribution analysis from the group of alloys at the top and bottom of the predicted creep values for further virtual and experimental validations. In conclusion, we anticipate that this workflow can be applied to screen desired conditions, such as chemistry and processing parameters, in high-dimensional space through physics-guided data analytics.

Alloy design↗
DOE OSTI2026

Microscopic constraints for the equation of state and structure of neutron stars: A Bayesian model mixing framework

Bayesian model mixing (BMM) is a statistical technique that can combine constraints from different regions of an input space in a principled way. Here we extend our BMM framework for the equation of state (EOS) of strongly interacting matter from symmetric nuclear matter to asymmetric matter, specifically focusing on zero-temperature, charge-neutral, 𝛽-equilibrated matter. We use Gaussian processes (GPs) to infer constraints on the neutron-star matter EOS at intermediate densities from two different microscopic theories: chiral effective-field theory (𝜒⁢EFT) at baryon densities around nuclear saturation, 𝑛 𝐵 ∼ 𝑛 0 , and perturbative QCD at asymptotically high baryon densities, 𝑛 𝐵 ⩾ 20⁢𝑛 0 . The uncertainties of the 𝜒⁢EFT and pQCD EOSs are obtained using the BUQEYE truncation error model. We demonstrate the flexibility of our framework through the use of two categories of GP kernels: conventional stationary kernels and a nonstationary changepoint kernel. We use the latter to explore potential constraints on the dense matter EOS by including exogenous data representing theory predictions and heavy-ion collision measurements at densities ⩾ 2⁢𝑛 0 . We also use our EOSs to obtain neutron-star mass-radius relations and their uncertainties. Finally, our framework, whose implementation will be available through a GitHub repository, provides a prior distribution for the EOS that can be used in large-scale neutron-star inference frameworks.

Bayesian methods↗