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,621 records · Page 90

Probabilistic Simulation of Combined Thermo-Mechanical Cyclic Fatigue in Composites

A methodology to compute probabilistically-combined thermo-mechanical fatigue life of polymer matrix laminated composites has been developed and is demonstrated. Matrix degradation effects caused by long-term environmental exposure and mechanical/thermal cyclic loads are accounted for in the simulation process. A unified time-temperature-stress-dependent multifactor-interaction relationship developed at NASA Glenn Research Center has been used to model the degradation/aging of material properties due to cyclic loads. The fast probability-integration method is used to compute probabilistic distribution of response. Sensitivities of fatigue life reliability to uncertainties in the primitive random variables (e.g., constituent properties, fiber volume ratio, void volume ratio, ply thickness, etc.) computed and their significance in the reliability-based design for maximum life is discussed. The effect of variation in the thermal cyclic loads on the fatigue reliability for a (0/+/-45/90)s graphite/epoxy laminate with a ply thickness of 0.127 mm, with respect to impending failure modes has been studied. The results show that, at low mechanical-cyclic loads and low thermal-cyclic amplitudes, fatigue life for 0.999 reliability is most sensitive to matrix compressive strength, matrix modulus, thermal expansion coefficient, and ply thickness. Whereas at high mechanical-cyclic loads and high thermal-cyclic amplitudes, fatigue life at 0.999 reliability is more sensitive to the shear strength of matrix, longitudinal fiber modulus, matrix modulus, and ply thickness.

Chamis, Christos C.↗

Probabilistic Simulation for Combined Cycle Fatigue in Composites

A methodology to compute probabilistic fatigue life of polymer matrix laminated composites has been developed and demonstrated. Matrix degradation effects caused by long term environmental exposure and mechanical/thermal cyclic loads are accounted for in the simulation process. A unified time-temperature-stress dependent multifactor interaction relationship developed at NASA Glenn Research Center has been used to model the degradation/aging of material properties due to cyclic loads. The fast probability integration method is used to compute probabilistic distribution of response. Sensitivities of fatigue life reliability to uncertainties in the primitive random variables (e.g., constituent properties, fiber volume ratio, void volume ratio, ply thickness, etc.) computed and their significance in the reliability-based design for maximum life is discussed. The effect of variation in the thermal cyclic loads on the fatigue reliability for a (0/+/- 45/90)s graphite/epoxy laminate with a ply thickness of 0.127 mm, with respect to impending failure modes has been studied. The results show that, at low mechanical cyclic loads and low thermal cyclic amplitudes, fatigue life for 0.999 reliability is most sensitive to matrix compressive strength, matrix modulus, thermal expansion coefficient, and ply thickness. Whereas at high mechanical cyclic loads and high thermal cyclic amplitudes, fatigue life at 0.999 reliability is more sensitive to the shear strength of matrix, longitudinal fiber modulus, matrix modulus, and ply thickness.

Chamis, Christos C.↗

Probabilistic Simulation of Combined Thermo-Mechanical Cyclic Fatigue in Composites

A methodology to compute probabilistically-combined thermo-mechanical fatigue life of polymer matrix laminated composites has been developed and is demonstrated. Matrix degradation effects caused by long-term environmental exposure and mechanical/thermal cyclic loads are accounted for in the simulation process. A unified time-temperature-stress-dependent multifactor-interaction relationship developed at NASA Glenn Research Center has been used to model the degradation/aging of material properties due to cyclic loads. The fast probability-integration method is used to compute probabilistic distribution of response. Sensitivities of fatigue life reliability to uncertainties in the primitive random variables (e.g., constituent properties, fiber volume ratio, void volume ratio, ply thickness, etc.) computed and their significance in the reliability-based design for maximum life is discussed. The effect of variation in the thermal cyclic loads on the fatigue reliability for a (0/+/-45/90)s graphite/epoxy laminate with a ply thickness of 0.127 mm, with respect to impending failure modes has been studied. The results show that, at low mechanical-cyclic loads and low thermal-cyclic amplitudes, fatigue life for 0.999 reliability is most sensitive to matrix compressive strength, matrix modulus, thermal expansion coefficient, and ply thickness. Whereas at high mechanical-cyclic loads and high thermal-cyclic amplitudes, fatigue life at 0.999 reliability is more sensitive to the shear strength of matrix, longitudinal fiber modulus, matrix modulus, and ply thickness.

Chamis, Christos C.↗

Smooth Pursuit of Flicker-Defined Motion

We examined the pursuit response to stimuli defined by space-variant flicker of a dense random dot carrier pattern. On each frame, every element of the pattern could change polarity, with a probability given by a two-dimensional Gaussian distribution. A normal distribution produces a circular region of twinkle, while inverting the distribution results in a spot of static texture in a twinkling surround. In this latter case, the carrier texture could be stationary, or could move with the twinkle modulator, thereby producing first-order motion in the region of the spot. While the twinkle-defined spot produces a strong sensation of motion, the complementary stimulus defined by the absence of twinkle does not, when viewed peripherally, it appears to move in steps even when the generating distribution moves smoothly. We examined pursuit responses to these stimuli using two techniques: 1) the eye movement correlogram, obtained by cross-correlating eye velocity with the velocity of a randomly-moving stimulus; and 2) delayed visual feedback, where transient stabilization of a target can produce spontaneous oscillations of the eye, with a period empirically observed to vary linearly with the applied delay. Both techniques provide an estimate of the internal processing time, which can be as short as 100 milliseconds for a first-order target. Assessed by the correlogram method, the response to flicker-defined motion is delayed by more than 100 milliseconds, and significantly weaker (especially in the vertical dimension). When initially presented in the delayed feedback condition, purely saccadic oscillation is observed. One subject eventually developed smooth oscillations (albeit with significant saccadic intrusions), showing a period-versus-delay slope similar to that observed for first-order targets. This result is somewhat surprising, given that we interpret the slope of the period-versus-delay-function as reflecting the balance between position- and velocity-sensitive inputs to pursuit.

eye movements↗

System-Level Prognostics for the National Airspace

In the National Airspace System (NAS), safety is assured through a set of rules, regulations, and procedures to respond to unsafe events. However, safety stands to benefit immensely from the introduction of tools and methodologies from Prognostics and Health Management (PHM). PHM will enable the NAS to stochastically predict unsafe states within the NAS, enabling a proactive preventative response strategy, as opposed to a reactive mitigative one. However, current PHM methods do not directly apply to the NAS for several reasons: they typically apply only at the component level, are implemented in a centralized manner, and are focused only on predicting remaining useful life. In this paper, we extend the model-based prognostics approach to PHM in order to provide a framework that can be applied to the NAS. We offer a system-level approach that supports a distributed implementation, and provide algorithms to predict the probability of an unsafe state, either at a specific time or within a time interval, and to predict the time of an unsafe state. Experimental results in simulation demonstrate the new approach.

prognostics↗

Smooth Pursuit of Flicker-Defined Motion

We examined the pursuit response to stimuli defined by space-variant flicker of a dense random dot carrier pattern. On each frame, every element of the pattern could change polarity, with a probability given by a two-dimensional Gaussian distribution. A normal distribution produces a circular region of twinkle, while inverting the distribution results in a spot of static texture in a twinkling surround. In this latter case, the carrier texture could be stationary, or could move with the twinkle modulator,thereby producing first-order motion in the region of the spot. While the twinkle-defined spot produces a strong sensation of motion, the complementary stimulus defined by the absence of twinkle does not;when viewed peripherally, it appears to move in steps even when the generating distribution moves smoothly. We examined pursuit responses to these stimuli using two techniques: 1) the eye movement correlogram, obtained by cross-correlating eye velocity with the velocity of a randomly-moving stimulus; and 2) delayed visual feedback, where transient stabilization of a target can produce spontaneous oscillations of the eye, with a period empirically observed to vary linearly with the applied delay. Both techniques provide an estimate of the internal processing time, which can be as short as100 milliseconds for a first-order target. Assessed by the correlogram method, the response to flicker defined motion is delayed by more than 100 milliseconds, and significantly weaker (especially in the vertical dimension). When initially presented in the delayed feedback condition, purely saccadic oscillation is observed. One subject eventually developed smooth oscillations (albeit with significant saccadic intrusions), showing a period-versus-delay slope similar to that observed for first-order targets. This result is somewhat surprising, given that we interpret the slope of the period-versus-delayfunction as reflecting the balance between position- and velocity-sensitive inputs to pursuit.

smooth pursuit↗

Improving the chi-squared approximation for bivariate normal tolerance regions

Let X be a two-dimensional random variable distributed according to N2(mu,Sigma) and let bar-X and S be the respective sample mean and covariance matrix calculated from N observations of X. Given a containment probability beta and a level of confidence gamma, we seek a number c, depending only on N, beta, and gamma such that the ellipsoid R = (x: (x - bar-X)'S(exp -1) (x - bar-X) less than or = c) is a tolerance region of content beta and level gamma; i.e., R has probability gamma of containing at least 100 beta percent of the distribution of X. Various approximations for c exist in the literature, but one of the simplest to compute -- a multiple of the ratio of certain chi-squared percentage points -- is badly biased for small N. For the bivariate normal case, most of the bias can be removed by simple adjustment using a factor A which depends on beta and gamma. This paper provides values of A for various beta and gamma so that the simple approximation for c can be made viable for any reasonable sample size. The methodology provides an illustrative example of how a combination of Monte-Carlo simulation and simple regression modelling can be used to improve an existing approximation.

Feiveson, Alan H.↗

Comparison of Three Multicylinder Icing Meters and Critique of Multicylinder Method

Three multicylinder cloud meters, fundamentally similar but differing in important details, were compared in use at the Mount Washington Observatory. Determinations of liquid water content were found to agree within the limits of the probable error, but the two instruments designed by the National Advisory Committee for Aeronautics indicated larger drop sizes than did the Observatory's instrument, apparently because of spurious ice catch on the rather rough surface of the larger cylinders. Comparisons of drop-size distribution were largely indeterminate., In a critique of the method, the probable error of determination of liquid water content was found to be +/-8 percent; of drop size, +/-6 percent; and of drop-size distribution, about +/-0.7 unit of the modulus of distribution. Of the systematic errors, run-off of unfrozen water is most important, blow-off and erosion seldom being hampering. Revision of collection-efficiency computations for cylinders in clouds with distributed drop sizes was found necessary and also revision of one of the correction-factor graphs heretofore used. The assumption of constant ice density in deriving cylinder size was found to be permissible for cylinders 1 inch or more in diameter.

Howell, Wallace E.↗

Guest Editorial: Planning and operation of resilient distribution system for integrated multi-energy

Resilience is the ability of power systems to prepare for and adapt to low-probability, high-impact incidents and withstand and recover rapidly from disruptions. With the ageing of electricity distribution infrastructure and increasing threats of weather-related incidents and natural disasters, the need to effectively enhance the resilience of the electricity distribution system has become urgent and has attracted worldwide attention. Although there are an increasing number of publications related to enhancing resilience strategies, resilience is an emerging concept in power systems. Existing practices are mostly focused on deploying distributed energy resources (DERs) and microgrids, hardening the existing infrastructures and building redundant capacities. However, from a broader perspective, resilience enhancement of power distribution system is a systematic engineering, involving long-term system planning and upgrading (e.g., deployment of smart grid technologies and intelligent switches), short-term proactive scheduling, real-time robust and resilient control of DERs, and post-event restoration and recovery strategies. Based on this point, this Special Issue in IET Energy System Integration focuses on soliciting the most recent and original technologies, scheduling and control strategies for improving the resilience of power distribution system. Eight papers are presented in this Special Issue, covering various aspects related to resilience enhancement of power distribution system, including fault-tolerant frequency measurement, robust scheduling of integrated electricity and district heating systems, robust control of DERs, efficient methods for safety verification as well as novel graph theory-based approach to restore the distribution systems after multiple simultaneous faults. A brief introduction of these 8 papers is provided below.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Cosmic gamma-ray lines

Recent observations of gamma-ray line emissions from gamma-ray bursts, the ISM, the Galactic center, and solar flares are reviewed, and the implications of these observations for high-energy processes in these sources are discussed. Line observations suggest that magnetized neutron stars are probably the best candidate objects for burst sources. Observations of the 1.809-MeV line from Al-26 decay provide evidence for ongoing nucleosynthesis in the Galaxy and information on the spatial distribution of nucleosynthetic sites. The compact 0.511-MeV line source is probably a black hole at or close to the Galactic center. Solar-flare studies have provided new information on the confinement and escape of charged particles at the sun and on multiple acceleration phases in solar flares.

Ramaty, R.↗

Combined Structural and Compositional Evolution of Planetary Rings Due to Micrometeoroid Impacts and Ballistic Transport

We introduce improved numerical techniques for simulating the structural and compositional evolution of planetary rings due to micrometeoroid bombardment and subsequent ballistic transport of impact ejecta. Our current, robust code is capable of modeling structural changes and pollution transport simultaneously over long times on both local and global scales. In this paper, we describe the methodology based on the original structural code of Durisen et al. (1989, Icarus 80, 136-166) and on the pollution transport code of Cuzzi and Estrada (1998, Icarus 132, 1-35). We provide demonstrative simulations to compare with, and extend upon previous work, as well as examples of how ballistic transport can maintain the observed structure in Saturn's rings using available Cassini occultation optical depth data. In particular, we explicitly verify the claim that the inner B (and presumably A) ring edge can be maintained over long periods of time due to an ejecta distribution that is heavily biased in the prograde direction through a balance between the sharpening effects of ballistic transport and the broadening effects of viscosity. We also see that a "ramp"-like feature forms over time just inside that edge. However, it does not remain linear for the duration of the runs presented here unless a less steep ejecta velocity distribution is adopted. We also model the C ring plateaus and find that their outer edges can be maintained at their observed sharpness for long periods due to ballistic transport. We hypothesize that the addition of a significant component of a retrograde-biased ejecta distribution may help explain the linearity of the ramp and is probably essential for maintaining the sharpness of C ring plateau inner edges. This component would arise for the subset of micrometeoroid impacts which are destructive rather than merely cratering. Such a distribution will be introduced in future work.

Cassini↗

Theoretical foundations for quantitative paleogenetics. III - The molecular divergence of nucleic acids and proteins for the case of genetic events of unequal probability

Theoretical equations are derived for molecular divergence with respect to gene and protein structure in the presence of genetic events with unequal probabilities: amino acid and base compositions, the frequencies of nucleotide replacements, the usage of degenerate codons, the distribution of fixed base replacements within codons and the distribution of fixed base replacements among codons. Results are presented in the form of tables relating the probabilities of given numbers of codon base changes with respect to the original codon for the alpha hemoglobin, beta hemoglobin, myoglobin, cytochrome c and parvalbumin group gene families. Application of the calculations to the rabbit alpha and beta hemoglobin mRNAs and proteins indicates that the genes are separated by about 425 fixed based replacements distributed over 114 codon sites, which is a factor of two greater than previous estimates. The theoretical results also suggest that many more base replacements are required to effect a given gene or protein structural change than previously believed.

Holmquist, R.↗

Parallel discrete event simulation: A shared memory approach

With traditional event list techniques, evaluating a detailed discrete event simulation model can often require hours or even days of computation time. Parallel simulation mimics the interacting servers and queues of a real system by assigning each simulated entity to a processor. By eliminating the event list and maintaining only sufficient synchronization to insure causality, parallel simulation can potentially provide speedups that are linear in the number of processors. A set of shared memory experiments is presented using the Chandy-Misra distributed simulation algorithm to simulate networks of queues. Parameters include queueing network topology and routing probabilities, number of processors, and assignment of network nodes to processors. These experiments show that Chandy-Misra distributed simulation is a questionable alternative to sequential simulation of most queueing network models.

Reed, Daniel A.↗

A Discussion of Time Management Concepts and Time Constraint Equations for Multi-Rate Federation Executions

The High Level Architecture (HLA) is a simulation interoperability standard developed by the Simulation Interoperability Standards Organization (SISO) and published as the international standard IEEE 1516-2010 by the Institute for Electrical and Electronics Engineers (IEEE). HLA is a widely used standard for the development and execution of collaborative distributed simulations. HLA provides a number of Management Services to simulation developers: Federation, Declaration, Object, Ownership, Data Distribution, and Time. Of those services, Time Management Services is probably one of the least understood and least used. However, Time Management Services are critical to technical simulations like those created for space systems using the Space Reference Federation Object Model (SpaceFOM). Time Management can be used to insure data coherence and execution repeatability in distributed simulations. When combined with real time execution policies, Time Management is being used to support real time execution of mixed software and hardware in the loop integration, verification, and validation simulations for active space systems development. This paper starts by providing an overview of the HLA Time Management Services. This provides the background to discuss the challenges associated with Time Management and its use, starting with simple common rate frame scheduled simulations, then simple multi-rate simulations, and ending with complex mixed rate simulations. The authors then formulate the significant time constraint relationships between identified frame scheduling parameters. The intent of the paper is to provide a concise discussion of how to use Time Management in both simple cases and in more complex mixed frame rate federation executions.

Simulation Interoperability↗

A Discussion of Time Management Concepts and Time Constraint Equations for Multi-Rate Federation Executions

The High Level Architecture (HLA) is a simulation interoperability standard developed by the Simulation Interoperability Standards Organization (SISO) and published as the international standard IEEE 1516-2010 by the Institute for Electrical and Electronics Engineers (IEEE). HLA is a widely used standard for the development and execution of collaborative distributed simulations. HLA provides a number of Management Services to simulation developers: Federation, Declaration, Object, Ownership, Data Distribution, and Time. Of those services, Time Management Services is probably one of the least understood and least used. However, Time Management Services are critical to technical simulations like those created for space systems using the Space Reference Federation Object Model (SpaceFOM). Time Management can be used to insure data coherence and execution repeatability in distributed simulations. When combined with real time execution policies, Time Management is being used to support real time execution of mixed software and hardware in the loop integration, verification, and validation simulations for active space systems development. This paper starts by providing an overview of the HLA Time Management Services. This provides the background to discuss the challenges associated with Time Management and its use, starting with simple common rate frame scheduled simulations, then simple multi-rate simulations, and ending with complex mixed rate simulations. The authors then formulate the significant time constraint relationships between identified frame scheduling parameters. The intent of the paper is to provide a concise discussion of how to use Time Management in both simple cases and in more complex mixed frame rate federation executions.

Simulation Interoperability↗

Improved microgrid resiliency through distributionally robust optimization under a policy-mode framework

Critical energy infrastructure are constantly under stress due to the ever increasing disruptions caused by wildfires, hurricanes, other weather related extreme events and cyber-attacks. Hence it becomes important to make critical infrastructure resilient to threats from such cyber-physical events. However, such events are hard to predict and numerous in nature and type and it becomes infeasible to make a system resilient to every possible such cyber-physical event. Such an approach can make the system operation overly conservative and impractical to operate. Furthermore, distributions of such events are hard to predict and historical data available on such events can be very sparse, making the problem even harder to solve. To deal with these issues, in this paper we present a policy-mode framework that enumerates and predicts the probability of various cyber-physical events and then a distributionally robust optimization (DRO) formulation that is robust to the sparsity of the available historical data. The proposed algorithm is illustrated on an islanded microgrid example: a modified IEEE 123-node feeder with distributed energy resources (DERs) and energy storage. Simulations are carried to validate the resiliency metrics under the sampled disruption events.

Nazir, Mohammad Nawaf↗

A Stochastic Intracellular Model of Anthrax Infection With Spore Germination Heterogeneity

We present a stochastic mathematical model of the intracellular infection dynamics of Bacillus anthracis in macrophages. Following inhalation of B. anthracis spores, these are ingested by alveolar phagocytes. Ingested spores then begin to germinate and divide intracellularly. This can lead to the eventual death of the host cell and the extracellular release of bacterial progeny. Some macrophages successfully eliminate the intracellular bacteria and will recover. Here, a stochastic birth-and-death process with catastrophe is proposed, which includes the mechanism of spore germination and maturation of B. anthracis . The resulting model is used to explore the potential for heterogeneity in the spore germination rate, with the consideration of two extreme cases for the rate distribution: continuous Gaussian and discrete Bernoulli. We make use of approximate Bayesian computation to calibrate our model using experimental measurements from in vitro infection of murine peritoneal macrophages with spores of the Sterne 34F2 strain of B. anthracis . The calibrated stochastic model allows us to compute the probability of rupture, mean time to rupture, and rupture size distribution, of a macrophage that has been infected with one spore. We also obtain the mean spore and bacterial loads over time for a population of cells, each assumed to be initially infected with a single spore. Our results support the existence of significant heterogeneity in the germination rate, with a subset of spores expected to germinate much later than the majority. Furthermore, in agreement with experimental evidence, our results suggest that most of the spores taken up by macrophages are likely to be eliminated by the host cell, but a few germinated spores may survive phagocytosis and lead to the death of the infected cell. Finally, we discuss how this stochastic modelling approach, together with dose-response data, allows us to quantify and predict individual infection risk following exposure.

59 BASIC BIOLOGICAL SCIENCES↗

The proximity of hotspots to convergent and divergent plate boundaries

An analysis of four different hotspot distributions, ranging from Morgan's (1972) original list of 19 to Vogt's (1981) list of 117 reveals that the hotspots are preferentially located near divergent plate boundaries. The probability of this proximity occurring by chance alone is quite remote, less than 0.01 for all four hotspot distributions. The same analysis also reveals that the hotspots are preferentially excluded from regions near convergent plate boundaries. The probability of this exclusion occurring by chance alone is 0.1 or less for three out of the four distributions examined. We interpret this behavior as being a consequence of the effects of large scale convective circulation on ascending mantle plumes. Mantle thermal plumes, the most probable source of hotspots, arise from instabilities in a basal thermal boundary layer. Plumes are suppressed from regions beneath convergent boundaries by descending flow and are entrained into the upwelling flow beneath spreading centers. Plate-scale convective circulation driven by subduction may also advect mantle thermal plumes toward spreading centers.

Weinstein, Stuart A.↗