Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation fidelity”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

High Fidelity Simulations of Air-Cooled Reactor Cavity Cooling System

High Temperature Gas Reactor (HTGR) designs incorporate passive safety systems (e.g., the Reactor Cavity Cooling System [RCCS]) that utilize natural principles to manage heat dissipation from the reactor pressure vessel (RPV) during accidents or routine shutdowns. The industry community is experiencing a pressing need for advanced simulation tools that can accurately assess the performance of these types of systems. In the literature, a knowledge gap exists concerning high-fidelity data for the RCCS, and this gap is one of the areas of focus of the present study. This research focuses on a specific RCCS designed for General Atomics' Modular High-Temperature Gas Reactor (GA-MHTGR). Experimental studies on a scaled version (can be seen in Figure 1) of the air-cooled RCCS used in GA-MHTGR were conducted by the University of Wisconsin-Madison (UW-Madison). This work contributes to a broader initiative aimed at establishing a numerical benchmark based on the UW-Madison experiments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

High-fidelity simulations of the run-in process for a pebble-bed reactor

Pebble-bed reactors (PBRs) rely on a continual feed of fuel pebbles being cycled through the core. As a result, they require a “run-in” period in order to reach an equilibrium state. The run-in period for a PBR is a complex, time-dependent problem that requires the injection of new fuel, different types of fuel, and power increases. This complexity in the run-in makes it important to capture the physical processes in order to generate an accurate representation. The present work details the creation of a high-fidelity Monte Carlo methodology for analyzing the run-in and subsequent approach to equilibrium for PBRs. The methodology entails a Python module wrapped around Serpent so as to perform neutronics calculations, move pebbles, refuel the core, and discharge pebbles, thereby modeling the explicit behavior of the PBR run-in. Further, three run-in simulations (a constant temperature profile, a linear temperature profile, and a constant temperature profile using control rods) were examined in order to identify the key physical phenomena present in the run-in process. Utilizing kugelpy, we found the inclusion of a temperature profile to be important for accurately capturing a discharge burnup (around 141 MWd/kg), a consistent k-eff (around 1.005), and an average pebble power (around 2.5 kW/pebble) that all fall within acceptable limits.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Data-driven learning of nonlocal models: from high-fidelity simulations to constitutive laws

We show that machine learning can improve the accuracy of simulations of stress waves in one-dimensional composite materials. We propose a data-driven technique to learn nonlocal constitutive laws for stress wave propagation models. The method is an optimization-based technique in which the nonlocal kernel function is approximated via Bernstein polynomials. The kernel, including both its functional form and parameters, is derived so that when used in a nonlocal solver, it generates solutions that closely match high-fidelity data. The optimal kernel therefore acts as a homogenized nonlocal continuum model that accurately reproduces wave motion in a smaller-scale, more detailed model that can include multiple materials. We apply this technique to wave propagation within a heterogeneous bar with a periodic microstructure. Several one-dimensional numerical tests illustrate the accuracy of our algorithm. The optimal kernel is demonstrated to reproduce high-fidelity data for a composite material in applications that are substantially different from the problems used as training data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

High-Fidelity Simulations of Gravity Currents Using a High-Order Finite-Difference Spectral Vanishing Viscosity Approach

This numerical work investigates the potential of a high-order finite-difference spectral vanishing viscosity approach to simulate gravity currents at high Reynolds numbers. The method introduces targeted numerical dissipation at small scales through altering the discretisation of the second derivatives of the viscous terms in the incompressible Navier-Stokes equations to mimic the spectral vanishing viscosity (SVV) operator, originally designed for the regularisation of spectral element method (SEM) solutions of pure advection problems. Using a sixth-order accurate finite-difference scheme, the adoption of the SVV method is straightforward and comes with a negligible additional computational cost. In order to assess the ability of this high-order finite-difference spectral vanishing viscosity approach, we performed large-eddy simulations (LES) of a gravity current in a channelised lock-exchange set-up with our SVV model and with the well-known explicit static and dynamic Smagorinsky sub-grid scale (SGS) models. The obtained data are compared with a direct numerical simulation (DNS) based on more than 800 million mesh nodes, and with experimental measurements. A framework for the energy budget is introduced to investigate the behaviour of the gravity current. First, it is found that the DNS is in good agreement with the experimental data for the evolution of the front location and velocity field as well as for the stirring and mixing inside the gravity current. Secondly, the LES performed with less than 0.4% of the total number of mesh nodes compared to the DNS, can reproduce the main features of the gravity currents, with the SVV model yielding slightly more accurate results. It is also found that the dynamic Smagorinsky model performs better than its static version. For the present study, the static and dynamic Smagorinsky models are 1.8 and 2.5 times more expensive than the SVV model, because the latter does not require the calculation of explicit SGS terms in the Navier-Stokes equations nor spatial filtering operations.

direct numerical simulation↗
DOE OSTI2023

Insights into the failure mode of the Chelyabinsk meteor from high-fidelity simulation

In 2013, an asteroid with a diameter of approximate 20 m broke up in the atmosphere above the Russian city of Chelyabinsk producing what is known as an airburst event. The event was recorded on video from multiple angles and a 500-kg fragment was recovered from Lake Chebarkul shortly thereafter. These aspects, coupled with the asteroid’s large size, distinguish the Chelyabinsk event as unique. In this study, we use the smoothed particle hydrodynamics method to examine the breakup mode of a Chelyabinsk-sized monolithic asteroid and probe whether a monolithic structure is consistent with observations from the Chelyabinsk event. We model the object as a homogeneous sphere with properties consistent with Chelyabinsk meteorites. We use planar-2D simulations to test the sensitivity of results to variations in the material model and to qualitatively inspect the cause of crack nucleation and fragmentation. We then compare results from a 3D simulation to observations from the event. Our 3D simulation predicts a burst height within 3 km of the observed peak luminosity suggesting the responsible asteroid may have been monolithic.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Experimentally validated high-fidelity simulations of a liquid jet in supersonic crossflow

Here, utilizing recent advancements in computational schemes for compressible, multiphase flows, this work features a parametric study of a pure liquid jet in supersonic crossflow that involves simulating the atomization process for four values of momentum-flux ratio. These simulations are validated against experimental results measured with high-speed X-ray imaging, which confirm the accuracy of the numerical approach. Also, the effect of numerical resolution on some flow behavior is investigated, revealing convergence of the jet shape and surface instability wavelength. Analysis of the resulting sprays includes statistical descriptions of the liquid distribution, liquid structures created through breakup, interfacial instabilities, and dominant flow features. As the flowrate increases, the spray penetrates further, it becomes more disperse, and less liquid impacts the wall, but the droplet size distribution changes little. The wavelength of instabilities on the windward side of the jet diverges from measured trends in subsonic crossflows. In a visualization of the time-averaged flow, counter-rotating vortices are observed along the jet core and in the wake, affecting the process of primary atomization and early droplet trajectories.

42 ENGINEERING↗
DOE OSTI2024

High-Fidelity Simulation of the Light-to-Dense Stratification Transient in the HiRJET Facility

Density stratification in a large enclosure is a crucial phenomenon to heat transfer and sustainable passive heat removal of a sodium fast reactor during reactivity transients. However, engineering turbulence models were identified to have unsatisfactory performance in predicting propagation of a stratified front. Yet, the scarcity of high-resolution data for stratification hampers the development of models. To explor e applications of leveraging direct numerical simulation (DNS) data to support turbulence model development, this work conducted DNS using NekRS to study a long stratification transient in the High-Resolution Jet (HiRJET) experimental facility. Here, this work considers an experiment run where light fluid is injected into a tank containing a denser fluid with a relative density difference of 1.5%. Formation of the stratified layer is identified as impingement of the buoyant jet promoting mixing of the two fluids. Based on the transient statistics, transport of the concentration can be characterized by regions with dominating effects of turbulent mixing, buoyant dissipation, and molecular diffusion, respectively, as moving away from the elevation of jet impingement. Concentration near the stratified front also exhibits oscillation at Brunt-Väisälä frequency. Preliminary validation of the simulation showed encouraging agreement of the concentration distribution with the reference experiment.

buoyant mixing↗
DOE OSTI2024

Hyperparameter Studies for Vision Transformers Trained on High-Fidelity Simulations

This library is a collection of python modules that define, train, and analyze vision-transformer (ViT) machine learning models. The code implements, with mild modifications, ViT models that have been made publicly available through publication and GitHub code. The training data for these models is hydrodynamic simulation output in the form of numpy arrays. This library contains code to train these ViT models on the hydrodynamic simulation output with a variety of hyperparameters, and to compare the results of such models. Furthermore, the library contains definitions of simple convolutional neural network (CNN) machine learning architectures which can be trained on the same hydrodynamic simulation output. These are included as a reference point to compare the ViT models to. Additionally, the library includes trained ViT and CNN models and example input data for demonstration purposes. The code is based on the PyTorch python library.

Callis, Skylar↗
DOE OSTI2019

Spray characterization for engine combustion network Spray G injector using high-fidelity simulation with detailed injector geometry

This article presents a computational fluid dynamics study of the engine combustion network Spray G, focusing on the transient characteristics of the spray during the start of injection and the impacts of nozzle geometry details derived from the manufacturing process. The large-eddy-simulation method, coupled with the volume-of-fluid method, was used to model the high-speed turbulent two-phase flow. A moving-needle boundary condition was applied to capture the internal flow boundary condition accurately. The injector geometry was measured with micron-level resolution using X-ray tomographic imaging at the Advanced Photon Source at Argonne National Laboratory, providing detailed machining tolerance and defects from manufacturing and a realistic rough surface. For comparison, a nominal geometry and a modified geometry incorporating measured large-scale geometric features but no surface details were also used in the simulations. Spray characteristics such as mass flow rate, injection velocity, and Sauter mean diameter were analyzed. Significantly distinct spray characteristics in terms of injection velocity, spray morphology, and primary breakup mechanism were predicted using different nozzle geometries, which is mainly attributable to the realistic surface finish and manufacturing defects. The measured high-resolution geometry predicts a lower injection velocity, a wider-spreading spray, and an overall slower breakup rate with evident injector tip wetting compared to the ideally smooth nozzle boundary. This result implies that the manufacturing details of the injector, which are usually ignored in fuel injection studies, have a significant impact on the spray development process and should be taken into account for design optimization.

Yue, Zongyu↗
DOE OSTI2019

Investigation of Marangoni convection with high-fidelity simulations for metal melt pool dynamics

The goal of this research is to investigate and characterize the dynamics of a molten metal pool through modeling and simulation. Two manufacturing processes are of interest, laser spot welding and selective laser melting. These are laser-based techniques for joining parts together and for building up parts. The laser power impacts the melt pool temperature and therefore the dynamics. Moreover, the metals used are often alloys, containing impurities, where the concentration of these impurities can have a large effect on surface tension and therefore influence the melt pool dynamics. The application space considered in this project is for laser spot welding as a case study, and is also relevant to 3D metal printing. The thermal evolution and fluid dynamics in a 2D metal block exposed to a Gaussian distributed laser source is numerically simulated using a high-order, fully-implicit code. Results are compared to published experimental data and other published numerical simulations for laser spot welding. Results confirm that unstable structures appear at higher laser powers with higher impurity concentrations compared to simulations with no impurities and lower laser powers.

42 ENGINEERING↗