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 91 records · Page 5

Estimates of Lightning NOx Production Based on OMI NO2 Observations Over the Gulf of Mexico

We evaluate nitrogen oxide (NO(sub x) NO + NO2) production from lightning over the Gulf of Mexico region using data from the Ozone Monitoring Instrument (OMI) aboard NASAs Aura satellite along with detection efficiency-adjusted lightning data from the World Wide Lightning Location Network (WWLLN). A special algorithm was developed to retrieve the lightning NOx [(LNO(sub x)] signal from OMI. The algorithm in its general form takes the total slant column NO2 from OMI and removes the stratospheric contribution and tropospheric background and includes an air mass factor appropriate for the profile of lightning NO(sub x) to convert the slant column LNO2 to a vertical column of LNO(sub x). WWLLN flashes are totaled over a period of 3 h prior to OMI overpass, which is the time an air parcel is expected to remain in a 1 deg. x 1 deg. grid box. The analysis is conducted for grid cells containing flash counts greater than a threshold value of 3000 flashes that yields an expected LNO(sub x) signal greater than the background. Pixels with cloud radiance fraction greater than a criterion value (0.9) indicative of highly reflective clouds are used. Results for the summer seasons during 2007-2011 yield mean LNO(sub x) production of approximately 80 +/- 45 mol per flash over the region for the two analysis methods after accounting for biases and uncertainties in the estimation method. These results are consistent with literature estimates and more robust than many prior estimates due to the large number of storms considered but are sensitive to several substantial sources of uncertainty.

Lightning↗

Methods for in-flight robustness evaluation

The goal of this program was to combine modern control concepts with new identification techniques to develop a comprehensive package for estimation of 'robust flutter boundaries' based on experimental data. The goal was to use flight data, combined with a fundamental physical understanding of flutter dynamics, to generate a prediction of flutter speed and an estimate of the accuracy of the prediction. This report is organized as follows: the specific contributions of this project will be listed first. Then, the problem under study will be stated and the general approach will be outlined. Third, the specific system under study (F- 18 SRA) will be described and a preliminary data analysis will be performed. Then, the various steps of the flutter boundary determination will be outlined and applied to tile F-18 SRA data and others.

Brenner, Marty↗

Adjoint-Based, Three-Dimensional Error Prediction and Grid Adaptation

Engineering computational fluid dynamics (CFD) analysis and design applications focus on output functions (e.g., lift, drag). Errors in these output functions are generally unknown and conservatively accurate solutions may be computed. Computable error estimates can offer the possibility to minimize computational work for a prescribed error tolerance. Such an estimate can be computed by solving the flow equations and the linear adjoint problem for the functional of interest. The computational mesh can be modified to minimize the uncertainty of a computed error estimate. This robust mesh-adaptation procedure automatically terminates when the simulation is within a user specified error tolerance. This procedure for estimating and adapting to error in a functional is demonstrated for three-dimensional Euler problems. An adaptive mesh procedure that links to a Computer Aided Design (CAD) surface representation is demonstrated for wing, wing-body, and extruded high lift airfoil configurations. The error estimation and adaptation procedure yielded corrected functions that are as accurate as functions calculated on uniformly refined grids with ten times as many grid points.

Park, Michael A.↗

Stochastic robustness of linear 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 non-Gaussian cases, including uncertain-but-bounded variations. Trivial extensions of the procedure admit 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 also can be estimated. Results are particularly amenable to graphical presentation.

Stengel, Robert F.↗

Robust Modal Filtering and Control of the X-56A Model with Simulated Fiber Optic Sensor Failures

The X-56A aircraft is a remotely-piloted aircraft with flutter modes intentionally designed into the flight envelope. The X-56A program must demonstrate flight control while suppressing all unstable modes. A previous X-56A model study demonstrated a distributed-sensing-based active shape and active flutter suppression controller. The controller relies on an estimator which is sensitive to bias. This estimator is improved herein, and a real-time robust estimator is derived and demonstrated on 1530 fiber optic sensors. It is shown in simulation that the estimator can simultaneously reject 230 worst-case fiber optic sensor failures automatically. These sensor failures include locations with high leverage (or importance). To reduce the impact of leverage outliers, concentration based on a Mahalanobis trim criterion is introduced. A redescending M-estimator with Tukey bisquare weights is used to improve location and dispersion estimates within each concentration step in the presence of asymmetry (or leverage). A dynamic simulation is used to compare the concentrated robust estimator to a state-of-the-art real-time robust multivariate estimator. The estimators support a previously-derived mu-optimal shape controller. It is found that during the failure scenario, the concentrated modal estimator keeps the system stable.

Modal Filtering↗

Robust Modal Filtering and Control of the X-56A Model with Simulated Fiber Optic Sensor Failures

The X-56A aircraft is a remotely-piloted aircraft with flutter modes intentionally designed into the flight envelope. The X-56A program must demonstrate flight control while suppressing all unstable modes. A previous X-56A model study demonstrated a distributed-sensing-based active shape and active flutter suppression controller. The controller relies on an estimator which is sensitive to bias. This estimator is improved herein, and a real-time robust estimator is derived and demonstrated on 1530 fiber optic sensors. It is shown in simulation that the estimator can simultaneously reject 230 worst-case fiber optic sensor failures automatically. These sensor failures include locations with high leverage (or importance). To reduce the impact of leverage outliers, concentration based on a Mahalanobis trim criterion is introduced. A redescending M-estimator with Tukey bisquare weights is used to improve location and dispersion estimates within each concentration step in the presence of asymmetry (or leverage). A dynamic simulation is used to compare the concentrated robust estimator to a state-of-the-art real-time robust multivariate estimator. The estimators support a previously-derived mu-optimal shape controller. It is found that during the failure scenario, the concentrated modal estimator keeps the system stable.

Modal Filtering↗

Educational Experiences of Embry-Riddle Students through NASA Research Collaboration

NASA's educational. programs benefit students and faculty while increasing the overall productivity of the organization. The NASA Graduate Student Research Program (GSRP) awards fellowships for graduate study leading to both masters and doctoral degrees in several technical fields. GSRP participants have the option to utilize NASA Centers andlor university research facilities. In addition, GSRP students can serve as mentors for undergrad students to provide a truly unique learning experience. NASA's Cooperative Education Program allows undergraduate students the chance to gain "real-world" work experience in the field. It also gives NASA a no risk capability to evaluate the true performance of a prospective new hire without relying solely on a "paper resume" while providing the students with a greater hiring potential upon graduation, at NASA or elsewhere. University faculty can also benefit by participating in the NASA Faculty Fellowship Program (NFFP). This program gives the faculty an opportunity to work with NASA peers. The Mission Analysis Branch of the Expendable Launch Vehicles Division at NASA Kennedy Space Center has utilized these two programs with students from Embry-Riddle Aeronautical University (ERAU) to conduct research in modeling and developing a parameter estimation method for spacecraft fuel slosh using simple pendulum analogs. Simple pendulum models are used to understand complicated spacecraft fuel slosh behavior. A robust parameter estimation process will help to identif' the parameters that will predict the response fairly accurately during the initial stages of design. These programs provide students with a unique opportunity to work on "real-world" aerospace problems, like spacecraft fuel slosh,. This in turn reinforces their problem solving abilities and their communication skills such as interviewing, resume writing, technical writing, and presentation. Faculty benefits by applying what they have learned to the classroom. Through university collaborations with NASA and industry help students to acquire skills that are vital for their success upon entering the workforce.

Schlee, Keith↗

Educational Experiences of Embry-Riddle Students through NASA Research Collaboration

NASA's educational programs benefit students while increasing the overall productivity of the organization. The NASA Graduate Student Research Program (GSRP) awards fellowships for graduate study leading to both masters and doctoral degrees in several technical fields, while the Cooperative Education program allows undergraduate and graduate students the chance to gain work experience in the field. The Mission Analysis Branch of the Expendable Launch Vehicles Division at NASA Kennedy Space Center has utilized these two programs with students from Embry-Riddle Aeronautical University to conduct research in modeling and developing a parameter estimation method for spacecraft fuel slosh using simple pendulum analogs. Simple pendulum models are used to understand complicated spacecraft fuel slosh behavior. A robust parameter estimation process will help to identiFy the parameters that will predict the response fairly accurately during the initial stages of design. NASA's Cooperative Education Program trains the next wave of new hires while allowing graduate and undergraduate college students to gain valuable "real-world" work experience. It gives NASA a no risk capability to evaluate the true performance of a prospective new hire without relying solely on a paper resume, while providing the students with a greater hiring potential upon graduation, at NASA or elsewhere. In addition, graduate students serve as mentors for undergrad students and provide a unique learning environment. Providing students with a unique opportunity to work on "real-world" aerospace problems ultimately reinforces their problem solving abilities and their communication skills (in terms of interviewing, resume writing, technical writing, presentation, and peer review) that are vital for the workforce to succeed.

Schlee, Keith↗

Robustness Analysis and Reliable Flight Regime Estimation of an Integrated Resilent Control System for a Transport Aircraft

Formal robustness analysis of aircraft control upset prevention and recovery systems could play an important role in their validation and ultimate certification. As a part of the validation process, this paper describes an analysis method for determining a reliable flight regime in the flight envelope within which an integrated resilent control system can achieve the desired performance of tracking command signals and detecting additive faults in the presence of parameter uncertainty and unmodeled dynamics. To calculate a reliable flight regime, a structured singular value analysis method is applied to analyze the closed-loop system over the entire flight envelope. To use the structured singular value analysis method, a linear fractional transform (LFT) model of a transport aircraft longitudinal dynamics is developed over the flight envelope by using a preliminary LFT modeling software tool developed at the NASA Langley Research Center, which utilizes a matrix-based computational approach. The developed LFT model can capture original nonlinear dynamics over the flight envelope with the ! block which contains key varying parameters: angle of attack and velocity, and real parameter uncertainty: aerodynamic coefficient uncertainty and moment of inertia uncertainty. Using the developed LFT model and a formal robustness analysis method, a reliable flight regime is calculated for a transport aircraft closed-loop system.

Shin, Jong-Yeob↗

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↗

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.↗

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↗

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↗

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.↗

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.↗

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↗

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.↗

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.↗