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

Studies of Single-Transverse Kinematic Variables for Neutrino Interactions on Argon

It is now well-known that a variety of theoretically-challenging nuclear effects play an important role in determining neutrino-nucleus interaction cross sections. Deficiencies in existing models of these effects represent a leading source of systematic uncertainty for neutrino oscillation experiments. In recent neutrino cross section measurements published by the T2K and MINERvA experiments, the study of kinematic observables called Single-Transverse (Kinematic) Variables (STVs) in the plane transverse to the neutrino beam have been demonstrated to provide a powerful tool on nuclear modeling uncertainties. In the near future, the SBN and DUNE experiments will both leverage the novel Liquid Argon Time Projection Chamber technology to pursue high-precision measurements of neutrino oscillations. As the interpretation of the results from these measurements will require a detailed understanding of neutrino-argon scattering cross sections and given the scarcity of available neutri no-argon data, a measurement of Single-Transverse Kinematic Variables distributions in neutrino-argon scattering is of high interest for the neutrino physics community. This thesis presents the first detailed comparison of model predictions for distributions of these observables obtained using four standard neutrino event simulation tools (GENIE, NuWro, GiBUU and NEUT). Results are shown for simulations of muon neutrinos from the Fermilab Booster Neutrino Beam interacting in the MicroBooNE LArTPC. Opportunities for discrimination between nuclear models with respect to future measurements is highlighted. Furthermore, first calculations of STVs from experimental argon data are shown. Data corresponding to 4.449x10^19 POT of exposure were collected using the MicroBooNE LArTPC in the Fermilab Booster Neutrino Beam.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Review and Modeling of Integrated Energy Systems with Nuclear Reactor Coupled Desalination and District Heating

Detailed reviews of a past advanced nuclear reactor based integrated energy system, as well as other nuclear reactor and fossil fuel based integrated energy systems have been performed for this work. Review of the utilization of heat from nuclear reactors for various applications and cogeneration has been done. The heat can be utilized by extraction of the steam from the turbine while the steam is still at a desired temperature. While use of nuclear process heat for district heating in countries like Finland, France, China, Poland, and elsewhere is discussed, more focus of the review has been given on nuclear desalination processes. Integrated energy systems (IES) where distinct types of reactors like PWR, BWR, sodium cooled fast reactor, heavy water reactor and other advanced reactors are coupled with various nuclear desalination processes like multi-effect distillation (MED), multi-stage flashing (MSF) and reverse osmosis (RO) methods have been discussed. The nuclear desalination plant at Aktau has been discussed in more detail due to its decades of successful operation. The IES of the Aktau plant coupled with 5-effect MED desalination plant has been taken as a reference for modeling the Open Modelica (OM) based IES of this work. Here, the OM IES model shows good agreement with the MED plant output of Aktau and can be extended for future applications of IES.

42 ENGINEERING↗
DOE OSTI2023

Multiscale modeling for high burnup structure formation and associated pulverization

This report summarizes the lower length scale modeling work performed in the fiscal year 2023 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to capture the microstrutural evolution and associated pulverization criteria in the high burnup regions. This year, the resolution mechanisms within the cluster dynamics code has been updated and it influence on the bubble growth in HBS regions at the mesoscale is studied. To improve the accuracy of phase-field models of fission gas bubble growth, a new Helmholtz free energy for high-density Xe gas is derived from a virial equation of state. Two previously used strategies for representing net vacancy production in phase-field models of fission gas bubble growth are compared with each other and with analytical models of bubble growth. The vacancy source only model is found to be more convenient to parameterize realistically compared with the vacancy source+sink model. Furthermore, the phase-field-fracture simulation have been performed using MD- informed failure stress values and realistic HBS structure obtained from the phase-field simulations. We also present the uncertainty bands on prediction of the critical stress to account for the effect of the lower length scale variabilities on the failure criteria at the mesoscale. It is observed that pulverization may occur in partially restructured regions with restructuring fraction as low as 17%. In light of this, an update to the BISON’s pulverization criteria is recommended.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1993

Nuclear Engine System Simulation (NESS). Version 2.0: Program user's guide

This Program User's Guide discusses the Nuclear Thermal Propulsion (NTP) engine system design features and capabilities modeled in the Nuclear Engine System Simulation (NESS): Version 2.0 program (referred to as NESS throughout the remainder of this document), as well as its operation. NESS was upgraded to include many new modeling capabilities not available in the original version delivered to NASA LeRC in Dec. 1991, NESS's new features include the following: (1) an improved input format; (2) an advanced solid-core NERVA-type reactor system model (ENABLER 2); (3) a bleed-cycle engine system option; (4) an axial-turbopump design option; (5) an automated pump-out turbopump assembly sizing option; (6) an off-design gas generator engine cycle design option; (7) updated hydrogen properties; (8) an improved output format; and (9) personal computer operation capability. Sample design cases are presented in the user's guide that demonstrate many of the new features associated with this upgraded version of NESS, as well as design modeling features associated with the original version of NESS.

Pelaccio, Dennis G.↗
DOE OSTI2024

SAM Code Improvements for Molten Salt Reactor Transient and Species Transport Modeling

The SAM code is under development and supported by DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. These advanced reactor concepts incorporate novel and improved approaches to achieve safety and economic feasibility. This report summarizes FY24 efforts in addressing the modeling gaps in SAM for molten salt reactor (MSR) applications. These efforts focused on the solver performance related to species transport, usability of the flowing decay heat and delayed neutron precursor capabilities, and modeling gaps related to fission product transport throughout the system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Efficient Multi-Scale, Thermo-Hydro-Mechanical Numerical Model for Simulation of Long-Term Stability of Rooms for Nuclear Waste Disposal

This project demonstrated the feasibility of developing a numerical modeling tool that provides an efficient and unique approach to analyze the geomechanical performance of a deep geological repository (DGR) for high-level nuclear waste. The numerical modeling tool is physics-based and models the thermo-hydro-mechanical (THM) processes during the required time (up to one million years) for a DGR. The unique capability of the developed code is explicit representation of fracturing and a discrete fracture network (DFN) in the model that includes the entire repository layout for simulated time of one million years. Therefore, it can potentially address one of the challenges of analysis of a DGR, which is efficient numerical simulation of a multi-scale, physics-based models.

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

X-Ray Energy Deposition Model for Simulating Asteroid Response to a Nuclear Planetary Defense Mitigation Mission

Abstract In the event of a potentially catastrophic asteroid impact, with sufficient warning time, deploying a nuclear device remains a powerful option for planetary defense if a kinetic impactor or other means of deflection proves insufficient. Predicting the effectiveness of a potential nuclear deflection or disruption mission depends on accurate multiphysics simulations of the device's X-ray energy deposition into the asteroid and the resulting material ablation. The relevant physics in these simulations span many orders of magnitude, require a variety of different complex physics packages, and are computationally expensive. Having an efficient and accurate way of modeling this system is necessary for exploring a mission's sensitivity to the asteroid's range of physical properties. To expedite future simulations, we present a completed X-ray energy deposition model developed using the radiation-hydrodynamics code Kull that can be used to initiate a nuclear mitigation mission calculation. The model spans a wide variety of possible mission initial conditions: four different asteroid-like materials at a range of porosities, two different source spectra, and a broad range of radiation fluences, source durations, and angles of incidence. Using blowoff momentum as the primary metric, the model-initiated simulation results match the full radiation-hydrodynamics results to within 10%.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Correlating Protein Dynamics and Catalytic Activity of a Model Hydrogenase Using Paramagnetic and Biological Nuclear Magnetic Resonance Spectroscopy

Rational catalyst design remains a significant challenge, with electronic structure, steric, and electrostatic effects known to contribute to activity. Recently, dynamics has been recognized as another factor that impacts catalysis, though identifying and predicting these effects has remained out of reach. Nickel-substituted rubredoxin (NiRd), a protein-based mimic of a hydrogenase enzyme, serves as a model catalytic system in which dynamics can be systematically investigated with respect to activity. While over 30 secondary-sphere mutants of NiRd have been shown to be catalytically active, no significant correlation was observed between the rates and catalytic overpotential or electronic structure, prompting questions about the protein-derived factors that modulate activity. Here, in this work, NMR spectroscopy was used to investigate the roles of substrate accessibility, protein dynamics, and protein stability in controlling catalysis. Significant paramagnetic effects from the nickel center (S = 1) isolate the methylene proton resonances of the metal-coordinating cysteine residues. The sensitivity of resonance positions and linewidths to local environment offers an opportunity to study dynamical molecular changes around the metal center with high resolution. Machine learning algorithms were employed to identify correlations between the catalytic activity and the paramagnetic NMR spectra. These analyses revealed spectroscopic features of specific cysteine protons that report on catalytic overpotential and increased turnover rates, which are further supported by the results obtained using high-field NMR techniques. Collectively, these studies indicate the potential for multifrequency NMR techniques to resolve key contributors to catalytic activity and highlight the importance of local and outer-sphere dynamics.

Protein Engineering↗
DOE OSTI2021

Numerical study of the chemo-mechanical behavior of FEBEX bentonite in nuclear waste disposal based on the Barcelona expansive model

Experimental studies show that compacted bentonite used as a backfill material for nuclear waste repository experiences strong coupling between chemical and mechanical processes. Here, we use a dual-structure expansive soil model, referred to as the Barcelona Expansive Model (BExM), to predict the behavior of bentonite buffer in subsurface emplacement tunnels for high-level radioactive waste. After emplacement, the bentonite is subjected to complex and coupled Thermal-Hydraulic-Mechanical-Chemical (THMC) processes. The BExM constitutive model is implemented in a multi-phase reactive transport and geomechanics simulator, TOUGHREACT-FLAC3D, and the model is verified with one benchmark test on FEBEX bentonite. We utilize a one-way Chemo-Mechanical (C-M) coupling approach, in which chemical changes affect the mechanical behavior of bentonite through the BExM linked with the evolution of mass fraction of smectite, exchangeable cation concentration, and ionic strength via osmotic suction. The parameters of the new coupled model for FEBEX bentonite are calibrated against a series of laboratory experiments with various salinity solutions. Finally, coupled THMC modeling is conducted for a generic argillite repository with bentonite buffer under high temperature, focusing on the long-term chemical change and its effect on the mechanical process.

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

Entanglement rearrangement in self-consistent nuclear structure calculations

Background: Entanglement plays a central role in a diverse array of increasingly important research areas, including quantum computation, simulation, measurement, sensing, and communication. Extensive suites of investigations have been performed to better understand entanglement in atomic and molecular quantum many-body systems, while the exploration of entanglement in the structure of nuclei and their reactions is presently in its infancy. Purpose: The goal of this work is to begin investigating the entanglement properties of nuclei from first-principles nuclear many-body calculations. We attempt to identify common features and emergent structures of entanglement that could ultimately lead to new and natural many-body schemes. With an eye toward quantum accelerators in future hybrid-supercomputers, criteria for partitioning nuclear many-body calculations into quantum and classical components may provide advantages in future large-scale computations. Along the way we look for explanations of the relative success of phenomenological models such as the nuclear shell model, and for better ways to match to low-energy nuclear effective field theories and lattice QCD calculations to nuclear many-body techniques that are based upon entanglement. Method: We explore the entanglement between single-particle states in 4 He and 6 He. The patterns of entanglement emerging from different single-particle bases are compared, and possible links with the convergence of observables are explored, in particular, ground-state energies. The nuclear wave functions are obtained by performing active-space no-core configuration-interaction calculations using a two-body nucleon-nucleon interaction derived from chiral effective field theory. Entanglement measures within single-particle bases exhibiting different degrees of complexity are determined, in particular, harmonic oscillator (HO), Hartree-Fock (HF), natural (NAT) and variational natural (VNAT) bases. Specifically, single-orbital entanglement entropy, two-orbital mutual information, and negativity are studied. Results: The entanglement structures in 4 He and 6 He are found to be more localized within NAT and VNAT bases than within a HO basis for the optimal HO parameters we have worked with. In particular a core-valence structure clearly emerges from the full no-core calculation of 6 He. Here, the two-nucleon mutual information shows that the VNAT basis, which typically exhibits good convergence properties, effectively decouples the active and inactive spaces. Conclusions: Measures of one- and two-nucleon entanglement are found to be useful in analyzing the structure of nuclear wave functions, in particular the efficacy of basis states, and may provide useful metrics toward developing more efficient schemes for ab initio computations of the structure and reactions of nuclei, and quantum many-body systems more generally.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A mechanistic model of a PWR-based nuclear power plant in response to external hazard-induced station blackout accidents

Natural hazard-induced nuclear accidents, such as the Fukushima Daiichi Accident that occurred in Japan in 2011, have significantly increased reactor safety studies in understanding nuclear power plant (NPP) responses to external hazard events such as earthquakes and floods. Natural hazards could cause the loss of offsite power in nuclear power plants, potentially leading to a Station Blackout (SBO) accident that significantly contributes to the overall risk of nuclear power plant accidents. Despite the fact that extensive research has been conducted on the station blackout accident for nuclear power plant, further understanding of these events is needed, particularly in the context of the dynamic nature of external hazards such as external flooding. This paper estimates the progression of station blackout events for a generic pressurized water reactor (PWR) in response to external flooding events. The original RELAP5-3D model of the Westinghouse four-loop design pressurized water reactor was adopted and modified to simulate the external flood-induced station blackout accident, including the short-term and long-term station blackout scenarios. A sensitivity analysis of long-term station blackout, examining reactor operation times and analyzing key parameters over time, was also conducted in this work. The results of the analyses, especially the critical timing parameters of key event sequences, provide useful insights about the time during the external flooding event, which is important for plant operators to make timely decisions to prevent potential core damage. This paper represents significant progress toward developing an integrated risk assessment framework for further identifying and assessing the effects of the critical sources of uncertainties of nuclear power plant under external hazard-induced events.

Liu, Tao↗
DOE OSTI2024

Nuclear mass radius and pressure in the Skyrme model

We compute the mass radius, scalar radius, tensor radius, baryon number radius, and mechanical radius of nuclei with baryon number B = 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 , 32, 108 in the Skyrme model. The relations between these radii and the nuclear gravitational form factors are investigated. We also compute the ‘pressure’ distribution and find that it is negative in the core region for all the nuclei with B > 1 . This suggests that the way mechanical stability is achieved in nuclei is qualitatively different than in the nucleon. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Improvement of mechanistic fuel-cladding chemical interaction modeling in BISON

This report describes work performed during FY2024 under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to inform and improve mechanistic models of fuel-cladding chemical interaction (FCCI) in metallic fuel. For fuel-side FCCI, atomistic simulations were performed to determine the diffusivity of iron (Fe) in the $\alpha$ and $\gamma$ phases of uranium (U). A model of liquid penetration of cladding due to melting of the fuel-side FCCI region was updated to account for the finite size of the FCCI region, and the model was validated through comparison with tests performed in the Fuel Behavior Test Apparatus (FBTA). For cladding wastage formation, a reduced-order model was improved by comparison with a multi-scale mechanistic model to better quantify the ROM parameters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Enhanced production of strange baryons in high-energy nuclear collisions from a multiphase transport model

Here, we introduce additional coalescence factors for the production of strange baryons in a multiphase transport (AMPT) model in order to describe the enhanced production of multistrange hadrons observed in Pb + Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV at the Large Hadron Collider (LHC) and Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV at Relativistic Heavy-Ion Collider (RHIC). This extended AMPT model is found to also give a reasonable description of the multiplicity dependence of the strangeness enhancement observed in high multiplicity events in $\textit{pp}$ collisions at $\sqrt{s}$ = 7 TeV and $\textit{p}$-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. We find that the coalescence factors depend on the system size but not much on whether the system is produced from A + A or $\textit{p}$ + A collisions. Here, the extended AMPT model thus provides a convenient way to model the mechanism underlying the observed strangeness enhancement in collisions of both small and large systems at RHIC and LHC energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2011

Statistical Models of Fracture Relevant to Nuclear-Grade Graphite: Review and Recommendations

The nuclear-grade (low-impurity) graphite needed for the fuel element and moderator material for next-generation (Gen IV) reactors displays large scatter in strength and a nonlinear stress-strain response from damage accumulation. This response can be characterized as quasi-brittle. In this expanded review, relevant statistical failure models for various brittle and quasi-brittle material systems are discussed with regard to strength distribution, size effect, multiaxial strength, and damage accumulation. This includes descriptions of the Weibull, Batdorf, and Burchell models as well as models that describe the strength response of composite materials, which involves distributed damage. Results from lattice simulations are included for a physics-based description of material breakdown. Consideration is given to the predicted transition between brittle and quasi-brittle damage behavior versus the density of damage (level of disorder) within the material system. The literature indicates that weakest-link-based failure modeling approaches appear to be reasonably robust in that they can be applied to materials that display distributed damage, provided that the level of disorder in the material is not too large. The Weibull distribution is argued to be the most appropriate statistical distribution to model the stochastic-strength response of graphite.

Nemeth, Noel N.↗
DOE OSTI2026

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗
DOE OSTI2022

Optimization of the deep neural network parameters for generating homogenized fuel assembly data for nodal codes

Homogenized fuel assembly (FA) data is a typical input data for nodal codes. Generating that data, however, could be time-consuming. One of promising ways to mitigate the computational burden of generating macroscopic cross-sections is to use trained artificial neural network (ANN) models for predicting nuclear data. However, there is a challenge to make the model support variable FA geometry. In this work, two most common types of FA were combined in one ANN model. Since there could be multiple ways of converting 2-dimensional FA data into 1-dimensional input vector for ANN, three different approaches of data flattening were evaluated. The input parameters included each fuel pin enrichment, fuel temperature, moderator temperature and boron concentration. The output parameters were 2-group macroscopic cross-sections (XS) and pin power distribution (HFF). A fully connected deep neural network (DNN) model was trained and tested using pre-generated data obtained with lattice physics code STREAM. The results of this study showed no statistically significant difference in the accuracy of XS and HFF generation for all 3 tested input vector orders. This means that fully connected DNN for XS generation demonstrated input sequence invariance. Results of comparing predicted XS data with reference solutions were found sufficiently close considering the reduction of computation time offered by ANN. Mean relative difference (MRD) for all output XS parameters was found below 0.7%, while HFF MRD was found higher compared to XS values, in some cases slightly exceeding 1%, mostly near guide tube locations. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗