Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “maximum likelihood 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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Small Body GN and C Research Report: G-SAMPLE - An In-Flight Dynamical Method for Identifying Sample Mass [External Release Version]

G-SAMPLE is an in-flight dynamical method for use by sample collection missions to identify the presence and quantity of collected sample material. The G-SAMPLE method implements a maximum-likelihood estimator to identify the collected sample mass, based on onboard force sensor measurements, thruster firings, and a dynamics model of the spacecraft. With G-SAMPLE, sample mass identification becomes a computation rather than an extra hardware requirement; the added cost of cameras or other sensors for sample mass detection is avoided. Realistic simulation examples are provided for a spacecraft configuration with a sample collection device mounted on the end of an extended boom. In one representative example, a 1000 gram sample mass is estimated to within 110 grams (95% confidence) under realistic assumptions of thruster profile error, spacecraft parameter uncertainty, and sensor noise. For convenience to future mission design, an overall sample-mass estimation error budget is developed to approximate the effect of model uncertainty, sensor noise, data rate, and thrust profile error on the expected estimate of collected sample mass.

in-flight sample-mass identification↗
NASA NTRS2012

Identification of Time-Varying Pilot Control Behavior in Multi-Axis Control Tasks

Recent developments in fly-by-wire control architectures for rotorcraft have introduced new interest in the identification of time-varying pilot control behavior in multi-axis control tasks. In this paper a maximum likelihood estimation method is used to estimate the parameters of a pilot model with time-dependent sigmoid functions to characterize time-varying human control behavior. An experiment was performed by 9 general aviation pilots who had to perform a simultaneous roll and pitch control task with time-varying aircraft dynamics. In 8 different conditions, the axis containing the time-varying dynamics and the growth factor of the dynamics were varied, allowing for an analysis of the performance of the estimation method when estimating time-dependent parameter functions. In addition, a detailed analysis of pilots adaptation to the time-varying aircraft dynamics in both the roll and pitch axes could be performed. Pilot control behavior in both axes was significantly affected by the time-varying aircraft dynamics in roll and pitch, and by the growth factor. The main effect was found in the axis that contained the time-varying dynamics. However, pilot control behavior also changed over time in the axis not containing the time-varying aircraft dynamics. This indicates that some cross coupling exists in the perception and control processes between the roll and pitch axes.

Zaal, Peter M. T.↗
NASA NTRS2014

How Well Does the Latest Anthropomorphic Test Device Mimic Human Impact Responses?

One of the goals of the NASA Occupant Protection Group is to understand the human tolerance to dynamic loading. This knowledge has to come through indirect approaches such as existing human response databases, anthropometric test devices (ATD), animal testing, post‐mortem human subjects, and models. This study investigated the biofidelity of the National Highway Traffic Safety Administration's ATD named the THOR (test device for human occupant restraint). If THOR responds comparably to humans, then it could potentially be used as a human surrogate to help validate space vehicle requirements for occupant protection. The THOR responses to frontal and spinal impacts (ranging from 8 to 12 G with rise times of 40, 70, and 100 ms) were measured and compared to human volunteer responses (95 trials in frontal and 58 in spinal) previously collected by the U. S. Air Force on the same horizontal impact accelerator. The impact acceleration profiles tested are within the expected range of multi‐purpose crew vehicle (MPCV) landing dynamics. A correlation score was calculated for each THOR to human comparison using CORA (CORrelation and Analysis) software. A two‐parameter beta distribution model fit was obtained for each dependent variable using maximum likelihood estimation. For frontal impacts, the THOR head x‐acceleration peak response correlated with the human response at 8 and 10‐G 100 ms but not 10‐G 70 ms. The phase lagged the human response. Head z‐acceleration was not correlated. Chest x‐acceleration was in phase, had a higher peak response, and was well correlated with lighter subjects (Cora = 0.8 for 46 kg vs. Cora = 0.4 for 126 kg). Head x‐displacement had a leading phase. Several subjects responded with the same peak displacement but the mean of the group was lower. The shoulder x‐displacement was in phase but had higher peaks than the human response. For spinal impacts, the THOR head x‐acceleration was not well correlated. Head and chest z‐acceleration was in phase but had a higher peak response. Chest z‐acceleration was highly correlated with heavier subjects at lower G pulses (Cora = 0.86 for 125 kg at 8 G). The human response was variable in shoulder z‐displacement but the THOR was in phase and was comparable to the mean peak response. Head xand z‐displacement was in phase but had higher peaks. Seat pan forces were well correlated, were in phase, but had a larger peak response than most subjects. The THOR does not respond to frontal and spinal impacts exactly the same way that a human does. Some responses are well matched and others are not. Understanding the strengths and weaknesses of this ATD is an important first step in determining its usefulness in occupant protection at NASA

Newby, N.↗
NASA NTRS2014

How Well Does the Latest Anthropomorphic Test Device Mimic Human Impact Responses?

One of the goals of the NASA Occupant Protection Group is to understand the human tolerance to dynamic loading. This knowledge has to come through indirect approaches such as existing human response databases, anthropometric test devices (ATD), animal testing, post-­‐mortem human subjects, and models. This study investigated the biofidelity of the National Highway Traffic Safety Administration's ATD named the THOR (test device for human occupant restraint). If THOR responds comparably to humans, then it could potentially be used as a human surrogate to help validate space vehicle requirements for occupant protection. The THOR responses to frontal and spinal impacts (ranging from 8 to 12 G with rise times of 40, 70, and 100 ms) were measured and compared to human volunteer responses (95 trials in frontal and 58 in spinal) previously collected by the U. S. Air Force on the same horizontal impact accelerator. The impact acceleration profiles tested are within the expected range of multi-­‐purpose crew vehicle (MPCV) landing dynamics. A correlation score was calculated for each THOR to human comparison using CORA (CORrelation and Analysis) software. A two-­‐parameter beta distribution model fit was obtained for each dependent variable using maximum likelihood estimation. For frontal impacts, the THOR head x-­‐acceleration peak response correlated with the human response at 8 and 10-­‐G 100 ms but not 10-­‐G 70 ms. The phase lagged the human response. Head z-­‐acceleration was not correlated. Chest x-­‐acceleration was in phase, had a higher peak response, and was well correlated with lighter subjects (Cora = 0.8 for 46 kg vs. Cora = 0.4 for 126 kg). Head x-­‐displacement had a leading phase. Several subjects responded with the same peak displacement but the mean of the group was lower. The shoulder x-­‐displacement was in phase but had higher peaks than the human response. For spinal impacts, the THOR head x-­‐acceleration was not well correlated. Head and chest z-­‐acceleration was in phase but had a higher peak response. Chest z-­‐acceleration was highly correlated with heavier subjects at lower G pulses (Cora = 0.86 for 125 kg at 8 G). The human response was variable in shoulder z-­‐displacement but the THOR was in phase and was comparable to the mean peak response. Head x-­‐ and z-­‐displacement was in phase but had higher peaks. Seat pan forces were well correlated, were in phase, but had a larger peak response than most subjects. The THOR does not respond to frontal and spinal impacts exactly the same way that a human does. Some responses are well matched and others are not. Understanding the strengths and weaknesses of this ATD is an important first step in determining its usefulness in occupant protection at NASA

Newby, Nate↗
NASA NTRS2015

Validation of Rain Rate Retrievals for the Airborne Hurricane Imaging Radiometer (HIRAD)

The NASA Hurricane and Severe Storm Sentinel (HS3) mission is an aircraft field measurements program using NASA's unmanned Global Hawk aircraft system for remote sensing and in situ observations of Atlantic and Caribbean Sea hurricanes. One of the principal microwave instruments is the Hurricane Imaging Radiometer (HIRAD), which measures surface wind speeds and rain rates. For validation of the HIRAD wind speed measurement in hurricanes, there exists a comprehensive set of comparisons with the Stepped Frequency Microwave Radiometer (SFMR) with in situ GPS dropwindsondes [1]. However, for rain rate measurements, there are only indirect correlations with rain imagery from other HS3 remote sensors (e.g., the dual-frequency Ka- & Ku-band doppler radar, HIWRAP), which is only qualitative in nature. However, this paper presents results from an unplanned rain rate measurement validation opportunity that occurred in 2013, when HIRAD flew over an intense tropical squall line that was simultaneously observed by the Tampa NEXRAD meteorological radar (Fig. 1). During this experiment, Global Hawk flying at an altitude of 18 km made 3 passes over the rapidly propagating thunderstorm, while the TAMPA NEXRAD perform volume scans on a 5-minute interval. Using the well-documented NEXRAD Z-R relationship, 2D images of rain rate (mm/hr) were obtained at two altitudes (3 km & 6 km), which serve as surface truth for the HIRAD rain rate retrievals. A preliminary comparison of HIRAD rain rate retrievals (image) for the first pass and the corresponding closest NEXRAD rain image is presented in Fig. 2 & 3. This paper describes the HIRAD instrument, which 1D synthetic-aperture thinned array radiometer (STAR) developed by NASA Marshall Space Flight Center [2]. The rain rate retrieval algorithm, developed by Amarin et al. [3], is based on the maximum likelihood estimation (MLE) technique, which compares the observed Tb's at the HIRAD operating frequencies of 4, 5, 6 and 6.6 GHz with corresponding theoretical Tb values from a forward radiative transfer model (RTM). The optimum solution is the integrated rain rate that minimizes the difference between RTM and observed values. Because the excess Tb from rain comes from the direct upwelling and the indirect reflected downwelling paths through the atmosphere, there are several assumptions made for the 2D rain distribution in the antenna incident plane (crosstrack to flight direction). The opportunity to knowing 2D rain surface truth from NEXRAD at two different altitudes will enable a comprehensive evaluation to be preformed and reported in this paper.

hurricane↗
NASA NTRS2017

Trajectory Dispersed Vehicle Process for Space Launch System

The Space Launch System (SLS) vehicle is part of NASA's deep space exploration plans that includes manned missions to Mars. Manufacturing uncertainties in design parameters are key considerations throughout SLS development as they have significant effects on focus parameters such as lift-off-thrust-to-weight, vehicle payload, maximum dynamic pressure, and compression loads. This presentation discusses how the SLS program captures these uncertainties by utilizing a 3 degree of freedom (DOF) process called Trajectory Dispersed (TD) analysis. This analysis biases nominal trajectories to identify extremes in the design parameters for various potential SLS configurations and missions. This process utilizes a Design of Experiments (DOE) and response surface methodologies (RSM) to statistically sample uncertainties, and develop resulting vehicles using a Maximum Likelihood Estimate (MLE) process for targeting uncertainties bias. These vehicles represent various missions and configurations which are used as key inputs into a variety of analyses in the SLS design process, including 6 DOF dispersions, separation clearances, and engine out failure studies.

Statham, Tamara↗
NASA NTRS2018

Modeling of Time-Dependent Strength Degradation of SiC/SiC Ceramic Matrix Composites via MAC/GMC Computer Code

Silicon carbide fiber reinforced silicon carbide (SiC/SiC) ceramic matrix composites (CMCs) display time-dependent strength degradation at intermediate temperatures (600 to 900 degrees Centigrade). This is generally believed to be an oxidation induced phenomenon. The understanding of the effect of temperature with environment (oxidation) is key towards development of SiC/SiC CMCs with a reliable load carrying capacity. Various theories have been proposed to explain the strength degradation. One suggests that the boron nitride (BN) coating deposited on the fibers oxidizes causing fusion of fibers. Another theory proposes that the SiC fibers are oxidized forming a silica scale leading to premature fiber failure. A more recent theory suggests that SiC fiber strength is intrinsically time-dependent due to slow crack growth of flaws in the fibers. An empirical model, termed as a "fiber classic model," which is based on a standard slow crack growth type power-law, has been implemented within NASA's micromechanics-based MAC/GMC computer code as a user routine. Model parameters for this "classic model" were calibrated from stress-rupture data of Hi-Nicalon (TM) monofilaments using the maximum likelihood estimation (MLE) technique. This new capability in the MAC/GMC computer code was then used to predict the stress-rupture behavior of Hi-Nicalon (TM) tows as well as 2-D SiC/SiC composites reinforced with Hi-Nicalon (TM) fibers. Results demonstrate that the MAC/GMC with this new capability successfully predicts the time-to-failure vs. applied stress within the intermediate temperature range at various scales as well as laminated composites in an oxidizing environment.

Environmental degradation↗
NASA NTRSDate not supplied

Reduced-Order Aerodynamic Modeling Based on CFD Frequency Responses from Multisine Inputs

A system identification analysis was performed to determine reduced-order models of a computational fluid dynamics (CFD) solver for linear aeroelastic analysis and control design. The application was to the FUN3D code and the flexible half-span wind tunnel test article, in transonic flow conditions, used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. Multiple inputs (structural mode displacements and control surface deflections) were simultaneously excited with orthogonal phase-optimized multisines and multiple outputs (generalized aerodynamic forces) were recorded, from which the matrix of frequency responses were computed using a single CFD run. A state-space model was then fit to the frequency response data using a maximum-likelihood estimator.

System identification↗
NASA NTRSDate not supplied

An Interactive MATLAB Program for Fitting Transfer Functions to Frequency Responses

A computer program called FRFit (Frequency Response Fitting) is described for fitting single-input single-output transfer function models to empirical frequency response data. The program is interactive in that the user specifies ``elementary factors'' (gain, delay, pure differentiators and integrators, and first- and second-order zeros and poles) by entering numerical values or moving sliders in a graphical user interface. A nonlinear optimization can then be performed to obtain maximum likelihood estimates of transfer function parameters and uncertainties to provide feedback on the modeling and refine estimates. Several examples are discussed, including the identification of aircraft pitch dynamics from simulation data and data reported in the literature, approximating Theodorsen's function of unsteady aerodynamics, and obtaining a reduced-order model of a computational fluid dynamics code describing the unsteady aerodynamics around an aeroelastic wing. The program has some usefulness as a teaching aid, and can be applied to model structure determination, reduced-order modeling, preliminary analysis, and simple system identification problems. The program was written in MATLAB and is planned for public release through the NASA Software Catalog.

System identification↗
NASA NTRSDate not supplied

An Interactive MATLAB Program for Fitting Transfer Functions to Frequency Responses

A computer program called FRFit (Frequency Response Fitting) for matching single-input single-output (SISO) transfer function models to empirical frequency response data is described. The program was written in MATLAB and has a graphical user interface (GUI). It is interactive in that the user manually builds the transfer function model using ``elementary factors'' (gain, delay, differentiators and integrators, and first- and second-order poles and zeros) and adjusts their values with sliders or entry fields. A nonlinear optimization can also be used to determine maximum-likelihood estimates of the transfer function parameters and their associated uncertainties. The program has some usefulness as a teaching aid, and can be applied to model structure determination, reduced-order modeling, preliminary analysis, and other system identification problems. FRFit is demonstrated using example problems, including the identification of aircraft transfer functions and rational function approximations of Theodorsen's function.

Frequency response↗
NASA NTRS2020

Synergistic Use of Hyperspectral UV-Visible OMI and Broadband Meteorological Imager MODIS Data for a Merged Aerosol Product

The retrieval of optimal aerosol datasets by the synergistic use of hyperspectral ultraviolet(UV)–visible and broadband meteorological imager (MI) techniques was investigated. The Aura Ozone Monitoring Instrument (OMI) Level 1B (L1B) was used as a proxy for hyperspectral UV–visible instrument data to which the Geostationary Environment Monitoring Spectrometer (GEMS) aerosol algorithm was applied. Moderate-Resolution Imaging Spectroradiometer (MODIS) L1B and dark target aerosol Level 2 (L2) data were used with a broadband MI to take advantage of the consistent time gap between the MODIS and the OMI. First, the use of cloud mask information from the MI infrared (IR) channel was tested for synergy. High-spatial-resolution and IR channels of the MI helped mask cirrus and sub-pixel cloud contamination of GEMS aerosol, as clearly seen in aerosol optical depth (AOD) validation with Aerosol Robotic Network (AERONET) data. Second, dust aerosols were distinguished in the GEMS aerosol-type classification algorithm by calculating the total dust confidence index (TDCI) from MODIS L1B IR channels. Statistical analysis indicates that the Probability of Correct Detection (POCD) between the forward and inversion aerosol dust models (DS) was increased from 72% to 94% by use of the TDCI for GEMS aerosol-type classification, and updated aerosol types were then applied to the GEMS algorithm. Use of the TDCI for DS type classification in the GEMS retrieval procedure gave improved single-scattering albedo (SSA) values for absorbing fine pollution particles (BC) and DS aerosols. Aerosol layer height (ALH) retrieved from GEMS was compared with Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data, which provides high-resolution vertical aerosol profile information. The CALIOP ALH was calculated from total attenuated backscatter data at 1064 nm, which is identical to the definition of GEMS ALH. Application of the TDCI value reduced the median bias of GEMS ALH data slightly. The GEMS ALH bias approximates zero, especially for GEMS AOD values of>~0.4 and GEMS SSA values of<~0.95.Finally, the AOD products from the GEMS algorithm and MI were used in aerosol merging with the maximum-likelihood estimation method, based on a weighting factor derived from the standard deviation of the original AOD products. With the advantage of the UV–visible channel in retrieving aerosol properties over bright surfaces, the combined AOD products demonstrated better spatial data availability than the original AOD products, with comparable accuracy. Furthermore, pixel-level error analysis of GEMS AOD data indicates improvement through MI synergy.

aerosol↗
NASA NTRS2001

Earth beacon acquisition and tracking under various illumination conditions

This paper develops an algorithm for determining the geocenter of the Earth regardless of the illumination by the sun using sub-pixel scanning and a simple thresholding technique. To complete this approach the acquisition algorithm is paired with a tracking technique based on maximum likelihood estimation.

beacon↗
NASA NTRSDate not supplied

Velocity and Temperature Measurements in High-speed Flows with Naturally Present Dust Particles Using Rayleigh and Mie Scattering

Dust particles and occasional moisture condensations are unavoidable reality of all wind tunnels. On the path to pursue a goal of velocity and temperature measurements in large transonic and supersonic wind tunnels we have created a tabletop, spectrally resolved, Rayleigh-Mie scattering setup around a small jet fed by ambient and lightly seeded air to determine the viability and accuracy of the technique. The other reality of a wind tunnel setup is the background scattering or the glare at laser frequency, which contaminates the Rayleigh-Mie scattered light. This is simulated by backgrounds with different reflectivity towards the collection optics. Light from a CW laser is delivered via an optical fiber and the scattered light is spectrally resolved using a stabilized Fabry-Perot interferometer, followed by imaging on an EMCCD camera. A model of the of the combined background glare, Mie scattering, and the Rayleigh spectrum was fitted to the camera image using maximum likelihood estimation. Since the background glare occurs at the known frequency of the incident light and the Mie scattering peak corresponds to the Doppler shift from the bulk velocity, both were easily identified, and provided a measure of flow velocity. Preliminary results show that the Rayleigh spectrum can also be resolved, which provides a measure of temperature. It is observed that a slight drift of the laser frequency during data collection affected fitting of the model function leading to larger error. A feedback loop-based stabilization system is on development to take advantage of slight tunability of the laser via a piezo-control. Preliminary results are presented in the abstract. More extensive data from a systematic survey will be presented in the final paper.

Rayleigh scattering↗
NASA NTRSDate not supplied

Velocity and Temperature Measurements in High-speed Flows with Naturally Present Dust Particles Using Rayleigh and Mie Scattering

Dust particles and occasional moisture condensations are unavoidable reality of all wind tunnels. On the path to pursue a goal of velocity and temperature measurements in high-speed wind tunnels we have created a tabletop, spectrally resolved, Rayleigh-Mie scattering setup around a small jet fed by ambient and air with different particle concentration to determine the accuracy of the technique. The other reality of a wind tunnel setup is the spurious reflection and scattering of the incident laser beam by solid surfaces (flare light), which contaminates the Rayleigh-Mie scattered light. This is simulated by backgrounds with different reflectivity towards the collection optics. Light from a CW laser is delivered via an optical fiber and the scattered light from a point on the laser path is spectrally resolved using a stabilized Fabry-Perot interferometer, followed by imaging on an EMCCD camera. A model of the combined background glare, Mie scattering, and the Rayleigh spectrum was fitted to the camera image using maximum likelihood estimation. It was observed that the present modeling approach can extract velocity and temperature with reasonable accuracy when the intensity of the Mie scattered light is less than or comparable to that of the Rayleigh scattered light; beyond that error in velocity measurement remains reasonable ±12m/s but the temperature measurement becomes progressively more inaccurate. Increasingly the flare light is found to cause a bias error in velocity. Similar trend is observed in the presence of both the flare light and the Mie scattered light. The setup has provided a set of data for the future improvement of the modeling procedure, and demonstrated that even in the case of large amount of dust particles and large flare light the present technique can provide measurement of velocity with reasonable accuracy.

Rayleigh scattering↗
DOE OSTI2022

Testing multiflavored ULDM models with SPARC

We perform maximum likelihood estimates (MLEs) for single and double flavor ultralight dark matter (ULDM) models using the Spitzer Photometry and Accurate Rotation Curves (SPARC) database. These estimates are compared to MLEs for several commonly used cold dark matter (CDM) models. By comparing various CDM models we find, in agreement with previous studies, that the Burkert and Einasto models tend to perform better than other commonly used CDM models. We focus on comparisons between the Einasto and ULDM models and analyze cases for which the ULDM particle masses are: free to vary; and fixed. For each of these analyses, we perform fits assuming the soliton and halo profiles are: summed together; and matched at a given radius. When we let the particle masses vary, we find a negligible preference for any particular range of particle masses, within $10^{-25}\,\text{eV}\leq m\leq10^{-19}\,\text{eV}$, when assuming the summed models. For the matched models, however, we find that almost all galaxies prefer particles masses in the range $10^{-23}\,\text{eV}\lesssim m\lesssim10^{-20}\,\text{eV}$. For both double flavor models we find that most galaxies prefer approximately equal particle masses. We find that the summed models give much larger variances with respect to the soliton-halo (SH) relation than the matched models. When the particle masses are fixed, the matched models give median and mean soliton and halo values that fall within the SH relation bounds, for most masses scanned. When the particle masses are fixed in the fitting procedure, we find the best fit results for the particle mass $m=10^{-20.5}\,\text{eV}$ (for the single flavor models) and $m_1=10^{-20.5}\,\text{eV}$, $m_2=10^{-20.2}\,\text{eV}$ for the double flavor, matched model. We discuss how our study will be furthered using a reinforcement learning algorithm.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1982

A comparison of minimum distance and maximum likelihood techniques for proportion estimation

The estimation of mixing proportions P sub 1, P sub 2,...P sub m in the mixture density f(x) = the sum of the series P sub i F sub i(X) with i = 1 to M is often encountered in agricultural remote sensing problems in which case the p sub i's usually represent crop proportions. In these remote sensing applications, component densities f sub i(x) have typically been assumed to be normally distributed, and parameter estimation has been accomplished using maximum likelihood (ML) techniques. Minimum distance (MD) estimation is examined as an alternative to ML where, in this investigation, both procedures are based upon normal components. Results indicate that ML techniques are superior to MD when component distributions actually are normal, while MD estimation provides better estimates than ML under symmetric departures from normality. When component distributions are not symmetric, however, it is seen that neither of these normal based techniques provides satisfactory results.

Woodward, W. A.↗