Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Models & methods for nuclear reactions”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

NRC Multiphysics Analysis Capability Deployment FY2020 - Part 4

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Sup- port for Advanced Non-Light Water Reactors.” The tasks completed for this report are: • Task 3a: This task demonstrates the effectiveness of Griffin’s methods for interpolating cross sections and SPH correction factors as a function of control drum rotation. The Monte Carlo code Serpent is used to generate reference results. The accuracy of Griffin-calculated eigenvalues and reaction rates are studied as a function of control drum rotation for cross section and SPH factor libraries of varying fidelity. Higher-fidelity libraries require less interpolation and are therefore inherently more accurate; however, they require more computational resources to create. Overall, Griffin is able to accurately model eigenvalue and reaction rates over the full range of drum rotation using a library generated from three to five discrete drum rotations for a particular microreactor design. • Task 3c: This preliminary task successfully demonstrates the performance of a coupled neutron-photon transport calculation with Griffin using the spherical harmonics approximation. The Griffin solutions are compared to reference solutions from the continuous energy Monte Carlo codes MCNP and Serpent. Furthermore, this work also demonstrates how to prepare the neutron and photon libraries using NJOY, MCNP, and Serpent with the latest ENDF/B-VIII data for 235 U and graphite. The Griffin total energy deposition in the active core region is calculated within 0.3% of the reference, but it degrades in the reflector regions with a maximum difference of -1.5%. Overall, the results look promising and future analysis should use an improved data stream and focus on determining the potential cancellation of error within each energy group.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A Predictive Model for Offgas Composition in Waste Glass Melters - Spring 2021 SULI General Audience Abstract

Approximately 56 million gallons of nuclear waste are currently stored in 177 underground storage tanks at the Hanford site in Washington state. A portion of the tank waste is targeted for immobilization by converting the nuclear waste to borosilicate glass through a process called vitrification. During vitrification, emissions are produced in addition to the intended glass product. These emissions, contained in a volume above the molten glass called the plenum, are composed of water vapor, air components, aerosols, particulates, and other gaseous species that result from evaporation, reactions, and air influx. The formation of these gasses causes the surface of the glass to foam. To reduce foaming, sucrose is added to the waste stream before the conversion to glass. Sucrose also reduces the formation of hazardous gasses in the plenum including nitrous oxides, NOx, and carboxides, COx. However, if excess sucrose is present, the reactions that produce the gasses in the plenum will be incomplete, and increased levels of NOx and COx can be observed. The Waste Treatment and Immobilization Plant at the Hanford site is equipped with the technology to treat the plenum emissions. Though many tests have been run to determine the presence of particulates in the plenum emissions, there is does not exist an established method to confidently predict the composition of the gaseous species in the plenum prior to treatment without physical testing. A model was created in the process modeling software Aspen Plus to predict the composition of NOx and COx in the plenum. The model employs a constant stirred tank reactor and user-defined reaction chemistry to determine the plenum composition. To verify the model, the offgas predicted by the model were compared to available offgas data for ten different feed compositions. It is demonstrated that the developed Aspen Plus model is capable of predicting the composition of NOx and COx in the plenum.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Direct and Indirect Measurements for a Better Understanding of the Primordial Nucleosynthesis

The Big Bang Nucleosynthesis (BBN) model is a great success of nuclear astrophysics due to the outstanding agreement between observational and predicted light elements abundances. One exception, however, is the so-called “lithium problem.” In this context, experimental efforts to measure the relevant reactions have been brought to an increased level of accuracy in measuring primordial abundances, and the introduction of indirect methods has allowed researchers to overcome the natural limitations of direct measurements in the energy range of interest for BBN. Here we review the results obtained from the application of the Trojan Horse Method to some of the most influential reactions of the standard network, such as 2 H(d,p) 3 H, 2 H(d,n) 3 He, 3 He(d,p) 4 He, 7 Li(p,α) 4 He, and 7 Be(n,α) 4 He. The relevant cross sections have been then used as new inputs to a classical BBN code, resulting in important constraints that make suggestions for a possible solution for the lithium problem outside of nuclear physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Fast emulation of quantum three-body scattering

Here, we develop a class of emulators for solving quantum three-body scattering problems. They are based on combining the variational method for scattering observables and the recently proposed eigenvector continuation concept. The emulators are first trained by the exact scattering solutions of the governing Hamiltonian at a small number of points in its parameter space, and then employed to make interpolations and extrapolations in that space. Through a schematic nuclear-physics model with finite-range two and three-body interactions, we demonstrate the emulators to be extremely accurate and efficient. The computing time for emulation is on the scale of milliseconds (on a laptop), with relative errors ranging from 10 –13 to 10 –4 depending on the case. The emulators also require little memory. We argue that these emulators can be generalized to even more challenging scattering problems. Furthermore, this general strategy may be applicable for building the same type of emulators in other fields, wherever variational methods can be developed for evaluating physical models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Validation of numerical models for seismic fluid-structure-interaction analysis of nuclear, safety-related equipment

Seismic design, qualification, and risk assessment of nuclear safety-related vessels filled with liquid will have to consider the interaction between the vessel (tank), the contained liquid, and submerged components, if any. Seismic fluid-structure-interaction (FSI) analysis of nuclear vessels will rely on numerical models, which are required to be verified and validated. This study validates previously-verified numerical models using test data generated from earthquake-simulator experiments involving a cylindrical tank. Two solvers in LS-DYNA for FSI simulations are used for the numerical analysis: Arbitrary-Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD). Numerical and test results are compared for responses critical to seismic design of advanced reactor vessels: hydrodynamic pressures on the tank wall, reactions at the support, and wave heights of the contained liquid. Analysis is performed for one-, two-, and three-directional seismic inputs with a range of intensity, and rocking motions. The accuracy of the numerical results and different methods for outputting wave heights are discussed. Recommendations for validation of seismic FSI numerical models of advanced reactors are provided. As a result, the validation exercise presented in this study is broadly applicable to cylindrical tanks, regardless of industry sector.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

Deciphering Fingerprints of Stellar Nucleosynthesis Through Nuclear Reaction Rate Measurements and Isotopic Analyses of Stardust

Microscopic stardust grains, found in trace amounts in primitive meteorites, afford a unique opportunity to study stellar nucleosynthesis in the laboratory by measuring their heavy element isotopic compositions. Applying LLNL’s state-of-the art material characterization capabilities to these micrometer-sized grains yields large new datasets with unprecedented precision. However, the astrophysical models necessary to interpret the isotopic record require improved nuclear physics data to constrain stellar nucleosynthesis conditions. For example, the branch point in the s-process path at 95 Zr controls the relative amounts of heavier nuclides but is poorly understood. We conducted a coordinated effort across disciplines to combine a new determination of the neutron capture cross section of 95 Zr with the laboratory analysis of hundreds of stardust grains to yield a view into stellar interiors. We developed new methods to quantify multiple isotopic systems simultaneously (e.g., Zr, Ba, and W) in stardust grains and improved the determination of the neutron capture cross section of unstable 95 Zr via the surrogate reaction method. This interdisciplinary approach is allowing us to fill gaps in our understanding of the environments in which elements are produced while enhancing LLNL’s experimental and theoretical capabilities relevant to nuclear threat reduction missions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Deciphering Fingerprints of Stellar Nucleosynthesis Through Nuclear Reaction Rate Measurements and Isotopic Analyses of Stardust

Microscopic stardust grains, found in trace amounts in primitive meteorites, afford a unique opportunity to study stellar nucleosynthesis in the laboratory by measuring their heavy element isotopic compositions. Applying LLNL’s state-of-the art material characterization capabilities to these micrometer-sized grains yields large new datasets with unprecedented precision. However, the astrophysical models necessary to interpret the isotopic record require improved nuclear physics data to constrain stellar nucleosynthesis conditions. For example, the branch point in the s-process path at 95 Zr controls the relative amounts of heavier nuclides but is poorly understood. We conducted a coordinated effort across disciplines to combine a new determination of the neutron capture cross section of 95 Zr with the laboratory analysis of hundreds of stardust grains to yield a view into stellar interiors. We developed new methods to quantify multiple isotopic systems simultaneously (e.g., Zr, Ba, and W) in stardust grains and improved the determination of the neutron capture cross section of unstable 95 Zr via the surrogate reaction method. This interdisciplinary approach is allowing us to fill gaps in our understanding of the environments in which elements are produced while enhancing LLNL’s experimental and theoretical capabilities relevant to nuclear threat reduction missions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Weighted active space protocol for multireference machine-learned potentials

Multireference methods such as multiconfiguration pair-density functional theory accurately capture electronic correlation in systems with strong multiconfigurational character, but their cost precludes direct use in molecular dynamics. Combining these methods with machine-learned interatomic potentials (MLPs) can extend their reach. However, the sensitivity of multireference calculations to the choice of the active space complicates the consistent evaluation of energies and gradients across structurally diverse nuclear configurations. To overcome this limitation, we introduce the weighted active space protocol (WASP), a systematic approach to assign a consistent active space for a given system across uncorrelated configurations. By integrating WASP with MLPs and enhanced sampling techniques, we propose a data-efficient active learning cycle that enables the training of an MLP on multireference data. We demonstrated the approach on the TiC + -catalyzed C–H activation of methane, a reaction that poses challenges for Kohn–Sham density functional theory due to its significant multireference character. This framework enables accurate and efficient modeling of catalytic dynamics, establishing a paradigm for simulating complex reactive processes beyond the limits of conventional electronic-structure methods.

enhanced sampling↗
NASA NTRS1983

Space radiation interaction mechanisms in materials

Models of charged-particle impact under conditions typical of the space environment are reported, with a focus on impact excitation and nuclear reactions, especially for heavy ions. Impact excitation is studied by using a global model for electronic excitation based on formal relations through the classical dielectric function to derive an approximation related to the local plasma (electron density distribution) within the atoms and molecules; and corrections to the model resulting from the nonfluid nature of this plasma are discussed. Nuclear reactions are studied by reducing quantum-mechanical treatment of this general N-body problem to an equivalent two-body problem that is solvable, and by comparing the results with experimental data. The equations for heavy-charged-particle transport are derived and solution techniques demonstrated. Finally, these methods of analysis are applied to study the change in the electrical properties of a GaAs semiconductor for photovoltaic applications; proton damage to GaAs crystals is found to arise from stable replacement defects and to be nonannealable, in contrast to electron-induced damage.

Wilson, J. W.↗
DOE OSTI2020

Bayesian analysis of the Zn 70 ( d , He 3 ) Cu 69 transfer reaction

Transfer reactions provide information about the single-particle nature of nuclear levels. In particular, the differential cross sections from these measurements are sensitive to the angular momentum of the transferred particle and the spectroscopic factor of the populated level. However, the process of extracting these properties is subject to uncertainties, both from experimental and theoretical sources. By integrating the distorted wave Born approximation into a Bayesian model, we propagate these uncertainties through to the spectroscopic factors and orbital angular momentum values. We use previously reported data of the proton pickup reaction Zn 70 ( d , He 3 ) Cu 69 as an example. By accounting for uncertainties in the experimental data, optical model parameters, and reaction mechanism, we find that the extracted spectroscopic factors for low-lying states of 69 Cu are subject to large, asymmetric uncertainties ranging from 35 to 108%. Additionally, Bayesian model comparison is employed to assign probabilities to each of the allowed angular momentum transfers. This method confirms the assignments for many states, but suggests that the data for a state lying at 3.70 MeV are better characterized by an ℓ =3 transfer, rather than the previously reported ℓ =2 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Measurement of the branching ratio of 16 N , 15 C , 12 B , and 13 B isotopes through the nuclear muon capture reaction in the Super-Kamiokande detector

The Super-Kamiokande detector has measured solar neutrinos for more than 25 years. The sensitivity for solar neutrino measurement is limited by the uncertainties of energy scale and background modeling. Decays of unstable isotopes with relatively long half-lives through nuclear muon capture, such as 16 N, 15 C, 12 B, and 13 B, are detected as background events for solar neutrino observations. Here, in this study, we developed a method to form a pair of stopping muon and decay candidate events and evaluated the production rates of such unstable isotopes. We then measured their branching ratios considering both their production rates and the estimated number of nuclear muon capture processes as Br⁡( 16 N) = (9.0 ± 0.1)%, Br⁡( 15 C) = (0.6 ± 0.1)%, Br⁡( 12 B) = (0.98 ± 0.18)%, Br⁡( 13 B) = (0.14 ± 0.12)%, respectively. The result for 16 N has world-leading precision at present and the results for 15 C, 12 B, and 13 B are the first branching ratio measurements for those isotopes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

SAM: A Modern System Code for Advanced Non-LWR Safety Analysis

The System Analysis Module (SAM), developed at Argonne National Laboratory and by collaborators at other organizations, is for advanced non–light water reactor safety analysis. SAM aims to provide fast-running, modest-fidelity, whole-plant transient analysis capabilities that are essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. To facilitate code development, SAM utilizes the MOOSE object-oriented application framework, its underlying finite element library, and linear and nonlinear solvers to leverage modern advanced software environments and numerical methods. SAM aims to solve tightly coupled physical phenomena, including fission reaction, heat transfer, fluid dynamics, and thermal-mechanical responses in advanced reactor structures, systems, and components with high accuracy and efficiency. Finally, this paper gives an overview of the SAM code development, including goals and functional requirements, physical models, current capabilities, verification and validation, software quality assurance, and examples of simulations for advanced nuclear reactor applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

How analytic choices can affect the extraction of electromagnetic form factors from elastic electron scattering cross section data

Scientists often try to incorporate prior knowledge into their regression algorithms, such as a particular analytic behavior or a known value at a kinematic endpoint. Unfortunately, there is often no unique way to make use of this prior knowledge, and thus, different analytic choices can lead to very different regression results from the same set of data. In this paper, to illustrate this point in the context of the proton electromagnetic form factors, we use the Mainz elastic data with its 1422 cross section points and 31 normalization parameters. Starting with a complex unbound non-linear regression, we will show how the addition of a single theory-motivated constraint removes an oscillation from the magnetic form factor and shifts the extracted proton charge radius. We then repeat both regressions using the same algorithm, but with a rebinned version of the Mainz dataset. These examples illustrate how analytic choices, such as the function that is being used or even the binning of the data, can dramatically affect the results of a complex regression. These results also demonstrate why it is critical when using regression algorithms to have either a physical model in mind or a firm mathematical basis

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Secondary neutron production from thick target deuteron breakup

Thick target deuteron breakup is a variable-energy accelerator-based source of high-energy neutrons, with applications in fundamental and applied nuclear science and engineering. However, the breakup mechanism remains poorly understood, and data on neutron yields from thick target breakup remain relatively scarce. In this work, the double-differential neutron yields from deuteron breakup have been measured on a thick beryllium target at ϵ d = 33 and 40 MeV, using both time-of-flight and activation techniques. We have also introduced a simple hybrid model for the double-differential deuteron breakup cross section, applicable in the ϵ d = 10-100 MeV energy range on light (Z ≤ 6) targets. This model features four empirical parameters that have been fit to reproduce experimental breakup measurements on beryllium targets, using the method of least squares. It was shown that these parameters extrapolate well to higher energies, and to other low-Z target materials. Here, we also include optimization of the parameters that modify the Kalbach systematics for compound and preequilibrium reactions, in order to better reproduce the experimental data for beryllium targets at large angles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

State densities of heavy nuclei in the static-path plus random-phase approximation

We report that nuclear state densities are important inputs to statistical models of compound-nucleus reactions. State densities are often calculated with self-consistent mean-field approximations that do not include important correlations and must be augmented with empirical collective enhancement factors. Here, we benchmark the static-path plus random-phase approximation (SPA + RPA) to the state density in a chain of samarium isotopes 148–155 Sm against exact results (up to statistical errors) obtained with the shell-model Monte Carlo (SMMC) method. The SPA + RPA method incorporates all static fluctuations beyond the mean field together with small-amplitude quantal fluctuations around each static fluctuation. Using a pairing plus quadrupole interaction, we show that the SPA + RPA state densities agree well with the exact SMMC densities for both the even- and odd-mass isotopes. For the even-mass isotopes, we also compare our results with mean-field state densities calculated with the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation. We find that the SPA + RPA repairs the deficiencies of the mean-field approximation associated with broken rotational symmetry in deformed nuclei and with the violation of particle-number conservation in the pairing condensate. In particular, in deformed nuclei the SPA + RPA reproduces the rotational enhancement of the state density relative to the mean-field state density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Development of HELIOS/BIPR/PARCS/MCNP6 Computation Route for WWER RPV Neutron Fluence Analysis and Validation Against Ex-Vessel Detector Measurement Data

Here, a computational route was developed for precise calculation of fast neutron fluence on a WWER-type reactor pressure vessel (RPV). The method is based on the transfer of neutronics data from HELIOS-2 lattice calculations and nodal diffusion neutronics data (power, density, and temperature) from BIPR7.1 and PARCS 3.36/PATHS core calculations into a three-dimensional (pinwise axially distributed) fixed neutron source for modeling of transport of fast neutrons from the reactor core to the outer surface of the RPV using MCNP6.2. Validation of the proposed computational method was carried out based on comparative analysis of MCNP6.2-predicted and neutron dosimetry–measured reaction rates [ 54 Fe(n,p) 54 Мn, 93 Nb(n,n') 93 mNb, and 58 Ni(n,р) 58 Со] on the outer surface of the Armenian Nuclear Power Plant (ANPP) Unit 2 RPV. Validation revealed that the MCNP6.2-predicted fast neutron fluence results are very sensitive to the ENDF-B neutron data. Particularly, MCNP6.2 with ENDF/B-VIII.0 significantly underpredicts (20% to 30%) fast neutron fluence while using ENDF/B-VII.1 data overpredicts it. Adding revised beta-released evaluations of 54 Fe, 56 Fe, 57 Fe, and 16 O from the International Nuclear Data Evaluation Network (INDEN) to ENDF/B-VIII.0 allows one to obtain reasonable agreement with measurement results for all types of measured reaction rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Spectroscopic Signatures of MQ-Resins in Silicone Elastomers

Polysiloxane elastomers have a large application space due to their versatile cross-linking chemistry and highly tunable physical and mechanical properties. One approach for improving the mechanical integrity of commercial polysiloxane “silicone” elastomers while maintaining their optical transparency is the addition of small, silicone-resin molecules to the network. However, both the poly(dimethylsiloxane) (PDMS) network and the silicone-resin particles have an amorphous structure and complex chemistry, which makes the characterization of their structural properties and segmental network dynamics difficult. Here, we report the synthesis and characterization of a series of model silicone networks modified with a specific class of silicone-resin known as MQ-resin using Raman and advanced nuclear magnetic resonance (NMR) spectroscopy methods. Raman spectroscopy was successfully used to quantify the contribution of the MQ-resin to the network, to determine the type of MQ-resin present in the network, and to investigate the completeness of the network cross-linking reaction. Furthermore, solid-state and 1H double-quantum (DQ) NMR spectroscopies were used not only as a detection method for the MQ-resins but also to quantify changes in the segmental dynamics of the network as a function of MQ-resin concentration. The combination of Raman and NMR spectroscopies describes a series of samples where the MQ-resin particles and PDMS chains maintain their independent segmental dynamics up to high concentrations of MQ-resins (40–50% MQ), where the physical properties of the resin dominate the physical properties of the overall network. The results from our spectroscopic analyses are consistent with the results from macroscopic characterization techniques such as solvent uptake and mechanical testing. The spectroscopic insights into the structure–property relationships of PDMS-MQ composites presented in this study are a valuable tool not only for the synthesis and reverse engineering of future generations of commercial silicone elastomers but also for understanding the mechanisms of aging and degradation over the material lifetime.

36 MATERIALS SCIENCE↗