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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Upgrades to Monteburns, version 3.0

Monteburns, a Monte Carlo burnup code which has the flexibility to model time-dependent isotopic changes for a variety of nuclear systems by linking the neutron transport code MCNP/X to a production/depletion code, has undergone several performance upgrades recently that have increased the code's capabilities. Once limited to a specific number of regions, enhancements have been implemented that afford a much greater number of burn materials, such that users will be more limited by the physical constraints of their computing environment as opposed to inherent limits built into the coding of Monteburns. In conjunction with the increase in the number of burn materials, parallel execution of a production/depletion code of choice has been implemented, such that users have the option of using CINDER90, 0RIGEN2, or the newly released version of ORIGEN. Finally, the recoverable energy per fission calculation was upgraded to include capture gamma energy deposited in all specified materials as a function of irradiation time. The sum of the prompt and delayed recoverable energies from fission was obtained as before. These upgrades were first tested on a rigorous 1/8 core model of a Pressurized Water Reactor with fresh, once- and twice-burned fuel. We can now model several orders of magnitude more materials using Monte Carlo techniques, which is a significant advance in the reactor modeling world. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

SCALE 6.2 Lattice Physics Performance Assessment

The US Nuclear Regulatory Commission relies on the lattice physics analysis capabilities of the SCALE code system to perform confirmatory licensing analyses. Either SCALE lattice physics code—TRITON/NEWT or Polaris—can be used to generate cross section data used by the PARCS nodal core simulator for full-core neutronics calculations. This report presents an assessment of the accuracy of SCALE lattice physics codes for preparation of lattice physics data that are used to support simulator codes such as the NRC’s PARCS, for UO 2 -mixed oxide (MOX)/Zr fueled light water reactor (LWR) analyses. Due to the nature of lattice physics calculations, critical reactor experiment benchmarks cannot be modeled in explicit detail in a lattice physics code. However, this limitation does not mean that these measurement data are not usable for lattice physics studies. Therefore, either geometry approximations or axial buckling must be implemented to determine the critical water height. These modeling limitations have led to development of a three-phase assessment strategy. In the first phase, selected critical experiment benchmarks are modeled using the SCALE 3D continuous-energy (CE) Monte Carlo (MC) code KENO, which is the most rigorous neutron transport method available in SCALE, with no approximations in the spatial, angle, or energy treatments. Biases and statistical uncertainties in quantities of interest such as k eff and pin power distributions are determined by comparing CE KENO results to experimental data. This first phase in the accuracy assessment is to demonstrate and establish the use of CE KENO as a reference solution for the second phase. In the second phase, 14 numerical test suites are used to compare SCALE lattice physics calculations with CE KENO as a reference solution. In the third phase, SCALE lattice physics depletion calculations are performed, and the spent fuel isotopic results are compared with available radiochemical assay measurements. Isotopic measurement comparisons provide quantitative assessment of isotopic density distribution predictions with the depletion models in SCALE lattice physics codes. This report documents results for all test suites. The assessment was performed using standard production techniques unless otherwise noted. Both TRITON/NEWT and Polaris exhibited acceptable accuracy for most test cases. For the few test cases in which acceptable accuracy criteria were not met, further code and data development are planned. The computer codes used in this assessment are as follows: SCALE 6.2 rev19189 (pre-release of 6.2.1) was used for CE KENO, TRITON/NEWT and Polaris (PWR only) calculations; SCALE 6.2.2 was used for Polaris calculations to address several updates in support of support boiling water reactor (BWR) lattice geometries. For this work, there is no difference between SCALE 6.2.1 and SCALE 6.2.2 beyond the inclusion of the Polaris BWR analysis capability; The ENDF/B-VII.1 continuous energy and 252 group libraries deployed in SCALE 6.2 were generated with AMPX 6.2 that is distributed with SCALE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

A Solid of Revolution Time Study using COG11.1 and MCNP6.1

The performance of two continuous energy Monte Carlo transport codes, COG11.1 and MCNP6.1, has been compared for the {sup 239}Pu Jezebel Benchmark using cross sections based on the ENDF/B-VII.1 evaluated nuclear data library. The purpose of this study is to determine the length of time it takes each code to complete a series of criticality calculations, using two different methods to model a solid of revolution. A solid of revolution is essential in many fields of engineering and mathematics, and has been a standard feature in computer-aided design and manufacturing since the 1970's. In the field of criticality safety, a solid of revolution can be used to reproduce and simplify complex engineering drawings and systems that would otherwise require excessive use of geometric primitives and computational resources to achieve a comparable level of fidelity. There are two basic methods to model a solid of revolution: Method 1 (COG) specifies a set of x-y pairs or polar coordinates and rotate the curve about the x-axis. Method 2 (COG and MCNP) creates a set of end-to-end conical frustums and/or cylinders. Method 1 is only available in COG and can be implemented using one surface and sector specification. Method 2 can be implemented in both COG and MCNP. This time study was initially conceived as a result of criticality safety analysis performed for a complex system with a high degree of rotational symmetry. Each COG input deck typically took 1-2 days to run in serial mode. COG input decks were then converted to MCNP format and run on the same computer cluster. Instead of expected 1-2 day runtimes, it took MCNP6.1 more than one week to complete each run. Since MCNP typically runs faster than many codes, these lengthy runtimes were considered unusual. This study attempts to investigate the runtime discrepancy between the two codes. The methods and benchmark in this time study were selected to provide a consistent, generalized approach to modeling an arbitrary solid of revolution. While there are many different ways to model the {sup 239}Pu Jezebel Benchmark, the focus of this study is on investigating the practical difference between the two basic methods to model a surface of revolution. Throughout this study, the only code that maintains reasonably low runtimes is COG11.1 (Method 1). Based on the observed performance of each code and method, it is likely that the divergent runtimes are the result of tracking particles across planes that lie between conical frustums, as opposed to defining the object boundary as a single surface. Slight fluctuations in graphed results were also observed, due to random sampling and CPU loading. This was confirmed by re-running input decks and comparing average runtimes to a preliminary curve fit. From the perspective of a criticality safety analyst, the results of this study are expected to contribute towards understanding one of the practical trade-offs between codes and methods as they apply to real-world analysis. Ultimately, it's about choosing the right tool for the job.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Sci-Sat AM: Radiation Dosimetry and Practical Therapy Solutions - 05: Not all geometries are equivalent for magnetic field Fano cavity tests

The coupling of MRI and radiation treatment systems for the application of magnetic resonance guided radiation therapy necessitates a reliable magnetic field capable Monte Carlo (MC) code. In addition to the influence of the magnetic field on dose distributions, the question of proper calibration has arisen due to the several percent variation of ion chamber and solid state detector responses in magnetic fields when compared to the 0 T case (Reynolds et al., Med Phys, 2013). In the absence of a magnetic field, EGSnrc has been shown to pass the Fano cavity test (a rigorous benchmarking tool of MC codes) at the 0.1 % level (Kawrakow, Med.Phys, 2000), and similar results should be required of magnetic field capable MC algorithms. To properly test such developing MC codes, the Fano cavity theorem has been adapted to function in a magnetic field (Bouchard et al., PMB, 2015). In this work, the Fano cavity test is applied in a slab and ion-chamber-like geometries to test the transport options of an implemented magnetic field algorithm in EGSnrc. Results show that the deviation of the MC dose from the expected Fano cavity theory value is highly sensitive to the choice of geometry, and the ion chamber geometry appears to pass the test more easily than larger slab geometries. As magnetic field MC codes begin to be used for dose simulations and correction factor calculations, care must be taken to apply the most rigorous Fano test geometries to ensure reliability of such algorithms.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Automatic treatment of the variance estimation bias in TRIPOLI-4 criticality calculations

The central limit (CLT) theorem States conditions under which the mean of a sufficiently large number of independent random variables, each with finite mean and variance, will be approximately normally distributed. The use of Monte Carlo transport codes, such as Tripoli4, relies on those conditions. While these are verified in protection applications (the cycles provide independent measurements of fluxes and related quantities), the hypothesis of independent estimates/cycles is broken in criticality mode. Indeed the power iteration technique used in this mode couples a generation to its progeny. Often, after what is called 'source convergence' this coupling almost disappears (the solution is closed to equilibrium) but for loosely coupled systems, such as for PWR or large nuclear cores, the equilibrium is never found, or at least may take time to reach, and the variance estimation such as allowed by the CLT is under-evaluated. In this paper we first propose, by the mean of two different methods, to evaluate the typical correlation length, as measured in cycles number, and then use this information to diagnose correlation problems and to provide an improved variance estimation. Those two methods are based on Fourier spectral decomposition and on the lag k autocorrelation calculation. A theoretical modeling of the autocorrelation function, based on Gauss-Markov stochastic processes, will also be presented. Tests will be performed with Tripoli4 on a PWR pin cell. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

New Features of TRIPOLI-4{sup R} Version 10

TRIPOLI-4{sup R} is a three-dimensional and continuous-energy Monte-Carlo transport code developed by the 'Service d'Etudes des Reacteurs et de Mathematiques Appliquees' (SERMA), at the CEA Paris-Saclay Center. It is devoted to shielding, reactor physics, criticality safety and nuclear instrumentation. TRIPOLI-4{sup R} has been developed from the mid 1990's in C++. A new version of the code is typically released every two years. In this paper, we present the new features of TRIPOLI- 4{sup R} version 10, to be soon released to NEA and RSICC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Analysis of Radiation Effects in Silicon using Kinetic Monte Carlo Methods

The transient degradation of semiconductor device performance under irradiation has long been an issue of concern. Neutron irradiation can instigate the formation of quasi-stable defect structures, thereby introducing new energy levels into the bandgap that alter carrier lifetimes and give rise to such phenomena as gain degradation in bipolar junction transistors. Normally, the initial defect formation phase is followed by a recovery phase in which defect-defect or defect-dopant interactions modify the characteristics of the damaged structure. A kinetic Monte Carlo (KMC) code has been developed to model both thermal and carrier injection annealing of initial defect structures in semiconductor materials. The code is employed to investigate annealing in electron-irradiated, p-type silicon as well as the recovery of base current in silicon transistors bombarded with neutrons at the Los Alamos Neutron Science Center (LANSCE) “Blue Room” facility. Our results reveal that KMC calculations agree well with these experiments once adjustments are made, within the appropriate uncertainty bounds, to some of the sensitive defect parameters.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Development of a BlueCRAB/MELCOR Framework for Supporting Realistic Mechanistic Source Term Calculations in Microreactors

Efforts are currently underway to deploy microreactor modeling and simulation tools to better support vendors and regulatory authorities in submitting and reviewing licensing applications. In particular, the Nuclear Regulatory Commission is expected to rely on the Comprehensive Reactor Analysis Bundle (BlueCRAB) multiphysics toolset in performing design- and beyond-design-basis accident analyses. In addition, the Nuclear Regulatory Commission has been using the MELCOR code to estimate mechanistic source terms during accidents. As MELCOR relies on isotopic inventory and reactor temperature/power evolution profiles during accident conditions—all of which can theoretically be obtained from BlueCRAB—the ultimate goal of this activity is to establish a common BlueCRAB-MELCOR framework. However, prior to the present research, BlueCRAB had never been used to calculate such quantities of interest at the full-core level. While there are many Monte Carlo (MC) codes capable of computing such quantities of interest, they are unable to readily account for multiphysics feedback. BlueCRAB allows for the coupling of different physics codes together to perform multiphysics-informed calculations. Therefore, the purpose of this fiscal year 2023 work is to investigate the feasibility and challenges of performing such calculations within BlueCRAB so as to generate the data that MELCOR relies on. To demonstrate the methodology, the proposed workflow was applied to a prototypical heat pipe-cooled microreactor model. To predict isotopic concentrations (taking into account the ac- cumulation of fission products during operation), the necessary microscopic cross sections were generated via OpenMC and tabulated with respect to temperature and burnup. Next, a recently developed capability in Griffin (the reactor physics application in BlueCRAB) was used to convert the OpenMC output format into the ISOXML format used by Griffin. A multiphysics microscopic depletion calculation that involved performing a coupled full-core, heterogeneous neutron trans- port and thermal calculation at each depletion step was conducted to deplete the core to end of life (EOL) conditions so as to provide both isotopics and the initial condition for the transient calculation. Following a brief null-transient to verify that the initial condition had been properly restarted and was indeed in thermal equilibrium, a heat pipe failure transient was simulated. Thus, the entire workflow of using BlueCRAB to generate MELCOR inputs, from cross-section generation to producing isotopic inventory and power/temperature evolution profiles during transients, is demonstrated. This report also details the identified gaps in the workflow and how they were (for the most part) addressed. Future work should focus on directly including MEL- COR into the workflow by performing a MELCOR calculation using the BlueCRAB-generated input data. In addition, the heat pipe reactor design should be improved so as to reflect more prototypical burnup characteristics at EOL.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

A Code of Radiation Damage by Gamma - Paper 62

In recent years, the study of radiation damage of materials has seen a thrust because of the expectation of the high rate of radiation damage in fusion reactors, accelerator driven sub critical systems and damage of a control system of space vehicles by the high energy neutrons and gamma showers. Embrittlement of reactor vessels, corrosion of steel and swelling of fuel rods are the well known empirical examples of severe radiation damage. It is also realized that damage by high energy gammas is in no way insignificant, and it is equivalent to the damage by the neutrons. Secondly, value addition of certain materials like enhancement of electrical resistance, mechanical strength, color enhancement of precious stones on irradiation by neutrons and gammas and polymerization in the presence of gamma radiation add to the necessity of a Monte Carlo (MC) code of radiation damage in a wide range of energy. The JAIPU code, for the study of radiation damage by neutrons up to 10 MeV energy, incorporates features like development of secondary atomic cascade and PKA energy corrected for the ionization lose and production of Frenkel pair. MC estimation of both the primary and secondary displacements has been compared with the Frenkel pairs calculated using the NRT model. In the present work, JAIPU code is extended for irradiation of materials by high energy gammas. Here cascade generated by the incident gamma produces electrons by way of photo-electric, Compton and pair production until it is absorbed. The recoiled Frenkel pairs produced in electron + atom collision have been considered for the purpose of atomic cascade. Standard cross section libraries for gamma and electron interaction are used while for atom + atom collision we have made use of the IOTA code. In addition, in the code, the lattice structure of the irradiated material is simulated to know the positions of the photon interactions with atoms and the energy of electron it adds an extra feature of the position and size of the voids and the clusters produced on irradiation. The code has the possibility of calculation of 'absorbed radiation dose' and 'KERMA' and the interatomic distances of the displaced atoms. It may be pointed out that the code uses only the MC simulation of the number of displacements unlike several other approaches where the NRT model with PKA energy is used for gross estimation of the number of the displaced atoms. Thus, the code keeps an account of the positions of displaced atoms for the in-depth modification of the irradiated material. 1 MeV gamma seems to damage Fe up to several millimeters of its depth and this is mostly by the recoiled Compton electrons than photo electrons. The code is in the developmental stage and it has the capability of applications in the developing field of Radiology besides the shielding and other studies of radiation damage. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2011

Calculation of k{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}} for several small detectors and for two linear accelerators using Monte Carlo simulations

Purpose: The scope of this study was to determine a complete set of correction factors for several detectors in static small photon fields for two linear accelerators (linacs) and for several detectors. Methods: Measurements for Monte Carlo (MC) commissioning were performed for two linacs, Siemens Primus and Elekta Synergy. After having determined the source parameters that best fit the measurements of field specific output factors, profiles, and tissue-phantom ratio, the generalized version of the classical beam quality correction factor for static small fields, k{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}}, were determined for several types of detectors by using the egs{sub c}hamber Monte Carlo user code which can accurately reproduce the geometry and the material composition of the detector. The influence of many parameters (energy and radial FWHM of the electron beam source, field dimensions, type of accelerator) on the value of k{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}} was evaluated. Moreover, a MC analysis of the parameters that influence the change of k{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}} as a function of field dimension was performed. A detailed analysis of uncertainties related to the measurements of the field specific output factor and to the Monte Carlo calculation of k{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}} was done. Results: The simulations demonstrated that the correction factor k{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}} can be considered independent from the quality beam factor Q in the range 0.68 {+-} 0.01 for all the detectors analyzed. The k{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}} of PTW 60012 and EDGE diodes can be assumed dependent only on the field size, for fields down to 0.5 x 0.5 cm{sup 2}. The microLion, and the microchambers, instead, must be used with some caution because they exhibit a slight dependence on the radial FWHM of the electron source, and therefore, a correction factor only dependent on field size can be used for fields {>=}0.75 x 0.75 and {>=}1.0 x 1.0 cm{sup 2}, respectively. The analysis of uncertainties gave an estimate of uncertainty for the 0.5 x 0.5 cm{sup 2} field of about 0.7% (1{sigma}) for k{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}} factor and of about 1.0% (1{sigma}) for the field output factor, {Omega}{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}}, of diodes, microchambers, and microLion. Conclusions: Stereotactic diodes with the appropriate k{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}} are recommended for determining {Omega}{sub Q{sub c{sub l{sub i{sub n,Q{sub m{sub s{sub r}{sup f{sub c}{sub l}{sub i}{sub n},f{sub m}{sub s}{sub r}}}}}}}}} of small photon beams.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Visualizing expanding warm dense matter heated by laser-generated ion beams

This PowerPoint presentation concluded with the following. We calculated the expected heating per atom and temperatures of various target materials using a Monte Carlo simulation code and SESAME EOS tables. We used aluminum ion beams to heat gold and diamond uniformly and isochorically. A streak camera imaged the expansion of warm dense gold (5.5 eV) and diamond (1.7 eV). GXI-X recorded all 16 x-ray images of the unheated gold bar targets proving that it could image the motion of the gold/diamond interface of the proposed target.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optimization of particle tracking methods for stochastic media

Random media emerge in several applications involving particle transport, encompassing e.g. photon propagation through Rayleigh-Taylor instabilities in fuel pellets for inertial confinement fusion, or neutron multiplication problems related to the assessment of re-criticality risk following severe accidents with fuel degradation. Reference calculations in such material configurations by means of Monte Carlo transport codes are particularly challenging, since high-density stochastic media might involve several hundreds of thousands of volumes and thus make particle tracking routines extremely cumbersome. In order to cope with these issues, two distinct strategies have been proposed so far: the use of neighbor maps, or the use of delta tracking. In this work we will compare these methods and illustrate their specific merits and drawbacks, as taken both alone and in combination with each other. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2012

Effects of adsorption and roughness upon the collision processes at the convertor surface of a plasma sputter negative ion source

Atomic collision processes associated with surface production of negative hydrogen ions (H{sup -}) by particle reflection at molybdenum surface immersed in hydrogen plasma have been investigated. To calculate sputtering yields of Cs, as well as energy spectra and angular distributions of reflected hydrogen atoms from molybdenum surface by H{sup +} ion and Cs{sup +} ion bombardments, a Monte Carlo simulation code ACAT (Atomic Collision in Amorphous Target) was run with the corresponding surface conditions. A fractal surface model has been developed and adopted to ACAT for evaluating the effect due to roughness of target material. The results obtained with ACAT have indicated that the retention of hydrogen atoms leads to the reduction in sputtering yields of Cs, and the surface roughness does largely affect the sputtering yields of Cs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2003

A new dynamical atmospheric ionizing radiation (AIR) model for epidemiological studies

A new Atmospheric Ionizing Radiation (AIR) model is currently being developed for use in radiation dose evaluation in epidemiological studies targeted to atmospheric flight personnel such as civilian airlines crewmembers. The model will allow computing values for biologically relevant parameters, e.g. dose equivalent and effective dose, for individual flights from 1945. Each flight is described by its actual three dimensional flight profile, i.e. geographic coordinates and altitudes varying with time. Solar modulated primary particles are filtered with a new analytical fully angular dependent geomagnetic cut off rigidity model, as a function of latitude, longitude, arrival direction, altitude and time. The particle transport results have been obtained with a technique based on the three-dimensional Monte Carlo transport code FLUKA, with a special procedure to deal with HZE particles. Particle fluxes are transformed into dose-related quantities and then integrated all along the flight path to obtain the overall flight dose. Preliminary validations of the particle transport technique using data from the AIR Project ER-2 flight campaign of measurements are encouraging. Future efforts will deal with modeling of the effects of the aircraft structure as well as inclusion of solar particle events. Published by Elsevier Ltd on behalf of COSPAR.

Aviation↗
DOE OSTI2015

Photon/Electron Benchmarks for Intercode Comparisons

The goal of this work was to improve accuracy and efficiency of two Monte-Carlo transport codes (MCNP and DIANE) with an emphasis on γ+electron physics. The approach involved intercode comparisons + measurements for gamma/e - energy deposition in a cylinder with a photon source and different materials (C, Pb) and the bombardment of 15-MeV electrons on thick targets (Al, Be, Pb). Comparisons of the codes DIANE and MCNP6 showed good agreement (differences < 3%) for gamma-electron energy deposition in a 2D cylinder, except for the first 0.1 μm of lead (difference < 10%). Comparisons with measurements showed generally good agreement, often better than 10%; best-performing codes/options are problem-dependent; and single-event discrepancies are in active use in reviewing electron elastic scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

IK-Frag: Frag data generator for the PHITS simulation with the inverse kinematic reaction producing a focused neutron beam

IK-Frag has been developed for the creation of the nuclear cross-section data format, which is named frag data and can be used in PHITS, a Monte Carlo simulation code. IK-Frag focuses on the inverse kinematic reactions between a lithium or beryllium ion and a proton target. These reactions achieve naturally collimated neutrons and potentially reduce the necessity of radiation shielding. IK-Frag enables PHITS users to conduct simulations for the inverse kinematic reactions. The present software aims to contribute to future development of the neutron source system using the inverse kinematic reactions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2013

Improved convergence of Monte Carlo generated multi-group scattering moments

This paper introduces an improved method of obtaining multi-group scattering moments from a Monte Carlo neutron transport code for use in deterministic transport solvers. The new method increases the information obtained from scattering events and therefore has more useful convergence characteristics than the currently used analog techniques. A prototype of the improved method was implemented in the OpenMC Monte Carlo transport code to compare the accuracy and convergence characteristics of the new method. The prototype showed that accuracy was retained (or improved) while increasing the figure-of-merit for the generation of multi-group scattering moments. (authors)

97 MATHEMATICS AND COMPUTING↗