Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Probability distributions”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

The Use of the Satellite Breakup Risk Assessment Model (SBRAM) to Characterize Collision Risk to Manned Spacecraft

NASA uses environment models such as ORDEM96 to characterize the long-term orbital debris collision hazard for spacecraft in LEO. Occasionally, however, there are breakups of satellites or rocket bodies that create enhanced collision hazard for a period of time. This enhanced collision hazard can pose increased risks to space operations - especially those involving manned missions where the tolerance for risk is very low. NASA has developed SBRAM to simulate the enhanced debris environment in the days and weeks that follow such a breakup. This simulation provides the kind of risk probabilities that can be used by mission planners to consider if changes are warranted for the mission. Announcements of breakups come to NASA from US Space Command as soon as they are identified. The pre-breakup orbit and time of breakup are used to determine the initial conditions of the explosion. SBRAM uses the latest explosion models developed at NASA to simulate a debris cloud for the breakup. The model uses a Monte Carlo technique to create a random debris cloud from the probability distributions in the breakup model. Each piece of debris randomly created in the cloud is propagated in a deterministic manner to include the effects of drag and other orbital perturbations. The detailed geometry of each simulated close approach to the target spacecraft is noted and logged and the collision probability is computed using an estimated probability density in down-range and cross-range positions of both the target spacecraft and debris object. The collision probability is computed from the overlap of these probability densities for each close-approach geometry and summed over all computed conjunctions. Cloud propagation runs over the desired time interval are then repeated until the scale of the collision risk can be estimated to a desired precision. This paper presents an overview of the SBRAM model and a number of examples, both real and hypothetical, to demonstrate its use. In addition, a number of different examples are shown how the data can be used by decision makers on issues such as spacecraft orientation and timing of EVAs.

Matney, Mark J.↗
NASA NTRS2012

Changes in the Intensity and Frequency of Atmospheric Blocking and Associated Heat Waves During Northern Summer Over Eurasia in the CMIP5 Model Simulations

The Russia heat wave and wild fires of the summer of 2010 was the most extreme weather event in the history of the country. Studies show that the root cause of the 2010 Russia heat wave/wild fires was an atmospheric blocking event which started to develop at the end of June and peaked around late July and early August. Atmospheric blocking in the summer of 2010 was anomalous in terms of the size, duration, and the location, which shifted to the east from the normal location. This and other similar continental scale severe summertime heat waves and blocking events in recent years have raised the question of whether such events are occurring more frequently and with higher intensity in a warmer climate induced by greenhouse gases. We studied the spatial and temporal distributions of the occurrence and intensity of atmospheric blocking and associated heat waves for northern summer over Eurasia based on CMIPS model simulations. To examine the global warming induced change of atmospheric blocking and heat waves, experiments for a high emissions scenario (RCP8.S) and a medium mitigation scenario (RCP4.S) are compared to the 20th century simulations (historical). Most models simulate the mean distributions of blockings reasonably well, including major blocking centers over Eurasia, northern Pacific, and northern Atlantic. However, the models tend to underestimate the number of blockings compared to MERRA and NCEPIDOE reanalysis, especially in western Siberia. Models also reproduced associated heat waves in terms of the shifting in the probability distribution function of near surface temperature. Seven out of eight models used in this study show that the frequency of atmospheric blocking over the Europe will likely decrease in a warmer climate, but slightly increase over the western Siberia. This spatial pattern resembles the blocking in the summer of 2010, indicating the possibility of more frequent occurrences of heat waves in western Siberia. In this talk, we will also discuss the potential effect of atmosphere-land feedback, particularly how the wetter spring affects the frequency and intensity of atmospheric blocking and heat wave during summer.

Kim, Kyu-Myong↗
DOE OSTI2015

Quantum theory as a description of robust experiments: Derivation of the Pauli equation

It is shown that the Pauli equation and the concept of spin naturally emerge from logical inference applied to experiments on a charged particle under the conditions that (i) space is homogeneous (ii) the observed events are logically independent, and (iii) the observed frequency distributions are robust with respect to small changes in the conditions under which the experiment is carried out. The derivation does not take recourse to concepts of quantum theory and is based on the same principles which have already been shown to lead to e.g. the Schrödinger equation and the probability distributions of pairs of particles in the singlet or triplet state. Application to Stern–Gerlach experiments with chargeless, magnetic particles, provides additional support for the thesis that quantum theory follows from logical inference applied to a well-defined class of experiments. - Highlights: • The Pauli equation is obtained through logical inference applied to robust experiments on a charged particle. • The concept of spin appears as an inference resulting from the treatment of two-valued data. • The same reasoning yields the quantum theoretical description of neutral magnetic particles. • Logical inference provides a framework to establish a bridge between objective knowledge gathered through experiments and their description in terms of concepts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Multi-variate model of T cell clonotype competition and homeostasis

Abstract Diversity of the naive T cell repertoire is maintained by competition for stimuli provided by self-peptides bound to major histocompatibility complexes (self-pMHCs). We extend an existing bi-variate competition model to a multi-variate model of the dynamics of multiple T cell clonotypes which share stimuli. In order to understand the late-time behaviour of the system, we analyse: (i) the dynamics until the extinction of the first clonotype, (ii) the time to the first extinction event, (iii) the probability of extinction of each clonotype, and (iv) the size of the surviving clonotypes when the first extinction event takes place. We also find the probability distribution of the number of cell divisions per clonotype before its extinction. The mean size of a new clonotype at quasi-steady state is an increasing function of the stimulus available to it, and a decreasing function of the fraction of stimuli it shares with other clonotypes. Thus, the probability of, and time to, extinction of a new clonotype entering the pool of T cell clonotypes is determined by the extent of competition for stimuli it experiences and by its initial number of cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Comparison of beam emission spectroscopy and gas puff imaging edge fluctuation measurements in National Spherical Torus Experiment

In this study, the close physical proximity of the Gas Puff Imaging (GPI) and Beam Emission Spectroscopy (BES) diagnostics on the National Spherical torus Experiment (NSTX) is leveraged to directly compare fluctuation measurements, and to study the local effects of the GPI neutral deuterium puff during H-mode plasmas without large Edge Localized Modes. The GPI and BES views on NSTX provide partially overlapping coverage of the edge and scrape-off layer (SOL) regions above the outboard midplane. The separation in the toroidal direction is 16°, and field lines passing through diagnostic views are separated by ∼20 cm in the direction perpendicular to the magnetic field. Strong cross-correlation is observed, and strong cross-coherence is seen for frequencies between 5 and 15 kHz. Also, probability distribution functions of fluctuations measured ∼3 cm inside the separatrix exhibit only minor deviations from a normal distribution for both diagnostics, and good agreement between correlation length estimates, decorrelation times, and structure velocities is found at the ±40% level. While the two instruments agree closely in many respects, some discrepancies are observed. Most notably, GPI normalized fluctuation levels exceed BES fluctuations by a factor of ∼9. BES mean intensity is found to be sensitive to the GPI neutral gas puff, and BES normalized fluctuation levels for frequencies between 1 and 10 kHz are observed to increase during the GPI puff.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

A Decision Theoretic Approach to Evaluate Radiation Detection Algorithms

There are a variety of sensor systems deployed at U.S. border crossings and ports of entry that scan for illicit nuclear material. In this work, we develop a framework for comparing the performance of detection algorithms that interpret the output of these scans and determine when secondary screening is needed. We optimize each algorithm to minimize its risk, or expected loss. We measure an algorithm’s risk by considering its performance over a sample, the probability distribution of threat sources, and the consequence of detection errors. While it is common to optimize algorithms by fixing one error rate and minimizing another, our framework allows one to simultaneously consider multiple types of detection errors. Our framework is flexible and easily adapted to many different assumptions regarding the probability of a vehicle containing illicit material, and the relative consequences of a false positive and false negative errors. Our methods can therefore inform decision makers of the algorithm family and parameter values which best reduce the threat from illicit nuclear material, given their understanding of the environment at any point in time. To illustrate the applicability of our methods, in this paper, we compare the risk from two families of detection algorithms and discuss the policy implications of our results.

Nobles, Mallory A.↗
NASA NTRS1975

Reliability-based optimum inspection and maintenance procedures

The development of reliability-based optimum inspection and maintenance schedules for engines needs an understanding of the fatigue behavior of the engines. Critical areas of the engine structure prone to fatigue damage are usually identified beforehand or after the fleet has been put into operation. In these areas, fatigue cracks initiate after several flight hours, and these cracks grow in length until failure takes place when these cracks attain the critical lengths. Crack initiation time and its growth rate are considered to be random variables. Usually, the inspection (fatigue) or test data from similar engines are used as prior distributions. The existing state-of-the-art is to ignore the different lengths of cracks obserbed at various inspections and to consider only the fact that a crack existed (or did not exist) at the time of inspection. In this paper, a procedure has been developed to obtain the probability of finding a crack of a given size at a certain time if the probability distributions for crack initiation and rates of growth are known. Application of the developed stochastic models to devise optimum procedures for inspection and maintenance are also discussed.

Nanagud, S.↗
NASA NTRS2006

Development of a Nonlinear Probability of Collision Tool for the Earth Observing System

The Earth Observing System (EOS) spacecraft Terra, Aqua, and Aura fly in constellation with several other spacecraft in 705-kilometer mean altitude sun-synchronous orbits. All three spacecraft are operated by the Earth Science Mission Operations (ESMO) Project at Goddard Space Flight Center (GSFC). In 2004, the ESMO project began assessing the probability of collision of the EOS spacecraft with other space objects. In addition to conjunctions with high relative velocities, the collision assessment method for the EOS spacecraft must address conjunctions with low relative velocities during potential collisions between constellation members. Probability of Collision algorithms that are based on assumptions of high relative velocities and linear relative trajectories are not suitable for these situations; therefore an algorithm for handling the nonlinear relative trajectories was developed. This paper describes this algorithm and presents results from its validation for operational use. The probability of collision is typically calculated by integrating a Gaussian probability distribution over the volume swept out by a sphere representing the size of the space objects involved in the conjunction. This sphere is defined as the Hard Body Radius. With the assumption of linear relative trajectories, this volume is a cylinder, which translates into simple limits of integration for the probability calculation. For the case of nonlinear relative trajectories, the volume becomes a complex geometry. However, with an appropriate choice of coordinate systems, the new algorithm breaks down the complex geometry into a series of simple cylinders that have simple limits of integration. This nonlinear algorithm will be discussed in detail in the paper. The nonlinear Probability of Collision algorithm was first verified by showing that, when used in high relative velocity cases, it yields similar answers to existing high relative velocity linear relative trajectory algorithms. The comparison with the existing high velocity/linear theory will also be used to determine at what relative velocity the analysis should use the new nonlinear theory in place of the existing linear theory. The nonlinear algorithm was also compared to a known exact solution for the probability of collision between two objects when the relative motion is strictly circular and the error covariance is spherically symmetric. Figure I shows preliminary results from this comparison by plotting the probabilities calculated from the new algorithm and those from the exact solution versus the Hard Body Radius to Covariance ratio. These results show about 5% error when the Hard Body Radius is equal to one half the spherical covariance magnitude. The algorithm was then combined with a high fidelity orbit state and error covariance propagator into a useful tool for analyzing low relative velocity nonlinear relative trajectories. The high fidelity propagator is capable of using atmospheric drag, central body gravitational, solar radiation, and third body forces to provide accurate prediction of the relative trajectories and covariance evolution. The covariance propagator also includes a process noise model to ensure realistic evolutions of the error covariance. This paper will describe the integration of the nonlinear probability algorithm and the propagators into a useful collision assessment tool. Finally, a hypothetical case study involving a low relative velocity conjunction between members of the Earth Observation System constellation will be presented.

McKinley, David P.↗
DOE OSTI2021

Energy sector portfolio analysis with uncertainty

Governments are dealing with the challenge of how to efficiently invest in research and development portfolios related to energy technologies. Research and development investment decisions in the energy space are especially difficult due to numerous risks and uncertainties, and due to the complexity of energy's interactions with the broad economy. Historically, much of the U.S. Department of Energy's in-depth research and development analyses focused on assessing the impact of a research and development activity in isolation from other available opportunities and did not substantially consider risk and uncertainty. Endeavoring to combine integrated energy-economy modeling with uncertainty analysis and technology-specific research and development activities, the U.S. Department of Energy commissioned the development of the Stochastic Energy Deployment System to support and improve public energy research and development decision-making. The Stochastic Energy Deployment System draws from expert-elicited probability distributions for research and development-driven improvements in technology cost and performance, and it uses Monte Carlo simulations to evaluate the likelihood of outcomes within a system dynamics energy-economy model. The framework estimates the uncertain benefits and costs of various research and development portfolios and provides insight into the probability of meeting national technology goals, while accounting for interactions with the larger economy and for interactions among research and development investments spanning many energy sectors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2017

A backward Monte-Carlo method for time-dependent runaway electron simulations

Kinetic descriptions of runaway electrons (REs) are usually based on Fokker-Planck models that determine the probability distribution function of REs in 2-dimensional momentum space. Despite the simplification involved, the Fokker-Planck equation can rarely be solved analytically and direct numerical approaches [e.g., continuum and particle-based Monte Carlo (MC)] can be time consuming, especially in the computation of asymptotic-type observables including the runaway probability, the slowing-down and runaway mean times, and the energy limit probability. In this study, we present a novel backward MC approach to these problems based on backward stochastic differential equations that describe the dynamics of the runaway probability by means of the Feynman-Kac theory. The key ingredient of the backward MC algorithm is to place all the particles in a runaway state and simulate them backward from the terminal time to the initial time. As such, our approach can provide much faster convergence than direct MC methods (by significantly reducing the number of particles required to achieve a prescribed accuracy) while at the same time maintaining the advantages of particle-based methods (compared to continuum approaches). The proposed algorithm is unconditionally stable and can be parallelized as easy as the direct MC method, and its extension to dimensions higher than two is straightforward, thus paving the way for conducting large-scale RE simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Projecting Future Energy Production from Operating Wind Farms in North America. Part I: Dynamical Downscaling

Abstract New simulations at 12-km grid spacing with the Weather and Research Forecasting (WRF) Model nested in the MPI Earth System Model (ESM) are used to quantify possible changes in wind power generation potential as a result of global warming. Annual capacity factors (CF; measures of electrical power production) computed by applying a power curve to hourly wind speeds at wind turbine hub height from this simulation are also used to illustrate the pitfalls in seeking to infer changes in wind power generation directly from low-spatial-resolution and time-averaged ESM output. WRF-derived CF are evaluated using observed daily CF from operating wind farms. The spatial correlation coefficient between modeled and observed mean CF is 0.65, and the root-mean-square error is 5.4 percentage points. Output from the MPI-WRF Model chain also captures some of the seasonal variability and the probability distribution of daily CF at operating wind farms. Projections of mean annual CF (CF A ) indicate no change to 2050 in the southern Great Plains and Northeast. Interannual variability of CF A increases in the Midwest, and CF A declines by up to 2 percentage points in the northern Great Plains. The probability of wind droughts (extended periods with anomalously low production) and wind bonus periods (high production) remains unchanged over most of the eastern United States. The probability of wind bonus periods exhibits some evidence of higher values over the Midwest in the 2040s, whereas the converse is true over the northern Great Plains. Significance Statement Wind energy is playing an increasingly important role in low-carbon-emission electricity generation. It is a “weather dependent” renewable energy source, and thus changes in the global atmosphere may cause changes in regional wind power production (PP) potential. We use PP data from operating wind farms to demonstrate that regional simulations exhibit skill in capturing actual power production. Projections to the middle of this century indicate that over most of North America east of the Rocky Mountains annual expected PP is largely unchanged, as is the probability of extended periods of anomalously high or low production. Any small declines in annual PP are of much smaller magnitude than changes due to technological innovation over the last two decades.

Meteorology & Atmospheric Sciences↗
DOE OSTI2018

Modeling Smart Microgrids for the Developing World with Probabilistic Supply and Demand Inputs

Developing world microgrids often balance insufficient supply with growing, unpredictable demand. Deterministic and probabilistic simulators exist to model these microgrids, and each focuses on different technical aspects. With the addition of smart meters into microgrids, monitoring and control is now available at high granularity, which enriches microgrid planning and operation. This research is designing a new simulator to model smart microgrids with discrete probability distributions as supply and demand inputs. In our model, smart meters allow real-time power clipping for demand side management, effectively smoothing the system load curve as needed. To compare clipping schemes for grid operation and generation mixes for planning, we aggregate inputs by convolution then compute expected energy sold and probability of avoiding power cuts. The simulator plots these values for different combinations of power clipping threshold and number of customers clipped.

13 HYDRO ENERGY↗
NASA NTRS1984

Development of a global model for atmospheric backscatter at CO2 wavelengths

The effect of aerosol microphysical processes on the backscatter from an aerosol plume undergoing long-range atmospheric transport was studied. A numerical model which examines the effects of coagulation and sedimentation on an aerosol size distribution is under development and the initial results for a single homogeneous layer were obtained. Use was made of the SAGE/SAM II data set to study the global variation of aerosol concentration and, hence, to predict the variation of Beta sub CO2. Computer programs were written to determine the mean, median, and the probability distribution of the measured aerosol extinction as a function of altitude, latitude and geographical conditions. The first data sets analyzed in this way are under study. Data was used to study aerosol behavior over the U.S.A. and the Pacific Ocean.

Kent, G. S.↗
DOE OSTI2014

Dynamical modeling of tidal streams

I present a new framework for modeling the dynamics of tidal streams. The framework consists of simple models for the initial action-angle distribution of tidal debris, which can be straightforwardly evolved forward in time. Taking advantage of the essentially one-dimensional nature of tidal streams, the transformation to position-velocity coordinates can be linearized and interpolated near a small number of points along the stream, thus allowing for efficient computations of a stream's properties in observable quantities. I illustrate how to calculate the stream's average location (its 'track') in different coordinate systems, how to quickly estimate the dispersion around its track, and how to draw mock stream data. As a generative model, this framework allows one to compute the full probability distribution function and marginalize over or condition it on certain phase-space dimensions as well as convolve it with observational uncertainties. This will be instrumental in proper data analysis of stream data. In addition to providing a computationally efficient practical tool for modeling the dynamics of tidal streams, the action-angle nature of the framework helps elucidate how the observed width of the stream relates to the velocity dispersion or mass of the progenitor, and how the progenitors of 'orphan' streams could be located. The practical usefulness of the proposed framework crucially depends on the ability to calculate action-angle variables for any orbit in any gravitational potential. A novel method for calculating actions, frequencies, and angles in any static potential using a single orbit integration is described in the Appendix.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Learning generative neural networks with physics knowledge

Deep generative neural networks have enabled modeling complex distributions, but incorporating physics knowledge into the neural networks is still challenging and is at the core of current physics-based machine learning research. To this end, we propose a physics generative neural network (PhysGNN), a new class of generative neural networks for learning unknown distributions in a physical system described by partial differential equations (PDE). PhysGNN couples PDE systems with generative neural networks. It is a fully differentiable model that allows back-propagation of gradients through both numerical PDE solvers and generative neural networks, and is trained by minimizing the discrete Wasserstein distance between generated and observed probability distributions of the PDE outputs using the stochastic gradient descent method. Moreover, PhysGNN does not require adversarial training like standard generative neural networks, which offers better stability than adversarial training. We show that PhysGNN can learn complex distributions in stochastic inverse problems, where conventional methods such as maximum likelihood estimation and momentum matching methods may be inapplicable when little knowledge is known about the form of unknown distributions or the physical model is too complex. Furthermore, our method allows physics-based generative neural network training for learning complex distributions in the context of differential equations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1984

An urnful of blinding functions

There is a fundamental connection between the mathematical theory of discrete probability distributions and the parametric curves and surfaces of computer-aided geometric design. It is no accident that the blending functions of Bezier curves and surfaces have an obvious probabilistic interpretation, nor is it a coincidence that the normalized uniform B-spline basis functions also model a simple stochastic process. The link between probability and geometry, and how to exploit simple probabilistic arguments to derive many of the classical geometric properties of the parametric curves and surfaces currently in vogue in computer-aided geometric design is discussed. This probabilistic approach is also used to introduce many new types of curves and surfaces, and it is demonstrated how probability theory can be used to simplify, unify, and generalize many well-known results.

Goldman, R. N.↗
DOE OSTI2018

The Absence of the Selfaveraging Property of the Entanglement Entropy of Disordered Free Fermions in One Dimension

We consider the macroscopic system of free lattice fermions in one dimension assuming that the one-body Hamiltonian of the system is the one dimensional discrete Schrödinger operator with independent identically distributed random potential. We show analytically and numerically that the variance of the entanglement entropy of the segment [−M,M] of the system is bounded away from zero as M→∞. This manifests the absence of the selfaveraging property of the entanglement entropy in our model, meaning that in the one-dimensional case the complete description of the entanglement entropy is provided by its whole probability distribution. This also may be contrasted the case of dimension two or more, where the variance of the entanglement entropy per unit surface area vanishes as M→∞ [8], thereby guaranteeing the representativity of its mean for large M in the multidimensional case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗