Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Code”

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

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

At least 1,153 records · Page 64

Methodology and Software for Gross Defect Detection of Spent Nuclear Fuel at the Atucha-I Reactor [Novel Methodology and Software for Spent Fuel Gross Defect Detection at the Atucha-I Reactor]

Here, fuel assemblies in the spent fuel pool are stored by suspending them in two vertically stacked layers at the Atucha Unit 1 nuclear power plant (Atucha-I). This introduces the unique problem of verifying the presence of fuel in either layer without physically moving the fuel assemblies. Given that the facility uses both natural uranium and slightly enriched uranium at 0.85 wt% 235 U and has been in operation since 1974, a wide range of burnups and cooling times can exist in any given pool. A gross defect detection tool, the spent fuel neutron counter (SFNC), has been used at the site to verify the presence of fuel up to burnups of 8000 MWd/t. At higher discharge burnups, the existing signal processing software of the tool was found to fail due to nonlinearity of the source term with burnup.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Strategies for thorium fuel cycle transition in the SD-TMSR

Liquid-fueled Molten Salt Reactor (MSR) systems represent advances in safety, economics, and sustainability. The MSR has been designed to operate with a Th/ 233 U fuel cycle with 233 U used as startup fissile material. Since 233 U does not exist in nature, we must examine other available fissile materials to start up these reactor concepts. This work investigates the fuel cycle and neutronics performance of the Single-fluid Double-zone Thorium-based Molten Salt Reactor (SD-TMSR) with different fissile material loadings at startup: High Assay Low Enriched Uranium (HALEU) (19.79%), Pu mixed with HALEU (19.79%), reactor-grade Pu (a mixture of Pu isotopes chemically extracted from Pressurized Water Reactor (PWR) spent nuclear fuel (SNF) with 33 GW d/tHM burnup), transuranic elements (TRU) from Light Water Reactor (LWR) SNF, and 233 U. The MSR burnup routine provided by SERPENT-2 is used to simulate the online reprocessing and refueling in the SD-TMSR. The effective multiplication factor, fuel salt composition evolution, and net production of 233 U are studied in the present work. Additionally, the neutron spectrum shift during the reactor operation is calculated. The results show that the continuous flow of reactorgrade Pu helps transition to the thorium fuel cycle within a relatively short time (≈ 4.5 years) compared to 26 years for 233 U startup fuel. Finally, using TRU as the initial fuel materials offers the possibility of operating the SD-TMSR for an extended period of time (≈ 40 years) without any external feed of 233 U.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Photonuclear production of 47 Ca for 47 Ca/ 47 Sc generator from natural CaCO 3 targets

Here this work investigated the indirect production of 47 Sc from natural Ca targets via 48 Ca(γ,n) 47 Ca → 47 Sc +β - +$\bar{\nu}_e$ with incident electron energies of 30, 35, and 40 MeV. The 47 Ca production yields were simulated using the PHITS Monte Carlo simulation code and compared to experimental data. The simulated production rates for all three irradiations are in good agreement with experimental data within uncertainties. As a demonstration of the 47 Ca/ 47 Sc generator system, one of the irradiated CaCO 3 targets was dissolved in nitric acid, and 47 Sc was isolated from the target material using commercially available Eichrom DGA resin. The 47 Sc was allowed to grow in, and the purification process was repeated with promising 47 Sc and Ca recovery yields.

07 ISOTOPE AND RADIATION SOURCES↗

ENDF/B-VII.1 Neutron Cross Section Data Testing with Critical Assembly Benchmarks and Reactor Experiments

The ENDF/B-VII.1 library is the latest revision to the United States Evaluated Nuclear Data File (ENDF). The ENDF library is currently in its seventh generation, with ENDF/B-VII.0 being released in 2006. This revision expands upon that library, including the addition of new evaluated files (was 393 neutron files previously, now 423 including replacement of elemental vanadium and zinc evaluations with isotopic evaluations) and extension or updating of many existing neutron data files. Complete details are provided in the companion paper [M. B. Chadwick et al., 'ENDF/B-VII.1 Nuclear Data for Science and Technology: Cross Sections, Covariances, Fission Product Yields and Decay Data,' Nuclear Data Sheets, 112, 2887 (2011)]. This paper focuses on how accurately application libraries may be expected to perform in criticality calculations with these data. Continuous energy cross section libraries, suitable for use with the MCNP Monte Carlo transport code, have been generated and applied to a suite of nearly one thousand critical benchmark assemblies defined in the International Criticality Safety Benchmark Evaluation Project's International Handbook of Evaluated Criticality Safety Benchmark Experiments. This suite covers uranium and plutonium fuel systems in a variety of forms such as metallic, oxide or solution, and under a variety of spectral conditions, including unmoderated (i.e., bare), metal reflected and water or other light element reflected. Assembly eigenvalues that were accurately predicted with ENDF/B-VII.0 cross sections such as unmoderated and uranium reflected {sup 235}U and {sup 239}Pu assemblies, HEU solution systems and LEU oxide lattice systems that mimic commercial PWR configurations continue to be accurately calculated with ENDF/B-VII.1 cross sections, and deficiencies in predicted eigenvalues for assemblies containing selected materials, including titanium, manganese, cadmium and tungsten are greatly reduced. Improvements are also confirmed for selected actinide reaction rates such as {sup 236}U, {sup 238,242}Pu and {sup 241,243}Am capture in fast systems. Other deficiencies, such as the overprediction of Pu solution system critical eigenvalues and a decreasing trend in calculated eigenvalue for {sup 233}U fueled systems as a function of Above-Thermal Fission Fraction remain. The comprehensive nature of this critical benchmark suite and the generally accurate calculated eigenvalues obtained with ENDF/B-VII.1 neutron cross sections support the conclusion that this is the most accurate general purpose ENDF/B cross section library yet released to the technical community.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Inventory Uncertainty Quantification using TENDL Covariance Data in Fispact-II

The new inventory code Fispact-II provides predictions of inventory, radiological quantities and their uncertainties using nuclear data covariance information. Central to the method is a novel fast pathways search algorithm using directed graphs. The pathways output provides (1) an aid to identifying important reactions, (2) fast estimates of uncertainties, (3) reduced models that retain important nuclides and reactions for use in the code's Monte Carlo sensitivity analysis module. Described are the methods that are being implemented for improving uncertainty predictions, quantification and propagation using the covariance data that the recent nuclear data libraries contain. In the TENDL library, above the upper energy of the resolved resonance range, a Monte Carlo method in which the covariance data come from uncertainties of the nuclear model calculations is used. The nuclear data files are read directly by FISPACT-II without any further intermediate processing. Variance and covariance data are processed and used by FISPACT-II to compute uncertainties in collapsed cross sections, and these are in turn used to predict uncertainties in inventories and all derived radiological data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurements and simulations of the cosmic-ray-induced neutron background

The cosmic-ray-induced neutron background at ground level has been measured and simulated in conjunction with EJ-309 organic liquid scintillators with an approximate deposited energy range of 0.5–6 MeV. Specifically, the pulse height distributions, net neutron count rates, and angular dependences were obtained. The simulations were carried out using the Monte Carlo transport code MCNPX-PoliMi combined with the (Cosmic-Ray Shower Generator) CRY source subroutine that returns secondary particles produced by cosmic rays. A scaling formula from literature was also implemented in the simulation. The angular dependence of the neutron count rate was measured by collimating the liquid scintillator with polyethylene to attain 18° angular resolution from 0° downwards to 72° horizontally. The neutron count rate was measured to be 23.10±1.69 h -1 sr -1 at 0°, and 7.20±0.78 h -1 sr -1 at 72°. The simulations and measurements compare well and show similar cosine anisotropy for the angular distribution. Finally, the study thus shows that the neutron background response in detector systems can be efficiently and accurately simulated using the procedures described.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SANDD: A directional antineutrino detector with segmented 6 Li-doped pulse-shape-sensitive plastic scintillator

We present a characterization of a small (9-liter) and mobile 0.1% 6 Li-doped pulse-shape-sensitive plastic scintillator antineutrino detector called SANDD (Segmented AntiNeutrino Directional Detector), constructed for the purpose of near-field reactor monitoring with sensitivity to antineutrino direction. SANDD comprises three different types of module. A detailed Monte Carlo simulation code was developed to match and validate the performance of each of the three modules. The combined model was then used to produce a prediction of the performance of the entire detector. Analysis cuts were established to isolate antineutrino inverse beta decay events while rejecting large fraction of backgrounds. The neutron and positron detection efficiencies are estimated to be 34.8% and 80.2%, respectively, while the coincidence detection efficiency is estimated to be 71.7%, resulting in inverse beta decay detection efficiency of 20.0% ± 0.2%(stat.) ± 2.1%(syst.). Finally, the predicted directional sensitivity of SANDD produces an uncertainty of 20° in the azimuthal direction per 100 detected antineutrino events.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Thick target neutron yields from Beryllium, Carbon, Tungsten, and Lead targets irradiated by 26.7 MeV/nucleon 4 H e ions

The angular and energy distributions of secondary neutrons produced from 26.7 MeV/nucleon 4 He ions, stopping in thick Be, C, W and Pb targets are measured by the time of flight method. The GEANT4, PHITS, and FLUKA Monte Carlo simulation codes with different physics models are employed to simulate the neutron yields and the results are compared with the experimental data. For the thick Be and C targets, GEANT4 with INCL++ reproduces the experimental data reasonably well. Whereas PHITS with JQMD-2.0 overestimates the neutron yields at intermediate energies and FLUKA with PEANUT underestimated the neutron yields at intermediate energies. For W and Pb targets, FLUKA with PEANUT physics model agree with experimental data best. PHITS with JQMD-2.0 overestimates the neutron yield at intermediate energies at 0°, and GEANT4 with INCL++ underestimates the neutron yield at intermediate energies at three angles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Shielding development for the spallation neutron source VENUS instrument

VENUS is a neutron imaging instrument that will use a broad range of neutron wavelengths, from epithermal to cold energies, and will include enhanced contrast mechanisms. It will offer novel energy-selective imaging techniques directly connecting complex engineering materials and systems’ structures, properties, and functions to reveal practical and fundamental answers about their real-world performance. Here, the instrument will be built at SNS beamline 10, facing the decoupled poisoned hydrogen moderator. The driving cost for the instrument is the beamline and instrument cave shielding. Final analyses were performed to evaluate the thickness and composition of shielding materials for the instrument cave and beamline to meet radiation safety criteria for the instrument to start up in 2024 after completing the SNS proton power upgrade.

43 PARTICLE ACCELERATORS↗

Low energy electrons and swift ion track structure in PADC

The current work aims at providing an accurate description of the ion track-structure in poly-allyl dyglycol carbonate (PADC) by using an up-to-date Monte-Carlo code-called TILDA-V (a French acronym for Transport d’Ions Lourds Dans l’Aqua & Vivo). In this simulation the ion track-structure in PADC is mainly described in terms of ejected electrons with a particular attention done to the Low Energy Electrons (LEEs). After a brief reminder of the most important channels through which LEEs are prone to break a chemical bond, we will report on the simulated energetic distributions of LEEs along an ion track in PADC for particular incident energies located on both sides of the Bragg-peak position. Lastly, based on the rare data dealing with LEEs interaction with polymers or organic molecules, we will emphasise the role played by the LEEs in the formation of a latent track in PADC, and more particularly the one played by the sub-ionization electrons.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Feasibility of a multigroup Boltzmann–Fokker–Planck solution for electron beam dose calculations

Legacy nuclear-reactor Boltzmann solvers start clinical deployment as an alternative to Monte Carlo (MC) codes and Fermi–Eyges semiemprical models in radiation oncology treatment planning. Today’s certified clinical solvers are limited to photon beams. In this paper, ELECTR, a state-of-the-art multigroup electron cross sections generation module in NJOY is presented and validated against Lockwood’s calorimetric measurements, EGS-nrc and GEANT-4 for 1–20 MeV unidirectional electron beams. The nuclear-reactor DRAGON-5 solver is upgraded to access the library and solve the Boltzmann–Fokker–Planck (BFP) equation. A variety of heterogeneous radiotherapy and radiosurgery phantom configurations were used for validation purpose. Case studies include a thorax benchmark, that of a typical breast Intra-Operative Radiotherapy and a high-heterogeneity patient-like benchmark. For all beams, 100% of the water voxels satisfied the American Association of Physicists in Medicine accuracy criterion for a BFP-MC dose error below 2%. At least, 97.0% of adipose, muscle, bone, lung, tumor and breast voxels satisfied the 2% criterion. The average BFP-MC relative error was about 0.56% for all voxels, beams and materials combined. By irradiating homogeneous slabs from Z = 1 (hydrogen) to Z = 99 (einsteinium), we reported performance and defects of the CEPXS mode [US. Sandia National Lab., SAND-89-1685] in ELECTR for the entire periodic table. For all Lockwood’s benchmarks, NJOY-DRAGON dose predictions are within the experimental data precision for 98% of voxels.

42 ENGINEERING↗

Characteristics of the Neutron Irradiation Facilities of the PSI Calibration Laboratory

The neutron radiation fields of the Calibration Laboratory at Paul Scherrer Institute (PSI) are traceable to the national standards of the Physikalisch-Technische Bundesanstalt (PTB) in Germany. A Berthold LB6411 neutron dose rate meter for neutron radiation is used as a secondary standard. Recently, a thorough characterization of the neutron irradiation fields of the {sup 241}Am-Be and {sup 252}Cf sources by means of reference measurements and a detailed MCNPX simulation of the irradiation facility has been initiated. In this work, the characteristics of the neutron radiation fields are summarized and presented together with model equations and an uncertainty analysis. MCNPX results are shown for the {sup 241}Am-Be source. A comparison of measured and simulated data shows an excellent agreement. From the simulation, valuable information about the neutron fields like the contribution of scattered neutrons in the fields and the energy spectra could be obtained.

07 ISOTOPE AND RADIATION SOURCES↗

Heavy reflector experiments in the IPEN/MB-01 reactor: Stainless steel, carbon steel and nickel

New experiments devoted to the measurements of physical parameters of a light water core surrounded by a heavy reflector were performed in the IPEN/MB-01 research reactor facility. These experiments comprise three sets of heavy reflector (SS-304, Carbon Steel, and Nickel) in a form of laminates around 3 mm thick. Each set was introduced individually in the west face of the core of the IPEN/MB-01 reactor. The aim here is to provide high quality experimental data for the interpretation and validation of the SS-304 heavy reflector calculation methods. The experiments of Carbon Steel, which is composed mainly of iron, and Nickel were performed to provide a consistent and an interpretative check for the SS-304 reflector experiment. The experimental results comprise critical control bank positions, temperatures and reactivities as a function of the number of the plates. Particularly to the case of Nickel, the experimental data are unique of its kind. The theoretical analysis was performed by MCNP-5 with the nuclear data library ENDF/B-VII.0. It was shown that this nuclear data library has a very good performance up to thirteen plates and overestimates the reactivity for higher number of plates independently of the type of the reflector.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Reactivity worth measurements at the IPEN/MB-01 nuclear reactor

Researches that aim to improve the performance of neutron transport codes and quality of nuclear cross section databases are very important to increase the accuracy of simulations and the quality of the analysis and prediction of phenomena in the nuclear field. In this context, relevant experimental data such as reactivity worth measurements are needed. The objective of this work was to perform a series of experiments of reactivity worth measurements, using a digital reactivity meter developed at IPEN. The experiments employed small metallic and ceramic samples inserted in the central region of the core of the experimental IPEN/MB-01 reactor. The theoretical analysis was performed by the MCNP-5 reactor physics code, developed and maintained by Los Alamos National Laboratory, and the ENDF/B-VII.0 nuclear data library.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Energy conversion and transfer for plasmas in a magnetic expansion configuration

A two-dimensional axisymmetric particle-in-cell code with Monte Carlo collision conditions has been used to study particle energy transfer in plasmas and conversion in applied magnetic and electric fields appropriate to coaxial acceleration. The research incorporates a computation scheme with: a model of single particle magnetic interactions; a model of single particle interactions in electric and magnetic fields; and a model of multi-particle collisional interactions in order to understand the energy transfer processes and conversion mechanisms of charged plasma particles. This approach predicts electron and ion motions along with their energy variations for physical conditions that occur in the related models; the results allow comparison with experimental data for magnetic field strengths of 0.01–0.05 T and electrode voltages of 22.0–32.0 V. With the incorporation of magnetic and electric field effects on charged particles, the multi-particle model includes electron-neutral ionization collisions, ion-neutral charge exchange collisions, and electron-ion Coulomb collisions. This research presents a new approach to achieve an underlying understanding of the plasma energy transfer and conversion in the external electric and magnetic fields that is not possible using magnetohydrodynamics continuum representations. Results indicate the following innovative conclusions: (1) Radial and azimuthal energies of magnetized electrons are converted into an axial electron energy component in the diverging magnetic field, and the azimuthal kinetic energy of unmagnetized ions is converted into axial and radial components. (2) In electric and magnetic fields, electric field energy is primarily converted into axial kinetic energy of magnetized electrons by the energy transformation effects of magnetic fields, and for unmagnetized ions, the radial kinetic energy component dominates in the conversion of electric field energy. (3) For the collisional plasma, electron kinetic energy tends to increase (or decrease) to a terminal value since electrons lose energy in collisions then gain energy again from the field. Ions acquire most energy directly from the electric field, although part of the electric field energy arrives to the ions by collisions. Further, the ion axial energy component dominates the total ion energy. The collision processes are found to be integral and essential for the conversion of the plasma non-directed energy gain to be converted into the resultant axial energy, the magnitudes of which are found to be in agreement with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Control rod calibration and reactivity effects at the IPEN/MB-01 reactor

Researches that aim to improve the performance of neutron transport codes and quality of nuclear cross section databases are very important to increase the accuracy of simulations and the quality of the analysis and prediction of phenomena in the nuclear field. In this context, relevant experimental data such as reactivity worth measurements are needed. Control rods may be made of several neutron absorbing materials that are used to adjust the reactivity of the core. For the reactor operation, these experimental data are also extremely important: with them it is possible to estimate the reactivity worth by the movement of the control rod, understand the reactor response at each rod position and to operate the reactor safely. This work presents a temperature correction approach for the control rod calibration problem. It is shown the control rod calibration data of the IPEN/MB-01 reactor, the integral and differential reactivity curves and a theoretical analysis, performed by the MCNP-5 reactor physics code, developed and maintained by Los Alamos National Laboratory, using the ENDF/B-VII.0 nuclear data library.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Design considerations for three-dimensional betavoltaics

Betavoltaic devices are suitable for delivering low-power over periods of years. Typically, their power density is on the order of nano to micro-Watts per cubic centimeter. In this work we evaluate the potential for using high-aspect ratio three-dimensional semiconductor structures to enhance the power and efficiency of these devices. The Monte Carlo transport code MCNP6 is used to provide realistic estimates of the theoretical levels of charge generation, which is in turn used to make predictions about the power output from three-dimensional betavoltaics. The focus of this work is on silicon and promethium-147, but other semiconductors and radioisotopes are considered as well. In the case of silicon diodes with three-dimensional features that are comparable to what is commercially available we estimate that power densities in the range of 20-25 mW/cm 3 can be achieved at efficiencies of 2.9-5.8% when coupled with promethium-147 oxide.

25 ENERGY STORAGE↗

Dependence of high- Z redeposition on the field-to-surface pitch angle and other sheath parameters in tokamaks

Accurately predicting redeposition is vital for high-Z plasma-facing component (PFC) survivability in magnetic confinement fusion. In this study, we categorize high-Z redeposition into three mechanisms: geometric-driven (prompt), sheath-driven (local), and scrape-off-layer-driven (far) redeposition. To investigate these mechanisms, we employ Monte Carlo transport codes to simulate azimuthally symmetric tungsten source erosion and redeposition in a tokamak. By iteratively analyzing critical parameters, we evaluate redeposition scaling for each mechanism. Specifically, we investigate the impact of magnetic-field-to-PFC pitch angle assumptions on PFC losses into the scrape-off layer. Our findings reveal significant pitch angle sensitivity due to an asymmetric prompt vs local redeposition trade-off. These results enhance our understanding of redeposition phenomena in fusion plasma environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗