TY - JOUR A1 - Grimm, Veronika A1 - Grübel, Julia A1 - Rückel, Bastian A1 - Sölch, Christian A1 - Zöttl, Gregor T1 - Storage investment and network expansion in distribution networks: The impact of regulatory frameworks JF - Applied Energy N2 - In this paper we propose a bi-level equilibrium model that allows to analyze the impact of different regulatory frameworks on storage and network investment in distribution networks. In our model, a regulated distribution system operator decides on network investment and operation while he anticipates the decisions of private agents on storage investment and operation. Since, especially in distribution networks, voltage stability and network losses have a decisive influence on network expansion and operation, we use a linearized AC power flow formulation to adequately account for these aspects. As adjustments of the current regulatory framework, we consider curtailment of renewable production, the introduction of a network fee based on the maximum renewable feed-in, and a subsidy scheme for storage investment. The performance of the different alternative frameworks is compared to the performance under rules that are commonly applied in various countries today, as well as to a system-optimal (first-best) benchmark. To illustrate the economic effects, we calibrate our model with data from the field project Smart Grid Solar. Our results reveal that curtailment and a redesign of network fees both have the potential to significantly reduce total system costs. On the contrary, investment subsidization of storage capacity has only a limited impact as long as the distribution system operator is not allowed to intervene in storage operation. KW - Computational Equilibrium Models; Electricity Markets; Investment Incentives; Distribution Network Expansion Planning; Storage Investment and Operation; Renewable Energy Production; (Self-) Consumption Y1 - 2019 U6 - https://doi.org/10.1016/j.apenergy.2019.114017 IS - 262 ER - TY - JOUR A1 - Runge, Philipp A1 - Sölch, Christian A1 - Albert, Jakob A1 - Wasserscheid, Peter A1 - Zöttl, Gregor A1 - Grimm, Veronika T1 - Economic comparison of different electric fuels for energy scenarios in 2035 JF - Applied Energy N2 - Electric fuels (e-fuels) enable CO2-neutral mobility and are therefore an alternative to battery-powered electric vehicles. This paper compares the cost-effectiveness of Fischer-Tropsch diesel, methanol and Liquid Organic Hydrogen Carriers. The production costs of those fuels are to a large part driven by the energy-intensive electrolytic hydrogen production. In this paper, we apply a multi-level electricity market model to calculate future hourly electricity prices for various electricity market designs in Germany for the year 2035. We then assess the economic efficiency of the different fuels under various future market conditions. In particular, we use the electricity price vectors derived from an electricity market model calibrated for 2035 as an input for a mathematical model of the entire process chain from hydrogen production and chemical bonding to the energetic utilization of the fuels in a vehicle. Within this model, we perform a sensitivity analysis, which quantifies the impact of various parameters on the fuel production cost. Most importantly, we consider prices resulting from own model calculations for different energy market designs, the investment cost for the electrolysis systems and the carbon dioxide purchase price. The results suggest that the use of hydrogen, which is temporarily bound to Liquid Organic Hydrogen Carriers, is a favorable alternative to the more widely discussed synthetic diesel and methanol. KW - Sector Coupling, Eletricity Markets, Electric fuels, Hydrogen Utilization, LOHC, Mobility Y1 - 2019 U6 - https://doi.org/10.1016/j.apenergy.2018.10.023 IS - 233-234 SP - 1078 EP - 1093 ER - TY - INPR A1 - Grimm, Veronika A1 - Martin, Alexander A1 - Sölch, Christian A1 - Weibelzahl, Martin A1 - Zöttl, Gregor T1 - Market-Based Redispatch May Result in Inefficient Dispatch N2 - In this paper we analyze a uniform price electricity spot market that is followed by redispatch in the case of network congestion. We assume that the transmission system operator is incentivized to minimize redispatch cost and compare a cost-based redispatch (CBR) to a market-based redispatch (MBR) mechanism. For networks with at least three nodes we show that in contrast to CBR, in the case of MBR the redispatch cost minimizing allocation may not be short-run efficient. As we demonstrate, in case of MBR the possibility of the transmission system operator to reduce redispatch cost at the expense of a reduced welfare may be driven by the electricity supply side or the electricity demand side. If, however, the transmission system operator is obliged to implement the welfare maximizing (instead of the redispatch cost minimizing) dispatch by regulation, this will result in an efficient dispatch also in case of MBR. KW - Electricity Markets, Redispatch, Congestion Management, Computational Equilibrium Models Y1 - 2019 ER - TY - JOUR A1 - Ambrosius, Mirjam A1 - Egerer, Jonas A1 - Grimm, Veronika A1 - van der Weijde, Adriaan H. T1 - Uncertain bidding zone configurations: The role of expectations for transmission and generation capacity expansion JF - European Journal of Operational Research N2 - Ongoing policy discussions on the reconfiguration of bidding zones in European electricity markets induce uncertainty about the future market design. This paper deals with the question of how this uncertainty affects market participants and their long-run investment decisions in generation and transmission capacity. Generalizing the literature on pro-active network expansion planning, we propose a stochastic multilevel model which incorporates generation capacity investment, network expansion, and market operation, taking into account uncertainty about the future bidding zone configuration. Using a stylized two-node network, we disentangle different effects that uncertainty has on market outcomes. If there is a possibility that future bidding zone configurations provide improved regional price signals, welfare gains materialize even if the change does not actually take place. As a consequence, welfare gains of an actual change of the bidding zone configuration are substantially lower due to those anticipatory effects. Additionally, we show substantial distributional effects in terms of both expected gains and risks, between producers and consumers and between different generation technologies. Y1 - 2019 U6 - https://doi.org/10.1016/j.ejor.2020.01.024 IS - 285 (1) SP - 343 EP - 359 ER - TY - JOUR A1 - Bichler, Martin A1 - Grimm, Veronika A1 - Kretschmer, Sandra A1 - Sutterer, Paul T1 - Market Design for Renewable Energy Auctions: An Analysis of Alternative Auction Formats JF - Energy Economics N2 - Auctions are widely used to determine the remuneration for renewable energies. They typically induce a high concentration of renewable energy plants at very productive sites far-off the main load centres, leading to an inefficient allocation as transmission line capacities are restricted but not considered in the allocation, resulting in an inefficient system configuration in the long run. To counteract these tendencies effectively, we propose a combinatorial auction design that allows to implement regional target capacities, provides a simple pricing rule and maintains a high level of competition between bidders by permitting package bids. By means of extensive numerical experiments we evaluate the combinatorial auction as compared to three further RES auction designs, the current German nationwide auction design, a simple nationwide auction, and regional auctions. We find that if bidders benefit from high enough economies of scale, the combinatorial auction design implements system-optimal target capacities without increasing the average remuneration per kWh as compared to the current German auction design. The prices resulting from the combinatorial auction are linear and anonymous for each region whenever possible, while minimal personalised markups on the linear prices are applied only when necessary. We show that realistic problem sizes can be solved in seconds, even though the problem is computationally hard. KW - renewable energy auction, combinatorial auction, market design Y1 - 2019 U6 - https://doi.org/10.1016/j.eneco.2020.104904 IS - 92 ER - TY - JOUR A1 - Grimm, Veronika A1 - Kretschmer, Sandra A1 - Mehl, Simon T1 - Green innovations: The organizational setup of pilot projects and its influence on consumer perceptions JF - Energy Policy N2 - Pilot-, test- and demonstration-projects (PTDs) are a prominent policy tool to promote the adoption of smart, green technologies. However, as technology adoption is heavily dependent on the individual attributes and beliefs of potential adopters, it is important to understand the influence of a PTD’s organizational setup on technology perception. By varying the information about a PTD’s organizational setup in a survey experiment among a selected sample of potential PTD-participants, we gather first experimental evidence for the effect of different setups on the perception of green technologies. We show that the organizational setup has a significant impact on a product’s perceived contribution to the energy transition, its establishment in the market, cost-reduction potential, innovativeness and environmental friendliness. In particular, full organizational cooperation between government, university and industry consistently improves perceptions compared to a partial setup. Regarding the willingness to participate in a PTD, we find that communication and support are the most imperative aspects and even more important than economic benefits. Our findings provide policy-makers with a more ample foundation on how PTDs should be designed to successfully transfer technologies to the market. KW - Pilot projects, Organizational setup, Adoption decision, Quadruple Helix, Innovation diffusion Y1 - 2019 U6 - https://doi.org/10.1016/j.enpol.2020.111474 IS - 142 ER - TY - INPR A1 - Veronika, Grimm A1 - Kretschmer, Sandra A1 - Mehl, Simon T1 - Home energy storage: A characterization of potential adopters in Germany N2 - Home energy storage (HES) poses a promising option to foster flexible electricity demand and grid-stabilising self-supply. Understanding the characteristics of potential storage adopters and their motives for adopting is essential to enable targeted measures to promote diffusion. This paper presents an in-depth analysis of motivational and psychological factors as well as product characteristics that affect the willingness to adopt domestic energy storage. Based on a selective sample of solar panel owners as the main energy storage target group, we find that potential adopters differ systematically in their evaluation of economic and non-economic characteristics of home energy storage. Using correlation-based average linkage clustering, we identify four distinct types of storage adopters: finance-oriented, independence-oriented, idealistic and multilaterally oriented households. Our findings emphasise that segment-specific policy strategies and support mechanisms are needed to stimulate HES adoption to pave the way for a sustainable energy system. KW - Energy storage, Consumer segmentation, Adoption decision, Innovation diffusion Y1 - 2019 ER - TY - JOUR A1 - Kleinert, Thomas A1 - Grimm, Veronika A1 - Schmidt, Martin T1 - Outer Approximation for Global Optimization of Mixed-Integer Quadratic Bilevel Problems JF - Mathematical Programming (Series B) N2 - Bilevel optimization problems have received a lot of attention in the last years and decades. Besides numerous theoretical developments there also evolved novel solution algorithms for mixed-integer linear bilevel problems and the most recent algorithms use branch-and-cut techniques from mixed-integer programming that are especially tailored for the bilevel context. In this paper, we consider MIQP-QP bilevel problems, i.e., models with a mixed-integer convex-quadratic upper level and a continuous convex-quadratic lower level. This setting allows for a strong-duality-based transformation of the lower level which yields, in general, an equivalent nonconvex single-level reformulation of the original bilevel problem. Under reasonable assumptions, we can derive both a multi- and a single-tree outer-approximation-based cutting-plane algorithm. We show finite termination and correctness of both methods and present extensive numerical results that illustrate the applicability of the approaches. It turns out that the proposed methods are capable of solving bilevel instances with several thousand variables and constraints and significantly outperform classical solution approaches. KW - Bilevel optimization KW - Outer approximation KW - Quadratic programming KW - Convex mixed-integer nonlinear optimization Y1 - 2019 ER - TY - JOUR A1 - Grimm, Veronika A1 - Kleinert, Thomas A1 - Liers, Frauke A1 - Schmidt, Martin A1 - Zöttl, Gregor T1 - Optimal Price Zones of Electricity Markets: A Mixed-Integer Multilevel Model and Global Solution Approaches JF - Optimization Methods and Software N2 - Mathematical modeling of market design issues in liberalized electricity markets often leads to mixed-integer nonlinear multilevel optimization problems for which no general-purpose solvers exist and which are intractable in general. In this work, we consider the problem of splitting a market area into a given number of price zones such that the resulting market design yields welfare-optimal outcomes. This problem leads to a challenging multilevel model that contains a graph-partitioning problem with multi-commodity flow connectivity constraints and nonlinearities due to proper economic modeling. Furthermore, it has highly symmetric solutions. We develop different problem-tailored solution approaches. In particular, we present an extended KKT transformation approach as well as a generalized Benders approach that both yield globally optimal solutions. These methods, enhanced with techniques such as symmetry breaking and primal heuristics, are evaluated in detail on academic as well as on realistic instances. It turns out that our approaches lead to effective solution methods for the difficult optimization tasks presented here, where the problem-specific generalized Benders approach performs considerably better than the methods based on KKT transformation. KW - Multilevel Optimization KW - Mixed-Integer Nonlinear Optimization KW - Graph Partitioning KW - Generalized Benders Decomposition KW - Electricity Market Design} Y1 - 2017 IS - 34(2) SP - 406 EP - 436 ER - TY - INPR A1 - Hiller, Benjamin A1 - Koch, Thorsten A1 - Schewe, Lars A1 - Schwarz, Robert A1 - Schweiger, Jonas T1 - A System to Evaluate Gas Network Capacities: Concepts and Implementation T2 - EJOR N2 - Since 2005, the gas market in the European Union is liberalized and the trading of natural gas is decoupled from its transport. The transport is done by so-called transmissions system operators (TSOs). The market model established by the European Union views the gas transmission network as a black box, providing shippers (gas traders and consumers) the opportunity to transport gas from any entry to any exit. TSOs are required to offer maximum independent capacities at each entry and exit such that the resulting gas flows can be realized by the network without compromising security of supply. Therefore, evaluating the available transport capacities is extremely important to the TSOs. This paper gives an overview of the toolset for evaluating gas network capacities that has been developed within the ForNe project, a joint research project of seven research partners initiated by Open Grid Europe, Germany's biggest TSO. While most of the relevant mathematics is described in the book "Evaluating Gas Network Capacities", this article sketches the system as a whole, describes some developments that have taken place recently, and gives some details about the current implementation. Y1 - 2017 ER - TY - JOUR A1 - Grimm, Veronika A1 - Schewe, Lars A1 - Schmidt, Martin A1 - Zöttl, Gregor T1 - A Multilevel Model of the European Entry-Exit Gas Market JF - Mathematical Methods of Operations Research N2 - In entry-exit gas markets as they are currently implemented in Europe, network constraints do not affect market interaction beyond the technical capacities determined by the TSO that restrict the quantities individual firms can trade at the market. It is an up to now unanswered question to what extent existing network capacity remains unused in an entry-exit design and to what extent feasible adjustments of the market design could alleviate inefficiencies. In this paper, we offer a four-level modeling framework that is capable of analyzing these issues and provide some first results on the model structure. In order to decouple gas trading from network congestion management, the TSO is required to determine technical capacities and corresponding booking fees at every entry and exit node up front. Firms book those capacities, which gives them the right to charge or discharge an amount of gas at a certain node up to this capacity in every scenario. Beyond these technical capacities and the resulting bookings, gas trade is unaffected by network constraints. The technical capacities have to ensure that transportation of traded quantities is always feasible. We assume that the TSO is regulated and determines technical capacities, fees, and transportation costs under a welfare objective. As a first step we moreover assume perfect competition among gas traders and show that the booking and nomination decisions can be analyzed in a single level. We prove that this aggregated model has a unique solution. We also show that the TSO's decisions can be subsumed in one level as well. If so, the model boils down to a mixed-integer nonlinear bilevel problem with robust aspects. In addition, we provide a first-best benchmark that allows to assess welfare losses that occur in an entry-exit system. Our approach provides a generic framework to analyze various aspects in the context of semi-liberalized gas markets. Therefore, we finally discuss and provide guidance on how to include several important aspects into the approach, such as network and production capacity investment, uncertain data, market power, and intra-day trading. KW - Entry-Exit System KW - Gas Market KW - Multilevel Modeling Y1 - 2017 IS - 89(2) SP - 223 EP - 255 ER - TY - INPR A1 - Aßmann, Denis A1 - Liers, Frauke A1 - Stingl, Michael A1 - Vera, Juan T1 - Deciding Robust Feasibility and Infeasibility Using a Set Containment Approach: An Application to Stationary Passive Gas Network Operations N2 - In this paper we study feasibility and infeasibility of nonlinear two-stage fully adjustable robust feasibility problems with an empty first stage. This is equivalent to deciding set containment of a projection of the feasible region and the uncertainty set. For answering this question, two very general approaches using methods from polynomial optimization are presented --- one for showing feasibility and one for showing infeasibility. The developed methods are approximated through sum of squares polynomials and solved using semidefinite programs. Deciding robust feasibility and infeasibility is important for gas network operations, which is a \nonconvex quadratic problem with absolute values functions. Concerning the gas network problem, different topologies are considered. It is shown that a tree structured network can be decided exactly using linear programming. Furthermore, a method is presented to reduce a tree network with one additional arc to a single cycle network. In this case, removing the absolute values and solving the problem can be decided with linearly many polynomial optimization problems. Lastly, the effectivity of the methods is tested on a variety of small cyclic networks. For instances where robust feasibility or infeasibility can be decided, level~2 or level~3 of the Lasserre relaxation hierarchy is typically sufficient. KW - robust optimization KW - polynomial optimization KW - stationary gas transport Y1 - 2017 ER - TY - THES A1 - Wollenberg, Tobias T1 - Two-Stage Stochastic Semidefinite Programming: Theory, Algorithms, and Application to AC Power Flow under Uncertainty N2 - In real life decision problems, one almost always is confronted with uncertainty and risk. For practical optimization problems this is manifested by unknown parameters within the input data, or, an inexact knowledge about the system description itself. In case the uncertain problem data is governed by a known probability distribution, stochastic programming offers a variety of models hedging against uncertainty and risk. Most widely employed are two-stage models, who admit a recourse structure: The first-stage decisions are taken before the random event occurs. After its outcome, a recourse (second-stage) action is made, often but not always understood as some "compensation". In the present thesis, the optimization problems that involve parameters which are not known with certainty are semidefinite programming problems. The constraint sets of these optimization problems are given by intersections of the cone of symmetric, positive semidefinite matrices with either affine or more general equations. Objective functions, formally, may be fairly general, although they often are linear as in the present thesis. We consider risk neutral and risk averse two-stage stochastic semidefinite programs with continuous and mixed-integer recourse, respectively. For these stochastic optimization problems we analyze their structure, derive solution methods relying on decomposition, and finally apply our results to unit commitment in alternating current (AC) power systems. Furthermore, deterministic unit commitment in AC power transmission systems is addressed. Beside traditional unit commitment constraints, the physics of power flow are included. To gain globally optimal solutions a recent semidefinite programming (SDP) approach is used which leads to large-scale semidefinite programs with discrete variables on top. As even the SDP relaxation of these programs is too large for being handled in an all-at-once manner by general SDP solvers, it requires an efficient and reliable method to tackle them. To this end, an algorithm based on Benders decomposition is proposed. With power demand (load) and in-feed from renewables serving as sources of uncertainty, two-stage stochastic programs are set up heading for unit commitment schedules which are both cost-effective and robust with respect to data perturbations. The impact of different, risk neutral and risk averse, stochastic criteria on the shapes of the optimal stochastic solutions will be examined. To tackle the resulting two-stage programs, we propose to approximate AC power flow by semidefinite relaxations. This leads to two-stage stochastic mixed-integer semidefinite programs having a special structure. To solve the latter, the L-shaped method and dual decomposition have been applied and compared. KW - Stochastic Programming KW - Semidefinite Programming KW - Decomposition Methods KW - Risk Aversion KW - Unit Commitment with AC Power Flow Y1 - 2017 ER - TY - JOUR A1 - Gugat, Martin A1 - Leugering, Günter A1 - Martin, Alexander A1 - Schmidt, Martin A1 - Sirvent, Mathias A1 - Wintergerst, David T1 - MIP-Based Instantaneous Control of Mixed-Integer PDE-Constrained Gas Transport Problems JF - Computational Optimization and Applications N2 - We study the transient optimization of gas transport networks including both discrete controls due to switching of controllable elements and nonlinear fluid dynamics described by the system of isothermal Euler equations, which are partial differential equations in time and 1-dimensional space. This combination leads to mixed-integer optimization problems subject to nonlinear hyperbolic partial differential equations on a graph. We propose an instantaneous control approach in which suitable Euler discretizations yield systems of ordinary differential equations on a graph. This networked system of ordinary differential equations is shown to be well-posed and affine-linear solutions of these systems are derived analytically. As a consequence, finite-dimensional mixed-integer linear optimization problems are obtained for every time step that can be solved to global optimality using general-purpose solvers. We illustrate our approach in practice by presenting numerical results on a realistic gas transport network. KW - Mixed-integer optimal control KW - Instantaneous control KW - Partial differential equations on graphs KW - Gas networks KW - Mixed-integer linear optimization Y1 - 2017 U6 - https://doi.org/10.1007/s10589-017-9970-1 VL - 70 IS - 1 SP - 267 EP - 294 ER - TY - JOUR A1 - Kramer, Anja A1 - Krebs, Vanessa A1 - Schmidt, Martin T1 - Strictly and Γ-Robust Counterparts of Electricity Market Models: Perfect Competition and Nash-Cournot Equilibria JF - Operations Research Perspectives N2 - This paper mainly studies two topics: linear complementarity problems for modeling electricity market equilibria and optimization under uncertainty. We consider both perfectly competitive and Nash–Cournot models of electricity markets and study their robustifications using strict robustness and the Γ-approach. For three out of the four combinations of economic competition and robustification, we derive algorithmically tractable convex optimization counterparts that have a clear-cut economic interpretation. In the case of perfect competition, this result corresponds to the two classical welfare theorems, which also apply in both considered robust cases that again yield convex robustified problems. Using the mentioned counterparts, we can also prove the existence and, in some cases, uniqueness of robust equilibria. Surprisingly, it turns out that there is no such economic sensible counterpart for the case of Γ-robustifications of Nash–Cournot models. Thus, an analogue of the welfare theorems does not hold in this case. Finally, we provide a computational case study that illustrates the different effects of the combination of economic competition and uncertainty modeling. KW - Robust optimization KW - Linear complementarity problems KW - Electricity market equilibrium models KW - Perfect competition KW - Nash-Cournot competition Y1 - 2018 IS - 89(2) SP - 100197 ER - TY - JOUR A1 - Hojny, Christopher A1 - Joormann, Imke A1 - Lüthen, Hendrik A1 - Schmidt, Martin T1 - Mixed-Integer Programming Techniques for the Connected Max-k-Cut Problem JF - Mathematical Programming Computation N2 - We consider an extended version of the classical Max-k-Cut problem in which we additionally require that the parts of the graph partition are connected. For this problem we study two alternative mixed-integer linear formulations and review existing as well as develop new branch-and-cut techniques like cuts, branching rules, propagation, primal heuristics, and symmetry breaking. The main focus of this paper is an extensive numerical study in which we analyze the impact of the different techniques for various test sets. It turns out that the techniques from the existing literature are not sufficient to solve an adequate fraction of the test sets. However, our novel techniques significantly outperform the existing ones both in terms of running times and the overall number of instances that can be solved. KW - Max-cut KW - Connectivity KW - Branch-and-cut KW - Mixed-integer programming Y1 - 2018 ER - TY - JOUR A1 - Hante, Falk T1 - Mixed-Integer Optimal Control for PDEs: Relaxation via Differential Inclusions and Applications to Gas Network Optimization JF - Mathematical Modelling, Optimization, Analytic and Numerical Solutions, Springer Series on Industrial and Applied Mathematics, 2020 N2 - We show that mixed-integer control problems for evolution type partial differential equations can be regarded as operator differential inclusions. This yields a relaxation result including a characterization of the optimal value for mixed-integer optimal control problems with control constraints. The theory is related to partial outer convexification and sum-up rounding methods. The results are applied to optimal valve switching control for gas pipeline operations. A numerical example illustrates the approach. Y1 - 2018 ER - TY - JOUR A1 - Schewe, Lars A1 - Schmidt, Martin A1 - Thürauf, Johannes T1 - Structural Properties of Feasible Bookings in the European Entry-Exit Gas Market System JF - 4OR N2 - In this work we analyze the structural properties of the set of feasible bookings in the European entry-exit gas market system. We present formal definitions of feasible bookings and then analyze properties that are important if one wants to optimize over them. Thus, we study whether the sets of feasible nominations and bookings are bounded, convex, connected, conic, and star-shaped. The results depend on the specific model of gas flow in a network. Here, we discuss a simple linear flow model with arc capacities as well as nonlinear and mixed-integer nonlinear models of passive and active networks, respectively. It turns out that the set of feasible bookings has some unintuitive properties. For instance, we show that the set is nonconvex even though only a simple linear flow model is used. KW - Gas networks KW - Booking KW - Entry-exit system KW - Convexity KW - Flow models Y1 - 2018 U6 - https://doi.org/10.1007/s10288-019-00411-3 IS - 18 SP - 197 EP - 218 ER - TY - JOUR A1 - Ambrosius, Mirjam A1 - Grimm, Veronika A1 - Kleinert, Thomas A1 - Liers, Frauke A1 - Schmidt, Martin A1 - Zöttl, Gregor T1 - Endogenous Price Zones and Investment Incentives in Electricity Markets: An Application of Multilevel Optimization with Graph Partitioning JF - Energy Economics N2 - In the course of the energy transition, load and supply centers are growing apart in electricity markets worldwide, rendering regional price signals even more important to provide adequate locational investment incentives. This paper focuses on electricity markets that operate under a zonal pricing market design. For a fixed number of zones, we endogenously derive the optimal configuration of price zones and available transfer capacities on a network in order to optimally govern investment and production decisions in the long run. In a multilevel mixed-integer nonlinear model that contains a graph partitioning problem on the first level, we determine welfare-maximizing price zones and available transfer capacities for a given electricity market and analyze their impact on market outcomes. Using a generalized Benders decomposition approach developed in Grimm et al. (2019) and a problem-tailored scenario clustering for reducing the input data size, we are able to solve the model to global optimality even for large instances. We apply the approach to the German electricity market as an example to examine the impact of optimal zoning on key performance indicators such as welfare, generation mix and locations, or electricity prices. It turns out that even for a small number of price zones, an optimal configuration of zones induces a welfare level that almost approaches the first best. KW - Electricity Markets KW - Price Zones KW - Investment Incentives KW - Multilevel Optimization KW - Graph Partitioning Y1 - 2018 IS - 92 ER - TY - JOUR A1 - Hante, Falk M. A1 - Schmidt, Martin T1 - Convergence of Finite-Dimensional Approximations for Mixed-Integer Optimization with Differential Equations JF - Control and Cybernetics N2 - We consider a direct approach to solve mixed-integer nonlinear optimization problems with constraints depending on initial and terminal conditions of an ordinary differential equation. In order to obtain a finite-dimensional problem, the dynamics are approximated using discretization methods. In the framework of general one-step methods, we provide sufficient conditions for the convergence of this approach in the sense of the corresponding optimal values. The results are obtained by considering the discretized problem as a parametric mixed-integer nonlinear optimization problem in finite dimensions, where the maximum step size for discretizing the dynamics is the parameter. In this setting, we prove the continuity of the optimal value function under a stability assumption for the integer feasible set and second-order conditions from nonlinear optimization. We address the necessity of the conditions on the example of pipe sizing problems for gas networks. KW - Optimization with differential equations KW - Optimal value function KW - Lipschitz continuity KW - Parametric optimization KW - Mixed-integer nonlinear programming Y1 - 2018 ER - TY - JOUR A1 - Labbé, Martine A1 - Plein, Fränk A1 - Schmidt, Martin T1 - Bookings in the European Gas Market: Characterisation of Feasibility and Computational Complexity Results JF - Optimization and Engineering N2 - As a consequence of the liberalisation of the European gas market in the last decades, gas trading and transport have been decoupled. At the core of this decoupling are so-called bookings and nominations. Bookings are special capacity right contracts that guarantee that a specified amount of gas can be supplied or withdrawn at certain entry or exit nodes of the network. These supplies and withdrawals are nominated at the day-ahead. The special property of bookings then is that they need to be feasible, i.e., every nomination that complies with the given bookings can be transported. While checking the feasibility of a nomination can typically be done by solving a mixed-integer nonlinear feasibility problem, the verification of feasibility of a set of bookings is much harder. The reason is the robust nature of feasibility of bookings - namely that for a set of bookings to be feasible, all compliant nominations, i.e., infinitely many, need to be checked for feasibility. In this paper, we consider the question of how to verify the feasibility of given bookings for a number of special cases. For our physics model we impose a steady-state potential-based flow model and disregard controllable network elements. For this case we derive a characterisation of feasible bookings, which is then used to show that the problem is in coNP for the general case but can be solved in polynomial time for linear potential-based flow models. Moreover, we present a dynamic programming approach for deciding the feasibility of a booking in tree-shaped networks even for nonlinear flow models. It turns out that the hardness of the problem mainly depends on the combination of the chosen physics model as well as the specific network structure under consideration. Thus, we give an overview over all settings for which the hardness of the problem is known and finally present a list of open problems. KW - Gas networks KW - Booking KW - Nomination KW - Computational complexity KW - Trees Y1 - 2018 IS - 21(1) SP - 305 EP - 334 ER - TY - JOUR A1 - Gugat, Martin A1 - Hante, Falk T1 - On the turnpike phenomenon for optimal boundary control problems with hyperbolic systems JF - SIAM Journal on Control and Optimization N2 - We study problems of optimal boundary control with systems governed by linear hyperbolic partial differential equations. The objective function is quadratic and given by an integral over the finite time interval (0,T) that depends on the boundary traces of the solution. If the time horizon T is sufficiently large, the solution of the dynamic optimal boundary control problem can be approximated by the solution of a steady state optimization problem. We show that for T to infinity the approximation error converges to zero in the sense of the norm in L^2(0,1) with the rate 1/T, if the time interval (0,T) is transformed to the fixed interval (0,1). Moreover, we show that also for optimal boundary control problems with integer constraints for the controls the turnpike phenomenon occurs. In this case the steady state optimization problem also has the integer constraints. If T is sufficiently large, the integer part of each solution of the dynamic optimal boundary control problem with integer constraints is equal to the integer part of a solution of the static problem. A numerical verification is given for a control problem in gas pipeline operations. Y1 - 2019 U6 - https://doi.org/10.1137/17M1134470 VL - 57 IS - 1 SP - 264 EP - 289 ER - TY - JOUR A1 - Krebs, Vanessa A1 - Schmidt, Martin T1 - Γ-Robust Linear Complementarity Problems JF - Optimization Methods and Software N2 - Complementarity problems are often used to compute equilibria made up of specifically coordinated solutions of different optimization problems. Specific examples are game-theoretic settings like the bimatrix game or energy market models like for electricity or natural gas. While optimization under uncertainties is rather well-developed, the field of equilibrium models represented by complementarity problems under uncertainty - especially using the concepts of robust optimization - is still in its infancy. In this paper, we extend the theory of strictly robust linear complementarity problems (LCPs) to Γ-robust settings, where existence of worst-case-hedged equilibria cannot be guaranteed. Thus, we study the minimization of the worst-case gap function of Γ-robust counterparts of LCPs. For box and l1-norm uncertainty sets we derive tractable convex counterparts for monotone LCPs and study their feasibility as well as the existence and uniqueness of solutions. To this end, we consider uncertainties in the vector and in the matrix defining the LCP. We additionally study so-called ρ-robust solutions, i.e., solutions of relaxed uncertain LCPs. Finally, we illustrate the Γ-robust concept applied to LCPs in the light of the above mentioned classical examples of bimatrix games and market equilibrium modeling. KW - Linear complementarity problems KW - Robust optimization KW - Optimization under uncertainty KW - Γ-robustness KW - Tractable counterparts Y1 - 2019 ER - TY - JOUR A1 - Farshbaf Shaker, Mohammad Hassan A1 - Gugat, Martin A1 - Heitsch, Holger A1 - Henrion, René T1 - Optimal Neumann boundary control of a vibrating string with uncertain initial data and probabilistic terminal constraints N2 - In optimal control problems, often initial data are required that are not known exactly in practice. In order to take into account this uncertainty, we consider optimal control problems for a system with an uncertain initial state. A finite terminal time is given. On account of the uncertainty of the initial state, it is not possible to prescribe an exact terminal state. Instead, we are looking for controls that steer the system into a given neighborhood of the desired terminal state with sufficiently high probability. This neighborhood is described in terms of an inequality for the terminal energy. The probabilistic constraint in the considered optimal control problem leads to optimal controls that are robust against the inevitable uncertainties of the initial state. We show the existence of such optimal controls. Numerical examples with optimal Neumann control of the wave equation are presented. KW - PDE constrained optimization, probabilistic constraints, uncertain initial data Y1 - U6 - https://doi.org/10.1137/19M1269944 ER - TY - JOUR A1 - Krug, Richard A1 - Mehrmann, Volker A1 - Schmidt, Martin T1 - Nonlinear Optimization of District Heating Networks JF - Optimization and Engineering N2 - We develop a complementarity-constrained nonlinear optimization model for the time-dependent control of district heating networks. The main physical aspects of water and heat flow in these networks are governed by nonlinear and hyperbolic 1d partial differential equations. In addition, a pooling-type mixing model is required at the nodes of the network to treat the mixing of different water temperatures. This mixing model can be recast using suitable complementarity constraints. The resulting problem is a mathematical program with complementarity constraints subject to nonlinear partial differential equations describing the physics. In order to obtain a tractable problem, we apply suitable discretizations in space and time, resulting in a finite-dimensional optimization problem with complementarity constraints for which we develop a suitable reformulation with improved constraint regularity. Moreover, we propose an instantaneous control approach for the discretized problem, discuss practically relevant penalty formulations, and present preprocessing techniques that are used to simplify the mixing model at the nodes of the network. Finally, we use all these techniques to solve realistic instances. Our numerical results show the applicability of our techniques in practice. KW - District heating networks KW - Nonlinear optimization KW - Euler equations KW - Differential-algebraic equations KW - Complementarity constraints Y1 - 2019 IS - 22(2) SP - 783 EP - 819 ER - TY - INPR A1 - Krebs, Vanessa A1 - Müller, Michael A1 - Schmidt, Martin T1 - Γ-Robust Linear Complementarity Problems with Ellipsoidal Uncertainty Sets T2 - International Transactions in Operational Research N2 - We study uncertain linear complementarity problems (LCPs), i.e., problems in which the LCP vector q or the LCP matrix M may contain uncertain parameters. To this end, we use the concept of Γ-robust optimization applied to the gap function formulation of the LCP. Thus, this work builds upon [16]. There, we studied Γ-robustified LCPs for l1- and box-uncertainty sets, whereas we now focus on ellipsoidal uncertainty set. For uncertainty in q or M, we derive conditions for the tractability of the robust counterparts. For these counterparts, we also give conditions for the existence and uniqueness of their solutions. Finally, a case study for the uncertain traffic equilibrium problem is considered, which illustrates the effects of the values of Γ on the feasibility and quality of the respective robustified solutions. KW - Robust optimization KW - Linear complementarity problems KW - Ellipsoidal uncertainty sets KW - Traffic equilibrium problems Y1 - 2019 IS - 29(1) SP - 417 EP - 441 ER - TY - INPR A1 - Schmidt, Martin A1 - Hiller, Benjamin A1 - Koch, Thorsten A1 - Pfetsch, Marc A1 - Geißler, Björn A1 - Henrion, René A1 - Joormann, Imke A1 - Martin, Alexander A1 - Morsi, Antonio A1 - Römisch, Werner A1 - Schewe, Lars A1 - Schultz, Rüdiger A1 - Steinbach, Marc C. T1 - Capacity Evaluation for Large-Scale Gas Networks N2 - Natural gas is important for the energy turnaround in many countries like in Germany, where it serves as a "bridging energy" towards a fossil-free energy supply in the future. About 20% of the total German energy demand is provided by natural gas, which is transported through a complex pipeline network with a total length of about 30000 km and the efficient use of the given transport infrastructure for natural gas is of political, economic, and societal importance. As a consequence of the liberalization of the European gas market in the last decades, gas trading and transport have been decoupled. This has led to new challenges for gas transport companies, and mathematical optimization is perfectly suited for tackling many of these challenges. However, the underlying mathematical problems are by far too hard to be solved by today's general-purpose software so that novel mathematical theory and algorithms are needed. The industrial research project "ForNe: Research Cooperation Network Optimization" has been initiated and funded by Open Grid Europe in 2009 and brought together experts in mathematical optimization from seven German universities and research institutes, which cover almost the entire range of mathematical optimization: integer and nonlinear optimization as well as optimization under uncertainty. The mathematical research results have been put together in a software package that has been delivered to Open Grid Europe at the end of the project. Moreover, the research is still continuing - e.g., in the Collaborative Research Center/Transregio 154 "Mathematical Modelling, Simulation and Optimization using the Example of Gas Networks" funded by the German Research Foundation. Y1 - 2019 ER - TY - JOUR A1 - Labbé, Martine A1 - Plein, Fränk A1 - Schmidt, Martin A1 - Thürauf, Johannes T1 - Deciding Feasibility of a Booking in the European Gas Market on a Cycle is in P for the Case of Passive Networks JF - Networks N2 - We show that the feasibility of a booking in the European entry-exit gas market can be decided in polynomial time on single-cycle networks that are passive, i.e., do not contain controllable elements. The feasibility of a booking can be characterized by solving polynomially many nonlinear potential-based flow models for computing so-called potential-difference maximizing load flow scenarios. We thus analyze the structure of these models and exploit both the cyclic graph structure as well as specific properties of potential-based flows. This enables us to solve the decision variant of the nonlinear potential-difference maximization by reducing it to a system of polynomials of constant dimension that is independent of the cycle's size. This system of fixed dimension can be handled with tools from real algebraic geometry to derive a polynomial-time algorithm. The characterization in terms of potential-difference maximizing load flow scenarios then leads to a polynomial-time algorithm for deciding the feasibility of a booking. Our theoretical results extend the existing knowledge about the complexity of deciding the feasibility of bookings from trees to single-cycle networks. KW - Gas networks KW - European entry-exit market KW - Bookings KW - Potential-based flows KW - Computational complexity Y1 - 2019 U6 - https://doi.org/10.1007/s00186-021-00752-y VL - 78 IS - 2 SP - 128 EP - 152 ER - TY - JOUR A1 - Egerer, Jonas A1 - Grimm, Veronika A1 - Kleinert, Thomas A1 - Schmidt, Martin A1 - Zöttl, Gregor T1 - The Impact of Neighboring Markets on Renewable Locations, Transmission Expansion, and Generation Investment JF - European Journal of Operational Research N2 - Many long-term investment planning models for liberalized electricity markets either optimize for the entire electricity system or focus on confined jurisdictions, abstracting from adjacent markets. In this paper, we provide models for analyzing the impact of the interdependencies between a core electricity market and its neighboring markets on key long-run decisions. This we do both for zonal and nodal pricing schemes. The identification of welfare optimal investments in transmission lines and renewable capacity within a core electricity market requires a spatially restricted objective function, which also accounts for benefits from cross-border electricity trading. This leads to mixed-integer nonlinear multilevel optimization problems with bilinear nonconvexities for which we adapt a Benders-like decomposition approach from the literature. In a case study, we use a stylized six-node network to disentangle different effects of optimal regional (as compared to supra-regional) investment planning. Regional planning alters investment in transmission and renewable capacity in the core region, which affects private investment in generation capacity also in adjacent regions and increases welfare in the core region at the cost of system welfare. Depending on the congestion-pricing scheme, the regulator of the core region follows different strategies to increase welfare causing distributional effects among stakeholders. KW - Neighboring Markets KW - Renewables KW - Network Expansion KW - Multilevel Optimization KW - Benders Decomposition Y1 - 2019 ER - TY - INPR A1 - Bohlayer, Markus A1 - Fleschutz, Markus A1 - Braun, Marco A1 - Zöttl, Gregor T1 - Energy-intense production-inventory planning with participation in sequential energy markets N2 - To support the uprise of demand response, especially in the context of industrial processes, we propose a new approach to integrally determine the production-inventory plan and the cost-minimizing bids to participate in sequential reserve and energy-only markets. In particular, our approach considers time-coupling constraints which occur in the context of a production-inventory planning problem. We extend this problem with a comprehensive bidding formulation, which allows evaluating revenues and potential cost from the market participation, considering price uncertainties and uncertain activations of committed reserve capacity. This results in a multistage stochastic mixed-integer linear program, which explicitly considers the stage-wise revelation of information in our setup. To illustrate the capabilities of our approach, we apply our model to a real-world case study in which we investigate the participation of a cement plant in the German energy-only and reserve markets. The results of our case study indicate significant revenues for flexible industrial processes when participating in German spot and reserve markets. KW - Demand side flexibility Load management Multi market bidding Stochastic programming Production planning Demand response Y1 - 2019 ER - TY - JOUR A1 - Grimm, Veronika A1 - Rückel, Bastian A1 - Sölch, Christian A1 - Zöttl, Gregor T1 - The impact of market design on transmission and generation investment in electricity markets JF - Energy Economics N2 - In this paper we propose an equilibrium model in order to analyze the impact of electricity market design on generation and transmission expansion in liberalized electricity markets. In a multi-level structure, our framework takes into account that generation investment and operation is decided by private investors, while network expansion and redispatch is decided by a regulated transmission system operator — as well as the different objectives of firms (profit maximization) and the regulator (welfare maximization). In order to illustrate the possibilities to quantify long term economic effects with our framework, we calibrate our model for the German electricity market. We consider various moderate adjustments of the market design: (i) the division of the market area into two price zones, (ii) the efficient curtailment of renewable production and (iii) a cost-benefit-driven balance between network expansion and network management measures. We then analyze the impact of these market designs on generation and transmission investment in case those design elements are anticipated upon network development planning. The resulting investment and production decisions are compared to a benchmark that reflects the current German electricity market design and to an overall optimal first-best benchmark. Our results reveal that price zones do have a significant impact on locational choice of generators and result in a reduced need for network expansion, but lead to only moderate annual welfare gains of approximately 0.9% of annual total system costs. Anticipation of optimal curtailment of renewables and a cost-benefit-driven use of redispatch operations upon network expansion planning, however, implies a welfare gain of over 4.9% of annual total system costs per year as compared to the existing market design, which equals 85% of the maximum possible welfare gain of the first-best benchmark. KW - Electricity Markets, Network Expansion, Generation Expansion, Investment Incentives, Market Design, Congestion Management, Computational Equilibrium Models Y1 - 2019 U6 - https://doi.org/10.1016/j.eneco.2020.104934 IS - 93 ER - TY - INPR A1 - Liers, Frauke A1 - Martin, Alexander A1 - Merkert, Maximilian A1 - Mertens, Nick A1 - Michaels, Dennis T1 - Solving Mixed-Integer Nonlinear Optimization Problems using Simultaneous Convexification - a Case Study for Gas Networks N2 - Solving mixed-integer nonlinear optimization problems (MINLPs) to global optimality is extremely challenging. An important step for enabling their solution consists in the design of convex relaxations of the feasible set. Known solution approaches based on spatial branch-and-bound become more effective the tighter the used relaxations are. Relaxations are commonly established by convex underestimators, where each constraint function is considered separately. Instead, a considerably tighter relaxation can be found via so-called simultaneous convexification, where convex underestimators are derived for more than one constraint function at a time. In this work, we present a global solution approach for solving mixed-integer nonlinear problems that uses simultaneous convexification. We introduce a separation method that relies on determining the convex envelope of linear combinations of the constraint functions and on solving a nonsmooth convex problem. In particular, we apply the method to quadratic absolute value functions and derive their convex envelopes. The practicality of the proposed solution approach is demonstrated on several test instances from gas network optimization, where the method outperforms standard approaches that use separate convex relaxations. KW - Mixed-Integer Nonlinear Programming KW - Simultaneous Convexification KW - Convex Envelope KW - Gas Network Optimization Y1 - 2020 ER - TY - JOUR A1 - Schewe, Lars A1 - Schmidt, Martin A1 - Thürauf, Johannes T1 - Computing Technical Capacities in the European Entry-Exit Gas Market is NP-Hard JF - Annals of Operations Research N2 - As a result of its liberalization, the European gas market is organized as an entry-exit system in order to decouple the trading and transport of natural gas. Roughly summarized, the gas market organization consists of four subsequent stages. First, the transmission system operator (TSO) is obliged to allocate so-called maximal technical capacities for the nodes of the network. Second, the TSO and the gas traders sign mid- to long-term capacity-right contracts, where the capacity is bounded above by the allocated technical capacities. These contracts are called bookings. Third, on a day-ahead basis, gas traders can nominate the amount of gas that they inject or withdraw from the network at entry and exit nodes, where the nominated amount is bounded above by the respective booking. Fourth and finally, the TSO has to operate the network such that the nominated amounts of gas can be transported. By signing the booking contract, the TSO guarantees that all possibly resulting nominations can indeed be transported. Consequently, maximal technical capacities have to satisfy that all nominations that comply with these technical capacities can be transported through the network. This leads to a highly challenging mathematical optimization problem. We consider the specific instantiations of this problem in which we assume capacitated linear as well as potential-based flow models. In this contribution, we formally introduce the problem of Computing Technical Capacities (CTC) and prove that it is NP-complete on trees and NP-hard in general. To this end, we first reduce the Subset Sum problem to CTC for the case of capacitated linear flows in trees. Afterward, we extend this result to CTC with potential-based flows and show that this problem is also NP-complete on trees by reducing it to the case of capacitated linear flow. Since the hardness results are obtained for the easiest case, i.e., on tree-shaped networks with capacitated linear as well as potential-based flows, this implies the hardness of CTC for more general graph classes. KW - European Entry-Exit Gas Market KW - Technical Capacities KW - Potential-Based Flows KW - Computational Complexity KW - NP-Hardness Y1 - 2020 U6 - https://doi.org/10.1007/s10479-020-03725-2 IS - 295 SP - 337 EP - 362 ER - TY - JOUR A1 - Çelebi, Emre A1 - Krebs, Vanessa A1 - Schmidt, Martin T1 - Γ-Robust Electricity Market Equilibrium Models with Transmission and Generation Investments JF - Energy Systems N2 - We consider uncertain robust electricity market equilibrium problems including transmission and generation investments. Electricity market equilibrium modeling has a long tradition but is, in most of the cases, applied in a deterministic setting in which all data of the model are known. Whereas there exist some literature on stochastic equilibrium problems, the field of robust equilibrium models is still in its infancy. We contribute to this new field of research by considering Γ-robust electricity market equilibrium models on lossless DC networks with transmission and generation investments. We state the nominal market equilibrium problem as a mixed complementarity problem as well as its variational inequality and welfare optimization counterparts. For the latter, we then derive a Γ-robust formulation and show that it is indeed the counterpart of a market equilibrium problem with robustified player problems. Finally, we present two case studies to gain insights into the general effects of robustification on electricity market models. In particular, our case studies reveal that the transmission system operator tends to act more risk-neutral in the robust setting, whereas generating firms clearly behave more risk-averse. KW - Robust optimization KW - Robust market equilibria KW - Electricity market equilibrium models KW - Transmission and generation investment KW - Perfect competition Y1 - 2020 ER - TY - INPR A1 - Aigner, Kevin-Martin A1 - Clarner, Jan-Patrick A1 - Liers, Frauke A1 - Martin, Alexander T1 - Robust Approximation of Chance Constrained DC Optimal Power Flow under Decision-Dependent Uncertainty N2 - We propose a mathematical optimization model and its solution for joint chance constrained DC Optimal Power Flow. In this application, it is particularly important that there is a high probability of transmission limits being satisfied, even in the case of uncertain or fluctuating feed-in from renewable energy sources. In critical network situations where the network risks overload, renewable energy feed-in has to be curtailed by the transmission system operator (TSO). The TSO can reduce the feed-in in discrete steps at each network node. The proposed optimization model minimizes curtailment while ensuring that there is a high probability of transmission limits being maintained. The latter is modeled via (joint) chance constraints that are computationally challenging. Thus, we propose a solution approach based on the robust safe approximation of these constraints. Hereby, probabilistic constraints are replaced by robust constraints with suitably defined uncertainty sets constructed from historical data. The uncertainty sets are calculated by encompassing randomly drawn scenarios using the scenario approach proposed by Margellos et al. (IEEE Transactions on Automatic Control, 59 (2014)). The ability to discretely control the power feed-in then leads to a robust optimization problem with decision-dependent uncertainties, i.e. the uncertainty sets depend on decision variables. We propose an equivalent mixed-integer linear reformulation for box uncertainties with the exact linearization of bilinear terms. Finally, we present numerical results for different test cases from the Nesta archive, as well as for a real network. We consider the discrete curtailment of solar feed-in, for which we use real-world weather and network data. The experimental tests demonstrate the effectiveness of this method and run times are very fast. Moreover, on average the calculated robust solutions lead only to a small increase in curtailment, when compared to nominal solutions. KW - OR in energy KW - optimal power flow KW - chance constrained programming KW - robust optimization KW - decision-dependent uncertainty Y1 - 2020 ER - TY - JOUR A1 - Adelhütte, Dennis A1 - Aßmann, Denis A1 - Gonzàlez Grandòn, Tatiana A1 - Gugat, Martin A1 - Heitsch, Holger A1 - Liers, Frauke A1 - Henrion, René A1 - Nitsche, Sabrina A1 - Schultz, Rüdiger A1 - Stingl, Michael A1 - Wintergerst, David T1 - Joint model of probabilistic/robust (probust) constraints applied to gas network optimization N2 - Optimization tasks under uncertain conditions abound in many real-life applications. Whereas solution approaches for probabilistic constraints are often developed in case the uncertainties can be assumed to follow a certain probability distribution, robust approaches are usually used in case solutions are sought that are feasible for all realizations of uncertainties within some pre-defined uncertainty set. As many applications contain different types of uncertainties that require robust as well as probabilistic treatments, we deal with a class of joint probabilistic/robust constraints as its appears in optimization problems under uncertainty. Focusing on complex uncertain gas network optimization problems, we show the relevance of this class of problems for the task of maximizing free booked capacities in an algebraic model for a stationary gas network. We furthermore present approaches for their solution. Finally, we study the problem of controlling a transient system that is governed by the wave equation. The task consists in determining controls such that a certain robustness measure remains below some given upper bound, with high probability. KW - robust optimization KW - chance constraints KW - optimal control KW - spheric-radial decomposition Y1 - 2017 U6 - https://doi.org/10.1007/s10013-020-00434-y ER - TY - JOUR A1 - Grimm, Veronika A1 - Grübel, Julia A1 - Schewe, Lars A1 - Schmidt, Martin A1 - Zöttl, Gregor T1 - Nonconvex Equilibrium Models for Gas Market Analysis: Failure of Standard Techniques and Alternative Modeling Approaches JF - European Journal on Operational Research N2 - This paper provides a first approach to assess gas market interaction on a network with nonconvex flow models. In the simplest possible setup that adequately reflects gas transport and market interaction, we elaborate on the relation of the solution of a simultaneous competitive gas market game, its corresponding mixed nonlinear complementarity problem (MNCP), and a first-best benchmark. We provide conditions under which the solution of the simultaneous game is also the solution of the corresponding MNCP. However, equilibria cannot be determined by the MNCP as the transmission system operator's (TSO’s) first-order conditions are insufficient, which goes back to nonconvexities of the gas flow model. This also implies that the welfare maximization problem may have multiple solutions that sometimes do not even coincide with any of the market equilibria. Our analysis shows that, even in the absence of strategic firms, market interaction fails to implement desirable outcomes from a welfare perspective due to the TSO’s incentive structure. We conclude that the technical environment calls for a market design that commits the TSO to a welfare objective through regulation and propose a design where the market solution corresponds to a welfare maximum and vice versa. KW - Natural Gas Markets KW - Nonconvex Equilibrium Models KW - Uniqueness KW - Multiplicity KW - Fundamental Welfare Theorems Y1 - 2017 IS - 273(3) SP - 1097 EP - 1108 ER - TY - JOUR A1 - Burlacu, Robert A1 - Egger, Herbert A1 - Groß, Martin A1 - Martin, Alexander A1 - Pfetsch, Marc A1 - Schewe, Lars A1 - Sirvent, Mathias A1 - Skutella, Martin T1 - Maximizing the storage capacity of gas networks: a global MINLP approach N2 - In this paper, we study the transient optimization of gas networks, focusing in particular on maximizing the storage capacity of the network. We include nonlinear gas physics and active elements such as valves and compressors, which due to their switching lead to discrete decisions. The former is described by a model derived from the Euler equations that is given by a coupled system of nonlinear parabolic partial differential equations (PDEs). We tackle the resulting mathematical optimization problem by a first-discretize-then-optimize approach. To this end, we introduce a new discretization of the underlying system of parabolic PDEs and prove well-posedness for the resulting nonlinear discretized system. Endowed with this discretization, we model the problem of maximizing the storage capacity as a non-convex mixed-integer nonlinear problem (MINLP). For the numerical solution of the MINLP, we algorithmically extend a well-known relaxation approach that has already been used very successfully in the field of stationary gas network optimization. This method allows us to solve the problem to global optimality by iteratively solving a series of mixed-integer problems (MIPs). Finally, we present two case studies that illustrate the applicability of our approach. KW - Mixed-Integer Nonlinear Programming KW - Transient Gas Transport Optimization KW - Storage Capacity Maximization KW - Power-to-Gas KW - First-Discretize-Then-Optimize Y1 - 2019 U6 - https://doi.org/10.1007/s11081-018-9414-5 VL - 20 SP - 543 EP - 573 ET - Optimization and Engineering ER - TY - JOUR A1 - Krebs, Vanessa A1 - Schmidt, Martin T1 - Uniqueness of Market Equilibria on Networks with Transport Costs JF - Operations Research Perspectives N2 - We study the existence and uniqueness of equilibria for perfectly competitive markets in capacitated transport networks. The model under consideration is rather general so that it captures basic aspects of related models in, e.g., gas or electricity networks. We formulate the market equilibrium model as a mixed complementarity problem and show the equivalence to a welfare maximization problem. Using the latter we prove uniqueness of the resulting equilibrium for piecewise linear and symmetric transport costs under additional mild assumptions. Moreover, we show the necessity of these assumptions by illustrating examples that possess multiple solutions if our assumptions are violated. KW - Market Equilibria KW - Networks KW - Transport Costs KW - Uniqueness KW - Perfect Competition Y1 - 2017 IS - 5 SP - 169 EP - 173 ER - TY - JOUR A1 - Mehrmann, Volker A1 - Schmidt, Martin A1 - Stolwijk, Jeroen J. T1 - Model and Discretization Error Adaptivity within Stationary Gas Transport Optimization JF - Vietnam Journal of Mathematics N2 - The minimization of operation costs for natural gas transport networks is studied. Based on a recently developed model hierarchy ranging from detailed models of instationary partial differential equations with temperature dependence to highly simplified algebraic equations, modeling and discretization error estimates are presented to control the overall error in an optimization method for stationary and isothermal gas flows. The error control is realized by switching to more detailed models or finer discretizations if necessary to guarantee that a prescribed model and discretization error tolerance is satisfied in the end. We prove convergence of the adaptively controlled optimization method and illustrate the new approach with numerical examples. KW - Gas network optimization KW - Isothermal stationary Euler equations KW - Model hierarchy KW - Adaptive error control KW - Marking strategy Y1 - 2017 IS - 46(4) SP - 779 EP - 801 ER - TY - JOUR A1 - Kleinert, Thomas A1 - Schmidt, Martin T1 - Computing Feasible Points of Bilevel Problems with a Penalty Alternating Direction Method JF - INFORMS Journal on Computing N2 - Bilevel problems are highly challenging optimization problems that appear in many applications of energy market design, critical infrastructure defense, transportation, pricing, etc. Often, these bilevel models are equipped with integer decisions, which makes the problems even harder to solve. Typically, in such a setting in mathematical optimization one develops primal heuristics in order to obtain feasible points of good quality quickly or to enhance the search process of exact global methods. However, there are comparably few heuristics for bilevel problems. In this paper, we develop such a primal heuristic for bilevel problems with mixed-integer linear or quadratic upper level and linear or quadratic lower level. The heuristic is based on a penalty alternating direction method, which allows for a theoretical analysis. We derive a convergence theory stating that the method converges to a stationary point of an equivalent single-level reformulation of the bilevel problem and extensively test the method on a test set of more than 2800 instances - which is one of the largest computational test sets ever used in bilevel programming. The study illustrates the very good performance of the proposed method, both in terms of running times and solution quality. This renders the method a suitable sub-routine in global bilevel solvers as well as a reasonable standalone approach. KW - Bilevel optimization KW - Mixed-integer bilevel optimization KW - Stationary points KW - Penalty methods KW - Alternating direction methods Y1 - 2019 U6 - https://doi.org/10.1287/ijoc.2019.0945 IS - 33 (1) SP - 198 EP - 215 ER - TY - JOUR A1 - Grübel, Julia A1 - Kleinert, Thomas A1 - Krebs, Vanessa A1 - Orlinskaya, Galina A1 - Schewe, Lars A1 - Schmidt, Martin A1 - Thürauf, Johannes T1 - On Electricity Market Equilibria with Storage: Modeling, Uniqueness, and a Distributed ADMM JF - Computers & Operations Research N2 - We consider spot-market trading of electricity including storage operators as additional agents besides producers and consumers. Storages allow for shifting produced electricity from one time period to a later one. Due to this, multiple market equilibria may occur even if classical uniqueness assumptions for the case without storages are satisfied. For models containing storage operators, we derive sufficient conditions that ensure uniqueness of generation and demand. We also prove uniqueness of the market equilibrium for the special case of a single storage operator. Nevertheless, in case of multiple storage operators, uniqueness fails to hold in general, which we show by illustrative examples. We conclude the theoretical discussion with a general ex-post condition for proving the uniqueness of a given solution. In contrast to classical settings without storages, the computation of market equilibria is much more challenging since storage operations couple all trading events over time. For this reason, we propose a tailored parallel and distributed alternating direction method of multipliers (ADMM) for efficiently computing spot-market equilibria over long time horizons. We first analyze the parallel performance of the method itself. Finally, we show that the parallel ADMM clearly outperforms solving the respective problems directly and that it is capable of solving instances with more than 42 million variables in less than 13 minutes. Y1 - 2019 U6 - https://doi.org/10.1016/j.cor.2019.104783 IS - 114 ER - TY - INPR A1 - Schewe, Lars A1 - Schmidt, Martin A1 - Weninger, Dieter T1 - A Decomposition Heuristic for Mixed-Integer Supply Chain Problems T2 - Operations Research Letters N2 - Mixed-integer supply chain models typically are very large but are also very sparse and can be decomposed into loosely coupled blocks. In this paper, we use general-purpose techniques to obtain a block decomposition of supply chain instances and apply a tailored penalty alternating direction method, which exploits the structural properties of the decomposed instances. We further describe problem-specific enhancements of the algorithm and present numerical results on real-world instances that illustrate the applicability of the approach. KW - Supply chain KW - Mixed-integer optimization KW - Decomposition KW - Penalty method KW - Alternating direction methods Y1 - 2019 IS - 48(3) SP - 225 EP - 232 ER - TY - INPR A1 - Kleinert, Thomas A1 - Labbé, Martine A1 - Plein, Fränk A1 - Schmidt, Martin T1 - There's No Free Lunch: On the Hardness of Choosing a Correct Big-M in Bilevel Optimization T2 - Operations Research N2 - One of the most frequently used approaches to solve linear bilevel optimization problems consists in replacing the lower-level problem with its Karush-Kuhn-Tucker (KKT) conditions and by reformulating the KKT complementarity conditions using techniques from mixed-integer linear optimization. The latter step requires to determine some big-M constant in order to bound the lower level's dual feasible set such that no bilevel-optimal solution is cut off. In practice, heuristics are often used to find a big-M although it is known that these approaches may fail. In this paper, we consider the hardness of two proxies for the above mentioned concept of a bilevel-correct big-M. First, we prove that verifying that a given big-M does not cut off any feasible vertex of the lower level's dual polyhedron cannot be done in polynomial time unless P=NP. Second, we show that verifying that a given big-M does not cut off any optimal point of the lower level's dual problem (for any point in the projection of the high-point relaxation onto the leader's decision space) is as hard as solving the original bilevel problem. KW - Bilevel optimization KW - Mathematical programs with complementarity constraints (MPCC) KW - Bounding polyhedra KW - Big-M KW - Hardness Y1 - 2019 IS - 68(6) SP - 1716 EP - 1721 ER - TY - JOUR A1 - Gugat, Martin A1 - Schultz, Rüdiger A1 - Schuster, Michael T1 - Convexity and Starshapedness of Feasible Sets in Stationary Flow Networks JF - Networks and Heterogeneous Media N2 - We deal with a stationary model for flow through a network. The flows are determined by the values at the boundary nodes of the network that we call the loads of the network. In the applications , the feasible loads must satisfy some box constraints. We analyze the structure of the set of feasible loads. Our analysis is motivated by gas pipeline flows, where the box constraints are pressure bounds. We present sufficient conditions to show, that the feasible set is star-shaped with respect to special points. For stronger conditions, we prove the convexity of the set of feasible loads. All the results are given for active and passive networks, i.e. networks with and without inner control. This analysis is motivated by the aim to use the spheric-radial decomposition for stochastic boundary data in this model. This paper can be used for simplifying the algorithmic use of the spheric-radial decomposition. KW - Validation of nominations, Stationary states, Isothermal Euler equations, Gas networks, Convexity, Star-Shapedness, Spheric-radial decomposition Y1 - 2020 U6 - https://doi.org/10.3934/nhm.2020008 SP - 171 EP - 195 ER - TY - JOUR A1 - Gugat, Martin A1 - Schuster, Michael ED - Weber, Gerhard-Wilhelm T1 - Stationary Gas Networks with Compressor Control and Random Loads: Optimization with Probabilistic Constraints JF - Mathematical Problems in Engineering N2 - We introduce a stationary model for gas flow based on simplified isothermal Euler equations in a non-cycled pipeline network. Especially the problem of the feasibility of a random load vector is analyzed. Feasibility in this context means the existence of a flow vectormeeting these loads, which satisfies the physical conservation laws with box constraints for the pressure. An important aspect of the model is the support of compressor stations, which counteract the pressure loss caused by friction in the pipes.The network is assumed to have only one influx node; all other nodes are efflux nodes.With these assumptions the set of feasible loads can be characterized analytically. In addition we show the existence of optimal solutions for some optimization problems with probabilistic constraints. A numerical example based on real data completes this paper. KW - Validation of nominations, Stationary states, Isothermal Euler equations, Ideal Gas, Gas network, Gas transport, Compressor control, Uncertainty, Chance Constraints, Spheric-radial decomposition Y1 - 2019 U6 - https://doi.org/https://doi.org/10.1155/2018/7984079 ER - TY - INPR A1 - Liers, Frauke A1 - Schewe, Lars A1 - Thürauf, Johannes T1 - Radius of Robust Feasibility for Mixed-Integer Problems N2 - For a mixed-integer linear problem (MIP) with uncertain constraints, the radius of robust feasibility (RRF) determines a value for the maximal “size” of the uncertainty set such that robust feasibility of the MIP can be guaranteed. The approaches for the RRF in the literature are restricted to continuous optimization problems. We first analyze relations between the RRF of a MIP and its continuous linear (LP) relaxation. In particular, we derive conditions under which a MIP and its LP relaxation have the same RRF. Afterward, we extend the notion of the RRF such that it can be applied to a large variety of optimization problems and uncertainty sets. In contrast to the setting commonly used in the literature, we consider for every constraint a potentially different uncertainty set that is not necessarily full-dimensional. Thus, we generalize the RRF to MIPs as well as to include “safe” variables and constraints, i.e., where uncertainties do not affect certain variables or constraints. In the extended setting, we again analyze relations between the RRF for a MIP and its LP relaxation. Afterward, we present methods for computing the RRF of LPs as well as of MIPs with safe variables and constraints. Finally, we show that the new methodologies can be successfully applied to the instances in the MIPLIB 2017 for computing the RRF. KW - Robust Optimization, Mixed-integer programming, Uncertainty sets, Robust feasibility Y1 - 2019 U6 - https://doi.org/10.1287/ijoc.2020.1030 PB - Informs Journal on Computing ER - TY - INPR A1 - Brodskyi, Yan A1 - Hante, Falk A1 - Seidel, Arno T1 - Stabilization of Partial Differential Equations by Sequential Action Control N2 - We extend the framework of sequential action control to systems of partial differential equations which can be posed as abstract linear control problems in a Hilbert space. We follow a late-lumping approach and show that the control action can be explicitly obtained from variational principles using adjoint information. Moreover, we analyze the closed-loop system obtained from the SAC feedback for quadratic stage costs. We apply this theory prototypically to an unstable heat equation and verify the results numerically. Y1 - 2019 ER - TY - JOUR A1 - Göttlich, Simone A1 - Hante, Falk A1 - Potschka, Andreas A1 - Schewe, Lars T1 - Penalty alternating direction methods for mixed-integer optimal control with combinatorial constraints JF - Mathematical Programming N2 - We consider mixed-integer optimal control problems with combinatorial constraints that couple over time such as minimum dwell times. We analyze a lifting and decomposition approach into a mixed-integer optimal control problem without combinatorial constraints and a mixed-integer problem for the combinatorial constraints in the control space. Both problems can be solved very efficiently with existing methods such as outer convexification with sum-up-rounding strategies and mixed-integer linear programming techniques. The coupling is handled using a penalty-approach. We provide an exactness result for the penalty which yields a solution approach that convergences to partial minima. We compare the quality of these dedicated points with those of other heuristics amongst an academic example and also for the optimization of electric transmission lines with switching of the network topology for flow reallocation in order to satisfy demands. Y1 - 2019 ER - TY - JOUR A1 - Gatzert, Nadine A1 - Martin, Alexander A1 - Schmidt, Martin A1 - Seith, Benjamin A1 - Vogl, Nikolai T1 - Portfolio Optimization with Irreversible Long-Term Investments in Renewable Energy under Policy Risk: A Mixed-Integer Multistage Stochastic Model and a Moving-Horizon Approach JF - European Journal of Operational Research N2 - Portfolio optimization is an ongoing hot topic of mathematical optimization and management science. Due to the current financial market environment with low interest rates and volatile stock markets, it is getting more and more important to extend portfolio optimization models by other types of investments than classical assets. In this paper, we present a mixed-integer multistage stochastic model that includes investment opportunities in irreversible and long-term infrastructure projects in the context of renewable energies, which are also subject to policy risk. On realistic time scales for investment problems of this type, the resulting instances are by far too large to be solved with today's most evolved optimization software. Thus, we present a tailored moving-horizon approach together with suitable approximations and simplifications of the model. We evaluate these approximations and simplifications in a computational sensitivity analysis and derive a final model that can be tackled on a realistic instance by our moving-horizon approach. KW - Mixed-integer optimization KW - Multistage stochastic optimization KW - Portfolio optimization KW - Illiquid investments KW - Policy risk Y1 - 2019 ER - TY - JOUR A1 - Grimm, Veronika A1 - Orlinskaya, Galina A1 - Schewe, Lars A1 - Schmidt, Martin A1 - Zöttl, Gregor T1 - Optimal Design of Retailer-Prosumer Electricity Tariffs Using Bilevel Optimization JF - Computers & Operations Research N2 - We compare various flexible tariffs that have been proposed to cost-effectively govern a prosumer's electricity management - in particular time-of-use (TOU), critical-peak-pricing (CPP), and a real-time-pricing tariff (RTP). As the outside option, we consider a fixed-price tariff (FP) that restricts the specific characteristics of TOU, CPP, and RTP, so that the flexible tariffs are at least as profitable for the prosumer as the fixed-price tariff. We propose bilevel models to determine the optimal interplay between the retailer's tariff design and the prosumer's decisions on using the storage, on consumption, and on electricity purchases from as well as electricity sales to the grid. The single-level reformulations of the considered bilevel models are computationally highly challenging optimization problems since they, e.g., combine bilinearities and mixed-integer aspects for modeling certain tariff structures. Based on a computational study using real-world data, we find that RTP increases retailer profits, however, leads to the largest price volatility for the prosumer. TOU and CPP only yield mild additional retailer profits and, due to the multiplicity of optimal plans on the part of the prosumer, imply uncertain revenues for the retailer. KW - Electricity tariffs KW - Pricing KW - Bilevel optimization KW - Mixed-integer optimization KW - Tariff design Y1 - 2019 IS - 114 ER - TY - INPR A1 - Hümbs, Lukas A1 - Martin, Alexander A1 - Schewe, Lars T1 - Exploiting complete linear descriptions for decentralized power market problems with integralities N2 - It is well known that linear prices supporting a competitive equilibrium exist in the case of convex markets, however, in the presence of integralities this is open and hard to decide in general. We present necessary and sufficient conditions for the existence of such prices for decentralized market problems where market participants have integral decision variables and their feasible sets are given in complete linear description. We utilize total unimodularity and the aforementioned conditions to show that such linear prices exist and present some applications. Furthermore, we compute competitive equilibria for two classes of decentralized market problems arising in energy markets and show that competitive equilibria may exist regardless of integralities. KW - Mixed-integer programming KW - Competitive equilibrium KW - Linear prices KW - Total unimodularity Y1 - 2019 N1 - Former title: "On the existence of solutions to market equilibrium problems with integralities" ER - TY - JOUR A1 - Robinius, Martin A1 - Schewe, Lars A1 - Schmidt, Martin A1 - Stolten, Detlef A1 - Thürauf, Johannes A1 - Welder, Lara T1 - Robust Optimal Discrete Arc Sizing for Tree-Shaped Potential Networks JF - Computational Optimization and Applications N2 - We consider the problem of discrete arc sizing for tree-shaped potential networks with respect to infinitely many demand scenarios. This means that the arc sizes need to be feasible for an infinite set of scenarios. The problem can be seen as a strictly robust counterpart of a single-scenario network design problem, which is shown to be NP-complete even on trees. In order to obtain a tractable problem, we introduce a method for generating a finite scenario set such that optimality of a sizing for this finite set implies the sizing's optimality for the originally given infinite set of scenarios. We further prove that the size of the finite scenario set is quadratically bounded above in the number of nodes of the underlying tree and that it can be computed in polynomial time. The resulting problem can then be solved as a standard mixed-integer linear optimization problem. Finally, we show the applicability of our theoretical results by computing globally optimal arc sizes for a realistic hydrogen transport network of Eastern Germany. KW - Discrete arc sizing KW - Mixed-integer linear optimization KW - Potential networks KW - Scenario generation KW - Robust optimization Y1 - 2018 U6 - https://doi.org/10.1007/s10589-019-00085-x IS - 73(3) SP - 791 EP - 819 ER - TY - JOUR A1 - Kleinert, Thomas A1 - Schmidt, Martin T1 - Global Optimization of Multilevel Electricity Market Models Including Network Design and Graph Partitioning JF - Discrete Optimization N2 - We consider the combination of a network design and graph partitioning model in a multilevel framework for determining the optimal network expansion and the optimal zonal configuration of zonal pricing electricity markets, which is an extension of the model discussed in [25] that does not include a network design problem. The two classical discrete optimization problems of network design and graph partitioning together with nonlinearities due to economic modeling yield extremely challenging mixed-integer nonlinear multilevel models for which we develop two problem-tailored solution techniques. The first approach relies on an equivalent bilevel formulation and a standard KKT transformation thereof including novel primal-dual bound tightening techniques, whereas the second is a tailored generalized Benders decomposition. For the latter, we strengthen the Benders cuts of [25] by using the structure of the newly introduced network design subproblem. We prove for both methods that they yield global optimal solutions. Afterward, we compare the approaches in a numerical study and show that the tailored Benders approach clearly outperforms the standard KKT transformation. Finally, we present a case study that illustrates the economic effects that are captured in our model. KW - Network design KW - Graph partitioning KW - Multilevel optimization KW - Mixed-integer optimization KW - Electricity market design Y1 - 2018 IS - 33 SP - 43 EP - 69 ER - TY - JOUR A1 - Schmidt, Martin A1 - Sirvent, Mathias A1 - Wollner, Winnifried T1 - The Cost of Not Knowing Enough: Mixed-Integer Optimization with Implicit Lipschitz Nonlinearities JF - Optimization Letters N2 - It is folklore knowledge that nonconvex mixed-integer nonlinear optimization problems can be notoriously hard to solve in practice. In this paper we go one step further and drop analytical properties that are usually taken for granted in mixed-integer nonlinear optimization. First, we only assume Lipschitz continuity of the nonlinear functions and additionally consider multivariate implicit constraint functions that cannot be solved for any parameter analytically. For this class of mixed-integer problems we propose a novel algorithm based on an approximation of the feasible set in the domain of the nonlinear function---in contrast to an approximation of the graph of the function considered in prior work. This method is shown to compute approximate global optimal solutions in finite time and we also provide a worst-case iteration bound. In some first numerical experiments we show that the ``cost of not knowing enough'' is rather high by comparing our approach with the open-source global solver SCIP. This reveals that a lot of work is still to be done for this highly challenging class of problems and we thus finally propose some possible directions of future research. KW - Mixed-Integer Nonlinear Optimization, Global Optimization, Lipschitz Optimization, Gas Networks Y1 - 2018 ER - TY - JOUR A1 - Gugat, Martin A1 - Rosier, Lionel A1 - Perrollaz, Vincent ED - Arendt, W. ED - Pierre, M. T1 - Boundary stabilization of quasilinear hyperbolic systems of balance laws: exponential decay for small source terms JF - Journal of Evolution Equations N2 - We investigate the long-time behaviour of solutions of quasilinear hyperbolic systems with transparent boundary conditions when small source terms are incorporated in the system. Even if the finite-time stability of the system is not preserved, it is shown here that an exponential convergence towards the steady state still holds with a decay rate which is proportional to the logarithm of the amplitude of the source term. The result is stated for a system with dynamical boundary conditions in order to deal with initial data that are free of any compatibility condition. The proof of the existence and uniqueness of a solution defined for all positive times is also provided in this paper. Y1 - 2018 U6 - https://doi.org/10.1007/s00028-018-0449-z ER - TY - JOUR A1 - Schewe, Lars A1 - Schmidt, Martin T1 - The Impact of Potential-Based Physics Models on Pricing in Energy Networks JF - Central European Journal of Operations Research N2 - Pricing of access to energy networks is an important issue in liberalized energy sectors because of the natural monopoly character of the underlying transport infrastructures. We introduce a general pricing framework for potential-based energy flows in arbitrarily structured transport networks. In different specifications of our general pricing model we discuss first- and second-best pricing results and compare different pricing outcomes of potential-free and potential-based energy flow models. Our results show that considering nonlinear laws of physics leads to significantly different pricing results on networks and that these differences can only be seen in sufficiently complex, e.g., cyclic, networks as they can be found in real-world situations. KW - Energy Networks KW - Pricing KW - Gas Networks KW - Electricity Networks Y1 - 2018 IS - 28(3) SP - 1003 EP - 1029 ER - TY - INPR A1 - Aßmann, Denis A1 - Liers, Frauke A1 - Stingl, Michael T1 - Decomposable Robust Two-Stage Optimization: An Application to Gas Network Operations Under Uncertainty N2 - We study gas network problems with compressors and control valves under uncertainty that can be formulated as two-stage robust optimization problems. Uncertain data are present in the physical parameters of the pipes as well as in the overall demand. We show how to exploit the special decomposable structure of the problem in order to reformulate the two-stage robust problem as a standard single-stage optimization problem. Since this structure is present in similar problems on e.g., water or direct current electricity networks, we investigate the consequences of the decomposable structure in an abstract setting: The right-hand side of the single-stage problem can be precomputed by solving a series of optimization problems and multiple elements of the right-hand side can be combined into one optimization task. In order to apply our results to gas network problems, we extend piecewise relaxations and preprocessing techniques to incorporate uncertain input data. The practical feasibility and effectiveness of our approach is demonstrated with benchmarks on realistic gas network instances. We observe large speedups due to the described aggregation method together with the developed preprocessing strategies. Furthermore, we are able to solve even comparably large gas network instances quickly for the price of slightly more conservative solutions. KW - robust optimization KW - gas networks KW - relaxations Y1 - 2017 ER - TY - INPR A1 - Gugat, Martin A1 - Schuster, Michael A1 - Sokolowski, Jan T1 - Location Problem for Compressor Stations in Pipeline Networks N2 - In the operation of pipeline networks, compressors play a crucial role in ensuring the network’s functionality for various scenarios. In this contribution we address the important question of finding the optimal location of the compressors. This problem is of a novel structure, since it is related with the gas dynamics that governs the network flow. That results in non-convex mixed integer stochastic optimization problems with probabilistic constraints. Using a steady state model for the gas flow in pipeline networks including compressor control and uncertain loads given by certain probability distributions, the problem of finding the optimal location for the control on the network, s.t. the control cost is minimal and the gas pressure stays within given bounds, is considered. In the deterministic setting, explicit bounds for the pipe length and the inlet pressure, s.t. a unique optimal compressor location with minimal control cost exists, are presented. In the probabilistic setting, an existence result for the optimal compressor location is presented and the uniqueness of the solution is discussed depending on the probability distribution. For Gaussian distributed loads a uniqueness result for the optimal compressor location is presented. Further the problem of finding the optimal compressor locations on networks including the number of compressor stations as variable is considered. Results for the existence of optimal locations on a graph in both, the deterministic and the probabilistic setting, are presented and the uniqueness of the solutions is discussed depending on probability distributions and graph topology. The paper concludes with an illustrative example demonstrating that the compressor locations determined using a steady state approach are also admissible in transient settings. KW - Gas Networks KW - Compressor Control KW - Weber Problem KW - Optimal Location KW - Uncertain Boundary Data Y1 - 2024 ER - TY - INPR A1 - Schuster, Michael T1 - On the Convergence of Optimization Problems with Kernel Density Estimated Probabilistic Constraints N2 - Uncertainty plays a significant role in applied mathematics and probabilistic constraints are widely used to model uncertainty in various fields, even if probabilistic constraints often demand computational challenges. Kernel density estimation (KDE) provides a data-driven approach for properly estimating probability density functions and efficiently evaluate corresponding probabilities. In this paper, we investigate optimization problems with probabilistic constraints, where the probabilities are approximated using a KDE approach. We establish sufficient conditions under which the solution of the KDE approximated optimization problem converges to the solution of the original problem as the sample size goes to infinity. The main results of this paper include three theorems: (1) For sufficiently large sample sizes, the solution of the original problem is also a solution of the approximated problem, if the probabilistic constraint is passive; (2) The limit of a convergent sequence of solutions of the approximated problems is a solution of the original problem, if the KDE uniformly converges; (3) We provide sufficient conditions for the existence of a convergent sequence of solutions of the approximated problems. KW - Probabilistic Constrained Optimization KW - Stochastic Optimization KW - Chance Constraints KW - Kernel Density Estimation KW - Convergence Analysis Y1 - 2024 ER - TY - INPR A1 - Lange, Christian T1 - Modeling and Optimal Control of the Flow of a Gas Mixture N2 - We consider the Euler equations for a pipeline flow of a mixture of two gases. An important application is hydrogen blending. Existence and uniqueness of semi-global solutions is shown and possible boundary conditions are analyzed. Secondly, we consider classes of associated optimal control problems and show existence of solutions. KW - Euler Equations KW - quasilinear hyperbolic system KW - gas transport modeling KW - hydrogen blending Y1 - 2024 ER - TY - JOUR A1 - Gugat, Martin A1 - Qian, Meizhi A1 - Sokolowski, Jan T1 - Network Design and Control: Shape and Topology Optimization for the Turnpike Property for the Wave Equation N2 - The optimal control problems for the wave equation are considered on networks. The turnpike property is shown for the state equation, the adjoint state equation as well as the optimal cost. The shape and topology optimization is performed for the network with the shape functional given by the optimality system of the control problem. The set of admissible shapes for the network is compact in finite dimensions, thus the use of turnpike property is straightforward. The topology optimization is analysed for an example of nucleation of a small cycle at the internal node of network. The topological derivative of the cost is introduced and evaluated in the framework of domain decomposition technique. Numerical examples are provided. KW - Turnpike Property KW - Wave Equation KW - Optimal Control KW - Spectral Method KW - Network Optimum Design Y1 - 2024 VL - J. Geom. Anal. IS - 34 ER - TY - JOUR A1 - Filipkovska, Maria T1 - Qualitative analysis of nonregular differential-algebraic equations and the dynamics of gas networks N2 - Conditions for the existence, uniqueness and boundedness of global solutions, as well as ultimate boundedness of solutions, and conditions for the blow-up of solutions of nonregular semilinear differential-algebraic equations have been obtained. An example demonstrating the application of the obtained results has been considered. Isothermal models of gas networks have been proposed as applications. KW - nonregular differential-algebraic equation KW - degenerate differential equation KW - singular pencil KW - gas network KW - global solvability Y1 - 2023 U6 - https://doi.org/https://doi.org/10.15407/mag19.04.719 VL - 19 IS - 4 SP - 719 EP - 765 ER - TY - JOUR A1 - Schuster, Michael A1 - Gugat, Martin A1 - Sokolowski, Jan T1 - The Location Problem for Compressor Stations in Pipeline Networks JF - Mathematics and Mechanics of Complex Systems N2 - In the operation of pipeline networks, compressors play a crucial role in ensuring the network’s functionality for various scenarios. In this contribution we address the important question of finding the optimal location of the compressors. This problem is of a novel structure, since it is related to the gas dynamics that governs the network flow. That results in nonconvex mixed integer stochastic optimization problems with probabilistic constraints. Using a steady state model for the gas flow in pipeline networks including compressor control and uncertain loads given by certain probability distributions, we consider the problem of finding the optimal location for the control on the network such that the control cost is minimal and the gas pressure stays within given bounds. In the deterministic setting, we present explicit bounds for the pipe length and the inlet pressure such that a unique optimal compressor location with minimal control cost exists. In the probabilistic setting, we give an existence result for the optimal compressor location and discuss the uniqueness of the solution depending on the probability distribution. For Gaussian distributed loads a uniqueness result for the optimal compressor location is presented. We further present the problem of finding optimal compressor locations on networks including the number of compressor stations as a variable. Results for the existence of optimal locations on a graph in both the deterministic and the probabilistic setting are presented, and the uniqueness of the solutions is discussed depending on probability distributions and graph topology. The paper concludes with an illustrative example on a diamond graph demonstrating that the minimal number of compressor stations is not necessarily equal to the optimal number of compressor stations. KW - gas network KW - compressor control KW - Weber problem KW - uncertain boundary data KW - non convex mixed integer stochastic problem Y1 - 2024 U6 - https://doi.org/10.2140 VL - 12 IS - 4 SP - 507 EP - 546 ER - TY - JOUR A1 - Gugat, Martin T1 - New Lyapunov functions for systems with source terms JF - Control and Cybernetics N2 - Lyapunov functions with exponential weights have been used successfully as a powerful tool for the stability analysis of hyperbolic systems of balance laws. In this paper we extend the class of weight functions to a family of hyperbolic functions and study the advantages in the analysis of 2 × 2 systems of balance laws. We present cases connected with the study of the limit of stabilizability, where the new weights provide Lyapunov functions that show exponential stability for a larger set of problem parameters than classical exponential weights. Moreover, we show that sufficiently large time-delays influence the limit of stabilizability in the sense that the parameter set, for which the system can be stabilized becomes substantially smaller. We also demonstrate that the hyperbolic weights are useful in the analysis of the boundary feedback stability of systems of balance laws that are governed by quasilinear hyperbolic partial differential equations. KW - Lyapunov Function KW - Exponential Weight KW - Hyperbolic Weight KW - Feedback law KW - Stabilization Y1 - 2025 U6 - https://doi.org/https://doi.org/10.2478/candc-2024-0008 VL - 53 IS - 1 SP - 163 EP - 187 ER - TY - JOUR A1 - Groß, Martin A1 - Marc E., Pfetsch A1 - Schewe, Lars A1 - Schmidt, Martin A1 - Skutella, Martin T1 - Algorithmic Results for Potential-Based Flows: Easy and Hard Cases N2 - 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. KW - Potential networks KW - Potential-based flows KW - Maximum flow problem KW - Series-parallel graphs KW - Network reduction Y1 - 2017 U6 - https://doi.org/10.1002/net.21865 VL - 73 IS - 3 SP - 303 EP - 324 ET - Networks ER - TY - INPR A1 - Wintergerst, David T1 - Application of Chance Constrained Optimization to Gas Networks N2 - We consider optimization problems with a joint probabilistic constraint under normally distributed uncertain parameters. The parametric constraints are replaced by one constraint stating that the probability of being feasible shall exceed or be equal to a prescribed threshold. In order to apply the concept to gas network optimization under uncertain boundary flows, which corresponds to the demand of customers, we derive an analytic gradient formula. The integral corresponding to the probability can be parameterized by spherical radial decomposition. For this parameterization gradient formulas are known under convexity assumptions of the parametric constraints in the parameter. For the application in gas networks that we have in mind, the convexity assumption of the parametric constraints is not satisfied. Therefore, we weaken it to convexity of the region of feasible parameters for a fixed optimization variable. We proceed to show that the assumptions needed for the gradient formula are met in the gas network optimization problem on a tree. For the numerical implementation we propose a multilevel sampling algorithm that uses a coarse approximation of the chance constraint to generate a warm start for the expensive approximation with fine sampling. The numerical results illustrate that this approach significantly reduces the computation time. KW - chance constraint KW - spherical radial decomposition KW - isothermal Euler equations KW - gas networks KW - multilevel Y1 - 2017 ER - TY - JOUR A1 - Schmidt, Martin A1 - Sirvent, Mathias A1 - Wollner, Winnifried T1 - A Decomposition Method for MINLPs with Lipschitz Continuous Nonlinearities JF - Mathematical Programming N2 - 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. KW - Mixed-Integer Nonlinear Optimization, Lipschitz Optimization, Inexact Function Evaluations, Decomposition Methods, Gas Networks Y1 - 2017 IS - 178(1) SP - 449 EP - 483 ER - TY - INPR A1 - Egger, Herbert A1 - Kugler, Thomas A1 - Liljegren-Sailer, Björn A1 - Marheineke, Nicole A1 - Mehrmann, Volker T1 - On structure preserving model reduction for damped wave propagation in transport networks N2 - We consider the discretization and subsequent model reduction of a system of partial differential-algebraic equations describing the propagation of pressure waves in a pipeline network. Important properties like conservation of mass, dissipation of energy, passivity, existence of steady states, and exponential stability can be preserved by an appropriate semi- discretization in space via a mixed finite element method and also during the further dimension reduction by structure preserving Galerkin projection which is the main focus of this paper. Krylov subspace methods are employed for the construction of the reduced models and we discuss modifications needed to satisfy certain algebraic compatibility conditions; these are required to ensure the well-posedness of the reduced models and the preservation of the key properties. Our analysis is based on the underlying infinite dimensional problem and its Galerkin approximations. The proposed algorithms therefore have a direct interpretation in function spaces; in principle, they are even applicable directly to the original system of partial differential-algebraic equations while the intermediate discretization by finite elements is only required for the actual computations. The performance of the proposed methods is illustrated with numerical tests and the necessity for the compatibility conditions is demonstrated by examples. KW - partial differential-algebraic equations KW - port-Hamiltonian systems KW - Galerkin projection Y1 - 2017 ER - TY - INPR A1 - Burlacu, Robert A1 - Geißler, Björn A1 - Schewe, Lars T1 - Solving Mixed-Integer Nonlinear Programs using Adaptively Refined Mixed-Integer Linear Programs N2 - We propose a method for solving mixed-integer nonlinear programs (MINLPs) to global optimality by discretization of occuring nonlinearities. The main idea is based on using piecewise linear functions to construct mixed-integer linear program (MIP) relaxations of the underlying MINLP. In order to find a global optimum of the given MINLP we develope an iterative algorithm which solves MIP relaxations that are adaptively refined. We are able to give convergence results for a wide range of MINLPs requiring only continuous nonlinearities with bounded domains and an oracle computing maxima of the nonlinearities on their domain. Moreover, the practicalness of our approach is shown numerically by an application from the field of gas network optimization. KW - Mixed-Integer Nonlinear Programming KW - Piecewise Linearization Y1 - 2017 ER - TY - JOUR A1 - Gugat, Martin A1 - Leugering, Günter A1 - Wang, Ke ED - Zhang, Xu T1 - Neumann boundary feedback stabilization for a nonlinear wave equation: A strict H2-Lyapunov function JF - Mathematical Control and Related Fields (MCRF) N2 - For a system that is governed by the isothermal Euler equations with friction for ideal gas, the corresponding field of characteristic curves is determined by the velocity of the flow. This velocity is determined by a second-order quasilinear hyperbolic equation. For the corresponding initial-boundary value problem with Neumann-boundary feedback, we consider non-stationary solutions locally around a stationary state on a finite time interval and discuss the well-posedness of this kind of problem. We introduce a strict H2-Lyapunov function and show that the boundary feedback constant can be chosen such that the H2-Lyapunov function and hence also the H2-norm of the difference between the non-stationary and the stationary state decays exponentially with time. KW - Boundary feedback control, feedback stabilization, exponential stability, isothermal Euler equations, second-order quasilinear equation, Lyapunov function, stationary state, non-stationary state, gas pipeline. Y1 - 2017 U6 - https://doi.org/10.3934/mcrf.2017015 VL - 7 IS - 3 SP - 419 EP - 448 ER - TY - JOUR A1 - Gugat, Martin A1 - Keimer, Alexander A1 - Leugering, Günter A1 - Wang, Zhiqiang ED - Piccoli, Benedetto T1 - Analysis of a system of nonlocal conservation laws for multi-commodity flow on networks JF -  Networks and Heterogeneous Media N2 - We consider a system of scalar nonlocal conservation laws on networks that model a highly re-entrant multi-commodity manufacturing system as encountered in semiconductor production. Every single commodity is mod-eled by a nonlocal conservation law, and the corresponding PDEs are coupled via a collective load, the work in progress. We illustrate the dynamics for two commodities. In the applications, directed acyclic networks naturally occur, therefore this type of networks is considered. On every edge of the network we have a system of coupled conservation laws with nonlocal velocity. At the junctions the right hand side boundary data of the foregoing edges is passed as left hand side boundary data to the following edges and PDEs. For distributing junctions, where we have more than one outgoing edge, we impose time dependent distribution functions that guarantee conservation of mass. We provide results of regularity, existence and well-posedness of the multi-commodity network model for L p-, BV-and W 1,p-data. Moreover, we define an L 2-tracking type objective and show the existence of minimizers that solve the corresponding optimal control problem. KW - conservation laws on network KW - nonlocal conservation laws KW - optimal nodal control KW - systems of hyperbolic pdes Y1 - 2016 U6 - https://doi.org/DOI: 10.3934/nhm.2015.10.749 VL - 10 IS - 4 SP - 749 EP - 785 ER - TY - JOUR A1 - Gugat, Martin A1 - Leugering, Günter ED - Zuazua, Enrique T1 - Time delay in optimal control loops for wave equations JF - ESAIM: COCV N2 - In optimal control loops delays can occur, for example through transmission via digital communication channels. Such delays influence the state that is generated by the implemented control. We study the effect of a delay in the implementation of L 2-norm minimal Neumann boundary controls for the wave equation. The optimal controls are computed as solutions of problems of exact optimal control, that is if they are implemented without delay, they steer the system to a position of rest in a given finite time T. We show that arbitrarily small delays δ > 0 can have a destabilizing effect in the sense that we can find initial states such that if the optimal control u is implemented in the form yx(t, 1) = u(t − δ) for t > δ, the energy of the system state at the terminal time T is almost twice as big as the initial energy. We also show that for more regular initial states, the effect of a delay in the implementation of the optimal control is bounded above in the sense that for initial positions with derivatives of BV-regularity and initial velocities with BV-regularity, the terminal energy is bounded above by the delay δ multiplied with a factor that depends on the BV-norm of the initial data. We show that for more general hyperbolic optimal exact control problems the situation is similar. For systems that have arbitrarily large eigenvalues, we can find terminal times T and arbitrarily small time delays δ, such that at the time T + δ, in the optimal control loop with delay the norm of the state is twice as large as the corresponding norm for the initial state. Moreover, if the initial state satisfies an additional regularity condition, there is an upper bound for the effect of time delay of the order of the delay with a constant that depends on the initial state only. KW - PDE constrained optimization KW - delay KW - wave equation KW - boundary control KW - hyperbolic system Y1 - 2016 U6 - https://doi.org/http://dx.doi.org/10.1051/cocv/2015038 ER - TY - JOUR A1 - Gugat, Martin A1 - Wintergerst, David A1 - Schultz, Rüdiger ED - Iske, Armin T1 - Networks of pipelines for gas with nonconstant compressibility factor: stationary states JF - Computational and Applied Mathematics N2 - For the management of gas transportation networks, it is essential to know how the stationary states of the system are determined by the boundary data. The isothermal Euler equations are an accurate pde-model for the gas flow through each pipe. A compressibility factor is used to model the nonlinear relationship between density and pressure that occurs in real gas in contrast to ideal gas. The gas flow through the nodes is governed by algebraic node conditions that require the conservation of mass and the continuity of the pressure. We examine networks that are described by arbitrary finite graphs and show that for suitably chosen boundary data, subsonic stationary states exist and are uniquely determined by the boundary data. Our construction of the stationary states is based upon explicit representations of the stationary states on each single pipe that can easily be evaluated numerically. We also use the monotonicity properties of these states as functions of the boundary data. Y1 - 2016 U6 - https://doi.org/10.1007/s40314-016-0383-z ER - TY - JOUR A1 - Bärmann, Andreas A1 - Liers, Frauke A1 - Martin, Alexander A1 - Merkert, Maximilian A1 - Thurner, Christoph A1 - Weninger, Dieter T1 - Solving network design problems via iterative aggregation JF - Mathematical Programming Computation N2 - In this work, we present an exact approach for solving network design problems that is based on an iterative graph aggregation procedure. The scheme allows existing preinstalled capacities. Starting with an initial aggregation, we solve a sequence of network design master problems over increasingly fine-grained representations of the original network. In each step, a subproblem is solved that either proves optimality of the solution or gives a directive where to refine the representation of the network in the subsequent iteration. The algorithm terminates with a globally optimal solution to the original problem. Our implementation uses a standard integer programming solver for solving the master problems as well as the subproblems. The computational results on random and realistic instances confirm the profitable use of the iterative aggregation technique. The computing time often reduces drastically when our method is compared to solving the original problem from scratch. KW - Aggregation KW - Network design KW - Combinatorial optimization KW - Mixed-integer programming KW - Branch-and-cut Y1 - 2015 U6 - https://doi.org/10.1007/s12532-015-0079-1 VL - 7 IS - 2 SP - 189 EP - 217 ER - TY - JOUR A1 - Liers, Frauke A1 - Merkert, Maximilian T1 - Structural Investigation of Piecewise Linearized Network Flow Problems N2 - In this work we study polyhedra in the context of network flow problems, where the flow value on each arc lies in one of several predefined intervals. This is motivated by nonlinear problems on transportation networks, where nonlinearities are handled by piecewise linear approximation or relaxation - a common and established approach in many applications. Several methods for modeling piecewise linear functions are known which provide a complete description for a single network arc. However, in general this property is lost when considering multiple arcs. We show how to strengthen the formulation for specific substructures consisting of multiple arcs by linear inequalities. For the case of paths of degree-two-nodes we give a complete description of the polyhedron projected to the integer variables. Our model is based on - but not limited to - the multiple choice method; we also show how to transfer our results to a formulation based on the incremental method. Computational results show that a state-of-the-art MIP-solver greatly benefits from using our cutting planes for random and realistic network topologies. KW - Combinatorial optimization KW - Complete description KW - Network flow problems KW - Piecewise linear functions Y1 - 2016 U6 - https://doi.org/10.1137/15M1006751 VL - 26 SP - 2863 EP - 2886 ER - TY - INPR A1 - Bärmann, Andreas A1 - Gellermann, Thorsten A1 - Merkert, Maximilian A1 - Schneider, Oskar T1 - Staircase Compatibility and its Applications in Scheduling and Piecewise Linearization N2 - We consider the clique problem with multiple-choice constraints (CPMC) and characterize a case where it is possible to give an efficient description of the convex hull of its feasible solutions. This case, which we call staircase compatibility, generalizes common properties in applications and allows for a linear description of the integer feasible solutions to (CPMC) with a totally unimodular constraint matrix of polynomial size. We derive two such totally unimodular reformulations for the problem: one that is obtained by a strengthening of the compatibility constraints and one that is based on a representation as a dual network flow problem. Furthermore, we show a natural way to derive integral solutions from fractional solutions to the problem by determining integral extreme points generating this fractional solution. We also evaluate our reformulations from a computational point of view by applying them to two different real-world applications. The first one is a problem in railway timetabling where we try to adapt a given timetable slightly such that energy costs from operating the trains are reduced. The second one is the piecewise linearization of non-linear flow problems on a gas network. In both cases, we are able to reduce the solution times significantly by passing to the theoretically stronger formulations of the problem. KW - Clique Problem KW - Multiple-Choice Constraints KW - Total Unimodularity KW - Scheduling KW - Piecewise Linearization Y1 - 2016 ER - TY - JOUR A1 - Hante, Falk A1 - Mommer, Mario A1 - Potschka, Andreas T1 - Newton-Picard preconditioners for time-periodic, parabolic optimal control problems JF - SIAM Journal on Numerical Analysis Y1 - 2016 U6 - https://doi.org/10.1137/140967969 VL - 53 IS - 5 SP - 2206 EP - 2225 ER - TY - JOUR A1 - Domschke, Pia A1 - Groß, Martin A1 - Hiller, Benjamin A1 - Hante, Falk A1 - Schewe, Lars A1 - Schmidt, Martin T1 - Mathematische Modellierung, Simulation und Optimierung von Gastransportnetzwerken JF - gwf-gas/erdgas Y1 - 2015 VL - 11 SP - 880 EP - 885 ER - TY - JOUR A1 - Schmidt, Martin A1 - Aßmann, Denis A1 - Burlacu, Robert A1 - Humpola, Jesco A1 - Joormann, Imke A1 - Kanelakis, Nikolaos A1 - Koch, Thorsten A1 - Oucherif, Djamal A1 - Pfetsch, Marc E. A1 - Schewe, Lars A1 - Schwarz, Robert A1 - Sirvent, Mathias T1 - GasLib – A Library of Gas Network Instances JF - Data N2 - 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. KW - Gas Transport KW - Networks KW - Problem Instances KW - Mixed-Integer Nonlinear Optimization KW - GasLib Y1 - 2017 U6 - https://doi.org/10.3390/data2040040 VL - 4 IS - 2 ER - TY - JOUR A1 - Hante, Falk A1 - Schmidt, Martin T1 - Complementarity-Based Nonlinear Programming Techniques for Optimal Mixing in Gas Networks JF - EURO Journal on Computational Optimization N2 - We consider nonlinear and nonsmooth mixing aspects in gas transport optimization problems. As mixed-integer reformulations of pooling-type mixing models already render small-size instances computationally intractable, we investigate the applicability of smooth nonlinear programming techniques for equivalent complementarity-based reformulations. Based on recent results for remodeling piecewise affine constraints using an inverse parametric quadratic programming approach, we show that classical stationarity concepts are meaningful for the resulting complementarity-based reformulation of the mixing equations. Further, we investigate in a numerical study the performance of this reformulation compared to a more compact complementarity-based one that does not feature such beneficial regularity properties. All computations are performed on publicly available data of real-world size problem instances from steady-state gas transport. KW - Gas transport networks KW - Mixing KW - Inverse parametric quadratic programming KW - Complementarity constraints KW - MPCC Y1 - 2017 IS - 7(3) SP - 299 EP - 323 ER - TY - CHAP A1 - Rüffler, Fabian A1 - Hante, Falk T1 - Optimality Conditions for Switching Operator Differential Equations N2 - We consider optimal switching of hybrid abstract evolution equations. The framework includes switching semilinear partial differential equations of parabolic or hyperbolic type, discontinuous state resets, switching costs and allows switching of the principle parts of the equations. We present adjoint-based formulae for the gradient of the cost functional with respect to position and number of switching time points that lead to first order necessary conditions. Moreover, we discuss an alternate-direction approach for implementing descent methods. As an application we consider optimal open/close-switching of valves and on/off-switching control of compressors in a gas network modelled by a graph with simplified euler equations on edges and suitable coupling conditions at nodes. Y1 - 2017 ER - TY - CHAP A1 - Gugat, Martin A1 - Herty, Michael A1 - Yu, Hui ED - Westdickenberg, Michael T1 - On the relaxation approximation for 2 × 2 hyperbolic balance laws N2 - The relaxation approximation for systems of conservation laws has been studied intensively for example by [17, 5, 19, 24]. In this paper the corresponding relaxation approximation for 2x2 systems of balance laws is studied. Our driving example is gas flow in pipelines described by the isothermal Euler equations. We are interested in the limiting behavior as the relaxation parameter tends to zero. We give conditions where the relaxation converges to the states of the original system and counterexamples for cases where the steady states depend on the space variable. Y1 - 2017 VL - x IS - x SP - x EP - x PB - Springer CY - Springer Proceedings in Mathematics and Statistics ER - TY - JFULL A1 - Leugering, Günter T1 - Domain Decomposition of an Optimal Control Problem for Semi-Linear Elliptic Equations on Metric Graphs with Application to Gas Networks N2 - We consider optimal control problems for the flow of gas in a pipe network. The equations of motions are taken to be represented by a semi-linear model derived from the fully nonlinear isothermal Euler gas equations. We formulate an optimal control problem on a given network and introduce a time discretization thereof. We then study the well-posedness of the corresponding time-discrete optimal control problem. In order to further reduce the complexity, we consider an instantaneous control strategy. The main part of the paper is concerned with a non-overlapping domain decomposition of the semi-linear elliptic optimal control problem on the graph into local problems on a small part of the network, ultimately on a single edge. KW - 28 KW - nonoverlapping omain decomposition, optimal control of semi-linear ellioptic systems on netowrks Y1 - 2017 U6 - https://doi.org/https://doi.org/10.4236/am.2017.88082 SN - 2152-7393 VL - 8 ER - TY - JOUR A1 - Krebs, Vanessa A1 - Schewe, Lars A1 - Schmidt, Martin T1 - Uniqueness and Multiplicity of Market Equilibria on DC Power Flow Networks JF - European Journal on Operations Research N2 - We consider uniqueness and multiplicity of market equilibria in a short-run setup where traded quantities of electricity are transported through a capacitated network in which power flows have to satisfy the classical lossless DC approximation. The firms face fluctuating demand and decide on their production, which is constrained by given capacities. Today, uniqueness of such market outcomes are especially important in more complicated multilevel models for measuring market (in)efficiency. Thus, our findings are important prerequisites for such studies. We show that market equilibria are unique on tree networks under mild assumptions and we also present a priori conditions under which equilibria are unique on cycle networks. On general networks, uniqueness fails to hold and we present simple examples for which multiple equilibria exist. However, we prove a posteriori criteria for the uniqueness of a given solution and characterize situations in which multiple solutions exist. KW - Market Equilibria KW - Uniqueness KW - Multiplicity KW - Networks KW - DC Power Flow Y1 - 2017 IS - 271(1) SP - 165 EP - 178 ER - TY - JOUR A1 - Gugat, Martin A1 - Ulbrich, Stefan ED - CHEN, GOONG T1 - The isothermal Euler equations for ideal gas with source term: Product solutions, flow reversal and no blow up JF - Journal of Mathematical Analysis and Applications KW - Global classical solutions Ideal gas Bi-directional flow Transsonic flow Y1 - 2017 U6 - https://doi.org/10.1016/j.jmaa.2017.04.064 VL - 454 IS - 1 SP - 439 EP - 452 ER - TY - JOUR A1 - Gugat, Martin A1 - Zuazua, Enrique ED - Grimble, Mike J. T1 - Exact penalization of terminal constraints for optimal control problems JF - OPTIMAL CONTROL APPLICATIONS AND METHODS N2 - We study optimal control problems for linear systems with prescribed initial and terminal states. We analyze the exact penalization of the terminal constraints. We show that for systems that are exactly controllable, the norm-minimal exact control can be computed as the solution of an optimization problem without terminal constraint but with a nonsmooth penalization of the end conditions in the objective function, if the penalty parameter is sufficiently large. We describe the application of the method for hyperbolic and parabolic systems of partial differential equations, considering the wave and heat equations as particular examples. Copyright © 2016 John Wiley & Sons, Ltd. Y1 - 2016 U6 - https://doi.org/10.1002/oca.2238 VL - 37 IS - 6 SP - 1329 EP - 1354 ER - TY - JOUR A1 - Hante, Falk A1 - Leugering, Günter A1 - Martin, Alexander A1 - Schewe, Lars A1 - Schmidt, Martin T1 - Challenges in optimal control problems for gas and fluid flow in networks of pipes and canals: From modeling to industrial applications N2 - 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. KW - Networks KW - pipes KW - optimal control KW - Euler and St. Venant equations KW - hierarchy of models Y1 - 2016 ER - TY - JOUR A1 - Schewe, Lars A1 - Schmidt, Martin T1 - Computing Feasible Points for Binary MINLPs with MPECs JF - Mathematical Programming Computation N2 - 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. KW - Mixed-Integer Nonlinear Optimization KW - MINLP KW - MPEC KW - Complementarity Constraints KW - Primal Heuristic Y1 - 2016 IS - 11(1) SP - 95 EP - 118 ER - TY - JOUR A1 - Gugat, Martin A1 - Hante, Falk T1 - Lipschitz Continuity of the Value Function in Mixed-Integer Optimal Control Problems JF - Mathematics of Control, Signals, and Systems Y1 - 2017 U6 - https://doi.org/10.1007/s00498-016-0183-4 VL - 29 IS - 1 ER - TY - INPR A1 - Gugat, Martin A1 - Perrollaz, Vincent A1 - Rosier, Lionel T1 - Boundary stabilization of quasilinear hyperbolic systems of balance laws: Exponential decay for small source terms N2 - We investigate the long-time behavior of solutions of quasilinear hyperbolic systems with transparent boundary conditions when small source terms are incorporated in the system. Even if the finite-time stability of the system is not preserved, it is shown here that an exponential convergence towards the steady state still holds with a decay rate which is proportional to the logarithm of the amplitude of the source term. The result is stated for a system with dynamical boundary conditions in order to deal with initial data that are free of any compatibility condition. Y1 - 2017 ER - TY - JOUR A1 - Leugering, Günter A1 - Martin, Alexander A1 - Schmidt, Martin A1 - Sirvent, Mathias T1 - Nonoverlapping Domain Decomposition for Optimal Control Problems governed by Semilinear Models for Gas Flow in Networks N2 - We consider optimal control problems for gas flow in pipeline networks. The equations of motion are taken to be represented by a first-order system of hyperbolic semilinear equations derived from the fully nonlinear isothermal Euler gas equations. We formulate an optimal control problem on a network and introduce a tailored time discretization thereof. In order to further reduce the complexity, we consider an instantaneous control strategy. The main part of the paper is concerned with a nonoverlapping domain decomposition of the optimal control problem on the graph into local problems on smaller sub-graphs - ultimately on single edges. We prove convergence of the domain decomposition method on networks and study the wellposedness of the corresponding time-discrete optimal control problems. The point of the paper is that we establish virtual control problems on the decomposed subgraphs such that the corresponding optimality systems are in fact equal to the systems obtained via the domain decomposition of the entire optimality system. KW - Optimal control, Gas networks, Euler's equation, Semilinear PDE, Nonoverlapping domain decomposition Y1 - 2017 VL - 46 IS - 3 SP - 191 EP - 225 PB - Control and Cybernetics ER - TY - JOUR A1 - Gugat, Martin ED - Khalique, Chaudry T1 - Exact Boundary Controllability for Free Traffic Flow with Lipschitz Continuous State JF - Mathematical Problems in Engineering N2 - We consider traffic flow governed by the LWR model. We show that a Lipschitz continuous initial density with free-flow and sufficiently small Lipschitz constant can be controlled exactly to an arbitrary constant free-flow density in finite time by a piecewise linear boundary control function that controls the density at the inflow boundary if the outflow boundary is absorbing. Moreover, this can be done in such a way that the generated state is Lipschitz continuous. Since the target states need not be close to the initial state, our result is a global exact controllability result. The Lipschitz constant of the generated state can be made arbitrarily small if the Lipschitz constant of the initial density is sufficiently small and the control time is sufficiently long. This is motivated by the idea that finite or even small Lipschitz constants are desirable in traffic flow since they might help to decrease the speed variation and lead to safer traffic. KW - Traffic flow Y1 - 2016 U6 - https://doi.org/doi:10.1155/2016/2743251 Creative Commons Attribution License VL - 2016 IS - 2016 SP - 1 EP - 11 ER - TY - JOUR A1 - Gugat, Martin ED - Bloch, Anthony M. T1 - Exponential Stabilization of the Wave Equation by Dirichlet Integral Feedback JF - SIAM Journal on Control and Optimization (SICON) N2 - We consider the problem of boundary feedback stabilization of a vibrating string that is fixed at one end and with control action at the other end. In contrast to previous studies that have required L 2-regularity for the initial position and H −1-regularity for the initial velocity, in this paper we allow for initial positions with L 1-regularity and initial velocities in W −1,1 on the space interval. It is well known that for a certain feedback parameter, for sufficiently regular initial states the classical energy of the closed-loop system with Neumann velocity feedback is controlled to zero after a finite time that is equal to the minimal time where exact controllability holds. In this paper, we present a Dirichlet boundary feedback that yields a well-defined closed-loop system in the (L 1 , W −1,1) framework and also has this property. Moreover, for all positive feedback parameters our feedback law leads to exponential decay of a suitably defined L 1-energy. For more regular initial states with (L 2 , H −1) regularity, the proposed feedback law leads to exponential decay of an energy that corresponds to this framework. If the initial states are even more regular with H 1-regularity of the initial position and L 2-regularity of the initial velocity, our feedback law also leads to exponential decay of the classical energy. KW - exponential stability KW - Dirichlet boundary control KW - energy decay KW - exact control Y1 - 2016 U6 - https://doi.org/DOI: 10.1137/140977023 VL - 53 IS - 1 SP - 526 EP - 546 ER - TY - JOUR A1 - Cacchiani, Valentina A1 - Jünger, Michael A1 - Liers, Frauke A1 - Lodi, Andrea A1 - Schmidt, Daniel T1 - Single-commodity robust network design with finite and Hose demand sets JF - Mathematical Programming N2 - We study a single-commodity robust network design problem (sRND) defined on an undirected graph. Our goal is to determine minimum cost capacities such that any traffic demand from a given uncertainty set can be satisfied by a feasible single-commodity flow. We consider two ways of representing the uncertainty set, either as a finite list of scenarios or as a polytope. We propose a branch-and-cut algorithm to derive optimal solutions to sRND, built on a capacity-based integer linear programming formulation. It is strengthened with valid inequalities derived as {0, 1/2}-Chvátal–Gomory cuts. Since the formulation contains exponentially many constraints, we provide practical separation algorithms. Extensive computational experiments show that our approach is effective, in comparison to existing approaches from the literature as well as to solving a flow based formulation by a general purpose solver. KW - Robust network design KW - Branch-and-cut KW - Cut-set inequalities KW - Polyhedral demand uncertainty KW - Separation under uncertainty Y1 - 2016 U6 - https://doi.org/10.1007/s10107-016-0991-9 VL - 157 IS - 1 SP - 297 EP - 342 ER - TY - JOUR A1 - Gottschalk, Corinna A1 - Koster, Arie M.C.A. A1 - Liers, Frauke A1 - Peis, Britta A1 - Schmand, Daniel A1 - Wierz, Andreas T1 - Robust Flows over Time: Models and Complexity Results N2 - We study dynamic network flows with uncertain input data under a robust optimization perspective. In the dynamic maximum flow problem, the goal is to maximize the flow reaching the sink within a given time horizon T, while flow requires a certain travel time to traverse an arc. In our setting, we account for uncertain travel times of flow. We investigate maximum flows over time under the assumption that at most Γ travel times may be prolonged simultaneously due to delay. We develop and study a mathematical model for this problem. As the dynamic robust flow problem generalizes the static version, it is NP-hard to compute an optimal flow. However, our dynamic version is considerably more complex than the static version. We show that it is NP-hard to verify feasibility of a given candidate solution. Furthermore, we investigate temporally repeated flows and show that in contrast to the non-robust case (i.e., without uncertainties) they no longer provide optimal solutions for the robust problem, but rather yield a worst case optimality gap of at least T. We finally show that for infinite delays, the optimality gap is at most O(k log T), where k is a newly introduced instance characteristic. The results obtained in this paper yield a first step towards understanding robust dynamic flow problems with uncertain travel times. KW - Dynamic Network Flows KW - Uncertain Travel Times KW - Dynamic Robust Flow Y1 - 2017 U6 - https://doi.org/10.1007/s10107-017-1170-3 ER - TY - JOUR A1 - Schweiger, Jonas A1 - Liers, Frauke T1 - A Decomposition Approach for Optimum Gas Network Extension with a Finite Set of Demand Scenarios N2 - Today's gas markets demand more flexibility from the network operators which in turn have to invest into their network infrastructure. As these investments are very cost-intensive and long-living, network extensions should not only focus on a single bottleneck scenario, but should increase the flexibility to fulfill different demand scenarios. In this work, we formulate a model for the network extension problem for multiple demand scenarios and propose a scenario decomposition in order to solve the arising challenging optimization tasks. In fact, euch subproblem consists of a mixed-integer nonlinear optimization problem (MINLP). Valid bounds are derived even without solving the subproblems to optimality. Furthermore, we develop heuristics that prove capable of improving the initial solutions substantially. Results of computational experiments on realistic network topologies are presented. It turns out that our method is able to solve these challenging instances to optimality within a reasonable amount of time. Y1 - 2016 ER - TY - JOUR A1 - Geißler, Björn A1 - Morsi, Antonio A1 - Schewe, Lars A1 - Schmidt, Martin T1 - Solving Highly Detailed Gas Transport MINLPs: Block Separability and Penalty Alternating Direction Methods JF - INFORMS Journal on Computing N2 - 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. Y1 - 2016 IS - 30(2) SP - 309 EP - 323 ER - TY - JOUR A1 - Geißler, Björn A1 - Morsi, Antonio A1 - Schewe, Lars A1 - Schmidt, Martin T1 - Penalty Alternating Direction Methods for Mixed-Integer Optimization: A New View on Feasibility Pumps JF - SIAM Journal on Optimization N2 - Feasibility pumps are highly effective primal heuristics for mixed-integer linear and nonlinear optimization. However, despite their success in practice there are only few works considering their theoretical properties. We show that feasibility pumps can be seen as alternating direction methods applied to special reformulations of the original problem, inheriting the convergence theory of these methods. Moreover, we propose a novel penalty framework that encompasses this alternating direction method, which allows us to refrain from random perturbations that are applied in standard versions of feasibility pumps in case of failure. We present a convergence theory for the new penalty based alternating direction method and compare the new variant of the feasibility pump with existing versions in an extensive numerical study for mixed-integer linear and nonlinear problems. Y1 - 2017 U6 - https://doi.org/10.1137/16M1069687 VL - 27 IS - 3 SP - 1611 EP - 1636 ER - TY - JOUR A1 - Grimm, Veronika A1 - Martin, Alexander A1 - Schmidt, Martin A1 - Weibelzahl, Martin A1 - Zöttl, Gregor T1 - Transmission and generation investment in electricity markets: The effects of market splitting and network fee regimes JF - European Journal of Operational Research N2 - 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. KW - Electricity market modeling KW - Mixed-integer nonlinear optimization KW - Multilevel programming KW - Network expansion KW - Transmission management Y1 - 2016 U6 - https://doi.org/10.1016/j.ejor.2016.03.044 VL - 254 IS - 2 SP - 493 EP - 509 ER -