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 - 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 - 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 - 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 - 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 - 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 - Heitsch, Holger A1 - Henrion, René T1 - An enumerative formula for the spherical cap discrepancy N2 - The spherical cap discrepancy is a widely used measure for how uniformly a sample of points on the sphere is distributed. Being hard to compute, this discrepancy measure is typically replaced by some lower or upper estimates when designing optimal sampling schemes for the uniform distribution on the sphere. In this paper, we provide a fully explicit, easy to implement enumerative formula for the spherical cap discrepancy. Not surprisingly, this formula is of combinatorial nature and, thus, its application is limited to spheres of small dimension and moderate sample sizes. Nonetheless, it may serve as a useful calibrating tool for testing the efficiency of sampling schemes and its explicit character might be useful also to establish necessary optimality conditions when minimizing the discrepancy with respect to a sample of given size. KW - spherical cap discrepancy KW - uniform distribution on sphere KW - optimality conditions Y1 - 2019 U6 - https://doi.org/10.1016/j.cam.2021.113409 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 - 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 - 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 - Grimm, Veronika A1 - Rückel, Bastian A1 - Sölch, Christian A1 - Zöttl, Gregor T1 - Regionally differentiated network fees to affect incentives for generation investment JF - Energy N2 - In this paper we propose an equilibrium model that allows to analyze subsidization schemes to affect locational choices for generation investment in electricity markets. Our framework takes into account generation investment decided by private investors and redispatch as well as network expansion decided by a regulated transmission system operator. In order to take into account the different objectives and decision variables of those agents, our approach uses a bi-level structure. We focus on the case of regionally differentiated network fees which have to be paid by generators (a so called g-component). The resulting investment and production decisions are compared to the outcome of an equilibrium model in the absence of such regionally differentiated investment incentives and to an overall optimal (first-best) benchmark. To illustrate possible economic effects, we calibrate our framework with data from the German electricity market. Our results reveal that while regionally differentiated network fees do have a significant impact on locational choice of generation capacities, we do not find significant effects on either welfare or network expansion. KW - Electricity Markets, Network Expansion, Generation Expansion, Investment Incentives, Computational Equilibrium Models Y1 - 2019 U6 - https://doi.org/10.1016/j.energy.2019.04.035 IS - 177 SP - 487 EP - 502 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 - 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 - 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 - JOUR A1 - Ulbrich, Stefan A1 - Schmitt, Johann Michael A1 - Schäfer Aguilar, Paloma A1 - Moos, Michael T1 - On the numerical discretization of optimal control problems for conservation laws N2 - We analyze the convergence of discretization schemes for the adjoint equation arising in the adjoint-based derivative computation for optimal control problems governed by entropy solutions of conservation laws. The difficulties arise from the fact that the correct adjoint state is the reversible solution of a transport equation with discontinuous coefficient and discontinuous end data. We derive the discrete adjoint scheme for monotone difference schemes in conservation form. It is known that convergence of the discrete adjoint can only be expected if the numerical scheme has viscosity of order O(h^\alpha) with appropriate 0 < \alpha < 1, which leads to quite viscous shock profiles. We show that by a slight modification of the end data of the discrete adjoint scheme convergence to the correct reversible solution can be obtained also for numerical schemes with viscosity of order O(h) and with sharp shock resolution. The theoretical findings are confirmed by numerical results. Y1 - 2019 ER - TY - INPR A1 - Klimm, Max A1 - Pfetsch, Marc E. A1 - Raber, Rico A1 - Skutella, Martin T1 - Packing under Convex Quadratic Constraints N2 - We consider a general class of binary packing problems with a convex quadratic knapsack constraint. We prove that these problems are APX-hard to approximate and present constant-factor approximation algorithms based upon three different algorithmic techniques: (1) a rounding technique tailored to a convex relaxation in conjunction with a non-convex relaxation whose approximation ratio equals the golden ratio; (2) a greedy strategy; (3) a randomized rounding method leading to an approximation algorithm for the more general case with multiple convex quadratic constraints. We further show that a combination of the first two strategies can be used to yield a monotone algorithm leading to a strategyproof mechanism for a game-theoretic variant of the problem. Finally, we present a computational study of the empirical approximation of the three algorithms for problem instances arising in the context of real-world gas transport networks. 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 - 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 - 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 - 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 - 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 - Heitsch, Holger T1 - On probabilistic capacity maximization in a stationary gas network JF - Optimization N2 - The question for the capacity of a given gas network, i.e., determining the maximal amount of gas that can be transported by a given network, appears as an essential question that network operators and political administrations are regularly faced with. In that context we present a novel mathematical approach in order to assist gas network operators in managing increasing uncertainty with respect to customers gas nominations and in exposing free network capacities while reliability of transmission and supply is taken into account. The approach is based on the rigorous examination of optimization problems with nonlinear probabilistic constraints. As consequence we deal with solving a problem belonging to the class of probabilistic/robust optimization problems, which can be formulated with some joint probabilistic constraint over an infinite system of random inequalities. We will show that the inequality system can be reduced to a finite one in the situation of considering a tree network topology. A detailed study of the problem of maximizing bookable capacities in a stationary gas network is presented that comes up with an algebraic model involving Kirchhoff's first and second laws. The focus will be on both the theoretical and numerical side. The analytical part consists in introducing and validating a generalized version of the known rank two constraint qualification implying the differentiability of the considered capacity problem. The results are important in order to solve the capacity problem numerically, where function and gradient evaluations of the probabilistic constraints are performed by an approach using spheric-radial decomposition applicable for multivariate Gaussian random variables and more general distributions. KW - stationary gas networks, booked capacities, probabilistic constraints, constraint qualification, spheric-radial decomposition Y1 - 2019 U6 - https://doi.org/10.1080/02331934.2019.1625353 ER - TY - JOUR A1 - Heitsch, Holger A1 - Strogies, Nikolai T1 - Consequences of uncertain friction for the transport of natural gas through passive networks of pipelines JF - Springer N2 - Assuming a pipe-wise constant structure of the friction coefficient in the modeling of natural gas transport through a passive network of pipes via semilinear systems of balance laws with associated linear coupling and boundary conditions, uncertainty in this parameter is quantified by a Markov chain Monte Carlo method. Information on the prior distribution is obtained from practitioners. The results are applied to the problem of validating technical feasibility under random exit demand in gas transport networks. The impact of quantified uncertainty to the probability level of technical feasible exit demand situations is studied by two example networks of small and medium size. The gas transport of the network is modeled by stationary solutions that are steady states of the time dependent semilinear problems. KW - uncertainty quantification, Markov chain Monte Carlo, reliability of gas networks, nomination validation, spheric-radial decomposition Y1 - 2019 SP - 211 EP - 238 ET - Topics in Applied Analysis and Optimisation 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 - Hauschild, Sarah-Alexa A1 - Marheineke, Nicole A1 - Mehrmann, Volker A1 - Mohring, Jan A1 - Badlyan, Arbi Moses A1 - Rein, Markus A1 - Schmidt, Martin T1 - Port-Hamiltonian modeling of district heating networks JF - Progress in Differential Algebraic Equations II (edited by Reis T., Grundel S., and Schöps S). Differential-Algebraic Equations Forum N2 - This paper provides a first contribution to port-Hamiltonian modeling of district heating networks. By introducing a model hierarchy of flow equations on the network, this work aims at a thermodynamically consistent port-Hamiltonian embedding of the partial differential-algebraic systems. We show that a spatially discretized network model describing the advection of the internal energy density with respect to an underlying incompressible stationary Euler-type hydrodynamics can be considered as a parameter-dependent finite-dimensional port-Hamiltonian system. Moreover, we present an infinite-dimensional port-Hamiltonian formulation for a compressible instationary thermodynamic fluid flow in a pipe. Based on these first promising results, we raise open questions and point out research perspectives concerning structure-preserving discretization, model reduction, and optimization. KW - Partial differential equations on networks KW - Port-Hamiltonian model framework KW - Energy-based formulation KW - District heating network KW - Thermodynamic fluid flow 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - Ullmann, Sebastian A1 - Müller, Christopher A1 - Lang, Jens T1 - Stochastic Galerkin Reduced Basis Methods for Parametrized Linear Convection-Diffusion-Reaction Equations N2 - We consider the estimation of parameter-dependent statistics of functional outputs of steady-state convection–diffusion–reaction equations with parametrized random and deterministic inputs in the framework of linear elliptic partial differential equations. For a given value of the deterministic parameter, a stochastic Galerkin finite element (SGFE) method can estimate the statistical moments of interest of a linear output at the cost of solving a single, large, block-structured linear system of equations. We propose a stochastic Galerkin reduced basis (SGRB) method as a means to lower the computational burden when statistical outputs are required for a large number of deterministic parameter queries. Our working assumption is that we have access to the computational resources necessary to set up such a reduced-order model for a spatial-stochastic weak formulation of the parameter-dependent model equations. In this scenario, the complexity of evaluating the SGRB model for a new value of the deterministic parameter only depends on the reduced dimension. To derive an SGRB model, we project the spatial-stochastic weak solution of a parameter-dependent SGFE model onto a reduced basis generated by a proper orthogonal decomposition (POD) of snapshots of SGFE solutions at representative values of the parameter. We propose residual-corrected estimates of the parameter-dependent expectation and variance of linear functional outputs and provide respective computable error bounds.We test the SGRB method numerically for a convection–diffusion–reaction problem, choosing the convective velocity as a deterministic parameter and the parametrized reactivity or diffusivity field as a random input. Compared to a standard reduced basis model embedded in a Monte Carlo sampling procedure, the SGRB model requires a similar number of reduced basis functions to meet a given tolerance requirement. However, only a single run of the SGRB model suffices to estimate a statistical output for a new deterministic parameter value, while the standard reduced basis model must be solved for each Monte Carlo sample. Y1 - 2018 U6 - https://doi.org/doi:10.3390/fluids6080263 VL - Fluids IS - 6 SP - 263 ER - TY - JOUR A1 - Mindt, Pascal A1 - Lang, Jens A1 - Domschke, Pia T1 - Entropy-Preserving Coupling of Hierarchical Gas Models N2 - This paper is concerned with coupling conditions at junctions for transport models which differ in their fidelity to describe transient flow in gas pipelines. It also includes the integration of compressors between two pipes with possibly different models. A hierarchy of three one-dimensional gas transport models is built through the 3 × 3 polytropic Euler equations, the 2 × 2 isentropic Euler equations and a simplified version of it for small velocities. To ensure entropy preservation, we make use of the novel entropy-preserving coupling conditions recently proposed by Lang and Mindt [Netw. Heterog. Media, 13:177-190, 2018] and require the equality of the total enthalpy at the junction and that the specific entropy for pipes with outgoing flow equals the convex combination of all entropies that belong to pipes with incoming flow. We prove the existence and uniqueness of solutions to generalised Riemann problems at a junction in the neighbourhood of constant coupling functions and stationary states which belong to the subsonic region. This provides the basis for the well-posedness of certain Cauchy problems for initial data with sufficiently small total variation. Y1 - 2018 U6 - https://doi.org/doi:10.1137/19M1240034 VL - SIAM Journal on Mathematical Analysis IS - 51 SP - 4754 EP - 4775 ER - TY - JOUR A1 - Schneider, Moritz A1 - Lang, Jens A1 - Weiner, Rüdiger T1 - Super-Convergent Implicit-Explicit Peer Methods with Variable Step Sizes N2 - Dynamical systems with sub-processes evolving on many different time scales are ubiquitous in applications. Their efficient solution is greatly enhanced by automatic time step variation. This paper is concerned with the theory, construction and application of IMEX-Peer methods that are super-convergent for variable step sizes and A-stable in the implicit part. IMEX schemes combine the necessary stability of implicit and low computational costs of explicit methods to efficiently solve systems of ordinary differential equations with both stiff and non-stiff parts included in the source term. To construct super-convergent IMEX-Peer methods which keep their higher order for variable step sizes and exhibit favourable linear stability properties, we derive necessary and sufficient conditions on the nodes and coefficient matrices and apply an extrapolation approach based on already computed stage values. New super-convergent IMEX-Peer methods of order s + 1 for s = 2, 3, 4 stages are given as result of additional order conditions which maintain the super-convergence property independent of step size changes. Numerical experiments and a comparison to other super-convergent IMEX-Peer methods show the potential of the new methods when applied with local error control. Y1 - 2019 U6 - https://doi.org/doi:10.1016/j.cam.2019.112501 VL - J. Comput. Appl. Math. IS - 387 SP - 112501 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 - 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 - 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 - 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 - Gonzalez Grandon, Tatiana A1 - Henrion, Rene A1 - Perez-Aros, Pedro T1 - Dynamic probabilistic constraints under continuous random distributions N2 - In this paper we address novel results on the theoretical structural analysis of dynamic joint probabilistic constraints under continuous random variables. This dynamic probabilistic function is important when decisions are time-dependent and when the modeler can react on past observations. We first study the continuity of dynamic probabilistic constraints and provide strong and weak semi-continuous results depending on whether the policies are supposed to be in the L^p or W^{1,p} spaces. Moreover, we prove the non-convexity of the feasible set of decisions induced by a dynamic probability function in the L^p space. Lastly, for a simple two-stage model, verifiable conditions for Lipschitz continuity and differentiability of this probability function are derived and endowed with explicit derivative formulae. Y1 - 2019 U6 - https://doi.org/10.1007/s10107-020-01593-z 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - Hajian, Soheil A1 - Hintermüller, Michael A1 - Schillings, Claudia A1 - Strogies, Nikolai T1 - A Bayesian approach to parameter identification in gas networks JF - Control and Cybernetics N2 - The inverse problem of identifying the friction coefficient in an isothermal semilinear Euler system is considered. Adopting a Bayesian approach, the goal is to identify the distribution of the quantity of interest based on a finite number of noisy measurements of the pressure at the boundaries of the domain. First well-posedness of the underlying non-linear PDE system is shown using semigroup theory, and then Lipschitz continuity of the solution operator with respect to the friction coefficient is established. Based on the Lipschitz property, well-posedness of the resulting Bayesian inverse problem for the identification of the friction coefficient is inferred. Numerical tests for scalar and distributed parameters are performed to validate the theoretical results. Y1 - 2018 VL - 48 SP - 377 EP - 402 ER - TY - JOUR A1 - Rüffler, Fabian A1 - Mehrmann, Volker A1 - Hante, Falk T1 - Optimal Model Switching for Gas Flow in Pipe Networks N2 - We consider model adaptivity for gas flow in pipeline networks. For each instant in time and for each pipe in the network a model for the gas flow is to be selected from a hierarchy of models in order to maximize a performance index that balances model accuracy and computational cost for a simulation of the entire network. This combinatorial problem involving partial differential equations is posed as an optimal switching control problem for abstract semilinear evolutions. We provide a theoretical and numerical framework for solving this problem using a two stage gradient descent approach based on switching time and mode insertion gradients. A numerical study demonstrates the practicability of the approach. Y1 - 2018 ER - TY - JOUR A1 - Gugat, Martin A1 - Steffensen, Sonja ED - Trelat, Emmanuel T1 - DYNAMIC BOUNDARY CONTROL GAMES WITH NETWORKS OF STRINGS JF - ESAIM: Control, Optimisation and Calculus of Variations (ESAIM: COCV) N2 - Consider a star-shaped network of strings. Each string is governed by the wave equation. At each boundary node of the network there is a player that performs Dirichlet boundary control action and in this way influences the system state. At the central node, the states are coupled by algebraic conditions in such a way that the energy is conserved. We consider the corresponding antagonistic game where each player minimizes a certain quadratic objective function that is given by the sum of a control cost and a tracking term for the final state. We prove that under suitable assumptions a unique Nash equilibrium exists and give an explicit representation of the equilibrium strategies. Y1 - 2018 U6 - https://doi.org/10.1051/cocv/2017082 ER - TY - JOUR A1 - Farshbaf-Shaker, M. Hassan A1 - Henrion, Rene A1 - Hömberg, Dietmar T1 - Properties of Chance Constraints in Infinite Dimensions with an Application to PDE Constrained Optimization JF - Set-Valued and Variational Analysis N2 - Chance constraints represent a popular tool for finding decisions that enforce the satisfaction of random inequality systems in terms of probability. They are widely used in optimization problems subject to uncertain parameters as they arise in many engineering applications. Most structural results of chance constraints (e.g., closedness, convexity, Lipschitz continuity, differentiability etc.) have been formulated in finite dimensions. The aim of this paper is to generalize some of these well-known semi-continuity and convexity properties as well as a stability result to an infinite dimensional setting. The abstract results are applied to a simple PDE constrained control problem subject to (uniform) state chance constraints. KW - Chance constraints KW - Probabilistic constraints KW - PDE constrained optimization Y1 - 2018 U6 - https://doi.org/doi:10.1007/s11228-017-0452-5 ER - TY - JOUR A1 - Hantoute, Abderrahim A1 - Henrion, Rene A1 - Perez-Aros, Pedro T1 - Subdifferential characterization of probability functions under Gaussian distribution N2 - Probability functions figure prominently in optimization problems of engineering. They may be nonsmooth even if all input data are smooth. This fact motivates the consideration of subdifferentials for such typically just continuous functions. The aim of this paper is to provide subdifferential formulae of such functions in the case of Gaussian distributions for possibly infinite-dimensional decision variables and nonsmooth (locally Lipschitzian) input data. These formulae are based on the spheric-radial decomposition of Gaussian random vectors on the one hand and on a cone of directions of moderate growth on the other. By successively adding additional hypotheses, conditions are satisfied under which the probability function is locally Lipschitzian or even differentiable. Y1 - 2018 U6 - https://doi.org/10.1007/s10107-018-1237-9 ER - TY - JOUR A1 - Hintermüller, Michael A1 - Strogies, Nikolai T1 - Identification of the friction function in a semilinear system for gas transport through a network JF - Optimization Methods and Software N2 - An identification problem for the friction parameter in a semilinear system of balance laws, describing the transport of gas through a passive network of pipelines, is considered. The existence of broad solutions to the state system is proven and sensitivity results for the corresponding solution operator are obtained. The existence of solutions to the output least squares formulation of the identification problem, based on noisy measurements over time at fixed spatial positions is established. Finally, numerical experiments validate the theoretical findings. Y1 - 2017 VL - 35 SP - 576 EP - 617 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 - Schneider, Moritz A1 - Lang, Jens A1 - Hundsdorfer, Willem T1 - Extrapolation-Based Super-Convergent Implicit-Explicit Peer Methods with A-stable Implicit Part N2 - In this paper we extend the implicit-explicit (IMEX) methods of Peer type recently developed in [Lang, Hundsdorfer, J. Comp. Phys., 337:203–215, 2017] to a broader class of two-step methods that allow the construction of super- convergent IMEX-Peer methods with A-stable implicit part. IMEX schemes combine the necessary stability of implicit and low computational costs of ex- plicit methods to efficiently solve systems of ordinary differential equations with both stiff and non-stiff parts included in the source term. To construct super- convergent IMEX-Peer methods with favourable stability properties, we derive necessary and sufficient conditions on the coefficient matrices and apply an extrapolation approach based on already computed stage values. Optimised super-convergent IMEX-Peer methods of order s + 1 for s = 2, 3, 4 stages are given as result of a search algorithm carefully designed to balance the size of the stability regions and the extrapolation errors. Numerical experiments and a comparison to other IMEX-Peer methods are included. Y1 - 2017 U6 - https://doi.org/10.1016/j.jcp.2018.04.006 VL - J. Comput. Physics IS - Vol. 367 SP - 121 EP - 133 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 - 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 - 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 - INPR A1 - Schmitt, Johann Michael A1 - Ulbrich, Stefan T1 - Optimal Boundary Control of Hyperbolic Balance Laws with State Constraints N2 - In this paper we analyze the optimal control of initial-boundary value problems for entropy solutions of scalar hyperbolic balance laws with pointwise state constraints. Hereby, we suppose that the initial and the boundary data switch between different C¹-functions at certain switching points, where the C¹ -functions and the switching points are considered as the control. For a class of cost functionals, we prove first order necessary optimality conditions for the corresponding optimal control problem with state constraints. Furthermore, we use a Moreau-Yosida type regularization to approximate the optimal control problem with state constraints. We derive optimality conditions for the regularized problems and finally prove convergence to the solution of the optimal control problem with state constraints. Y1 - 2021 ER - TY - RPRT A1 - Hiller, Benjamin A1 - Walther, Tom T1 - Improving branching for disjunctive polyhedral models using approximate convex decompositions N2 - Disjunctive sets arise in a variety of optimization models and much esearch has been devoted to obtain strong relaxations for them. This paper focuses on the evaluation of the relaxation during the branch-and-bound search process. We argue that the branching possibilities (\ie binary variables) of the usual formulations are unsuitable to obtain strong bounds early in the search process as they do not capture the overall shape of the the entire disjunctive set. To analyze and exploit the shape of the disjunctive set we propose to compute a hierarchy of approximate convex decompositions and show how to extend the known formulations to obtain improved branching behavior. KW - disjunctive set KW - convex hull reformulation KW - branch-and-bound KW - approximate convex decomposition Y1 - 2017 ER - TY - INPR A1 - Egger, Herbert A1 - Kugler, Thomas A1 - Wollner, Winnifried T1 - Numerical optimal control of instationary gas transport with control and state constraints N2 - We consider the optimal control of a nonlinear hyperbolic system of balance laws on a one-dimensional network which arises in the context of gas transport in pipeline systems. State constraints, which are required for the safe operation of the system, are incorporated by a barrier method. We discuss the well-posedness of the governing system of partial differential-algebraic equations and investigate the existence of minimizers. For the numerical solution, we then consider the approximation of the state equation by mixed finite elements in space and a particular linear implicit time integration scheme that can be interpreted as a discontinuous Galerkin approximation. We establish well- posedness of this discretization scheme and prove the existence of minimizers for the corresponding discretized optimal control problem and discuss its numerical solution by a projected Gauß-Newton method. The efficient realization of the Jacobian and Hessian of the quadratic approximations that have to be minimized in every iteration of the Gauß-Newton method can be obtained via the solution of discretized sensitivity and adjoint equations. These are obtained by formal differentiation and transposition of the Galerkin methods employed for the discretization of the state equations. All approximations obtained after discretization can thus be interpreted as functions on the continuous level and, since the functional analytic setting is not changed by the Galerkin discretization, we observe mesh independence of the resulting fully discrete methods. For illustration of our theoretical results and to demonstrate the efficiency of the proposed method, we present numerical results for two test problems that model typical situations that may arise in the daily operation of gas networks. Y1 - 2017 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 - INPR A1 - Egger, Herbert A1 - Böttcher, Anke T1 - Energy stable discretization of Allen-Cahn type problems modeling the motion of phase boundaries N2 - We study the systematic numerical approximation of a class of Allen-Cahn type problems modeling the motion of phase interfaces. The common feature of these models is an underlying gradient flow structure which gives rise to a decay of an associated energy functional along solution trajectories. We first study the discretization in space by a conforming Galerkin approximation of a variational principle which characterizes smooth solutions of the problem. Well-posedness of the resulting semi-discretization is established and the energy decay along discrete solution trajectories is proven. A problem adapted implicit time-stepping scheme is then proposed and we establish its well-posed and decay of the free energy for the fully discrete scheme. Some details about the numerical realization by finite elements are discussed, in particular the iterative solution of the nonlinear problems arising in every time-step. The theoretical results are illustrated by numerical tests which also provide further evidence for asymptotic expansions of the interface velocities derived by Alber et al. Y1 - 2017 ER - TY - INPR A1 - Egger, Herbert A1 - Radu, Bogdan T1 - Super-convergence and post-processing for mixed finite element approximations of the wave equation N2 - We consider the numerical approximation of acoustic wave propagation problems by mixed BDM(k+1)-P(k) finite elements on unstructured meshes. Optimal convergence of the discrete velocity and super-convergence of the pressure by one order are established. Based on these results, we propose a post-processing strategy that allows us to construct an improved pressure approximation from the numerical solution. Corresponding results are well-known for mixed finite element approximations of elliptic problems and we extend these analyses here to the hyperbolic problem under consideration. We also consider the subsequent time discretization by the Crank-Nicolson method and show that the analysis and the post-processing strategy can be generalized to the fully discrete schemes. Our proofs do not rely on duality arguments or inverse inequalities and the results therefore apply also for non-convex domains and non-uniform meshes. Y1 - 2017 ER - TY - INPR A1 - Egger, Herbert T1 - A mixed variational discretization for non-isothermal compressible flow in pipelines N2 - We consider the non-isothermal flow of a compressible fluid through pipes. Starting from the full set of Euler equations, we propose a variational characterization of solutions that encodes the conservation of mass, energy, and entropy in a very direct manner. This variational principle is suitable for a conforming Galerkin approximation in space which automatically inherits the basic physical conservation laws. Three different spaces are used for approximation of density, mass flux, and temperature, and we consider a mixed finite element method as one possible choice of suitable approximation spaces. We also investigate the subsequent discretization in time by a problem adapted implicit time stepping scheme for which exact conservation of mass as well as a slight dissipation of energy and increase of entropy are proven which are due to the numerical dissipation of the implicit time discretization. The main arguments of our analysis are rather general and allow us to extend the approach with minor modification to more general boundary conditions and flow models taking into account friction, viscosity, heat conduction, and heat exchange with the surrounding medium. Y1 - 2017 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 - Egger, Herbert A1 - Pietschmann, Jan-Frederik A1 - Schlottbom, Matthias T1 - On the uniqueness of nonlinear diffusion coefficients in the presence of lower order terms N2 - We consider the identification of nonlinear diffusion coefficients of the form a(t,u) or a(u) in quasi-linear parabolic and elliptic equations. Uniqueness for this inverse problem is established under very general assumptions using partial knowledge of the Dirichlet-to-Neumann map. The proof of our main result relies on the construction of a series of appropriate Dirichlet data and test functions with a particular singular behavior at the boundary. This allows us to localize the analysis and to separate the principal part of the equation from the remaining terms. We therefore do not require specific knowledge of lower order terms or initial data which allows to apply our results to a variety of applications. This is illustrated by discussing some typical examples in detail. Y1 - 2017 ER - TY - INPR A1 - Egger, Herbert T1 - Enhancement of flow measurements using fluid dynamic constraints N2 - Novel experimental modalities acquire spatially resolved velocity measurements for steady state and transient flows which are of interest for engineering and biological applications. One of the drawbacks of such high resolution velocity data is their susceptibility to measurement errors. In this paper, we propose a novel filtering strategy that allows enhancement of noisy measurements to obtain reconstruction of smooth divergence free velocity and corresponding pressure fields, which together approximately comply to a prescribed flow model. The main step in our approach consists of the appropriate use of the velocity measurements in the design of a linearized flow model which can be shown to be well-posed and consistent with the true velocity and pressure fields up to measurement and modeling errors. The reconstruction procedure is formulated as a linear quadratic optimal control problem and the resulting filter has analyzable smoothing and approximation properties. We also discuss briefly the discretization of our approach by finite element methods and comment on the efficient solution of the linear optimality system by iterative solvers. The capability of the proposed method to significantly reduce data noise is demonstrated by numerical tests in which we also compare to other methods like smoothing and solenoidal filtering. Y1 - 2017 ER - TY - INPR A1 - Egger, Herbert A1 - Kugler, Thomas T1 - An asymptotic preserving mixed finite element method for wave propagation in pipelines N2 - We consider a parameter dependent family of damped hyperbolic equations with interesting limit behavior: the system approaches steady states exponentially fast and for parameter to zero the solutions converge to that of a parabolic limit problem. We establish sharp estimates and elaborate their dependence on the model parameters. For the numerical approximation we then consider a mixed finite element method in space together with a Runge-Kutta method in time. Due to the variational and dissipative nature of this approximation, the limit behavior of the infinite dimensional level is inherited almost automatically by the discrete problems. The resulting numerical method thus is asymptotic preserving in the parabolic limit and uniformly exponentially stable. These results are further shown to be independent of the discretization parameters. Numerical tests are presented for a simple model problem which illustrate that the derived estimates are sharp in general. Y1 - 2017 ER - TY - INPR A1 - Egger, Herbert T1 - A class of Galerkin schemes for time-dependent radiative transfer N2 - The numerical solution of time-dependent radiative transfer problems is challenging, both, due to the high dimension as well as the anisotropic structure of the underlying integro-partial differential equation. In this paper we propose a general framework for designing numerical methods for time-dependent radiative transfer based on a Galerkin discretization in space and angle combined with appropriate time stepping schemes. This allows us to systematically incorporate boundary conditions and to preserve basic properties like exponential stability and decay to equilibrium also on the discrete level. We present the basic a-priori error analysis and provide abstract error estimates that cover a wide class of methods. The starting point for our considerations is to rewrite the radiative transfer problem as a system of evolution equations which has a similar structure like first order hyperbolic systems in acoustics or electrodynamics. This analogy allows us to generalize the main arguments of the numerical analysis for such applications to the radiative transfer problem under investigation. We also discuss a particular discretization scheme based on a truncated spherical harmonic expansion in angle, a finite element discretization in space, and the implicit Euler method in time. The performance of the resulting mixed PN-finite element time stepping scheme is demonstrated by computational results. Y1 - 2017 ER - TY - INPR A1 - Beckers, Susanne A1 - Behrens, Jörn A1 - Wollner, Winnifried T1 - Duality Based Error Estimation in the Presence of Discontinuities N2 - Goal-oriented mesh adaptation, in particular using the dual-weighted residual (DWR) method, is known in many cases to produce very efficient meshes. For obtaining such meshes the (numerical) solution of an adjoint problem is needed to weight the residuals appropriately with respect to their relevance for the overall error. For hyperbolic problems already the weak primal problem requires in general an additional entropy condition to assert uniqueness of solutions; this difficulty is also reflected when considering adjoints to hyperbolic problems involving discontinuities where again an additional requirement (reversibility) is needed to select appropriate solutions. Within this article, an approach to the DWR method for hyperbolic problems based on an artificial viscosity approximation is proposed. It is discussed why the proposed method provides a well-posed dual problem, while a direct, formal, application of the dual problem does not. Moreover, we will discuss a further, novel, approach in which the forward problem need not be modified, thus allowing for an unchanged forward solution. The latter procedure introduces an additional residual term in the error estimation, accounting for the inconsistency between primal and dual problem. Finally, the effectivity of the extended error estimator, assessing the global error by a suitable functional of interest, is tested numerically; and the advantage over a formal estimator approach is demonstrated. KW - dual weighted residual, hyperbolic problems, discontinuous Galerkin, artificial viscosity Y1 - 2017 ER - TY - JOUR A1 - Hanke, Michael A1 - März, Roswitha A1 - Tischendorf, Caren T1 - Least-Squares Collocation for Higher-Index Linear Differential-Algebraic Equations: Estimating the Instability Threshold JF - Math. Comp. N2 - Differential-algebraic equations with higher index give rise to essentially ill-posed problems. The least-squares collocation by discretizing the pre-image space is not much more computationally expensive than standard collocation methods used in the numerical solution of ordinary differential equations and index-1 differential-algebraic equations. This approach has displayed excellent convergence properties in numerical experiments, however, theoretically, till now convergence could be established merely for regular linear differential-algebraic equations with constant coefficients. We present now an estimate of the instability threshold which serves as the basic key for proving convergence for general regular linear DAEs. KW - differential-algebraic equation, higher index, essentially ill-posed problem, collocation, boundary value problem, initial value problem Y1 - 2017 U6 - https://doi.org/10.1090/mcom/3393 VL - 88 SP - 1647 EP - 1683 ER - TY - RPRT A1 - Walther, Tom A1 - Hiller, Benjamin T1 - Modelling compressor stations in gas networks N2 - Gas networks are an important application area for optimization. When considering long-range transmission, compressor stations play a crucial role in these applications. The purpose of this report is to collect and systematize the models used for compressor stations in the literature. The emphasis is on recent work on simple yet accurate polyhedral models that may replace more simplified traditional models without increasing model complexity. The report also describes an extension of the compressor station data available in GasLib (http://gaslib.zib.de/) with the parameters of these models. KW - gas compressor modeling KW - compressor station KW - polyhedral model KW - MINLP model KW - nonlinear model Y1 - 2017 ER - TY - INPR A1 - Walther, Tom A1 - Hiller, Benjamin A1 - Saitenmacher, René T1 - Polyhedral 3D Models for compressors in gas networks N2 - Compressor machines are crucial elements in a gas transmission network, required to compensate for the pressure loss caused by friction in the pipes. Modelling all physical and technical details of a compressor machine involves a large amount of nonlinearity, which makes it hard to use such models in the optimization of large-scale gas networks. In this paper, we are going to describe a modelling approach for the operating range of a compressor machine, starting from a physical reference model and resulting in a polyhedral representation in the 3D space of mass flow throughput as well as in- and outlet pressure. KW - gas compressor modeling KW - polyhedral model KW - accuracy KW - z factor 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 - INPR A1 - Hiller, Benjamin A1 - Saitenmacher, René A1 - Walther, Tom T1 - Analysis of operating modes of complex compressor stations N2 - We consider the modeling of operation modes for complex compressor stations (i.e., ones with several in- or outlets) in gas networks. In particular, we propose a refined model that allows to precompute tighter relaxations for each operation mode. These relaxations may be used to strengthen the compressor station submodels in gas network optimization problems. We provide a procedure to obtain the refined model from the input data for the original model. This procedure is based on a nontrivial reduction of the graph representing the gas flow through the compressor station in an operation mode. KW - model reformulation KW - bound tightening KW - multi-way compressor station KW - flow patterns Y1 - 2017 ER - TY - JOUR A1 - Lang, Jens A1 - Mindt, Pascal T1 - Entropy-Preserving Coupling Conditions for One-dimensional Euler Systems at Junctions N2 - This paper is concerned with a set of novel coupling conditions for the 3x3 one-dimensional Euler system with source terms at a junction of pipes with possibly different cross-sectional areas. Beside conservation of mass, we require the equality of the total enthalpy at the junction and that the specific entropy for pipes with outgoing flow equals the convex combination of all entropies that belong to pipes with incoming flow. These conditions ensure energy as well as entropy conservation at the junction. We prove the existence and uniqueness of solutions to the generalised Riemann problem at a junction in the neighbourhood of constant stationary states which belong to the subsonic region. This provides the basis for the well-posedness of the homogeneous and inhomogeneous Cauchy problems for initial data with sufficiently small total variation. Y1 - 2017 U6 - https://doi.org/10.3934/nhm.2018008 VL - Networks and Heterogeneous Media IS - Vol. 13 SP - 177 EP - 190 ER - TY - GEN A1 - Lang, Jens A1 - Leugering, Günter A1 - Martin, Alexander A1 - Tischendorf, Caren T1 - Gasnetzwerke: Mathematische Modellierung, Simulation und Optimierung N2 - Im Mai 2014 wurde seitens der DFG der Transregio 154 Mathematische Modellierung, Simulation und Optimierung am Beispiel von Gasnetzwerken bewilligt. Die Forschungsarbeiten an den beteiligten Standorten, der Friedrich-Alexander-Universität Erlangen-Nürnberg (Sprecheruniversität; Sprecher: Alexander Martin), der Technischen Universität Darmstadt (stellvertretender Sprecher: Jens Lang), der Technischen Universität Berlin, der Humboldt Universität (stellvertretende Sprecherin: Caren Tischendorf) sowie den Partnerinstitutionen Weierstraß-Institut (Berlin), Konrad-Zuse-Zentrum (Berlin) und Universität Duisburg-Essen haben im Oktober 2014 begonnen. Y1 - 2015 U6 - https://doi.org/10.1515/dmvm-2015-0013 ER - TY - JOUR A1 - Lang, Jens A1 - Hundsdorfer, Willem T1 - Extrapolation-based implicit–explicit Peer methods with optimised stability regions N2 - In this paper we investigate a new class of implicit–explicit (IMEX) two-step methods of Peer type for systems of ordinary differential equations with both non-stiff and stiff parts included in the source term. An extrapolation approach based on already computed stage values is applied to construct IMEX methods with favourable stability properties. Optimised IMEX-Peer methods of order p=2,3,4, are given as result of a search algorithm carefully designed to balance the size of the stability regions and the extrapolation errors. Numerical experiments and a comparison to other implicit–explicit methods are included. Y1 - 2017 U6 - https://doi.org/doi:10.1016/j.jcp.2017.02.034 VL - J. Comput. Phys. IS - 337 SP - 203 EP - 215 ER - TY - JOUR A1 - Ullmann, Sebastian A1 - Rotkvic, Marko A1 - Lang, Jens T1 - POD-Galerkin reduced-order modeling with adaptive finite element snapshots N2 - We consider model order reduction by proper orthogonal decomposition (POD) for parametrized partial differential equations, where the underlying snapshots are computed with adaptive finite elements. We address computational and theoretical issues arising from the fact that the snapshots are members of different finite element spaces. We propose a method to create a POD-Galerkin model without interpolating the snapshots onto their common finite element mesh. The error of the reduced-order solution is not necessarily Galerkin orthogonal to the reduced space created from space-adapted snapshot. We analyze how this influences the error assessment for POD-Galerkin models of linear elliptic boundary value problems. As a numerical example we consider a two-dimensional convection–diffusion equation with a parametrized convective direction. To illustrate the applicability of our techniques to non-linear time-dependent problems, we present a test case of a two-dimensional viscous Burgers equation with parametrized initial data. Y1 - 2016 U6 - https://doi.org/10.1016/j.jcp.2016.08.018 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 - CHAP A1 - Leugering, Günter A1 - Mophou, Gisèle ED - Schulz, Volker T1 - Instantaneous optimal control of friction dominated flow in a gas-network T2 - DFG-AIMS-Workshop, in Mbour, Senegal, 13-16. March 2017 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 nonlinear model derived from a semi-linear approximation of 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. This involves a p-Laplace-type problem on the graph with p =3/2. We prove well-posedness, existence of optimal controls and derive a first order optimality condition. Y1 - 2017 VL - International Series for Numerical Mathematics PB - Birkhäuser ER - TY - INPR A1 - Domschke, Pia A1 - Hiller, Benjamin A1 - Lang, Jens A1 - Tischendorf, Caren T1 - Modellierung von Gasnetzwerken: Eine Übersicht N2 - Mit dieser Übersicht wollen wir eine Zusammenstellung von unterschiedlichen Modellen zur Beschreibung des Gasflusses in Netzwerken bereitstellen, um den Einstieg in das Thema zu erleichtern. Besonderes Augenmerk wird dabei auf die der Modellierung inneliegende hierarchische Struktur und die detaillierte Beschreibung einzelner Bauteile wie Ventile und Kompressoren gelegt. Daneben finden sich ebenfalls Netzmodellklassen, die auf rein algebraische Relationen aufbauen. Am Ende geben wir einen kurzen Überblick über grundlegende numerische Verfahren und Konzepte zur Behandlung von hyperbolischen Bilanzgleichungen. Wir erheben keinen Anspruch auf Vollständigkeit und verweisen an vielen Stellen auf die bestehende Literatur. Die Idee eines Modellkataloges ist uns im Rahmen der Antragstellung zum SFB/Transregio 154 „Mathematische Modellierung, Simulation und Optimierung am Beispiel von Gasnetzwerken“ gekommen. Wir möchten an dieser Stelle die Förderung durch die DFG dankend erwähnen. KW - Euler-Gleichungen, isotherme Euler-Gleichungen, Modellhierarchie, Netzelemente Y1 - 2017 VL - 2717 ER - TY - JOUR A1 - Hintermüller, Michael A1 - Rautenberg, Carlos A1 - Strogies, Nikolai T1 - Dissipative and Non-dissipative Evolutionary Quasi-variational Inequalities with Gradient Constraints JF - Set-Valued and Variational Analysis N2 - Evolutionary quasi-variational inequality (QVI) problems of dissipative and non-dissipative nature with pointwise constraints on the gradient are studied. A semi-discretization in time is employed for the study of the problems and the derivation of a numerical solution scheme, respectively. Convergence of the discretization procedure is proven and properties of the original infinite dimensional problem, such as existence, extra regularity and non-decrease in time, are derived. The proposed numerical solver reduces to a finite number of gradient-constrained convex optimization problems which can be solved rather efficiently. The paper ends with a report on numerical tests obtained by a variable splitting algorithm involving different nonlinearities and types of constraints. Y1 - 2017 VL - 27 SP - 433 EP - 468 ER - TY - JOUR A1 - Hintermüller, Michael A1 - Strogies, Nikolai T1 - On the consistency of Runge--Kutta methods up to order three applied to the optimal control of scalar conservation laws JF - Numerical Analysis and Optimization N2 - Higher-order Runge-Kutta (RK) time discretization methods for the optimal control of scalar conservation laws are analyzed and numerically tested. The hyperbolic nature of the state system introduces specific requirements on discretization schemes such that the discrete adjoint states associated with the control problem converge as well. Moreover, conditions on the RK-coefficients are derived that coincide with those characterizing strong stability preserving Runge-Kutta methods. As a consequence, the optimal order for the adjoint state is limited, e.g., to two even in the case where the conservation law is discretized by a third-order method. Finally, numerical tests for controlling Burgers equation validate the theoretical results. Y1 - 2017 U6 - https://doi.org/10.20347/WIAS.PREPRINT.2442 SP - 119 EP - 154 ER - TY - JOUR A1 - Gugat, Martin A1 - Herty, Michael A1 - Müller, Siegfried ED - Piccoli, Benedetto T1 - Coupling conditions for the transition from supersonic to subsonic fluid states JF - Networks and Heterogeneous Media (NHM) N2 - We discuss coupling conditions for the p-system in case of a transition from supersonic states to subsonic states. A single junction with adjacent pipes is considered where on each pipe the gas ow is governed by a general p-system. By extending the notion of demand and supply known from traffic fiow analysis we obtain a constructive existence result of solutions compatible with the introduced conditions. KW - p-system KW - coupling conditions KW - network Y1 - 2017 U6 - https://doi.org/10.3934/nhm.2017016 VL - 12 IS - 3 SP - 371 EP - 380 ER - TY - JOUR A1 - Habeck, Oliver A1 - Pfetsch, Marc E. A1 - Ulbrich, Stefan T1 - Global optimization of mixed-integer ODE constrained network problems using the example of stationary gas transport N2 - In this paper we propose a new approach for finding global solutions of mixed-integer nonlinear optimization problems with ordinary differential equation constraints on networks. Instead of using a first discretize then optimize approach, we combine spatial and variable branching with appropriate discretizations of the differential equations to derive relaxations of the original problem. To construct the relaxations we derive convex under- and concave over-estimators for the ODE solution operators using numerical discretization schemes. Thereby, we make use of the underlying network structure, where the solutions of the ODEs only need to be known at a finite number of points. This property enables us to adaptively refine the discretization and relaxation without introducing new variables. The incorporation into a spatial branch-and-bound process allows to compute global epsilon-optimal solutions or decide infeasibility. We prove that this algorithm terminates finitely under some natural assumptions. We then show how this approach works for the example of stationary gas transport and provide some illustrative computational examples. KW - Global Optimization KW - Mixed-Integer Nonlinear Optimization Y1 - 2017 U6 - https://doi.org/10.1137/17M1152668 VL - 29 IS - 4 SP - 2949 EP - 2985 ET - SIAM Journal of Optimization ER - TY - CHAP A1 - Benner, Peter A1 - Grundel, Sarah A1 - Himpe, Christian A1 - Huck, Christoph A1 - Streubel, Tom A1 - Tischendorf, Caren T1 - Gas Network Benchmark Models T2 - Applications of Differential-Algebraic Equations: Examples and Benchmarks N2 - The simulation of gas transportation networks becomes increasingly more important as its use-cases broaden to more complex applications. Classically, the purpose of the gas network was the transportation of predominantly natural gas from a supplier to the consumer for long-term scheduled volumes. With the rise of renewable energy sources, gas-fired power plants are often chosen to compensate for the fluctuating nature of the renewables, due to their on-demand power generation capability. Such an only short-term plannable supply and demand setting requires sophisticated simulations of the gas network prior to the dispatch to ensure the supply of all customers for a range of possible scenarios and to prevent damages to the gas network. In this work we describe the modeling of gas networks and present benchmark systems to test implementations and compare new or extended models. Y1 - 2017 U6 - https://doi.org/10.1007/11221_2018_5 VL - DAE-FORUM SP - 171 EP - 197 PB - Springer CY - Cham ER - TY - CHAP A1 - Huck, Christoph A1 - Tischendorf, Caren T1 - Transient Modeling and Simulation of Gas Pipe Networks with Characteristic Diagram Models for Compressors T2 - Proceedings in Applied Mathematics and Mechanics N2 - One challenge for the simulation and optimization of real gas pipe networks is the treatment of compressors. Their behavior is usually described by characteristic diagrams reflecting the connection of the volumetric flow and the enthalpy change or shaft torque. Such models are commonly used for an optimal control of compressors and compressor stations [4, 7] using stationary models for the gas flow through the pipes. For transient simulations of gas networks, simplified compressor models have been studied in [1–3]. Here, we present a transient simulation of gas pipe networks with characteristic diagram models of compressors using a stable network formulation as (partial) differential-algebraic system. Y1 - 2017 U6 - https://doi.org/10.1002/pamm.201710322 VL - 17 IS - 1 SP - 707 EP - 708 ER - TY - INPR A1 - Huck, Christoph A1 - Tischendorf, Caren T1 - Perturbation analysis of hyperbolic PDAEs describing flow networks 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 - 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 - 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 -