Implicit-Explicit Runge-Kutta schemes for numerical discretization of optimal control problems

Implicit-explicit (IMEX) Runge-Kutta methods play a major rule in the numerical treatment of differential systems governed by stiff and non-stiff terms. This paper discusses order conditions and symplecticity properties of a class of IMEX Runge-Kutta methods in the context of optimal control problems. The analysis of the schemes is based on the continuous optimality system. Using suitable transformations of the adjoint equation, order conditions up to order three are proven as well as the relation between adjoint schemes obtained through different transformations is investigated. Conditions for the IMEX Runge-Kutta methods to be symplectic are also derived. A numerical example illustrating the theoretical properties is presented

Comparative model accuracy of a data-fitted generalized Aw-Rascle-Zhang model

The Aw-Rascle-Zhang (ARZ) model can be interpreted as a generalization of the Lighthill-Whitham-Richards (LWR) model, possessing a family of fundamental diagram curves, each of which represents a class of drivers with a different empty road velocity. A weakness of this approach is that different drivers possess vastly different densities at which traffic flow stagnates. This drawback can be overcome by modifying the pressure relation in the ARZ model, leading to the generalized Aw-Rascle-Zhang (GARZ) model. We present an approach to determine the parameter functions of the GARZ model from fundamental diagram measurement data. The predictive accuracy of the resulting data-fitted GARZ model is compared to other traffic models by means of a three-detector test setup, employing two types of data: vehicle trajectory data, and sensor data. This work also considers the extension of the ARZ and the GARZ models to models with a relaxation term, and conducts an investigation of the optimal relaxation time.Comment: 30 pages, 10 figures, 3 table

Modeling and analysis of pooled stepped chutes

We consider an application of pooled stepped chutes where the transport in each pooled step is described by the shallow--water equations. Such systems can be found for example at large dams in order to release overflowing water. We analyze the mathematical conditions coupling the flows between different chutes taken from the engineering literature. We present the solution to a Riemann problem in the large and also a well--posedness result for the coupled problem. We finally report on some numerical experiments.Comment: 17 pages, 31 figure

Kinetic description of optimal control problems and applications to opinion consensus

In this paper an optimal control problem for a large system of interacting agents is considered using a kinetic perspective. As a prototype model we analyze a microscopic model of opinion formation under constraints. For this problem a Boltzmann-type equation based on a model predictive control formulation is introduced and discussed. In particular, the receding horizon strategy permits to embed the minimization of suitable cost functional into binary particle interactions. The corresponding Fokker-Planck asymptotic limit is also derived and explicit expressions of stationary solutions are given. Several numerical results showing the robustness of the present approach are finally reported.Comment: 25 pages, 18 figure

Linear multistep methods for optimal control problems and applications to hyperbolic relaxation systems

We are interested in high-order linear multistep schemes for time discretization of adjoint equations arising within optimal control problems. First we consider optimal control problems for ordinary differential equations and show loss of accuracy for Adams-Moulton and Adams-Bashford methods, whereas BDF methods preserve high--order accuracy. Subsequently we extend these results to semi--lagrangian discretizations of hyperbolic relaxation systems. Computational results illustrate theoretical findings

A two-dimensional data-driven model for traffic flow on highways

Based on experimental traffic data obtained from German and US highways, we propose a novel two-dimensional first-order macroscopic traffic flow model. The goal is to reproduce a detailed description of traffic dynamics for the real road geometry. In our approach both the dynamic along the road and across the lanes is continuous. The closure relations, being necessary to complete the hydrodynamic equation, are obtained by regression on fundamental diagram data. Comparison with prediction of one-dimensional models shows the improvement in performance of the novel model.Comment: 27 page

Optimal treatment planning governed by kinetic equations

In this paper we study a problem in radiotherapy treatment planning, where the evolution of the radiation field is governed by a deterministic Boltzmann transport equation. We show existence, uniqueness and regularity of solutions to an optimal dose distribution problem constrained by the Boltzmann Continuous Slowing-Down equation in an appropriate function space. The main new difficulty is the treatment of the stopping power term. Furthermore, we characterize optimal controls for problems governed by this transport equation.Comment: 15 pages, 1 figur

Hybrid stochastic kinetic description of two-dimensional traffic dynamics

In this work we present a two-dimensional kinetic traffic model which takes into account speed changes both when vehicles interact along the road lanes and when they change lane. Assuming that lane changes are less frequent than interactions along the same lane and considering that their mathematical description can be done up to some uncertainty in the model parameters, we derive a hybrid stochastic Fokker-Planck-Boltzmann equation in the quasi-invariant interaction limit. By means of suitable numerical methods, precisely structure preserving and direct Monte Carlo schemes, we use this equation to compute theoretical speed-density diagrams of traffic both along and across the lanes, including estimates of the data dispersion, and validate them against real data