Search

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.

We develop an explicit second order staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge domain boundary values that are significantly more complex than in the case of a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

We present a computational method for open-loop minimum-norm control synthesis for fixed-endpoint transfer of bilinear ensemble systems that are indexed by two continuously varying parameters. We suppose that one ensemble parameter scales the homogeneous, linear part of the dynamics, and the second parameter scales the effect of the applied control inputs on the inhomogeneous, bilinear dynamics. This class of dynamical systems is motivated by robust quantum control pulse synthesis, where the ensemble parameters correspond to uncertainty in the free Hamiltonian and inhomogeneity in the control Hamiltonian, respectively. Our computational method is based on polynomial approximation of the ensemble state in parameter space and discretization of the evolution equations in the time domain using a product of matrix exponentials corresponding to zero-order hold controls over the time intervals. The dynamics are successively linearized about control and trajectory iterates to formulate a sequence of quadratic programs for computing perturbations to the control that successively improve the objective until the iteration converges. We use a two-stage computation to first ensure transfer to the desired terminal state, and then minimize the norm of the control function. The method is demonstrated for the canonical uniform transfer problem for the Bloch system that appears in nuclear magnetic resonance, as well as the matter-wave splitting problem for the Raman-Nath system that appears in ultra-cold atom interferometry.

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.

We propose a predictor-corrector adaptive method for the simulation of hyperbolic partial differential equations (PDEs) on networks under general uncertainty in parameters, initial conditions, or boundary conditions. The approach is based on the stochastic finite volume (SFV) framework that circumvents sampling schemes or simulation ensembles while also preserving fundamental properties, in particular hyperbolicity of the resulting systems and conservation of the discrete solutions. The initial boundary value problem (IBVP) on a set of network-connected one-dimensional domains that represent a pipeline is represented using active discretization of the physical and stochastic spaces, and we evaluate the propagation of uncertainty through network nodes by solving a junction Riemann problem. The adaptivity of our method in refining discretization based on error metrics enables computationally tractable evaluation of intertemporal uncertainty in order to support decisions about timing and quantity of pipeline operations to maximize delivery under transient and uncertain conditions. We illustrate our computational method using simulations for a representative network.

We propose a predictor-corrector adaptive method for the study of hyperbolic partial differential equations (PDEs) under uncertainty. Constructed around the framework of stochastic finite volume (SFV) methods, our approach circumvents sampling schemes or simulation ensembles while also preserving fundamental properties, in particular hyperbolicity of the resulting systems and conservation of the discrete solutions. Furthermore, we augment the existing SFV theory with a priori convergence results for statistical quantities, in particular push-forward densities, which we demonstrate through numerical experiments. By linking refinement indicators to regions of the physical and stochastic spaces, we drive anisotropic refinements of the discretizations, introducing new degrees of freedom (DoFs) where deemed profitable. To illustrate our proposed method, we consider a series of numerical examples for non-linear hyperbolic PDEs based on Burgers' and Euler's equations.

It is known that periodic forcing of nonlinear flows can result in a chaotic response under certain conditions. Such non-periodic and chaotic solutions have been observed in simulations of heterogeneous gas flow in a pipeline with periodic, time-varying boundary conditions. In this paper, we examine a proportional feedback law for boundary control of a parabolic partial differential equation system that represents the flow of two gases through a pipe. We demonstrate that periodic variation of the mass fraction of the lighter gas at the pipe inlet can result in the chaotic propagation of gas pressure waves, and show that appropriate flow control can suppress this response. We examine phase space solutions for the single pipe system subject to boundary control, and use numerical experiments to characterize conditions for the controller gain to suppress chaos.

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.

We present a framework for designing optimal optical pulses for the matter-wave splitting of a Bose-Einstein Condensate (BEC) under the influence of experimental inhomogeneities, so that the sample is transferred from an initial rest position into a singular higher diffraction order. To represent the evolution of the population of atoms, the Schroedinger's equation is reinterpreted as a parameterized ensemble of dynamical units that are disparately impacted by the beam light-shift potential in a continuous manner. The derived infinite-dimensional coupled Raman-Nath equations are truncated to a finite system of diffraction levels, and we suppose that the parameter that defines the inhomogeneity in the control applied to the ensemble system is restricted to a compact interval. We first design baseline square pulse sequences for the excitation of BEC beam-splitter states following a previous study, subject to dynamic constraints for either a nominal system assuming no inhomogeneity or for several samples of the uncertain parameter. We then approximate the continuum state-space of the ensemble of dynamics using a spectral approach based on Legendre moments, which is truncated at a finite order. Control functions that steer the BEC system from an equivalent rest position to a desired final excitation are designed using a constrained optimal control approach developed for handling nonlinear dynamics. This representation results in a minimal dimension of the computational problem and is shown to be highly robust to inhomogeneity in comparison to the baseline approach. Our method accomplishes the BEC-splitting state transfer for each subsystem in the ensemble, and is promising for precise excitation in experimental settings where robustness to environmental and intrinsic noise is paramount.

We derive a linear system of ordinary differential equations (ODEs) to approximate the dynamics of natural gas in pipeline networks. Although a closed-form expression of the eigenvalues of the state matrix does not generally exist, the poles of an irrational transfer function corresponding to the linearized partial differential equations are used to approximate the eigenvalues of the ODE system. Our analysis qualitatively demonstrates that the eigenvalues of the state matrix of the entire network system are "pipeline separable" in the sense that the eigenvalues are dominated by the individual pipeline parameters and not the incidence connectivity of the network graph. The linear system is used as the dynamic constraints of a linear optimal control problem (OCP) to design the control actions of compressor units to minimize the energy that they expend. The motivation of this work is to reduce the computational complexity of optimizing gas dynamics in large networks to meet the unpredictable and highly variable demand from electric generators. The linear and corresponding nonlinear OCPs are discretized in time to obtain linear and nonlinear optimization problems, which are demonstrated on a test network to illustrate the validity of linear programming. Moreover, an analytical bound on the error between the solutions of the linear and nonlinear flow dynamics is presented using Lyapunov functions and verified computationally by plotting the error against the size of the flow variation around the steady-state solution.

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.

The blending of hydrogen generated using clean energy into natural gas pipeline networks is proposed in order to utilize existing energy systems for their planned lifetimes while reducing their reliance on fossil fuels. We formulate a system of partial differential equations (PDEs) that govern the flow dynamics of mixtures of gases in pipeline networks under the influence of time-varying compressor and regulator control actions. The formulation is derived for general gas networks that can inject or withdraw arbitrary time-varying mixtures of gases into or from the network at arbitrarily specified nodes. The PDE formulation is discretized in space to form a nonlinear control system that is used to prove that homogeneous mixtures are well-behaved and heterogeneous mixtures may be ill-behaved in the sense of monotone-ordering of solutions. We use numerical simulations to compute interfaces in the parameter region of sinusoidal boundary conditions that delimit monotonic, periodic, and chaotic system responses. The interfaces suggest that any solution in the monotonic response region is not chaotic and will eventually approach a periodic orbit. The results are demonstrated using examples for a single pipeline and a small test network.

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., Comment: 12 pages

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., Comment: In Proceedings of the 2023 American Control Conference

We formulate a control system model for the distributed flow of mixtures of highly heterogeneous gases through large-scale pipeline networks with time-varying injections of constituents, withdrawals, and control actions of compressors. This study is motivated by the proposed blending of clean hydrogen into natural gas pipelines as an interim means to reducing end use carbon emissions while utilizing existing infrastructure for its planned lifetime. We reformulate the partial differential equations for gas dynamics on pipelines and balance conditions at junctions using lumped elements to a sparse nonlinear differential algebraic equation system. Our key advance is modeling the mixing of constituents in time throughout the network, which requires doubling the state space needed for a single gas and increases numerical ill-conditioning. The reduced model is shown to be a consistent approximation of the original system, which we use as the dynamic constraints in a model-predictive optimal control problem for minimizing the energy expended by applying time-varying compressor operating profiles to guarantee time-varying delivery profiles subject to system pressure limits. The optimal control problem is implemented after time discretization using a nonlinear program, with validation of the results done using a transient simulation. We demonstrate the methodology for a small test network, and discuss scalability and potential applications., Comment: In Proceedings of the 2023 American Control Conference

We develop a weakly intrusive framework to simulate the propagation of uncertainty in solutions of generic hyperbolic partial differential equation systems on graph-connected domains with nodal coupling and boundary conditions. The method is based on the Stochastic Finite Volume (SFV) approach, and can be applied for uncertainty quantification (UQ) of the dynamical state of fluid flow over actuated transport networks. The numerical scheme has specific advantages for modeling intertemporal uncertainty in time-varying boundary parameters, which cannot be characterized by strict upper and lower (interval) bounds. We describe the scheme for a single pipe, and then formulate the controlled junction Riemann problem (JRP) that enables the extension to general network structures. We demonstrate the method's capabilities and performance characteristics using a standard benchmark test network.

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., Comment: 12 pages

The capabilities of natural neural systems have inspired new generations of machine learning algorithms as well as neuromorphic very large-scale integrated (VLSI) circuits capable of fast, low-power information processing. However, it has been argued that most modern machine learning algorithms are not neurophysiologically plausible. In particular, the workhorse of modern deep learning, the backpropagation algorithm, has proven difficult to translate to neuromorphic hardware. In this study, we present a neuromorphic, spiking backpropagation algorithm based on synfire-gated dynamical information coordination and processing, implemented on Intel's Loihi neuromorphic research processor. We demonstrate a proof-of-principle three-layer circuit that learns to classify digits from the MNIST dataset. To our knowledge, this is the first work to show a Spiking Neural Network (SNN) implementation of the backpropagation algorithm that is fully on-chip, without a computer in the loop. It is competitive in accuracy with off-chip trained SNNs and achieves an energy-delay product suitable for edge computing. This implementation shows a path for using in-memory, massively parallel neuromorphic processors for low-power, low-latency implementation of modern deep learning applications., Comment: 21 pages, 5 figures, Changes v1->v2: minor changes of text and formatting, correction of total power in supplementary Table III

We propose a general method for optimizing periodic input waveforms for global entrainment of weakly forced limit-cycle oscillators based on phase reduction and nonlinear programming. We derive averaged phase dynamics from the mathematical model of a limit-cycle oscillator driven by a weak periodic input and optimize the Fourier coefficients of the input waveform to maximize prescribed objective functions. In contrast to the optimization methods that rely on the calculus of variations, the proposed method can be applied to a wider class of optimization problems including global entrainment objectives. As an illustration, we consider two optimization problems, one for achieving fast global convergence of the oscillator to the entrained state and the other for realizing prescribed global phase distributions in a population of identical uncoupled noisy oscillators. We show that the proposed method can successfully yield optimal input waveforms to realize the desired states in both cases., Comment: 17pages, 8 figures

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., Comment: 13 pages

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., Comment: 8 pages

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.

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., Comment: Corrected alignment of histograms in Figure 2

Systematic discovery of reduced-order closure models for multi-scale processes remains an important open problem in complex dynamical systems. Even when an effective lower-dimensional representation exists, reduced models are difficult to obtain using solely analytical methods. Rigorous methodologies for finding such coarse-grained representations of multi-scale phenomena would enable accelerated computational simulations and provide fundamental insights into the complex dynamics of interest. We focus on a heterogeneous population of oscillators of Kuramoto type as a canonical model of complex dynamics, and develop a data-driven approach for inferring its coarse-grained description. Our method is based on a numerical optimization of the coefficients in a general equation of motion informed by analytical derivations in the thermodynamic limit. We show that certain assumptions are required to obtain an autonomous coarse-grained equation of motion. However, optimizing coefficient values enables coarse-grained models with conceptually disparate functional forms, yet comparable quality of representation, to provide accurate reduced-order descriptions of the underlying system., Comment: 13 pages, 10 figures

We derive conditions for monotonicity properties that characterize general flows of a commodity over a network, where the flow is described by potential and flow dynamics on the edges, as well as potential continuity and Kirchhoff-Neumann mass balance requirements at nodes. The transported commodity may be injected or withdrawn at any of the network nodes, and its movement throughout the network is controlled by nodal actuators. For a class of dissipative nonlinear parabolic partial differential equation (PDE) systems on networks, we derive conditions for monotonicity properties in steady-state flow, as well as for propagation of monotone ordering of states with respect to time-varying boundary condition parameters. In the latter case, initial conditions, as well as time-varying parameters in the coupling conditions at vertices, provide an initial boundary value problem (IBVP). We prove that ordering properties of the solution to the IBVP are preserved when the initial conditions and the parameters of the time-varying coupling law are appropriately ordered. Then, we prove that when monotone ordering is not preserved, the first crossing of solutions occurs at a network node. We consider the implications for robust optimization and optimal control formulations and real-time monitoring of uncertain dynamic flows on networks, and discuss application to subsonic compressible fluid flow with energy dissipation on physical networks. The main result and monitoring policy are demonstrated for gas pipeline test networks and a case study using data corresponding to a real working system. We propose applications of this general result to the control and monitoring of natural gas transmission networks., Comment: arXiv admin note: text overlap with arXiv:1601.05102

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.

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., Comment: In Proc. 3rd IEEE Conference on Control Technology and Applications. arXiv admin note: substantial text overlap with arXiv:1803.07156

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., Comment: To appear in Proc. 58th IEEE Conference on Decision and Control

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.

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., Comment: Published in IEEE Transactions on Control Systems Technology

A significant portion of electricity consumed worldwide is used to power thermostatically controlled loads (TCLs) such as air conditioners, refrigerators, and water heaters. Because the short-term timing of operation of such systems is inconsequential as long as their long-run average power consumption is maintained, they are increasingly used in demand response (DR) programs to balance supply and demand on the power grid. Here, we present an \textit{ab initio} phase model for general TCLs, and use the concept to develop a continuous oscillator model of a TCL and compute its phase response to changes in temperature and applied power. This yields a simple control system model that can be used to evaluate control policies for modulating the power consumption of aggregated loads with parameter heterogeneity and stochastic drift. We demonstrate this concept by comparing simulations of ensembles of heterogeneous loads using the continuous state model and an established hybrid state model. The developed phase model approach is a novel means of evaluating DR provision using TCLs, and is instrumental in estimating the capacity of ancillary services or DR on different time scales. We further propose a novel phase response based open-loop control policy that effectively modulates the aggregate power of a heterogeneous TCL population while maintaining load diversity and minimizing power overshoots. This is demonstrated by low-error tracking of a regulation signal by filtering it into frequency bands and using TCL sub-ensembles with duty cycles in corresponding ranges. Control policies that can maintain a uniform distribution of power consumption by aggregated heterogeneous loads will enable distribution system management (DSM) approaches that maintain stability as well as power quality, and further allow more integration of renewable energy sources., Comment: 17 pages, 16 figures

We present an explicit second order staggered finite difference (FD) discretization scheme for forward simulation of natural gas transport in pipeline networks. By construction, this discretization approach guarantees that the conservation of mass condition is satisfied exactly. The mathematical model is formulated in terms of density, pressure, and mass flux variables, and as a result permits the use of a general equation of state to define the relation between the gas density and pressure for a given temperature. In a single pipe, the model represents the dynamics of the density by propagation of a non-linear wave according to a variable wave speed. We derive compatibility conditions for linking domain boundary values to enable efficient, explicit simulation of gas flows propagating through a network with pressure changes created by gas compressors. We compare Kiuchi's implicit method and an explicit operator splitting method with our staggered grid method, and perform numerical experiments to validate the convergence order of the new method. In addition, we perform several computations to investigate the influence of non-ideal equation of state models and temperature effects into pipeline simulations with boundary conditions over various time and space scales., Comment: 28 pages, 11 figures, 5 tables

We examine the design of the entrainment process for an uncountably infinite collection of coupled phase oscillators that are all subject to the same periodic driving signal. In the absence of coupling, an appropriately designed input can result in each oscillator attaining the frequency of the driving signal, with a phase offset determined by its natural frequency. We consider a special case of interacting oscillators in which the coupling tends to destabilize the phase configuration to which the driving signal would send the collection in the absence of coupling. In this setting we derive stability results that characterize the trade-off between the effects of driving and coupling, and compare these results to the well-known Kuramoto model of a collection of free-running coupled oscillators.

We develop an operator splitting method to simulate flows of isothermal compressible natural gas over transmission pipelines. The method solves a system of nonlinear hyperbolic partial differential equations (PDEs) of hydrodynamic type for mass flow and pressure on a metric graph, where turbulent losses of momentum are modeled by phenomenological Darcy-Weisbach friction. Mass flow balance is maintained through the boundary conditions at the network nodes, where natural gas is injected or withdrawn from the system. Gas flow through the network is controlled by compressors boosting pressure at the inlet of the adjoint pipe. Our operator splitting numerical scheme is unconditionally stable and it is second order accurate in space and time. The scheme is explicit, and it is formulated to work with general networks with loops. We test the scheme over range of regimes and network configurations, also comparing its performance with performance of two other state of the art implicit schemes., Comment: 21 pages, 11 figures, submitted to SIAM Journal of Scientific Computing

We derive conditions for the propagation of monotone ordering properties for a class of nonlinear parabolic partial differential equation (PDE) systems on metric graphs. For such systems, PDE equations with a general nonlinear dissipation term define evolution on each edge, and balance laws create Kirchhoff-Neumann boundary conditions at the vertices. Initial conditions, as well as time-varying parameters in the coupling conditions at vertices, provide an initial value problem (IVP). We first prove that ordering properties of the solution to the IVP are preserved when the initial conditions and time-varying coupling law parameters at vertices are appropriately ordered. In addition, we prove that when monotone ordering is not preserved, the first crossing of solutions occurs at a graph vertex. We consider the implications for robust optimal control formulations and real-time monitoring involving uncertain dynamic flows on networks, and discuss application to subsonic compressible fluid flow with energy dissipation on physical networks.

We derive a monotonicity property for general, transient flows of a commodity transferred throughout a network, where the flow is characterized by density and mass flux dynamics on the edges with density continuity and mass balance conditions at the nodes. The dynamics on each edge are represented by a general system of partial differential equations that approximates subsonic compressible fluid flow with energy dissipation. The transferred commodity may be injected or withdrawn at any of the nodes, and is propelled throughout the network by nodally located compressors. These compressors are controllable actuators that provide a means to manipulate flows through the network, which we therefore consider as a control system. A canonical problem requires compressor control protocols to be chosen such that time-varying nodal commodity withdrawal profiles are delivered and the density remains within strict limits while an economic or operational cost objective is optimized. In this manuscript, we consider the situation where each nodal commodity withdrawal profile is uncertain, but is bounded within known maximum and minimum time-dependent limits. We introduce the monotone parameterized control system property, and prove that general dynamic dissipative network flows possess this characteristic under certain conditions. This property facilitates very efficient formulation of optimal control problems for such systems in which the solutions must be robust with respect to commodity withdrawal uncertainty. We discuss several applications in which such control problems arise and where monotonicity enables simplified characterization of system behavior.

The growing reliance of electric power systems on gas-fired generation to balance intermittent sources of renewable energy has increased the variation and volume of flows through natural gas transmission pipelines. Adapting pipeline operations to maintain efficiency and security under these new conditions requires optimization methods that account for transients and that can quickly compute solutions in reaction to generator re-dispatch. This paper presents an efficient scheme to minimize compression costs under dynamic conditions where deliveries to customers are described by time-dependent mass flow. The optimization scheme relies on a compact representation of gas flow physics, a trapezoidal discretization in time and space, and a two-stage approach to minimize energy costs and maximize smoothness. The resulting large-scale nonlinear programs are solved using a modern interior-point method. The proposed optimization scheme is validated against an integration of dynamic equations with adaptive time-stepping, as well as a recently proposed state-of-the-art optimal control method. The comparison shows that the solutions are feasible for the continuous problem and also practical from an operational standpoint. The results also indicate that our scheme provides at least an order of magnitude reduction in computation time relative to the state-of-the-art and scales to large gas transmission networks with more than 6000 kilometers of total pipeline.

We outline a new control system model for the distributed dynamics of compressible gas flow through large-scale pipeline networks with time-varying injections, withdrawals, and control actions of compressors and regulators. The gas dynamics PDE equations over the pipelines, together with boundary conditions at junctions, are reduced using lumped elements to a sparse nonlinear ODE system expressed in vector-matrix form using graph theoretic notation. This system, which we call the reduced network flow (RNF) model, is a consistent discretization of the PDE equations for gas flow. The RNF forms the dynamic constraints for optimal control problems for pipeline systems with known time-varying withdrawals and injections and gas pressure limits throughout the network. The objectives include economic transient compression (ETC) and minimum load shedding (MLS), which involve minimizing compression costs or, if that is infeasible, minimizing the unfulfilled deliveries, respectively. These continuous functional optimization problems are approximated using the Legendre-Gauss-Lobatto (LGL) pseudospectral collocation scheme to yield a family of nonlinear programs, whose solutions approach the optima with finer discretization. Simulation and optimization of time-varying scenarios on an example natural gas transmission network demonstrate the gains in security and efficiency over methods that assume steady-state behavior.

For many natural and engineered systems, a central function or design goal is the synchronization of one or more rhythmic or oscillating processes to an external forcing signal, which may be periodic on a different time-scale from the actuated process. Such subharmonic synchrony, which is dynamically established when N control cycles occur for every M cycles of a forced oscillator, is referred to as N:M entrainment. In many applications, entrainment must be established in an optimal manner, for example by minimizing control energy or the transient time to phase locking. We present a theory for deriving inputs that establish subharmonic N:M entrainment of general nonlinear oscillators, or of collections of rhythmic dynamical units, while optimizing such objectives. Ordinary differential equation models of oscillating systems are reduced to phase variable representations, each of which consists of a natural frequency and phase response curve. Formal averaging and the calculus of variations are then applied to such reduced models in order to derive optimal subharmonic entrainment waveforms. The optimal entrainment of a canonical model for a spiking neuron is used to illustrate this approach, which is readily extended to arbitrary oscillating systems.