MCNP6.3 Executions in Parallel on Snow
Explore the source record for details and available documents.
MCNP, Unstructured Mesh↗
SEARCH · Search NASA
Results for “MCNP Unstructured Mesh”
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.
Evaluation of Attila and MCNP computational methods for dose and exposure estimation
Radiation transport calculations are often used to estimate dose or exposure to components and personnel surrounding a radiation source. The sources for these calculations are decaying radionuclides within various nuclear materials. Historically, dose calculations use MCNP (Monte Carlo N-Particle) transport code as the primary particle transport tool without a secondary computational tool to validate the results from the MCNP simulations [1]. The goal of this study is to make an independent check of the Monte Carlo solution from MCNP6 Version 6.2.1 with the discrete ordinates solution from Attila 10.2.0 Beta 3. As an example problem for this study, water-filled, stainless-steel vessels, modeled with an unstructured mesh (UM) with both MCNP and Attila [2], are exposed to 252Cf and 60Co point sources. This report also includes a discussion of the limitations of unstructured mesh in a MCNP calculation.
61 RADIATION PROTECTION AND DOSIMETRY↗
Niowave Neutron Source Converter: Lead-Bismuth-Eutectic (LBE) Windowless Target Design and Evaluation
Los Alamos National Laboratory (LANL) is working with Niowave on the design and evaluation of their lead-bismuth-eutectic (LBE) windowless target (i.e., neutron source converter). Niowave plans to use 200 kW electron beam at 40 MeV beam energy to produce neutrons by photonuclear reaction with LBE. Then, the neutrons undergo fission at the surrounding uranium target assembly (UTA) to produce Molybdenum 99 (Mo-99) as a fission product, which eventually decays to Technetium-99 (Tc-99m). Tc-99m is one of important radioisotopes that is used for medical diagnostics. LANL conducted 3D multiphysics analysis for the Niowave neutron converter design and provided design assessment in thermohydraulic aspects. LANL conducted radiation transport calculations using Monte-Carlo N-Particle (MCNP) code with unstructured meshing scheme. The 3D volumetric heating profiles in the LBE and Stainless-Steel (SS) housing were imported into multiphase computational fluid dynamics (CFD) to obtain 3D temperature profiles of LBE and SS through conjugate heat transfer (CHT) analysis. The key findings are: LBE film thickness at the center of the beam is approximately 1.6 cm with a maximum LBE velocity of approximately 1.8 m/s, which is below a 2 m/s limit to avoid erosion issues on supporting structures; Heat deposition in the LBE peaks at ~1 cm depth from the LBE free-surface because of the forward interactions of electron, photon, and neutron with LBE; LBE maximum temperature is ~360 °C, which is below LBE evaporation initiative temperature, ~450 °C; LBE-SS interface temperature is ~350 °C, which is below the safety thermal limit to prevent severe corrosion on SS. The results indicate that Niowave’s neutron converter design satisfies both hydraulic and thermal criteria for safe operation. By virtue of such computational analysis, Niowave can move toward establishing an experimental setup to experimentally test their LBE neutron converter. The following sections describe the detailed work done by LANL.
43 PARTICLE ACCELERATORS↗
Processing MCNP Elemental Edit Outputs
The Los Alamos National Laboratory's (LANL) Monte Carlo N-Particle® (MCNP®) transport code version 6 has the capability for tracking particles on unstructured mesh (UM) geometry models. The MCNP UM feature was originally designed for the UM models created by the Abaqus/CAE software product. The MCNP version 6.2.0 and later recognizes the UM models read from Abaqus INP files or MCNPUM files converted from Abaqus INP files. Other finite element analysis software packages may generate UM models and then convert these models into Abaqus INP formats.
97 MATHEMATICS AND COMPUTING↗
MCNP6.3 Unstructured Mesh Verification: GodivR and CANDU Models
A geometric cell of the Monte Carlo N-Particle (MCNP)1 transport code is traditionally created by using Boolean operators on defined surfaces. This constructive solid geometry (CSG) capability has been available in the MCNP code since its beginning. However, a CSG model approach is limited when it comes to constructing a representative geometry for a complex model in its ability to capture a correct model representation. Starting with the version 6.0, the MCNP code has the ability of embedding an unstructured mesh (UM) model into a CSG cell to create a hybrid geometry [1]. The MCNP UM feature provides the flexibility of defining very complex geometries because computer aided design (CAD) and mesh generation software packages can be utilized to construct UM models for MCNP simulations.
97 MATHEMATICS AND COMPUTING↗
Oktavian Modeling and Analysis with MCNP6.3 [Poster]
Project Goals: Apply MCNP6 variance reduction (VR) techniques to improve the Oktavian benchmark experiment calculations; Verify VR methods and electron transport with MCNP6.3 unstructured mesh (UM) geometry. All MCNP calculations in this poster are for mode n, p, e.
42 ENGINEERING↗
MCNP6 Element-Wise Densities and Temperatures
Explore the source record for details and available documents.
97 MATHEMATICS AND COMPUTING↗
Coupling of MCNP and ABAQUS for Reactor Multiphysics (MARM) [Slides]
Summer 2021 Internship Objectives: Learn MCNP (Monte Carlo N-Particle) transport code - Unstructured Mesh & Burnup Features, and Implement Python Code for MARM framework - MARM = MCNP and Abaqus based Reactor Multiphysics. Slides provide background and details for achieving learning objectives, as well as ideas for future work in the area.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Athena-I CUBIT Journal Files
The Monte Carlo N-Particle (MCNP)1 transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into con structive solid geometry (CSG) cells. This feature has been developed for performing calculations of complex geometry models because manually creating CSG models is a time-consuming and error prone process as the complexities of geometries increase. The UM capability was originally designed to work with UM models created with the Abaqus software and ASCII input files that it generates. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. Starting with version 6.3, MCNP can also process HDF5 mesh input files. External codes must be used to generate Abaqus input files for MCNP UM calculations. CUBIT, the Sandia National Laboratory automated mesh generation toolkit, can generate a UM model formatted as an Abaqus input file. However, the Abaqus input files exported from CUBIT cannot be used for MCNP simulations because it lacks the proper syntax. A Python script was developed to convert an Abaqus file created by CUBIT to an Abaqus file format that MCNP can process. Creating UM models for complex geometries is not an easy task. The process of creating UM models in CUBIT for MCNP simulations is detailed in. CUBIT provides several user interface options including a graphical user interface (GUI) and a command line interface. A GUI provides an easy way to use CUBIT without learning the CUBIT command syntax. When using CUBIT with either interface option, command lines are written into an ASCII file known as a journal file; this journal file can be edited and archived so that it can be played back in CUBIT to automatically generate a UM model. This report describes the CUBIT journal files of the UM models developed for Athena-I. The Athena platform, an energy-tuning assembly, was developed to spectrally shape the National Ignition Facility (NIF) deuterium-tritium fusion neutron source to a thermonuclear (fusion) plus prompt fission neutron spectrum with capability to act as a short-pulse neutron source. MCNP6 was used for the Athena experiment design analysis.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Verifying Burnup Calculations on Unstructured Mesh
Explore the source record for details and available documents.
MCNP UM ACTIUM AARE GOVIDA↗
Activation Analysis in Preparation for a Tungsten Irradiation Experiment at LANSCE
To organize the safe handling of activated material, knowing the residual dose rates is crucial. Here, in this work, we present the pre-experiment activation analysis for an experiment in which tungsten blocks are irradiated by 800-MeV protons. In this analysis, we use the Monte Carlo N-Particle (MCNP) code for radiation transport, Attila4MC for unstructured mesh generation, and Activation in Accelerator Radiation Environments (AARE), including CINDER2008, for activation analysis. If the tungsten blocks must be removed within a day after the experiment, then exposure to personnel entering the room must be reduced. One exposure-reduction strategy is to add carbon steel shielding around the tungsten blocks, efficiently reducing the dose from the activated tungsten. However, the shielding becomes activated itself during irradiation: 56Mn is the dominant contributor for short decay times. The actual schedule at the time of the experiment allowed sufficient cool-off time for the tungsten in the room so that additional shielding was not necessary. A less rigorous comparison of the calculated values with the post-experiment measurements showed reasonable agreement.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Optimization of the Second Target Station cold source moderators using an automated workflow
The Second Target Station (STS) at the US Department of Energy’s Oak Ridge National Laboratory is designed to become the world’s highest peak-brightness spallation source of cold neutrons. Successful completion of the STS, which is currently in the preliminary design phase, will provide transformative new capabilities to examine novel materials for future technologies. At STS, neutrons will be generated by spallation reactions in a solid tungsten target. They will be moderated and thermalized in two cold (20 K) para-hydrogen moderators. Careful optimization of these moderators is essential to the project’s success. To find optimal moderator designs, an advanced optimization workflow integrates high-fidelity neutronics calculations using the Monte Carlo N-Particle (MCNP) transport code MCNP6.2 with state-of-the-art optimization algorithms in the Dakota optimization toolkit. For each design iteration, a parametrized solid CAD geometry is generated in Creo and automatically converted into an unstructured mesh geometry by Attila 4MC for the neutronics calculation with MCNP. Iterations repeat until optimal designs are found. Herein this paper presents the results of a sensitivity and optimization study for the cylindrical and tube moderators. Both moderators can be optimized for maximum peak brightness, maximum time-integrated brightness, or any combination between these extremes. Maximum peak brightness is achieved by using smaller optimal dimensions of the moderators, whereas maximum time-integrated brightness is achieved by using larger dimensions. A Pareto front details the designs that optimally balance both brightness metrics. The Pareto front can be found in only 40–110 iterations with 4–5 design parameters when using the efficient global and Pareto-set optimization algorithms in Dakota. Additionally, important engineering constraints can be taken into account, such as the coupling between the cylindrical moderator radius and aluminum vessel wall thicknesses required to ensure structural integrity of the vessels. This interaction has a significant impact on the resulting optimal designs. Our new, highly efficient, fully automated optimization workflow will be used to optimize additional STS components in the future and can be adopted for design and optimization studies at other experimental neutron and accelerator facilities.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
An unstructured mesh based neutronics optimization workflow
We have developed a fully automated workflow to optimize the neutronics performance of the Second Target Station (STS) at the Oak Ridge National Laboratory’s Spallation Neutron Source. The optimization workflow starts with the parametrized solid CAD engineering models and converts them into the unstructured mesh (UM) models for the neutronics calculations with MCNP6.2. Calculations are executed and their results are loaded into the Dakota optimization toolkit. Dakota analyzes the results and proposes new geometry parameters for the next design iteration. The cycle repeats until the optimal parameters are found. The automated CAD to MCNP conversion, the use of high-fidelity UM models, and the use of modern optimizer are the key elements that advance the entire optimization workflow in comparison with the original workflow. The original workflow was based on a simplified constructive solid geometry (CSG) modeling with MCNPX, mcnp_pstudy tool, and an in-house optimizer. Herein to demonstrate the new workflow, we present a case of neutronics optimization of the moderator–reflector assembly (MRA). Apart from the MRA, the workflow can optimize other major STS components, such as the spallation target, neutron beamlines, radiation shielding, and various accelerator components. Importantly, the new workflow opens the door to the advanced multi-physics multi-parameter optimization and has the potential for use in other nuclear physics and accelerator applications.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Processing MCNP Elemental Edit Outputs
The Monte Carlo N-Particle (MCNP) transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. A UM geometry is a collection of elements representing a solid geometry. The first step of MCNP UM modeling is using other software packages to create a finite element mesh representation of a solid 3D geometry. Computer-aided design (CAD) or computer-aided manufacturing (CAM) software is typically used to create a solid geometry model, which is later imported into mesh generation software to create a UM model. The MCNP UM feature was originally designed for models generated by the Abaqus/CAE software. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. MCNP can process a UM model consisting of several different element types including linear tetrahedral or hexahedral elements and calculate quantities of interest such as flux and energy deposition at elements. An MCNP UM simulation provides high-fidelity elemental edit (i.e., tally) outputs, which can be further used in multiphysics calculations. The MCNP UM feature was used for multiphysics simulations where quantities of interest calculated by MCNP are used as inputs for heat transfer calculations in Abaqus. MCNP6.3 can produce two types of elemental edit output (EEOUT) file formats: ASCII and HDF5. An EEOUT file type must be requested on an EMBED card while output type (flux or energy deposition) must be requested on an EMBEE card. We wrote Python3 scripts to extract energy deposition values in an ASCII or HDF5 EEOUT file and compute a heat flux profile for an Abaqus heat transfer calculation.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Advancing the STS Neutron Moderator Design with an Automated Optimization Workflow and Unstructured Mesh Modeling
With the Second Target Station approaching its final design phase, a detailed neutronics evaluation of its critical components is necessary. Optimizing the dimensions of the two cold-source moderators that are at the heart of this facility presents a multi-objective optimization problem for which an accurate geometric description is crucial. We have applied a fully automated optimization workflow in which a detailed unstructured mesh geometry is automatically generated with Attila4MC, starting from a parametrized CREO geometry followed by preprocessing with SpaceClaim. With this geometry, a MCNP run is performed to calculate the brightness metrics, which are subsequently provided to the optimization algorithm in DAKOTA that provides new parameters and drives the optimization loop until convergence. In this paper, we show the results of the analysis that are used for the final design of the cylindrical and tube moderator. The optimization simulations provide a refinement to and confirmation of the conclusions of the previous design iteration. Additional to the optimization, a sensitivity study is performed to study the effect of minor geometry changes, which is important for the final engineering design. In conclusion, with these studies, we demonstrate that the automated workflow and high-fidelity unstructured mesh modeling are efficient tools for a thorough design evaluation.
DAKOTA↗
Athena-I Modeling with MCNP6.3 [Poster]
This objective of this project is to develop a model of the Athena-I experiment by creating unstructured mesh (UM) geometries in CUBIT and to perform calculations of the models in Monte Carlo N-Particle (MCNP) transport code version 6.3 (MCNP6.3).
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Two MCNP Models for Computational Performance Benchmarking
The purpose of this report is to describe two non-trivial MCNP models and execution configurations that are suitable to characterizing computational performance. To that end, this report uses a constructive solid geometry (CSG) representation of the Oak Ridge National Laboratory (ORNL) Pool Critical Assembly (PCA) [1–3] to perform a k-eigenvalue calculation and an unstructured mesh (UM) representation of the International Commission on Radiological Protection (ICRP) publication 145 (ICRP145) male human phantom [4, 5] to perform a fixed-source calculation assuming that a 1 MeV photon source is distributed throughout the phantom’s liver.
97 MATHEMATICS AND COMPUTING↗
$χ$-$MeRA$: Computationally efficient adaptive mesh refinement of Monte Carlo mesh based tallies
Here, the reactor physics community is always focused on reducing the computational time and memory required for simulations. $χ$-$MeRA$, which stands for flux-based-($χ$)-Mesh tally Refinement Adaptively, was built to reduce the computational time and memory required to solve the neutronics side of a multiphysics problem when compared to traditional methods for mesh based tallies in Monte Carlo (MC) simulations. $χ$-$MeRA$ couples a MC code with an adaptive mesh refinement (AMR) algorithm to take advantage of the accuracy of a MC code and the efficiency of an AMR algorithm. Also developed within $χ$-$MeRA$ was a set of metrics to assess the effects of the refinement on various parameters in the simulation space. For a plutonium sphere, $χ$-$MeRA$ shows a reduction in memory usage and computation time when compared to a fully refined mesh by a factor of 14.7 and 6.7, respectively. When compared to an unstructured mesh, improvement of 1.3 and 4.8 was achieved for memory usage and computation time. The development of $χ$-$MeRA$ helps solve the neutronics side of a multiphysics problem in a faster, more computationally efficient manner than traditional methods, and the final mesh created contains accurate results that can be passed onto the next physics code.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗