Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo integration”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Comparisons of Integrated Radiation Transport Models with Microdosimetry Data in Spaceflight

Astronauts are exposed to galactic cosmic rays (GCR), trapped protons, and possible solar particle events (SPE) during spaceflight. For such complicated mixtures of radiation types and kinetic energies, tissue equivalent proportional counters (TEPC's) represent a simple time-dependent approach for radiation monitoring. Of interest in radiation protection is the average quality factor of a radiation field defined as a function of linear energy transfer, LET, Q(sub ave)(LET). However TEPC's measure the average quality factors as a function of lineal energy (y), Q(sub ave)(y) defined as the average energy deposition in a volume divided by the average chord length of the volume. Lineal energy, y deviates from LET due to energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector. Using integrated space radiation models that includes the transport code HZETRN/BRYNTRN, the quantum nuclear interaction model, QMSFRG, and results from Monte-Carlo track simulations of TEPC's response to ions, we consider comparisons of model calculations to TEPC results from NASA missions in low Earth orbit and make predictions for lunar and Mars missions. Good agreement between the model and measured spectra from past NASA missions is found. A finding of this work is that TEPC's values for trapped or solar protons of Q(sub ave)(y) range from 1.9-2.5, overestimating Q(sub ave)(LET), which ranges from 1.4-1.6 with both quantities increasing with shielding depth due to nuclear secondaries Comparisons for the complete GCR spectra show that Q(sub ave)(LET) for GCR is approximately 3.5-4.5, while TEPC's measure 2.9-3.4 for Q(sub ave)(y) with the GCR values decreasing with depth as heavy ions are absorbed in shielding material. Our results support the use of TEPC's for space radiation environmental monitoring when computational analysis is used for proper data interpretation.

Cucinotta, Francis A.↗
DOE OSTI2023

Multi-fidelity microstructure-induced uncertainty quantification by advanced Monte Carlo methods

Quantifying uncertainty associated with the microstructure variation of a material can be a computationally daunting task, especially when dealing with advanced constitutive models and fine mesh resolutions in the crystal plasticity finite element method (CPFEM). Numerous studies have been conducted regarding the sensitivity of material properties and performance to the mesh resolution and choice of constitutive model. However, a unified approach that accounts for various fidelity parameters, such as mesh resolutions, integration time-steps and constitutive models simultaneously is currently lacking. Further, this paper proposes a novel uncertainty quantification (UQ) approach for computing the properties and performance of homogenized materials using CPFEM, that exploits a hierarchy of approximations with different levels of fidelity. In particular, we illustrate how multi-level sampling methods, such as multi-level Monte Carlo (MLMC) and multi-index Monte Carlo (MIMC), can be applied to assess the impact of variations in the microstructure of polycrystalline materials on the predictions of homogenized materials properties. We show that by adaptively exploiting the fidelity hierarchy, we can significantly reduce the number of microstructures required to reach a certain prescribed accuracy. Finally, we show how our approach can be extended to a multi-fidelity framework, where we allow the underlying constitutive model to be chosen from either a phenomenological plasticity model or a dislocation-density-based model.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, in this work, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Status of the OECD/NEA Watts Bar unit 1 benchmark and calculations of local reactor core power data of the Zero Power physics tests

The paper aims to provide an update on the status of the Organization for Economic Cooperation and Development (OECD) / the Nuclear Energy Agency (NEA) Tennessee Valley Authority (TVA) Watts Bar 1 (WB1) Multi-Physics Multi-Cycle depletion benchmark and the results associated with Exercise 1. The benchmark relies on the set of benchmark progression problems developed by the Department of Energy (DOE) Consortium for the Advanced Simulation of Light Water Reactors (CASL) for the Virtual Environment for Reactor (VERA), which are based on real plant design and operational data. The OECD/NEA TVA WB1 benchmark is designed for validation of both traditional and novel high-fidelity multi-physics codes to analyze Pressurized Water Reactors (PWR) depletion cycles. The activities are conducted under the Expert Group on Reactor Systems Multi-Physics (EGMUP) at NEA/OECD. In this work, we analyze Exercise 1 of the benchmark: stand-alone Three-Dimensional (3D) neutronics at Start-up Zero Power Physics Test (ZPPT) at Hot Zero Power Conditions (HZP) and the resulting power maps. The code used to model the exercise is the continuous energy Monte Carlo code Serpent 2.1.31 along with the nuclear data library ENDF/B-VII.1. Serpent results are compared with the high-fidelity deterministic code MPACT, which is part of VERA. The paper presents three selected results from power calculations required output. The results compared are the normalized axially integrated radial core power maps, the normalized axial averaged core power shapes, and the normalized core hottest and coldest assemblies radial power maps. The Root Mean Square Deviation (RMSD) between Serpent and MPACT is 0.51 % for the axially integrated radial core power map, 1.41 % for the axial averaged core power shape, 2.42 % for the hottest assembly axial power shape, 0.88 % for the coldest assembly axial integrated hottest assembly radial pin power map, and 0.17 % for the axially integrated coldest assembly radial pin power map. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Preliminary verification of the MCNP perturbation and fixed-source tally sensitivity tools

Integral benchmark experiments are vital in the adjustment and validation of the nuclear data that govern predictive simulations across the nuclear community. The nuclear data sensitivity capabilities of the Monte Carlo N-Particle (MCNP ®) transport code are currently limited; however, expanding sensitivity capabilities will allow benchmark experiments to be designed to resolve compensating errors and adjust nuclear data where previously prohibitively difficult. This paper provides details of a preliminary verification for the use of (i.) the recently revised perturbation and (ii.) developmental fixed-source sensitivity tools within MCNP to calculate sensitivities of tallied responses (such as current integrated over a surface, F1, and flux averaged over a cell, F4) to nuclear data in fixed-source simulations. Energy-binned and energy-integrated sensitivities calculated with these tools are compared against sensitivities calculated using a central-difference approximation. Here, the verification is completed for four configurations of a benchmarked system using a 4.5-kg plutonium sphere surrounded by varying amounts of copper and/or polyethylene. The results show that sensitivities calculated with the perturbation and fixed-source sensitivity tools agree with the central-difference-based approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Hybrid Parameter Search and Dynamic Model Selection for Mixed-Variable Bayesian Optimization

Herein this article presents a new type of hybrid model for Bayesian optimization (BO) adept at managing mixed variables, encompassing both quantitative (continuous and integer) and qualitative (categorical) types. Our proposed new hybrid models (named hybridM) merge the Monte Carlo Tree Search structure (MCTS) for categorical variables with Gaussian Processes (GP) for continuous ones. hybridM leverages the upper confidence bound tree search (UCTS) for MCTS strategy, showcasing the tree architecture’s integration into Bayesian optimization. Our innovations, including dynamic online kernel selection in the surrogate modeling phase and a unique UCTS search strategy, position our hybrid models as an advancement in mixed-variable surrogate models. Numerical experiments underscore the superiority of hybrid models, highlighting their potential in Bayesian optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Preliminary design of control rods in the single-fluid double-zone thorium molten salt reactor (SD-TMSR)

Recent studies on Molten Salt Reactors (MSRs) showed that the excess reactivity at the beginning of the operation is large for many fueling strategies and must be compensated by a reactivity control system. The current work introduces a reliable safety system based on control rods in addition to the online feed system reactivity control in the Single-fluid Double-zone Thorium-based Molten Salt Reactor (SD-TMSR). Three different initial fissile loadings are considered: 233 U, reactor-grade Pu, and transuranic (TRU) elements as a startup fuel. We applied six different absorbing materials to investigate the main operational and safety parameters in the SD-TMSR: natural B 4 C, enriched B 4 C with 90% 10 B, HfB 2 , HfH 1.62 , Eu 2 O 3 , and Gd 2 O 3 . The present work focuses on control rod design, integral and differential control rod worth, shutdown margin, and shadowing effects at steady-state. We employed the SERPENT-2 Monte-Carlo code to calculate the reactivity worth and analyze the performance of the reactivity control system. We showed that 233 U and reactor-grade Pu startup cores maintain adequate shutdown margin with all considered absorbers. Finally, this report proposes a design of control rod clusters that compensate the excess reactivity of the SD-TMSR loaded with different initial fissile material.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

A dissolution model of alite coupling surface topography and ions transport under different hydrodynamics conditions at microscale

Portland cement is the most produced material in the world. The hydration process of cement consists of a group of complex chemical reactions. In order to investigate the mechanism of cement hydration, it is vital to study the hydration of each phase separately. An integrated model is proposed in this paper to simulate the dissolution of alite under different hydrodynamic conditions at microscale, coupling Kinetic Monte Carlo model (KMC), Lattice Boltzmann method (LBM) and diffusion boundary layer (DBL). The dissolution of alite is initialised with KMC. Two Multiple-relaxation-time (MRT) LB models are used to simulate the fluid flow and transport of ions, respectively. For solid-liquid interface, DBL is adapted to calculate the concentration gradient and dissolution flux. The model is validated with experiment from literature. The simulation results show good agreements with the results published in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2020

BioSTEAM-LCA: An Integrated Modeling Framework for Agile Life Cycle Assessment of Biorefineries under Uncertainty

Biorefineries will play a critical role in sustainable bioeconomies, but projections of their environmental impacts vary widely. A core challenge with life cycle assessments (LCAs) of biorefineries is that they are often disconnected from biorefinery design, simulation, and techno-economic analysis (TEA). This lack of integration is a barrier to early stage technology and process evaluations, reducing consistency and transparency across sustainability indicators while limiting our understanding of the relative importance of individual factors (e.g., design decisions, greenhouse gas emission accounting procedures), how these factors interact, and trade-offs or synergies with process economics. In this study, we propose a new agile LCA framework, BioSTEAMLCA, which layers onto BioSTEAM (Biorefinery Simulation and Techno-Economic Analysis Modules, which automates biorefinery design, simulation, and TEA) to characterize the environmental impacts of biorefineries across a landscape of designs, technology performance assumptions, and contexts. Inventory databases and impact assessment methods are integrated to enable flexible user defined LCA system models, and the implications of uncertainties throughout the production system are characterized via Monte Carlo simulation. To demonstrate the capabilities of BioSTEAMLCA, we present a case study for sugarcane ethanol production. Overall, BioSTEAM-LCA enables computationally efficient, agile gate-to-gate LCA to evaluate biorefinery processes, the production of candidate biofuels and bioproducts, and trade-offs among productivity, economics, and environmental impacts under uncertainty.

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

CoVVVR: Control Variates & Vegas Variance Reduction

This package is a wrapper over the vegas integration package. The control variate variance reduction method is applied to the function when integrated along with the techniques applied in vegas such as importance sampling. To understand control variates, lets first look at how Monte Carlo works.

Scott, JacobL. [Univ. of Kansas, Lawrence, KS (Uni↗
DOE OSTI2023

Bayesian Cloud Property Retrievals from ARM Active and Passive Measurements

The optimum use of the continuous measurements of thermodynamics, radiation, aerosols, clouds and precipitation from the DOE Atmospheric Radiation Measurement (ARM) program is key to achieve the DOE Atmospheric System Research (ASR)’s objectives. One of the key mission requirements is to retrieve cloud and precipitation properties, as well as vertical motion parameters, along the vertical cross- section defined by the profiling active sensors. Such retrievals are challenging to perform continuously in the entire spectrum of cloud and precipitation conditions due to the large natural microphysical and dynamical variability, the often-limited information content in the measurements, and the lack of proper characterization of measurement quality and uncertainty. Today, the acquisition of new remote and in-situ sensors by the ARM program creates opportunities to address the microphysical retrieval problem by exploiting new, more robust retrieval techniques and integrating various scattered advancements in both sensor techniques and retrieval algorithms. During this project, we constructed a robust Bayesian Markov chain Monte Carlo (MCMC) cloud property retrieval algorithm that includes a state of the art radar forward model. Our MCMC-based retrieval produces both the best estimate of height-resolved cloud and precipitation properties in the radar profile, as well as an estimate of the in-cloud vertical motion and turbulence. In addition, the MCMC algorithm automatically produces robust and flexible estimates of retrieval uncertainty. We tested the algorithm on several synthetic cloud profiles obtained from large eddy simulation (LES) models with bin-resolved microphysics.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Space Radiation Transport Methods Development

Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the progress towards a full three-dimensional and computationally efficient deterministic code for which the current HZETRN evaluates the lowest order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard Finite Element Method (FEM) geometry common to engineering design practice enabling development of integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 milliseconds and severely limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of reconfigurable computing and could be utilized in the final design as verification of the deterministic method optimized design.

Wilson, J. W.↗
NASA NTRS2004

A space radiation transport method development

Improved spacecraft shield design requires early entry of radiation constraints into the design process to maximize performance and minimize costs. As a result, we have been investigating high-speed computational procedures to allow shield analysis from the preliminary design concepts to the final design. In particular, we will discuss the progress towards a full three-dimensional and computationally efficient deterministic code for which the current HZETRN evaluates the lowest-order asymptotic term. HZETRN is the first deterministic solution to the Boltzmann equation allowing field mapping within the International Space Station (ISS) in tens of minutes using standard finite element method (FEM) geometry common to engineering design practice enabling development of integrated multidisciplinary design optimization methods. A single ray trace in ISS FEM geometry requires 14 ms and severely limits application of Monte Carlo methods to such engineering models. A potential means of improving the Monte Carlo efficiency in coupling to spacecraft geometry is given in terms of re-configurable computing and could be utilized in the final design as verification of the deterministic method optimized design. Published by Elsevier Ltd on behalf of COSPAR.

NASA Center JSC↗
DOE OSTI2020

Constrained Bayesian Optimization of Criticality Experiments at LLNL [Slides]

The design of criticality experiments is typically an iterative process that employs a Monte Carlo transport code. The goal is to find a design that optimizes some variable, like the sensitivity of a response to a cross section, while simultaneously ensuring criticality. The high fidelity of the Monte Carlo code is a great asset, but it makes exploring the design space computationally expensive. Herein, we present how a constrained Bayesian optimization algorithm can be used to efficiently design a criticality experiment. It uses Gaussian processes as a surrogate model to probe the design space and to reduce the number of code executions that are needed to find the optimum. We demonstrate constrained Bayesian optimization with a Pu-239/polyethylene solution system and a TEX experiment that is designed for critical ity safety validation of a nuclear waste model at the Hanford Site. For both systems, a global optimum was found within 75 Monte Carlo simulations.

42 ENGINEERING↗
DOE OSTI2021

Measurement of the production cross section of pairs of isolated photons in pp collisions at 13 TeV with the ATLAS detector

A measurement of prompt photon-pair production in proton-proton collisions at $\sqrt{s}$ = 13 TeV is presented. The data were recorded by the ATLAS detector at the LHC with an integrated luminosity of 139 fb -1 . Events with two photons in the well-instrumented region of the detector are selected. The photons are required to be isolated and have a transverse momentum of p T, γ1(2) > 40 (30) GeV for the leading (sub-leading) photon. The differential cross sections as functions of several observables for the diphoton system are measured and compared with theoretical predictions from state-of-the-art Monte Carlo and fixed-order calculations. The QCD predictions from next-to-next-to-leading-order calculations and multi-leg merged calculations are able to describe the measured integrated and differential cross sections within uncertainties, whereas lower-order calculations show significant deviations, demonstrating that higher-order perturbative QCD corrections are crucial for this process. The resummed predictions with parton showers additionally provide an excellent description of the low transverse-momentum regime of the diphoton system.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2013

An Approach for the Assessment of System Upset Resilience

This report describes an approach for the assessment of upset resilience that is applicable to systems in general, including safety-critical, real-time systems. For this work, resilience is defined as the ability to preserve and restore service availability and integrity under stated conditions of configuration, functional inputs and environmental conditions. To enable a quantitative approach, we define novel system service degradation metrics and propose a new mathematical definition of resilience. These behavioral-level metrics are based on the fundamental service classification criteria of correctness, detectability, symmetry and persistence. This approach consists of a Monte-Carlo-based stimulus injection experiment, on a physical implementation or an error-propagation model of a system, to generate a system response set that can be characterized in terms of dimensional error metrics and integrated to form an overall measure of resilience. We expect this approach to be helpful in gaining insight into the error containment and repair capabilities of systems for a wide range of conditions.

Torres-Pomales, Wilfredo↗
NASA NTRS2014

Optimizing Medical Kits for Spaceflight

The Integrated Medical Model (IMM) is a probabilistic model that estimates medical event occurrences and mission outcomes for different mission profiles. IMM simulation outcomes describing the impact of medical events on the mission may be used to optimize the allocation of resources in medical kits. Efficient allocation of medical resources, subject to certain mass and volume constraints, is crucial to ensuring the best outcomes of in-flight medical events. We implement a new approach to this medical kit optimization problem. METHODS We frame medical kit optimization as a modified knapsack problem and implement an algorithm utilizing a dynamic programming technique. Using this algorithm, optimized medical kits were generated for 3 different mission scenarios with the goal of minimizing the probability of evacuation and maximizing the Crew Health Index (CHI) for each mission subject to mass and volume constraints. Simulation outcomes using these kits were also compared to outcomes using kits optimized..RESULTS The optimized medical kits generated by the algorithm described here resulted in predicted mission outcomes more closely approached the unlimited-resource scenario for Crew Health Index (CHI) than the implementation in under all optimization priorities. Furthermore, the approach described here improves upon in reducing evacuation when the optimization priority is minimizing the probability of evacuation. CONCLUSIONS This algorithm provides an efficient, effective means to objectively allocate medical resources for spaceflight missions using the Integrated Medical Model.

Optimization↗
DOE OSTI2024

Control And Optimization Modular Modeling Application For Nuclear Deployment

The purpose of the COMMAND code is to provide a flexible, scalable tool for use in developing, integrating, and testing the technologies necessary for achieving autonomous operations of advanced nuclear reactors. The code enables users to efficiently implement custom simulations and experiments by combining key methods from different software modules. These modules are focused on: modeling and simulation tools, such as nuclear simulation tools used for high-fidelity modeling (e.g., Reactor Excursion and Leak Analysis Program [RELAP5-3D] and Monte Carlo N-Particle [MCNP]); machine learning and optimization tools (e.g., anomaly detection and data-driven modeling techniques); advanced control in its digital, high-performance, and supervisory control forms (e.g., proportional integral derivative (PID) control and model predictive control (MPC); and integration with hardware through industrial communication protocols. To ensure flexibility and scalability, COMMAND was designed to be both modular—the software “pieces” all inherit from generic building blocks and can be combined and connected to create complicated simulations—and high performing—designed for parallel processing, enabling simulations and experiments to take advantage of multi-core computers, servers, and nodes. The code is written in the Python programming language due to the language's popularity, active community, and open-source and cross-platform nature. Maintaining consistency with other simulation tools used within the nuclear energy community, users implement simulations and experiments through text input files, which define components, parameters, connections, etc., through lines of text. Given that COMMAND is written in Python, these input files are native Python scripts, and so use the standard Python structure and formatting. This also enables users to take advantage of Python's extensive package library to develop custom capabilities for their specific use cases.

Faber, Jacob [Idaho National Laboratory (INL), Ida↗