Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo code”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Hexagonal Geometries in MPACT

The MPACT code is a high-fidelity light-water reactor analysis code using whole-core pin-resolved neutron transport calculations on modern parallel-computing hardware. MPACT uses the 2D/1D method to solve 3D neutron transport problems by decomposing the problem into a stack of 2D slices, each of which is solved independently using the method of characteristics (MOC). The slices are then coupled axially using the P3 nodal expansion method (NEM-P3) for the 1D axial calculations. MPACT also employs the coarse mesh finite difference (CMFD) method to accelerate calculations. This manuscript details work supporting advanced reactor designs using hexagonal pins and hexagonal assemblies such as the VVER-1000. If performed correctly, MOC is geometry agnostic. However, MPACT previously had optimizations in place for Cartesian geometries, specifically in the modularization and current calculations. Sections 2 and 3 detail the changes made to MPACT to support MOC and CMFD calculations on hexagonal geometries. Section 4 reports results demonstrating solution consistency for problems run with and without CMFD acceleration, results demonstrating solution consistency when run in serial and parallel, and pincell results using the Monte Carlo code, McCard’s benchmark results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Radiological Characterisation of VLLW Produced at High-Energy Particle Accelerators - 19118

In this paper we describe the characterisation method of very-low-level radioactive waste (VLLW) produced at the European Organization for Nuclear Research (CERN). The method is based on the calculation of a preliminary radionuclide inventory using either analytical or Monte Carlo codes, on the measurement of γ-emitters and on the quantification of pure-α, pure-β and low-energy X-ray emitters, called difficult-to-measure (DTM) radionuclides, using the so-called scaling factor (SF) technique. The first stage of the characterisation process is the calculation of the radionuclide inventory, which is a list of radionuclides that can potentially be produced via the interaction of primary and secondary particles with the surrounding media. Once the preliminary radionuclide inventory is obtained, the γ-emitting radionuclides are measured via γ-ray spectrometry on each package of the waste population. The major γemitter, called key nuclide (KN), is also identified. The scaling factor method estimates the activity of DTM radionuclides by checking for a consistent and repeated relationship between the key nuclide and the activity of the difficult to measure radionuclides from samples. Between 2016 and 2018, several hundred samples were collected at CERN and over 1000 radiochemical analyses were performed to estimate robust statistics to quantify scaling factors for the most common DTM radionuclides, such as H-3, Fe-55 and Ni-63. This paper illustrates the sample collection strategy, the tests performed to choose the most appropriate statistical models describing the behaviour of scaling factors and a complete quantification of the uncertainties associated with the characterisation process. Finally, we present the results obtained during the characterisation of over 2000 m{sup 3} radioactive metals and cables disposed of at the French National repository for VLLW. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

Development of Critical Experiments to Provide Validation Data for Multiphysics Coupling Methods (Final Report)

Computational simulation has become one of the most important tools in the modern toolbox for nuclear reactor design. Reactor simulations, like the reactors themselves, involve multiple physical processes including neutronics, fluid flow and heat transfer (thermal-hydraulics), thermo-mechanics, fuel behavior, chemistry, and balance of plant. Each process can affect the behavior of the others, and advanced computational tools have been developed that couple separate physics models together to simulate the feedback effects. The objective of this work is to generate experimental data to support the validation of multi-physics reactor codes. It specifically targets the SHARP toolset from the NEAMS Reactor Product Line, but it is intended to be useful in any set of coded that couples thermal-hydraulics with neutronics simulations. The project takes advantage of the inherent flexibility of the RPI Reactor Critical Facility (RCF), a low power open-pool research reactor. The standard configuration of the RCF is used to explore the effects of moderator temperature on neutronic behavior and reactivity. Because the low power of the reactor generates no appreciable temperature change in the moderator, the moderator temperature is artificially changed with two 18 kW electric immersion heaters. The temperature feedback effects can be observed by recording the change in the reactor power as the reactor goes from a supercritical state, to critical, to a subcritical state. To expand the temperature range and temperature rate of change that can be investigated, an experimental apparatus consisting of a heated water loop was designed and installed in the RCF. In this setup, water from a heated reservoir tank is pumped through a pipe placed through the center of the reactor core. The water in the reservoir can be heated to around 70 °C and pumped through the system. The pump has a variable flow rate and electronic control, so a number of scenarios can be developed. A slug of water can be pushed through the system to generate a reasonably fast change in reactivity. The pump speed can be increased or decreased over the course of a measurement, or even oscillated at different periods, to observe the response of the reactor. To complement the experimental data generated by the RCF, a set of computational models have been developed. They are intended to illustrate the usefulness of the experimental data generated in validating coupled multi-physics simulations. A model previously constructed in MCNP has been used in conjunction with safety analysis at the RCF for many years and was modified to support the design of the heated water loop experiment. This model was translated into Serpent 2, a similar Monte Carlo code, which provided the added advantage of being able to produce collapsed cross-section data files for use in the deterministic neutron transport codes that are part of the SHARP toolset. A complete reconstruction of the RCF model, both with and without the heated fluid loop experiment, were necessary for use with the SHARP codes, PROTEUS (neutronics) and Nek5000 (thermal hydraulics), as these require an unstructured mesh format to operate. The computational requirements for fully coupled three dimensional simulations are significant, and challenges to full implementation and comparison with RCF experiments remain. A wide variety of experiments have been conducted for eventual benchmarking and validation of multi-physics reactor codes. A complete database of detailed time- and temperature-dependent measurements have been generated and stored. Several examples of the experimental results are presented here, with refinement and analysis of experiments ongoing. While these experiments still only address a finite range of operating conditions, in particular being limited to water in the sub-boiling regime, they represent a unique (to date) look at the feedback between temperature and neutronic behavior suited for validation of the coupling routines in multi-physics codes. It is hoped that these will be the foundation of future benchmarks that will provide confidence in the tools that will be used to develop the next generation of nuclear reactor systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

The importance of delayed neutron precursors in gamma dose calculations for activated primary heat exchanger components in the Molten Salt Breeder Reactor

The Molten Salt Reactor (MSR) Multiphysics Applications technical area in the Nuclear Energy Advanced Modeling and Simulation program has supported the development of 3D Monte Carlo models of the Molten Salt Breeder Reactor (MSBR) over the last couple of years. This MSBR model was previously run with the Shift Monte Carlo code to perform radiation shielding calculations in the reactor cell area. The MSBR is a 2250 MWth (1000 MWe) liquid-fueled molten salt reactor design developed at Oak Ridge National Laboratory in the 1970s. Determining the source terms from activated primary heat exchanger (HX) components is important because delayed and prompt neutron fluxes incident on these components affect the dose rate in the primary HX maintenance areas. This information can be used in the development of remote handling procedures required during shutdown for maintenance. A methodology has been developed and is proposed in this paper to quantify the activated source term from the primary HX components as a result of the movement of the delayed neutron precursors in flowing primary fuel salt through the primary HXs in the MSBR. The goal of this research is to evaluate the gamma dose rates in the maintenance hatches above the primary HX using the activated HX source terms. The study showed that the gamma dose rates are approximately two orders of magnitudes higher when accounting for the neutron activation from the movement of delayed neutron precursors through the HXs than when flowing fuel is not considered. Thus, the movement of delayed neutron precursors must be taken into account for accurately predicting the neutron activation of primary loop components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Baryon Anticorrelations in PYTHIA

Abstract We present a computational investigation of a problem of hadron collisions from recent years, that of baryon anticorrelations. This is an experimental dearth of baryons near other baryons in phase space, not seen upon examining numerical Monte Carlo simulations. We have addressed one of the best known Monte Carlo codes, , to see what baryon (anti)correlations it produces, where they are originated at the string-fragmentation level in the underlying Lund model, and what simple modifications could lead to better agreement with data. We propose two ad-hoc alterations of the fragmentation code, a “one-baryon” and an “always-baryon” policies that qualitatively reproduce the data behaviour, i.e anticorrelation, and suggest that lacking Pauli-principle induced corrections at the quark level could be the culprit behind the current disagreement between computations and experiment.

Physics↗
DOE OSTI2018

Validation of statistical uncertainties in subcritical benchmark measurements: Part II – Measured data

Part I of this work presented an uncertainty approach which incorporates the singles and doubles counting rates determined by the Hage-Cifarelli formalism of the Feynman Variance-to-Mean method. This moments-based approach utilizes time correlations between prompt neutrons detections to assess the system multiplication. Additionally, it presented a validation which utilized simulated data generated using a 0-D point-kinetics Monte Carlo code. In Part II, this same method and validation approach is applied to measured data of a subcritical benchmark experiment. These measurements were performed with a 4.5 kg sphere of alpha-phase weapons grade plutonium reflected by copper and/or polyethylene. Finally, the results of 17 configurations are shown including detailed validation for four of the configurations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Validation of Statistical Uncertainties in Subcritical Benchmark Measurements: Part I - Theory and Simulations

Subcritical neutron noise measurements are often used to determine multiplication of a system with Special Nuclear Material. This work presents an uncertainty approach which incorporates the singles and doubles counting rates determined by the Hage-Cifarelli formalism. This is a moments-based approach which utilizes time correlations between prompt neutrons in order to assess the system multiplication. After the method is described, a validation is presented which utilizes simulated data generated using a 0-D point-kinetics Monte Carlo code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Utilization of ACE nuclear data file toolkit ACEtk to calculate relative sensitivity coefficients of point-kinetics parameters

Sensitivity and uncertainty methods are quintessential for nuclear criticality safety and experiment design. This type of analysis relies on calculations of sensitivity coefficients; sensitivity coefficients of the effective neutron multiplication factor with respect to some nuclear data are predominantly calculated and used. As a part of the Laboratory Directed Research & Development project EUCLID (Experiments Underpinned by Computational Learning for Improvements in nuclear Data) at Los Alamos National Laboratory, sensitivity coefficients of many radiation detector measurement responses with respect to nuclear data were investigated. Specifically, this paper outlines a method to calculate point-kinetics parameters relative sensitivity coefficients with respect to nuclear data. Point-kinetics parameters such as the prompt neutron decay constant, effective delayed neutron fraction, and neutron generation time are especially important to experimenters and reactor operators designing systems with dynamic neutron populations. This method couples capabilities of the ACE (A Compact ENDF) nuclear data file toolkit, ACEtk, with the ability to load cross sections into the radiation transport code Monte Carlo N-Particle (MCNP). In conclusion, key aspects of optimizing this method for a particular application and sensitivity profiles of the Jezebel criticality experiment are examined and discussed.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2020

Accelerated version of NUBEAM capabilities in DIII-D using neural networks

A neural network model of the effects of neutral beam injection on DIII-D has been developed. The training and testing data used by the model have been generated by the NUBEAM module of TRANSP for experimental discharges from the 2018 DIII-D campaign. Using a principle component analysis to reduce the dimensionality of profile data, the model has been shown to reproduce the results of the Monte Carlo code NUBEAM with a high level of accuracy and an execution time orders of magnitude faster than the execution time of NUBEAM. Furthermore, this makes the neural network model uniquely suited to applications in model-based scenario planning (off-line) and active control (on-line), where a large number of simulation runs are required by the associated optimization tasks that need to be performed before and during the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Analysis of fast-ion losses measured in MAST-U via infrared thermography and a Fast Ion Loss Detector

Fast-ion losses need to be monitored to avoid damage to plasma facing components. In existing experimental devices, the scintillator-based fast-ion loss detector (FILD) is the most advanced diagnostic for measuring fast-ion losses. However, FILDs provide only local information about the losses. Infrared (IR) thermography can be used as a complementary tool for more global monitoring of the deposition of fast-ion losses on the wall, at the expense of no velocity-space resolution. IR cameras measure the temperature of the plasma facing components. This measurement, determined by a combined effect of the thermal plasma, radiation, neutrons and fast-ion losses, can be decomposed to infer the fast-ion load on the tokamak wall. In this manuscript, a workflow to estimate fast-ion losses via IR thermography is applied to the MAST-U spherical tokamak, using a 1D approximation to extract the experimental heat flux on the FILD front face from IR data. To numerically estimate the different contributions to this total heat flux, the field-line tracing environment SMITER is used to calculate the thermal plasma contribution, the orbit-following Monte-Carlo code ASCOT to estimate the fast-ion losses, and bolometry measurements for the radiation. To validate the workflow, two discharges, L-mode plasmas with low MHD activity, were executed using on and off-axis beams, respectively. The experimentally and numerically estimated heat flux are of the same order of magnitude for the on-axis heated scenario, with a strong dependence of the estimated fast-ion losses contribution on the fit to the kinetic profiles used as input. This is also true for the off-axis heated scenario, where the total numerically estimated heat flux is 2.1 or 1.3 times higher than the maximum experimentally estimated heat flux, depending on the ASCOT input used.

FILD↗
DOE OSTI2026

Preliminary neutronic and thermal analysis of Inertial Fusion Energy systems with a focus on thick liquid wall concepts

This paper provides an assessment of the neutronic behavior of Inertial Fusion Energy systems using the Serpent 2 Monte Carlo code. The study explores different types of reactor concepts, including dry-wall and thick liquid wall systems. In particular, it characterizes the tritium breeding and energy multiplication performance, as well as the effect of long-term irradiation on structural materials. For this purpose, the displacement per atom (DPA) and the gas production rate in structural materials are calculated. A simplified geometry with varying geometrical parameters is used as a representation for the different kinds of devices. Additionally, the case of the HYLIFE-II concept is used to assess modeling inaccuracies arising from certain approximations. A comparison between the simplified geometry and a realistic one is performed, thanks to Serpent’s capability to run transport calculations on unstructured meshes. The study then addresses the impact of these discrepancies on the temperature distribution in the first wall using the GeN-Foam code for conjugate heat transfer analysis. The results highlight the necessity to consider the heterogeneous nature of the liquid wall, as simplified models tend to underestimate radiation damage metrics and heat loads.

Conjugate heat transfer↗
DOE OSTI2021

Evaluation of gamma-ray transmission through rectangular collimator slits for application in nuclear fuel spectrometry

Gamma-ray spectrometry is widely applied in several science fields, and in particular in non-destructive gamma scanning and gamma emission tomography of irradiated nuclear fuel. Usually, a collimator is used in the experimental setup, to selectively interrogate a region of interest in the fuel. For the optimization of instrument design, as well as for planning measurement campaigns, predictive models for the transmitted gamma-ray intensity through the collimator are needed. Commonly, Monte Carlo Radiation Transport tools are used for accurate prediction of gamma-ray transport, however, the long computation time requirements when used in low-efficiency experimental setups present challenges. Here, the full-energy peak intensity transmitted through a rectangular collimator slit was examined. A uniform planar surface source emitting isotropically was considered, and the rate of photons reaching an ideal counter plane on the opposite side of the collimator was evaluated by analytical integration. To find a closed-form primitive function, some idealizations were required, and thereby parametric models were obtained for the optical field of view, dependent on slit dimensions (length, height and width) and source-to-collimator distance. For contributions from outside the optical field of view, where a closed-form expression cannot be found, fast numerical integral methods were instead used. The results were validated using the Monte Carlo code MCNP6 and show an agreement within three percent for the numerical method. For the analytical method, deviations up to tens of percent were obtained, which is deemed to still be sufficient for instrument design and measurement planning, where often the order of magnitude of the count rate is not a priori known. The method is planned for use in iterative optimization routines in the design of Gamma Emission Tomography devices, as well as for the prediction of gamma spectra obtained in the planning of fuel inspections. An application of the proposed method was demonstrated in spectrum prediction for a short cooling-time fuel rod test from the Halden reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Evaluation of the spatial self-shielding impact for TRISO-based nuclear fuel depletion

Reactor physics analyses of nuclear cores with nuclear fuel concepts containing tristructural isotropic (TRISO) particles, such as pebbles or compact fuel elements, rely on various degrees of simplification to keep these highly heterogeneous problems computationally tractable. One such limitation regards the level of spatial discretization employed during burnup calculations, where traditionally only a limited number of spatial zones are modeled at the full core level and assume that the spectrum is constant within the fuel elements and TRISO particles in this depletion zone. This type of assumption neglects the impact of spatial self-shielding effect within the kernels (microscale level) as well as within the compact or pebbles (mesoscale level). Furthermore, the Monte Carlo code Serpent 2 contains many relevant features for efficiently modeling this type of geometry, including a collision-based domain decomposition intended for very large burnup calculations, which we leveraged for this work to quantify the impact of capturing neutron flux variations occurring at the micro- and mesoscale level on a series of high-temperature gas-cooled reactor fuel element depletion problems. While spatial self-shielding is observed at both scales, with differences from a volume-averaged burnup of ±7% within the kernels and ±2% between TRISO particles within the fuel element, the conjugated effect on nuclide inventories and multiplication factor are negligible, hence confirming that assuming a single average spectrum value may be sufficient for most applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Preliminary analysis of TREAT free-field experiments using openmc

This work analyses activation calculations for dosimetry materials during a steady-state irradiation in the Transient Reactor Test (TREAT) reactor core. Hence, we developed a workflow based on the Monte Carlo code OpenMC alongside a custom depletion solver. The irradiation-induced activity as a function of time is computed, and several sensitivity studies are performed to evaluate uncertainty. This study has shown activity computations are sensitive to flux amplitude, irradiation time, atoms quantity and microscopic cross sections. Stochastic uncertainties have been propagated to evaluate the activity uncertainty for each dosimetry material. Most uncertainties are below our target of 3%, which demonstrates OpenMC as a powerful predictive and analysis tool. The precise results obtained through this newly developed computation scheme will be used in future experiments to characterize quantities of interest when operating the TREAT reactor in new configurations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Preservation of kinetics parameters generated by Monte Carlo calculations in two-step deterministic calculations

The generation of accurate kinetic parameters such as mean generation time Λ and effective delayed neutron fraction β eff via Monte Carlo codes is established. Employing these in downstream deterministic codes warrants another step to ensure no additional error is introduced by the low-order transport operator when computing forward and adjoint fluxes for bilinear weighting of these parameters. Another complexity stems from applying superhomogenization (SPH) equivalence in non-fundamental mode approximations, where reference and low-order calculations rely on a 3D full core model. In these cases, SPH factors can optionally be computed for only part of the geometry while preserving reaction rates and K-effective, but the impact of such approximations on kinetics parameters has not been thoroughly studied. This paper aims at studying the preservation of bilinearly-weighted quantities in the Serpent–Griffin calculation procedure. Diffusion and transport evaluations of IPEN/MB-01, Godiva, and Flattop were carried out with the Griffin reactor physics code, testing available modeling options using Serpent-generated multigroup cross sections and equivalence data. Verifying Griffin against Serpent indicates sensitivities to multigroup energy grid selection and regional application of SPH equivalence, introducing significant errors; these were demonstrated to be reduced through the use of a transport method together with a finer energy grid.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Development of broadband x-ray radiography for diagnosing magnetically driven cylindrically compressed matter

Experiments and modeling of x-ray radiography of millimeter diameter solid Al wires with laser-produced broadband x rays are reported. Experiments were performed using the 50-TW Leopard short-pulse laser in a laser and pulsed power chamber at the Nevada Terawatt Facility. To characterize broadband x rays and demonstrate a radiographic capability, bremsstrahlung, escaping electrons, and radiograph images of Al wires were simultaneously measured. The angularly resolved x-ray spectra are modeled by comparing measured bremsstrahlung signals in the range between 10 and ~500 keV with hybrid particle-in-cell simulations. Transmission of Al wires from the radiograph images is further simulated with a Monte Carlo code. The measured transmission profiles of Al wires with three different diameters agree with calculations when a simulated x-ray spectrum composed of line emissions and bremsstrahlung is used with a source size of 600 ± 200 μm. Transmission calculations with only 22 keV Ag Kα or an exponential x-ray spectrum do not reproduce the measurement, suggesting that the accurate determination of an x-ray source spectrum, as well as the inclusion of the photon sensitivity of the detector, is critical in transmission calculations to infer the density of an object. As a result, the laser-based broadband x-ray radiography that was developed has been successfully implemented in a pulsed power chamber for future laser-pulsed-power coupled experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Characterization of a 1D-imaging high-energy x-ray backlighter driven by the National Ignition Facility Advanced Radiographic Capability laser

Plastic deformation of samples compressed to Mbar pressures at high strain rates at the National Ignition Facility (NIF) forms the basis of ongoing material strength experiments in conditions relevant to meteor impacts, geophysics, armor development, and inertial confinement fusion. Hard x-ray radiography is the primary means of measuring the evolution of these samples, typically employing a slit-collimated high-Z microdot driven by the NIF laser to generate >40 keV x rays [E. Gumbrell et al., Rev. Sci. Instrum. 89, 10G118 (2018) and C. M. Huntington et al., Rev. Sci. Instrum. 89, 10G121 (2018)]. Alternatively, a dysprosium “micro-flag” target driven by the Advanced Radiographic Capability laser (~2 kJ, 10 ps) can deliver significantly higher spatiotemporal resolution [M. P. Hill et al., Rev. Sci. Instrum. 92, 033535 (2021)], especially in high-opacity samples. Here, initial experiments revealed problematic brightness and spectral gradients from this source, but by radiographing a set of diamond-turned, 105 µm-thick Pb test objects and supported by simulations using the 3D Monte Carlo code GEANT4, these geometry-dependent gradients across the field of view are quantified and mitigation strategies are assessed. In addition to significantly enhancing the modulation transfer function compared to the existing system, image stacking from multiple layers of image plate is shown to almost double the signal to noise ratio that will reduce uncertainties in future dynamic strength experiments.

47 OTHER INSTRUMENTATION↗