Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,675 records · Page 93

Accelerating detector simulations with Celeritas: Profiling and performance optimizations

Celeritas is a GPU-optimized Monte Carlo (MC) particle transport code designed to meet the growing computational demands of next-generation high energy physics (HEP) experiments. It provides efficient simulation of electromagnetic (EM) physics processes in complex geometries with magnetic fields, detector hit scoring, and seamless integration into Geant4-driven applications to offload EM physics to GPUs. Recent efforts have focused on performance optimizations and expanding profiling capabilities. This paper presents some key advancements, including the integration of the Perfetto system profiling tool for detailed performance analysis and the development of track-sorting methods to improve computational efficiency.

Lund, Amanda [Argonne National Laboratory (ANL)]↗

Comparative studies of actinide and sub-actinide fission cross section calculation from MCNP6 and TALYS

Comparative studies of actinide and sub-actinide fission cross section calculation from MCNP6 and TALYS have been conducted. In this work, fission cross section resulted from MCNP6 prediction will be compared with result from TALYS calculation. MCNP6 with its event generator CEM03.03 and LAQGSM03.03 have been validated and verified for several intermediate and heavy nuclides fission reaction data and also has a good agreement with experimental data for fission reaction that induced by photons, pions, and nucleons at energy from several ten of MeV to about 1 TeV. The calculation that induced within TALYS will be focused mainly to several hundred MeV for actinide and sub-actinide nuclides and will be compared with MCNP6 code and several experimental data from other evaluator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gamma-ray dose from an overhead plume

Standard plume models can underestimate the gamma-ray dose when most of the radioactive material is above the heads of the receptors. Typically, a model is used to calculate the air concentration at the height of the receptor, and the dose is calculated by multiplying the air concentration by a concentration-to-dose conversion factor. Models indicate that if the plume is emitted from a stack during stable atmospheric conditions, the lower edges of the plume may not reach the ground, in which case both the ground-level concentration and the dose are usually reported as zero. However, in such cases, the dose from overhead gamma-emitting radionuclides may be substantial. Such underestimates could impact decision making in emergency situations. The Monte Carlo N-Particle code, MCNP, was used to calculate the overhead shine dose and to compare with standard plume models. At long distances and during unstable atmospheric conditions, the MCNP results agree with the standard models. As a result, at short distances, where many models calculate zero, the true dose (as modeled by MCNP) can be estimated with simple equations.

61 RADIATION PROTECTION AND DOSIMETRY↗

Induced Radioactivity in Lead Shielding at the National Synchrotron Light Source

The National Synchrotron Light Source (NSLS) at Brookhaven National Laboratory was shut down in September 2014. Lead bricks used as radiological shadow shielding within the accelerator were exposed to stray radiation fields during normal operations. The FLUKA code, a fully integrated Monte Carlo simulation package for the interaction and transport of particles and nuclei in matter, was used to estimate induced radioactivity in this shielding and stainless steel beam pipe from known beam losses. The FLUKA output was processed using MICROSHIELD® to estimate on-contact exposure rates with individually exposed bricks to help design and optimize the radiological survey process. This entire process can be modeled using FLUKA, but use of MICROSHIELD® as a secondary method was chosen because of the project’s resource constraints. Due to the compressed schedule and lack of shielding configuration data, simple FLUKA models were developed in this paper. FLUKA activity estimates for stainless steel were compared with sampling data to validate results, which show that simple FLUKA models and irradiation geometries can be used to predict radioactivity inventories accurately in exposed materials. During decommissioning 0.1% of the lead bricks were found to have measurable levels of induced radioactivity. Finally, post-processing with MICROSHIELD® provides an acceptable secondary method of estimating residual exposure rates.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Production of Th 229 for medical applications: Excitation functions of low-energy protons on Th 232 targets

As a part of a general program to evaluate production routes for 229 Th, we studied production of 229 Th via proton-induced reactions on 232 Th targets bombarded with low-energy protons, E p ≤ 40MeV. The reported excitation functions include those for proton-induced reactions on natural thorium yielding to 228,229,230&232 Pa isotopes; 232 Th( p,xn ) reactions, where x=1, 3, 4, and 5, at proton energy ranges of 12–40 MeV. Although the data for 232 Th(p,n) 228 Pa, 232 Th(p,3n) 230 Pa, and 232 Th(p,5n) 232 Pa reactions were deduced by direct analysis of the thorium foils after irradiation, the data for 232 Th(p,4n) 229 Pa were obtained by radiochemical techniques. The half-life of 229 Pa was evaluated and determined to be 1.55 ± 0.01 d. Further, the α-branching ratio, α/(α + EC) of 229 Pa was evaluated to be 0.53 ± 0.10% by allowing 229 Pa to decay for ~7d, then chemically extracting and quantifying the 225 Ac (t 1/2 = 10.0±0.1d) from 229 Pa samples. In addition, we report the effective production cross section of 229 Th in a thick 232 Th target in the proton energy range of 23–33 MeV. The peak of the excitation function for the 232 Th(p,4n) 229 Pa reaction occurs at 162 ± 14 mb and E p =29.7±0.5MeV. This is only slightly larger than the effective cross section for the 232 Th(p,x) 229 Th reaction (obtained from a thick target experiment). This data indicates that the 232 Th(p,4n) 229 Pa reaction is the major reaction pathway for the cumulative 232 Th(p,x) 229 Th reaction cross section in this energy range. The measured cross sections were compared with theoretical cross sections using the simulation codes Particle and Heavy Ion Transport code System (PHITS) and Monte Carlo Neutral Particle 6 (MCNP6). Furthermore, at proton energy ranges of 12–33 MeV, the cumulative excitation function predicted by PHITS for the reactions leading to 229 Th was in close agreement with the experimental function, whereas the function predicted by MCNP6 was a factor of two higher at the peak of the excitation function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

First Industrial Tests of a Matrix Monitor Correction for the Differential Die-away Technique of Historical Waste Drums

The fissile mass in radioactive waste drums filled with compacted metallic residues (spent fuel hulls and nozzles) produced at AREVA NC La Hague reprocessing plant is measured by neutron interrogation with the Differential Die-away measurement Technique (DDT). In the next years, old hulls and nozzles mixed with Ion-Exchange Resins will be measured. The ion-exchange resins increase neutron moderation in the matrix, compared to the waste measured in the current process. In this context, the Nuclear Measurement Laboratory (LMN) of CEA Cadarache has studied a matrix effect correction method, based on a drum monitor, namely a 3He proportional counter located inside the measurement cavity. After feasibility studies performed with LMN's PROMETHEE 6 laboratory measurement cell and with MCNPX simulations, this paper presents first experimental tests performed on the industrial ACC (hulls and nozzles compaction facility) measurement system. A calculation vs. experiment benchmark has been carried out by performing dedicated calibration measurements with a representative drum and {sup 235}U samples. The comparison between calculation and experiment shows a satisfactory agreement for the drum monitor. The final objective of this work is to confirm the reliability of the modeling approach and the industrial feasibility of the method, which will be implemented on the industrial station for the measurement of historical wastes. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evaluation of brachytherapy lung implant dose distributions from photon-emitting sources due to tissue heterogeneities

Purpose: Photon-emitting brachytherapy sources are used for permanent implantation to treat lung cancer. However, the current brachytherapy dose calculation formalism assumes a homogeneous water medium without considering the influence of radiation scatter or tissue heterogeneities. The purpose of this study was to determine the dosimetric effects of tissue heterogeneities for permanent lung brachytherapy. Methods: The MCNP5 v1.40 radiation transport code was used for Monte Carlo (MC) simulations. Point sources with energies of 0.02, 0.03, 0.05, 0.1, 0.2, and 0.4 MeV were simulated to cover the range of pertinent brachytherapy energies and to glean dosimetric trends independent of specific radionuclide emissions. Source positions from postimplant CT scans of five patient implants were used for source coordinates, with dose normalized to 200 Gy at the center of each implant. With the presence of fibrosis (around the implant), cortical bone, lung, and healthy tissues, dose distributions and {sub PTV}DVH were calculated using the MCNP *FMESH4 tally and the NIST mass-energy absorption coefficients. This process was repeated upon replacing all tissues with water. For all photon energies, 10{sup 9} histories were simulated to achieve statistical errors (k = 1) typically of 1%. Results: The mean PTV doses calculated using tissue heterogeneities for all five patients changed (compared to dose to water) by only a few percent over the examined photon energy range, as did PTV dose at the implant center. The {sub PTV}V{sub 100} values were 81.2%, 90.0% (as normalized), 94.3%, 93.9%, 92.7%, and 92.2% for 0.02, 0.03, 0.05, 0.1, 0.2, and 0.4 MeV source photons, respectively. Relative to water, the maximum bone doses were higher by factors of 3.7, 5.1, 5.2, 2.4, 1.2, and 1.0 The maximum lung doses were about 0.98, 0.94, 0.91, 0.94, 0.97, and 0.99. Relative to water, the maximum healthy tissue doses at the mediastinal position were higher by factors of 9.8, 2.2, 1.3, 1.1, 1.1, and 1.1. However, the maximum doses to these healthy tissues were only 3.1, 7.2, 11.3, 10.9, 9.0, and 8.1 Gy while maximum bone doses were 66, 177, 236, 106, 49, and 39 Gy, respectively. Similarly, maximum lung doses were 55, 66, 73, 74, 73, and 73 Gy, respectively. Conclusions: The current brachytherapy dose calculation formalism overestimates PTV dose and significantly underestimates doses to bone and healthy tissue. Further investigation using specific brachytherapy source models and patient-based CT datasets as MC input may indicate whether the observed trends can be generalized for low-energy lung brachytherapy dosimetry.

61 RADIATION PROTECTION AND DOSIMETRY↗

Precision Description of Processes at Colliders in the SANC System

A review of the SANC computer system intended for theoretical support of experiments at colliders is presented. The system allows the precision predictions to be obtained semiautomatically for a wide class of particle interaction processes within the Standard Model. In this case, one considers total contributions of one-loop radiative corrections. Basic elements of the computer system are described, and examples of using the system for the high-precision description of certain processes of the Standard Model particle interactions are given. The MCSANC code developed using the Monte Carlo simulation is presented, which is used to calculate cross sections of a number of processes, investigated at the Large Hadron Collider. The prospects of research development in this field for theoretical support at future hadron and lepton colliders are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

lanl/jaybenne

Monte Carlo thermal radiation transport code

Ryan, Ben↗

Fully Implicit Conjugate Heat Transfer Analysis of the ARC-Class Vacuum Vessel

The coupled simulation of fusion reactor blankets including neutronics, thermal-hydraulics and thermo-mechanics is expected to speed up the design cycle of fusion reactor design concepts. In this work we demonstrate tight implicit coupling of conjugate heat transfer using the open-source Computational Fluid Dynamics software OpenFOAM for thermo-fluid mechanics and Diablo for thermo-solid mechanics. The heat transfer analysis is augmented by volumetric energy deposition from neutronic calculations using the Monte Carlo N-particle code on both solid and fluid parts of the vacuum vessel. An additional heat flux is imposed on the first wall estimated from the design power of the reactor. The tight coupling is realized through the open-source coupling library, preCICE, and tested on the vacuum vessel of the affordable, robust, compact reactor design by Commonwealth Fusion Systems. The features of the coupling and the influence of different coupling parameters such as coupling schemes, acceleration techniques and convergence criterion are discussed. The coupled simulation results are compared to a thermal-hydraulics simulation which includes only the fluid domains (the liquid immersion molten salt blanket and cooling channel) to demonstrate usefulness of a coupled simulation. Further analysis is performed to identify regions of hot spots for subsequent design improvement. This introduces the outline for integrating conjugate electromagnetics and fluid/solid mechanics (e.g., allow for deformation of the cooling channel walls) with our present approach for future analysis.

Sircar, Arpan↗

Sensitivity-based Similarity Metrics for New Experiment Design Optimization

The nuclear data used in advanced reactor simulations requires validation. Data from nuclear criticality experiments can provide this validation. New nuclear criticality experiment design requires extensive knowledge and expert judgement such that the experimental design parameters are selected in such a way to keep the experiment subcritical. To aide in this experimental design process, professionals can utilize sensitivity and uncertainty analysis. Sensitivity and uncertainty analysis relies on matching new application experiments with currently existing benchmark experiments. Currently, there is functionality in the Whisper 1.1 software package to calculate a similarity metric based on neutron multiplication factor sensitivity coefficients between a new application designed by the user and existing International Criticality Safety Benchmark Experiment Project (ICSBEP) benchmarks. The Whisper 1.1 software package is included in Monte Carlo N-Particle ® Code Version 6.21 (MCNP ® 6.2). This work is geared toward expanding this capability to new similarity metrics based on beta-effective sensitivity coefficients and reactivity coefficient sensitivity coefficients. While the investigation of these sensitivity coefficients is presented in detail in separate works at this same conference, this work will be primarily focused on studying the similarity metrics in more detail. These similarity metrics will then be incorporated into the optimization algorithms used for experiment design in EUCLID (Experiments Underpinned by Computational Learning for Improvements in nuclear Data), which is a Los Alamos National Laboratory (LANL) project designed to constrain nuclear data of interest, such that adjustments can be made to possible inaccuracies. A more detailed optimization can be subsequently performed by breaking down these similarity metrics by isotope, reaction, and energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

2012 SARA Students Technical Report

The Service Academy Research Associates (SARA) program provides an opportunity for Midshipmen and Cadets from US Service Academies to participate in research at Los Alamos National Laboratory (LANL), Lawrence Livermore National Laboratory (LLNL), and Sandia National Laboratory for several weeks during the summer as part of their summer training assignments. During the summer of 2012, three Midshipmen were assigned to work with the XCP Division at LANL for approximately 5-6 weeks. As one of the nation's top national security science laboratories, LANL stretches across 36 square miles, has over 2,100 facilities, and employs over 9,000 individuals including a significant number of students and postdocs. LANL's mission is to 'apply science and technology to: ensure the safety, security, and reliability of the US nuclear deterrent, reduce global threats, and solve other emerging national security challenges.' While LANL officially operates under the US Department of Energy (DoE), fulfilling this mission requires mutual cooperation with the US Department of Defense (DoD) as well. LANL's high concentration of knowledge and experience provides interns a chance to perform research in many disciplines, and its connection with the DoD in both operation and personnel gives SARA students insight to career possibilities both during and after military service. SARA students have plenty of opportunity to enjoy hiking, camping, the Los Alamos YMCA, and many other outdoor activities in New Mexico while staying at the Buffalo Thunder Resort, located 20 miles east of the lab. XCP Division is the Computational Physics division of LANL's Weapons Department. Working with XCP Division requires individuals to be Q cleared by the DoE. This means it is significantly more convenient for SARA students to be assigned to XCP Division than their civilian counterparts as the DoD CNWDI clearance held by SARA students is easily transferred to the lab prior to the students arriving at the start of the summer. SARA students working with XCP Division were given a comprehensive introduction into nuclear engineering and physics, nuclear weapons, and radiation transport and detection via texts and lectures at various classification levels. Students also attended tours of several prominent facilities at LANL including TA-41 Ice House, TA-55 PF-4 plutonium facility, the Nicholas C. Metropolis Center for Modeling and Simulation, also known as the Secure Computing Center (SCC), and the Dual-Axis Radiological Hydro Test (DARHT) facility; in addition, SARA students accompanied by LANL staff traveled to Minot AFB in North Dakota for tours of the 5th Bomb Wing and 91st Missile Wing facilities. Students participated in a week long class on the Monte Carlo N Particle (MCNP) code to supplement their understanding of radiation transport simulations. SARA students were then tasked with using this knowledge to model radiation detectors and use MCNP to compare their models to experimental data and previously accepted models.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Characterization of Delayed-Particle Emission Signatures for Pyroprocessing, Part 1: Fuel Assembly

An initial study was conducted with nondestructive assay analysis using the MCNP6 Monte Carlo radiation-transport code to calculate delayed-neutron and delayed- gamma emission signature monitors to characterize burnup and isotopic composition for pyroprocessing. The study consists of five parts. In Part 1, a generic Westinghouse fuel assembly model was modified to create models with 3% and 5% 235 U enrichments, 20, 30, 40, and 50 GWd/MTU burnups, and 3-, 5-, 10-, 20-, and 30-year cooling times. MCNP6 was used to create detail-rich delayed-particle source terms for fuel assembly. Delayed-gamma activity ratios were then developed for conventional 134 Cs/ 137 Cs, 134 Cs/ 154 Eu, and 154 Eu/ 137 Cs markers as a function of burnup and of actinide mass for elemental uranium, neptunium, plutonium, americium, and curium. In addition, data mining using photon tally data was employed to discover potential new delayed-gamma activity ratios. The data mining effort was streamlined by the creation and use of several utility codes to extract and process data. The study confirms behavior for burnup dependence of 134 Cs/ 137 Cs, 134 Cs/ 154 Eu, and 154 Eu/ 137 Cs markers. A process for the discovery of potential new delayed-gamma markers has been developed, and results are presented. The actinide and nonactinide inventories for the fuel assembly models are used as source materials for the electrorefining models. Parts 2–5 contain analysis of a pyrochemical separations model that includes an MCNP6 mockup of the INL Mark IV electrorefiner (see companion documents). Part 2 focuses on delayed-particle emissions from the electrorefiner with loaded fuel dissolution baskets before electrorefining. Part 3 investigates delayed-particle emissions from the electrorefiner following electrorefining. Part 4 presents delayed-particle emissions from other pyroprocessing units that contain material developed by the electrorefining process in Part 3. Part 5 parallels Part 4, except for a 10 wt% plutonium reduction.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Delayed-Neutron and Delayed-Photon Dose-Rate Estimations for Chemical Processing of Spent Fuel Using MCNP6

A nondestructive assay (NDA) computational study is conducted using the MCNP6 Monte Carlo radiation-transport code to calculate dose estimates caused by delayed-neutron (DN) and delayed-gamma (DG) emission for spent nuclear fuel (SNF) that is immersed in nitric acid. The SNF isotopic inventories and delayed-particle (DP) sources are prepared using a generic Westinghouse fuel assembly model for three sets of parameters.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

MCNP6 Class

This report provides a Monte Carlo method for simulating radiation transport and provides a common sense approach for MCNP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Impact of the High Flux Isotope Reactor HEU to LEU Fuel Conversion on Cold Source Nuclear Heat Generation Rates

Under the sponsorship of the US Department of Energy National Nuclear Security Administration, staff members at the Oak Ridge National Laboratory have been conducting studies to determine whether the High Flux Isotope Reactor (HFIR) can be converted from high enriched uranium (HEU) fuel to low enriched uranium (LEU) fuel. As part of these ongoing studies, an assessment of the impact that the HEU to LEU fuel conversion has on the nuclear heat generation rates in regions of the HFIR cold source system and its moderator vessel was performed and is documented in this report. Silicon production rates in the cold source aluminum regions and few-group neutron fluxes in the cold source moderator were also estimated. Neutronics calculations were performed with the Monte Carlo N-Particle code to determine the nuclear heat generation rates in regions of the HFIR cold source and its vessel for the HEU core operating at a full reactor power (FP) of 85 MW(t) and the reference LEU core operating at an FP of 100 MW(t). Calculations were performed with beginning-of-cycle (BOC) and end-of-cycle (EOC) conditions to bound typical irradiation conditions. Average specific BOC heat generation rates of 12.76 and 12.92 W/g, respectively, were calculated for the hemispherical region of the cold source liquid hydrogen (LH2) for the HEU and LEU cores, and EOC heat generation rates of 13.25 and 12.86 W/g, respectively, were calculated for the HEU and LEU cores. Thus, the greatest heat generation rates were calculated for the EOC HEU core, and it is concluded that the conversion from HEU to LEU fuel and the resulting increase of FP from 85 MW to 100 MW will not impact the ability of the heat removal equipment to remove the heat deposited in the cold source system. Silicon production rates in the cold source aluminum regions are estimated to be about 12.0% greater at BOC and 2.7% greater at EOC for the LEU core in comparison to the HEU core. Silicon is aluminum s major transmutation product and affects mechanical properties of aluminum including density, neutron irradiation hardening, swelling, and loss of ductility. Because slightly greater quantities of silicon will be produced in the cold source moderator vessel for the LEU core, these effects will be slightly greater for the LEU core than for the HEU core. Three-group (thermal, epithermal, and fast) neutron flux results tallied in the cold source LH2 hemisphere show greater values for the LEU core under both BOC and EOC conditions. The thermal neutron flux in the LH2 hemisphere for the LEU core is about 12.4% greater at BOC and 2.7% greater at EOC than for the HEU core. Therefore, cold neutron scattering will not be adversely affected and the 4 12 neutrons conveyed to the cold neutron guide hall for research applications will be enhanced.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Characterization of Delayed-Particle Emission Signatures for Pyroprocessing. Part 1: ABTR Fuel Assembly

A three-part study is conducted using the MCNP6 Monte Carlo radiation-transport code to calculate delayed-neutron (DN) and delayed-gamma (DG) emission signatures for nondestructive assay (NDA) metal-fuel pyroprocessing. In Part 1, MCNP6 is used to produce irradiation-induced used nuclear fuel (UNF) isotopic inventories for an Argonne National Laboratory (ANL) Advanced Burner Test Reactor (ABTR) preconceptual design fuel assembly (FA) model. The initial fuel inventory consists of uranium mixed with light-water-reactor transuranic (TRU) waste and 10 wt% zirconium (U-LWR-SFTRU-10%Zr). To facilitate understanding, parametric evaluation is done using models for 3% and 5% initial 235 U a% enrichments, burnups of 5, 10, 15, 20, 30, …, 120 GWd/MTIHM, and 3-, 5-, 10-, 20-, and 30- year cooling times. Detailed delayed-particle radioisotope source terms for the irradiate FA are created using BAMF-DRT and SOURCES3A. Using simulation tallies, DG activity ratios (DGARs) are developed for 134 Cs/ 137 Cs 134 Cs/ 154 Eu, and 154 Eu/ 137 Cs markers as a function of (1) burnup and (2) actinide mass, including elemental uranium, neptunium, plutonium, americium, and curium. Spectral-integrated DN emission is also tallied. The study reveals a rich assortment of DGAR behavior as a function of DGAR type, enrichment, burnup, and cooling time. Similarly, DN emission plots show variation as a function of burnup and of actinide mass. Sensitivity of DGAR and DN signatures to initial 235 U enrichment, burnup, and cooling time is evident. Comparisons of the ABTR radiation signatures and radiation signatures previously reported for a generic Westinghouse oxide-fuel assembly indicate that there are pronounced differences in the ABTR and Westinghouse oxide-fuel DN and DG signatures. These differences are largely attributable to the initial TRU inventory in the ABTR fuel. The actinide and nonactinide inventories for the FA models serve as source materials for the pre- and postelectrorefining models to be reported in Parts 2 and 3.

07 ISOTOPE AND RADIATION SOURCES↗

Simulation of Hohlraum Wall Texture for Improved Performance in Hohlraums

The performance of hohlraums for use in NIF (National Ignition Facility) is explored using Cassio: a LANL produced radiation hydrodynamics code that implements implicit Monte Carlo radiation transport. We examine the effect on blowoff of adding texture to the inside of the hohlraum wall. These new designs are compared in simulation against current designs and the data are analyzed for the possible use of such hohlraums in future high energy density physics experiments.

42 ENGINEERING↗