Search NASA⌕ Search

SEARCH · Search NASA

Results for “Advanced Modeling and Simulation”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Mesoscale Modeling for Restructuring and Fragmentation in High Burnup UO 2

This report summarizes the mesoscale modeling conducted in fiscal year 2025 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, focusing on the microstructural evolution and restructuring in high burnup UO 2 nuclear fuel and its impact on fuel fragmentation. We developed a pioneering phase-field model to simulate restructuring behavior across different regions of high burnup fuel, including the dark zone and rim region. A grand-potential-based phase-field model is employed to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of temperature and burnup rate were studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to initiate around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. The restructuring occurs faster with higher burnup rate. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. The predicted microstructures are consistent with experimental observations of the restructured regions. Based on the mesoscale simulations, a mechanistic model for restructuring and grain size evolution was implemented in BISON. Thus, this work provides a first-of-its-kind restructuring model for different regions of high-burnup fuel to BISON. It was shown that the model predicts appropriate grain size evolution along fuel radius as those observed in experiments. This model enables BISON to predict the effect of restructuring on the fission gas release. Finally, the phase-field fracture simulation with dark zone specific microstructures were presented to provide the fragmentation criteria for the dark zone. This work introduces a first-of-its-kind restructuring model for high burnup fuel, significantly enhancing BISON's predictive capabilities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Effect of Surrogate Fission Products on Thermal Conductivity and Viscosity of NaCl-UCl 3

Thermophysical properties of salt systems relevant to molten salt reactors (MSRs) are experimentally measured in support of the US Department of Energy Office of Nuclear Energy (DOE-NE) MSR campaign within the Advanced Reactor Technology (ART) Program. These property measurements also support the development of the thermophysical arm of the Molten Salt Database (MSD) within DOE-NE’s Nu clear Energy Advanced Modeling and Simulation (NEAMS) Program. Several US Department of Energy (DOE) laboratories, including Oak Ridge National Laboratory (ORNL), have been capitalizing on modern methodologies and sample characterization techniques to conduct thermophysical property measurements to support these programs and to provide MSR developers and modelers with more recent and higher-quality data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Lattice Design and Advanced Modeling to Guide the Design of High-Performance Lightweight Structural Materials

Lightweight structural materials are required to increase the mobility of fission batteries. The materials must feature a robust combination of mechanical properties to demonstrate structural resilience. The primary objective of this project is to produce lightweight structural materials whose strength-to-weight ratios exceed those of the current widely used structural materials such as 316L stainless steels (316L SS). To achieve this, advanced modeling and simulation tools were employed to design lattice structures with different lattice parameters and different lattice types. A process was successfully developed for transforming lattice-structures models into Multiphysics Object Oriented Simulation Environment (MOOSE) inputs. Finite element modeling (FEM) was used to simulate the uniaxial tensile testing of the lattice-structured parts to investigate the stress distribution at a given displacement. The preliminary results showed that the lattice-structured sample displayed a lower Young’s modulus in comparison with the solid material and that the unit cell size of the lattice had a minimal effect. The novelty here is to apply up-front modeling to determine the best structure for the application before actually producing the sample. The approach of using modeling as a guiding tool for preliminary material design can significantly save time and cost for material development.

36 MATERIALS SCIENCE↗
NASA NTRS2016

Development and Use of Engineering Standards for Computational Fluid Dynamics for Complex Aerospace Systems

Computational fluid dynamics (CFD) and other advanced modeling and simulation (M&S) methods are increasingly relied on for predictive performance, reliability and safety of engineering systems. Analysts, designers, decision makers, and project managers, who must depend on simulation, need practical techniques and methods for assessing simulation credibility. The AIAA Guide for Verification and Validation of Computational Fluid Dynamics Simulations (AIAA G-077-1998 (2002)), originally published in 1998, was the first engineering standards document available to the engineering community for verification and validation (V&V) of simulations. Much progress has been made in these areas since 1998. The AIAA Committee on Standards for CFD is currently updating this Guide to incorporate in it the important developments that have taken place in V&V concepts, methods, and practices, particularly with regard to the broader context of predictive capability and uncertainty quantification (UQ) methods and approaches. This paper will provide an overview of the changes and extensions currently underway to update the AIAA Guide. Specifically, a framework for predictive capability will be described for incorporating a wide range of error and uncertainty sources identified during the modeling, verification, and validation processes, with the goal of estimating the total prediction uncertainty of the simulation. The Guide's goal is to provide a foundation for understanding and addressing major issues and concepts in predictive CFD. However, this Guide will not recommend specific approaches in these areas as the field is rapidly evolving. It is hoped that the guidelines provided in this paper, and explained in more detail in the Guide, will aid in the research, development, and use of CFD in engineering decision-making.

Lee, Hyung B.↗
NASA NTRS2019

Sensor Analysis, Modeling, and Test for Robust Propulsion System Autonomy

An approach is presented supporting analysis, modeling, and test validation of operational flight instrumentation (OFI) that facilitates critical functions for the Space Launch System (SLS) main propulsion system (MPS). Certain types of OFI sensors were shown to exhibit highly nonlinear and non-gaussian noise characteristics during acceptance testing, motivating the development of advanced modeling and simulation (M&S) capability to support algorithm verification and flight certification. Hardware model and algorithm simulation fidelity was informed by a risk scoring metric; redesign of high-risk algorithms using test-validated sensor models significantly improved their expected performance as evaluated using Monte Carlo acceptance sampling methods. Autonomous functions include closed-loop ullage pressure regulation, pressurant leak detection, and fault isolation for automated safing and crew caution and warning (C&W).

Orr, Jeb S.↗
DOE OSTI2018

FEASIBILITY ASSESSMENT FOR DEVELOPING AN INTEGRAL LOCA TESTING CAPABILITY AT THE TRANSIENT RESEARCH TEST (TREAT) REACTOR

The study of Light Water Reactor (LWR) fuel performance in Loss of Coolant Accidents (LOCAs) has been a significant and evolving field since the beginning of the commercialization of LWR technology. The following paper presents the results of feasibility assessment to develop an in-pile, integral LOCA testing capability using the recently restarted TREAT reactor located at the Idaho National Laboratory (INL). Critical performance objectives for the testing platform are derived from a review of known LWR fuel degradation phenomena in LOCA scenarios. Advanced modeling and simulation tools were then used to gain additional insight into both the system performance and fuel performance characteristics of LOCA testing at TREAT. Results show that TREAT is well suited to conduct in-pile integral LOCA testing.

42 ENGINEERING↗
DOE OSTI2019

The Future Landscape of High-Redshift Galaxy Cluster Science

Modern galaxy cluster science is a multi-wavelength endeavor with cornerstones provided by X-ray, optical/IR, mm, and radio measurements. In combination, these observations enable the construction of large, clean, complete cluster catalogs, and provide precise redshifts and robust mass calibration. The complementary nature of these multi-wavelength data dramatically reduces the impact of systematic effects that limit the utility of measurements made in any single waveband. The future of multi-wavelength cluster science is compelling, with cluster catalogs set to expand by orders of magnitude in size, and extend, for the first time, into the high-redshift regime where massive, virialized structures first formed. Unlocking astrophysical and cosmological insight from the coming catalogs will require new observing facilities that combine high spatial and spectral resolution with large collecting areas, as well as concurrent advances in simulation modeling campaigns. Together, future multi-wavelength observations will resolve the thermodynamic structure in and around the first groups and clusters, distinguishing the signals from active and star-forming galaxies, and unveiling the interrelated stories of galaxy evolution and structure formation during the epoch of peak cosmic activity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Surrogate Models for TREAT Transient Calculations

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 DOE-NE 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 [1]. 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 here 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) [2] method is employed here for approximating the responses of interest in terms of functions of lower dimensions. This method has been implemented in RAVEN [3] 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 the MAMMOTH model of the TREAT startup testing program Transient Test #15 core [4], 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 #15 model, as well as root mean square error metric.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Refinements of a treat model for transient analysis

In order to prepare for the TREAT reactor being taken off standby, the Gateway for Advanced Modeling and Simulation Lab has used a set of KENO models to simulate certain TREAT temperature- limited transients using the Transient Reactor excursion simulator (T-ReX) code, a hybrid deterministic-stochastic code. These transients are from the M8CAL experiment series, a set of experiments done to calibrate the reactor for the planned M8 series of fuel tests. Recently, the model has been updated in media definitions and geometry. Models were examined and calculations were performed to low uncertainty values to find the k{sub eff} and reactivity insertions. Alterations to modeled values of shutdown and initial transient rod insertion depth were made to create models which had experimental values of reactivity insertion, pre-transient k{sub eff}, or both. It was found that two of the altered media inputs, fuel and Zircaloy-3 cladding, had significant effect on the k{sub eff}. Surprisingly, the Zircaloy-3 composition had a larger effect on k{sub eff}. The reactivity insertions in this model's 2855, 2856 and 2857 temperature-limited transients were found to be above experimental values by a factor of 5.5% to 7%, and that pre-transient k{sub eff} was too high, at 1.002. It was found that increasing shutdown rod insertion by 2 to 2.5 cm decreased pre-transient k{sub eff} to 1.0, and that increasing shutdown rod insertion by 3 to 5 cm and decreasing initial transient rod insertion by 1 to 2 cm gave perfect pre-transient k{sub eff} and total reactivity insertion. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
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↗
DOE OSTI2016

Orientation-Dependent Surface Faceting in UO{sub 2}

The surface behavior of UO{sub 2} plays an important role throughout the lifetime of this widely-used commercial nuclear fuel, from manufacturing to waste storage. UO{sub 2} surface properties affect fuel swelling behavior by determining the morphology, size and distribution of fission gas bubbles. Additionally, UO{sub 2} surface characteristics influence fracture behavior, which is controlled by the thermodynamics of fracture surfaces. Therefore, detailed understanding of the fundamental surface behavior will help advance prediction of the fuel performance of UO{sub 2}. Due to significant anisotropy in surface tension, faceting features form after annealing in both UO{sub 2} and many other crystalline materials. Knowledge of the correlation between crystallographic orientation and surface faceting features will expand the fundamental knowledge of UO{sub 2} surface characteristics, leading to a better understanding of fission gas swelling and fracture mechanisms. Inherent characterization technique limited previous studies on the crystallographic orientation of UO{sub 2} surface facets to low Miller index surfaces in single crystal specimens. In this study, synchrotron Laue microdiffraction was utilized to examine the crystallographic orientation of grains on the coarse surface of a polycrystalline UO{sub 2} specimen. The results were then compared with the surface morphology, obtained from scanning electron microscopy (SEM), so that a general correlation between the surface faceting features and crystallographic orientation in UO{sub 2} can be established. Both the fundamental surface characteristics and the crystallographic orientation measurements obtained in this study help advance the experiment validation of MARMOT, a mesoscale multi-physics code that predicts the coevolution of microstructure and properties within reactor fuel during its lifetime in the reactor. The MARMOT code is an indispensable component of the MOOSE-BISON-MARMOT (MBM) code suite that has been developed by Idaho National Laboratory to enable next generation fuel performance modeling capability as a part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program Fuels Product Line (FPL). In this study, the surface-faceting features of UO{sub 2} were found to be dominated by the triple-plane structure and its variants. With this knowledge, by measuring the directions of the <110> edges, the crystallographic orientation of a surface UO{sub 2} grain can be precisely deduced in most regions of the standard stereogram triangle (SST) of the surface normal. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Comparison of Surrogate Models with Physical Models for Dynamic Probabilistic Risk Analysis Using the RAVEN Code

Dynamic probabilistic risk analysis (DPRA) has emerged in the nuclear industry as an effective tool to identify points of weakness in nuclear reactors and create a comprehensive approach to quantifying risk. While DPRA allows for detailed analysis, it usually requires thousands of runs due to the necessity of exploring low probability, high consequence events, which can be computationally prohibitively difficult to assess. Risk Analysis Virtual Environment (RAVEN) is a software package being developed under the Nuclear Energy Advanced Modeling and Simulation program at the Idaho National Laboratory. It is able to manage complex control logic to drive a RELAP-7 simulation and to handle sampling strategies including Monte Carlo, Latin Hypercube, Grid and Adaptive samplers. By performing a thorough investigation of the possibility space of an accident scenario with RELAP-7, RAVEN is able to identify consequences of accident scenarios that have not been encountered. The most computational resource intensive step of a RAVEN analysis is running the RELAP-7 code modeling the physical system behavior. Using surrogate models to represent the outcomes of more complex simulations rather than the full model can greatly reduce the computation effort required. In situations where the uncertainty distributions of a model are unknown, or when there are multiple sets of uncertainty distributions, DPRA generally requires rerunning the physical model for each set of uncertainty distributions used. By taking an unbiased sample of the physical model and creating a surrogate model from the results, the surrogate model can be sampled under a different set of uncertainty distributions to provide equivalent results. This paper presents an example of this approach. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Verification of High-Fidelity Neutronics Code PROTEUS for C5G7 Benchmark Problems

Under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the 3D whole core transport code PROTEUS has been developed at Argonne National Laboratory to perform the detailed 3D transport calculation for various reactor types including irregular or complex geometric reactors. PROTEUS includes two distinct transport solvers: SN and MOC. Among them, the MOC solver based on the extruded geometry, thus named MOCEX, can produce the faithful 3-D transport solutions by employing the rigorous and robust 2D/3D MOC formulation. This formulation overcomes the drawbacks of a full 3D MOC method requiring enormous memory and computational resources and resolves the instability issues of the 2D/1D method as well. The solution of MOCEX converges to that of a Monte Carlo code as the discretization of angle, space and energy is refined, unlike the 2D/1D method. As a verification and validation (V and V) effort of PROTEUS-MOCEX, this paper presents verification test results using the well-known 3D C5G7 benchmark problems by solving and comparing eigenvalues and pin powers with MCNP solutions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Surrogate Models for TREAT Transient Calculations

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.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Implementation of CMFD Acceleration Scheme in PROTEUS-MOC

The PROTEUS-MOC code is a three-dimensional (3D) transport code, which is being developed at Argonne National Laboratory under the DOE-NE's Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. It provides a fully consistent 3D discretization of the transport equation by combining the two-dimensional (2D) method of characteristics (MOC) for the planar direction with the discontinuous Galerkin method for the treatment of the axial variable. PROTEUS-MOC uses the unstructured finite element method (FEM) and thus allows the explicit representation of any extruded geometries. It also has a capability of massively parallel computation in energy group, angle, plane, and space using a top-level GMRES linear algebra solver. For parallel computation in space, it uses an unstructured domain decomposition by adapting METIS. The computational burden for direct heterogeneous core calculations based on unstructured geometries is excessive for practical reactor problems, and hence it is essential to implement an efficient acceleration scheme in PROTEUS-MOC. Recently, the coarse mesh finite difference (CMFD) method, which was originally applied to the nodal diffusion methods, has been successfully applied to the MOC transport solvers, resulting in a remarkable reduction in the number of fission source iterations. In addition, a modified CMFD method has been developed to overcome the limitation of the conventional CMFD that is applicable only to neutronics solvers with a conservative discretization. This CMFD method was derived to be consistent with the FEM, which is based on a nonconservative discretization. In this method, the neutron balance in each coarse mesh is enforced by introducing a pseudo absorption cross section (PAXS), and the well-established alternating solution process of CMFD and transport calculations is employed to accelerate source convergence. This modified CMFD acceleration scheme was adopted in PROTEUS-MOC to accelerate the fission source iteration. Pin-cell and assembly based CMFD schemes have been implemented by defining individual pin cells or assemblies as a coarse mesh. For this, an algorithm to map the finite elements to coarse meshes was developed, and the routines were developed to determine the coupling coefficients and PAXS of the CMFD equations from the transport solution. Preliminary performance tests were carried out using the C5G7 benchmark problem. This paper presents the main features of this CMFD acceleration scheme and the preliminary test results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗