Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cross Section generation”

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 91 records · Page 5

Status of PyGriffin Development for Integration into the NEAMS Workbench

The integration of Griffin into Workbench, by way of the PyGriffin code package, was initiated by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY-2022. PyGriffin was developed as a Python wrapper to Griffin to streamline the complex workflow involving mesh, cross section, and Griffin input generation, code execution, and results postprocessing, by improving user experience with high fidelity neutronics analysis through the Workbench GUI interface. PyGriffin can be used as a standalone application or through the PyARC code suite. FY-2023 saw the continued development of PyGriffin to expand its capabilities including advanced post-processing of simulation results and initial development of the Monte Carlo (MC) cross section generation workflow with the Shift MC code. PyGriffin was also approved for open-source software status in FY-2023, initiating plans to move the code repository into an open location. In addition to PyGriffin development, there were many improvements and developments made in PyARC leading to several software releases (latest version is v2.3.0), but the focus of this report is PyGriffin.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Advancements in Multiphysics Microdepletion Analysis of an eVinci TM -like Microreactor Leveraging OpenMC-CRAB Workflow

Nuclear microreactors (MRs) are a class of nuclear reactor technology, characterized by reduced dimensions, modular design, and reduced power output in contrast to conventional Light Water Reactors (LWRs). MRs are proposed for supplying electricity and eventual process heat to remote locations, such as military installations and disaster-affected areas. Current research work sponsored by the US Department of Energy Microreactor Program (MRP) is devoted to the development of novel modeling and simulation tools to better support MR vendors and regulatory bodies. Notably, the NRC is projected to utilize the CRAB multiphysics software driver for executing both design and beyond-design-basis accident analyses. Furthermore, the NRC has been utilizing the MELCOR code to calculate mechanistic source terms during accidents. Since MELCOR relies on isotopic inventory and reactor temperature/power profiles under accident conditions, which theoretically can be derived from CRAB, the goal is to establish a comprehensive CRAB-MELCOR computational framework. Past work was focused on testing and demonstrating CRAB's capability to generate results that can be used to inform mechanistic source term calculations in MELCOR. In particular, a computational workflow leveraging OpenMC-generated microscopic cross sections and CRAB was first applied to perform multiphysics microscopic depletion calculation followed by an accident scenario for a stylized microreactor problem. In fiscal year 2024, the research work has been focused on applying the OpenMC-CRAB workflow, which was first tested in fiscal year 2023, to a realistic 3D heat-pipe cooled MR problem representative of the eVinci TM design. The latter computational problem was developed with inputs from WEC to conserve selected neutronic and thermal characteristics of the eVinci TM design without releasing proprietary data. The results of this simulation, encompassing isotopic inventory, power density distribution, and kinetic parameters, will inform both MELCOR and the WEC-developed FATE code for mechanistic source terms calculations. The results from the two codes will then be compared for code verification purposes. This report contains the design characteristics of the realist heat pipe cooled microreactor developed as a use-case for the verification exercise, and the current results for the multiphysics microscopic depletion performed with the OpenMC-CRAB workflow. The results include eigenvalue as a function of time, power distribution at EOL, in addition to nuclides inventory's time evolution and spatial distribution. Finally, we report improvements to the workflow efficiency achieved through a collaboration with the NEAMS programs. Through this collaborative effort, we were able to strongly decrease the computational time for the multiphysics microdepletion calculation (i.e., from 17.4 hours to 5.7 hours on 280 processors) in addition to simplifying the interface to generate isotopics spatial distribution utilizable by FATE and MELCOR. Future work, including the improvement of the current microscopic cross-sections' library and the simulation of an accident scenario at EOL, is also discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

HTGR Simulation Methods & International Collaborations

ART-GCR “Methods” activity is split between Experimental Validation data from the ANL NSTF and OSU HTTF (next three presentations). HTGR core simulation (this presentation). International collaboration within OECD Generation-IV (Gen-IV) and USA/Japan bi-lateral agreements (this presentation) HTGR Simulation Methods No new NE-52 funding for HTGR Methods support in FY20; ~$200K FY19 carry-over funds only. Consists of international code-to-code benchmarks (IAEA CRP on HTGR UAM and OECD/NEA MHTGR-350) and refinement of a few-group Pebble Bed Reactor (PBR) cross section (XS) generation methodology. Funding will be requested in FY21 to produce the final reports for the two benchmarks and continue the development of the PBR XS generation methodology. Additional (non-ART) HTGR-related support work at INL NEAMS: HTR-Application work package at INL Create a benchmark for the pebble shuffling and depletion algorithms being developed for NEAMS Griffin code. iFOA award with X-Energy: Develop independent Monte Carlo model of Xe-100 design. Independent design confirmatory analysis of Xe-100 design using NEAMS tools Griffin and Pronghorn. Support X-Energy design team to use their own legacy design tools (VSOP99 and MGT). Support NEAMS Griffin and Pronghorn development team for the iFOA needs (received $50K additional funding for required development).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Recurrent and convolutional neural networks for sequential multispectral optoacoustic tomography ( MSOT ) imaging

Abstract Multispectral optoacoustic tomography (MSOT) is a beneficial technique for diagnosing and analyzing biological samples since it provides meticulous details in anatomy and physiology. However, acquiring high through‐plane resolution volumetric MSOT is time‐consuming. Here, we propose a deep learning model based on hybrid recurrent and convolutional neural networks to generate sequential cross‐sectional images for an MSOT system. This system provides three modalities (MSOT, ultrasound, and optoacoustic imaging of a specific exogenous contrast agent) in a single scan. This study used ICG‐conjugated nanoworms particles (NWs‐ICG) as the contrast agent. Instead of acquiring seven images with a step size of 0.1 mm, we can receive two images with a step size of 0.6 mm as input for the proposed deep learning model. The deep learning model can generate five other images with a step size of 0.1 mm between these two input images meaning we can reduce acquisition time by approximately 71%.

Juhong, Aniwat↗

Preservation of kinetics parameters generated by Monte Carlo calculations in two-step deterministic calculations

The generation of accurate kinetic parameters such as mean generation time Λ and effective delayed neutron fraction β eff via Monte Carlo codes is established. Employing these in downstream deterministic codes warrants another step to ensure no additional error is introduced by the low-order transport operator when computing forward and adjoint fluxes for bilinear weighting of these parameters. Another complexity stems from applying superhomogenization (SPH) equivalence in non-fundamental mode approximations, where reference and low-order calculations rely on a 3D full core model. In these cases, SPH factors can optionally be computed for only part of the geometry while preserving reaction rates and K-effective, but the impact of such approximations on kinetics parameters has not been thoroughly studied. This paper aims at studying the preservation of bilinearly-weighted quantities in the Serpent–Griffin calculation procedure. Diffusion and transport evaluations of IPEN/MB-01, Godiva, and Flattop were carried out with the Griffin reactor physics code, testing available modeling options using Serpent-generated multigroup cross sections and equivalence data. Verifying Griffin against Serpent indicates sensitivities to multigroup energy grid selection and regional application of SPH equivalence, introducing significant errors; these were demonstrated to be reduced through the use of a transport method together with a finer energy grid.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Ion Mobility Spectrometry-Mass Spectrometry for High-Throughput Analysis

Ion mobility spectrometry is a widely used analytical technique providing gas phase separation of molecules. It has received increasing attention in the recent years with the advancement in technology development and the availability of commercial instruments. In this chapter, we introduced the ion mobility fundamental theory and provided examples of IMS applications, especially for isomer separation, collision cross section database generation, high throughput analysis workflows, software tools for IMS data analysis, and ongoing high resolution SLIM IMS development. While IMS is not yet routinely utilized in drug discovery and pharmaceutical industry, there has been increased interest in high throughput library screening and antibody characterization. With all the ongoing development in IMS technology and informatics, we foresee more and more exciting applications of high throughput IMS analysis in different fields including omics studies, drug discovery and clinical applications in the near future.

Ross, Dylan H.↗

Generation and optimization of cross-sections for electron-C 4 F 7 N collisions

Abstract A complete and consistent set of electron-neutral collision cross-sections for the novel insulating gas C 4 F 7 N is reported. The set is composed of a combination of cross-sections previously reported in literature, optimized via a genetic algorithm in conjunction with a multi-term Boltzmann equation solver, and calculatedab initiousing the R-matrix code Quantemol-EC. The finalized set accurately reproduces reported macroscopic rate and transport coefficients as well as Townsend coefficients and critical electric field strengths in C 4 F 7 N and its mixtures with nitrogen, carbon dioxide, and argon.

Physics↗

SEFOR Core I-E

This repository contains the VTB model developed for the Southwest Experimental Fast Oxide Reactor (SEFOR) core configuration I-E using the NEAMS tools. The MOOSE Reactor Module was employed to set up the geometry mesh, MC2-3 was applied to generate multigroup cross sections at various temperatures, and Griffin was used to calculate the k-eff values of Core I-E across these temperature conditions. To validate the deterministic results from Griffin, reference solutions were obtained with the Monte Carlo code Shift. The reactivity feedback was then derived and compared against SEFOR experimental data from isothermal tests in order to assess the performance of Griffin. Reproducing these results requires licensed access to several code systems, including Shift, Griffin, and MC2-3. Instructions for obtaining and accessing these codes can be found at the following locations: https://www.ornl.gov/scale/releases for SCALE/Shift, https://mooseframework.inl.gov/ncrc/applications/ncrc_root_griffin.html for Griffin, and https://www.anl.gov/nse/mc23 for MC2-3.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Nuclear Data: What is the big deal? [Slides]

This presentation begins with a look at the practical application: differential and integral data. It also includes a nuclear data evaluation: cross section formalisms, data processing and cross section library generation, computer tools, and criticality safety assessment: nuclear data. The presentation ends with some concluding remarks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A COUPLED DETERMINISTIC TRANSPORT CALCULATION IN COMSOL USING PROPER ORTHOGONAL DECOMPOSITION

A coupled physics problem consisting of radiation transport and heat transfer was modeled using COMSOL Multiphysics to simulate a criticality accident. Reduced order models (ROM) were created to accelerate the radiation transport calculation through Proper Orthogonal Decomposition (POD). Additionally, on-the-fly neutron cross sections were generated from a non-linear function through a feed forward neural network. The Livelink for MATLAB module allowed for data transfer between an otherwise independent solver and in-house modules. A data driven approach to POD was developed to act as a surrogate model to couple with heat transfer, with the POD module written in MATLAB coupled to an isolated heat transfer solve in COMSOL. A coupled transient solution resulted in a good approximation of the flux and temperature, with the maximum error being 3E-2 and 2.5E-3 respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Neutronics Modeling of the pulsed plasma rocket reactor using rattlesnake

In the pulsed plasma rocket (PPR) reactor, neutron induced fission processes are utilized to implement a series of pulsed micro-explosions of very high power and eject plasma as a propellant. More specifically, projectiles (bullets) composed of moderated uranium are sent through the chamber of an unmoderated uranium barrel. By inducing rapid fission within the bullets, a plasma can be generated with appropriate delivery of neutrons. This study aims to assess the neutronics performance of the PPR reactor during normal operations using the MOOSE-based Rattlesnake code through the evaluation of the impact of the movement of the fuel bullet and the rotation of control drums on the criticality of the system. The Monte Carlo (MC) code Serpent 2 was employed to generate material-based cross sections for use in Rattlesnake and the reference neutronics solution. Cubit was used to generate the mesh for the Rattlesnake model. Parametric studies were conducted to evaluate the best approaches for cross section and mesh generation to ensure accurate results from Rattlesnake. As part of the verification process, the eigenvalue results of the system at various fuel bullet positions were obtained using Rattlesnake and compared with the reference solutions. The acceptable differences show that the Rattlesnake model with appropriate cross section and mesh generation procedure is a sufficiently accurate approximation of the continuous energy (CE) MC model for micro-sized reactors like the PPR reactor. Next, the impact of the rotation of the control drums on the system criticality was evaluated using the verified model. It was shown that the high-fidelity simulation using the deterministic code Rattlesnake can produce sufficiently accurate results for the evaluation of the reactor's neutronics performance in different phases of the normal operation with acceptable computational cost.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Methodology for Generating Covariance Data of Thermal Neutron Scattering Cross Sections

This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Generation of Enrichment-Dependent Thermal Neutron Scattering Data

This work details the generation of enrichment-dependent thermal neutron scattering cross sections for several crucial uranium fuel compounds. The evaluations of the thermal scattering law (TSL) and associated cross sections for uranium dioxide (UO 2 ), uranium carbide (UC), and uranium nitride (UN) were performed using standard ab initio lattice dynamics (AILD) methods. The data for uranium metal was produced using a novel hybrid approach of molecular dynamics combined with lattice dynamics methods. 235 U enrichments of 5%, 10% (LEU+), 19.75% (HALEU), 93% (HEU), and 100% were considered, in addition to natural uranium. The enrichment-dependent masses and free atom cross sections were used in the generation of elastic and inelastic thermal neutron scattering cross sections, while the calculation of the phonon density of states (DOS) and resulting TSL considered only the natural isotopic composition of uranium. The use of an identical DOS for all enrichments is expected to have minimal impact on the final data, as the small change in uranium mass should not significantly affect lattice vibrations. The cross sections are shown to exhibit significant dependence on 235 U enrichment. The submission of this data to the National Nuclear Data Center (NNDC) for release in the ENDF/B-VIII.1 database should support the design of advanced reactor concepts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

From inclusive to semi-inclusive one-nucleon knockout in neutrino event generators

In neutrino event generators, for models for neutrino and electron scattering only inclusive cross sections are implemented. When these models are used to describe a semi-inclusive cross section, the event generator attaches the hadron variables based on some assumptions. In this work we compare the nucleon kinematics given by the method used in the GENIE event generator, e.g. in the implementation of the SuSAv2 model, to a fully unfactorized calculation using the relativistic distorted wave impulse approximation (RDWIA). We focus on kinematics relevant to the $e4\nu$ analysis and show that observables obtained with RDWIA differ significantly from those of the approximate method used in GENIE, the latter should be considered unrealistic.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Tutorial for Generating Correlated Random Samples in the Context of Replica Cross Section Data Used in the Propagation of Uncertainty (Second Edition)

The following sample problem write-ups are designed as a tutorial for generating correlated random samples (e.g., replica multi-group cross section data) for use in propagation of uncertainty problems. These problems were set up and solved in MATLAB, but any programming environment with basic statistical functions can be used to generate similar results. Since linear sample problems were chosen, helpful comparisons to the “sandwich” rule propagation of uncertainty are available and were used.

97 MATHEMATICS AND COMPUTING↗

Inclusive and differential cross section measurements of $ \textrm{t}\overline{\textrm{t}}\textrm{b}\overline{\textrm{b}} $ production in the lepton+jets channel at $ \sqrt{s} $ = 13 TeV

Measurements of inclusive and normalized differential cross sections of the associated production of top quark-antiquark and bottom quark-antiquark pairs, $ \textrm{t}\overline{\textrm{t}}\textrm{b}\overline{\textrm{b}} $, are presented. The results are based on data from proton-proton collisions collected by the CMS detector at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{−1}$. The cross sections are measured in the lepton+jets decay channel of the top quark pair, using events containing exactly one isolated electron or muon and at least five jets. Measurements are made in four fiducial phase space regions, targeting different aspects of the $ \textrm{t}\overline{\textrm{t}}\textrm{b}\overline{\textrm{b}} $ process. Distributions are unfolded to the particle level through maximum likelihood fits, and compared with predictions from several event generators. The inclusive cross section measurements of this process in the fiducial phase space regions are the most precise to date. In most cases, the measured inclusive cross sections exceed the predictions with the chosen generator settings. The only exception is when using a particular choice of dynamic renormalization scale, $ {\mu}_{\textrm{R}}=\frac{1}{2}{\prod}_{i=\textrm{t},\overline{\textrm{t}},\textrm{b},\overline{\textrm{b}}}{m}_{\textrm{T},i}^{1/4} $, where $ {m}_{\textrm{T},i}^2={m}_i^2+{p}_{\textrm{T},i}^2 $ are the transverse masses of top and bottom quarks. The differential cross sections show varying degrees of compatibility with the theoretical predictions, and none of the tested generators with the chosen settings simultaneously describe all the measured distributions.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Validation of the single-event method for low-energy electron transport via stopping power calculations with $\mathrm{MCNP}$

Monte Carlo simulations of low-energy ( <50 keV) electron transport in matter are essential for a broad range of application fields. Several Monte Carlo codes have developed specialized treatments for this case, but a comprehensive validation of low-energy electron transport for general-purpose simulations remains lacking in the literature. One approach to accomplish this validation is calculation of stopping power using low-energy electron transport physics, as stopping power is a fundamental radiation transport quantity which must be simulated accurately for nearly any application. Here in this work, we use the Monte Carlo N-Particle (MCNP) radiation transport code with the single-event method for electron transport to calculate stopping powers of low-energy electrons (50 eV to 30 keV) in 41 elemental solids, 14 compound solids, and five rare gas solids, comparing simulation results to published semi-empirical stopping power calculations from optical measurements. In general, the simulations give good agreement (typically within ±10%) with semi-empirical stopping power values at higher energies: 300 eV and above for most elemental solids, 1 keV and above for compound solids, and 400 eV and above for rare gas solids. Agreement between MCNP and semi-empirical values is worse below these energies. The most significant source of error is the EPRDATA14 cross section data, which does not account for changes in electronic structure due to solid-state bonding, particularly in compound materials. The simplistic model of atomic excitation used to generate the EPRDATA14 cross sections is another key source of error. Additionally, the breakdown of the continuous slowing-down approximation introduces significant uncertainty at low energies, although this is a limitation of the calculation method and not of the simulation procedure. Accounting for these and other uncertainty sources, the single event method in MCNP is robust and able to give good accuracy for a variety of low-energy electron transport problems through diverse kinds of materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗