Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo codes”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Application of MORSE to radiation analysis of nuclear flight propulsion modules

Several modifications and additions were made to the multigroup Monte Carlo code (MORSE) to implement its use in a computational procedure for performing radiation analyses of NERVA nuclear flight propulsion modules. These changes include the incorporation of a new general geometry module; the inclusion of an expectation tracklength estimator; and the option to obtain source information from two-dimensional discrete ordinates calculations. Computations comparing MORSE and a point cross section Monte Carlo code, COHORT, were made in which a coupled discrete ordinates/Monte Carlo procedure was used to calculate the gamma dose rate at tank top locations of a typical propulsion module. The dose rates obtained from the MORSE computation agreed with the dose rates obtained from the COHORT computation to within the limits of the statistical accuracy of the calculations.

Woolson, W. A.↗
NASA NTRS1997

CGRO Guest Investigator Program

The following are highlights from the research supported by this grant: (1) Theory of gamma-ray blazars: We studied the theory of gamma-ray blazars, being among the first investigators to propose that the GeV emission arises from Comptonization of diffuse radiation surrounding the jet, rather than from the synchrotron-self-Compton mechanism. In related work, we uncovered possible connections between the mechanisms of gamma-ray blazars and those of intraday radio variability, and have conducted a general study of the role of Compton radiation drag on the dynamics of relativistic jets. (2) A Nonlinear Monte Carlo code for gamma-ray spectrum formation: We developed, tested, and applied the first Nonlinear Monte Carlo (NLMC) code for simulating gamma-ray production and transfer under much more general (and realistic) conditions than are accessible with other techniques. The present version of the code is designed to simulate conditions thought to be present in active galactic nuclei and certain types of X-ray binaries, and includes the physics needed to model thermal and nonthermal electron-positron pair cascades. Unlike traditional Monte-Carlo techniques, our method can accurately handle highly non-linear systems in which the radiation and particle backgrounds must be determined self-consistently and in which the particle energies span many orders of magnitude. Unlike models based on kinetic equations, our code can handle arbitrary source geometries and relativistic kinematic effects In its first important application following testing, we showed that popular semi-analytic accretion disk corona models for Seyfert spectra are seriously in error, and demonstrated how the spectra can be simulated if the disk is sparsely covered by localized 'flares'.

Begelman, Mitchell C.↗
DOE OSTI2025

Initial steady-state core simulation capability or thermal and pool-type molten salt reactors, coupling reactor physics, thermal-hydraulics, and evolving chemistry

This report presents the development and validation of an initial steady-state multiphysics capability for molten salt reactors (MSRs) under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in Fiscal Year 2025. The framework integrates neutronics, thermal-hydraulics, species transport, and thermochemistry to capture the coupled dynamics of liquid-fueled systems. Implementation and testing were performed on two representative designs: the Molten Salt Reactor Experiment (MSRE), a thermal-spectrum, channeled-flow reactor, and the Lotus Molten Salt Reactor (L-MSR), a fast-spectrum, pool-type reactor. The modeling suite employs Griffin for reactor physics and depletion, Pronghorn and SAM for thermal-hydraulics, Thermochimica for chemistry, and Saline for thermophysical properties, with benchmarking and validation carried out against historical MSRE data, experimental flow-loop measurements, and reference depletion calculations from Monte Carlo codes. The framework demonstrated the ability to reproduce key reactor behaviors including temperature feedback, reactivity losses, delayed neutron precursor transport, xenon poisoning, and redox potential evolution. The results confirm the feasibility and accuracy of the coupled models in predicting steady-state and selected transient MSR behaviors. This latter ones are used in this report as a proxy indicating that the steady-state models from which the transient starts are accurate. For MSRE, validation showed good agreement with pump start-up and natural circulation tests, while for the L-MSR, benchmarking confirmed hydraulic calibration and consistency of neutronics–thermal coupling. The tools also provided new insights into species transport, noble metal deposition, and salt solidification dynamics. On the Xenon transport front, the code is validated against the steady state Xenon poisoining measurement and showed good agreement with the experimental value. Identified areas for future work include advanced void transport modeling, three-dimensional simulations, improved alloy corrosion models, and tighter integration with high-fidelity Monte Carlo codes. These developments provide a foundation for high-fidelity MSR simulations that can support reactor design optimization, safety assessments, and long-term operational strategies.

42 - ENGINEERING↗
NASA NTRS1972

Common radiation analysis model for 75,000 pound thrust NERVA engine (1137400E)

The mathematical model and sources of radiation used for the radiation analysis and shielding activities in support of the design of the 1137400E version of the 75,000 lbs thrust NERVA engine are presented. The nuclear subsystem (NSS) and non-nuclear components are discussed. The geometrical model for the NSS is two dimensional as required for the DOT discrete ordinates computer code or for an azimuthally symetrical three dimensional Point Kernel or Monte Carlo code. The geometrical model for the non-nuclear components is three dimensional in the FASTER geometry format. This geometry routine is inherent in the ANSC versions of the QAD and GGG Point Kernal programs and the COHORT Monte Carlo program. Data are included pertaining to a pressure vessel surface radiation source data tape which has been used as the basis for starting ANSC analyses with the DASH code to bridge into the COHORT Monte Carlo code using the WANL supplied DOT angular flux leakage data. In addition to the model descriptions and sources of radiation, the methods of analyses are briefly described.

Warman, E. A.↗
DOE OSTI2020

NRC Multiphysics Analysis Capability Deployment FY2020 - Part 4

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Sup- port for Advanced Non-Light Water Reactors.” The tasks completed for this report are: • Task 3a: This task demonstrates the effectiveness of Griffin’s methods for interpolating cross sections and SPH correction factors as a function of control drum rotation. The Monte Carlo code Serpent is used to generate reference results. The accuracy of Griffin-calculated eigenvalues and reaction rates are studied as a function of control drum rotation for cross section and SPH factor libraries of varying fidelity. Higher-fidelity libraries require less interpolation and are therefore inherently more accurate; however, they require more computational resources to create. Overall, Griffin is able to accurately model eigenvalue and reaction rates over the full range of drum rotation using a library generated from three to five discrete drum rotations for a particular microreactor design. • Task 3c: This preliminary task successfully demonstrates the performance of a coupled neutron-photon transport calculation with Griffin using the spherical harmonics approximation. The Griffin solutions are compared to reference solutions from the continuous energy Monte Carlo codes MCNP and Serpent. Furthermore, this work also demonstrates how to prepare the neutron and photon libraries using NJOY, MCNP, and Serpent with the latest ENDF/B-VIII data for 235 U and graphite. The Griffin total energy deposition in the active core region is calculated within 0.3% of the reference, but it degrades in the reflector regions with a maximum difference of -1.5%. Overall, the results look promising and future analysis should use an improved data stream and focus on determining the potential cancellation of error within each energy group.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Status of GPU capabilities within the Shift Monte Carlo radiation transport code

Shift is a general-purpose Monte Carlo (MC) radiation transport code for fission, fusion, and national security applications. Shift has been adapted to efficiently run on GPUs in order to leverage leadership-class supercomputers. This work presents Shift’s current GPU capabilities. These include core radiation transport capabilities for eigenvalue and fixed-source simulations, and support for non-uniform domain decomposition, Doppler broadening, free-gas elastic scattering, general-purpose geometry, hybrid MC/deterministic transport, and depletion. Transport results demonstrate a 2–5× GPU-to-CPU speedup on a per-node basis for an eigenvalue problem on the Frontier supercomputer and a 28× speedup for a fixed-source problem on the Summit supercomputer.

Biondo, Elliott [ORNL] (ORCID:0000000290881360)↗
DOE OSTI2021

Evaluation of Polyethylene and Blank Pulsed Sphere Experiments Using Deuteron Transport Feature in COG

From the late 1960s to about 1985, LLNL undertook a series of ‘pulsed sphere’ experiments for 31 materials involving 148 different experiments using 75 different spheres. The purpose was to measure the neutron leakage spectra for different materials from 14-MeV neutrons generated by 3 H(d,n) 4 He reactions induced by an incident D+ beam from the Insulated Core Transfer (ICT) accelerator at LLNL. This deuteron beam was focused to impinge upon tritium loaded onto a titanium substrate within a low mass target assembly placed within the center of a spherical shell constructed from the material of interest. COG is a three-dimensional, continuous energy, Monte Carlo code, developed by Lawrence Livermore National Laboratory (LLNL). Recent additions to COG include alpha and deuteron particle transport in order to directly calculate (α,n) and 3 H(d,n) 4 He reactions, and resulting neutron emission spectra. Results of the direct deuteron transport calculation for the 0.8 MFP (mean free path) polyethylene sphere with PILOT B detector was previously published. The purpose of this paper is to present expanded COG results showing the neutron peak from deuterium buildup ( 2 H(d,n) 3 He) due to tritium decay in the target. Cross section data libraries used are ENDF/B-VIII.0 for neutrons and ENDF/B-VII.1 for deuterons, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Comparison of Temperature-Dependent Cross Section Treatments Within the Shift Monte Carlo Radiation Transport Code

Multiple calculation methods have been developed to account for Doppler broadening of cross sections in Monte Carlo radiation transport simulations. This work compares the accuracy and performance of four different broadening methods, including an on-the-fly multipole treatment, applied to three different geometries of increasing complexity, including a nuclear reactor assembly with a realistic temperature profile. The comparisons reveal how the combination of physics implementations in the Shift Monte Carlo code result in slightly different temperature-dependent behavior, and different levels of performance, for a variety of materials and situations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 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 ↗
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.↗