TY - JOUR A1 - Lang, Jens A1 - Schmitt, Bernhard A. T1 - A Stiff MOL Boundary Control Problem for the 1D Heat Equation with Exact Discrete Solution N2 - Method-of-lines discretizations are demanding test problems for stiff inte- gration methods. However, for PDE problems with known analytic solution the presence of space discretization errors or the need to use codes to compute reference solutions may limit the validity of numerical test results. To over- come these drawbacks we present in this short note a simple test problem with boundary control, a situation where one-step methods may suffer from order reduction. We derive exact formulas for the solution of an optimal boundary control problem governed by a one-dimensional discrete heat equation and an objective function that measures the distance of the final state from the target and the control costs. This analytical setting is used to compare the numeri- cally observed convergence orders for selected implicit Runge-Kutta and Peer two-step methods of classical order four which are suitable for optimal control problems. Y1 - U6 - https://doi.org/https://doi.org/10.1007/s10957-022-02154-4 VL - Journal of Optimization Theory and Applications IS - Vol. 196 SP - 1106 EP - 1118 ER - TY - JOUR A1 - Strelow, Erik Laurin A1 - Gerisch, Alf A1 - Lang, Jens A1 - Pfetsch, Marc E. T1 - Physics-Informed Neural Networks: A Case Study for Gas Transport Problems N2 - Physics informed neural networks have been recently proposed and offer a new promising method to solve differential equations. They have been adapted to many more scenarios and different variations of the original method have been proposed. In this case study we review many of these variations. We focus on variants that can compensate for imbalances in the loss function and perform a comprehensive numerical comparison of these variants with application to gas transport problems. Our case study includes different formulations of the loss function, different algorithmic loss balancing methods, different optimization schemes and different numbers of parameters and sampling points. We conclude that the original PINN approach with specifically chosen constant weights in the loss function gives the best results in our tests. These weights have been obtained by a computationally expensive random-search scheme. We further conclude for our test case that loss balancing methods which were developed for other differential equations have no benefit for gas transport problems, that the control volume physics informed formulation has no benefit against the initial formulation and that the best optimization strategy is the L-BFGS method. Y1 - VL - Journal of Computational Physics IS - Vol. 481 SP - 112041 ER - TY - JOUR A1 - Lang, Jens A1 - Schmitt, Bernhard A. T1 - Implicit A-Stable Peer Triplets for ODE Constrained Optimal Control Problems N2 - This paper is concerned with the construction and convergence analysis of novel implicit Peer triplets of two-step nature with four stages for nonlinear ODE constrained optimal control problems. We combine the property of superconvergence of some standard Peer method for inner grid points with carefully designed starting and end methods to achieve order four for the state variables and order three for the adjoint variables in a first-discretize-then-optimize approach together with A-stability. The notion triplets emphasizes that these three different Peer methods have to satisfy additional matching conditions. Four such Peer triplets of practical interest are constructed. Also as a benchmark method, the well-known backward differentiation formula BDF4, which is only A(73.35)-stable, is extended to a special Peer triplet to supply an adjoint consistent method of higher order and BDF type with equidistant nodes. Within the class of Peer triplets, we found a diagonally implicit A(84)-stable method with nodes symmetric in [0,1] to a common center that performs equally well. Numerical tests with three well established optimal control problems confirm the theoretical findings also concerning A-stability. Y1 - U6 - https://doi.org/https://doi.org/10.3390/a15090310 VL - Algorithms IS - Vol. 15 ER - TY - INPR A1 - Domschke, Pia A1 - Giesselmann, Jan A1 - Lang, Jens A1 - Breiten, Tobias A1 - Mehrmann, Volker A1 - Morandin, Riccardo A1 - Hiller, Benjamin A1 - Tischendorf, Caren T1 - Gas Network Modeling: An Overview (Extended English Version) N2 - With this overview we want to provide a compilation of different models for the description of gas flow in networks in order to facilitate the introduction to the topic. Special attention is paid to the hierarchical structure inherent to the modeling, and the detailed description of individual components such as valves and compressors. Also included are network model classes based on purely algebraic relations, and energy-based port-Hamiltonian models. A short overview of basic numerical methods and concepts for the treatment of hyperbolic balance equations is also given. We do not claim completeness and refer in many places to the existing literature. Y1 - 2023 ER - TY - INPR A1 - Lang, Jens A1 - Schmitt, Bernhard A. T1 - Implicit Peer Triplets in Gradient-Based Solution Algorithms for ODE Constrained Optimal Control N2 - It is common practice to apply gradient-based optimization algorithms to numerically solve large-scale ODE constrained optimal control problems. Gradients of the objective function are most efficiently computed by approximate adjoint variables. High accuracy with moderate computing time can be achieved by such time integration methods that satisfy a sufficiently large number of adjoint order conditions and supply gradients with higher orders of consistency. In this paper, we upgrade our former implicit two-step Peer triplets constructed in [Algorithms, 15:310, 2022] to meet those new requirements. Since Peer methods use several stages of the same high stage order, a decisive advantage is their lack of order reduction as for semi-discretized PDE problems with boundary control. Additional order conditions for the control and certain positivity requirements now intensify the demands on the Peer triplet. We discuss the construction of 4-stage methods with order pairs (4,3) and (3,3) in detail and provide three Peer triplets of practical interest. We prove convergence for s-stage methods, for instance, order s for the state variables even if the adjoint method and the control satisfy the conditions for order s-1, only. Numerical tests show the expected order of convergence for the new Peer triplets. Y1 - 2023 VL - http://arxiv.org/abs/2303.18180 ER - TY - INPR A1 - Graser, Gertrud A1 - Kreimeier, Timo A1 - Walther, Andrea T1 - Solving Linear Generalized Nash Games Using an Active Signature Method N2 - We propose a method to solve linear generalized Nash equilibrium problems (LGNEPs). For this purpose, a reformulation of the LGNEPs as piecewise linear problems is considered. This requires the calculation of all vertices for a special kind of unbounded convex polyhedra. Then the active signature method for constrained abs-linear problems can be used to determine the Nash equilibria. We analyse the computational effort for the resulting solution procedure. This includes also the verification of suitable optimality conditions. Finally, we present and analyse numerical results for some test problems. Y1 - 2024 ER - TY - INPR A1 - Alldredge, Graham A1 - Frank, Martin A1 - Giesselmann, Jan T1 - On the convergence of the regularized entropy-based moment method for kinetic equations N2 - The entropy-based moment method is a well-known discretization for the velocity variable in kinetic equations which has many desirable theoretical properties but is difficult to implement with high-order numerical methods. The regularized entropy-based moment method was recently introduced to remove one of the main challenges in the implementation of the entropy-based moment method, namely the requirement of the realizability of the numerical solution. In this work we use the method of relative entropy to prove the convergence of the regularized method to the original method as the regularization parameter goes to zero and give convergence rates. Our main assumptions are the boundedness of the velocity domain and that the original moment solution is Lipschitz continuous in space and bounded away from the boundary of realizability. We provide results from numerical simulations showing that the convergence rates we prove are optimal. Y1 - 2023 U6 - https://doi.org/https://doi.org/10.5802/smai-jcm.93 VL - 9 ER - TY - INPR A1 - Giesselmann, Jan A1 - Kolbe, Niklas T1 - A posteriori error analysis of a positivity preserving scheme for the power-law diffusion Keller-Segel model N2 - We study a finite volume scheme approximating a parabolic-elliptic Keller-Segel system with power law diffusion with exponent γ∈[1,3] and periodic boundary conditions. We derive conditional a posteriori bounds for the error measured in the L∞(0,T;H1(Ω)) norm for the chemoattractant and by a quasi-norm-like quantity for the density. These results are based on stability estimates and suitable conforming reconstructions of the numerical solution. We perform numerical experiments showing that our error bounds are linear in mesh width and elucidating the behaviour of the error estimator under changes of γ. KW - Keller-Segel KW - chemotaxis; KW - nonlinear diffusion KW - finite volume scheme KW - a posteriori error analysis Y1 - 2023 ER - TY - INPR A1 - Giesselmann, Jan A1 - Kwon, Kiwoong T1 - A posteriori error control for a Discontinuous Galerkin approximation of a Keller-Segel model N2 - We provide a posteriori error estimates for a discontinuous Galerkin scheme for the parabolic-elliptic Keller-Segel system in 2 or 3 space dimensions. The estimates are conditional, in the sense that an a posteriori computable quantity needs to be small enough - which can be ensured by mesh refinement - and optimal in the sense that the error estimator decays with the same order as the error under mesh refinement. A specific feature of our error estimator is that it can be used to prove existence of a weak solution up to a certain time based on numerical results. KW - Keller-Segel KW - chemotaxis KW - nonlinear diffusion KW - discontinuous Galerkin scheme KW - a posteriori error analysis Y1 - 2023 ER - TY - INPR A1 - Giesselmann, Jan A1 - Krupa, Sam T1 - Theory of shifts, shocks, and the intimate connections to L2-type a posteriori error analysis of numerical schemes for hyperbolic problems N2 - In this paper, we develop reliable a posteriori error estimates for numerical approximations of scalar hyperbolic conservation laws in one space dimension. Our methods have no inherent small-data limitations and are a step towards error control of numerical schemes for systems. We are careful not to appeal to the Kruzhkov theory for scalar conservation laws. Instead, we derive novel quantitative stability estimates that extend the theory of shifts, and in particular, the framework for proving stability first developed by the second author and Vasseur. This is the first time this methodology has been used for quantitative estimates. We work entirely within the context of the theory of shifts and a-contraction, techniques which adapt well to systems. In fact, the stability framework by the second author and Vasseur has itself recently been pushed to systems [Chen-Krupa-Vasseur. Uniqueness and weak-BV stability for 2×2 conservation laws. Arch. Ration. Mech. Anal., 246(1):299--332, 2022]. Our theoretical findings are complemented by a numerical implementation in MATLAB and numerical experiments. KW - Conservation laws KW - entropy conditions KW - entropy solutions KW - shocks, KW - a posteriori error estimates Y1 - 2023 ER - TY - INPR A1 - Egger, Herbert A1 - Giesselmann, Jan T1 - Regularity and long time behavior of a doubly nonlinear parabolic problem and its discretization N2 - We study a doubly nonlinear parabolic problem arising in the modeling of gas transport in pipelines. Using convexity arguments and relative entropy estimates we show uniform bounds and exponential stability of discrete approximations obtained by a finite element method and implicit time stepping. Due to convergence of the approximations to weak solutions of the problem, our results also imply regularity, uniqueness, and long time stability of weak solutions of the continuous problem. KW - gas transport KW - doubly nonlinear parabolic problems KW - relative entropy estimates KW - exponential stability KW - structure preserving discretization Y1 - 2023 ER - TY - INPR A1 - Gugat, Martin A1 - Giesselmann, Jan T1 - An Observer for pipeline flow with hydrogen blending in gas networks: exponential synchronization N2 - We consider a state estimation problem for gas flows in pipeline networks where hydrogen is blended into the natural gas. The flow is modeled by the quasi-linear isothermal Euler equations coupled to an advection equation on a graph. The flow through the vertices where the pipes are connected is governed by algebraic node conditions. The state is approximated by an observer system that uses nodal measurements. We prove that the state of the observer system converges to the original system state exponentially fast in the L2-norm if the measurements are exact. If measurement errors are present we show that the observer state approximates the original system state up to an error that is proportional to the maximal measurement error. The proof of the synchronization result uses Lyapunov functions with exponential weights. Y1 - 2023 ER - TY - INPR A1 - Bernhard, Daniela A1 - Liers, Frauke A1 - Stingl, Michael A1 - Uihlein, Andrian T1 - A Gradient-Based Method for Joint Chance-Constrained Optimization with Continuous Distributions N2 - The input parameters of an optimization problem are often affected by uncertainties. Chance constraints are a common way to model stochastic uncertainties in the constraints. Typically, algorithms for solving chance-constrained problems require convex functions or discrete probability distributions. In this work, we go one step further and allow non-convexities as well as continuous distributions. We propose a gradient-based approach to approximately solve joint chance-constrained models. We approximate the original problem by smoothing indicator functions. Then, the smoothed chance constraints are relaxed by penalizing their violation in the objective function. The approximation problem is solved with the Continuous Stochastic Gradient method that is an enhanced version of the stochastic gradient descent and has recently been introduced in the literature. We present a convergence theory for the smoothing and penalty approximations. Under very mild assumptions, our approach is applicable to a wide range of chance-constrained optimization problems. As an example, we illustrate its computational efficiency on difficult practical problems arising in the operation of gas networks. The numerical experiments demonstrate that the approach quickly finds nearly feasible solutions for joint chance-constrained problems with non-convex constraint functions and continuous distributions, even for realistically-sized instances. Y1 - ER - TY - JOUR A1 - Geiersbach, Caroline A1 - Henrion, René T1 - Optimality conditions in control problems with random state constraints in probabilistic or almost-sure form N2 - In this paper, we discuss optimality conditions for optimization problems {involving} random state constraints, which are modeled in probabilistic or almost sure form. While the latter can be understood as the limiting case of the former, the derivation of optimality conditions requires substantially different approaches. We apply them to a linear elliptic partial differential equation (PDE) with random inputs. In the probabilistic case, we rely on the spherical-radial decomposition of Gaussian random vectors in order to formulate fully explicit optimality conditions involving a spherical integral. In the almost sure case, we derive optimality conditions and compare them to a model based on robust constraints with respect to the (compact) support of the given distribution. Y1 - 2023 ER - TY - JOUR A1 - Geiersbach, Caroline A1 - Henrion, René A1 - Pérez-Aros, Pedro T1 - Numerical solution of an optimal control problem with probabilistic and almost sure state constraints N2 - We consider the optimal control of a PDE with random source term subject to probabilistic or almost sure state constraints. In the main theoretical result, we provide an exact formula for the Clarke subdifferential of the probability function without a restrictive assumption made in an earlier paper. The focus of the paper is on numerical solution algorithms. As for probabilistic constraints, we apply the method of spherical radial decomposition. Almost sure constraints are dealt with a Moreau--Yosida smoothing of the constraint function accompanied by Monte Carlo sampling of the given distribution or its support or even just the boundary of its support. Moreover, one can understand the almost sure constraint as a probabilistic constraint with safety level one which offers yet another perspective. Finally, robust optimization can be applied efficiently when the support is sufficiently simple. A comparative study of these five different methodologies is carried out and illustrated. Y1 - 2023 ER - TY - INPR A1 - Hante, Falk A1 - Kuchler, Christian T1 - Indirect methods for optimal control of parabolic hybrid PDE-dynamical / switching systems using relaxation N2 - We propose a novel algorithmic approach to computationally solve optimal control problems governed by linear parabolic partial differential equations (PDEs) including a state-dependent control-regime switching mechanism. We state an equivalent mixed-integer formulation featuring vanishing constraints (VCs) arising from methods of disjunctive programming. We embed the problem into the class of equilibrium constraints (ECs) by introduction of an additional slack variable. Based on theoretical results associated with Sum-Up-Rounding (SUR) strategies, we proceed with the solution of the related relaxed formulation by an indirect approach. In order to obtain a computationally tractable optimality system, we apply a Moreau-Yosida type penalty approach for the VCs. After a theoretical discussion, we introduce and exert the algorithmic framework founded on a semismooth Newton method. Finally, we communicate computational experiments based on the proposed approach. Y1 - ER - TY - INPR A1 - Grimm, Veronika A1 - Grübel, Julia A1 - Schmidt, Martin A1 - Schwartz, Alexandra A1 - Wiertz, Ann-Kathrin A1 - Zöttl, Gregor T1 - On a Tractable Single-Level Reformulation of a Multilevel Model of the European Entry-Exit Gas Market with Market Power N2 - We propose a framework that allows to quantitatively analyze the interplay of the different agents involved in gas trade and transport in the context of the European entry-exit system. Previous contributions have focused on the case of perfectly competitive buyers and sellers of gas, which allows to replace the respective market equilibrium problem by a single welfare maximization problem. Our novel framework considers the mathematically more challenging case of a monopolistic and thus strategic gas seller. In this framework, the objective functions of the gas sellers and buyers cannot be aggregated into a common objective function, which is why a multilevel formulation is necessary to accurately capture the sequential nature of the decisions taken. For this setup, we derive sufficient conditions that allow for reformulating the challenging four-level model as a computationally tractable single-level reformulation. We prove the correctness of this reformulation and use it for solving several test instances to illustrate the applicability of our approach. KW - Multilevel optimization KW - Reformulations KW - Gas markets KW - Market power Y1 - 2023 ER - TY - JOUR A1 - Shyshkanova, Ganna A1 - Zaytseva, Tetyana A1 - Zhushman, V A1 - Levchenko, Ntaliia A1 - Korotunova, Olena T1 - Solving three-dimensional contact problems for foundation design in green building N2 - Design of foundations on an elastic base is carried out using the solution of three-dimensional problems of contact interaction. Improving the accuracy of engineering calculations is necessary to ensure economic efficiency and increase energy savings in green building. The problems of indentation of punches with a flat base bounded by doubly connected close to polygonal contact areas are researched in the present work. Small parameter method is used to obtain explicit analytical expressions for the contact pressure distribution and the punch displacement dependence in a simplified form, which is convenient for engineering practice. The found load-displacement dependence satisfies the known inequalities that are valid for an arbitrary contact domain. Also a numerical-analytical method is in consideration. It uses the simple layer potential expansion and successive approximations for the problems accounting roughness of the elastic half-space. Roughness coefficient is considered as a parameter of regularization of the integral equation for the smooth contact problem. The results of both methods coincide with sufficient accuracy. Y1 - 2023 U6 - https://doi.org/10.1088/1742-6596/2609/1/012001 ER - TY - THES A1 - Krug, Richard T1 - Decomposition Methods for Time-Dependent Mixed-Integer Nonlinear Optimization Problems on Graphs N2 - Decomposition can be the method of choice to deal with optimization problems that contain hard to solve model structures or that are of large scale. The main idea is to decompose the problematic aspects of the problem into multiple smaller blocks that can be solved more easily. Here, the challenge is to combine the single pieces to a solution that is not only feasible but maybe even optimal for the original problem. In many cases, this can be done by introducing an iteration that eventually converges to a desired solution. In this cumulative dissertation, we present several iterative decomposition methods that are tailored to different types of optimization models and use distinct approaches to split up the problems. Our main motivation for this originates from the optimization of gas transport networks, where we encounter partial differential equations as well as discrete control decisions. Additionally, we engage in the related field of district heating network optimization to study the challenges arising from large-scale and fully discretized systems as well as undesirable model features such as, e.g., complementarity constraints. Here, we introduce two temperature mixing models that are well suited for optimization and a number of techniques to speed up the solution process, which are applied in numerical experiments. As a next step, we develop an iterative time-domain decomposition method that is applied to optimal control problems subject to semilinear hyperbolic systems of partial differential equations. For this, we derive first-order optimality conditions that are then split using a non-overlapping decomposition of the time horizon. We exploit the fact that the resulting systems have a primal interpretation as so-called virtual control problems. We prove the convergence of the iterative method and develop a posteriori error estimates. Later, we extend the scheme to systems of ordinary differential equations with mixed- integer controls by using Pontryagin’s maximum principle. We again show the convergence and conduct a numerical case study. Moreover, we use a consensus-based version of the classic penalty alternating direction method to solve tailored reformulations of transient gas network problems that allow us to minimize the number of coupling constraints between sub-problems. Here, we utilize the quasi-separable structure of the network to decompose it into sub-networks with more desirable properties. We also discuss different decomposition strategies and test them in a numerical case study. Finally, we present a successive linear relaxation method for mixed-integer nonlinear problems with multivariate Lipschitz continuous nonlinearities. The distinguishing feature of this algorithm is that it exploits no properties of the nonlinearities besides the Lipschitz constants. Therefore, the method is applicable for problems with non-convex or even non-differentiable constraints. The nonlinearities do not even need to be given in a closed form, which allows us to integrate black-box constraints into the model. We prove that the algorithm converges to an approximate global optimum and we provide a worst-case estimate for the number of iterations. The iterative method is applied to stationary gas transport problems, where implicitly given solutions of the differential equations are modeled via black-box constraints. Y1 - 2023 ER - TY - JOUR A1 - Ouanes, Nesrine A1 - González Grandón, Tatiana A1 - Heitsch, Holger A1 - Henrion, René T1 - Optimizing the economic dispatch of weakly-connected mini-grids under uncertainty using joint chance constraints N2 - In this paper, we deal with a renewable-powered mini-grid, connected to an unreliable main grid, in a Joint Chance Constrained (JCC) programming setting. In many countries with low energy access rates, grid-connected mini-grid system operators contend with four different types of uncertainties: stochastic solar power and demand forecast errors; absolute uncertain national grid outage onset times; and outages duration subjected to statistical analysis. These uncertainties pose new challenges to the classical power system’s operation tasks. Two alternatives to the JCC problem are presented. In particular, we present an Individual Chance Constraint (ICC) and a purely deterministic dispatch model. The JCC model has the capability to address all four uncertainties, while the ICC covers only three of them, overlooking the uncertainty about the outage duration. In contrast, the purely deterministic model completely ignores any uncertain parameters. We illustrate the three models through a comparison of outcomes attained from a real mini-grid in Lake Victoria, Tanzania. Results show how the dispatch is modified across the models to plan the battery and diesel reserves in the chance-constrained models, with the reserves in the JCC being larger than in the ICC model. In comparison between all models, we prove that the JCC model offers the most robust results, since it can handle uncertainties about forecasting errors, on the one hand, and grid outages, on the other. The results also show that the decrease in profits due to the hedging with reserves kept in the MG is significantly small compared to the high level of reliability reached and the potential load shedding that could be avoided in the case of an outage. Y1 - 2023 ER - TY - JOUR A1 - Shyshkanova, Ganna A1 - Walther, Andrea T1 - Optimization of a punch shape with a doubly connected contact domain N2 - The objective is to optimize the pressure distribution under a rigid punch having a doubly connected contact domain close to a circular ring and interacting with an elastic half-space. The required design variable is the punch shape. The functional to be minimized is the root-mean-square deviation of the pressure distribution from some given distribution. An analytical technique is developed for solving the problem for the punches with doubly connected shape, by reducing to a sequence of similar problems for the circular ring punches using expansions of the simple layer potential. The method of expansion in terms of a small parameter is used. The simple layer potential expansion is proposed when mapping a doubly connected integration domain onto a circular ring by transforming the integration variables and transforming the coordinates of the pole of the kernel. As a result, a sequence of similar problems was obtained for a circular ring to determine the functions characterizing the distribution of normal pressure under the punch in the form of a non-circular ring, as well as the normal displacements, from where the optimal punch shape is determined. Y1 - 2023 ER - TY - INPR A1 - Goerigk, Marc A1 - Kurtz, Jannis A1 - Schmidt, Martin A1 - Thürauf, Johannes T1 - Connections between Robust and Bilevel Optimization N2 - Robust and bilevel optimization share the common feature that they involve a certain multilevel structure. Hence, although they model something rather different when used in practice, they seem to have a similar mathematical structure. In this paper, we analyze the connections between different types of robust problems (static robust problems with and without decision-dependence of their uncertainty sets, worst-case regret problems, and two-stage robust problems) as well as of bilevel problems (optimistic problems, pessimistic problems, and robust bilevel problems). It turns out that bilevel optimization seems to be more general in the sense that for most types of robust problems, one can find proper reformulations as bilevel problems but not necessarily the other way around. We hope that these results pave the way for a stronger connection between the two fields - in particular to use both theory and algorithms from one field in the other and vice versa. KW - Bilevel optimization KW - Robust optimization KW - Reformulations Y1 - 2023 ER - TY - INPR A1 - Gugat, Martin A1 - Qian, Meizhi A1 - Sokolowski, Jan T1 - Topological derivative method for control of wave equation on networks N2 - The dynamical, boundary optimal control problems on networks are considered. The domain of definition for the distributed parameter system is given by a graph G. The optimal cost function for control problem is further optimized with respect to the shape and topology of the graph Ω. The small cycle is introduced and the topological derivative of the cost with respect to the size of the cycle is determined. In this way, the singular perturbations of the graph can be analyzed in order to change the topology Ω. The topological derivative method in shape and topology optimization is a new tool which can be used to minimize the shape functionals under the Partial Differential Equations (PDEs) constraints. The topological derivative is used as well for solution of optimum design problems for graphs. In optimal control problems the topological derivative is used for optimum design of the domain of integration of the state equation. As an example, optimal control problems are considered on a cross with a small cycle. The state equation is the wave equation on the graph. The boundary control problem by Neumann conditions at a boundary vertex is solved for a tracking cost function. The shape functional is given by the optimal value of the control cost. The topological derivative of the shape functional is determined for the steady state model with the size of a cycle ε → 0. Numerical results for a model problem are presented. KW - distributed parameter system KW - optimal control KW - shape optimization KW - topological derivative KW - network modelling Y1 - 2023 ER - TY - INPR A1 - Kannan, Aswin A1 - Kreimeier, Timo A1 - Walther, Andrea T1 - On Solving Nonsmooth Retail Portfolio Maximization Problems Using Active Signature Methods Y1 - 2023 ER - TY - INPR A1 - Hante, Falk A1 - Kuchler, Christian T1 - An Algorithmic Framework for Optimal Control of Hybrid Dynamical System with Parabolic PDEs N2 - We present an algorithmic approach for the computational solution of optimal control problems with hybrid nature governed by linear parabolic PDEs featuring implicit switches. We propose a stepwise reformulation of the original formulation into a more tractable setting via application of methods from disjunctive programming and a time transformation method. After removal of the implicit switching rule at the cost of the introduction of explicit switching variables and vanishing constraints, the connection of the resulting formulation to problems with equilibrium constraints is established and studied. The previous steps in combination with smoothening and a Moreau-Yosida type penalty approach allow the derivation of necessary first order optimality conditions to characterize candidates for optimality to the original system. Following the discussion of each individual reformulation step, we introduce the algorithmic framework founded on a semismooth Newton method. Finally, we report on computational of the proposed framework. Y1 - 2023 ER - TY - INPR A1 - Göß, Adrian A1 - Martin, Alexander A1 - Pokutta, Sebastian A1 - Sharma, Kartikey T1 - Norm-induced Cuts: Optimization with Lipschitzian Black-box Functions N2 - Optimal control problems usually involve constraints which model physical states and their possible transitions. These are represented by ordinary or partial differential equations (ODEs/PDEs) which add a component of infinite dimension to the problem. In recent literature, one method to simulate such ODEs/PDEs are physics-informed neural networks. Typically, neural networks are highly non-linear which makes their addition to optimization problems challenging. Hence, we leverage their often available Lipschitz property on a compact domain. The respective Lipschitz constants have to be computed only once and are accessible thereafter. We present a method that, based on this property, iteratively adds cuts involving the violation of the constraints by the current incumbent and the Lipschitz constant. Hereby, the “shape” of a cut depends on the norm used. We prove the correctness of the method by showing that it either returns an optimal solution when terminating or creates a sequence with optimal accumulation points. This is complemented by a discussion about the termination in the infeasible case, as well as an analysis of the problem complexity. For the analysis, we show that the lower and upper iteration bound asymptotically coincide when the relative approximation error goes to zero. In the end, we visualize the method on a small example based on a two-dimensional non-convex optimization problem, as well as stress the necessity of having a globally optimal oracle for the sub-problems by another example. KW - Global Optimization KW - Lipschitz Optimization KW - Black-box Optimization KW - Derivative-free Optimization Y1 - ER - TY - INPR A1 - Bongarti, Marcelo A1 - Hintermüller, T1 - Optimal boundary control of the isothermal semilinear Euler equation for gas dynamics on a network N2 - The analysis and boundary optimal control of the nonlinear transport of gas on a network of pipelines is considered. The evolution of the gas distribution on a given pipe is modeled by an isothermal semilinear compressible Euler system in one space dimension. On the network, solutions satisfying (at nodes) the Kirchhoff flux continuity conditions are shown to exist in a neighborhood of an equilibrium state. The associated nonlinear optimization problem then aims at steering such dynamics to a given target distribution by means of suitable (network) boundary controls while keeping the distribution within given (state) constraints. The existence of local optimal controls is established and a corresponding Karush-Kuhn-Tucker (KKT) stationarity system with an almost surely non-singular Lagrange multiplier is derived. KW - optimal boundary control KW - gas dynamics KW - gas networks KW - isothermal Euler equation KW - compressible fluid dynamics Y1 - 2023 ER - TY - CHAP A1 - Gugat, Martin A1 - Schuster, Michael T1 - Max-p optimal boundary control of gas flow N2 - In the transition to renewable energy sources, hydrogen will potentially play an important role for energy storage. The efficient transport of this gas is possible via pipelines. An understanding of the possibilities to control the gas flow in pipelines is one of the main building blocks towards the optimal use of gas. For the operation of gas transport networks it is important to take into account the randomness of the consumers’ demand, where often information on the probability distribution is available. Hence in an efficient optimal control model the corresponding probability should be included and the optimal control should be such that the state that is generated by the optimal control satisfies given state constraints with large probability. We comment on the modelling of gas pipeline flow and the problems of optimal nodal control with random demand, where the aim of the optimization is to determine controls that generate states that satisfy given pressure bounds with large probability. We include the H2 norm of the control as control cost, since this avoids large pressure fluctuations which are harmful in the transport of hydrogen since they can cause embrittlement of the pipeline metal. KW - gas pipeline flow KW - nodal control KW - hyperbolic differential equation KW - random demand KW - state constraints Y1 - 2022 U6 - https://doi.org/https://doi.org/10.15495/EPub_UBT_00006809 VL - Extended Abstracts of the 25th International Symposium on Mathematical Theory of Networks and Systems Bayreuth, Germany, 12-16 September 2022 ER - TY - JOUR A1 - Gugat, Martin A1 - Lazar, Martin T1 - Turnpike Properties for Partially Uncontrollable Systems N2 - We analyse the turnpike properties for a general, infinite dimensional, linear-quadratic (LQ) optimal control problem, both in the deterministic and in the stochastic case. The novelty of the paper is twofold. Firstly, it obtains positive turnpike results for systems that are (partially) uncontrollable. Secondly, it provides turnpike results for averaged control associated to a family of problems that depend on a random parameter, which is the first turnpike type result in the averaged controllability framework. KW - Measure Turnpike KW - Averaged Control KW - LQ optimal control problem KW - Infinite-time admissibility KW - Turnpike phenomenon Y1 - 2023 VL - Automatica IS - 149 ER - TY - INPR A1 - Schuster, Michael A1 - Sakamoto, Noboru T1 - A Turnpike Result for Optimal Boundary Control Problems with the Transport Equation under Uncertainty N2 - In this paper we analyze the turnpike phenomenon for optimal boundary control problems with a linear transport equation with source term. The convex objective function depends on the boundary traces of the transport equation and is strictly convex with respect to the boundary control. We show an integral turnpike result for an optimal Dirichlet boundary control problem in the sense that if the time horizon goes to infinity, then the dynamic optimal control converges to the corresponding steady state optimal control. The novelty of this work is two-sided. On the one hand, even if turnpike results for this kind of optimal boundary control problem already exist, we present a new direct proof without using adjoint calculus that leads to sharper estimates. On the other hand we consider uncertainty in the initial data and/or in the source term. We show that the integral turnpike result also holds considering uncertainty. Throughout the paper we use numerical examples to illustrate the results. KW - Turnpike KW - Boundary Control KW - Transport Equation KW - Random Boundary Data Y1 - 2023 ER - TY - JOUR A1 - Giesselmann, Jan A1 - Gugat, Martin A1 - Kunkel, Teresa T1 - Observer-based data assimilation for barotropic gas transport using distributed measurements N2 - We consider a state estimation problem for gas pipeline flow modeled by the one-dimensional barotropic Euler equations. In order to reconstruct the system state, we construct an observer system of Luenberger type based on distributed measurements of one state variable. First, we show the existence of Lipschitz-continuous semi-global solutions of the observer system and of the original system for initial and boundary data satisfying smallness and compatibility conditions for a single pipe and for general networks. Second, based on an extension of the relative energy method we prove that the state of the observer system converges exponentially in the long time limit towards the original system state. We show this for a single pipe and for star-shaped networks. Y1 - 2023 U6 - https://doi.org/10.4310/CMS.240918203214 VL - 22 SP - 2271 EP - 2309 PB - Communications in Mathematical Sciences ER - TY - JOUR A1 - Egerer, Jonas A1 - Grimm, Veronika A1 - Niazmand, Kiana A1 - Runge, Philipp T1 - The economics of global green ammonia trade – "Shipping Australian wind and sunshine to Germany" JF - Applied Energy N2 - This paper contributes to understanding the transformation of global energy trade to green energy carriers, focusing on green ammonia as the foreseeable first green hydrogen carrier. We provide a comprehensive overview of today's ammonia trade and assess scaling options for the trade of green ammonia. To that aim, we develop an optimization model for the integrated assessment of the green ammonia value chain that covers all steps from green ammonia production in an exporting country, up to delivery to a harbor in an importing country. The model endogenously chooses among different technology options and determines cost minimal operation. In a case study, we apply the model to the large-scale import of ammonia from Australia to Germany in a scenario for 2030. The results show that green ammonia can reach cost parity with gray ammonia even for moderate gas prices (but not necessarily with blue ammonia) if CO2 prices are high enough. We also provide a sensitivity analysis with respect to the interest rate and other key technical and economic parameters and show that cracking ammonia to provide pure hydrogen comes at a 45 % cost markup per MWh at the destination. KW - green ammonia KW - ammonia trade KW - optimization model KW - case study Australia to Germany Y1 - 2023 U6 - https://doi.org/https://doi.org/10.1016/j.apenergy.2023.120662 VL - 334 ER - TY - INPR A1 - Hante, Falk M. A1 - Schmidt, Martin T1 - Gas Transport Network Optimization: Mixed-Integer Nonlinear Models N2 - Although modern societies strive towards energy systems that are entirely based on renewable energy carriers, natural gas is still one of the most important energy sources. This became even more obvious in Europe with Russia's 2022 war against the Ukraine and the resulting stop of gas supplies from Russia. Besides that it is very important to use this scarce resource efficiently. To this end, it is also of significant relevance that its transport is organized in the most efficient, i.e., cost- or energy-efficient, way. The corresponding mathematical optimization models have gained a lot of attention in the last decades in different optimization communities. These models are highly nonlinear mixed-integer problems that are constrained by algebraic constraints and partial differential equations (PDEs), which usually leads to models that are not tractable. Hence, simplifications have to be made and in this chapter, we present a commonly accepted finite-dimensional stationary model, i.e., a model in which the steady-state solutions of the PDEs are approximated with algebraic constraints. For more details about the involved PDEs and the treatment of transient descriptions we refer to Hante and Schmidt (2023). The presented finite-dimensional as well as mixed-integer nonlinear and nonconvex model is still highly challenging if it needs to be solved for real-world gas transport networks. Hence, we also review some classic solution approaches from the literature. KW - Gas networks KW - Mixed-integer nonlinear optimization KW - Mixed-integer linear optimization KW - Nonlinear optimization Y1 - 2023 ER - TY - INPR A1 - Hante, Falk M. A1 - Schmidt, Martin T1 - Gas Transport Network Optimization: PDE-Constrained Models N2 - The optimal control of gas transport networks was and still is a very important topic for modern economies and societies. Accordingly, a lot of research has been carried out on this topic during the last years and decades. Besides mixed-integer aspects in gas transport network optimization, one of the main challenges is that a physically and technically detailed modeling of transient gas dynamics leads to theoretically and computationally highly demanding models involving nonlinear partial differential equations (PDEs). For further background on the application, historical notes and a detailed discussion of mixed-integer aspects for stationary descriptions we refer to Hante and Schmidt (2023). In this chapter, we focus on the most common modeling approaches concerning transient descriptions, point out the challenges, and summarize important contributions concerning the optimization of the most relevant control parameters for this particular class of problems. KW - Gas networks KW - Partial differential equations KW - Optimal control KW - PDE-constrained optimization KW - Modeling Y1 - 2023 ER - TY - INPR A1 - Krug, Richard A1 - Leugering, Günter A1 - Martin, Alexander A1 - Schmidt, Martin A1 - Weninger, Dieter T1 - A Consensus-Based Alternating Direction Method for Mixed-Integer and PDE-Constrained Gas Transport Problems N2 - We consider dynamic gas transport optimization problems, which lead to large-scale and nonconvex mixed-integer nonlinear optimization problems (MINLPs) on graphs. Usually, the resulting instances are too challenging to be solved by state-of-the-art MINLP solvers. In this paper, we use graph decompositions to obtain multiple optimization problems on smaller blocks, which can be solved in parallel and which may result in simpler classes of optimization problems since not every block necessarily contains mixed-integer or nonlinear aspects. For achieving feasibility at the interfaces of the several blocks, we employ a tailored consensus-based penalty alternating direction method. Our numerical results show that such decomposition techniques can outperform the baseline approach of just solving the overall MINLP from scratch. However, a complete answer to the question of how to decompose MINLPs on graphs in dependence of the given model is still an open topic for future research. KW - Gas transport networks KW - Mixed-integer nonlinear optimization KW - Alternating direction methods KW - Graph decomposition KW - Penalty methods Y1 - 2022 ER - TY - INPR A1 - Schmidt, Martin A1 - Thürauf, Johannes T1 - An Exact Method for Nonlinear Network Flow Interdiction Problems N2 - We study network flow interdiction problems with nonlinear and nonconvex flow models. The resulting model is a max-min bilevel optimization problem in which the follower's problem is nonlinear and nonconvex. In this game, the leader attacks a limited number of arcs with the goal to maximize the load shed and the follower aims at minimizing the load shed by solving a transport problem in the interdicted network. We develop an exact algorithm consisting of lower and upper bounding schemes that computes an optimal interdiction under the assumption that the interdicted network remains weakly connected. The main challenge consists of computing valid upper bounds for the maximal load shed, whereas lower bounds can directly be derived from the follower's problem. To compute an upper bound, we propose solving a specific bilevel problem, which is derived from restricting the flexibility of the follower when adjusting the load flow. This bilevel problem still has a nonlinear and nonconvex follower's problem, for which we then prove necessary and sufficient optimality conditions. Consequently, we obtain equivalent single-level reformulations of the specific bilevel model to compute upper bounds. Our numerical results show the applicability of this exact approach using the example of gas networks. KW - Interdiction Games KW - Bilevel Optimization KW - Potential-Based Flows KW - Mixed-Integer Nonlinear Optimization Y1 - 2022 ER - TY - INPR A1 - Shyshkanova, Ganna A1 - Walther, Andrea T1 - Contact Pressure over Doubly Connected Rectangular Domains and Punch Shape Optimization N2 - Contact problems arise in a variety of industrial processes, engineering and biomechanical systems. 3-D contact problem for a rigid punch with a doubly connected base bounded by the lines close to rectangles is in consideration. An analytic-numerical technique is developed for its solving. The problem contains Fredholm integral equations of the first kind, which are transformed into the second kind by means of regularization. Using the simple layer potential expansion, the kernels of the integrals are presented in the form of expansions in the powers of the polar radius. The difference between the values of the desired function at different points and the subsequent interpolation of the terms are proposed to smooth the kernels and eliminate singularities. The integral equations are reduced to one-dimension and then solved using quadrature formulas. Subsequently a punch shape is taken as a desired function, and as a minimizing functional is considered the root-mean-square deviation of the pressure distribution arising under the punch from some optimal distribution. In this case, the values of the total forces and moments applied to the punch are assumed to be given, which leads to restrictions imposed on the distributions by the equilibrium conditions. The normal displacements are determined which arising under the action of the found contact pressure on the elastic half-space. The desired punch shape is found using the simple layer potential. A solution to the problem is obtained for the punch with the doubly connected base bounded by lines close to rectangles. KW - contact problem KW - shape optimization KW - simple layer potential KW - integral equations Y1 - 2022 ER - TY - INPR A1 - Kreimeier, Timo A1 - Pokutta, Sebastian A1 - Walther, Andrea A1 - Woodstock, Zev T1 - On a Frank-Wolfe Approach for Abs-smooth Functions N2 - We propose an algorithm which appears to be the first bridge between the fields of conditional gradient methods and abs-smooth optimization. Our nonsmooth nonconvex problem setting is motivated by machine learning, since the broad class of abs-smooth functions includes, for instance, the squared $\ell_2$-error of a neural network with ReLU or hinge Loss activation. To overcome the nonsmoothness in our problem, we propose a generalization to the traditional Frank-Wolfe gap and prove that first-order minimality is achieved when it vanishes. We derive a convergence rate for our algorithm which is identical to the smooth case. Although our algorithm necessitates the solution of a subproblem which is more challenging than the smooth case, we provide an efficient numerical method for its partial solution, and we identify several applications where our approach fully solves the subproblem. Numerical and theoretical convergence is demonstrated, yielding several conjectures. KW - Frank-Wolfe algorithm KW - Active Signature Method KW - abs-smooth functions KW - nonsmooth optimization KW - convergence rate Y1 - 2022 ER - TY - INPR A1 - Grübel, Julia A1 - Krug, Richard A1 - Schmidt, Martin A1 - Wollner, Winnifried T1 - A Successive Linear Relaxation Method for MINLPs with Multivariate Lipschitz Continuous Nonlinearities N2 - We present a novel method for mixed-integer optimization problems with multivariate and Lipschitz continuous nonlinearities. In particular, we do not assume that the nonlinear constraints are explicitly given but that we can only evaluate them and that we know their global Lipschitz constants. The algorithm is a successive linear relaxation method in which we alternate between solving a master problem, which is a mixed-integer linear relaxation of the original problem, and a subproblem, which is designed to tighten the linear relaxation of the next master problem by using the Lipschitz information about the respective functions. By doing so, we follow the ideas of Schmidt et al. (2018, 2021) and improve the tackling of multivariate constraints. Although multivariate nonlinearities obviously increase modeling capabilities, their incorporation also significantly increases the computational burden of the proposed algorithm. We prove the correctness of our method and also derive a worst-case iteration bound. Finally, we show the generality of the addressed problem class and the proposed method by illustrating that both bilevel optimization problems with nonconvex and quadratic lower levels as well as nonlinear and mixed-integer models of gas transport can be tackled by our method. We provide the necessary theory for both applications and briefly illustrate the outcomes of the new method when applied to these two problems. KW - Mixed-Integer Nonlinear Optimization KW - Global Optimization KW - Lipschitz Optimization KW - Bilevel Optimization KW - Gas Networks Y1 - 2022 ER - TY - INPR A1 - Horländer, Andreas A1 - Schmidt, Martin T1 - A Penalty Branch-and-Bound Method for Mixed-Integer Quadratic Bilevel Problems N2 - We propose an algorithm for solving bilevel problems with mixed-integer convex-quadratic upper level as well as convex-quadratic and continuous lower level. The method is based on a classic branch-and-bound procedure, where branching is performed on the integer constraints and on the complementarity constraints resulting from the KKT reformulation of the lower-level problem. However, instead of branching on constraints as usual, suitably chosen penalty terms are added to the objective function in order to create new subproblems in the tree. We prove the correctness of the method and present its applicability by some first numerical results. KW - Bilevel optimization KW - Branch-and-bound KW - Penalty methods KW - Mixed-integer optimization Y1 - 2022 ER - TY - INPR A1 - Aigner, Kevin-Martin A1 - Bärmann, Andreas A1 - Braun, Kristin A1 - Liers, Frauke A1 - Pokutta, Sebastian A1 - Schneider, Oskar A1 - Sharma, Kartikey A1 - Tschuppik, Sebastian T1 - Data-driven Distributionally Robust Optimization over Time N2 - Stochastic Optimization (SO) is a classical approach for optimization under uncertainty that typically requires knowledge about the probability distribution of uncertain parameters. As the latter is often unknown, Distributionally Robust Optimization (DRO) provides a strong alternative that determines the best guaranteed solution over a set of distributions (ambiguity set). In this work, we present an approach for DRO over time that uses online learning and scenario observations arriving as a data stream to learn more about the uncertainty. Our robust solutions adapt over time and reduce the cost of protection with shrinking ambiguity. For various kinds of ambiguity sets, the robust solutions converge to the SO solution. Our algorithm achieves the optimization and learning goals without solving the DRO problem exactly at any step. We also provide a regret bound for the quality of the online strategy which converges at a rate of $ O(\log T / \sqrt{T})$, where $T$ is the number of iterations. Furthermore, we illustrate the effectiveness of our procedure by numerical experiments on mixed-integer optimization instances from popular benchmark libraries and give practical examples stemming from telecommunications and routing. Our algorithm is able to solve the DRO over time problem significantly faster than standard reformulations. KW - distributionally robust optimization KW - learning over time KW - gradient descent Y1 - 2023 ER - TY - INPR A1 - Molan, Ioana A1 - Schmidt, Martin T1 - Using Neural Networks to Solve Linear Bilevel Problems with Unknown Lower Level N2 - Bilevel problems are used to model the interaction between two decision makers in which the lower-level problem, the so-called follower's problem, appears as a constraint in the upper-level problem of the so-called leader. One issue in many practical situations is that the follower's problem is not explicitly known by the leader. For such bilevel problems with unknown lower-level model we propose the use of neural networks to learn the follower's optimal response for given decisions of the leader based on available historical data of pairs of leader and follower decisions. Integrating the resulting neural network in a single-level reformulation of the bilevel problem leads to a challenging model with a black-box constraint. We exploit Lipschitz optimization techniques from the literature to solve this reformulation and illustrate the applicability of the proposed method with some preliminary case studies using academic and linear bilevel instances. KW - Bilevel optimization KW - Unknown follower problems KW - Neural networks KW - Lipschitz optimization Y1 - 2022 ER - TY - INPR A1 - Gutina, Daria A1 - Bärmann, Andreas A1 - Roeder, Georg A1 - Schellenberger, Martin A1 - Liers, Frauke T1 - Optimisation over Decision Trees – A Case Study for the Design of Stable Direct-Current Electricity Networks N2 - In many real-world mixed-integer optimisation problems from engineering, the side constraints can be subdivided into two categories: constraints which describe a certain logic to model a feasible allocation of resources (such as a maximal number of available assets, working time requirements, maintenance requirements, contractual obligations, etc.), and constraints which model physical processes and the related quantities (such as current, pressure, temperature, etc.). While the first type of constraints can often easily be stated in terms of a mixed-integer program (MIP), the second part may involve the incorporation of complex non-linearities, partial differential equations or even a black-box simulation of the involved physical process. In this work, we propose the integration of a trained tree-based classifier – a decision-tree or a random forest, into a mixed-integer optimization model as a possible remedy. We assume that the classifier has been trained on data points produced by a detailed simulation of a given complex process to represent the functional relationship between the involved physical quantities. We then derive MIP-representable reformulations of the trained classifier such that the resulting model can be solved using state-of-the-art solvers. At the hand of several use cases in terms of possible optimisation goals, we show the broad applicability of our framework that is easily extendable to other tasks beyond engineering. In a detailed real-world computational study for the design of stable direct- current power networks, we demonstrate that our approach yields high-quality solutions in reasonable computation times. KW - decision trees KW - random forests KW - mixed-integer programming KW - power networks Y1 - 2022 ER - TY - JOUR A1 - Gugat, Martin A1 - Schuster, Michael A1 - Steffensen, Sonja T1 - A Dynamic Multilevel Model of the European Gas Market N2 - The European gas market is governed by rules that are agreed on by the European Union. We present a mathematical market model that takes into account this structure, where the technical system operator (TSO) offers certain transportation capacities that can be booked and later nominated within the previously chosen bookings. The TSO also fixes booking fees and defines an operational control of the gas pipeline system in order to deliver the gas according to the nominations. Since the gas flow is governed by a system of partial differential equations, to realize this control structure partial differential equations (PDEs) should be involved in the model. While the four level gas market model has been discussed previously, in this paper we take into account the flow model by PDEs in the discussion of the model and in the reduction to a single level problem, where we also state the corresponding necessary optimality conditions. KW - Gas Dynamics KW - Gas Market KW - Nodal Control KW - Isothermal Euler Equations Y1 - 2022 UR - http://cot.mathres.org/archives/1488 IS - Volume 2023 SP - 1 EP - 26 PB - Communications in Optimization Theory ER - TY - THES A1 - Grübel, Julia T1 - Existence, Uniqueness, and Algorithms for Equilibria in Competitive Energy Markets N2 - Due to the transition towards climate neutrality, energy markets are rapidly evolving. New technologies are developed that allow electricity from renewable energy sources to be stored or to be converted into other energy commodities. As a consequence, new players enter the markets and existing players gain more importance. Market equilibrium problems are capable of capturing these changes and therefore enable us to answer contemporary research questions with regard to energy market design and climate policy. This cumulative dissertation is devoted to the study of different market equilibrium problems that address such emerging aspects in liberalized energy markets. In the first part, we review a well-studied competitive equilibrium model for energy commodity markets and extend this model by sector coupling, by temporal coupling, and by a more detailed representation of physical laws and technical requirements. Moreover, we summarize our main contributions of the last years with respect to analyzing the market equilibria of the resulting equilibrium problems. For the extension regarding sector coupling, we derive sufficient conditions for ensuring uniqueness of the short-run equilibrium a priori and for verifying uniqueness of the long-run equilibrium a posteriori. Furthermore, we present illustrative examples that each of the derived conditions is indeed necessary to guarantee uniqueness in general. For the extension regarding temporal coupling, we provide sufficient conditions for ensuring uniqueness of demand and production a priori. These conditions also imply uniqueness of the short-run equilibrium in case of a single storage operator. However, in case of multiple storage operators, examples illustrate that charging and discharging decisions are not unique in general. We conclude the equilibrium analysis with an a posteriori criterion for verifying uniqueness of a given short-run equilibrium. Since the computation of equilibria is much more challenging due to the temporal coupling, we shortly review why a tailored parallel and distributed alternating direction method of multipliers enables to efficiently compute market equilibria. For the extension regarding physical laws and technical requirements, we show that, in nonconvex settings, existence of an equilibrium is not guaranteed and that the fundamental welfare theorems therefore fail to hold. In addition, we argue that the welfare theorems can be re-established in a market design in which the system operator is committed to a welfare objective. For the case of a profit-maximizing system operator, we propose an algorithm that indicates existence of an equilibrium and that computes an equilibrium in the case of existence. Based on well-known instances from the literature on the gas and electricity sector, we demonstrate the broad applicability of our algorithm. Our computational results suggest that an equilibrium often exists for an application involving nonconvex but continuous stationary gas physics. In turn, integralities introduced due to the switchability of DC lines in DC electricity networks lead to many instances without an equilibrium. Finally, we state sufficient conditions under which the gas application has a unique equilibrium and the line switching application has finitely many. In the second part, all preprints belonging to this cumulative dissertation are provided. These preprints, as well as two journal articles to which the author of this thesis contributed, are referenced within the extended summary in the first part and contain more details. KW - Energy markets KW - Equilibrium computation KW - Existence KW - Perfect competition KW - Uniqueness Y1 - U6 - https://doi.org/10.25353/ubtr-xxxx-4c49-7f53 ER - TY - INPR A1 - Beck, Yasmine A1 - Ljubic, Ivana A1 - Schmidt, Martin T1 - A Survey on Bilevel Optimization Under Uncertainty N2 - Bilevel optimization is a very active field of applied mathematics. The main reason is that bilevel optimization problems can serve as a powerful tool for modeling hierarchical decision making processes. This ability, however, also makes the resulting problems challenging to solve - both in theory and practice. Fortunately, there have been significant algorithmic advances in the field of bilevel optimization so that we can solve much larger and also more complicated problems today compared to what was possible to solve two decades ago. This results in more and more challenging bilevel problems that researchers try to solve today. This survey gives a detailed overview of one of these more challenging classes of bilevel problems: bilevel optimization under uncertainty. We review the classic ways of addressing uncertainties in bilevel optimization using stochastic or robust techniques. Moreover, we highlight that the sources of uncertainty in bilevel optimization are much richer than for usual, i.e., single-level, problems since not only the problem's data can be uncertain but also the (observation of the) decisions of the two players can be subject to uncertainty. We thus also review the field of bilevel optimization under limited observability, the area of problems considering only near-optimal decisions, and discuss intermediate solution concepts between the optimistic and pessimistic cases. Finally, we also review the rich literature on applications studied using uncertain bilevel problems such as in energy, for interdiction games and security applications, in management sciences, and networks. KW - Bilevel optimization KW - Optimization under uncertainty KW - Bounded rationality KW - Survey Y1 - 2022 ER - TY - JOUR A1 - Hernandez, Martin A1 - Lecaros, Rodrigo A1 - Zamorano, Sebastian T1 - Averaged turnpike property for differential equations with random constant coefficients N2 - This paper studies the integral turnpike and turnpike in average for a class of random or- dinary differential equations. We prove that, under suitable assumptions on the matrices that define the system, the optimal solutions for an optimal distributed control tracking problem remain, in an averaged sense, sufficiently close to the associated random stationary optimal solution for the majority of the time horizon Y1 - 2022 ER - TY - CHAP A1 - Leugering, Günter T1 - Nonoverlapping Domain Decomposition in Space and Time for Optimal Control Problems on Metric Graphs by the Example of Gas Flow in Pipe Networks N2 - We consider non-overlapping domain decomposition methods for ordinary and partial differential equations and corresponding optimal control problems on metric graphs. As an exemplary context, we chose a semilinear approximation of the Euler system and a doubly nonlinear parabolic model that has come to be known as friction dominated flow in gas pipe networks. By this choice, we encounter hyperbolic, parabolic and elliptic linear and nonlinear problems in a single highly motivating application. We depart from the classical domain decomposition methods described by P.L. Lions and J.L. Lions and O. Pironneau and extend those to problems on metric graphs. The choice of methods is determined by the desire to use a control concept that has come to be known as virtual controls which, in turn, possibly lead to a fully parallel decomposition of the corresponding optimality systems. In a second step, we extend the methods to p-Laplace problems on networks. The analysis, due to space limitations, will appear in a forthcoming publication. See however J.E. Lagnese and G. Leugering. Furthermore, we then describe methods for space and time domain decomposition or optimal control problems in the spirit of J.E. Lagnese and G. Leugering. We finally provide some comments on PINN-based approximations of the methods described before. We provide numerical evidence for all algorithms discussed. KW - Optimal control KW - PDEs on graphs KW - p-Laplace problem on a graph KW - p-parabolic problems KW - instantaneous control Y1 - 2022 ER - TY - INPR A1 - Beck, Yasmine A1 - Schmidt, Martin A1 - Thürauf, Johannes A1 - Bienstock, Daniel T1 - On a Computationally Ill-Behaved Bilevel Problem with a Continuous and Nonconvex Lower Level N2 - It is well known that bilevel optimization problems are hard to solve both in theory and practice. In this paper, we highlight a further computational difficulty when it comes to solving bilevel problems with continuous but nonconvex lower levels. Even if the lower-level problem is solved to ɛ-feasibility regarding its nonlinear constraints for an arbitrarily small but positive ɛ, the obtained bilevel solution as well as its objective value may be arbitrarily far away from the actual bilevel solution and its actual objective value. This result even holds for bilevel problems for which the nonconvex lower level is uniquely solvable, for which the strict complementarity condition holds, for which the feasible set is convex, and for which Slater's constraint qualification is satisfied for all feasible upper-level decisions. Since the consideration of ɛ-feasibility cannot be avoided when solving nonconvex problems to global optimality, our result shows that computational bilevel optimization with continuous and nonconvex lower levels needs to be done with great care. Finally, we illustrate that the nonlinearities in the lower level are the key reason for the observed bad behavior by showing that linear bilevel problems behave much better at least on the level of feasible solutions. KW - Bilevel optimization KW - Nonconvex lower levels KW - Approximate feasibility KW - Global optimization Y1 - 2022 ER - TY - INPR A1 - Hannes, Dänschel A1 - Volker, Mehrmann A1 - Roland, Marius A1 - Schmidt, Martin T1 - Adaptive Nonlinear Optimization of District Heating Networks Based on Model and Discretization Catalogs N2 - We propose an adaptive optimization algorithm for operating district heating networks in a stationary regime. The behavior of hot water flow in the pipe network is modeled using the incompressible Euler equations and a suitably chosen energy equation. By applying different simplifications to these equations, we derive a catalog of models. Our algorithm is based on this catalog and adaptively controls where in the network which model is used. Moreover, the granularity of the applied discretization is controlled in a similar adaptive manner. By doing so, we are able to obtain optimal solutions at low computational costs that satisfy a prescribed tolerance w.r.t. the most accurate modeling level. To adaptively control the switching between different levels and the adaptation of the discretization grids, we derive error measure formulas and a posteriori error measure estimators. Under reasonable assumptions we prove that the adaptive algorithm terminates after finitely many iterations. Our numerical results show that the algorithm is able to produce solutions for problem instances that have not been solvable before. KW - District heating networks KW - Adaptive methods KW - Nonlinear optimization Y1 - 2022 ER -