Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear power flow”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling active control of resistive wall mode with power saturation and sensor noise on HL-2M

The resistive wall mode (RWM) control on the HL-2M tokamak is simulated with the MARS-F code (Liu et al 2000 Phys. Plasmas 7 3681), aiming at quantifying control current and voltage requirements when more realistic issues are taken into account, i.e. the control power saturation and the sensor signal noise. The fluid model predicts a narrow stability region for the n = 1 RWM without magnetic feedback, in the 2D parameter space of the plasma pressure versus the toroidal flow speed. Magnetic feedback can fully stabilize the RWM on HL-2M. Without considering the voltage limitation and the sensor signal noise, it is found that plasma flow helps active control of the mode, by reducing the required critical feedback gain for both flux-to-current and flux-to-voltage control schemes. In the absence of the sensor signal noise, the lowest control voltage saturation level, below which the RWM control is lost, is found to roughly satisfy a linear relation to the plasma flow frequency, indicating that subsonic plasma flow is effective in relaxing the control power requirement for the RWM feedback stabilization. The presence of the sensor signal noise substantially modifies the feedback results. A statistical study finds that the sensor signal noise, with the standard deviation of 0.1 G on HL-2M, roughly doubles the required control voltage for successful mode control. The synergistic stabilization effect due to plasma flow is somewhat weakened by the presence of the sensor signal noise. At a given rotation, the tolerable voltage limit generally increases with increasing feedback gain due to the sensor signal noise.

Physics↗
DOE OSTI2021

AGGREGATE: dAta-driven modelinG preservinG contRollable dEr for outaGe mAnagemenT and rEsiliency (Final Report)

The AGGREGATE project team successfully developed and validated various modules for outage management. Brief summaries of each module are provided to showcase their strength for outage management and restoration for a distribution system with a high penetration of connected distribution energy resources (DERs). In recent years, inverter-based DERs have been widely deployed in distribution system. A most of behind-the-meter (BTM) solar power generation is not visible to the utility. The data-driven DER and load estimation modules are using machine learning (ML) and artificial intelligence (AI) to manage this issue, which provides an opportunity for distribution system operators (DSOs) to operate systems and make decisions in real-time for a distribution system with a high penetration of DERs deployed. Also, the estimated DER and true load can be further leveraged in network aggregation and cold-load pick up estimation for reducing the computing complexity and providing for fast restoration. After load demand and DER power generations have been estimated, the information will support topology and state estimation (SE). The topology estimation module demonstrated the viability of mixed integer linear programming (MILP) formulation to estimate the most likely operational radial topology and outage sections using power flow measurements, historical/estimated load and DERs data and smart meter ping measurements. Formulation includes continuous (power flow, load and DERs data) and binary measurements (smart meter ping measurements) in a single formulation. Errors in continuous data and binary data are modeled as normal distribution and Bernoulli distribution, respectively. In the future distribution grid, the power injection from controllable DERs will be essential for efficient and resilient grid operation. However, determining the optimal DER injections and restoration actions is dependent on knowledge of the system states. State estimation (SE), already the cornerstone of transmission energy management systems, will become commonplace in distribution management systems as more measurements become available from deployment of automated metering infrastructure (AMI). Observability analysis is the first step in SE, as it determines the sufficiency of the available measurements for accurately estimating the current system states. A new type of pseudo-measurement called a Correlational Measurement (CM) is introduced in this module, to enhance the observability of the system to enable more accurate SE. CMs encapsulate knowledge of correlation between demand patterns for similar classes of loads as well as injection patterns for same-technology renewable DERs. During grid contingency scenarios, DERs have been traditionally disconnected, without any fault ride-through capabilities. However, with new regulations and better technology, it is feasible for these resources to contribute to the grid’s restoration after an adverse event and hence enhance resilience. The controllability module proposes a two-step restoration scheme for the power system restoration process by leveraging additional degrees of freedom in power electronics interfaced DERs for mitigating voltage problems. In a resilience mode without the utility system, the distribution grid relies on DERs to serve critical load. In such a severe event with multiple faults on the distribution feeders, actuation of various protective devices (PDs) divides the distribution system into electrical islands. The undetected actuated PDs due to fault current contributions from DERs can delay the restoration process, thereby reducing the system resilience. The Advanced Outage Management (AOM) and the Advanced Feeder Restoration (AFR) modules developed in this project provide improved system resilience with multiple DERs. AOM identifies the faulted sections and actuated PDs in a distribution system with DERs by incorporating smart meter data. The most credible outage scenario including fault locations, PD actuations, and fault indicator (FI) failures is identified by a set of binary integer linear programming incorporating hypotheses. The AFR module serves to restore a distribution system with available energy resources taking into consideration the availability of utility sources and DERs. By partitioning the system into islands, critical load will be served with the available generation resources within islands based on the solution of a MILP. When the utility systems become available, the optimal path will be determined by a spanning tree search algorithm that reconnects these islands back to substations and restores the remaining load. The transmission and distribution (T&D) co-simulation module was used to validate the effect of a control action performed on the distribution side assets as it propagates to the transmission side. This ensures that the control action performed results in a feasible operating point on both the transmission and the distribution system. In addition to validation, the team used the T&D co-simulation module to demonstrate how distribution system assets can be used to mitigate issues on the transmission system. Specifically, the team demonstrated that appropriate switching operations on the distribution side can alleviate the line overload condition on the transmission side without causing new operational constraint violations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

HyKKT

HyKKT (pronounced as "hiked") is a package for solving systems of linear equations of Karush-Kuhn-Tucker (KKT) form, which typically arise in optimization problems, such as optimal power flow analysis. HyKKT uses Cholesky instead of LDL^T factorization and solves the general KKT system to a desired numerical precision via block reduction and conjugate gradient on the Schur complement. Such implementation is more suitable for implementation on graphic processing units (GPUs).

Regev, Shaked↗
DOE OSTI2020

Novel Technique for Developing Linearized Convex System Models from Experimentally Derived Data

This paper presents a novel technique for generating a convex system model from an experimentally derived data set which features variance among repeated measurements. The convex system model developed as a test case characterizes the dynamic system losses of a vanadium redox flow battery as a function of the active power output and the battery state of charge. The technique hinges on a pre-cleaning via clustering procedure which precedes the formation of a planar convex hull comprised of triangular simplices. The clustering procedure efficiently reduces the experimental data set while mitigating variance among repeated measurements and removes outliers. Ultimately, the lower evelope of the planar convex hull serves as the desired convex system model. The proposed technique reduces systematic model error which is otherwise present when directly developing a planar convex hull model based on an unreduced data set.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Toward a scalable robust security-constrained optimal power flow using a proximal projection bundle method

Robust security-constrained optimal power flow (rSCOPF) aims to find the worst-case contingencies of alternating current optimal power flow (ACOPF) in power systems. With the rise of GPU architectures on the upcoming supercomputer architectures, optimization algorithms that rely on sparse linear algebra and indefinite linear systems are becoming increasingly hard to solve efficiently (e.g. interior-point method). To address this we revisit a maximin optimization formulation of the rSCOPF and the single-level mixed-integer semidefinite programming (MISDP) reformulation, which is obtained by taking the Lagrangian relaxation of the inner minimization ACOPF problem. In this paper, we focus on the development of a proximal projection bundle method (PPBM) for solving continuous relaxation node subproblems of the MISDP problem, based primarily on the well-known alternating direction method of multipliers. Cutting planes reminiscent of bundle method ideas are also applied in coordination with updates of the proximal parameter. The cutting-plane method can generate a large number of linear inequalities, leading to a large scale but decomposable quadratic programming (QP) subproblem that is amenable to GPUs. We present the numerical results on the IEEE 30, 57, 118, and 300-bus systems by using our PBMM method. We discuss the main computational bottleneck of our method, which is the time taken to solve each iteration of a QP subproblem instance of the PPBM, and how GPU architectures can accelerate this solution process.

bundle method↗
DOE OSTI2020

Hydrokinetic Power Conversion Using Vortex-Induced Oscillation with Cubic Restoring Force

A cubic-spring restoring function with high-deformation stiffening is introduced to passively improve the harnessed marine hydrokinetic power by using flow-induced oscillations/vibrations (FIO/V) of a cylinder. In these FIO/V experiments, a smooth, rigid, single-cylinder on elastic end-supports is tested at Reynolds numbers ranging from 24,000 < Re < 120,000. The parameters of the tested current energy converter (CEC) are cubic stiffness and linear damping. Using the second generation of digital virtual spring-damping (Vck) controller developed by the Marine Renewable Energy Laboratory (MRELab), the cubic modeling of the oscillator stiffness is tested. Experimental results show the influence of the parameter variation on the amplitude, frequency, energy conversion, energy efficiency, and power of the converter. All experiments are conducted in the low turbulence-free surface water (LTFSW) channel of the MRELab of the University of Michigan. The main conclusions are: (1) The nonlinearity in the cubic oscillator is an effective way to extend the vortex-induced vibration (VIV) upper branch, which results in higher harnessing power and efficiency compared to the linear stiffness cylinder converter. (2) Compared to the linear converter, the overall power increase is substantial. The nonlinear power optimum, occurring at the end of the VIV upper branch, is 63% higher than its linear counterpart. (3) The cubic stiffness converter with low harnessing damping achieves consistently good performance in all the VIV regions because of the hardening restoring force, especially at higher flow velocity.

13 HYDRO ENERGY↗
DOE OSTI2021

Dissipation in the magnetic turbulence of reversed field pinch plasmas

Reversed field pinch plasmas are subject to tearing instability that creates a broad spectrum of magnetic fluctuations. The dominant fluctuations have poloidal and toroidal mode numbers (m,n) = (1,6 - 10) and can grow to 2 - 3% of the mean magnetic field. Through non-linear coupling this growth culminates in a strong reconnection event and broadening of the magnetic spectrum extending to the ion gyroradius scale. Multiple developments occur during the reconnection stage: ions and electrons are energized, magnetic fluctuation amplitudes increase, plasma flow is halted, and the toroidal magnetic flux increases in a sawtooth-like fashion as the RFP dynamo becomes stronger. Magnetic fluctuations are measured in the plasma edge at multiple radial locations from r/a = 0.75 to 0.96 to assess and characterize the magnetic turbulence. The measured spectrum perpendicular to the mean field, S(k ⊥ ), can be fit to a model spectrum consisting of power law and exponential component with one free parameter that characterizes dissipation. The measured dissipation is much larger than estimated from classical viscous or resistive dissipation, but it is consistent with a flow damping measurement of anomalous viscosity. The measurements show an evolution of the spectrum during which fluctuation power builds up in the smallest wavenumbers and cascades to the larger wavenumber due to the non-linear coupling between the linear (m,n) = (1, > 6) and the non-linear (m,n) = (0, 1) tearing mo

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Power capture and power take-off load of a self-balanced dual-flap oscillating surge wave energy converter

Wave energy converters are an important part of future renewable energy infrastructure. Predicting their power matrix, capture width ratio, and power take-off loads at a targeted site is required for performance assessment before deployment. Because their testing is very expensive, numerical modeling and simulations play a significant role in those assessments. Linear potential flow theory has limited accuracy under large amplitude wave forcing. More accurate predictions can be obtained by using higher-fidelity models, which are computationally expensive. We present a framework for multi-fidelity numerical simulations to determine the hydrodynamic response, wave capture capability, and power take-off load of a full-scale dual-flap oscillating surge wave energy converter. This design exploits out-of-phase motion by setting the distance between the flaps to half the wavelength of the most occurring wave. The simulations are validated using a 1:10 model experiments in a wave tank. Based on these validations, it was determined that Euler simulations provide an acceptable prediction with 90% reduction in computational time with only 11% error. Utilizing Euler simulations at full-scale, the results demonstrate that the annual electrical energy output is 1.79 GWh under regular wave conditions. Here, one significant improvement over single-flap designs is the capture width ratio which exceeds unity.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Egg yolk as a model for gelation: From rheometry to flow physics

Egg yolks are an excellent model for studying sol-gel transitions, particularly the power law viscoelasticity that defines the critical point of gelation. However, prior studies lack comprehensive datasets and fail to visualize flow behavior linked to temperature and time-dependent linear and nonlinear rheology. Here, we present a detailed dataset characterizing egg yolk viscoelasticity across temperature, time, and forcing amplitude using oscillatory shear, step strain, step stress, and constant high strain rate. Novel protorheology visualizations link rheological properties with observable flow behavior. Our findings highlight the nuanced determination of the critical gel point, emphasizing observation timescale dependencies. We compare methods to identify critical temperatures for gelation, including power law viscoelasticity, moduli crossover, diverging zero-shear viscosity, and emerging equilibrium elastic modulus, while visualizing flow consequences near these transitions. Egg yolk is an accessible, realistic, and nontoxic material relevant to the physicist and the chef alike, making it ideal for understanding the rheology of critical gels. By integrating protorheology photos and videos with rigorous rheometric data, we deepen the understanding of critical gels, with broader impacts for studying other materials with sol-gel transitions.

Marsh, Maxwell C. [Department of Mechanical Scienc↗
DOE OSTI2023

Towards Efficient Alternating Current Optimal Power Flow Analysis on Graphical Processing Units

We present a solution of sparse ACOPF analysis on GPU. In particular, we discuss the performance bottlenecks and detail our efforts to accelerate the linear solver, a core component of ACOPF that dominates the computational time. ACOPF solutions of two large-scale systems, synthetic Northeast (25,000 buses) and Eastern (70,000 buses) \cite{birchfield2017tamu-cases} on GPU show promising speed-up compared to CPU based solution using a state-of-the-art solver. To our knowledge, this is the first result demonstrating acceleration of sparse ACOPF on GPUs.

Power grid analysis, GPU↗
DOE OSTI2021

2D measurements of parallel counter-streaming flows in the W7-X scrape-off layer for attached and detached plasmas

Abstract Investigations of particle parallel flow velocities have been carried out for the scrape-off layer (SOL) of the Wendelstein 7-X (W7-X) stellarator, in order to gain insights on the SOL transport properties during attached and detached plasma scenarios. The experimental evidence is based on the coherence imaging spectroscopy (CIS) diagnostic, able to measure 2D impurity emission intensity and flow velocity. The impurity monitored by CIS is C 2+ , characterized by a line-emission intensity observed to be linearly proportional to the total plasma radiated power in both attached and detached plasmas. The related C 2+ velocity shows a strong dependence on the line-averaged electron density while remaining insensitive to the input power. During attached plasmas, the velocity increases with increasing line-averaged density. The tendency reverses in the transition to and during detachment, in which the velocity decreases by at least a factor of 2. The sharp drop in velocity, together with a rise in line-emission intensity, is reliably correlated to the detachment transition and can therefore be used as one of its signatures. The impurity flow velocity appears to be well coupled with the main ions’ one, thus implying the dominant role of impurity-main ion friction in the parallelimpurity transport dynamics. In view of this SOL impurity transport regime, the CIS measurement results are here interpreted with the help of EMC3-Eirene simulations, and their major trends are already explainable with a simple 1D fluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Vehicle-to-Grid planning framework incorporating electric vehicle user equilibrium and distribution network flexibility enhancement

The rapid surge in electric vehicle (EV) adoption, coupled with advancements in charging technologies, emphasizes the critical necessity for expanding EV recharging infrastructure. Simultaneously, the Distribution Network (DN) encounters escalating challenges in meeting charging demand during peak traffic periods. Consequently, there is a mounting demand for the deployment of innovative Vehicle-to-Grid (V2G) technologies to augment the DN’s flexibility in power dispatch and alleviate travel costs for EV users. Hence, this paper proposes an EV-user-equilibrium-(UE)-constrained V2G planning framework that enhances flexibility in the DN. The framework aims to ascertain the optimal placement and capacity of EV charging stations (EVCSs) and V2G charging piles within the Transportation Network (TN). It takes into account the equilibrium condition stemming from competitive EV charging and routing behaviors alongside the optimal expansion of DN energy resources to accommodate the electricity supplied by the V2G piles. This study commences by analyzing EV drivers’ travel decisions, considering the influence of charging and V2G pile locations and sizes. Subsequently, we tackle the Traffic Assignment Problem with User Equilibrium (TAP-UE) model to characterize the steady-state traffic flow distribution of EVs. Following this, we formulate the optimization model for the Coordinated Power and Transportation Network (CPTN), which encompasses the optimal expansion of DN facilities and traffic flow regulation under UE conditions. To mitigate the computational complexity associated with the V2G planning model, we introduce a series of linearization methods to obtain a manageable Mixed-Integer Linear Programming (MILP) solution. Finally, to validate the efficacy of our proposed planning framework, we apply it to two test systems, including a real-world case study. Through these case studies, we explore the necessity and potential benefits of V2G technologies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Magnetohydrodynamics in free surface liquid metal flow relevant to plasma-facing components

Abstract While flowing Liquid Metal (LM) Plasma-Facing Components (PFCs) represent a potentially transformative technology to enable long-pulse operation with high-power exhaust for fusion reactors, Magnetohydrodynamic (MHD) drag in the conducting LM will reduce the flow speed. Experiments have been completed in the linear open-channel LMX-U device [Hvasta et al 2018 Nucl. Fusion 58 01602] for validation of MHD drag calculations with either insulating or conducting walls, with codes similar to those used to design flowing LM PFCs for a Fusion Nuclear Science Facility [Kessel et al 2019 Fusion Sci. Technol . 75 886]. We observe that the average channel flow speed decreased with the use of conducting walls and the strength of the applied transverse magnetic field. The MHD drag from the retarding Lorentz force resulted in an increase of the LM depth in the channel that ‘piled up’ near the inlet, but not the outlet. As reproduced by OpenFOAM and ANSYS CFX calculations, the magnitude and characteristics of the pileup in the flow direction increased with the applied traverse magnetic field by up to 120%, as compared to the case without an applied magnetic field, corresponding to an average velocity reduction of ∼45%. Particle tracking measurements confirmed a predicted shear in the flow speed, with the surface velocity increasing by 300%, despite the 45% drop in the average bulk speed. The MHD effect makes the bulk flow laminarized but keeps surface waves aligned along the magnetic field lines due to the anisotropy of MHD drag. The 3D fringe field and high surface velocity generate ripples around the outlet region. It was also confirmed that the MHD drag strongly depends on the conductivity of the channel walls, magnetic field, and volumetric flow rate, in agreement with the simulations and a developed analytical model. These validated models are now available to begin to determine the conditions under which the ideal LM channel design of a constant flow speed and fluid depth could be attained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fuel Bonding and its Impact on Axial Gas Communication Behavior in Light-Water Reactor Fuel Rods

Axial gas communication concerns the flow along the axial axis of nuclear fuel rods during ramp and loss of coolant accident (LOCA) conditions. During power ramps, the higher linear heat generation rate may cause fuel-to-clad gap closure that may prevent transport of released fission gases to the plenum. Upon reduction in power the gas then can communicate to the plenum. This phenomenon has been experimentally observed by short power dips during ramp experiments completed at the Risø reactor. At higher burnups it is observed that the UO2 fuel and Zircaloy cladding forms a chemical bond. This bond results in complete closure of the gap. When these high burnup rods are subjected to a LOCA, the bond has implications on both the mechanical response (i.e., ballooning) of the cladding and subsequent fuel relocation and axial gas communication. In the LOCA scenario, gas communication is of interest in two different regimes: 1) pre-rupture communication from the plenum towards the lower pressure ballooning area and 2) the post-rupture depressurization of the plenum to the external system pressure. In both regimes the presence of a fuel-to-cladding bond will impact the rate of depressurization. In this work we present a fuel-to-clad bonding model that is coupled to an existing axial gas communication model framework in the BISON fuel performance code. The effect of considering the bond on fuel performance modeling predictions is presented through comparisons to existing experimental data. Experiments considered include several rods from the Halden IFA-650 test series. An evaluation on a full-length rod that explores the combined effect of plenum size and bonding status on axial gas communication behavior is also presented.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

SAM™ (System Advisor Model™) [SWR-16-02, SWR-10-13]

See the SAM™ website to build a desktop version of the National Laboratory of the Rockies' (NLR's) System Advisor Model™ (SAM). https://sam.nlr.gov/ The System Advisor Model™ (SAM™) is a free techno-economic software model that facilitates decision-making for people in the renewable energy industry: -Project managers and engineers -Policy analysts -Technology developers -Researchers SAM can model many types of renewable energy systems: -Photovoltaic systems, from small residential rooftop to large utility-scale systems -Battery storage with Lithium ion, lead acid, or flow batteries for front-of-meter or behind-the-meter applications -Concentrating Solar Power systems for electric power generation, including parabolic trough, power tower, and linear Fresnel -Industrial process heat from parabolic trough and linear Fresnel systems -Wind power, from individual turbines to large wind farms -Marine energy wave and tidal systems -Solar water heating -Fuel cells -Geothermal power generation -Biomass combustion for power generation -High concentration photovoltaic systems SAM's financial models are for the following types of projects: -Residential and commercial projects where the renewable energy system is on the customer side of the electric utility meter (behind the meter), and power from the system is used to reduce the customer's electricity bill. -Power purchase agreement (PPA) projects where the system is connected to the grid at an interconnection point, and the project earns revenue through power sales. The project may be owned and operated by a single owner or by a partnership involving a flip or leaseback arrangement. -Third party ownership where the system is installed on the customer's (host) property and owned by a separate entity (developer), and the host is compensated for power generated by the system through either a PPA or lease agreement. For a more detailed description of SAM, see Blair et al. (2018), System Advisor Model (SAM) General Description (Version 2017.9.5), NREL/TP-6A20-70414. https://www.nrel.gov/docs/fy18osti/70414.pdf

Ryberg, David↗
DOE OSTI2020

Optimal Network Topology for Node-Breaker Representations With AC Power Flow Constraints

It has been demonstrated that network topology optimization (NTO) may change the topology of power system networks, and consequently, provide additional flexibility to reduce network congestion and violations. Most NTO problems are formulated based on the bus-branch model in which it is challenging to represent a realistic picture of all substation configurations. In this paper, we explore advantages of substation reconfiguration modeling based on node-breaker representations for NTO problem with full nonlinear alternating current power flow. It also proposes a tailored solution algorithm to solve this nonconvex mixed-integer nonlinear programming through the outer approximation method. The proposed solution approach iterates between a mixed-integer linear programming and a nonlinear subproblem. Additional enhancements to further accelerate the iteration process are illustrated. Numerical case studies demonstrate the relative economic and operational impact of optimal network topology with node-breaker representations.

42 ENGINEERING↗
DOE OSTI2021

Additively manufactured electrodes for plasma and power-flow studies in high-power transmission lines on the 1-MA MAIZE facility

Power-flow studies on the 30-MA, 100-ns Z facility at Sandia National Laboratories have shown that plasmas in the facility’s magnetically insulated transmission lines (MITLs) and double post-hole convolute can result in a loss of current delivered to the load. To study power-flow physics on the 1-MA, 100-ns MAIZE facility at the University of Michigan, planar MITL loads and planar post-hole convolute loads have been developed that extend into the lines of sight for various imaging diagnostics on MAIZE. These loads use 3D-printed dielectric support structures lined with thin foils of either aluminum or stainless steel. Here, the metal foils serve as the current-carrying power-flow surfaces, which generate plasma during the current pulse. The foil thickness (50 µm) and widths (11.5–16 mm) are selected to ensure a sufficient linear current density (0.5–0.7 MA/cm) for plasma formation. Laser backlighting (532 nm) and visible-light self-emission imaging capture the overall plasma evolution in the anode–cathode gaps, including the gap closure velocities (1–4 cm/μs).

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Optimal PMU Restoration for Power System Observability Recovery After Massive Attacks

Cyber-physically resilient power system operation requires rapid recovery of situational awareness after disastrous scenarios such as massive cyber attacks. In this article, the concept of optimal sequential restoration of Phasor Measurement Units is introduced, aiming at rapid recovery of power system observability in post-attack scenarios. Additionally, the observability recovery problem is formulated as a Mixed Integer Linear Programming problem to maximize the cumulative observability level gained over time during the restoration process. Three alternative objective functions-maximizing the number of observable buses, the number of observable branches, and the total amount of observable power flow-are considered. The effectiveness of the proposed optimal strategy is verified by comparing it with heuristic approaches on the IEEE 57-bus system.

24 POWER TRANSMISSION AND DISTRIBUTION↗