Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Grid disturbance”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Cyber-Physical Anomaly Detection for Wide-Area Protection Using Machine Learning

Wide-area protection scheme (WAPS) provides system-wide protection by detecting and mitigating small and large-scale disturbances that are difficult to resolve using local protection schemes. As this protection scheme is evolving from a substation-based distributed remedial action scheme (DRAS) to the control center-based centralized RAS (CRAS), it presents severe challenges to their cybersecurity because of its heavy reliance on an insecure grid communication, and its compromise would lead to system failure. This article presents an architecture and methodology for developing a cyber-physical anomaly detection system (CPADS) that utilizes synchrophasor measurements and properties of network packets to detect data integrity and communication failure attacks on measurement and control signals in CRAS. The proposed machine leaning-based methodology applies a rules-based approach to select relevant input features, utilizes variational mode decomposition (VMD) and decision tree (DT) algorithms to develop multiple classification models, and performs final event identification using a rules-based decision logic. Here, we have evaluated the proposed methodology of CPADS using the IEEE 39 bus system for several performance measures (accuracy, recall, precision, and F-measure) in a cyber-physical testbed environment. Furthermore, our experimental results reveal that the proposed algorithm (VMD-DT) of CPADS outperforms the existing machine learning classifiers during noisy and noise-free measurements while incurring an acceptable processing overhead.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

An Analysis of Post Attack Impacts and Effects of Learning Parameters on Vulnerability Assessment of Power Grid

Due to the increasing number of heterogeneous devices connected to electric power grid, the attack surface increases the threat actors. Game theory and machine learning are being used to study the power system failures caused by external manipulation. Most of existing works in the literature focus on one-shot process of attacks and fail to show the dynamic evolution of the defense strategy. In this paper, we focus on an adversarial multistage sequential game between the adversaries of the smart electric power transmission and distribution system. We study the impact of exploration rate and convergence of the attack strategies (sequences of action that creates large scale blackout based on the system capacity) based on the reinforcement learning approach. We also illustrate how the learned attack actions disrupt the normal operation of the grid by creating transmission line outages, bus voltage violations, and generation loss. This simulation studies are conducted on IEEE 9 and 39 bus systems. The results show the improvement of the defense strategy through the learning process. The results also prove the feasibility of the learned attack actions by replicating the disturbances created in simulated power system.

attack impacts↗
NASA NTRS1975

Implicit approximate-factorization schemes for the low-frequency transonic equation

Two- and three-level implicit finite-difference algorithms for the low-frequency transonic small disturbance-equation are constructed using approximate factorization techniques. The schemes are unconditionally stable for the model linear problem. For nonlinear mixed flows, the schemes maintain stability by the use of conservatively switched difference operators for which stability is maintained only if shock propagation is restricted to be less than one spatial grid point per time step. The shock-capturing properties of the schemes were studied for various shock motions that might be encountered in problems of engineering interest. Computed results for a model airfoil problem that produces a flow field similar to that about a helicopter rotor in forward flight show the development of a shock wave and its subsequent propagation upstream off the front of the airfoil.

Ballhaus, W. F.↗
NASA NTRS1990

A model of stratospheric chemistry and transport on an isentropic surface

This paper presents a new photochemical transport model designed to simulate the behavior of stratospheric trace species in the middle stratosphere. The model has an Eulerian grid with the latitude and longitude coordinates on a single isentropic surface (hemispheric or global), in which both the dynamical and the photochemical processes can be accurately represented. The model is intgegrated for 12 days with winds and temperatures supplied by three-dimensional integration of an idealized wavenumber-one disturbance. The results for the long-lived tracers such as N2O showed excellent correlation with the potential vorticity distribution, validating the transport scheme. Calculations with zonally averaged wind and temperature fields showed that discrepancies in the calculation of the zonal mean were less than 10 percent for O3 and HNO3, compared with the zonal mean of the previous results.

Austin, John↗
NASA NTRS1993

Ionospheric plasma flow about a system of electrically biased flat plates

The steady state interaction of two electrically biased parallel plates immersed in a flowing plasma characteristic of low earth orbit is studied numerically. Fluid equations are developed to describe the motion of the cold positively charged plasma ions, and are solved using finite-differences in two dimensions on a Cartesian grid. The behavior of the plasma electrons is assumed to be described by the Maxwell-Boltzmann distribution. Results are compared to an analytical and a particle simulation technique for a simplified flow geometry consisting of a single semi-infinite negatively biased plate. Comparison of the extent of the electrical disturbance into the flowing plasma and the magnitude of the current collected by the plate is very good. The interaction of two equally biased parallel plates is studied as a function of applied potential. The separation distance at which the current collected by either plate decreases by five and twenty percent is determined as a function of applied potential. The percent decreases were based on a non-interacting case. The decrease in overall current is caused by a decrease in ionic density in the region between the plates. As the separation between the plates decreases, the plates collect the ions at a faster rate than they are supplied to the middle region by the oncoming plasma flow. The docking of spacecraft in orbit is simulated by moving two plates of unequal potential toward one another in a quasi-static manner. One plate is held at a large negative potential while the other floats electrically in the resulting potential field. It is found that the floating plate does not charge continuously negative as it approaches the other more negatively biased plate. Instead, it charges more and then less negative as ionic current decreases and then increases respectively upon approach. When the two plates come into contact, it is expected that the electrically floating plate will charge rapidly negative to a potential near that of the other plate.

Herr, Joel L.↗
DOE OSTI2025

Cyberattack Detection and Mitigation on Central Volt‐VAr Using Circuit Law and Machine Learning

ABSTRACT In a distribution grid, voltage is maintained within a nominal range through a Volt‐VAr function that controls capacitor banks, reactive power of distributed energy resources (DER), and on‐load tap changers (OLTC). Availability of communications helps with the implementation of central Volt‐VAr control; however, it also opens the system to cyberattacks, causing voltage disturbances. Previous work has shown the adverse impacts of false data injection (FDI) on the central Volt‐VAr control; however, very few works have studied methods to detect and mitigate FDI on Volt‐VAr control. This paper addresses gaps in the detection and mitigation of FDI on the measurement packets of a central Volt‐VAr control. This work uses a two‐stage algorithm for cyberattack detection since the accuracy of a single‐stage machine learning (ML)–based detection method decreases while dealing with unseen data. The first stage is based on the verification of measurements against circuit laws, and the second stage utilizes a tree search algorithm and an ML method to detect the falsified data. This paper compares long short‐term memory (LSTM) and bidirectional LSTM (BiLSTM) as the employed ML algorithms. Finally, the mitigation algorithm replaces the falsified data with the estimated output of the ML algorithm. The effectiveness of the proposed method is tested for several cases using the IEEE 13‐bus test system in PSCAD software.

Beikbabaei, Milad [Bradley Department of Electrica↗
NASA NTRS1991

Reynolds number effects on supersonic asymmetrical flows over a cone at high angle of attack

The supersonic viscous flow over a 5-degree half-angle cone at an angle of attack of four times the cone half-angle is studied computationally using both the conical and the three-dimensional Navier-Stokes equations. The numerical solutions were obtained with an implicit, upwind-biased algorithm. Asymmetrical flowfields of the absolute-instability type are found using the conical-flow equations which agree with published results. However, the absolute instabilities of the originally symmetric flow found with the conical equations do not occur in the three-dimensional simulations, although spurious asymmetric three-dimensional flows for symmetric bodies arise if the grid resolution is insufficient in the nose region. The asymmetric flows computed with the three-dimensional equations are convective instabilities and are possible if the local Reynolds number exceeds a critical value and a fixed geometric asymmetry is imposed. A continuous range of asymmetries can be developed, depending on the size of the disturbance and the Reynolds number. As the Reynolds number is increased, the asymmetries demonstrate a bistable behavior at levels of side force consistent with those predicted using the conical equations. Below a certain critical Reynolds number, any flow asymmetries arising from a geometrical asymmetry are damped with increasing distance downstream from the geometrical asymmetry.

Thomas, J. L.↗
NASA NTRS1986

Applications of potential theory computations to transonic aeroelasticity

Unsteady aerodynamic and aeroelastic stability calculations based upon transonic small disturbance (TSD) potential theory are presented. Results from the two-dimensional XTRAN2L code and the three-dimensional XTRAN3S code are compared with experiment to demonstrate the ability of TSD codes to treat transonic effects. The necessity of nonisentropic corrections to transonic potential theory is demonstrated. Dynamic computational effects resulting from the choice of grid and boundary conditions are illustrated. Unsteady airloads for a number of parameter variations including airfoil shape and thickness, Mach number, frequency, and amplitude are given. Finally, samples of transonic aeroelastic calculations are given. A key observation is the extent to which unsteady transonic airloads calculated by inviscid potential theory may be treated in a locally linear manner.

Edwards, J. W.↗
DOE OSTI2025

An Open-Access Repository of Synchrophasor Data Quality Examples: Curation and Example Applications

Synchrophasor measurements are critical in providing wide-area situational awareness to power system operators. However, data artifacts may be introduced due to various issues such as loss of communication, loss of GPS signal, internal clock error, and vendor-specific implementation of phasor estimation algorithms. Tools designed to provide actionable insights from synchrophasor data, hence, must be designed to be robust to these data quality issues. In this work, two years of synchrophasor data sourced from multiple electric utilities in the United States were analyzed to identify examples of data quality problems. These examples were then labeled and published in the Grid Event Signature Library, a publicly available repository of power system measurements hosted by the Oak Ridge National Laboratory. This paper describes the data curation process, and illustrates two application use cases where the dataset can be valuable to the research community. In the first use case, a random forest classifier is trained to distinguish power system disturbance signatures from data anomalies introduced in synchrophasor measurements due to clock errors. The second use case studies the impact of data quality issues on an example synchrophasor application (specifically, event start time determination). The choice of data quality problems investigated is informed by the examples in the repository curated in this work.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2001

Development of a Hybrid RANS/LES Method for Turbulent Mixing Layers

Significant research has been underway for several years in NASA Glenn Research Center's nozzle branch to develop advanced computational methods for simulating turbulent flows in exhaust nozzles. The primary efforts of this research have concentrated on improving our ability to calculate the turbulent mixing layers that dominate flows both in the exhaust systems of modern-day aircraft and in those of hypersonic vehicles under development. As part of these efforts, a hybrid numerical method was recently developed to simulate such turbulent mixing layers. The method developed here is intended for configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. Interest in Large Eddy Simulation (LES) methods have increased in recent years, but applying an LES method to calculate the wide range of turbulent scales from small eddies in the wall-bounded regions to large eddies in the mixing region is not yet possible with current computers. As a result, the hybrid method developed here uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall-bounded regions entering a mixing section and uses a LES procedure to calculate the mixing-dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. With this technique, closure for the RANS equations is obtained by using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The LES equations are closed using the Smagorinsky subgrid scale model. Although the function of the Cebeci-Smith model to replace all of the turbulent stresses is quite different from that of the Smagorinsky subgrid model, which only replaces the small subgrid turbulent stresses, both are eddy viscosity models and both are derived at least in part from mixing-length theory. The similar formulation of these two models enables the RANS and LES equations to be solved with a single solution scheme and computational grid. The hybrid RANS-LES method has been applied to a benchmark compressible mixing layer experiment in which two isolated supersonic streams, separated by a splitter plate, provide the flows to a constant-area mixing section. Although the configuration is largely two dimensional in nature, three-dimensional calculations were found to be necessary to enable disturbances to develop in three spatial directions and to transition to turbulence. The flow in the initial part of the mixing section consists of a periodic vortex shedding downstream of the splitter plate trailing edge. This organized vortex shedding then rapidly transitions to a turbulent structure, which is very similar to the flow development observed in the experiments. Although the qualitative nature of the large-scale turbulent development in the entire mixing section is captured well by the LES part of the current hybrid method, further efforts are planned to directly calculate a greater portion of the turbulence spectrum and to limit the subgrid scale modeling to only the very small scales. This will be accomplished by the use of higher accuracy solution schemes and more powerful computers, measured both in speed and memory capabilities.

Georgiadis, Nicholas J.↗
NASA NTRS1991

A finite-difference, frequency-domain numerical scheme for the solution of the linearized unsteady Euler equations

A numerical method is developed for solving periodic, three-dimensional, vortical flows around lifting airfoils in subsonic flow. The first-order method, that is presented, fully accounts for the distortion effects of the nonuniform mean flow on the convected upstream vortical disturbances. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Using an elliptic coordinate transformation, the unsteady boundary value problem is solved in the frequency domain on grids which are determined as a function of the Mach number and reduced frequency. Extensive comparisons are made with known solutions to unsteady vortical flow problems, and it is seen that the agreement is generally very good for reduced frequencies ranging from 0 up to 4.

Scott, James R.↗
NASA NTRS2012

Multi-Angle Implementation of Atmospheric Correction for MODIS (MAIAC). Part 3: Atmospheric Correction

This paper describes the atmospheric correction (AC) component of the Multi-Angle Implementation of Atmospheric Correction algorithm (MAIAC) which introduces a new way to compute parameters of the Ross-Thick Li-Sparse (RTLS) Bi-directional reflectance distribution function (BRDF), spectral surface albedo and bidirectional reflectance factors (BRF) from satellite measurements obtained by the Moderate Resolution Imaging Spectroradiometer (MODIS). MAIAC uses a time series and spatial analysis for cloud detection, aerosol retrievals and atmospheric correction. It implements a moving window of up to 16 days of MODIS data gridded to 1 km resolution in a selected projection. The RTLS parameters are computed directly by fitting the cloud-free MODIS top of atmosphere (TOA) reflectance data stored in the processing queue. The RTLS retrieval is applied when the land surface is stable or changes slowly. In case of rapid or large magnitude change (as for instance caused by disturbance), MAIAC follows the MODIS operational BRDF/albedo algorithm and uses a scaling approach where the BRDF shape is assumed stable but its magnitude is adjusted based on the latest single measurement. To assess the stability of the surface, MAIAC features a change detection algorithm which analyzes relative change of reflectance in the Red and NIR bands during the accumulation period. To adjust for the reflectance variability with the sun-observer geometry and allow comparison among different days (view geometries), the BRFs are normalized to the fixed view geometry using the RTLS model. An empirical analysis of MODIS data suggests that the RTLS inversion remains robust when the relative change of geometry-normalized reflectance stays below 15%. This first of two papers introduces the algorithm, a second, companion paper illustrates its potential by analyzing MODIS data over a tropical rainforest and assessing errors and uncertainties of MAIAC compared to conventional MODIS products.

Lyapustin, A.↗
DOE OSTI2020

Mitigating Cascading Outages in Severe Weather Using Simulation-Based Optimization

Severe weather events can trigger cascading power outages and lead to significant losses. In this work, we investigate cascading outage mitigation under severe weather conditions. Given day-ahead weather forecasts and component failure models, we aim to identify a set of power lines that can be hardened to minimize the expected impact of potential cascading outages. Since the expected load shedding cannot be expressed as an explicit function of line hardening decisions and system states, we developed a cascading outage simulator to estimate the expected value of load shedding under various initial weather-related disruption scenarios generated using a weather forecast. To avoid massive enumeration of all possible combinations of line hardening decisions and reduce the simulation efforts, we employed an efficient simulation-based optimization approach that quickly identifies the (near) optimal line hardening decisions in the presence of both large simulation noises due to the highly variable initial disturbances and system states, and significant randomness in the subsequent cascades. Furthermore, the algorithm is also able to utilize parallel computing to dramatically reduce computation time to support decision making in preparation for severe weather conditions. We performed a case study on the Northeast Power Coordinating Council (NPCC) 140-bus system model to demonstrate that our approach can significantly improve power grid resilience to adverse weather events.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2002

Verification Assessment of Flow Boundary Conditions for CFD Analysis of Supersonic Inlet Flows

Boundary conditions for subsonic inflow, bleed, and subsonic outflow as implemented into the WIND CFD code are assessed with respect to verification for steady and unsteady flows associated with supersonic inlets. Verification procedures include grid convergence studies and comparisons to analytical data. The objective is to examine errors, limitations, capabilities, and behavior of the boundary conditions. Computational studies were performed on configurations derived from a "parameterized" supersonic inlet. These include steady supersonic flows with normal and oblique shocks, steady subsonic flow in a diffuser, and unsteady flow with the propagation and reflection of an acoustic disturbance.

Slater, John W.↗
DOE OSTI2025

Large-signal Stability Analysis of Grid-forming Inverters with Equivalent-circuit Models

Here, this paper proposes an energy function-based direct method for large-signal stability assessment of grid-forming (GFM) inverters leveraging an equivalent-circuit representation of all involved control- and physical-layer dynamics. Three different primary controls, a standard inner-current outer-voltage cascaded-control architecture, output LCL filter, and reference-current saturation limiting are featured in the modeling and analysis framework. A composite energy function for the GFM inverter is obtained by summing up individual energy contributions gleaned from the circuit representation. The approach can readily be generalized to different primary controls, output-filter arrangements, and current limiters since it is based on a circuit-theoretic foundation. Numerical simulations validate the efficacy of the approach in estimating the critical clearing time following a large-signal disturbance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Root Cause Analysis of AC Overcurrent in July 2020 San Fernando Disturbance

On July 7 2020, approximately 1,000 MW reduction in solar photovoltaic (PV) output was experienced by the bulk power system in Southern California after a three-phase fault. Disturbance analysis indicates that instantaneous AC current caused inverter tripping. The projection of the root cause of AC overcurrent from the NERC report is that inverter current control is not tight. This letter provides a quantitative analysis on ac overcurrent and addresses the following questions: (i) what type of current controls may lead to ac overcurrent; and (ii) how controller parameters and grid strength influence ac overcurrent. Here, the letter provides frequency-domain analysis and time-domain simulation results of a grid-integrated inverter to illustrate the effect of control parameters and grid strength on ac overcurrent.

14 SOLAR ENERGY↗
DOE OSTI2025

Comparative Analysis of Inter-Area Oscillations in the US Eastern and Western Interconnections

This paper presents a comparative analysis of interarea oscillations in the US Eastern and Western Interconnections using frequency disturbance data collected from the advanced wide-area Frequency Monitoring Network (FNET/GridEye), enabling us to investigate and compare the oscillation characteristics of both regions. The study analyzes the statistical data from the two interconnections, including total oscillation events, average dominant frequencies, damping ratios, and maximum amplitudes. We also explore the impact of seasonal and daily variations on oscillation occurrences and the influence of different grid topologies and operational practices. The results provide insights into both interconnections' stability and control characteristics, offering valuable information for power system operators to enhance grid stability and oscillation suppression measures.

Fu, Hao [University of Tennessee, Knoxville (UTK)]↗
DOE OSTI2022

NGEE Arctic LANL Overview [Slides]

The Next-Generation Ecosystem Experiments (NGEE Arctic) project has a goal to "deliver a process-rich ecosystem model, extending from bedrock to the top of the vegetative canopy/atmospheric interface, in which the evolution of Arctic ecosystems in a changing climate can be modeled at the scale of a high resolution Earth System Model (ESM) grid cell." LANL works across multiple NGEE Arctic science questions, such as: Q1. How does the structure and organization of the landscape control permafrost evolution and associated carbon and nutrient fluxes in a changing climate? Q5. Where, when, and why will the Arctic become wetter or drier, and what are the implications for climate forcing? and Q6. What controls the vulnerability of Arctic ecosystems to disturbance, and how do disturbances alter the structure and function of these ecosystems? This report features LANL's research summaries.

54 ENVIRONMENTAL SCIENCES↗