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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Enhancing Scalability for FUN3D Rotorcraft Simulations with Yoga: an Overset Grid Assembler

FUN3D, an unstructured grid Navier-Stokes CFD code is capable of overset grid simu- lations, but does not have an internal method for assembling an overset grid system from a group of component grids. FUN3D currently relies on the third party codes Suggar++ and DiRTlib to perform domain assembly and provide intergrid connectivity for overset simulations. Rotorcraft simulations with moving, deforming blades require domain assembly and mesh deformation at each time step. For these simulations, the three primary drivers of computational cost for each time step are: deforming the mesh, performing domain assembly, and performing subiterations of the flow solver. FUN3D exhibits strong and weak scalability for the flow solver subiterations and the mesh deformation. However, FUN3D is currently hardwired directly to the serial version of Suggar++, which has a fixed cost for a given mesh system. Therefore, domain assembly begins to dominate the total cost of each time step as grid systems become larger. An integrated method for parallel domain assembly is presented that addresses scalability for large grid systems. Restructuring within FUN3D to accommodate integrated domain assembly is also discussed, which could enable use of the parallel Suggar++ library.

Cameron T Druyor↗

NASA NTRS2013

Hybrid Large-Eddy/Reynolds-Averaged Simulation of a Supersonic Cavity Using VULCAN

Simulations of a supersonic recessed-cavity flow are performed using a hybrid large-eddy/Reynolds-averaged simulation approach utilizing an inflow turbulence recycling procedure and hybridized inviscid flux scheme. Calorically perfect air enters a three-dimensional domain at a free stream Mach number of 2.92. Simulations are performed to assess grid sensitivity of the solution, efficacy of the turbulence recycling, and the effect of the shock sensor used with the hybridized inviscid flux scheme. Analysis of the turbulent boundary layer upstream of the rearward-facing step for each case indicates excellent agreement with theoretical predictions. Mean velocity and pressure results are compared to Reynolds-averaged simulations and experimental data for each case and indicate good agreement on the finest grid. Simulations are repeated on a coarsened grid, and results indicate strong grid density sensitivity. Simulations are performed with and without inflow turbulence recycling on the coarse grid to isolate the effect of the recycling procedure, which is demonstrably critical to capturing the relevant shear layer dynamics. Shock sensor formulations of Ducros and Larsson are found to predict mean flow statistics equally well.

Quinlan, Jesse↗

NASA NTRS2011

Evaluation of the Monotonic Lagrangian Grid and Lat-Long Grid for Air Traffic Management

The Air Traffic Monotonic Lagrangian Grid (ATMLG) is used to simulate a 24 hour period of air traffic flow in the National Airspace System (NAS). During this time period, there are 41,594 flights over the United States, and the flight plan information (departure and arrival airports and times, and waypoints along the way) are obtained from an Federal Aviation Administration (FAA) Enhanced Traffic Management System (ETMS) dataset. Two simulation procedures are tested and compared: one based on the Monotonic Lagrangian Grid (MLG), and the other based on the stationary Latitude-Longitude (Lat- Long) grid. Simulating one full day of air traffic over the United States required the following amounts of CPU time on a single processor of an SGI Altix: 88 s for the MLG method, and 163 s for the Lat-Long grid method. We present a discussion of the amount of CPU time required for each of the simulation processes (updating aircraft trajectories, sorting, conflict detection and resolution, etc.), and show that the main advantage of the MLG method is that it is a general sorting algorithm that can sort on multiple properties. We discuss how many MLG neighbors must be considered in the separation assurance procedure in order to ensure a five-mile separation buffer between aircraft, and we investigate the effect of removing waypoints from aircraft trajectories. When aircraft choose their own trajectory, there are more flights with shorter duration times and fewer CD&R maneuvers, resulting in significant fuel savings.

Kaplan, Carolyn↗

DOE OSTI2022

Transmission-Distribution Dynamic Co-Simulation of Electric Vehicles Providing Grid Frequency Response

This paper investigates the impacts of electric vehicles (EVs) on power system frequency regulation based on an open-source transmission-and-distribution (T&D) dynamic co-simulation framework. The development of a EV dynamic model based on appropriate WECC dynamic model is introduced first, then the T&D dynamic co-simulation platform is described. The advantage of the overall platform is that DERs such as distributed photovoltaic (DPV) and EV are modeled explicitly in both transmission and distribution simulators for frequency dynamics and voltage, respectively. In the case studies, EVs frequency responses (i.e., primary and/or secondary) are simulated after the system is exposed to N-1 contingency such as a generation trip. Various EV frequency regulation participation strategies are investigated as well to study their impacts on system frequency response. The studies shows that EVs have potential capability to provide effective frequency regulation services.

electric vehicle↗

DOE OSTI2022

Transmission-Distribution Dynamic Co-Simulation of Electric Vehicles Providing Grid Frequency Response: Preprint

electric vehicle↗

NASA NTRS2013

Understanding the Impact of Field-Emitter Characteristics on Electron Beam Focusing in the VAPoR Time-of-Fight Mass Spectrometer

Simulations of field emission of electrons from an electron gun are used to determine the angular distribution of the emitted electron beam and the percentage of charge transmitted through the grid. The simulations are a first step towards understanding the spherical aberration present after focusing the electron beam. The effect of offset of the cathode with respect to the grid and the separation between cathode and grid on the angular distributions of emitted electrons and transmission of the grid are explored.

electron gun↗

NASA NTRS1999

Large Eddy Simulation of Cirrus Clouds

The Regional Atmospheric Modeling System (RAMS) with mesoscale interactive nested-grids and a Large-Eddy Simulation (LES) version of RAMS, coupled to two-moment microphysics and a new two-stream radiative code were used to investigate the dynamic, microphysical, and radiative aspects of the November 26, 1991 cirrus event. Wu (1998) describes the results of that research in full detail and is enclosed as Appendix 1. The mesoscale nested grid simulation successfully reproduced the large scale circulation as compared to the Mesoscale Analysis and Prediction System's (MAPS) analyses and other observations. Three cloud bands which match nicely to the three cloud lines identified in an observational study (Mace et al., 1995) are predicted on Grid #2 of the nested grids, even though the mesoscale simulation predicts a larger west-east cloud width than what was observed. Large-eddy simulations (LES) were performed to study the dynamical, microphysical, and radiative processes in the 26 November 1991 FIRE 11 cirrus event. The LES model is based on the RAMS version 3b developed at Colorado State University. It includes a new radiation scheme developed by Harrington (1997) and a new subgrid scale model developed by Kosovic (1996). The LES model simulated a single cloud layer for Case 1 and a two-layer cloud structure for Case 2. The simulations demonstrated that latent heat release can play a significant role in the formation and development of cirrus clouds. For the thin cirrus in Case 1, the latent heat release was insufficient for the cirrus clouds to become positively buoyant. However, in some special cases such as Case 2, positively buoyant cells can be embedded within the cirrus layers. These cells were so active that the rising updraft induced its own pressure perturbations that affected the cloud evolution. Vertical profiles of the total radiative and latent heating rates indicated that for well developed, deep, and active cirrus clouds, radiative cooling and latent heating could be comparable in magnitude in the cloudy layer. This implies that latent heating cannot be neglected in the construction of a cirrus cloud model. The probability density function (PDF) of w was analyzed to assist in the parameterization of cloud-scale velocities in large-scale models. For the more radiatively-driven, thin cirrus case, the PDFs are approximately Gaussian. However, in the interior of the deep, convectively unstable case, the PDFs of w are multi-modal and very broad, indicating that parameterizing cloud-scale motions for such clouds can be very challenging. The results of this research are described in detail in a paper submitted to the Journal of Atmospheric Science (Wu and Cotton, 1999), which is enclosed as Appendix 2. Using soundings extracted from a mesoscale simulation of the November 26, 1991 cirrus event, the radiative effects on vapor deposition/sublimation of ice crystals was studied using a two-dimensional cloud-resolving model (CRM) version of RAMS, coupled to an explicit bin-resolving microphysics. The CRM simulations of the November 26, 1991 cirrus event demonstrate that the radiative impact on the diffusional growth (or sublimation) of ice crystals is significant. In this case, the ice particles experienced radiative warming. Model results show that radiative feedbacks in the diffusional growth of ice particles can be very complex. Radiative warming of an ice particle will restrict the particle's diffusional growth. In the case of radiative warming, ice particles larger than a certain size will experience so much radiative warming that surface ice saturation vapor pressures become large enough to cause sublimation of the larger crystals, while smaller crystals are growing by vapor deposition. However, ice mass production can be enhanced in the case of radiative cooling of an ice particle. For the November 26, 1991 cirrus event, radiative feedback results in significant reduction in the total ice mass, especially in the production of large ice crystals, and consequently, both radiative and dynamic properties of the cirrus cloud are significantly affected. A complete description of this research has been submitted as a paper to the Journal of Atmospheric Science (Wu et al., 1999), and included as Appendix 3.

Wu, Ting↗

DOE OSTI2022

Experimental Setup for Grid Control Device Software Updates in Supply Chain Cyber-Security

Supply chain cyberattacks that exploit insecure third-party software are a growing concern for the security of the electric power grid. These attacks seek to deploy malicious software in grid control devices during the fabrication, shipment, installation, and maintenance stages, or as part of routine software updates. Malicious software on grid control devices may inject bad data or execute bad commands, which can cause blackouts and damage power equipment. This paper describes an experimental setup to simulate the software update process of a commercial power relay as part of a hardware-in-the-loop simulation for grid supply chain cyber-security assessment. The laboratory setup was successfully utilized to study three supply chain cyber-security use cases.

Keller, Joseph↗

DOE OSTI2024

Positive-Sequence Modeling of Droop-Controlled Grid-Forming Inverters for Transient Stability Simulation of Transmission Systems

Here, this paper describes a positive-sequence model to represent two widely reported droop-controlled grid-forming inverters for bulk power systems simulation study. Methods of how to develop the equivalent voltage source behind impedance to represent inverters with and without inner control loops, modeling of P-f and Q-V droop controls, active and reactive power limiting modeling, and modeling of fault current limiting controls are described in detail. The model is implemented in commercially available positive-sequence simulation tools and compared with detailed electromagnetic transient models in both a single-grid-forming-inverter infinite-bus system and a modified IEEE 39-bus system. Finally, the model is tested on the full US Western Electricity Coordinating Council (WECC) system. Study results show that the model has a good level of accuracy compared to the detailed EMT models, and at the same time also achieves a high computational efficiency, which is suitable for the bulk power system simulation study.

42 ENGINEERING↗

NASA NTRS2015

DNS/LES Simulations of Separated Flows at High Reynolds Numbers

Direct numerical simulations (DNS) and large-eddy simulations (LES) simulations of flow through a periodic channel with a constriction are performed using the dynamic Smagorinsky model at two Reynolds numbers of 2800 and 10595. The LES equations are solved using higher order compact schemes. DNS are performed for the lower Reynolds number case using a fine grid and the data are used to validate the LES results obtained with a coarse and a medium size grid. LES simulations are also performed for the higher Reynolds number case using a coarse and a medium size grid. The results are compared with an existing reference data set. The DNS and LES results agreed well with the reference data. Reynolds stresses, sub-grid eddy viscosity, and the budgets for the turbulent kinetic energy are also presented. It is found that the turbulent fluctuations in the normal and spanwise directions have the same magnitude. The turbulent kinetic energy budget shows that the production peaks near the separation point region and the production to dissipation ratio is very high on the order of five in this region. It is also observed that the production is balanced by the advection, diffusion, and dissipation in the shear layer region. The dominant term is the turbulent diffusion that is about two times the molecular dissipation.

Balakumar, P.↗

NASA NTRS1982

Numerical simulation of a conducting disk partially covered by an insulator interacting with a plasma

A numerical simulation has been made of the 'pinhole effect' which produces the characteristic 'S-shaped' current-voltage curve. A disk-shaped conducting probe immersed in a plasma is modeled using a particle-in-cell (PIC) code. A probe partially covered by a very thin insulating layer is considered, as well as a probe mounted on an insulating disk. The simulation uses a cylindrical particle mover and allows for a variable number of particles in the system. The simulation space grid uses three different mesh sizes, the coarsest being away from the probe and the finest near the probe, in order to accurately calculate the trajectories of the simulation particles contributing the current to the probe and to the surface charge density on the dielectric. The calculation of the electrostatic potential is done self-consistently using successive over-relaxation (SOR). Backscattering and secondary electron emission are included for the case of positive probe voltage.

Chaky, R. C.↗

DOE OSTI2024

Investigating the effects of cooperative transmission expansion planning on grid performance during heat waves with varying spatial scales

There is growing recognition of the advantages of interregional transmission capacity to decarbonize electricity grids. A less explored benefit is potential performance improvements during extreme weather events. This study examines the impacts of cooperative transmission expansion planning using an advanced modeling chain to simulate power grid operations of the United States Western Interconnection in 2019 and 2059 under different levels of collaboration between transmission planning regions. Two historical heat waves in 2019 with varying geographical coverage are replayed under future climate change in 2059 to assess the transmission cooperation benefits during grid stress. The results show that cooperative transmission planning yields the best outcomes in terms of reducing wholesale electricity prices and minimizing energy outages both for the whole interconnection and individual transmission planning regions. Compared to individual planning, cooperative planning reduces wholesale electricity prices by 64.3 % and interconnection-wide total costs (transmission investments + grid operations) by 34.6 % in 2059. It also helps decrease greenhouse gas emissions by increasing renewable energy utilization. However, the benefits of cooperation diminish during the widespread heat wave when all regions face extreme electricity demand due to higher space cooling needs. Despite this, cooperative transmission planning remains advantageous, particularly for California Independent System Operator with significant diurnal solar generation capacity. This study suggests that cooperation in transmission planning is crucial for reducing costs and increasing reliability both during normal periods and extreme weather events. It highlights the importance of optimizing the strategic investments to mitigate challenges posed by wider-scale extreme weather events of the future.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

PIAFS: A 2D nonlinear hydrodynamics code to model gaseous optics

The survivability of final optics is expected to be a major challenge for all future inertial fusion energy concepts. Due to their higher damage threshold, gaseous optics have been identified as a promising solution to this problem. Gaseous optics can be created through the photoabsorption of spatially modulated UV light, which induces various chemical processes that heat the gas. This heating leads to a pressure perturbation, which in turn launches a density perturbation that can imprint a refractive index modulation such as a grating. In this article, we introduce a parallel C/C++ code to simulate gaseous optics. PIAFS2D is a high-order conservative finite-difference code to solve the compressible Navier–Stokes equations along with the photochemical heating sources on Cartesian grids. The simulations are validated by the linear theory derived in a previous paper [Michel et al., Phys. Rev. Appl. 22, 024014 (2024)]. For larger perturbations, the behavior of the system—particularly the evolution of the generated acoustic wave—demonstrates strong nonlinearity. PIAFS2D allows the study of nonlinear behaviors and can be used for the design of high-efficiency gaseous optics elements in realistic experimental conditions.

Oudin, A. [Lawrence Livermore National Laboratory ↗

DOE OSTI2024

SNNPG: Using Spiking Neural Networks to Detect Attacks in the Power Grid

We explore the potential of Spiking Neural Networks (SNN) to enhance the security of power grid operations by detecting False Data Injection (FDI) attacks. These attacks manipulate PMU readings, leading to erroneous control decisions and grid disruptions. We develop a method to convert Phase Measurement Unit (PMU) data into spike trains, capturing both temporal and spatial dimensions. Using an SNN model, we conduct evaluations with simulated power grid data, showcasing accuracy in detecting FDI attacks. SNN models rapidly identify anomalies in real-time PMU data, safeguarding grid operations by alerting operators to irregular readings and preventing incorrect decisions.

artificial intelligence↗

DOE OSTI2022

Quantifying the effects of varietal types × management on the spatial variability of sorghum biomass across US environments

Abstract Regional‐scale estimations of sorghum biomass production allow identification of optimum genotype×environment×management (G×E×M) combinations for bioenergy generation. The objective of this study was to determine the degree of contributions of G, E, and M toward variability in sorghum biomass in the United States. Using the Agricultural Production Systems sIMulator in a grid computing platform, biomass was simulated for irrigated and rainfed conditions for 30 years across the United States for four sorghum varietal types (grain—GS, sudangrass—SS, photosensitive—PS, and photo‐insensitive—PI). Simulated biomass was assessed by environments clustered using the sum of intercepted solar radiation ( ir ), mean of temperature stress factor ( tp ) and water stress factor ( sw ). Simulated biomass ranged from 5.8 t ha −1 (GS‐rainfed) to 27.5 t ha −1 (PI‐irrigated). Under high‐temperature environments (mean annual temperature = 25°C), rainfed biomass between 40 and 80 days after planting (DAP) was strongly correlated with sw ( r = 0.64–0.86) and irrigated biomass with ir ( r = 0.68–0.81). Under low‐temperature environments (mean annual temperature = 18°C) after 40 DAP, tp and ir had greater effects than sw ( r = 0.55–0.82). Biomass variance was mainly explained by varietal type (50%–76%) in all environments×irrigation combinations, except in the high‐ and mid‐temperature environments under rainfed conditions where rainfall had the major effect (25%–45%). However, when mean temperature during the growing season decreased from 25°C (high environments) to 18°C (low environments), the contribution of mean temperature to biomass variance increased from 7% to 34% (rainfed) and from 4% to 36% (irrigated). Varietal type had the larger interactions with other factors independently of the environment and irrigation. We demonstrated a need to quantify (i) the main G×E×M drivers of biomass variability based on environmental stress factors and (ii) the variance contribution of these drivers on sorghum biomass. Our regional‐scale estimations are key inputs for future robust biomass projections of energy sorghum genotypes integrating G×E×M under climate change scenarios.

Ojeda, Jonathan J.↗

NASA NTRS2018

Air-Sea Interactions in a High-Resolution Ocean-Atmosphere Simulation

During the past few years the Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have produced, respectively, global atmosphere-only and ocean-only simulations with km-scale grid spacing. These simulations have proved invaluable for process studies and for the development of satellite and in-situ sampling strategies. Nevertheless, a key limitation of these "nature" simulations is the lack of interaction between the ocean and the atmosphere, which limits their usefulness for studying air-sea interactions and for designing observing missions to study these interactions. We present here results from a coupled GEOS-MIT "nature run" simulation, wherein we have coupled a cubed-sphere-720 (~ 1/8) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere simulation to equivalent atmosphere-only and ocean-only simulations. A particular focus of the comparisons is the coupled versus uncoupled differences in interactions between Sea Surface Temperature (SST) and ocean surface wind. We discuss, in particular, a several-day mode of temporal variability in the SST-wind cycle and how it is represented in the different model simulations and in observationally-based products. A mechanism for the cycle, which is driven by SST-wind feedback, is proposed.

Strobach, Ehud↗

NASA NTRS2019

An Ocean-Atmosphere Simulation for Studying Air-Sea Interactions

During the past few years the Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have produced, respectively, global atmosphere-only and ocean-only simulations with km-scale grid spacing. These simulations have proved invaluable for process studies and for the development of satellite and in-situ sampling strategies. Nevertheless, a key limitation of these "nature" simulations is the lack of interaction between the ocean and the atmosphere, which limits their usefulness for studying air-sea interactions and for designing observing missions to study these interactions. To remove this limitation, we aim to perform a coupled simulation using the km-scale GEOS atmosphere and the km-scale MIT ocean models. The initial attempt at the km-scale coupled simulation resulted in computational issues which will be presented here. As a preliminary step towards the km-scale objective, we present results from a high resolution but not yet km-scale simulation, wherein we have coupled a cubed-sphere-720 (~ 1/8) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere set-up to equivalent atmosphere-only and ocean-only simulations. A particular focus of the comparisons is the differences in interactions between Sea Surface Temperature (SST) and ocean surface wind for the coupled and uncoupled simulations. We discuss observed and modeled high temporal variability (~days) SST-wind cycle and how it is represented in the different systems. A mechanism for the cycle, which is driven by SST-wind feedback, is proposed.

Strobach, Ehud↗

NASA NTRS2019

An Ocean-Atmosphere Simulation for Studying Air-Sea Interactions

Strobach, Ehud↗