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 163 records · Page 9

Study of space charge compensation phenomena in charged particle beams

The propagation of a charged particle beam is accompanied by the production of secondary particles created in the interaction of the beam itself with the background gas flowing in the accelerator tube. In the drift region, where the electric field of the electrodes is negligible, secondary particles may accumulate giving a plasma which shields the self-induced potential of the charged beam. This phenomenon, known as space charge compensation is a typical issue in accelerator physics, where it is usually addressed by means of 1D radial transport codes or Monte Carlo codes. The present paper describes some theoretical studies on this phenomenon, presenting a Particle in Cell-Monte Carlo (PIC-MC) Code developed ad hoc where both radial and axial confinements of secondary particles are calculated. The features of the model, offering a new insight on the problem, are described and some results discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

An Overview of the Monte Carlo Methods, Codes, & Applications Group

This report sketches the work of the Group to deliver first-principle Monte Carlo methods, production quality codes, and radiation transport-based computational and experimental assessments using the codes MCNP and MCATK for such applications as criticality safety, non-proliferation, nuclear energy, nuclear threat reduction and response, radiation detection and measurement, radiation health protection, and stockpile stewardship.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Verification of Monte Carlo transport codes against measured small angle p-, d-, and t-emission in carbon fragmentation at 600 MeV/nucleon

Momentum spectra of hydrogen isotopes have been measured at 3.5° from 12 C fragmentation on a Be target. Momentum spectra cover both the region of fragmentation maximum and the cumulative region. Differential cross sections span five orders of magnitude. The data are compared to predictions of four Monte Carlo codes: QMD, LAQGSM, BC, and INCL++. There are large differences between the data and predictions of some models in the high momentum region. The INCL++ code gives the best and almost perfect description of the data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Generation of the Cross Section Library for PROTEUS

The cross sections for a high-fidelity neutron transport code PROTEUS can be generated in multiple ways. Multigroup cross sections in the ISOTXS format can be generated using the external codes such as MC 2 -3 and Monte Carlo codes such as Serpent and OpenMC. For the Monte Carlo codes, their multigroup cross section outputs should be converted to the ISOTXS format using a post processor named GenISOTXS developed by ANL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Results for Phase I of the IAEA Coordinated Research Program on HTGR Uncertainties

The quantification of uncertainties in design and safety analysis of reactors is today not only broadly accepted, but in many cases became the preferred way to replace traditional conservative analysis for safety and licensing analysis. The use of a more fundamental methodology is also consistent with the reliable high fidelity physics models and robust, efficient, and accurate codes available today. To facilitate uncertainty analysis applications a comprehensive approach and methodology must be developed and applied. High Temperature Gas-cooled Reactors (HTGR) has its own peculiarities, coated particle design, large graphite quantities, different materials and high temperatures that also require other simulation requirements. The IAEA has therefore launched a Coordinated Research Project (CRP) on the HTGR Uncertainty Analysis in Modeling (UAM) in 2013 to study uncertainty propagation specifically in the HTGR analysis chain. Two benchmark problems are defined, with the prismatic design represented by the General Atomics (GA) MHTGR-350 and a 250 MW modular pebble bed design similar to the HTR-PM (INET, China). This report summarizes the contributions of the HTGR Methods Simulation group at Idaho National Laboratory (INL) up to this point of the CRP. The activities at INL have been focused so far on creating the problem specifications for the prismatic design, as well as providing reference solutions for the exercises defined for Phase I. An overview is provided of the HTGR UAM objectives and scope, and the detailed specifications for Exercises I-1, I-2, I-3 and I-4 are also included here for completeness. The main focus of the report is the compilation and discussion of reference results for Phase I (i.e. for input parameters at their nominal or best-estimate values), which is defined as the first step of the uncertainty quantification process. These reference results can be used by other CRP participants for comparison with other codes or their own reference results. The status on the Monte Carlo modeling of the experimental VHTRC facility is also discussed. Reference results were obtained for the neutronics stand-alone cases (Ex. I-1 and Ex. I-2) using the (relatively new) Monte Carlo code Serpent, and comparisons were performed with the more established Monte Carlo codes MCNP and KENO-VI. For the thermal-fluids stand-alone cases (Ex. I-3 and I-4) the commercial CFD code CFX was utilized to obtain reference results that can be compared with lower fidelity tools.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Study of applied magnetic field magnetoplasmadynamic thrusters with particle-in-cell code with Monte Carlo collision. I. Computation methods and physical processes

A two-dimensional axisymmetric electromagnetic particle-in-cell code with Monte Carlo collision conditions has been developed for an applied-field magnetoplasmadynamic thruster simulation. This theoretical approach establishes a particle acceleration model to investigate the microscopic and macroscopic characteristics of particles. This new simulation code was used to study the physical processes associated with applied magnetic fields. In this paper (I), detail of the computation procedure and results of predictions of local plasma and field properties are presented. The numerical model was applied to the configuration of a NASA Lewis Research Center 100-kW magnetoplasmadynamic thruster which has well documented experimental results. The applied magnetic field strength was varied from 0 to 0.12 T, and the effects on thrust were calculated as a basis for verification of the theoretical approach. With this confirmation, the changes in the distributions of ion density, velocity, and temperature throughout the acceleration region related to the applied magnetic fields were investigated. Using these results, the effects of applied field on physical processes in the thruster discharge region could be represented in detail, and those results are reported.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Implementation and Validation of Reference Sensitivity Profile Calculations in TRIPOLI 4 - 14247

In current neutron transport calculations, particularly on continuous-energy Monte-Carlo codes, the major source of uncertainty is coming from nuclear data. In order to improve nuclear data evaluations, integral experiments are used. Feedback is then precisely obtained by re-estimating cross-section levels and correlations according to their sensitivity. k{sub eff}-sensitivity coefficients are traditionally obtained using deterministic codes since hands-on methods were not yet implemented in continuous energy Monte Carlo codes. In this paper, an Iterated-Fission-Probability-based method newly developed in the continuous energy Monte Carlo code TRIPOLI-4 will be used to calculate sensitivity coefficient breakdown on isotopes and reactions and results will be compared to those of ERANOS deterministic code. Objective is twofold: to show the reliability of the new TRIPOLI-4 method and to identify where and why profiles are slightly different. Sensitivity profiles, forward and adjoint fluxes are thus meticulously compared using a fine 1968-group energy structure. An OECD international benchmark (Uncertainty Analysis for Criticality Safety Assessment - UACSA - Phase III) is taken as working example to enforce validation and to compare to independent results (e.g. MCNP-6). This benchmark is constituted of two bare spheres filled with homogeneous mixture of uranium and polyethylene (respectively 2% and 50% U-235 enrichments) and gives an overview of issues encountered in thermal reactor cores. Several energy-group structures (1968G, 172G, 33G) are also tested in ERANOS to identify multi-group and implicit effects. Lessons learnt are that continuous Monte Carlo methods seem to be mature and compare well with deterministic codes. Only damper the U-238 elastic - and on certain conditions the U-238 capture - sensitivity coefficient for which difference can reach 30% between TRIPOLI-4 and ERANOS. Indeed concerning U-238, implicit effects are emphasized and multi-group theory remains limited by hypothesis; even a 1968-group structure is not refined enough to properly calculate U-238 elastic scattering k{sub eff}-sensitivity. This reference 3D continuous-energy sensitivity capability in TRIPOLI4-IFP will be extensively used in the uncertainty automated tool RIB of the French Criticality-Safety package CRISTAL. (authors)

42 ENGINEERING↗

TART2012 An Overview of A Coupled Neutron Photon 3-D, Combinatorial Geometry Time Dependent Monte Carlo Transport Code

This document is designed to provide an overview of TART2012, a coupled neutron photon, 3 Dimensional, combinatorial geometry, time dependent Monte Carlo radiation transport code. This code can run on any modern computer. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART2012 is also incredibly FAST; if you have used similar codes, you will be amazed at how fast this code is compared to other similar codes. Use of the entire system can save you a great deal of time and energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Monte Carlo simulations of {sup 3}He ion physical characteristics in a water phantom and evaluation of radiobiological effectiveness

Purpose: {sup 3}He ions may hold great potential for clinical therapy because of both their physical and biological properties. In this study, the authors investigated the physical properties, i.e., the depth-dose curves from primary and secondary particles, and the energy distributions of helium ({sup 3}He) ions. A relative biological effectiveness (RBE) model was applied to assess the biological effectiveness on survival of multiple cell lines. Methods: In light of the lack of experimental measurements and cross sections, the authors used Monte Carlo methods to study the energy deposition of {sup 3}He ions. The transport of {sup 3}He ions in water was simulated by using three Monte Carlo codes—FLUKA, GEANT4, and MCNPX—for incident beams with Gaussian energy distributions with average energies of 527 and 699 MeV and a full width at half maximum of 3.3 MeV in both cases. The RBE of each was evaluated by using the repair-misrepair-fixation model. In all of the simulations with each of the three Monte Carlo codes, the same geometry and primary beam parameters were used. Results: Energy deposition as a function of depth and energy spectra with high resolution was calculated on the central axis of the beam. Secondary proton dose from the primary {sup 3}He beams was predicted quite differently by the three Monte Carlo systems. The predictions differed by as much as a factor of 2. Microdosimetric parameters such as dose mean lineal energy (y{sub D}), frequency mean lineal energy (y{sub F}), and frequency mean specific energy (z{sub F}) were used to characterize the radiation beam quality at four depths of the Bragg curve. Calculated RBE values were close to 1 at the entrance, reached on average 1.8 and 1.6 for prostate and head and neck cancer cell lines at the Bragg peak for both energies, but showed some variations between the different Monte Carlo codes. Conclusions: Although the Monte Carlo codes provided different results in energy deposition and especially in secondary particle production (most of the differences between the three codes were observed close to the Bragg peak, where the energy spectrum broadens), the results in terms of RBE were generally similar.

61 RADIATION PROTECTION AND DOSIMETRY↗

SUPREM-DSMC: A New Scalable, Parallel, Reacting, Multidimensional Direct Simulation Monte Carlo Flow Code

An AFRL/NRL team has recently been selected to develop a scalable, parallel, reacting, multidimensional (SUPREM) Direct Simulation Monte Carlo (DSMC) code for the DoD user community under the High Performance Computing Modernization Office (HPCMO) Common High Performance Computing Software Support Initiative (CHSSI). This paper will introduce the JANNAF Exhaust Plume community to this three-year development effort and present the overall goals, schedule, and current status of this new code.

Campbell, David↗

Methodology comparison for gamma-heating calculations in material-testing reactors

The Jules Horowitz Reactor (JHR) is a Material-Testing Reactor (MTR) under construction in the south of France at CEA Cadarache (French Alternative Energies and Atomic Energy Commission). It will typically host about 20 simultaneous irradiation experiments in the core and in the beryllium reflector. These experiments will help us better understand the complex phenomena occurring during the accelerated ageing of materials and the irradiation of nuclear fuels. Gamma heating, i.e. photon energy deposition, is mainly responsible for temperature rise in non-fuelled zones of nuclear reactors, including JHR internal structures and irradiation devices. As temperature is a key parameter for physical models describing the behavior of material, accurate control of temperature, and hence gamma heating, is required in irradiation devices and samples in order to perform an advanced suitable analysis of future experimental results. From a broader point of view, JHR global attractiveness as a MTR depends on its ability to monitor experimental parameters with high accuracy, including gamma heating. Strict control of temperature levels is also necessary in terms of safety. As JHR structures are warmed up by gamma heating, they must be appropriately cooled down to prevent creep deformation or melting. Cooling-power sizing is based on calculated levels of gamma heating in the JHR. Due to these safety concerns, accurate calculation of gamma heating with well-controlled bias and associated uncertainty as low as possible is all the more important. There are two main kinds of calculation bias: bias coming from nuclear data on the one hand and bias coming from physical approximations assumed by computer codes and by general calculation route on the other hand. The former must be determined by comparison between calculation and experimental data; the latter by calculation comparisons between codes and between methodologies. In this presentation, we focus on this latter kind of bias. Nuclear heating is represented by the physical quantity called absorbed dose (energy deposition induced by particle-matter interactions, divided by mass). Its calculation with Monte Carlo codes is possible but computationally expensive as it requires transport simulation of charged particles, along with neutrons and photons. For that reason, the calculation of another physical quantity, called KERMA, is often preferred, as KERMA calculation with Monte Carlo codes only requires transport of neutral particles. However, KERMA is only an estimator of the absorbed dose and many conditions must be fulfilled for KERMA to be equal to absorbed dose, including so-called condition of electronic equilibrium. Also, Monte Carlo computations of absorbed dose still present some physical approximations, even though there is only a limited number of them. Some of these approximations are linked to the way how Monte Carlo codes apprehend the transport simulation of charged particles and the productive and destructive interactions between photons, electrons and positrons. There exists a huge variety of electromagnetic shower models which tackle this topic. Differences in the implementation of these models can lead to discrepancies in calculated values of absorbed dose between different Monte Carlo codes. The magnitude of order of such potential discrepancies should be quantified for JHR gamma-heating calculations. We consequently present a two-pronged plan. In a first phase, we intend to perform compared absorbed dose / KERMA Monte Carlo calculations in the JHR. This way, we will study the presence or absence of electronic equilibrium in the different JHR structures and experimental devices and we will give recommendations for the choice of KERMA or absorbed dose when calculating gamma heating in the JHR. In a second phase, we intend to perform compared TRIPOLI4 / MCNP absorbed dose calculations in a simplified JHR-representative geometry. For this comparison, we will use the same nuclear data library for both codes (the European library JEFF3.1.1 and photon library EPDL97) so as to isolate the effects from electromagnetic shower models on absorbed dose calculation. This way, we hope to get insightful feedback on these models and their implementation in Monte Carlo codes. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Spent fuel utilization in a compact traveling wave reactor

In recent years, several innovative designs of nuclear reactors are proposed. One of them is Traveling Wave Reactor (TWR). The unique characteristic of a TWR is the capability of breeding its own fuel in the reactor. The reactor is fueled by mostly depleted, natural uranium or spent nuclear fuel and a small amount of enriched uranium to initiate the fission process. Later on in the core, the reactor gradually converts the non-fissile material into the fissile in a process like a traveling wave. In this work, a TWR with spent nuclear fuel blanket was studied. Several parameters such as reactivity coefficients, delayed neutron fraction, prompt neutron generation lifetime, and fission power, were analyzed. The discharge burnup composition was also analyzed. The calculation is performed by a continuous energy Monte Carlo code McCARD.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Nuclear fragmentation measurements for hadrontherapy and space radiation protection

Nuclear fragmentation measurements are necessary in hadrontherapy and space radiation protection, to predict the effects of the ion nuclear interactions within the human body. Nowadays, a very limited set of carbon fragmentation cross sections has been measured and in particular, to our knowledge, no double differential fragmentation cross sections at intermediate energies are available in literature. We have measured the double differential cross sections and the angular distributions of the secondary fragments produced in the {sup 12}C fragmentation at 62 AMeV on a thin carbon target. The experimental data have been also used to benchmark the prediction capability of the Geant4 Monte Carlo code at intermediate energies, where it was never tested before.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Spectra of photons and neutrons generated in a heterogeneous head of a 15 MV LINAC at differents field sizes

Spectra of photons and neutrons were calculated, using the Monte Carlo code MCNP-5 using the e/p/n mode. A heterogeneous model was used to define the linac head where the collimators were modeled to produce five different treatment fields at the isocenter. Photon and neutron spectra were estimated in several points along two directions from the isocenter. The total photon fluence beyond 60 cm behaves according to 1/r{sup 2} rule, while total neutron fluence, beyond 80 cm, can be described by diffusion theory using an infinite plane as a neutron source.

43 PARTICLE ACCELERATORS↗

Performance of a Low Activity Beta-Sensitive SR{sup 90} Water Monitor for Fukushima

There are large volumes of contaminated water from the stabilization efforts at the damaged Fukushima Nuclear Power Plants. This water is being processed to remove radioactivity for eventual release to the environment. An on-line continuously operating system to confirm that the clean-up system is working properly, and to provide prompt feedback of the results is required by the system operator. While gamma emitting nuclides allow for the straight forward approach of gamma spectroscopy to identify and quantify radioactivity in water, pure beta emitting nuclides such as Sr{sup 90} pose a challenging problem. The relatively short range of beta radiation in water requires optimization of the measurement geometry in terms of the source-detector distance and source-detector interface while retaining a background sensitivity low enough to meet the Minimum Detectable Concentration (MDC) of 10 Bq/kg in 180 minutes. This issue is complicated by the continuum nature of the beta spectrum which does not allow for simple nuclide identification. The use of the Monte-Carlo code MCNP to estimate system performance before prototyping vastly increases the success of the end product. Various parameters such as detector size and thickness, water chamber size, water chamber construction materials were evaluated to help choose the optimum geometry. The final design was a system consisting of two large-area (16 x 35 cm) and thin (0.15 mm) plastic scintillators placed very close to a sealed a water chamber. The size of the chamber was optimized to obtain the maximum efficiency for the nuclide being measured (Sr/Y{sup 90}) but to minimize the efficiency for possible interferences (Ru/Rh{sup 106}, Cs{sup 137}). A thin carbon fiber window was selected with adequate material and thickness to contain the water under pressure, but also thin enough (0.5 mm) to allow enough beta radiation to pass through to the active detector volume. The entire measurement geometry is then housed in a thick lead shield to reduce contributions from external sources to an acceptable level. Data acquisition is accomplished through customized application-specific software that allows for long counting times to attain a low MDC, but also simultaneously provides alarms on short averaging times to achieve a fast response to sudden changes in activity concentration. Multiple monitors are then linked to supervisory software where real time data and alarms are available for analysis in remote locations. The system also allows for remote operation of the unit; check sources, background checks, systems settings and more may be accessed remotely. Testing of the production devices has shown that we can achieve the 10 Bk/kg MDC requirement for Sr{sup 90} in equilibrium with Y{sup 90} with a count time of approximately 20 minutes. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Application of electron paramagnetic resonance imaging to the characterization of the Ultem(R) exposed to 1 MeV electrons. Correlation of radical density data to tiger code calculations

A major long-term goal of the Materials Division at the NASA Langley Research Center is the characterization of new high-performance materials that have potential applications in the aircraft industry, and in space. The materials used for space applications are often subjected to a harsh and potentially damaging radiation environment. The present study constitutes the application of a novel technique to obtain reliable data for ascertaining the molecular basis for the resilience and durability of materials that have been exposed to simulated space radiations. The radiations of greatest concern are energetic electrons and protons, as well as galactic cosmic rays. Presently, the effects of such radiation on matter are not understood in their entirety. It is clear however, that electron radiation causes ionization and homolytic bond rupture, resulting in the formation of paramagnetic spin centers in the polymer matrices of the structural materials. Since the detection and structure elucidation of paramagnetic species are most readily accomplished using Electron Paramagnetic Resonance (EPR) Spectroscopy, the NASA LaRC EPR system was brought back on-line during the 1991 ASEE term. The subsequent 1992 ASEE term was devoted to the adaptation of the EPR core system to meet the requirements for EPR Imaging (EPRI), which provides detailed information on the spatial distribution of paramagnetic species in bulk media. The present (1994) ASEE term was devoted to the calibration of this EPR Imaging system, as well as to the application of this technology to study the effects of electron irradiation on Ultem(exp R), a high performance polymer which is a candidate for applications in aerospace. The Ultem was exposed to a dose of 2.4 x 10(exp 9) Rads (1-MeV energy/electron) at the LaRC electron accelerator facility. Subsequently, the exposed specimens were stored in liquid nitrogen, until immediately prior to analyses by EPRI. The intensity and dimensions of the EPR Images that were generated for the irradiated specimens showed that the electrons penetrated the material to a depth of approximately 0.125 inch. These data show a very high degree of correlation to the energy deposition profile as predicted by the Tiger Code, a Monte Carlo code that provides guidelines for the transport of electrons in matter. Subsequent efforts will focus on delineating the transport properties of energetic protons in Ultem(R).

Suleman, Naushadalli K.↗

Impact Study on the Methodology used for Photon- Heating Determination in Material-Testing Reactors

Determination of photon heating by calculation is an important issue for the Jules Horowitz Reactor (JHR), the next international Material-Testing Reactor (MTR) under construction at CEA Cadarache research center in the south of France. Accurate knowledge of photon heating in structure materials and irradiation devices is necessary for JHR design and safety studies. In this paper, we quantify the impact of different photon-heating calculation routes by comparing absorbed dose and KERMA calculations (Kinetic Energy Released per Mass) from two different Monte Carlo codes, TRIPOLI-4.9{sup R} and MCNP (Monte Carlo N-Particle transport code). These calculations are carried out in JHR-representative geometries with the nuclear-data library JEFF3.1.1 and the photon-data library EPDL97. Discrepancies amounting to up to 18% between absorbed dose and KERMA are found in JHR irradiation devices and are linked to charged-particle transport effects taking place in heterogeneous materials of small dimensions. In a JHR assembly cell, discrepancies of about 1% on photon KERMA and of about 3% on absorbed dose are highlighted between the two Monte Carlo codes. These latter discrepancies are small compared to typical sources of uncertainty for Monte Carlo calculation (for instance, nuclear data uncertainty) and are supposed to be due to differences in the processing of gamma-production data by neutron interactions and to differences in electromagnetic-shower models and implementation between the two codes. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗