Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo code”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Development and application of a reverse Monte Carlo radiative transfer code for rocket plume base heating

A reverse Monte Carlo radiative transfer code to predict rocket plume base heating is presented. In this technique rays representing the radiation propagation are traced backwards in time from the receiving surface to the point of emission in the plume. This increases the computational efficiency relative to the forward Monte Carlo technique when calculating the radiation reaching a specific point, as only the rays that strike the receiving point are considered.

Everson, John↗
DOE OSTI2020

Electron transport in gaseous detectors with a Python-based Monte Carlo simulation code

Understanding electron drift and diffusion in gases and gas mixtures is a topic of central importance for the development of modern particle detection instrumentation. The industry-standard MagBoltz code has become an invaluable tool during its 20 years of development, providing capability to solve for electron transport (‘swarm’) properties based on a growing encyclopedia of built-in collision cross sections. We have made a refactorization of this code from FORTRAN into Cython, and studied a range of gas mixtures of interest in high energy and nuclear physics. The results from the new open source PyBoltz package match the outputs from the original MagBoltz code, with comparable simulation speed. An extension to the capabilities of the original code is demonstrated, in implementation of a new Modified Effective Range Theory interface. We hope that the versatility afforded by the new Python code-base will encourage continued use and development of the MagBoltz tools by the particle physics community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

ITER ECH Transmission Line System Design and Status

The electron cyclotron (EC) heating & current drive (H&CD) system on ITER provides plasma heating by generating, transmitting, and launching high-intensity, high-frequency (170 GHz) electromagnetic wave energy steerable across the plasma cross-section. The transmission line (TL) subsystem connects the Matching Optics Unit (MOU) on each of the 24 gyrotrons to the 32 feed points in the four upper launchers and the 24 feed points in the equatorial launcher. Each TL must be able to operate at up to 1.2 MW of input power for up to 1 hour pulse lengths. The TL system contains 50 mm water-cooled corrugated waveguide, 90o miter bends, 140o miter bends, polarizer miter bend pairs, switches, expansion units, pumpouts, DC breaks, MOU-TL adapters, Radio Frequency (RF) loads, and isolation shutter valves. A detailed finite element analysis has been used to verify the thermo-mechanical performance of each component. The microwave performance has been analyzed using a 2-D electromagnetic code combined with a Monte Carlo code. This approach allows the impact of manufacturing and installation tolerances to be assessed and optimized to provide a high probability of achieving the system performance requirements. Prototypes of the waveguide and TL components have been fabricated and tested at high-power. Production contracts are now being issued for fabrication and delivery of the waveguide and components to ITER.

Hanson, Gregory↗
DOE OSTI2022

Reactor physics benchmark of Westinghouse PWR core design suite for high burnup/high enrichment fuel - Part 1: Pin and lattice

Westinghouse has performed a comprehensive set of code-to-code benchmark comparisons to corroborate application of its state-of-the-art NEXUS-based code suite to high burnup and high enrichment (HB/HE) fueled PWR core design and reload analysis. This paper focuses on analysis of the predictions for pin and lattice designs representative of HB/HE core reloads from the Westinghouse advanced lattice code PARAGON2, the Monte-Carlo code SERPENT2 and MPACT, the deterministic neutronic component of VERA. A companion paper presented at this conference analyzes results obtained with the Westinghouse ANC9 core simulator, using lattice data provided by PARAGON2 with the NEXUS cross-section methodology, and MPACT for 2D and 3D core configurations incorporating HB/HE fuel reloads to achieve cycle length extension to 24-month in high-power density PWRs. The results of the pin and lattice benchmark presented in this paper show that the neutronic predictions from PARAGON2 for HB/HE fueled PWRs are in remarkable agreement with predictions from both SERPENT2 and MPACT including for depletion to high fuel burnup values, and employing a set of challenging lattices which feature combined use of burnable absorbers to obtain the required reactivity hold-down to support PWR transition to HB/HE fuel and 24-month operational cycles. As corroborated by this analysis and the analysis contained in the companion paper, the Westinghouse NEXUS-based PWR core analysis package can be confidently applied by industry for design analysis of PWR core reloads adopting HB/HE fuel and properly support the industry commercialization plans for in-reactor deployment. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

(U) LNK3DNT Capability in SENSMG

The SENSMG multigroup neutron sensitivity code has been extended to read materials and geometry from redoin and lnk3dnt files written by the PARTISN multigroup discrete ordinates code or the MCNP6 Monte Carlo code. This report describes the user interface for this capability and the code verification that was done. Not all of the capability of SENSMG is available for redoin/lnk3dnt inputs. The only responses allowed are k eff or α eigenvalues. Fixed-source problems are not allowed, and there is no mechanism to input reaction rates or reaction-rate ratios. The next section of this report describes the user interface. Section III describes the verification of the capability. Section IV is a summary and conclusions. The input files used for verification are listed in the appendix.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

(U) LNK3DNT Capability in SENSMG Extended

The SENSMG multigroup neutron sensitivity code can read materials and geometry from redoin and lnk3dnt files written by the PARTISN multigroup discrete ordinates code or the MCNP6 Monte Carlo code. This report describes the user interface for this capability and the code verification that was done. Reference 5 described an initial version of this capability. Two important extensions have been made: 1) Reaction-rate ratios are now allowed as responses with redoin/lnk3dnt inputs; and 2) Sensitivities of reaction-rate ratios to mass densities and interface locations are now calculated (for all problems, not just redoin/lnk3dnt problems).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

(U) LNK3DNT Capability in SENSMG Tested with Transport Correction

The SENSMG multigroup neutron sensitivity code can read materials and geometry from redoin and lnk3dnt files written by the PARTISN multigroup discrete ordinates code or the MCNP6 Monte Carlo code. Reference 5 described an initial version of this capability. Reference 6 introduced two important extensions: 1) Reaction-rate ratios are now allowed as responses with redoin/lnk3dnt inputs; and 2) Sensitivities of reaction-rate ratios to mass densities and interface locations are now calculated (for all problems, not just redoin/lnk3dnt problems). In this report, we document the testing of the redoin/lnk3dnt capability when used with transport correction. PARTISN’s transport corrections are described in Ref. 3, and Refs. 7 and 8 discuss the calculation of cross-section sensitivities when transport corrections are applied. References 7 and 8 verify SENSMG for problems with transport correction using central differences. This report verifies the capability with redoin/lnk3dnt inputs by comparison with regular SENSMG inputs. Unlike Refs. 5 and 6, in this report we use a spherical test problem because the sensitivity of the response with respect to the highest-order scattering is zero in cylinders due to the use of flux moments in the inner products.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Modeling the Ex-core Detector Response into a Production Option in MPACT

In the past few years, the Consortium for Advanced Simulation of Light Water Reactors (CASL) has overseen a successful activity to extend the radiation transport capability from the reactor in-core geometry to the ex-core geometry by coupling the deterministic neutronics code MPACT, with the Monte Carlo code Shift. However, transient problems such as RIAs require an efficient prediction of the ex-core detector response to determine the timeframe for a reactor trip. Therefore, in addition to the high-fidelity MPACT-Shift coupled calculation, separate efforts have been made to explore simplified methods for ex-core detector response with MPACT standalone calculation. In a previous milestone L3_RTM.MCH.P17.06, several simplified methods for ex-core detector response were developed in MPACT. In particular, the double kernel method is based on the default MPACT core geometry, so extensive modeling or transport calculation with full reflector and ex-core details are not needed in this method. Since the 1-D diffusion and neutron streaming kernels are used for extended ex-core calculations, the extra computational efforts are very small. The work implemented in the current milestone provides a fully integrated option for MPACT to approximately calculate the ex-core detector response using the double kernel method. To overcome the difficulty of obtaining the consistent diffusion coefficients (these were obtained through separate CMFD calculation and edits), the 1-D diffusion solver is replaced by a 1-D Sn transport solver to compute the neutron flux at the vessel outer surface. The new VERA inputs for ex-core detector are processed in MPACT. The detector signals are linked to the SCRAM logic as an option to trip the core during transient calculation. Two sets of 3-D quarter core problems are run with MPACT standalone and VeraShift to verify the radial and axial effects on detector response by varying the moderator density and control rod position. For the cases with moderator density changes, the detector responses from MPACT are within 2 sigma error of VeraShift. For the cases at different control rod positions, the axial offsets computed from MPACT are within 2% error of VeraShift results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1984

A highly optimized vectorized code for Monte Carlo simulations of SU(3) lattice gauge theories

New methods are introduced for improving the performance of the vectorized Monte Carlo SU(3) lattice gauge theory algorithm using the CDC CYBER 205. Structure, algorithm and programming considerations are discussed. The performance achieved for a 16(4) lattice on a 2-pipe system may be phrased in terms of the link update time or overall MFLOPS rates. For 32-bit arithmetic, it is 36.3 microsecond/link for 8 hits per iteration (40.9 microsecond for 10 hits) or 101.5 MFLOPS.

Barkai, D.↗
DOE OSTI2021

TNSL Overview

Thermal neutron scattering law (TNSL) data describe low-energy neutrons scattering off of bound materials, and can have a significant impact on modeling any system with slow neutrons, including nuclear reactors. Previous work to introduce TNSL data to neutron transport codes at LLNL focused on COG and TART [1], with the limitation that these codes require highly specialized data processing and formatting. We have recently increased efforts to process TNSL data with the central LLNL nuclear data processing code FUDGE, to be stored in the generalized nuclear database structure (GNDS) for use in any general transport code with the ability to read GNDS data. The first step in this effort is to verify that the TNSL processing with FUDGE yields results comparable to results obtained using the LANL nuclear data processing code NJOY. The next step is to verify the transport of thermal neutrons in Mercury (a Monte Carlo code) and Ardra (a deterministic code) against one another, as well as against the LANL Monte Carlo neutron transport code MCNP. This verification step has not been completed, due to a number of discrepancies between results obtained using differently processed data. There is ongoing effort to understand differences between FUDGE and NJOY. Finally, we map out our current capability to validate TNSL data against benchmark systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

The History of Monte Carlo and MCNP at Los Alamos [Slides]

The Monte Carlo method for radiation particle transport has its origins at LANL dating back to the 1940’s. The creators of these methods were Drs. Stanislaw Ulam, John von Neumann, Robert Richtmyer, and Nicholas Metropolis. Monte Carlo methods for particle transport have been driving computational developments since the beginning of modern computers; this continues today. In the 1950’s and 1960’s, these new methods were organized into a series of special-purpose Monte Carlo codes, including MCS, MCN, MCP, and MCG. These codes were able to transport neutrons and photons for specialized LANL applications. In 1977, these separate codes were combined to create the LANL Monte Carlo N-Particle (MCNP) radiation particle transport code. In 1983, MCNP3 was released for public distribution to the Radiation Safety Information Computational Center (RSICC). MCNP6.3 will be released late in 2021. Each year, LANL has ~100 MCNP new users and RSICC ~1250 new licenses distributed. This talk will review the Los Alamos history of the development of the modern Monte Carlo method and MCNP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗