SEARCH · Engineering Papers

Results for “PROBABILITY”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,639 records · Page 91

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Optimum space shuttle launch times relative to natural environment

Three sets of meteorological criteria were analyzed to determine the probabilities of favorable launch and landing conditions. Probabilities were computed for every 3 hours on a yearly basis using 14 years of weather data. These temporal probability distributions, applicable to the three sets of weather criteria encompassing benign, moderate and severe weather conditions, were computed for both Kennedy Space Center (KSC) and Edwards Air Force Base. In addition, conditional probabilities were computed for unfavorable weather conditions occurring after a delay which may or may not be due to weather conditions. Also, for KSC, the probabilities of favorable landing conditions at various times after favorable launch conditions have prevailed have been computed so that mission probabilities may be more accurately computed for those time periods when persistence strongly correlates weather conditions. Moreover, the probabilities and conditional probabilities of the occurrence of both favorable and unfavorable events for each individual criterion were computed to indicate the significance of each weather element to the overall result.

King, R. L.↗

NASA NTRS1960

Measurement of the Errors of Service Altimeter Installations During Landing-Approach and Take-Off Operations

The overall errors of the service altimeter installations of a variety of civil transport, military, and general-aviation airplanes have been experimentally determined during normal landing-approach and take-off operations. The average height above the runway at which the data were obtained was about 280 feet for the landings and about 440 feet for the take-offs. An analysis of the data obtained from 196 airplanes during 415 landing approaches and from 70 airplanes during 152 take-offs showed that: 1. The overall error of the altimeter installations in the landing- approach condition had a probable value (50 percent probability) of +/- 36 feet and a maximum probable value (99.7 percent probability) of +/- 159 feet with a bias of +10 feet. 2. The overall error in the take-off condition had a probable value of +/- 47 feet and a maximum probable value of +/- 207 feet with a bias of -33 feet. 3. The overall errors of the military airplanes were generally larger than those of the civil transports in both the landing-approach and take-off conditions. In the landing-approach condition the probable error and the maximum probable error of the military airplanes were +/- 43 and +/- 189 feet, respectively, with a bias of +15 feet, whereas those for the civil transports were +/- 22 and +/- 96 feet, respectively, with a bias of +1 foot. 4. The bias values of the error distributions (+10 feet for the landings and -33 feet for the take-offs) appear to represent a measure of the hysteresis characteristics (after effect and recovery) and friction of the instrument and the pressure lag of the tubing-instrument system.

Gracey, William↗

NASA NTRS2010

Tool for Forecasting Cool-Season Peak Winds Across Kennedy Space Center and Cape Canaveral Air Force Station

The expected peak wind speed for the day is an important element in the daily morning forecast for ground and space launch operations at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). The 45th Weather Squadron (45 WS) must issue forecast advisories for KSC/CCAFS when they expect peak gusts for >= 25, >= 35, and >= 50 kt thresholds at any level from the surface to 300 ft. In Phase I of this task, the 45 WS tasked the Applied Meteorology Unit (AMU) to develop a cool-season (October - April) tool to help forecast the non-convective peak wind from the surface to 300 ft at KSC/CCAFS. During the warm season, these wind speeds are rarely exceeded except during convective winds or under the influence of tropical cyclones, for which other techniques are already in use. The tool used single and multiple linear regression equations to predict the peak wind from the morning sounding. The forecaster manually entered several observed sounding parameters into a Microsoft Excel graphical user interface (GUI), and then the tool displayed the forecast peak wind speed, average wind speed at the time of the peak wind, the timing of the peak wind and the probability the peak wind will meet or exceed 35, 50 and 60 kt. The 45 WS customers later dropped the requirement for >= 60 kt wind warnings. During Phase II of this task, the AMU expanded the period of record (POR) by six years to increase the number of observations used to create the forecast equations. A large number of possible predictors were evaluated from archived soundings, including inversion depth and strength, low-level wind shear, mixing height, temperature lapse rate and winds from the surface to 3000 ft. Each day in the POR was stratified in a number of ways, such as by low-level wind direction, synoptic weather pattern, precipitation and Bulk Richardson number. The most accurate Phase II equations were then selected for an independent verification. The Phase I and II forecast methods were compared using an independent verification data set. The two methods were compared to climatology, wind warnings and advisories issued by the 45 WS, and North American Mesoscale (NAM) model (MesoNAM) forecast winds. The performance of the Phase I and II methods were similar with respect to mean absolute error. Since the Phase I data were not stratified by precipitation, this method's peak wind forecasts had a large negative bias on days with precipitation and a small positive bias on days with no precipitation. Overall, the climatology methods performed the worst while the MesoNAM performed the best. Since the MesoNAM winds were the most accurate in the comparison, the final version of the tool was based on the MesoNAM winds. The probability the peak wind will meet or exceed the warning thresholds were based on the one standard deviation error bars from the linear regression. For example, the linear regression might forecast the most likely peak speed to be 35 kt and the error bars used to calculate that the probability of >= 25 kt = 76%, the probability of >= 35 kt = 50%, and the probability of >= 50 kt = 19%. The authors have not seen this application of linear regression error bars in any other meteorological applications. Although probability forecast tools should usually be developed with logistic regression, this technique could be easily generalized to any linear regression forecast tool to estimate the probability of exceeding any desired threshold . This could be useful for previously developed linear regression forecast tools or new forecast applications where statistical analysis software to perform logistic regression is not available. The tool was delivered in two formats - a Microsoft Excel GUI and a Tool Command Language/Tool Kit (Tcl/Tk) GUI in the Meteorological Interactive Data Display System (MIDDS). The Microsoft Excel GUI reads a MesoNAM text file containing hourly forecasts from 0 to 84 hours, from one model run (00 or 12 UTC). The GUI then displays e peak wind speed, average wind speed, and the probability the peak wind will meet or exceed the 25-, 35- and 50-kt thresholds. The user can display the Day-1 through Day-3 peak wind forecasts, and separate forecasts are made for precipitation and non-precipitation days. The MIDDS GUI uses data from the NAM and Global Forecast System (GFS), instead of the MesoNAM. It can display Day-1 and Day-2 forecasts using NAM data, and Day-1 through Day-5 forecasts using GFS data. The timing of the peak wind is not displayed, since the independent verification showed that none of the forecast methods performed significantly better than climatology. The forecaster should use the climatological timing of the peak wind (2248 UTC) as a first guess and then adjust it based on the movement of weather features.

Barrett, Joe H., III↗

NASA NTRS2011

Tracking Object Existence From an Autonomous Patrol Vehicle

An autonomous vehicle patrols a large region, during which an algorithm receives measurements of detected potential objects within its sensor range. The goal of the algorithm is to track all objects in the region over time. This problem differs from traditional multi-target tracking scenarios because the region of interest is much larger than the sensor range and relies on the movement of the sensor through this region for coverage. The goal is to know whether anything has changed between visits to the same location. In particular, two kinds of alert conditions must be detected: (1) a previously detected object has disappeared and (2) a new object has appeared in a location already checked. For the time an object is within sensor range, the object can be assumed to remain stationary, changing position only between visits. The problem is difficult because the upstream object detection processing is likely to make many errors, resulting in heavy clutter (false positives) and missed detections (false negatives), and because only noisy, bearings-only measurements are available. This work has three main goals: (1) Associate incoming measurements with known objects or mark them as new objects or false positives, as appropriate. For this, a multiple hypothesis tracker was adapted to this scenario. (2) Localize the objects using multiple bearings-only measurements to provide estimates of global position (e.g., latitude and longitude). A nonlinear Kalman filter extension provides these 2D position estimates using the 1D measurements. (3) Calculate the probability that a suspected object truly exists (in the estimated position), and determine whether alert conditions have been triggered (for new objects or disappeared objects). The concept of a probability of existence was created, and a new Bayesian method for updating this probability at each time step was developed. A probabilistic multiple hypothesis approach is chosen because of its superiority in handling the uncertainty arising from errors in sensors and upstream processes. However, traditional target tracking methods typically assume a stationary detection volume of interest, whereas in this case, one must make adjustments for being able to see only a small portion of the region of interest and understand when an alert situation has occurred. To track object existence inside and outside the vehicle's sensor range, a probability of existence was defined for each hypothesized object, and this value was updated at every time step in a Bayesian manner based on expected characteristics of the sensor and object and whether that object has been detected in the most recent time step. Then, this value feeds into a sequential probability ratio test (SPRT) to determine the status of the object (suspected, confirmed, or deleted). Alerts are sent upon selected status transitions. Additionally, in order to track objects that move in and out of sensor range and update the probability of existence appropriately a variable probability detection has been defined and the hypothesis probability equations have been re-derived to accommodate this change. Unsupervised object tracking is a pervasive issue in automated perception systems. This work could apply to any mobile platform (ground vehicle, sea vessel, air vehicle, or orbiter) that intermittently revisits regions of interest and needs to determine whether anything interesting has changed.

Wolf, Michael↗

DOE OSTI2022

Fiber Uncertainty Visualization for Bivariate Data With Parametric and Nonparametric Noise Models

Visualization and analysis of multivariate data and their uncertainty are top research challenges in data visualization. Constructing fiber surfaces is a popular technique for multivariate data visualization that generalizes the idea of level-set visualization for univariate data to multivariate data. Here, in this paper, we present a statistical framework to quantify positional probabilities of fibers extracted from uncertain bivariate fields. Specifically, we extend the state-of-the-art Gaussian models of uncertainty for bivariate data to other parametric distributions (e.g., uniform and Epanechnikov) and more general nonparametric probability distributions (e.g., histograms and kernel density estimation) and derive corresponding spatial probabilities of fibers. In our proposed framework, we leverage Green's theorem for closed-form computation of fiber probabilities when bivariate data are assumed to have independent parametric and nonparametric noise. Additionally, we present a nonparametric approach combined with numerical integration to study the positional probability of fibers when bivariate data are assumed to have correlated noise. For uncertainty analysis, we visualize the derived probability volumes for fibers via volume rendering and extracting level sets based on probability thresholds. We present the utility of our proposed techniques via experiments on synthetic and simulation datasets.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2014

Quantifying the 'naturalness' of the curvaton model

We investigate the probability of obtaining an observable curvature perturbation, using as an example the minimal curvaton-higgs (MCH) model. We determine ''probably observable'' and ''probably excluded'' regions of parameter space assuming generic initial conditions and applying a stochastic approach for the curvaton's evolution during inflation. Inflation is assumed to last longer than the N{sub obs} ≅ 55 observable e-folds, and the total number of e-folds of inflation determines the particular ranges of parameters that are probable. For the MCH model, these ''probably observable'' regions always lie within the range 8 × 10{sup 4} GeV ≤ m{sub σ} ≤ 2 × 10{sup 7} GeV, where m{sub σ} is the curvaton mass, and the Hubble scale at horizon exit is chosen as H{sub *} = 10{sup 10} GeV. Because the ''probably observable'' region depends on the total duration of inflation, information on parameters in the Lagrangian from particle physics and from precision CMB observations can therefore provide information about the total duration of inflation, not just the last N{sub obs} e-folds. This method could also be applied to any model that contains additional scalar fields to determine the probability that these scalar fields contribute to the curvature perturbation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2018

Trinal Nature of the Excitation of Bulk Plasmons by a Fast Charged Particle at Grazing Incidence on the Interface between Vacuum and a Metal with Strong Spatial Dispersion

The excitation and propagation of bulk and surface (surface waves of transition radiation of plasmons at frequencies above the plasma frequency) plasma waves by incident electrons moving both in vacuum toward the surface of a metal and inside the metal whose boundary elastically and spectacularly reflects internal nonequilibrium electrons have been analyzed. In contrast to work [B. N. Libenson, J. Exp. Theor. Phys. 113, 553 (2011)], attention in this work is focused on the influence of surface effects on bulk-plasmon excitation by incident electrons. The probabilities and spectra of single characteristic energy loss of an intermediate- energy electron (50–500 eV) moving at an angle to the surface of the medium in three regions: in vacuum, in the medium, and again in vacuum after the electron leaves the medium are calculated. The kinetic approximation is used for the dielectric function, where the entire range of the plasmon spectrum is taken into account correctly for the problem under consideration. In the indicated energy range of incident electrons, surface effects, on the one hand, significantly reduces the probability of excitation of bulk plasma waves in the medium with strong spatial dispersion, in particular, as compared to the results obtained in [B.N. Libenson, J. Exp. Theor. Phys. 113, 553 (2011)], where surface effects were disregarded and the probability of bulkplasmon excitation by a 200-eV electron incident and emitted perpendicularly to the boundary is about one third of that in the unbounded medium. On the other hand, at grazing incidence from vacuum, the probability of transition radiation of bulk plasmons increases significantly and can lead to a change in the character of the angular dependence of the intensity of bulk plasma energy loss. Thus, the main result of this work is that a decrease in the glancing angle of the fast electron with respect to the vacuum–metal interface is accompanied both by an increase in the contribution from the transition mechanism to the probability of bulk-plasmon excitation in the vacuum region and by a decrease in the contribution from Cherenkov and bremsstrahlung mechanisms of excitation in the medium. The probability of bulk-plasmon excitation in the vacuum region exceeds the probability of excitation at the further motion of the electron in metallic aluminum at angles of incidence larger than 65°, 70°, and 75° at the energies E = 200, 350, and 500 eV, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Reply to “Comment on ‘The Reduced Detection Rate of Signals That Are Hidden by Earthquakes: Case Studies with Spotlight Detectors That Operate at Seismic Arrays,’ by Joshua D. Carmichael, Brent G. Delbridge, and Richard Alfaro-Diaz” by Paul G. Richards

We reply to a comment that Paul G. Richards (Richards, 2026) directed to an article by Carmichael et al. (2025). Richards expressed concern over two issues. The first issue relates to the size of the explosions discussed in the article. The second relates to the probability that earthquakes and explosions sourced near the same location coincide in time, by chance. We do not address the first issue because the original article does not assign significance to the size of the explosions from the point of view of an experimental party. We do respond to Richards’ comment on the second issue and claim that it can be explained through a clarification about conditional versus joint probabilities. Our present discussion, therefore, describes the difference between the conditional probability that a correlator fails, given that an explosion and earthquake coincidently occur in the same place, and the joint probability that a correlator fails while an earthquake and explosion also coincidentally occur in the same place. This latter, joint probability appears to be rare. We clarify that our original article did not treat the joint probability, only the conditional probability.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2013

Probabilistic Hazard Assessment for Tornadoes, Straight-line Wind, and Extreme Precipitation at the Savannah River Site

Recent data sets for three meteorological phenomena with the potential to inflict damage on SRS facilities - tornadoes, straight winds, and heavy precipitation - are analyzed using appropriate statistical techniques to estimate occurrence probabilities for these events in the future. Summaries of the results for DOE-mandated return periods and comparisons to similar calculations performed in 1998 by Weber, et al., are given. Using tornado statistics for the states of Georgia and South Carolina, we calculated the probability per year of any location within a 2⁰ square area surrounding SRS being struck by a tornado (the ‘strike’ probability) and the probability that any point will experience winds above set thresholds. The strike probability was calculated to be 1.15E-3 (1 chance in 870) per year and wind speeds for DOE mandated return periods of 50,000 years, 125,000 years, and 1E+7 years (USDOE, 2012) were estimated to be 136 mph, 151 mph and 221 mph, respectively. In 1998 the strike probability for SRS was estimated to be 3.53 E-4 and the return period wind speeds were 148 mph every 50,000 years and 180 mph every 125,000 years. A 1E+7 year tornado wind speed was not calculated in 1998; however a 3E+6 year wind speed was 260 mph. The lower wind speeds resulting from this most recent analysis are largely due to new data since 1998, and to a lesser degree differences in the models used. By contrast, default tornado wind speeds taken from ANSI/ANS-2.3-2011 are somewhat higher: 161 mph for return periods of 50,000 years, 173 mph every 125,000 years, and 230 mph every 1E+7 years (ANS, 2011). Although the ANS model and the SRS models are very similar, the region defined in ANS 2.3 that encompasses the SRS also includes areas of the Great Plains and lower Midwest, regions with much higher occurrence frequencies of strong tornadoes. The SRS straight wind values associated with various return periods were calculated by fitting existing wind data to a Gumbel distribution, and extrapolating the values for any return period from the tail of that function. For the DOE mandated return periods, we expect straight winds of 123 mph every 2500 years, and 132mph every 6250 years at any point within the SRS. These values are similar to those from the W98 report (which also used the Gumbel distribution for wind speeds) which gave wind speeds of 115mph and 122 mph for return periods of 2500 years and 6250 years, respectively. For extreme precipitation accumulation periods, we compared the fits of three different theoretical extreme-value distributions, and in the end decided to maintain the use of the Gumbel distribution for each period. The DOE mandated 6-hr accumulated rainfall for return periods of 2500 years and 6250 years was estimated as 7.8 inches and 8.4 inches, respectively. For the 24-hr rainfall return periods of 10,000 years and 25,000 years, total rainfall estimates were 10.4 inches and 11.1 inches, respectively. These values are substantially lower than comparable values provided in the W98 report. This is largely a consequence of the W98 use of a different extreme value distribution with its corresponding higher extreme probabilities.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Compare $\overline{ν}$ parameterizations for 235 U(n,f) with fitted fission barriers

In this document, we compare observables computed with CGMF to some available data for 235 U(n,f), to investigate the multi-chance fission probabilities (fission barriers from CoH were refit) and two TKE parametrization. The four parametrizations include i) default CGMF multi-chance fission probabilities with bends in TKE for all fission chances, ii) default CGMF multi-chance fission probabilities with a TKE bend in only first-chance fission iii) refitted multi-chance fission probabilities with bends in TKE for all fission chances, and iv) refitted multi-chance fission probabilities with a bend in TKE for only first-chance fission. In all figures where the four parametrizations are on the same plot, the ones with bends in TKE for all chance fission are in black, with only first-chance bend in red, the default multi-chance fission probabilities are dashed and the new multi-chance fission probabilities are the solid lines.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗

DOE OSTI2018

Cyber–Physical Correlation Effects in Defense Games for Large Discrete Infrastructures

In certain critical infrastructures, correlations between cyber and physical components can be exploited to launch strategic attacks, so that disruptions to one component may affect others and possibly the entire infrastructure. Such correlations must be explicitly taken into account in ensuring the survival of the infrastructure. For large discrete infrastructures characterized by the number of cyber and physical components, we characterize the cyber–physical interactions at two levels: (i) the cyber–physical failure correlation function specifies the conditional survival probability of the cyber sub-infrastructure given that of the physical sub-infrastructure (both specified by their marginal probabilities), and (ii) individual survival probabilities of both sub-infrastructures are characterized by first-order differential conditions expressed in terms of their multiplier functions. We formulate an abstract problem of ensuring the survival probability of a cyber–physical infrastructure with discrete components as a game between the provider and attacker, whose utility functions are composed of infrastructure survival probability terms and cost terms, both expressed in terms of the number of components attacked and reinforced. We derive Nash equilibrium conditions and sensitivity functions that highlight the dependence of infrastructure survival probability on cost terms, correlation functions, multiplier functions, and sub-infrastructure survival probabilities. We apply these analytical results to characterize the defense postures of simplified models of metro systems, cloud computing infrastructures, and smart power grids.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1990

Frame synchronization for a channel with different data rates

In some communication links, different data rates are used at different times on the same channel and each is assigned a separate frame synchronization (sync) marker. The frame search is completed within a finite number of frame intervals once both phase-lock and bit-sync are achieved because the incoming data rate can change. The probability of failure to frame-lock, and an upper bound of probability of false frame-lock in a frame search mode are derived for a communication link with a 10 to the -5th bit error probability. The probability of failure to detect data-rate change and the probability of false detection of data-rate change in a frame-lock mode are discussed. In a frame-lock mode, the probability of failure to detect data-rate change and probability detection of data-rate change an equal to 4.96 x 10 to the -12th. A frame-lock device can be used for the Space Station Freedom Ku-band return link.

Kwon, Hyuck M.↗

NASA NTRS1992

Inverse sequential detection of parameter changes in developing time series

Progressive values of two probabilities are obtained for parameter estimates derived from an existing set of values and from the same set enlarged by one or more new values, respectively. One probability is that of erroneously preferring the second of these estimates for the existing data ('type 1 error'), while the second probability is that of erroneously accepting their estimates for the enlarged test ('type 2 error'). A more stable combined 'no change' probability which always falls between 0.5 and 0 is derived from the (logarithmic) width of the uncertainty region of an equivalent 'inverted' sequential probability ratio test (SPRT, Wald 1945) in which the error probabilities are calculated rather than prescribed. A parameter change is indicated when the compound probability undergoes a progressive decrease. The test is explicitly formulated and exemplified for Gaussian samples.

Radok, Uwe↗

NASA NTRS1999

Kennedy Space Center Coronary Heart Disease Risk Screening Program

The number one cause of death in the U.S. is coronary heart disease (CHD). It is probably a major cause of death and disability in the lives of employees at Kennedy Space Center (KSC) as well. The KSC Biomedical Office used a multifactorial mathematical formula from the Framingham Heart Study to calculate CHD risk probabilities for individuals in a segment of the KSC population that required medical evaluation for job certification. Those assessed to be high-risk probabilities will be targeted for intervention. Every year, several thousand KSC employees require medical evaluations for job related certifications. Most medical information for these evaluations is gathered on-site at one of the KSC or Cape Canaveral Air Station (CCAS) medical clinics. The formula used in the Framingham Heart Study allows calculation of a person's probability of acquiring CHD within 10 years. The formula contains the following variables: Age, Diabetes, Smoking, Left Ventricular Hypertrophy, Blood Pressure (Systolic or Diastolic), Cholesterol, and HDL cholesterol. The formula is also gender specific. It was used to calculate the 10-year probabilities of CHD in KSC employees who required medical evaluations for job certifications during a one-year time frame. This KSC population was profiled and CHD risk reduction interventions could be targeted to those at high risk. Population risk could also be periodically reevaluated to determine the effectiveness of intervention. A 10-year CHD risk probability can be calculated for an individual quite easily while gathering routine medical information. An employee population's CHD risk probability can be profiled graphically revealing high risk segments of the population which can be targeted for risk reduction intervention. The small audience of NASA/contractor physicians, nurses and exercise/fitness professionals at the breakout session received the lecture very well. Approximately one third indicated by a show of hands that they would be interested in implementing a similar program at their NASA Center. Questions were asked pertaining to standardization for age, the validity of using the idealized male values also for the female population, and indications of the screening test's sensitivity and specificity.

Tipton, David A.↗

NASA NTRS2005

Bidirectional Classical Stochastic Processes with Measurements and Feedback

A measurement on a quantum system is said to cause the "collapse" of the quantum state vector or density matrix. An analogous collapse occurs with measurements on a classical stochastic process. This paper addresses the question of describing the response of a classical stochastic process when there is feedback from the output of a measurement to the input, and is intended to give a model for quantum-mechanical processes that occur along a space-like reaction coordinate. The classical system can be thought of in physical terms as two counterflowing probability streams, which stochastically exchange probability currents in a way that the net probability current, and hence the overall probability, suitably interpreted, is conserved. The proposed formalism extends the . mathematics of those stochastic processes describable with linear, single-step, unidirectional transition probabilities, known as Markov chains and stochastic matrices. It is shown that a certain rearrangement and combination of the input and output of two stochastic matrices of the same order yields another matrix of the same type. Each measurement causes the partial collapse of the probability current distribution in the midst of such a process, giving rise to calculable, but non-Markov, values for the ensuing modification of the system's output probability distribution. The paper concludes with an analysis of a classical probabilistic version of the so-called grandfather paradox.

Hahne, G. E.↗

NASA NTRS2001

Spacecraft attitude control systems with dynamic methods and structures for processing star tracker signals

Methods are provided for dynamically processing successively-generated star tracker data frames and associated valid flags to generate processed star tracker signals that have reduced noise and a probability greater than a selected probability P.sub.slctd of being valid. These methods maintain accurate spacecraft attitude control in the presence of spurious inputs (e.g., impinging protons) that corrupt collected charges in spacecraft star trackers. The methods of the invention enhance the probability of generating valid star tracker signals because they respond to a current frame probability P.sub.frm by dynamically selecting the largest valid frame combination whose combination probability P.sub.cmb satisfies a selected probability P.sub.slctd. Noise is thus reduced while the probability of finding a valid frame combination is enhanced. Spacecraft structures are also provided for practicing the methods of the invention.

Liu, Yong↗

NASA NTRS2008

Completion of the Edward Air Force Base Statistical Guidance Wind Tool

The goal of this task was to develop a GUI using EAFB wind tower data similar to the KSC SLF peak wind tool that is already in operations at SMG. In 2004, MSFC personnel began work to replicate the KSC SLF tool using several wind towers at EAFB. They completed the analysis and QC of the data, but due to higher priority work did not start development of the GUI. MSFC personnel calculated wind climatologies and probabilities of 10-minute peak wind occurrence based on the 2-minute average wind speed for several EAFB wind towers. Once the data were QC'ed and analyzed the climatologies were calculated following the methodology outlined in Lambert (2003). The climatologies were calculated for each tower and month, and then were stratified by hour, direction (10" sectors), and direction (45" sectors)/hour. For all climatologies, MSFC calculated the mean, standard deviation and observation counts of the Zminute average and 10-minute peak wind speeds. MSFC personnel also calculated empirical and modeled probabilities of meeting or exceeding specific 10- minute peak wind speeds using PDFs. The empirical PDFs were asymmetrical and bounded on the left by the 2- minute average wind speed. They calculated the parametric PDFs by fitting the GEV distribution to the empirical distributions. Parametric PDFs were calculated in order to smooth and interpolate over variations in the observed values due to possible under-sampling of certain peak winds and to estimate probabilities associated with average winds outside the observed range. MSFC calculated the individual probabilities of meeting or exceeding specific 10- minute peak wind speeds by integrating the area under each curve. The probabilities assist SMG forecasters in assessing the shuttle FR for various Zminute average wind speeds. The A M ' obtained the processed EAFB data from Dr. Lee Bums of MSFC and reformatted them for input to Excel PivotTables, which allow users to display different values with point-click-drag techniques. The GUI was created from the PivotTables using VBA code. It is run through a macro within Excel and allows forecasters to quickly display and interpret peak wind climatology and probabilities in a fast-paced operational environment. The GUI was designed to look and operate exactly the same as the KSC SLF tool since SMG forecasters were already familiar with that product. SMG feedback was continually incorporated into the GUI ensuring the end product met their needs. The final version of the GUI along with all climatologies, PDFs, and probabilities has been delivered to SMG and will be put into operational use.

Dreher, Joseph G.↗