Refine
Year of publication
Keywords
We consider optimal control problems for the flow of gas or fresh water in pipe networks as well as drainage or sewer systems in open canals. The equations of motion are taken to be represented by the nonlinear isothermal Euler gas equations, the water hammer equations, or the St.~Venant equations for flow. We formulate model hierarchies and derive an abstract model for such network flow problems including pipes, junctions, and controllable elements such as valves, weirs, pumps, as well as compressors. We use the abstract model to give an overview of the known results and challenges concerning equilibria, well-posedness, controllability, and optimal control. A major challenge concerning the optimization is to deal with switching on-off states that are inherent to controllable devices in such applications combined with
continuous simulation and optimization of the gas flow. We formulate the corresponding mixed-integer nonlinear optimal control problems and outline a decomposition approach as a solution technique.
Nonconvex mixed-binary nonlinear optimization problems frequently appear in practice and are typically extremely hard to solve. In this paper we discuss a class of primal heuristics that are based on a reformulation of the problem as a mathematical program with equilibrium constraints. We then use different regularization schemes for this class of problems and use an iterative solution procedure for solving series of regularized problems. In the case of success, these procedures result in a feasible solution of the original mixed-binary nonlinear problem. Since we rely on local nonlinear programming solvers the resulting method is fast and we further improve its reliability by additional algorithmic techniques. We show the strength of our method by an extensive computational study on 662 MINLPLib2 instances, where our methods are able to produce feasible solutions for 60% of all instances in at most 10s.
Detailed modeling of gas transport problems leads to nonlinear
and nonconvex mixed-integer optimization or feasibility models
(MINLPs) because both the incorporation of discrete controls of the
network as well as accurate physical and technical modeling is
required in order to achieve practical solutions. Hence, ignoring
certain parts of the physics model is not valid for practice. In the
present contribution we extend an approach based on linear relaxations
of the underlying nonlinearities by tailored model reformulation
techniques yielding block-separable MINLPs. This combination of
techniques allows us to apply a penalty alternating direction method
and thus to solve highly detailed MINLPs for large-scale real-world
instances. The practical strength of the proposed method is
demonstrated by a computational study in which we apply the method to
instances from steady-state gas transport
including both pooling effects with respect to the mixing of gases of
different composition and a highly detailed compressor station model.
We propose an equilibrium model that allows to analyze the long-run impact of the electricity market design on transmission line expansion by the regulator and investment in generation capacity by private firms in liberalized electricity markets. The model incorporates investment decisions of the transmission system operator and private firms in expectation of an energy-only market and cost-based redispatch. In different specifications we consider the cases of one vs. multiple price zones (market splitting) and analyze different approaches to recover network cost—in particular lump sum, generation capacity based, and energy based fees. In order to compare the outcomes of our multilevel market model with a first best benchmark, we also solve the corresponding integrated planner problem. Using two test networks we illustrate that energy-only markets can lead to suboptimal locational decisions for generation capacity and thus imply excessive network expansion. Market splitting heals these problems only partially. These results are valid for all considered types of network tariffs, although investment slightly differs across those regimes.
We present a solution algorithm for problems from
steady-state gas transport optimization.
Due to nonlinear and nonconvex physics and engineering models as
well as discrete controllability of active network devices, these
problems lead to difficult nonconvex mixed-integer nonlinear optimization
models.
The proposed method is based on mixed-integer linear techniques using
piecewise linear relaxations of the nonlinearities and a tailored
alternating direction method.
Most other publications in the field of gas transport optimization only consider
pressure and flow as main physical quantities. In this work, we additionally
incorporate heat power supplies and demands as well as a mixing model for
different gas qualities.
We demonstrate the capabilities of our method on Germany's largest
transport networks and hereby present numerical results on the largest
instances that were ever reported in the literature for this problem
class.
Potential-based flows are an extension of classical network flows in which the flow on an arc is determined by the difference of the potentials of its incident nodes. Such flows are unique and arise, for example, in energy networks. Two important algorithmic problems are to determine whether there exists a feasible flow and to maximize the flow between two designated nodes. We show that these problems can be solved for the single source and sink case by reducing the network to a single arc. However, if we additionally consider switches that allow to force the flow to 0 and decouple the potentials, these problems are NP-hard. Nevertheless, for particular series-parallel networks, one can use algorithms for the subset sum problem. Moreover, applying network presolving based on generalized series-parallel structures allows to significantly reduce the size of realistic energy networks.
In the
Steiner Forest
problem, we are given a graph and a collection of source-sink
pairs, and the
goal is to find a subgraph of minimum total length such that all
pairs are connected. The problem is
APX-Hard and can be
2
-approximated by, e.g., the elegant primal-dual algorithm
of Agrawal, Klein, and
Ravi from 1995.
We give a local-search-based constant-factor approximati
on for the problem. Local search brings in
new techniques to an area that has for long not seen any improv
ements and might be a step towards
a combinatorial algorithm for the more general survivable n
etwork design problem. Moreover, local
search was an essential tool to tackle the dynamic MST/Stein
er Tree problem, whereas dynamic Steiner
Forest is still wide open.
It is easy to see that any constant factor local search algori
thm requires steps that add/drop many edges
together. We propose natural local moves which, at each step
, either (a) add a shortest path in the current
graph and then drop a bunch of inessential edges, or (b) add a s
et of edges to the current solution. This
second type of moves is motivated by the potential function w
e use to measure progress, combining the
cost of the solution with a penalty for each connected compon
ent. Our carefully-chosen local moves and
potential function work in tandem to eliminate bad local min
ima that arise when using more traditional
local moves.
Our analysis first considers the case where the local optimum
is a single tree, and shows optimality w.r.t.
moves that add a single edge (and drop a set of edges) is enough
to bound the locality gap. For the
general case, we show how to “project” the optimal solution o
nto the different trees of the local optimum
without incurring too much cost (and this argument uses opti
mality w.r.t. both kinds of moves), followed
by a tree-by-tree argument. We hope both the potential funct
ion, and our analysis techniques will be
useful to develop and analyze local-search algorithms in ot
her contexts.
Many mixed-integer optimization problems are constrained by nonlinear functions that do not possess desirable analytical properties like convexity or factorability or cannot even be evaluated exactly. This is, e.g., the case for problems constrained by differential equations or for models that rely on black-box simulation runs. For these problem classes, we present, analyze, and test algorithms that solve mixed-integer problems with only Lipschitz continuous nonlinearities. Our theoretical results depend on the assumptions made on the (in)exactness of function evaluations and on the knowledge of Lipschitz constants. If Lipschitz constants are known, we prove finite termination at approximate globally optimal points both for the case of exact and inexact function evaluations. If only approximate Lipschitz constants are known, we prove finite termination and derive additional conditions under which infeasibility can be detected. A computational study for gas transport problems and an academic case study show the applicability of our algorithms to real-world problems and how different assumptions on the constraint functions up- or downgrade the practical performance of the methods.
The development of mathematical simulation and optimization models and algorithms for solving gas transport problems is an active field of research. In order to test and compare these models and algorithms, gas network instances together with demand data are needed. The goal of GasLib is to provide a set of publicly available gas network instances that can be used by researchers in the field of gas transport. The advantages are that researchers save time by using these instances and that different models and algorithms can be compared on the same specified test sets. The library instances are encoded in an XML format. In this paper, we explain this format and present the instances that are available in the library.