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 145 records · Page 8

New Developments and Capabilities Within WEC-Sim: Preprint

WEC-Sim is an open-source software for simulating wave energy converters, which has been actively developed and applied since its initial release in 2014 to simulate a wide variety of device archetypes. WEC-Sim is developed jointly by the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (Sandia) within the MATLAB/SIMULINK environment. A general wave-to-wire model begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of wave energy converters (WEC), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modeling the entire wave-to-wire is difficult and encompasses multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze, and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM)-based frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). With this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Additional WEC-Sim functionalities include pre-built Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems. The developers of WEC-Sim continue to release new versions of the software, at least annually, with our latest release in September 2022. These releases include bug fixes, updates to software documentation, as well as new features to expand WEC-Sim's capabilities to model a wide range of WEC concepts. This publication will highlight the new features added to WEC-Sim between versions 4.1.0 to 5.0.1 which spans over a two year period from June 2020 to September 2022. New features to be described will include topics such as continuous integration checks, revised Morison Element and nonlinear hydro implementations, run directly from Simulink (required for hardware-in-the-loop execution), BEMIO updates to import Capytaine BEM hydrodynamics, addition of cable blocks, and new wave visualization features.

TIDAL AND WAVE POWER↗

Finite time step and spatial grid effects in δf simulation of warm plasmas

This paper introduces a technique for analyzing time integration methods used with the particle weight equations in δf method particle-in-cell (PIC) schemes. The analysis applies to the simulation of warm, uniform, periodic or infinite plasmas in the linear regime and considers the collective behavior similar to the analysis performed by Langdon for full-f PIC schemes [1,2]. We perform both a time integration analysis and spatial grid analysis for a kinetic ion, adiabatic electron model of ion acoustic waves. An implicit time integration scheme is studied in detail for δf simulations using our weight equation analysis and for full-f simulations using the method of Langdon. It is found that the δf method exhibits a CFL-like stability condition for low temperature ions, which is independent of the parameter characterizing the implicitness of the scheme. The accuracy of the real frequency and damping rate due to the discrete time and spatial schemes is also derived using a perturbative method. The theoretical analysis of numerical error presented here may be useful for the verification of simulations and for providing intuition for the design of new implicit time integration schemes for the δf method, as well as understanding differences between δf and full-f approaches to plasma simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hardware-in-the-Loop Co-simulation of Distribution Grid for Demand Response

In modern power systems, co-simulation is proposed as an enabler for analyzing the interactions between disparate systems. This paper introduces the co-simulation platform Virtual Grid Integration Laboratory (VirGIL) including Hardware-in-the-Loop testing, and demonstrates its potential to assess demand response strategies. VirGIL is based on a modular architecture using the Functional Mock-up Interface industrial standard to integrate new simulators. VirGIL combines state-of-the-art simulators in power systems, communications, buildings, and control. In this work, VirGIL is extended with a Hardware-in-the-Loop component to control the ventilation system of a real 12-story building in Denmark. VirGIL capabilities are illustrated in three scenarios: load following, primary reserves and load following aggregation. Experimental results show that the system can track one minute changing signals and it can provide primary reserves for up-regulation. Furthermore, the potential of aggregating several ventilation systems is evaluated considering the impact at distribution grid level and the communications protocol effect.

Rotger-Griful, Sergi↗

Recent Developments in the WEC-Sim Open-Source Design Tool: Preprint

WEC-Sim (Wave Energy Converter SIMulator) is an open-source code for simulating wave energy converters, which has been actively developed and applied to simulate a wide variety of device archetypes and has become a popular tool since its initial release in 2014. WEC-Sim is developed jointly by the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (SNL) within the MATLAB/SIMULINK environment. Figure 1 illustrates a general wave-to-wire model which begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of a single WEC (or array), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modelling the entire wave-to-wire is difficult and encompass multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM) based-frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). Within this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Table 1 lists additional WEC-Sim functionalities, which are created using prebuilt Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems.

hydrodynamics modeling↗

ParNCL and ParGAL: Data-parallel Tools for Postprocessing of Large-scale Earth Science Data

Earth science high-performance applications often require extensive analysis of their output in order to complete the scien- tific goals or produce a visual image or animation. Often this analysis cannot be done in situ because it requires calculating time-series statistics from state sampled over the entire length of the run or analyzing the relationship between similar time series from previous simulations or observations. Many of the tools used for this postprocessing are not themselves high- performance applications, but the new Parallel Gridded Analysis Library (ParGAL) provides high-performance data-parallel versions of several common analysis algorithms for data from a structured or unstructured grid simulation. The library builds on several scalable systems, including the Mesh Oriented DataBase (MOAB), a library for representing mesh data that sup- ports structured, unstructured finite element, and polyhedral grids; the Parallel-NetCDF (PNetCDF) library; and Intrepid, an extensible library for computing operators (such as gradient, curl, and divergence) acting on discretized fields. We have used ParGAL to implement a parallel version of the NCAR Command Language (NCL) a scripting language widely used in the climate community for analysis and visualization. The data-parallel algorithms in ParGAL/ParNCL are both higher performing and more flexible than their serial counterparts.

97 MATHEMATICS AND COMPUTING↗

Hardware-in-the-Loop (HIL) Simulations for Smart Grid Impact Studies

Hardware-in-the-loop (HIL) simulations are increasingly employed in power engineering as more inverter-based generation and smart appliances are connected to the electric grid. HIL techniques allow for the co-simulation of analytical models with actual devices whose complex behavior is computationally inefficient or difficult to model. They also allow for testing the behavior of these devices under adverse conditions that rarely occur in the field but are important to evaluate. This paper provides an overview of how HIL simulations have been used to date and proposes that HIL simulations should play an important role in evaluating new control strategies, especially at the distribution level, that are being proposed to allow for continued affordable and reliable operation of the electric grid. We present a new capability that was developed to evaluate the interactions between residential loads and the smart grid: smart home hardware-in-the-loop. The paper includes results from an HIL experiment that incorporates multiple technologies and controls.

24 POWER TRANSMISSION AND DISTRIBUTION↗

PV Inverter Testing for Momentary Cessation and Rate-of-Change-of-Frequency Events

To understand the power system stability and develop better electromagnetic transient (EMT) models of field deployed photovoltaic (PV) inverters, it is important to characterize inverters' response to abnormal voltage and frequency scenarios. Because EMT models are not typically available for small distribution-connected PV inverters, and because inerterconnection standards historically did not specify desired ride-through behaviors, we tested two such inverters in the lab to characterize their responses to severe undervoltage events and high rate-of-change-of-frequency (ROCOF) conditions. The inverters tested were pre-IEEE 1547-2018 residential PV inverters widely used in the Hawaiian Electric territory and many other areas. The testing results for undervoltage scenarios showed that the inverter from one vendor exhibited momentary cessation while the inverter from the other vendor did not exhibit momentary cessation behavior or tripping for most of the events below the 120 ms undervoltage trip threshold duration set by IEEE 1547-2003. The testing results for ROCOF scenarios showed that the inverter from one vendor temporarily lost synchronization during ROCOF conditions while the inverter from the other vendor did not lose synchronization or cease generation for any ROCOF conditions. Both the inverters were also tested for EMT-simulated grid events with severe changes in frequency and voltage. The observed responses of the inverters were different from the simulated response of PV inverters represented the best available assumptions from pre-existing information. The results from these experiments can be used to update the inverter models used in bulk power system studies.

aggregates↗

PV Inverter Testing for Momentary Cessation and Rate-of-Change-of-Frequency Events: Preprint

To understand the power system stability and develop better electromagnetic transient (EMT) models of field deployed photovoltaic (PV) inverters, it is important to characterize inverters' response to abnormal voltage and frequency scenarios. Because EMT models are not typically available for small distribution-connected PV inverters, and because inerterconnection standards historically did not specify desired ride-through behaviors, we tested two such inverters in the lab to characterize their responses to severe undervoltage events and high rate-of-change-of-frequency (ROCOF) conditions. The inverters tested were pre-IEEE 1547-2018 residential PV inverters widely used in the Hawaiian Electric territory and many other areas. The testing results for undervoltage scenarios showed that the inverter from one vendor exhibited momentary cessation while the inverter from the other vendor did not exhibit momentary cessation behavior or tripping for most of the events below the 120 ms undervoltage trip threshold duration set by IEEE 1547-2003. The testing results for ROCOF scenarios showed that the inverter from one vendor temporarily lost synchronization during ROCOF conditions while the inverter from the other vendor did not lose synchronization or cease generation for any ROCOF conditions. Both the inverters were also tested for EMT-simulated grid events with severe changes in frequency and voltage. The observed responses of the inverters were different from the simulated response of PV inverters represented the best available assumptions from pre-existing information. The results from these experiments can be used to update the inverter models used in bulk power system studies.

IEEE 1547 2018↗

Hardware-in-the-Loop (HIL) Simulations for Smart Grid Impact Studies: Preprint

Hardware-in-the-loop (HIL) simulation is increasingly employed in power engineering as more inverter-based generation and smart appliances are connected to the electric grid. HIL techniques allow for co-simulation of analytical models with actual devices whose complex behavior is computationally inefficient or difficult to capture. It also allows for testing of the behavior of these devices under adverse conditions that only occur rarely in the field, but that are important to evaluate. An overview of how HIL simulation has been used to date is provided. This paper further proposes that HIL simulation should play an important role in evaluating new control strategies, especially at the distribution level, that are being proposed to allow for continued affordable and reliable operation of the electric grid with the introduction of these new technologies. It presents a new capability that was developed to enable evaluation of interactions between residential loads and the smart grid: smart home hardware-in-the-loop. The paper includes results from an HIL simulation that incorporates multiple technologies and controls.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Analysis of Transformer Magnetizing Reactance Using Describing Function in Direct Power Control of Back-to-Back Modular Multilevel Converter with Advanced Grid Support

This paper presents a Port Controlled Hamiltonian based direct power control architecture for a back-to-back modular multilevel converter system connecting two ac sources at different frequencies. The system features advanced grid support functionalities based on IEEE 1547-2018, implemented on the inverter side. The rectifier side controller ensures reference following for active and/or reactive powers and maintains the commanded dc bus voltage. The inverter side controller ensures power command following for active and reactive powers.The proposed control architecture is designed to suppress second harmonic oscillations in powers during unbalanced grid voltage sags by dynamically adjusting the currents on each ac side. This also ensures effective elimination of any second harmonic oscillations in the equivalent dc bus voltage. Validation is performed on an OPAL-RT real-time platform with case studies on unbalanced and balanced sags, demonstrating the controller's effectiveness during real time implementation. A reduced-scale laboratory prototype further verifies these case studies, with experimental results for balanced sags due to grid simulator limitations. The results confirm the robustness and efficiency of the proposed control strategy in ensuring stable and reliable operation under various grid conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗

EMT-TS Hybrid Simulation for Large Power Grids Considering IBR-Driven Dynamics

The escalating integration of inverter-based resources (IBRs) poses new challenges to power systems by introducing fast dynamics with higher frequencies, which may need to be simulated by an electromagnetic transient (EMT) program. As an alternative to conducting EMT simulations for the entire system, which is typically time consuming, hybrid simulation between EMT and phasor-domain transient stability (TS) can greatly reduce the computational burden while preserving the detailed fast dynamics in the EMT zone. This paper establishes an EMT-TS hybrid simulation platform using open-source tools, specifically ParaEMT, GridPACK, and HELICS, which are the EMT simulator, TS simulator, and interface framework, respectively. Case studies on the 240-bus Western Electricity Coordinating Council (WECC) system demonstrate that the developed ParaEMT-HELICS-GridPACK hybrid simulator can accurately capture both slow electromechanical and fast IBR-driven dynamics with a 2.4x speedup.

electromagnetic transient↗

Grid Impacts of Hyperloop Technology in the United States

This report provides an overview of Hyperloop technologies including the current technology developments and a brief discussion on key components of the Hyperloop system that determine the electric power requirements as a function of time. At-scale Hyperloop systems do not exist yet so a model was used to generate a set of load profiles for conceptual systems at four locations in the United States: two in California, one in Colorado, and one in Ohio. In all four cases, the modeled Hyperloop load profiles showed a pulsating characteristic in both active and reactive power. Grid modeling was performed by Pacific Northwest National Laboratory to estimate the grid impacts and to discuss grid integration challenges. This report discusses the outcomes of the grid simulations for four study cases and offers some perspective of how existing large industrial loads with similar load characteristics are treated to accommodate current grid planning guidelines. The report concludes with a discussion on how energy storage system might be used to eliminate the pulsating load characteristics or significantly reduce it.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Hestia-SWIFL: hourly anthropogenic fossil fuel CO2 and heat on the 2km WRF grid, version 1.1

The Hestia-SWIFL version 1.1 anthropogenic heat (AH) and fossil fuel CO2 (FFCO2) emissions data product represent emissions due to the combustion of fossil fuel and cement production within the state of Arizona from 2019 to 2022. This product was developed as part of the Southwest Urban Corridor Integrated Field Laboratory (SW-IFL) project, which aims to provide new knowledge and tools that address extreme heat, air quality, climate change and related urban environmental issues by integrating high-resolution observations, modeling, and civic engagement. The emissions are generated using a bottom-up/engineering approach and are tied to results generated by the Vulcan Project version 4, an effort to quantify space/time-resolved FFCO2 & AH emissions for the entire United States landscape. A large number of data sources are combined to best estimate the emissions at fine scales such as air quality emissions data, traffic flow data, building information, sociodemographic information, and fuel statistics. The AH product provides emissions for two emissions sources (transportation and point source emissions) in units of Watts per hour per square meter (W/m2) per year (annual files) or per hour (hourly files). The FFCO2 product provides emissions from nine individual emission sectors as well as the total, and in units of tons of carbon (tC) per grid cell per year or per hour. The output made available here places the native spatial resolution of the Hestia FFCO2 & AH emissions data product (points, lines, and polygons) into a regularized 2km x 2km grid at hourly and annual temporal resolutions, and stored in netCDF files. The exact spatial extent is defined by the ASU Weather Research Forecast (WRF) simulation grid. All data are processed using R/Python pm high-performance computing system. 2-27-2026 updates: Bugs in airport hourly profile (both AH and FFCO2) and building spatial patterns (FFCO2 only) were fixed. Hourly emissions are reprocessed for all years to reflect those changes.

54 ENVIRONMENTAL SCIENCES↗

Progress in Unsteady Turbopump Flow Simulations Using Overset Grid Systems

This viewgraph presentation provides information on unsteady flow simulations for the Second Generation RLV (Reusable Launch Vehicle) baseline turbopump. Three impeller rotations were simulated by using a 34.3 million grid points model. MPI/OpenMP hybrid parallelism and MLP shared memory parallelism has been implemented and benchmarked in INS3D, an incompressible Navier-Stokes solver. For RLV turbopump simulations a speed up of more than 30 times has been obtained. Moving boundary capability is obtained by using the DCF module. Scripting capability from CAD geometry to solution is developed. Unsteady flow simulations for advanced consortium impeller/diffuser by using a 39 million grid points model are currently underway. 1.2 impeller rotations are completed. The fluid/structure coupling is initiated.

Kiris, Cetin C.↗

VULCAN-CFD User Manual: Ver. 7.2.0

VULCAN-CFD offers a comprehensive set of capabilities to enable the simulation of continuum flowfields from subsonic to hypersonic conditions. The governing equations that are employed include allowances for both chemical and thermal nonequilibrium processes, coupled with a wide variety of turbulence models for both Reynolds-averaged and large eddy simulations. The software package can simulate two-dimensional, axisymmetric, or three-dimensional problems on structured multiblock meshes or fully unstructured meshes. A parabolic (i.e., space-marching) treatment can also be used for any subset of a structured mesh that can accommodate this solution strategy. The flow solver provides a significant level of geometric flexibility for structured grid simulations by allowing for arbitrary face-to-face C(0) continuous and non-C(0) continuous block interface connectivities. The unstructured grid paradigm allows for mixed element unstructured meshes that contain any combination of tetrahedral, prismatic, pyramidal, and hexahedral cell elements. The flow solver is also fully parallelized using MPI (Message Passing Interface) libraries in a data-parallel fashion, allowing for efficient simulations on modern High Performance Computing (HPC) systems. This document provides information related to the installation and execution of the VULCAN-CFD software package. A detailed description of the physical and numerical models available in the software are provided in the VULCAN-CFD Theory Manual.

VULCAN-CFD User Manual↗

Simulating Virtual Terminal Area Weather Data Bases for Use in the Wake Vortex Avoidance System (Wake VAS) Prediction Algorithm

During the research project, sounding datasets were generated for the region surrounding 9 major airports, including Dallas, TX, Boston, MA, New York, NY, Chicago, IL, St. Louis, MO, Atlanta, GA, Miami, FL, San Francico, CA, and Los Angeles, CA. The numerical simulation of winter and summer environments during which no instrument flight rule impact was occurring at these 9 terminals was performed using the most contemporary version of the Terminal Area PBL Prediction System (TAPPS) model nested from 36 km to 6 km to 1 km horizontal resolution and very detailed vertical resolution in the planetary boundary layer. The soundings from the 1 km model were archived at 30 minute time intervals for a 24 hour period and the vertical dependent variables as well as derived quantities, i.e., 3-dimensional wind components, temperatures, pressures, mixing ratios, turbulence kinetic energy and eddy dissipation rates were then interpolated to 5 m vertical resolution up to 1000 m elevation above ground level. After partial validation against field experiment datasets for Dallas as well as larger scale and much coarser resolution observations at the other 8 airports, these sounding datasets were sent to NASA for use in the Virtual Air Space and Modeling program. The application of these datasets being to determine representative airport weather environments to diagnose the response of simulated wake vortices to realistic atmospheric environments. These virtual datasets are based on large scale observed atmospheric initial conditions that are dynamically interpolated in space and time. The 1 km nested-grid simulated datasets providing a very coarse and highly smoothed representation of airport environment meteorological conditions. Details concerning the airport surface forcing are virtually absent from these simulated datasets although the observed background atmospheric processes have been compared to the simulated fields and the fields were found to accurately replicate the flows surrounding the airport where coarse verification data were available as well as where airport scale datasets were available.

Kaplan, Michael L.↗

Impedance Analysis and PHIL Demonstration of Reactive Power Oscillations in a Wind Power Plant Using a 4-MW Wind Turbine

This paper presents impedance-based analysis, mitigation, and power-hardware-in-the-loop (PHIL) demonstration of reactive power oscillations in a wind power plant using a 4-MW Type III wind turbine drivetrain. Because such low-frequency oscillations result from interactions among slower control loops of wind turbines regulating phasor quantities—active and reactive power output of the wind turbine and the magnitude of voltages at the point of interconnection (POI)—a new type of admittance is defined in terms of phasor quantities for their analysis. The so-called power-domain admittance of a wind turbine is defined as the transfer function from the frequency and magnitude of voltages at the POI to the active and reactive power output of the turbine. The power-domain admittance responses of the 4-MW wind turbine are measured using a 7-MVA grid simulator to identify the source of the reactive power oscillations. Power-domain impedance analysis and PHIL experiments are performed to explain how a resonant mode manifests as turbine-to-turbine and plant-to-grid reactive power oscillations. It is discovered that weaker grids exhibiting high inductive impedance mitigate oscillations in the reactive power output of wind power plants; however, a higher grid impedance does not help in damping turbine-to-turbine reactive power oscillations. This paper presents a simple droop-based solution to eliminate both turbine-to-turbine and plant-to-grid reactive power oscillations.

17 WIND ENERGY↗