Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid simulator testing”

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 703 records · Page 39

Parallel Implementation of the Discontinuous Galerkin Method

This paper describes a parallel implementation of the discontinuous Galerkin method. Discontinuous Galerkin is a spatially compact method that retains its accuracy and robustness on non-smooth unstructured grids and is well suited for time dependent simulations. Several parallelization approaches are studied and evaluated. The most natural and symmetric of the approaches has been implemented in all object-oriented code used to simulate aeroacoustic scattering. The parallel implementation is MPI-based and has been tested on various parallel platforms such as the SGI Origin, IBM SP2, and clusters of SGI and Sun workstations. The scalability results presented for the SGI Origin show slightly superlinear speedup on a fixed-size problem due to cache effects.

Baggag, Abdalkader↗

Reinforcement Learning Applied to Cognitive Space Communications

The future of space exploration depends on robust, reliable communication systems. As the number of such communication systems increase, automation is fast becoming a requirement to achieve this goal. A reinforcement learning solution can be employed as a possible automation method for such systems. The goal of this study is to build a reinforcement learning algorithm which optimizes data throughput of a single actor. A training environment was created to simulate a link within the NASA Space Communication and Navigation (SCaN) infrastructure, using state of the art simulation tools developed by the SCaN Center for Engineering, Networks, Integration, and Communications (SCENIC) laboratory at NASA Glenn Research Center to obtain the closest possible representation of the real operating environment. Reinforcement learning was then used to train an agent inside this environment to maximize data throughput. The simulation environment contained a single actor in low earth orbit capable of communicating with twenty-five ground stations that compose the Near-Earth Network (NEN). Initial experiments showed promising training results, so additional complexity was added by augmenting simulation data with link fading profiles obtained from real communication events with the International Space Station. A grid search was performed to find the optimal hyperparameters and model architecture for the agent. Using the results of the grid search, an agent was trained on the augmented training data. Testing shows that the agent performs well inside the training environment and can be used as a foundation for future studies with added complexity and eventually tested in the real space environment.

Schubert, Carson D.↗

Uncertainty and Sensitivity Analysis of TREAT Transient Test no. 15

The U.S. Department of Energy (DOE) is preparing for the resumption of transient testing using the Transient Reactor Test Facility (TREAT) located at the Idaho National Laboratory (INL). In order to support the analysis of TREAT operation and experiment planning, the MAMMOTH reactor physics team has recently begun development of uncertainty and sensitivity analysis for the TREAT core simulated via MAMMOTH. MAMMOTH is a high-fidelity three-dimensional transient multi-physics code developed at INL under the DOE-NE Advanced Modeling and Simulation program (NEAMS). This research is part of the development presenting some initial uncertainty and sensitivity analysis of TREAT Transient Test no. 15. In this summary, we employ the Sobol sampling strategy implemented in RAVEN to efficiently perform the sensitivity and uncertainty analysis. The idea is using Sobol functional decomposition with generalized polynomial chaos (gPC) expansions to estimate the responses of interest. The gPC coefficients are then computed via the sparse grid spectral projection method. In addition, this strategy includes evaluation of the second or higher order sensitivity indices, allowing the study of interactions between several input variables. The principal focus of this analysis is to quantify the uncertainties in the total power and maximum fuel temperatures with respect to uncertainties in the initial transient conditions, such as rated power, fuel temperature, graphite fuel density and control rod boron concentration. In addition, both the first order and second order time-dependent sensitivity indices are computed and analyzed. The results show that the interactions between initial control rod boron concentration and the rated power will contribute at most 1% to the total uncertainty of the power peak.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Sensitivity Analysis of the Wind Forcing Effect on the Accuracy of Large-Wave Hindcasting

Deployment of wave energy converters (WECs) relies on consistent and accurate wave resource characterization, which is typically achieved through numerical modeling using deterministic wave models. The accurate predictions of large-wave events are critical to the success of wave resource characterization because of the risk on WEC installation, maintenance, and damage caused by extreme sea states. Because wind forcing is the primary driver of wave models, the quality of wind data plays an important role in the accuracy of wave predictions. This study evaluates the sensitivity of large-wave prediction to different wind-forcing products, and identifies a feasible approach to improve wave model results through improved wind forcing. Using a multi-level nested-grid modeling approach, we perform a series of sensitivity tests at four representative National Data Buoy Center buoy locations on the U.S. East and West Coasts. The selected wind-forcing products include the Climate Forecast System Reanalysis global wind product and North American Regional Reanalysis regional wind product as well as the observed wind at the buoys. Sensitivity test results indicate a consistent improvement in model predictions for the large-wave events (e.g., >90th percentile of significant wave height) at all buoys when observed-wind data were used to drive the wave model simulations.

16 TIDAL AND WAVE POWER↗

Modeling Framework for Data Center

This chapter highlights the critical need for advanced modeling of data centers due to their rapidly increasing energy consumption and impact on grid reliability. Driven by the demand for AI applications, data centers are projected to consume a significant portion of US energy by 2028, putting stress on an already challenged power grid. The chapter emphasizes the importance of "fast" time-scale models to understand the dynamic interactions between data centers and the grid, especially given the rapid power fluctuations of AI workloads. It outlines a modeling framework that includes both offline and real-time EMT domain simulations, detailing the necessary representations for various components like utility interfaces, transformers, IT loads, UPS, cooling loads, Battery Energy Storage Systems (BESS), generators, protection systems, and higher-level control systems. While standard simulation tools like PSCAD offer basic models, custom development is often required to accurately capture the unique and fast-changing behaviors of modern data centers. The chapter also discusses key metrics and test cases for validating these models, focusing on transient load responses, protection relay coordination, and demand flexibility. Finally, it addresses the challenges of modeling large-scale data centers, such as computational complexity and the trade-off between model fidelity and practicality, suggesting hybrid modeling approaches as a solution. The overarching goal is to create a robust framework that helps assess data center impacts on grid stability, identify vulnerabilities, and inform the development of standards for reliable integration of these large loads into the bulk power system.

25 ENERGY STORAGE↗

Data-driven deconvolution for large eddy simulations of Kraichnan turbulence

In this article, we demonstrate the use of artificial neural networks as optimal maps which are utilized for convolution and deconvolution of coarse-grained fields to account for sub-grid scale turbulence effects. Here, we demonstrate that an effective eddy-viscosity is predicted by our purely data-driven large eddy simulation framework without explicit utilization of phenomenological arguments. In addition, our data-driven framework precludes the knowledge of true sub-grid stress information during the training phase due to its focus on estimating an effective filter and its inverse so that grid-resolved variables may be related to direct numerical simulation data statistically. The proposed predictive framework is also combined with a statistical truncation mechanism for ensuring numerical realizability in an explicit formulation. Through this, we seek to unite structural and functional modeling strategies for modeling non-linear partial differential equations using reduced degrees of freedom. Both a priori and a posteriori results are shown for a two-dimensional decaying turbulence case in addition to a detailed description of validation and testing. A hyperparameter sensitivity study also shows that the proposed dual network framework simplifies learning complexity and is viable with exceedingly simple network architectures. Our findings indicate that the proposed framework approximates a robust and stable sub-grid closure which compares favorably to the Smagorinsky and Leith hypotheses for capturing the theoretical k -3 scaling in Kraichnan turbulence.

42 ENGINEERING↗

FUN3D Analyses of the 5th AIAA Propulsion Aerodynamics Workshop Inlet Test Case

This work presents the results of FUN3D analyses that were performed for the 5th AIAA Propulsion Aerodynamics Workshop. The workshop was separated into two sections that focused on inlet and nozzle flows. This paper focuses on the inlet section of the workshop, which considered an axisymmetric inlet in close-proximity to the ground and subjected to a crosswind. Four levels of grid refinement were provided by the workshop committee. Participants were asked to perform simulations for three crosswind velocities on a minimum of three levels of grid refinement, for a total of nine required cases. This work employed an updated set of grids provided by the committee after completion of the workshop based on feedback from participants. The updated grids featured a larger computational domain and increased refinement near the ground vortex region. This work employed the three coarsest grid levels from the updated grid set. Note that an additional extension downstream of the inlet was added to the provided grids for this work to mitigate observed convergence issues. Additionally, a limited turbulence model study was performed for the highest value of crosswind velocity and the finest grid. The results show that grid convergence was only achieved for the case with the lowest value of crosswind velocity. However, the results were generally observed to improve, relative to the experimental data, with increasing grid refinement. The FUN3D predictions were observed to agree qualitatively with the provided experimental data. However, FUN3D exhibited a tendency to underpredict the minimum value of the radially averaged total pressure at the Aerodynamic Interface Plane, which is consistent with the general findings from the workshop. Overall, FUN3D was found to perform well relative to its peers for the PAW5 workshop inlet test case.

PAI↗

FUN3D Analyses of the 5th AIAA Propulsion Aerodynamics Workshop Inlet Test Case

This work presents the results of FUN3D analyses that were performed for the 5th AIAA Propulsion Aerodynamics Workshop. The workshop was separated into two sections that focused on inlet and nozzle flows. This paper focuses on the inlet section of the workshop, which considered an axisymmetric inlet in close-proximity to the ground and subjected to a crosswind. Four levels of grid refinement were provided by the workshop committee. Participants were asked to perform simulations for three crosswind velocities on a minimum of three levels of grid refinement, for a total of nine required cases. This work employed an updated set of grids provided by the committee after completion of the workshop based on feedback from participants. The updated grids featured a larger computational domain and increased refinement near the ground vortex region. This work employed the three coarsest grid levels from the updated grid set. Note that an additional extension downstream of the inlet was added to the provided grids for this work to mitigate observed convergence issues. Additionally, a limited turbulence model study was performed for the highest value of crosswind velocity and the finest grid. The results show that grid convergence was only achieved for the case with the lowest value of crosswind velocity. However, the results were generally observed to improve, relative to the experimental data, with increasing grid refinement. The FUN3D predictions were observed to agree qualitatively with the provided experimental data. However, FUN3D exhibited a tendency to underpredict the minimum value of the radially averaged total pressure at the Aerodynamic Interface Plane, which is consistent with the general findings from the workshop. Overall, FUN3D was found to perform well relative to its peers for the PAW5 workshop inlet test case.

PAI↗

Analysis of iron-chromium-aluminum samples exposed to accident conditions followed by quench in the QUENCH-19 experiment

The QUENCH-19 experiment was a first-of-its-kind full-bundle test simulating accident conditions followed by water quench on accident-tolerant fuel (ATF) cladding. Here, a type of FeCrAl(Y) alloy, B136Y3, was developed at Oak Ridge National Laboratory and tested at the Karlsruhe Institute of Technology using Kanthal APM corner rods, a shroud, and Kanthal AF spacer grids. Testing conditions were similar to those in QUENCH-15—which tested ZIRLO cladding behavior—so that B136Y3 and ZIRLO cladding could be compared. QUENCH-19 consisted of an initial pre-oxidation heating followed by a transient. Then, a maximum power hold, which was not present in QUENCH-15, was executed to extend the heating period for the FeCrAl(Y) rods. Finally, a rapid water quench was executed that was similar to emergency core coolant system (ECCS) actuation. Compared with the ZIRLO rods in QUENCH-15, the bundle in QUENCH-19 released significantly less H 2 (9.2 g vs. 47.6 g) and achieved a much lower maximum temperature (1455°C vs. 1880°C). Furthermore, no breakaway oxidation was observed in QUENCH-19. Metallographic mounts revealed that despite the symmetry of the setup, at elevations near the maximum temperature, cladding and thermocouples were heavily damaged, substantial melting and oxidation occurred, and the cladding underwent chemical interaction with the thermocouple sheaths. Additionally, the ZrO 2 spacers detrimentally interacted with the cladding, leading to mixed oxide debris and the full destruction of some rods. Additional failure was found in certain cooler rods that may have risen due to the high thermal expansion coefficient of FeCrAl alloys. This paper presents an analysis of this work, which suggests that FeCrAl cladding can chemically survive anticipated loss-of-coolant accident events followed by rapid ECCS quench if the correct geometry and core design are present.

36 MATERIALS SCIENCE↗

Modeling Distributed Computing Infrastructures for HEP Applications

Predicting the performance of various infrastructure design options in complex federated infrastructures with computing sites distributed over a wide area network that support a plethora of users and workflows, such as the Worldwide LHC Computing Grid (WLCG), is not trivial. Due to the complexity and size of these infrastructures, it is not feasible to deploy experimental test-beds at large scales merely for the purpose of comparing and evaluating alternate designs. An alternative is to study the behaviours of these systems using simulation. This approach has been used successfully in the past to identify efficient and practical infrastructure designs for High Energy Physics (HEP). A prominent example is the Monarc simulation framework, which was used to study the initial structure of the WLCG. New simulation capabilities are needed to simulate large-scale heterogeneous computing systems with complex networks, data access and caching patterns. A modern tool to simulate HEP workloads that execute on distributed computing infrastructures based on the SimGrid and WRENCH simulation frameworks is outlined. Studies of its accuracy and scalability are presented using HEP as a case-study. Hypothetical adjustments to prevailing computing architectures in HEP are studied providing insights into the dynamics of a part of the WLCG and candidates for improvements.

Horzela, Maximilian↗

TREAT HODOSCOPE – A Predictive Simulation Model

The Transient Reactor Test Facility (TREAT) reactor is being restarted after over two decades of standby operation. TREAT is designed to test the time-dependent behavior of nuclear fuel during simulated reactor accidents. One of the principle data collection systems of TREAT is the fast neutron hodoscope which is used to observe the motion of the fuel test specimen as it is being damaged, deformed, or melted over the course of the power transient. The hodoscope is a grid of ex-core collimator channels paired with a bank of fast neutron sensitive radiation detectors. These detectors have historically been calibrated using empirical methods utilizing measurement data from smaller power pulses. Given modern computational capabilities this work sets out to establish a predictive computational capability of the fast neutron flux field exiting the hodoscope to support these calibrations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Direct Numerical Simulation of Hypersonic Turbulent Boundary Layer Flow using SPARC: Initial Evaluation

This report documents the initial testing of the Sandia Parallel Aerodynamics and Reentry Code (SPARC) to directly simulate hypersonic, turbulent boundary layer flow over a sharp 7- degree half-angle cone. This type of computation involves a tremendously large range of scales both in time and space, requiring a large number of grid cells and the efficient utilization of a large pool of resources. The goal of the simulation is to mimic and verify a wind tunnel experiment that seeks to measure the turbulent surface pressure fluctuations. These data are necessary for building a model to predict random vibration loading in the reentry flight environment. A low-dissipation flux scheme in SPARC is used on a 2.7 billion cell mesh to capture the turbulent fluctuations in the boundary layer flow. The grid is divided into 115200 partitions and simulated using the Knight's Landings (KNL) partition of the Trinity system. The parallel performance of SPARC is explored on the Trinity system, as well as some of the other new architectures. Extracting data from the simulation shows good agreement with the experiment as well as a colleague's simulation. The data provide a guide for which a new model can be built for better prediction of the reentry random vibration loads.

42 ENGINEERING↗

Domain connectivity among systems of overset grids

This progress report describes research performed to achieve algorithmic improvements in the use of overset grids for the computation of unsteady flowfields about geometrically complex and moving component configurations. A hole-map computational method is detailed that reduces computational demand and corresponding demand on human resources. Results indicate that the method is both efficient and flexible for use in complex simulation cases. Tests were performed on shuttle orbiter/external tank, wing/store, and V-22 tilt-rotor configurations. The number of simulated points and the corresponding time on a SGI 4D-210 workstation are presented.

Meakin, Robert L.↗

Study of the Protection Improvements for a Weak Grid Area With High Inverter-Based Resources (IBRs)

This project designs enhanced protection scheme for the real-world weak grid area with a high penetration of IBRs. As the existing protection schemes are originally designed for traditional synchronous machines, we first evaluate if the protection scheme will continue to operate reliably in systems with high levels of IBRs. Hardware relays are tested using a controller-hardware-in-the-loop setup. PSCAD electromagnetic transient simulation with IBR original equipment manufacturer black-box models is used to perform fault studies and generate COMTRADE data, which are replayed by a real-time digital simulator (RTDS) to feed input to the hardware relays. Three scenarios are analyzed: normal operation, an N-1 contingency, and an IBR-only scenario. The evaluation results reveal the following: 1) the protection scheme remains reliable under normal conditions and N-1 contingencies and 2) in IBR-only scenarios, differential protection (87L) continues to operate reliably, whereas local protection elements, such as distance and directional elements, fail because of the lack of regulated negative sequence current contributed by IBRs. Enhanced protection is designed to address the challenge of lack of negative sequence current from IBRs, including increased restraining factors a2 and k2 to block 32Q or using V instead QV ORDER for ground faults, enhanced mho distance element with voltage and phase angle supervision for L-L faults. The efficacy of enhanced protection logic is validated and proven to work reliably. Additionally, IEEE Std. 2800-2022 negative sequence current compliant GFL and GFM IBRs from another vendor are tested and proven to work reliably without need for enhanced logic. Therefore, this work provides valuable decision-making for utilities facing protection system challenges due to IBRs, either designing enhanced protection scheme or requesting their IBRs being IEEE Std. 2800-2022 compliant to produce regulated negative sequence current for protection relay to make correct decision.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Multiscale Modeling Framework Using Element‐Based Galerkin Methods for Moist Atmospheric Limited‐Area Simulations

This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. Numerical results from the MMF showed enhanced representation of cloud processes compared to the coarse model.

Kang, Soonpil [Naval Postgraduate School, Monterey↗

Numerical simulation of a hovering rotor using embedded grids

The flow field for a rotor blade in hover was computed by numerically solving the compressible thin-layer Navier-Stokes equations on embedded grids. In this work, three embedded grids were used to discretize the flow field - one for the rotor blade and two to convect the rotor wake. The computations were performed at two hovering test conditions, for a two-bladed rectangular rotor of aspect ratio six. The results compare fairly with experiment and illustrates the use of embedded grids in solving helicopter type flow fields.

Duque, Earl-Peter N.↗

Opportunities for Breakthroughs in Large-Scale Computational Simulation and Design

Opportunities for breakthroughs in the large-scale computational simulation and design of aerospace vehicles are presented. Computational fluid dynamics tools to be used within multidisciplinary analysis and design methods are emphasized. The opportunities stem from speedups and robustness improvements in the underlying unit operations associated with simulation (geometry modeling, grid generation, physical modeling, analysis, etc.). Further, an improved programming environment can synergistically integrate these unit operations to leverage the gains. The speedups result from reducing the problem setup time through geometry modeling and grid generation operations, and reducing the solution time through the operation counts associated with solving the discretized equations to a sufficient accuracy. The opportunities are addressed only at a general level here, but an extensive list of references containing further details is included. The opportunities discussed are being addressed through the Fast Adaptive Aerospace Tools (FAAST) element of the Advanced Systems Concept to Test (ASCoT) and the third Generation Reusable Launch Vehicles (RLV) projects at NASA Langley Research Center. The overall goal is to enable greater inroads into the design process with large-scale simulations.

Alexandrov, Natalia↗

HTR solver: An open-source exascale-oriented task-based multi-GPU high-order code for hypersonic aerothermodynamics

In this study, the open-source Hypersonics Task-based Research (HTR) solver for hypersonic aerothermodynamics is described. The physical formulation of the code includes thermochemical effects induced by high temperatures (vibrational excitation and chemical dissociation). The HTR solver uses high-order TENO-based spatial discretization on structured grids and efficient time integrators for stiff systems, is highly scalable in GPU-based supercomputers as a result of its implementation in the Regent/Legion stack, and is designed for direct numerical simulations of canonical hypersonic flows at high Reynolds numbers. Additionally, the performance of the HTR solver is tested with benchmark cases including inviscid vortex advection, low- and high-speed laminar boundary layers, inviscid one-dimensional compressible flows in shock tubes, supersonic turbulent channel flows, and hypersonic transitional boundary layers of both calorically perfect gases and dissociating air.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗