Verification and Validation of Nuclear Data Processing by NJOY Using ACEtk and MCNP [Slides]
Abstract not provided.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
SEARCH · Engineering Papers
Results for “NJOY”
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.
Processing files with NJOY - Reading and manipulating ACE files with ACEtk
Explore the source record for details and available documents.
NJOY, MCNP, ACE, ACEtk↗
Overview of NJOY work for NCSP FY22
Explore the source record for details and available documents.
NJOY, NCSP↗
Investigation of the impact of difference between FRENDY and NJOY2016 on neutronics calculations
A nuclear data library is used as a starting input for all subsequent neutronics calculations. NJOY has been used worldwide as a nuclear data processing code to create cross section libraries for a long time. For the verification of NJOY method and for providing an alternative nuclear data processing tool, JAEA has been developed the new nuclear library processing code FRENDY. In this paper, nuclear calculations were performed using the ACE files and the multigroup libraries created by both FRENDY and NJOY, and the impacts on the neutronics characteristics due to nuclear data processing were investigated using those libraries. MCNP was used to compare the ACE files by calculating many benchmark problems including ICSBEP and it was confirmed that the k-eff values generally agree with each other within the range of statistical errors. The multigroup cross sections are verified by the BWR design codes LANCR/AETNA through calculation of a commercial BWR5 equilibrium core loaded with 9*9 fuels. It was confirmed that fuel assembly and core characteristics are consistent with each other. From the above investigations, it was confirmed that FRENDY can provide comparable continuous/multi-group neutron cross sections with NJOY. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
ENDFtk: A robust tool for reading and writing ENDF-formatted nuclear data
ENDFtk is a recently developed C++ and Python interface to interact with ENDF-6 formatted nuclear data files. It provides a robust and complete interface, allowing the reading and writing of all formats currently part of the ENDF-6 formats manual, as well as some non-ENDF formats used by the NJOY processing code. It provides an interface that mimics the names in the ENDF-6 formats manual as well as an equivalent interface using human-readable attribute names. It is robust and powerful enogh for nuclear data experts to develop complex applications, while also simple enough to be used non-experts to retrieve and manipulate evaluated nuclear data. ENDFtk offers the ability to easily interrogate and manipulate data either in large-scale code projects or in simple Python scripts. Here, in this paper, a brief overview of the interface is given, as well as more substantial examples demonstrating plotting simple data, interacting with more complex data, and writing new data to files. ENDFtk is open source and available for download via GitHub (https://github.com/njoy/ENDFtk).
97 MATHEMATICS AND COMPUTING↗
Advanced Modeling and Simulation Methods for Evaluation of Thermal Neutron Scattering Materials
With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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↗
Extending deterministic transport capabilities for very-high and ultra-high energy electron beams
Focused Very-High Energy Electron (VHEE, 50–300 MeV) and Ultra-High Energy Electron (UHEE, > 300 MeV) beams can accurately target both large and deeply seated human tumors with high sparing properties, while avoiding the spatial requirements and cost of proton and heavy ion facilities. Advanced testing phases are underway at the CLEAR facilities at CERN (Switzerland), NLCTA at Stanford (USA), and SPARC at INFN (Italy), aiming to accelerate the transition to clinical application. Currently, Monte Carlo (MC) transport is the sole paradigm supporting preclinical trials and imminent clinical deployment. In this paper, we propose an alternative: the first extension of the nuclear-reactor deterministic chain Njoy-Dragon for VHEE and UHEE applications. We have extended the Boltzmann-Fokker-Planck (BFP) multigroup formalism and validated it using standard radio-oncology benchmarks, complex assemblies with a wide range of atomic numbers, and comprehensive irradiation of the entire periodic table. We report that 99% of water voxels exhibit a BFP-MC deviation below 2% for electron energies under 1.5 GeV. Additionally, we demonstrate that at least 97% of voxels of bone, lung, adipose tissue, muscle, soft tissue, tumor, steel, and aluminum meet the same criterion between 50 MeV and 1.5 GeV. For water, the thorax, and the breast intra-operative benchmark, typical average BFP-MC deviations of 0.3% and 0.4% were observed at 300 MeV and 1 GeV, respectively. By irradiating the entire periodic table, we observed similar performance between lithium (Z = 3) and cerium (Z = 58). Deficiencies observed between praseodymium (Z = 59) and einsteinium (Z = 99) have been reported, analyzed, and quantified, offering critical insights for the ongoing development of the Evaluated Nuclear Data File mode in NJOY.
62 RADIOLOGY AND NUCLEAR MEDICINE↗
IRMA
IRMA (In)elastic Representation of Materials As S(α,β) evaluations IRMA turns one phonon model into three outputs that usually require three separate tool chains: an evaluated nuclear-data file, predicted neutron-scattering spectra, and scattering kernels for Monte Carlo transport. The three outputs draw on a single, consistent description of the material, so the evaluation, the spectroscopy that can validate it, and the transport that uses it always agree about the physics. Nuclear data. IRMA writes ENDF-6 File 7 thermal scattering evaluations on automatically constructed (α, β) grids. This part reimplements and generalizes NJOY's LEAPR: the classic kernels reproduce freshly generated NJOY2016 tapes digit for digit and published reference tapes to about 1e-4, and the generalized paths add the exact coherent one-phonon term, anisotropic Debye-Waller tensors, coherent elastic for arbitrary crystals, and a per-species partition for polyatomic materials. The tapes feed NJOY, AMPX, FUDGE, and every transport code downstream of them. Neutron spectroscopy. The irma.spectra forward model projects the same physics onto an instrument's kinematics and resolution: INS spectra for VISION and generic indirect geometries, and 2-D S(Q,E) powder maps for direct-geometry spectrometers, from a phonopy model or straight from a phonon DOS. It can be used to predict a proposed measurement before beam time; in analysis, it supplies the calculated single-scattering counterpart of a measured spectrum, from the same material description the evaluation was built from. Monte Carlo transport. The irma.ncrystal exporter writes per-temperature scattering kernels for the companion NCrystal plugin, so McStas, OpenMC, and other NCrystal-aware codes sample the same physics. The exported kernels carry the per-site anisotropic Debye-Waller tensors, keeping directional coherent-elastic physics that NCrystal's standard scalar treatment does not represent. With the same physics inside a transport code, an entire beamline becomes a virtual experiment: IRMA's end-to-end validation ran a custom McStas implementation of the ARCS spectrometer, assembled from the existing McVine and McStas models, against measured data. From a bare crystal structure. The irma mlip front end builds the phonon model itself: a structure file and a choice of potential are enough. Nine pretrained machine-learned interatomic potentials are supported, on a laptop CPU, with no first-principles calculation; an approximate phonon model for a new material costs minutes, not a DFT campaign, and the build emits prefilled inputs for all three outputs. The result is a good starting point rather than a finished evaluation: survey-quality physics with every parameter exposed for review. A converged atomistic calculation enters the same way, as a phonopy model, when higher fidelity is needed.
Ramic, Kemal [Oak Ridge National Laboratory (ORNL)↗
(U) Average Scattering Angle Cosine Sensitivities
The Evaluated Nuclear Data File (ENDF) provides covariances of the average scattering angle cosine $\overline{μ}$ for 117 nuclides, and the NJOY nuclear data processing code can process these covariances. To propagate them to the uncertainty of a response due to uncertainties in $\overline{μ}$ requires sensitivities of the response with respect to $\overline{μ}$. This report uses NJOY’s definition of $\overline{μ}$ and derives a relationship between the sensitivity of a response to $\overline{μ}$ and the sensitivity of the response to the anisotropic P 1 elastic scattering cross section. This relationship has been implemented in the multigroup neutron sensitivity code SENSMG, which relies on the discrete ordinates code PARTISN for the neutron transport calculations. Kodeli has also calculated sensitivities with respect to $\overline{μ}$ and associated uncertainties. However, the formal derivation of this paper shows that the definition of the sensitivity with respect to $\overline{μ}$ is ambiguous because there are multiple ways to change $\overline{μ}$ by changing scattering cross sections. Here we derive two such ways. Both are implemented in SENSMG. This report is organized as follows. Section II discusses $\overline{μ}$ and its significance. Section III presents the equations used for the sensitivity of a response with respect to $\overline{μ}$, and Section IV discusses how to compute these sensitivities with SENSMG and MCNP6.3. Section V is a fast-spectrum test problem, and Sec. VI is a thermal-spectrum problem. Section VII is a summary. The example problem inputs are listed in the appendix.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
A Primer on Nuclear "Recoil" Data
For many years the ENDF-6 format has existed to contain evaluated nuclear data. One primary driver for the format and the data is neutron transport calculations. Evaluated data can be processed by a code like NJOY into either continuous energy form (ACE format) for Monte Carlo codes like MCNP or into multi-group form (NDI tables) for deterministic codes like Partisn. The reaction cross sections are found in the MF 3 section of the ENDF-6 format. If the reaction produces one or more neutrons as outputs, then secondary neutron data must also be given in MF 4,5 or 6 sections of the format. Energy and angular distributions for the output neutrons must be given in one of several available formats. The most general formats are found in MF 6. MF 4 is for angular distributions, MF 5 is for Energy distributions, and MF 6 contains both. In recent years, more interest has developed in the other output particles (i.e., the “recoil” particles) from neutron induced reactions. This has been driven by interest in charged particle transport and in more specialized partial kermas. ( A separate total neutron kerma has been available for a long time.) Partial kermas are a breakdown of the Kinetic Energy Released into the MAterial by output particle or by neutron reaction. The sum of the partial kermas should be equal to the total kerma on a group-wise basis. The ENDF-6 format is general enough for these new data requirements. Ideally, evaluated data would exist for every output particle (including the secondary neutrons) from every neutron-induced reaction. In the MF 6 format section, data for multi-particle outputs may be entered using the LAW =1 option. This option explicitly allows energy and angular output distributions for each particle produced in the reaction. Such information should preserve the balance between partial kermas and the total kerma at the groupwise level as well as the individual particle averaged energies. For simpler 2-body reactions, LAW = 2 is available. This allows the evaluation to specify only enough data to specify the 2-body reaction fully. Such a simplification does not exist for 3 (or more) body breakups. A simplified form of LAW 2, i.e., LAW 3, also exists to generate approximate recoil output distributions in the absence of full data. However, evaluated data does not generally exist at this fine granularity for reactions involving more than 2 output particles ( e.g., the 3-body break-up reaction). When the multi-body detailed output distribution data is not available, LAW=6 may be employed in NJOY. LAW 6 produces approximate output distributions for all the particles produced in the reaction. Just like in the 2-body case, the smaller particles will generally carry off more energy. At any particular energy, the balance between the sum of the partial kermas and the total kerma will not necessarily be preserved. However, the average energy of each output particle from a reaction is preserved.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Advanced Modeling and Simulations for Evaluation of Thermal Neutron Scattering Materials [Slides]
In this presentation, researchers detail tradeoffs between two simulation suites for thermal neutron scattering data. OXLIMAX + MCViNE provide a more accurate elastic peak shape and NJOY + MCNP have more accurate spectra between phonon mode peaks (most likely due to multiple neutron scattering). Additionally, some outstanding issues with MCViNE are noted, including multiple neutron scattering and non-uniform (Q,E) grid for data.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Multigroup Cross-section Generation in MCNP6.3 [Slides]
This presentation states that in comparison to the NJOY-produced multigroup cross sections, the MCNP-produced multigroup cross sections are generally consistent. Statistical uncertainties, however, are challenging and the unresolved resonance region may be looked at in the future. It also discusses how the SPM and LCS options were compared to each other for internal consistency. Additionally, some reactor pin-cell-like problems were used to compare to multigroup capabilities in other Monte Carlo codes (e.g., Serpent, OpenMC).
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Evaluation of δ-Phase ZrH1.4 to ZrH1.7 Thermal Neutron Scattering Laws Using Ab Initio Molecular Dynamics Simulations
Zirconium hydride is commonly used for next-generation reactor designs due to its excellent hydrogen retention capacity at temperatures below 1000 K. These types of reactors operate at thermal neutron energies and require accurate representation of thermal scattering laws (TSLs) to optimize moderator performance and evaluate the safety indicators for reactor design. In this work, we present an atomic-scale representation of sub-stoichiometric ZrH2−x(0.3≤x≤0.6), which relies on ab initio molecular dynamics (AIMD) in tandem with velocity auto-correlation (VAC) analysis to generate phonon density of states (DOS) for TSL development. The novel NJOY+NCrystal tool, developed by the European Spallation Source community, was utilized to generate the TSL formulations in the A Compact ENDF (ACE) format for its utility in neutron transport software. First, stoichiometric zirconium hydride cross sections were benchmarked with experiments. Then sub-stoichiometric zirconium hydride TSLs were developed. Significant deviations were observed between the new δ-phase ZrH2−x TSLs and the TSLs in the current ENDF release. It was also observed that varying the hydrogen vacancy defect concentration and sites did not cause as significant a change in the TSLs (e.g., ZrH1.4 vs. ZrH1.7) as was caused by the lattice transformation from ϵ- to δ-phase.
42 ENGINEERING↗
A proposed method for addressing large unphysical uncertainties in mubar
Mubar uncertainties in Section MF 34 MT 2 of ENDF/B-VIII.0 data are sometimes too large. Physical limitations of the bounds of mubar limit the maximum uncertainty to be < 1.0 and usually $\ll$ 1.0. Two ad hoc methods are proposed to artificially constrain the randomly sampled values of mubar to allowable values. The application of NJOY’s mubar covariance matrix to the “sandwich rule” is also discussed and a consistent way to apply sensitivity coefficients is described. Mubar uncertainties for the Jezebel k eff were found to be 156 pcm – a value comparable to other evaluated cross section uncertainties.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
A Study of the Covariance Data in ENDF/B VIII.0 for Low Z Isotopes
Thirty group covariance data have been produced from ENDF/B VIII.0 MF 33 for H1, H2, He4, Li6, Li7, Be9, B10, B11, C12, C13, N15, O16, the JENDL O16 and F19. Multi-group cross section covariance data was produced by NJOY routine ERRORR. Negative eigenvalues were found only for F19, but the covariance constraints required for the condition of the sum of partial cross sections being equal to the resultant total cross section were found to be incorrect for H2, Li6, B10, C12, C13, and O16. An easy correction is suggested for Li6 and C13. In closing, an expansion of the ENDF-6 formats manual is suggested for the MF 33 chapter and a new proposed covariance normalization scheme is proposed.
07 ISOTOPE AND RADIATION SOURCES↗
(U) SENSMG: First-Order Sensitivities of Neutron Reaction Rates, Reaction-Rate Ratios, Leakage, k eff , α , and Subcritical Multiplication Using PARTISN
SENSMG is a tool for computing first-order sensitivities of neutron reaction rates, reaction-rate ratios, leakage, k eff , α, and subcritical multiplication using the PARTISN multigroup discrete-ordinates code. SENSMG computes sensitivities to all of the transport cross sections and data (total, fission, a nu, chi, and all scattering moments), two edit cross sections (absorption and capture), and the density for every nuclide and energy group. It also computes sensitivities to the mass density for every material and derivatives with respect to all interface locations and outer boundaries. It computes sensitivities to user specified reactions whose cross sections are available in a user-supplied NJOY output file. The tool can be used for one-dimensional spherical and slab (r) and two-dimensional cylindrical (r-z) geometries. The tool can be used for fixed-source and eigenvalue problems. For most responses, the tool implements Generalized Perturbation Theory (GPT) as discussed by Williams and Stacey. The tool is thus limited to computing sensitivities only for GPT-allowable responses. For subcritical multiplication, the tool implements sensitivities derived by O’Brien and Clark. SENSMG has a similar role as the old SWANLAKE (Ref. 8), FORSS (Ref. 9), and SENSIT (Ref. 10) codes. It has capabilities similar to those of SUSD3D (Refs. 11 and 12), which also uses PARTISN. Section II of this report describes the theory behind adjoint-based sensitivities, gives the equations that SENSMG solves, and defines the sensitivities that are output. Section III describes the user interface, including the input file and command line options. Section IV describes the output. Section V gives some notes about the coding that may be of interest. Section VI presents some sample problems and discusses verification, which is ongoing. Section VII lists needs and ideas for future work. Appendix A lists most of the input files whose results are presented in Sec. VI. Appendix B provides some useful details on one of the cross-section libraries that SENSMG supports.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Contributions for ENDF/B-VIII.1 Paper
It was discussed at WANDA2020 that some issues were found by the community and users in ENDF/B-VIII.0 covariances after their release. It was mentioned that one issue in testing ENDF/B-VIII.0 covariances was that they were provided by evaluators too late in the β-release phases that insufficient time for thorough testing was available. Therefore, the executive committee requested that first nuclear data covariances should already be delivered for ENDF/B-VIII.1β 1 in February 2023—about half a year earlier than was the case for ENDF/B-VIII.0. This additional time allowed for more stringent testing with LANL’s ENDFtk and CovVal codes, and ORNL’s AMPX code during the release phases. The LLNL covariance testing employed a combination of NJOY, empy tools from EMPIRE, and x4i.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗