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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A Subgrid Approach for Modeling Microtopography Effects on Overland Flow

Microtopography, or heterogeneities in the elevation across scales much smaller than the domain of interest, plays a critical role in surface water retention, surface/subsurface interactions, and runoff. Resolving microtopographic influences on flow requires high–resolution simulations, which are computationally demanding even when considering the surface system in isolation and even more so when surface flow is one component in integrated simulations that couple surface flow with unsaturated subsurface flow. There is thus significant motivation for models that allow the effects of subgrid microtopography to be better represented. Subgrid models modify coarsened models to capture the microtopography–induced nonlinear effects on hydrologic processes. We present a subgrid model that alters the water storage and flow terms in the diffusion wave equation for surface flow. Stochastically generated microtopography with strongly contrasting spatial structure, high–resolution digital elevation maps from a polygonal tundra site on the North Slope of Alaska and a hummocky microtopography from a field site in Northern Minnesota are used to assess the accuracy and applicability of the subgrid model to disparate landscapes. Approaches for determining subgrid model parameters are tested and simulation results using the subgrid model are compared to benchmark fine–scale simulations and to coarse simulations that ignore microtopography. Our findings confirm that a properly parameterized subgrid model greatly improves the coarse–scale representation of hydrographs and total water content in the system. Furthermore, using the polygonal tundra example, we propose and test a strategy for moving to application–relevant spatial scales by combining microtopography classification and a few fine–scale simulations on small subdomains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Ten questions concerning occupant behavior in buildings: The big picture

Occupant behavior has significant impacts on building energy performance and occupant comfort. However, occupant behavior is not well understood and is often oversimplified in the building life cycle, due to its stochastic, diverse, complex, and interdisciplinary nature. The use of simplified methods or tools to quantify the impacts of occupant behavior in building performance simulations significantly contributes to performance gaps between simulated models and actual building energy consumption. Therefore, it is crucial to understand occupant behavior in a comprehensive way, integrating qualitative approaches and data- and model-driven quantitative approaches, and employing appropriate tools to guide the design and operation of low-energy residential and commercial buildings that integrate technological and human dimensions. This paper presents ten questions, highlighting some of the most important issues regarding concepts, applications, and methodologies in occupant behavior research. The proposed questions and answers aim to provide insights into occupant behavior for current and future researchers, designers, and policy makers, and most importantly, to inspire innovative research and applications to increase energy efficiency and reduce energy use in buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2016

Solving the master equation without kinetic Monte Carlo: Tensor train approximations for a CO oxidation model

In multiscale modeling of heterogeneous catalytic processes, one crucial point is the solution of a Markovian master equation describing the stochastic reaction kinetics. Usually, this is too high-dimensional to be solved with standard numerical techniques and one has to rely on sampling approaches based on the kinetic Monte Carlo method. In this study we break the curse of dimensionality for the direct solution of the Markovian master equation by exploiting the Tensor Train Format for this purpose. The performance of the approach is demonstrated on a first principles based, reduced model for the CO oxidation on the RuO{sub 2}(110) surface. We investigate the complexity for increasing system size and for various reaction conditions. The advantage over the stochastic simulation approach is illustrated by a problem with increased stiffness.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Multi-fidelity Gaussian process regression for prediction of random fields

We propose a new multi-fidelity Gaussian process regression (GPR) approach for prediction of random fields based on observations of surrogate models or hierarchies of surrogate models. Our method builds upon recent work on recursive Bayesian techniques, in particular recursive co-kriging, and extends it to vector-valued fields and various types of covariances, including separable and non-separable ones. The framework we propose is general and can be used to perform uncertainty propagation and quantification in model-based simulations, multi-fidelity data fusion, and surrogate-based optimization. We demonstrate the effectiveness of the proposed recursive GPR techniques through various examples. Specifically, we study the stochastic Burgers equation and the stochastic Oberbeck–Boussinesq equations describing natural convection within a square enclosure. In both cases we find that the standard deviation of the Gaussian predictors as well as the absolute errors relative to benchmark stochastic solutions are very small, suggesting that the proposed multi-fidelity GPR approaches can yield highly accurate results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Swelling and Diffusion during Methanol Sorption into Hydrated Nafion

Diffusion within polymer electrolyte membranes is often coincident with time-dependent processes such as swelling and polymer relaxation, which are factors that limit their ability to block molecular crossover during use. The solution-diffusion model of membrane permeation, which is the accepted theory for dense polymers, applies only to steady-state processes and does not address dynamic internal structural changes that can accompany permeation. To begin discovery of how such changes can be coupled to the permeation process, we have constructed here a stochastic multiscale reaction-diffusion model that examines time-dependent methanol uptake into and swelling of hydrated Nafion. Several potential mechanisms of diffusion and polymer response are tested. The simulation predictions are compared to real-time Fourier transform infrared attenuated total reflectance spectroscopy (FTIR-ATR) absorbance reported in the literature [ Hallinan , D. T. , Jr. ; Elabd , Y. A. J. Phys. Chem. B 2007 , 111 , 13221 - 13230 ]. Of the proposed polymer response mechanisms, only one, a reaction-limited, local response to increasing methanol concentration that takes the entire experimental time frame of 600 s, produces simulated FTIR-ATR data consistent with experiment. The simulations show that water diffusion out of the membrane is minimal during methanol sorption and that changes in the measured infrared absorbances are due primarily to the increase in methanol concentration accompanied by dilution of water during swelling. Swelling involves densification of the polymer structure even as there is an overall volume expansion of the film. Potential connections between the polymer densification and molecular-level structural changes of Nafion in methanol are discussed. These results indicate that the interaction between methanol and Nafion serves to increase Nafion's capacity to accommodate large volumes of methanol-water solutions, facilitating increased permeation across the membrane relative to pure water.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of diffusion constant on the morphology of dendrite growth in lithium metal batteries

Lithium dendrites can lead to a short circuit and battery failure, and developing strategies for their suppression is of considerable importance. In this work, we study the growth of dendrites in a simple model system where the solvent is a continuum and the lithium ions are hard spheres that can deposit by sticking to existing spheres or the electrode surface. Using stochastic dynamics simulations, we investigate the effect of applied voltage and diffusion constant on the growth of dendrites. We find that the diffusion constant is the most significant factor, and the inhomogeneity of the electric field does not play a significant role. Additionally, the growth is most pronounced when the applied voltage and diffusion constant are both low. We observe a structural change from broccoli to cauliflower shape as the diffusion constant is increased. The simulations suggest that a control of electrolyte parameters that impact lithium diffusion might be an attractive route to controlling dendrite growth.

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

Stochastic Parameterization: Toward a New View of Weather and Climate Models

The last decade has seen the success of stochastic parameterizations in short-term, medium-range, and seasonal forecasts: operational weather centers now routinely use stochastic parameterization schemes to represent model inadequacy better and to improve the quantification of forecast uncertainty. Developed initially for numerical weather prediction, the inclusion of stochastic parameterizations not only provides better estimates of uncertainty, but it is also extremely promising for reducing long-standing climate biases and is relevant for determining the climate response to external forcing. This article highlights recent developments from different research groups that show that the stochastic representation of unresolved processes in the atmosphere, oceans, land surface, and cryosphere of comprehensive weather and climate models 1) gives rise to more reliable probabilistic forecasts of weather and climate and 2) reduces systematic model bias. We make a case that the use of mathematically stringent methods for the derivation of stochastic dynamic equations will lead to substantial improvements in our ability to accurately simulate weather and climate at all scales. Recent work in mathematics, statistical mechanics, and turbulence is reviewed; its relevance for the climate problem is demonstrated; and future research directions are outlined

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Fast Computational Algorithms for Partial Differential Equations and Uncertainty Quantifications

This project concerned the construction, testing and analysis of computational algorithms for solving parameterized and stochastic partial differential equations. The study and understanding of equations of this type is of fundamental importance in numerous engineering and scientific applications. Examples include simulation of plasma dynamics in models of electric propulsion and nuclear fusion, simulation of multiphase flows, such as the flow of water, gas and oil in reservoirs, and structural analysis of the dependence of structures on materials. Parametrization is used in such settings when properties of the models such as viscosity of fluids or electric resistivity of materials are not precisely understood and instead are treated as random variables. The resulting solutions are themselves random, and having such solutions will enable engineers to use probabilistic methods to assess the likelihood of events, for example, whether a pollutant in a liquid will exceed a limit, and to use such analyses to develop ways to ensure positive outcomes. Construction of accurate (high resolution) computational solutions is expensive, requiring significant computer time and computational resources, and there is need to reduce computational cost to make simulation useful and effective. The aim of the project was to construct algorithms to efficiently compute surrogate solutions to parameterized problems to allow for efficient and accurate simulation. The technical approach used focused on two related strategies, based on rank-reduction methods and reduced-order models. These methods construct surrogate solutions of parameter-dependent models by projection or interpolation into low-dimensional approximation spaces. Cost savings are achieved if the low-dimensional spaces can be identified and constructed efficiently and if the resulting low-dimensional algebraic systems can be solved cheaply. Accomplishments include: Theoretical and empirical demonstration of the effectiveness of fast multigrid solution strategies for computing low-rank representations of parameter-dependent solutions to discrete partial differential equations, including the first proof establishing so-called textbook convergence properties for low-rank methods. Development of efficient solution algorithms for solving nonlinear parameter-dependent partial differential equations used in models of fluid dynamics. Developent of efficient algorithms for low-rank representation of solutions of time-dependent simulations of fluid dynamics using multi-dimensional tensor representations of solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Towards Unifying the Planetary Boundary Layer and Shallow Convection in CAM5 with the Eddy-Diffusivity/Mass-Flux Approach

The modular structure of the boundary layer and convection parameterizations in atmospheric models have long been affecting the numerical representation of subgrid-scale motions and their mutual interactions. A promising alternative, the eddy-diffusivity/mass-flux approach (EDMF), has the potential for unifying the existing formulations into a consistent scheme and improving some of the long-standing issues. This study documents a step towards developing such a unified approach by implementing a stochastic multi-plume EDMF scheme into the Community Atmosphere Model (Version 5.0). Its performance in single-column mode is evaluated against the control parameterization and large-eddy simulation (LES) for two benchmark cases: marine and continental shallow convection. Overall, the results for the two parameterizations agree well with each other and with LES in terms of mean profiles of moist conserved variables and their vertical fluxes, as well as the updraft properties. However, systematic differences between the two schemes, especially for transient continental convection, are also documented. Using EDMF helps improve some of the parameterized features of shallow convection. In particular, for the highest tested vertical resolution, the EDMF cloud base and top heights and the vertical fluxes of energy and water are remarkably close to LES.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2012

Signatures of the Rayleigh-Plateau Instability Revealed by Imposing Synthetic Perturbations on Nanometer Liquid Metals on Substrates

Fluid jets destabilize into droplets according to the well known Rayleigh Plateau (RP) instability. The applicability of the RP model to a nanoscale liquid metal supported on a substrate (a so called rivulet) was investigated using molecular dynamics (MD) simulations as it has important implications for nanoscale materials assembly. Stochastic fluctuations were found to affect rivulet stability and dynamics when compared with microscopic predictions of rivulet break up; here a broad and delayed distribution of break up times was observed. Yet, hallmark characteristics of the RP theory were observed; (1) a critical wavelength demarcating stability, (2) stable and unstable growth modes and (3) neck formation induced by axial pressure gradients. Notably, MD simulations implementing synthetic perturbations, yielded significantly reduced dispersion in the final nanodroplet spacing and clearly revealed the underlying physics driving rivulet break up.

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

Application of Optimal Production Control theory for Home Energy Management in a Micro Grid

We consider the optimal stochastic control problem for home energy systems with solar and energy storage devices when the demand is realized from the grid. The demand is subject to Brownian motions with both drift and variance parameters modulated by a continuous-time Markov chain that represents the regime of electricity price. We model the systems as pure stochastic differential equation models, and then we follow the completing square technique to solve the stochastic home energy management problem. The effectiveness of the efficiency of the proposed approach is validated through a simulation example. For practical situations with constraints consistent to those studied here, our results imply the proposed framework could reduce the electricity cost from short-term purchase in peak hour market.

Malikopoulos, Andreas↗
DOE OSTI2015

Optimization of a dry, mixed nuclear fuel storage array for nuclear criticality safety - 14339

A dry storage array of used nuclear fuel at the Idaho National Laboratory contains a mixture of more than twenty different types of research and test reactor fuels. Significant neutronic interaction occurs between fuel canisters within this large storage arrangement of 636 fuel storage canisters. Nuclear criticality safety is maintained primarily by moderation control and through control of neutron interaction and leakage provided by the credited geometry. Over the nearly forty years of operation the collected fuels are now arranged in a configuration that does not minimize the system neutron multiplication factor (k{sub eff}). This creates a concern since the entire remaining facility storage capacity cannot be utilized without exceeding the system minimum subcritical state (k{sub safe}) for ensuring nuclear criticality safety. This work determines a more optimal arrangement for the stored fuels with a goal to minimize the system k{sub eff}, but also with a minimum of potential fuel canister relocation movements. The ultimate goal is to establish a solution for this multiple-objective optimization problem by finding the Pareto Front for maximum k{sub eff} reduction versus the number of canister relocations. The process first required additional analysis of the array by Monte-Carlo neutron transport codes (specifically, MCNP) to gain a greater understanding of the interaction characteristics of the disparate fuel types. Some stored fuels are composed of primarily metallic and oxide components that provide very little inherent moderation, while other fuel types contain graphite or zirconium-hydride which can contribute a significant amount of moderation. These factors cause a complex relationship within a fuel arrangement with respect to thermal utilization and axial leakage probability. The new MCNP analysis was utilized to develop a simplified spreadsheet-based empirical model which simulates the array k{sub eff} result. This empirical model was benchmarked against additional MCNP calculations for various arrangements and then employed in the optimization analysis. A Tabu metaheuristic and stochastic approximation optimization was applied to the empirical model with the optimization function to minimize the array k{sub eff}. The optimization algorithm interrogated the empirical model to determine the optimal rearrangement for up to sixty canister movements. The Pareto Front was found, which shows that a limited number of canister relocation movements could allow the existing thirty-one empty array positions to be filled with common research reactor fuels while also yielding a reduction in the overall array k{sub eff} relative to the current arrangement. This will impart an improvement in the facility utilization while also offering an enhancement in the safety margin. The reduction in array reactivity has been confirmed by a final set of MCNP calculations with the proposed new arrangements. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

A Continuum Dislocation Dynamics Crystal Plasticity Approach to Irradiated Body-Centered Cubic α-Iron

Radiation-induced embrittlement of reactor pressure vessel (RPV) steels can potentially limit the operating life of nuclear power plants. Over extended exposure to radiation doses, these body-centered cubic (BCC) irons demonstrate irradiation damage. Here, we present a continuum dislocation density (CDD) crystal plasticity model to capture the interaction among dislocations and self-interstitial atom (SIA) loops in α-iron. We demonstrate the importance of modeling cross slip using a combined stochastic Monte Carlo approach and the role of slip system strength anisotropy in capturing stochastic cross slip interactions. Through these captured interactions, the CDD crystal plasticity model can capture both the stress response and the physical evolution of dislocations on different slip system planes. Single-crystal verification experiments are used to calibrate the CDD crystal plasticity model, and a set of simplified polycrystalline simulations demonstrates the model’s ability to capture the stress response from tensile experiments on α-iron.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Crystalline islands of semiconductor films

Silicon carbide (SiC) nucleation in the form of powder in a discharge plasma and the formation of thin film islands on a Si(100) substrate in the course of gas-phase epitaxy are simulated numerically. Models of plasma-like media and nonequilibrium processes accompanying phase transitions of the first kind (such as condensation and crystallization) in the initial fast (fluctuation) stage are described. The nonstationary evolution of nuclei size distribution functions is modeled by solving kinetic equations in partial derivatives and stochastic Ito-Stratonovich analog equations. This makes it possible to refine the formation mechanisms of microcrystalline state polytypes and calculate the nucleation rate and the initial roughness of a SiC coating.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Random Process Simulation for stochastic fatigue analysis

A simulation technique is described which directly synthesizes the extrema of a random process and is more efficient than the Gaussian simulation method. Such a technique is particularly useful in stochastic fatigue analysis because the required stress range moment E(R sup m), is a function only of the extrema of the random stress process. The family of autoregressive moving average (ARMA) models is reviewed and an autoregressive model is presented for modeling the extrema of any random process which has a unimodal power spectral density (psd). The proposed autoregressive technique is found to produce rainflow stress range moments which compare favorably with those computed by the Gaussian technique and to average 11.7 times faster than the Gaussian technique. The autoregressive technique is also adapted for processes having bimodal psd's. The adaptation involves using two autoregressive processes to simulate the extrema due to each mode and the superposition of these two extrema sequences. The proposed autoregressive superposition technique is 9 to 13 times faster than the Gaussian technique and produces comparable values for E(R sup m) for bimodal psd's having the frequency of one mode at least 2.5 times that of the other mode.

Larsen, Curtis E.↗
DOE OSTI2018

The Use of Hydrogen as a Fuel for Engines in the Energy Cycle of Remote Production Facilities

The approach to using hydrogen as fuel, which ensures the smooth operation of autonomous power systems that use renewable energy sources (wind or solar power installations) with the stochastic mode of power generation, has been presented. The fundamental possibility of implementing the nondetonation combustion of hydrogen via the addition of ecologically clean components or a small percentage of methane has been demonstrated by methods of mathematical modeling.

09 BIOMASS FUELS↗
DOE OSTI2023

High-fidelity multiphysics load following and accidental transient modeling of microreactors using NEAMS tools: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The feasibility of modeling microreactors using high-fidelity models with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tools is investigated in this report. Three overarching questions guided this research: can NEAMS tools readily be applied for high-fidelity multiphysics modeling of different types of transients in microreactor designs; how accurate are the results obtained; and are improvements needed in accuracy or user experience of NEAMS tools, especially considering newly developed capabilities? This work builds upon FY-2022 work, and two microreactor concepts considering heat pipe (HP-MR) and gas-cooled (GC-MR) technologies were further analyzed using high-fidelity multiphysics simulations. The NEAMS tools considered and coupled within the MultiApp environment are Griffin for neutronics, BISON for thermo-mechanics, Sockeye for heat pipe modeling (in HP-MR), SAM for 1D Fluid – 3D solid modeling of coolant channels and system modeling of balance of plant components (in GC-MR), and the SWIFT code for hydrogen redistribution in hydride moderator. The Heat Pipe MicroReactor (HP-MR) concept was further analyzed in FY-2023 to demonstrate the stochastic TRISO failure modeling capability in BISON to check operational limits of the TRISO fuel. A new full-core Gas-Cooled MicroReactor (GC-MR) model was developed based on the initial assembly-model used in Y-2022 and used for steady-state and accidental depressurization transient simulations. Accuracy of the simulations performed was assessed through 1) verification analyses completed on the different physics with code-to-code comparison, and 2) validation of the multiphysics simulations based on modeling of the Kilopower Reactor Using Stirling Technology (KRUSTY) experiment. In FY-2023, the mesh and model of KRUSTY was updated to closely match publicly available data, and the neutronic model was verified and validated against experimental control rod worth measurements. The multiphysics model of KRUSTY was developed and used for steady-state analysis and for modeling reactivity insertion transient. The calculated power increase and stabilization agrees well with experimental data following adjustment in fuel thermal expansion coefficient. As an important component of this project, the ANL team gathered experience with a wide range of NEAMS tools: the MOOSE Mesh System, Griffin, BISON, SWIFT, Sockeye, SAM, Workbench, and the MOOSE MultiApp System, and provided assessment of new capabilities. Noteworthy are the user assessment of the “vapor-only” flow model in Sockeye and development of a multiphysics startup transient in HP-MR unit cell for use as tutorial in Sockeye. The full-core GC-MR model was used for assessment of SAM for balance of plant modeling and for demonstrating the SWIFT code capability for hydrogen redistribution modeling in multiphysics transient analyses. In this process, several bugs/issues were identified and reported to developers. Finally, the assembly GC-MR model developed in FY-2022 coupling Griffin, BISON and SAM through flow blockage and rod ejection transients was published to the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB). The Heat Pipe MicroReactor (HP-MR) concept high-fidelity multiphysics coupling of Griffin/BISON/Sockeye in load-following and heat pipe failure transients was also published on the VTB. Those submissions are enabling thorough review of these models as well as wide distribution to industry, regulator, and university users. In this analysis, several new research questions were uncovered, and follow-up analyses are recommended to further improve some models, consider additional transients, and continue development of VTB models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗