Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Probability distribution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,171 records · Page 65

Numerical Flow Visualization in Basic- and Hyper-Cluster Spheres

Packed spherical particle beds have wide applications throughout the process industry and are usually analyzed using an appropriate combination of laminar and turbulent flows with empirically derived coefficients of which the Ergun (1952) relation is probably the best known. The 3-D complexity of the void distribution within the beds has precluded detailed studies of sphere clusters. Numerical modeling and flow vector visualization for the basic tetra- and hexa-sphere clusters and two hyper-sphere clusters are presented at two Reynolds numbers, 400 and 1200. Cutting planes are used to enable visualization of the complex flows generated within the sphere clusters and are discussed herein. The boundary conditions and flow fields for the simple clusters are also compared to the hyper-clusters with larger variations found for hexa-clusters.

Hendricks, R. C.↗

Predictive Game Theory

Probability theory governs the outcome of a game; there is a distribution over mixed strat.'s, not a single "equilibrium". To predict a single mixed strategy must use our loss function (external to the game's players. Provides a quantification of any strategy's rationality. Prove rationality falls as cost of computation rises (for players who have not previously interacted). All extends to games with varying numbers of players.

Wolpert, David H.↗

Analytic Development of a Reference Profile for the First Entry in a Skip Atmospheric Entry

This note shows that a feasible reference drag profile for the first entry portion of a skip entry can be generated as a polynomial expression of the velocity. The coefficients of that polynomial are found through the resolution of a system composed of m + 1 equations, where m is the degree of the drag polynomial. It has been shown that a minimum of five equations (m = 4) are required to establish the range and the initial and final conditions on velocity and flight path angle. It has been shown that at least one constraint on the trajectory can be imposed through the addition of one extra equation in the system, which must be accompanied by the increase in the degree of the drag polynomial. In order to simplify the resolution of the system of equations, the drag was considered as being a probability density function of the velocity, with the velocity as a distribution function of the drag. Combining this notion with the introduction of empirically derived constants, it has been shown that the system of equations required to generate the drag profile can be successfully reduced to a system of linear algebraic equations. For completeness, the resulting drag profiles have been flown using the feedback linearization method of differential geometric control as a guidance law with the error dynamics of a second order homogeneous equation in the form of a damped oscillator. Satisfactory results were achieved when the gains in the error dynamics were changed at a certain point along the trajectory that is dependent on the velocity and the curvature of the drag as a function of the velocity. Future work should study the capacity to update the drag profile in flight when dispersions are introduced. Also, future studies should attempt to link the first entry, as presented and controlled in this note, with a more standard control concept for the second entry, such as the Apollo entry guidance, to try to assess the overall skip entry performance. A guidance law that includes an integral feedback term, as is the case in the actual Space Shuttle entry guidance and as is proposed in Ref 29, could be tried in future studies to assess whether its use results in an improvement of the tracking performance, and to evaluate the design needs when determining the control gains.

Garcia-Llama, Eduardo↗

Second Cancers After Fractionated Radiotherapy: Stochastic Population Dynamics Effects

When ionizing radiation is used in cancer therapy it can induce second cancers in nearby organs. Mainly due to longer patient survival times, these second cancers have become of increasing concern. Estimating the risk of solid second cancers involves modeling: because of long latency times, available data is usually for older, obsolescent treatment regimens. Moreover, modeling second cancers gives unique insights into human carcinogenesis, since the therapy involves administering well characterized doses of a well studied carcinogen, followed by long-term monitoring. In addition to putative radiation initiation that produces pre-malignant cells, inactivation (i.e. cell killing), and subsequent cell repopulation by proliferation can be important at the doses relevant to second cancer situations. A recent initiation/inactivation/proliferation (IIP) model characterized quantitatively the observed occurrence of second breast and lung cancers, using a deterministic cell population dynamics approach. To analyze ifradiation-initiated pre-malignant clones become extinct before full repopulation can occur, we here give a stochastic version of this I I model. Combining Monte Carlo simulations with standard solutions for time-inhomogeneous birth-death equations, we show that repeated cycles of inactivation and repopulation, as occur during fractionated radiation therapy, can lead to distributions of pre-malignant cells per patient with variance >> mean, even when pre-malignant clones are Poisson-distributed. Thus fewer patients would be affected, but with a higher probability, than a deterministic model, tracking average pre-malignant cell numbers, would predict. Our results are applied to data on breast cancers after radiotherapy for Hodgkin disease. The stochastic IIP analysis, unlike the deterministic one, indicates: a) initiated, pre-malignant cells can have a growth advantage during repopulation, not just during the longer tumor latency period that follows; b) weekend treatment gaps during radiotherapy, apart from decreasing the probability of eradicating the primary cancer, substantially increase the risk of later second cancers.

Sachs, Rainer K.↗

Tidal Disruptions of Main Sequence Stars - I. Observable Quantities and their Dependence on Stellar and Black Hole Mass

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses(MM0.1510)and six black hole masses( MM105 105BH7). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of Må for Må/Me3. Independent of Må, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of Må; at MBH/Me;5×106 this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of0.35−2, depending on Må and MBH, implying a corresponding change(∝(ΔE)−3/2)in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star’s angular momentum and the fraction of Må remaining. Within the “empty loss-cone” regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions.

Taeho Ryu↗

Tidal Disruptions of Main Sequence Stars - II. Simulation Methodology and Stellar Mass Dependence of the Character of Full Tidal Disruptions

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses (0.15 < M_/M<10)and six black hole masses( ́10<MM105 105BH7). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of Må for Må/Me3. Independent of Må, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of ∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of Må; at MBH/Me5×106this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of 0.35−2, depending on Må and MBH, implying a corresponding change(∝(ΔE)−3/2)in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star’s angular momentum and the fraction of Må remaining. Within the “empty loss-cone” regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions

Taeho Ryu↗

The use of prior probabilities in maximum likelihood classification of remotely sensed data

Possibilities for the improvement of classification accuracies by the use of prior information about the expected distribution of classes in the maximum likelihood classification of remote sensing data are examined. The modification of the maximum likelihood decision rule to take into account one or several sets of probabilities for the occurrence of classes which probabilities are based on independent knowledge of the area surveyed is demonstrated. It is then shown that the use of prior probabilities is sufficiently versatile so as to allow the prior weighting of output classes based on their anticipated sizes as well as the merging of continuously varying measurements with discrete collateral information data sets and the construction of time-sequential classification systems in which an earlier classification modifies the outcome of a latter one.

Strahler, A. H.↗

The statistics of gravitational lenses - The distributions of image angular separations and lens redshifts

An analytical calculation is presented which determines the probability that the light from a distant quasar can be lensed by objects in a different redshift range. Three types of information are used to make the calculation: (1) the physical properties of the intervening scattering lenses; (2) the intrinsic properties of the quasars and (3) the selection procedures by which the quasar samples are acquired. A series of differential equations is given in which the probability is derived, and a number of recommendations are offered with respect to VLA and Space Telescope observations of quasars lensed in this way.

Turner, E. L.↗

Probabilistic Prediction of Lifetimes of Ceramic Parts

ANSYS/CARES/PDS is a software system that combines the ANSYS Probabilistic Design System (PDS) software with a modified version of the Ceramics Analysis and Reliability Evaluation of Structures Life (CARES/Life) Version 6.0 software. [A prior version of CARES/Life was reported in Program for Evaluation of Reliability of Ceramic Parts (LEW-16018), NASA Tech Briefs, Vol. 20, No. 3 (March 1996), page 28.] CARES/Life models effects of stochastic strength, slow crack growth, and stress distribution on the overall reliability of a ceramic component. The essence of the enhancement in CARES/Life 6.0 is the capability to predict the probability of failure using results from transient finite-element analysis. ANSYS PDS models the effects of uncertainty in material properties, dimensions, and loading on the stress distribution and deformation. ANSYS/CARES/PDS accounts for the effects of probabilistic strength, probabilistic loads, probabilistic material properties, and probabilistic tolerances on the lifetime and reliability of the component. Even failure probability becomes a stochastic quantity that can be tracked as a response variable. ANSYS/CARES/PDS enables tracking of all stochastic quantities in the design space, thereby enabling more precise probabilistic prediction of lifetimes of ceramic components.

Nemeth, Noel N.↗

Public Risk Assessment Program

The Public Entry Risk Assessment (PERA) program addresses risk to the public from shuttle or other spacecraft re-entry trajectories. Managing public risk to acceptable levels is a major component of safe spacecraft operation. PERA is given scenario inputs of vehicle trajectory, probability of failure along that trajectory, the resulting debris characteristics, and field size and distribution, and returns risk metrics that quantify the individual and collective risk posed by that scenario. Due to the large volume of data required to perform such a risk analysis, PERA was designed to streamline the analysis process by using innovative mathematical analysis of the risk assessment equations. Real-time analysis in the event of a shuttle contingency operation, such as damage to the Orbiter, is possible because PERA allows for a change to the probability of failure models, therefore providing a much quicker estimation of public risk. PERA also provides the ability to generate movie files showing how the entry risk changes as the entry develops. PERA was designed to streamline the computation of the enormous amounts of data needed for this type of risk assessment by using an average distribution of debris on the ground, rather than pinpointing the impact point of every piece of debris. This has reduced the amount of computational time significantly without reducing the accuracy of the results. PERA was written in MATLAB; a compiled version can run from a DOS or UNIX prompt.

Mendeck, Gavin↗

Radar Cross Section of Orbital Debris Objects

This discussion is concerned with the radar-data analysis and usage involved in the building of model orbital debris (OD) populations in the near-Earth environment, focusing on radar cross section (RCS). While varying with radar wavelength, physical dimension, material composition, overall shape and structure, the RCS of an irregular object is also strongly dependent on its spatial orientation. The historical records of observed RCSs for cataloged OD objects in the Space Surveillance Network are usually distributed over an RCS range, forming respective characteristic patterns. The National Aeronautics and Space Administration (NASA) Size Estimation Model provides an empirical probability-density function of RCS as a function of “effective diameter” (or “characteristic length”), which makes it feasible to predict possible RCS distributions for a given model OD population and to link data with model from a statistical perspective. The discussion also includes application of the widely used method of moments (MoM) and the Generalized Multi-particle Mie-solution (GMM) in the prediction of the RCS of arbitrarily shaped objects. Theoretical calculation results for an aluminum cube are compared with corresponding experimental measurements.

Xu, Yu-Lin↗

Risk Assessment of Bone Fracture During Space Exploration Missions to the Moon and Mars

The possibility of a traumatic bone fracture in space is a concern due to the observed decrease in astronaut bone mineral density (BMD) during spaceflight and because of the physical demands of the mission. The Bone Fracture Risk Module (BFxRM) was developed to quantify the probability of fracture at the femoral neck and lumbar spine during space exploration missions. The BFxRM is scenario-based, providing predictions for specific activities or events during a particular space mission. The key elements of the BFxRM are the mission parameters, the biomechanical loading models, the bone loss and fracture models and the incidence rate of the activity or event. Uncertainties in the model parameters arise due to variations within the population and unknowns associated with the effects of the space environment. Consequently, parameter distributions were used in Monte Carlo simulations to obtain an estimate of fracture probability under real mission scenarios. The model predicts an increase in the probability of fracture as the mission length increases and fracture is more likely in the higher gravitational field of Mars than on the moon. The resulting probability predictions and sensitivity analyses of the BFxRM can be used as an engineering tool for mission operation and resource planning in order to mitigate the risk of bone fracture in space.

Lewandowski, Beth E.↗

Risk Assessment of Bone Fracture During Space Exploration Missions to the Moon and Mars

The possibility of a traumatic bone fracture in space is a concern due to the observed decrease in astronaut bone mineral density (BMD) during spaceflight and because of the physical demands of the mission. The Bone Fracture Risk Module (BFxRM) was developed to quantify the probability of fracture at the femoral neck and lumbar spine during space exploration missions. The BFxRM is scenario-based, providing predictions for specific activities or events during a particular space mission. The key elements of the BFxRM are the mission parameters, the biomechanical loading models, the bone loss and fracture models and the incidence rate of the activity or event. Uncertainties in the model parameters arise due to variations within the population and unknowns associated with the effects of the space environment. Consequently, parameter distributions were used in Monte Carlo simulations to obtain an estimate of fracture probability under real mission scenarios. The model predicts an increase in the probability of fracture as the mission length increases and fracture is more likely in the higher gravitational field of Mars than on the moon. The resulting probability predictions and sensitivity analyses of the BFxRM can be used as an engineering tool for mission operation and resource planning in order to mitigate the risk of bone fracture in space.

Lewandowski, Beth E.↗

Photoionization of the diffuse interstellar medium and galactic halo by OB associtations

Assuming smoothly varying H I distributions in te Galactic disk, we have calculated the geometry of diffuse II regions due to OB associations in the Galactic plane. Near the solar circle, OB associations with a Lyman continuum (Lyc) photon luminosity Psi(sub Lyc) = 3.3 x 10(exp 7) cm(exp -2) s(exp -1), produce H II regions that are density bounded in the vertical direction (H II chimneys) allowing Lyc to escape the gaseous disk and penetrate into the Galactic halo. We provide analytic formulae for the Lyc escape fraction as functions of S(sub 0) O-star catalog of Garmany and a new Lyc stellar stellar Lyc stellar flux calibration, we find a production rate of Lyc photons by OB associations within 2.5 kpc of Psi(sub Lyc) = 3.3 x 10(exp 7) cm(exp -2) s(exp -1). Integrating the fraction of Lyc photons that escape the disk over our adopted luminosity function of OB associations, we estimate that approximately 7% of the ionizing photons, or Phi(sub Lyc) = 2.3 x 10(exp 6) cm(exp -2) s(exp -1), escape each side of the H I disk layer and penetrate the diffuse ionized medium ('Reynolds layer'). This flux is sufficient to explain the potoionization of this, although we have not constructed a model for the observed H-alpha emission and pulsar dispersion measures that is fully consistent with the absorption rate of Lyc in the H II layer. Since our quiescent model does not account for the effects of dynamic chimneys and superbubbles, which should enhance Lyc escape, we conclude the O stars are the probable source of ionizing radiation for the Reynolds layer. For a random distribution of OB associations throughout the disk, the Lyc flux is nearly uniform for heights Z is greater than approximately 0.8 kpc above the midplane.

Dove, James B.↗

Boost-phase discrimination research

The final report describes the combined work of the Computational Chemistry and Aerothermodynamics branches within the Thermosciences Division at NASA Ames Research Center directed at understanding the signatures of shock-heated air. Considerable progress was made in determining accurate transition probabilities for the important band systems of NO that account for much of the emission in the ultraviolet region. Research carried out under this project showed that in order to reproduce the observed radiation from the bow shock region of missiles in their boost phase it is necessary to include the Burnett terms in the constituent equation, account for the non-Boltzmann energy distribution, correctly model the NO formation and rotational excitation process, and use accurate transition probabilities for the NO band systems. This work resulted in significant improvements in the computer code NEQAIR that models both the radiation and fluid dynamics in the shock region.

Langhoff, Stephen R.↗

Plasma turbulence in the downstream ionosheath of Venus

Observations made by the Pioneer Venus Orbiter plasma analyzer and the plasma wave instrument in the Venus ionosheath are compared. Large increases in plasma wave turbulence levels appear to be connected with changing plasma distributions and interpenetrating plasma beams. Some of these plasma waves are identified as Doppler - shifted ion acoustic waves due to beam/beam interactions, but it is noted that different forms of instabilities are probably also operative. The changes in the temperature, intensity and energy of the peak in the PVO plasma distributions are similar to those observed by Venera 10 closer to the planet and appear to be evidence for rarefaction and compression in the downstream ionosheath. Some of the changes in the PVO plasma distributions may be related to the presence of a second ion population or the acceleration of protons.

Intriligator, D. S.↗

Redundancy: How Many Unreliable Spares are Needed for High Reliability and Confidence on a Time Limited Mission?

This paper investigates the number of redundant units needed to achieve high reliability with high confidence. The approach applies to the case where the unit failure rate is too high for a single unit to provide the required reliability over the mission duration. To achieve high reliability, the design then uses N redundant units, one operating unit and N – 1 spares. If the unit failure rate is f, the mission length is L, and f * L is small (not the case assumed here), the unit failure probability over the mission duration is F1 = f * L << 1. In this case, the probability that all N units will fail is FN = F1N, and the needed N = LN(FN)/LN(F1). For the case of large f * L assumed here, F1 = f * L > 1, and F1 is the expected number of failures during the mission. The needed redundancy, N, to achieve the specified N unit reliability, FN, can be computed using the cumulative Poisson distribution with mean equal to F1. The number of spares, N - 1, is increased until the probability - that the total number of failures will be less than N -1 - achieves the required reliability. The confidence that this reliability can be achieved can be computed using the cumulative Poisson distribution or the chi-square distribution. Since the measured unit failure rate, f, has some uncertainty, the confidence that the rate is not lower than the actual failure rate and the required reliability is not overestimated is about 50%. Adding more redundant units increases the confidence that the required reliability, FN, will be achieved. For a fixed number of redundant units, the expected reliability and confidence can be traded off, since lower reliability goals have higher confidence in being achieved. Both the required reliability and confidence can be specified initially and the needed number of redundant units computed using the measured failure rate. The unit failure rate is determined by initial reliability growth testing to remove design errors and to better estimate the final constant failure rate. Reducing the failure rate and reducing its variance both reduce the number of redundant units needed for the required reliability and confidence. Since the total cost is the sum of the costs of the units and of the testing, there is an optimum test time that produces minimum cost.

Harry W. Jones↗

Monte Carlo analysis of lobular gas-surface scattering in tubes applied to thermal transpiration

A model of rarefied gas flow in tubes was developed which combines a lobular distribution with diffuse reflection at the wall. The model with Monte Carlo techniques was used to explain previously observed deviations in the free molecular thermal transpiration ratio which suggest molecules can have a greater tube transmission probability in a hot-to-cold direction than in a cold-to-hot direction. The model yields correct magnitudes of transmission probability ratios for helium in Pyrex tubing (1.09 to 1.14), and some effects of wall-temperature distribution, tube surface roughness, tube dimensions, gas temperature, and gas molecular mass.

Smith, J. D.↗