TY - JOUR A1 - Krug, Richard A1 - Leugering, Günter A1 - Martin, Alexander A1 - Schmidt, Martin A1 - Weninger, Dieter T1 - Time-Domain Decomposition for Optimal Control Problems Governed by Semilinear Hyperbolic Systems with Mixed Two-Point Boundary Conditions JF - Control and Cybernetics N2 - In this article, we continue our work (Krug et al., 2021) on time-domain decomposition of optimal control problems for systems of semilinear hyperbolic equations in that we now consider mixed two-point boundary value problems and provide an in-depth well-posedness analysis. The more general boundary conditions significantly enlarge the scope of applications, e.g., to hyperbolic problems on metric graphs with cycles. We design an iterative method based on the optimality systems that can be interpreted as a decomposition method for the original optimal control problem into virtual control problems on smaller time domains. KW - Time-domain decomposition KW - Optimal control KW - Semilinear hyperbolic systems KW - Convergence Y1 - 2021 ER - TY - JOUR A1 - Krug, Richard A1 - Leugering, Günter A1 - Martin, Alexander A1 - Schmidt, Martin A1 - Weninger, Dieter T1 - Time-Domain Decomposition for Optimal Control Problems Governed by Semilinear Hyperbolic Systems JF - SIAM Journal on Control and Optimization N2 - In this article, we extend the time-domain decomposition method described by Lagnese and Leugering (2003) to semilinear optimal control problems for hyperbolic balance laws with spatio-temporal varying coefficients. We provide the design of the iterative method applied to the global first-order optimality system, prove its convergence, and derive an a posteriori error estimate. The analysis is done entirely on the continuous level. A distinguishing feature of the method is that the decomposed optimality system can be interpreted as an optimality system of a local "virtual" optimal control problem. Thus, the iterative time-domain decomposition of the optimality system can be interpreted as an iterative parallel scheme for virtual optimal control problems on the subintervals. A typical example and further comments are given to show the range of potential applications. Moreover, we provide some numerical experiments to give a first interpretation of the role of the parameters involved in the iterative process. KW - Time-domain decomposition KW - Optimal control KW - Semilinear hyperbolic systems KW - Convergence KW - A posteriori error estimates Y1 - 2020 ER - TY - INPR A1 - Leugering, Günter T1 - Space-time-domain decomposition for optimal control problems governed by linear hyperbolic systems N2 - In this article, we combine a domain decomposition method in space and time for optimal control problems with PDE-constraints described by Lagnese and Leugering to a simultaneous space-time decomposition applied to optimal control problems for systems of linear hyperbolic equations with distributed control. We thereby extend the recent work by Krug et al. and answer a long standing open question as to whether the combination of time- and space domain decomposition for the method under consideration can be put into one single convergent iteration procedure. The algorithm is designed for a semi-elliptic system of equations obtained from the hyperbolic optimality system by the way of reduction to the adjoint state. The focus is on the relation to the classical procedure introduced by Lions for elliptic problems. KW - Space- and time-domain decomposition KW - Optimal control KW - linear hyperbolic systems KW - Convergence KW - A posteriori error estimates Y1 - 2021 ER - TY - INPR A1 - Leugering, Günter T1 - Space-Time-Domain Decomposition for Optimal Control Problems Governed by Linear Hyperbolic Systems N2 - In this article, we combine a domain decomposition method in space and time for optimal control problems with PDE-constraints described by Lagnese and Leugering to a simultaneous space-time decomposition applied to optimal control problems for systems of linear hyperbolic equations with distributed control. We thereby extend the recent work by Krug et al. and answer a long standing open question as to whether the combination of time- and space domain decomposition for the method under consideration can be put into one single convergent iteration procedure. The algorithm is designed for a semi-elliptic system of equations obtained from the hyperbolic optimality system by the way of reduction to the adjoint state. The focus is on the relation to the classical procedure introduced by Lions for elliptic problems. KW - Space- and time-domain decomposition KW - Optimal control KW - linear hyperbolic systems KW - Convergence KW - A posteriori error estimates Y1 - 2021 ER - TY - JOUR A1 - Leugering, Günter T1 - Industrial Applications of Optimal Control for Partial Differential Equations on Networks: Reduction and Decomposition Methods Applied to the Discrete–Continuous Control of Gas Flow in Complex Pipe Systems N2 - This chapter provides an exemplary road map—in a nutshell—from a given industrial application, the control of gas networks, which is far too complex for a direct approach, to a problem that can be actually handled using well-known methods in control theory. It also provides an iterative non-overlapping domain decomposition that can be interpreted as an Uzawa method. The chapter outline two strategies. The first one can be seen as a Jacobi-type approach. In the second approach, fix the integer controls s and decompose the corresponding optimality system for the entire graph into the subgraphs Gk by a another, but very similar, non-overlapping domain decomposition. The problem is the intrinsic coupling of integer controls, continuous controls, and nonlinear dynamics on a metric graph. The idea is to introduce a virtual control that aims at controlling classical in homogeneous Neumann condition including the iteration history at the interface as inhomogeneity to the Robin-type condition that appears in the decomposition. Y1 - 2021 VL - Computational Science and Its Applications IS - 1st Edition SP - 25 EP - 40 ER - TY - JOUR A1 - Gugat, Martin A1 - Leugering, Günter A1 - Martin, Alexander A1 - Schmidt, Martin A1 - Sirvent, Mathias A1 - Wintergerst, David T1 - MIP-Based Instantaneous Control of Mixed-Integer PDE-Constrained Gas Transport Problems JF - Computational Optimization and Applications N2 - We study the transient optimization of gas transport networks including both discrete controls due to switching of controllable elements and nonlinear fluid dynamics described by the system of isothermal Euler equations, which are partial differential equations in time and 1-dimensional space. This combination leads to mixed-integer optimization problems subject to nonlinear hyperbolic partial differential equations on a graph. We propose an instantaneous control approach in which suitable Euler discretizations yield systems of ordinary differential equations on a graph. This networked system of ordinary differential equations is shown to be well-posed and affine-linear solutions of these systems are derived analytically. As a consequence, finite-dimensional mixed-integer linear optimization problems are obtained for every time step that can be solved to global optimality using general-purpose solvers. We illustrate our approach in practice by presenting numerical results on a realistic gas transport network. KW - Mixed-integer optimal control KW - Instantaneous control KW - Partial differential equations on graphs KW - Gas networks KW - Mixed-integer linear optimization Y1 - 2017 U6 - https://doi.org/10.1007/s10589-017-9970-1 VL - 70 IS - 1 SP - 267 EP - 294 ER - TY - JOUR A1 - 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 - 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 - JOUR A1 - Gugat, Martin A1 - Keimer, Alexander A1 - Leugering, Günter A1 - Wang, Zhiqiang ED - Piccoli, Benedetto T1 - Analysis of a system of nonlocal conservation laws for multi-commodity flow on networks JF -  Networks and Heterogeneous Media N2 - We consider a system of scalar nonlocal conservation laws on networks that model a highly re-entrant multi-commodity manufacturing system as encountered in semiconductor production. Every single commodity is mod-eled by a nonlocal conservation law, and the corresponding PDEs are coupled via a collective load, the work in progress. We illustrate the dynamics for two commodities. In the applications, directed acyclic networks naturally occur, therefore this type of networks is considered. On every edge of the network we have a system of coupled conservation laws with nonlocal velocity. At the junctions the right hand side boundary data of the foregoing edges is passed as left hand side boundary data to the following edges and PDEs. For distributing junctions, where we have more than one outgoing edge, we impose time dependent distribution functions that guarantee conservation of mass. We provide results of regularity, existence and well-posedness of the multi-commodity network model for L p-, BV-and W 1,p-data. Moreover, we define an L 2-tracking type objective and show the existence of minimizers that solve the corresponding optimal control problem. KW - conservation laws on network KW - nonlocal conservation laws KW - optimal nodal control KW - systems of hyperbolic pdes Y1 - 2016 U6 - https://doi.org/DOI: 10.3934/nhm.2015.10.749 VL - 10 IS - 4 SP - 749 EP - 785 ER - TY - JOUR A1 - Gugat, Martin A1 - Leugering, Günter ED - Zuazua, Enrique T1 - Time delay in optimal control loops for wave equations JF - ESAIM: COCV N2 - In optimal control loops delays can occur, for example through transmission via digital communication channels. Such delays influence the state that is generated by the implemented control. We study the effect of a delay in the implementation of L 2-norm minimal Neumann boundary controls for the wave equation. The optimal controls are computed as solutions of problems of exact optimal control, that is if they are implemented without delay, they steer the system to a position of rest in a given finite time T. We show that arbitrarily small delays δ > 0 can have a destabilizing effect in the sense that we can find initial states such that if the optimal control u is implemented in the form yx(t, 1) = u(t − δ) for t > δ, the energy of the system state at the terminal time T is almost twice as big as the initial energy. We also show that for more regular initial states, the effect of a delay in the implementation of the optimal control is bounded above in the sense that for initial positions with derivatives of BV-regularity and initial velocities with BV-regularity, the terminal energy is bounded above by the delay δ multiplied with a factor that depends on the BV-norm of the initial data. We show that for more general hyperbolic optimal exact control problems the situation is similar. For systems that have arbitrarily large eigenvalues, we can find terminal times T and arbitrarily small time delays δ, such that at the time T + δ, in the optimal control loop with delay the norm of the state is twice as large as the corresponding norm for the initial state. Moreover, if the initial state satisfies an additional regularity condition, there is an upper bound for the effect of time delay of the order of the delay with a constant that depends on the initial state only. KW - PDE constrained optimization KW - delay KW - wave equation KW - boundary control KW - hyperbolic system Y1 - 2016 U6 - https://doi.org/http://dx.doi.org/10.1051/cocv/2015038 ER - TY - 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 - Hante, Falk A1 - Leugering, Günter A1 - Martin, Alexander A1 - Schewe, Lars A1 - Schmidt, Martin T1 - Challenges in optimal control problems for gas and fluid flow in networks of pipes and canals: From modeling to industrial applications N2 - We consider optimal control problems for the flow of gas or fresh water in pipe networks as well as drainage or sewer systems in open canals. The equations of motion are taken to be represented by the nonlinear isothermal Euler gas equations, the water hammer equations, or the St.~Venant equations for flow. We formulate model hierarchies and derive an abstract model for such network flow problems including pipes, junctions, and controllable elements such as valves, weirs, pumps, as well as compressors. We use the abstract model to give an overview of the known results and challenges concerning equilibria, well-posedness, controllability, and optimal control. A major challenge concerning the optimization is to deal with switching on-off states that are inherent to controllable devices in such applications combined with continuous simulation and optimization of the gas flow. We formulate the corresponding mixed-integer nonlinear optimal control problems and outline a decomposition approach as a solution technique. KW - Networks KW - pipes KW - optimal control KW - Euler and St. Venant equations KW - hierarchy of models Y1 - 2016 ER - TY - 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 - Leugering, Günter A1 - Martin, Alexander A1 - Schmidt, Martin A1 - Sirvent, Mathias T1 - Nonoverlapping Domain Decomposition for Optimal Control Problems governed by Semilinear Models for Gas Flow in Networks N2 - We consider optimal control problems for gas flow in pipeline networks. The equations of motion are taken to be represented by a first-order system of hyperbolic semilinear equations derived from the fully nonlinear isothermal Euler gas equations. We formulate an optimal control problem on a network and introduce a tailored time discretization thereof. In order to further reduce the complexity, we consider an instantaneous control strategy. The main part of the paper is concerned with a nonoverlapping domain decomposition of the optimal control problem on the graph into local problems on smaller sub-graphs - ultimately on single edges. We prove convergence of the domain decomposition method on networks and study the wellposedness of the corresponding time-discrete optimal control problems. The point of the paper is that we establish virtual control problems on the decomposed subgraphs such that the corresponding optimality systems are in fact equal to the systems obtained via the domain decomposition of the entire optimality system. KW - Optimal control, Gas networks, Euler's equation, Semilinear PDE, Nonoverlapping domain decomposition Y1 - 2017 VL - 46 IS - 3 SP - 191 EP - 225 PB - Control and Cybernetics ER - TY - JOUR A1 - Gugat, Martin A1 - Leugering, Günter A1 - Martin, Alexander A1 - Schmidt, Martin A1 - Sirvent, Mathias A1 - Wintergerst, David T1 - Towards Simulation Based Mixed-Integer Optimization with Differential Equations JF - Networks N2 - We propose a decomposition based method for solving mixed-integer nonlinear optimization problems with “black-box” nonlinearities, where the latter, e.g., may arise due to differential equations or expensive simulation runs. The method alternatingly solves a mixed-integer linear master problem and a separation problem for iteratively refining the mixed-integer linear relaxation of the nonlinear equalities. The latter yield nonconvex feasible sets for the optimization model but we have to restrict ourselves to convex and monotone constraint functions. Under these assumptions, we prove that our algorithm finitely terminates with an approximate feasible global optimal solution of the mixed integer nonlinear problem. Additionally, we show the applicability of our approach for three applications from optimal control with integer variables, from the field of pressurized flows in pipes with elastic walls, and from steady-state gas transport. For the latter we also present promising numerical results of our method applied to real-world instances that particularly show the effectiveness of our method for problems defined on networks. KW - Mixed-Integer Optimization KW - Simulation Based Optimization KW - Optimization with Differential Equations KW - Decomposition Method KW - Gas Transport Networks Y1 - 2018 U6 - https://doi.org/10.1002/net.21812 ER - TY - JOUR A1 - Gugat, Martin A1 - Leugering, Günter A1 - Hante, Falk ED - Piccoli, Benedetto T1 - Stationary States in Gas Networks JF - Networks and Heterogeneous Media N2 - Pipeline networks for gas transportation often contain circles. For such networks it is more difficult to determine the stationary states than for networks without circles. We present a method that allows to compute the stationary states for subsonic pipe flow governed by the isothermal Euler equations for certain pipeline networks that contain circles. We also show that suitably chosen boundary data determine the stationary states uniquely. The construction is based upon novel explicit representations of the stationary states on single pipes for the cases with zero slope and with nonzero slope. In the case with zero slope, the state can be represented using the Lambert-W function. KW - Network Y1 - 2016 U6 - https://doi.org/doi:10.3934/nhm.2015.10.295 VL - 10 IS - 2 SP - 295 EP - 320 ER - TY - INPR A1 - Leugering, Günter T1 - Nonoverlapping Domain Decomposition for Instantaneous Optimal Control of Friction Dominated Flow in a Gas-Network N2 - We consider a non-overlapping domain decomposition method for an optimal control problem related to the flow of gas in a pipe network. The equations of motions are taken to be represented by a friction dominated model derived from a semi-linear approximation of the fully nonlinear isothermal Euler gas equations. This involves a p-Laplace-type problem on the graph with p = 3/2. We continue the work by Leugering and Mophou where such a problem has been discussed in the context of an instantaneous control strategy. We provide a non-overlapping domain decomposition in the spirit of P.L. Lions for elliptic problems and extend the method to the first order optimality system. KW - Optimal control KW - Gas networks KW - p-Laplace problem on a graph KW - Optimality system KW - Domain decomposition 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 - 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 -