Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Modeling”

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 415 records · Page 23

Model for Independent Particle Motion

Independent particle model in nuclear physics assumes that the nucleon in the nucleus moves in the average (mean field) potential generated by all other nucleons. This chapter gives a short overview of basic features of the independent particle motion in atomic nuclei and its theoretical realization in the framework of shell models for spherical, deformed, and rotating nuclei as well as in more sophisticated approaches such as microscopic+macroscopic model and density functional theories. Independent particle motion of nucleons leads to global and single-particle consequences. The global ones manifest themselves in the shell structure, and its consequences for the global structure of nuclear landscape, the existence of superheavy nuclei, and the superdeformation at high spin are briefly reviewed. Finally, the latter shows itself in the single-particle properties such as energies, alignments, and densities; their manifestations are illustrated on specific examples.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Earth Filter Mathematica Toolbox V4 User’s Guide: Simplifying, Visualizing and Quick Modeling of All Fission Products from 60 Nuclear Sources

The Earth Filter Mathematica toolbox is a set of user-friendly routines that accurately calculate the ingrowth and decay of any of the 908 fission products for an infinite combination of sources over a large dynamic range and large time scales. This toolbox was developed by Dr. W. P. Brug at Los Alamos National Laboratory and eliminates the need for the user to define and solve differential equations by automatically scanning the ENDF-349 dataset and creating a complete set of differential equations (including delayed neutron branches to both the ground and metastable states) for every fission product. Additionally, this code allows the user to choose the data source for half-lives from several key references (i.e., TORI, IAEA, ENDF or NUDAT.) The complete set of ENDF-349 independent and cumulative fission yields for all 60 isotopic targets and 4 different neutron spectrums are included. This allows nearly any nuclear fission source to be accurately modeled as a linear combination of endpoint sources. The Earth Filter Toolbox has no limit as to the length of split list, and split lists of 5 or more endpoint sources can be quite common.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multiscale and Multimodal Characterization of 2D Titanium Carbonitride MXene

Here, a comprehensive study on the prototype solid solution phase carbonitride MXene Ti 3 CN is conducted using nuclear magnetic resonance, electron spin resonance, total and quasi-elastic neutron scattering, combined with density functional theory-based electronic structure and molecular dynamic calculations. The combination of experiment and theory lead toward rational atomic structural models of Ti 3 CN. The remnant Al ions from the etching process significantly tune the interlayer spacing, distinct from the more typical MXene, Ti 3 C 2 , prepared similarly. Neutron scattering indicates the surface terminations of Ti 3 CN display high oxygen and fluorine concentrations and rather low hydroxyl and hydrogen concentrations. Calculations show that the structure including both the residual Al ions and mixed surface terminations give the best agreement with the measurements. The water molecules in Ti 3 CN are highly immobile, in strong contrast to those in Ti 3 C 2 . The analysis of the electronic structure suggests that the nitride MXene displays higher conductivity than the carbides. Finally, the absence of hydroxyl groups in terminations, the solid-solution in the anion sites, the remnants within layers, and immobile water altogether make the carbonitrides a unique series in the MXene family, implying a further exploration of their exotic properties and applications in energy storage.

36 MATERIALS SCIENCE↗

Constrained Bayesian optimization of criticality experiments

The design of criticality experiments is typically an iterative process that employs a Monte Carlo transport code. The goal is to find a design that optimizes some variable, like the sensitivity of a response to a cross section, while simultaneously ensuring criticality. The high fidelity of the Monte Carlo code is a great asset, but it makes exploring the design space computationally expensive. Herein, we present how a constrained Bayesian optimization algorithm can be used to efficiently design a criticality experiment. It uses Gaussian processes as a surrogate model to probe the design space and to reduce the number of code executions that are needed to find the optimum. Furthermore, we demonstrate constrained Bayesian optimization with a Pu-239/polyethylene solution system and a TEX experiment that is designed for criticality safety validation of a nuclear waste model at the Hanford Site. For both systems, a global optimum was found within 75 Monte Carlo simulations.

42 ENGINEERING↗

Irradiation Vehicles for Evaluating SiC/SiC Cladding Bowing Under Neutron Flux Gradients

Silicon carbide fiber–reinforced silicon carbide matrix (SiC/SiC) composites are among the most promising candidates for long term accident-tolerant nuclear fuel cladding. A key challenge related to their deployment is lateral bowing caused by differential radiation-induced swelling under dose or temperature gradients, which could obstruct coolant flow or interfere with control rod/blade movements. Although previous modeling efforts have predicted bowing behavior in light-water reactor (LWR) environments, experimental validation remains limited, especially at prototypic LWR temperatures. This study addresses that gap by irradiating six reduced-length SiC/SiC cladding tubes (~600 mm) in the High Flux Isotope Reactor (HFIR) at ~300°C, which is representative of LWR conditions. The tubes will be housed in a sealed vessel with an inert gas gap to maintain target temperatures and prevent direct coolant contact. Arranged in three pairs, each set will receive a different radiation dose (low, medium, high), with the central pair receiving ~0.1 displacements per atom (dpa)—the expected dose for peak bowing. The experiment will determine the dose-dependent bowing behavior and validate predictive finite element models. In this work, the tubes are freely suspended from pins to allow for unconstrained bowing; however, we present a concept for introducing localized constraints to represent grid spacer effects. Post-irradiation examination will include swelling measurements and profilometry to assess bowing and compare the results with model predictions. This work aims to confirm the conditions under which maximum bowing occurs so as to improve the reliability of SiC/SiC performance models in nuclear applications.

Russell, Nick [ORNL] (ORCID:0000000296099820)↗

Demonstrating Computational Equivalence Between Continuous and Discrete Adjoint Methods by Calculating Time-Dependent Adjoint Solutions with Neutron Diffusion Models

The continuous adjoint method and the discrete adjoint method are two alternative approaches used to calculate adjoint solutions for adjoint systems. The continuous adjoint method derives adjoint equations analytically from continuous forward equations and then solves the adjoint equations either analytically or numerically in a discretized form whereas the discrete adjoint method calculates the adjoint solutions directly from the discretized forward equations. With regard to the methodology development and calculation procedure, distinct differences are well recognized between the two methods. For certain reasons, both methods are exclusively preferred and commonly used by different computational communities, but limited studies clarify the connections between the two adjoint methods from either of the communities. Herein, this paper demonstrates the computational equivalence between the continuous and discrete adjoint methods by investigating time-dependent adjoint solutions to the two-group neutron diffusion model in nuclear reactor analysis problems using both methods. Adjoint solutions can be used to estimate system parameters for reactor safety analysis. Appropriate final state conditions for the adjoint systems are specified in both of the methods, and the conditions are clarified with proper physical explanations. With the help of an event-based case study on neutron diffusion models, the accuracy of the time-dependent adjoint fluxes obtained from both methods is verified, and the pros and cons of both adjoint methods are examined. More importantly, the computational equivalence of both methods is demonstrated when they are applied to multigroup neutron diffusion systems. The advantage of calculating time-dependent adjoint fluxes by directly solving time-dependent adjoint systems rather than taking steady-state approximations as in common practice is also demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

From chiral effective field theory to perturbative QCD: A Bayesian model mixing approach to symmetric nuclear matter

Constraining the equation of state (EOS) of strongly interacting, dense matter is the focus of intense experimental, observational, and theoretical effort. Chiral effective field theory (𝜒⁢EFT ) can describe the EOS between the typical densities of nuclei and those in the outer cores of neutron stars, while perturbative QCD (pQCD) can be applied to properties of deconfined quark matter, both with quantified theoretical uncertainties. However, describing the full range of densities in between with a single EOS that has well-quantified uncertainties is a challenging problem. Bayesian multimodel inference from 𝜒⁢EFT and pQCD can help bridge the gap between the two theories. In this work, we introduce a correlated Bayesian model mixing framework that uses a Gaussian process (GP) to assimilate different information into a single QCD EOS for symmetric nuclear matter. The present implementation uses a stationary GP to infer this mixed EOS solely from the EOSs of 𝜒⁢EFT and pQCD while accounting for the truncation errors of each theory. The GP is trained on the pressure as a function of number density in the low- and high-density regions where 𝜒⁢EFT and pQCD are, respectively, valid. We impose priors on the GP kernel hyperparameters to suppress unphysical correlations between these regimes. This, together with the assumption of stationarity, results in smooth 𝜒⁢EFT-to-pQCD curves for both the pressure and the speed of sound. We show that using uncorrelated mixing requires uncontrolled extrapolation of at least one of 𝜒⁢EFT or pQCD into regions where the perturbative series breaks down and leads to an acausal EOS. Here, we also discuss extensions of this framework to nonstationary and less differentiable GP kernels, its future application to neutron-star matter, and the incorporation of additional constraints from nuclear theory, experiment, and multimessenger astronomy.

Bayesian methods↗

Modeling and Optimization of a Nuclear Integrated Energy System for the Remote Microgrid on El Hierro

Nuclear microreactors are a potential technology to provide heat and electricity for remote microgrids. There is potential for the microgrid on the island of El Hierro to use a microreactor, within an integrated energy system (IES), to generate electricity and provide desalinated water. This work proposes a workflow for optimizing and analyzing IESs for microgrids. In this study, an IES incorporating a microreactor, thermal energy storage (TES) system, combined heat and power plant, and a thermal desalination plant was designed, optimized, and analyzed using Idaho National Laboratory’s Framework for Optimization of Resources and Economics (FORCE) toolset. The optimization tool, Holistic Energy Resource Optimization Network (HERON), was used to determine the optimal capacity sizes and dispatch for the reactor and thermal energy storage systems to meet demand. The optimized reactor and TES sizes were found to be 11.61 MWth and 58.47 MWhth, respectively, when optimizing the IES to replace 95% of the island’s existing diesel generation needs. A dynamic model of the system was created in the Modelica language, using models from the HYBRID repository, to analyze and verify the dispatch from the optimizer. The dynamic model was able to meet the ramp rates while maintaining reactor power with minimal control adjustments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Insights on Nuclear Energy from Long-Term System Models: Inter-Model Comparison Results

A variety of US electric power sector capacity expansion models are used by decision makers and analysts to evaluate competition among generation, transmission, and storage technologies to meeting the demands of the system. CEMs use least-cost optimization to identify optimal portfolios of investments capable of satisfying all specified requirements. While CEMs are a useful tool to inform pathways to meet future needs, projections can differ significantly between tools for apparently similar scenario assumptions. Differences in model structure, scope, and input assumption contribute to this issue. This work compares model response with harmonization between four modeling teams on issues significant to the representation and development of nuclear energy.

capacity expansion model↗

Applying Constrained Bayesian Optimization to the Design of Critical Experiments

Often when planning a criticality experiment, many design configurations are iteratively investigated with a Monte Carlo transport code. The goal is that the experiment will be optimal with respect to some variable, like the fraction of fissions occurring at a certain energy range, while simultaneously being critical. Unfortunately, the Monte Carlo transport simulations are expensive, which can ultimately limit the number of configurations that can be explored. In this work, we present how Gaussian processes (GPs) can be used as a reduced-order model in a constrained Bayesian optimization (CBO) algorithm to design a criticality experiment. The GPs replace the Monte Carlo transport simulations that explore the design space. The CBO algorithm efficiently identifies new points in the design space to run the Monte Carlo transport code while respecting the criticality constraint. It does so in a manner that both improves the accuracy of the GP and finds the approximate global optimum. We demonstrate the performance of CBO with the design of a Thermal Epithermal eXperiment (TEX) for the criticality safety validation of nuclear waste models of the Hanford Tank Farm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling of the High Temperature Test Facility Using RELAP-7

The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is an electrically heated, helium cooled, experimental facility. The HTTF was modeled using RELAP-7 in both 1-D and 3-D for validation. Due to core symmetry, 1/6th of the HTTF core was modeled in 3-D using the coupled heat conduction and forced convection capabilities of RELAP-7. This served as a validation case for RELAP-7 and its capabilities to model advanced nuclear technologies such as high temperature gas-cooled reactors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Viscosity Measurements of the HLW APPS Glasses

The U.S. Department of Energy is responsible for building the Hanford Waste Treatment and Immobilization Plant (WTP) at the Hanford site in Washington to remediate 56 million gallons of radioactive waste historically stored in 177 underground tanks. The Office of River Protection has requested that the Savannah River National Laboratory (SRNL) contribute in areas of recognized capabilities and expertise for glass waste form development to support successful startup of the WTP. SRNL support of this work is defined in the Task Technical and Quality Assurance Plan.1 This report provides results from viscosity measurements on a series of simulated nuclear waste glasses designed and fabricated at Pacific Northwest National Laboratory (PNNL). The glasses were designated the high-level waste – Aspen Process Performance Simulation (HLW-APPS) study glasses. The data provided in this report is to be used in the development, validation, and implementation of enhanced property/composition models for nuclear waste glasses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Survey of Relevant Data from the MSRP to Guide Development of MSR Chemistry Modeling Benchmarks

The Multiphysics Applications technical area of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is tasked with assessing, demonstrating, and applying NEAMS tools in solving multiphysics problems of nuclear reactors, such as liquid-fueled molten salt reactors (MSRs), which are the focus of this report. MSRs are actively being pursued as a potential candidate for power and heat generation by the nuclear industry. However, the inherent multiphysics nature of liquid-fueled MSRs, stemming from the strong interrelationship of neutronics, thermal fluids, and chemistry phenomena, provides unique challenges in modeling and simulation (mod/sim). Therefore, it will be important to develop mod/sim tools with varying types of multiphysics coupling that depend on the problem. The current work is focused on the initiation and development of MSR chemistry modeling benchmarks useful for validating current and potential future NEAMS tools. Development of such benchmarks include the following actions: 1) Summarize available chemistry data from the Oak Ridge National Laboratory (ORNL) MSR Program (MSRP) including operation of the Molten Salt Reactor Experiment (MSRE) and design studies for the Molten Salt Breeder Reactor (MSBR) concept; 2) Recommend simulation problems in MSR chemistry mod/sim based on (1); 3) Assess the current state of NEAMS tools that may support (2); 4) Demonstrate and validate the capabilities of the NEAMS tools in (3) while providing iterative feedback on future code development activities. The objective of this report is to initiate this effort by completing actions (1) and (2), which are discussed in Section 2. An overview of the relevant available MSRE experimental data is provided with examples of how this data may be useful in chemistry mod/sim problems, with the caveat that most of this data is over 50 years old, therefore some problem details as well as uncertainty estimates are often not provided. In Section 3, the current state of NEAMS tools is assessed for potential use in these mod/sim problems, with considerations for future code development activities, supporting action (3). Future work supporting this project under NEAMS may include a deeper dive into the specific phenomena discussed here including tasks such as the compilation of additional available data, updates in code development activities, and ultimately the demonstration and validation of these tools, supporting action (4).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Pair production as a probe for the dynamics of nuclear fission and α decay

Electron-positron pairs can be produced via the Schwinger mechanism in the presence of strong electric fields. In particular, the fields involved in α decay and nuclear fission are strong enough to produce them. The energy of the e⁺⁢e⁻ pair is related to the relative distance and velocity of the daughter nuclei. Thus, the energy distribution of the produced pairs can give information about the dynamics of the fission and α-decay processes. A neck model of nuclear fission is used to illustrate how the pairs can be used as a probe of the dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Virtual Test Bed (VTB) repository: a library of multiphysics reference reactor models using NEAMS tools

With the next generation of nuclear reactors under development, modeling and simulation (MS) tools are being developed by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in order to support their design, licensing, and future operation. Mirroring the physical test beds currently under construction (i.e., EBR-II and ZPPR), the Virtual Test Bed (VTB) was launched by the National Reactor Innovation Center (NRIC) in collaboration with NEAMS to support the advanced reactor community. This collaborative effort, which involves multiple teams at both Idaho National Laboratory and Argonne National Laboratory aims to use NEAMS tools to model a wide range of reactor designs. Those models are automatically tested to ensure their continued functionality as the tools are further developed. Examples are extensively documented, each acting as a tutorial for applying the relevant NEAMS tools to that reactor design. Currently, five advanced reactor types (with a total of eight specific design variants) are simulated by a variety of different models. These models range from steady-state, core multiphysics simulations to integrated plant analysis during loss-of flow transients. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Peridynamic modeling of cementitious materials for nuclear waste management

Radioactive and hazardous waste generated from fuel processing plants, nuclear reactors, and hospitals, requires effective management strategies. Cementitious materials are widely applied for these needs, serving as structural materials, reactive barriers, or waste forms. In these applications cracking poses a significant risk to performance, driven by inconsistent shrinkage of the components and the varying strength and permeability of their interfaces. Traditional modeling approaches face challenges in representing the complex fracture behavior of cementitious materials due to the reliance on spatial derivatives and difficulties with mesh generation. Here, to overcome these limitations, we employ peridynamics, a novel continuum mechanics formulation that uses integrals to describe mechanical equilibrium, avoiding discontinuities associated with traditional methods. Through incorporation of a bi-linear softening model and quasistatics, an experimentally validated model for Portland cement concrete samples was created. Mechanical parameters, including compressive strength and elastic modulus were validated and variation due to aggregate packing was evaluated. Additionally, sensitivity analysis of the peridynamic parameters to the bulk material properties was established. The results lay the groundwork for evaluating the impact of unique conditions of cementitious waste forms that can be assessed to improve the reliability of waste management strategies.

Aggregates↗