Source author record

Benedetto Piccoli

Benedetto Piccoli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2026arXiv

Supervised and Unsupervised Neural Network Solver for First Order Hyperbolic Nonlinear PDEs

We present a neural network-based method for learning scalar hyperbolic conservation laws. Our method replaces the traditional numerical flux in finite volume schemes with a trainable neural network while preserving the conservative structure of the scheme. The model can be trained both in a supervised setting with efficiently generated synthetic data or in an unsupervised manner, leveraging the weak formulation of the partial differential equation. We provide theoretical results that our model can perform arbitrarily well, and provide associated upper bounds on neural network size. Extensive experiments demonstrate that our method often outperforms efficient schemes such as Godunov's scheme, WENO, and Discontinuous Galerkin for comparable computational budgets. Finally, we demonstrate the effectiveness of our method on a traffic prediction task, leveraging field experimental highway data from the Berkeley DeepDrive drone dataset.

preprint2022arXiv

A measure model for the spread of viral infections with mutations

Genetic variations in the COVID-19 virus are one of the main causes of the COVID-19 pandemic outbreak in 2020 and 2021. In this article, we aim to introduce a new type of model, a system coupled with ordinary differential equations (ODEs), and measure differential equation (MDE), stemming from the classical SIR model for the variants distribution. Specifically, we model the evolution of susceptible $S$ and removed $R$ populations by ODEs and the infected $I$ population by an MDE comprised of a probability vector field (PVF) and a source term. In addition, the ODEs for $S$ and $R$ contain terms that are related to the measure $I$. We establish analytically the well-posedness of the coupled ODE-MDE system by using generalized Wasserstein distance. We give two examples to show that the proposed ODE-MDE model coincides with the classical SIR model in the case of constant or time-dependent

preprint2022arXiv

A rigorous multi-population multi-lane hybrid traffic model and its mean-field limit for dissipation of waves via autonomous vehicles

In this paper, a multi-lane multi-population microscopic model, which presents stop and go waves, is proposed to simulate traffic on a ring-road. Vehicles are divided between human-driven and autonomous vehicles (AV). Control strategies are designed with the ultimate goal of using a small number of AVs (less than 5\% penetration rate) to represent Lagrangian control actuators that can smooth the multilane traffic flow and dissipate the stop-and-go waves. This in turn may reduce fuel consumption and emissions. The lane-changing mechanism is based on three components that we treat as parameters in the model: safety, incentive and cool-down time. The choice of these parameters in the lane-change mechanism is critical to modeling traffic accurately, because different parameter values can lead to drastically different traffic behaviors. In particular, the number of lane-changes and the speed variance are highly affected by the choice of parameters. Despite this modeling issue, when using sufficiently simple and robust controllers for AVs, the stabilization of uniform flow steady-state is effective for any realistic value of the parameters, and ultimately bypasses the observed modeling issue. Our approach is based on accurate and rigorous mathematical models, which allows a limit procedure that is termed, in gas dynamic terminology, mean-field. In simple words, from increasing the human-driven population to infinity, a system of coupled ordinary and partial differential equations are obtained. Moreover, control problems also pass to the limit, allowing the design to be tackled at different scales.

preprint2022arXiv

Emissions minimization on road networks via Generic Second Order Models

In this paper we consider the problem of estimating emissions due to vehicular traffic on complex networks, and minimizing their effect by regulating traffic at junctions. For the traffic evolution, we consider a Generic Second Order Model, which encompasses the majority of two-equations (i.e. second-order) models available in the literature, and extend it to road networks with merge and diverge junctions. The dynamics on the whole network is determined by selecting a solution to the Riemann Problems at junctions, i.e. the Cauchy problems with constant initial data on each incident road. The latter are solved assuming the maximization of the flow and assigning a traffic distribution coefficient for outgoing roads of diverges, and a priority rule for incoming roads of merges. A general emission model is considered and its parameters are tuned to the NOx emission rate. The minimization of emissions is then formulated in terms of the traffic distribution and priority parameters, taking into account travel times. A comparison is provided between roundabouts with optimized parameters and traffic lights, which correspond to time-varying traffic priorities. Our approach can be adapted to manage traffic in complex networks in order to reduce emissions while keeping travel time at acceptable levels.

preprint2022arXiv

Limitations and Improvements of the Intelligent Driver Model (IDM)

This contribution analyzes the widely used and well-known "intelligent driver model (briefly IDM), which is a second order car-following model governed by a system of ordinary differential equations. Although this model was intensively studied in recent years for properly capturing traffic phenomena and driver braking behavior, a rigorous study of the well-posedness has, to our knowledge, never been performed. First it is shown that, for a specific class of initial data, the vehicles' velocities become negative or even diverge to $-\infty$ in finite time, both undesirable properties for a car-following model. Various modifications of the IDM are then proposed in order to avoid such ill-posedness. The theoretical remediation of the model, rather than post facto by ad-hoc modification of code implementations, allows a more sound numerical implementation and preservation of the model features. Indeed, to avoid inconsistencies and ensure dynamics close to the one of the original model, one may need to inspect and clean large input data, which may result in practically impossible scenarios for large-scale simulations. Although well-posedness issues occur only for specific initial data, this may happen frequently when different traffic scenarios are analyzed, and especially in presence of lane-changing, on ramps and other network components as it is the case for most commonly used micro-simulators. On the other side, it is shown that well-posedness can be guaranteed by straight-forward improvements, such as those obtained by slightly changing the acceleration to prevent the velocity from becoming negative.

preprint2022arXiv

Optimization of vaccination for COVID-19 in the midst of a pandemic

During the Covid-19 pandemic a key role is played by vaccination to combat the virus. There are many possible policies for prioritizing vaccines, and different criteria for optimization: minimize death, time to herd immunity, functioning of the health system. Using an age-structured population compartmental finite-dimensional optimal control model, our results suggest that the eldest to youngest vaccination policy is optimal to minimize deaths. Our model includes the possible infection of vaccinated populations. We apply our model to real-life data from the US Census for New Jersey and Florida, which have a significantly different population structure. We also provide various estimates of the number of lives saved by optimizing the vaccine schedule and compared to no vaccination.

preprint2022arXiv

Stability of multi-population traffic flows

Traffic waves, known also as stop-and-go waves or phantom hams, appear naturally as traffic instabilities, also in confined environments as a ring-road. A multi-population traffic is studied on a ring-road, comprised of drivers with stable and unstable behavior. There exists a critical penetration rate of stable vehicles above which the system is stable, and under which the system is unstable. In the latter case, stop-and-go waves appear, provided enough cars are on the road. The critical penetration rate is explicitly computable, and, in reasonable situations, a small minority of aggressive drivers is enough to destabilize an otherwise very stable flow. This is a source of instability that a single population model would not be able to explain. Also, the multi-population system can be stable below the critical penetration rate if the number of cars is sufficiently small. Instability emerges as the number of cars increases, even if the traffic density remains the same (i.e. number of cars and road size increase similarly). This shows that small experiments could lead to deducing imprecise stability conditions.

preprint2021arXiv

Generalized solutions to bounded-confidence models

Bounded-confidence models in social dynamics describe multi-agent systems, where each individual interacts only locally with others. Several models are written as systems of ordinary differential equations with discontinuous right-hand side: this is a direct consequence of restricting interactions to a bounded region with non-vanishing strength at the boundary. Various works in the literature analyzed properties of solutions, such as barycenter invariance and clustering. On the other side, the problem of giving a precise definition of solution, from an analytical point of view, was often overlooked. However, a rich literature proposing different concepts of solution to discontinuous differential equations is available. Using several concepts of solution, we show how existence is granted under general assumptions, while uniqueness may fail even in dimension one, but holds for almost every initial conditions. Consequently, various properties of solutions depend on the used definition and initial conditions.

preprint2020arXiv

A multiscale model for traffic regulation via autonomous vehicles

Autonomous vehicles (AVs) allow new ways of regulating the traffic flow on road networks. Most of available results in this direction are based on microscopic approaches, where ODEs describe the evolution of regular cars and AVs. In this paper, we propose a multiscale approach, based on recently developed models for moving bottlenecks. Our main result is the proof of existence of solutions for open-loop controls with bounded variation.

preprint2020arXiv

Metabolic graphs, LIFE method and the modeling of drug action on Mycobacterium tuberculosis

This paper serves as a framework for designing advanced models for drug action on metabolism. Drug treatment may affect metabolism by either enhancing or inhibiting metabolic reactions comprising a metabolic network. We introduce the concept of \textit{metabolic graphs}, a generalization of hypergraphs having specialized features common to metabolic networks. Linear-in-flux-expression (briefly LIFE) is a methodology for analyzing metabolic networks and simulating virtual patients. We extend LIFE dynamics to be compatible with metabolic graphs, including the more complex interactions of enhancer and inhibitor molecules that affect biochemical reactions. We discuss results considering network structure required for existence and uniqueness of equilibria on metabolic graphs and show simulations of drug action on \textit{Mycobacterium tuberculosis} (briefly MTB)

preprint2019arXiv

Model-based assessment of the impact of driver-assist vehicles using kinetic theory

In this paper we consider a kinetic description of follow-the-leader traffic models, which we use to study the effect of vehicle-wise driver-assist control strategies at various scales, from that of the local traffic up to that of the macroscopic stream of vehicles. We provide a theoretical evidence of the fact that some typical control strategies, such as the alignment of the speeds and the optimisation of the time headways, impact on the local traffic features (for instance, the speed and headway dispersion responsible for local traffic instabilities) but have virtually no effect on the observable macroscopic traffic trends (for instance, the flux/throughput of vehicles). This unobvious conclusion, which is in very nice agreement with recent field studies on autonomous vehicles, suggests that the kinetic approach may be a valid tool for an organic multiscale investigation and possibly design of driver-assist algorithms.

preprint2016arXiv

Continuity of the path delay operator for dynamic network loading with spillback

This paper establishes the continuity of the path delay operators for dynamic network loading (DNL) problems based on the Lighthill-Whitham-Richards model, which explicitly capture vehicle spillback. The DNL describes and predicts the spatial-temporal evolution of traffic flow and congestion on a network that is consistent with established route and departure time choices of travelers. The LWR-based DNL model is first formulated as a system of partial differential algebraic equations (PDAEs). We then investigate the continuous dependence of merge and diverge junction models with respect to their initial/boundary conditions, which leads to the continuity of the path delay operator through the wave-front tracking methodology and the generalized tangent vector technique. As part of our analysis leading up to the main continuity result, we also provide an estimation of the minimum network supply without resort to any numerical computation. In particular, it is shown that gridlock can never occur in a finite time horizon in the DNL model.

preprint2016arXiv

Continuous-time link-based kinematic wave model: formulation, solution existence, and well-posedness

We present a continuous-time link-based kinematic wave model (LKWM) for dynamic traffic networks based on the scalar conservation law model. Derivation of the LKWM involves the variational principle for the Hamilton-Jacobi equation and junction models defined via the notions of demand and supply. We show that the proposed LKWM can be formulated as a system of differential algebraic equations (DAEs), which captures shock formation and propagation, as well as queue spillback. The DAE system, as we show in this paper, is the continuous-time counterpart of the link transmission model. In addition, we present a solution existence theory for the continuous-time network model and investigate continuous dependence of the solution on the initial data, a property known as well-posedness. We test the DAE system extensively on several small and large networks and demonstrate its numerical efficiency.

preprint2016arXiv

Control of reaction-diffusion equations on time-evolving manifolds

Among the main actors of organism development there are morphogens, which are signaling molecules diffusing in the developing organism and acting on cells to produce local responses. Growth is thus determined by the distribution of such signal. Meanwhile, the diffusion of the signal is itself affected by the changes in shape and size of the organism. In other words, there is a complete coupling between the diffusion of the signal and the change of the shapes. In this paper, we introduce a mathematical model to investigate such coupling. The shape is given by a manifold, that varies in time as the result of a deformation given by a transport equation. The signal is represented by a density, diffusing on the manifold via a diffusion equation. We show the non-commutativity of the transport and diffusion evolution by introducing a new concept of Lie bracket between the diffusion and the transport operator. We also provide numerical simulations showing this phenomenon.

preprint2016arXiv

Interaction Network, State Space and Control in Social Dynamics

In the present chapter we study the emergence of global patterns in large groups in first and second-order multi-agent systems, focusing on two ingredients that influence the dynamics: the interaction network and the state space. The state space determines the types of equilibrium that can be reached by the system. Meanwhile, convergence to specific equilibria depends on the connectivity of the interaction network and on the interaction potential. When the system does not satisfy the necessary conditions for convergence to the desired equilibrium, control can be exerted, both on finite-dimensional systems and on their mean-field limit.

preprint2016arXiv

On the Continuum Approximation of the On-and-off Signal Control on Dynamic Traffic Networks

In the modeling of traffic networks, a signalized junction is typically treated using a binary variable to model the on-and-off nature of signal operation. While accurate, the use of binary variables can cause problems when studying large networks with many intersections. Instead, the signal control can be approximated through a continuum approach where the on-and-off control variable is replaced by a priority parameter. Advantages of such approximation include elimination of the need for binary variables, lower time resolution requirements, and more flexibility and robustness in a decision environment. It also resolves the issue of discontinuous travel time functions arising from the context of dynamic traffic assignment. Despite these advantages in application, it is not clear from a theoretical point of view how accurate is such continuum approach; i.e., to what extent is this a valid approximation for the on-and-off case. The goal of this paper is to answer these basic research questions and provide further guidance for the application of such continuum signal model. In particular, by employing the Lighthill-Whitham-Richards model (Lighthill and Whitham, 1955; Richards, 1956) on a traffic network, we investigate the convergence of the on-and-off signal model to the continuum model in regimes of diminishing signal cycles. We also provide numerical analyses on the continuum approximation error when the signal cycles are not infinitesimal. As we explain, such convergence results and error estimates depend on the type of fundamental diagram assumed and whether or not vehicle spillback occurs in a network. Finally, a traffic signal optimization problem is presented and solved which illustrates the unique advantages of applying the continuum signal model instead of the on-and-off one.

preprint2016arXiv

Priority-based Riemann solver for traffic flow on networks

In this article we introduce a new Riemann solver for traffic flow on networks. The Priority Riemann solver (PRS) provides a solution at junctions by taking into consideration priorities for the incoming roads and maximization of through flux. We prove existence of solutions for the solver for junctions with up to two incoming and two outgoing roads and show numerically the comparison with previous Riemann solvers. Additionally, we introduce a second version of the solver that considers the priorities as softer constraints and illustrate numerically the differences between the two solvers.

preprint2016arXiv

Regularization of chattering phenomena via bounded variation control

In control theory, the term chattering is used to refer to strong oscillations of controls, such as an infinite number of switchings over a compact interval of times. In this paper we focus on three typical occurences of chattering: the Fuller phenomenon, referring to situations where an optimal control switches an infinite number of times over a compact set; the Robbins phenomenon, concerning optimal control problems with state constraints, meaning that the optimal trajectory touches the boundary of the constraint set an infinite number of times over a compact time interval; the Zeno phenomenon, referring as well to an infinite number of switchings over a compact set, for hybrid optimal control problems. From the practical point of view, when trying to compute an optimal trajectory, for instance by means of a shooting method, chattering may be a serious obstacle to convergence. In this paper we propose a general regularization procedure, by adding an appropriate penalization of the total variation. This produces a quasi-optimal control, and we prove that the family of quasi-optimal solutions converges to the optimal solution of the initial problem as the penalization tends to zero. Under additional assumptions, we also quantify the quasi-optimality property by determining a speed of convergence of the costs.

preprint2016arXiv

Second-order models and traffic data from mobile sensors

Mobile sensing enabled by GPS or smart phones has become an increasingly important source of traffic data. For sufficient coverage of the traffic stream, it is important to maintain a reasonable penetration rate of probe vehicles. From the standpoint of capturing higher-order traffic quantities such as acceleration/deceleration, emission and fuel consumption rates, it is desirable to examine the impact on the estimation accuracy of sampling frequency on vehicle position. Of the two issues raised above, the latter is rarely studied in the literature. This paper addresses the impact of both sampling frequency and penetration rate on mobile sensing of highway traffic. To capture inhomogeneous driving conditions and deviation of traffic from the equilibrium state, we employ the second-order phase transition model (PTM). Several data fusion schemes that incorporate vehicle trajectory data into the PTM are proposed. And, a case study of the NGSIM dataset is presented which shows the estimation results of various Eulerian and Lagrangian traffic quantities. The findings show that while first-order traffic quantities can be accurately estimated even with a low sampling frequency, higher-order traffic quantities, such as acceleration, deviation, and emission rate, tend to be misinterpreted due to insufficiently sampled vehicle locations. We also show that a correction factor approach has the potential to reduce the sensing error arising from low sampling frequency and penetration rate, making the estimation of higher-order quantities more robust against insufficient data coverage of the highway traffic.

preprint2016arXiv

Traffic regulation via controlled speed limit

We study an optimal control problem for traffic regulation via variable speed limit. The traffic flow dynamics is described with the Lighthill-Whitham-Richards (LWR) model with Newell-Daganzo flux function. We aim at minimizing the $L^2$ quadratic error to a desired outflow, given an inflow on a single road. We first provide existence of a minimizer and compute analytically the cost functional variations due to needle-like variation in the control policy. Then, we compare three strategies: instantaneous policy; random exploration of control space; steepest descent using numerical expression of gradient. We show that the gradient technique is able to achieve a cost within 10% of random exploration minimum with better computational performances.

preprint2015arXiv

Developmental Partial Differential Equations

In this paper, we introduce the concept of Developmental Partial Differential Equation (DPDE), which consists of a Partial Differential Equation (PDE) on a time-varying manifold with complete coupling between the PDE and the manifold's evolution. In other words, the manifold's evolution depends on the solution to the PDE, and vice versa the differential operator of the PDE depends on the manifold's geometry. DPDE is used to study a diffusion equation with source on a growing surface whose growth depends on the intensity of the diffused quantity. The surface may, for instance, represent the membrane of an egg chamber and the diffused quantity a protein activating a signaling pathway leading to growth. Our main objective is to show controllability of the surface shape using a fixed source with variable intensity for the diffusion. More specifically, we look for a control driving a symmetric manifold shape to any other symmetric shape in a given time interval. For the diffusion we take directly the Laplace-Beltrami operator of the surface, while the surface growth is assumed to be equal to the value of the diffused quantity. We introduce a theoretical framework, provide approximate controllability and show numerical results. Future applications include a specific model for the oogenesis of Drosophila melanogaster.

preprint2015arXiv

Optimal control of a bioreactor for biofuel production

Dynamic flux balance analysis of a bioreactor is based on the coupling between a dynamic problem, which models the evolution of biomass, feeding substrates and metabolites, and a linear program, which encodes the metabolic activity inside cells. We cast the problem in the language of optimal control and propose a hybrid formulation to model the full coupling between macroscopic and microscopic level. On a given location of the hybrid system we analyze necessary conditions given by the Pontryagin Maximum Principle and discuss the presence of singular arcs. In particular, for the single-input case we prove that optimal controls are bang-bang. For the multi-input case, under suitable assumptions, we prove that generically with respect to initial conditions optimal controls are bang-bang.

preprint2015arXiv

Optimal Control of a Collective Migration Model

Collective migration of animals in a cohesive group is rendered possible by a strategic distribution of tasks among members: some track the travel route, which is time and energy-consuming, while the others follow the group by interacting among themselves. In this paper, we study a social dynamics system modeling collective migration. We consider a group of agents able to align their velocities to a global target velocity, or to follow the group via interaction with the other agents. The balance between these two attractive forces is our control for each agent, as we aim to drive the group to consensus at the target velocity. We show that the optimal control strategies in the case of final and integral costs consist of controlling the agents whose velocities are the furthest from the target one: these agents sense only the target velocity and become leaders, while the uncontrolled ones sense only the group, and become followers. Moreover, in the case of final cost, we prove an "Inactivation" principle: there exist initial conditions such that the optimal control strategy consists of letting the system evolve freely for an initial period of time, before acting with full control on the agent furthest from the target velocity.

preprint2014arXiv

Control to flocking of the kinetic Cucker-Smale model

The well-known Cucker-Smale model is a macroscopic system reflecting flocking, i.e. the alignment of velocities in a group of autonomous agents having mutual interactions. In the present paper, we consider the mean-field limit of that model, called the kinetic Cucker-Smale model, which is a transport partial differential equation involving nonlocal terms. It is known that flocking is reached asymptotically whenever the initial conditions of the group of agents are in a favorable configuration. For other initial configurations, it is natural to investigate whether flocking can be enforced by means of an appropriate external force, applied to an adequate time-varying subdomain. In this paper we prove that we can drive to flocking any group of agents governed by the kinetic Cucker-Smale model, by means of a sparse centralized control strategy, and this, for any initial configuration of the crowd. Here, "sparse control" means that the action at each time is limited over an arbitrary proportion of the crowd, or, as a variant, of the space of configurations; "centralized" means that the strategy is computed by an external agent knowing the configuration of all agents. We stress that we do not only design a control function (in a sampled feedback form), but also a time-varying control domain on which the action is applied. The sparsity constraint reflects the fact that one cannot act on the whole crowd at every instant of time. Our approach is based on geometric considerations on the velocity field of the kinetic Cucker-Smale PDE, and in particular on the analysis of the particle flow generated by this vector field. The control domain and the control functions are designed to satisfy appropriate constraints, and such that, for any initial configuration, the velocity part of the support of the measure solution asymptotically shrinks to a singleton, which means flocking.

preprint2014arXiv

Mean-Field Sparse Optimal Control

We introduce the rigorous limit process connecting finite dimensional sparse optimal control problems with ODE constraints, modeling parsimonious interventions on the dynamics of a moving population divided into leaders and followers, to an infinite dimensional optimal control problem with a constraint given by a system of ODE for the leaders coupled with a PDE of Vlasov-type, governing the dynamics of the probability distribution of the followers. In the classical mean-field theory one studies the behavior of a large number of small individuals freely interacting with each other, by simplifying the effect of all the other individuals on any given individual by a single averaged effect. In this paper we address instead the situation where the leaders are actually influenced also by an external policy maker, and we propagate its effect for the number $N$ of followers going to infinity. The technical derivation of the sparse mean-field optimal control is realized by the simultaneous development of the mean-field limit of the equations governing the followers dynamics together with the $Γ$-limit of the finite dimensional sparse optimal control problems.

preprint2014arXiv

On properties of the Generalized Wasserstein distance

The Wasserstein distances $W_p$ ($p\geq 1$), defined in terms of solution to the Monge-Kantorovich problem, are known to be a useful tool to investigate transport equations. In particular, the Benamou-Brenier formula characterizes the square of the Wasserstein distance $W_2$ as the infimum of the kinetic energy, or action functional, of all vector fields moving one measure to the other. Another important property of the Wasserstein distances is the Kantorovich-Rubinstein duality stating the equality between the distance $W_1$ and the supremum of the integrals of Lipschitz continuous functions with Lipschitz constant bounded by one. An intrinsic limitation of Wasserstein distances is the fact that they are defined only between measures having the same mass. To overcome such limitation, we recently introduced the generalized Wasserstein distances $W_p^{a,b}$, defined in terms of both the classical Wasserstein distance $W_p$ and the total variation (or $L^1$) distance. Here $p$ plays the same role as for the classic Wasserstein distance, while $a$ and $b$ are weights for the transport and the total variation term. In this paper we prove two important properties of the generalized Wasserstein distances: 1) a generalized Benamou-Brenier formula providing the equality between $W_2^{a,b}$ and the supremum of an action functional, which includes a transport term (kinetic energy) and a source term. 2) a duality à la Kantorovich-Rubinstein establishing the equality between $W_1^{1,1}$ and the flat metric.

preprint2014arXiv

Runge-Kutta Discontinuous Galerkin Method for Traffic Flow Model on Networks

We propose a bound-preserving Runge-Kutta (RK) discontinuous Galerkin (DG) method as an efficient, effective and compact numerical approach for numerical simulation of traffic flow problems on networks, with arbitrary high order accuracy. Road networks are modeled by graphs, composed of a finite number of roads that meet at junctions. On each road, a scalar conservation law describes the dynamics, while coupling conditions are specified at junctions to define flow separation or convergence at the points where roads meet. We incorporate such coupling conditions in the RK DG framework, and apply an arbitrary high order bound preserving limiter to the RK DG method to preserve the physical bounds on the network solutions (car density). We showcase the proposed algorithm on several benchmark test cases from the literature, as well as several new challenging examples with rich solution structures. Modeling and simulation of Cauchy problems for traffic flows on networks is notorious for lack of uniqueness or (Lipschitz) continuous dependence. The discontinuous Galerkin method proposed here deals elegantly with these problems, and is perhaps the only realistic and efficient high-order method for network problems.

preprint2014arXiv

Sparse Stabilization and Control of Alignment Models

From a mathematical point of view self-organization can be described as patterns to which certain dynamical systems modeling social dynamics tend spontaneously to be attracted. In this paper we explore situations beyond self-organization, in particular how to externally control such dynamical systems in order to eventually enforce pattern formation also in those situations where this wished phenomenon does not result from spontaneous convergence. Our focus is on dynamical systems of Cucker-Smale type, modeling consensus emergence, and we question the existence of stabilization and optimal control strategies which require the minimal amount of external intervention for nevertheless inducing consensus in a group of interacting agents. We provide a variational criterion to explicitly design feedback controls that are componentwise sparse, i.e. with at most one nonzero component at every instant of time. Controls sharing this sparsity feature are very realistic and convenient for practical issues. Moreover, the maximally sparse ones are instantaneously optimal in terms of the decay rate of a suitably designed Lyapunov functional, measuring the distance from consensus. As a consequence we provide a mathematical justification to the general principle according to which "sparse is better" in the sense that a policy maker, who is not allowed to predict future developments, should always consider more favorable to intervene with stronger action on the fewest possible instantaneous optimal leaders rather than trying to control more agents with minor strength in order to achieve group consensus. We then establish local and global sparse controllability properties to consensus and, finally, we analyze the sparsity of solutions of the finite time optimal control problem where the minimization criterion is a combination of the distance from consensus and of the l1-norm of the control.

preprint2012arXiv

Generalized Wasserstein distance and its application to transport equations with source

In this article, we generalize the Wasserstein distance to measures with different masses. We study the properties of such distance. In particular, we show that it metrizes weak convergence for tight sequences. We use this generalized Wasserstein distance to study a transport equation with source, in which both the vector field and the source depend on the measure itself. We prove existence and uniqueness of the solution to the Cauchy problem when the vector field and the source are Lipschitzian with respect to the generalized Wasserstein distance.

preprint2012arXiv

How can macroscopic models reveal self-organization in traffic flow?

In this paper we propose a new modeling technique for vehicular traffic flow, designed for capturing at a macroscopic level some effects, due to the microscopic granularity of the flow of cars, which would be lost with a purely continuous approach. The starting point is a multiscale method for pedestrian modeling, recently introduced in Cristiani et al., Multiscale Model. Simul., 2011, in which measure-theoretic tools are used to manage the microscopic and the macroscopic scales under a unique framework. In the resulting coupled model the two scales coexist and share information, in the sense that the same system is simultaneously described from both a discrete (microscopic) and a continuous (macroscopic) perspective. This way it is possible to perform numerical simulations in which the single trajectories and the average density of the moving agents affect each other. Such a method is here revisited in order to deal with multi-population traffic flow on networks. For illustrative purposes, we focus on the simple case of the intersection of two roads. By exploiting one of the main features of the multiscale method, namely its dimension-independence, we treat one-dimensional roads and two-dimensional junctions in a natural way, without referring to classical network theory. Furthermore, thanks to the coupling between the microscopic and the macroscopic scales, we model the continuous flow of cars without losing the right amount of granularity, which characterizes the real physical system and triggers self-organization effects, such as, for example, the oscillatory patterns visible at jammed uncontrolled crossroads.

preprint2012arXiv

Numerical schemes for the optimal input flow of a supply-chain

An innovative numerical technique is presented to adjust the inflow to a supply chain in order to achieve a desired outflow, reducing the costs of inventory, or the goods timing in warehouses. The supply chain is modelled by a conservation law for the density of processed parts coupled to an ODE for the queue buffer occupancy. The control problem is stated as the minimization of a cost functional J measuring the queue size and the quadratic difference between the outflow and the expected one. The main novelty is the extensive use of generalized tangent vectors to a piecewise constant control, which represent time shifts of discontinuity points. Such method allows convergence results and error estimates for an Upwind- Euler steepest descent algorithm, which is also tested by numerical simulations.

preprint2012arXiv

Transport equation with nonlocal velocity in Wasserstein spaces: convergence of numerical schemes

Motivated by pedestrian modelling, we study evolution of measures in the Wasserstein space. In particular, we consider the Cauchy problem for a transport equation, where the velocity field depends on the measure itself. We deal with numerical schemes for this problem and prove convergence of a Lagrangian scheme to the solution, when the discretization parameters approach zero. We also prove convergence of an Eulerian scheme, under more strict hypotheses. Both schemes are discretizations of the push-forward formula defined by the transport equation. As a by-product, we obtain existence and uniqueness of the solution. All the results of convergence are proved with respect to the Wasserstein distance. We also show that $L^1$ spaces are not natural for such equations, since we lose uniqueness of the solution.

preprint2011arXiv

Left invertibility of discrete-time output-quantized systems: the linear case with finite inputs

This paper studies left invertibility of discrete-time linear output-quantized systems. Quantized outputs are generated according to a given partition of the state-space, while inputs are sequences on a finite alphabet. Left invertibility, i.e. injectivity of I/O map, is reduced to left D-invertibility, under suitable conditions. While left invertibility takes into account membership to sets of a given partition, left D-invertibility considers only membership to a single set, and is much easier to detect. The condition under which left invertibility and left D-invertibility are equivalent is that the elements of the dynamic matrix of the system form an algebraically independent set. Our main result is a method to compute left D-invertibility for all linear systems with no eigenvalue of modulus one. Therefore we are able to check left invertibility of output-quantized linear systems for a full measure set of matrices. Some examples are presented to show the application of the proposed method.

preprint2010arXiv

Effects of anisotropic interactions on the structure of animal groups

This paper proposes an agent-based model which reproduces different structures of animal groups. The shape and structure of the group is the effect of simple interaction rules among individuals: each animal deploys itself depending on the position of a limited number of close group mates. The proposed model is shown to produce clustered formations, as well as lines and V-like formations. The key factors which trigger the onset of different patterns are argued to be the relative strength of attraction and repulsion forces and, most important, the anisotropy in their application.

preprint2010arXiv

Multiscale modeling of granular flows with application to crowd dynamics

In this paper a new multiscale modeling technique is proposed. It relies on a recently introduced measure-theoretic approach, which allows to manage the microscopic and the macroscopic scale under a unique framework. In the resulting coupled model the two scales coexist and share information. This allows to perform numerical simulations in which the trajectories and the density of the particles affect each other. Crowd dynamics is the motivating application throughout the paper.

preprint2010arXiv

Sensor Deployment for Network-like Environments

This paper considers the problem of optimally deploying omnidirectional sensors, with potentially limited sensing radius, in a network-like environment. This model provides a compact and effective description of complex environments as well as a proper representation of road or river networks. We present a two-step procedure based on a discrete-time gradient ascent algorithm to find a local optimum for this problem. The first step performs a coarse optimization where sensors are allowed to move in the plane, to vary their sensing radius and to make use of a reduced model of the environment called collapsed network. It is made up of a finite discrete set of points, barycenters, produced by collapsing network edges. Sensors can be also clustered to reduce the complexity of this phase. The sensors' positions found in the first step are then projected on the network and used in the second finer optimization, where sensors are constrained to move only on the network. The second step can be performed on-line, in a distributed fashion, by sensors moving in the real environment, and can make use of the full network as well as of the collapsed one. The adoption of a less constrained initial optimization has the merit of reducing the negative impact of the presence of a large number of local optima. The effectiveness of the presented procedure is illustrated by a simulated deployment problem in an airport environment.

preprint2010arXiv

Time-evolving measures and macroscopic modeling of pedestrian flow

This paper deals with the early results of a new model of pedestrian flow, conceived within a measure-theoretical framework. The modeling approach consists in a discrete-time Eulerian macroscopic representation of the system via a family of measures which, pushed forward by some motion mappings, provide an estimate of the space occupancy by pedestrians at successive time steps. From the modeling point of view, this setting is particularly suitable to treat nonlocal interactions among pedestrians, obstacles, and wall boundary conditions. In addition, analysis and numerical approximation of the resulting mathematical structures, which is the main target of this work, follow more easily and straightforwardly than in case of standard hyperbolic conservation laws, also used in the specialized literature by some Authors to address analogous problems.

preprint2009arXiv

Detection of Gaussian signals via hexagonal sensor networks

This paper considers a special case of the problem of identifying a static scalar signal, depending on the location, using a planar network of sensors in a distributed fashion. Motivated by the application to monitoring wild-fires spreading and pollutants dispersion, we assume the signal to be Gaussian in space. Using a network of sensors positioned to form a regular hexagonal tessellation, we prove that each node can estimate the parameters of the Gaussian from local measurements. Moreover, we study the sensitivity of these estimates to additive errors affecting the measurements. Finally, we show how a consensus algorithm can be designed to fuse the local estimates into a shared global estimate, effectively compensating the measurement errors.

preprint2009arXiv

Modeling self-organization in pedestrians and animal groups from macroscopic and microscopic viewpoints

This paper is concerned with mathematical modeling of intelligent systems, such as human crowds and animal groups. In particular, the focus is on the emergence of different self-organized patterns from non-locality and anisotropy of the interactions among individuals. A mathematical technique by time-evolving measures is introduced to deal with both macroscopic and microscopic scales within a unified modeling framework. Then self-organization issues are investigated and numerically reproduced at the proper scale, according to the kind of agents under consideration.