Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grid 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 865 records · Page 48

Simulating Transonic Buffet Aerodynamics for the Boeing Transonic Truss-Braced Wing Aircraft

The flow past the Boeing Transonic Truss-Braced Wing (TTBW) aircraft is simulated with a Hybrid Reynolds-Averaged Navier Stokes Large-Eddy Simulations (HRLES) turbulence modeling approach to study transonic buffet onset and evaluate the predictive capability of the numerical approach. All simulations used structured overset curvilinear grids for the free-air configuration, with all simulations performed using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) computational framework. Aerodynamic loads and surface pressure obtained from the HRLES solutions are compared to results from unsteady RANS (URANS) and experiments conducted at NASA Ames 11-by 11-Foot Transonic Wind Tunnel. Overall good agreement is obtained in the predicted loads and surface pressure. Unsteady pressure data from HRLES and URANS simulations is utilized to compute power spectral density (PSD) to predict transonic buffet onset and compared to the Kulite data obtained from the experiment. Comparisons of the PSD spectrum also show reasonable agreement with the wind tunnel data, with HRLES and URANS predicting buffet onset at a slightly earlier angle-of-attack than the experiment.

AATT↗

Impact of surface coupling grids on tropical cyclone extremes in high-resolution atmospheric simulations

This article discusses the sensitivity of tropical cyclone climatology to surface coupling strategy in high-resolution configurations of the Community Earth System Model. Using two supported model setups, we demonstrate that the choice of grid on which the lowest model level wind stress and surface fluxes are computed may lead to differences in cyclone strength in multi-decadal climate simulations, particularly for the most intense cyclones. Using a deterministic framework, we show that when these surface quantities are calculated on an ocean grid that is coarser than the atmosphere, the computed frictional stress is misaligned with wind vectors in individual atmospheric grid cells. This reduces the effective surface drag, and results in more intense cyclones when compared to a model configuration where the ocean and atmosphere are of equivalent resolution. Our results demonstrate that the choice of computation grid for atmosphere–ocean interactions is non-negligible when considering climate extremes at high horizontal resolution, especially when model components are on highly disparate grids.

54 ENVIRONMENTAL SCIENCES↗

Simulation Of Flow In A Turbine Cascade

Report discusses application of lower-upper (LU) implicit numerical-integration scheme, together with interactive grid-generation method, to simulation of inviscid, compressible flow in two-dimensional model of turbine cascade. Study part of continuing effort to develop techniques for computation of complicated three-dimensional flows in turbomachinery.

Choo, Yung K.↗

Optimizing the Resolution of Hydrodynamic Simulations for MCRaT Radiative Transfer Calculations

Despite their discovery about half a century ago, the Gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydro dynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affect the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

Jose Arita-Escalante↗

Sub-grid drag model for immersed vertical cylinders in fluidized beds

Immersed vertical cylinders are often used as heat exchanger in gas-solid fluidized beds. Computational Fluid Dynamics (CFD) simulations are computationally expensive for large scale systems with bundles of cylinders. Therefore sub-grid models are required to facilitate simulations on a coarse grid, where internal cylinders are treated as a porous medium. The influence of cylinders on the gas-solid flow tends to enhance segregation and affect the gas-solid drag. A correction to gas-solid drag must be modeled using a suitable sub-grid constitutive relationship. In the past, Sarkar et al. have developed a sub-grid drag model for horizontal cylinder arrays based on 2D simulations. However, the effect of a vertical cylinder arrangement was not considered due to computational complexities. In this study, highly resolved 3D simulations with vertical cylinders were performed in small periodic domains. These simulations were filtered to construct a sub-grid drag model which can then be implemented in coarse-grid simulations. Gas-solid drag was filtered for different solids fractions and a significant reduction in drag was identified when compared with simulation without cylinders and simulation with horizontal cylinders. Slip velocities significantly increase when vertical cylinders are present. Lastly, vertical suspension drag due to vertical cylinders is insignificant however substantial horizontal suspension drag is observed which is consistent to the finding for horizontal cylinders.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analysis of subgrid-scale eddy viscosity with use of results from direct numerical simulations

Without resort to any modeling, subgrid-scale eddy viscosity is computed from the results of high-resolution (64 cubed and 128 cubed grid points) direct numerical simulations of three-dimensional homogeneous isotropic decaying turbulence. In these simulations the eddy viscosity peaks sharply at the cutoff wavenumber, in rough agreement with the results of Kraichnan (1976). In addition, in the low-wavenumber range the eddy viscosity may be negative, contrary to the generally accepted concept of subgrid-scale eddy viscosity. Some possible explanations of this behavior are discussed.

Domaradzki, J. Andrzej↗

Analyzing at-scale distribution grid response to extreme temperatures

Threats against power grids continue to increase, as extreme weather conditions and natural disasters (extreme events) become more frequent. Hence, there is a need for the simulation and modeling of power grids to reflect realistic conditions during extreme events conditions, especially distribution systems. Herein, this paper presents a modeling and simulation platform for electric distribution grids which can estimate overall power demand during extreme weather conditions. The presented platform's efficacy is shown by demonstrating estimation of electrical demand for 1) Electricity Reliability Council of Texas (ERCOT) during winter storm Uri in 2021, and 2) alternative hypothetical scenarios of integrating Distributed Energy Resources (DERs), weatherization, and load electrification. In comparing to the actual demand served by ERCOT during the winter storm Uri of 2021, the proposed platform estimates approximately 34 GW of peak capacity deficit. These numbers are consistent with state-of-the-art prediction results published in the literature. For the case of the future electrification of heating loads, peak capacity of 78 GW (124% increase) is estimated, which would be reduced to 47 GW (38% increase) with the adoption of efficient heating appliances and improved thermal insulation. Integrating distributed solar PV and storage into the grid causes improvement in the local energy utilization and hence reduces the potential unmet energy by 31% and 40%, respectively.

24 POWER TRANSMISSION AND DISTRIBUTION↗

System-of-Systems Approach for Integrated Energy Systems Modeling and Simulation: Preprint

Today’s electricity grid is the most complex system ever built—and the future grid is likely to be even more complex because it will incorporate distributed energy resources (DERs) such as wind, solar, and various other sources of generation and energy storage. The complexity is further augmented by the possible evolution to new retail market structures that provide incentives to owners of DERs to support the grid. To understand and test new retail market structures and technologies such as DERs, demand-response equipment, and energy management systems while providing reliable electricity to all customers, an Integrated Energy System Model (IESM) is being developed at NREL. The IESM is composed of a power flow simulator (GridLAB-D), home energy management systems implemented using GAMS/Pyomo, a market layer, and hardware-in-the-loop simulation (testing appliances such as HVAC, dishwasher, etc.). The IESM is a system-of-systems (SoS) simulator wherein the constituent systems are brought together in a virtual testbed. We will describe an SoS approach for developing a distributed simulation environment. We will elaborate on the methodology and the control mechanisms used in the co-simulation illustrated by a case study.

97 MATHEMATICS AND COMPUTING↗

Evaluating Large-Eddy Simulations of Boundary-Layer Cloud Response to Climate Change using MAGIC Observations

Over the cool waters off California, a persistent and extensive deck of low-lying clouds forms atop a shallow, turbulent layer of moist air extending 500-2000 meters above the ocean, above which the air is much warmer and drier. Its west edge varies with weather and season and is typically half way to Hawaii. This cloud deck and similar cloud decks elsewhere reflect sunlight back to space, cooling the underlying ocean and indeed the entire planet. They must be accurately represented by a climate model to correctly simulate our current climate or future global warming. Most climate models suggest these cloud decks will thin and break up in a warmer CO 2 -enriched climate, so extra sunlight would reach the ocean surface and amplify warming worldwide. But these models have difficulty accurately simulating the cloud thickness and extent because the turbulent cloud-forming air currents are much too small to be simulated by their coarse computational grid. In 2012-2013, the DOE Atmospheric Radiation Measurement MAGIC program instrumented a container ship transiting between California and Hawaii to observe clouds and rain, small particles (aerosols), and weather conditions. Data were obtained from 20 4-6 day cruises in each direction. This study tested whether a sophisticated ‘large-eddy simulation’ (LES) of the turbulent air motions that create the clouds in this region accurately predicts the variations in cloud properties seen across the different cruises. If so, that model is more believable for predicting cloud changes in this region in a future warmer climate. Our LES used a numerical mesh over a 6.4x6.4 km region with 50 m horizontal and 5 m vertical grid spacing to simulate the cloud-forming eddies and the sharp overlying transition between moist and dry air. To better compare with MAGIC observations, we configured our model to simulate a region that moved along with the ship. We used a weather forecast model to specify effects of the large-scale air flow. We compared our results with balloon ascents, a radar and other instruments pointed upward at the clouds, measurements of downward sunlight and infrared radiation, aerosols, and weather variables for all suitable cruise periods. Our LES was remarkably skillful in predicting the cloudiness variations measured by the ship. MAGIC also sampled mesoscale cells, regions of cloud thickening and thinning over distances of 10-50 km. An LES run over a larger region 100 km on a side shows these cells form spontaneously, enhanced by strong radiative heat loss at the cloud tops.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the uncertainty introduced by discretization and time-averaging in two-fluid model predictions

The two-fluid model (TFM) has become a tool for the design and troubleshooting of industrial fluidized bed reactors. To use TFM for scale up with confidence, the uncertainty in its predictions must be quantified. Here, we study two sources of uncertainty: discretization and time-averaging. First, we show that successive grid refinement may not yield grid-independent transient quantities, including cross-section–averaged quantities. Successive grid refinement would yield grid-independent time-averaged quantities on sufficiently fine grids. A Richardson extrapolation can then be used to estimate the discretization error, and the grid convergence index gives an estimate of the uncertainty. Richardson extrapolation may not work for industrial-scale simulations that use coarse grids. We present an alternative method for coarse grids and assess its ability to estimate the discretization error. Second, we assess two methods (autocorrelation and binning) and find that the autocorrelation method is more reliable for estimating the uncertainty introduced by time-averaging TFM data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High-Fidelity Simulation of Turbulent Flow Past a Gaussian Bump

A spanwise-periodic computation of a turbulent flow past a Gaussian bump is performed in the form of a hybrid direct numerical simulation and wall-resolved large-eddy simulation. A fourth-order spatially-accurate flow solver is employed to perform the simulation, using 10.2 billion grid points for a Reynolds number of 170000 based on the bump height. The key findings from the simulation are reported in the acceleration and deceleration flow regions associated with the bump shape. Significant anisotropy in the normal Reynolds stresses, along both the wall-normal and streamwise directions, is observed within the acceleration region. The ratio between the Reynolds shear stress and turbulent kinetic energy in that region also experiences significant deviations from the norms of a zero pressure gradient turbulent boundary layer. The chosen Reynolds number generates strong flow separation in the adverse pressure gradient region, which is in contrast with a previous simulation at half the Reynolds number that only indicated incipient separation. An internal layer generated in the acceleration region evolves into a free shear layer that develops in the deceleration region and separates. Proper modeling of this inner layer appears crucial to predict the flow separation. Surface curvature effects on the attached flow development are also discussed.

Turbulent Boundary Layer↗

Simulation of Internal Combustion Engines with High-Performance Computing Tools

Traditional Lagrangian spray modeling approaches for internal combustion engines are highly griddependent due to insufficient resolution in the near nozzle region. This is primarily because of inherent restrictions of volume fraction with the Lagrangian assumption together with high computational costs associated with small grid sizes. A state-of-the-art grid-convergent spray modeling approach was developed and implemented by Senecal et al. (ASME-ICEF2012-92043) in the CONVERGE software. The key features of the methodology include Adaptive Mesh Refinement (AMR), advanced liquid-gas momentum coupling, and improved distribution of the liquid phase, which enables use of cell sizes smaller than the nozzle diameter. This modeling approach was rigorously validated against nonevaporating, evaporating, and reacting data from the literature. The current numerical study focusses on further demonstration of the grid-convergent modeling approach for simulating a single-cylinder Cat® compression ignition engine. The simulated injector is characterized with a nominal nozzle exit diameter of 259 μm. Simulations using various minimum grid sizes (ranging from 125 μm to 1000 μm) are compared for engine performance and emissions parameters of interest such as pressure, heat release rate, ignition delay, NOx, HC, and soot emissions. The peak cell-count for the highest resolution simulation was on the order of 34 million. These computationally expensive simulations were facilitated at a high-performance computing facility at Argonne National Laboratory. METIS load-balancing algorithm was developed and implemented in Converge code for the simulations. Scaling studies were also performed. The validity of previously recommended grid settings (ASME-ICEF2012-92043) for accuracy/runtime trade-off is further assessed. Efficacy of a simplified combustion model is also compared against a detailed chemical kinetics based combustion model.

33 ADVANCED PROPULSION SYSTEMS↗

Similar Estimates of Temperature Impacts on Global Wheat Yield by Three Independent Methods

The potential impact of global temperature change on global crop yield has recently been assessed with different methods. Here we show that grid-based and point-based simulations and statistical regressions (from historic records), without deliberate adaptation or CO2 fertilization effects, produce similar estimates of temperature impact on wheat yields at global and national scales. With a 1 C global temperature increase, global wheat yield is projected to decline between 4.1% and 6.4%. Projected relative temperature impacts from different methods were similar for major wheat-producing countries China, India, USA and France, but less so for Russia. Point-based and grid-based simulations, and to some extent the statistical regressions, were consistent in projecting that warmer regions are likely to suffer more yield loss with increasing temperature than cooler regions. By forming a multi-method ensemble, it was possible to quantify 'method uncertainty' in addition to model uncertainty. This significantly improves confidence in estimates of climate impacts on global food security.

Climate impacts↗

Fast Dynamic Simulation-Based Small Signal Stability Assessment and Control

Power grid planning and operation decisions are made based on simulation of the dynamic behavior of the system. Enabling substantial energy savings while increasing the reliability of the aging North American power grid through improved utilization of existing transmission assets hinges on the adoption of wide-area measurement systems (WAMS) for power system stabilization. However, adoption of WAMS alone will not suffice if the power system is to reach its full entitlement in stability and reliability. It is necessary to enhance predictability with "faster than real-time" dynamic simulations that will enable the dynamic stability margins, proactive real-time control, and improve grid resiliency to fast time-scale phenomena such as cascading network failures. Present-day dynamic simulations are performed only during offline planning studies, considering only worst case conditions such as summer peak, winter peak days, etc. With widespread deployment of renewable generation, controllable loads, energy storage devices and plug-in hybrid electric vehicles expected in the near future and greater integration of cyber infrastructure (communications, computation and control), monitoring and controlling the dynamic performance of the grid in real-time would become increasingly important. The state-of-the-art dynamic simulation tools have limited computational speed and are not suitable for real-time applications, given the large set of contingency conditions to be evaluated. These tools are optimized for best performance of single-processor computers, but the simulation is still several times slower than real-time due to its computational complexity. With recent significant advances in numerical methods and computational hardware, the expectations have been rising towards more efficient and faster techniques to be implemented in power system simulators. This is a natural expectation, given that the core solution algorithms of most commercial simulators were developed decades ago, when High Performance Computing (HPC) resources were not commonly available.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Computation Of Flow About A Helicopter Rotor

Improved method developed to simulate numerically flow about helicopter rotor with multiple blades. Computes vortical wake beneath rotor without using simplified ad hoc mathematical models of effects in wake. Patched grids used in numerical simulation of flow about rotor blades. Cylindrical configuration simplifies interpolation at interface between rotating and stationary grids.

Chen, C. L.↗

Open-Source PSCAD Grid-Following and Grid-Forming Inverters and a Benchmark for Zero-Inertia Power System Simulations: Preprint

This paper presents open-source, flexible, and easily-scalable models of grid following and grid forming inverters for the PSCAD software platform. The models are intended for system integration studies, particularly stability analyses of power systems with high penetration of inverter-based generation. To verify the model functionality, they are implemented in a IEEE9-bus system in a zero-inertia operational scenario of 100% inverter-based generation. The models have been made available open source at the PyPSCAD NREL GitHub page.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Characterizing uncertainties in ground “truth” of precipitation over complex terrain through high-resolution numerical modeling

Ground observation of precipitation over complex terrain is subject to large uncertainties due to inadequate sampling. This study explores a method that combines limited gauge data and a high-resolution numerical simulation to quantify the precipitation uncertainties in central Himalaya. Specifically, the Coefficient of spatial Variability (CV) of precipitation and the minimum Number of Required Stations (NRS) to obtain areal-mean precipitation ground truth values within a 0.25-degree area are investigated using fine-scale meteorological simulation at 1.5 km grid spacing. Evaluation over a densely-gauged area demonstrates comparable CV and NRS values between station observations and simulations. The simulated CV and NRS values in a larger area show a strong and positive dependence on each other and an expected positive (negative) correlation with topographic complexity (temporal scale). The proposed method sheds lights on evaluating precipitation products and holds promise for informing the layout of rain gauge networks in complex terrain.

Ouyang, Lin↗