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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects-Based Monitoring of Geomagnetically-Induced Current Using a Convolutional Neural Network

Geomagnetically-induced current (GIC) due to space weather can flow in the power grid causing undesirable effects such as transformer overheating, misoperation of protection devices, and potential blackouts. It is therefore important to monitor GIC in the power grid to improve online situational awareness and decision-making of system operators during a geomagnetic disturbance. To avoid the costly installation of GIC monitors at transformers’ neutrals, it is desirable to find correlations between GIC and already-monitored parameters. Hence, this work proposed the use of a convolutional neural network (CNN) to compute GIC amplitudes from learned patterns in the time-series data of transformer even harmonic currents. Using an electromagnetic transient program, GIC injection simulations were performed for a modeled Dominion Energy Virginia (DEV) substation with two 504 MVA, 500/230 kV transformers. Data collected from these offline simulations were used to train the CNN to provide online GIC monitoring. Testing the CNN performance involved using real GIC measurements from published literature and from a physical GIC monitor in the DEV area. Finally, the results showed that the proposed method was able to provide GIC readings with a root mean squared error of 1.56 A/phase (equivalent to an average accuracy of 94%) for these real GIC waveforms.

42 ENGINEERING↗
DOE OSTI2025

Unified Universal Control and Coordination of Inverter-Based Resources, and Validation for a PV + Battery Hybrid Plant

As renewable energy deployment grows, hybrid power plants (HPPs) combining photovoltaic (PV) and battery systems must evolve to offer both energy and grid stability services. These systems typically include a mix of grid-following (GFL) and grid-forming (GFM) inverters, presenting unique coordination and control challenges. This Department of Energy–funded project developed and validated a Unified Universal Control and Coordination (UUCC) framework for such PV + battery hybrid plants, enabling seamless and stable operation, including ultrafast black start, autonomous synchronization, and robust frequency and voltage regulation, under different grid conditions. The project significantly advanced the understanding of inverter-based resource (IBR) control by developing and validating three complementary system-level approaches for hybrid GFL/GFM operation: 1. A combined Virtual Resistance (VR)-based GFL and Virtual Oscillator Control (VOC)-based GFM method, where each inverter type is governed by a specialized control strategy. Together, these achieve stable, fast-response coordination, eliminating inrush current and enabling smooth black start and grid synchronization across a wide range of grid strengths. 2. A Deadbeat-based UUCC strategy, which uses discrete-time, switching-cycle-level control for both GFL and GFM inverters. This approach replaces traditional PI/PLL control with a control parameter-free, high-bandwidth framework that supports stable LVRT and instantaneous synchronization under all conditions. 3. A benchmark comparison with Siemens’ commercial GFM microgrid controller, which provided a fast baseline platform. The commercial approach decoupled v & f control was implemented on a commercial microgrid controller.The baseline commercial benchmark helped highlight superior transient response and black start performance offered by the deadbeat and VOC approaches. These technical contributions offer substantial improvements over conventional inverter control schemes, which often rely on slow phase-locked loop (PLL)-based synchronization, require careful control parameters tuning, and prone to unstable in weak grids with GFL inverters and in stiff grid with GFM inverters therefore challenging for hybrid GFL+GFM under all grid conditions. The deadbeat-based UUCC framework enables simpler, faster, and more robust operation of hybrid IBR systems using wide-bandgap (WBG) devices such as SiC power semiconductors. The rapid expansion of hybrid distributed energy resources (DERs), including residential and commercial PV-BESS installations such as Tesla Powerwall, PV with vehicle-to-grid (V2G) capability, and other integrated configurations, presents complex operational challenges for medium-voltage radial distribution feeders. These networks are subject to frequent disturbances such as faults, switching operations, rapid reclosing sequences, and feeder reconfigurations, all of which introduce dynamic stress on IBRs. In addition, planned feeder segmentation and deliberate islanding for resilience will require DERs that can autonomously perform blackstart, establish voltage and frequency references, and resynchronize with the main grid. The advanced deadbeat-based UUCC control and blackstart functionalities developed in this project directly address these requirements, enabling decentralized and autonomous operation of inverter-dominated DERs in distribution systems under a wide range of fault and reconfiguration scenarios. From a public benefit perspective, these innovations enable more reliable and cost-effective integration of renewable energy into distribution networks. The ability to autonomously black start and stabilize grids under varying grid conditions support accelerates recovery from outages and support decentralized resilient energy systems. By reducing system complexity and improving performance, this project lays critical groundwork for future inverter-dominated power grids that are clean, reliable, and accessible to all.

14 SOLAR ENERGY↗
NASA NTRS1985

Approximate-factorization schemes for solving the transonic full-potential equation

The present paper provides a general discussion of approximate-factorization techniques applied to the transonic full-potential equation. Giving particular attention to the AF2 approximate-factorization scheme. This scheme was first introduced by Ballhaus and Steger (1975) for solving the low-frequency (unsteady), transonic small-disturbance equation. The full-potential equation algorithm is examined, taking into account the governing equations, grid generation, the artificial density scheme (spatial differencing), the alternating direction implicit scheme, the AF2 iteration scheme, temporal damping, and boundary conditions. Computed results are also presented. It is shown that fast, fully-implicit algorithms of the approximate-factorization variety are both efficient and reliable for solving the conservative full-potential equation.

Holst, T. L.↗
DOE OSTI2022

Precise ROCOF estimation algorithm for low inertia power grids

The precise estimation of Rate Of Change Of Frequencies (ROCOFs) in a generation trip event can be helpful on power system inertia estimation, fast system response, and accurate event capturing. However, the ROCOF estimations from the existing Synchronized Measurement Devices (SMDs) are usually simply calculated by the finite difference between two adjacent frequency measurement points. The ineluctable noises, disturbances, and spikes from real-world frequency measurements can bring large dynamics to the ROCOF estimation and thus can result in an inaccurate estimation of the initial ROCOF. This issue becomes more serious when the target power grid has low inertia, and a large amount of distributed energy sources are deployed. To address this issue, a precise ROCOF estimation algorithm is designed based on least square method with flexible window size. In addition, a median filter is also designed and applied on the frequency measurements before using the proposed algorithm. The window size and thresholds in the proposed algorithm are determined with historical event data analysis. The proposed algorithm is deployed in a low cost, flexible, and distribution level universal grid analyzer (UGA) platform. In conclusion, multiple experiments are conducted in both a laboratory and the Hawaiian Islands to verify the effectiveness of the proposed algorithm.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Bulk Grid Frequency Support Using Electrolyzers

Increasing renewable penetration in the bulk power systems has resulted in declining inertial response of the generation during system disturbances. This results in larger frequency deviations and oscillations during such events. This paper proposes a strategy for controlled dispatch of electrolyzers as loads to improve the system’s frequency response. A droop based controller with constraints based on electrolyzer operating limits has been developed to regulate its load. Consecutively, a localized frequency based control strategy has been developed to support the bulk system during disturbances. The impact of the electrolyzer support has been demonstrated using the 23-bus SAVNW system, and the 240-bus WECC system for baseline, 25% and 50% renewable penetration. The impact of changing droop settings and the available controllable electrolyzer capacities on the system’s frequency response has also been presented. Overall, the dispatchable electrolyzers helps improve (reduce) both the maximum frequency deviation and the settling times during common disturbances like loss of load/generation, and line faults. Systems with lower inertia were seen to benefit more from the controlled electrolyzer dispatch.

39 EE - Hydrogen and Fuel Cell Technologies (EE-3F↗
DOE OSTI2025

Distributed Control and Optimization of Virtual Inertia in Power Systems

The integration of renewable energy sources into power grids through inverter-based resources (IBRs) reduces system inertia, leading to faster frequency dynamics and potential grid instability. This paper proposes a distributed, consensus-based approach for the real-time control and optimization of inertia sources during system disturbances, enhancing both system stability and economic performance. Furthermore, we integrate this approach with the traditional inertia placement problem, demonstrating its suitability for existing frameworks. We validate its effectiveness in a 4-bus test system and a two-area 8-bus test system, showcasing its potential for scalable, dynamic inertia management in low-inertia grids.

Yadav, Ajay [ORNL] (ORCID:000000016111881X)↗
NASA NTRS1976

The multigrid method: Fast relaxation

A multi-level grid method was studied as a possible means of accelerating convergence in relaxation calculations for transonic flows. The method employs a hierarchy of grids, ranging from very coarse (e.g. 4 x 2 mesh cells) to fine (e.g. 64 x 32); the coarser grids are used to diminish the magnitude of the smooth part of the residuals, hopefully with far less total work than would be required with optimal iterations on the finest grid. To date the method was applied quite successfully to the solution of the transonic small-disturbance equation for the velocity potential in conservation form. Nonlifting transonic flow past a parabolic arc airfoil is the example studied, with meshes of both constant and variable step size.

South, J. C., Jr.↗
NASA NTRS1982

A Numerical simulation of transition in plane channel flow

A numerical simulation of the final stages of transition to turbulence in plane channel flow at a Reynolds number of 7500 is described. Three dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time evolution of two and three dimensional finite amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32 by 33 by 32 grid. Solutions indicate the existence of structures similar to those observed in the laboratory and which are characteristic of various stages of transition that lead to final breakdown. Details of the resulting flow field after breakdown indicate the evolution of streak-like formations found in turbulent flows. Although the flow field does approach a steady state (turbulent channel flow), implementation of subgrid-scale terms are necessary to obtain proper turbulent statistics.

Goglia, G.↗
NASA NTRS1983

A numerical simulation of transition in plane channel flow

This paper involves a numerical simulation of the final stages of transition to turbulence in plane channel flow at a Reynolds number of 7500. Three-dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time-evolution of two- and three-dimensional finite-amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32 x 33 x 32 grid. Solutions indicate the existence of structures similar to those observed in the laboratory and which are characteristic of various stages of transition that lead to final breakdown. Details of the resulting flow field after breakdown indicate the evolution of streak-like formations found in turbulent flows. Although the flow field does approach a steady-state (turbulent channel flow), implementation of subgrid-scale terms are necesary to obtain proper turbulent statistics.

Biringen, S.↗
NASA NTRS1984

Final stages of transition to turbulence in plane channel flow

This paper involves a numerical simulation of the final stages of transition to turbulence in plane channel flow at a Reynolds number of 1500. Three-dimensional incompressible Navier-Stokes equations are numerically integrated to obtain the time evolution of two- and three-dimensional finite-amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32 x 51 x 32 grid. Solutions indicate the existence of structures similar to those observed in the laboratory and characteristics of the various stages of transition that lead to final breakdown. In particular, evidence points to the formation of a upside-down-V-shaped vortex and the subsequent system of horseshoe vortices inclined to the main flow direction as the primary elements of transition. Details of the resulting flow field after breakdown indicate the evolution of streaklike formations found in turbulent flows. Although the flow field does approach a steady state (turbulent channel flow), the introduction of subgrid-scale terms seems necessary to obtain fully developed turbulence statistics.

Biringen, S.↗
NASA NTRS1992

Direct simulation of high-speed mixing layers

A computational study of a nonreacting high-speed mixing layer is performed. A higher order algorithm with sufficient grid points is used to resolve all relevant scales. In all cases, a temporal free-stream disturbance is introduced. The resulting flow is time-sampled to generate a statistical cross section of the flow properties. The studies are conducted at two convective Mach numbers, three free-stream turbulence intensities, three Reynolds numbers, and two types of initial profiles-hyperbolic tangent (tanh) and boundary layer. The boundary-layer profile leads to more realistic predictions of the transition processes. The predicted transition Reynolds number of 0.18 x 10(exp 6) compares well with experimental data. Normalized vortex spacings for the boundary-layer case are about 3.5 and compare favorably with the 1.5 to 2.5 found in experimental measurements. The tanh profile produces spacings of about 10. The growth rate of the layer is shown to be moderately affected by the initial disturbance field, but comparison with experimental data shows moderate agreement. For the boundary-layer case, it is shown that noise at the Strouhal number of 0.007 is selectively amplified and shows little Reynolds number dependence.

Mukunda, H. S.↗
NASA NTRS1993

Large-Eddy Simulation of Boundary Layer Transition on Swept Wings

The large-eddy simulation of the spatial evolution of a stationary crossflow vortex packet in a three-dimensional boundary layer was performed. Although a coarse grid was used (compared to that required by a direct numerical simulation) the essential features of the disturbance evolution, such as the spanwise disturbance spreading and the vortex rollover, were captured accurately. The eddy viscosity became significant only in the late nonlinear stages of the simulation.

Huai, Xiaoli↗
NASA NTRS1986

Acoustic propagation using computational fluid dynamics

The propagation characteristics of several helicopter airfoil profiles have been investigated using the transonic small disturbance equation. A test case was performed to generate a moving shock that propagated off the airfoil. Various grids were then examined to determine their ability to accurately capture these propagating shock waves. Finally, the case of airfoil-vortex interactions was thoroughly studied over a wide range of Mach numbers and airfoil shapes with particular emphasis on the transonic regime; this results in a highly conplicated fluctuation of lift, drag, and pitching moment. The calculated acoustic intensity levels, along with the details of the computational flow field, provide new insights into the understanding of transonic airfoil-vortex interactions.

Baeder, J. D.↗
NASA NTRS1985

Three-dimensional inviscid flow analysis of turbofan forced mixers

A three-dimensional potential analysis has been formulated and applied to the inviscid flow over a turbofan forced mixer. The method uses a unique small disturbance formulation to analytically uncouple the circumferential flow from the radial and axial flow problem, thereby reducing the analysis to the solution of a series of axisymmetric problems. These equations are discretized using a flux volume formulation along a Cartesian grid. The method extends earlier applications of the Cartesian method to complex cambered geometries. The effects of power addition are also included within the potential formulation. Good agreement is obtained with an alternate small disturbance analysis for a symmetric mixer in a planar duct. In addition calculations showing pressure distributions and induced secondary vorticity fields are presented for practical turbofan mixer configurations, and where possible, comparison has been made with available experimental data.

Barber, T. J.↗
NASA NTRS1988

Turbofan forced mixer lobe flow modeling. 2: Three-dimensional inviscid mixer analysis (FLOMIX)

A three-dimensional potential analysis (FLOMIX) was formulated and applied to the inviscid flow over a turbofan foced mixer. The method uses a small disturbance formulation to analytically uncouple the circumferential flow from the radial and axial flow problem, thereby reducing the analysis to the solution of a series of axisymmetric problems. These equations are discretized using a flux volume formulation along a Cartesian grid. The method extends earlier applications of the Cartesian method to complex cambered geometries. The effects of power addition are also included within the potential formulation. Good agreement is obtained with an alternate small disturbance analysis for a high penetration symmetric mixer in a planar duct. In addition, calculations showing pressure distributions and induced secondary vorticity fields are presented for practical trubofan mixer configurations, and where possible, comparison was made with available experimental data. A detailed description of the required data input and coordinate definition is presented along with a sample data set for a practical forced mixer configuration. A brief description of the program structure and subroutines is also provided.

Barber, T.↗
NASA NTRS1978

Analysis of Satellite-to-Satellite Tracking (SST) and altimetry data from GEOS-C

Radar altimetry and satellite-to-satellite (SST) range and range rate tracking measurements were used to infer the exterior gravitational field of the earth and the structure of the geoid from GEOS-C metric data. Under the SST analysis, a direct point-by-point estimate of gravity disturbance by means of a recursive filter with backward smoothing was attempted but had to be forsaken because of poor convergence. The adopted representation consists of a more or less uniform grid of discrete masses at a depth of approximately 400 km from the earth's surface. The layer is superimposed on a spherical harmonics model. The procedure for smoothing the altimetry and inferring the fine-structured gravity field over the Atlantic test area is described. The local disturbances are represented by means of a density layer. The altimeter height biases were first estimated by a least squares adjustment at orbital crossover points. After taking out the bias, long wavelength contributions from GEM-6 as well as a calibration correction were subtracted. The residual heights were then represented by a mass distribution beneath the earth's surface.

Wong, L.↗
NASA NTRS1981

A potential-flow/boundary-layer method for calculating subsonic and transonic airfoil flow with trailing-edge separation

The development of a potential-flow/boundary-layer method for calculating subsonic and transonic turbulent flow past airfoils with trailing-edge separation is reported. A moment-of-momentum integral boundary-layer method is used which employs the law-of-the-wall/law-of-the-wake velocity profile and a two-layer eddy-viscosity model and ignores the laminar sublayer. All integrals across the boundary layer are obtained in closed form. Separation is assumed to occur when the shearing-stress velocity vanishes. A closed-form solution is derived for separated-flow regions where the shearing stress is negligible. In the potential-flow method, the exact form of the airfoil boundary condition is used, but it is applied at the chord line rather than the airfoil surface. This allows the accurate computation of flow about airfoils at large angles of attack but permits the use of body-oriented Cartesian computational grids. The governing equation for the perturbation velocity potential contains several terms in addition to the classical small-disturbance terms.

Barnwell, R. W.↗