Gabriella Puppo

Gabriella Puppo, Matteo Semplice, Andrea Tosin et al.

The purpose of this paper is to study the properties of kinetic models for traffic flow described by a Boltzmann-type approach and based on a continuous space of microscopic velocities. In our models, the particular structure of the collision kernel allows one to find the analytical expression of a class of steady-state distributions, which are characterized by being supported on a quantized space of microscopic speeds. The number of these velocities is determined by a physical parameter describing the typical acceleration of a vehicle and the uniqueness of this class of solutions is supported by numerical investigations. This shows that it is possible to have the full richness of a kinetic approach with the simplicity of a space of microscopic velocities characterized by a small number of modes. Moreover, the explicit expression of the asymptotic distribution paves the way to deriving new macroscopic equations using the closure provided by the kinetic model.

Gabriella Puppo, Matteo Semplice, Andrea Tosin et al.

Experimental studies on vehicular traffic provide data on quantities like density, flux, and mean speed of the vehicles. However, the diagrams relating these variables (the fundamental and speed diagrams) show some peculiarities not yet fully reproduced nor explained by mathematical models. In this paper, resting on the methods of kinetic theory, we introduce a new traffic model which takes into account the heterogeneous nature of the flow of vehicles along a road. In more detail, the model considers traffic as a mixture of two or more populations of vehicles (e.g., cars and trucks) with different microscopic characteristics, in particular different lengths and/or maximum speeds. With this approach we gain some insights into the scattering of the data in the regime of congested traffic clearly shown by actual measurements.

Isabella Cravero, Gabriella Puppo, Matteo Semplice et al.

This work is dedicated to the development and comparison of WENO-type reconstructions for hyperbolic systems of balance laws. We are particularly interested in high order shock capturing non-oscillatory schemes with uniform accuracy within each cell and low spurious effects. We need therefore to develop a tool to measure the artifacts introduced by a numerical scheme. To this end, we study the deformation of a single Fourier mode and introduce the notion of distorsive errors, which measure the amplitude of the spurious modes created by a discrete derivative operator. Further we refine this notion with the idea of temperature, in which the amplitude of the spurious modes is weighted with its distance in frequency space from the exact mode. Following this approach linear schemes have zero temperature, but to prevent oscillations it is necessary to introduce nonlinearities in the scheme, thus increasing their temperature. However it is important to heat the linear scheme just enough to prevent spurious oscillations. With several tests we show that the newly introduced CWENOZ schemes are cooler than other existing WENO-type operators, while maintaining good non-oscillatory properties.

Gabriella Puppo, Matteo Semplice, Andrea Tosin et al.

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.

Gabriella Puppo

This minicourse contains a description of recent results on the modelling of rarefied gases in weakly non equilibrium regimes, and the numerical methods used to approximate the resulting equations. Therefore this work focuses on BGK type approximations, rather than on full Boltzmann models. Within this framework, models for polyatomic gases and for mixtures will be considered. We will also address numerical issues characteristic of the difficulties one encounters when integrating kinetic equations. In particular, we will consider asymptotic preserving schemes, which are designed to approximate equilibrium solutions, without resolving the fast scales of the approach to equilibrium.

Giuseppe Visconti, Michael Herty, Gabriella Puppo et al.

Starting from interaction rules based on two levels of stochasticity we study the influence of the microscopic dynamics on the macroscopic properties of vehicular flow. In particular, we study the qualitative structure of the resulting flux-density and speed-density diagrams for different choices of the desired speeds. We are able to recover multivalued diagrams as a result of the existence of a one-parameter family of stationary distributions, whose expression is analytically found by means of a Fokker-Planck approximation of the initial Boltzmann-type model.

Fausto Cavalli, Giovanni Naldi, Gabriella Puppo et al.

In this work we compare two semidiscrete schemes for the solution of hyperbolic conservation laws, namely the relaxation and the Kurganov Tadmor central scheme. We are particularly interested in their behavior under small time steps, in view of future applications to convection diffusion problems. The schemes are tested on two benchmark problems, with one space variable.

Fausto Cavalli, Giovanni Naldi, Gabriella Puppo et al.

The application of Runge-Kutta schemes designed to enjoy a large region of absolute stability can significantly increase the efficiency of numerical methods for PDEs based on a method of lines approach. In this work we investigate the improvement in the efficiency of the time integration of relaxation schemes for degenerate diffusion problems, using SSP Runge-Kutta schemes and computing the maximal CFL coefficients. This technique can be extended to other PDEs, linear and nonlinear, provided the space operator has eigenvalues with a non-zero real part.

Fausto Cavalli, Giovanni Naldi, Gabriella Puppo et al.

Several relaxation approximations to partial differential equations have been recently proposed. Examples include conservation laws, Hamilton-Jacobi equations, convection-diffusion problems, gas dynamics problems. The present paper focuses onto diffusive relaxation schemes for the numerical approximation of nonlinear parabolic equations. These schemes are based on suitable semilinear hyperbolic system with relaxation terms. High order methods are obtained by coupling ENO and WENO schemes for space discretization with IMEX schemes for time integration. Error estimates and convergence analysis are developed for semidiscrete schemes with numerical analysis for fully discrete relaxed schemes. Various numerical results in one and two dimension illustrate the high accuracy and good properties of the proposed numerical schemes. These schemes can be easily implemented for parallel computer and applied to more general system of nonlinear parabolic equations in two- and three-dimensional cases.