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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Two-Step Procedure for Liquid-Salt-Cooled-Reactor Analysis

Here, this paper presents the development of a lattice physics–to–core simulator two-step procedure for the rapid analysis of the Advanced High Temperature Reactor (AHTR). Lattice physics, reflector, and control blade models were developed from which cross-section libraries could be generated for a nodal core simulator. Few-group structures for the core simulator were also generated to account for the neutronic characteristics of AHTR. After developing the AHTR two-step procedure, cross-section libraries were generated using the SERPENT continuous-energy Monte Carlo code. These libraries were then used in the core simulator NESTLE to perform full-core calculations, which were in turn benchmarked against reference SERPENT full-core models. Benchmarking results showed reasonable accuracy of the developed two-step procedure but revealed an inherent inadequacy in the one-dimensional radial reflector model and showed a likely need for a greater number of energy groups than were used in this study.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Modeling of the Molten Salt Reactor Experiment with SCALE

A SCALE model was developed for the Molten Salt Reactor Experiment (MSRE) benchmark that was recently added to the International Handbook of Evaluated Reactor Physics Benchmark Experiments. This SCALE model served as a basis for criticality calculations and nuclear data sensitivity and uncertainty analyses with the Monte Carlo code Shift and the TSUNAMI computational capabilities in the SCALE code system. The focus of this work is the assessment of the impact of nuclear data on the calculated eigenvalue results in support of the discussion of differences between the calculated and the experimental eigenvalue result. The differences in the eigenvalues obtained using the ENDF/B-VII.0, ENDF/B-VII.1, and ENDF/B-VIII.0 nuclear data libraries cover a relatively small range of ~230 pcm. Since eigenvalue sensitivity of the MSRE is dominated by the neutron multiplicity and neutron capture of 235 U and elastic scattering in graphite, relevant changes in the ENDF/B libraries for nuclear reactions (such as carbon capture) that caused large differences in other graphite-moderated systems did not have a significant impact. Propagation of nuclear data uncertainty results in an eigenvalue uncertainty of ~700 pcm with the major contributors being 235 U neutron multiplicity, graphite elastic scattering, and 7Li neutron capture. All calculations resulted in large differences of ~2000 pcm in eigenvalue compared to the benchmark experimental value. Several potential contributors to this difference—including uncertainties and gaps in the knowledge of the material, geometry, and nuclear data—were identified. Simplified models of the full MSRE core were developed, and similarity assessments were conduced with the full MSRE core model. It was found that simplified models can serve as adequate surrogates of the full-core model such that they can be used for performing selected nuclear data performance assessments with a lower computational burden.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Benchmark of Neutron Thermalization in Graphite Using a Pulsed Slowing-Down-Time Experiment

A benchmark has been developed using a pulsed slowing-down-time experiment to isolate the thermalization process in graphite. The experiment was conducted at the Oak Ridge Electron Linear Accelerator facility at Oak Ridge National Laboratory, and it measured the time spectrum of neutrons leaking from a graphite pile during slowing down and thermalization within graphite. Simulations of the benchmark experiment were performed using the MCNP6.1 Monte Carlo code and the ENDF/B-VII.1 and ENDF/B-VIII.0 cross-section databases. The benchmark provides a time spectrum (i.e., time-dependent counts in a detector) that allows for validation of the graphite thermal scattering libraries (TSLs). The impact on the simulations using a suite of graphite TSLs was compared with the experimental results. Given the density of nuclear graphite, the TSL corresponding to graphite with 30% porosity, as implemented in ENDF/B-VIII.0, was found to most accurately represent the measured time spectrum corresponding to the thermal energy range with an average deviation of ±1.7%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Kinetic simulation of a low-pressure helium discharge with comparison to experimental measurements

Modern computational validation efforts rely on comparison of known experimental quantities such as current, voltage, particle densities, and other plasma properties with the same values determined through simulation. A discrete photon approach for radiation transport was recently incorporated into a particle-in-cell/direct simulation Monte Carlo code. As a result, spatially and temporally resolved synthetic spectra may be generated even for non-equilibrium plasmas. The generation of this synthetic spectra lends itself to potentially new validation opportunities. In this work, initial comparisons of synthetic spectra are made with experimentally gathered optical emission spectroscopy. A custom test apparatus was constructed that contains a 0.5 cm gap distance parallel plane discharge in ultra high purity helium gas (99.9999%) at a pressure of 75 Torr. Plasma generation is initiated with the application of a fast rise-time, 100 ns full-width half maximum, 2.0 kV voltage pulse. Transient electrical diagnostics are captured along with time-resolved emission spectra. A one-dimensional simulation is run under the same conditions and compared against the experiment to determine if sufficient physics are included to model the discharge. To sync the current measurements from experiment and simulation, significant effort was undertaken to understand the kinetic scheme required to reproduce the observed features. Additionally, the role of the helium molecule excimer emission and atomic helium resonance emission on photocurrent from the cathode are studied to understand which effect dominates photo-feedback processes. Results indicate that during discharge development, atomic helium resonance emission dominates the photo-flux at the cathode even though it is strongly self-absorbed. A comparison between the experiment and simulation demonstrates that the simulation reproduces observed features in the experimental discharge current waveform. Here, the synthesized spectra from the kinetic method produces more favorable agreement with the experimental data than a simple local thermodynamic equilibrium calculation and is a first step towards using spectra generated from a kinetic method in validation procedures. The results of this study produced a detailed compilation of important helium plasma chemistry reactions for simulating transient helium plasma discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Determination of the mean energy of fast electron losses and anisotropies through thick-target emission on WEST

Abstract A new method to obtain the mean energy of fast electron losses in fusion plasmas using a versatile multi-energy hard x-ray (HXR) detector is presented. The method is based on measuring the thick-target emission of tungsten in the divertor region produced by fast electron losses interacting with the target and modeling the tungsten spectra by a Monte Carlo code which simulates the interaction between a beam of electrons and a solid target. The mean energy of the fast electron losses is determined through the comparison between the experimental and synthetic emission. The results show that fast electron losses during lower hybrid current drive discharges at WEST have a mean energy of 90–140 keV and represent only 2% of the total heat flux at the target. Additionally, anisotropic HXR emission has been detected for the first time at the WEST core and edge plasma, with opposite directions. It is due to the forward-peak emission of two distinctive populations of fast electrons: co-current fast electrons in the core and counter-current fast electron losses at the inner strike point. In view of future experiments like ITER where electron cyclotron current drive will generate a fast electron population, this technique could serve as a real-time monitor of fast electron losses and eventually feed an actuator on the current drive generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Validation of a synthetic fast ion loss detector model for Wendelstein 7-X

Abstract We present the first validated synthetic diagnostic for fast ion loss detectors (FILDs) in the Wendelstein 7-X (W7-X) stellarator. This model has been developed on, and validated against experimental data from, a FILD provided by the National Institute for Fusion Science (NIFS-FILD), with potential future applicability to the existing Faraday Cup FILD (FC-FILD) on W7-X as well as the scintillating FILD (S-FILD) currently under development. A workflow combining Monte Carlo codes BEAMS3D and ASCOT5 is used to track fast ions produced by neutral beam injection from the moment of ionization until they are thermalized or lost from the last closed flux surface, and from there to a virtual plane which serves as a projection of the entrance aperture to the FILD. Simulations in ASCOT5 are analyzed via a geometric method to determine the probability of transmission through the FILD aperture and onto the detector as a function of normalized momentum, pitch angle, gyrophase, and position at the virtual plane. This probability is then applied to the simulated ions arriving from the plasma, producing a simulated signal from a computationally tractable number of simulated fast ions. Simulated signals are presented for two W7-X experiments with neutral beam injection and quantitatively compared with experimental measurements from the NIFS-FILD diagnostic. An estimate of the frequency of charge-exchange with neutral particles in the edge is performed, and it is found that this process may have a significant impact on the measured signals.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)↗
DOE OSTI2025

Measurement of the mean excitation energy of liquid argon

The mean excitation energy (I-value) of liquid argon is a critical input for energy estimation in neutrino oscillation experiments. It is measured to be 199 +4 -5 eV using the range of 402.2 MeV protons from the Fermilab Linac. This compares to the author's recent evaluation of (197 ± 7) eV based on a combination of an oscillator strength distribution analysis, gaseous argon range measurements, sparse stopping power data on solid argon, and an extrapolation of data on the effect of phase from other substances. Using all sources of information, we recommend a value of (198 ± 4) eV for liquid argon, which is significantly higher than 188 eV, from ICRU-37's gaseous argon evaluation, commonly used in Monte Carlo codes such as Geant4.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A low-luminosity core-collapse supernova very similar to SN 2005cs

ABSTRACT We present observations and analysis of PSN J17292918+7542390, a low-luminosity Type II-P supernova (LL SN IIP). The observed sample of such events is still low, and their nature is still under debate. Such SNe are similar to SN 2005cs, a well-observed LL Type II-P event, having low expansion velocities, and small ejected 56Ni mass. We have developed a robust and relatively fast Monte Carlo code that fits semi-analytic models to light curves of core-collapse SNe. This allows the estimation of the most important physical parameters, like the radius of the progenitor star, the mass of the ejected envelope, the mass of the radioactive nickel synthesized during the explosion, among others. PSN J17292918+7542390 has $R_0 = 91_{-70}^{+119} \times 10^{11} \, \text{cm}$, $M_\text{ej} = 9.89_{-1.00}^{+2.10} \, \mathrm{ M}_{\odot }$, $E_{\mbox{kin}} = 0.65_{-0.18}^{+0.19} \, \text{foe}$, and $v_{\mbox{exp}} = 3332_{-347}^{+216}$ km s−1, for its progenitor radius, ejecta mass, kinetic energy, and expansion velocity, respectively. The initial nickel mass of the PSN J17292918+7542390 turned out to be $1.55_{-0.70}^{+0.75} \times 10^{-3} \, \mathrm{M}_{\odot }$. The measured photospheric velocity at the earliest observed phase is 7000 km s−1. As far as we can tell based on the small population of observed LL SNe IIP, the determined values are typical for these events.

Jäger, Jr, Zoltán↗
DOE OSTI2020

Correlations Between Fission Fragment and Neutron Anisotropies in Neutron-Induced Fission

Several sources of angular anisotropy for fission fragments and prompt neutrons have been studied in neutron-induced fission reactions. These include kinematic recoils of the target from the incident neutron beam and the fragments from the emission of the prompt neutrons, preferential directions of the emission of the fission fragments with respect to the beam axis due to the population of particular transition states at the fission barrier, and forward-peaked angular distributions of pre-equilibrium neutrons which are emitted before the formation of a compound nucleus. In addition, there are several potential sources of angular anisotropies that are more difficult to disentangle: the angular distributions of prompt neutrons from fully accelerated fragments or from scission neutrons, and the emission of neutrons from fission fragments that are not fully accelerated. Here, we study the effects of the first group of anisotropy sources, particularly exploring the correlations between the fission-fragment anisotropy and the resulting neutron anisotropy. While kinematic effects were already accounted for in our Hauser-Feshbach Monte Carlo code cgmf, anisotropic angular distributions for the fission fragments and pre-equilibrium neutrons resulting from neutron-induced fission on 233,234,235,238 U, 239,241 Pu, and 237 Np have been introduced for the first time. The effects of these sources of anisotropy are examined over a range of incident neutron energies, from thermal to 20 MeV, and compared with experimental data from the Chi-Nu liquid scintillator array. The anisotropy of the fission fragments is reflected in the anisotropy of the prompt neutrons, especially as the outgoing energy of the prompt neutrons increases, allowing for an extraction of the fission-fragment anisotropy to be made from a measurement of the neutrons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

JETSCAPE framework: p+p results

The JETSCAPE framework is a modular and versatile Monte Carlo software package for the simulation of high energy nuclear collisions. In this work we present a new tune of JETSCAPE, called PP19, and validate it by comparison to jet-based measurements in p + p collisions, including inclusive single jet cross sections, jet shape observables, fragmentation functions, charged hadron cross sections, and dijet mass cross sections. These observables in p+p collisions provide the baseline for their counterparts in nuclear collisions. Quantifying the level of agreement of JETSCAPE results with p+p data is thus necessary for meaningful applications of JETSCAPE to A+A collisions. The calculations use the JETSCAPE PP19 tune, defined in this paper, based on version 1.0 of the JETSCAPE framework. For the observables discussed in this work calculations using JETSCAPE PP19 agree with data over a wide range of collision energies at a level comparable to standard Monte Carlo codes. These results demonstrate the physics capabilities of the JETSCAPE framework and provide benchmarks for JETSCAPE users.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Design and validation of a fission-product source for radioactive beams

A neutron-generator-based source of fission products (nuCARIBU) is being developed to replace the existing source at Argonne National Laboratory. Neutrons generated by reaction are thermalized in a compact polyethylene [(moderator and impinge on a thin 235 U foil to produce fission products. Experiments have been conducted to verify the nuCARIBU concept and simulation results using a 4 MeV, 0.16 µ⁢A proton beam from the ATLAS accelerator irradiating a LiF neutron production target and two different simplified versions of the nuCARIBU (C 2 ⁢H 4 ) n moderator with thicknesses of 3.81 cm and 1.91 cm separating it from the enriched 235 U fission foil. Here, the emitted γ-ray peaks of 1383.93 keV ( 92 Sr), 358.0 keV ( 104 Tc), 565.992 keV ( 134 Te), and 258.411 keV ( 138 Xe) obtained from fission products are used to determine the neutron-induced fission rates in the fission foil. The measured fission rates are in good agreement with those simulated using the PHITS Monte Carlo code with user-defined neutron production cross sections based on well-known cross section data. These experimental results, consistent with the simulations, show that the current 252 Cf source (CARIBU) can be replaced by the neutron induced fission source (nuCARIBU).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Simulations of multicomponent relativistic thermalization

Multicomponent relativistic fluids have been studied for decades. However, simulating the dynamics of the particles and fluids in such a mixture has been a challenge due to the fact that such simulations are computationally expensive in three spatial dimensions. Here, we report on the development and application of a multidimensional relativistic Monte Carlo code to explore the thermalization process in a relativistic multicomponent environment in a computationally inexpensive way. As an illustration we simulate the fully relativistic three-dimensional Brownian-motion-like solution to the thermalization of a high-mass particle (proton) in a bath of relativistic low-mass particles (electrons). Here, we follow the thermalization and ultimate equilibrium distribution of the Brownian-like particle as can happen in the cosmic plasma during big-bang nucleosynthesis. We also simulate the thermalization of energetic particles injected into the plasma as can occur, for example, by the decay of massive unstable particles during the big bang.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Conceptual Design of a Scintillator-Based Fast-Ion Loss Detector for the Wendelstein 7-X Stellarator

Here, a conceptual design of a scintillator-based fast-ion loss detector (FILD) has been developed for the Wendelstein 7-X stellarator (W7-X). Simulations using the Monte Carlo codes ASCOT5 and BEAMS3D have been performed to calculate the expected flux of neutral beam injection (NBI)-generated fast hydrogen ions onto the conceptual detector probe head. These fast-ion loss fluxes have been calculated for several magnetic field configurations as well as probe insertion positions. At the maximum insertion position, both co-and counter-going losses with high incident pitch angles are observed; however, at retracted positions, only co-going fast ions reach the probe head. The FILDSIM code has been used to optimize the geometry of the detector entrance and collimating elements to achieve a wide velocity space coverage as well as a high velocity–space resolution. A synthetic FILD signal is calculated for the expected loss distribution via forward modeling using the instrument response function. The synthetic signal is found to largely retain the velocity space features of the loss distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Ion beam sputtering of silicon: Energy distributions of sputtered and scattered ions

The properties of sputtered and scattered ions are studied for ion beam sputtering of Si by bombardment with noble gas ions. The energy distributions in dependence on ion beam parameters (ion energy: 0.5-1 keV; ion species: Ne, Ar, Xe) and geometrical parameters (ion incidence angle, polar emission angle, and scattering angle) are measured by means of energy-selective mass spectrometry. The presence of anisotropic effects due to direct sputtering and scattering is discussed and correlated with process parameters. The experimental results are compared to calculations based on a simple elastic binary collision model and to simulations using the Monte-Carlo code SDTRIMSP . The influence of the contribution of implanted primary ions on energy distributions of sputtered and scattered particles is studied in simulations. It is found that a 10% variation of the target composition leads to detectable but small differences in the energy distributions of scattered ions. Comparison with previously reported data for other ion/target configurations confirms the presence of similar trends and anisotropic effects: The number of high-energy sputtered ions increases with increasing energy of incident ions and decreasing scattering angle. The effect of the ion/target mass ratio is additionally investigated. Small differences are observed with the change of the primary ion species: The closer the mass ratio to unity, the higher the average energy of sputtered ions. Finally, the presence of peaks, assigned to different mechanisms of direct scattering, strongly depends on the ion/target mass ratio.

42 ENGINEERING↗
DOE OSTI2024

Radiative corrections: from medium to high energy experiments

Abstract Radiative corrections are crucial for modern high-precision physics experiments, and are an area of active research in the experimental and theoretical community. Here we provide an overview of the state of the field of radiative corrections with a focus on several topics: lepton–proton scattering, QED corrections in deep-inelastic scattering, and in radiative light-hadron decays. Particular emphasis is placed on the two-photon exchange, believed to be responsible for the proton form-factor discrepancy, and associated Monte-Carlo codes. We encourage the community to continue developing theoretical techniques to treat radiative corrections, and perform experimental tests of these corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Challenges and opportunities in verification and validation of low temperature plasma simulations and experiments

This paper describes the verification and validation (V&V) framework developed for the stochastic Particle-in-Cell, Direct Simulation Monte Carlo code Aleph. An ideal framework for V&V from the viewpoint of the authors is described where a physics problem is defined, and relevant physics models and parameters to the defined problem are assessed and captured in a Phenomena Identification and Ranking Table (PIRT). Furthermore, numerous V&V examples guided by the PIRT for a simple gas discharge are shown to demonstrate the V&V process applied to a real-world simulation tool with the overall goal to demonstrably increase the confidence in the results for the simulation tool and its predictive capability. Although many examples are provided here to demonstrate elements of the framework, the primary goal of this work is to introduce this framework and not to provide a fully complete implementation, which would be a much longer document. Comparisons and contrasts are made to more usual approaches to V&V, and techniques new to the low-temperature plasma community are introduced. Specific challenges relating to the sufficiency of available data (e.g., cross sections), the limits of ad hoc validation approaches, the additional difficulty of utilizing a stochastic simulation tool, and the extreme cost of formal validation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Developments on two lithium vapor-box linear test-stand experiments

The lithium vapor-box divertor is a possible fusion power exhaust solution.It uses condensation pumping to create a gradient of vapor density in a divertor slot; this should allow a stable detachment front without active feedback.As initial explorations of the concept, two test stands which take the form of three connected cylindrical stainless steel boxes are being developed: one without plasma at PPPL, to test models of lithium evaporation and flow; and one for the linear plasma device Magnum-PSI (at DIFFER in Eindhoven, The Netherlands) to test the ability of a lithium vapor cloud to induce volumetric detachment and redistribute the plasma power.The first experiment uses boxes with diameters of 6 cm, joined by apertures with diameters of 2.2 cm. Up to 1 g of Li is placed in one box, which is heated to up to 600 degrees C. The Li evaporates, then flows to and condenses in the two other, cooler boxes over several minutes. The quantity of Li transported is assessed by weighing the boxes before and after the heating cycle, and is compared to the quantity predicted to flow for the box at its measured temperature using a Direct Simulation Monte Carlo code, SPARTA. With good experimental conditions, the two values agree to within 15%.The experiment on Magnum-PSI is in the conceptual design stage.The design is assessed by simulations using the code B2.5-Eunomia.They show that when the hydrogen-ion plasma beam, with n_e = 4e20 per cubic meter, T_e = 1.5 eV, and r = 1 cm, is passed through a 16 cm long, 12 Pa, 625 degree C Li vapor cloud, the plasma heat flux and pressure on the target are significantly reduced compared to the case without Li.With the Li present, the plasma is cooled by excitation of Li neutrals followed by radiation until it volumetrically recombines, lowering the heat flux from 3.7 MW/m^2 to 0.13 MW/m^2, and the pressure is reduced by 93%, largely by collisions of hydrogen ions with neutral Li.

divertor↗