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 109 records · Page 6

Adaptive Neuro Fuzzy Inference System for Cyber-Intrusion Detection in a Smart Grid

The evolution of the power grid has brought increasing deployment of advance metering infrastructure, penetration of intelligent electronic devices, and integration of physical power system components with information and communications technologies. With the fast-expanding connectivity, cyber vulnerabilities arise due to the use of internet-based communication systems. These systems are targets of cyber-intrusions which attempt to disturb the normal power system functions. Traditional intrusion detection algorithms have been developed without an explicit model of the cyber components. In this paper, an algorithm to detect false data injections in the power system is proposed considering both cyber and physical models of the power system. The algorithm is based on an Adaptive Neuro Fuzzy Inference System (ANFIS) which collects information from state variables of the cyber-physical system to meet the performance requirements of the grid. Simulations of the proposed approach using the IEEE 13-bus test system validate the effectiveness of this artificial intelligence-based algorithm.

Adaptive Neuro Fuzzy Inference System (ANFIS), cyb↗

Provision of Grid Services by PV Plants with Integrated Battery Energy Storage System: Preprint

Battery energy storage systems (BESS), due to their tremendous range of uses and configurations, may assist PV integration in any number of ways by increasing power system flexibility. In this paper we describe results of a research project conducted by NREL and First Solar for developing controls and demonstrating many use cases for PV-BESS systems including (1) matching generation to loads through time shifting, reducing PV curtailment,; (2) promoting higher levels of PV penetrations by balancing the grid through ancillary services; (3) use of PV-BESS systems to provide wide-area stability services in the form of oscillations damping controls; and (r) use of gird forming BESS for enabling black start and islanded applications for PV-BESS systems. In the paper we will also describe the unique multi-MW test and validation platform developed at NREL for conducting research on fielded utility-scale PV-BESS systems. The platform consists of First Solar 0.43 MW PV power plant, NREL’s 1MW/1MW BESS, 2MW PV power supplies and 7MVA power electronic grid simulator interfaced with industrial plant controllers and power-hardware-in-the-loop (PHIL) capability.

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

Frame invariant neural network closures for Kraichnan turbulence

Numerical simulations of geophysical and atmospheric flows have to rely on parameterizations of subgrid scale processes due to their limited spatial resolution. Despite substantial progress in developing parameterization (or closure) models for subgrid scale (SGS) processes using physical insights and mathematical approximations, they remain imperfect and can lead to inaccurate predictions. In recent years, machine learning has been successful in extracting complex patterns from high-resolution spatio-temporal data, leading to improved parameterization models, and ultimately better coarse grid prediction. However, the inability to satisfy known physics and poor generalization hinders the application of these models for real-world problems. In this work, we put forth a frame invariant closure approach to improve the accuracy and generalizability of deep learning-based subgrid scale closure models by embedding physical symmetries directly into the structure of the neural network. Specifically, we utilized specialized layers within the convolutional neural network in such a way that desired constraints are theoretically guaranteed without the need for any regularization terms. We demonstrate our framework for a two-dimensional decaying turbulence test case mostly characterized by the forward enstrophy cascade. We show that our frame invariant SGS model (i) accurately predicts the subgrid scale source term, (ii) respects the physical symmetries such as translation, Galilean, and rotation invariance, and (iii) is numerically stable when implemented in coarse-grid simulation with generalization to different initial conditions and Reynolds number. This work opens up a possibility of connecting physics-based theories and data-driven modeling paradigms, and thus represents a promising step towards the development of physically consistent data-driven turbulence closure models.

42 ENGINEERING↗

Grid Modeling and Hosting Capacity Analysis

An interface with the OpenDSS distribution grid simulator, facilitating users in calibrating and validating models. Additionally, three distinct tools have been developed: PV Hosting Capacity: This tool allows users to determine the additional amount of photovoltaic (PV) generation that can be integrated into the system without breaching operational constraints. EV Hosting Capacity: This tool focuses on identifying the capacity for incorporating extra electric vehicle (EV) load into the system without surpassing operational limitations. Project Impact Analysis Tool: This tool is designed to assess and report all potential grid violations associated with a specific project, providing valuable insights into its impact on the distribution grid.

Poudel, Shiva↗

Near real-time traffic routing

A near real-time physical transportation network routing system comprising: a traffic simulation computing grid and a dynamic traffic routing service computing grid. The traffic simulator produces traffic network travel time predictions for a physical transportation network using a traffic simulation model and common input data. The physical transportation network is divided into a multiple sections. Each section has a primary zone and a buffer zone. The traffic simulation computing grid includes multiple of traffic simulation computing nodes. The common input data includes static network characteristics, an origin-destination data table, dynamic traffic information data and historical traffic data. The dynamic traffic routing service computing grid includes multiple dynamic traffic routing computing nodes and generates traffic route(s) using the traffic network travel time predictions.

Yang, Chaowei↗

IM3 GO WEST Parameter Search Dataset

GO WEST is an open-source power grid modeling framework for U.S. Western Interconnection, which allows users to tailor the model depending on their research study and science questions. It is developed to address weather and water dynamics, and associated vulnerabilities in this bulk power system. It covers 28 balancing authorities (BA) and 12 states in U.S. Western Interconnection. GO WEST allows users to select different number of nodes and come up with a simplified network by utilizing 10,000 nodal topology of U.S. Western Interconnection created by Texas A&M University. Users can try and select different number of nodes, mathematical formulations (linear programming vs. mixed-integer linear programming), transmission line limit scaling factors, and hurdle rate scaling factors. GO WEST offers a unit commitment and economic dispatch (UC/ED) module to simulate grid operations on an hourly scale. In this sense, users can calibrate and validate their model versions by comparing model outputs to historical datasets. Therefore, GO WEST can help researchers to strike a balance between model fidelity (i.e. accuracy) and computational complexity (i.e. runtime). This dataset includes model inputs and outputs from 600 model versions for each 2019, 2020, and 2021. The folder naming convention is as follows: Exp{Number of Nodes}_{Mathematical Formulation}_{Transmission Line Limit Scaling Factor in MW}_{Hurdle Rate Scaling Factor in %}_{Year}. Linear programming is designated with "simple" label whereas mixed-integer linear programming is designated with "coal" label. For example, "Exp100_simple_1000_50_2019" folder contains inputs and outputs from 2019 model version with 100 nodes, linear programming, +1000 MW transmission line limit scaling factor, and +50% hurdle rate scaling factor. GO WEST GitHub repository hosts all raw datasets, processing scripts, and model scripts. Please refer to the README file for a detailed description of the included files.

Economics↗

A Methodology for Evaluating Operator Usage of Machine Learning Recommendations for Power Grid Contingency Analysis

This work presents the application of a methodology to measure domain expert trust and workload, elicit feedback, and understand the technological usability and impact when a machine learning assistant is introduced into contingency analysis for real-time power grid simulation. The goal of this framework is to rapidly collect and analyze a broad variety of human factors data in order to accelerate the development and evaluation loop for deploying machine learning applications. We describe our methodology and analysis, and we discuss insights gained from a pilot participant about the current usability state of an early technology readiness level (TRL) artificial neural network (ANN) recommender.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗

Novel Representation of Leaf Phenology Improves Simulation of Amazonian Evergreen Forest Photosynthesis in a Land Surface Model

Leaf phenology in the humid tropics largely regulates the seasonality of forest carbon and water exchange. However, it is inadequately represented in most global land surface models due to limited understanding of its controls. Based on intensive field studies at four Amazonian evergreen forests, we propose a novel, quantitative representation of tropical forest leaf phenology, which links multiple environmental variables with the seasonality of new leaf production and old leaf litterfall. The new phenology simulates higher rates of leaf turnover (new leaves replacing old leaves) in dry seasons with more sunlight, which is then implemented in ORCHIDEE, together with recent findings of ontogeny-associated photosynthetic capacity, and is evaluated against ground-based measurements of leaf phenology (canopy leaf area index and litterfall), eddy covariance fluxes (photosynthesis and latent heat), and carbon allocations from field observations. Results show the periodical cycles of solar radiation and vapor pressure deficit are the two most important environmental variables that are empirically related to new leaf production and old leaf abscission in tropical evergreen forests. The model with new representation of leaf phenology captures the seasonality of canopy photosynthesis at three out of four sites, as well as the seasonality of litterfall, latent heat, and light use efficiency of photosynthesis at all tested sites, and improves the seasonality of carbon allocations to leaves, roots, and sapwoods. This study advances understanding of the environmental controls on tropical leaf phenology and offers an improved modeling tool for gridded simulations of interannual CO 2 and water fluxes in the tropics.

54 ENVIRONMENTAL SCIENCES↗

Higher Order Chemistry Models in the CFD Simulation of Laser-Ablated Carbon Plumes

Production of single-walled carbon nanotubes (SWNT) has taken place for a number of years and by a variety of methods such-as laser ablation, chemical vapor deposition, and arc-jet ablation. Yet, little is actually understood about the exact chemical kinetics and processes that occur in SWNT formation. In recent time, NASA Johnson Space Center has devoted a considerable effort to the experimental evaluation of the laser ablation production process for SWNT originally developed at Rice University. To fully understand the nature of the laser ablation process it is necessary to understand the development of the carbon plume dynamics within the laser ablation oven. The present work is a continuation of previous studies into the efforts to model plume dynamics using computational fluid dynamics (CFD). The ultimate goal of the work is to improve understanding of the laser ablation process, and through that improved understanding, refine the laser ablation production of SWNT. Fig. 1 shows a basic schematic of the laser-ablation oven at NASA-JSC. Construction of the facility is simple in concept. Two concentric quartz tubes of 1.5 mm thickness form the inner and outer tubes with inside diameters of 2.2 and 5.08 cm respectively. At one end of the inner tube are located two 60 Hz pulsed lasers operating at 1064 nm and 532 nm wavelength with beam diameters of 5 mm aligned coaxially with the longitudinal axis of the inner quartz tube. For standard nanotube production runs, a 10 ns 532 nm pulse is followed 50 ns later by a 10 ns 1064 nm pulse. Each pulse is of 300 mJ energy. A target of carbon graphite with approximately 1% nickel and cobalt catalysts is located at the other end of the inner quartz tube. In the ordinary processing of SWNT, a base flow of 100 sccm of argon is maintained from the laser location and exits past the carbon target at a pressure of 66.7 kPa. These conditions yield a baseline mass flow through the chamber of 2.723x10(exp -6)kg/s of argon. The whole oven facility is heated to a temperature of 1473 K prior to nanotube production runs. Upon laser irradiation, part of the carbon target ablates immediately and forms a carbon vapor plume that penetrates into the argon base flow towards the laser initially at supersonic velocities. In the previous studies either a single carbon species, C3, was used to model the plume development, or a simplified 11 species carbon model reduced down from the Krestinin and Moravsky model for full fullerene chemistry was used. While both of these options yielded significant results, it was felt that the actual chemistry occ&g in the carbon plume might have a greater affect on the plume than assumed. Indeed, in the earlier attempts at modeling the carbon plume, several thermophysical characteristics could never be matched to experimental observations of plume development - mainly the propagation distance of the plume itself. In the present study, two additional chemistry models will be used to duplicate the previous studies simulations of the carbon plume. The first chemistry model used in this study is again a reduced form of the Krestinin and Moravsky rates. However the highest order carbon species allowed has been increased from C6 to C30 - therefore allowing the simulation of up to a half of the standard C60 fullerene. The second chemistry model investigated is a reduced form of a full carbon nanotube model developed at NASA-JSC. The C30 studies have already been accomplished at the present time, and the reduced SWNT model studies are currently underway. To pursue the current study, one sacrifice had to be made in that the simulation grid spacing had to be increased from 0.5 mm spacing to 1 mm spacing for the sake of computational efficiency since computational effort is proportional to the square of the number of grid points multiplied by the number of species considered. propagation that is far more in line with the experimental results observed by Puretzky et al as shown in Fig. 2. e C6 studies had yielded a far greater propagation in previous studies. In addition, chemical species development with the C30 model indicates that many higher order carbon species are produced outside of the plume proper (indicated by plotting contours of the background argon concentrations in Fig. 3) - this result was not observed in previous studies. In fact, some species primarily occurred outside of the plume itself - as shown for C27 in Fig. 4 when it is compared to Fig. 3. It could be asserted that this has occurred because all of the C27 in the plume had already been consumed in the formation of C30, but this does not seem to be indicated over time. Several other factors that arose in the previous studies have also been made more clear by the use of the higher order chemical models - one being that the use of c6 as an indicator species was mistaken. C6 is the only carbon species in the previous studies that was not injected into the flowfield as a boundary condition; it was therefore hoped that this species would provide insight into the formation of higher order carbon species for comparison to full SWNT production. But, when the plot of total mass in the plume is examined on a species by species basis in Fig. 5, it is seen that Cg was a fairly insignificant contributor to the total carbon mass in the plume and would not provide information on higher order carbon formation. the thermophysical characteristics of the carbon plume as well as simulate the carbon plume using the reduced SWNT model to provide an even better simulation of full chemistry effects upon plume propagation.

Scott, C. D.↗

Accelerating geostatistical modeling using geostatistics-informed machine Learning

Ordinary Kriging (OK) is a popular geostatistical algorithm for spatial interpolation and estimation. The computational complexity of OK changes quadratically and cubically for memory and speed, respectively, given the number of data. Therefore, it is computationally intensive and also challenging to process a large set of data, especially in three-dimensional (3D) cases. This paper develops a geostatistics-informed machine learning (GIML) model to improve the efficiency of OK by reducing the number of points required to be estimated using OK. Specifically, only a very few of the unknown points are estimated by OK to get the weights and estimations, which are used as the training dataset. Moreover, the governing equations of OK are used to guide our proposed machine learning to better reproduce the spatial distributions. Our results show that the proposed GIML can reduce the computational time of OK by at least one order of magnitude. The effectiveness of the GIML is evaluated and compared using a 2D case. Furthermore, we demonstrate its efficiency and robustness by considering a different number of training samples on various 3D simulation grids.

58 GEOSCIENCES↗

Wave resource characterization at regional and nearshore scales for the U.S. Alaska coast based on a 32-year high-resolution hindcast

A wave resource characterization was performed for the southern coast of Alaska based on a 32-year hindcast covering the period from 1979 to 2010. The characterization closely followed International Electrotechnical Commission (IEC) Technical Specifications. An unstructured-grid Simulating WAves Nearshore (SWAN) model, which had an approximate spatial resolution of 300 m within 30 km from the nearest shoreline, was developed. Additionally, extensive model validation and error characterization was performed. The model was found to perform well with an average absolute percent error of 8.6% in significant wave height, averaged over 18 buoys. Statistics for the six IEC wave resource parameters were calculated and aggregated at 20 km from shore to quantify the incident wave power and its variability at a regional scale. The southern coast of the Aleutian Archipelago was found to have the most available wave energy in the region. A nearshore resource assessment was performed by evaluating resource hotspots located 1 km from shore. More than 900 nearshore stations had an Optimum Hotspot Identifier value of 5 (kW/m) at diverse water depths, thereby positioning Alaska as a promising location for wave energy development.

54 ENVIRONMENTAL SCIENCES↗

Source Contributions to Fine Particulate Matter and Attributable Mortality in India and the Surrounding Region

Fine particulate matter (PM 2.5 ) is a leading mortality risk factor in India and the surrounding region of South Asia. This study evaluates the contribution of emission sectors and fuels to PM 2.5 mass for 29 states in India and 6 surrounding countries (Pakistan, Bangladesh, Nepal, Bhutan, Sri Lanka, and Myanmar) by combining source-specific emission estimates, stretched grid simulations from a chemical transport model, high resolution hybrid PM 2.5 , and disease-specific mortality estimates. We find that 1.02 (95% Confidence Interval (CI): 0.78-1.26) million deaths in South Asia attributable to ambient PM 2.5 in 2019 were primarily from three leading sectors: residential combustion (28%), industry (15%), and power generation (12%). Solid biofuel is the leading combustible fuel contributing to the PM 2.5 -attributable mortality (31%), followed by coal (17%), and oil and gas (14%). State-level analyses reveal higher residential combustion contributions (35%-39%) in states (Delhi, Uttar-Pradesh, Haryana) with high ambient PM 2.5 (>95 µg/m 3 ). The combined mortality burden associated with residential combustion (ambient) and household air pollution (HAP) in India is 0.72 million (95% CI:0.54-0.89) (68% attributable to HAP, 32% attributable to residential combustion). Our results illustrate the potential to reduce PM 2.5 mass and improve population health by reducing emissions from traditional energy sources across multiple sources in South Asia.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Selective bond formation triggered by short optical pulses: quantum dynamics of a four-center ring closure

Here we report bond formation induced by an ultrashort UV pulse. The photochemical process is described by quantum dynamics as coherent electronic and nuclear motions during the ultrashort pulse induced ring closure of norbornadiene to quadricyclane. Norbornadiene consists of two ethylene moieties connected by a rigid (CH 2 ) 3 bridge. Upon photoexcitation, two new sigma bonds are formed, resulting in the closure of a four-atom ring. As a medium-sized polyatomic molecule, norbornadiene exhibits a high density of strongly coupled electronic states from about 6 eV above the ground state. We report on inducing the formation of the new bonds using a short femtosecond UV pulse to pump a non-equilibrium electronic density in the open form that evolves towards the closed ring form. As the coherent electronic-nuclear coupled dynamics unfold, the excited states change character through non-adiabatic interactions and become valence states for the two new C–C bonds of quadricyclane. Our three-dimensional fully quantum dynamical grid simulations during the first 200 fs show that short UV pulses of different polarization initiate markedly different initial non-equilibrium electronic densities that follow different dynamical paths to the S 0 /S 1 conical intersection. They lead to different initial relative yields of quadricyclane, thereby opening the way to controlling bond-making with attopulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Joint Matrix Completion and Compressed Sensing for State Estimation in Low-observable Distribution System

Limited measurement availability at the distribution grid presents challenges for state estimation and situational awareness. This paper combines the advantages of two sparsity-based state estimation approaches (matrix completion and compressive sensing) that have been proposed recently to address the challenge of unobservability. The proposed approach exploits both the low rank structure and a suitable transform domain representation to leverage the correlation structure of the spatio-temporal data matrix while incorporating the powerflow constraints of the distribution grid. Simulations are carried out on three phase unbalanced IEEE 37 test system to verify the effectiveness of the proposed approach. The performance results reveal - (1) the superiority over traditional matrix completion and (2) very low state estimation errors for high compression ratios representing very low observability.

Dahale, Shweta↗

Cyber Framework for Steering and Measurements Collection Over Instrument-Computing Ecosystems

We propose a framework to develop cyber solutions to support the remote steering of science instruments and measurements collection over instrument-computing ecosystems. It is based on provisioning separate data and control connections at the network level, and developing software modules consisting of Python wrappers for instrument commands and Pyro server-client codes that make them available across the ecosystem network. We demonstrate automated measurement transfers and remote steering operations in a microscopy use case for materials research over an ecosystem of Nion microscopes and computing platforms connected over site networks. The proposed framework is currently under further refinement and being adopted to science workflows with automated remote experiments steering for autonomous chemistry laboratories and smart energy grid simulations.

Al Najjar, Anees↗

Sequence Impedance Measurement of Utility-Scale Wind Turbines and Inverters – Reference Frame, Frequency Coupling, and MIMO/SISO Forms

Sequence impedance responses with or without considering frequency coupling in both MIMO and SISO forms are increasingly used for the stability analysis of three-phase power electronic systems; however, many aspects of sequence impedance measurement are not fully explored. It is not clear if the sequence impedance has a reference frame similar to the dq impedance. If so, the role of the grid voltage angle estimation in aligning the sequence impedance reference frame has not been discussed. Additionally, existing methods for measuring the sequence impedance with frequency coupling are complicated, are not feasible for large wind turbines and inverters, and provide the sequence impedance responses in only either MIMO or SISO form. This paper presents a sequence impedance measurement method that considers the frequency coupling, performs reference frame alignment, demonstrates the impact of the grid voltage angle estimation, and obtains the sequence impedance response in both MIMO and SISO forms. This paper demonstrates the proposed method and practical problems associated with the sequence impedance measurement of utility-scale wind turbines and inverters on a 1.9-MW Type III wind turbine and a 2.2-MVA inverter using an impedance measurement system built around a 7-MW/13.8-kV grid simulator and a 5-MW dynamometer.

17 WIND ENERGY↗

Electrification Futures Study: Operational Analysis of U.S. Power Systems with Increased Electrification and Demand-Side Flexibility

This report provides a high-resolution U.S. national-scale power systems assessment of scenarios with widespread electrification using detailed grid simulations. We find that high electrification scenarios envisioned in the Electrification Futures Study can be operated at the hourly level to meet future increased levels of electrified demand and, under the most aggressive cases, with up to 66% of generation from wind and solar technologies. The analysis also examines the value of flexibility from demand by supporting system operability under high electrification and high renewable systems. Demand-side flexibility—especially flexibility from newly electrified loads—is found to enhance the operational efficiency by reducing VRE curtailment and increasing utilization of generators with lower operating costs. It can also enhance the ability of electrification to decarbonize the energy sector.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗