Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled 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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Overview of New Tools to Perform Safety Analysis: BWR Station Black Out Test Case

Dynamic Probabilistic Risk Assessment (DPRA) methodologies couple system simulator codes (e.g., RELAP, MELCOR) with simulation controller codes (e.g., RAVEN, ADAPT). While system simulator codes accurately model system dynamics deterministically, simulation controller codes introduce both deterministic (e.g., system control logic, operating procedures) and stochastic (e.g., component failures, parameter uncertainties) elements into the simulation. Typically, a DPRA is performed by: 1) sampling values of a set of parameters from the uncertainty space of interest (using the simulation controller codes), and 2) simulating the system behavior for that specific set of parameter values (using the system simulator codes). For complex systems, one of the major challenges in using DPRA methodologies is to analyze the large amount of information (i.e., large number of scenarios ) generated, where clustering techniques are typically employed to allow users to better organize and interpret the data. In this paper, we focus on the analysis of a nuclear simulation dataset that is part of the Risk Informed Safety Margin Characterization (RISMC) Boiling Water Reactor (BWR) station blackout (SBO) case study. We apply a software tool that provides the domain experts with an interactive analysis and visualization environment for understanding the structures of such high-dimensional nuclear simulation datasets. Our tool encodes traditional and topology-based clustering techniques, where the latter partitions the data points into clusters based on their uniform gradient flow behavior. We demonstrate through our case study that both types of clustering techniques complement each other in bringing enhanced structural understanding of the data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Development of Transient Coupled Chemical-Thermal-Fluids Simulation for Unsealed Aluminum-Clad Spent Nuclear Fuel Storage Canister

This study examines a preliminary three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model of the unsealed canisters containing aluminum-clad spent nuclear fuel from INL’s ATR in the irradiated fuel storage facility. This 3D CFD model considers fully resolved fuel plates in an unsealed canister with small mass exchange with surrounding air as time evolves in the simulation and contains a network of chemical reactions to consider due to several of the side reactions occurring in the presence of water vapor, nitrogen, oxygen and carbon dioxide rather than in helium backfilled DOE standard canisters. Prior to the CFD modeling, the full set of chemical reactions under these conditions was analyzed, reduced, and refit to a more manageable number of equations for the computational efficiency of the 3D CFD simulation. To deal with the long time-scale of interest for tracking the possible buildup of trace species such as nitric acid, NOx, and hydrogen, a segregated time stepping algorithm is used, wherein the velocity field is only updated periodically, and assume to be in a pseudo steady state. This study first examines the theory of the model, workflow to establish radiolysis reaction network, and initial simulations of the evolutions of thermal fields, hydrogen gas and nitric acid concentrations within unsealed canisters in conditions relevant to the INL INTEC CPP-603 facility over a long period of time.

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

Development of Transient Coupled Chemical-Thermal-Fluids Simulation for Unsealed Aluminum-Clad Spent Nuclear Fuel Storage Canister

This study examines a preliminary three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model of the unsealed canisters containing aluminum-clad spent nuclear fuel from INL’s ATR in the irradiated fuel storage facility. This 3D CFD model considers fully resolved fuel plates in an unsealed canister with small mass exchange with surrounding air as time evolves in the simulation and contains a network of chemical reactions to consider due to several of the side reactions occurring in the presence of water vapor, nitrogen, oxygen and carbon dioxide rather than in helium backfilled DOE standard canisters. Prior to the CFD modeling, the full set of chemical reactions under these conditions was analyzed, reduced, and refit to a more manageable number of equations for the computational efficiency of the 3D CFD simulation. To deal with the long time-scale of interest for tracking the possible buildup of trace species such as nitric acid, NOx, and hydrogen, a segregated time stepping algorithm is used, wherein the velocity field is only updated periodically, and assume to be in a pseudo steady state. This study first examines the theory of the model, workflow to establish radiolysis reaction network, and initial simulations of the evolutions of thermal fields, hydrogen gas and nitric acid concentrations within unsealed canisters in conditions relevant to the INL INTEC CPP-603 facility over a long period of time.

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

Modeling FLEX Human Actions Using the EMRALD Dynamic Risk Assessment Tool

Most current efforts have modified and applied existing human reliability analysis (HRA) methods to treat human actions related to beyond-design-basis external events in which diverse and flexible coping strategies (FLEX) equipment would be deployed. However, many questions remain regarding the suitability of legacy HRA methods to address FLEX human actions, sparking the need for a new method of reasonably evaluating them. In this context, Idaho National Laboratory (INL) researched a relatively new approach to treating FLEX human actions via Event Modeling Risk Assessment Using Linked Diagrams (EMRALD) software. EMRALD was developed to support the increasing need for dynamic probabilistic risk assessment (PRA) models that can respond to evolving plant conditions during simulations. A couple benefits were identified when analyzing FLEX human actions through this software. In general, it is especially useful for evaluating a strategy’s feasibility, including FLEX human actions that require a relatively long time to perform. In this paper, we suggest how FLEX human actions can be modeled using the EMRALD software. Two different HRA modeling approaches using EMRALD are introduced: (1) procedure-based modeling and (2) PRA/HRA-based modeling. The former approach was introduced in the authors’ previous paper, whereas this paper mainly discusses how the latter approach works for an extended loss of AC power (ELAP) scenario with relevant procedures and PRA models. A hybrid method combining the two modeling approaches is introduced at the paper’s conclusion.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2019

Development of Transient Coupled Chemical-Thermal-Fluids Simulation for Unsealed Aluminum-Clad Spent Nuclear Fuel Storage Canister - 19273

This study examines a preliminary three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model of the unsealed canisters containing aluminum-clad spent nuclear fuel from INL's ATR in the irradiated fuel storage facility. This 3D CFD model considers fully resolved fuel plates in an unsealed canister with small mass exchange with surrounding air as time evolves in the simulation and contains a network of chemical reactions to consider due to several of the side reactions occurring in the presence of water vapor, nitrogen, oxygen and carbon dioxide rather than in helium backfilled DOE standard canisters. Prior to the CFD modeling, the full set of chemical reactions under these conditions was analyzed, reduced, and refit to a more manageable number of equations for the computational efficiency of the 3D CFD simulation. To deal with the long time-scale of interest for tracking the possible buildup of trace species such as nitric acid, NOx, and hydrogen, a segregated time stepping algorithm is used, wherein the velocity field is only updated periodically, and assumed to be in a pseudo steady state. This study first examines the theory of the model, workflow to establish radiolysis reaction network, and initial simulations of the evolutions of thermal fields, hydrogen gas and nitric acid concentrations within unsealed canisters in conditions relevant to the INL INTEC CPP-603 facility over a long period of time. (authors)

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

Connect microstructure evolution to property degradation with validated simulation

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500$^o$C to 1500$^o$C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 $\frac{D}{m^2s}$. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Realistic simulations of spin squeezing and cooperative coupling effects in large ensembles of interacting two-level systems

We describe an efficient numerical method for simulating the dynamics of interacting spin ensembles in the presence of dephasing and decay. The method builds on the discrete truncated Wigner approximation for isolated systems, which combines the mean-field dynamics of a spin ensemble with a Monte Carlo sampling of discrete initial spin values to account for quantum correlations. Here we show how this approach can be generalized for dissipative spin systems by replacing the deterministic mean-field evolution by a stochastic process, which describes the decay of coherences and populations while preserving the length of each spin. We demonstrate the application of this technique for simulating nonclassical spin-squeezing effects or the dynamics and steady states of cavity QED models with 10 5 interacting two-level systems. This opens up the possibility to perform accurate real-scale simulations of a diverse range of experiments in quantum optics or with solid-state spin ensembles under realistic laboratory conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Coupling Mesoscale Budget Components to Large-Eddy Simulations for Wind-Energy Applications

To simulate the airflow through a wind farm across a wide range of atmospheric conditions, microscale models (e.g., large-eddy simulation, LES, models) have to be coupled with mesoscale models, because microscale models lack the atmospheric physical processes to represent time-varying local forcing. Here we couple mesoscale model outputs to a LES solver by applying mesoscale momentum- and temperature-budget components from the Weather Research and Forecasting model to the governing equations of the Simulator fOr Wind Farm Applications model. We test whether averaging the budget components affects the LES results with regard to quantities of interest to wind energy. Our study focuses on flat terrain during a quiescent diurnal cycle. The simulation results are compared with observations from a 200-m tall meteorological tower and a wind-profiling radar, by analyzing time series, profiles, rotor-averaged quantities, and spectra. However, while results show that averaging reduces the spatio-temporal variability of the mesoscale momentum-budget components, when coupled with the LES model, the mesoscale bias (in comparison with observations of wind speed and direction, and potential temperature) is not reduced. In contrast, the LES technique can correct for shear and veer. In both cases, however, averaging the budget components shows no significant impact on the mean flow quantities in the microscale and is not necessary when coupling mesocale budget components to the LES model.

17 WIND ENERGY↗
DOE OSTI2025

NEML2: An efficient and modular multiphysics constitutive modeling library for hybrid computing environments

This paper presents NEML2, an open-source, high-performance library developed for constitutive material modeling, designed to support the flexible and modular development of models for complex material behavior. Building on the foundational structure of its predecessor, NEML, the NEML2 library introduces significant improvements, including enhanced vectorization, automatic differentiation, and seamless integration with PyTorch, facilitating the application of machine learning techniques in material simulations. NEML2 provides a C++ backend with Python bindings, enabling users to create custom material models that can be executed efficiently on both CPU and GPU platforms. The library also supports coupling with Multiphysics simulation frameworks like MOOSE, making it suitable for realistic simulations involving coupled physical processes. Rigorous quality assurance through unit and regression testing ensures the reliability of results, while the extensible, user-friendly design encourages collaboration and reproducibility across the scientific community. This paper provides an overview of NEML2’s architecture, core features, and applications, highlighting its impact on accelerating material qualification and advancing computational methods in materials science.

GPU↗
DOE OSTI2023

Kinetic simulation of a 50 mTorr capacitively coupled argon discharge over a range of frequencies and comparison to experiments

The effect of driving frequency in the range of 13.56–73 MHz on electron energy distribution and electron heating modes in a 50 mTorr capacitively coupled argon plasma discharge is studied using 1D-3V particle-in-cell simulations. Calculated electron energy probability functions exhibit three distinct “temperatures” for low-, mid-, and high-energy electrons at all the studied driving frequencies. When compared to published experimental data, the calculated probability functions show a reasonable agreement for the energy range resolved in the measurements (about 2–10 eV). Discrepancies due to limitations in experimental energy resolution outside this range lead to differences between computational and experimental values of the electron number density determined from the distribution functions, and the predicted effective electron temperature is within 25% of experimental values. The impedance of the discharge is interpreted in terms of a homogeneous equivalent circuit model, and the driving frequency dependence of the inferred combined sheath thickness is found to obey a known, theoretically derived, power law. The average power transferred from the field to the electrons (electron heating) is computed, and a region of negative heating near the sheath edge, particularly at higher driving frequencies, is identified. Analysis of the electron momentum equation shows that electron inertia, which on temporal averaging would be zero in a linear regime, is responsible for negative values of power deposition near the sheath edge at high driving frequencies due to the highly nonlinear behavior of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗