Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Engineering”

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 1,441 records · Page 80

Fission neutron activations of stainless steel and its corresponding elemental components

Characterizing neutron activation of stainless steel is a critical task for many applications, and yet little activation data is available for this alloy. To address the need, a series of neutron activation experiments have been performed using the Flattop critical assembly at the National Criticality Experiments Research Center. Several experimental conditions were probed for complete understanding of the neutron environment, including foil composition and position within the assembly. The experimental results reported here will aid in the production of stainless steel benchmarks for a range of nuclear data applications including those used in nuclear forensics analysis.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Film condensation with high heat fluxes and scaled experiments using pure steam for reactor containment cooling

Condensation tests were performed using a newly developed test facility for scaling the passive containment cooling system (PCCS) to a small modular reactor (SMR). The PCCS of the SMR plays a pivotal role in ensuring greater safety, reliability, and compactness than what is afforded by traditional reactors. Therefore, a well-designed PCCS is essential to SMRs. However, previous studies and test data were unsuitable for scaling, due to high variation in the test geometry and operating conditions. This study intends to close this research gap by using a novel designed scaled test facility consisting of vertical condensing test sections featuring 1-, 2-, and 4-inch-diameter condensing tubes with annular water cooling, and by applying superheated and saturated steam with different steam mass flow ranges of 5–25 g/s. Further, the primary test data, including axial temperatures, mass flow rates, and pressures, were used in conjunction with a standard data reduction method to estimate critical parameters such as heat fluxes, heat transfer coefficients, and condensation rates. These scaled test data would support improving empirical correlations and validating condensation models to identify scaling distortion for SMR PCCSs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

CFD modeling of turbulent air flow in self-heated gyroid TPMS structures: Thermal-hydraulic performance and validation

The application of mathematically derived geometries, such as triply periodic minimal surface (TPMS) lattices, has garnered significant interest across various fields, including the nuclear sector, due to their superior thermal-hydraulic characteristics for heat transfer compared to traditional plain or finned tubes. Here, this study validates a computational fluid dynamics (CFD) model, evaluates different turbulence models and CFD model settings, and performs uncertainty quantification to provide a comprehensive analysis. Despite extensive research on CFD modeling of TPMS lattices, such as gyroid and diamond geometries, there is a notable lack of publicly available literature providing comprehensive details on numerical analysis aspects, including convergence and methodological best practices. This study embarks on a benchmark analysis of a gyroid geometry to evaluate its thermal-hydraulic performance under turbulent flow conditions and scrutinize various CFD model configurations. The main contributions of this work include validating the CFD model, assessing and comparing different turbulence models, and enhancing pressure drop and temperature prediction capabilities. The results aim to support the development of methodologies needed to benchmark and enhance numerical analysis techniques for TPMS lattices. This work seeks to complement the existing body of knowledge, support the development of TPMS reactor concepts, and improve best practices for CFD modeling of TPMS lattices, ultimately advancing methodologies to support future applications in this domain.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

CFD modeling of natural circulation in LiCl-KCl molten salt closed loop

Characterizing flow within a molten salt closed-loop system is crucial for assessing system requirements, evaluating performance, and identifying potential flaws. Direct flow measurement using instrumentation is challenging due to extreme environmental conditions and the limitations associated with measuring molten salt flow under natural convection. Here, this study aims to provide comprehensive insights into the thermal-hydraulic behavior of a closed loop, with a particular focus on temperature distribution and velocity prediction. The Computational Fluid Dynamics (CFD) model demonstrated the capability to effectively simulate and predict both temperature distributions and flow velocities within the molten salt loop. The CFD model's predictive capability was validated by its ability to replicate temperature measurements under varying boundary conditions. The analysis revealed that the CFD model tends to underpredict temperatures in the cold leg and overpredict them in the hot leg, highlighting the need for continuous model refinement and acknowledging the limitations of using a steady-state approach. Furthermore, the potential of using external temperature measurements to estimate internal molten salt temperatures and predict flow velocity was explored, revealing that this approach could introduce up to a 5.5% error in flow velocity calculations. Line probes mapping temperature distributions across the tube's cross-section and molten salt provided valuable insights into temperature gradients, emphasizing the need for a thermal conductivity equation for molten salt with lower uncertainty to achieve more accurate temperature predictions of the system.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

A methodology for domain overlapping coupling of thermal-hydraulic systems

Multi-scale coupling has increasingly drawn attention as a promising approach for modeling thermal systems. Traditional system codes provide fast and robust predictions at the plant scale, while high-fidelity computational fluid dynamics (CFD)-based tools resolve localized flow and heat transfer phenomena with greater accuracy. By combining these complementary strengths, co-simulations enable multi-scale analysis that would otherwise be computationally prohibitive for a standalone CFD code. Here, this work introduces a robust and problem-agnostic domain overlapping (DO) coupling between the system thermal-hydraulic (STH) code System Analysis Module (SAM) and the coarse-mesh CFD code Pronghorn. Both applications belong to the Comprehensive Reactor Analysis Bundle (BlueCRAB) code suite, a code suite in active development at the Idaho National Laboratory (INL), tailored for multi-physics analysis of advanced reactors. Unlike previous approaches, BlueCRAB supports an agnostic interface between codes based on different fidelity, while its coupling formulation can address arbitrary flow geometries with multiple inlets and outlets in coupled components. The implemented method leads to consistent pressure drops, enthalpies, and scalar concentrations between coupled SAM and Pronghorn simulations. The methodology is demonstrated through two verification tests, which ensure the numerical consistency and conservation across the codes, and through one validation test against experimental data. The proposed problems explore different physical aspects inherent to thermal systems, with particular attention given to nuclear reactor analysis. These include buoyancy-driven flows, complex flow patterns, and setups with multiple inlets and outlets, representing challenges in advanced reactor applications.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Overcoming challenges to support US resumption of high specific activity cobalt-60

Domestic production of high specific activity 60 Co was halted after a target rupture in 2012 at the Advanced Test Reactor (ATR). The Isotope Program (IP) within the US Department of Energy (DOE) Office of Science tasked a multilaboratory team of researchers and managers from Oak Ridge and Idaho National Laboratories with the redesign the radioisotope capsule. The objective of this effort was to create a more robust and reliable design, compared to the pre-2012 target. The team successfully completed this task to produce the DOE-IP cobalt (Co) production capsule design. Furthermore, 66 capsules were successfully fabricated by Oak Ridge National Laboratory (ORNL) and delivered to Idaho National Laboratory (INL) for irradiation in the ATR between January 2014 and October 2016. This paper describes the efforts of the team to prepare and disposition the two initial DOE-IP Co production capsules that were processed in March 2020. These efforts include performing accurate production predictions, experimentally validating predictions with assay measurements, shipping with the Orano-furnished Battelle Energy Alliance Research Reactor shipping package, and disassembling capsules at the isotope vendor site.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Damage development in neutron-irradiated concrete in a test reactor: Hygro-thermal and mechanical simulations

This paper reports the development of a 3D mesoscale hygro-thermal-mechanical simulation approach to predict damage in concrete irradiated in a test reactor. This framework, developed in MOOSE, considers the effects of elevated temperature, moisture content, and high neutron fluence (energy threshold, E > 0.1 MeV) on the mortar and aggregates separately. The first-stage simulation implements hygro-thermal analysis to determine the temperature and RH inside the specimen as a function of imposed radiation energy. These are used as inputs to the second stage, which considers radiation-induced volumetric expansion (RIVE) of aggregates, and creep, shrinkage, and stress-strain response of mortar to predict the expansion, stresses, and damage in specimens made using different coarse aggregates and subjected to different irradiation times. The irradiation time-dependent damage in the mortar is expressed using an isotropic damage parameter. This multi-physics model serves as a predictive tool for damage quantification in concrete due to neutron irradiation.

36 MATERIALS SCIENCE↗

The MOOSE fluid properties module

The Fluid Properties module within the Multiphysics Object-Oriented Simulation Environment (MOOSE) is used to compute fluid properties for numerous applications, ranging from nuclear reactor thermal hydraulics to geothermal energy. Those applications drove the development of the module to enable numerous different fluid equations of states, property lookups with primitive and conserved flow variable to cater to pressure and density-driven solvers, and an object-oriented design facilitating expansion and maintenance. Each fluid property is implemented in its own class but inherits capabilities such as automatic differentiation, automated out-of-bounds handling or variable conversion capabilities. Here, this paper presents the module, its design, its user and developer interface, its content in terms of fluids and properties, and several of its applications showing its major role in the MOOSE simulation ecosystem.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The influence of cooling rate on condensation of iron, aluminum, and uranium oxide nanoparticles

Fundamental observations of particle size distributions are needed to develop models that predict the fate and transport of radioactive materials in the atmosphere following a nuclear incident. The extent of material transport is influenced by the time scales of particle formation processes (e.g., condensation, coagulation). In this study, we investigated the influence of cooling time scales on size distributions of uranium, aluminum, and iron oxide particles that are synthesized separately under identical run conditions inside the controlled environment of an argon plasma flow reactor. Two distinct temperature distributions are imposed along the flow reactor by varying the argon flow rate downstream of the plasma torch. The vaporized reactants of uranium, aluminum, and iron are cooled from about 5000K to 1000K before they are collected on silicon wafers for ex situ scanning electron microscope analysis. The microscope images show that the sizes of the largest aluminum and iron oxide particles heavily depend on the cooling time scales, whereas significant size variation with cooling rate is not observed for uranium oxide particles. In addition, the size distribution of aluminum oxide particles exhibits the broadest range among all three metal oxides studied. We performed simulations of particle size distributions using a kinetic model that couples gas phase oxidation chemistry with particle formation processes, including nucleation, condensation, and coagulation. The model results demonstrate the strong sensitivity of particle size distribution to different cooling histories (i.e., temperature vs residence time) along the flow reactor. In conclusion, the kinetic model also helps identify directions for future research to improve the predictions.

36 MATERIALS SCIENCE↗

Meshless discretization of the discrete-ordinates transport equation with integration based on Voronoi cells

The time-dependent, gray, linear radiation transport equation is discretized using the meshless local Petrov-Galerkin method with reproducing kernels. The integration is performed using a Voronoi tessellation, which creates a partition of unity that only depends on the position and extent of the kernels. The resolution of the integration automatically follows the particles and requires no manual adjustment. The discretization includes streamline-upwind Petrov-Galerkin stabilization to prevent oscillations and improve numerical conditioning. The angular quadrature is selectively refineable to increase angular resolution in chosen directions. The time discretization is done using backward Euler. The transport solve for each direction and the solve for the scattering source are both done using Krylov iterative methods. The results indicate first-order convergence in time and second-order convergence in space for linear reproducing kernels.

97 MATHEMATICS AND COMPUTING↗

Iron-rich microstructure records of high temperature multi-component silicate melt behavior in nuclear fallout

Above-ground nuclear explosions that interact with the surface of the earth entrain materials from the surrounding environment, influencing the resulting physical and chemical evolution of the fireball, which can affect the final chemical phase and mobility of hazardous radionuclides that are dispersed in the environment as fallout particles. The interaction of iron with a nuclear explosion is of specific interest due to the potential for iron to act as a redox buffer and because of the likelihood of significant masses of metals to be present in urban environments. We investigated fallout from a historic surface interacting nuclear explosion conducted on a steel tower and report the discovery of widespread and diverse iron-rich micro-structures preserved within the samples, including crystalline dendrites and micron-scale iron-rich spheres with liquid immiscibility textures. We assert these micro-structures reflect local redox conditions and cooling rates and can inform interpretation of high temperature events, enabling new insights into fireball condensation physics and chemistry when metals from the local environment (i.e. structural steel) are vaporized or entrained. Finally, these observations also significantly expand the availability of silicate immiscibility datasets applicable to rapidly quenched systems such as meteorite impact melt glass.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Multi-material additive manufacturing of aluminum 6061-T6 alloy with stainless steel 304: Suppression of intermetallic compounds and interface growth mechanism

Fabricating integrated aluminum-steel structures is difficult because of metallurgical incompatibilities that promote intermetallic compound (IMC) formation at their interfaces. This work introduces a practical route for additively manufacturing 6061-T6 aluminum alloy (AA6061-T6) directly onto 304 stainless steel (SS304) using a high shear strain rate-assisted interfacial bonding mechanism. Through the friction extrusion deposition-based additive manufacturing method, both single- and multi-layer deposits of AA6061-T6 were successfully fabricated on SS304, yielding a robust hybrid multi-material structure. Comprehensive analyses of deposition quality, interface porosity, and bonding performance showed that the interface achieved a tensile strength exceeding 154 MPa under quasi-static loading. Microstructural observations revealed that the deposited aluminum layers experienced severe plastic deformation, resulting in pronounced grain refinement. Importantly, the interfacial zone was found to be free of brittle IMCs, instead containing a thin amorphous layer that evolved through a non-linear reaction-zone growth law. This solid-state additive strategy establishes a promising pathway for lightweight structural systems, nuclear energy component cladding, and multifunctional engineering components, redefining how dissimilar metals can be integrated for advanced applications.

Aluminum alloy↗

Gas-solid reaction pathway for chlorination of rare earth and actinide metals using hydrogen and chlorine gas

Hydriding and chlorination of buttons of metallic cerium have been tested as a head-end step to separating metal impurities from actinides. The objective here was to achieve complete conversion to chlorides while maintaining high surface area via small particles and open porosity. Hydriding by reaction with 100% H 2 gas at 573 K was successful at reducing cm-scale metal buttons into small particles/fragments. This step was followed by reaction with flowing 99.5% anhydrous Cl 2 gas at temperatures ranging from 423 to 573 K for 60 min. In experiments reported here, up to 92% conversion to cerium chloride was achieved from reacting with the anhydrous chlorine gas. The final chloride product has high porosity coupled with small particle sizes (microns to 1600 μm), which makes it suitable for subsequent purification processes. Higher chlorine concentrations in the gas favor higher conversion to chloride. The optimal temperature to achieve high chlorination is 523 K. Oxygen is a minor impurity in the chlorine gas which appears to prevent complete chlorination from being achieved.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Hydrogen trapping energetics at BCC iron-helium interfaces

Density functional theory (DFT) calculations have been performed to assess the trapping and segregation strength of hydrogen (H) to helium (He) interfaces in BCC (Fe) as a function of He density and surface orientation. The He density ranges from 0 to 2 He/V which corresponds to equilibrium bubbles that are expected to form in the structural components in a fusion reactor environment. The DFT calculations consist of a slab of BCC metal and an initial lattice of FCC He. Here, the report on the binding energy of H to these interfaces as well as the migration barriers into and away from the surface, which together provide information on the de-trapping energy for H at He bubbles. It was found that the binding energy of H to He-Fe interfaces decreases with increasing He density. Migration barriers are also modified due to the He density, which are sensitive to the surface orientation as well as the diffusion pathway. These results provide valuable first-principles energetics that are necessary for higher-order, mesoscale models, and provide insight into the extent to which He bubbles may trap tritium in fusion structural materials.

42 ENGINEERING↗

Verification of Bison fission product species conservation under TRISO reactor conditions

When assessing the reliability and predictive capabilities of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. During this process, numerical results in a discrete solution are compared to the analytical solution of the mathematical model. Here, in this paper, the code verification process is applied to one-dimensional spatiotemporal problems that exercise partial differential equation governing the conservation of fission product species (or mass diffusion). Numerical experiments were performed in the Bison fuel performance code to evaluate its predictive capability under various TRISO reactor conditions such as base irradiation and safety heating test conditions for either short- or long-lived fission product species, as well as a case concerning evaporation from the outer surface of a particle. The code predictions were compared with the expected exact results obtained from the analytical expressions, and the fact that they demonstrate the correct analytical behavior provides strong evidence of proper numerical algorithm implementation.

07 ISOTOPE AND RADIATION SOURCES↗

Determination of a surrogate for plutonium electrorefining

Conducting research experiments on plutonium electrorefining is difficult due to the significant hazards and regulations associated with nuclear materials. Finding a surrogate for plutonium electrorefining studies would enable more fundamental research to be conducted. Potential surrogates were first identified by determining the physical properties required to conduct electrorefining at the same conditions commonly used in plutonium electrorefining, a molten metal and molten CaCl 2 at 1123 K. Ce-CeCl 3 , In-InCl 3 , and Pb-PbCl 2 were the only potential surrogates identified using these constraints. Sn-SnCl 2 was also tested at these same conditions. More potential surrogates were identified by changing the matrix salt and operating temperature. This expanded the potential surrogate list to also include Zn-ZnCl 2 , Sn-SnCl 2 , and Bi-BiCl 3 . Zn-ZnCl 2 was used with the LiCl-CaCl 2 (65:35 mol%) eutectic at 773 K. Sn-SnCl2 and Bi-BiCl 3 were used with the LiCl-KCl-CaCl 2 (50.5:44.2:5.3 mol%) eutectic at 673–773 K. Ce electrorefining in molten CaCl 2 resulted in a difficult to separate colloid mixture of Ce, Ca and Cl. Electrorefining rates for In in molten CaCl 2 were too slow due to InCl 3 volatilizing out of the molten salt. Only trace amounts of SnCl 2 was retained in the CaCl 2 at 1123 K resulting in impractical electrorefining rates. Zn metal product was successfully collected in the LiCl-CaCl 2 eutectic molten salt, but the metal obtained did not coalesce into one piece. Sn and Bi were successfully electrorefined in the LiCl-KCl-CaCl 2 eutectic molten salt and coalesced into product rings with high yields and coulombic efficiencies. Finally, while a surrogate could not be identified using the same conditions as plutonium electrorefining, two possible surrogates, Sn-SnCl 2 and Bi-BiCl 3 , were found that could imitate the physical configuration (i.e., molten salt on top of molten metal) of plutonium electrorefining at a reduced temperature using the eutectic LiCl-KCl-CaCl 2 salt at 673–773 K in place of CaCl 2 at 1123 K.

36 MATERIALS SCIENCE↗