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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
NASA NTRS2019

SO2 Emission Estimates Using OMI SO2 Retrievals for 2005-2017

SO2 column densities from OMI provide important information on emission trends and missing sources, but there are discrepancies between different retrieval products. We employ three OMI SO2 retrieval products (NASA standard (SP), NASA prototype, and BIRA) to study the magnitude and trend of SO2 emissions. SO2 column densities from these retrievals are most consistent when viewing angles and solar zenith angles are small, suggesting more robust emission estimates in summer and at low latitudes. We then apply a hybrid 4DVar/mass balance emission inversion to derive monthly SO2 emissions from the NASA SP and BIRA products. Compared to HTAPv2 emissions in 2010, both posterior emission estimates are lower in US, India and Southeast China, but show different changes of emissions in North China Plain. The discrepancies between monthly NASA and BIRA posterior emissions in 2010 are less than or equal to 17% in China and 34% in India. SO2 emissions increase from 2005 to 2016 by 35% (NASA) 48% (BIRA) in India, but decrease in China by 23% (NASA) 33% (BIRA) since 2008. Compared to insitu measurements, the posterior GEOSChem surface SO2 concentrations have reduced NMB in China, the US, and India but not in South Korea in 2010. BIRA posteriors have better consistency with the annual growth rate of surface SO2 measurement in China and spatial variability of SO2 concentration in China, South Korea and India, whereas NASA SP posteriors have better seasonality. These evaluations demonstrate the capability to recover SO2 emissions using OMI observations.

Data assimilation↗
NASA NTRS2007

Using Remotely-Sensed Estimates of Soil Moisture to Infer Soil Texture and Hydraulic Properties across a Semi-arid Watershed

Near-surface soil moisture is a critical component of land surface energy and water balance studies encompassing a wide range of disciplines. However, the processes of infiltration, runoff, and evapotranspiration in the vadose zone of the soil are not easy to quantify or predict because of the difficulty in accurately representing soil texture and hydraulic properties in land surface models. This study approaches the problem of parameterizing soils from a unique perspective based on components originally developed for operational estimation of soil moisture for mobility assessments. Estimates of near-surface soil moisture derived from passive (L-band) microwave remote sensing were acquired on six dates during the Monsoon '90 experiment in southeastern Arizona, and used to calibrate hydraulic properties in an offline land surface model and infer information on the soil conditions of the region. Specifically, a robust parameter estimation tool (PEST) was used to calibrate the Noah land surface model and run at very high spatial resolution across the Walnut Gulch Experimental Watershed. Errors in simulated versus observed soil moisture were minimized by adjusting the soil texture, which in turn controls the hydraulic properties through the use of pedotransfer functions. By estimating a continuous range of widely applicable soil properties such as sand, silt, and clay percentages rather than applying rigid soil texture classes, lookup tables, or large parameter sets as in previous studies, the physical accuracy and consistency of the resulting soils could then be assessed. In addition, the sensitivity of this calibration method to the number and timing of microwave retrievals is determined in relation to the temporal patterns in precipitation and soil drying. The resultant soil properties were applied to an extended time period demonstrating the improvement in simulated soil moisture over that using default or county-level soil parameters. The methodology is also applied to an independent case at Walnut Gulch using a new soil moisture product from active (C-band) radar imagery with much lower spatial and temporal resolution. Overall, results demonstrate the potential to gain physically meaningful soils information using simple parameter estimation with few but appropriately timed remote sensing retrievals.

Santanello, Joseph A.↗
DOE OSTI2020

Stochasticity and robustness in spiking neural networks

Despite drawing inspiration from biological systems which are inherently noisy and variable, artificial neural networks have been shown to require precise weights to carry out the task which they are trained to accomplish. This creates a challenge when adapting these artificial networks to specialized execution platforms which may encode weights in a manner which restricts their accuracy and/or precision.Reflecting back on the non-idealities which are observed in biological systems, we investigated the effect these properties have on the robustness of spiking neural networks under perturbations to weights. First, we examined techniques extant in conventional neural networks which resemble noisy processes, and postulated they may produce similar beneficial effects in spiking neural networks. Second, we evolved a set of spiking neural networks utilizing biological non-idealities to solve a pole-balancing task, and estimated their robustness. We showed it is higher in networks using noisy neurons, and demonstrated that one of these networks can perform well under the variance expected when a hafnium-oxide based resistive memory is used to encode synaptic weights. Lastly, we trained a series of networks using a surrogate gradient method on the MNIST classification task. We confirmed that these networks demonstrate similar trends in robustness to the evolved networks. We discuss these results and argue that they display empirical evidence supporting the role of noise as a regularizer which can increase network robustness.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2021

Estimating Bare Earth in Sparse Boreal Forests With WorldView Stereo Imagery

Circumboreal forests are currently experiencing rapid climate warming which is altering their structure, productivity, and status as a carbon sink. Very-high resolution (VHR, <2 m) stereo derived digital surface models are available to monitor these forests, but, a similar resolution digital terrain model (DTM) is required to extract information about tree height, which is often used to estimate carbon content. To our knowledge, no openly available VHR DTM currently exists. To address this need, we developed approaches to extract DTMs by filtering VHR stereo point clouds (PCs) in sparse canopies of Alaska. Our evaluation consisted of two stereo processing methods with three PC search radii, at six different tree canopy cover (TCC) intervals. We found VHR DTMs were robust for estimating bare ground at TCC intervals less than 40% with vertical errors <1.6 m using airborne small footprint light detection and ranging (LiDAR) as reference.

Christopher S R Neigh↗
NASA NTRS2021

Estimating Bare Earth in Sparse Boreal Forests With WorldView Stereo Imagery

Circumboreal forests are currently experiencing rapid climate warming which is altering their structure, productivity, and status as a carbon sink. Very-high resolution (VHR, <2 m) stereo derived digital surface models are available to monitor these forests, but, a similar resolution digital terrain model (DTM) is required to extract information about tree height, which is often used to estimate carbon content. To our knowledge, no openly available VHR DTM currently exists. To address this need, we developed approaches to extract DTMs by filtering VHR stereo point clouds (PCs) in sparse canopies of Alaska. Our evaluation consisted of two stereo processing methods with three PC search radii, at six different tree canopy cover (TCC) intervals. We found VHR DTMs were robust for estimating bare ground at TCC intervals less than 40% with vertical errors <1.6 m using airborne small footprint light detection and ranging (LiDAR) as reference.

Christopher S. R. Neigh↗
NASA NTRS1991

Stochastic robustness of linear time-invariant control systems

A simple numerical procedure for estimating the stochastic robustness of a linear, time-invariant system is described. Monte Carlo evaluation of the system's eigenvalues allows the probability of instability and the related stochastic root locus to be estimated. This definition of robustness is an alternative to existing deterministic definitions that address both structured and unstructured parameter variations directly. This analysis approach treats not only Gaussian parameter uncertainties but also non-Gaussian cases, including uncertain-but-bounded variations. Trivial extensions of the procedure permit alternate discriminants to be considered. Thus the probabilities that stipulated degrees of instability will be exceeded or that closed-loop roots will leave desirable regions can also be estimated. Results are particularly amenable to graphical presentation.

Stengel, Robert F.↗
NASA NTRS1997

Exponential Approximations Using Fourier Series Partial Sums

The problem of accurately reconstructing a piece-wise smooth, 2(pi)-periodic function f and its first few derivatives, given only a truncated Fourier series representation of f, is studied and solved. The reconstruction process is divided into two steps. In the first step, the first 2N + 1 Fourier coefficients of f are used to approximate the locations and magnitudes of the discontinuities in f and its first M derivatives. This is accomplished by first finding initial estimates of these quantities based on certain properties of Gibbs phenomenon, and then refining these estimates by fitting the asymptotic form of the Fourier coefficients to the given coefficients using a least-squares approach. It is conjectured that the locations of the singularities are approximated to within O(N(sup -M-2), and the associated jump of the k(sup th) derivative of f is approximated to within O(N(sup -M-l+k), as N approaches infinity, and the method is robust. These estimates are then used with a class of singular basis functions, which have certain 'built-in' singularities, to construct a new sequence of approximations to f. Each of these new approximations is the sum of a piecewise smooth function and a new Fourier series partial sum. When N is proportional to M, it is shown that these new approximations, and their derivatives, converge exponentially in the maximum norm to f, and its corresponding derivatives, except in the union of a finite number of small open intervals containing the points of singularity of f. The total measure of these intervals decreases exponentially to zero as M approaches infinity. The technique is illustrated with several examples.

Banerjee, Nana S.↗
NASA NTRS1995

Attitude determination using an adaptive multiple model filtering Scheme

Attitude determination has been considered as a permanent topic of active research and perhaps remaining as a forever-lasting interest for spacecraft system designers. Its role is to provide a reference for controls such as pointing the directional antennas or solar panels, stabilizing the spacecraft or maneuvering the spacecraft to a new orbit. Least Square Estimation (LSE) technique was utilized to provide attitude determination for the Nimbus 6 and G. Despite its poor performance (estimation accuracy consideration), LSE was considered as an effective and practical approach to meet the urgent need and requirement back in the 70's. One reason for this poor performance associated with the LSE scheme is the lack of dynamic filtering or 'compensation'. In other words, the scheme is based totally on the measurements and no attempts were made to model the dynamic equations of motion of the spacecraft. We propose an adaptive filtering approach which employs a bank of Kalman filters to perform robust attitude estimation. The proposed approach, whose architecture is depicted, is essentially based on the latest proof on the interactive multiple model design framework to handle the unknown of the system noise characteristics or statistics. The concept fundamentally employs a bank of Kalman filter or submodel, instead of using fixed values for the system noise statistics for each submodel (per operating condition) as the traditional multiple model approach does, we use an on-line dynamic system noise identifier to 'identify' the system noise level (statistics) and update the filter noise statistics using 'live' information from the sensor model. The advanced noise identifier, whose architecture is also shown, is implemented using an advanced system identifier. To insure the robust performance for the proposed advanced system identifier, it is also further reinforced by a learning system which is implemented (in the outer loop) using neural networks to identify other unknown quantities such as spacecraft dynamics parameters, gyro biases, dynamic disturbances, or environment variations.

Lam, Quang↗
DOE OSTI2022

Active learning with multifidelity modeling for efficient rare event simulation

Here, while multifidelity modeling provides a cost-effective way to conduct uncertainty quantification with computationally expensive models, much greater efficiency can be achieved by adaptively deciding the number of required high-fidelity (HF) simulations, depending on the type and complexity of the problem and the desired accuracy in the results. We propose a framework for active learning with multifidelity modeling emphasizing the efficient estimation of rare events. Our framework works by fusing a low-fidelity (LF) prediction with an HF-inferred correction, filtering the corrected LF prediction to decide whether to call the high-fidelity model, and for enhanced subsequent accuracy, adapting the correction for the LF prediction after every HF model call. The framework does not make any assumptions as to the LF model type or its correlations with the HF model. In addition, for improved robustness when estimating smaller failure probabilities, we propose using dynamic active learning functions that decide when to call the HF model. We demonstrate our framework using several academic case studies (including some high-dimensional problems) and two finite element model case studies: estimating Navier-Stokes velocities using the Stokes approximation and estimating stresses in a transversely isotropic model subjected to displacements via a coarsely meshed isotropic model. Across these case studies, not only did the proposed framework estimate the failure probabilities accurately, but compared with either Monte Carlo or a standard variance reduction method, it also required only a small fraction of the calls to the HF model.

42 ENGINEERING↗
NASA NTRS2020

Dynamic Medical Risk Assessment Supported by Inference Networks

The Human Research Program's next generation risk model, the Medical Extensible Dynamic Probabilistic Risk Assessment Tool (MEDPRAT), provides researchers with estimates of astronaut medical and health risk. MEDPRATs Susceptibility Inference Network (SIN) facilitates the dynamic component of this tool. The SIN provides MEDPRAT with a "memory" that allows the system to use knowledge of what simulation events have occurred to alter the representative probability that future simulation events will occur during a given trial. The SIN allows for the medical events being simulated to be related to, and influence the likelihood of one another, providing a more robust risk estimate. We present initial work of our efforts to mathematically quantify and represent these dependent relationships between medical events.

McIntyre, L.↗
NASA NTRS2002

Simulation of TRMM Microwave Imager Brightness Temperature using Precipitation Radar Reflectivity for Convective and Stratiform Rain Areas over Land

Rain is highly variable in space and time. In order to measure rainfall over global land with satellites, we need observations with very high spatial resolution and frequency in time. On board the Tropical Rainfall Measuring Mission (TRMM) satellite, the Precipitation Radar (PR) and Microwave Imager (TMI) are flown together for the purpose of estimating rain rate. The basic method to estimate rain from PR has been developed over the past several decades. On the other hand, the TMI method of rain estimation is still in the state development, particularly over land. The objective of this technical memorandum is to develop a theoretical framework that helps relate the observations made by these two instruments. The principle result of this study is that in order to match the PR observations with the TMI observations in convective rain areas, a mixed layer of graupel and supercooled water drops above the freezing level is needed. On the other hand, to match these observations in the stratiform region, a layer of snowflakes with appropriate densities above the freezing level, and a melting layer below the freezing level, are needed. This understanding can lead to a robust rainfall estimation technique from the microwave radiometer observations.

Prabhakara, C.↗
DOE OSTI2021

Gap-filling eddy covariance methane fluxes: Comparison of machine learning model predictions and uncertainties at FLUXNET-CH4 wetlands

Time series of methane fluxes measured by eddy-covariance require gap-filling to estimate annual emissions. Gap-filling methane fluxes is challenging because of high variability and complex responses to multiple drivers. To date, there is no widely established gap-filling standard for methane, with regards both to the best model algorithms and predictors. In this study, we address the need for standardization by synthesizing results of gap-filling methods applied at 17 wetland sites spanning boreal to tropical regions including all major wetlands classes and two rice paddies. We introduce new procedures for: 1) creating realistic artificial gap scenarios, 2) training and evaluating gap-filling models without overstating performance, and 3) predicting half-hourly methane fluxes and annual emissions with robust uncertainty estimates. We tested a conventional method (marginal distribution sampling) and four machine learning algorithms - penalized linear regression, artificial neural networks, random forests, and boosted decision trees - and four predictor sets, including temporal, meteorological, ecosystem carbon and energy flux, and soil predictors. We find that the conventional method can achieve similar median performance to the machine learning models but is worse than the best machine learning models and relatively insensitive to predictor choices. Of the machine learning models, decision tree algorithms performed the best in cross-validation experiments, even with a baseline predictor set, and artificial neural networks showed comparable performance when using all predictors. Soil temperature was frequently the most important predictor whilst water table depth was important at sites with substantial water table fluctuations, highlighting the value of data on soil conditions. Raw gap-filling uncertainties from the machine learning models were underestimated and we propose a method to calibrate uncertainties to observations. Finally, we gap-fill and provide summary evaluation metrics for all 81 sites in the FLUXNET-CH4 community dataset and publicly release the python code for model development, evaluation, and uncertainty estimation.

42 ENGINEERING↗
DOE OSTI2023

Daily and Seasonal Variations of Shallow Seismic Velocities in Southern California From Joint Analysis of H/V Ratios and Autocorrelations of Seismic Waveforms

Establishing a baseline of ongoing secular velocity variations at the subsurface can improve the accuracy of detecting and interpreting short-term velocity changes, and advance the understanding of observed seismic motions and the behavior of subsurface materials. Toward these goals, we develop and apply a deconvolved autocorrelation (DA) method to estimate regional daily and seasonal changes of seismic velocities in southern California. The DA method combines advantages of traditional autocorrelation and Horizontal-to-Vertical Spectral Ratio, and is used to analyze over 10 years of data recorded by 50 stations. The results indicate widespread daily and seasonal changes of up to 10% and 4%, respectively, in the top tens of meters of the crust. The thickness of the surface layer, distance from the coast, and topographic variations are important factors controlling the amplitudes of the resolved velocity variations. The results suggest that changes of soil moisture and thermoelastic strain are likely dominant factors affecting the daily and seasonal variations, respectively. The developed DA method can improve the accuracy and robustness of estimated changes of subsurface materials at other locations.

58 GEOSCIENCES↗
DOE OSTI2022

Transportable absolute yields of underground nuclear explosions. Application to the North Korean nuclear tests

The method developed and presented here improves our nuclear explosion monitoring capability by providing a new independent and transportable method to determine the explosive yields of seismically recorded events associated with foreign weapons development. Here we report new independent absolute yield and depth of burial estimates of the North Korean nuclear tests which are commensurate with previously determined source parameters. The key advantage of this method over other methods which estimate absolute explosive yields is that this method does not require any a priori calibration, and can be immediately applied to any region of interest. This method uses the source information obtained from the spectral ratios of envelopes of measured seismic coda waves to invert for model-derived yield and depth of burial. The only requirements of this method are that the region of interest must contain several (> 3) seismic sources which have been well recorded (SNR > 2) by a set of shared stations. For regions with only a few events (~ 3) and/or extremely band-limited observations for which the depth of burial is unknown, an earthquake with a known magnitude may be required to obtain an accurate and robust yield estimate.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS2010

Assessment of Current Estimates of Global and Regional Mean Sea Level from the TOPEX/Poseidon, Jason-1, and OSTM 17-Year Record

The science value of satellite altimeter observations has grown dramatically over time as enabling models and technologies have increased the value of data acquired on both past and present missions. With the prospect of an observational time series extending into several decades from TOPEX/Poseidon through Jason-1 and the Ocean Surface Topography Mission (OSTM), and further in time with a future set of operational altimeters, researchers are pushing the bounds of current technology and modeling capability in order to monitor global sea level rate at an accuracy of a few tenths of a mm/yr. The measurement of mean sea-level change from satellite altimetry requires an extreme stability of the altimeter measurement system since the signal being measured is at the level of a few mm/yr. This means that the orbit and reference frame within which the altimeter measurements are situated, and the associated altimeter corrections, must be stable and accurate enough to permit a robust MSL estimate. Foremost, orbit quality and consistency are critical to satellite altimeter measurement accuracy. The orbit defines the altimeter reference frame, and orbit error directly affects the altimeter measurement. Orbit error remains a major component in the error budget of all past and present altimeter missions. For example, inconsistencies in the International Terrestrial Reference Frame (ITRF) used to produce the precision orbits at different times cause systematic inconsistencies to appear in the multimission time-frame between TOPEX and Jason-1, and can affect the intermission calibration of these data. In an effort to adhere to cross mission consistency, we have generated the full time series of orbits for TOPEX/Poseidon (TP), Jason-1, and OSTM based on recent improvements in the satellite force models, reference systems, and modeling strategies. The recent release of the entire revised Jason-1 Geophysical Data Records, and recalibration of the microwave radiometer correction also require the further re-examination of inter-mission consistency issues. Here we present an assessment of these recent improvements to the accuracy of the 17 -year sea surface height time series, and evaluate the subsequent impact on global and regional mean sea level estimates.

Beckley, Brian D.↗
DOE OSTI2023

Indoor Occupancy Sensing via Networked Nodes (2012–2022): A Review

In the past decade, different sensing mechanisms and algorithms have been developed to detect or estimate indoor occupancy. One of the most recent advancements is using networked sensor nodes to create a more comprehensive occupancy detection system where multiple sensors can identify human presence within more expansive areas while delivering enhanced accuracy compared to a system that relies on stand-alone sensor nodes. The present work reviews the studies from 2012 to 2022 that use networked sensor nodes to detect indoor occupancy, focusing on PIR-based sensors. Methods are compared based on pivotal ADPs that play a significant role in selecting an occupancy detection system for applications such as Health and Safety or occupant comfort. These parameters include accuracy, information requirement, maximum sensor failure and minimum observation rate, and feasible detection area. We briefly describe the overview of occupancy detection criteria used by each study and introduce a metric called “sensor node deployment density” through our analysis. This metric captures the strength of network-level data filtering and fusion algorithms found in the literature. It is hinged on the fact that a robust occupancy estimation algorithm requires a minimal number of nodes to estimate occupancy. This review only focuses on the occupancy estimation models for networked sensor nodes. It thus provides a standardized insight into networked nodes’ occupancy sensing pipelines, which employ data fusion strategies, network-level machine learning algorithms, and occupancy estimation algorithms. This review thus helps determine the suitability of the reviewed methods to a standard set of application areas by analyzing their gaps.

Emad-Ud-Din, Muhammad (ORCID:0000000279515538)↗
DOE OSTI2024

Specifying angular reference for robust three-phase state estimation

While the phase angle of any of the bus voltages can be chosen as the angular reference in the state estimation formulation of positive sequence networks, the same approach does not readily extend to three-phase network state estimation problem. It is commonly assumed that there is at least one bus where the bus voltages are perfectly balanced with phase angles displaced ±120° and these balanced three phase voltages are used as the three-phase reference in solving the three-phase state estimation problem. This assumption may be quite realistic in transmission networks, and for distribution networks with a strong transmission system connection. However, it might not be realistic to assume existence of a perfectly balanced reference bus in today’s distribution systems with ever increasing penetration of renewable sources or for isolated operation of microgrids. In this paper, a novel state estimation formulation will be presented which facilitates correct solution irrespective of the existence of buses with perfectly balanced voltages. The new formulation is general, and lends itself to bad data processing. It yields accurate results in any three-phase power system irrespective of its operating conditions (balanced or highly unbalanced), configuration (isolated microgrid, connected to transmission system, etc.) and whether or not it contains any synchronous generators. The performance of the method is validated using the IEEE 123 bus three-phase system.

24 POWER TRANSMISSION AND DISTRIBUTION↗