SEARCH · Search NASA
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.
Nontargeted molecules of alpha-immunoconjugate (AIC) intravenously injected in clinical trials of targeted alpha therapy (TAT) could be transported by convection and diffusion along with blood or lymphatic circulation.
With ongoing advances in fusion and advanced fission reactors, quantifying irradiation effects in materials is critical. Transmutation-induced helium in cladding and structural materials can drive swelling and embrittlement, thereby reducing these components’ lifespans. Yet most studies ignore the considerable uncertainties in predicting helium generation rates. In this work, we created a code wrapper, F-SCATTER, that automatically performs simulations in FISPACT-II. We used this tool to investigate potential variance in helium generation rate, or He/dpa, calculations based on deviations in alloy composition, irradiating neutron flux spectrum, computational methodology, and nuclear data sources. We used 12 wt% Cr HT9 steel as the reference case and observed a 6.5%–98.3% He/dpa spread based on compositional variation within a single chemical specification, a 1.8%–11.5% He/dpa variation upon the incorporation of a 15% artificial uncertainty in flux at each energy, and a He/dpa difference as high as 231% when using ENDF/B-VIII.0 versus TENDL-2021 data libraries. Similar results were found for other prominent iron-based alloys, including Grade 91, castable nano-structured alloy, and 316H—where additional variations exist based on reactor type (e.g. thermal, fast, or fusion) and alloying elements such as carbon, nitrogen, and nickel. Based on the simulated results, we conclude that a significant part of the heat-to-heat variability in swelling responses of Fe-based alloys can be driven by impurity content in alloy compositions, and, therefore, chemical control should be a key element in supply chain design for advanced nuclear energy systems. Furthermore, we provide critical recommendations on best practices for evaluating and reporting helium production and lattice damage rates when computing predictions with multiphysics programs such as FISPACT-II.
A computational conjugate heat transfer methodology was developed, as a first step towards an efficient and accurate multiphysics, thermo-fluid computational methodology to predict environments for hypothetical solid-core, nuclear thermal engine thrust chamber and components. A solid conduction heat transfer procedure was implemented onto a pressure-based, multidimensional, finite-volume, turbulent, chemically reacting, thermally radiating, and unstructured grid computational fluid dynamics formulation. The conjugate heat transfer of a cylindrical material specimen heated by an impinging hot hydrogen jet inside an enclosed test fixture was simulated and analyzed. The solid conduction heat transfer procedure was anchored with a standard solid heat transfer code. Transient analyses were then performed with ,variable thermal conductivities representing three composites of a material utilized as flow element in a legacy engine test. It was found that material thermal conductivity strongly influences the transient heat conduction characteristics. In addition, it was observed that high thermal gradient occur inside the cylindrical specimen during an impulsive or a 10 s ramp start sequence, but not during steady-state operations.
An in-house code is developed to simulate the inductively coupled plasma (ICP). The model comprises the fluid, electromagnetic and transformer submodels. Fluid equations are solved to evaluate the plasma flow parameters, including the plasma and neutral densities, ion and neutral velocities, electron flux, electron temperature, and electric potential. The model relies on the ambipolar approximation and offers the evaluation of plasma parameters without solving the sheath region. The electromagnetic model handles the calculation of the electric and magnetic fields using the magnetic vector potential. The transformer model captures the effect of the matching circuit utilized in laboratory experiments for RF power deposition. The continuity and momentum equations are solved using finite volume method. The energy, electric potential, and magnetic vector potential equations are solved using finite difference method. The resulting linear systems of equations are solved with iterative solvers including Jacobi and GMRES. The code is written using the C++ programming language, it works in parallel and has graphical user interface. The model is applied to study ICP characteristics of a plasma confined within a cylindrical chamber with dielectric walls for two different power deposition cases. The results obtained from the developed model are verified using the plasma module of COMSOL Multiphysics. The model is also applied to a plasma source configuration, and it is demonstrated that there is an overall increase in the plasma potential when current is extracted from ICP with a biased wall electrode.
Computational fluid dynamics simulations have been performed to study the applicability of engineering CFD methods for thermals striping analysis. Thermal striping is the fluctuating temperature profile in a solid caused by fluctuating fluid temperature, and the resulting fluctuating thermal stresses cause high cycle fatigue and eventual material failure. Here this study presents the methodology for thermal striping analysis that including a transient conjugate heat transfer RANS model and benchmark the results against a sodium triple jet experiment. The multiphysics calculations are implemented to analyze the thermal stresses in the solid domain caused by the coupled heat transfer between the fluid and solid domains. Validation data from a liquid sodium triple jet experiment include time-averaged temperature measurements and power spectra of the temperature signal. The numerical results agree well with these experimental measures, demonstrating key features such as the dominant frequency of temperature fluctuations. The applicability of a low-cost wall treatment method is demonstrated, enabling key computational savings. Finally, the performance of two finite element stress analysis software packages is compared, and the validity of the lower-cost method is confirmed. These results demonstrate the applicability of engineering methods for computational thermal striping calculations, enabling thermal striping estimations in large fluid systems such as the core of a Sodium-cooled nuclear reactor.
The design and development process of a Nuclear Thermal Propulsion (NTP) system requires extensive multiphysics modeling to couple the neutron physics and thermal feedback effects to determine the reactor’s power shape. Propulsion system performance codes utilize this power shape to determine NTP key performance parameters. While the power shape is heavily dependent on the temperature profile and geometry of the reactor, many analyses either assume a constant power shape, or use neutronics analysis to determine a power shape for a specific reactor configuration. The development of a coupling interface for a propulsion system performance code and a Monte Carlo neutron transport code (OpenMC) allows for the reactor power shape to be calculated in an iteration loop. The interface utilizes a file share system to transfer geometry dimensions, temperatures, and material identifiers to OpenMC, which is used to perform a neutron transport simulation of a design like the government Testing Reference Design reactor. The interface is then able to post-process the results from OpenMC and use the same file share system to share a power shape and other important neutron transport parameters to the system performance code. Initial results show that neglecting the changes to power shape when comparing reactor configurations can yield inaccurate results. Furthermore, utilizing propellants other than hydrogen gas can cause significant changes to the power shape, and thus, the thermal performance of a specific reactor design. This methodology is being expanded to allow for multiple families of NTP reactors to be analyzed, including block moderator, particle bed, and NERVA-derived reactors.
The design and development process of a Nuclear Thermal Propulsion (NTP) system requires extensive multiphysics modeling to couple the neutron physics and thermal feedback effects to determine the reactor’s power shape. Propulsion system performance codes utilize this power shape to determine NTP key performance parameters. While the power shape is heavily dependent on the temperature profile and geometry of the reactor, many analyses either assume a constant power shape, or use neutronics analysis to determine a power shape for a specific reactor configuration. The development of a coupling interface for a propulsion system performance code and a Monte Carlo neutron transport code (OpenMC) allows for the reactor power shape to be calculated in an iteration loop. The interface utilizes a file share system to transfer geometry dimensions, temperatures, and material identifiers to OpenMC, which is used to perform a neutron transport simulation of a design like the government Testing Reference Design reactor. The interface is then able to post-process the results from OpenMC and use the same file share system to share a power shape and other important neutron transport parameters to the system performance code. Initial results show that neglecting the changes to power shape when comparing reactor configurations can yield inaccurate results. Furthermore, utilizing propellants other than hydrogen gas can cause significant changes to the power shape, and thus, the thermal performance of a specific reactor design. This methodology is being expanded to allow for multiple families of NTP reactors to be analyzed, including block moderator, particle bed, and NERVA-derived reactors.
RELAP-7 is a nuclear systems safety analysis code being developed at the Idaho National Laboratory, and is the next generation tool in the RELAP reactor safety/systems analysis application series. RELAP-7 development began in 2011 to support the Risk Informed Safety Margins Characterization (RISMC) Pathway of the Light Water Reactor Sustainability (LWRS) program. The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical methods, and physical models in order to provide capabilities needed for the RISMC methodology and to support nuclear power safety analysis. The code is being developed based on Idaho National Laboratory’s modern scientific software development framework – MOOSE (the Multi-Physics Object-Oriented Simulation Environment). The initial development goal of the RELAP-7 approach focused primarily on the development of an implicit algorithm capable of strong (nonlinear) coupling of the dependent hydrodynamic variables contained in the 1-D/2-D flow models with the various 0-D system reactor components that compose various boiling water reactor (BWR) and pressurized water reactor nuclear power plants (NPPs). During Fiscal Year (FY) 2015, the RELAP-7 code has been further improved with expanded capability to support boiling water reactor (BWR) and pressurized water reactor NPPs analysis. The accumulator model has been developed. The code has also been coupled with other MOOSE-based applications such as neutronics code RattleSnake and fuel performance code BISON to perform multiphysics analysis. A major design requirement for the implicit algorithm in RELAP-7 is that it is capable of second-order discretization accuracy in both space and time, which eliminates the traditional first-order approximation errors. The second-order temporal is achieved by a second-order backward temporal difference, and the one-dimensional second-order accurate spatial discretization is achieved with the Galerkin approximation of Lagrange finite elements. During FY-2015, we have done numerical verification work to verify that the RELAP-7 code indeed achieves 2nd-order accuracy in both time and space for single phase models at the system level.
The Transient Reactor Test Facility (TREAT) at Idaho National Laboratory (INL) is currently preparing for resumption of operation after more than 20 years in a long-term standby state. In parallel with preparation for operations and a new transient testing program, the DOENE Advanced Modeling and Simulation program (NEAMS) has been supporting the development of three-dimensional transient multi-physics modeling capabilities within the reactor physics application MAMMOTH. Recently, the MAMMOTH reactor physics team has begun development of uncertainty and sensitivity analysis for the TREAT core. This research will be used to support the analysis of reactor operation and experiment planning. Traditional approaches, i.e., Monte Carlo simulations, for such tasks require a large number of model evaluations under both steady state and transient conditions. However, exhaustive sets of sampling runs are simply not practical for TREAT since evaluations of each model can take a substantial amount of time. A solution to this problem is to build a surrogate model employing a smaller number of model evaluations. The surrogate model is designed to provide an accurate and fast approximation to the detailed TREAT model. In this manner, the surrogate model can be used in place of the TREAT transient model for uncertainty and sensitivity analysis at a fraction of the computational cost. The high dimensional model representation (HDMR) method is employed here for approximating the responses of interest in terms of functions of lower dimensions. This method has been implemented in RAVEN and all TREAT transient calculations are performed using the multiphysics code MAMMOTH. In addition, a code interface has also been developed in order to couple RAVEN with MAMMOTH. The surrogate model is then tested and compared for TREAT startup testing program Transient Test no. 15 core, and its accuracy is measured by its ability to estimate the mean, coefficient of variance, skewness, excess kurtosis and probability density function of the outputs of Transient Test no. 15 model, as well as root mean square error metric.
Here, a computational framework for assisting in the development of novel textiles is presented. Electronic textiles are key in the rapidly growing field of wearable electronics for both consumer and military uses. There are two main challenges to the modeling of electronic textiles: the discretization of the textile microstructure and the interaction between electromagnetic and mechanical fields. A director-based beam formulation with an assumed electrical current is used to discretize the fabric at the level of individual fibrils. The open-source package FEniCS was used to implement the finite element model. Contact integrals were added into the FEniCS framework so that multiphysics contact laws can be incorporated in the same framework, leveraging the code generation and automated differentiation capabilities of FEniCS to produce the tangents needed by the implicit solution method. The computational model is used to construct and determine the mechanical, thermal, and electrical properties of a representative volume elements of a plain woven textile. Dynamic relaxation to solve the mechanical fields and the electrical and thermal fields is solved statically for a given mechanical state. The simulated electrical responses are fit to a simplified Kirchhoff network model to determine effective resistances of the textile.
U.S. industry is expected to deploy tristructural isotropic (TRISO) particle fuel technologies for commercial reactors within the next decade. In our previous work, we defined a preliminary transient design space for TRISO fuels, identified potential gaps in the available data, and began to develop multiphysics modeling tools that could be applied to design targeted Transient Reactor Test Facility (TREAT) experiments to fill these gaps. This work builds on that foundation by (1) updating BISON fuel performance and Griffin reactor physics models to reflect the current TREAT experiment tube and capsule designs,(2) coupling the codes to improve the accuracy and usability of the transient design analyses, and (3) demonstrating their use over an expanded design space that includes fuel burnup. The simulated mechanical responses of the TRISO particles were complex functions of fission product accumulation, fission gas release, and irradiation-induced dimensional change in the pyrolytic carbon layers. Predicted tangential stresses in particle silicon carbide layers were least compressive for preheated tests involving fresh fuels but remained compressive throughout the ranges of temperature, heat rate, and burnup considered in this work. Finally, comparisons between the potential TREAT transients and historical test reactor irradiations showed that the TREAT tests would produce significantly lower average energy deposition rates, yielding less severe transients with greater relevance to near-term commercial applications. Use of these predictive capabilities has the potential to increase the value of each test, improving the overall efficiency and cost effectiveness of transient testing for TRISO and other advanced fuels.
Fusion devices are anticipated to produce and consume several kilograms of tritium per year. This rare fuel resource is both highly mobile and radioactive, making tracking inventories a priority for operation and safety. The fusion safety program at the Idaho National Laboratory has been developing the Tritium Migration and Analysis Program (TMAP), of which the latest version is a MOOSE-based application. TMAP8 is verified against its predecessors and possesses additional multi-dimensional tritium migration modeling capabilities. As we extend its capabilities towards both whole device (in multiple dimensions) and whole plant (with multiple components) simulations, the syntax of inputs must become compact, descriptive, compatible with quality assurance processes, and as error-proof as achievable. The new Physics system developed MOOSE can set up equations and instantiating them on plant components. The system permits the automatic definition of complex discretization with a consistency between object parameters achieved programmatically. The Physics system can currently instantiate the equations for heat conduction and Navier Stokes weakly compressible flow. In MOOSE-terms, it automates the definition of kernels, boundary conditions, and several core and helper materials and fields. As part of this effort, Physics classes were developed for tritium migration, trapping and advection within either a multi-dimensional Navier Stokes fluid dynamics simulation, or a 1D thermal hydraulics piping system. In this presentation, we will showcase the new syntax, its application to several verification and validation cases which were already studied using the classical TMAP8 syntax, and a demonstration of the new coupling capabilities for the migration of tritium into blanket coolant channels and the subsequent advection into the coolant loop.
The U.S. nuclear industry is expected to deploy tristructural isotropic (TRISO) particle fuel technologies for commercial reactors within the next decade. In previous work, we defined a preliminary transient design space for TRISO fuels, identified potential gaps in the available data, and began to develop multiphysics modeling tools that could be applied to design targeted Transient Reactor Test Facility (TREAT) experiments to fill these gaps. Here, this work builds on that foundation by (1) updating BISON fuel performance and Griffin reactor physics models to reflect the current TREAT experiment tube and capsule designs, (2) coupling the codes to improve the accuracy and usability of the transient design analyses, and (3) demonstrating their use over an expanded design space that includes fuel burnup. The simulated mechanical responses of the TRISO particles were complex functions of fission product accumulation, fission gas release, and irradiation-induced dimensional change in the pyrolytic carbon layers. The predicted tangential stresses in the particles' silicon carbide layers were least compressive for preheated tests involving fresh fuels but remained compressive throughout the ranges of temperature, heat rate, and burnup considered in this work. Finally, comparisons between the potential TREAT transients and historical test reactor irradiations showed that the TREAT tests would produce significantly lower average energy deposition rates, yielding less severe transients with greater relevance to near-term commercial applications. The use of these predictive capabilities has the potential to increase the value of each test, improving the overall efficiency and cost-effectiveness of transient testing for TRISO and other advanced fuels.
The design of a power electronics module is a multiphysics problem and involves electrical, thermal, and reliability targets and constraints. Although the electrical and thermal design consider multiple aspects such as materials, geometry, and energy losses, reliability design is often limited to the selection of the attachment materials and electrical interconnect types. In this paper, we incorporate reliability metrics in the design phase by investigating the impact of package geometry on the thermomechanical behavior of the die-attachment material. To this end, we conducted thermal and power cycling simulations of a commercial six-pack power module with silicon carbide devices to compute junction temperature and strain energy density per cycle, respectively. We performed multiple simulations with different geometric dimensions and established a correlation between input features and output variables using subspace-based dimension reduction. The machine learning-based dimension reduction method serves as a surrogate model, which can be employed to identify the optimal module design from a thermal and reliability standpoint.
Welding processes used in the production of pressure vessels impart residual stresses in the manufactured component. Computational modeling is critical to predicting these residual stress fields and understanding how they interact with notches and flaws to impact pressure vessel durability. Here, in this work, we present a finite element model for a resistance forge weld and validate it using laboratory measurements. Extensive microstructural changes, near-melt temperatures, and large localized deformations along the weld interface pose significant challenges to Lagrangian finite element modeling. The proposed modeling approach overcomes these roadblocks in order to provide a high-fidelity simulation that can predict the residual stress state in the manufactured pressure vessel; a rich microstructural constitutive model accounts for material recrystallization dynamics, a frictional-to-tied contact model is coordinated with the constitutive model to represent interfacial bonding, and adaptive remeshing is employed to alleviate severe mesh distortion. An interrupted-weld approach is applied to the simulation to facilitate comparison to displacement measures. Several techniques are employed for residual stress measurement in order to validate the finite element model: neutron diffraction, the contour method, and the slitting method. Model-measurement comparisons are supplemented with detailed simulations that reflect the configurations of the residual-stress measurement processes themselves. The model results show general agreement with experimental measurements, and we observe some similarities in the features around the weld region. Factors that contribute to model-measurement differences are identified. Finally, we conclude with some discussion of the model development and residual stress measurement strategies, including how to best leverage the efforts put forth here for other weld problems.
Core radial expansion in liquid-metal cooled fast reactor systems is a well-known phenomenon that produces strong reactivity feedback effects. It is crucial from a safety standpoint to design a reactor with the proper core restraint system to guide expansion of the fuel assemblies into a formation that produces negative reactivity feedback in accident conditions. However, the modeling of core radial expansion and its associated reactivity feedback is extremely challenging. Liquid metal-cooled fast reactor designers (both industry and commercial) have pointed to improved modeling of core radial expansion as a key modeling challenge for these types of reactors. The Department of Energy Nuclear Energy Advanced Modeling and Simulation (DOE-NEAMS) program has an important role to play in helping to develop a robust, predictive tool for this phenomenon. A literature review has been performed to identify current simulation capabilities applicable to the prediction of core radial expansion in liquid metal-cooled fast reactors and the resultant reactivity feedback effects. An overview as well as a detailed description of the multi-physics phenomena underlying core radial expansion is given in the first two chapters. The objective of this report is to identify existing software that has been or could be applied to the core radial expansion problem, and then to identify possible paths to develop a robust, coupled Multiphysics tool for advanced analysis, in which all three physics (radiation transport, structural mechanics, thermal fluids) are taken into account considering feedback effects from each physics to the other. Several computer codes and workflows have been developed in the United States and abroad to perform individual aspects of radial expansion calculations or to perform loosely coupled analysis. No code system currently exists which tightly and robustly couples the neutronics, thermal hydraulics, and thermal mechanical physics with enough detail to fully resolve the complex core radial expansion reactivity feedback effects. The future availability of such a code system is of vital importance to fully understanding the reactivity feedback effects that occur due to radial expansion, and consequently to optimizing the design of the core restraint system. A potential code development path forward using MOOSE-based tools is suggested in order to leverage the natural tight coupling and robustness of MOOSE-based applications. Additionally, potential gaps for modeling and common assumptions are listed which require further investigation to ensure accuracy for prediction of deformations.
Cardinal is a MOOSE application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics to the MOOSE framework, closing the neutronics and thermal-fluid gaps in conducting tightly-coupled, high-resolution multiscale and multiphysics analyses. By leveraging MOOSE's interfaces for wrapping external codes, Cardinal overcomes many challenges encountered in earlier multiphysics coupling works, such as file-based I/O or overly-restrictive geometry mapping requirements. In this work, we leverage a subset of the multiphysics interfaces in Cardinal to perform coupling of OpenMC neutron transport, MOOSE heat conduction, and THM thermal-fluids for steady-state modeling of a prismatic gas reactor fuel assembly.
This study presents a novel multiphysics phase-field fracture model to analyze high-burnup uranium dioxide (UO2) fuel behavior under transient reactor conditions. Fracture is treated as a stochastic phase transition, which inherently accounts for the random microstructural effects that lead to variations in the value of fracture strength. Moreover, the model takes into consideration the effects of temperature and burnup on thermal conductivity. Therefore, the model is able to predict crack initiation, propagation, and complex morphologies in response to thermal gradients and stress distributions. Several simulations were conducted to investigate the effects of operational and transient conditions on fracture behavior and the resulting cracking patterns. High-burnup fuels exhibit reduced thermal conductivity, elevating temperature gradients and resulting in extensive radial and circumferential cracks. Transient heating rates and temperatures significantly affect fracture patterns, with higher heating rates generating steeper gradients and more irregular crack trajectories. This approach provides critical insights into fuel integrity during accident scenarios and supports the safety evaluation of extended burnup limits.