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 919 records · Page 51

Reliability analysis in the Office of Safety, Environmental, and Mission Assurance (OSEMA)

The technical personnel in the SEMA office are working to provide the highest degree of value-added activities to their support of the NASA Langley Research Center mission. Management perceives that reliability analysis tools and an understanding of a comprehensive systems approach to reliability will be a foundation of this change process. Since the office is involved in a broad range of activities supporting space mission projects and operating activities (such as wind tunnels and facilities), it was not clear what reliability tools the office should be familiar with and how these tools could serve as a flexible knowledge base for organizational growth. Interviews and discussions with the office personnel (both technicians and engineers) revealed that job responsibilities ranged from incoming inspection to component or system analysis to safety and risk. It was apparent that a broad base in applied probability and reliability along with tools for practical application was required by the office. A series of ten class sessions with a duration of two hours each was organized and scheduled. Hand-out materials were developed and practical examples based on the type of work performed by the office personnel were included. Topics covered were: Reliability Systems - a broad system oriented approach to reliability; Probability Distributions - discrete and continuous distributions; Sampling and Confidence Intervals - random sampling and sampling plans; Data Analysis and Estimation - Model selection and parameter estimates; and Reliability Tools - block diagrams, fault trees, event trees, FMEA. In the future, this information will be used to review and assess existing equipment and processes from a reliability system perspective. An analysis of incoming materials sampling plans was also completed. This study looked at the issues associated with Mil Std 105 and changes for a zero defect acceptance sampling plan.

Kauffmann, Paul J.↗

Advanced Multi-Phase Flow CFD Model Development for Solid Rocket Motor Flowfield Analysis

It is known that the simulations of solid rocket motor internal flow field with AL-based propellants require complex multi-phase turbulent flow model. The objective of this study is to develop an advanced particulate multi-phase flow model which includes the effects of particle dynamics, chemical reaction and hot gas flow turbulence. The inclusion of particle agglomeration, particle/gas reaction and mass transfer, particle collision, coalescence and breakup mechanisms in modeling the particle dynamics will allow the proposed model to realistically simulate the flowfield inside a solid rocket motor. The Finite Difference Navier-Stokes numerical code FDNS is used to simulate the steady-state multi-phase particulate flow field for a 3-zone 2-D axisymmetric ASRM model and a 6-zone 3-D ASRM model at launch conditions. The 2-D model includes aft-end cavity and submerged nozzle. The 3-D model represents the whole ASRM geometry, including additional grain port area in the gas cavity and two inhibitors. FDNS is a pressure based finite difference Navier-Stokes flow solver with time-accurate adaptive second-order upwind schemes, standard and extended k-epsilon models with compressibility corrections, multi zone body-fitted formulations, and turbulence particle interaction model. Eulerian/Lagrangian multi-phase solution method is applied for multi-zone mesh. To simulate the chemical reaction, penalty function corrected efficient finite-rate chemistry integration method is used in FDNS. For the AL particle combustion rate, the Hermsen correlation is employed. To simulate the turbulent dispersion of particles, the Gaussian probability distribution with standard deviation equal to (2k/3)(exp 1/2) is used for the random turbulent velocity components. The computational results reveal that the flow field near the juncture of aft-end cavity and the submerged nozzle is very complex. The effects of the turbulent particles affect the flow field significantly and provide better prediction of the ASRM performance. The multi-phase flow analysis using the FDNS code in the present research can be used as a design tool for solid rocket motor applications.

Liaw, Paul↗

Probabilistic evaluation of fuselage-type composite structures

A methodology is developed to computationally simulate the uncertain behavior of composite structures. The uncertain behavior includes buckling loads, natural frequencies, displacements, stress/strain etc., which are the consequences of the random variation (scatter) of the primitive (independent random) variables in the constituent, ply, laminate and structural levels. This methodology is implemented in the IPACS (Integrated Probabilistic Assessment of Composite Structures) computer code. A fuselage-type composite structure is analyzed to demonstrate the code's capability. The probability distribution functions of the buckling loads, natural frequency, displacement, strain and stress are computed. The sensitivity of each primitive (independent random) variable to a given structural response is also identified from the analyses.

Shiao, Michael C.↗

Comparison of large-scale structures and velocities in the local universe

Comparison of the large-scale density and velocity fields in the local universe shows detailed agreement, strengthening the standard paradigm of the gravitational origin of these structures. Quantitative analysis can determine the cosmological density parameter, Omega, and biasing factor, b; there is virtually no sensitivity in any local analyses to the cosmologial constant, lambda. Comparison of the dipole anisotropy of the cosmic microwave background with the acceleration due to the Infrared Astronomy Satellite (IRAS) galaxies puts the linear growth factor in the range beta approximately equals Omega (exp 0.6)/b = 0.6(+0.7/-0.3) (95% confidence). A direct comparison of the density and velocity fields of nearby galaxies gives beta = 1.3 (+0.7/-0.6), and from nonlinear analysis the weaker limit (Omega greater than 0.45 for b greater than 0.5 (again 95% confidence). A tighter limit (Omega greater than 0.3 (4-6 sigma)), is obtained by a reconstruction of the probability distribution function of the initial fluctuations from which the structures observed today arose. The last two methods depend critically on the smooth velocity field determined from the observed velocities of nearby galaxies by the POTENT method. A new analysis of these velocities, with more than three times the data used to obtain the above quoted results, is now underway and promises to tighten the uncertainties considerably, as well as reduce systematic bias.

Yahil, Amos↗

The albedo of fractal stratocumulus clouds

An increase in the planetary albedo of the earth-atmosphere system by only 10% can decrease the equilibrium surface temperature to that of the last ice age. Nevertheless, albedo biases of 10% or greater would be introduced into large regions of current climate models if clouds were given their observed liquid water amounts, because of the treatment of clouds as plane parallel. The focus on marine stratocumulus clouds is due to their important role in cloud radiative forcing and also that, of the wide variety of earth's cloud types, they are most nearly plane parallel, so that they have the least albedo bias. The fractal model employed here reproduces both the probability distribution and the wavenumber spectrum of the stratocumulus liquid water path, as observed during the First ISCCP Regional Experiment (FIRE). A single new fractal parameter 0 less than or equal to f less than or equal to 1, is introduced and determined empirically by the variance of the logarithm of the vertically integrated liquid water. The reduced reflectivity of fractal stratocumulus clouds is approximately given by the plane-parallel reflectivity evaluated at a reduced 'effective optical thickness,' which when f = 0.5 is tau(sub eff) approximately equal to 10. Study of the diurnal cycle of stratocumulus liquid water during FIRE leads to a key unexpected result: the plane-parallel albedo bias is largest when the cloud fraction reaches 100%, that is, when any bias associated with the cloud fraction vanishes. This is primarily due to the variability increase with cloud fraction. Thus, the within-cloud fractal structure of stratocumulus has a more significant impact on estimates of its mesoscale-average albedo than does the cloud fraction.

Cahalan, Robert F.↗

Effects of a polar stratosphere cloud parameterization on ozone depletion due to stratospheric aircraft in a two-dimensional model

A parameterization of Type 1 and 2 polar stratospheric cloud (PSC) formation is presented which is appropriate for use in two-dimensional (2-D) photochemical models of the stratosphere. The calculations of PSC frequency of occurrence and surface area density uses climatological temperature probability distributions obtained from National Meteorological Center data to avoid using zonal mean temperatures, which are not good predictors of PSC behavior. The parameterization does not attempt to model the microphysics of PSCs. The parameterization predicts changes in PSC formation and heterogeneous processing due to perturbations of stratospheric trace constituents. It is therefore useful in assessing the potential effects of a fleet of stratospheric aircraft (high speed civil transports, or HSCTs) on stratospheric composition. the model calculated frequency of PSC occurrence agrees well with a climatology based on stratospheric aerosol measurement (SAM) 2 observations. PSCs are predicted to occur in the tropics. Their vertical range is narrow, however, and their impact on model O3 fields is small. When PSC and sulfate aerosol heterogeneous processes are included in the model calculations, the O3 change for 1980 - 1990 is in substantially better agreement with the total ozone mapping spectrometer (TOMS)-derived O3 trend than otherwise. The overall changes in model O3 response to standard HSCT perturbation scenarios produced by the parameterization are small and tend to decrease the model sensitivity to the HSCT perturbation. However, in the southern hemisphere spring a significant increase in O3 sensitivity to HSCT perturbations is found. At this location and time, increased PSC formation leads to increased levels of active chlorine, which produce the O3 decreases.

Considine, David B.↗

An efficient algorithm for choosing scattering directions in Monte Carlo work with arbitrary phase functions

This paper describes an efficient Monte Carlo algorithm for choosing a new direction of a photon after a scattering interaction. The algorithm chooses a scattering angle by linear interpolation in a table of the inverse cumulative scattering probability. A Legendre expansion of the phase function makes it easy to apply Clenshaw's algorithm to build the interpolation table. The points in the table are close enough together that linear interpolation is accurate. With a table of 100,000 entries, we can keep the absolute and relative errors in matching the probability distribution below 10(exp -5).

Barkstrom, Bruce R.↗

Pseudo-wavelet analysis of turbulence patterns in three vegetation layers

Ramp patterns in scalar traces such as temperature are the signature of coherent structures. A pseudo-wavelet analysis technique was developed in which ideal saw-tooth patterns of varying size were used as basis functions and fitted to temperature and velocity data. Data recorded from three very different vegetation stands were examined in this study. It was found that the most probable structure duration for the forest canopy was in the range 35-40 s, for the orchard canopy it was 20-25 s and for the maize it was 15-20 s. When expressed in non-dimensional form, the structure duration probability distribution for the maize canopy was about a decade larger than for the forest canopy, with the orchard canopy intermediate. The mean eddy duration versus wind shear relation falls on a narrow band for all three canopies, indicating that wind shear at the canopy top is the determining factor for the scale of the coherent eddies. The inverse of duration and intermittency of coherent structures exhibits a tendency of independence from wind shear at higher wind shear values. Coherent structures transport heat in a more efficient way than do smaller scale, less coherent motions. In all the canopies, the heat flux fractions associated with coherent structures are at least 10% higher than the corresponding time fraction.

Qiu, Jie↗

System-Availability And Resource-Allocation Program

ACARA analyzes availability, life-cycle cost, and scheduling of resources. Uses statistical Monte Carlo method to simulate capacity states of system as well as failure and repair of components. Failures of components modeled mathematically by use of combination of exponential and Weibull probability distributions. Schedules replacement of components to optimize performance of system. Made to comply with any constraints on production of components, capacities of resupply vehicles, spares kept on site, crews, and/or equipment. Written in APL2.

Viterna, L. A.↗

Electromagnetic Compatibility Testing Studies

This report discusses the results on analytical models and measurement and simulation of statistical properties from a study of microwave reverberation (mode-stirred) chambers performed at Texas Tech University. Two analytical models of power transfer vs. frequency in a chamber, one for antenna-to-antenna transfer and the other for antenna to D-dot sensor, were experimentally validated in our chamber. Two examples are presented of the measurement and calculation of chamber Q, one for each of the models. Measurements of EM power density validate a theoretical probability distribution on and away from the chamber walls and also yield a distribution with larger standard deviation at frequencies below the range of validity of the theory. Measurements of EM power density at pairs of points which validate a theoretical spatial correlation function on the chamber walls and also yield a correlation function with larger correlation length, R(sub corr), at frequencies below the range of validity of the theory. A numerical simulation, employing a rectangular cavity with a moving wall shows agreement with the measurements. The determination that the lowest frequency at which the theoretical spatial correlation function is valid in our chamber is considerably higher than the lowest frequency recommended by current guidelines for utilizing reverberation chambers in EMC testing. Two suggestions have been made for future studies related to EMC testing.

Trost, Thomas F.↗

Accurate Predictions of Mean Geomagnetic Dipole Excursion and Reversal Frequencies, Mean Paleomagnetic Field Intensity, and the Radius of Earth's Core Using McLeod's Rule

The geomagnetic spatial power spectrum R(sub n)(r) is the mean square magnetic induction represented by degree n spherical harmonic coefficients of the internal scalar potential averaged over the geocentric sphere of radius r. McLeod's Rule for the magnetic field generated by Earth's core geodynamo says that the expected core surface power spectrum (R(sub nc)(c)) is inversely proportional to (2n + 1) for 1 less than n less than or equal to N(sub E). McLeod's Rule is verified by locating Earth's core with main field models of Magsat data; the estimated core radius of 3485 kn is close to the seismologic value for c of 3480 km. McLeod's Rule and similar forms are then calibrated with the model values of R(sub n) for 3 less than or = n less than or = 12. Extrapolation to the degree 1 dipole predicts the expectation value of Earth's dipole moment to be about 5.89 x 10(exp 22) Am(exp 2)rms (74.5% of the 1980 value) and the expected geomagnetic intensity to be about 35.6 (mu)T rms at Earth's surface. Archeo- and paleomagnetic field intensity data show these and related predictions to be reasonably accurate. The probability distribution chi(exp 2) with 2n+1 degrees of freedom is assigned to (2n + 1)R(sub nc)/(R(sub nc). Extending this to the dipole implies that an exceptionally weak absolute dipole moment (less than or = 20% of the 1980 value) will exist during 2.5% of geologic time. The mean duration for such major geomagnetic dipole power excursions, one quarter of which feature durable axial dipole reversal, is estimated from the modern dipole power time-scale and the statistical model of excursions. The resulting mean excursion duration of 2767 years forces us to predict an average of 9.04 excursions per million years, 2.26 axial dipole reversals per million years, and a mean reversal duration of 5533 years. Paleomagnetic data show these predictions to be quite accurate. McLeod's Rule led to accurate predictions of Earth's core radius, mean paleomagnetic field intensity, and mean geomagnetic dipole power excursion and axial dipole reversal frequencies. We conclude that McLeod's Rule helps unify geo-paleomagnetism, correctly relates theoretically predictable statistical properties of the core geodynamo to magnetic observation, and provides a priori information required for stochastic inversion of paleo-, archeo-, and/or historical geomagnetic measurements.

Voorhies, Coerte V.↗

Bit-Wise Arithmetic Coding For Compression Of Data

Bit-wise arithmetic coding is data-compression scheme intended especially for use with uniformly quantized data from source with Gaussian, Laplacian, or similar probability distribution function. Code words of fixed length, and bits treated as being independent. Scheme serves as means of progressive transmission or of overcoming buffer-overflow or rate constraint limitations sometimes arising when data compression used.

Kiely, Aaron↗

A Human Factors Analysis of EVA Time Requirements

Human Factors Engineering (HFE) is a discipline whose goal is to engineer a safer, more efficient interface between humans and machines. HFE makes use of a wide range of tools and techniques to fulfill this goal. One of these tools is known as motion and time study, a technique used to develop time standards for given tasks. During the summer of 1995, a human factors motion and time study was initiated with the goals of developing a database of EVA task times and developing a method of utilizing the database to predict how long an EVA should take. Initial development relied on the EVA activities performed during the STS-61 (Hubble) mission. The first step of the study was to become familiar with EVA's, the previous task-time studies, and documents produced on EVA's. After reviewing these documents, an initial set of task primitives and task-time modifiers was developed. Data was collected from videotaped footage of two entire STS-61 EVA missions and portions of several others, each with two EVA astronauts. Feedback from the analysis of the data was used to further refine the primitives and modifiers used. The project was continued during the summer of 1996, during which data on human errors was also collected and analyzed. Additional data from the STS-71 mission was also collected. Analysis of variance techniques for categorical data was used to determine which factors may affect the primitive times and how much of an effect they have. Probability distributions for the various task were also generated. Further analysis of the modifiers and interactions is planned.

Pate, Dennis W.↗

Demonstration of a Probabilistic Technique for the Determination of Economic Viability of Very Large Transport Configurations

Over the past few years, modem aircraft design has experienced a paradigm shift from designing for performance to designing for affordability. This report contains a probabilistic approach that will allow traditional deterministic design methods to be extended to account for disciplinary, economic, and technological uncertainty. The probabilistic approach was facilitated by the Fast Probability Integration (FPI) technique; a technique which allows the designer to gather valuable information about the vehicle's behavior in the design space. This technique is efficient for assessing multi-attribute, multi-constraint problems in a more realistic fashion. For implementation purposes, this technique is applied to illustrate how both economic and technological uncertainty associated with a Very Large Transport aircraft concept may be assessed. The assessment is evaluated with the FPI technique to determine the cumulative probability distributions of the design space, as bound by economic objectives and performance constraints. These distributions were compared to established targets for a comparable large capacity aircraft, similar in size to the Boeing 747-400. The conventional baseline configuration design space was determined to be unfeasible and marginally viable, motivating the infusion of advanced technologies, including reductions in drag, specific fuel consumption, wing weight, and Research, Development, Testing, and Evaluation costs. The resulting system design space was qualitatively assessed with technology metric "k" factors. The infusion of technologies shifted the VLT design into regions of feasibility and greater viability. The study also demonstrated a method and relationship by which the impact of new technologies may be assessed in a more system focused approach.

Mavris, Dimitri N.↗

Probabilistic Modeling of Ceramic Matrix Composite Strength

Uncertainties associated with the primitive random variables such as manufacturing process (processing temperature, fiber volume ratio, void volume ratio), constituent properties (fiber, matrix and interface), and geometric parameters (ply thickness, interphase thickness) have been simulated to quantify the scatter in the first matrix cracking strength (FMCS) and the ultimate tensile strength of SCS-6/RBSN (SiC fiber (SCS-6) reinforced reaction-bonded silicon nitride composite) ceramic matrix composite laminate at room temperature. Cumulative probability distribution function for the FMCS and ultimate tensile strength at room temperature (RT) of (0)(sub 8), (0(sub 2)/90(sub 2), and (+/-45(sub 2))(sub S) laminates have been simulated and the sensitivity of primitive variables to the respective strengths have been quantified. Computationally predicted scatter of the strengths for a uniaxial laminate have been compared with those from limited experimental data. Also the experimental procedure used in the tests has been described briefly. Results show a very good agreement between the computational simulation and the experimental data. Dominating failure modes in (0)(sub 8), (0/90)(sub s) and (+/-45)(sub S) laminates have been identified. Results indicate that the first matrix cracking strength for the (0)(sub S), and (0/90)(sub S) laminates is sensitive to the thermal properties, modulus and strengths of both the fiber and matrix whereas the ultimate tensile strength is sensitive to the fiber strength and the fiber volume ratio. In the case of a (+/-45)(sub S), laminate, both the FMCS and the ultimate tensile strengths have a small scatter range and are sensitive to the fiber tensile strength as well as the fiber volume ratio.

Shan, Ashwin R.↗

Design Protocols and Analytical Strategies that Incorporate Structural Reliability Models

The general goal of this project is to establish design protocols that enable the engineer to analyze and predict certain types of behavior in ceramic composites. Sections of the final report addresses the following: Description of the Problem that Motivated the Technology Development, Description of the New Technology that was Developed, Unique and Novel Features of the Technology and Results/Benefits of Application (year by year accomplishments), and Utilization of New Technology in Non-Aerospace Applications. Activities for this reporting period included the development of a design analysis as part of a cooperative agreement with general Electric Aircraft Engines. The effort focused on modifying the Toughened Ceramics Analysis and Reliability Evaluation of Structures (TCARES) algorithm for use in the design of engine components fabricated from NiAl. Other activities related to the development of an ASTM standard practice for estimating Weibull parameters. The standard focuses on the evaluation and reporting of uniaxial strength data, and the estimation of probability distribution parameters for ceramics which fail in a brittle fashion.

Duffy, Stephen F.↗

On Geomagnetism and Paleomagnetism

A statistical description of Earth's broad scale, core-source magnetic field has been developed and tested. The description features an expected, or mean, spatial magnetic power spectrum that is neither "flat" nor "while" at any depth, but is akin to spectra advanced by Stevenson and McLeod. This multipole spectrum describes the magnetic energy range; it is not steep enough for Gubbins' magnetic dissipation range. Natural variations of core multipole powers about their mean values are to be expected over geologic time and are described via trial probability distribution functions that neither require nor prohibit magnetic isotropy. The description is thus applicable to core-source dipole and low degree non-dipole fields despite axial dipole anisotropy. The description is combined with main field models of modem satellite and surface geomagnetic measurements to make testable predictions of: (1) the radius of Earth's core, (2) mean paleomagnetic field intensity, and (3) the mean rates and durations of both dipole power excursions and durable axial dipole reversals. The predicted core radius is 0.7% above the 3480 km seismologic value. The predicted root mean square paleointensity (35.6 mu T) and mean Virtual Axial Dipole Moment (about 6.2 lx 1022 Am(exp 2)) are within the range of various mean paleointensity estimates. The predicted mean rate of dipole power excursions, as defined by an absolute dipole moment <20% of the 1980 value, is 9.04/Myr and 14% less than obtained by analysis of a 4 Myr paleointensity record. The predicted mean rate of durable axial dipole reversals (2.26/Myr) is 2.3% more than established by the polarity time-scale for the past 84 Myr. The predicted mean duration of axial dipole reversals (5533 yr) is indistinguishable from an observational value. The accuracy of these predictions demonstrates the power and utility of the description, which is thought to merit further development and testing. It is suggested that strong stable stratification of Earth's uppermost outer core leads to a geologically long interval of no dipole reversals and a very nearly axisymmetric field outside the core. Statistical descriptions of other planetary magnetic fields are outlined.

Voorhies, Coerte V.↗

Aircraft NO(x) had no Unique Fingerprint on Sonex; Lightning Dominated Fresh NO(x) Sources

Key questions to which SONEX was directed were the following: Can aircraft corridors be detected? Is there a unique tracer for aircraft NO(x)? Can a "background" NO(x) (or NO(y) be defined? What fraction of NO(x) measured during SONEX was from aircraft? How representative was SONEX of the North Atlantic in 1997 and how typical of other years? We attempt to answer these questions through species-species correlations, probability distribution functions (PDFs), and meteorological history. There is not a unique aircraft tracer, largely due to the high variability of air mass origins and tracer ratios, which render "average" quantities meaningless. The greatest NO and NO(y) signals were associated with lightning and convective NO sources. Well-defined background CO, NO(y) and NO(y)/ozone ratio appear in subsets of two cross-track flights with subtropical origins and five flights with predominantly mid-latitude air. Forty percent of the observations on these 7 flights showed NO(y)/ozone to be above background, evidently due to unreacted NO(x). This NO(x) is a combination of aircraft, lightning and surface pollution injected by convection. The strongly subtropical signatures in SONEX observations, confirmed by pv (potential vorticity) values along flight tracks, argues for most of the unreacted NO(x) originating from lightning. Potential vorticity statistics along SONEX flight tracks in 1992-1998, and for the North Atlantic as a whole, show the SONEX meteorological environment to be representative of the North Atlantic flight corridor in the October-November period.

Thompson, A.↗