Showing 1–50 of 60 results for author: Sundar, K

$A^*$ for Graphs of Convex Sets

Authors: Kaarthik Sundar, Sivakumar Rathinam

Abstract: We present a novel algorithm that fuses the existing convex-programming based approach with heuristic information to find optimality guarantees and near-optimal paths for the Shortest Path Problem in the Graph of Convex Sets (SPP-GCS). Our method, inspired by $A^*$, initiates a best-first-like procedure from a designated subset of vertices and iteratively expands it until further growth is neither… ▽ More We present a novel algorithm that fuses the existing convex-programming based approach with heuristic information to find optimality guarantees and near-optimal paths for the Shortest Path Problem in the Graph of Convex Sets (SPP-GCS). Our method, inspired by $A^*$, initiates a best-first-like procedure from a designated subset of vertices and iteratively expands it until further growth is neither possible nor beneficial. Traditionally, obtaining solutions with bounds for an optimization problem involves solving a relaxation, modifying the relaxed solution to a feasible one, and then comparing the two solutions to establish bounds. However, for SPP-GCS, we demonstrate that reversing this process can be more advantageous, especially with Euclidean travel costs. In other words, we initially employ $A^*$ to find a feasible solution for SPP-GCS, then solve a convex relaxation restricted to the vertices explored by $A^*$ to obtain a relaxed solution, and finally, compare the solutions to derive bounds. We present numerical results to highlight the advantages of our algorithm over the existing approach in terms of the sizes of the convex programs solved and computation time. △ Less

Submitted 24 July, 2024; v1 submitted 24 July, 2024; originally announced July 2024.

Optimization of District Heating Network Parameters in Steady-State Operation

Authors: Sai Krishna K. Hari, Anatoly Zlotnik, Shriram Srinivasan, Kaarthik Sundar, Mary Ewers

Abstract: We examine the modeling, simulation, and optimization of district heating systems, which are widely used for thermal transport using steam or hot water as a carrier. We propose a generalizable framework to specify network models and scenario parameters, and develop an optimization method for evaluating system states including pressures, fluid flow rates, and temperatures throughout the network. Th… ▽ More We examine the modeling, simulation, and optimization of district heating systems, which are widely used for thermal transport using steam or hot water as a carrier. We propose a generalizable framework to specify network models and scenario parameters, and develop an optimization method for evaluating system states including pressures, fluid flow rates, and temperatures throughout the network. The network modeling includes pipes, thermal plants, pumps, and passive or controllable loads as system components. We propose basic models for thermodynamic fluid transport and enforce the balance of physical quantities in steady-state flow over co-located outgoing and return networks. We formulate an optimization problem with steam and hot water as the outgoing and return carriers, as in legacy 20th century systems. The physical laws and engineering limitations are specified for each component type, and the thermal network flow optimization (TNFO) problem is formulated and solved for a realistic test network under several scenarios. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Report number: LA-UR-24-23145 MSC Class: 76N25; 90C26; 34B45 ACM Class: J.2

Equitable Routing -- Rethinking the Multiple Traveling Salesman Problem

Authors: Abhay Singh Bhadoriya, Deepjyoti Deka, Kaarthik Sundar

Abstract: The Multiple Traveling Salesman Problem (MTSP) with a single depot is a generalization of the well-known Traveling Salesman Problem (TSP) that involves an additional parameter, namely, the number of salesmen. In the MTSP, several salesmen at the depot need to visit a set of interconnected targets, such that each target is visited precisely once by at most one salesman while minimizing the total le… ▽ More The Multiple Traveling Salesman Problem (MTSP) with a single depot is a generalization of the well-known Traveling Salesman Problem (TSP) that involves an additional parameter, namely, the number of salesmen. In the MTSP, several salesmen at the depot need to visit a set of interconnected targets, such that each target is visited precisely once by at most one salesman while minimizing the total length of their tours. An equally important variant of the MTSP, the min-max MTSP, aims to distribute the workload (length of the individual tours) among salesmen by requiring the longest tour of all the salesmen to be as short as possible, i.e., minimizing the maximum tour length among all salesmen. The min-max MTSP appears in real-life applications to ensure a good balance of workloads for the salesmen. It is known in the literature that the min-max MTSP is notoriously difficult to solve to optimality due to the poor lower bounds its linear relaxations provide. In this paper, we formulate two novel parametric variants of the MTSP called the "fair-MTSP". One variant is formulated as a Mixed-Integer Second Order Cone Program (MISOCP), and the other as a Mixed Integer Linear Program (MILP). Both focus on enforcing the workloads for the salesmen to be equitable, i.e., the distribution of tour lengths for the salesmen to be fair while minimizing the total cost of their tours. We present algorithms to solve the two variants of the fair-MTSP to global optimality and computational results on benchmark and real-world test instances that make a case for fair-MTSP as a viable alternative to the min-max MTSP. △ Less

Submitted 16 April, 2024; v1 submitted 11 April, 2024; originally announced April 2024.

Comments: 10 pages

Report number: LA-UR-24-23196 MSC Class: 90B06; 90C27; 90C57; 90C90

Deep Reinforcement Learning-Based Approach for a Single Vehicle Persistent Surveillance Problem with Fuel Constraints

Authors: Manav Mishra, Hritik Bana, Saswata Sarkar, Sujeevraja Sanjeevi, PB Sujit, Kaarthik Sundar

Abstract: This article presents a deep reinforcement learning-based approach to tackle a persistent surveillance mission requiring a single unmanned aerial vehicle initially stationed at a depot with fuel or time-of-flight constraints to repeatedly visit a set of targets with equal priority. Owing to the vehicle's fuel or time-of-flight constraints, the vehicle must be regularly refueled, or its battery mus… ▽ More This article presents a deep reinforcement learning-based approach to tackle a persistent surveillance mission requiring a single unmanned aerial vehicle initially stationed at a depot with fuel or time-of-flight constraints to repeatedly visit a set of targets with equal priority. Owing to the vehicle's fuel or time-of-flight constraints, the vehicle must be regularly refueled, or its battery must be recharged at the depot. The objective of the problem is to determine an optimal sequence of visits to the targets that minimizes the maximum time elapsed between successive visits to any target while ensuring that the vehicle never runs out of fuel or charge. We present a deep reinforcement learning algorithm to solve this problem and present the results of numerical experiments that corroborate the effectiveness of this approach in comparison with common-sense greedy heuristics. △ Less

Submitted 2 May, 2024; v1 submitted 9 April, 2024; originally announced April 2024.

Comments: 6 pages

Report number: LA-UR-24-23186

Stochastic Finite Volume Method for Uncertainty Management in Gas Pipeline Network Flows

Authors: Saif R. Kazi, Sidhant Misra, Svetlana Tokareva, Kaarthik Sundar, Anatoly Zlotnik

Abstract: Natural gas consumption by users of pipeline networks is subject to increasing uncertainty that originates from the intermittent nature of electric power loads serviced by gas-fired generators. To enable computationally efficient optimization of gas network flows subject to uncertainty, we develop a finite volume representation of stochastic solutions of hyperbolic partial differential equation (P… ▽ More Natural gas consumption by users of pipeline networks is subject to increasing uncertainty that originates from the intermittent nature of electric power loads serviced by gas-fired generators. To enable computationally efficient optimization of gas network flows subject to uncertainty, we develop a finite volume representation of stochastic solutions of hyperbolic partial differential equation (PDE) systems on graph-connected domains with nodal coupling and boundary conditions. The representation is used to express the physical constraints in stochastic optimization problems for gas flow allocation subject to uncertain parameters. The method is based on the stochastic finite volume approach that was recently developed for uncertainty quantification in transient flows represented by hyperbolic PDEs on graphs. In this study, we develop optimization formulations for steady-state gas flow over actuated transport networks subject to probabilistic constraints. In addition to the distributions for the physical solutions, we examine the dual variables that are produced by way of the optimization, and interpret them as price distributions that quantify the financial volatility that arises through demand uncertainty modeled in an optimization-driven gas market mechanism. We demonstrate the computation and distributional analysis using a single-pipe example and a small test network. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Report number: LA-UR-24-22647 MSC Class: 65Kxx; 90C35; 90C15

Exact and Heuristic Approaches for the Stochastic N-k Interdiction in Power Grids

Authors: Kaarthik Sundar, Andrew Mastin, Manuel Garcia, Russell Bent, Jean-Paul Watson

Abstract: The article introduces the stochastic N-k interdiction problem for power grid operations and planning that aims to identify a subset of k components (out of N components) that maximizes the expected damage, measured in terms of load shed. Uncertainty is modeled through a fixed set of outage scenarios, where each scenario represents a subset of components removed from the grid. We formulate the sto… ▽ More The article introduces the stochastic N-k interdiction problem for power grid operations and planning that aims to identify a subset of k components (out of N components) that maximizes the expected damage, measured in terms of load shed. Uncertainty is modeled through a fixed set of outage scenarios, where each scenario represents a subset of components removed from the grid. We formulate the stochastic N-k interdiction problem as a bi-level optimization problem and propose two algorithmic solutions. The first approach reformulates the bi-level stochastic optimization problem to a single level, mixed-integer linear program (MILP) by dualizing the inner problem and solving the resulting problem directly using a MILP solver to global optimality. The second is a heuristic cutting-plane approach, which is exact under certain assumptions. We compare these approaches in terms of computation time and solution quality using the IEEE-Reliability Test System and present avenues for future research. △ Less

Submitted 31 January, 2024; originally announced February 2024.

Comments: 6 pages, 2 figures

Report number: LA-UR-24-20722

Locational Marginal Pricing of Energy in Pipeline Transport of Natural Gas and Hydrogen with Carbon Offset Incentives

Authors: Mo Sodwatana, Saif R. Kazi, Kaarthik Sundar, Adam Brandt, Anatoly Zlotnik

Abstract: We propose an optimization formulation for locational pricing of energy transported through a pipeline network that carries mixtures of natural gas and hydrogen from distributed sources to consumers. The objective includes the economic value provided by the pipeline to consumers of energy and suppliers of natural gas and green hydrogen, as well as incentives to lower carbon emissions by consuming… ▽ More We propose an optimization formulation for locational pricing of energy transported through a pipeline network that carries mixtures of natural gas and hydrogen from distributed sources to consumers. The objective includes the economic value provided by the pipeline to consumers of energy and suppliers of natural gas and green hydrogen, as well as incentives to lower carbon emissions by consuming the latter instead of the former. The optimization is subject to the physics of gas flow and mixing in the pipeline network as well as engineering limits. In addition to formulating this mathematical program, we synthesize the Lagrangian and derive analytical expressions for the dual variables. We propose that the dual solution can be used to derive locational marginal prices of natural gas, hydrogen, and energy, as well as the decarbonization premium paid by consumers that receive hydrogen. We derive several properties of solutions obtained using the proposed market mechanism, and demonstrate them using case studies for standard 8-node and 40-node pipeline test networks. Finally, we show that optimization-based analysis of the type proposed here is critical for making sound decisions about economic policy and infrastructure expansion for blending green hydrogen into existing natural gas pipelines. △ Less

Submitted 16 May, 2024; v1 submitted 19 October, 2023; originally announced October 2023.

Report number: LA-UR-23-31810 MSC Class: 90C46; 90C35

Fairly Extreme: Minimizing Outages Equitably

Authors: Kaarthik Sundar, Deepjyoti Deka, Russell Bent

Abstract: This paper focuses on the problem of minimizing the outages due to extreme events on the power grid equitably among all customers of the grid. The paper presents two ways of incorporating fairness into the existing formulations that seek to minimize the total outage in the power grid. The first method is motivated by existing literature on incorporating fairness in optimization problems and this i… ▽ More This paper focuses on the problem of minimizing the outages due to extreme events on the power grid equitably among all customers of the grid. The paper presents two ways of incorporating fairness into the existing formulations that seek to minimize the total outage in the power grid. The first method is motivated by existing literature on incorporating fairness in optimization problems and this is done by modifying the problem's objective function. The second method introduces a novel notion of fairness, termed $\varepsilon$-fairness, that can be incorporated into existing problem formulations through a single second-order cone constraint. Both these methods are very general and can be used to incorporate fairness in existing planning and operational optimization problems in the power grid and beyond. Extensive computational case studies that examine the effectiveness of both these methods to characterize fairness is presented followed by conclusions and ways forward. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: 8 pages

Report number: LA-UR-23-31124

On the Existence of Steady-State Solutions to the Equations Governing Fluid Flow in Networks

Authors: Shriram Srinivasan, Nishant Panda, Kaarthik Sundar

Abstract: The steady-state solution of fluid flow in pipeline infrastructure networks driven by junction/node potentials is a crucial ingredient in various decision-support tools for system design and operation. While the nonlinear system is known to have a unique solution (when one exists), the absence of a definite result on the existence of solutions hobbles the development of computational algorithms, f… ▽ More The steady-state solution of fluid flow in pipeline infrastructure networks driven by junction/node potentials is a crucial ingredient in various decision-support tools for system design and operation. While the nonlinear system is known to have a unique solution (when one exists), the absence of a definite result on the existence of solutions hobbles the development of computational algorithms, for it is not possible to distinguish between algorithm failure and non-existence of a solution. In this letter, we show that for any fluid whose equation of state is a scaled monomial, a unique solution exists for such nonlinear systems if the term solution is interpreted in terms of potentials and flows rather than pressures and flows. However, for gases following the CNGA equation of state, while the question of existence remains open, we construct an alternative system that always has a unique solution and show that the solution to this system is a good approximant of the true solution. The existence result for flow of natural gas in networks also applies to other fluid flow networks such as water distribution networks or networks that transport carbon dioxide in carbon capture and sequestration. Most importantly, our result enables correct diagnosis of algorithmic failure, problem stiffness, and non-convergence in computational algorithms. △ Less

Submitted 8 April, 2024; v1 submitted 6 September, 2023; originally announced September 2023.

Comments: 6 pages, 3 figures

Report number: LA-UR-23-30213

Relaxations of the Steady Optimal Gas Flow Problem for a Non-Ideal Gas

Authors: Sai Krishna Kanth Hari, Kaarthik Sundar, Shriram Srinivasan, Russell Bent

Abstract: Natural gas ranks second in consumption among primary energy sources in the United States. The majority of production sites are in remote locations, hence natural gas needs to be transported through a pipeline network equipped with a variety of physical components such as compressors, valves, etc. Thus, from the point of view of both economics and reliability, it is desirable to achieve optimal tr… ▽ More Natural gas ranks second in consumption among primary energy sources in the United States. The majority of production sites are in remote locations, hence natural gas needs to be transported through a pipeline network equipped with a variety of physical components such as compressors, valves, etc. Thus, from the point of view of both economics and reliability, it is desirable to achieve optimal transportation of natural gas using these pipeline networks. The physics that governs the flow of natural gas through various components in a pipeline network is governed by nonlinear and non-convex equality and inequality constraints and the most general steady-flow operations problem takes the form of a Mixed Integer Nonlinear Program. In this paper, we consider one example of steady-flow operations -- the Optimal Gas Flow (OGF) problem for a natural gas pipeline network that minimizes the production cost subject to the physics of steady-flow of natural gas. The ability to quickly determine global optimal solution and a lower bound to the objective value of the OGF for different demand profiles plays a key role in efficient day-to-day operations. One strategy to accomplish this relies on tight relaxations to the nonlinear constraints of the OGF. Currently, many nonlinear constraints that arise due to modeling the non-ideal equation of state either do not have relaxations or have relaxations that scale poorly for realistic network sizes. In this work, we combine recent advancements in the development of polyhedral relaxations for univariate functions to obtain tight relaxations that can be solved within a few seconds on a standard laptop. We demonstrate the quality of these relaxations through extensive numerical experiments on very large scale test networks available in the literature and find that the proposed relaxation is able to prove optimality in 92% of the instances. △ Less

Submitted 24 August, 2023; originally announced August 2023.

Comments: 28 pages

Dynamic Optimization and Optimal Control of Hydrogen Blending Operations in Natural Gas Networks

Authors: Saif R. Kazi, Kaarthik Sundar, Anatoly Zlotnik

Abstract: We present a dynamic model for the optimal control problem (OCP) of hydrogen blending into natural gas pipeline networks subject to inequality constraints. The dynamic model is derived using the first principles partial differential equations (PDEs) for the transport of heterogeneous gas mixtures through long distance pipes. Hydrogen concentration is tracked together with the pressure and mass flo… ▽ More We present a dynamic model for the optimal control problem (OCP) of hydrogen blending into natural gas pipeline networks subject to inequality constraints. The dynamic model is derived using the first principles partial differential equations (PDEs) for the transport of heterogeneous gas mixtures through long distance pipes. Hydrogen concentration is tracked together with the pressure and mass flow dynamics within the pipelines, as well as mixing and compatibility conditions at nodes, actuation by compressors, and injection of hydrogen or natural gas into the system or withdrawal of the mixture from the network. We implement a lumped parameter approximation to reduce the full PDE model to a differential algebraic equation (DAE) system that can be easily discretized and solved using nonlinear optimization or programming (NLP) solvers. We examine a temporal discretization that is advantageous for time-periodic boundary conditions, parameters, and inequality constraint bound values. The method is applied to solve case studies for a single pipe and a multi-pipe network with time-varying parameters in order to explore how mixing of heterogeneous gases affects pipeline transient optimization. △ Less

Submitted 23 February, 2024; v1 submitted 5 April, 2023; originally announced April 2023.

Report number: LA-UR-23-23439

Optimal Power System Topology Control Under Uncertain Wildfire Risk

Authors: Yuqi Zhou, Kaarthik Sundar, Deepjyoti Deka, Hao Zhu

Abstract: Wildfires pose a significant threat to the safe and reliable operation of electric power systems. They can quickly spread and cause severe damage to power infrastructure. To reduce the risk, public safety power shutoffs are often used to restore power balance and prevent widespread blackouts. However, the unpredictability of wildfires makes it challenging to implement effective counter-measures in… ▽ More Wildfires pose a significant threat to the safe and reliable operation of electric power systems. They can quickly spread and cause severe damage to power infrastructure. To reduce the risk, public safety power shutoffs are often used to restore power balance and prevent widespread blackouts. However, the unpredictability of wildfires makes it challenging to implement effective counter-measures in a timely manner. This paper proposes an optimization-based topology control problem as a solution to mitigate the impact of wildfires. The goal is to find the optimal network topology that minimizes total operating costs and prevents load shedding during power shutoffs under uncertain line shutoff scenarios caused by uncertain spreading wildfires. The solution involves solving two-stage stochastic mixed-integer linear programs, with preventive and corrective control actions taken when the risk of a wildfire and corresponding outage line scenarios are known with low and high confidence, respectively. The Progressive Hedging algorithm is used to solve both problems. The effectiveness of the proposed approach is demonstrated using data from the RTS-GMLC system that is artificially geo-located in Southern California, including actual wildfire incidents and risk maps. Our work provides a crucial study of the comparative benefits due to accurate risk forecast and corresponding preventive control over real-time corrective control that may not be realistic. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 10 pages, 6 figures

Modeling and Optimization of Steady Flow of Natural Gas and Hydrogen Mixtures in Pipeline Networks

Authors: Saif R. Kazi, Kaarthik Sundar, Shriram Srinivasan, Anatoly Zlotnik

Abstract: We extend the canonical problems of simulation and optimization of steady-state gas flows in pipeline networks with compressors to the transport of mixtures of highly heterogeneous gases injected throughout a network. Our study is motivated by proposed projects to blend hydrogen generated using clean energy into existing natural gas pipeline systems as part of efforts to reduce the reliance of ene… ▽ More We extend the canonical problems of simulation and optimization of steady-state gas flows in pipeline networks with compressors to the transport of mixtures of highly heterogeneous gases injected throughout a network. Our study is motivated by proposed projects to blend hydrogen generated using clean energy into existing natural gas pipeline systems as part of efforts to reduce the reliance of energy systems on fossil fuels. Flow in a pipe is related to endpoint pressures by a basic Weymouth equation model, with an ideal gas equation of state, where the wave speed depends on the hydrogen concentration. At vertices, in addition to mass balance, we also consider mixing of incoming flows of varying hydrogen concentrations. The problems of interest are the heterogeneous gas flow simulation (HGFS), which determines system pressures and flows given fixed boundary conditions and compressor settings, as well as the heterogeneous gas flow optimization (HGFO), which extremizes an objective by determining optimal boundary conditions and compressor settings. We examine conditions for uniqueness of solutions to the HGFS, as well as compare and contrast mixed-integer and continuous nonlinear programming formulations for the HGFO. We develop computational methods to solve both problems, and examine their performance using four test networks of increasing complexity. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Comments: 12 pages

Report number: LA-UR-2232518

Optimization of Hydrogen Blending in Natural Gas Networks for Carbon Emissions Reduction

Authors: Mo Sodwatana, Saif R. Kazi, Kaarthik Sundar, Anatoly Zlotnik

Abstract: We present an economic optimization problem for allocating the flow of natural gas and hydrogen blends through a large-scale transportation pipeline network. Physical flow of the gas mixture is modeled using a steady-state relation between pressure decrease and flow rate, which depends on mass concentration of the constituents as it varies by location in the network. The objective reflects the eco… ▽ More We present an economic optimization problem for allocating the flow of natural gas and hydrogen blends through a large-scale transportation pipeline network. Physical flow of the gas mixture is modeled using a steady-state relation between pressure decrease and flow rate, which depends on mass concentration of the constituents as it varies by location in the network. The objective reflects the economic value provided by the system, accounting for delivered energy in withdrawn flows, the cost of natural gas and hydrogen injections, and avoided carbon emissions. The problem is solved subject to physical flow equations, nodal balance and mixing laws, and engineering inequality constraints. The desired energy delivery rate and minimum hydrogen concentration can be specified as upper and lower bound values, respectively, of inequality constraints, and we examine the sensitivity of the physical pressure and flow solution to these parameters for two test networks. The results confirm that increasing hydrogen concentration requires greater energy expended for compression to deliver the same energy content, and the formulation could be used for valuation of the resulting mitigation of carbon emissions. △ Less

Submitted 22 May, 2023; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: In Proceedings of the 2023 American Control Conference

Report number: LA-UR-22-30127

Heuristics for Multi-Vehicle Routing Problem Considering Human-Robot Interactions

Authors: Venkata Sirimuvva Chirala, Kaarthik Sundar, Saravanan Venkatachalam, Jonathon M. Smereka, Sam Kassoumeh

Abstract: Unmanned ground vehicles (UGVs) are being used extensively in civilian and military applications for applications such as underground mining, nuclear plant operations, planetary exploration, intelligence, surveillance and reconnaissance (ISR) missions and manned-unmanned teaming. We consider a multi-objective, multiple-vehicle routing problem in which teams of manned ground vehicles (MGVs) and UGV… ▽ More Unmanned ground vehicles (UGVs) are being used extensively in civilian and military applications for applications such as underground mining, nuclear plant operations, planetary exploration, intelligence, surveillance and reconnaissance (ISR) missions and manned-unmanned teaming. We consider a multi-objective, multiple-vehicle routing problem in which teams of manned ground vehicles (MGVs) and UGVs are deployed respectively in a leader-follower framework to execute missions with differing requirements for MGVs and UGVs while considering human-robot interactions (HRI). HRI studies highlight the costs of managing a team of follower UGVs by a leader MGV. This paper aims to compute feasible paths, replenishments, team compositions and number of MGV-UGV teams deployed such that the requirements for MGVs and UGVs for the missions are met and the path, replenishment, HRI and team deployment costs are at minimum. The problem is first modeled as a a mixed-integer linear program (MILP) that can be solved to optimality by off-the-shelf commercial solvers for small-sized instances. For larger instances, a variable neighborhood search algorithm is offered to compute near optimal solutions and address the challenges that arise when solving the combinatorial multi-objective routing optimization problem. Finally, computational experiments that corroborate the effectiveness of the proposed algorithms are presented. △ Less

Submitted 20 August, 2022; originally announced August 2022.

Comments: 10 pages

Numerical Solution of the Steady-State Network Flow Equations for a Non-Ideal Gas

Authors: Shriram Srinivasan, Kaarthik Sundar, Vitaliy Gyrya, Anatoly Zlotnik

Abstract: We formulate a steady-state network flow problem for non-ideal gas that relates injection rates and nodal pressures in the network to flows in pipes. For this problem, we present and prove a theorem on uniqueness of generalized solution for a broad class of non-ideal pressure-density relations that satisfy a monotonicity property. Further, we develop a Newton-Raphson algorithm for numerical soluti… ▽ More We formulate a steady-state network flow problem for non-ideal gas that relates injection rates and nodal pressures in the network to flows in pipes. For this problem, we present and prove a theorem on uniqueness of generalized solution for a broad class of non-ideal pressure-density relations that satisfy a monotonicity property. Further, we develop a Newton-Raphson algorithm for numerical solution of the steady-state problem, which is made possible by a systematic non-dimensionalization of the equations. The developed algorithm has been extensively tested on benchmark instances and shown to converge robustly to a generalized solution. Previous results indicate that the steady-state network flow equations for an ideal gas are difficult to solve by the Newton-Raphson method because of its extreme sensitivity to the initial guess. In contrast, we find that non-dimensionalization of the steady-state problem is key to robust convergence of the Newton-Raphson method. We identify criteria based on the uniqueness of solutions under which the existence of a non-physical generalized solution found by a non-linear solver implies non-existence of a physical solution, i.e., infeasibility of the problem. Finally, we compare pressure and flow solutions based on ideal and non-ideal equations of state to demonstrate the need to apply the latter in practice. The solver developed in this article is open-source and is made available for both the academic and research communities as well as the industry. △ Less

Submitted 29 November, 2022; v1 submitted 31 March, 2022; originally announced April 2022.

Comments: 12 pages

Report number: LA-UR-22-21789

Routing Problem for Unmanned Aerial Vehicle Patrolling Missions -- A Progressive Hedging Algorithm

Authors: Sudarshan Rajan, Kaarthik Sundar, Natarajan Gautam

Abstract: The paper presents a two-stage stochastic program to model a routing problem involving an Unmanned Aerial Vehicle (UAV) in the context of patrolling missions. In particular, given a set of targets and a set of supplemental targets corresponding to each target, the first stage decisions involve finding the sequence in which the vehicle has to visit the set of targets. Upon reaching each target, the… ▽ More The paper presents a two-stage stochastic program to model a routing problem involving an Unmanned Aerial Vehicle (UAV) in the context of patrolling missions. In particular, given a set of targets and a set of supplemental targets corresponding to each target, the first stage decisions involve finding the sequence in which the vehicle has to visit the set of targets. Upon reaching each target, the UAV collects information and if the operator of the UAV deems that the information collected is not of sufficient fidelity, then the UAV has to visit all the supplemental targets corresponding to that target to collect additional information before proceeding to visit the next target. The problem is solved using a progressive hedging algorithm and extensive computational results corroborating the effectiveness of the proposed model and the solution methodology is presented. △ Less

Submitted 29 June, 2021; v1 submitted 15 June, 2021; originally announced June 2021.

Comments: 16 pages

Report number: LA-UR-21-24769

Demand Response Analogues for Residential Loads in Natural Gas Networks

Authors: Fuyu Hu, Kaarthik Sundar, Shriram Srinivasan, Russell Bent

Abstract: Demand response for electrical power networks is a mature field that has yielded numerous efficiency and resilience benefits like managing the peaks and valleys of electricity usage, reduction in peak electricity usage, to name a few. However, only recently has the study of the counterpart to demand response for electric power in natural gas started to receive similar levels of attention. Natural… ▽ More Demand response for electrical power networks is a mature field that has yielded numerous efficiency and resilience benefits like managing the peaks and valleys of electricity usage, reduction in peak electricity usage, to name a few. However, only recently has the study of the counterpart to demand response for electric power in natural gas started to receive similar levels of attention. Natural gas systems are increasingly operating at or near capacity, which challenges these systems to meet all the needs for gas, especially during severe winter weather. However, unlike demand response programs in electrical networks, demand response in gas networks cannot shift peak usage. Here, we develop analogues to demand response that can help improve the resilience of natural gas systems by reducing peak consumption and thereby limiting potential disruptions such events can cause. This paper develops a mathematical formulation to support a residential-level demand response analogue for natural gas based on current and anticipated smart thermostat technologies. The mathematical formulation takes the form of an optimal control problem (OCP) that leverages physics constraints to model temperature changes, balance equitable service, and optimize the gas consumption for collections of houses. On test problems, the formulation demonstrates significant benefits, including the ability to cut peak demand by 15% while still ensuring equitable service to customers. △ Less

Submitted 1 December, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: 27 pages

Report number: LA-UR-21-23328

Operation of Natural Gas Pipeline Networks With Storage Under Transient Flow Conditions

Authors: Sai Krishna Kanth Hari, Kaarthik Sundar, Shriram Srinivasan, Anatoly Zlotnik, Russell Bent

Abstract: We formulate a nonlinear optimal control problem for intra-day operation of a natural gas pipeline network that includes storage reservoirs. The dynamics of compressible gas flow through pipes, compressors, reservoirs, and wells are considered. In particular, a reservoir is modeled as a rigid, hollow container that stores gas under isothermal conditions and uniform density, and a well is modeled a… ▽ More We formulate a nonlinear optimal control problem for intra-day operation of a natural gas pipeline network that includes storage reservoirs. The dynamics of compressible gas flow through pipes, compressors, reservoirs, and wells are considered. In particular, a reservoir is modeled as a rigid, hollow container that stores gas under isothermal conditions and uniform density, and a well is modeled as a vertical pipe. For each pipe, flow dynamics are described by a coupled partial differential equation (PDE) system in density and mass flux variables, with momentum dissipation modeled using the Darcy-Wiesbach friction approximation. Compressors are modeled as scaling up the pressure of gas between inlet and outlet. The governing equations for all network components are spatially discretized and assembled into a nonlinear differential-algebraic equation (DAE) system, which synthesizes above-ground pipeline and subsurface reservoir dynamics into a single reduced-order model. We seek to maximize an objective function that quantifies economic profit and network efficiency subject to the flow equations and inequalities that represent operating limitations. The problem is solved using a primal-dual interior point solver and the solutions are validated in computational experiments and simulations on several pipeline test networks to demonstrate the effectiveness of the proposed methodology. △ Less

Submitted 4 May, 2021; v1 submitted 3 March, 2021; originally announced March 2021.

Comments: 13 pages

Report number: LA-UR-20-27442

Robust Gas Pipeline Network Expansion Planning to Support Power System Reliability

Authors: Kaarthik Sundar, Sidhant Misra, Anatoly Zlotnik, Russell Bent

Abstract: We examine the problem of optimal transport capacity expansion planning for a gas pipeline network to service the growing demand of gas-fired power plants that are increasingly used to provide base load, flexibility, and reserve generation for bulk electric system. The aim is to determine the minimal cost set of additional pipes and gas compressors that can be added to the network to provide the a… ▽ More We examine the problem of optimal transport capacity expansion planning for a gas pipeline network to service the growing demand of gas-fired power plants that are increasingly used to provide base load, flexibility, and reserve generation for bulk electric system. The aim is to determine the minimal cost set of additional pipes and gas compressors that can be added to the network to provide the additional capacity to service future loads. This combinatorial optimization problem is initially formulated as a mixed-integer nonlinear program, which we then extend to account for the variability that is inherent to the demands of gas-fired electricity production and uncertainty in expected future loads. We consider here steady-state flow modeling while ensuring that the solution is feasible for all possible values of interval uncertainty in loads, which results in a challenging semi-infinite problem. We apply previously derived monotonicity properties that enable simplification of the problem to require constraint satisfaction in the two extremal scenarios only, and then formulate the robust gas pipeline network expansion planning problem using a mixed-integer second order cone formulation. We consider case studies on the Belgian network test case to examine the performance of the proposed approach. △ Less

Submitted 25 January, 2021; originally announced January 2021.

Comments: 8 pages

Report number: LA-UR-19-30191

Deployable, Data-Driven Unmanned Vehicle Navigation System in GPS-Denied, Feature-Deficient Environments

Authors: Sohum Misra, Kaarthik Sundar, Rajnikant Sharma, Kevin Brink

Abstract: This paper presents a novel data-driven navigation system to navigate an Unmanned Vehicle (UV) in GPS-denied, feature-deficient environments such as tunnels, or mines. The method utilizes landmarks that vehicle can deploy and measure range from to enable localization as the vehicle traverses its pre-defined path through the tunnel. A key question that arises in such scenario is to estimate and red… ▽ More This paper presents a novel data-driven navigation system to navigate an Unmanned Vehicle (UV) in GPS-denied, feature-deficient environments such as tunnels, or mines. The method utilizes landmarks that vehicle can deploy and measure range from to enable localization as the vehicle traverses its pre-defined path through the tunnel. A key question that arises in such scenario is to estimate and reduce the number of landmarks that needs to be deployed for localization before the start of the mission, given some information about the environment. The main focus is to keep the maximum position uncertainty at a desired value. In this article, we develop a novel vehicle navigation system in GPS-denied, feature-deficient environment by combining techniques from estimation, machine learning, and mixed-integer convex optimization. This article develops a novel, systematic method to perform localization and navigate the UV through the environment with minimum number of landmarks while maintaining desired localization accuracy. We also present extensive simulation experiments on different scenarios that corroborate the effectiveness of the proposed navigation system. △ Less

Submitted 2 November, 2021; v1 submitted 24 January, 2021; originally announced January 2021.

Comments: 37 pages

Report number: LA-UR-21-20249

Optimal Economic Operation of Liquid Petroleum Products Pipeline Systems

Authors: Elena Khlebnikova, Kaarthik Sundar, Anatoly Zlotnik, Russell Bent, Mary Ewers, Byron Tasseff

Abstract: The majority of overland transport needs for crude petroleum and refined petroleum products are met using pipelines. Numerous studies have developed optimization methods for design of these systems in order to minimize construction costs while meeting capacity requirements. Here, we formulate problems to optimize the operations of existing single liquid commodity pipeline systems subject to physic… ▽ More The majority of overland transport needs for crude petroleum and refined petroleum products are met using pipelines. Numerous studies have developed optimization methods for design of these systems in order to minimize construction costs while meeting capacity requirements. Here, we formulate problems to optimize the operations of existing single liquid commodity pipeline systems subject to physical flow and pump engineering constraints. The objectives are to maximize the economic value created for users of the system and to minimize operating costs. We present a general computational method for this class of continuous, non-convex nonlinear programs, and examine the use of pump operating settings and flow allocations as decision variables. The approach is applied to compute optimal operating regimes and perform engineering economic sensitivity analyses for a case study of a crude oil pipeline developed using publicly available data. △ Less

Submitted 21 December, 2020; originally announced December 2020.

MSC Class: 76B75; 90C35

Graph-Based Modeling and Decomposition of Energy Infrastructures

Authors: Sungho Shin, Carleton Coffrin, Kaarthik Sundar, Victor M. Zavala

Abstract: Nonlinear optimization problems are found at the heart of real-time operations of critical infrastructures. These problems are computationally challenging because they embed complex physical models that exhibit space-time dynamics. We propose modeling these problems as graph-structured optimization problems, and illustrate how their structure can be exploited at the modeling level (for parallelizi… ▽ More Nonlinear optimization problems are found at the heart of real-time operations of critical infrastructures. These problems are computationally challenging because they embed complex physical models that exhibit space-time dynamics. We propose modeling these problems as graph-structured optimization problems, and illustrate how their structure can be exploited at the modeling level (for parallelizing function/derivative computations) and at the solver level (for parallelizing linear algebra operations). Specifically, we present a restricted additive Schwarz scheme that enables flexible decomposition of complex graph structures within an interior-point algorithm. The proposed approach is implemented as a general-purpose nonlinear programming solver that we call MadNLP.jl; this Julia-based solver is interfaced to the graph-based modeling package Plasmo.jl. The efficiency of this framework is demonstrated via problems arising in transient gas network optimization and multi-period AC optimal power flow. We show that our framework accelerates the solution (compared to off-the-shelf tools) by over 300%; specifically, solution times are reduced from 72.36 sec to 23.84 sec for the gas problem and from 515.81 sec to 149.45 sec for the power flow problem. △ Less

Submitted 2 March, 2021; v1 submitted 5 October, 2020; originally announced October 2020.

Natural Gas Maximal Load Delivery for Multi-contingency Analysis

Authors: Byron Tasseff, Carleton Coffrin, Russell Bent, Kaarthik Sundar, Anatoly Zlotnik

Abstract: As the use of renewable generation has increased, electric power systems have become increasingly reliant on natural gas-fired power plants as fast ramping sources for meeting fluctuating bulk power demands. This dependence has introduced new vulnerabilities to the power grid, including disruptions to gas transmission networks from natural and man-made disasters. To address the operational challen… ▽ More As the use of renewable generation has increased, electric power systems have become increasingly reliant on natural gas-fired power plants as fast ramping sources for meeting fluctuating bulk power demands. This dependence has introduced new vulnerabilities to the power grid, including disruptions to gas transmission networks from natural and man-made disasters. To address the operational challenges arising from these disruptions, we consider the task of determining a feasible steady-state operating point for a damaged gas pipeline network while ensuring the maximal delivery of load. We formulate the mixed-integer nonconvex maximal load delivery (MLD) problem, which proves difficult to solve on large-scale networks. To address this challenge, we present a mixed-integer convex relaxation of the MLD problem and use it to determine bounds on the transport capacity of a gas pipeline system. To demonstrate the effectiveness of the relaxation, the exact and relaxed formulations are compared across a large number of randomized damage scenarios on nine natural gas pipeline network models ranging in size from 11 to 4197 junctions. A proof of concept application, which assumes network damage from a set of synthetically generated earthquakes, is also presented to demonstrate the utility of the proposed optimization-based capacity evaluation in the context of risk assessment for natural disasters. For all but the largest network, the relaxation-based method is found to be suitable for use in evaluating the impacts of multi-contingency network disruptions, often converging to the optimal solution of the relaxed formulation in less than ten seconds. △ Less

Submitted 1 October, 2020; v1 submitted 30 September, 2020; originally announced September 2020.

Comments: Corrected alignment of histograms in Figure 2

Report number: LA-UR-20-27610

An Uncertainty Management Framework for Integrated Gas-Electric Energy Systems

Authors: Line A. Roald, Kaarthik Sundar, Anatoly Zlotnik, Sidhant Misra, Göran Andersson

Abstract: In many parts of the world, electric power systems have seen a significant shift towards generation from renewable energy and natural gas. Because of their ability to flexibly adjust power generation in real time, gas-fired power plants are frequently seen as the perfect partner for variable renewable generation. However, this reliance on gas generation increases interdependence and propagates unc… ▽ More In many parts of the world, electric power systems have seen a significant shift towards generation from renewable energy and natural gas. Because of their ability to flexibly adjust power generation in real time, gas-fired power plants are frequently seen as the perfect partner for variable renewable generation. However, this reliance on gas generation increases interdependence and propagates uncertainty between power grids and gas pipelines, and brings coordination and uncertainty management challenges. To address these issues, we propose an uncertainty management framework for uncertain, but bounded gas consumption by gas-fired power plants. The admissible ranges are computed based on a joint optimization problem for the combined gas and electricity networks, which involves chance-constrained scheduling for the electric grid and a novel robust optimization formulation for the natural gas network. This formulation ensures feasibility of the integrated system with a high probability, while providing a tractable numerical formulation. A key advance with respect to existing methods is that our method is based on a physically accurate, validated model for transient gas pipeline flows. Our case study benchmarks our proposed formulation against methods that ignore how reserve activation impacts the fuel use of gas power plants, and only consider predetermined gas consumption. The results demonstrate the importance of considering uncertainty to avoid operating constraint violations and curtailment of gas to the generators. △ Less

Submitted 25 June, 2020; originally announced June 2020.

Sequence of Polyhedral Relaxations for Nonlinear Univariate Functions

Authors: Kaarthik Sundar, Sujeevraja Sanjeevi, Harsha Nagarajan

Abstract: Given a nonlinear, univariate, bounded, and differentiable function $f(x)$, this article develops a sequence of Mixed Integer Linear Programming (MILP) and Linear Programming (LP) relaxations that converge to the graph of $f(x)$ and its convex hull, respectively. Theoretical convergence of the sequence of relaxations to the graph of the function and its convex hull is established. For nonlinear no… ▽ More Given a nonlinear, univariate, bounded, and differentiable function $f(x)$, this article develops a sequence of Mixed Integer Linear Programming (MILP) and Linear Programming (LP) relaxations that converge to the graph of $f(x)$ and its convex hull, respectively. Theoretical convergence of the sequence of relaxations to the graph of the function and its convex hull is established. For nonlinear non-convex optimization problems, the relaxations presented in this article can be used to construct tight MILP and LP relaxations. These MILP and the LP relaxations can also be used with MILP-based and spatial branch-and-bound based global optimization algorithms, respectively. △ Less

Submitted 23 February, 2021; v1 submitted 27 May, 2020; originally announced May 2020.

Comments: 15 pages

Report number: LA-UR-20-23980

Out-of-Distribution Detection in Multi-Label Datasets using Latent Space of $β$-VAE

Authors: Vijaya Kumar Sundar, Shreyas Ramakrishna, Zahra Rahiminasab, Arvind Easwaran, Abhishek Dubey

Abstract: Learning Enabled Components (LECs) are widely being used in a variety of perception based autonomy tasks like image segmentation, object detection, end-to-end driving, etc. These components are trained with large image datasets with multimodal factors like weather conditions, time-of-day, traffic-density, etc. The LECs learn from these factors during training, and while testing if there is variati… ▽ More Learning Enabled Components (LECs) are widely being used in a variety of perception based autonomy tasks like image segmentation, object detection, end-to-end driving, etc. These components are trained with large image datasets with multimodal factors like weather conditions, time-of-day, traffic-density, etc. The LECs learn from these factors during training, and while testing if there is variation in any of these factors, the components get confused resulting in low confidence predictions. The images with factors not seen during training is commonly referred to as Out-of-Distribution (OOD). For safe autonomy it is important to identify the OOD images, so that a suitable mitigation strategy can be performed. Classical one-class classifiers like SVM and SVDD are used to perform OOD detection. However, the multiple labels attached to the images in these datasets, restricts the direct application of these techniques. We address this problem using the latent space of the $β$-Variational Autoencoder ($β$-VAE). We use the fact that compact latent space generated by an appropriately selected $β$-VAE will encode the information about these factors in a few latent variables, and that can be used for computationally inexpensive detection. We evaluate our approach on the nuScenes dataset, and our results shows the latent space of $β$-VAE is sensitive to encode changes in the values of the generative factor. △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: Workshop on Assured Autonomy (WAAS) -2020

Piecewise Polyhedral Formulations for a Multilinear Term

Authors: Kaarthik Sundar, Harsha Nagarajan, Jeff Linderoth, Site Wang, Russell Bent

Abstract: In this paper, we present a mixed-integer linear programming (MILP) formulation of a piecewise, polyhedral relaxation (PPR) of a multilinear term using its convex hull representation. Based on the solution of the PPR, we also present a MILP formulation whose solutions are feasible for nonconvex, multilinear equations. We then present computational results showing the effectiveness of proposed form… ▽ More In this paper, we present a mixed-integer linear programming (MILP) formulation of a piecewise, polyhedral relaxation (PPR) of a multilinear term using its convex hull representation. Based on the solution of the PPR, we also present a MILP formulation whose solutions are feasible for nonconvex, multilinear equations. We then present computational results showing the effectiveness of proposed formulations on instances of standard benchmarks of nonlinear programs (NLPs) with multilinear terms and compare the proposed formulation with a traditional formulation that is built by recursively relaxing bilinear groupings of multilinear terms. △ Less

Submitted 3 December, 2020; v1 submitted 2 January, 2020; originally announced January 2020.

Dynamic State and Parameter Estimation for Natural Gas Networks using Real Pipeline System Data

Authors: Kaarthik Sundar, Anatoly Zlotnik

Abstract: We present a method for joint state and parameter estimation for natural gas networks where gas pressures and flows through a network of pipes depend on time-varying injections, withdrawals, and compression controls. The estimation is posed as an optimal control problem constrained by coupled partial differential equations on each pipe that describe space- and time-dependent density and mass flux.… ▽ More We present a method for joint state and parameter estimation for natural gas networks where gas pressures and flows through a network of pipes depend on time-varying injections, withdrawals, and compression controls. The estimation is posed as an optimal control problem constrained by coupled partial differential equations on each pipe that describe space- and time-dependent density and mass flux. These are discretized and combined with nodal conditions to give dynamic constraints for posing the estimation as nonlinear least squares problems. We develop a rapid, scalable computational method for performing the estimation in the presence of measurement and process noise. Finally, we evaluate its effectiveness using a data set from a capacity planning model for an actual pipeline system and a month of time-series data from its supervisory control and data acquisition (SCADA) system. △ Less

Submitted 5 December, 2019; originally announced December 2019.

Comments: In Proc. 3rd IEEE Conference on Control Technology and Applications. arXiv admin note: substantial text overlap with arXiv:1803.07156

Report number: LA-UR-17-27462

A Branch-and-Price Algorithm for a Team Orienteering Problem with Fixed-Wing Drones

Authors: Kaarthik Sundar, Sujeevraja Sanjeevi, Christopher Montez

Abstract: This paper formulates a team orienteering problem with multiple fixed-wing drones and develops a branch-and-price algorithm to solve the problem to optimality. Fixed-wing drones, unlike rotary drones, have kinematic constraints associated with them, thereby preventing them to make on-the-spot turns and restricting them to a minimum turn radius. This paper presents the implications of these constra… ▽ More This paper formulates a team orienteering problem with multiple fixed-wing drones and develops a branch-and-price algorithm to solve the problem to optimality. Fixed-wing drones, unlike rotary drones, have kinematic constraints associated with them, thereby preventing them to make on-the-spot turns and restricting them to a minimum turn radius. This paper presents the implications of these constraints on the drone routing problem formulation and proposes a systematic technique to address them in the context of the team orienteering problem. Furthermore, a novel branch-and-price algorithm with branching techniques specific to the constraints imposed due to fixed-wing drones are proposed. Extensive computational experiments on benchmark instances corroborating the effectiveness of the algorithms are also presented. △ Less

Submitted 1 December, 2021; v1 submitted 9 December, 2019; originally announced December 2019.

Comments: 24 pages including references

Report number: LA-UR-19-32483

Optimal Control for Scheduling and Pricing Intra-day Natural Gas Transport on Pipeline Networks

Authors: Anatoly Zlotnik, Kaarthik Sundar, Aleksandr M. Rudkevich, Aleksandr Beylin, Xindi Li

Abstract: We formulate an economic optimal control problem for transport of natural gas over a large-scale transmission pipeline network under transient flow conditions. The objective is to maximize economic welfare for users of the pipeline system, who provide time-dependent price and quantity bids to purchase or supply gas at metered locations on a system with time-varying injections, withdrawals, and con… ▽ More We formulate an economic optimal control problem for transport of natural gas over a large-scale transmission pipeline network under transient flow conditions. The objective is to maximize economic welfare for users of the pipeline system, who provide time-dependent price and quantity bids to purchase or supply gas at metered locations on a system with time-varying injections, withdrawals, and control actions of compressors and regulators. Our formulation ensures that pipeline hydraulic limitations, compressor station constraints, operational factors, and pre-existing contracts for gas transport are satisfied. A pipeline is modeled as a metric graph with gas dynamics partial differential equations on edges and coupling conditions at the nodes. These dynamic constraints are reduced using lumped elements to a sparse nonlinear differential algebraic equation system. A highly efficient temporal discretization scheme for time-periodic formulations is introduced, which we extend to develop a rolling-horizon model-predictive control scheme. We apply the computational methodology to a pipeline system test network case study. In addition to the physical flow and compressor control solution, the optimization yields dual functions that we interpret as the time-dependent economic values of gas at each location in the network. △ Less

Submitted 5 December, 2019; originally announced December 2019.

Comments: To appear in Proc. 58th IEEE Conference on Decision and Control

Report number: LA-UR-19-28750

Proving Global Optimality of ACOPF Solutions

Authors: S. Gopinath, H. L. Hijazi, T. Weisser, H. Nagarajan, M. Yetkin, K. Sundar, R. W. Bent

Abstract: We present our latest contributions in terms of mathematical modeling and algorithm development for the global optimization of the ACOPF problem. These contributions allow us to close the optimality gap on a number of open instances in the PGLIB and NESTA benchmark libraries. This is achieved by combining valid cut generation with semidefinite programming-based bound tightening. The mathematical f… ▽ More We present our latest contributions in terms of mathematical modeling and algorithm development for the global optimization of the ACOPF problem. These contributions allow us to close the optimality gap on a number of open instances in the PGLIB and NESTA benchmark libraries. This is achieved by combining valid cut generation with semidefinite programming-based bound tightening. The mathematical formulations along with the solution algorithms are implemented in the modeling framework Gravity (www.gravityopt.com), an open-source platform for reproducible numerical experiments. △ Less

Submitted 19 February, 2020; v1 submitted 8 October, 2019; originally announced October 2019.

An Adversarial Model for Attack Vector Vulnerability Analysis on Power and Gas Delivery Operations

Authors: Ignacio Losada Carreno, Anna Scaglione, Anatoly Zlotnik, Deepjyoti Deka, Kaarthik Sundar

Abstract: Power systems often rely on natural gas pipeline networks to supply fuel for gas-fired generation. Market inefficiencies and a lack of formal coordination between the wholesale power and gas delivery infrastructures may magnify the broader impact of a cyber-attack on a natural gas pipeline. In this study we present a model that can be used to quantify the impact of cyber-attacks on electricity and… ▽ More Power systems often rely on natural gas pipeline networks to supply fuel for gas-fired generation. Market inefficiencies and a lack of formal coordination between the wholesale power and gas delivery infrastructures may magnify the broader impact of a cyber-attack on a natural gas pipeline. In this study we present a model that can be used to quantify the impact of cyber-attacks on electricity and gas delivery operations. We model activation of cyber-attack vectors that attempt to gain access to pipeline gas compressor controls using a continuous-time Markov chain over a state space based on the gas operator Industrial Control System firewall zone partition. Our approach evaluates the operating states and decision-making in the networks using physically realistic and operationally representative models. We summarize these models, the sequence of analyses used to quantify the impacts of a cyber-incident, and propose a Monte Carlo simulation approach to quantify the resulting effect on the reliability of the bulk power system by the increase in operational cost. The methodology is applied to a case study of interacting power, gas, and cyber test networks. △ Less

Submitted 8 October, 2019; originally announced October 2019.

Credible Interdiction for Transmission Systems

Authors: Kaarthik Sundar, Sidhant Misra, Russell Bent, Feng Pan

Abstract: This paper presents novel formulations and algorithms for $N$-$k$ interdiction problem in transmission networks. In particular, it formulates spatial and topological resource constraints on attackers for $N$-$k$ interdiction problems and illustrates the formulation with two new classes of $N$-$k$ attacks: (i) Spatial $N$-$k$ attacks where the attack is constrained by geographic distance of a bus c… ▽ More This paper presents novel formulations and algorithms for $N$-$k$ interdiction problem in transmission networks. In particular, it formulates spatial and topological resource constraints on attackers for $N$-$k$ interdiction problems and illustrates the formulation with two new classes of $N$-$k$ attacks: (i) Spatial $N$-$k$ attacks where the attack is constrained by geographic distance of a bus chosen by an attacker and (ii) Topological $N$-$k$ attacks where the attack is constrained to connected components. These two specific types of $N$-$k$ attacks compute interdiction plans designed to better model localized attacks, such as those induced by natural disasters or physical attacks. We then formulate these two resource-constrained interdiction problems as bilevel, max-min optimization problems and present a novel constraint generation algorithm to solve these formulations. Detailed case studies analyzing the behavior of spatially and topologically resource-constrained problems and comparing them to the traditional $N$-$k$ interdiction problem are also presented. △ Less

Submitted 30 November, 2020; v1 submitted 17 April, 2019; originally announced April 2019.

Report number: LA-UR-19-23362

Routing Multiple Unmanned Vehicles in GPS-Denied Environments

Authors: Bingyu Wang, Sohum Misra, Sivakumar Rathinam, Rajnikant Sharma, Kaarthik Sundar

Abstract: This article aims to develop novel path planning algorithms required to deploy multiple unmanned vehicles in Global Positioning System (GPS) denied environments. Unmanned vehicles (ground or aerial) are ideal platforms for executing monitoring and data gathering tasks in civil infrastructure management, agriculture, public safety, law enforcement, disaster relief and transportation. Significant ad… ▽ More This article aims to develop novel path planning algorithms required to deploy multiple unmanned vehicles in Global Positioning System (GPS) denied environments. Unmanned vehicles (ground or aerial) are ideal platforms for executing monitoring and data gathering tasks in civil infrastructure management, agriculture, public safety, law enforcement, disaster relief and transportation. Significant advancement in the area of path planning for unmanned vehicles over the last decade has resulted in a suite of algorithms that can handle heterogeneity, motion and other on-board resource constraints for these vehicles. However, most of these routing and path planning algorithms rely on the availability of the GPS information. Unintentional and intentional interference and design errors can cause GPS service outages, which in turn, can crucially affect all the systems that depend on GPS information. This article addresses a multiple vehicle path planning problem that arises while deploying a team of unmanned vehicles for monitoring applications in GPS-denied environments and presents a mathematical formulation and algorithms for solving the problem. Simulation results are also presented to corroborate the performance of the proposed algorithms. △ Less

Submitted 31 December, 2018; originally announced January 2019.

Comments: 10 pages. arXiv admin note: text overlap with arXiv:1708.03269

Optimization-Based Bound Tightening using a Strengthened QC-Relaxation of the Optimal Power Flow Problem

Authors: Kaarthik Sundar, Harsha Nagarajan, Sidhant Misra, Mowen Lu, Carleton Coffrin, Russell Bent

Abstract: This article develops a strengthened convex quadratic convex (QC) relaxation of the AC Optimal Power Flow (AC-OPF) problem and presents an optimization-based bound-tightening (OBBT) algorithm to compute tight, feasible bounds on the voltage magnitude variables for each bus and the phase angle difference variables for each branch in the network. Theoretical properties of the strengthened QC relaxat… ▽ More This article develops a strengthened convex quadratic convex (QC) relaxation of the AC Optimal Power Flow (AC-OPF) problem and presents an optimization-based bound-tightening (OBBT) algorithm to compute tight, feasible bounds on the voltage magnitude variables for each bus and the phase angle difference variables for each branch in the network. Theoretical properties of the strengthened QC relaxation that show its dominance over the other variants of the QC relaxation studied in the literature are also derived. The effectiveness of the strengthened QC relaxation is corroborated via extensive numerical results on benchmark AC-OPF test networks. In particular, the results demonstrate that the proposed relaxation consistently provides the tightest variable bounds and optimality gaps with negligible impacts on runtime performance. △ Less

Submitted 29 January, 2019; v1 submitted 12 September, 2018; originally announced September 2018.

Estimation of Location and Orientation for Underwater Vehicles from Range Measurements

Authors: Sai Krishna Kanth Hari, Kaarthik Sundar, Jose Braga, Joao Teixeira, Joao Sousa, Swaroop Darbha

Abstract: Localization is an important required task for enabling vehicle autonomy for underwater vehicles. Localization entails the determination of position of the center of mass and orientation of a vehicle from the available measurements. In this paper, we focus on localization by using One-Way Travel Time (OWTT) measurements available to a vehicle from the communication of its multiple on-board receive… ▽ More Localization is an important required task for enabling vehicle autonomy for underwater vehicles. Localization entails the determination of position of the center of mass and orientation of a vehicle from the available measurements. In this paper, we focus on localization by using One-Way Travel Time (OWTT) measurements available to a vehicle from the communication of its multiple on-board receivers with acoustic beacons, more specifically, long baseline (LBL) beacons. Range can be inferred by multiplying the OWTT with speed of sound; however, water conditions can change spatially and temporally resulting in uncertainty in range measurement. The farther a beacon is from a receiver, the larger is the uncertainty. The proposed method for localization accounts captures this uncertainty by bounding the true distance with an increasing (calibrating) function of the range measurement. Determination of this calibration function is formulated as polynomial optimization problem and is a crucial step for localization. The proposed two-step procedure for localization is as follows: based on the range measurements specific to a receiver from the beacons, a convex optimization problem is proposed to estimate the location of the receiver. The estimate is essentially a center of the set of possible locations of the receiver. In the second step, the location estimates of the vehicle are corrected using rigid body motion constraints and the orientation of the rigid body is thus determined. Numerical examples and experimental results provided at the end corroborate the procedures developed in this paper. △ Less

Submitted 1 July, 2019; v1 submitted 9 August, 2018; originally announced August 2018.

Comments: 15 pages

Convex Hull Formulations for Mixed-Integer Multilinear Functions

Authors: Harsha Nagarajan, Kaarthik Sundar, Hassan Hijazi, Russell Bent

Abstract: In this paper, we present convex hull formulations for a mixed-integer, multilinear term/function (MIMF) that features products of multiple continuous and binary variables. We develop two equivalent convex relaxations of an MIMF and study their polyhedral properties in their corresponding higher-dimensional spaces. We numerically observe that the proposed formulations consistently perform better t… ▽ More In this paper, we present convex hull formulations for a mixed-integer, multilinear term/function (MIMF) that features products of multiple continuous and binary variables. We develop two equivalent convex relaxations of an MIMF and study their polyhedral properties in their corresponding higher-dimensional spaces. We numerically observe that the proposed formulations consistently perform better than state-of-the-art relaxation approaches. △ Less

Submitted 29 July, 2018; originally announced July 2018.

Comparing Respiratory Monitoring Performance of Commercial Wireless Devices

Authors: Peter Hillyard, Anh Luong, Alemayehu Solomon Abrar, Neal Patwari, Krishna Sundar, Robert Farney, Jason Burch, Christina A. Porucznik, Sarah Hatch Pollard

Abstract: This paper addresses the performance of systems which use commercial wireless devices to make bistatic RF channel measurements for non-contact respiration sensing. Published research has typically presented results from short controlled experiments on one system. In this paper, we deploy an extensive real-world comparative human subject study. We observe twenty patients during their overnight slee… ▽ More This paper addresses the performance of systems which use commercial wireless devices to make bistatic RF channel measurements for non-contact respiration sensing. Published research has typically presented results from short controlled experiments on one system. In this paper, we deploy an extensive real-world comparative human subject study. We observe twenty patients during their overnight sleep (a total of 160 hours), during which contact sensors record ground-truth breathing data, patient position is recorded, and four different RF breathing monitoring systems simultaneously record measurements. We evaluate published methods and algorithms. We find that WiFi channel state information measurements provide the most robust respiratory rate estimates of the four RF systems tested. However, all four RF systems have periods during which RF-based breathing estimates are not reliable. △ Less

Submitted 17 July, 2018; originally announced July 2018.

Cooperative Planning for Fuel-constrained Aerial Vehicles and Ground-based Refueling Vehicles for Large-Scale Coverage

Authors: Parikshit Maini, Kaarthik Sundar, Sivakumar Rathinam, PB Sujit

Abstract: Low cost Unmanned Aerial Vehicles (UAVs) need multiple refuels to accomplish large area coverage. The number of refueling stations and their placement plays a vital role in determining coverage efficiency. In this paper, we propose the use of a ground-based refueling vehicle (RV) to increase the operational range of a UAV in both spatial and temporal domains. Determining optimal routes for the UAV… ▽ More Low cost Unmanned Aerial Vehicles (UAVs) need multiple refuels to accomplish large area coverage. The number of refueling stations and their placement plays a vital role in determining coverage efficiency. In this paper, we propose the use of a ground-based refueling vehicle (RV) to increase the operational range of a UAV in both spatial and temporal domains. Determining optimal routes for the UAV and RV, and selecting optimized locations for refueling to aid in minimizing coverage time is a challenging problem due to different vehicle speeds, coupling between refueling location placement, and the coverage area at each location. We develop a two-stage strategy for coupled route planning for UAV and RV to perform a coverage mission. The first stage computes a minimal set of refueling sites that permit a feasible UAV route. In the second stage, multiple Mixed-Integer Linear Programming (MILP) formulations are developed to plan optimal routes for the UAV and the refueling vehicle taking into account the feasible set of refueling sites generated in stage one. The performance of different formulations is compared empirically. In addition, computationally efficient heuristics are developed to solve the routing problem. Extensive simulations are conducted to corroborate the effectiveness of proposed approaches. △ Less

Submitted 11 May, 2018; originally announced May 2018.

Cooperative Routing for an Air-Ground Vehicle Team -- Exact Algorithm, Transformation Method, and Heuristics

Authors: Satyanarayana G. Manyam, Kaarthik Sundar, David W. Casbeer

Abstract: This article considers a cooperative vehicle routing problem for an intelligence, surveillance, and reconnaissance mission in the presence of communication constraints between the vehicles. The proposed framework uses a ground vehicle and an Unmanned Aerial Vehicle (UAV) that travel cooperatively and visit a set of targets while satisfying the communication constraints. The problem is formulated a… ▽ More This article considers a cooperative vehicle routing problem for an intelligence, surveillance, and reconnaissance mission in the presence of communication constraints between the vehicles. The proposed framework uses a ground vehicle and an Unmanned Aerial Vehicle (UAV) that travel cooperatively and visit a set of targets while satisfying the communication constraints. The problem is formulated as a mixed-integer linear program, and a branch-and-cut algorithm is developed to solve the problem to optimality. Furthermore, a transformation method and a heuristic are also developed for the problem. The effectiveness of all the algorithms is corroborated through extensive computational experiments on several randomly generated instances. △ Less

Submitted 28 February, 2019; v1 submitted 23 April, 2018; originally announced April 2018.

Comments: 12 pages. arXiv admin note: text overlap with arXiv:1605.09739

State and Parameter Estimation for Natural Gas Pipeline Networks using Transient State Data

Authors: Kaarthik Sundar, Anatoly Zlotnik

Abstract: We formulate two estimation problems for pipeline systems in which measurements of compressible gas flow through a network of pipes is affected by time-varying injections, withdrawals, and compression. We consider a state estimation problem that is then extended to a joint state and parameter estimation problem that can be used for data assimilation. In both formulations, the flow dynamics are des… ▽ More We formulate two estimation problems for pipeline systems in which measurements of compressible gas flow through a network of pipes is affected by time-varying injections, withdrawals, and compression. We consider a state estimation problem that is then extended to a joint state and parameter estimation problem that can be used for data assimilation. In both formulations, the flow dynamics are described on each pipe by space- and time-dependent density and mass flux that evolve according to a system of coupled partial differential equations, in which momentum dissipation is modelled using the Darcy-Wiesbach friction approximation. These dynamics are first spatially discretized to obtain a system of nonlinear ordinary differential equations on which state and parameter estimation formulations are given as nonlinear least squares problems. A rapid, scalable computational method for performing a nonlinear least squares estimation is developed. Extensive simulations and computational experiments on multiple pipeline test networks demonstrate the effectiveness of the formulations in obtaining state and parameter estimates in the presence of measurement and process noise. △ Less

Submitted 5 July, 2018; v1 submitted 19 March, 2018; originally announced March 2018.

Comments: Published in IEEE Transactions on Control Systems Technology

Hierarchical Predictive Control Algorithms for Optimal Design and Operation of Microgrids

Authors: Sai Krishna Kanth Hari, Kaarthik Sundar, Harsha Nagarajan, Russell Bent, Scott Backhaus

Abstract: In recent years, microgrids, i.e., disconnected distribution systems, have received increasing interest from power system utilities to support the economic and resiliency posture of their systems. The economics of long distance transmission lines prevent many remote communities from connecting to bulk transmission systems and these communities rely on off-grid microgrid technology. Furthermore, co… ▽ More In recent years, microgrids, i.e., disconnected distribution systems, have received increasing interest from power system utilities to support the economic and resiliency posture of their systems. The economics of long distance transmission lines prevent many remote communities from connecting to bulk transmission systems and these communities rely on off-grid microgrid technology. Furthermore, communities that are connected to the bulk transmission system are investigating microgrid technologies that will support their ability to disconnect and operate independently during extreme events. In each of these cases, it is important to develop methodologies that support the capability to design and operate microgrids in the absence of transmission over long periods of time. Unfortunately, such planning problems tend to be computationally difficult to solve and those that are straightforward to solve often lack the modeling fidelity that inspires confidence in the results. To address these issues, we first develop a high fidelity model for design and operations of a microgrid that include component efficiencies, component operating limits, battery modeling, unit commitment, capacity expansion, and power flow physics; the resulting model is a mixed-integer quadratically-constrained quadratic program (MIQCQP). We then develop an iterative algorithm, referred to as the Model Predictive Control (MPC) algorithm, that allows us to solve the resulting MIQCQP. We show, through extensive computational experiments, that the MPC-based method can scale to problems that have a very long planning horizon and provide high quality solutions that lie within 5\% of optimal. △ Less

Submitted 18 March, 2018; originally announced March 2018.

Comments: To appear in "Power Systems Computation Conference", Dublin, Ireland

Landmark Placement for Localization in a GPS-denied Environment

Authors: Kaarthik Sundar, Shriram Srinivasan, Sohum Misra, Sivakumar Rathinam, Rajnikant Sharma

Abstract: Path planning algorithms for unmanned aerial or ground vehicles, in many surveillance applications, rely on Global Positioning System (GPS) information for localization. However, disruption of GPS signals, by intention or otherwise, can render these plans and algorithms ineffective. This article provides a way of addressing this issue by utilizing stationary landmarks to aid localization in such G… ▽ More Path planning algorithms for unmanned aerial or ground vehicles, in many surveillance applications, rely on Global Positioning System (GPS) information for localization. However, disruption of GPS signals, by intention or otherwise, can render these plans and algorithms ineffective. This article provides a way of addressing this issue by utilizing stationary landmarks to aid localization in such GPS-disrupted or GPS-denied environment. In particular, given the vehicle's path, we formulate a landmark-placement problem and present algorithms to place the minimum number of landmarks while satisfying the localization, sensing, and collision-avoidance constraints. The performance of such a placement is also evaluated via extensive simulations on ground robots. △ Less

Submitted 10 July, 2018; v1 submitted 21 February, 2018; originally announced February 2018.

Comments: 7 pages, accepted to American Control Conference 2018

A Two-Stage Approach for Routing Multiple Unmanned Aerial Vehicles with Stochastic Fuel Consumption

Authors: Saravanan Venkatachalam, Kaarthik Sundar, Sivakumar Rathinam

Abstract: The past decade has seen a substantial increase in the use of small unmanned aerial vehicles (UAVs) in both civil and military applications. This article addresses an important aspect of refueling in the context of routing multiple small UAVs to complete a surveillance or data collection mission. Specifically, this article formulates a multiple-UAV routing problem with the refueling constraint of… ▽ More The past decade has seen a substantial increase in the use of small unmanned aerial vehicles (UAVs) in both civil and military applications. This article addresses an important aspect of refueling in the context of routing multiple small UAVs to complete a surveillance or data collection mission. Specifically, this article formulates a multiple-UAV routing problem with the refueling constraint of minimizing the overall fuel consumption for all of the vehicles as a two-stage stochastic optimization problem with uncertainty associated with the fuel consumption of each vehicle. The two-stage model allows for the application of sample average approximation (SAA). Although the SAA solution asymptotically converges to the optimal solution for the two-stage model, the SAA run time can be prohibitive for medium- and large-scale test instances. Hence, we develop a tabu-search-based heuristic that exploits the model structure while considering the uncertainty in fuel consumption. Extensive computational experiments corroborate the benefits of the two-stage model compared to a deterministic model and the effectiveness of the heuristic for obtaining high-quality solutions. △ Less

Submitted 24 October, 2018; v1 submitted 13 November, 2017; originally announced November 2017.

Comments: 18 pages

PowerModels.jl: An Open-Source Framework for Exploring Power Flow Formulations

Authors: Carleton Coffrin, Russell Bent, Kaarthik Sundar, Yeesian Ng, Miles Lubin

Abstract: In recent years, the power system research community has seen an explosion of novel methods for formulating and solving power network optimization problems. These emerging methods range from new power flow approximations, which go beyond the traditional DC power flow by capturing reactive power, to convex relaxations, which provide solution quality and runtime performance guarantees. Unfortunately… ▽ More In recent years, the power system research community has seen an explosion of novel methods for formulating and solving power network optimization problems. These emerging methods range from new power flow approximations, which go beyond the traditional DC power flow by capturing reactive power, to convex relaxations, which provide solution quality and runtime performance guarantees. Unfortunately, the sophistication of these emerging methods often presents a significant barrier to evaluating them on a wide variety of power system optimization applications. To address this issue, this work proposes PowerModels, an open-source platform for comparing power flow formulations. From its inception, PowerModels was designed to streamline the process of evaluating different power flow formulations on shared optimization problem specifications. This work provides a brief introduction to the design of PowerModels, validates its implementation, and demonstrates its effectiveness with a proof-of-concept study analyzing five different formulations of the Optimal Power Flow problem. △ Less

Submitted 12 March, 2018; v1 submitted 5 November, 2017; originally announced November 2017.

Report number: LA-UR-17-29326

Relaxations of AC Maximal Load Delivery for Severe Contingency Analysis

Authors: Carleton Coffrin, Russell Bent, Byron Tasseff, Kaarthik Sundar, Scott Backhaus

Abstract: This work considers the task of finding an AC-feasible operating point of a severely damaged transmission network while ensuring that a maximal amount of active power loads can be delivered. This AC Maximal Load Delivery (AC-MLD) task is a nonconvex nonlinear optimization problem that is incredibly challenging to solve on large-scale transmission system datasets. This work demonstrates that convex… ▽ More This work considers the task of finding an AC-feasible operating point of a severely damaged transmission network while ensuring that a maximal amount of active power loads can be delivered. This AC Maximal Load Delivery (AC-MLD) task is a nonconvex nonlinear optimization problem that is incredibly challenging to solve on large-scale transmission system datasets. This work demonstrates that convex relaxations of the AC-MLD problem provide a reliable and scalable method for finding high-quality bounds on the amount of active power that can be delivered in the AC-MLD problem. To demonstrate their effectiveness, the solution methods proposed in this work are rigorously evaluated on 1000 N-k scenarios on seven power networks ranging in size from 70 to 6000 buses. The most effective relaxation of the AC-MLD problem converges in less than 20 seconds on commodity computing hardware for all 7000 of the scenarios considered. △ Less

Submitted 14 October, 2018; v1 submitted 21 October, 2017; originally announced October 2017.

Comments: The problem considered in this work is also known as AC Minimal Load-Shedding

Report number: LA-UR-17-29379, LA-UR-18-29573

Routing Unmanned Vehicles in GPS-Denied Environments

Authors: Kaarthik Sundar, Sohum Misra, Sivakumar Rathinam, Rajnikant Sharma

Abstract: Most of the routing algorithms for unmanned vehicles, that arise in data gathering and monitoring applications in the literature, rely on the Global Positioning System (GPS) information for localization. However, disruption of GPS signals either intentionally or unintentionally could potentially render these algorithms not applicable. In this article, we present a novel method to address this diff… ▽ More Most of the routing algorithms for unmanned vehicles, that arise in data gathering and monitoring applications in the literature, rely on the Global Positioning System (GPS) information for localization. However, disruption of GPS signals either intentionally or unintentionally could potentially render these algorithms not applicable. In this article, we present a novel method to address this difficulty by combining methods from cooperative localization and routing. In particular, the article formulates a fundamental combinatorial optimization problem to plan routes for an unmanned vehicle in a GPS-restricted environment while enabling localization for the vehicle. We also develop algorithms to compute optimal paths for the vehicle using the proposed formulation. Extensive simulation results are also presented to corroborate the effectiveness and performance of the proposed formulation and algorithms. △ Less

Submitted 19 December, 2017; v1 submitted 10 August, 2017; originally announced August 2017.

Comments: Publised in International Conference on Umanned Aerial Systems

An Adaptive, Multivariate Partitioning Algorithm for Global Optimization of Nonconvex Programs

Authors: Harsha Nagarajan, Mowen Lu, Site Wang, Russell Bent, Kaarthik Sundar

Abstract: In this work, we develop an adaptive, multivariate partitioning algorithm for solving mixed-integer nonlinear programs (MINLP) with multi-linear terms to global optimality. This iterative algorithm primarily exploits the advantages of piecewise polyhedral relaxation approaches via disjunctive formulations to solve MINLPs to global optimality in contrast to the conventional spatial branch-and-bound… ▽ More In this work, we develop an adaptive, multivariate partitioning algorithm for solving mixed-integer nonlinear programs (MINLP) with multi-linear terms to global optimality. This iterative algorithm primarily exploits the advantages of piecewise polyhedral relaxation approaches via disjunctive formulations to solve MINLPs to global optimality in contrast to the conventional spatial branch-and-bound approaches. In order to maintain relatively small-scale mixed-integer linear programs at every iteration of the algorithm, we adaptively partition the variable domains appearing in the multi-linear terms. We also provide proofs on convergence guarantees of the proposed algorithm to a global solution. Further, we discuss a few algorithmic enhancements based on the sequential bound-tightening procedure as a presolve step, where we observe the importance of solving piecewise relaxations compared to basic convex relaxations to speed-up the convergence of the algorithm to global optimality. We demonstrate the effectiveness of our disjunctive formulations and the algorithm on well-known benchmark problems (including Pooling and Blending instances) from MINLPLib and compare with state-of-the-art global optimization solvers. With this novel approach, we solve several large-scale instances which are, in some cases, intractable by the global optimization solver. We also shrink the best known optimality gap for one of the hard, generalized pooling problem instance. △ Less

Submitted 1 February, 2019; v1 submitted 8 July, 2017; originally announced July 2017.

Probabilistic $N$-$k$ Failure-Identification for Power Systems

Authors: Kaarthik Sundar, Carleton Coffrin, Harsha Nagarajan, Russell Bent

Abstract: This paper considers a probabilistic generalization of the $N$-$k$ failure-identification problem in power transmission networks, where the probability of failure of each component in the network is known a priori and the goal of the problem is to find a set of $k$ components that maximizes disruption to the system loads weighted by the probability of simultaneous failure of the $k$ components. Th… ▽ More This paper considers a probabilistic generalization of the $N$-$k$ failure-identification problem in power transmission networks, where the probability of failure of each component in the network is known a priori and the goal of the problem is to find a set of $k$ components that maximizes disruption to the system loads weighted by the probability of simultaneous failure of the $k$ components. The resulting problem is formulated as a bilevel mixed-integer nonlinear program. Convex relaxations, linear approximations, and heuristics are developed to obtain feasible solutions that are close to the optimum. A general cutting-plane algorithm is proposed to solve the convex relaxation and linear approximations of the $N$-$k$ problem. Extensive numerical results corroborate the effectiveness of the proposed algorithms on small-, medium-, and large-scale test instances, the test instances include the IEEE 14-bus system, the IEEE single-area and three-area RTS96 systems, the IEEE 118-bus system, the WECC 240-bus test system, the 1354-bus PEGASE system, and the 2383-bus Polish winter-peak test system. △ Less

Submitted 10 July, 2018; v1 submitted 18 April, 2017; originally announced April 2017.

Comments: 17 pages, Networks, 2018