Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “robust estimation”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Wetter California Projected by CMIP6 Models With Observational Constraints Under a High GHG Emission Scenario

As the world's fifth-largest economy entity, California (CA) is vulnerable to climate changes especially the magnitude of winter (wet season) precipitation that is closely linked to regional drought severity, vegetation growth, and wildfire activities. However, the fate of CA precipitation in the future remains highly uncertain, for example, the state-of-the-art Earth system models (ESMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6) projected –3% to +42% and –27% to +63% precipitation changes in northern and central-southern CA, respectively, under the high greenhouse gas (GHG) emission scenario (SSP585). In this work, we applied the Pareto optimality concept and used observed teleconnection patterns (causal relationships between ocean regions and CA precipitation) to mechanistically constrain CMIP6 projected CA precipitations. We estimated that precipitation will robustly increase by 0.4–1.3 mm d –1 (10–34%) over northern CA and increase by 0.1–0.5 mm d –1 (7–32%) over central-southern CA by the end of the 21st century compared with present-day. Up to 71% of ESM projection uncertainties were reduced mainly due to the strong and consistent causal relationship between North American west coast sea level pressure and CA precipitation in both observations and CMIP6 models. Our results suggest that teleconnection patterns are powerful mechanistic constraints that can help explain and reduce uncertainties in ESM projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Dust dominates high-altitude snow darkening and melt over high-mountain Asia

Westerly-driven, long-range transportation of dust particles in elevated aerosol layers (EALs) is a persistent phenomenon during spring and summer over the Indian subcontinent. During the snow accumulation season, EALs transport substantial dust to the snow-covered slopes of high-mountain Asia (HMA). Here we use multiple satellite-based estimates to demonstrate a robust physical association between the EALs and dust-induced snow darkening over HMA. Results from a fully-coupled atmosphere-chemistry-snow model support these observations, revealing a signature of increasing dust-induced snow darkening with surface elevation over HMA that peaks near 4500 m. Moreover, the influence of dust on snow darkening is greater than that of black carbon above 4000 m. Our findings suggest a discernable role of dust in the observed spatial heterogeneity of snowmelt and snowline trends over HMA and highlight an increasing contribution of dust to snowmelt as the snowline rises with warming.

Elevated aerosol layers (EAL), Light absorbing par↗
DOE OSTI2023

First Moments of a Polyhedron Clipped by a Paraboloid

We provide closed-form expressions for the first moments (i.e., the volume and volume-weighted centroid) of a polyhedron clipped by a paraboloid, that is, of a polyhedron intersected with the subset of the three-dimensional real space located on one side of a paraboloid. These closed-form expressions are derived following successive applications of the divergence theorem and the judicious parametrization of the intersection of the polyhedron’s faces with the paraboloid. Here, we provide means for identifying ambiguous discrete intersection topologies, and propose a corrective procedure for preventing their occurence. Finally, we put our proposed closed-form expressions and numerical approach to the test with millions of random and manually engineered polyhedron/paraboloid intersection configurations. The results of these tests show that we are able to provide robust machine-accurate estimates of the first moments at a computational cost that is within one order of magnitude of that of state-of-the-art half-space clipping algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Novel estimation of aerosol processes with particle size distribution measurements: a case study with the TOMAS algorithm v1.0.0

Atmospheric aerosol microphysical processes are a significant source of uncertainty in predicting climate change. Specifically, aerosol nucleation, emissions, and growth rates, which are simulated in chemical transport models to predict the particle size distribution, are not understood well. However, long-term size distribution measurements made at several ground-based sites across Europe implicitly contain information about the processes that created those size distributions. This work aims to extract that information by developing and applying an inverse technique to constrain aerosol emissions as well as nucleation and growth rates based on hourly size distribution measurements. We developed an inverse method based upon process control theory into an online estimation technique to scale aerosol nucleation, emissions, and growth so that the model–measurement bias in three measured aerosol properties exponentially decays. The properties, which are calculated from the measured and predicted size distributions, used to constrain aerosol nucleation, emission, and growth rates are the number of particles with a diameter between 3 and 6 nm, the number with a diameter greater than 10 nm, and the total dry volume of aerosol ( N 3–6 , N 10 , V dry ), respectively. In this paper, we focus on developing and applying the estimation methodology in a zero-dimensional “box” model as a proof of concept before applying it to a three-dimensional simulation in subsequent work. The methodology is first tested on a dataset of synthetic and perfect measurements that span diverse environments in which the true particle emissions, growth, and nucleation rates are known. The inverse technique accurately estimates the aerosol microphysical process rates with an average and maximum error of 2 % and 13 %, respectively. Next, we investigate the effect that measurement noise has on the estimated rates. The method is robust to typical instrument noise in the aerosol properties as there is a negligible increase in the bias of the estimated process rates. Finally, the methodology is applied to long-term datasets of in situ size distribution measurements in western Europe from May 2006 through June 2007. At Melpitz, Germany, and Hyytiälä, Finland, the average diurnal profiles of estimated 3 nm particle formation rates are reasonable, having peaks near noon local time with average peak values of 1 and 0.15 cm –3 s –1 , respectively. The normalized absolute error in estimated N 3–6 , N 10 , and V dry at three European measurement sites is less than 15 %, showing that the estimation framework developed here has potential to decrease model–measurement bias while constraining uncertain aerosol microphysical processes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1983

Estimating time and time-lag in time-of-flight velocimetry

Estimating time and time-lag in time-of-flight velocimeters is investigated. Statistics of a filtered Poisson point process is given. A Maximum Likelihood Estimator is compared with suboptimum estimators in terms of robustness. For a dominating background combined spatial and temporal processing can improve the robustness compared with purely temporal processing. Schemes for the spatial filters are given.

Lading, L.↗
NASA NTRS2004

Filtering in Hybrid Dynamic Bayesian Networks

We demonstrate experimentally that inference in a complex hybrid Dynamic Bayesian Network (DBN) is possible using the 2 - T i e Slice DBN (2T-DBN) from [Koller & Lerner, 20001 to model fault detection in a watertank system. In [Koller & Lerner, 20001 a generic Particle Filter (PF) is used for inference. We extend the experiment and perform approximate inference using The Extended Kalman Filter (EKF) and the Unscented Kalman Filter (UKF). Furthermore, we combine these techniques in a 'non-strict' Rao-Blackwellisation framework and apply it to the watertank system. We show that UKF and UKF in a PF framework outperfom the generic PF, EKF and EKF in a PF framework with respect to accuracy and robustness in terms of estimation RMSE. Especially we demonstrate the superiority of UKF in a PF framework when our beliefs of how data was generated are wrong. We also show that the choice of network structure is very important for the performance of the generic PF and the EKF algorithms, but not for the UKF algorithms. Furthermore, we investigate the influence of data noise in the water[ank simulation. Theory and implementation is based on the theory presented.

Andersen, Morten Nonboe↗
NASA NTRS2002

The SIM astronmetric grid

The Space Interferometry Mission (SIM) is fundamentally a one-dimensional instrument with a 15-degree field-of-regard. Mission objectives require a global reference grid of thousands of well-understood stars with positions known to 4 microarcseconds which will be used to establish the instrument baseline vector during scientific observations. This accuracy will be achieved by frequently observing a set of stars throughout the mission and performing a global fit of the observations to determine position, proper motion and parallax for each star. Each star will be observed approximately 200 times with about 6.5 stars per single instrument field on the sky. We describe the nature of the reference grid, the candidate objects, and the results of simulations demonstrating grid performance, including estimates of the grid robustness when including effects such as instrument drift and possible contamination of the grid star sample by undetected binaries.

SIM grid astrometry Space Interferometry Mission S↗
NASA NTRS2005

An Analysis of Spacecraft Localization from Descent Image Data for Pinpoint Landing on Mars and other Cratered Bodies Data Acquisition

A pinpoint landing capability will be a critical component for many planned NASA missions to Mars and beyond. Implicit in the requirement is the ability to accurately localize the spacecraft with respect to the terrain during descent. In this paper, we present evidence that a vision-based solution using craters as landmarks is both practical and will meet the requirements of next generation missions. Our emphasis in this paper is on the feasibility of such a system in terms of (a) localization accuracy and (b) applicability to Martian terrain. We show that accuracy of well under 100 meters can be expected under suitable conditions. We also present a sensitivity analysis that makes an explicit connection between input data and robustness of our pose estimate. In addition, we present an analysis of the susceptibility of our technique to inherently ambiguous configurations of craters. We show that probability of failure due to such ambiguity is becoming increasingly small.

landforms↗
NASA NTRS2012

A Robust Initialization Scheme for a Lateral Trajectory Optimization Problem with Time of Arrival Windows

We present a robust initialization scheme that estimates parameter values for the numerical solution of a two-point boundary value problem. The two-point boundary value problem formulation stems from the optimization of a cost functional subject to the dynamics of a simplified lateral aircraft model and other constraints. Leveraging regular perturbation methods, initial parameter estimates are analytically determined and used to initialize a gradient descent optimization routine which is shown to rapidly converge over a range of initial aircraft positions and heading angles. Additionally, the velocity of the aircraft is optimized to ensure the trajectory of the aircraft terminates within a desired region in both time and space.

lateral trajectory optimization↗
NASA NTRS2002

The StarLight interferometer architecture and operational concepts

The StarLight flight project was designed to demonstrate the key technologies of spacebome long-baseline stellar interferometry and precision formation flying for potential use on the Terrestrial Planet Finder (TPF) and other future astrophysics missions. Interferometer performance validation could be achieved over a 6-12 month period by obtaining several hundred fringe visibility amplitude measurements for stars in the band 600-1000 nm for a variety of stellar visibilities (0.2-1.0), magnitudes (Mv = 2-5), and baselines (B = 30-125 meters). Interferometry could be performed both in a 1 meter fixed-baseline combiner-only mode and in a two-spacecraft formation mode. In formation mode, the combiner spacecraft would remain at the focus of a virtual parabola, while the collector spacecraft assumed various positions along the parabola such that the two arms of the interferometer remained equal over a variety of separations and bearing angles. Challenges to be encountered in flight include high-bandwidth inter-spacecraft stellar and metrology pointing control, alignment and shear correction, delay and delay-rate estimation, visibility calibration, and robust fringe tracking in the presence of local and inter-spacecraft dynamics. This paper is based on the StarLight project design-capture of March 2002 and will describe the StarLight Interferometer System architecture and selected operational concepts (both of which have relevance to the on-going TPF Technology Program).

stellar↗
NASA NTRS2017

A new concept of stability in orbit propagation, useful for quantifying numerical errors

We present the concept of topological stability in the numerical propagation of orbits, and show how it results in a useful new method for measuring the global numerical error of an orbit propagation. The concept applies to any problem in orbital dynamics. Moreover, it can be extended to any three-dimensional system of di erential equations of second order. In order to assess the topological stability of a given integration a special metric is introduced, which can be used to estimate the numerical errors robustly. The method is particularly well suited for dealing with strongly perturbed and chaotic systems. The construction is based on the constraint imposed by the Hopf map that supports the Kustaanheimo-Stiefel transformation. Generic concepts of stability are translated to KS space.

Pelaez, Jesus↗
DOE OSTI2021

A Scalable Gradient Free Method for Bayesian Experimental Design with Implicit Models

Bayesian experimental design (BED) is to answer the question that how to choose designs that maximize the information gathering. For implicit models, where the likelihood is intractable but sampling is possible, conventional BED methods have difficulties in efficiently estimating the posterior distribution and maximizing the mutual information (MI) between data and parameters. Recent work proposed the use of gradient ascent to maximize a lower bound on MI to deal with these issues. However, the approach requires a sampling path to compute the pathwise gradient of the MI lower bound with respect to the design variables, and such a pathwise gradient is usually inaccessible for implicit models. In this paper, we propose a novel approach that leverages recent advances in stochastic approximate gradient ascent incorporated with a smoothed variational MI estimator for efficient and robust BED. Without the necessity of pathwise gradients, our approach allows the design process to be achieved through a unified procedure with an approximate gradient for implicit models. Several experiments show that our approach outperforms baseline methods, and significantly improves the scalability of BED in high-dimensional problems.

Zhang, Jiaxin↗
DOE OSTI2021

Robust verification of stochastic simulation codes

We introduce a robust verification tool for computational codes, which we call Stochastic Robust Extrapolation based Error Quantification (StREEQ). Unlike the prevalent Grid Convergence Index (GCI) [1] method, our approach is suitable for both stochastic and deterministic computational codes and is generalizable to any number of discretization variables. Building on ideas introduced in the Robust Verification [2] approach, we estimate the converged solution and orders of convergence with uncertainty using multiple fits of a discretization error model. In contrast to Robust Verification, we perform these fits to many bootstrap samples yielding a larger set of predictions with smoother statistics. Here, bootstrap resampling is performed on the lack-of-fit errors for deterministic code responses, and directly on the noisy data set for stochastic responses. This approach lends a degree of robustness to the overall results, capable of yielding precise verification results for sufficiently resolved data sets, and appropriately expanding the uncertainty when the data set does not support a precise result. For stochastic responses, a credibility assessment is also performed to give the analyst an indication of the trustworthiness of the results. Furthermore, this approach is suitable for both code and solution verification, and is particularly useful for solution verification of high-consequence simulations..

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Optimal 1D Ly α forest power spectrum estimation – I. DESI-lite spectra

ABSTRACT The 1D Ly α forest flux power spectrum P1D is sensitive to scales smaller than a typical galaxy survey, and hence ties to the intergalactic medium’s thermal state, suppression from neutrino masses, and new dark matter models. It has emerged as a competitive framework to study new physics, but also has come with various challenges and systematic errors in analysis. In this work, we revisit the optimal quadratic estimator for P1D, which is robust against the relevant problems such as pixel masking, time evolution within spectrum, and quasar continuum errors. We further improve the estimator by introducing a fiducial power spectrum, which enables us to extract more information by alleviating the discreteness of band powers. We meticulously apply our method to synthetic Dark Energy Spectroscopic Instrument (DESI) spectra and demonstrate how the estimator overcomes each challenge. We further apply an optimization scheme that approximates the Fisher matrix to three elements per row and reduces computation time by 60 per cent. We show that we can achieve per cent precision in P1D with 5-yr DESI data in the absence of systematics and provide forecasts for different spectral qualities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Evidential Deep Learning: Enhancing Predictive Uncertainty Estimation for Earth System Science Applications

Abstract Robust quantification of predictive uncertainty is a critical addition needed for machine learning applied to weather and climate problems to improve the understanding of what is driving prediction sensitivity. Ensembles of machine learning models provide predictive uncertainty estimates in a conceptually simple way but require multiple models for training and prediction, increasing computational cost and latency. Parametric deep learning can estimate uncertainty with one model by predicting the parameters of a probability distribution but does not account for epistemic uncertainty. Evidential deep learning, a technique that extends parametric deep learning to higher-order distributions, can account for both aleatoric and epistemic uncertainties with one model. This study compares the uncertainty derived from evidential neural networks to that obtained from ensembles. Through applications of the classification of winter precipitation type and regression of surface-layer fluxes, we show evidential deep learning models attaining predictive accuracy rivaling standard methods while robustly quantifying both sources of uncertainty. We evaluate the uncertainty in terms of how well the predictions are calibrated and how well the uncertainty correlates with prediction error. Analyses of uncertainty in the context of the inputs reveal sensitivities to underlying meteorological processes, facilitating interpretation of the models. The conceptual simplicity, interpretability, and computational efficiency of evidential neural networks make them highly extensible, offering a promising approach for reliable and practical uncertainty quantification in Earth system science modeling. To encourage broader adoption of evidential deep learning, we have developed a new Python package, Machine Integration and Learning for Earth Systems (MILES) group Generalized Uncertainty for Earth System Science (GUESS) (MILES-GUESS) ( https://github.com/ai2es/miles-guess ), that enables users to train and evaluate both evidential and ensemble deep learning. Significance Statement This study demonstrates a new technique, evidential deep learning, for robust and computationally efficient uncertainty quantification in modeling the Earth system. The method integrates probabilistic principles into deep neural networks, enabling the estimation of both aleatoric uncertainty from noisy data and epistemic uncertainty from model limitations using a single model. Our analyses reveal how decomposing these uncertainties provides valuable insights into reliability, accuracy, and model shortcomings. We show that the approach can rival standard methods in classification and regression tasks within atmospheric science while offering practical advantages such as computational efficiency. With further advances, evidential networks have the potential to enhance risk assessment and decision-making across meteorology by improving uncertainty quantification, a longstanding challenge. This work establishes a strong foundation and motivation for the broader adoption of evidential learning, where properly quantifying uncertainties is critical yet lacking.

Schreck, John S.↗
DOE OSTI2021

Novel 2D velocity estimation method for large transient events in plasmas

Dynamics of fast transient events are challenging to be analyzed with high time resolution. Such events can occur in fusion plasmas such as the filaments during edge-localized modes (ELMs). Here, we present a robust method—the spatial displacement estimation—for estimating the displacements of structures with fast dynamics from high spatial and time resolution imaging diagnostics [e.g., gas-puff imaging (GPI)] with sampling time temporal resolution. First, a background suppression method is shown, which suppresses the slowly time-evolving and spatially non-uniform background in the signal. In the second step, a two-dimensional polynomial trend subtraction method is presented to tackle the remaining polynomial order trend in the signal. After performing these pre-processing steps, the spatial displacement of the propagating structure is estimated from the two-dimensional spatial cross-correlation coefficient function calculated between consecutive frames. The method is tested for its robustness and accuracy by simulated Gaussian events and spatially displaced random noise. An example application of the method is presented on propagating ELM filaments measured by the GPI system on the National Spherical Torus Experiment spherical tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Stochastic robustness synthesis for a benchmark problem

Stochastic robustness analysis guides the synthesis of robust linear quadratic Gaussian (LQG) regulators for a benchmark control problem. Probabilities of exceeding allowable design limits, including stability, setting time, and control usage, are estimated by Monte Carlo evaluation. Robust, low-gain compensators that fulfill objectives are designed by numerically minimizing quadratic functions of these probabilities. The method is straightforward and makes use of uncomplicated design principles.

Stengel, Robert F.↗