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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Benchmark Analysis of Pion Contribution from Galactic Cosmic Rays

Shielding strategies for extended stays in space must include a comprehensive resolution of the secondary radiation environment inside the spacecraft induced by the primary, external radiation. The distribution of absorbed dose and dose equivalent is a function of the type, energy and population of these secondary products. A systematic verification and validation effort is underway for HZETRN, which is a space radiation transport code currently used by NASA. It performs neutron, proton and heavy ion transport explicitly, but it does not take into account the production and transport of mesons, photons and leptons. The question naturally arises as to what is the contribution of these particles to space radiation. The pion has a production kinetic energy threshold of about 280 MeV. The Galactic cosmic ray (GCR) spectra, coincidentally, reaches flux maxima in the hundreds of MeV range, corresponding to the pion production threshold. We present results from the Monte Carlo code MCNPX, showing the effect of lepton and meson physics when produced and transported explicitly in a GCR environment.

Aghara, Sukesh K.↗
DOE OSTI2014

First principles transport coefficients and reaction rates of Ar{sub 2}{sup +} ions in argon for cold plasma jet modeling

Momentum-transfer collision cross-sections and integral collision cross-sections for the collision-induced dissociation are calculated for collisions of ionized argon dimers with argon atoms using a nonadiabatic semiclassical method with the electronic Hamiltonian calculated on the fly via a diatomics-in-molecules semiempirical model as well as inverse-method modeling based on simple isotropic rigid-core potential. The collision cross-sections are then used in an optimized Monte Carlo code for evaluations of the Ar{sub 2}{sup +} mobility in argon gas, longitudinal diffusion coefficient, and collision-induced dissociation rates. A thorough comparison of various theoretical calculations as well as with available experimental data on the Ar{sub 2}{sup +} mobility and collision cross-sections is performed. Good agreement is found between both theoretical approaches and the experiment. Analysis of the role of inelastic processes in Ar{sub 2}{sup +}/Ar collisions is also provided.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2012

Benchmark of Atucha-2 PHWR RELAP5-3D control rod model by Monte Carlo MCNP5 core calculation

Atucha-2 is a Siemens-designed PHWR reactor under construction in the Republic of Argentina. Its geometrical complexity and peculiarities require the adoption of advanced Monte Carlo codes for performing realistic neutronic simulations. Therefore core models of Atucha-2 PHWR were developed using MCNP5. In this work a methodology was set up to collect the flux in the hexagonal mesh by which the Atucha-2 core is represented. The scope of this activity is to evaluate the effect of obliquely inserted control rod on neutron flux in order to validate the RELAP5-3D{sup C}/NESTLE three dimensional neutron kinetic coupled thermal-hydraulic model, applied by GRNSPG/UNIPI for performing selected transients of Chapter 15 FSAR of Atucha-2. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

Temperature effects of resonance scattering for epithermal neutrons in MCNP

Epithermal neutron elastic scattering can be significantly affected by the thermal motion of target nuclides. Since the 1950's continuous-energy Monte Carlo codes have generally accounted for the target motion using a free gas scattering model, with the assumption that the scattering cross-section is constant in energy. Recent work has shown the importance of resonance scattering, and several methods for an improved free-gas treatment have been developed. We have implemented a rejection-based sampling scheme in the MCNP free-gas treatment to account for cross-section variation. The modified MCNP code was used to investigate a number of practical concerns: results for an LWR Doppler defect benchmark; computational costs; and energy limits for the free-gas treatment. Additionally, the impact on a suite of ICSBEP criticality benchmark problems (at room temperature) was determined to be negligible, an important result since such problems are used extensively in testing and evaluating revisions to ENDF/B-VII nuclear data. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

High energy x-ray phase contrast CT using glancing-angle grating interferometers

Purpose: The authors present initial progress toward a clinically compatible x-ray phase contrast CT system, using glancing-angle x-ray grating interferometry to provide high contrast soft tissue images at estimated by computer simulation dose levels comparable to conventional absorption based CT. Methods: DPC-CT scans of a joint phantom and of soft tissues were performed in order to answer several important questions from a clinical setup point of view. A comparison between high and low fringe visibility systems is presented. The standard phase stepping method was compared with sliding window interlaced scanning. Using estimated dose values obtained with a Monte-Carlo code the authors studied the dependence of the phase image contrast on exposure time and dose. Results: Using a glancing angle interferometer at high x-ray energy (∼45 keV mean value) in combination with a conventional x-ray tube the authors achieved fringe visibility values of nearly 50%, never reported before. High fringe visibility is shown to be an indispensable parameter for a potential clinical scanner. Sliding window interlaced scanning proved to have higher SNRs and CNRs in a region of interest and to also be a crucial part of a low dose CT system. DPC-CT images of a soft tissue phantom at exposures in the range typical for absorption based CT of musculoskeletal extremities were obtained. Assuming a human knee as the CT target, good soft tissue phase contrast could be obtained at an estimated absorbed dose level around 8 mGy, similar to conventional CT. Conclusions: DPC-CT with glancing-angle interferometers provides improved soft tissue contrast over absorption CT even at clinically compatible dose levels (estimated by a Monte-Carlo computer simulation). Further steps in image processing, data reconstruction, and spectral matching could make the technique fully clinically compatible. Nevertheless, due to its increased scan time and complexity the technique should be thought of not as replacing, but as complimentary to conventional CT, to be used in specific applications.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1989

Computer simulation of dust grain evolution

The latest results are reported from a Monte Carlo code that is being developed at NASA Ames. The goal of this program, is to derive from the observed and presumed properties of the interstellar medium (ISM) the following information: (1) the size spectrum of interstellar dust; (2) the chemical structure of interstellar dust; (3) interstellar abundances; and (4) the lifetime of a dust grain in the ISM. Presently this study is restricted to refractory interstellar material, i.e., the formation and destruction of ices are not included in the program. The program is embedded in an analytic solution for the bulk chemical evolution of a two-phase interstellar medium in which stars are born in molecular clouds, but new nucleosynthesis products and stellar return are entered into a complementary intercloud medium. The well-mixed matter of each interstellar phase is repeatedly cycled stochastically through the complementary phase and back. Refractory dust is created by thermal condensation as stellar matter flows away from sites of nucleosynthesis such as novae and supernovae and/or from the matter returned from evolved intermediate stars. The history of each particle is traced by standard Monte Carlo techniques as it is sputtered and fragmented by supernova shock waves in the intercloud medium. It also accretes an amorphous mantle of gaseous refractory atoms when its local medium joins with the molecular cloud medium. Finally it encounters the possibility of astration (destruction by star formation) within the molecular clouds.

Liffman, K.↗
DOE OSTI2015

Fast Biological Modeling for Voxel-based Heavy Ion Treatment Planning Using the Mechanistic Repair-Misrepair-Fixation Model and Nuclear Fragment Spectra

Purpose: The physical and biological differences between heavy ions and photons have not been fully exploited and could improve treatment outcomes. In carbon ion therapy, treatment planning must account for physical properties, such as the absorbed dose and nuclear fragmentation, and for differences in the relative biological effectiveness (RBE) of ions compared with photons. We combined the mechanistic repair-misrepair-fixation (RMF) model with Monte Carlo-generated fragmentation spectra for biological optimization of carbon ion treatment plans. Methods and Materials: Relative changes in double-strand break yields and radiosensitivity parameters with particle type and energy were determined using the independently benchmarked Monte Carlo damage simulation and the RMF model to estimate the RBE values for primary carbon ions and secondary fragments. Depth-dependent energy spectra were generated with the Monte Carlo code FLUKA for clinically relevant initial carbon ion energies. The predicted trends in RBE were compared with the published experimental data. Biological optimization for carbon ions was implemented in a 3-dimensional research treatment planning tool. Results: We compared the RBE and RBE-weighted dose (RWD) distributions of different carbon ion treatment scenarios with and without nuclear fragments. The inclusion of fragments in the simulations led to smaller RBE predictions. A validation of RMF against measured cell survival data reported in published studies showed reasonable agreement. We calculated and optimized the RWD distributions on patient data and compared the RMF predictions with those from other biological models. The RBE values in an astrocytoma tumor ranged from 2.2 to 4.9 (mean 2.8) for a RWD of 3 Gy(RBE) assuming (α/β){sub X} = 2 Gy. Conclusions: These studies provide new information to quantify and assess uncertainties in the clinically relevant RBE values for carbon ion therapy based on biophysical mechanisms. We present results from the first biological optimization of carbon ion radiation therapy beams on patient data using a combined RMF and Monte Carlo damage simulation modeling approach. The presented method is advantageous for fast biological optimization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

Temperature Effects of Resonance Scattering for Epithermal Neutrons in MCNP

Epithermal neutron elastic scattering can be significantly affected by the thermal motion of target nuclides. Since the 1950s continuous-energy Monte Carlo codes have generally accounted for the target motion using a free gas scattering model, with the assumption that the scattering cross section is constant in energy. Recent work has shown the importance of resonance scattering, and several methods for an improved free-gas treatment have been developed. We have implemented a rejection-based sampling scheme in the MCNP free-gas treatment to account for cross-section variation. The modified MCNP code was used to investigate a number of practical concerns: results for an LWR Doppler defect benchmark; computational costs; and energy limits for the free gas treatment. Additionally, the impact on a suite of ICSBEP criticality benchmark problems (at room temperature) was determined to be negligible, an important result since such problems are used extensively in testing and evaluating revisions to ENDFIB-VII nuclear data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Derivation and verification of the direct-sampling method for simulating Monte Carlo flight paths in tetrahedral meshes with linear finite-element cross sections

This paper provides a derivation of a direct-sampling approach for modeling continuously varying cross sections in tetrahedral-mesh-based Monte Carlo codes. Specifically, cross sections are spatially approximated using linear nodal finite elements. A linearization strategy is provided for non-linearly varying cross sections. The method is verified against seven analytical pure-absorber test problems. These test problems also highlight the benefit of using linear finite elements over element-wise-constant cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Intergalactic magnetic fields and gamma-ray observations of extreme TeV blazars

The intergalactic magnetic field (IGMF) in cosmic voids can be indirectly probed through its effect on electromagnetic cascades initiated by a source of teraelectronvolt (TeV) gamma-rays, such as active galactic nuclei (AGNs). AGNs that are sufficiently luminous at TeV energies, 'extreme TeV blazars', can produce detectable levels of secondary radiation from inverse Compton scattering of the electrons in the cascade, provided that the IGMF is not too large. We review recent work in the literature that utilizes this idea to derive constraints on the IGMF for three TeV-detected blazars, 1ES 0229+200, 1ES 1218+304, and RGB J0710+591, and we also investigate four other hard-spectrum TeV blazars in the same framework. Through a recently developed, detailed, three-dimensional particle-tracking Monte Carlo code, incorporating all major effects of QED and cosmological expansion, we research the effects of major uncertainties, such as the spectral properties of the source, uncertainty in the ultraviolet and far-infrared extragalactic background light, undersampled very high energy (energy ≥100 GeV) coverage, past history of gamma-ray emission, source versus observer geometry, and the jet AGN Doppler factor. The implications of these effects on the recently reported lower limits of the IGMF are thoroughly examined to conclude that the presently available data are compatible with a zero-IGMF hypothesis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SU-E-T-560: Monte Carlo Simulation of the Neutron Radiation Field Around a Medical 18 MV Linac

Purpose: Today the majority of radiation therapy treatments are performed at medical electron linear accelerators (linacs). The accelerated electrons are used for the generation of bremsstrahlung photons. The use of higher electron respectively photon energies has some advantages over lower energies such as the longer dose build-up. However photons with energies higher than ∼7 MeV can additionally to the interaction with bound electrons undergo inelastic reactions with nuclei. These photonuclear reactions lead to the emission of fast neutrons which contaminate the primary photon field. The neutrons might penetrate through the collimators and deliver out-of-field dose to the patient. Furthermore the materials inside the linac head as well as the air inside the treatment room get activated which might deliver dose to the medical employees even when the linac is not in operation. A detailed knowledge of these effects is essential for adequate radiation protection of the employees and an optimal patient treatment. Methods: It is a common method to study the radiation fields of such linacs by means of Monte Carlo simulations. For the investigation of the effects caused by photonuclear reactions a typical linac in high energy mode (Varian Clinac 18 MV-X) as well as the surrounding bunker were modelled and simulated using the Monte Carlo code FLUKA which includes extensive nuclear reaction and neutron transport models additional to electron-photon transport as well as capabilities for a detailed study of effective dose distributions and activation yields. Results: Neutron spectra as well as neutron effective dose distributions within the bunker were obtained, reaching up to some mSv/Gy in the patient’s plane. The results are normalized per Gy in the depth dose maximum at 10×10 cm{sup 2} field size. Therefore an absolute interpretation is possible. Conclusion: The obtained data gives a better understanding of the photonuclear reaction caused effects.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

Increasing efficiency of BEAMnrc-simulated Co-60 beams using directional source biasing

Purpose: This study describes the implementation of a directional source biasing (DSB) scheme for efficiently simulating Cobalt-60 treatment heads using the BEAMnrc Monte Carlo code. Previous simulation of Co-60 beams with BEAMnrc was impractical because of the time required to track photons not directed into the treatment field and to simulate secondary charged particles. Methods: In DSB, efficiency is increased by splitting each photon emitted by the Co-60 source a user-defined number of times. Only those split primary photons directed into a user-defined splitting field (encompassing the treatment field) are sampled, yielding many low-weight photons directed into the field. Efficiency can be further increased by taking advantage of radial symmetry at the top of the treatment head to reduce the number of split primary photons tracked in this portion. There is also an option to generate contaminant electrons in DSB. Results: The DSB scheme in BEAMnrc increases the photon fluence calculation efficiency in a 10 × 10 cm{sup 2} Co-60 beam by a factor of 1800 with a concurrent increase in contaminant electron fluence calculation efficiency by a factor of 1200. Implementation of DSB in beampp, a C++ code for accelerator simulations based on EGSnrc and the C++ class library, egspp, increases photon fluence efficiency by a factor of 2800 and contaminant electron fluence efficiency by a factor of 1600. Optimum splitting numbers are in the range of 20 000–40 000. For dose calculations in a water phantom (0.5 × 0.5 × 0.5 cm{sup 3} voxels) this translates into a factor of ∼400 increase in dose calculation efficiency (all doses > 0.5 × D{sub max}). An example calculation of the ratio of dose to water to dose to chamber (the basis of the beam quality correction factor) to within 0.2% in a realistic chamber using a full simulation of a Co-60 treatment head as a source indicates the practicality of Co-60 simulations with DSB. Conclusions: The efficiency improvement resulting from DSB makes Monte Carlo commissioning of Co-60 beams and calculation of beam quality correction factors feasible.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Implementation of Triply Periodic Minimal Surfaces (TPMS) as surface objects in OpenMC

Triply Periodic Minimal Surfaces (TPMS) represent a promising geometry for future fuel designs due to their significant surface-to-volume ratio, which facilitates efficient cooling of nuclear fuel, a crucial factor for safety and efficiency. Demonstrating the remarkable capabilities of TPMS fuel requires initial modeling and simulation. This paper presents an implementation of TPMS in the Monte Carlo code OpenMC, enabling reactor physics modeling of TPMS. Here, the primary advantages over traditional methods using CAD files include reduced memory requirements for computations and high-fidelity implementation. This implementation has been tested against CAD files loaded in Serpent2, yielding promising results with low biases in the $k_{\textrm{eff}}$, comparable to biases in the material balance sheet. The implementation presented in this work will be used in future reactor physics computations related to new reactor designs involving TPMS-based fuels.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Exascale Multiphysics Nuclear Reactor Simulations for Advanced Designs

ENRICO is a coupled application developed under the U.S. Department of Energy's Exascale Computing Project (ECP) targeting the modeling of advanced nuclear reactors. It couples radiation transport with heat and fluid simulation, including the high-fidelity, highresolution Monte-Carlo code Shift and the Computational fluid dynamics code NekRS. NekRS is a highly-performant open-source code for simulation of incompressible and low-Mach fluid flow, heat transfer, and combustion with a particular focus on turbulent flows in complex domains. It is based on rapidly convergent high-order spectral element discretizations that feature minimal numerical dissipation and dispersion. State-of-the-art multilevel preconditioners, efficient high-order time-splitting methods, and runtime-adaptive communication strategies are built on a fast OCCA-based kernel library, libParanumal, to provide scalability and portability across the spectrum of current and future high-performance computing platforms. On Frontier, Nek5000/RS has recently achieved an unprecedented milestone in breaching over 1 billion spectral elements and 350 billion degrees of freedom. Shift has demonstrated the capability to transport upwards of 1 billion particles per second in full core nuclear reactor simulations featuring complete temperature-dependent, continuous-energy physics on Frontier. Shift achieved a weak-scaling efficiency of 97.8% on 8192 nodes of Frontier and calculated 6 reactions in 214,896 fuel pin regions below 1% statistical error yielding first-of-a-kind resolution for a Monte Carlo transport application.

Hamilton, Steven P.↗
DOE OSTI2013

Fusion alpha-particle losses in a high-beta rippled tokamak

In tokamak plasmas, the confinement of energetic ions depends on the magnetic field structure. If the plasma pressure is finite, the equilibrium current (i.e., the Pfirsch-Schlüter current and diamagnetic current) flows in the plasma to maintain the magnetohydrodynamic (MHD) equilibrium. These plasma currents generate poloidal and toroidal magnetic field and alter the field structure. Moreover, if we consider the non-axisymmetry of magnetic field structures such as toroidal field (TF) ripples, the non-axisymmetric component of the equilibrium current can alter TF ripples themselves. When the plasma beta becomes high, the changes in the field structure due to the equilibrium current might affect the confinement of energetic ions significantly. We intend to clarify how these currents alter the field structure and affect the confinement of alpha particles in high-beta plasma. The MHD equilibrium is calculated using VMEC and the orbits of fusion alpha particles are followed by using the fully three-dimensional magnetic field orbit-following Monte Carlo code. In relatively low-beta plasma (e.g., the volume-averaged beta value <β>≤2%), the changes in the magnetic field component due to the plasma current negligibly affect the confinement of alpha particles except for the Shafranov shift effect. However, for <β>≥3%, the diamagnetic effect reduces the magnetic field strength and significantly increases alpha-particle losses. In these high-beta cases, the non-axisymmetric field component generated by the equilibrium current also increases these losses, but not as effectively as compared to the diamagnetic effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Fission matrix capability for MCNP, Part II - Applications

This paper describes the initial experience and results from implementing a fission matrix capability into the MCNP Monte Carlo code. The fission matrix is obtained at essentially no cost during the normal simulation for criticality calculations. It can be used to provide estimates of the fundamental mode power distribution, the reactor dominance ratio, the eigenvalue spectrum, and higher mode spatial eigenfunctions. It can also be used to accelerate the convergence of the power method iterations. Past difficulties and limitations of the fission matrix approach are overcome with a new sparse representation of the matrix, permitting much larger and more accurate fission matrix representations. Numerous examples are presented. A companion paper (Part I - Theory) describes the theoretical basis for the fission matrix method. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Multigroup cross section generation capability in GRIFFIN

GRIFFIN is an advanced reactor multiphysics application built on the object-oriented simulation environment (MOOSE) and is jointly developed by Idaho National Laboratory and Argonne National Laboratory. The cross section application programming interface, originally developed for the PROTEUS code, has been integrated into GRIFFIN to prepare cross sections for thermal reactor applications with heterogeneous geometries. Additional improvements have been made by implementing an on-the-fly slowing down method, a double heterogeneity treatment capability, and updating the procedure to generate the fine multigroup library. The cross section preparation capability in GRIFFIN was verified for graphite-moderated TRISO fuel-based reactor benchmark problems: unit-cell problems of VHTR and EMPIRE micro reactor and HTTR assembly problems. Eigenvalues and multigroup cross sections of GRIFFIN agreed very well with those of the continuous-energy Monte Carlo code Serpent2 within 200 pcm in eigenvalue and 2% in cross sections. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

SPECTROPOLARIMETRY OF SUPERLUMINOUS SUPERNOVAE: INSIGHT INTO THEIR GEOMETRY

We present the first spectropolarimetric observations of a hydrogen-free superluminous supernova (SLSN) at z = 0.1136, namely SN 2015bn. The transient shows significant polarization at both of the observed epochs: one 24 days before maximum light in the rest-frame, and the other at 27 days after peak luminosity. Analysis of the Q – U plane suggests the presence of a dominant axis and no physical departure from the main axis at either epoch. The polarization spectrum along the dominant axis is characterized by a strong wavelength dependence and an increase in the signal from the first to the second epoch. We use a Monte Carlo code to demonstrate that these properties are consistent with a simple toy model that adopts an axisymmetric ellipsoidal configuration for the ejecta. We find that the wavelength dependence of the polarization is possibly due to a strong wavelength dependence in the line opacity, while the higher level of polarization at the second epoch is a consequence of the increase in the asphericity of the inner layers of the ejecta or the fact that the photosphere recedes into less spherical layers. The geometry of the SLSN is similar to that of stripped-envelope core-collapse SNe connected to GRB, while the overall evolution of the ejecta shape could be consistent with a central engine.

79 ASTRONOMY AND ASTROPHYSICS↗