Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonparametric”

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.

196 records · Page 11

Confidence Intervals for Laboratory Sonic Boom Annoyance Tests

Commercial supersonic flight is currently forbidden over land because sonic booms have historically caused unacceptable annoyance levels in overflown communities. NASA is providing data and expertise to noise regulators as they consider relaxing the ban for future quiet supersonic aircraft. One deliverable NASA will provide is a predictive model for indoor annoyance to aid in setting an acceptable quiet sonic boom threshold. A laboratory study was conducted to determine how indoor vibrations caused by sonic booms affect annoyance judgments. The test method required finding the point of subjective equality (PSE) between sonic boom signals that cause vibrations and signals not causing vibrations played at various amplitudes. This presentation focuses on a few statistical techniques for estimating the interval around the PSE. The techniques examined are the Delta Method, Parametric and Nonparametric Bootstrapping, and Bayesian Posterior Estimation.

Rathsam, Jonathan↗

Asymmetric Expansion of the Youngest Galactic Supernova Remnant G1.9+0.3

The youngest Galactic supernova remnant (SNR) G1.9+0.3, produced by a (probable) SN Ia that exploded approximately 1900 CE, is strongly asymmetric at radio wavelengths, much brighter in the north, but bilaterally symmetric in X-rays. We present the results of X-ray expansion measurements that illuminate the origin of the radio asymmetry. We confirm the mean expansion rate (2011-2015) of 0.58% per yr, but large spatial variations are present. Using the nonparametric 'Demons' method, we measure the velocity field throughout the entire SNR, finding that motions vary by a factor of 5, from 0.''09 to 0.''44 per yr. The slowest shocks are at the outer boundary of the bright northern radio rim, with velocities v(sub s) as low as 3600 km per sec (for an assumed distance of 8.5 kpc), much less than v(sub s) = 12,000-13,000 km per sec along the X-ray-bright major axis. Such strong deceleration of the northern blast wave most likely arises from the collision of SN ejecta with a much denser than average ambient medium there. This asymmetric ambient medium naturally explains the radio asymmetry. In several locations, significant morphological changes and strongly nonradial motions are apparent. The spatially integrated X-ray flux continues to increase with time. Based on Chandra observations spanning 8.3 yr, we measure its increase at 1.3% +/- 0.8% per yr. The SN ejecta are likely colliding with the asymmetric circumstellar medium ejected by the SN progenitor prior to its explosion.

ISM: individual objects (G1.9+0.3) ISM: individual↗

A Comparative Analysis of Micrometeorological Determinants of Evapotranspiration Rates Within a Heterogeneous Urban Environment

Variability in micrometeorological conditions and their influence on estimated reference evapotranspiration (RET) rates were evaluated across a heterogeneous urban environment. Micrometeorological data sets (incoming solar radiation, air temperature, relative humidity and wind speed) were collected over a one-year period at six weather stations in New York City, NY (USA). Weather stations are located at four new urban green space monitoring sites and two airports. Reference evapotranspiration (RET) rates were estimated from the micrometeorological data sets for a short reference surface at a daily time-step using the ASCE Standardized Reference Evapotranspiration Equation, a Penman-Monteith based combination equation. Nonparametric comparative statistical analyses (Kruskal-Wallis) revealed statistically significant differences (at significance level α = 0.05) in micrometeorological conditions and estimated RET rates between the six sites. On a cumulative annual basis, estimated RET varied by up to 40 percent between the sites. A new technique for adjusting weather data collected at one location (e.g. regional airports) for use at another location (e.g. interior engineered urban green spaces) was evaluated. The study highlights the importance, for accurate estimation of ET, of onsite micrometeorological data sets, but concludes that additional research is needed to more thoroughly characterize micrometeorological variability across heterogeneous urban environments, and also to evaluate the influence of non-meteorological determinants, e.g. vegetation type, soil/media type, media moisture conditions and anthropogenic heat fluxes, on urban ET.

Urban environment↗

Practical Aspects of the Frequency Domain Approach for Aircraft System Identification

Practical aspects of the frequency-domain approach for aircraft system identification are explained and demonstrated. Topics related to experiment design, flight data analysis, and dynamic modeling are included. For demonstration purposes, simulated time series data and simulated flight data from an F-16 nonlinear simulation with realistic noise are used. This approach enables detailed evaluations of the techniques and results, because the true characteristics of the data and aircraft dynamics are known for the simulated data. Analytical techniques and practical considerations are examined for the finite Fourier transform, nonparametric frequency response estimation, parametric modeling in the frequency domain, experiment design for frequency-domain modeling, data analysis and modeling in the frequency domain, and real-time calculations. Flight data from a subscale jet transport aircraft are used to demonstrate some of the techniques and technical issues.

Morelli, Eugene A.↗

NCA-LDAS: Overview and Analysis of Hydrologic Trends for the National Climate Assessment

Terrestrial hydrologic trends over the conterminous United States are estimated for 1980-47 2015 using the National Climate Assessment-Land Data Assimilation System (NCA-LDAS) reanalysis. NCA-LDAS employs the uncoupled Noah Version 3.3 land surface model at 0.125 degreesx0.125 degrees forced with NLDAS-2 meteorology, rescaled Climate Prediction Center precipitation, and assimilated satellite-based soil moisture, snow depth and irrigation products. Mean annual trends are reported using the nonparametric Mann-Kendall test at p<0.1 significance. Results illustrate the interrelationship between regional gradients in forcing trends, and trends in other land energy and water stores and fluxes. Mean precipitation trends range from +3 to +9 mm/yr in the Upper Great Plains and Northeast to -1 to -9 in the West and South; net radiation flux trends range from +0.05 to +0.20 W/m(exp2)/yr in the East to -0.05 to -0.20 in the West; U.S.-wide temperature trends average about +0.03K/yr. Trends in soil moisture, snow cover, latent and sensible heat fluxes and runoff are consistent with forcings, contributing to increasing evaporative fraction trends from West to East. Evaluation of NCA-LDAS trends compared to independent data indicates mixed results. The RMSE of U.S.-wide trends in number of snow cover days improved from 3.13 to 2.89 days/yr while trend detection increased 11%; Trends in latent heat flux were hardly affected, RMSE decreasing only 0.17 to 0.16 W/m(exp 2)/yr, while trend detection increased 2%; NCA-LDAS runoff trends degraded significantly from 2.6 to 16.1 mm/yr while trend detection was unaffected. Analysis also indicated that NCA-LDAS exhibits relatively more skill in low precipitation station density areas, suggesting there are limits to the effectiveness of satellite data assimilation in densely gaged regions. Overall, NCA-LDAS demonstrates capability for quantifying physically consistent, U.S. hydrologic climate trends over the satellite era.

Snow depth↗

New Approach for Temporal Stability Evaluation of Pseudo-Invariant Calibration Sites (PICS)

Pseudo-Invariant Calibration Sites (PICS) are one of the most popular methods for in-flight vicarious radiometric calibration of Earth remote sensing satellites. The fundamental question of PICS temporal stability has not been adequately addressed. However, the main purpose of this work is to evaluate the temporal stability of a few PICS using a new approach. The analysis was performed over six PICS (Libya 1, Libya 4, Niger 1, Niger 2, Egypt 1 and Sudan 1). The concept of a “Virtual Constellation” was developed to provide greater temporal coverage and also to overcome the dependence limitation of any specific characteristic derived from one particular sensor. TOA reflectance data from four sensors consistently demonstrating “stable” calibration to within 5%—the Landsat 7 ETM+ (Enhanced Thematic Mapper Plus), Landsat 8 OLI (Operational Land Imager), Terra MODIS (Moderate Resolution Imaging Spectroradiometer) and Sentinel-2A MSI (Multispectral Instrument)–were merged into a seamless dataset. Instead of using the traditional method of trend analysis (Student’s T test), a nonparametric Seasonal Mann-Kendall test was used for determining the PICS stability. The analysis results indicate that Libya 4 and Egypt 1 do not exhibit any monotonic trend in six reflective solar bands common to all of the studied sensors, indicating temporal stability. A decreasing monotonic trend was statistically detected in all bands, except SWIR 2, for Sudan 1 and the Green and Red bands for Niger 1. An increasing trend was detected in the Blue band for Niger 2 and the NIR band for Libya 1. These results do not suggest abandoning PICS as a viable calibration source. Rather, they indicate that PICS temporal stability cannot be assumed and should be regularly monitored as part of the sensor calibration process.

Tuli, Fatima Tuz Zafrin↗

Reliability Analysis in the Presence of Aleatory Uncertainty

This paper proposes a method for modeling a system’s response using data. In contrast to approaches that identify a limit state function, we focus on the case in which not all uncertain parameters affecting the response are observable and the measured response is corrupted by noise. To this end, the system response is not characterized by a limit state function but instead by a Random Predictor Model (RPM) having a nonparametric structure. Consequently, the resulting failure probability is not a scalar but a random variable. This variable accounts for the aleatory contributions of the model-form uncertainty and the measurement noise into the response. Furthermore, we propose a framework that enables trading off the predicted range of failure probabilities resulting from such an analysis with a measure of risk. In this context, risk is the percentage of all predicted outcomes the analyst is willing to ignore. The reliability analysis of an aeroelastic structure subject to flutter is used to illustrate the ideas proposed.

Crespo, L. G.↗

Aircraft System Identification from Multisine Inputs and Frequency Responses

The identification of aircraft flight dynamics is often performed using frequency responses, which are nonparametric models that quantify the steady-state magnitude and phase of a dynamic system response to sinusoidal inputs, as a function of frequency. Frequency responses are computed from measured input and output data, and then model parameters, such as stability and control derivatives, are estimated to best fit a parametric model to the empirical frequency response data. The utility of this approach is due to the familiarity of engineers with frequency responses, a number of theoretical and practical advantages under specific conditions, the availability of software packages, and many other reasons.

System identification↗

Threats to Resiliency of Redundant Systems Due to Destructive SEE

Destructive SEE pose serious challenges for the reliable use of COTS devices in space systems. We used system-level modeling to determine SEL rates that would likely compromise system reliability, resilience and capabilities. We then assembled a representative dataset of COTS CMOS parts and used nonparametric statistical techniques to assess the threat posed to redundant systems by destructive SEE.

single-event effects↗

NCA-LDAS: Overview and Analysis of Hydrologic Trends for the National Climate Assessment

Terrestrial hydrologic trends over the conterminous United States are estimated for 1980–2015 using the National Climate Assessment Land Data Assimilation System (NCA-LDAS) reanalysis. NCA-LDAS employs the uncoupled Noah version 3.3 land surface model at 0.125° 3 0.125° forced with NLDAS-2 meteorology, rescaled Climate Prediction Center precipitation, and assimilated satellite-based soil moisture, snow depth, and irrigation products. Mean annual trends are reported using the nonparametric Mann–Kendall test at p < 0.1 significance. Results illustrate the interrelationship between regional gradients in forcing trends and trends in other land energy and water stores and fluxes. Mean precipitation trends range from +3 to +9 mm/yr in the upper Great Plains and Northeast to -1 to -9 mm/yr in the West and South, net radiation flux trends range from +0.05 to +0.20 W/sq. m yr in the East to -0.05 to -0.20 W/sq. m yr in the West, and U.S.-wide temperature trends average about +0.03K/yr. Trends in soil moisture, snow cover, latent and sensible heat fluxes, and runoff are consistent with forcings, contributing to increasing evaporative fraction trends from west to east. Evaluation of NCA-LDAS trends compared to independent data indicates mixed results. The RMSE of U.S.-wide trends in number of snow cover days improved from 3.13 to 2.89 days/yr while trend detection increased 11%. Trends in latent heat flux were hardly affected, with RMSE decreasing only from 0.17 to 0.16 W/sq. m yr, while trend detection increased 2%. NCA-LDAS runoff trends degraded significantly from 2.6 to 16.1 mm/yr while trend detection was unaffected. Analysis also indicated that NCA-LDAS exhibits relatively more skill in low precipitation station density areas, suggesting there are limits to the effectiveness of satellite data assimilation in densely gauged regions. Overall, NCA-LDAS demonstrates capability for quantifying physically consistent, U.S. hydrologic climate trends over the satellite era.

Michael F. Jasinski↗

NASA Engineering and Safety Center Technical Bulletin No. 22-04: Uncertainty Quantification of Reduced Order Structural Dynamic Models

Uncertainty quantification (UQ) provides statistical bounds on prediction accuracy based on finite element model (FEM) uncertainty. An alternate method for UQ, called the Hybrid Parametric Variation (HPV) combines a parametric variation of the Hurty/Craig-Bampton (HCB) fixed-interface (FI) modal frequencies with a nonparametric variation (NPV) method. This provides a UQ method that can be traced to test data, which can be updated as additional data and improved correlated models become available.

Uncertainty Quantification↗

Uncertainty Models for the Hybrid Parametric Variation Method of Uncertainty Quantification; Analysis

There is some level of uncertainty in every finite element model (FEM), which flows to a level of uncertainty in predicted results. The purpose of uncertainty quantification (UQ) is to provide statistical bounds on prediction accuracy based on model uncertainty. This is distinct from model updating, which attempts to modify models to improve their accuracy. UQ does not improve the accuracy of models, but accepts that the models are inherently inaccurate and attempts to quantify the impact of that inaccuracy on predicted results. Previously, an alternate method for UQ, called the Hybrid Parametric Variation (HPV) method, was applied to Space Launch System (SLS) Hurty/Craig-Bampton (HCB) components to predict system-level statistics for launch vehicle attitude control transfer functions and core stage section loads due to buffet. The HPV method combines a parametric variation of the HCB fixed-interface (FI) modal frequencies with a nonparametric variation (NPV) method that randomly varies the HCB mass and stiffness matrices as Wishart random matrix distributions using random matrix theory (RMT). Alternatively, the most common method for modeling uncertainty in the structural dynamics community is a parametric approach, which varies physical parameters in the model. However, there are several disadvantages associated with the parametric method. Determining a reduced set of parameters that have a significant impact on the system response can be time consuming, and the selected parameter probability distributions are rarely reliably known. Therefore, in practice, the parameters are surrogates for the actual errors, and the link to parameter uncertainty is unknown. Another major drawback is that the uncertainty that can be represented is limited to the form of the nominal FEM. It is the experience of the authors that based on numerous aerospace programs, almost all FEM errors are in form rather than parameter values. This hypothesis is supported by the observation of the authors that it is almost never possible to ‘tune’ a FEM to match modal test results by only modifying model parameters. Model-form uncertainty cannot be directly represented by FEM input parameters nor included in a parametric approach. However, model-form uncertainty can be modeled using RMT, where a probability distribution is developed for the matrix ensemble of interest. The major advantage of the NPV method is that it covers errors in model form. The HPV method anchors uncertainty at the HCB component level to component modal test results by matching the HCB and test modes based on mode descriptions or other methods, and then applying differing levels of frequency variation. The specific variations depend on the confidence to which a component FEM has been validated through modal testing. The NPV method is layered on the frequency variation to match modal test self-orthogonality and cross-orthogonality (XO) results. Once the component uncertainty models are identified, they are assembled, and the uncertainty is propagated to the system level using a Monte Carlo (MC) analysis approach that generates statistics for system-level predictions This provides a UQ method that can be traced to test data, which can be updated as additional data and improved correlated models become available. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The second example is an application to a component that is design specific to the SLS. Based on this work and other assessments, the HPV method provides another tool to the toolset used for complex system UQ analysis. From experience gathered to date using the HPV method, additional design specific applications must be investigated to provide further confidence in the validity of the HPV method of UQ analysis.

Uncertainty quantification↗

Autoregressive distributed lag-based dynamic uniformity modeling and monitoring approaches for superconductor manufacturing

High-temperature superconductors (HTS), known for their high efficiency and low energy loss, have found profound applications across various fields, driving the demand for long, uniformly performing tapes. However, ensuring uniform performance over extended lengths of HTS tapes, often characterized by the consistency of critical current, remains challenging due to fluctuations in growth conditions during manufacturing. Here, to elucidate the mechanisms underlying variations in tape uniformity and enable real-time monitoring of associated parameters, we propose an Autoregressive Distributed Lag (ADL)-based Dynamic Uniformity Modeling and Monitoring (ADUM2) approach. This method integrates uniformity measurement, the identification of critical process parameters and real-time monitoring within the manufacturing process. The ADUM2 approach is applied to the advanced metal organic chemical vapor deposition (A-MOCVD) process, a pilot-scale method for superconductor manufacturing. Our model demonstrates superior performance compared to benchmark methods, accounting for over 80% of the total variance in the data and identifying 13 key process parameters influencing the uniformity of HTS tapes. This study offers significant insights into the high-temperature superconductor manufacturing process and holds the potential to facilitate the production of cost-effective, uniformly performing long superconducting tapes in the future.

autoregressive distributed lag analysis↗

Multiclass Classification Using Bayesian Multivariate Adaptive Regression Splines

We present a new Bayesian model for the problem of multiclass classification. In this model, the probabilities of class membership of a given observation are determined by the mean of a latent Gaussian distribution. The mean functions of this latent distribution consist of combinations of highly flexible basis functions of the inputs: multivariate adaptive regression splines (MARS), first developed for multiple regression. We use reversible jump Markov chain Monte Carlo to make inference on the classification model, including the number of basis functions. We compare the probabilistic classification performance of our proposed approach to existing methods on simulated and benchmark data, and compare uncertainty estimates on simulated data. Our proposed method compares favorably with existing Bayesian and frequentist multiclass classification methods in out-of-sample probabilistic classification, and uncertainty estimation of these probabilistic classifications. We examine the fit of the proposed method to a data set of hurricane storm surge levels near Delaware Bay, US, and conclude that sea level rise is a key contributor to damage delivered by storm surge.

97 MATHEMATICS AND COMPUTING↗

Lumbar spine disc heights and curvature: upright posture vs. supine compression harness

INTRODUCTION: Spinal lengthening in microgravity is thought to cause back pain in astronauts. A spinal compression harness can compress the spine to eliminate lengthening but the loading condition with harness is different than physiologic conditions. Our purpose was to compare the effect of spine compression with a harness in supine position on disk height and spinal curvature in the lumbar spine to that of upright position as measured using a vertically open magnetic resonance imaging system. METHODS: Fifteen healthy subjects volunteered. On day 1, each subject lay supine for an hour and a baseline scan of the lumbar spine was performed. After applying a load of fifty percent of body weight with the harness for thirty minutes, the lumbar spine was scanned again. On day 2, after a baseline scan, a follow up scan was performed after kneeling for thirty minutes within the gap between two vertically oriented magnetic coils. Anterior and posterior disk heights, posterior disk bulging, and spinal curvature were measured from the baseline and follow up scans. RESULTS: Anterior disk heights increased and posterior disk heights decreased compared with baseline scans both after spinal compression with harness and upright posture. The spinal curvature increased by both loading conditions of the spine. DISCUSSION: The spinal compression with specially designed harness has the same effect as the physiologic loading of the spine in the kneeling upright position. The harness shows some promise as a tool to increase the diagnostic capabilities of a conventional MR system.

Non-NASA Center↗

Angiogenesis and osteogenesis in an orthopedically expanded suture

The purpose of this study was to examine the angiogenic and the subsequent osteogenic responses during a 96-hour time-course after sutural expansion. Fifty rats were divided into: (1) a control group that received only angiogenic induction through injection of 5 ng/gm recombinant human endothelial cell growth factor (rhECGF); (2) an experimental group that received orthopedic expansion and rhECGF; (3) a sham group that received expansion and sodium chloride (NaCl) injection; and (4) a baseline group that received no expansion or injection. All rats were injected with 3H-thymidine (1.0 microCi/gm) 1 hour before death to label the DNA of S-phase cells. Demineralized sections (4 microm thick) were stained with hematoxylin and eosin. Angiogenesis and cell migration were analyzed with a previously established cell kinetics model. Analysis of variance was used to test the hypothesis that enhancement of angiogenesis stimulates reestablishment of osteogenic capability. Blood vessel number, area, and endothelial cell-labeled index significantly increased in experimental groups, but no difference was found between control and baseline groups. Labeled-pericyte index and activated pericyte numbers in the experimental group were also higher than in the sham groups. These results show that supplemental rhECGF enhances angiogenesis in expanded sutures but not in nonexpanded sutures. Data also suggest that pericytes are the source of osteoblasts in an orthopedically expanded suture.

NASA Discipline Musculoskeletal↗