SEARCH · Engineering Papers
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.
We develop single-event risk Prior probability distributions based on historical and heritage data. The Priors can be used to bound single-event effects risk for testing, part selection and design.
We investigate changes in daily precipitation extremes using Tropical Rainfall Measuring Mission (TRMM) data (1998-2013), which coincides with the "global warming hiatus." Results show a change in probability distribution functions of local precipitation events (LPEs) during this period consistent with previous global warming studies, indicating increasing contrast between wet and dry extremes, with more intense LPE, less moderate LPE, and more dry (no rain) days globally. Analyses for land and ocean separately reveal more complex and nuanced changes over land, characterized by a strong positive trend (+12.0% per decade, 99% confidence level (c.l.)) in frequency of extreme LPEs over the Northern Hemisphere extratropics during the wet season but a negative global trend (−6.6% per decade, 95% c.l.) during the dry season. A significant global drying trend (3.2% per decade, 99% c.l.) over land is also found during the dry season. Regions of pronounced increased dry events include western and central U.S., northeastern Asia, and Southern Europe/Mediterranean.
The mountain regions of the northeastern United States are a critical socioeconomic resource for Vermont, New York State, New Hampshire, Maine, and southern Quebec. While global climate models (GCMs) are important tools for climate change risk assessment at regional scales, even the increased spatial resolution of statistically downscaled GCMs (commonly approximately 1/ 8 deg) is not sufficient for hydrologic, ecologic, and land-use modeling of small watersheds within the mountainous Northeast. To address this limitation, an ensemble of topographically downscaled, high-resolution (30"), daily 2-m maximum air temperature; 2-m minimum air temperature; and precipitation simulations are developed for the mountainous Northeast by applying an additional level of downscaling to intermediately downscaled (1/ 8 deg) data using high-resolution topography and station observations. First, observed relationships between 2-m air temperature and elevation and between precipitation and elevation are derived. Then, these relationships are combined with spatial interpolation to enhance the resolution of intermediately downscaled GCM simulations. The resulting topographically downscaled dataset is analyzed for its ability to reproduce station observations. Topographic downscaling adds value to intermediately downscaled maximum and minimum 2-m air temperature at high-elevation stations, as well as moderately improves domain-averaged maximum and minimum 2-m air temperature. Topographic downscaling also improves mean precipitation but not daily probability distributions of precipitation. Overall, the utility of topographic downscaling is dependent on the initial bias of the intermediately downscaled product and the magnitude of the elevation adjustment. As the initial bias or elevation adjustment increases, more value is added to the topographically downscaled product.
This presentation documents Kennedy Space Center's Independent Assessment work completed on three assessments for the Ground Systems Development and Operations (GSDO) Program to assist the Chief Safety and Mission Assurance Officer during key programmatic reviews and provided the GSDO Program with analyses of how egress time affects the likelihood of astronaut and ground worker survival during an emergency. For each assessment, a team developed probability distributions for hazard scenarios to address statistical uncertainty, resulting in survivability plots over time. The first assessment developed a mathematical model of probabilistic survivability versus time to reach a safe location using an ideal Emergency Egress System at Launch Complex 39B (LC-39B); the second used the first model to evaluate and compare various egress systems under consideration at LC-39B. The third used a modified LC-39B model to determine if a specific hazard decreased survivability more rapidly than other events during flight hardware processing in Kennedy's Vehicle Assembly Building.
The lidar on the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission, makes robust measurements of dust and has generated a length of record that is significant both seasonally and inter-annually. We exploit this record to determine a multi-year climatology of the properties of Asian and Saharan dust, in particular seasonal optical depths, layer frequencies, and layer heights of dust gridded in accordance with the Level 3 data products protocol, between 2006-2015. The data are screened using standard CALIPSO quality assurance flags, cloud aerosol discrimination (CAD) scores, overlying features and layer properties. To evaluate the effects of transport on the morphology, vertical extent and size of the dust layers, we compare probability distribution functions of the layer integrated volume depolarization ratios, geometric depths and integrated attenuated color ratios near the source to the same distributions in the far field or transport region. CALIPSO is collaboration between NASA and Centre National D'études Spatiales (CNES), was launched in April 2006 to provide vertically resolved measurements of cloud and aerosol distributions. The primary instrument on the CALIPSO satellite is the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP), a near-nadir viewing two-wavelength polarization-sensitive instrument. The unique nature of CALIOP measurements make it quite challenging to validate backscatter profiles, aerosol type, and cloud phase, all of which are used to retrieve extinction and optical depth. To evaluate the uncertainty in the lidar ratios, we compare the values computed from dust layers overlying opaque water clouds, considered nominal, with the constant lidar ratio value used in the CALIOP algorithms for dust. We also explore the effects of noise on the CALIOP retrievals at daytime by comparing the distributions of the properties at daytime to the nighttime distributions.
Multi-person teams are sometimes responsible for critical tasks, such as flying an airliner. Here we present a method using gaze tracking data to assess shared visual attention, a term we use to describe the situation where team members are attending to a common set of elements in the environment. Gaze data are quantized with respect to a set of N areas of interest (AOIs); these are then used to construct a time series of N dimensional vectors, with each vector component representing one of the AOIs, all set to 0 except for the component corresponding to the currently fixated AOI, which is set to 1. The resulting sequence of vectors can be averaged in time, with the result that each vector component represents the proportion of time that the corresponding AOI was fixated within the given time interval. We present two methods for comparing sequences of this sort, one based on computing the time-varying correlation of the averaged vectors, and another based on a chi-square test testing the hypothesis that the observed gaze proportions are drawn from identical probability distributions. We have evaluated the method using synthetic data sets, in which the behavior was modeled as a series of "activities," each of which was modeled as a first-order Markov process. By tabulating distributions for pairs of identical and disparate activities, we are able to perform a receiver operating characteristic (ROC) analysis, allowing us to choose appropriate criteria and estimate error rates. We have applied the methods to data from airline crews, collected in a high-fidelity flight simulator (Haslbeck, Gontar & Schubert, 2014). We conclude by considering the problem of automatic (blind) discovery of activities, using methods developed for text analysis.
With a proliferation of new and unconventional vehicles and operations expected in the future, the ab initio airspace design will require new approaches to trajectory prediction for separation assurance and other air traffic management functions. This paper presents an approach to probabilistic modeling of the trajectory of an aircraft when its intent is unknown. The approach uses a set of feature functions to constrain a maximum entropy probability distribution based on a set of observed aircraft trajectories. This model can be used to sample new aircraft trajectories to form an ensemble reflecting the variability in an aircraft's intent. The model learning process ensures that the variability in this ensemble reflects the behavior observed in the original data set. Computational examples are presented.
Deep convective clouds (DCCs) are ideal visible calibration targets because they are bright nearly isotropic solar reflectors located over the tropics and they can be easily identified using a simple infrared threshold. Because all satellites view DCCs, DCCs provide the opportunity to uniformly monitor the stability of all operational sensors, both historical and present. A collective DCC anisotropically corrected radiance calibration approach is used to construct monthly probability distribution functions (PDFs) to monitor sensor stability. The DCC calibration targets were stable to within 0.5% and 0.3% per decade when the selection criteria were optimized based on Aqua MODerate Resolution Imaging Spectroradiometer 0.65-micrometer-band radiances. The Tropical Western Pacific (TWP), African, and South American regions were identified as the dominant DCC domains. For the 0.65-micrometer band, the PDF mode statistic is preferable, providing 0.3%regional consistency and 1%temporal uncertainty over land regions. It was found that the DCC within the TWP had the lowest radiometric response and DCC over land did not necessarily have the highest radiometric response. For wavelengths greater than 1 micrometer, the mean statistic is preferred, and land regions provided a regional variability of 0.7%with a temporal uncertainty of 1.1% where the DCC land response was higher than the response over ocean. Unlike stratus and cirrus clouds, the DCC spectra were not affected by water vapor absorption.
Changes in urine chemistry, during and post flight, potentially increases the risk of renal stones in astronauts. Although much is known about the effects of space flight on urine chemistry, no inflight incidence of renal stones in US astronauts exists and the question "How much does this risk change with space flight?" remains difficult to accurately quantify. In this discussion, we tackle this question utilizing a combination of deterministic and probabilistic modeling that implements the physics behind free stone growth and agglomeration, speciation of urine chemistry and published observations of population renal stone incidences to estimate changes in the rate of renal stone presentation. The modeling process utilizes a Population Balance Equation based model developed in the companion IWS abstract by Kassemi et al. (2016) to evaluate the maximum growth and agglomeration potential from a specified set of urine chemistry values. Changes in renal stone occurrence rates are obtained from this model in a probabilistic simulation that interrogates the range of possible urine chemistries using Monte Carlo techniques. Subsequently, each randomly sampled urine chemistry undergoes speciation analysis using the well-established Joint Expert Speciation System (JESS) code to calculate critical values, such as ionic strength and relative supersaturation. The Kassemi model utilizes this information to predict the mean and maximum stone size. We close the assessment loop by using a transfer function that estimates the rate of stone formation from combining the relative supersaturation and both the mean and maximum free stone growth sizes. The transfer function is established by a simulation analysis which combines population stone formation rates and Poisson regression. Training this transfer function requires using the output of the aforementioned assessment steps with inputs from known non-stone-former and known stone-former urine chemistries. Established in a Monte Carlo system, the entire renal stone analysis model produces a probability distribution of the stone formation rate and an expected uncertainty in the estimate. The utility of this analysis will be demonstrated by showing the change in renal stone occurrence predicted by this method using urine chemistry distributions published in Whitson et al. 2009. A comparison to the model predictions to previous assessments of renal stone risk will be used to illustrate initial validation of the model.
This work presents a computationally-efficient inverse approach to probabilistic damage diagnosis. Given strain data at a limited number of measurement locations, Bayesian inference and Markov Chain Monte Carlo (MCMC) sampling are used to estimate probability distributions of the unknown location, size, and orientation of damage. Substantial computational speedup is obtained by replacing a three-dimensional finite element (FE) model with an efficient surrogate model. The approach is experimentally validated on cracked test specimens where full field strains are determined using digital image correlation (DIC). Access to full field DIC data allows for testing of different hypothetical sensor arrangements, facilitating the study of strain-based diagnosis effectiveness as the distance between damage and measurement locations increases. The ability of the framework to effectively perform both probabilistic damage localization and characterization in cracked plates is demonstrated and the impact of measurement location on uncertainty in the predictions is shown. Furthermore, the analysis time to produce these predictions is orders of magnitude less than a baseline Bayesian approach with the FE method by utilizing surrogate modeling and effective numerical sampling approaches.
A human-in-the-loop experiment was conducted with 15 retired air traffic controllers to investigate two research questions: (a) what procedures are appropriate for the use of unmanned aircraft system (UAS) detect-and-avoid systems, and (b) how long in advance of a predicted close encounter should pilots request or execute a separation maneuver. The controller participants managed a busy Oakland air route traffic control sector with mixed commercial/general aviation and manned/UAS traffic, providing separation services, miles-in-trail restrictions and issuing traffic advisories. Controllers filled out post-scenario and post-simulation questionnaires, and metrics were collected on the acceptability of procedural options and temporal thresholds. The states of aircraft were also recorded when controllers issued traffic advisories. Subjective feedback indicated a strong preference for pilots to request maneuvers to remain well clear from intruder aircraft rather than deviate from their IFR clearance. Controllers also reported that maneuvering at 120 seconds until closest point of approach (CPA) was too early; maneuvers executed with less than 90 seconds until CPA were more acceptable. The magnitudes of the requested maneuvers were frequently judged to be too large, indicating a possible discrepancy between the quantitative UAS well clear standard and the one employed subjectively by manned pilots. The ranges between pairs of aircraft and the times to CPA at which traffic advisories were issued were used to construct empirical probability distributions of those metrics. Given these distributions, we propose that UAS pilots wait until an intruder aircraft is approximately 80 seconds to CPA or 6 nmi away before requesting a maneuver, and maneuver immediately if the intruder is within 60 seconds and 4 nmi. These thresholds should make the use of UAS detect and avoid systems compatible with current airspace procedures and controller expectations.
The performance and facility effect characterization tests of NASA's 12.5-kW Hall Effect Rocket with Magnetic Shielding had been completed. As a part of these tests, three plasma oscillation characterization studies were performed to help determine operation settings and quantify margins. The studies included the magnetic field strength variation study, background pressure effect study, and cathode flow fraction study. Separate high-speed videos of the thruster including the cathode and of only the cathode were recorded. Breathing mode at 10-15 kHz and cathode gradient-driven mode at 60-75 kHz were observed. An additional high frequency (40-70 kHz) global oscillation mode with sinusoidal probability distribution function was identified.
Physical characteristics of NEAs are an essential input to modeling behavior during atmospheric entry and to assess the risk of impact but determining these properties requires a non-trivial investment of time and resources. The characteristics relevant to these models include size, density, strength and ablation coefficient. Some of these characteristics cannot be directly measured, but rather must be inferred from related measurements of asteroids and/or meteorites. Furthermore, for the majority of NEAs, only the basic measurements exist so often properties must be inferred from statistics of the population of more completely characterized objects. The Asteroid Threat Assessment Project at NASA Ames Research Center has developed a probabilistic asteroid impact risk (PAIR) model in order to assess the risk of asteroid impact. Our PAIR model and its use to develop probability distributions of impact risk are discussed in other contributions to PDC 2017 (e.g., Mathias et al.). Here we utilize PAIR to investigate which NEA characteristics are important for assessing the impact threat by investigating how changes in these characteristics alter the damage predicted by PAIR. We will also provide an assessment of the current state of knowledge of the NEA characteristics of importance for asteroid threat assessment. The relative importance of different properties as identified using PAIR will be combined with our assessment of the current state of knowledge to identify potential high impact investigations. In addition, we will discuss an ongoing effort to collate the existing measurements of NEA properties of interest to the planetary defense community into a readily accessible database.
It is hypothesized that the interactions between many of the most massive derelicts in low Earth orbit are more frequent than modeled by the traditional combination of kinetic theory of gases and Poisson probability distribution function. This is suggested by the fact that there are clusters of derelicts where members' inclinations are nearly identical and their apogees/perigees overlap significantly resulting in periodic synchronization of the objects' orbits. In order to address this proposition, an experiment was designed and conducted over the last two years. Results from this monitoring and characterization experiment are presented with implications for proposed debris remediation strategies. Four separate clusters of massive derelicts were examined that are centered around 775km, 850km, 975km, and 1500km, respectively. In aggregate, the constituents of these clusters contain around 500 objects and about 800,000kg of mass; this equates to a third of all derelict mass in LEO. Preliminary analysis indicates that encounter rates over this time period for these objects are greater than is estimated by traditional techniques. Hypothesized dependencies between latitude of encounter, relative velocity, frequency of encounters, inclination, and differential semi-major axis were established and verified. This experiment also identified specific repeatable cluster dynamics that may reduce the cost/risk and enhance the effectiveness of debris remediation activities and also enable new operational debris remediation options.
The lack of well-developed techniques for modeling changing statistical moments in our observations has stymied the application of stochastic process theory for many scientific and engineering applications. Non linear effects of the observation methodology is one of the most perplexing aspects to modeling non stationary processes. This perplexing problem was encountered when modeling the effect of non stationary receiver fluctuations on the performance of radiometer calibration architectures. Existing modeling approaches were found not applicable; particularly problematic is modeling processes across scales over which they begin to exhibit non stationary behavior within the time interval of the calibration algorithm. Alternatively, the radiometer output is modeled as samples from a sequence random variables; the random variables are treated using a conditional probability distribution function conditioned on the use of the variable in the calibration algorithm. This approach of treating a process as a sequence of random variables with non stationary stochastic moments produce sensible predictions of temporal effects of calibration algorithms. To test these model predictions, an experiment using the Millimeter wave Imaging Radiometer (MIR) was conducted. The MIR with its two black body calibration references was configured in a laboratory setting to observe a third ultra-stable reference (CryoTarget). The MIR was programmed to sequentially sample each of the three references in approximately a 1 second cycle. Data were collected over a six-hour interval. The sequence of reference measurements form an ensemble sample set comprised of a series of three reference measurements. Two references are required to estimate the receiver response. A third reference is used to estimate the uncertainty in the estimate. Typically, calibration algorithms are designed to suppress the non stationary effects of receiver fluctuations. By treating the data sequence as an ensemble collection, it is possible to apply temporal algorithms which exacerbate the non stationary effects. By varying the algorithm, information about the properties of the non stationary receiver fluctuations is obtained. Comparisons of analytical calculations and statistical analysis of data demonstrate impressive agreement.
MODIS reflective solar bands are calibrated on-orbit using a solar diffuser and near-monthly lunar observations. To monitor the performance and effectiveness of the on-orbit calibrations, pseudo-invariant targets such as deep convective clouds (DCCs), Libya-4, and Dome-C are used to track the long-term stability of MODIS Level 1B product. However, the current MODIS operational DCC technique (DCCT) simply uses the criteria set for the 0.65- m band. We optimize several critical DCCT parameters including the 11- micrometer IR-band Brightness Temperature (BT11) threshold for DCC identification, DCC core size and uniformity to help locate DCCs at convection centers, data collection time interval, and probability distribution function (PDF) bin increment for each channel. The mode reflectances corresponding to the PDF peaks are utilized as the DCC reflectances. Results show that the BT11 threshold and time interval are most critical for the Short Wave Infrared (SWIR) bands. The Bidirectional Reflectance Distribution Function model is most effective in reducing the DCC anisotropy for the visible channels. The uniformity filters and PDF bin size have minimal impacts on the visible channels and a larger impact on the SWIR bands. The newly optimized DCCT will be used for future evaluation of MODIS on-orbit calibration by MODIS Characterization Support Team.
Stochastic Reduced Order Models with Python (SROMPy) is a software package developed to enable user-friendly utilization of the stochastic reduced order model (SROM) approach for uncertainty quantification. A SROM is a low dimensional, discrete approximation to a random quantity that enables efficient and non-intrusive stochastic computations. With SROMPy, a user can easily generate a SROM to approximate a random variable or vector described by several different types of probability distributions using the Python programming language. Once a SROM is constructed, the software can be used to propagate uncertainty through a user-defined computational model to estimate statistics of a given quantity of interest. This report is meant to introduce the SROMPy module and brie y demonstrate its capabilities. A simple example of a spring-mass system with a random input is included to illustrate the practicality of the SROM approach to uncertainty quantification and relative ease of applying it with SROMPy. The example includes a comparison with a solution obtained using classical Monte Carlo simulation, demonstrating the similarities and advantages of using the SROM approach.
Advances in computing capabilities and scientific development have come together to evolve general circulation models into multi-scale Earth system modeling tools. The Goddard Earth Observing System (GEOS) model is one such example of this evolution. The GEOS model is driven by the finite-volume cubed-sphere (FV3) non-hydrostatic dynamical core. Surrounding FV3 is a scale-aware physics package and data assimilation capability permitting multi-scale application of GEOS for sub-seasonal to seasonal climate prediction, medium range weather prediction, and global mesoscale modeling at convection allowing resolutions. GEOS also includes a comprehensive chemistry package representing a range of capabilities from basic chemistry and interactve aerosols and gaseous species, to carbon emissions and uptake, and complex ozone photochemistry. In this study, we apply the GEOS model to study the fidelity of these processes with increasing horizontal resolution and scale-aware processes in GEOS on extreme events. The GEOS model is run for 40-days beginning in August 2016 at three uniform global resolutions of 13-km (c768), 6-km (c1536) and 3-km (c3072) with 72 vertical levels up to 0.01mb. The model physics use the Grell-Freitas scale-aware convection scheme to dynamically reduce the role of parameterized deep convection as resolved scale processes in the model take over at higher resolutions. We will include high-resolution global emissions and fluxes of aerosols and carbon downscaled from recent satellite observations. We will compare these simulations with reanalyses and observations, focusing on rainfall, clouds and radiative forcing at hourly to monthly timescales. We will closely examine the probability distribution of precipitation intensities and radiative properties of clouds on a daily time scale. In addition, we will focus on extreme events, in particular the diurnal cycle of convection over the US and the frequency and physical nature of organized convection and heavy rain events across the globe.