Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stochastic modelling”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Blob velocities and sizes in the Alcator C-Mod scrape-off layer for ohmic and high confinement mode plasmas

An improved time delay estimation method is used to calculate the velocity of cross-field blob motion in the scrape-off layer of Alcator C-Mod for an ohmic and two high confinement (H-mode) plasmas; an edge localized mode free and an enhanced D-alpha H-mode. The gas puff imaging data analysis results are interpreted in the framework of a stochastic model that describes the fluctuations as a super-position of uncorrelated blob-like structures. In all confinement modes investigated, the scrape-off layer is dominated by large amplitude, blob-like filaments moving radially outwards with velocities in the range from 400 to 1000 m/s. Blobs in high confinement mode plasmas have similar velocities and sizes as in ohmic plasma, which is consistent with the close similarity of conditionally averaged burst shapes and frequency spectra for the confinement modes investigated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Phase transitions in Schloegl's second model for autocatalysis on a Bethe lattice

Schloegl's second model (also known as the quadratic contact process) on a lattice involves spontaneous particle annihilation at rate p and autocatalytic particle creation at empty sites with n ≥ 2 occupied neighbors. The particle creation rate for exactly n occupied neighbors is selected here as n(n - 1)/[z(z - 1)] for lattice coordination number z. We analyze this model on a Bethe lattice. Precise behavior for stochastic models on regular periodic infinite lattices is usually surmised from kinetic Monte Carlo simulation on a finite lattice with periodic boundary conditions. However, the persistence of boundary effects for a Bethe lattice complicates this process, e.g., by inducing spatially heterogenous states. This motivates the exploration of various boundary conditions and unconventional simulation ensembles on the Bethe lattice to predict behavior for infinite size. Here, we focus on z = 3, and predict a discontinuous transition to the vacuum state on the infinite lattice when p exceeds a threshold value of around 0.053.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2009

Development of a GCR Event-based Risk Model

A goal at NASA is to develop event-based systems biology models of space radiation risks that will replace the current dose-based empirical models. Complex and varied biochemical signaling processes transmit the initial DNA and oxidative damage from space radiation into cellular and tissue responses. Mis-repaired damage or aberrant signals can lead to genomic instability, persistent oxidative stress or inflammation, which are causative of cancer and CNS risks. Protective signaling through adaptive responses or cell repopulation is also possible. We are developing a computational simulation approach to galactic cosmic ray (GCR) effects that is based on biological events rather than average quantities such as dose, fluence, or dose equivalent. The goal of the GCR Event-based Risk Model (GERMcode) is to provide a simulation tool to describe and integrate physical and biological events into stochastic models of space radiation risks. We used the quantum multiple scattering model of heavy ion fragmentation (QMSFRG) and well known energy loss processes to develop a stochastic Monte-Carlo based model of GCR transport in spacecraft shielding and tissue. We validated the accuracy of the model by comparing to physical data from the NASA Space Radiation Laboratory (NSRL). Our simulation approach allows us to time-tag each GCR proton or heavy ion interaction in tissue including correlated secondary ions often of high multiplicity. Conventional space radiation risk assessment employs average quantities, and assumes linearity and additivity of responses over the complete range of GCR charge and energies. To investigate possible deviations from these assumptions, we studied several biological response pathway models of varying induction and relaxation times including the ATM, TGF -Smad, and WNT signaling pathways. We then considered small volumes of interacting cells and the time-dependent biophysical events that the GCR would produce within these tissue volumes to estimate how GCR event rates mapped to biological signaling induction and relaxation times. We considered several hypotheses related to signaling and cancer risk, and then performed simulations for conditions where aberrant or adaptive signaling would occur on long-duration space mission. Our results do not support the conventional assumptions of dose, linearity and additivity. A discussion on how event-based systems biology models, which focus on biological signaling as the mechanism to propagate damage or adaptation, can be further developed for cancer and CNS space radiation risk projections is given.

Cucinotta, Francis A.↗
DOE OSTI2025

Resolving turbulent magnetohydrodynamics: a hybrid operator-diffusion framework

We present a hybrid machine learning framework that combines physics-informed neural operators (PINOs) with score-based generative diffusion models to simulate the full spatio-temporal evolution of two-dimensional, incompressible, resistive magnetohydrodynamic turbulence across a broad range of Reynolds numbers (Re). The framework leverages the equation-constrained generalization capabilities of PINOs to predict coherent, low-frequency dynamics, while a conditional diffusion model stochastically corrects high-frequency residuals, enabling accurate modeling of fully developed turbulence. Trained on a comprehensive ensemble of high-fidelity simulations with Re ϵ {100, 250, 500, 750, 1000, 3000, 10000}, the approach achieves state-of-the-art accuracy in regimes previously inaccessible to deterministic surrogates. At Re = 1000 and 3000, the model faithfully reconstructs the full spectral energy distributions of both velocity and magnetic fields late into the simulation, capturing non-Gaussian statistics, intermittent structures, and cross-field correlations with high fidelity. At extreme turbulence levels (Re = 10 000), it remains the first surrogate capable of recovering the high-wavenumber evolution of the magnetic field, preserving large-scale morphology and enabling statistically meaningful predictions.

Diffusion-Integrated Neural Operators↗
NASA NTRS1978

Analytical techniques for the study of some parameters of multispectral scanner systems for remote sensing

The author has identified the following significant results. The concept of average mutual information in the received spectral random process about the spectral scene was developed. Techniques amenable to implementation on a digital computer were also developed to make the required average mutual information calculations. These techniques required identification of models for the spectral response process of scenes. Stochastic modeling techniques were adapted for use. These techniques were demonstrated on empirical data from wheat and vegetation scenes.

Wiswell, E. R.↗
NASA NTRS1978

Simplified model of statistically stationary spacecraft rotation and associated induced gravity environments

A stochastic model of spacecraft motion was developed based on the assumption that the net torque vector due to crew activity and rocket thruster firings is a statistically stationary Gaussian vector process. The process had zero ensemble mean value, and the components of the torque vector were mutually stochastically independent. The linearized rigid-body equations of motion were used to derive the autospectral density functions of the components of the spacecraft rotation vector. The cross-spectral density functions of the components of the rotation vector vanish for all frequencies so that the components of rotation were mutually stochastically independent. The autospectral and cross-spectral density functions of the induced gravity environment imparted to scientific apparatus rigidly attached to the spacecraft were calculated from the rotation rate spectral density functions via linearized inertial frame to body-fixed principal axis frame transformation formulae. The induced gravity process was a Gaussian one with zero mean value. Transformation formulae were used to rotate the principal axis body-fixed frame to which the rotation rate and induced gravity vector were referred to a body-fixed frame in which the components of the induced gravity vector were stochastically independent. Rice's theory of exceedances was used to calculate expected exceedance rates of the components of the rotation and induced gravity vector processes.

Fichtl, G. H.↗
DOE OSTI2026

Thermodynamically consistent incorporation of the Langmuir adsorption model into compressible fluctuating hydrodynamics

For a gas–solid interfacial system where chemical species undergo reversible adsorption, we develop a mesoscopic stochastic modeling method that simulates both gas-phase hydrodynamics and surface coverage dynamics by coupling the Langmuir adsorption model with compressible fluctuating hydrodynamics. To this end, we derive a thermodynamically consistent mass–energy update scheme that accounts for how the mass and energy variables in the gas and surface subsystems should be updated according to the changes in the number of molecules of each species in each subsystem due to adsorption and desorption events. By performing a stochastic analysis for the ideal Langmuir model and the full hydrodynamic system, we analytically confirm that our mass–energy update scheme captures thermodynamic equilibrium predicted by equilibrium statistical mechanics. We find that an internal energy correction term is needed, which is attributed to the difference in the mean kinetic energy of gas molecules colliding with the surface from that computed from the Maxwell–Boltzmann distribution. By performing an equilibrium simulation study for an ideal gas mixture of CO and Ar, with CO undergoing reversible adsorption, we validate our overall simulation method and implementation.

Adsorption↗
NASA NTRS1984

Data and results of a laboratory investigation of microprocessor upset caused by simulated lightning-induced analog transients

Advanced composite aircraft designs include fault-tolerant computer-based digital control systems with thigh reliability requirements for adverse as well as optimum operating environments. Since aircraft penetrate intense electromagnetic fields during thunderstorms, onboard computer systems maya be subjected to field-induced transient voltages and currents resulting in functional error modes which are collectively referred to as digital system upset. A methodology was developed for assessing the upset susceptibility of a computer system onboard an aircraft flying through a lightning environment. Upset error modes in a general-purpose microprocessor were studied via tests which involved the random input of analog transients which model lightning-induced signals onto interface lines of an 8080-based microcomputer from which upset error data were recorded. The application of Markov modeling to upset susceptibility estimation is discussed and a stochastic model development.

Belcastro, C. M.↗
NASA NTRS2001

Cloud Structure and Cloud-Radiation Interaction: From Smallest Observable Scales (Less than 10 m) to GCM Grid Scales (More than 100 km)

In this presentation we review the fractal nature of internal cloud structure from cm- to km-scales as captured by in-situ probes during long horizontal penetrations by aircraft. We uncover the non-Poissonian spatial distribution of cloud droplets at submeter scales and confirm scale-invariant behavior for large scales. Based on these structural characteristics, we generate simple fractal cloud models that reproduce statistical scaling properties of real clouds. These stochastic models represent a link between nonlinear science, in general, and cloud-radiation interaction, in particular. Next we run three-dimensional radiative transfer computations on these synthetic fractal clouds and compare the structure of the resulting radiation fields with the known structure of the cloud model and with satellite images of real clouds. The different behaviors observed for small and large-scale variabilities will be discussed in detail. We find that while the large-scale fluctuations of the resulting radiation fields resemble those in the original scale-invariant cloud structure, the radiation at small scales is much smoother than its cloud liquid water counterpart. This violates scale-invariance and produces a scale-break at 0.2-0.5 km that is clearly observed in high-resolution satellite data such as from Landsat. Finally, we show how radiative transfer Green function theory in the photon diffusion limit explains (and predicts) the above phenomena of "radiative smoothing."

Marshak, Alexander↗
NASA NTRS1985

Measurements and predictions of the structure of evaporating sprays

In the present experimental and theoretical study of turbulent, evaporating sprays, round, Freon-11 sprays produced by an air-atomizing injector that is directed vertically downward in still air were subjected to structure measurements for mean and fluctuating gas velocities, total Freon-11 concentration, drop size and velocity contributions, mean gas temperature, and liquid flux distributions. An evaluation was then conducted of three spray models: (1) locally homogeneous flow, (2) deterministic separated flow, and (3) stochastic separated flow. The first two of these were found to perform poorly; the stochastic model yielded the best agreement between predictions and measurements.

Solomon, A. S. P.↗
NASA NTRS1996

Two-stage model of radon-induced malignant lung tumors in rats: effects of cell killing

A two-stage stochastic model of carcinogenesis is used to analyze lung tumor incidence in 3750 rats exposed to varying regimens of radon carried on a constant-concentration uranium ore dust aerosol. New to this analysis is the parameterization of the model such that cell killing by the alpha particles could be included. The model contains parameters characterizing the rate of the first mutation, the net proliferation rate of initiated cells, the ratio of the rates of cell loss (cell killing plus differentiation) and cell division, and the lag time between the appearance of the first malignant cell and the tumor. Data analysis was by standard maximum likelihood estimation techniques. Results indicate that the rate of the first mutation is dependent on radon and consistent with in vitro rates measured experimentally, and that the rate of the second mutation is not dependent on radon. An initial sharp rise in the net proliferation rate of initiated cell was found with increasing exposure rate (denoted model I), which leads to an unrealistically high cell-killing coefficient. A second model (model II) was studied, in which the initial rise was attributed to promotion via a step function, implying that it is due not to radon but to the uranium ore dust. This model resulted in values for the cell-killing coefficient consistent with those found for in vitro cells. An "inverse dose-rate" effect is seen, i.e. an increase in the lifetime probability of tumor with a decrease in exposure rate. This is attributed in large part to promotion of intermediate lesions. Since model II is preferable on biological grounds (it yields a plausible cell-killing coefficient), such as uranium ore dust. This analysis presents evidence that a two-stage model describes the data adequately and generates hypotheses regarding the mechanism of radon-induced carcinogenesis.

NASA Program Radiation Health↗
NASA NTRS2006

Towards a Theory of Sampled-Data Piecewise-Deterministic Markov Processes

The analysis and design of practical control systems requires that stochastic models be employed. Analysis and design tools have been developed, for example, for Markovian jump linear continuous and discrete-time systems, piecewise-deterministic processes (PDP's), and general stochastic hybrid systems (GSHS's). These model classes have been used in many applications, including fault tolerant control and networked control systems. This paper presents initial results on the analysis of a sampled-data PDP representation of a nonlinear sampled-data system with a jump linear controller. In particular, it is shown that the state of the sampled-data PDP satisfies the strong Markov property. In addition, a relation between the invariant measures of a sampled-data system driven by a stochastic process and its associated discrete-time representation are presented. As an application, when the plant is linear with no external input, a sufficient testable condition for the convergence in distribution to the invariant delta Dirac measure is given.

Herencia-Zapana, Heber↗
NASA NTRS1995

Failure models for textile composites

The goals of this investigation were to: (1) identify mechanisms of failure and determine how the architecture of reinforcing fibers in 3D woven composites controlled stiffness, strength, strain to failure, work of fracture, notch sensitivity, and fatigue life; and (2) to model composite stiffness, strength, and fatigue life. A total of 11 different angle and orthogonal interlock woven composites were examined. Composite properties depended on the weave architecture, the tow size, and the spatial distributions and strength of geometrical flaws. Simple models were developed for elastic properties, strength, and fatigue life. A more complicated stochastic model, called the 'Binary Model,' was developed for damage tolerance and ultimate failure. These 3D woven composites possessed an extraordinary combination of strength, damage tolerance, and notch insensitivity.

Cox, Brian↗
NASA NTRS2010

Fiber Breakage Model for Carbon Composite Stress Rupture Phenomenon: Theoretical Development and Applications

Stress rupture failure of Carbon Composite Overwrapped Pressure Vessels (COPVs) is of serious concern to Science Mission and Constellation programs since there are a number of COPVs on board space vehicles with stored gases under high pressure for long durations of time. It has become customary to establish the reliability of these vessels using the so called classic models. The classical models are based on Weibull statistics fitted to observed stress rupture data. These stochastic models cannot account for any additional damage due to the complex pressure-time histories characteristic of COPVs being supplied for NASA missions. In particular, it is suspected that the effects of proof test could significantly reduce the stress rupture lifetime of COPVs. The focus of this paper is to present an analytical appraisal of a model that incorporates damage due to proof test. The model examined in the current paper is based on physical mechanisms such as micromechanics based load sharing concepts coupled with creep rupture and Weibull statistics. For example, the classic model cannot accommodate for damage due to proof testing which every flight vessel undergoes. The paper compares current model to the classic model with a number of examples. In addition, several applications of the model to current ISS and Constellation program issues are also examined.

Murthy, Pappu L. N.↗
DOE OSTI2023

Subordinated Gaussian Processes for Solar Irradiance

Traditionally the power grid has been a one-way street with power flowing from large transmission-connected generators through the distribution network to consumers. This paradigm is changing with the introduction of distributed renewable energy resources (DERs), and with it, the way the grid is managed. There is currently a dearth of high fidelity solar irradiance datasets available to help grid researchers understand how expansion of DERs could affect future power system operations. Realistic simulations of by-the-second solar irradiances are needed to study how DER variability affects the grid. Irradiance data are highly non-stationary and non-Gaussian, and even modern time series models are challenged by their distributional properties. We develop a subordinated non-Gaussian stochastic model whose simulations realistically capture the distribution and dependence structure in measured irradiance. We illustrate our approach on a fine resolution dataset from Hawaii, where our approach outperforms standard nonlinear time series models.

MATHEMATICS AND COMPUTING,SOLAR ENERGY↗
NASA NTRS1977

Models and techniques for evaluating the effectiveness of aircraft computing systems

Models, measures and techniques were developed for evaluating the effectiveness of aircraft computing systems. The concept of effectiveness involves aspects of system performance, reliability and worth. Specifically done was a detailed development of model hierarchy at mission, functional task, and computational task levels. An appropriate class of stochastic models was investigated which served as bottom level models in the hierarchial scheme. A unified measure of effectiveness called 'performability' was defined and formulated.

Meyer, J. F.↗
NASA NTRS2004

Designing to Mitigate Food Growing Failures in Space

Future space life support systems may use crop plants to grow most of the crew s food. A harvest failure can reduce the food available for future consumption. If the previously stored food is insufficient to reach the next harvest, the crew may go hungry. This paper considers how the overall food supply system should be modified to cope with food production failures. The food supply concept for a mission will use grown food, or stored food, cIr both. The optimum food supply mix depends on the costs and failure probabilities of stored and grown food. A simple food system model assumes that either we obtain the nominal harvest or a failure occurs and no food is harvested. Given the probability that any particular harvest fails, it is easy to compute the expected number of failures and the total food shortfall over a mission. If some food is grown and the probability of harvest failure is high, a non-redundant system has an unacceptable likelihood that the crew will have no food for a full harvest period. Food supply reliability must be increased either by supplying more food initially or by increasing food production capacity. We can obtain a very reliable food supply even when the harvest failure rate is high. If the cost of growing food is much less than the cost of providing stored food, it is better to provide redundant food growing capacity than to increase initial storage. A more realistic biomass production failure model allows the harvest amount or time to vary around the nominal values, using stochastic modeling with repeated Monte Carlo simulation, but such failures have minor impact compared to a complete harvest failure.

Jones, Harry↗
NASA NTRS1993

Narrow band noise as a model of time-dependent accelerations - Study of the stability of a fluid surface in a microgravity environment

We introduce a stochastic model to analyze in quantitative detail the effect of the high frequency components of the residual accelerations onboard spacecraft (often called g-jitter) on fluid motion. The residual acceleration field is modeled as a narrow band noise characterized by three independent parameters: its intensity G squared, a dominant frequency Omega, and a characteristic spectral width tau exp -1. The white noise limit corresponds to Omega tau goes to O, with G squared tau finite, and the limit of a periodic g-jitter (or deterministic limit) can be recovered for Omega tau goes to infinity, G squared finite. The analysis of the response of a fluid surface subjected to a fluctuating gravitational field leads to the stochastic Mathieu equation driven by both additive and multiplicative noise. We discuss the stability of the solutions of this equation in the two limits of white noise and deterministic forcing, and in the general case of narrow band noise. The results are then applied to typical microgravity conditions.

Casademunt, Jaume↗