Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “probability generating function”

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 217 records · Page 12

A REDUCED ORDER MODELING APPROACH TO PROBABILISTIC CREEP-DAMAGE PREDICTIONS IN FINITE ELEMENT ANALYSIS

This paper introduces a computationally efficient Reduced Order Modeling (ROM) approach for the probabilistic prediction of creep-damage failure. Component-level probabilistic simulations are needed to assess the reliability and safety of high-temperature components. Full-scale probabilistic creep-damage modeling in finite element (FE) approach is computationally expensive requiring many hundreds of simulations to replicate the uncertainty of component failure. To that end, ROM is proposed to minimize the elevated computational cost while controlling the loss of accuracy. It is proposed that full-scale probabilistic simulations can be completed in 1D at a reduced cost, the extremum conditions extracted, and those conditions applied for lower-cost 2D/3D probabilistic simulations of components that capture the mean and uncertainty of failure. The probabilistic Sine-hyperbolic (Sinh) model is selected which in previous work was calibrated to alloy 304 stainless steel. The Sinh model includes probability density functions (pdfs) for test condition (stress and temperature), initial damage (i.e. microstructure), and material properties uncertainty. The Sinh model is programmed into ANSYS finite element software using the USERCREEP.F material subroutine. First, the Sinh model and FE code are subject to verification and validation to ensure the accuracy of the simulations. Numerous Monte Carlo simulations are executed in a 1D model to generate probabilistic creep deformation, damage, and rupture data. This data is analyzed and the probabilistic parameters corresponding to extreme creep response are extracted. The ROM concept is applied where only the extreme conditions are applied in the 2D probabilistic prediction of a component. The probabilistic predictions between the 1D and 2D geometry is compared to assess ROM for creep. The accuracy of the probabilistic prediction employing the ROM approach will potentially reduce the time and cost of simulating complex engineering systems. Future studies will introduce multi-stage Sinh, stochasticity, and spatial uncertainty for improved prediction.

36 MATERIALS SCIENCE↗

Statistical Analysis of Fading Power Vectors for Real-Time Atmospheric Channel Emulation

The design and testing of a free-space optical communication system requires assessment of the impact of random fluctuations in received power from a laser beam transmitted over an atmospheric channel. A number of methods for generating fading power vectors for in-lab emulation of an atmospheric channel have previously been reported. These techniques include spectral shaping and filtering of a signal from a normally distributed pseudo-random number generator, full wave optics simulations with random phase screens, and pre-recorded measurements from experimental free-space links. In this work, a statistical analysis of atmospheric fading is presented with the goal of producing a practical engineering model for the fading channel suitable for generating synthetic fade vectors in real-time for long-duration receiver testing with channel interleaving. Specifically, a parametric model is developed for turbulence-induced fade on space-to-ground links with large-aperture receivers, including aperture-averaging and the effects of aperture size on the instantaneous coupling efficiency for mode- limited receivers. In particular, we analyze the probability density function and temporal power spectrum for fluctuations of the coupling efficiency for few-mode fibers in a range of turbulence conditions.

free-space optical communication↗

Stochastic decision analysis

Small space flight project design at NASA Langley Research Center goes through a multi-phase process from preliminary analysis to flight operations. The process insures that each system achieves its technical objectives with demonstrated quality and within planned budgets and schedules. A key technical component of early phases is decision analysis, which is a structure procedure for determining the best of a number of feasible concepts based upon project objectives. Feasible system concepts are generated by the designers and analyzed for schedule, cost, risk, and technical measures. Each performance measure value is normalized between the best and worst values and a weighted average score of all measures is calculated for each concept. The concept(s) with the highest scores are retained, while others are eliminated from further analysis. This project automated and enhanced the decision analysis process. Automation of the decision analysis process was done by creating a user-friendly, menu-driven, spreadsheet macro based decision analysis software program. The program contains data entry dialog boxes, automated data and output report generation, and automated output chart generation. The enhancements to the decision analysis process permit stochastic data entry and analysis. Rather than enter single measure values, the designers enter the range and most likely value for each measure and concept. The data can be entered at the system or subsystem level. System level data can be calculated as either sum, maximum, or product functions of the subsystem data. For each concept, the probability distributions are approximated for each measure and the total score for each concept as either constant, triangular, normal, or log-normal distributions. Based on these distributions, formulas are derived for the probability that the concept meets any given constraint, the probability that the concept meets all constraints, and the probability that the concept is within a given amount of the best score. Formulas are also derived for the probability that one concept's total score is within a given amount of a second concept's total score. These probabilistic calculations provide more realistic data entry and output information for designers, enabling designers to better determine which concepts to eliminate and which concepts to retain at the decision points of each design phase.

Lacksonen, Thomas A.↗

Statistical Analysis of Fading Power Vectors for Real-Time Atmospheric Channel Emulation

The design and testing of a free-space optical communication system requires assessment of the impact of random fluctuations in received power from a laser beam transmitted over an atmospheric channel. A number of methods for generating fading power vectors for in-lab emulation of an atmospheric channel have previously been reported. These techniques include spectral shaping and filtering of a signal from a normally distributed pseudo-random number generator, full wave optics simulations with random phase screens, and pre-recorded measurements from experimental free-space links. In this work, a statistical analysis of atmospheric fading is presented with the goal of producing a practical engineering model for the fading channel suitable for generating synthetic fade vectors in real-time for long-duration receiver testing with channel interleaving. Specifically, a parametric model is developed for turbulence-induced fade on space-to-ground links with large-aperture receivers, including aperture-averaging and the effects of aperture size on the instantaneous coupling efficiency for mode-limited receivers. In particular, we analyze the probability density function and temporal power spectrum for fluctuations of the coupling efficiency for few-mode fibers in a range of turbulence conditions.

free-space optical communication↗

Space shuttle solid rocket booster recovery system definition. Volume 2: SRB water impact Monte Carlo computer program, user's manual

The HD 220 program was created as part of the space shuttle solid rocket booster recovery system definition. The model was generated to investigate the damage to SRB components under water impact loads. The random nature of environmental parameters, such as ocean waves and wind conditions, necessitates estimation of the relative frequency of occurrence for these parameters. The nondeterministic nature of component strengths also lends itself to probabilistic simulation. The Monte Carlo technique allows the simultaneous perturbation of multiple independent parameters and provides outputs describing the probability distribution functions of the dependent parameters. This allows the user to determine the required statistics for each output parameter.

Source record↗

Computation of Domain-Averaged Shortwave Irradiance by a One-Dimensional Algorithm Incorporating Correlations between Optical Thickness and Direct Incident Radiation

A one-dimensional radiative transfer algorithm that accounts for correlations between the optical thickness and the incident direct solar radiation is developed to compute the domain-averaged shortwave irradiance profile. It divides the direct irradiance into four components and treats the direct irradiance in two separate, clear and cloudy columns to account for the fact that clouds attenuate the direct irradiance more than clear-sky. The horizontal inhomogeneity of clouds in the cloudy column is treated by the gamma weighted two-stream approximation, which assumes that the optical thickness of clouds follows a gamma distribution. The algorithm inputs the cloud fraction, cumulative cloud fraction as a function of height, and a parameter expressing the shape of the probability density function of the cloud optical thickness distribution in addition to inputs required for a two-stream radiative transfer model. These cloud property inputs can be obtained using ground- and satellite-based instruments. Therefore, the algorithm can treat realistic cloud overlap features and horizontal inhomogeneity of clouds in a framework of one- dimensional radiative transfer. Heating rates computed by the algorithm using cloud fields generated by cloud resolving models agree with those computed with a Monte Carlo model. If optical properties in computational layers that divide a vertically extensive cloud are correlated, the irradiance profile computed by the algorithm further improves.

Kato, S.↗

New method for estimating low-earth-orbit collision probabilities

An unconventional but general method is described for estimating the probability of collision between an earth-orbiting spacecraft and orbital debris. This method uses a Monte Caralo simulation of the orbital motion of the target spacecraft and each discrete debris object to generate an empirical set of distances, each distance representing the separation between the spacecraft and the nearest debris object at random times. Using concepts from the asymptotic theory of extreme order statistics, an analytical density function is fitted to this set of minimum distances. From this function, it is possible to generate realistic collision estimates for the spacecraft.

Vedder, John D.↗

Award-Winning CARES/Life Ceramics Durability Evaluation Software Is Making Advanced Technology Accessible

Products made from advanced ceramics show great promise for revolutionizing aerospace and terrestrial propulsion and power generation. However, ceramic components are difficult to design because brittle materials in general have widely varying strength values. The CARES/Life software developed at the NASA Lewis Research Center eases this by providing a tool that uses probabilistic reliability analysis techniques to optimize the design and manufacture of brittle material components. CARES/Life is an integrated package that predicts the probability of a monolithic ceramic component's failure as a function of its time in service. It couples commercial finite element programs--which resolve a component's temperature and stress distribution - with reliability evaluation and fracture mechanics routines for modeling strength - limiting defects. These routines are based on calculations of the probabilistic nature of the brittle material's strength.

Source record↗

Wavelet flow for extragalactic foreground simulations

Extragalactic foregrounds in cosmic microwave background (CMB) observations are both a source of cosmological and astrophysical information and a nuisance to the CMB. Effective field-level modeling that captures their non-Gaussian statistical distributions is increasingly important for optimal information extraction, particularly given the low-noise observations from current and upcoming experiments. Here, we explore the use of Wavelet Flow (WF) models to tackle the novel task of modeling the field-level probability distributions of multi-component CMB secondaries and foregrounds. Specifically, we jointly train correlated CMB lensing convergence (κ) and cosmic infrared background (CIB) maps with a WF model and obtain a network that statistically recovers the input to high accuracy — the trained network generates samples of κ and CIB fields whose average power spectra are within a few percent of the inputs across all scales, and whose Minkowski functionals are similarly accurate compared to the inputs. Leveraging the multiscale architecture of these models, we fine-tune both the model parameters and the priors at each scale independently, optimizing performance across different resolutions. These results demonstrate that WF models can accurately simulate correlated components of CMB secondaries, supporting improved analysis of cosmological data. Our code and trained models can be found on this GitHub repo.

cosmological simulations↗

Microwave backscattering theory and active remote sensing of the ocean surface

The status is reviewed of electromagnetic scattering theory relative to the interpretation of microwave remote sensing data acquired from spaceborne platforms over the ocean surface. Particular emphasis is given to the assumptions which are either implicit or explicit in the theory. The multiple scale scattering theory developed during this investigation is extended to non-Gaussian surface statistics. It is shown that the important statistic for the case is the probability density function of the small scale heights conditioned on the large scale slopes; this dependence may explain the anisotropic scattering measurements recently obtained with the AAFE Radscat. It is noted that present surface measurements are inadequate to verify or reject the existing scattering theories. Surface measurements are recommended for qualifying sensor data from radar altimeters and scatterometers. Additional scattering investigations are suggested for imaging type radars employing synthetically generated apertures.

Brown, G. S.↗

Magnetospheric Structures: Uranus and Neptune

Magnetospheric structures that might be encountered at Uranus and Neptune are described. Statistics indicate a sufficiently high probability to warrant consideration of their likely properties in advance of the Voyager encounters. Because the spin axis of Uranus lies nearly in the ecliptic and presently points approximately sunward, Voyager is likely to encounter the unique pole on configuration that has special theoretical significance. Corotation in the magnetospheres of Uranus and Neptune would probably exclude solar wind drive convection as an important driver of global magnetospheric dynamics, as it does at Jupiter and Saturn. The magnetospheres of Uranus and Neptune probably lack sufficient internal sources of plasma to produce significant levels of rotationally driven convection. The reported observation of auroral emission from Uranus has therefore motivated the development of an alternative model in which solar wind motion is coupled directly to the rotation of the ionosphere to establish a dynamo circuit which generates Birkeland currents and polar cap aurora. This model predicts the strength and configuration of the aurora as functions of the magnitude and polarity, respectively, of the planetary magnetic moment.

T W Hill↗

Transition probability of the Al II 2669 intersystem line

Time-resolved observations of the spin-changing, or 'intersystem' emission at 2669.157 A obtained by the ion storage technique are used to measure the transition probability of the 3s2 1S0 - 3s3p 3P1 exp 0 line in Al II. A laser-generated plasma was used as the source of the metastable Al(+) ions. The A-value result obtained for the intersystem transition is 3.33 + or - 0.23 x 10 to the 3rd/sec at the 90-percent confidence level; this value is used to derive two line-intensity ratios which involve the intersystem line as a function of electron density and temperature.

Johnson, B. C.↗

A Privacy-Preserving Distributed Control of Optimal Power Flow

Here, we consider a distributed optimal power flow formulated as an optimization problem that maximizes a nondifferentiable concave function. Solving such a problem by the existing distributed algorithms can lead to data privacy issues because the solution information exchanged within the algorithms can be utilized by an adversary to infer the data. To preserve data privacy, in this paper we propose a differentially private projected subgradient (DP-PS) algorithm that includes a solution encryption step. We show that a sequence generated by DP-PS converges in expectation, in probability, and with probability 1. Moreover, we show that the rate of convergence in expectation is affected by a target privacy level of DP-PS chosen by the user. We conduct numerical experiments that demonstrate the convergence and data privacy preservation of DP-PS.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Magnetospheric energization by interaction between planetary spin and the solar wind

If the solar wind is capable of driving magnetospheric convection, then solar-wind flow past any spinning, magnetized planet with a conducting ionosphere must cause the magnetic field lines in the outer part of its magnetospheric tail to be twisted into a helix. Such a magnetic field configuration requires magnetically field-aligned (Birkeland) currents in the tail that flow in and near the magnetopause and close by driving Pedersen currents through the planetary ionosphere. The strength of the Birkeland currents (and, by current continuity, the Pedersen currents) is, to first order, independent of the angle between the planetary-spin vector and the solar-wind velocity vector. Rather, the total current is a function of the magnetic moment of the planet, the radius of the tail, the angular velocity of planetary spin, the conductivity of the ionosphere, and the solar wind speed. For Jupiter, Saturn, Uranus, and perhaps Neptune, the power these currents deliver to the ionosphere is significant with regard to magnetospheric dynamics, such as the production of aurora and the generation of low-frequency radio emissions. For Mercury, Venus, earth, Mars, and probably Pluto, these currents are relatively small, although observable effects may be marginally detectable for the case of the earth's magnetosphere.

Isbell, J.↗

Clustering and Cliques in Preferential Attachment Random Graphs with Edge Insertion

In this paper, we investigate the global clustering coefficient (a.k.a transitivity) and clique number of graphs generated by a preferential attachment random graph model with an additional feature of allowing edge connections between existing vertices. Specifically, at each time step t, either a new vertex is added with probability f(t), or an edge is added between two existing vertices with probability 1 – f(t). We establish concentration inequalities for the global clustering and clique number of the resulting graphs under the assumption that f(t) is a regularly varying function at infinity with index of regular variation –$\gamma$, where $\gamma$ $\in$ [0, 1). Finally, we also demonstrate an inverse relation between these two statistics: the clique number is essentially the reciprocal of the global clustering coefficient.

97 MATHEMATICS AND COMPUTING↗

Generation of random geological models using multi-randomization for machine learning

Generating high-fidelity geological models is essential for advancing machine learning (ML) methods in automated seismic interpretation. For instance, seismic images paired with corresponding fault labels are foundational for ML-based fault detection from seismic migration sections. While several open-access datasets of random geological models exist, open-source tools specifically designed to produce large volumes of such models for ML applications remain scarce. To address this gap, we present RGM (Random Geological Model), an open-source software package for efficiently generating 2D and 3D synthetic geological models tailored for ML workflows. RGM supports the creation of diverse model components, including medium property distributions (P-/S-wave velocities and density), seismic reflectivity images (i.e., synthetic migration sections), relative geological time, and discrete fault attributes such as probability, dip, strike, rake, and displacement. It also accommodates the creation of complex geological features such as salt bodies and unconformities. The model generation algorithm employs a multi-randomization strategy, yielding an effectively infinite-dimensional model space that encompasses a wide range of geological scenarios and associated seismic features. Furthermore, RGM incorporates a method to generate synthetic elastic migration images using analytical elastic reflection coefficients combined with frequency-dependent scaling. This functionality enables the creation of training datasets for ML models that leverage elastic seismic images. RGM is implemented in modern object-oriented Fortran, allowing users to flexibly control statistical parameters governing model variability. We demonstrate the capability, performance, and geological realism of the package through comprehensive 2D and 3D examples.

58 GEOSCIENCES↗

Analysis of ground tests of a microwave, earth-occultation, pressure-reference-level system.

A two-satellite, microwave, earth-occultation system can supplement an infrared sensor by providing an accurate altitude reference that will serve to fix, as a function of height, the derived temperature profile of the infrared sounder. The results of ground tests made in Hawaii to estimate the likely effects of scintillation and fading on an occultation system are described. It was found that the microwave signal suffered periods of intense fading; extensive computer analyses of the data were performed in which aircraft-generated refractivity profiles were subjected to ray tracing. Results of these analyses indicate that the probable cause of the observed fading was multipath, a low-altitude phenomenon usually attributed to water vapor inhomogeneities. It is maintained that multipath will therefore have minimal effect on the pressure-reference microwave occultation system, which would operate at a relatively high closest-approach altitude (about 8 km).

Ungar, S. G.↗

Directed Design of Experiments for Validating Probability of Detection Capability of NDE Systems (DOEPOD)

The capability of an inspection system is established by applications of various methodologies to determine the probability of detection (POD). One accepted metric of an adequate inspection system is that there is 95% confidence that the POD is greater than 90% (90/95 POD). Directed design of experiments for probability of detection (DOEPOD) has been developed to provide an efficient and accurate methodology that yields observed POD and confidence bounds for both Hit- Miss or signal amplitude testing. Specifically, DOEPOD demands utilization of observance of occurrences. Directed DOEPOD does not assume prescribed POD logarithmic or similar functions with assumed adequacy over a wide range of flaw sizes and inspection system technologies, so that multi-parameter curve fitting or model optimization approaches to generate a POD curve are not required.

Generazio, Edward R.↗