Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “input 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 325 records · Page 18

Real-Time Estimation of Bare-Airframe Frequency Responses from Closed-Loop Data and Multisine Inputs

A method is presented for computing frequency responses of multiple-input multiple-output bare-airframe dynamics from flight test data containing feedback control and/or mixing of control effectors. Orthogonal phase-optimized multisines are used to simultaneously excite each input with unique harmonic frequencies, at which frequency responses are computed as ratios of output-to-input Fourier transform data. The confounding effects of feedback and mixing for frequency response estimation are resolved by interpolating the frequency responses among all the harmonic frequencies. The method can be run in batch for post-flight analysis, or in real time as the aircraft is flying. The effectiveness of the method was verified using closed-loop simulations of the subscale NASA T-2 generic transport airplane. The method was also demonstrated using flight test data from the X-56A MUTT aeroelastic airplane, which was flown with feedback control and mixing.

Jared A Grauer↗

Vibration-based health characterization of multiple IMUs in UAV applications

Inertial measurement units (IMUs) are vital in UAV navigation for vehicle attitude and position estimation, especially in the frequently-occurring case of temporary loss of GPS capabilities. Degradation in these units could yield inaccurate position estimates, resulting in incorrect input for vehicle control systems. A wide variety of these units exist to meet the needs of UAV operators, general aviation, and military vehicles. In recent decades, lightweight and inexpensive MEMS sensors have grown in popularity with respect to UAVs. Studies into the failure modes of MEMS devices suggest that vibration is one of the leading causes of MEMS degradation. MEMS sensors are generally encased in their own hermetic packaging, as their small size and sensitivity cause them to be particularly vulnerable to air damping and particulate contamination when exposed. Exposure to excessive vibration can cause cracking in this packaging as well as damage to connecting electronics, with recorded results that cause sensors to act outside of their specifications, potentially providing unreliable input to controllers. A prior study from the Politecnico di Torino looked into the degradation of MEMS IMUs, specifically the AXIS-AIS402, for aerospace and vibrating environments. In the study, the authors found that after exposing MEMS IMUs to simulated vibrations from working conditions of aerospace applications, the sensors showed increased noise and increased bias instability. At NASA Ames, several widely-used MEMS IMUs were selected for a new IMU degradation study, comparing accelerometer performance before and after exposure to similar levels of vibration simulating aerospace conditions. The IMUs were selected from the Pixhawk autopilot systems, hobby sUAS, or prior experimental studies at NASA. Additionally, the industrial grade VectorNav-100 IMU was selected as a ground truth reference point to which the performance of the test sensors can be compared. In the study, each sensor is exposed to functional levels of UAV vibration as replicated on a single-frequency vibratory table in the NASA Ames SHARP Laboratory, and each accelerometer is characterized for bias, drift, and noise characteristics both before and after vibratory exposure. The focus of this study is to monitor the health of the IMU sensors and determine if sensor degradation can be detected from output error and frequency spectrum analysis, rather than from microscopic examination of the sensors. In this case, sensor degradation and faulty data can be more easily detectable by non-experts hoping to use these sensors as an important piece of an integrated package. In addition, specific sensor degradation modes potentially leading to in-flight hazards could be identified, such as accelerated levels of drift and altered temperature response. The overall goal of this work is to increase levels of safety for future users of these sensors and to provide a clear analysis of sensor limitations currently lacking in most sensor specification documentation.

IMUs↗

Winter Arctic Sea Ice Thickness from ICESat-2 Freeboards

NASA’s ICESat-2 mission was launched in September 2018 with the primary goal of monitoring our rapidly changing polar regions. The sole instrument onboard, the Advanced Topographic Laser Altimeter System (ATLAS), is now providing routine, very high-resolution, surface elevation data across the globe, including the Arctic and Southern Oceans. In this study we demonstrate our new processing chain for converting the along-track ICESat-2 sea ice freeboard product (ATL10) into sea ice thickness, focusing our initial efforts on the Arctic Ocean. For this conversion we primarily make use of snow depth and density data from the NASA Eulerian Snow on Sea Ice Model (NESOSIM). The coarse resolution (~100 km) snow data are redistributed onto the high resolution (~30 - 100 m) ATL10 freeboards using relationships obtained from snow depth and freeboard data collected by NASA’s Operation IceBridge mission. We present regional sea ice thickness distributions and highlight their seasonal evolution through our first winter-season of data collection. We include ice thickness uncertainty estimates, while also acknowledging the limitations of these estimates. We generate a gridded monthly thickness product and compare this with various monthly sea ice thickness estimates obtained from ESA’s CryoSat-2 satellite mission, with ICE 33 Sat-2 showing consistently lower thicknesses. Finally, we compare our February/March 2019 thickness estimates to ICESat February/March (19th February – 21st March) 2008 ice thickness estimates using the same input assumptions, which show a ~0.40 m or ~21% thinning across an inner Arctic Ocean domain in this 11-year time period.

Alek A Petty↗

A Quantum Algorithm to Simulate Open Quantum Systems

Given the advent of quantum algorithms for a wide array of problems in linear algebra and machine learning, it is important to develop general methods for the simulation of arbitrary (ie non-unitary) operators on quantum hardware. In this talk, we present a novel quantum algorithm based on the quantum singular value transformation (QSVT) to apply an arbitrary operator K to some input state and subsequently estimate the expectation value of some observable. Our construction then immediately yields a route to estimating observables of states undergoing open quantum dynamics, whose effect is captured by a set of non-unitary Kraus operators. Our algorithm succeeds deterministically given the Sz-Nagy dilation, and we provide details on the algorithm's query and gate complexity, numerical verification, and comparisons with prior methods.

Quantum computing↗

The Containment Assurance Risk Framework of the Mars Sample Return Program

The Mars Sample Return campaign aims at bringing rock and atmospheric samples from Mars to Earth through a series of robotic missions. These missions would collect the samples being cached and deposited on Martian soil by the Perseverance rover, place them in a container, and launch them into Martian orbit for subsequent capture by an orbiter that would bring them back. Given there exists a non-zero probability that the samples contain biological material, precautions are being taken to design systems that would break the chain of contact between Mars and Earth. These include techniques such as sterilization of Martian particles, redundant containment vessels, and a robust reentry capsule capable of accurate landings without a parachute. Requirements exist that the probability of containment not assured of Martian-contaminated material into Earth’s biosphere be less than one in a million. To demonstrate compliance with this strict requirement, a statistical framework was developed to assess the likelihood of containment loss during each sample return phase and make a statement about the total combined mission probability of containment not assured. The work presented here describes this framework, which considers failure modes or fault conditions that can initiate failure sequences ultimately leading to containment not assured. Reliability estimates are generated from databases, design heritage, component specifications, or expert opinion in the form of probability density functions or point estimates and provided as inputs to the mathematical models that simulate the different failure sequences. The probabilistic outputs are then combined following the logic of several fault trees to compute the ultimate probability of containment not assured. Given the multidisciplinary nature of the problem and the different types of mathematical models used, the statistical tools needed for analysis are required to be computationally efficient. While standard Monte Carlo approaches are used for fast models, a multi-fidelity approach to rare event probabilities is proposed for expensive models. In this paradigm, inexpensive low-fidelity models are developed for computational acceleration purposes while the expensive high-fidelity model is kept in the loop to retain accuracy in the results. This work presents an example of end-to-end application of this framework highlighting the computational benefits of a multi-fidelity approach.

Giuseppe Cataldo↗

Probabilistic Hydrological Estimation of LandSlides (PHELS): Global Ensemble Landslide Hazard Modelling

In this study we present a model for the global Probabilistic Hydrological Estimation of LandSlides (PHELS). PHELS estimates the daily hazard of hydrologically triggered landslides at a coarse spatial resolution of 36 km, by combining landslide susceptibility (LSS) and (percentiles of) hydrological variable(s). The latter include daily rainfall, a 7-day antecedent rainfall index (ARI7) or root-zone soil moisture content (rzmc) as hydrological predictor variables, or the combination of rainfall and rzmc. The hazard estimates with any of these predictor variables have areas under the receiver operating characteristic curve (AUC) above 0.68. The best performance was found with combined rainfall and rzmc predictors (AUC = 0.79), which resulted in the lowest number of missed alarms (especially during spring) and false alarms. Furthermore, PHELS provides hazard uncertainty estimates by generating ensemble simulations based on repeated sampling of LSS and the hydrological predictor variables. The estimated hazard uncertainty follows the behaviour of the input variable uncertainties, is about 13.6 % of the estimated hazard value on average across the globe and in time and is smallest for very low and very high hazard values.

Anne Felsberg↗

A Performance Forecasting Framework for Concentrating Solar Power Systems

This presentation will summarize a methodology used to characterize the uncertainty in performance of a contrived CSP project. As well as develop a case study in which we backcast the projected and actual annual performance of an existing plant. Finally will highlight the importance of several assumptions and the use of time series inputs in obtaining realistic estimates.

concentrating solar power↗

Assessment of Antarctic Ice-Sheet Mass Balance Estimates: 1992 - 2009

Published mass balance estimates for the Antarctic Ice Sheet (AIS) lie between approximately +50 to -250 Gt/year for 1992 to 2009, which span a range equivalent to 15% of the annual mass input and 0.8 mm/year Sea Level Equivalent (SLE). Two estimates from radar-altimeter measurements of elevation change by European Remote-sensing Satellites (ERS) (+28 and -31 Gt/year) lie in the upper part, whereas estimates from the Input-minus-Output Method (IOM) and the Gravity Recovery and Climate Experiment (GRACE) lie in the lower part (-40 to -246 Gt/year). We compare the various estimates, discuss the methodology used, and critically assess the results. Although recent reports of large and accelerating rates of mass loss from GRACE=based studies cite agreement with IOM results, our evaluation does not support that conclusion. We find that the extrapolation used in the published IOM estimates for the 15 % of the periphery for which discharge velocities are not observed gives twice the rate of discharge per unit of associated ice-sheet area than the 85% faster-moving parts. Our calculations show that the published extrapolation overestimates the ice discharge by 282 Gt/yr compared to our assumption that the slower moving areas have 70% as much discharge per area as the faster moving parts. Also, published data on the time-series of discharge velocities and accumulation/precipitation do not support mass output increases or input decreases with time, respectively. Our modified IOM estimate, using the 70% discharge assumption and substituting input from a field-data compilation for input from an atmospheric model over 6% of area, gives a loss of only 13 Gt/year (versus 136 Gt/year) for the period around 2000. Two ERS-based estimates, our modified IOM, and a GRACE-based estimate for observations within 1992 to 2005 lie in a narrowed range of +27 to - 40 Gt/year, which is about 3% of the annual mass input and only 0.2 mm/year SLE. Our preferred estimate for 1992-2001 is - 47 Gt/year for West Antarctica, + 16 Gt/year for East Antarctica, and -31 Gt/year overall (+0.1 mm/year SLE), not including part of the Antarctic Peninsula (1.07 % of the AIS area)

Zwally, H. Jay↗

Mg II emission from shock waves in Mira variables

In order to investigate the atmospheric shock dynamics of Miras, a program to monitor the ultraviolet Mg II flux, the visual magnitude (using the IUE FES-magnitude) and whenever possible the hydrogen Balmer emission for several Miras was initiated. Results suggest that it should be possible to derive the density structure of the atmosphere and hence estimate quite accurately the kinetic energy input into the stellar atmosphere. The atmospheric shock models suggest that at phase equal 1.0 most of the energy which enters the region where the shock is located ends up in the wind. This is after the strongest Mg II emission is seen. From an observed velocity at infinity and an estimate of the energy input derived from the Mg II flux, it will be possible to determine the mass loss rate for these objects.

Brugel, Edward W.↗

SOC Microstructural Property Estimator

This pre-trained ML model is a tool that uses basic compositional parameters for porous solid oxide cell (SOC) electrodes - the phase fractions and mean particle/pore diameters – as inputs and uses them to estimate additional electrochemical performance parameters: active (i.e., connected) TPB density, all tortuosity factors, and phase pair specific interfacial areas. The electrode is assumed to be composed of two solid phases and a pore phase. The property calculations are performed using neural network regression models trained on a large bank of synthetic electrode microstructural data that NETL has generated using the program DREAM3D (that bank is also hosted on EDX: https://edx.netl.doe.gov/dataset/soc-synthetic-microstructure-bank). This means the generated parameters are based on training from actual measured properties from 3D microstructures, not estimated from geometric simplifications. This tool was developed and is intended to replace percolation theory calculations in models that use hypothetical electrode properties. An example use case would be running SOC performance simulations across a parametric sweep of electrode designs (e.g., varying phase fractions and particle sizes) and assessing how it impacts the electrochemical performance of the SOC. Within the parameter space of the training data (statistics of that parameter space is provided in the readme file), this model achieves sub-5% mean absolute percent errors, an order of magnitude less error than percolation theory across the same parameter space. However, be aware that this tool was developed with parametric simulations in mind, and users are encouraged to assess accuracy for their own specific use case rather than taking accuracy metrics at face value. More info, including a usage guide, is in the included readme file. This tool should be cited with the DOI number provided.

Electrode Microstructure↗

Space radiation shielding strategies and requirements for deep space missions

Methods for estimating crew exposure to radiation and for evaluating shield requirements for spacecraft equipment and crew are examined. The development status for deterministic space radiation transport computer codes and models of their nuclear interaction inputs, which are useful for estimating the composition and thickness of shield materials, is discussed. The relation between shield thickness and exposures is studied. Estimates of deep-space shield requirements are proposed. It is noted that an assessment of input parameter uncertainties reveals that a factor of 2 uncertainty in predicted dose equivalent could increase shield thickness by an order of magnitude of more.

Townsend, Lawrence W.↗

How Does the Use of Different Soil Mineralogical Atlases Impact Soluble Iron Deposition Estimates?

Atmospheric deposition is the primary input of soluble iron (Fe) in remote regions of the open ocean. Fe constitutes a fundamental micro-nutrient for marine biota, affecting biogeochemical cycles and ultimately carbon uptake. Fe is emitted to the atmosphere primarily from dust sources. Fe solubility at emission is origin-dependent, being mostly insoluble when associated with dust minerals. Earth System Models (ESMs) commonly assume that dust aerosols have a globally uniform composition, neglecting known regional variations in soil mineralogy. This work assesses the implications of soil composition uncertainties on the bio-available Fe delivery to the ocean by using a state-of-the-art ESM with a detailed atmospheric Fe cycle (EC-Earth3-Iron). We run two 1-year simulations considering two different Soil Mineralogy Atlases (SMAs). Our results reveal a non-negligible impact of soil mineralogy uncertainties on soluble Fe estimates. By only changing the source of information for the soil composition and updating the Fe content in minerals, we find differences of 31.5% in mineral Fe emissions and 26.7% in soluble Fe deposition budgets.

soluble iron↗

TLSLEAF: Automatic Leaf Angle Estimates From Single-Scan Terrestrial Laser Scanning

Leaf angle distribution (LAD) in forest canopies affects estimates of leaf area, light interception, and global-scale photosynthesis, but is often simplified to a single theoretical value. Here, we present TLSLeAF (Terrestrial Laser Scanning Leaf Angle Function), an automated open-source method of deriving LADs from terrestrial laser scanning. TLSLeAF produces canopy-scale leaf angle and LADs by relying on gridded laser scanning data. The approach increases processing speed, improves angle estimates, and requires minimal user input. Key features are automation, leaf–wood classification, beta parameter output, and implementation in R to increase accessibility for the ecology community. TLSLeAF precisely estimates leaf angle with minimal distance effects on angular estimates while rapidly producing LADs on a consumer-grade machine. We challenge the popular spherical LAD assumption, showing sensitivity to ecosystem type in plant area index and foliage profile estimates that translate to c. 25% and c. 11% increases in canopy net photosynthesis (c. 25%) and solar-induced chlorophyll fluorescence (c. 11%). TLSLeAF can now be applied to the vast catalog of laser scanning data already available from ecosystems around the globe. The ease of use will enable widespread adoption of the method outside of remote-sensing experts, allowing greater accessibility for addressing ecological hypotheses and large-scale ecosystem modeling efforts.

3D↗

A dual-input nonlinear system analysis of autonomic modulation of heart rate

Linear analyses of fluctuations in heart rate and other hemodynamic variables have been used to elucidate cardiovascular regulatory mechanisms. The role of nonlinear contributions to fluctuations in hemodynamic variables has not been fully explored. This paper presents a nonlinear system analysis of the effect of fluctuations in instantaneous lung volume (ILV) and arterial blood pressure (ABP) on heart rate (HR) fluctuations. To successfully employ a nonlinear analysis based on the Laguerre expansion technique (LET), we introduce an efficient procedure for broadening the spectral content of the ILV and ABP inputs to the model by adding white noise. Results from computer simulations demonstrate the effectiveness of broadening the spectral band of input signals to obtain consistent and stable kernel estimates with the use of the LET. Without broadening the band of the ILV and ABP inputs, the LET did not provide stable kernel estimates. Moreover, we extend the LET to the case of multiple inputs in order to accommodate the analysis of the combined effect of ILV and ABP effect on heart rate. Analyzes of data based on the second-order Volterra-Wiener model reveal an important contribution of the second-order kernels to the description of the effect of lung volume and arterial blood pressure on heart rate. Furthermore, physiological effects of the autonomic blocking agents propranolol and atropine on changes in the first- and second-order kernels are also discussed.

NASA Discipline Regulatory Physiology↗

Using Bayesian Estimation to Quantify the Risks of Spaceflight

Complex PRA is the best method and current gold standard, as it allows us to reflect changes (improvements) in individual spacecraft system risks over time. Yet, Bayesian methods – using simple inputs – allowed us to get estimates that are in line with high-complexity PRA. Here Bayesian estimation acts as a “gut check” on PRA (supposing a “reasonable” prior is used). Frequentist calculations can overestimate P(failure), especially when total number of trials are few (e.g., after Challenger). In general, Bayesian probability estimation can be used in any risk assessment situation to integrate what we think we know ahead of time with what we observe over time.

Reynolds, Robert J.↗

Technology Benefit Estimator (T/BEST): User's Manual

The Technology Benefit Estimator (T/BEST) system is a formal method to assess advanced technologies and quantify the benefit contributions for prioritization. T/BEST may be used to provide guidelines to identify and prioritize high payoff research areas, help manage research and limited resources, show the link between advanced concepts and the bottom line, i.e., accrued benefit and value, and to communicate credibly the benefits of research. The T/BEST software computer program is specifically designed to estimating benefits, and benefit sensitivities, of introducing new technologies into existing propulsion systems. Key engine cycle, structural, fluid, mission and cost analysis modules are used to provide a framework for interfacing with advanced technologies. An open-ended, modular approach is used to allow for modification and addition of both key and advanced technology modules. T/BEST has a hierarchical framework that yields varying levels of benefit estimation accuracy that are dependent on the degree of input detail available. This hierarchical feature permits rapid estimation of technology benefits even when the technology is at the conceptual stage. As knowledge of the technology details increases the accuracy of the benefit analysis increases. Included in T/BEST's framework are correlations developed from a statistical data base that is relied upon if there is insufficient information given in a particular area, e.g., fuel capacity or aircraft landing weight. Statistical predictions are not required if these data are specified in the mission requirements. The engine cycle, structural fluid, cost, noise, and emissions analyses interact with the default or user material and component libraries to yield estimates of specific global benefits: range, speed, thrust, capacity, component life, noise, emissions, specific fuel consumption, component and engine weights, pre-certification test, mission performance engine cost, direct operating cost, life cycle cost, manufacturing cost, development cost, risk, and development time. Currently, T/BEST operates on stand-alone or networked workstations, and uses a UNIX shell or script to control the operation of interfaced FORTRAN based analyses. T/BEST's interface structure works equally well with non-FORTRAN or mixed software analysis. This interface structure is designed to maintain the integrity of the expert's analyses by interfacing with expert's existing input and output files. Parameter input and output data (e.g., number of blades, hub diameters, etc.) are passed via T/BEST's neutral file, while copious data (e.g., finite element models, profiles, etc.) are passed via file pointers that point to the expert's analyses output files. In order to make the communications between the T/BEST's neutral file and attached analyses codes simple, only two software commands, PUT and GET, are required. This simplicity permits easy access to all input and output variables contained within the neutral file. Both public domain and proprietary analyses codes may be attached with a minimal amount of effort, while maintaining full data and analysis integrity, and security. T/BESt's sotware framework, status, beginner-to-expert operation, interface architecture, analysis module addition, and key analysis modules are discussed. Representative examples of T/BEST benefit analyses are shown.

Generazio, Edward R.↗

Estimating the Acute Health Impacts of Fire-Originated PM 2.5 Exposure During the 2017 California Wildfires: Sensitivity to Choices of Inputs

Exposure to wildfire smoke increases the risk of respiratory and cardiovascular hospital admissions. Health impact assessments, used to inform decision-making processes, characterize the health impacts of environmental exposures by combining preexisting epidemiological concentration–response functions (CRFs) with estimates of exposure. These two key inputs influence the magnitude and uncertainty of the health impacts estimated, but for wildfire-related impact assessments the extent of their impact is largely unknown. We first estimated the number of respiratory, cardiovascular, and asthma hospital admissions attributable to fire-originated PM 2.5 exposure in central California during the October 2017 wildfires, using Monte Carlo simulations to quantify uncertainty with respect to the exposure and epidemiological inputs. We next conducted sensitivity analyses, comparing four estimates of fire-originated PM 2.5 and two CRFs, wildfire and nonwildfire specific, to understand their impact on the estimation of excess admissions and sources of uncertainty. We estimate the fires accounted for an excess 240 (95% CI: 114, 404) respiratory, 68 (95% CI: -10, 159) cardiovascular, and 45 (95% CI: 18, 81) asthma hospital admissions, with 56% of admissions occurring in the Bay Area. Although differences between impact assessment methods are not statistically significant, the admissions estimates' magnitude is particularly sensitive to the CRF specified while the uncertainty is most sensitive to estimates of fire-originated PM 2.5 . Not accounting for the exposure surface's uncertainty leads to an underestimation of the uncertainty of the health impacts estimated. Employing context-specific CRFs and using accurate exposure estimates that combine multiple data sets generates more certain estimates of the acute health impacts of wildfires.

58 GEOSCIENCES↗

High-dimensional and permutation invariant anomaly detection

Methods for anomaly detection of new physics processes are often limited to low-dimensional spaces due to the difficulty of learning high-dimensional probability densities. Particularly at the constituent level, incorporating desirable properties such as permutation invariance and variable-length inputs becomes difficult within popular density estimation methods. In this work, we introduce a permutation-invariant density estimator for particle physics data based on diffusion models, specifically designed to handle variable-length inputs. We demonstrate the efficacy of our methodology by utilizing the learned density as a permutation-invariant anomaly detection score, effectively identifying jets with low likelihood under the background-only hypothesis. To validate our density estimation method, we investigate the ratio of learned densities and compare to those obtained by a supervised classification algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗