Engineering Papers⌕ Search

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.

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

At least 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Evaluation of Behavior and Survival of Fish Exposed to an Axial-Flow Hydrokinetic Turbine

Abstract Previous studies have evaluated fish injury and mortality at hydrokinetic (HK) turbines, but because these studies focused on the impacts of these turbines in situ they were unable to evaluate fish responses to controlled environmental characteristics (e.g., current velocity and light or dark conditions). In this study, we used juvenile hybrid Striped Bass (HSB; Striped Bass Morone saxatilis × White Bass M. chrysops; N = 620), Rainbow Trout Oncorhynchus mykiss (N = 3,719), and White Sturgeon Acipenser transmontanus (N = 294) in a series of laboratory experiments to (1) evaluate the ability of fish to avoid entrainment through an axial-flow HK turbine, (2) evaluate fish injury and survival associated with turbine entrainment, and (3) compare the effects of different HK turbines on fish. We found that the probability of turbine entrainment was species dependent and highest for HSB. Across species, current velocity influenced entrainment probability. Among entrained fish, observed survival rates were generally >0.95. The probability of injury for surviving entrained fish only differed from that for nonentrained fish for Rainbow Trout and in general was not >0.20. The probability of injury following entrainment was greater only for HSB, although there were no differences in injury rates between fish that were turbine entrained and those that were not, suggesting that injuries were not turbine related. Taking turbine entrainment, survival, and injury estimates together allowed us to estimate the probability of a randomly selected fish in a population proximate to an HK turbine surviving passage or remaining uninjured after passage. For species and current velocities for which there was a significant effect due to entrainment, we estimated, for instance, that HSB had a survival probability of 0.95 and that Rainbow Trout and White Sturgeon had a >0.99 probability of survival. Similarly, by combining these estimates with those from previous studies, we derived total passage survival probabilities >0.90 but generally approaching 1.00 across different HK turbine types, fish species, and fish lengths. Received May 14, 2014; accepted October 23, 2014

Amaral, Stephen V.↗
NASA NTRS2006

Convective Weather Forecast Quality Metrics for Air Traffic Management Decision-Making

Since numerical weather prediction models are unable to accurately forecast the severity and the location of the storm cells several hours into the future when compared with observation data, there has been a growing interest in probabilistic description of convective weather. The classical approach for generating uncertainty bounds consists of integrating the state equations and covariance propagation equations forward in time. This step is readily recognized as the process update step of the Kalman Filter algorithm. The second well known method, known as the Monte Carlo method, consists of generating output samples by driving the forecast algorithm with input samples selected from distributions. The statistical properties of the distributions of the output samples are then used for defining the uncertainty bounds of the output variables. This method is computationally expensive for a complex model compared to the covariance propagation method. The main advantage of the Monte Carlo method is that a complex non-linear model can be easily handled. Recently, a few different methods for probabilistic forecasting have appeared in the literature. A method for computing probability of convection in a region using forecast data is described in Ref. 5. Probability at a grid location is computed as the fraction of grid points, within a box of specified dimensions around the grid location, with forecast convection precipitation exceeding a specified threshold. The main limitation of this method is that the results are dependent on the chosen dimensions of the box. The examples presented Ref. 5 show that this process is equivalent to low-pass filtering of the forecast data with a finite support spatial filter. References 6 and 7 describe the technique for computing percentage coverage within a 92 x 92 square-kilometer box and assigning the value to the center 4 x 4 square-kilometer box. This technique is same as that described in Ref. 5. Characterizing the forecast, following the process described in Refs. 5 through 7, in terms of percentage coverage or confidence level is notionally sound compared to characterizing in terms of probabilities because the probability of the forecast being correct can only be determined using actual observations. References 5 through 7 only use the forecast data and not the observations. The method for computing the probability of detection, false alarm ratio and several forecast quality metrics (Skill Scores) using both the forecast and observation data are given in Ref. 2. This paper extends the statistical verification method in Ref. 2 to determine co-occurrence probabilities. The method consists of computing the probability that a severe weather cell (grid location) is detected in the observation data in the neighborhood of the severe weather cell in the forecast data. Probabilities of occurrence at the grid location and in its neighborhood with higher severity, and with lower severity in the observation data compared to that in the forecast data are examined. The method proposed in Refs. 5 through 7 is used for computing the probability that a certain number of cells in the neighborhood of severe weather cells in the forecast data are seen as severe weather cells in the observation data. Finally, the probability of existence of gaps in the observation data in the neighborhood of severe weather cells in forecast data is computed. Gaps are defined as openings between severe weather cells through which an aircraft can safely fly to its intended destination. The rest of the paper is organized as follows. Section II summarizes the statistical verification method described in Ref. 2. The extension of this method for computing the co-occurrence probabilities in discussed in Section HI. Numerical examples using NCWF forecast data and NCWD observation data are presented in Section III to elucidate the characteristics of the co-occurrence probabilities. This section also discusses the procedure for computing throbabilities that the severity of convection in the observation data will be higher or lower in the neighborhood of grid locations compared to that indicated at the grid locations in the forecast data. The probability of coverage of neighborhood grid cells is also described via examples in this section. Section IV discusses the gap detection algorithm and presents a numerical example to illustrate the method. The locations of the detected gaps in the observation data are used along with the locations of convective weather cells in the forecast data to determine the probability of existence of gaps in the neighborhood of these cells. Finally, the paper is concluded in Section V.

Chatterji, Gano B.↗
NASA NTRS2012

Communicating the Threat of a Tropical Cyclone to the Eastern Range

The 45th Weather Squadron (45 WS) has developed a tool to help visualize the Wind Speed Probability product from the National Hurricane Center (NHC) and to help communicate that information to space launch customers and decision makers at the 45th Space Wing (45 SW) and Kennedy Space Center (KSC) located in east central Florida. This paper reviews previous work and presents the new visualization tool, including initial feedback as well as the pros and cons. The NHC began issuing their Wind Speed Probability product for tropical cyclones publicly in 2006. The 45 WS uses this product to provide a threat assessment to 45 SW and KSC leadership for risk evaluations with an approaching tropical cyclone. Although the wind speed probabilities convey the uncertainty of a tropical cyclone well, communicating this information to customers is a challenge. The 45 WS continually strives to provide the wind speed probability information to customers in a context which clearly communicates the threat of a tropical cyclone. First, an intern from the Florida Institute of Technology (FIT) Atmospheric Sciences department, sponsored by Scitor Corporation, independently evaluated the NHC wind speed probability product. This work was later extended into a M.S. thesis at FIT, partially funded by Scitor Corporation and KSC. A second thesis at FIT further extended the evaluation partially funded by KSC. Using this analysis, the 45 WS categorized the probabilities into five probability interpretation categories: Very Low, Low, Moderate, High, and Very High. These probability interpretation categories convert the forecast probability and forecast interval into easily understood categories that are consistent across all ranges of probabilities and forecast intervals. As a follow-on project, KSC funded a summer intern to evaluate the human factors of the probability interpretation categories, which ultimately refined some of the thresholds. The 45 WS created a visualization tool to express the timing and risk for multiple locations in a single graphic. Preliminary results on an on-going project by FIT will be included in this paper. This project is developing a new method of assigning the probability interpretation categories and updating the evaluation of the performance of the NHC wind speed probability analysis.

Winters, Katherine A.↗
DOE OSTI2016

Probabilistic Risk Assessment of Spent Fuel Pool Cooling System Failure Considering Safety and Security Initiating Events

In order to quantify the probable risk associated with a safety system, nuclear safety risk evaluation process is applied which is generally performed in three steps. The first step is to calculate the damage frequency of a critical system. This step is followed by the radioactive source term estimation due to the potential system damage. The final step is the assessment of off-site consequences. The safety risk formula is shown by Eq. (1). P{sub F} is the probability of safety system failure (the probability used describes a random, independent event), and C is the consequence value. R{sub SAF} = P{sub F} * C (1) The security risk associated with an attack on critical infrastructure and/or special nuclear materials theft is usually evaluated using three- term risk equation (probability of attack, probability of adversary success, and consequence). A security threat and the related initiating events are not considered in a pure safety risk analysis. The vulnerability of the physical protection system is performed independently even though the consequence resulting from a security initiating event can be similar to a safety related initiating event and thus gets double counted. The security risk can be calculated using Eq. (2) R{sub SEC} = P{sub A} * P{sub S} * C (2) PA is the probability of attack, P{sub S} is the probability of adversary success (P{sub S} = 1 - P{sub I} P{sub N}). P{sub I} being the probability of interruption which accounts for the probabilities of detection of all timely detection elements along an adversary path. PN being the probability of neutralization: the result of engagement of the response force with the adversaries after the interruption has occurred. C is the consequence magnitude from a successful attack. The two security probabilities do not describe random, independent events. As a result, because of the differences in the randomness and independence of the respective safety and security probabilities, the safety and the security risks are currently calculated separately. The consequence term in nuclear safety and security risk evaluation have the same definition. An effort is made here to study the effect on the overall risk by considering security initiating events along with the subsequent failure of the intermediate safety element. The target selected is a spent fuel storage pool for analyzing the public risk arising from the interplay between nuclear security and safety initiating events. EASI (Estimate of Adversary Sequence Interruption) model is used to calculate the probability of interruption related to a security system that leads to the calculation of the security initiating event frequency. Further progression of this initiating security event to a top event is analyzed using a Probabilistic Safety Assessment (PSA) tool, SAPHIRE (version 8) to obtain a joint pool damage frequency. For the progression of events leading to a top event, the failure of the intermediate safety elements is assumed to be external random events such as heavy load drop and fire events. The objectives of this study is to (1) determine the double accounting of frequency of failure, risk and allocated resources in the scenario of a non-random security initiating event leading to a safety random event chain (2) compare the damage frequency from pure security with the frequency of combined security-safety failure event and pure safety failure events. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

WE-E-BRE-11: New Method to Simulate DNA Damage Using Ionization Cross-Sections and a Geometrical Nucleosome Model

Purpose: To obtain probability distributions of various DNA damage types as a function of the incident electron kinetic energy. Methods: Using Geant4-DNA electron ionization cross-sections, we calculated path length distributions for electrons of energies between 10 eV and 1 MeV, defined as the length between two subsequent interactions. These path lengths were then convolved with probability distributions for the creation of same-strand damage, opposite-strand damage, clustered damage, isolated damage, and same DNA strand target damage. These probability distributions of DNA damage were obtained by a Monte Carlo routine calculating probabilities of interaction in DNA targets inside a nucleosome geometrical model. Results represent the probability of a secondary electron, initially created inside a DNA strand target, of undergoing its next interaction: (1) in the opposite strand (DSB), (2) in the same strand (SSB+), (3) in either the opposite or same-strand (clustered), (4) in the same DNA target (multiple-hit) or (5) more than 10 base pairs away (isolated). Results: Electrons with kinetic energy between 50 and 250 eV have a maximal probability of creating DSB, SSB+, clustered damage and multiple-hits in the same target The probabilities for these damage patterns have values of 2.5%, 4.3%, 6.7% and 5.4%, respectively. Isolated damage is most probable between 700 eV to 900 eV with a probability of 0.2%. Conclusion: We obtained DNA damage probability distributions as a function of electron incident energy. We showed that electrons with kinetic energies between 50 and 250 eV have the highest probability of producing complex forms of DNA damage (DSB, SSB+). We also showed that a double ionization within the same DNA target is the most frequent outcome occurring 5% of the time. It is expected that electron slowing down spectra can be convolved with our formalism to calculate source specific DNA damage patterns. Research grants from governments of Canada and Quebec. PP acknowledges partial support by the CREATE Medical Physics Research Training Network grant of the Natural Sciences and Engineering Research Council (Grant number: 432290)

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

ESTIMATING THE PROPERTIES OF HARD X-RAY SOLAR FLARES BY CONSTRAINING MODEL PARAMETERS

We wish to better constrain the properties of solar flares by exploring how parameterized models of solar flares interact with uncertainty estimation methods. We compare four different methods of calculating uncertainty estimates in fitting parameterized models to Ramaty High Energy Solar Spectroscopic Imager X-ray spectra, considering only statistical sources of error. Three of the four methods are based on estimating the scale-size of the minimum in a hypersurface formed by the weighted sum of the squares of the differences between the model fit and the data as a function of the fit parameters, and are implemented as commonly practiced. The fourth method is also based on the difference between the data and the model, but instead uses Bayesian data analysis and Markov chain Monte Carlo (MCMC) techniques to calculate an uncertainty estimate. Two flare spectra are modeled: one from the Geostationary Operational Environmental Satellite X1.3 class flare of 2005 January 19, and the other from the X4.8 flare of 2002 July 23. We find that the four methods give approximately the same uncertainty estimates for the 2005 January 19 spectral fit parameters, but lead to very different uncertainty estimates for the 2002 July 23 spectral fit. This is because each method implements different analyses of the hypersurface, yielding method-dependent results that can differ greatly depending on the shape of the hypersurface. The hypersurface arising from the 2005 January 19 analysis is consistent with a normal distribution; therefore, the assumptions behind the three non-Bayesian uncertainty estimation methods are satisfied and similar estimates are found. The 2002 July 23 analysis shows that the hypersurface is not consistent with a normal distribution, indicating that the assumptions behind the three non-Bayesian uncertainty estimation methods are not satisfied, leading to differing estimates of the uncertainty. We find that the shape of the hypersurface is crucial in understanding the output from each uncertainty estimation technique, and that a crucial factor determining the shape of hypersurface is the location of the low-energy cutoff relative to energies where the thermal emission dominates. The Bayesian/MCMC approach also allows us to provide detailed information on probable values of the low-energy cutoff, E{sub c} , a crucial parameter in defining the energy content of the flare-accelerated electrons. We show that for the 2002 July 23 flare data, there is a 95% probability that E{sub c} lies below approximately 40 keV, and a 68% probability that it lies in the range 7-36 keV. Further, the low-energy cutoff is more likely to be in the range 25-35 keV than in any other 10 keV wide energy range. The low-energy cutoff for the 2005 January 19 flare is more tightly constrained to 107 {+-} 4 keV with 68% probability. Using the Bayesian/MCMC approach, we also estimate for the first time probability density functions for the total number of flare-accelerated electrons and the energy they carry for each flare studied. For the 2002 July 23 event, these probability density functions are asymmetric with long tails orders of magnitude higher than the most probable value, caused by the poorly constrained value of the low-energy cutoff. The most probable electron power is estimated at 10{sup 28.1} erg s{sup -1}, with a 68% credible interval estimated at 10{sup 28.1}-10{sup 29.0} erg s{sup -1}, and a 95% credible interval estimated at 10{sup 28.0}-10{sup 30.2} erg s{sup -1}. For the 2005 January 19 flare spectrum, the probability density functions for the total number of flare-accelerated electrons and their energy are much more symmetric and narrow: the most probable electron power is estimated at 10{sup 27.66{+-}0.01} erg s{sup -1} (68% credible intervals). However, in this case the uncertainty due to systematic sources of error is estimated to dominate the uncertainty due to statistical sources of error.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A Poisson method application to the assessment of the earthquake hazard in the North Anatolian Fault Zone, Turkey

North Anatolian Fault (NAF) is one from the most important strike-slip fault zones in the world and located among regions in the highest seismic activity. The NAFZ observed very large earthquakes from the past to present. The aim of this study; the important parameters of Gutenberg-Richter relationship (a and b values) estimated and this parameters taking into account, earthquakes were examined in the between years 1900-2015 for 10 different seismic source regions in the NAFZ. After that estimated occurrence probabilities and return periods of occurring earthquakes in fault zone in the next years, and is being assessed with Poisson method the earthquake hazard of the NAFZ. The Region 2 were observed the largest earthquakes for the only historical period and hasn’t been observed large earthquake for the instrumental period in this region. Two historical earthquakes (1766, M{sub S}=7.3 and 1897, M{sub S}=7.0) are included for Region 2 (Marmara Region) where a large earthquake is expected in the next years. The 10 different seismic source regions are determined the relationships between the cumulative number-magnitude which estimated a and b parameters with the equation of LogN=a-bM in the Gutenberg-Richter. A homogenous earthquake catalog for M{sub S} magnitude which is equal or larger than 4.0 is used for the time period between 1900 and 2015. The database of catalog used in the study has been created from International Seismological Center (ISC) and Boğazici University Kandilli observation and earthquake research institute (KOERI). The earthquake data were obtained until from 1900 to 1974 from KOERI and ISC until from 1974 to 2015 from KOERI. The probabilities of the earthquake occurring are estimated for the next 10, 20, 30, 40, 50, 60, 70, 80, 90 and 100 years in the 10 different seismic source regions. The highest earthquake occur probabilities in 10 different seismic source regions in the next years estimated that the region Tokat-Erzincan (Region 9) %99 with an earthquake occur probability for magnitude 6.5 which the return period 24.7 year, %92 with an earthquake occur probability for magnitude 7 which the return period 39.1 year, %80 with an earthquake occur probability for magnitude 7.5 which the return period 62.1 year, %64 with an earthquake occur probability for magnitude 8 which the return period 98.5 year. For the Marmara Region (Region 2) in the next 100 year estimated that %89 with an earthquake occur probability for magnitude 6 which the return period 44.9 year, %45 with an earthquake occur probability for magnitude 6.5 which the return period 87 year, %45 with an earthquake occur probability for magnitude 7 which the return period 168.6 year.

58 GEOSCIENCES↗
DOE OSTI2013

Environmental Impacts of Wind Power Development on the Population Biology of Greater Prairie-Chickens

1. We investigated the impacts of wind power development on the demography, movements, and population genetics of Greater Prairie-Chickens (Tympanuchus cupido) at three sites in northcentral and eastern Kansas for a 7-year period. Only 1 of 3 sites was developed for wind power, the 201MW Meridan Way Wind Power Facility at the Smoky Hills site in northcentral Kansas. Our project report is based on population data for prairie chickens collected during a 2-year preconstruction period (2007-2008), a 3-year postconstruction period (2009-2011) and one final year of lek surveys (2012). Where relevant, we present preconstruction data from our field studies at reference sites in the northern Flint Hills (2007-2009) and southern Flint Hills (2006-2008). 2. We addressed seven potential impacts of wind power development on prairie chickens: lek attendance, mating behavior, use of breeding habitat, fecundity rates, natal dispersal, survival rates, and population numbers. Our analyses of pre- and postconstruction impacts are based on an analysis of covariance design where we modeled population performance as a function of treatment period, distance to eventual or actual site of the nearest wind turbine, and the interaction of these factors. Our demographic and movement data from the 6-year study period at the Smoky Hills site included 23 lek sites, 251 radio-marked females monitored for 287 bird-years, and 264 nesting attempts. Our genetic data were based on genotypes of 1,760 females, males and chicks that were screened with a set of 27 microsatellite markers that were optimized in the lab. 3. In our analyses of lek attendance, the annual probability of lek persistence during the preconstruction period was ~0.9. During the postconstruction period, distance to nearest turbine did not have a significant effect on the probability of lek persistence. However, the probability of lek persistence increased from 0.69 at 0 m to 0.89 at 30 km from turbines, and most abandoned lek sites were located <5 km from turbines. Probability of lek persistence was significantly related to habitat and number of males. Leks had a higher probability of persistence in grasslands than agricultural fields, and increased from ~0.2 for leks of 5 males, to >0.9 for leks of 10 or more males. Large leks in grasslands should be a higher priority for conservation. Overall, wind power development had a weak effect on the annual probability of lek persistence. 3. We used molecular methods to investigate the mating behavior of prairie chickens. The prevailing view for lek-mating grouse is that females mate once to fertilize the clutch and that conspecific nest parasitism is rare. We found evidence that females mate multiple times to fertilize the clutch (8-18% of broods, 4-38% of chicks) and will parasitize nests of other females during egg-laying (~17% of nests). Variable rates of parentage were highest in the fragmented landscapes at the Smoky Hills field site, and were lower at the Flint Hills field site. Comparisons of the pre- and postconstruction periods showed that wind energy development did not affect the mating behaviors of prairie chickens. 4. We examined use of breeding habitats by radio-marked females and conducted separate analyses for nest site selection, and movements of females not attending nests or broods. The landscape was a mix of native prairie and agricultural habitats, and nest site selection was not random because females preferred to nest in grasslands. Nests tended to be closer to turbines during the postconstruction period and there was no evidence of behavioral avoidance of turbines by females during nest site selection. Movements of females not attending nests or broods showed that females crossed the site of the wind power development at higher rates during the preconstruction period (20%) than the postconstruction period (11%), and that movements away from turbines were more frequent during the postconstruction period. Thus, wind power development appears to affect movements in breeding habitats but not nest site selection of female prairie chickens during the breeding season. 5. We tested the effects of wind power development on five components of female fecundity: timing of clutch initiation, clutch size of first nests and renests, nest survival, and hatchability of eggs. Average date of clutch initiation was 26 April, clutch size was 12.7 and 10.6 eggs for first nests and renests, probability of nest survival was low at 0.18, but egg hatchability was high at 0.79. Wind power development had no impact on reproductive effort or nesting success, and all five components of fecundity were not related to treatment period or distance to turbine. Nest survival was the main factor limiting reproductive output of female prairie chickens and most losses were due to predation. Daily nest survival was strongly related to vegetative cover at the nest. Changes to rangeland management practices that would double nesting cover from 2.5 to 5 dm would triple the probability of nest survival from 0.17 to 0.52. Grass and forb cover had weak positive effects on daily nest survival whereas shrub cover, proximity to woodlands, and recent rainfall had negative effects. Reproductive performance of prairie chickens is low in managed rangelands in northcentral Kansas and efforts to improve range conditions and reduce predator activity would aid recovery of prairie chicken populations. 6. We used molecular methods to investigate patterns of natal dispersal in prairie chickens. High rates of nest failure limited the number of young that we could sample. Direct detections of natal dispersal were limited because survival of newly hatched chicks to become adults were low and because we were unable to detect dispersal distances outside of our study area. Direct observations of natal movements were limited and were inadequate to make conclusions about the potential impacts of wind energy development on natal dispersal. Spatial correlograms of genetic distance among males at leks were a more sensitive measure of population structure, and indicated a weak effect of wind energy development on the spatial genetic structure of prairie chickens. 7. We tested the effects of wind power development on female survival with time-to-event models, and on residual body mass of males with analyses of covariance. Distance to turbine and the interaction of distance and treatment period had no effect on female survival. Contrary to predictions of negative impacts of wind power development, the probability of female was lowest during the preconstruction period (0.274) and increased significantly during the postconstruction period (0.543). Inspection of hazard functions indicated that the difference in annual survival could be attributed to a higher risk of mortality during the lekking season in the preconstruction period. We suggest that wind power development may have improved ecological conditions for prairie chickens by disrupting the foraging behavior of diurnal raptors that kill prairie chickens at lek sites. In support of this idea, raptor kills tended to be farther from turbines during the postconstruction period whereas mammalian kills were closer. Analyses of the major causes of mortality did not support our hypothesis because the odds of raptor predation were greater after development was completed. Most mortality losses of radio-marked females were due to predation, and losses to collision mortality or harvest were rare events. Low rates of natural mortality during fall and winter imply that harvest is likely to be additive mortality in prairie chickens. Wind energy development reduced the residual body mass of male Greater Prairie-Chickens at lek sites near turbines. Low values of residual body mass could have a negative impact on individual survival or fecundity rates, or may be related to predation risk and flight performance of males displaying at open lek sites. 8. We tested for impacts of wind power development on population numbers of prairie chickens with monitoring of male numbers at leks, and with genetic measures of population structure. Lek counts indicated that wind power development did not affect the population size of prairie chickens. Peak counts of males at leks were recorded the first year after construction was completed and the highest rates of population change were observed during the interval when the wind power facility was constructed. Population numbers of prairie chickens near and distant from turbines appeared to covary in parallel, probably because bird numbers were linked by dispersal movements. Estimates of population viability based on genetic diversity, effective population size and rates of population exchange did not show annual changes and were unaffected by wind development during our study. Estimates of relatedness among males at the same and different leks suggested that wind power development has either reduced dispersal rates or changed settlement patterns, leading to higher rates of relatedness among males displaying at the same lek site. 9. Greater Prairie-Chickens were not strongly affected by wind power development in Kansas. Negative impacts of wind power development included a trend for reductions in lek persistence near turbines, behavioral avoidance of turbines by females during their breeding season movements, and changes in the genetic structure of males at leks that were consistent with reduced dispersal or recruitment rates. We found no impacts of wind power development on nest site selection, female reproductive effort or nesting success, or population numbers. Positive impacts of wind power development included an increase in female survival rates. We hypothesized that the unexpected increase in female survival was related to changes in trophic interactions and disruption of the foraging behavior of raptors that kill prairie chickens at lek site. 10. Research funding for this project included a grant from the 20% Wind by 2030 Program of the Department of Energy (this final report), and grants from the Kansas Department of Wildlife, Parks and Tourism, the National Fish and Wildlife Federation, and initial funding from the National Wind Coordinating Collaborative. Research products from data collected during the preconstruction period have included five peer-reviewed research articles and two PhD dissertations at Kansas State University. Additional manuscripts are in review for possible publication in 2013-2014.

17 WIND ENERGY↗
NASA NTRS1999

Structural Reliability Analysis and Optimization: Use of Approximations

This report is intended for the demonstration of function approximation concepts and their applicability in reliability analysis and design. Particularly, approximations in the calculation of the safety index, failure probability and structural optimization (modification of design variables) are developed. With this scope in mind, extensive details on probability theory are avoided. Definitions relevant to the stated objectives have been taken from standard text books. The idea of function approximations is to minimize the repetitive use of computationally intensive calculations by replacing them with simpler closed-form equations, which could be nonlinear. Typically, the approximations provide good accuracy around the points where they are constructed, and they need to be periodically updated to extend their utility. There are approximations in calculating the failure probability of a limit state function. The first one, which is most commonly discussed, is how the limit state is approximated at the design point. Most of the time this could be a first-order Taylor series expansion, also known as the First Order Reliability Method (FORM), or a second-order Taylor series expansion (paraboloid), also known as the Second Order Reliability Method (SORM). From the computational procedure point of view, this step comes after the design point identification; however, the order of approximation for the probability of failure calculation is discussed first, and it is denoted by either FORM or SORM. The other approximation of interest is how the design point, or the most probable failure point (MPP), is identified. For iteratively finding this point, again the limit state is approximated. The accuracy and efficiency of the approximations make the search process quite practical for analysis intensive approaches such as the finite element methods; therefore, the crux of this research is to develop excellent approximations for MPP identification and also different approximations including the higher-order reliability methods (HORM) for representing the failure surface. This report is divided into several parts to emphasize different segments of the structural reliability analysis and design. Broadly, it consists of mathematical foundations, methods and applications. Chapter I discusses the fundamental definitions of the probability theory, which are mostly available in standard text books. Probability density function descriptions relevant to this work are addressed. In Chapter 2, the concept and utility of function approximation are discussed for a general application in engineering analysis. Various forms of function representations and the latest developments in nonlinear adaptive approximations are presented with comparison studies. Research work accomplished in reliability analysis is presented in Chapter 3. First, the definition of safety index and most probable point of failure are introduced. Efficient ways of computing the safety index with a fewer number of iterations is emphasized. In chapter 4, the probability of failure prediction is presented using first-order, second-order and higher-order methods. System reliability methods are discussed in chapter 5. Chapter 6 presents optimization techniques for the modification and redistribution of structural sizes for improving the structural reliability. The report also contains several appendices on probability parameters.

Grandhi, Ramana V.↗
DOE OSTI2014

Analytical solutions for quantum walks on 1D chain with different shift operators

In this paper, we study the discrete-time quantum walks on 1D Chain with the moving and swapping shift operators. We derive analytical solutions for the eigenvalues and eigenstates of the evolution operator U{sup -hat} using the Chebyshev polynomial technique, and calculate the long-time averaged probabilities for the two different shift operators respectively. It is found that the probability distributions for the moving and swapping shift operators display completely different characteristics. For the moving shift operator, the probability distribution exhibits high symmetry where the probabilities at mirror positions are equal. The probabilities are inversely proportional to the system size N and approach to zero as N→∞. On the contrary, for the swapping shift operator, the probability distribution is not symmetric, the probability distribution approaches to a power-law stationary distribution as N→∞ under certain coin parameter condition. We show that such power-law stationary distribution is determined by the eigenstates of the eigenvalues ±1 and calculate the intrinsic probability for different starting positions. Our findings suggest that the eigenstates corresponding to eigenvalues ±1 play an important role for the swapping shift operator. - Highlights: • QWs on 1D chain with the moving and swapping operators are studied for the first time. • We derive analytical results for the probability distribution for the two operators. •We compare the dynamics of QWs with two different shift operators. • We find the particular eigenvalues ±1 play an important role for the dynamics. • We use the Chebyshev technique to treat the problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Role of the site of synaptic competition and the balance of learning forces for Hebbian encoding of probabilistic Markov sequences

The majority of distinct sensory and motor events occur as temporally ordered sequences with rich probabilistic structure. Sequences can be characterized by the probability of transitioning from the current state to upcoming states (forward probability), as well as the probability of having transitioned to the current state from previous states (backward probability). Despite the prevalence of probabilistic sequencing of both sensory and motor events, the Hebbian mechanisms that mold synapses to reflect the statistics of experienced probabilistic sequences are not well understood. Here, we show through analytic calculations and numerical simulations that Hebbian plasticity (correlation, covariance, and STDP) with pre-synaptic competition can develop synaptic weights equal to the conditional forward transition probabilities present in the input sequence. In contrast, post-synaptic competition can develop synaptic weights proportional to the conditional backward probabilities of the same input sequence. We demonstrate that to stably reflect the conditional probability of a neuron’s inputs and outputs, local Hebbian plasticity requires balance between competitive learning forces that promote synaptic differentiation and homogenizing learning forces that promote synaptic stabilization. The balance between these forces dictates a prior over the distribution of learned synaptic weights, strongly influencing both the rate at which structure emerges and the entropy of the final distribution of synaptic weights. Together, these results demonstrate a simple correspondence between the biophysical organization of neurons, the site of synaptic competition, and the temporal flow of information encoded in synaptic weights by Hebbian plasticity while highlighting the utility of balancing learning forces to accurately encode probability distributions, and prior expectations over such probability distributions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Probabilistic inference in very large universes

Our current favored cosmological theories allow for the striking and controversial possibility that the observable universe is just a small part of a much larger universe in which parameters that describe the effective, low-energy laws of physics vary from one region to another. The controversy is largely driven by the fact that such a “very large universe” is mostly observationally inaccessible to us, so the issue arises of how we can reasonably assess a theory that describes such a universe. In this paper, we propose a Bayesian method for theory assessment based on theory-generated probability distributions for our observations. We focus on the principles that define this method, leaving aside concerns about how, in practice, one would carry out the required calculations. (One important issue that we set aside is the measure problem.) We argue that cosmological theories can be tested by the standard method of Bayesian updating, but we need to use theoretical predictions for “first-person” probabilities—that is, probabilities that we should use for our observations, taking into account all relevant selection effects. These selection effects can vary from one observer to another and can vary with time, so, in principle, first-person probabilities are defined for each observer instant—an observer at a specific instant of time. Calculations of first-person probabilities should take into account everything that the observer believes about herself and her surroundings, which we refer to as her subjective state. If the universe is very large, a theory might predict that there are many observer instants in the same subjective state; we argue that first-person probabilities should be calculated using a principle of self-locating indifference (PSLI), the assumption that any real observer should make predictions for her future as if she were chosen randomly and uniformly from the theoretically predicted observer instants that share her subjective state. We believe the PSLI is intuitively very reasonable, but we also argue that, if the theory is correct, the use of this principle maximizes the expected fraction of observers who will make correct predictions. A further complication is that cosmological theories are not expected to fully predict the detailed properties of the universe, but rather will predict a set of possible universes, each with a probability. Different possible universes will generically have different numbers of observers. We argue that, in the calculation of first-person probabilities, the probability for each possible universe should be weighted by the number of observer instants in the specified subjective state that it contains. These issues have been controversial in the literature, so we also provide a rebuttal to the claim that principles like the PSLI involve a “selection fallacy”; a rebuttal to what we dub the principle of required certainty; an argument rejecting theories that predict a preponderance of Boltzmann brains; a rebuttal to a parable about humans and Jovians used by Hartle and Srednicki to argue that assumptions of typicality can lead to absurd consequences; and, finally, a discussion about how the use of “old evidence” can be fit into a Bayesian mold.

Azhar, Feraz [University of Notre Dame, IN (United↗
DOE OSTI2014

Factors Affecting Route Selection and Survival of Steelhead Kelts at Snake River Dams in 2012 and 2013

In 2012 and 2013, Pacific Northwest National Laboratory conducted a study that summarized the passage proportions and route-specific survival rates of steelhead kelts that passed through Federal Columbia River Power System (FCRPS) dams. To accomplish this, a total of 811 steelhead kelts were tagged with Juvenile Salmon Acoustic Telemetry System (JSATS) transmitters. Acoustic receivers, both autonomous and cabled, were deployed throughout the FCRPS to monitor the downstream movements of tagged-kelts. Kelts were also tagged with Passive Integrated Transponder tags to monitor passage through juvenile bypass systems and detect returning fish. The current study evaluated data collected in 2012 and 2013 to identify individual, behavioral, environmental and dam operation variables that were related to passage and survival of steelhead kelts that passed through FCRPS dams. Bayesian model averaging of multivariable logistic regression models was used to identify the environmental, temporal, operational, individual, and behavioral variables that had the highest probability of influencing the route of passage and the route-specific survival probabilities for kelts that passed Lower Granite (LGR), Little Goose (LGS), and Lower Monumental (LMN) dams in 2012 and 2013. The posterior probabilities of the best models for predicting route of passage ranged from 0.106 for traditional spill at LMN to 0.720 for turbine passage at LGS. Generally, the behavior (depth and near-dam searching activity) of kelts in the forebay appeared to have the greatest influence on their route of passage. Shallower-migrating kelts had a higher probability of passing via the weir and deeper-migrating kelts had a higher probability of passing via the JBS and turbines than other routes. Kelts that displayed a higher level of near-dam searching activity had a higher probability of passing via the spillway weir and those that did less near-dam searching had a higher probability of passing via the JBS and turbines. The side of the river in which kelts approached the dam and dam operations also affected route of passage. Dam operations and the size and condition of kelts were found to have the greatest effect on route-specific survival probabilities for fish that passed via the spillway at LGS. That is, longer kelts and those in fair condition had a lower probability of survival for fish that passed via the spillway weir. The survival of spillway weir- and deep-spill passed kelts was positively correlated with the percent of the total discharge that passed through turbine unit 4. Too few kelts passed through the traditional spill, JBS, and turbine units to evaluate survival through these routes. The information gathered in this study describes Snake River steelhead kelt passage behavior, rates, and distributions through the FCRPS as well as provide information to biologists and engineers about the dam operations and abiotic conditions that are related to passage and survival of steelhead kelts.

13 HYDRO ENERGY↗
NASA NTRS1999

Lunar Meteorites and Implications for Compositional Remote Sensing of the Lunar Surface

Lunar meteorites (LMs) are rocks found on Earth that were ejected from the Moon by impact of an asteroidal meteoroid. Three factors make the LMs important to remote-sensing studies: (1) Most are breccias composed of regolith or fragmental material; (2) all are rocks that resided (or breccias composed of material that resided) in the upper few meters of the Moon prior to launch and (3) most apparently come from areas distant from the Apollo sites. How Many Lunar Locations? At this writing (June 1999), there are 18 known lunar meteorite specimens. When unambiguous cases of terrestrial pairing are considered, the number of actual LMs reduces to 13. (Terrestrial pairing is when a single piece of lunar rock entered Earth's atmosphere, but multiple fragments were produced because the meteoroid broke apart on entry, upon hitting the ground or ice, or while being transported through the ice.) We have no reason to believe that LMs preferentially derive from any specific region(s) of the Moon; i.e., we believe that they are samples from random locations. However, we do not know how many different locations are represented by the LMs; mathematically, it could be as few as 1 or as many as 13. The actual maximum is < 13 because in some cases a single impact appears to have yielded more than one LM. Yamato 793169 and Asuka 881757 are considered "source-crater paired" or "launch paired" because they are compositionally and petrographically similar to each other and distinct from the others, and both have similar cosmic-ray exposure (CRE) histories. The same can be said of QUE 94281 and Y 793274. Thus the 13 meteorites probably represent a maximum of 11 locations on the Moon. The minimum number of likely source craters is debated and in flux as new data for different isotopic systems are obtained. Conservatively, considering CRE data only, a minimum of about 5 impacts is required. Compositional and petrographic data offer only probabilistic constraints. An extreme, but not unreasonable viewpoint, is that such data offer no constraint. For example, if one were to cut up the Apollo 17 landing site (which was selected for its diversity) into softball-sized pieces, some of those pieces (e.g., sample 70135) would be crystalline mare basalts like Y 793169 whereas others (e.g., sample 73131 would be feldspathic regolith breccias like MAC 88104/ 88105. However, nature is not so devious. Warren argues that LMs come from craters of only a few kilometers in diameter. If so, even though CRE data allow, for example, that ALHA 81005 and Y 791197) were launched simultaneously from the same crater, the probability is nevertheless low because the two meteorites are compositionally and mineralogically distinct. Thus, within the allowed range (5-11) for the number of locations represented by the LMs, values at the high end of the range are probably more likely. Mare Meteorites: Three LMs consist almost entirely of mare basalt. Two, Y 793169 and Asuka 881757, are unbrecciated, low-Ti, crystalline rocks that are compositionally and mineralogically similar (but not identical) to each other; they probably derive from a single lunar-mare location. The third, EET 87521/96008, is a fragmental breccia consisting predominantly of VLT mare basalt. Thus, these LMs probably represent only two lunar mare locations. The basaltic LMs have mineral and bulk compositions distinct from Apollo mare basalts. The petrography of Calcalong Creek has not been described in detail, but compositionally it is unique in that it corresponds to a mixture (breccia) of about one-half feldspathic material (i.e., the mean composition of the feldspathic lunar meteorites, below), one-fourth KREEP norite, one-fourth VLT mare basalt (like EET 87521), and 1% CI chondrite. With 4 micro g/g Th and correspondingly high concentrations of other incompatible elements, it is the only lunar meteorite that is likely to have come from within the Procellarum KREEP Terrane (PKT). Yamato 793274 and QUE 94281 are together distinct in being fragmental breccias containing subequal parts of feldspathic highland material and VLT mare basalt. Jolliff et al. estimate a mare to highland ratio of 54:46 for QUE 94281 and 62:38 for Y 793274; this difference is well within the range observed for soils collected only centimeters apart (in cores) at interface site like Apollo 15 and 17 [11]. Although the two meteorites were found on opposite sides of Antarctica, they are probably launch-paired. The strongest evidence is that the pyroclastic glass spherules that occur in both are of two compositional groups and the two groups are essentially the same in both meteorites. Yamato 791197 is nominally a feldspathic lunar meteorite (below), but among FLMs, it probably contains the highest abundance of clasts and glasses of mare derivation. As a consequence, its composition is at the high-Fe, low-Mg end of the range for FLMs and is not included in the FLM average of Table 1. Its composition is consistent with about 10% mare-derived material. Similarly, the two small (Y 82) pieces of Y 82192/82193186032 are more mafic than the large (Y 86) piece, probably as a result of about 7% mare-derived material. All Apollo missions went to areas in or near the PKT, and, consequently, all Apollo regolith samples are contaminated with Th-rich material from the PKT. At the nominally "typical" highland site, Apollo 16, about 30% of the regolith (<1-mm fines) is Th-rich ejecta from the Imbrium impact and about 6% is mare material probably derived from mare basins. Thus Apollo 16 regolith is not typical of the highlands. Among Apollo rocks, the compositions of the FLMs correspond most closely to the feldspathic granulitic breccias of Apollo 16 and 17. (Additional information is contained in original)

Korotev, R. L.↗
NASA NTRS2017

Sensitivity Analysis of the Bone Fracture Risk Model

Introduction: The probability of bone fracture during and after spaceflight is quantified to aid in mission planning, to determine required astronaut fitness standards and training requirements and to inform countermeasure research and design. Probability is quantified with a probabilistic modeling approach where distributions of model parameter values, instead of single deterministic values, capture the parameter variability within the astronaut population and fracture predictions are probability distributions with a mean value and an associated uncertainty. Because of this uncertainty, the model in its current state cannot discern an effect of countermeasures on fracture probability, for example between use and non-use of bisphosphonates or between spaceflight exercise performed with the Advanced Resistive Exercise Device (ARED) or on devices prior to installation of ARED on the International Space Station. This is thought to be due to the inability to measure key contributors to bone strength, for example, geometry and volumetric distributions of bone mass, with areal bone mineral density (BMD) measurement techniques. To further the applicability of model, we performed a parameter sensitivity study aimed at identifying those parameter uncertainties that most effect the model forecasts in order to determine what areas of the model needed enhancements for reducing uncertainty. Methods: The bone fracture risk model (BFxRM), originally published in (Nelson et al) is a probabilistic model that can assess the risk of astronaut bone fracture. This is accomplished by utilizing biomechanical models to assess the applied loads; utilizing models of spaceflight BMD loss in at-risk skeletal locations; quantifying bone strength through a relationship between areal BMD and bone failure load; and relating fracture risk index (FRI), the ratio of applied load to bone strength, to fracture probability. There are many factors associated with these calculations including environmental factors, factors associated with the fall event, mass and anthropometric values of the astronaut, BMD characteristics, characteristics of the relationship between BMD and bone strength and bone fracture characteristics. The uncertainty in these factors is captured through the use of parameter distributions and the fracture predictions are probability distributions with a mean value and an associated uncertainty. To determine parameter sensitivity, a correlation coefficient is found between the sample set of each model parameter and the calculated fracture probabilities. Each parameters contribution to the variance is found by squaring the correlation coefficients, dividing by the sum of the squared correlation coefficients, and multiplying by 100. Results: Sensitivity analyses of BFxRM simulations of preflight, 0 days post-flight and 365 days post-flight falls onto the hip revealed a subset of the twelve factors within the model which cause the most variation in the fracture predictions. These factors include the spring constant used in the hip biomechanical model, the midpoint FRI parameter within the equation used to convert FRI to fracture probability and preflight BMD values. Future work: Plans are underway to update the BFxRM by incorporating bone strength information from finite element models (FEM) into the bone strength portion of the BFxRM. Also, FEM bone strength information along with fracture outcome data will be incorporated into the FRI to fracture probability.

mathematical models↗
DOE OSTI2013

Preliminary Safety Analysis of the Gorleben Site: Safety Concept and Application to Scenario Development Based on a Site-Specific Features, Events and Processes (FEP) Database - 13304

Based upon the German safety criteria, released in 2010 by the Federal Ministry of the Environment (BMU), a safety concept and a safety assessment concept for the disposal of heat-generating high-level waste have both been developed in the framework of the preliminary safety case for the Gorleben site (Project VSG). The main objective of the disposal is to contain the radioactive waste inside a defined rock zone, which is called containment-providing rock zone. The radionuclides shall remain essentially at the emplacement site, and at the most, a small defined quantity of material shall be able to leave this rock zone. This shall be accomplished by the geological barrier and a technical barrier system, which is required to seal the inevitable penetration of the geological barrier by the construction of the mine. The safe containment has to be demonstrated for probable and less probable evolutions of the site, while evolutions with very low probability (less than 1 % over the demonstration period of 1 million years) need not to be considered. Owing to the uncertainty in predicting the real evolution of the site, plausible scenarios have been derived in a systematic manner. Therefore, a comprehensive site-specific features, events and processes (FEP) data base for the Gorleben site has been developed. The safety concept was directly taken into account, e.g. by identification of FEP with direct influence on the barriers that provide the containment. No effort was spared to identify the interactions of the FEP, their probabilities of occurrence, and their characteristics (values). The information stored in the data base provided the basis for the development of scenarios. The scenario development methodology is based on FEP related to an impairment of the functionality of a subset of barriers, called initial barriers. By taking these FEP into account in their probable characteristics the reference scenario is derived. Thus, the reference scenario describes a comprehensive set of probable future evolutions of the repository site. By stepwise consideration of less probable FEP or less probable characteristics of FEP within this process, alternative scenarios are also developed, which are characterized by a lower probability of occurrence. An important methodological aspect is that some assumptions had to be made for the scenario development. They allow, on the one hand, to deal systematically with incomplete knowledge regarding the geological situation below ground owing to restricted site investigations, and, on the other hand, to structure the scenario development process. The consideration of alternative assumptions may result in additional alternative scenarios. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

Single-shot, multi-signature remote detection of uranium by filament-induced breakdown spectroscopy

Filament-induced breakdown spectroscopy (FIBS) is an attractive remote detection approach that is potentially applicable to nuclear materials, such as uranium (U). Simultaneous detection of atomic U and uranium monoxide features is demonstrated in a single laser shot with detection probabilities of 52.2% and 39.0%, respectively, and with a 1% false-alarm probability. It is shown that a nearly 100% detection probability can be reached with 1% false alarm probability in .1 second when using an 80-Hz laser and data acquisition system. Atomic U and UO signal features are identifiable in the remote measurement with comparable signal-to-background ratios, suggesting that, despite being broader than the U I emission line, the UO band is also suitable for remote detection of uranium by FIBS.

(020.3260) Isotope shifts, (140.3440) Laser-induce↗
DOE OSTI2022

Cumulative effects of piscivorous colonial waterbirds on juvenile salmonids: A multi predator-prey species evaluation

We investigated the cumulative effects of predation by piscivorous colonial waterbirds on the survival of multiple salmonid ( Oncorhynchus spp.) populations listed under the U.S. Endangered Species Act (ESA) and determined what proportion of all sources of fish mortality (1 –survival) were due to birds in the Columbia River basin, USA. Anadromous juvenile salmonids (smolts) were exposed to predation by Caspian terns ( Hydroprogne caspia ), double-crested cormorants ( Nannopterum auritum ), California gulls ( Larus californicus ), and ring-billed gulls ( L . delawarensis ), birds known to consume both live and dead fish. Avian consumption and survival probabilities (proportion of available fish consumed or alive) were estimated for steelhead trout ( O . mykiss ), yearling Chinook salmon ( O . tshawytscha ), sub-yearling Chinook salmon, and sockeye salmon ( O . nerka ) during out-migration from the lower Snake River to the Pacific Ocean during an 11-year study period (2008–2018). Results indicated that probabilities of avian consumption varied greatly across salmonid populations, bird species, colony location, river reach, and year. Cumulative consumption probabilities (consumption by birds from all colonies combined) were consistently the highest for steelhead, with annual estimates ranging from 0.22 (95% credible interval = 0.20–0.26) to 0.51 (0.43–0.60) of available smolts. The cumulative effects of avian consumption were significantly lower for yearling and sub-yearling Chinook salmon, with consumption probabilities ranging annually from 0.04 (0.02–0.07) to 0.10 (0.07–0.15) and from 0.06 (0.3–0.09) to 0.15 (0.10–0.23), respectively. Avian consumption probabilities for sockeye salmon smolts was generally higher than for Chinook salmon smolts, but lower than for steelhead smolts, ranging annually from 0.08 (0.03–0.22) to 0.25 (0.14–0.44). Although annual consumption probabilities for birds from certain colonies were more than 0.20 of available smolts, probabilities from other colonies were less than 0.01 of available smolts, indicating that not all colonies of birds posed a substantial risk to smolt mortality. Consumption probabilities were lowest for small colonies and for colonies located a considerable distance from the Snake and Columbia rivers. Total mortality attributed to avian consumption was relatively small for Chinook salmon (less than 10%) but was the single greatest source of mortality for steelhead (greater than 50%) in all years evaluated. Results suggest that the potential benefits to salmonid populations of managing birds to reduce smolt mortality would vary widely depending on the salmonid population, the species of bird, and the size and location of the breeding colony.

Evans, Allen F.↗