Mattia Zanella

Lorenzo Pareschi, Mattia Zanella

In this paper we focus on the construction of numerical schemes for nonlinear Fokker-Planck equations that preserve the structural properties, like non negativity of the solution, entropy dissipation and large time behavior. The methods here developed are second order accurate, they do not require any restriction on the mesh size and are capable to capture the asymptotic steady states with arbitrary accuracy. These properties are essential for a correct description of the underlying physical problem. Applications of the schemes to several nonlinear Fokker-Planck equations with nonlocal terms describing emerging collective behavior in socio-economic and life sciences are presented.

Giacomo Albi, Lorenzo Pareschi, Mattia Zanella

In this paper we consider the modeling of opinion dynamics over time dependent large scale networks. A kinetic description of the agents' distribution over the evolving network is considered which combines an opinion update based on binary interactions between agents with a dynamic creation and removal process of new connections. The number of connections of each agent influences the spreading of opinions in the network but also the way connections are created is influenced by the agents' opinion. The evolution of the network of connections is studied by showing that its asymptotic behavior is consistent both with Poisson distributions and truncated power-laws. In order to study the large time behavior of the opinion dynamics a mean field description is derived which allows to compute exact stationary solutions in some simplified situations. Numerical methods which are capable to describe correctly the large time behavior of the system are also introduced and discussed. Finally, several numerical examples showing the influence of the agents' number of connections in the opinion dynamics are reported.

Giacomo Albi, Lorenzo Pareschi, Mattia Zanella

In this paper the optimal control of flocking models with random inputs is investigated from a numerical point of view. The effect of uncertainty in the interaction parameters is studied for a Cucker-Smale type model using a generalized polynomial chaos (gPC) approach. Numerical evidence of threshold effects in the alignment dynamic due to the random parameters is given. The use of a selective model predictive control permits to steer the system towards the desired state even in unstable regimes.

Lorenzo Pareschi, Pierluigi Vellucci, Mattia Zanella

We introduce and discuss kinetic models describing the influence of the competence in the evolution of decisions in a multi-agent system. The original exchange mechanism, which is based on the human tendency to compromise and change opinion through self-thinking, is here modified to include the role of the agents' competence. In particular, we take into account the agents' tendency to behave in the same way as if they were as good, or as bad, as their partner: the so-called equality bias. This occurred in a situation where a wide gap separated the competence of group members. We discuss the main properties of the kinetic models and numerically investigate some examples of collective decision under the influence of the equality bias. The results confirm that the equality bias leads the group to suboptimal decisions.

Andrea Tosin, Mattia Zanella

We develop a hierarchical description of traffic flow control by means of driver-assist vehicles aimed at the mitigation of speed-dependent road risk factors. Microscopic feedback control strategies are designed at the level of vehicle-to-vehicle interactions and then upscaled to the global flow via a kinetic approach based on a Boltzmann-type equation. Then first and second order hydrodynamic traffic models, which naturally embed the microscopic control strategies, are consistently derived from the kinetic-controlled framework via suitable closure methods. Several numerical examples illustrate the effectiveness of such a hierarchical approach at the various scales.

Michael Herty, Andrea Tosin, Giuseppe Visconti et al.

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.

Lorenzo Pareschi, Mattia Zanella

In this paper we consider the development of numerical schemes for mean-field equations describing the collective behavior of a large group of interacting agents. The schemes are based on a generalization of the classical Chang-Cooper approach and are capable to preserve the main structural properties of the systems, namely nonnegativity of the solution, physical conservation laws, entropy dissipation and stationary solutions. In particular, the methods here derived are second order accurate in transient regimes whereas they can reach arbitrary accuracy asymptotically for large times. Several examples are reported to show the generality of the approach.

Horacio Tettamanti, Giulia Guicciardi, Mattia Zanella

In this work, we present a multiscale kinetic framework for consensus-based image segmentation. By interpreting an image as a system of interacting particles, each pixel is characterised by its spatial position and an internal feature encoding color information. We introduce a coupled interaction scheme governing the evolution of particles in both position and feature spaces, from which we derive a kinetic formulation for the particle density in the space-feature domain combining transport, aggregation, and diffusion effects. Furthermore, through a suitable scaling, we obtain a first-order macroscopic model describing the evolution of the fraction of pixels carrying information on the fraction of pixels having a certain feature. Based on this reduced-complexity model, we present a data-oriented approach where we make use of particle-based optimisation techniques for the accurate segmentation of images. Numerical tests show the effectiveness of the proposed framework and its robustness under different noise conditions.

Raffaella Fiamma Cabini, Anna Pichiecchio, Alessandro Lascialfari et al.

In this work, we apply a kinetic version of a bounded confidence consensus model to biomedical segmentation problems. In the presented approach, time-dependent information on the microscopic state of each particle/pixel includes its space position and a feature representing a static characteristic of the system, i.e. the gray level of each pixel. From the introduced microscopic model we derive a kinetic formulation of the model. The large time behavior of the system is then computed with the aid of a surrogate Fokker-Planck approach that can be obtained in the quasi-invariant scaling. We exploit the computational efficiency of direct simulation Monte Carlo methods for the obtained Boltzmann-type description of the problem for parameter identification tasks. Based on a suitable loss function measuring the distance between the ground truth segmentation mask and the evaluated mask, we minimize the introduced segmentation metric for a relevant set of 2D gray-scale images. Applications to biomedical segmentation concentrate on different imaging research contexts.

Monica Caloi, Mattia Zanella

Inspired by recently developed Fokker--Planck models for Bose--Einstein statistics, we study a consensus formation model with condensation effects driven by a polynomial diffusion coefficient vanishing at the domain boundaries. For the underlying kinetic model, given by a nonlinear Fokker--Planck equation with superlinear drift, it was shown that if the initial mass exceeds a critical threshold, the solution may exhibit finite-time concentration in certain parameter regimes. Here, we show that this supercritical mass phenomenon persists for a broader class of diffusion functions and provide estimates of the critical mass required to induce finite-time loss of regularity.

Giacomo Dimarco, Lorenzo Pareschi, Mattia Zanella

Kinetic equations play a major rule in modeling large systems of interacting particles. Recently the legacy of classical kinetic theory found novel applications in socio-economic and life sciences, where processes characterized by large groups of agents exhibit spontaneous emergence of social structures. Well-known examples are the formation of clusters in opinion dynamics, the appearance of inequalities in wealth distributions, flocking and milling behaviors in swarming models, synchronization phenomena in biological systems and lane formation in pedestrian traffic. The construction of kinetic models describing the above processes, however, has to face the difficulty of the lack of fundamental principles since physical forces are replaced by empirical social forces. These empirical forces are typically constructed with the aim to reproduce qualitatively the observed system behaviors, like the emergence of social structures, and are at best known in terms of statistical information of the modeling parameters. For this reason the presence of random inputs characterizing the parameters uncertainty should be considered as an essential feature in the modeling process. In this survey we introduce several examples of such kinetic models, that are mathematically described by nonlinear Vlasov and Fokker--Planck equations, and present different numerical approaches for uncertainty quantification which preserve the main features of the kinetic solution.