Search NASA⌕ Search

SEARCH · Search NASA

Results for “Engineering (General)”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

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↗
DOE OSTI2025

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↗
DOE OSTI2025

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↗
DOE OSTI2020

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↗
DOE OSTI2020

Modelling dislocation-graphene interactions in a BCC Fe matrix by molecular dynamics simulations and gradient plasticity theory

We report graphene nanosheets (GNS) can enhance the strength and ductility of metal-based composites as they can obstruct the propagation of dislocations. The present article employs Molecular Dynamics (MD) simulations to investigate dislocation-GNS interaction mechanisms and possible influencing factors, including the number of GNS layers, the thickness of the metallic amorphous layer and the C - C bond strength. The results indicated that the shear strength of the metal/GNS interface and the bending stiffness of GNS determined the ability of GNS to block dislocation transmission. A physically based phenomenological parameter that can capture such dislocation-GNS interactions is the mechanical interface energy that has been put forth within gradient plasticity. By fitting the theoretical expressions to the simulation data, it was possible to obtain estimates for the mechanical interface energy for the GNS. It was found that increasing the GNS layers and adding an amorphous layer resulted in a strengthening in the stress–strain response and increased the value of this interfacial parameter. This indicates that the mechanical interfacial energy can be a unified measure for capturing and tuning the strength of various interfaces such as grain boundaries, GNS, amorphous-crystalline interface and bimetal interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2021

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↗
DOE OSTI2021

An integrated numerical model for coupled poro-hydro-mechanics and fracture propagation using embedded meshes

Integrated models for fluid-driven fracture propagation and general multiphase flow in porous media are valuable to the study and engineering of several systems, including hydraulic fracturing, underground disposal of waste, and geohazard mitigation across such applications. This work extends the coupled model multiphase flow and poromechanical model of Ren et al. (2018) to admit fracture propagation (FP). The coupled XFEM-EDFM scheme utilizes a separate fracture mesh that is embedded on a static background mesh. The onset and dynamics of fracture propagation are governed by the equivalent stress intensity factor (SIF) criterion. A domain-integral method (J integral) is applied to compute this information. Additionally, an adaptive time-marching scheme is proposed to rapidly restrict and grow temporal resolution to match the underlying time-scales. The proposed model is verified with analytical solutions, and shows the capability to accurately and adaptively co-simulate fluid transport and deformation as well as the propagation of multiple fractures.

42 ENGINEERING↗
DOE OSTI2021

Investigation of grid-based vorticity-velocity large eddy simulation off-body solvers for application to overset CFD

Accurately predicting unsteady wakes and vortex-dominated flows is essential to a wide range of engineering applications, including aircraft, rotorcraft, shipboard operations, bio-inspired unsteady flight and propulsion, wind turbines, and urban flows. While current CFD software can model the complete flow field and wake system, the computational costs incurred in high Reynolds number unsteady turbulent flow simulations often remain prohibitive for routine engineering use, particularly for applications involving moving components. Prior work has demonstrated that by adopting a vorticity-velocity formulation in a grid-based off-body flow solver (VorTran-M and VorTran-M2) one can lower these costs by several orders of magnitude when compared to conventional approaches. This paper describes the extensions made to VorTran-M2 to support turbulent flows, and associated benchmarking activity to assess its performance for problems involving strong stretching and diffusion processes, whose competing contributions to the vorticity field are core drivers of turbulent flow evolution. Predictions are presented for: (i) the Kida-Pelz problem whose inviscid form is of mathematical interest due to its apparent formation of singular flow in finite time; and (ii) the Taylor Green vortex arrangement, which has been extensively studied as a fundamental simulation challenge in the turbulent modeling community. Here, the results are used to evaluate the overall predictive ability and performance of two sub-grid scale models incorporated into VorTran-M2. Results indicate that the computational cost savings seen previously for inviscid and convection dominated problems extend to turbulent flow simulations supporting the viability of VorTran-M2 as a low cost means for accurately modeling the far-field and background flow, particularly when long duration vorticity evolution is of interest.

42 ENGINEERING↗
DOE OSTI2021

Transverse impact by RCCs on S-glass and Kevlar® FRC strips

This study aims to isolate interactions between plies and tows and reveal fundamental physics involved in transverse impact on fiber-reinforced composite (FRC) structures. Composite strips were sectioned from large panels and characterized by optical photography, three-dimensional synchrotron X-ray computed tomography, and scanning electron microscopy (SEM). Each strip was impacted perpendicularly by a right circular cylinder (RCC) projectile at a velocity ranging from ~ 150 to 600 m/s. The global strip behavior, as well as localized deformation and failure of the strip near the projectile corner, were both captured by high-speed optical imaging. S-Glass FRC strips were observed to fail in tension ahead of the RCC projectiles’ flat surfaces while Kevlar® FRC strips fractured at the projectile corners. The concept of critical velocity region previously used for impact on yarns was introduced to define different failure modes of each composite strip type. The strip damage extent was found to increase with the impact velocity and reach the maximum at the upper limit of the critical velocity region. Above the critical velocity region, the damage extent decreased with impact velocity. Wave propagations and load histories in the composite strips during impact were quantified and compared with Smith’s theory. Finally, critical velocities of single fibers, yarns, and composite strips and ballistic limits of single-ply and multi-ply composite panels were compared to provide insight into the design of impact-resistant fabrics and composites.

42 ENGINEERING↗
DOE OSTI2024

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↗
DOE OSTI2026

The high explosives & affected targets (HEAT) dataset

Artificial Intelligence (AI) surrogate models offer a computationally efficient alternative to full-physics simulations, yet no existing datasets are publicly available for training, testing, and validation of machine learning models of the dynamics of high-explosive driven shocks through multiple materials. Shock propagation through materials is a computationally challenging problem because simulations must include material-specific equations of state (EOS) along with descriptions of other physical processes such as plastic deformation, phase change, damage processes, fluid instabilities, and multi-material interactions. Shocks are typically initiated by high-velocity impacts or explosive loading. The latter case necessitates the addition of models of reactive materials to represent high-explosive (HE) detonation. Here, to address the lack of an expansive dataset for multi-material shock propagation in the AI/ML community, we present the High-Explosives and Affected Targets (HEAT) Dataset. HEAT is a physics-rich collection of two-dimensional, cylindrically symmetric, simulations generated using an Eulerian, multi-material, shock-propagation code developed at Los Alamos National Laboratory. The dataset includes two partitions: (1) the expanding shock-cylinder (CYL) simulations, Figs. 1, and (2) the Perturbed Layered Interface (PLI) simulations, Fig. 2. Entries in both partitions consist of time series of arrays of thermodynamic fields (pressure, density, and temperature), kinematic fields (position and velocity), and additional fields that depend on thermodynamic and/or kinematic fields (e.g., material stress). Materials in the CYL partition include solids (aluminium, copper, depleted uranium, stainless steel, tantalum, and a generic polymer), a liquid (water), gases (air, nitrogen), and a generic detonating material (high explosive, HE). The PLI partition spans a highly varying geometry but consists of fixed materials across entries: Copper, aluminium, stainless steel, generic polymer, and generic HE. HEAT captures critical phenomena such as momentum transfer, shock propagation, plastic deformation, and thermal effects, making HEAT a valuable benchmark for development of AI/ML emulation of multi-material shock propagation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dispatch analysis of flexible power operation with multi-unit small modular reactors

The relevance of nuclear power plant flexible power operation (FPO) is rising due to increased penetration of variable renewables. Small Modular Reactors (SMRs) such as NuScale are envisaged to have multiple units on one site. This provides opportunities for coordinating FPO between units, but also introduces challenges as common services (e.g., refueling equipment) must be shared between units and fuel loadings may be standardized. It is therefore important to quantify whether FPO with SMRs is economically beneficial. Here, this paper quantifies the economics of FPO on different timescales to analyze refueling outages of such multi-unit SMRs. A 24-h price-taker profit maximization of SMR operation is first solved considering revenues from wholesale power and ancillary service (AS) markets. The results are then used to find the most profitable operation strategy over several years for multiple units at the site accounting for the physics and materials limits on NPP operation, including refueling outages. Results show a small but appreciable participation of nuclear into AS, contributing 4% of revenues. Due to the longer fuel cycle, over a decade of operation the units’ refueling outages drifted by 3 months, ultimately leading to refueling during summer.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

The simulation of 3D hypervelocity spallation using a hydrocode PAGOSA with FLIP+MPM

n this work, a hydrocode PAGOSA with FLIP+MPM is presented and exercised to investigate the fracture in ductile material. The merit of PAGOSA with FLIP+MPM to solve the advection problem is first illustrated by a solid piston periodically moving in a sealed tube with air. Furthermore, the ability of PAGOSA with FLIP+MPM to capture the fracture in material is shown by a simple stretching fracture in ductile material. In both of two benchmark problems, the PAGOSA results and analytical solutions are also presented for comparison. Then PAGOSA with FLIP+MPM is used to model complex spall in ductile material, which is a challenging problem in engineering applications. The convergences of PAGOSA with FLIP+MPM—based on both the mesh size and the marker density—are investigated by monitoring free surface velocity. To further show the ability of PAGOSA with FLIP+MPM to predict fracture in ductile material, the numerical results are compared with the experimental results and other published numerical results. Moreover, the effect of the spall parameter in PAGOSA with FLIP+MPM on the numerical simulation is also analyzed by investigating the free surface velocity. Finally, the effects of the peak compressive stress, the tensile strain rate, and the loading rate on the spallation are further investigated. The numerical results show that PAGOSA with FLIP+MPM can improve PAGOSA's performance when applied to predicting fracture in ductile materials, and is robust enough to accurately predict spall fracture in ductile materials.

42 ENGINEERING↗
DOE OSTI2023

The dynamics of suspensions of prolate spheroidal particles—Effects of volume fraction

Here, the effect of volume fraction on the dynamics of a gravity driven suspension of prolate spheroidal solid particles in a fully periodic domain is examined by fully resolved numerical simulations, for relatively modest Reynolds numbers (around 20, depending on the volume fraction). Three systems are examined, for volume fractions of 2.2% (20 particles), 5.5% (50 particles), and 9.9% (90 particles). The results show a transition from flow at low volume fractions dominated by the hydrodynamic interactions between the particles and the fluid, although modified by collisions, to flow dominated by collisions at higher volume fractions. The distribution of the particles with respect to each other, as measured by the probability distribution of nearest distances, is nearly random at low volume fractions but at high volume fractions, there is more clustering than for a random distribution. At lower volume fractions most of the particles fall broadside on whereas at the highest volume fraction, their orientation is essentially random.

42 ENGINEERING↗
DOE OSTI2020

Catalytic consequences of hydrogen addition events and solvent-adsorbate interactions during guaiacol-H 2 reactions at the H 2 O-Ru(0 0 0 1) interface

Catalytic reactions of biomass-derived phenolics and H 2 occur on transition metal surfaces via competitive C–O cleavage and ring saturation pathways, with both requiring multiple hydrogen addition events before forming their respective rate limiting transition states. These events are markedly affected by solvent chemical identity, with polar protic solvents ionizing hydrogen adatoms (H*) to interfacial protons (H + ) and opening up new catalytic routes. Here, we establish the reaction coordinate space for guaiacol-H 2 reactions on Ru(0 0 0 1) using density functional theory and describe the atomic-scale effect of a polar protic solvent, H 2 O. Coupled H + and H* attack leads to quasi-equilibrated enol and keto intermediates as the precursors for C–O cleavage and ring saturation, respectively. For C–O cleavage, H 2 O solvent enables a lower energy pathway via concomitant transfer of the hydroxyl H + to the methoxy oxygen during C–OCH 3 cleavage, forming a charge separated [Ru(s)–(C 6 H 5 O – )…(H + )…OCH 3 ] transition state and reducing the barrier by up to 0.8 eV as compared to unassisted C–OCH 3 cleavage. For ring saturation, H* attack onto an unsaturated meta carbon is rate limiting with no direct solvent participation, suggesting that protic polar solvents selectively promote the C–O cleavage pathway. Taken together, we show that activating guaiacol for either C–O bond cleavage or ring saturation product formation depends on the reactive hydrogen identity (H* or H + ), enol/keto isomerization equilibrium, and accessibility of the proton assisted Car–OCH 3 cleavage transition state. Here, all such factors are tunable via changes to the solvent or metal identity.

09 BIOMASS FUELS↗
DOE OSTI2020

Crystal and magnetic structure of polar oxide HoCrWO 6

In this study, polar magnetic oxide HoCrWO 6 is synthesized and its crystal structure, magnetic structure, and thermodynamic properties are investigated. HoCrWO 6 forms the polar crystal structure (space group Pna 2 1 (#33)) due to the cation ordering of W 6+ and Cr 3+ . There is an antiferromagnetic transition at T N = 24.5 K along with the magnetic entropy change (~5 J.Kg. –1 K –1 at 70 kOe). Neutron diffraction measurement indicates that both Cr and Ho sublattices are ordered with the moment of 2.32(5)μ B and 8.7(4)μ B at 2 K, respectively. While Cr forms A-type collinear antiferromagnetic (AFM) structure with magnetic moment along the b axis, Ho sublattice orders in a non-coplanar AFM arrangement. A comparison with isostructural DyFeWO 6 and DyCrWO 6 indicates that the magnetic structure of this family of compounds is controlled by the presence or absence of e g electrons in the transition metal sublattice.

42 ENGINEERING↗