Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Probability distribution”

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 631 records · Page 35

Comparing Evaporative Sources of Terrestrial Precipitation and Their Extremes in MERRA Using Relative Entropy

A quasi-isentropic back trajectory scheme is applied to output from the Modern Era Retrospective-analysis for Research and Applications and a land-only replay with corrected precipitation to estimate surface evaporative sources of moisture supplying precipitation over every ice-free land location for the period 1979-2005. The evaporative source patterns for any location and time period are effectively two dimensional probability distributions. As such, the evaporative sources for extreme situations like droughts or wet intervals can be compared to the corresponding climatological distributions using the method of relative entropy. Significant differences are found to be common and widespread for droughts, but not wet periods, when monthly data are examined. At pentad temporal resolution, which is more able to isolate floods and situations of atmospheric rivers, values of relative entropy over North America are typically 50-400 larger than at monthly time scales. Significant differences suggest that moisture transport may be the key to precipitation extremes. Where evaporative sources do not change significantly, it implies other local causes may underlie the extreme events.

Climatological Distributions↗

Analytical Assessment of Simultaneous Parallel Approach Feasibility from Total System Error

In a simultaneous paired approach to closely-spaced parallel runways, a pair of aircraft flies in close proximity on parallel approach paths. The aircraft pair must maintain a longitudinal separation within a range that avoids wake encounters and, if one of the aircraft blunders, avoids collision. Wake avoidance defines the rear gate of the longitudinal separation. The lead aircraft generates a wake vortex that, with the aid of crosswinds, can travel laterally onto the path of the trail aircraft. As runway separation decreases, the wake has less distance to traverse to reach the path of the trail aircraft. The total system error of each aircraft further reduces this distance. The total system error is often modeled as a probability distribution function. Therefore, Monte-Carlo simulations are a favored tool for assessing a "safe" rear-gate. However, safety for paired approaches typically requires that a catastrophic wake encounter be a rare one-in-a-billion event during normal operation. Using a Monte-Carlo simulation to assert this event rarity with confidence requires a massive number of runs. Such large runs do not lend themselves to rapid turn-around during the early stages of investigation when the goal is to eliminate the infeasible regions of the solution space and to perform trades among the independent variables in the operational concept. One can employ statistical analysis using simplified models more efficiently to narrow the solution space and identify promising trades for more in-depth investigation using Monte-Carlo simulations. These simple, analytical models not only have to address the uncertainty of the total system error but also the uncertainty in navigation sources used to alert an abort of the procedure. This paper presents a method for integrating total system error, procedure abort rates, avionics failures, and surveillance errors into a statistical analysis that identifies the likely feasible runway separations for simultaneous paired approaches.

Madden, Michael M.↗

On a Solar Origin for the Cosmogenic Nuclide Event of 775 A.D.

We explore requirements for a solar particle event (SPE) and flare capable of producing the cosmogenic nuclide event of 775 A.D., and review solar circumstances at that time. A solar source for 775 would require a greater than 1 GV spectrum approximately 45 times stronger than that of the intense high-energy SPE of 1956 February 23. This implies a greater than 30 MeV proton fluence (F(sub 30)) of approximately 8 × 10(exp 10) proton cm(exp −2), approximately 10 times larger than that of the strongest 3 month interval of SPE activity in the modern era. This inferred F(sub 30) value for the 775 SPE is inconsistent with the occurrence probability distribution for greater than 30 MeV solar proton events. The best guess value for the soft X-ray classification (total energy) of an associated flare is approximately X230 (approximately 9 × 10(exp 33) erg). For comparison, the flares on 2003 November 4 and 1859 September 1 had observed/inferred values of approximately X35 (approximately 10(exp 33) erg) and approximately X45 (approximately 2 × 10(exp 33) erg), respectively. The estimated size of the source active region for a approximately 10(exp 34) erg flare is approximately 2.5 times that of the largest region yet recorded. The 775 event occurred during a period of relatively low solar activity, with a peak smoothed amplitude about half that of the second half of the 20th century. The approximately 1945-1995 interval, the most active of the last approximately 2000 yr, failed to witness a SPE comparable to that required for the proposed solar event in 775. These considerations challenge a recent suggestion that the 775 event is likely of solar origin.

solar particles↗

NASA Instrument Cost/Schedule Model

NASA's Office of Independent Program and Cost Evaluation (IPCE) has established a number of initiatives to improve its cost and schedule estimating capabilities. 12One of these initiatives has resulted in the JPL developed NASA Instrument Cost Model. NICM is a cost and schedule estimator that contains: A system level cost estimation tool; a subsystem level cost estimation tool; a database of cost and technical parameters of over 140 previously flown remote sensing and in-situ instruments; a schedule estimator; a set of rules to estimate cost and schedule by life cycle phases (B/C/D); and a novel tool for developing joint probability distributions for cost and schedule risk (Joint Confidence Level (JCL)). This paper describes the development and use of NICM, including the data normalization processes, data mining methods (cluster analysis, principal components analysis, regression analysis and bootstrap cross validation), the estimating equations themselves and a demonstration of the NICM tool suite.

JCL↗

The Gemini NICI Planet-Finding Campaign: The Offset Ring of HR 4796 A

We present J, H, CH4 short (1.578 micrometers), CH4 long (1.652 micrometers) and K(sub s)-band images of the dust ring around the 10 Myr old star HR 4796 A obtained using the Near Infrared Coronagraphic Imager (NICI) on the Gemini-South 8.1 m Telescope. Our images clearly show for the first time the position of the star relative to its circumstellar ring thanks to NICI's translucent focal plane occulting mask. We employ a Bayesian Markov chain Monte Carlo method to constrain the offset vector between the two. The resulting probability distribution shows that the ring center is offset from the star by 16.7 +/- 1.3 milliarcseconds along a position angle of 26 +/- 3deg, along the PA of the ring, 26.47 +/- 0.04deg. We find that the size of this offset is not large enough to explain the brightness asymmetry of the ring. The ring is measured to have mostly red reflectivity across the JHKs filters, which seems to indicate micron-sized grains. Just like Neptune's 3:2 and 2:1 mean-motion resonances delineate the inner and outer edges of the classical Kuiper belt, we find that the radial extent of the HR 4796 A and the Fomalhaut rings could correspond to the 3:2 and 2:1 mean-motion resonances of hypothetical planets at 54.7 AU and 97.7 AU in the two systems, respectively. A planet orbiting HR 4796 A at 54.7 AU would have to be less massive than 1.6 Jup mass so as not to widen the ring too much by stirring.

Near Infrared Coronagraphic Imager (NICI)↗

High-Resolution NU-WRF Simulations of a Deep Convective-Precipitation System During MC3E: Comparisons Between Goddard Microphysics Schemes and Observations - Part 1

The Goddard microphysics scheme was recently improved by adding a 4th ice class (frozen dropshail). This new 4ICE scheme was implemented and tested in the Goddard Cumulus Ensemble model (GCE) for an intense continental squall line and a moderate,less-organized continental case. Simulated peak radar reflectivity profiles were improved both in intensity and shape for both cases as were the overall reflectivity probability distributions versus observations. In this study, the new Goddard 4ICE scheme is implemented into the regional-scale NASA Unified - Weather Research and Forecasting model (NU-WRF) and tested on an intense mesoscale convective system that occurred during the Midlatitude Continental Convective Clouds Experiment (MC3E). The NU42WRF simulated radar reflectivities, rainfall intensities, and vertical and horizontal structure using the new 4ICE scheme agree as well as or significantly better with observations than when using previous versions of the Goddard 3ICE (graupel or hail) schemes. In the 4ICE scheme, the bin microphysics-based rain evaporation correction produces more erect convective cores, while modification of the unrealistic collection of ice by dry hail produces narrow and intense cores, allowing more slow-falling snow to be transported rearward. Together with a revised snow size mapping, the 4ICE scheme produces a more horizontally stratified trailing stratiform region with a broad, more coherent light rain area. In addition, the NU-WRF 4ICE simulated radar reflectivity distributions are consistent with and generally superior to those using the GCE due to the less restrictive open lateral boundaries

Regional Scale Model↗

Probabilistic Survivability Versus Time Modeling

This technical paper documents Kennedy Space Centers Independent Assessment team work completed on three assessments for the Ground Systems Development and Operations (GSDO) Program to assist the Chief Safety and Mission Assurance Officer (CSO) and GSDO management during key programmatic reviews. The assessments provided the GSDO Program with an analysis of how egress time affects the likelihood of astronaut and worker survival during an emergency. For each assessment, the 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 Kennedys Vehicle Assembly Building (VAB).Based on the composite survivability versus time graphs from the first two assessments, there was a soft knee in the Figure of Merit graphs at eight minutes (ten minutes after egress ordered). Thus, the graphs illustrated to the decision makers that the final emergency egress design selected should have the capability of transporting the flight crew from the top of LC 39B to a safe location in eight minutes or less. Results for the third assessment were dominated by hazards that were classified as instantaneous in nature (e.g. stacking mishaps) and therefore had no effect on survivability vs time to egress the VAB. VAB emergency scenarios that degraded over time (e.g. fire) produced survivability vs time graphs that were line with aerospace industry norms.

Joyner, James J., Sr.↗

L'CO/LFIR Relations with CO Rotational Ladders of Galaxies Across the Herschel SPIRE Archive

We present a catalog of all CO (carbon monoxide) (J (total angular momentum quantum number) equals 1-0 through J equals 13-12), [CI], and [NII] lines available from extragalactic spectra from the Herschel SPIRE (Spectral and Photometric Imaging Receiver) Fourier Transform Spectrometer (FTS) archive combined with observations of the low-J CO lines from the literature and from the Arizona Radio Observatory. This work examines the relationships between L (sub FIR (Far Infra Red)), L prime (sub CO), and L (sub CO) / L (sub CO, 1-0). We also present a new method for estimating probability distribution functions (PDFs) from marginal signal-to-noise ratio Herschel FTS spectra, which takes into account the instrumental "ringing" and the resulting highly correlated nature of the spectra. The slopes of Log (L (sub FIR)) versus Log (L prime (sub CO)) are linear for all mid- to high-J CO lines. The mid- to high-J CO luminosity relative to CO J equals 1-0 increases with increasing L (sub FIR), indicating higher excitement of the molecular gas, though these ratios do not exceed approximately 100. The luminosities relative to CO J equals 1-0 remain relatively at from J equals 6-5 through J equals13-2, across many orders of magnitude of L (sub FIR). Qualitative comparisons to current theoretical models do not match these flat SLED (spectral line energy distributions) shapes, indicating the need for more comprehensive modeling of the excitation processes of warm molecular gas in nearby galaxies.

Submillimeter Data↗

Detecting Climate Signals in Precipitation Extremes from TRMM (1998-2013) - Increasing Contrast Between Wet and Dry Extremes During the "Global Warming Hiatus"

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.

Precipitation↗

Development and Evaluation of High-Resolution Climate Simulations Over the Mountainous Northeastern United States

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.

climatology↗

Probabilistic Survivability Versus Time Modeling

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.

Operational Hazards↗

Nearly a Decade of CALIPSO Observations of Asian and Saharan Dust Properties Near Source and Transport Regions

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.

Omar, Ali H.↗

Measuring and Modeling Shared Visual Attention

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.

attention↗

Probabilistic Modeling of Aircraft Trajectories for Dynamic Separation Volumes

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.

Lewis, Timothy A.↗

The Characterization of Deep Convective Clouds as an Invariant Calibration Target and as a Visible Calibration Technique

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.

Doelling, David R.↗

Estimating the Rate of Occurrence of Renal Stones in Astronauts

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.

gravitational physiology↗

A Computationally-Efficient Inverse Approach to Probabilistic Strain-Based Damage Diagnosis

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.

Warner, James E.↗