Analysis of a heterogeneous kinetic model for traffic flow

In this work we extend a recent kinetic traffic model to the case of more than one class of vehicles, each of which is characterized by few different microscopic features. We consider a Boltzmann-like framework with only binary interactions, which take place among vehicles belonging to the various classes. Our approach differs from the multi-population kinetic model based on a lattice of speeds because here we assume continuous velocity spaces and we introduce a parameter describing the physical velocity jump performed by a vehicle that increases its speed after an interaction. The model is discretized in order to investigate numerically the structure of the resulting fundamental diagrams and the system of equations is analyzed by studying well posedness. Moreover, we compute the equilibria of the discretized model and we show that the exact asymptotic kinetic distributions can be obtained with a small number of velocities in the grid. Finally, we introduce a new probability law in order to attenuate the sharp capacity drop occurring in the diagrams of traffic.Comment: 31 page

Preconditioned fully implicit PDE solvers for monument conservation

Mathematical models for the description, in a quantitative way, of the damages induced on the monuments by the action of specific pollutants are often systems of nonlinear, possibly degenerate, parabolic equations. Although some the asymptotic properties of the solutions are known, for a short window of time, one needs a numerical approximation scheme in order to have a quantitative forecast at any time of interest. In this paper a fully implicit numerical method is proposed, analyzed and numerically tested for parabolic equations of porous media type and on a systems of two PDEs that models the sulfation of marble in monuments. Due to the nonlinear nature of the underlying mathematical model, the use of a fixed point scheme is required and every step implies the solution of large, locally structured, linear systems. A special effort is devoted to the spectral analysis of the relevant matrices and to the design of appropriate iterative or multi-iterative solvers, with special attention to preconditioned Krylov methods and to multigrid procedures. Numerical experiments for the validation of the analysis complement this contribution.Comment: 26 pages, 13 figure