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 - 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 - INPR A1 - Gugat, Martin A1 - Krug, Richard A1 - Martin, Alexander T1 - Transient gas pipeline flow: Analytical examples, numerical simulation and a comparison to the quasi-static approach N2 - The operation of gas pipeline flow with high pressure and small Mach numbers allows to model the flow by a semilinear hyperbolic system of partial differential equations. In this paper we present a number of transient and stationary analytical solutions of this model. They are used to discuss and clarify why a pde model is necessary to handle certain dynamic situations in the operation of gas transportation networks. We show that adequate numerical discretizations can capture the dynamical behavior sufficiently accurate. We also present examples that show that in certain cases an optimization approach that is based upon multi-period optimization of steady states does not lead to approximations that converge to the optimal state. Y1 - 2021 U6 - https://doi.org/10.1007/s11081-021-09690-4 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 - Gugat, Martin A1 - Habermann, Jens A1 - Hintermüller, Michael A1 - Huber, Olivier T1 - Constrained exact boundary controllability of a semilinear model for pipeline gas flow N2 - While the quasilinear isothermal Euler equations are an excellent model for gas pipeline flow, the operation of the pipeline flow with high pressure and small Mach numbers allows us to obtain approximate solutions by a simpler semilinear model. We provide a derivation of the semilinear model that shows that the semilinear model is valid for sufficiently low Mach numbers and sufficiently high pressures. We prove an existence result for continuous solutions of the semilinear model that takes into account lower and upper bounds for the pressure and an upper bound for the magnitude of the Mach number of the gas flow. These state constraints are important both in the operation of gas pipelines and to guarantee that the solution remains in the set where the model is physically valid. We show the constrained exact boundary controllability of the system with the same pressure and Mach number constraints. Y1 - 2021 ER - TY - JOUR A1 - Martin, Gugat A1 - Giesselmann, Jan A1 - Kunkel, Teresa T1 - Exponential synchronization of a nodal observer for a semilinear model for the flow in gas networks N2 - The flow of gas through networks of pipes can be modeled by coupling hyperbolic systems of partial differential equations that describe the flow through the pipes that form the edges of the graph of the network by algebraic node conditions that model the flow through the vertices of the graph. In the network, measurements of the state are available at certain points in space.Based upon these nodal observations, the complete system state can be approximated using an observer system. In this paper we present a nodal observer, and prove that the state of the observer system converges to the original state exponentially fast. Numerical experiments confirm the theoretical findings. Y1 - 2021 U6 - https://doi.org/10.1093/imamci/dnab029 CY - IMA Journal of Mathematical Control and Information ER - TY - JOUR A1 - Gugat, Martin A1 - Hante, Falk A1 - Jin, Li T1 - Closed loop control of gas flow in a pipe: Stability for a transient model JF - at - Automatisierungstechnik N2 - This contribution focuses on the analysis and control of friction-dominated flow of gas in pipes. The pressure in the gas flow is governed by a partial differential equation that is a doubly nonlinear parabolic equation of p-Laplace type, where p=2/3. Such equations exhibit positive solutions, finite speed of propagation and satisfy a maximum principle. The pressure is fixed on one end (upstream), and the flow is specified on the other end (downstream). These boundary conditions determine a unique steady equilibrium flow. We present a boundary feedback flow control scheme, that ensures local exponential stability of the equilibrium in an L2-sense. The analysis is done both for the pde system and an ode system that is obtained by a suitable spatial semi-discretization. The proofs are based upon suitably chosen Lyapunov functions. Y1 - 2020 ER - TY - INPR A1 - Gugat, Martin A1 - Herty, Michael T1 - Modeling, Control and Numerics of Gas Networks N2 - In this article we survey recent progress on mathematical results on gas flow in pipe networks with a special focus on questions of control and stabilization. We briefly present the modeling of gas flow and coupling conditions for flow through vertices of a network. Our main focus is on gas models for spatially one-dimensional flow governed by hyperbolic balance laws. We survey results on classical solutions as well as weak solutions. We present results on well–posedness, controllability, feedback stabilization, the inclusion of uncertainty in the models and numerical methods. KW - Hyperbolic Balance Laws, Stabilization, Exact Controllability, Modeling of Gas Flow, Finite-Volume Schemes, Optimal control, Uncertainty Y1 - 2020 ER - TY - JOUR A1 - Schuster, Michael A1 - Strauch, Elisa A1 - Gugat, Martin A1 - Lang, Jens T1 - Probabilistic Constrained Optimization on Flow Networks N2 - Uncertainty often plays an important role in dynamic flow problems. In this paper, we consider both, a stationary and a dynamic flow model with uncertain boundary data on networks. We introduce two different ways how to compute the probability for random boundary data to be feasible, discussing their advantages and disadvantages. In this context, feasible means, that the flow corresponding to the random boundary data meets some box constraints at the network junctions. The first method is the spheric radial decomposition and the second method is a kernel density estimation. In both settings, we consider certain optimization problems and we compute derivatives of the probabilistic constraint using the kernel density estimator. Moreover, we derive necessary optimality conditions for the stationary and the dynamic case. Throughout the paper, we use numerical examples to illustrate our results by comparing them with a classical Monte Carlo approach to compute the desired probability. KW - Probabilistic Constraints KW - Flow Networks KW - Gas Networks KW - Spheric Radial Decomposition KW - Kernel Density Estimator Y1 - 2020 U6 - https://doi.org/https://doi.org/10.1007/s11081-021-09619-x VL - Optimization and Engineering ER - TY - JOUR A1 - Gugat, Martin A1 - Schuster, Michael A1 - Zuazua, Enrique T1 - M. Gugat, M. Schuster, E. Zuazua. The Finite-Time Turnpike Phenomenon for Optimal Control Problems: Stabilization by Non-Smooth Tracking Terms, in “Stabilization of Distributed Parameter Systems: Design Methods and Applications”. Grigory Sklyar Alexander Zuyev Eds., ICIAM 2019 SEMA SIMAI Springer Series 2, p. 17-42. ISSN 2199-3041 N2 - In this paper, problems of optimal control are considered where in the objective function, in addition to the control cost, there is a tracking term that measures the distance to a desired stationary state. The tracking term is given by some norm, and therefore it is in general not differentiable. In the optimal control problem, the initial state is prescribed. We assume that the system is either exactly controllable in the classical sense or nodal profile controllable. We show that both for systems that are governed by ordinary differential equations and for infinite-dimensional systems, for example, for boundary control systems governed by the wave equation, under certain assumptions, the optimal system state is steered exactly to the desired state after finite time. Y1 - 2020 VL - SEMA SIMAI Springer Series 2 SP - 17 EP - 42 PB - Springer International Publishing ET - Grigory Sklyar Alexander Zuyev Eds., ICIAM 2019 ER - TY - JOUR A1 - Gugat, Martin A1 - Giesselmann, Jan T1 - Boundary feedback stabilization of a semilinear model for the flow in star-shaped gas networks N2 - The flow of gas through a pipeline network can be modelled by a coupled system of 1-d quasilinear hyperbolic equations. In this system, the influence of certain source terms that model friction effects is essential. Often for the solution of control problems it is convenient to replace the quasilinear model by a simpler semilinear model. In this paper, we analyze the behavior of such a semilinear model on a star-shaped network. The model is derived from the diagonal form of the quasilinear model by replacing the eigenvalues by the sound speed multiplied by 1 or -1 respectively. Thus in the corresponding eigenvalues the influence of the gas velocity is neglected, which is justified in the applications since it is much smaller than the sound speed in the gas. For a star-shaped network of horizontal pipes for suitable coupling conditions we present boundary feedback laws that stabilize the system state exponentially fast to a position of rest for sufficiently small initial data. We show the exponential decay of the $H^1$-norm for arbitrarily long pipes. This is remarkable since in general even for linear systems, for certain source terms the system can become exponentially unstable if the space interval is too long. Our proofs are based upon observability inequalities for the $L^2$ and the $H^1$-norm. Y1 - 2020 U6 - https://doi.org/10.1051/cocv/2021061 CY - ESAIM:COCV ER - TY - INPR A1 - Gugat, Martin A1 - Herty, Michael T1 - Limits of stabilizabilizy for a semilinear model for gas pipeline flow N2 - We present a positive and a negative stabilization result for a semilinear model of gas flow in pipelines. For feedback boundary conditions we obtain an unconditional stabilization result in the absence and conditional instability in the presence of the source term. We also obtain unconditional instability for the corresponding quasilinear model given by the isothermal Euler equations Y1 - 2020 ER - TY - INPR A1 - Gugat, Martin A1 - Sokolowski, Jan T1 - On Problems of Dynamic Optimal Nodal control for Gas Networks N2 - We consider a dynamic ptimal control problem for gas pipeline systems. The flow is governed by a quasilinear hyperbolic model. Since in the operation of the gas networks regular solutions without shocks are desirable, we impose appropriate state and control constraint in order to guarantee that a classical solution is generated. Due to a W^{2;inf}-regularization term in the objective function, we can show the existence of an optimal control. Moreover, we give conditions that guarantee that the control becomes constant a the end of the control time interval if the weight of the regularization term is suffciently large. KW - optimal nodal control KW - gas network KW - turnpike property KW - quasilinear hyperbolic problem KW - dynamic control Y1 - 2021 ER - TY - JOUR A1 - Gugat, Martin T1 - On the turnpike property with interior decay for optimal control problems JF - Mathematics of Control, Signals, and Systems N2 - In this paper the turnpike phenomenon is studied for problems of optimal control where both pointwise-in-time state and control constraints can appear. We assume that in the objective function, a tracking term appears that is given as an integral over the time-interval [0, T] and measures the distance to a desired stationary state. In the optimal control problem, both the initial and the desired terminal state are prescribed. We assume that the system is exactly controllable in an abstract sense if the time horizon is long enough. We show that that the corresponding optimal control problems on the time intervals [0, T] give rise to a turnpike structure in the sense that for natural numbers n if T is su� ciently large, the contribution of the objective function from subintervals of [0, T] of the form [t - t/2^n, t + (T-t)/2^n] is of the order 1/min{t^n, (T-t)^n}. We also show that a similar result holds for epsilon-optimal solutions of the optimal control problems if epsilon > 0 is chosen suffciently small. At the end of the paper we present both systems that are governed by ordinary differential equations and systems governed by partial differential equations where the results can be applied. KW - Optimal control KW - Turnpike KW - Control Constraint KW - State constraint KW - Exact controllability Y1 - 2021 U6 - https://doi.org/https://doi.org/10.1007/s00498-021-00280-4 ER - TY - INPR A1 - Gugat, Martin T1 - Optimal boundary control of the wave equation: The finite-time turnpike phenomenon N2 - It is well-known that vibrating strings can be steered to a position of rest in finite time by suitably defined boundary control functions, if the time horizon is suffciently long. In optimal control problems, the desired terminal state is often enforced by terminal conditions, that add an additional diffculty to the optimal control problem. In this paper we present an optimal control problem for the wave equation with a time-dependent weight in the objective function such that for a suffciently long time horizon, the optimal state reaches a position of rest in finite time without prescribing a terminal constraint. This situation can be seen as a realization of the finite-time turnpike phenomenon that has been studied recently in [1]. Y1 - 2021 ER - TY - JOUR A1 - Gugat, Martin A1 - Ulbrich, Stefan T1 - On Lipschitz Solutions of Initial Boundary Value Problems for Balance Laws JF - Mathematical Models and Methods in Applied Sciences N2 - The flow of gas through networks of pipes can be modeled by the isothermal Euler equations and algebraic node conditions that model the flow through the vertices of the network graph. We prove the well-posedness of the system for gas with nonconstant compressibility factor that is given by an affine linear function. We consider initial data and control functions that are Lipschitz continuous and compatible with the node and boundary conditions. We show the existence of semi--global Lipschitz continuous solutions of the initial boundary value problem. The construction of the solution is based upon a fixed point iteration along the characteristic curves. The solutions of the intial boundary value problem on arbitrary networks satisfy a maximum principle in terms of the Riemann invariants that states that the maximum of the absolute values is attained for the initial or the boundary data. Y1 - 2017 ER - TY - JOUR A1 - Gugat, Martin A1 - Leugering, Günter A1 - Wang, Ke ED - Zhang, Xu T1 - Neumann boundary feedback stabilization for a nonlinear wave equation: A strict H2-Lyapunov function JF - Mathematical Control and Related Fields (MCRF) N2 - For a system that is governed by the isothermal Euler equations with friction for ideal gas, the corresponding field of characteristic curves is determined by the velocity of the flow. This velocity is determined by a second-order quasilinear hyperbolic equation. For the corresponding initial-boundary value problem with Neumann-boundary feedback, we consider non-stationary solutions locally around a stationary state on a finite time interval and discuss the well-posedness of this kind of problem. We introduce a strict H2-Lyapunov function and show that the boundary feedback constant can be chosen such that the H2-Lyapunov function and hence also the H2-norm of the difference between the non-stationary and the stationary state decays exponentially with time. KW - Boundary feedback control, feedback stabilization, exponential stability, isothermal Euler equations, second-order quasilinear equation, Lyapunov function, stationary state, non-stationary state, gas pipeline. Y1 - 2017 U6 - https://doi.org/10.3934/mcrf.2017015 VL - 7 IS - 3 SP - 419 EP - 448 ER - TY - JOUR A1 - Martin, Gugat A1 - Herty, Michael A1 - Müller, Siegfried ED - Piccoli, Benedetto T1 - Coupling conditions for the transition from supersonic to subsonic fluid states JF - Networks and Heterogeneous Media (NHM) N2 - We discuss coupling conditions for the p-system in case of a transition from supersonic states to subsonic states. A single junction with adjacent pipes is considered where on each pipe the gas ow is governed by a general p-system. By extending the notion of demand and supply known from traffic fiow analysis we obtain a constructive existence result of solutions compatible with the introduced conditions. KW - p-system KW - coupling conditions KW - network Y1 - 2017 U6 - https://doi.org/10.3934/nhm.2017016 VL - 12 IS - 3 SP - 371 EP - 380 ER - TY - JOUR A1 - Gugat, Martin A1 - Zuazua, Enrique ED - Grimble, Mike J. T1 - Exact penalization of terminal constraints for optimal control problems JF - OPTIMAL CONTROL APPLICATIONS AND METHODS N2 - We study optimal control problems for linear systems with prescribed initial and terminal states. We analyze the exact penalization of the terminal constraints. We show that for systems that are exactly controllable, the norm-minimal exact control can be computed as the solution of an optimization problem without terminal constraint but with a nonsmooth penalization of the end conditions in the objective function, if the penalty parameter is sufficiently large. We describe the application of the method for hyperbolic and parabolic systems of partial differential equations, considering the wave and heat equations as particular examples. Copyright © 2016 John Wiley & Sons, Ltd. Y1 - 2016 U6 - https://doi.org/10.1002/oca.2238 VL - 37 IS - 6 SP - 1329 EP - 1354 ER - TY - JOUR A1 - Gugat, Martin A1 - Hante, Falk T1 - Lipschitz Continuity of the Value Function in Mixed-Integer Optimal Control Problems JF - Mathematics of Control, Signals, and Systems Y1 - 2017 U6 - https://doi.org/10.1007/s00498-016-0183-4 VL - 29 IS - 1 ER - TY - INPR A1 - Gugat, Martin A1 - Schultz, Rüdiger T1 - Boundary feedback stabilization of the isothermal Euler-equations with uncertain boundary data N2 - In a gas transport system, the customer behavior is uncertain. Motivated by this situation, we consider a boundary stabilization problem for the flow through a gas pipeline, where the outflow at one end of the pipe %that is governed by the customer's behavior is uncertain. The control action is located at the other end of the pipe. The feedback law is a classical Neumann velocity feedback with a feedback parameter $k>0$. We show that as long as the $H^1$-norm of the function that describes the noise in the customer's behavior decays exponentially with a rate that is sufficiently large, the velocity of the gas can be stabilized exponentially fast in the sense that a suitably chosen Lyapunov function decays exponentially. For the exponential stability it is sufficient that the feedback parameter $k$ is sufficiently large and the stationary state to which the system is stabilized is sufficiently small. The stability result is local, that is it holds for initial states that are sufficiently close to the stationary state. This result is an example for the exponential boundary feedback stabilization of a quasilinear hyperbolic system with uncertain boundary data. The analysis is based upon the choice of a suitably Lyapunov function. The decay of this Lyapunov function implies that also the $L^2$-norm of the difference of the system state and the stationary state decays exponentially. Y1 - 2017 ER - TY - CHAP A1 - Gugat, Martin A1 - Herty, Michael A1 - Yu, Hui ED - Westdickenberg, Michael T1 - On the relaxation approximation for 2 × 2 hyperbolic balance laws N2 - The relaxation approximation for systems of conservation laws has been studied intensively for example by [17, 5, 19, 24]. In this paper the corresponding relaxation approximation for 2x2 systems of balance laws is studied. Our driving example is gas flow in pipelines described by the isothermal Euler equations. We are interested in the limiting behavior as the relaxation parameter tends to zero. We give conditions where the relaxation converges to the states of the original system and counterexamples for cases where the steady states depend on the space variable. Y1 - 2017 VL - x IS - x SP - x EP - x PB - Springer CY - Springer Proceedings in Mathematics and Statistics ER - TY - JOUR A1 - Gugat, Martin A1 - Ulbrich, Stefan ED - CHEN, GOONG T1 - The isothermal Euler equations for ideal gas with source term: Product solutions, flow reversal and no blow up JF - Journal of Mathematical Analysis and Applications KW - Global classical solutions Ideal gas Bi-directional flow Transsonic flow Y1 - 2017 U6 - https://doi.org/10.1016/j.jmaa.2017.04.064 VL - 454 IS - 1 SP - 439 EP - 452 ER - TY - JOUR A1 - Gugat, Martin A1 - Keimer, Alexander A1 - Leugering, Günter A1 - Wang, Zhiqiang ED - Piccoli, Benedetto T1 - Analysis of a system of nonlocal conservation laws for multi-commodity flow on networks JF -  Networks and Heterogeneous Media N2 - We consider a system of scalar nonlocal conservation laws on networks that model a highly re-entrant multi-commodity manufacturing system as encountered in semiconductor production. Every single commodity is mod-eled by a nonlocal conservation law, and the corresponding PDEs are coupled via a collective load, the work in progress. We illustrate the dynamics for two commodities. In the applications, directed acyclic networks naturally occur, therefore this type of networks is considered. On every edge of the network we have a system of coupled conservation laws with nonlocal velocity. At the junctions the right hand side boundary data of the foregoing edges is passed as left hand side boundary data to the following edges and PDEs. For distributing junctions, where we have more than one outgoing edge, we impose time dependent distribution functions that guarantee conservation of mass. We provide results of regularity, existence and well-posedness of the multi-commodity network model for L p-, BV-and W 1,p-data. Moreover, we define an L 2-tracking type objective and show the existence of minimizers that solve the corresponding optimal control problem. KW - conservation laws on network KW - nonlocal conservation laws KW - optimal nodal control KW - systems of hyperbolic pdes Y1 - 2016 U6 - https://doi.org/DOI: 10.3934/nhm.2015.10.749 VL - 10 IS - 4 SP - 749 EP - 785 ER - TY - JOUR A1 - Gugat, Martin A1 - Leugering, Günter ED - Zuazua, Enrique T1 - Time delay in optimal control loops for wave equations JF - ESAIM: COCV N2 - In optimal control loops delays can occur, for example through transmission via digital communication channels. Such delays influence the state that is generated by the implemented control. We study the effect of a delay in the implementation of L 2-norm minimal Neumann boundary controls for the wave equation. The optimal controls are computed as solutions of problems of exact optimal control, that is if they are implemented without delay, they steer the system to a position of rest in a given finite time T. We show that arbitrarily small delays δ > 0 can have a destabilizing effect in the sense that we can find initial states such that if the optimal control u is implemented in the form yx(t, 1) = u(t − δ) for t > δ, the energy of the system state at the terminal time T is almost twice as big as the initial energy. We also show that for more regular initial states, the effect of a delay in the implementation of the optimal control is bounded above in the sense that for initial positions with derivatives of BV-regularity and initial velocities with BV-regularity, the terminal energy is bounded above by the delay δ multiplied with a factor that depends on the BV-norm of the initial data. We show that for more general hyperbolic optimal exact control problems the situation is similar. For systems that have arbitrarily large eigenvalues, we can find terminal times T and arbitrarily small time delays δ, such that at the time T + δ, in the optimal control loop with delay the norm of the state is twice as large as the corresponding norm for the initial state. Moreover, if the initial state satisfies an additional regularity condition, there is an upper bound for the effect of time delay of the order of the delay with a constant that depends on the initial state only. KW - PDE constrained optimization KW - delay KW - wave equation KW - boundary control KW - hyperbolic system Y1 - 2016 U6 - https://doi.org/http://dx.doi.org/10.1051/cocv/2015038 ER - TY - JOUR A1 - Gugat, Martin A1 - Wintergerst, David A1 - Schultz, Rüdiger ED - Iske, Armin T1 - Networks of pipelines for gas with nonconstant compressibility factor: stationary states JF - Computational and Applied Mathematics N2 - For the management of gas transportation networks, it is essential to know how the stationary states of the system are determined by the boundary data. The isothermal Euler equations are an accurate pde-model for the gas flow through each pipe. A compressibility factor is used to model the nonlinear relationship between density and pressure that occurs in real gas in contrast to ideal gas. The gas flow through the nodes is governed by algebraic node conditions that require the conservation of mass and the continuity of the pressure. We examine networks that are described by arbitrary finite graphs and show that for suitably chosen boundary data, subsonic stationary states exist and are uniquely determined by the boundary data. Our construction of the stationary states is based upon explicit representations of the stationary states on each single pipe that can easily be evaluated numerically. We also use the monotonicity properties of these states as functions of the boundary data. Y1 - 2016 U6 - https://doi.org/10.1007/s40314-016-0383-z ER - TY - INPR A1 - Gugat, Martin A1 - Perrollaz, Vincent A1 - Rosier, Lionel T1 - Boundary stabilization of quasilinear hyperbolic systems of balance laws: Exponential decay for small source terms N2 - We investigate the long-time behavior of solutions of quasilinear hyperbolic systems with transparent boundary conditions when small source terms are incorporated in the system. Even if the finite-time stability of the system is not preserved, it is shown here that an exponential convergence towards the steady state still holds with a decay rate which is proportional to the logarithm of the amplitude of the source term. The result is stated for a system with dynamical boundary conditions in order to deal with initial data that are free of any compatibility condition. Y1 - 2017 ER - TY - JOUR A1 - Gugat, Martin ED - Khalique, Chaudry T1 - Exact Boundary Controllability for Free Traffic Flow with Lipschitz Continuous State JF - Mathematical Problems in Engineering N2 - We consider traffic flow governed by the LWR model. We show that a Lipschitz continuous initial density with free-flow and sufficiently small Lipschitz constant can be controlled exactly to an arbitrary constant free-flow density in finite time by a piecewise linear boundary control function that controls the density at the inflow boundary if the outflow boundary is absorbing. Moreover, this can be done in such a way that the generated state is Lipschitz continuous. Since the target states need not be close to the initial state, our result is a global exact controllability result. The Lipschitz constant of the generated state can be made arbitrarily small if the Lipschitz constant of the initial density is sufficiently small and the control time is sufficiently long. This is motivated by the idea that finite or even small Lipschitz constants are desirable in traffic flow since they might help to decrease the speed variation and lead to safer traffic. KW - Traffic flow Y1 - 2016 U6 - https://doi.org/doi:10.1155/2016/2743251 Creative Commons Attribution License VL - 2016 IS - 2016 SP - 1 EP - 11 ER - TY - JOUR A1 - Gugat, Martin ED - Bloch, Anthony M. T1 - Exponential Stabilization of the Wave Equation by Dirichlet Integral Feedback JF - SIAM Journal on Control and Optimization (SICON) N2 - We consider the problem of boundary feedback stabilization of a vibrating string that is fixed at one end and with control action at the other end. In contrast to previous studies that have required L 2-regularity for the initial position and H −1-regularity for the initial velocity, in this paper we allow for initial positions with L 1-regularity and initial velocities in W −1,1 on the space interval. It is well known that for a certain feedback parameter, for sufficiently regular initial states the classical energy of the closed-loop system with Neumann velocity feedback is controlled to zero after a finite time that is equal to the minimal time where exact controllability holds. In this paper, we present a Dirichlet boundary feedback that yields a well-defined closed-loop system in the (L 1 , W −1,1) framework and also has this property. Moreover, for all positive feedback parameters our feedback law leads to exponential decay of a suitably defined L 1-energy. For more regular initial states with (L 2 , H −1) regularity, the proposed feedback law leads to exponential decay of an energy that corresponds to this framework. If the initial states are even more regular with H 1-regularity of the initial position and L 2-regularity of the initial velocity, our feedback law also leads to exponential decay of the classical energy. KW - exponential stability KW - Dirichlet boundary control KW - energy decay KW - exact control Y1 - 2016 U6 - https://doi.org/DOI: 10.1137/140977023 VL - 53 IS - 1 SP - 526 EP - 546 ER - TY - JOUR A1 - Gugat, Martin A1 - Trelat, Emmanuel A1 - Zuazua, Enrique ED - Sepulchre, Rodolphe T1 - Optimal Neumann control for the 1D wave equation: Finite horizon, infinite horizon, boundary tracking terms and the turnpike property JF - Systems & Control Letters N2 - We consider a vibrating string that is fixed at one end with Neumann control action at the other end. We investigate the optimal control problem of steering this system from given initial data to rest, in time TT, by minimizing an objective functional that is the convex sum of the L2L2-norm of the control and of a boundary Neumann tracking term. We provide an explicit solution of this optimal control problem, showing that if the weight of the tracking term is positive, then the optimal control action is concentrated at the beginning and at the end of the time interval, and in-between it decays exponentially. We show that the optimal control can actually be written in that case as the sum of an exponentially decaying term and of an exponentially increasing term. This implies that, if the time TT is large, then the optimal trajectory approximately consists of three arcs, where the first and the third short-time arcs are transient arcs, and in the middle arc the optimal control and the corresponding state are exponentially close to 00. This is an example of a turnpike phenomenon for a problem of optimal boundary control. If T=+∞T=+∞ (infinite time horizon problem), then only the exponentially decaying component of the control remains, and the norms of the optimal control action and of the optimal state decay exponentially in time. In contrast to this situation, if the weight of the tracking term is zero and only the control cost is minimized, then the optimal control is distributed uniformly along the whole interval [0,T][0,T] and coincides with the control given by the Hilbert Uniqueness Method. In addition, we establish a similarity theorem stating that, for every T>0T>0, there exists an appropriate weight λ<1λ<1 for which the optimal solutions of the corresponding finite horizon optimal control problem and of the infinite horizon optimal control problem coincide along the first part of the time interval [0,2][0,2]. We also discuss the turnpike phenomenon from the perspective of a general framework with a strongly continuous semi-group. KW - Neumann boundary control KW - Turnpike phenomenon KW - Exact control Y1 - 2016 U6 - https://doi.org/10.1016/j.sysconle.2016.02.001 VL - 90 SP - 61 EP - 70 ER - TY - INPR A1 - Wintergerst, David A1 - Gugat, Martin T1 - Finite Time Blow-up of Traveling Wave Solutions for the Flow of Real Gas through Pipeline Networks N2 - In the context of gas transportation, analytical solutions are essential for the understanding of the underlying dynamics described by a system of partial differential equations. We derive traveling wave solutions for the 1-d isothermal Euler equations. A non-constant compressibility factor is used to describe the correlation between density and pressure. The blow-up of the traveling wave solution in � finite time is proven. We then extend our analysis to networks under appropriate coupling conditions and derive compatibility conditions to fulfill these coupling conditions. KW - isothermal Euler equations KW - real gas KW - finite time blow-up KW - traveling waves KW - networks Y1 - 2016 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 - JOUR A1 - Gugat, Martin A1 - Weiland, Sven T1 - Nodal Stabilization of the Flow in a Network with a Cycle JF - Journal of Optimization, Differential Equations and their Applications N2 - In this paper we discuss an approach to the stability analysis for classical solutions of closed loop systems that is based upon the tracing of the evolution of the Riemann invariants along the characteristics. We consider a network where several edges are coupled through node conditions that govern the evolution of the Riemann invariants through the nodes of the network. The analysis of the decay of the Riemann invariants requires to follow backwards all the characteristics that enter such a node and contribute to the evolution. This means that with each nodal reflection/crossing the number of characteristics that contribute to the evolution increases. We show how for simple networks with a suffcient number of damping nodal controlers it is possible to keep track of this family of characteristics and use this approach to analyze the exponential stability of the system. The analysis is based on an adapted version of Gronwall's lemma that allows us to take into account the possible increase of the Riemann invariants when the characteristic curves cross a node of the network. Our example is motivated by applications in the control of gas pipeline flow, where the graphs of the networks often contain many cycles. KW - Nodal Stabilization KW - Classical Solutions KW - Cycled Networks KW - Gas Pipelines KW - Nodal Control Y1 - 2021 VL - 29 IS - 2 SP - 1 EP - 24 ER - TY - JOUR A1 - Gugat, Martin A1 - Hante, Falk T1 - On the turnpike phenomenon for optimal boundary control problems with hyperbolic systems JF - SIAM Journal on Control and Optimization N2 - We study problems of optimal boundary control with systems governed by linear hyperbolic partial differential equations. The objective function is quadratic and given by an integral over the finite time interval (0,T) that depends on the boundary traces of the solution. If the time horizon T is sufficiently large, the solution of the dynamic optimal boundary control problem can be approximated by the solution of a steady state optimization problem. We show that for T to infinity the approximation error converges to zero in the sense of the norm in L^2(0,1) with the rate 1/T, if the time interval (0,T) is transformed to the fixed interval (0,1). Moreover, we show that also for optimal boundary control problems with integer constraints for the controls the turnpike phenomenon occurs. In this case the steady state optimization problem also has the integer constraints. If T is sufficiently large, the integer part of each solution of the dynamic optimal boundary control problem with integer constraints is equal to the integer part of a solution of the static problem. A numerical verification is given for a control problem in gas pipeline operations. Y1 - 2019 U6 - https://doi.org/10.1137/17M1134470 VL - 57 IS - 1 SP - 264 EP - 289 ER - TY - JOUR A1 - Farshbaf Shaker, Mohammad Hassan A1 - Gugat, Martin A1 - Heitsch, Holger A1 - Henrion, René T1 - Optimal Neumann boundary control of a vibrating string with uncertain initial data and probabilistic terminal constraints N2 - In optimal control problems, often initial data are required that are not known exactly in practice. In order to take into account this uncertainty, we consider optimal control problems for a system with an uncertain initial state. A finite terminal time is given. On account of the uncertainty of the initial state, it is not possible to prescribe an exact terminal state. Instead, we are looking for controls that steer the system into a given neighborhood of the desired terminal state with sufficiently high probability. This neighborhood is described in terms of an inequality for the terminal energy. The probabilistic constraint in the considered optimal control problem leads to optimal controls that are robust against the inevitable uncertainties of the initial state. We show the existence of such optimal controls. Numerical examples with optimal Neumann control of the wave equation are presented. KW - PDE constrained optimization, probabilistic constraints, uncertain initial data Y1 - U6 - https://doi.org/10.1137/19M1269944 ER - TY - JOUR A1 - Adelhütte, Dennis A1 - Aßmann, Denis A1 - Gonzàlez Grandòn, Tatiana A1 - Gugat, Martin A1 - Heitsch, Holger A1 - Liers, Frauke A1 - Henrion, René A1 - Nitsche, Sabrina A1 - Schultz, Rüdiger A1 - Stingl, Michael A1 - Wintergerst, David T1 - Joint model of probabilistic/robust (probust) constraints applied to gas network optimization N2 - Optimization tasks under uncertain conditions abound in many real-life applications. Whereas solution approaches for probabilistic constraints are often developed in case the uncertainties can be assumed to follow a certain probability distribution, robust approaches are usually used in case solutions are sought that are feasible for all realizations of uncertainties within some pre-defined uncertainty set. As many applications contain different types of uncertainties that require robust as well as probabilistic treatments, we deal with a class of joint probabilistic/robust constraints as its appears in optimization problems under uncertainty. Focusing on complex uncertain gas network optimization problems, we show the relevance of this class of problems for the task of maximizing free booked capacities in an algebraic model for a stationary gas network. We furthermore present approaches for their solution. Finally, we study the problem of controlling a transient system that is governed by the wave equation. The task consists in determining controls such that a certain robustness measure remains below some given upper bound, with high probability. KW - robust optimization KW - chance constraints KW - optimal control KW - spheric-radial decomposition Y1 - 2017 U6 - https://doi.org/10.1007/s10013-020-00434-y ER - TY - JOUR A1 - Gugat, Martin A1 - Schultz, Rüdiger A1 - Schuster, Michael T1 - Convexity and Starshapedness of Feasible Sets in Stationary Flow Networks JF - Networks and Heterogeneous Media N2 - We deal with a stationary model for flow through a network. The flows are determined by the values at the boundary nodes of the network that we call the loads of the network. In the applications , the feasible loads must satisfy some box constraints. We analyze the structure of the set of feasible loads. Our analysis is motivated by gas pipeline flows, where the box constraints are pressure bounds. We present sufficient conditions to show, that the feasible set is star-shaped with respect to special points. For stronger conditions, we prove the convexity of the set of feasible loads. All the results are given for active and passive networks, i.e. networks with and without inner control. This analysis is motivated by the aim to use the spheric-radial decomposition for stochastic boundary data in this model. This paper can be used for simplifying the algorithmic use of the spheric-radial decomposition. KW - Validation of nominations, Stationary states, Isothermal Euler equations, Gas networks, Convexity, Star-Shapedness, Spheric-radial decomposition Y1 - 2020 U6 - https://doi.org/10.3934/nhm.2020008 SP - 171 EP - 195 ER - TY - JOUR A1 - Gugat, Martin A1 - Schuster, Michael ED - Weber, Gerhard-Wilhelm T1 - Stationary Gas Networks with Compressor Control and Random Loads: Optimization with Probabilistic Constraints JF - Mathematical Problems in Engineering N2 - We introduce a stationary model for gas flow based on simplified isothermal Euler equations in a non-cycled pipeline network. Especially the problem of the feasibility of a random load vector is analyzed. Feasibility in this context means the existence of a flow vectormeeting these loads, which satisfies the physical conservation laws with box constraints for the pressure. An important aspect of the model is the support of compressor stations, which counteract the pressure loss caused by friction in the pipes.The network is assumed to have only one influx node; all other nodes are efflux nodes.With these assumptions the set of feasible loads can be characterized analytically. In addition we show the existence of optimal solutions for some optimization problems with probabilistic constraints. A numerical example based on real data completes this paper. KW - Validation of nominations, Stationary states, Isothermal Euler equations, Ideal Gas, Gas network, Gas transport, Compressor control, Uncertainty, Chance Constraints, Spheric-radial decomposition Y1 - 2019 U6 - https://doi.org/https://doi.org/10.1155/2018/7984079 ER - TY - JOUR A1 - Gugat, Martin A1 - Steffensen, Sonja ED - Trelat, Emmanuel T1 - DYNAMIC BOUNDARY CONTROL GAMES WITH NETWORKS OF STRINGS JF - ESAIM: Control, Optimisation and Calculus of Variations (ESAIM: COCV) N2 - Consider a star-shaped network of strings. Each string is governed by the wave equation. At each boundary node of the network there is a player that performs Dirichlet boundary control action and in this way influences the system state. At the central node, the states are coupled by algebraic conditions in such a way that the energy is conserved. We consider the corresponding antagonistic game where each player minimizes a certain quadratic objective function that is given by the sum of a control cost and a tracking term for the final state. We prove that under suitable assumptions a unique Nash equilibrium exists and give an explicit representation of the equilibrium strategies. Y1 - 2018 U6 - https://doi.org/10.1051/cocv/2017082 ER - TY - JOUR A1 - Gugat, Martin A1 - Rosier, Lionel A1 - Perrollaz, Vincent ED - Arendt, W. ED - Pierre, M. T1 - Boundary stabilization of quasilinear hyperbolic systems of balance laws: exponential decay for small source terms JF - Journal of Evolution Equations N2 - We investigate the long-time behaviour of solutions of quasilinear hyperbolic systems with transparent boundary conditions when small source terms are incorporated in the system. Even if the finite-time stability of the system is not preserved, it is shown here that an exponential convergence towards the steady state still holds with a decay rate which is proportional to the logarithm of the amplitude of the source term. The result is stated for a system with dynamical boundary conditions in order to deal with initial data that are free of any compatibility condition. The proof of the existence and uniqueness of a solution defined for all positive times is also provided in this paper. Y1 - 2018 U6 - https://doi.org/10.1007/s00028-018-0449-z ER - TY - INPR 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 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 - JOUR A1 - Gugat, Martin A1 - Henrion, René A1 - Heitsch, Holger T1 - A turnpike property for optimal control problems with dynamic probabilistic constraints JF - Journal of Convex Analysis N2 - In this paper we consider systems that are governed by linear time-discrete dynamics with an initial condition and a terminal condition for the expected values. We study optimal control problems where in the objective function a term of tracking type for the expected values and a control cost appear. In addition, the feasible states have to satisfy a conservative probabilistic constraint that requires that the probability that the trajectories remain in a given set F is greater than or equal to a given lower bound. An application are optimal control problems related to storage management systems with uncertain in- and output. We give suffcient conditions that imply that the optimal expected trajectories remain close to a certain state that can be characterized as the solution of an optimal control problem without prescribed initial- and terminal condition. Hence we contribute to the study of the turnpike phenomenon that is well-known in mathematical economics. KW - Probabilistic Constraints KW - Probabilistic Robustness KW - here-and-now decision KW - Turnpike phenomenon KW - Measure turnpike Y1 - 2021 VL - 30 IS - 3 SP - 1025 EP - 1052 PB - Heldermann Verlag 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 - 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 - 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 - Gugat, Martin A1 - Schuster, Michael A1 - Sokolowski, Jan T1 - Location Problem for Compressor Stations in Pipeline Networks N2 - In the operation of pipeline networks, compressors play a crucial role in ensuring the network’s functionality for various scenarios. In this contribution we address the important question of finding the optimal location of the compressors. This problem is of a novel structure, since it is related with the gas dynamics that governs the network flow. That results in non-convex mixed integer stochastic optimization problems with probabilistic constraints. Using a steady state model for the gas flow in pipeline networks including compressor control and uncertain loads given by certain probability distributions, the problem of finding the optimal location for the control on the network, s.t. the control cost is minimal and the gas pressure stays within given bounds, is considered. In the deterministic setting, explicit bounds for the pipe length and the inlet pressure, s.t. a unique optimal compressor location with minimal control cost exists, are presented. In the probabilistic setting, an existence result for the optimal compressor location is presented and the uniqueness of the solution is discussed depending on the probability distribution. For Gaussian distributed loads a uniqueness result for the optimal compressor location is presented. Further the problem of finding the optimal compressor locations on networks including the number of compressor stations as variable is considered. Results for the existence of optimal locations on a graph in both, the deterministic and the probabilistic setting, are presented and the uniqueness of the solutions is discussed depending on probability distributions and graph topology. The paper concludes with an illustrative example demonstrating that the compressor locations determined using a steady state approach are also admissible in transient settings. KW - Gas Networks KW - Compressor Control KW - Weber Problem KW - Optimal Location KW - Uncertain Boundary Data Y1 - 2024 ER - TY - JOUR A1 - Gugat, Martin A1 - Qian, Meizhi A1 - Sokolowski, Jan T1 - Network Design and Control: Shape and Topology Optimization for the Turnpike Property for the Wave Equation N2 - The optimal control problems for the wave equation are considered on networks. The turnpike property is shown for the state equation, the adjoint state equation as well as the optimal cost. The shape and topology optimization is performed for the network with the shape functional given by the optimality system of the control problem. The set of admissible shapes for the network is compact in finite dimensions, thus the use of turnpike property is straightforward. The topology optimization is analysed for an example of nucleation of a small cycle at the internal node of network. The topological derivative of the cost is introduced and evaluated in the framework of domain decomposition technique. Numerical examples are provided. KW - Turnpike Property KW - Wave Equation KW - Optimal Control KW - Spectral Method KW - Network Optimum Design Y1 - 2024 VL - J. Geom. Anal. IS - 34 ER -