Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stochastic Modeling and Simulation”

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 OSTI2017

Diffusion of two-dimensional epitaxial clusters on metal (100) surfaces: Facile versus nucleation-mediated behavior and their merging for larger sizes

For diffusion of two-dimensional homoepitaxial clusters of N atoms on metal(100) surfaces mediated by edge atom hopping, macroscale continuum theory suggests that the diffusion coefficient scales like DN ~ N -β with β = 3/2. However, we find quite different and diverse behavior in multiple size regimes. These include: (i) facile diffusion for small sizes N < 9; (ii) slow nucleation-mediated diffusion with small β < 1 for “perfect” sizes N = N p = L 2 or L(L+1), for L = 3, 4,… having unique ground state shapes, for moderate sizes 9 ≤ N ≤ O(10 2 ); the same also applies for N = N p +3, N p + 4,… (iii) facile diffusion but with large β > 2 for N = Np + 1 and N p + 2 also for moderate sizes 9 ≤ N ≤ O(10 2 ); (iv) merging of the above distinct branches and subsequent anomalous scaling with 1 ≲ β < 3/2, reflecting the quasi-facetted structure of clusters, for larger N = O(10 2 ) to N = O(10 3 ); and (v) classic scaling with β = 3/2 for very large N = O(103) and above. The specified size ranges apply for typical model parameters. We focus on the moderate size regime where show that diffusivity cycles quasi-periodically from the slowest branch for N p + 3 (not Np) to the fastest branch for Np + 1. Behavior is quantified by Kinetic Monte Carlo simulation of an appropriate stochastic lattice-gas model. However, precise analysis must account for a strong enhancement of diffusivity for short time increments due to back-correlation in the cluster motion. Further understanding of this enhancement, of anomalous size scaling behavior, and of the merging of various branches, is facilitated by combinatorial analysis of the number of the ground state and low-lying excited state cluster configurations, and also of kink populations.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Numerical analysis of randomly forced glycolitic oscillations

Randomly forced glycolytic oscillations in Higgins model are studied both numerically and analytically. Numerical analysis is based on the direct simulation of the solutions of stochastic system. Non-uniformity of the stochastic bundle along the deterministic cycle is shown. For the analytical investigation of the randomly forced Higgins model, the stochastic sensitivity function technique and confidence domains method are applied. Results of the influence of additive noise on the cycle of this model are given.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Learning stochastic dynamics with statistics-informed neural network

We introduce a machine-learning framework named statistics-informed neural network (SINN) for learning stochastic dynamics from data. This new architecture was theoretically inspired by a universal approximation theorem for stochastic systems, which we introduce in this paper, and the projection-operator formalism for stochastic modeling. Here, we devise mechanisms for training the neural network model to reproduce the correct statistical behavior of a target stochastic process. Numerical simulation results demonstrate that a well-trained SINN can reliably approximate both Markovian and non-Markovian stochastic dynamics. We demonstrate the applicability of SINN to coarse-graining problems and the modeling of transition dynamics. Furthermore, we show that the obtained reduced-order model can be trained on temporally coarse-grained data and hence is well suited for rare-event simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Determining Reduced Order Models for Optimal Stochastic Reduced Order Models

The use of parameterized reduced order models(PROMs) within the stochastic reduced order model (SROM) framework is a logical progression for both methods. In this report, five different parameterized reduced order models are selected and critiqued against the other models along with truth model for the example of the Brake-Reuss beam. The models are: a Taylor series using finite difference, a proper orthogonal decomposition of the the output, a Craig-Bampton representation of the model, a method that uses Hyper-Dual numbers to determine the sensitivities, and a Meta-Model method that uses the Hyper-Dual results and constructs a polynomial curve to better represent the output data. The methods are compared against a parameter sweep and a distribution propagation where the first four statistical moments are used as a comparison. Each method produces very accurate results with the Craig-Bampton reduction having the least accurate results. The models are also compared based on time requirements for the evaluation of each model where the Meta- Model requires the least amount of time for computation by a significant amount. Each of the five models provided accurate results in a reasonable time frame. The determination of which model to use is dependent on the availability of the high-fidelity model and how many evaluations can be performed. Analysis of the output distribution is examined by using a large Monte-Carlo simulation along with a reduced simulation using Latin Hypercube and the stochastic reduced order model sampling technique. Both techniques produced accurate results. The stochastic reduced order modeling technique produced less error when compared to an exhaustive sampling for the majority of methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Stochastic growth logistic model with aftereffect for batch fermentation process

In this paper, the stochastic growth logistic model with aftereffect for the cell growth of C. acetobutylicum P262 and Luedeking-Piret equations for solvent production in batch fermentation system is introduced. The parameters values of the mathematical models are estimated via Levenberg-Marquardt optimization method of non-linear least squares. We apply Milstein scheme for solving the stochastic models numerically. The effciency of mathematical models is measured by comparing the simulated result and the experimental data of the microbial growth and solvent production in batch system. Low values of Root Mean-Square Error (RMSE) of stochastic models with aftereffect indicate good fits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Reservoir Modeling by Data Integration via Intermediate Spaces and Artificial Intelligence Tools in MPS Simulation Frameworks

Conditioning stochastic simulations are very important in many geostatistical applications that call for the introduction of nonlinear and multiple-point data in reservoir modeling. Here, a new methodology is proposed for the incorporation of different data types into multiple-point statistics (MPS) simulation frameworks. Unlike the previous techniques that call for an approximate forward model (filter) for integration of secondary data into geologically constructed models, the proposed approach develops an intermediate space where all the primary and secondary data are easily mapped onto. Definition of the intermediate space, as may be achieved via application of artificial intelligence tools like neural networks and fuzzy inference systems, eliminates the need for using filters as in previous techniques. The applicability of the proposed approach in conditioning MPS simulations to static and geologic data is verified by modeling a real example of discrete fracture networks using conventional well-log data. The training patterns are well reproduced in the realizations, while the model is also consistent with the map of secondary data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2014

Simulations of DSB Yields and Radiation-induced Chromosomal Aberrations in Human Cells Based on the Stochastic Track Structure Induced by HZE Particles

The formation of double‐strand breaks (DSBs) and chromosomal aberrations (CAs) is of great importance in radiation research and, specifically, in space applications. We are presenting a new particle track and DNA damage model, in which the particle stochastic track structure is combined with the random walk (RW) structure of chromosomes in a cell nucleus. The motivation for this effort stems from the fact that the model with the RW chromosomes, NASARTI (NASA radiation track image) previously relied on amorphous track structure, while the stochastic track structure model RITRACKS (Relativistic Ion Tracks) was focused on more microscopic targets than the entire genome. We have combined chromosomes simulated by RWs with stochastic track structure, which uses nanoscopic dose calculations performed with the Monte‐Carlo simulation by RITRACKS in a voxelized space. The new simulations produce the number of DSBs as function of dose and particle fluence for high‐energy particles, including iron, carbon and protons, using voxels of 20 nm dimension. The combined model also calculates yields of radiation‐induced CAs and unrejoined chromosome breaks in normal and repair deficient cells. The joined computational model is calibrated using the relative frequencies and distributions of chromosomal aberrations reported in the literature. The model considers fractionated deposition of energy to approximate dose rates of the space flight environment. The joined model also predicts of the yields and sizes of translocations, dicentrics, rings, and more complex‐type aberrations formed in the G0/G1 cell cycle phase during the first cell division after irradiation. We found that the main advantage of the joined model is our ability to simulate small doses: 0.05‐0.5 Gy. At such low doses, the stochastic track structure proved to be indispensable, as the action of individual delta‐rays becomes more important.

Ponomarev, Artem↗
NASA NTRS2014

Simulations of DSB Yields and Radiation-induced Chromosomal Aberrations in Human Cells Based on the Stochastic Track Structure iIduced by HZE Particles

The formation of double‐strand breaks (DSBs) and chromosomal aberrations (CAs) is of great importance in radiation research and, specifically, in space applications. We are presenting a new particle track and DNA damage model, in which the particle stochastic track structure is combined with the random walk (RW) structure of chromosomes in a cell nucleus. The motivation for this effort stems from the fact that the model with the RW chromosomes, NASARTI (NASA radiation track image) previously relied on amorphous track structure, while the stochastic track structure model RITRACKS (Relativistic Ion Tracks) was focused on more microscopic targets than the entire genome. We have combined chromosomes simulated by RWs with stochastic track structure, which uses nanoscopic dose calculations performed with the Monte‐Carlo simulation by RITRACKS in a voxelized space. The new simulations produce the number of DSBs as function of dose and particle fluence for high‐energy particles, including iron, carbon and protons, using voxels of 20 nm dimension. The combined model also calculates yields of radiation‐induced CAs and unrejoined chromosome breaks in normal and repair deficient cells. The joined computational model is calibrated using the relative frequencies and distributions of chromosomal aberrations reported in the literature. The model considers fractionated deposition of energy to approximate dose rates of the space flight environment. The joined model also predicts of the yields and sizes of translocations, dicentrics, rings, and more complex‐type aberrations formed in the G0/G1 cell cycle phase during the first cell division after irradiation. We found that the main advantage of the joined model is our ability to simulate small doses: 0.05‐0.5 Gy. At such low doses, the stochastic track structure proved to be indispensable, as the action of individual delta‐rays becomes more important.

Ponomarev, Artem↗
DOE OSTI2024

Deep Generative Models for Fast Photon Shower Simulation in ATLAS

The need for large-scale production of highly accurate simulated event samples for the extensive physics programme of the ATLAS experiment at the Large Hadron Collider motivates the development of new simulation techniques. Building on the recent success of deep learning algorithms, variational autoencoders and generative adversarial networks are investigated for modelling the response of the central region of the ATLAS electromagnetic calorimeter to photons of various energies. The properties of synthesised showers are compared with showers from a full detector simulation using GEANT4 . Both variational autoencoders and generative adversarial networks are capable of quickly simulating electromagnetic showers with correct total energies and stochasticity, though the modelling of some shower shape distributions requires more refinement. This feasibility study demonstrates the potential of using such algorithms for ATLAS fast calorimeter simulation in the future and shows a possible way to complement current simulation techniques.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Modeling of Diffusivity for 2D Vacancy Nanopits and Comparison with 2D Adatom Nanoislands on Metal(100) Surfaces Including Analysis for Ag(100)

Diffusion coefficients, D N , for 2D vacancy nanopits are compared with those for 2D homoepitaxial adatom nanoislands on metal(100) surfaces, focusing on the variation of DN with size, N. Here, N is measured in missing atoms for pits and adatoms for islands. Analysis of D N is based on kinetic Monte Carlo simulations of a tailored stochastic lattice-gas model, where pit and island diffusion are mediated by periphery diffusion, i.e., by edge atom hopping. Precise determination of D N versus N for typical parameters reveals a cyclical variation with an overall decrease in magnitude for increasing moderate O(10 2 ) ≤ N ≤ O(10 3 ). Monotonic decay, D N ~ N -β , is found for N ≥ O(10 2 ) with effective exponents, β = β eff , for both pits and islands, both well below the macroscopic value of βmacro = 3/2. D N values for vacancy pits are significantly lower (higher) than for adatom islands for moderate N in the case of low (high) kink rounding barrier. However, D N values for pits and islands slowly merge, and β eff → 3/2 for sufficiently large N. The latter feature is expected from continuum Langevin formulations appropriate for large sizes. Finally, we compare predictions from our model incorporating appropriate energetic parameters for Ag(100) with different sets of experimental data for diffusivity at 300 K, including assessment of β eff , for experimentally observed sizes N from ~100 to ~1000.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Effect of Reaction Kinetics on Graphene-Edge Morphology and Composition

New simulations of graphene growth in flame environments are presented. The simulations employed a kinetic Monte Carlo (KMC) algorithm coupled to molecular mechanics (MM) geometry optimization to track individual graphenic species as they evolve. Focus was given to incorporation of five-member rings and resulting curvature and edge defects. The model code was re-written to be more computationally efficient enabling a larger set of simulations to be run, decreasing stochastic fluctuations in the averaged results. The model also included updated rate coefficients for graphene edge reactions recently published in the literature. Sensitivity analysis, enabled by the more efficient coding, confirmed the importance of the updated rates as well as providing further insight on key reactions controlling layer growth and morphology. The new KMC code was applied to simulate graphene-layer growth in the environments of laminar premixed flames. Finally, in these flame simulations, C-H ratios of evolving structures were computed and compared to those from experiment and from an alternate model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Anisotropic physics-regularized interpretable machine learning of microstructure evolution

Anisotropic Physics-Regularized Interpretable Machine Learning Microstructure Evolution (APRIMME) is a general-purpose machine learning solution for grain growth simulations. In prior work, PRIMME employed a deep neural network to predict site-specific migration as a function of its neighboring sites to model normal, isotropic, grain growth behavior. This work aims to extend this method by incorporating grain boundary misorientation-based grain growth behavior. APRIMME is trained on anisotropic simulations created using the Monte Carlo-Potts (MCP) model. Furthermore, the results of this work are compared statistically using grain radius, number of sides per grain, mean neighborhood misorientations, and the standard deviation of triple junction dihedral angles, and are found to match in most cases. The exceptions are small and seem to be related to two causes: (1) the deterministic model of APRIMME is learning from the stochastic simulations of MCP, which seems to accentuate triple junction behaviors; and, (2) a bias against very small grains is made evident in a quicker decrease in grains than expected at the beginning of an APRIMME simulation. APRIMME is also evaluated for its general ability to capture anisotropic grain growth behavior by first investigating different test case initial conditions, including a circle grain, three grain, and hexagonal grain microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Detection of Historical and Future Precipitation Variations and Extremes Over the Continental United States

Problem: The overall goal of this proposal is to detect observed seasonal-mean precipitation variations and extreme event occurrences over the United States. Detection, e.g. the process of demonstrating that an observed change in climate is unusual, first requires some means of estimating the range of internal variability absent any external drivers. Ideally, the internal variability would be derived from the observations themselves, however generally the observed variability is a confluence of both internal variability and variability in response to external drivers. Further, numerical climate models—the standard tool for detection studies—have their own estimates of intrinsic variability, which may differ substantially from that found in the observed system as well as other model systems. These problems are further compounded for weather and climate extremes, which as singular events are particularly ill-suited for detection studies because of their infrequent occurrence, limited spatial range, and underestimation within global and even regional numerical models. Rationale: As a basis for this research we will show how stochastic daily-precipitation models—models in which the simulated interannual-to-multidecadal precipitation variance is purely the result of the random evolution of daily precipitation events within a given time period—can be used to address many of these issues simultaneously. Through the novel application of these well-established models, we can first estimate the changes/trends in various means and extremes that can occur even with fixed daily-precipitation characteristics, e.g. that can occur simply as a result of the stochastic evolution of daily weather events within a given climate. Detection of a change in the observed climate—either naturally or anthropogenically forced—can then be defined as any change relative to this stochastic variability, e.g. as changes/trends in the means and extremes that could only have occurred through a change in the underlying climate. As such, this method is capable of detecting “hot spot” regions—as well as “flare ups” within the hot spot regions—that have experienced interannual to multi-decadal scale variations and trends in seasonal-mean precipitation and extreme events. Further by applying the same methods to numerical climate models we can discern the fidelity of the current-generation climate models in representing detectability within the observed climate system. In this way, we can objectively determine the utility of these model systems for performing detection studies of historical and future climate change.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

The cost of conservative synchronization in parallel discrete event simulations

The performance of a synchronous conservative parallel discrete-event simulation protocol is analyzed. The class of simulation models considered is oriented around a physical domain and possesses a limited ability to predict future behavior. A stochastic model is used to show that as the volume of simulation activity in the model increases relative to a fixed architecture, the complexity of the average per-event overhead due to synchronization, event list manipulation, lookahead calculations, and processor idle time approach the complexity of the average per-event overhead of a serial simulation. The method is therefore within a constant factor of optimal. The analysis demonstrates that on large problems--those for which parallel processing is ideally suited--there is often enough parallel workload so that processors are not usually idle. The viability of the method is also demonstrated empirically, showing how good performance is achieved on large problems using a thirty-two node Intel iPSC/2 distributed memory multiprocessor.

Nicol, David M.↗
DOE OSTI2012

Simulation Of Microtearing Turbulence In NSTX

Thermal energy confinement times in NSTX dimensionless parameter scans increase with decreasing collisionality. While ion thermal transport is neoclassical, the source of anomalous electron thermal transport in these discharges remains unclear, leading to considerable uncertainty when extrapolating to future ST devices at much lower collisionality. Linear gyrokinetic simulations find microtearing modes to be unstable in high collisionality discharges. First non-linear gyrokinetic simulations of microtearing turbulence in NSTX show they can yield experimental levels of transport. Magnetic flutter is responsible for almost all the transport (~98%), perturbed field line trajectories are globally stochastic, and a test particle stochastic transport model agrees to within 25% of the simulated transport. Most significantly, microtearing transport is predicted to increase with electron collisionality, consistent with the observed NSTX confinement scaling. While this suggests microtearing modes may be the source of electron thermal transport, the predictions are also very sensitive to electron temperature gradient, indicating the scaling of the instability threshold is important. In addition, microtearing turbulence is susceptible to suppression via sheared E-B flows as experimental values of E-B shear (comparable to the linear growth rates) dramatically reduce the transport below experimental values. Refinements in numerical resolution and physics model assumptions are expected to minimize the apparent discrepancy. In cases where the predicted transport is strong, calculations suggest that a proposed polarimetry diagnostic may be sensitive to the magnetic perturbations associated with the unique structure of microtearing turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

System Modeling Techniques for PRA P-200

Course Objectives: Build PRA modeling and analysis skills – Event tree and fault tree model development – Dependent failures and common cause modeling – Component failure mechanisms, Improve understanding of quantification process, Improve ability to extract key results from a PRA, Greater familiarity with PRA goals and process, and Aleatory (stochastic) versus Epistemic (state of knowledge) uncertainty

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Design of experiments to spectroscopically characterize radiation flow in stochastic media

Precise characterization of experimental radiation flow is required to validate the high energy density physics models, numerical methods, and codes that are used to simulate radiation-hydrodynamics phenomena such as thermal radiation transport in stochastic media. The Cassio code is used to simulate thermal radiation flow through inhomogeneous, stochastic-media-foam configurations containing optically thick clumps dispersed within an optically thin background aerogel. Cassio can model small inhomogeneous problems directly, but most problems require approximations to meet computer limitations on run-times and memory usage. Various examples of these approximations are methods that produce, in one calculation, an ensemble-averaged solution and associated standard deviation; reduced spatial dimensionality with approximate geometries; and full material homogenization with no geometric detail. Cassio simulations are used to design experiments at the OMEGA-60 Laser Facility that can measure the radiation flow using the spatially resolved COAX absorption spectroscopy diagnostic. The experimental platforms flow radiation through foam targets ranging from a background-only aerogel, to a single configuration of a specified stochastic medium, to a fully homogenized foam of the background and clump materials. Under constant total clump mass, larger clumps (here, larger than 10 μm diameter) will mix more slowly with the background such that the bulk radiation flow is faster than it would be in a fully homogenized material. The COAX platform can be used to infer temperature and density profiles in both the background material and clumps, simultaneously, and therefore to differentiate radiation flow in a range of stochastic and homogeneous media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

A stochastic model for particle impingements on orbiting spacecraft

A general methodology for simulating particle impingements on orbiting spacecraft is developed. Major steps in the modeling process are presented as (1) modeling objective, (2) construction of the spacecraft geometrical model, (3) simulation of the particles in the space environment, (4) particle impact and subsequent events of interest, and (5) results of the simulation. A simulation of the expected meteoroid impingements on the Hubble Space Telescope and the resulting angular momentum transfers which can cause telescope pointing disturbances is given to illustrate these methods.

Howell, L. W., Jr.↗