Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo N Particle”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Monte Carlo N-Particle forward modeling for density reconstruction of double shell capsule radiographs

In the Double Shell Inertial Confinement Fusion concept, characterizing the shape asymmetry of imploding metal shells is vital for understanding energy-efficient compression and radiative losses of the thermonuclear fuel. The Monte Carlo N-Particle MCNP® code forward models radiography of Double Shell capsule implosions using the Advanced Radiographic Capability at the National Ignition Facility. A procedure is developed for using MCNP to reconstruct density profiles from the radiograph image intensity. For a given Double Shell imploding target geometry, MCNP radiographs predict image contrast, which can help guide experimental design. Finally, in future work, the calculated MCNP synthetic radiographs will be compared with experimental radiographs to determine the radial and azimuthal density profiles of the Double Shell capsules.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Monte Carlo N-Particle Tracking of Ultrafine Particle Flow in Bent Micro-Tubes

The problem of large pressure-differential driven laminar convective-diffusive ultrafine aerosol flow through bent micro-tubes is of interest in several contemporary research areas including; release of contents from pressurized containment vessels, aerosol sampling equipment, advanced scientific instruments, gas-phase micro-heat exchangers, and microfluidic devices. In each of these areas, the predominant problem is the determination of the fraction of particles entering the micro-tube that is deposited within the tube and the fraction that is transmitted through. Due to the extensive parameter restrictions of this class of problems, a Lagrangian particle tracking method making use of the coupling of the analytical stream line solutions of Dean and the simplified Langevin equation is quite a useful tool in problem characterization. This method is a direct analog to the Monte Carlo N-Particle method of particle transport extensively used in nuclear physics and engineering. In this work, 10 nm diameter particles with a density of 1 g/cm3 are tracked within micro-tubes with toroidal bends with pressure differentials ranging between 0.2175 and 0.87 atmospheres. The tubes have radii of 25 microns and 50 microns and the radius of curvature is between 1 m and 0.3183 cm. The carrier gas is helium, and temperatures of 298 K and 558 K are considered. Numerical convergence is considered as a function of time step size and of the number of particles per simulation. Particle transmission rates and deposition patterns within the bent micro-tubes are calculated.

Casella, Andrew M.↗
DOE OSTI2025

Monte Carlo N-Particle Transport Performance of Predicting Digital Radiographic IQI Inspection

The identification of porosity, geometric noncompliance, and other defect types are critical to the qualification of materials and components. X-ray radiographic nondestructive testing is a common industrial inspection method for process quality control and component qualification and certification. Digital radiography provides a quick and efficient alternative when compared to traditional film-based inspection. The quality of radiographic inspection is dependent on equipment specifications, such as the source spot size and detector pixel size, and the specific parameters selected for use for the radiographic technique. To evaluate if an x-ray system and technique is sufficient for a given requirement, a radiographic image quality indicator (IQI) can be used. Radiographic IQIs in hard to machine materials or hard to manufacture defects can be time consuming and expensive to manufacture. This study was conducted to evaluate current Savannah River National Laboratory (SRNL) x-ray imaging systems with a custom tantalum IQI and using Monte Carlo simulations to predict the performance of future systems. The tantalum IQI was tested using a Siefert Isovolt 420 keV x-ray tube with a Perkin Elmer XRD 1611 flat panel with 100-micron pixels. Using the Monte Carlo N-Particle transport software, the radiographic tally was used to simulate the photon flux through an identical tantalum IQI. These simulations provided a benchmark as to the best theoretical identification on a given system using our tantalum IQI. The simulations were refined to match SRNL’s current systems’ noise levels, leading to confidence in their ability to predict the performance of other systems that may be purchased and deployed in the future at the Savannah River Site. Future studies will be conducted to prove this research can be extended to artificially evaluate the ability for systems to identify critical defect sizes through x-ray radiographic inspection, drastically reducing the cost and time burdens of producing high-fidelity radiographic test articles.

digital X-ray radiography↗
DOE OSTI2024

Montepy

Montepy is an object-oriented Python interface for MCNP input files. It excels at reading, editing, creating, and writing MCNP input files. The Monte Carlo N Particle (MCNP) code is a Monte Carlo transport solver for simulating many different types of radiation problems and is developed by Los Alamos National Laboratory. Montepy creates an object-oriented interface for the MCNP input file to allow for efficient automated modifications to large, complex simulations.

Gale, Micah [Idaho National Laboratory (INL), Idah↗
DOE OSTI2010

An introduction to TR-X: a simplified tool for standardized analysis

TR-X is a multi-platform program that provides a graphical interface to Monte Carlo N-Particle transport (MCNP) and Monte Carlo N-Particle transport eXtended (MCNPX) codes. Included in this interface are tools to reduce the tedium of input file creation, provide standardization of model creation and analysis, and expedite the execution of the created models. TR-X provides tools to make the rapid testing of multiple permutations of these models easier, while also building in standardization that allows multiple solutions to be compared.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Shielding Benchmark Comparison - MCNP6.2

This report documents the calculations performed for a shielding code comparison between various Department of Energy (DOE) sites. These shielding calculations were performed using an experiment drawn from the International Handbook of Evaluated Criticality Safety Benchmark Experiments, published by the Organisation for Economic Cooperation and Development/Nuclear Energy Agency (OECD/NEA). The benchmark selected for comparison is ALARM-CF-AIR-LAB-001 (“Neutron Fields in the Three-Section Concrete Labyrinth from Cf-252 Source, Benchmark ALARM CF AIR LAB-001”). The Y-12 submission for this code comparison was performed with Monte Carlo N Particle (MCNP) Transport Code System, Version 6.2 and Automated Variance Reduction Generator (ADVANTG).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

Initial Neutronics Analyses for HEU to LEU Fuel Conversion of the Transient Reactor Test Facility (TREAT) at the Idaho National Laboratory

The purpose of the TREAT reactor is to generate large transient neutron pulses in test samples without over-heating the core to simulate fuel assembly accident conditions. The power transients in the present HEU core are inherently self-limiting such that the core prevents itself from overheating even in the event of a reactivity insertion accident. The objective of this study was to support the assessment of the feasibility of the TREAT core conversion based on the present reactor performance metrics and the technical specifications of the HEU core. The LEU fuel assembly studied had the same overall design, materials (UO 2 particles finely dispersed in graphite) and impurities content as the HEU fuel assembly. The Monte Carlo N–Particle code (MCNP) and the point kinetics code TREKIN were used in the analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Uncertainty Analysis for a Criticality Benchmark [Slides]

This presentation discusses how KRUSTY was analyzed with LANL’s MCNP neutron transport code, Monte Carlo N-Particle. In the Monte Carlo method, neutrons are born in fission events and fly around the geometry having events (e.g., scatter, capture, fission, escape) according to the natural probabilities. “Natural probabilities” are given as neutron cross sections , which are measured and evaluated nuclear data. One of the goals of the benchmark project is to improve the nuclear cross sections (and the codes that use them). In this work, the response of interest is called k eff , which describes the neutron multiplication in a near critical system. Like most quantities, k eff is not measured directly; it is only inferred. In summation, MCNP results for k eff matched the measurements extremely well. The evaluated uncertainties are asymmetric and smaller than 0.00100 (e.g., + 0.00080 / – 0.00062). The largest uncertainties are due to alignment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Uncertainty Analysis for a Criticality Benchmark [Slides]

This presentation discusses how KRUSTY was analyzed with LANL’s MCNP neutron transport code, Monte Carlo N-Particle. In the Monte Carlo method, neutrons are born in fission events and fly around the geometry having events (e.g., scatter, capture, fission, escape) according to the natural probabilities. “Natural probabilities” are given as neutron cross sections , which are measured and evaluated nuclear data. One of the goals of the benchmark project is to improve the nuclear cross sections (and the codes that use them). In this work, the response of interest is called k eff , which describes the neutron multiplication in a near critical system. Like most quantities, k eff is not measured directly; it is only inferred. In summation, MCNP results for k eff matched the measurements extremely well. The evaluated uncertainties are asymmetric and smaller than 0.00100 (e.g., + 0.00080 / – 0.00062). The largest uncertainties are due to alignment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2012

Verification of analytic energy moments for the one-dimensional energy dependent neutron diffusion equation with MCNP5 and Attila-7.1.0

The energy dependent neutron diffusion equation (EDNDE) is converted into a moment equation which is solved analytically for the 1-D problem of a bare sphere of pure 235U. The normalized moments 0–5 generated analytically are compared to normalized energy moments, from Monte Carlo N Particle 5 version 1.40 (MCNP5) and Attila-7.1.0-beta version (Attila). The analytic normalized neutron energy moments, fall between the results from MCNP5 (lower bound) and Attila (upper bound) and are accurate compared to MCNP5 neutron energy moments when error in this Monte Carlo simulation are considered. The error range is from 0% to 14%. The Attila moments are less accurate when compared to MCNP5 than the analytical moments derived in this work. The method of moments is shown to be a fast reliable method, compared to either Monte Carlo methods (MCNP5) or 30 multi-energy group methods (Attila).

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Comparison of Attila Transport Code with MCNP in the Analysis of Small Modular Reactor Shielding

Of great importance in the design of a nuclear power plant is the evaluation of shielding from the ionizing radiation produced in the core of the nuclear reactor. This requires simulation software to predict the dose rates in the areas of interest. These problems are referred to as deep penetration problems, as they involve the reduction of a radiation field by two or three orders of magnitude or more. This makes accurate analyses difficult for Monte Carlo based solutions that rely on enough particles being sampled in an area of interest. Comparing the simulation of the same modeled system with more than one computer application, especially one that uses a different algorithmic approach, is a way to gain confidence in the results. In this research, two models were created, one in Monte Carlo N Particle (MCNP) which uses a stochastic algorithmic approach, and the other in Attila, which uses a deterministic algorithmic approach to test the validity of the model and to estimate the dose rates in the area of interest. Successful comparison of the two models is a way we can validate the model used for that estimation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2012

Validation of energy moments from the one-dimensional energy dependent neutron diffusion equation, MCNP5 and Attila-7.1.0 with the GODIVA experiment

Normalized neutron energy moments (moments) from the one-dimensional energy dependent neutron diffusion equation (EDNDE), Monte Carlo N Particle 5 version 1.40 (MCNP5) and Attila-7.1.0-beta version (Attila) are validated with the GODIVA experiment (GODIVA). Energy moments 0–5 for all three methods are compared to GODIVA moments. GODIVA moments are measured with two methods. The 1st method is a time of flight (T-O-F) measurement of the average energy (moment 1) of the leaking neutrons from the surface of GODIVA and the 2nd method is from back calculating moments from foil activation analysis of various metal foils at the center of GODIVA. The error range of the EDNDE normalized moments compared to GODIVA is from 0% to 24%. The MCNP5 error range compared to GODIVA is 0–12% and the Attila error range is 0–79%. The method of moments is shown to be a fast reliable method, compared to either Monte Carlo methods (MCNP5) or 30 multi-energy group methods (Attila) with regard to the GODIVA experiment.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2013

Verification and validation of the maximum entropy method for reconstructing neutron flux, with MCNP5, Attila-7.1.0 and the GODIVA experiment

Verification and validation of reconstructed neutron flux based on the maximum entropy method is presented in this paper. The verification is carried out by comparing the neutron flux spectrum from the maximum entropy method with Monte Carlo N Particle 5 version 1.40 (MCNP5) and Attila-7.1.0-beta (Attila). A spherical 100% 235U critical assembly is modeled as the test case to compare the three methods. The verification error range for the maximum entropy method is 15–21% where MCNP5 is taken to be the comparison standard. Attila relative error for the critical assembly is 20–35%. Validation is accomplished by comparing a neutron flux spectrum that is back calculated from foil activation measurements performed in the GODIVA experiment (GODIVA). The error range of the reconstructed flux compared to GODIVA is 0–10%. The error range of the neutron flux spectrum from MCNP5 compared to GODIVA is 0–20% and the Attila error range compared to the GODIVA is 0–35%. The maximum entropy method is shown to be a fast reliable method, compared to either Monte Carlo methods (MCNP5) or 30 multienergy group methods (Attila) and with respect to the GODIVA experiment.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2014

SU-E-J-100: Reconstruction of Prompt Gamma Ray Three Dimensional SPECT Image From Boron Neutron Capture Therapy(BNCT)

Purpose: Purpose of paper is to confirm the feasibility of acquisition of three dimensional single photon emission computed tomography (SPECT) image from boron neutron capture therapy (BNCT) using Monte Carlo simulation. Methods: In case of simulation, the pixelated SPECT detector, collimator and phantom were simulated using Monte Carlo n particle extended (MCNPX) simulation tool. A thermal neutron source (<1 eV) was used to react with the boron uptake region (BUR) in the phantom. Each geometry had a spherical pattern, and three different BURs (A, B and C region, density: 2.08 g/cm3) were located in the middle of the brain phantom. The data from 128 projections for each sorting process were used to achieve image reconstruction. The ordered subset expectation maximization (OSEM) reconstruction algorithm was used to obtain a tomographic image with eight subsets and five iterations. The receiver operating characteristic (ROC) curve analysis was used to evaluate the geometric accuracy of reconstructed image. Results: The OSEM image was compared with the original phantom pattern image. The area under the curve (AUC) was calculated as the gross area under each ROC curve. The three calculated AUC values were 0.738 (A region), 0.623 (B region), and 0.817 (C region). The differences between length of centers of two boron regions and distance of maximum count points were 0.3 cm, 1.6 cm and 1.4 cm. Conclusion: The possibility of extracting a 3D BNCT SPECT image was confirmed using the Monte Carlo simulation and OSEM algorithm. The prospects for obtaining an actual BNCT SPECT image were estimated from the quality of the simulated image and the simulation conditions. When multiple tumor region should be treated using the BNCT, a reasonable model to determine how many useful images can be obtained from the SPECT could be provided to the BNCT facilities. This research was supported by the Leading Foreign Research Institute Recruitment Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, Information and Communication Technologies (ICT) and Future Planning (MSIP)(Grant No.200900420) and the Radiation Technology Research and Development program (Grant No.2013043498), Republic of Korea.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2013

Characterization of Filters Loaded With Reactor Strontium Carbonate - 13203

A collection of three highly radioactive filters containing reactor strontium carbonate were being prepared for disposal. All three filters were approximately characterized at the time of manufacture by gravimetric methods. The first filter had been partially emptied, and the quantity of residual activity was uncertain. Dose rate to activity modeling using the Monte-Carlo N Particle (MCNP) code was selected to confirm the gravimetric characterization of the full filters, and to fully characterize the partially emptied filter. Although dose rate to activity modeling using MCNP is a common technique, it is not often used for Bremsstrahlung-dominant materials such as reactor strontium. As a result, different MCNP modeling options were compared to determine the optimum approach. This comparison indicated that the accuracy of the results were heavily dependent on the MCNP modeling details and the location of the dose rate measurement point. The optimum model utilized a photon spectrum generated by the Oak Ridge Isotope Generation and Depletion (ORIGEN) code and dose rates measured at 30 cm. Results from the optimum model agreed with the gravimetric estimates within 15%. It was demonstrated that dose rate to activity modeling can be successful for Bremsstrahlung-dominant radioactive materials. However, the degree of success is heavily dependent on the choice of modeling techniques. (authors)

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

SU-C-16A-03: Direction Modulated Brachytherapy for HDR Treatment of Cervical Cancer

Purpose: To investigate a new Directional Modulated Brachytherapy (DMBT) intra-uterine tandem using various 192-Ir after-loaders. Methods: Dose distributions from the 192-Ir sources were modulated using a 6.3mm diameter tungsten shield (18.0g/cm3). The source moved along 6 longitudinal grooves, each 1.3mm in diameter, evenly spaced along periphery of the shield, The tungsten rod was enclosqed by 0.5mm thick Delrin (1.41g/cc). Monte Carlo N particle (MCNPX) was used to calculate dose distributions. 51million particles were calculated on 504 cores of a supercomputer. Fifteen different patients originally treated with a traditional tandem-and-ovoid applicator, with 5 fractions each, (15 patients X 5 fxs = 75 plans) were re-planned with the DMBT applicator combined with traditional ovoids, on an in-house developed HDR brachytherapy planning platform, which used intensity modulated planning capabilities using a constrained gradient optimization algorithm. For all plans the prescription dose was 6 Gy and they were normalized to match the clinical treated V100. Results: Generally, the DMBT plan quality was a remarkable improvement from conventional T and O plans because of the anisotropic dose distribution of DMBT. The largest difference was to the bladder which had a 0.59±0.87 Gy (8.5±28.7%) reduction in dose. This was because of the the horseshoe shape (U-shape) of the bladder. The dose reduction to rectum and sigmoid were 0.48±0.55 Gy (21.1±27.2%) and 0.10±0.38 Gy (40.6±214.9%), respectively. The D90 to the HRCTV was 6.55±0.96 Gy (conventional T and O) and 6.59±1.06 Gy (DMBT). Conclusion: For image guided adaptive brachytherapy, greater flexibility of radiation intensity is essential and DMBT can be the solution.

60 APPLIED LIFE SCIENCES↗