Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mcnp”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Integration of the Full Tokamak Reference Model with the Complex Model for ITER Neutronic Analysis

The ITER International Organization has developed a number of reference Monte Carlo N-Particle (MCNP) models including the tokamak machine C-model, the Tokamak Complex model, and the neutral beam injection (NBI) systems model. The Tokamak Complex model primarily describes building structures beyond the bioshield. Representation of the tokamak and its systems are not included in this model. The Oak Ridge National Laboratory Radiation Transport Group has conducted two ITER neutronic analysis model integrations: (1) integration of the tokamak C-model with the Tokamak Complex model for shutdown dose rate characterization in Port Cell 16 at level B1, and (2) integration of the NBI model with the Tokamak Complex model for estimating the spatial distribution of biological dose rate at levels L1, L2, and L3 of the Tokamak Complex. The integrated models were further extended to include models of system components that are essential to the neutronic analyses. Furthermore, this paper presents the approach and computer tools used to integrate existing reference models, describes the additional design details implemented in the integrated models, and provides representative neutronic calculations based on the extended models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Algorithmic Improvements to MCNP5 for High-Resolution Fusion Neutronics Analyses

Neutronics analyses of the ITER experimental fusion reactor rely on increasingly complex geometry models and estimates of energy-dependent neutron flux and radiation dose-rate distributions generated at ever higher resolutions. There are significant practical challenges with applying the Monte Carlo N-Particle (MCNP) continuous-energy transport code to high-resolution analyses. For models consisting of more than 100 000 surfaces and cells, geometry initialization can take several hours, thus slowing down model integration and transport analysis efforts. In multithreaded simulations, the amount of memory consumed by superimposed mesh tally data increases in proportion to the number of threads. This behavior limits either the tally resolution or the number of processor cores that can be utilized in the simulation. This paper describes algorithmic improvements that were implemented in a modified version of MCNP5 to overcome these limitations. These improvements are referred to as the Oak Ridge National Laboratory Transformative Neutronics (ORNL-TN) upgrade. Here, a comparison of the performance and memory usage of both MCNP5 and ORNL-TN on several relevant fusion neutronics models is presented. In these tests and in actual high-resolution neutronics analyses, ORNL-TN reduces geometry processing times from hours to a few seconds and increases in-memory mesh tally capacity from the order of 10 8 to 10 10 space-energy bins.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Application of the Denovo Discrete Ordinates Radiation Transport Code to Large-Scale Fusion Neutronics

Fusion energy systems pose unique challenges to the modeling and simulation community. These challenges must be met to ensure the success of the ITER experimental fusion reactor. ITER’s complex systems require detailed modeling that goes beyond the scale of comparable simulations to date. In this work, the Denovo radiation transport code was used to calculate neutron fluence and kerma for the JET streaming benchmark. This work was performed on the Titan supercomputer at the Oak Ridge Leadership Computing Facility. Denovo is a novel three-dimensional discrete ordinates transport code designed to be highly scalable. Sensitivity studies have been completed to examine the impact of several deterministic parameters. Furthermore, results were compared against experiment as well as the MCNP and Shift Monte Carlo codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

MHD Analysis of Dual-Coolant Lead-Lithium Blanket for Spherical Tokamak Advanced Reactor

The tritium breeding blanket is vital for future fusion power plants, with the Spherical Tokamak Advanced Reactor (STAR) project highlighting the dual-coolant lead-lithium (DCLL) design. The DCLL blanket performs shielding, energy exhaust, and tritium breeding using a lead-lithium alloy , with lithium as the breeder and lead as the neutron multiplier. It also serves as the primary coolant, with helium providing supplemental cooling. Reduced-activation ferritic/martensitic steel is used for the blanket structure. Magnetohydrodynamic (MHD) phenomena influence the liquid metal flow in a magnetic field, affecting heat transfer in the breeder affected by energetic neutrons. Understanding key flow parameters in such conditions is critical for efficient DCLL design. This study uses three-dimensional thermofluid MHD analysis with ANSYS CFX software, modified at Princeton Plasma Physics Laboratory, to simulate high Hartmann flows. The neutronics code MCNP, coupled with plasma equilibrium, provides heat source distribution. In conclusion, we examine electromagnetic interactions in adjacent fluid domains and analyze the magnetic field’s impact on flow distribution in the inboard and outboard blanket layout, using detailed mesh generation for accurate results.

DCLL↗
DOE OSTI2024

Transfer functions for Q A /Q B international regulatory limits for the safe transport of radioactive materials

This paper presents a proposed revision of the International Atomic Energy Agency transport regulations, related to the A 1 and A 2 limit values used to determine the radioactive transport classification. Based on the 'Q system', a novel methodology was introduced to derive Q A and Q B values related to scenarios involving external exposure from a distant source. These values are key parameters that respectively represent the total effective dose and total equivalent dose to the skin, from all primary and secondary particles contributing to radiation exposure. The International Working Group (WG A 1 /A 2 ) is established and associated with the TRANSSC Technical Expert Group on Radiation Protection. A review of the A 1 and A 2 values is performed in response to identified limitations within the existing Q system. The followed approach is based on Monte Carlo simulations that enabled the development of transfer functions aimed at reducing computational time and increasing the flexibility of dose evaluations for any radionuclide with known particle emission spectra. This method allows updating the Q A and Q B values to account for future data evolutions (decay data, fluence-to-dose conversion coefficients) and standardizing the calculation of regulation limits across all referenced radionuclides and scenarios related to external exposure. The transfer functions are established using three Monte Carlo simulation codes—FLUKA, Geant4, and MCNP—and address the previous limitations of the 'Q system', reflecting the latest International Commission for Radiation Protection recommendations and improvements in calculation techniques. The results of the WG show consistent agreement across the codes, with minor discrepancies observed at low primary energies due to statistical uncertainties and different handling of stopping power for electrons/positrons in the codes. This revised approach aligns with current standards and recommendations, ensuring that the radiological consequences of transport accidents are acceptable for the new A 1 and A 2 limits from a radiological protection perspective.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2019

Multiphysics approach to plasma neutron source modelling at the JET tokamak

A novel multiphysics methodology for the computation of realistic plasma neutron sources has been developed. The approach is based on state-of-the-art plasma transport and neutron spectrum calculations, coupled with a Monte Carlo neutron transport code, bridging the gap between plasma physics and neutronics. In the paper two JET neutronics tokamak models are used to show the application of the developed plasma neutron sources and validate them. Diagnostic data for the record JET D discharge 92436 are used as input for the TRANSP code, modelling neutron emission in two external plasma heating scenarios, namely using only neutral beam injection and a combination of the latter and ion cyclotron resonance heating. Neutron spectra, based on plasma transport results, are computed using the DRESS code. The developed PLANET code package is employed to generate plasma neutron source descriptions and couple them with the MCNP code. The effects of using the developed sources in neutron transport calculations on the response of JET neutron diagnostic systems is studied and compared to the results obtained with a generic plasma neutron source. It is shown that, although there are significant differences in the emissivity profiles, spectra shape and anisotropy between the neutron sources, the integral response of the time-resolved ex-vessel neutron detectors is largely insensitive to source changes, with major relative deviations of up to several percent. However it is calculated that, because of the broadening of neutron spectra as a consequence of external plasma heating, larger differences may occur in activation of materials which have threshold reactions located at DD neutron peak energies. The PLANET plasma neutron source computational methodology is demonstrated to be suitable for detailed neutron source effect studies on JET during DT experiments and can be applied to ITER analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental validation of an integrated modelling approach to neutron emission studies at JET

An integrated modelling methodology for the calculation of realistic plasma neutron sources for the JET tokamak has been developed. The computational chain comprises TRANSP plasma transport and DRESS neutron spectrum calculations, and their coupling to the MCNP neutron transport code, bridging plasma physics and neutronics. Here, we apply the developed methodology to the analysis of neutron emission properties of deuterium and helium plasmas at JET, and validate individual modelling steps against neutron diagnostic measurements. Two types of JET discharges are modelled—baseline-like and three-ion radio-frequency scenarios—due to their diversity in plasma heating, characteristics of the induced fast ion population, and the imprint of these on neutron emission properties. The neutron emission modelling results are quantitatively compared to the total neutron yield from fission chambers, neutron emissivity profiles from the neutron camera, neutron spectra from the time-of-flight spectrometer, and neutron activation measurements. The agreement between measured and calculated quantities is found to be satisfactory for all four diagnostic systems within the estimated experimental and computational uncertainties. Additionally, the effect of neutrons not originating from the dominating D(D, n) 3 He reactions is studied through modelling of triton burnup DT neutrons, and, in mixed D- 3 He plasmas, neutrons produced in the 9 Be(D, nγ) 10 B reaction on impurities. It is found that these reactions can contribute up to several percent to the total neutron yield and dominate the neutron activation of samples. Finally, the effect of MeV-range fast ions on the neutron activation of 115 In and 27 Al samples is measured and computationally validated.

integrated modelling↗
DOE OSTI2022

AR reactor neutronics multi-code validation*

The affordable, robust, compact (ARC) reactor is a tokamak fusion reactor concept currently under development by Commonwealth Fusion Systems and Massachusetts Institute of Technology. There are three important neutronics considerations for the operation of the ARC reactor: (1) breeding of enough tritium in the blanket to sustain the D–T reaction in the plasma; (2) ensuring low fluence on the superconducting toroidal field coils; and (3) assessing neutron volumetric heating in structural components. Here, this work aims to perform a validation of the neutronics analysis approach by code-to-code comparison. State-of-the-art software stacks are employed for the neutronics analysis of the ARC reactor, and a computer-aided design (CAD) model is used directly for Monte Carlo (MC) neutron transport calculations. Three software stacks, Attila-MCNP, OpenMC-DAGMC, and Shift-DAGMC, are used to perform neutronic analyses of a 90° sector CAD model of the ARC reactor. In conclusion, results show that the flux tallies calculated by the three software stacks are very close. Also, the volumetric heating and tritium breeding values have less than 0.6% relative difference between codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation

Abstract Stand-off runaway electron termination by injected tungsten particulates offers a plausible option in the toolbox of disruption mitigation. Tungsten is an attractive material choice for this application due to large electron stopping power and high melting point. To assess the feasibility of this scheme, we simulate runaway collisions with tungsten particulates using the MCNP program for incident runaway energies ranging from 1 to 10 MeV. We assess runaway termination from energetics and collisional kinematics perspectives. Energetically, the simulations show that 99% of runaway beam energy is removed by tungsten particulates on a timescale of 4–9 µ s. Kinematically, the simulations show that 99% of runaways are terminated by absorption or backscattering on a timescale of 3–4 µ s. By either metric, the runaway beam is effectively terminated before the onset of particulate melting. Furthermore, the simulations show that secondary radiation emission by tungsten particulates does not significantly impact the runaway termination efficacy of this scheme. Secondary radiation is emitted at lower particle energies than the incident runaways and with a broad angular distribution such that the majority of secondary electrons emitted will not experience efficient runaway re-acceleration. Overall, the stand-off runaway termination scheme is a promising concept as a last line of defense against runaway damage in ITER, SPARC, and other future burning-plasma tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Neutronics Analyses for the Spallation Neutron Source Second Target Station

The Spallation Neutron Source (SNS) and the High Flux Isotope Reactor (HFIR), two of the word-class neutron scattering facilities, are located at the Oak Ridge National Laboratory. The SNS and HFIR are funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Science, and are operated as user facilities, available to researchers from all over the world. Currently there are thirteen neutron scattering instruments in operation at the HFIR and twenty at the SNS First Target Station.Producing the cutting edge science requires continuous improvements and development of the facilities and instruments. The SNS was designed from the beginning to allow addition of the Second Target Station (STS), and an upgrade of the accelerator power. At this time both advances: the accelerator upgrade and the construction of the STS are in preparation.Current baseline design for the STS calls for a rotating compact tungsten target. The target station is driven with short (<1 μs) proton pulses at 10 Hz repetition rate and 467 kW proton beam power, and is optimized for high intensity and high resolution long wavelength neutron applications. The proton beam footprint as small as acceptable from the mechanical and heat removal aspects is planned to generate a compact-volume neutron production zone in the target, which is essential for tight coupling of the target and the moderators and for achieving high-intensity peak thermal and cold neutron fluxes. The STS will allow operation of approximately 22 beamlines and will expand and complement the current national neutron scattering capabilities.This paper will present an update of the status of the neutronics analyses performed for the STS and will discuss the performance of the moderators, heating rates, activation analyses, and shielding calculations which provide input in the engineering design.

Remec, Igor↗
DOE OSTI2016

Reducing statistical uncertainties in simulated organ doses of phantoms immersed in water

In this study, methods are addressed to reduce the computational time to compute organ-dose rate coefficients using Monte Carlo techniques. Several variance reduction techniques are compared including the reciprocity method, importance sampling, weight windows and the use of the ADVANTG software package. For low-energy photons, the runtime was reduced by a factor of 10 5 when using the reciprocity method for kerma computation for immersion of a phantom in contaminated water. This is particularly significant since impractically long simulation times are required to achieve reasonable statistical uncertainties in organ dose for low-energy photons in this source medium and geometry. Although the MCNP Monte Carlo code is used in this paper, the reciprocity technique can be used equally well with other Monte Carlo codes.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

Organ and effective dose coefficients for cranial and caudal irradiation geometries: Neutrons

Dose coefficients based on the recommendations of International Commission on Radiological Protection (ICRP) Publication 103 were reported in ICRP Publication 116, the revision of ICRP Publication 74 and ICRU Publication 57 for the six reference irradiation geometries: anterior–posterior, posterior–anterior, right and left lateral, rotational and isotropic. In this work, dose coefficients for neutron irradiation of the body with parallel beams directed upward from below the feet (caudal) and downward from above the head (cranial) using the ICRP 103 methodology were computed using the MCNP 6.1 radiation transport code. The dose coefficients were determined for neutrons ranging in energy from 10 –9 MeV to 10 GeV. Here, at energies below about 500 MeV, the cranial and caudal dose coefficients are less than those for the six reference geometries reported in ICRP Publication 116.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

Implementation of a Portable HPGe for Field Contamination Assay

In this work, using MCNP to construct a detector model based initially on x-ray images of a portable high purity germanium (HPGe) detector followed by normalizing covering material values to also agree with check source responses, a validation of the model was attained. By calibrating the detector parameters using large count spectra, rigorous reproducibility is attained for high activity measurements but does not prevent deviations from normality in error distributions at the very low count events where spectral peaks are not always identifiable. The resulting model was created to allow operational assay of contamination over large areal distributions that could not otherwise be measured, such as the exhaust shaft at the Waste Isolation Pilot Plant (WIPP). Results indicate that contamination levels of activity in the exhaust shaft can be assayed to within a factor of 2. Lastly, detection limits are evaluated to be well below the contamination levels, which would constitute a legal environmental release if unfiltered ventilation of the underground facility were used.

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

Alpha Air Sample Counting Efficiency Versus Dust Loading: Evaluation of a Large Data Set

Dust loading on air sample filters is known to cause a loss of efficiency for direct counting of alpha activity on the filters, but the amount of dust loading and the correction factor needed to account for attenuated alpha particles is difficult to assess. In this paper, correction factors are developed by statistical analysis of a large database of air sample results for a uranium and plutonium processing facility at the Savannah River Site. As is typically the case, dust-loading data is not directly available, but sample volume is found to be a reasonable proxy measure; the amount of dust loading is inferred by a combination of the derived correction factors and a Monte Carlo model. The technique compares the distribution of activity ratios [beta/(beta + alpha)] by volume and applies a range of correction factors on the raw alpha count rate. The best-fit results with this method are compared with MCNP modeling of activity uniformly deposited in the dust and analytical laboratory results of digested filters. Finally, a linear fit is proposed to evenly-deposited alpha activity collected on filters with dust loading over a range of about 2 mg cm -2 to 1,000 mg cm -2 .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Trabecula: A Random Generator of Computational Phantoms for Bone Marrow Dosimetry

This study was motivated by the efforts on evaluation of radiation risk for leukemia incidence in the Techa River Cohort, where the main bone marrow dose contributors were 89,90Sr (bone-seeking beta emitters). Energy deposition in bone marrow targets was evaluated by simulating radiation particle transport using computational phantoms. The present paper describes the computer program “Trabecula” implementing an algorithm for parametric generation of computational phantoms, which serve as the basis for calculating bone marrow doses. “Trabecula” is a user-friendly tool that automatically converts analytical models into voxelized representations that are directly compatible as input to MCNP.

Shishkina, Elena A.↗
DOE OSTI2022

Applying ALARA Principles in the Design of New Radiological Facilities

The application of ALARA (As Low As Reasonably Achievable) principles to the design of new radiological facilities at Argonne National Laboratory provides a consistent radiation safety basis for future facility operations. The Department of Energy Rule 10 CFR 835 specifies design objectives to be adopted during the design of new facilities for controlling personal radiation exposure. One is to keep exposure levels below 20% of the applicable standards in 10 CFR 835.202. For a radiation worker 20% of the standard corresponds to 5 µSv hr -1 for a 2,000 hour work year. For a member of the public the DOE design objective corresponds to 0.2 mSv in a calendar year. For a new facility worst case radionuclides and their source strengths are chosen. Local shielding is specified to reduce exposure rates to less than 50 µSv hr -1 at 30 cm from the shielding. The current version MCNP6 of the Los Alamos radiation shielding computer program MCNP is then used to calculate the exposure rates elsewhere. Design modifications are made to meet the criteria. The calculations and resulting facility design modifications are discussed for two new radiological facilities.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Neutron-Neutron Correlations in the Photofission of 238U

Background: In the fission of actinides, the nearly back-to-back motion of the fission fragments has a strong effect on the kinematics of fission neutrons. This leads to a favoring of opening angles near 0? and 180? in the neutron-neutron (n-n) opening angle distributions of correlated neutron pairs from the same fission event. Purpose: To measure the n-n opening angle and energy correlations in the photofission of 238U. As of this writing, measurements of correlated n-n opening angle distributions have been reported only for the spontaneous and neutron-induced fission of actinides. This work is the first to report such a measurement using photofission and will provide useful experimental input for photofission models used in codes such as MCNP and FREYA. Method: Fission is induced using bremsstrahlung photons produced via a low duty factor, pulsed, linear electron accelerator. The bremsstrahlung photon beam impinges upon a 238U target that is surrounded by a large neutron scintillation detection system capable of measuring particle position and time of flight, from which n-n opening angle and energy are measured. Neutron-neutron angular correlations are determined by taking the ratio between a correlated neutron distribution and an uncorrelated neutron distribution formed by the pairing of neutrons produced during different beam pulses. This analysis technique greatly diminishes effects due to detector efficiencies, acceptance, and experimental drifts. Results: The angular correlation of neutrons from the photofission of 238U shows a high dependence on neutron energy as well as a dependence on the angle of the emitted neutrons with respect to the incoming photon beam. Angular correlations were also measured using neutrons from the spontaneous fission of 252Cf, showing good agreement with past measurements. Conclusions: The measured angular correlations reflect the underlying back-to-back nature of the fission fragments. An anomalous decline in n-n yield was observed for opening angles near 180? for 238U.

Burggraf, Jeff↗
DOE OSTI2020

Cross section of neutrons from the H 2 ( n , 2 n ) reaction at E n = 15 MeV

In this work, the double-differential cross section of the deuteron breakup reaction 2 H(n, 2n) has been studied experimentally with a neutron beam energy of 15 MeV. Special attention has been devoted to estimation of background condition and multiple scattering effect in the scattering sample. Experimental data have been compared with models based on phase-space approximation used in the ENDF/B-VIII.0 data library and in the MCNP neutron transport code, as well as with rigorous model based on Faddeev equations used for cross section evaluations in JENDL data library. It was found that experimental data are better reproduced by Faddeev model, however, the model overestimates data in the low-energy region of the neutron spectrum (<4 MeV)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗