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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Progress on Demonstration of a MOOSE-Based Coupled Capability for Hot Channel Factors in Fast Reactors

Hot channel factors (HCFs) are computed values that account for the impact on predicted peak fuel, cladding, and coolant temperatures due to uncertainties in the as-built reactor’s material properties and geometry as well as uncertainties due to modeling approximations. Reduction in computed HCF values via reduction or elimination of modeling approximations may translate to significant economic savings if the reactor power can be raised due to the extra temperature margin gained. While limited historical datasets exist for sodium-cooled fast reactors (SFRs), there are no available HCF data for lead-cooled fast reactors (LFRs) outside of work generated previously within NEAMS. The computation of HCFs involves insights from reactor physics, thermal fluids and heat conduction calculations to determine how the peak temperatures respond to various uncertainties in the design. Due to the significant advantages for multi-physics coupling offered by the MOOSE framework, Griffin (MOOSE-based reactor physics code), MOOSE Heat Conduction Module, and Cardinal (MOOSE-wrapped multi-physics application which includes the NekRS thermal fluids code) are being coupled together using the MOOSE MultiApp System to develop a highfidelity multi-physics modeling capability for HCF simulations. This high-fidelity coupling workflow may also be beneficial for other fast reactor applications in the future. In previous work, Griffin and NekRS were individually assessed to ensure the necessary capabilities were in place. This work describes initial efforts to couple the codes (including folding in the MOOSE Heat Conduction Module) and determining the workflow for the perturbed calculations which will leverage the Stochastic Tools Module (STM). To our knowledge, this is the first coupling of Griffin and NekRS as well as the first exploratory use of Stochastic Tools Module for Cardinal. In this report, the neutronics code Griffin, the heat conduction solver in MOOSE, and the MOOSE-wrapped application containing NekRS (Cardinal) are linked together to demonstrate the coupled capability. Griffin and Cardinal are linked dynamically by specifying shared libraries. Different coupling hierarchies are tested for selecting the most appropriate coupling strategy. A coupling scheme is selected based on the efficiency of calculation and ease of data communication. Multiple tests are performed to choose suitable mesh structure, model configurations, scheme setup and boundary conditions to avoid loss of energy due to data interpolation between different modules or weak imposition of fluxes in finite element codes. Computational experiments are performed to study the tolerance control of each type of iteration to avoid false convergence. The coupled capability is demonstrated in both single pin and 7-pin models based on LFR materials and geometry. The study finds that the use of too large a time step size in the heat conduction module can lead to temperature oscillation even though the heat conduction equation does not have a time-derivative kernel, but only the time-dependent boundary condition. A 7-pin model without duct region achieved good convergence in the coupled calculation while a 7 pin model with duct region experienced data communication issues which need to be resolved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2001

Free Vibration of Uncertain Unsymmetrically Laminated Beams

Monte Carlo Simulation and Stochastic FEA are used to predict randomness in the free vibration response of thin unsymmetrically laminated beams. For the present study, it is assumed that randomness in the response is only caused by uncertainties in the ply orientations. The ply orientations may become random or uncertain during the manufacturing process. A new 16-dof beam element, based on the first-order shear deformation beam theory, is used to study the stochastic nature of the natural frequencies. Using variational principles, the element stiffness matrix and mass matrix are obtained through analytical integration. Using a random sequence a large data set is generated, containing possible random ply-orientations. This data is assumed to be symmetric. The stochastic-based finite element model for free vibrations predicts the relation between the randomness in fundamental natural frequencies and the randomness in ply-orientation. The sensitivity derivatives are calculated numerically through an exact formulation. The squared fundamental natural frequencies are expressed in terms of deterministic and probabilistic quantities, allowing to determine how sensitive they are to variations in ply angles. The predicted mean-valued fundamental natural frequency squared and the variance of the present model are in good agreement with Monte Carlo Simulation. Results, also, show that variations between plus or minus 5 degrees in ply-angles can affect free vibration response of unsymmetrically and symmetrically laminated beams.

Kapania, Rakesh K.↗
DOE OSTI2019

Pilot Point Optimization of Mining Boundaries for Lateritic Metal Deposits: Finding the Trade-off Between Dilution and Ore Loss

Geological contacts in lateritic metal deposits (footwall topographies) often delineate the orebody boundaries. Spatial variations seen in such contacts are frequently higher than those for the metal grades of the deposit. Therefore, borehole spacing chosen based on the grade variations cannot adequately capture the geological contact variability. Consequently, models created using such boreholes cause high volumetric uncertainties in the actual and targeted ore extraction volumes, which, in turn, lead to high unplanned dilution and ore losses. In this paper, a method to design optimum ore/mining boundaries for lateritic metal deposits is presented. The proposed approach minimizes the dilution/ore losses and comprises two main steps. First, the uncertainty on the orebody boundary is represented using a set of stochastic realizations generated with a multiple-point statistics algorithm. Then, the optimal orebody boundary is determined using an optimization technique inspired by a model calibration method called Pilot Points. The pilot points represent synthetic elevation values, and they are used to construct smooth mining boundaries using the multilevel B-spline technique. The performance of a generated surface is evaluated using the expected sum of losses in each of the stochastic realizations. The simulated annealing algorithm is used to iteratively determine the pilot point values which minimize the expected losses. The results show a significant reduction in the dilution volume as compared to those obtained from the actual mining operation.

58 GEOSCIENCES↗
DOE OSTI2018

Analytic Analysis for Dynamic System Frequency in Power Systems Under Uncertain Variability

In this paper, a new systematic method is proposed to analyze the system frequency in power systems under uncertain variability (UV). UV gradually increases in power systems, which originates from the renewable energy generation, random loads, frequency measurement, and communication channels, exerting noteworthy impacts on dynamic system frequency. The UV is modeled as a stochastic process in this paper. Then, the stochastic differential equations are used to describe the dynamic system frequency response (SFR) of a power system under UV, which is based on the SFR model. To assess system frequency dynamics under UV, an index of intra-range probability is put forward. Based on stochastic analysis theory, an analytic method is proposed to analyze the system frequency under UV, by which the intra-range probability can be analytically solved. Compared with Monte Carlo simulation in a typical SFR case and the Iceland electricity transmission network, the proposed method shows almost the same results using much less computation resource. Finally, insights for improving the SFR under UV are provided. The proposed method can be used to quickly verify whether the system frequency could withstand the UV and stay in the security range.

42 ENGINEERING↗
DOE OSTI2012

Spectral analysis of fundamental signal and noise performances in photoconductors for mammography

Purpose: This study investigates the fundamental signal and noise performance limitations imposed by the stochastic nature of x-ray interactions in selected photoconductor materials, such as Si, a-Se, CdZnTe, HgI{sub 2}, PbI{sub 2}, PbO, and TlBr, for x-ray spectra typically used in mammography. Methods: It is shown how Monte Carlo simulations can be combined with a cascaded model to determine the absorbed energy distribution for each combination of photoconductor and x-ray spectrum. The model is used to determine the quantum efficiency, mean energy absorption per interaction, Swank noise factor, secondary quantum noise, and zero-frequency detective quantum efficiency (DQE). Results: The quantum efficiency of materials with higher atomic number and density demonstrates a larger dependence on convertor thickness than those with lower atomic number and density with the exception of a-Se. The mean deposited energy increases with increasing average energy of the incident x-ray spectrum. HgI{sub 2}, PbI{sub 2}, and CdZnTe demonstrate the largest increase in deposited energy with increasing mass loading and a-Se and Si the smallest. The best DQE performances are achieved with PbO and TlBr. For mass loading greater than 100 mg cm{sup -2}, a-Se, HgI{sub 2}, and PbI{sub 2} provide similar DQE values to PbO and TlBr. Conclusions: The quantum absorption efficiency, average deposited energy per interacting x-ray, Swank noise factor, and detective quantum efficiency are tabulated by means of graphs which may help with the design and selection of materials for photoconductor-based mammography detectors. Neglecting the electrical characteristics of photoconductor materials and taking into account only x-ray interactions, it is concluded that PbO shows the strongest signal-to-noise ratio performance of the materials investigated in this study.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Influence of grain charge gradients on the dynamics of macroparticles in an electrostatic trap

An analytical model of anomalous heating of charged dust grains (macroparticles) caused by their stochastic motion in a bounded plasma volume is proposed. Analytical expressions allowing one to describe the pumping (heating) of interacting grains with additional stochastic energy due to grain charge gradients are derived. The analytical results are verified by numerical simulation of the problem. It is shown that spatial variations in the charges of dust grains can lead to their anomalous heating in laboratory plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Data From: Influence of Hydrological Perturbations and Riverbed Sediment Characteristics on Hyporheic Zone Respiration of CO2 and N-2, Journal of Geophysical Research-Biogeosciences

This data package contains pumping data (.txt), parameter matrices, and R code (.R, .RData) to perform bootstrapping for parameter selection for the bioclogging model development. The pumping data were collected from the Russian River Riverbank Filtration site located in Sonoma County, California from 2010-2017 from three riverbank collection wells located alongside the study site. The pumping data is directly correlated with water table oscillations, so the code performs these correlations and simulates stochastic versions of water table oscillations. See Metadata Description.pdf for full details on dataset production. This dataset must be used with the R programming language. This dataset and R code is associated with the publication "Influence of Hydrological Perturbations and Riverbed Sediment Characteristics on Hyporheic Zone Respiration of CO2 and N-2"This research was supported by the Jane Lewis Fellowship from the University of California, Berkeley, the Sonoma County Water Agency (SCWA), the Roy G. Post Foundation Scholarship, the U.S. Department of Energy, Office of Science Graduate Student Research (SCGSR) Program, U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research under award DE-AC02-05CH11231, and the UFZ-Helmholtz Centre for Environmental Research, Leipzig, Germany.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1991

A Scheduling Algorithm for Replicated Real-Time Tasks

We present an algorithm for scheduling real-time periodic tasks on a multiprocessor system under fault-tolerant requirement. Our approach incorporates both the redundancy and masking technique and the imprecise computation model. Since the tasks in hard real-time systems have stringent timing constraints, the redundancy and masking technique are more appropriate than the rollback techniques which usually require extra time for error recovery. The imprecise computation model provides flexible functionality by trading off the quality of the result produced by a task with the amount of processing time required to produce it. It therefore permits the performance of a real-time system to degrade gracefully. We evaluate the algorithm by stochastic analysis and Monte Carlo simulations. The results show that the algorithm is resilient under hardware failures.

Yu, Albert C.↗
NASA NTRS2017

The Mechanism for the Energy Buildup Driving Solar Eruptive Events

The underlying origin of solar eruptive events (SEEs), ranging from giant coronal mass ejections to small coronalhole jets, is that the lowest-lying magnetic flux in the Sun's corona undergoes continual buildup of stress and free energy. This magnetic stress has long been observed as the phenomenon of "filament channels:" strongly sheared magnetic field localized around photospheric polarity inversion lines. However, the mechanism for the stress buildup-formation of filament channels-is still debated. We present magnetohydrodynamic simulations of a coronal volume that is driven by transient, cellular boundary flows designed to model the processes by which the photosphere drives the corona. The key feature of our simulations is that they accurately preserve magnetic helicity, the topological quantity that is conserved even in the presence of ubiquitous magnetic reconnection. Although small-scale random stress is injected everywhere at the photosphere, driving stochastic reconnection throughout the corona, the net result of the magnetic evolution is a coherent shearing of the lowest-lying field lines. This highly counterintuitive result-magnetic stress builds up locally rather than spreading out to attain a minimum energy state-explains the formation of filament channels and is the fundamental mechanism underlying SEEs. Furthermore, this process is likely to be relevant to other astrophysical and laboratory plasmas.

Knizhnik, K. J.↗
DOE OSTI2015

A Multiobjective Optimization Framework for Online Stochastic Optimal Control in Hybrid Electric Vehicles

The increasing urgency to extract additional efficiency from hybrid propulsion systems has led to the development of advanced power management control algorithms. In this paper we address the problem of online optimization of the supervisory power management control in parallel hybrid electric vehicles (HEVs). We model HEV operation as a controlled Markov chain and we show that the control policy yielding the Pareto optimal solution minimizes online the long-run expected average cost per unit time criterion. The effectiveness of the proposed solution is validated through simulation and compared to the solution derived with dynamic programming using the average cost criterion. Both solutions achieved the same cumulative fuel consumption demonstrating that the online Pareto control policy is an optimal control policy.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Letter: Modeling reactive shock waves in heterogeneous solids at the continuum level with stochastic differential equations

In this paper, a new paradigm is introduced for modeling reactive shock waves in heterogeneous solids at the continuum level. Inspired by the probability density function methods from turbulent reactive flows, it is hypothesized that the unreacted material microstructures lead to a distribution of heat release rates from chemical reaction. Fluctuations in heat release, rather than velocity, are coupled to the reactive Euler equations which are then solved via the Riemann problem. Finally, a numerically efficient, one-dimensional hydrocode is used to demonstrate this new approach, and simulation results of a representative impact calculation (inert flyer into explosive target) are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Enabling site-specific well leakage risk estimation during geologic carbon sequestration using a modular deep-learning-based wellbore leakage model

Geologic carbon sequestration (GCS) is a promising technology for mitigating net carbon emissions and growing climate concern by storing CO 2 in reservoirs. Oil and gas brownfields are an attractive option for CO 2 storage, but these sites have many historical wellbores from petroleum production and can be a potential leakage pathway for CO 2 or formation brine. Therefore, risk management of GCS operations requires an assessment of potential well leakage. Due to the high uncertainty of the system, stochastic approaches are ideal for quantifying the range of risk behaviors, but they must be computationally efficient in the face of complex physics. Here, we develop a new physics-centric deep learning wellbore model to predict the leakage of CO 2 and brine through leaky wellbores. Multi-physics numerical simulations were used to generate data sets, and physics-informed features were introduced. Neural networks were optimized with an automated searching algorithm. Feature analysis quantifies the impact of each feature on model prediction and confirms the role of physics-inspired parameters. The model shows high predictive performance across a wide range of geologic and injection conditions and well attributes. In conclusion, a case study illustrates how the model is applied to assess well leakage in GCS operations.

58 GEOSCIENCES↗
DOE OSTI2024

Improving Robustness of Spectrogram Classifiers with Neural Stochastic Differential Equations

Signal analysis and classification is fraught with high levels of noise and perturbation. Computer-vision-based deep learning models applied to spectrograms have proven useful in the field of signal classification and detection; however, these methods aren't designed to handle the low signal-to-noise ratios inherent within non-vision signal processing tasks. While they are powerful, they are currently not the method of choice in the inherently noisy and dynamic critical infrastructure domain, such as smart-grid sensing, anomaly detection, and non-intrusive load monitoring. Currently, these models can be brittle, which makes them susceptible to noisy input. This also means they have sub-optimal stability of explanation outputs. Experts and technicians using these models to make decisions in real world scenarios need assurance that a model is performing as it is supposed to. The classification or prediction outputs it generates should be sound and grounded, not likely to change in the presence of shifting noise landscapes. In this work, we explore the idea of Neural Stochastic Differential Equations (NSDE's) to improve the robustness of models trained to classify time series data and the effect of NSDE's on the explainability of outputs. We then test the effectiveness of these approaches by applying them to a non-intrusive load monitoring (NILM) dataset that consists of simulated harmonic signals injected into a real building.

Brogan, Joel↗
DOE OSTI2016

Towards Next Generation Lithium-Sulfur Batteries: Non-Conventional Carbon Compartments/Sulfur Electrodes and Multi-Scale Analysis

In this work, a novel heterofunctional, bimodal-porous carbon morphology, termed the carbon compartment (CC), is utilized as a sulfur host as a lithium-sulfur battery cathode. A multi-scale model explores the physics and chemistry of the lithium-sulfur battery cathode. The CCs are synthesized by a rapid, low cost process to improve electrode-electrolyte interfacial contact and accommodate volumetric expansion associated with sulfide formation. The CCs demonstrate high sulfur loading (47 %-wt. S) and ca. 700 mAh g -1 reversible capacity with high coulombic efficiency due to their unique structures. Density functional theory and ab initio Molecular Dynamics characterize the interface between the C/S composite and electrolyte during the sulfur reduction mechanism. Stochastic realizations of 3D electrode microstructures are reconstructed based on representative SEM images to study the influence of solid sulfur loading and lithium sulfide precipitation on microstructural and electrochemical properties. A macroscale electrochemical performance model is developed to analyze the performance of lithium-sulfur batteries. The combined multi-scale simulation studies explain key fundamentals of sulfur reduction and its relation to the polysulfide shuttle mechanism: how the process is affected due to the presence of carbon substrate, thermodynamics of lithium sulfide formation and deposition on carbon, and microstructural effects on the overall cell performance.

25 ENERGY STORAGE↗
NASA NTRS1987

Stochastic acceleration in the transrelativistic region and pion production in solar flares

The stochastic Fermi acceleration spectrum in the transrelativistic region obtained from a Monte-Carlo simulation for an energy-independent alpha(T) is much harder than the extension of the nonrelativistic analytic spectrum to this energy range for the same alpha(T). The latter, with alpha(T) = 0.043, was used to model the pion and nuclear line emissions for the impulsive phase of the 3 Jun. 1982 flare, as well as the 2.223 MeV emission from this flare. We find that the ratios of these three emissions for the Monte-Carlo spectrum with alpha(T) = 0.028 are essentially the same as those for the analytical spectrum with alpha(T) = 0.043. We also find that the acceleration time from approximately 30 MeV to approximately 1 GeV is less than or approximately = 10 s, consistent with the observations of the 3 Jun. 1982 flare.

Miller, James A.↗
NASA NTRS2019

The Role of Magnetic Helicity in Coronal Heating

One of the greatest challenges in solar physics is understanding the heating of the Sun's corona. Most theories for coronal heating postulate that free energy in the form of magnetic twist/stress is injected by the photosphere into the corona where the free energy is converted into heat either through reconnection or wave dissipation. The magnetic helicity associated with the twist/stress, however, is expected to be conserved and appear in the corona. In previous works, we showed that the helicity associated with the small-scale twists undergoes an inverse cascade via stochastic reconnection in the corona and ends up as the observed large-scale shear of filament channels. Our "helicity condensation" model accounts for both the formation of filament channels and the observed smooth, laminar structure of coronal loops. In this paper, we demonstrate, using helicity- and energy-conserving numerical simulations of a coronal system driven by photospheric motions, that the model also provides a natural mechanism for heating the corona. We show that the heat generated by the reconnection responsible for the helicity condensation process is sufficient to account for the observed coronal heating. We study the role that helicity injection plays in determining coronal heating and find that, crucially, the heating rate is only weakly dependent on the net helicity preference of the photospheric driving. Our calculations demonstrate that motions with 100% helicity preference are least efficient at heating the corona; those with 0% preference are most efficient. We discuss the physical origins of this result and its implications for the observed corona.

Knizhnik, K. J.↗
DOE OSTI2017

Flux in a 1D-Reactor Cell by Response Matrix

Legacy computational reactor physics software tools and protocols currently used in support of the Advanced Test Reactor (ATR) core fuel management and safety assurance, and to some extent, experiment management, are inconsistent with the state of modern nuclear engineering practice. Thus, they are difficult, if not impossible, to verify and validate (V and V) according to modern standards. Furthermore, the legacy staff knowledge required for effective application of their application and protocols from the 1960's and 1970's is rapidly being lost due to staff turnover and retirements. In late 2009, the Idaho National Laboratory (INL) initiated a focused effort, the ATR Core Modeling Update Project, to address this situation through the introduction of modern high-fidelity transport theory based computational software and models to replace the legacy methods, with appropriate (V and V) protocols according to applicable national standards. The new computational tools include: - The 2-dimensional HELIOS deterministic transport code for integrated ATR core fuel and experiment management and Core Safety Analysis Package (CSAP) physics data generation. The MCNP5 3-dimensional stochastic simulation (Monte Carlo) code for verification and validation applications and for specialized ATR core analysis along with its longstanding traditional role in experiment support - The 2-dimensional deterministic NEWT transport code and the 3-dimensional KENO6 stochastic simulation code, with their SCALE-based support infrastructure for use in V and V applications and as supporting modules for sensitivity-uncertainty analysis via the TSUNAMI system. The ATTILA 3-dimensional discrete-ordinates transport code for specialized experiment support and possibly additional applications. Regardless of application, every neutron transport code, new or legacy, requires verification. The aim of this presentation is to provide a neutron transport benchmark for a simple 1D-reactor cell in plane geometry. The benchmark finds use as a standard of comparison for the codes listed above intended for the ATR upgrade. The method developed here is in contrast to the finite Fourier transform solution found in Ref. 1, which is limited to isotropic scattering and particularly difficult to numerically implement unless the total cross section is constant over the cell.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Development of random geometry capability in RMC code for explicit analysis of polytype particle/pebble - 18493

Monte Carlo method plays an important role in accurate simulation of High Temperature Gas-Cooled Reactors (HTR) and molten salt reactors, owing to its advantages of the flexible geometry modeling and the use of continuous-energy nuclear cross sections. For the innovative design of pebble bed, prism HTRs and molten salt reactors, there may be more than one kind of particles or pebbles in the system, such as the burnable poison particles, graphite pebbles or fuel pebbles with different enrichments. Different treatments have been used to deal with the double-heterogeneous geometry of these reactors. In some researches, the particle distributions in stochastic media are assumed the regular lattice, while in the others, the explicit packing methods such as Random Sequential Addition (RSA) and Discrete Element Method (DEM) are used, by generating the random coordinates of the particles/pebbles. However, traditional geometry with repeated structure has difficulties in explicit modeling polytype particle/pebble, due to the arbitrary packing fractions of different kinds of particles or pebbles. Although, packing methods such as Random Sequential Addition (RSA) and Discrete Element Method (DEM) could solve the polytype particle/pebble problems to some extent, but these methods are troublesome in geometry modeling and time consuming. In this paper, a new Random Universe Geometry method was developed and implemented in Monte Carlo code RMC, for explicit analysis of polytype particle/pebble in double heterogeneous geometry systems, such as HTRs and molten salt reactors. This method was applied to constructed mixed particles/pebbles problems and also simplified HTR-10 full core problem to test the feasibility, and the results show that the random geometry modeling can provide explicit packing fraction of different kinds of particles or pebbles. The results are also compared with that calculated by RSA method and regular structure, and agree quite well. This work would provide an effective tool for the innovative design of pebble bed, prism HTRs and molten salt reactors with polytype particles or pebbles using Monte Carlo method. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗