skip to main content


Title: Natural Gas Flow Equations: Uniqueness and an MI-SOCP Solver
The critical role of gas fired-plants to compensate renewable generation has increased the operational variability in natural gas networks (GN). Towards developing more reliable and efficient computational tools for GN monitoring, control, and planning, this work considers the task of solving the nonlinear equations governing steady-state flows and pressures in GNs. It is first shown that if the gas flow equations are feasible, they enjoy a unique solution. To the best of our knowledge, this is the first result proving uniqueness of the steady-state gas flow solution over the entire feasible domain of gas injections. To find this solution, we put forth a mixed-integer second-order cone program (MI-SOCP)-based solver relying on a relaxation of the gas flow equations. This relaxation is provably exact under specific network topologies. Unlike existing alternatives, the devised solver does not need proper initialization or knowing the gas flow directions beforehand, and can handle gas networks with compressors. Numerical tests on tree and meshed networks indicate that the relaxation is exact even when the derived conditions are not met.  more » « less
Award ID(s):
1711587
PAR ID:
10179047
Author(s) / Creator(s):
;
Date Published:
Journal Name:
2019 American Control Conference (ACC)
Page Range / eLocation ID:
2114 to 2120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In this work, we consider two-stage quadratic optimization problems under ellipsoidal uncertainty. In the first stage, one needs to decide upon the values of a subset of optimization variables (control variables). In the second stage, the uncertainty is revealed, and the rest of the optimization variables (state variables) are set up as a solution to a known system of possibly nonlinear equations. This type of problem occurs, for instance, in optimization for dynamical systems, such as electric power systems as well as gas and water networks. We propose a convergent iterative algorithm to build a sequence of approximately robustly feasible solutions with an improving objective value. At each iteration, the algorithm optimizes over a subset of the feasible set and uses affine approximations of the second-stage equations while preserving the nonlinearity of other constraints. We implement our approach and demonstrate its performance on Matpower instances of AC optimal power flow. Although this paper focuses on quadratic problems, the approach is suitable for more general setups. 
    more » « less
  2. This paper focuses on the problem of finding multiple paths with relay nodes to maximize throughput for ultra-high-rate millimeter wave (mmWave) backhaul networks in urban environments. Relays are selected between a pair of source and destination base stations to form multiple interference-free paths. We first formulate the problem of feasibility of multi-path construction as a constraint satisfaction problem that includes constraints on intra-path and inter-path interference and several other constraints that arise from the problem setting. Based on the derived equations, we transform the multiple paths construction problem into a Boolean satisfiability problem. This problem can then be solved through use of a satisfiability (SAT) solver, which however results in a very high running time for realistic problem sizes. To address this, we propose a heuristic algorithm that runs in a fraction of the time of the SAT solver and finds multiple interference-free paths using a modification of a maximum flow algorithm. Simulation results based on 3-D models of a section of downtown Atlanta show that the heuristic algorithm finds multiple paths in almost all the feasible cases (those where the SAT solver succeeds in finding a solution) and produces paths with higher average throughput than the SAT solver. Furthermore, the heuristic increases throughput by 50-100% in typical cases compared to a single-path solution. 
    more » « less
  3. The behaviour of low-viscosity, pressure-driven compressible pore fluid flows in viscously deformable porous media is studied here with specific application to gas flow in lava domes. The combined flow of gas and lava is shown to be governed by a two-equation set of nonlinear mixed hyperbolic–parabolic type partial differential equations describing the evolution of gas pore pressure and lava porosity. Steady state solution of this system is achieved when the gas pore pressure is magmastatic and the porosity profile accommodates the magmastatic pressure condition by increased compaction of the medium with depth. A one-dimensional (vertical) numerical linear stability analysis (LSA) is presented here. As a consequence of the pore-fluid compressibility and the presence of gravitation compaction, the gradients present in the steady-state solution cause variable coefficients in the linearized equations which generate instability in the LSA despite the diffusion-like and dissipative terms in the original system. The onset of this instability is shown to be strongly controlled by the thickness of the flow and the maximum porosity, itself a function of the mass flow rate of gas. Numerical solutions of the fully nonlinear system are also presented and exhibit nonlinear wave propagation features such as shock formation. As applied to gas flow within lava domes, the details of this dynamics help explain observations of cyclic lava dome extrusion and explosion episodes. Because the instability is stronger in thicker flows, the continued extrusion and thickening of a lava dome constitutes an increasing likelihood of instability onset, pressure wave growth and ultimately explosion. 
    more » « less
  4. Emerging distribution systems with a proliferation of distributed energy resources are facing with new challenges, such as voltage collapse and power flow congestion in unsymmetrical network configurations. As a fundamental tool that could help quantify these new challenges and further mitigate their impacts on the secure and economic operation of distribution systems, effective AC optimal power flow (ACOPF) models and solution approaches are in urgent need. This study focuses on ACOPF of three‐phase four‐conductor configured distribution systems, in which neutral conductors and ground resistances are modelled explicitly to reflect practical situation. In addition, by leveraging the Kirchhoff's current law (KCL) theorem and the effect of zero injections, voltage variables of neutrals and zero‐injection phases can be effectively eliminated. The ACOPF problem is formulated as a convex semidefinite programming (SDP) relaxation model in complex domain. In recognising possible solution inexactness of SDP relaxation model, a Karush–Kuhn–Tucker condition based process is further proposed to effectively recover feasible solutions to the original ACOPF problem by calculating a set of computational‐inexpensive non‐linear equations. Numerical studies on a modified IEEE 123‐bus system show the effectiveness of the proposed SDP relaxation model with variable reductions and the feasible solution recovery process for three‐phase four‐conductor configured distribution systems.

     
    more » « less
  5. As the number of distributed energy resources participating in power networks increases, it becomes increasingly more important to actively manage network constraints to ensure safe operation. One proposed method that has gained significant attention and implementation, particularly in Australia, is the use of dynamic operating envelopes. Operating envelopes represent net export limits set by the system operator on every node in the distribution network that change as system conditions change. They are calculated using an optimal power flow problem, frequently using a linearization or relaxation of the nonlinear power flow equations. This paper presents two case studies and some numerical analysis to explain why a second order cone relaxation of the power flow equations will lead to ineffective operating envelopes. A modification to the objective function which allows the second order cone relaxation to nearly recover the solution to the nonlinear formulation is also presented. 
    more » « less