G. Dimarco, G. Toscani, M. Zanella
Preprint arXiv, 2021.
The spread of COVID-19 has been thwarted in most countries through non-pharmaceutical interventions. In particular, the most effective measures in this direction have been the stay-at-home and closure strategies of businesses and schools. However, population-wide lockdowns are far from being optimal carrying heavy economic consequences. Therefore, there is nowadays a strong interest in designing more efficient restrictions. In this work, starting from a recent kinetic-type model which takes into account the heterogeneity described by the social contact of individuals, we analyze the effects of introducing an optimal control strategy into the system, to limit selectively the mean number of contacts and reduce consequently the number of infected cases. Thanks to a data-driven approach, we show that this new mathematical model permits to assess
the effects of the social limitations. Finally, using the model introduced here and starting from the available data, we show the effectivity of the proposed selective measures to dampen the epidemic trends
G. Toscani, M. Zanella
Atti Accad. Naz. Lincei Rend. Lincei Mat. Appl., to appear. Preprint arXiv.
We study the relaxation to equilibrium for a class linear one-dimensional Fokker–Planck equations characterized by a particular subcritical confinement potential.
An interesting feature of this class of Fokker–Planck equations is that, for any given probability density $e(x)$, the diffusion coefficient can be built to have $e(x)$ as steady state. This representation of the equilibrium density can be fruitfully used to obtain one-dimensional Wirtinger-type inequalities and to recover, for a sufficiently regular density $e(x) $, a polynomial rate of convergence to equilibrium.Numerical results then confirm the theoretical analysis, and allow to conjecture that convergence to equilibrium with positive rate still holds for steady states characterized by a very slow polynomial decay at infinity.
M. Azzi, C. Bardelli, S. Deandrea, G. Dimarco, S. Figini, P. Perotti, G. Toscani, M. Zanella
Math. Mod. Meth. Appl. Sci., to appear. Preprint arXiv.
In this work, using a detailed dataset furnished by National Health Authorities concerning the Province of Pavia (Lombardy, Italy), we propose to determine the essential features of the ongoing COVID-19 pandemic in term of contact dynamics. Our contribution is devoted to provide a possible planning of the needs of medical infrastructures in the Pavia Province and to suggest different scenarios about the vaccination campaign which possibly help in reducing the fatalities and/or reducing the number of infected in the population.
The proposed research combines a new mathematical description of the spread of an infectious diseases which takes into account both age and average daily social contacts with a detailed analysis of the dataset of all traced infected individuals in the Province of Pavia. These information are used to develop a data-driven model in which calibration and feeding of the model are extensively used. The epidemiological evolution is obtained by relying on an approach based on statical mechanics. This leads to study the evolution over time of a system of probability distributions characterizing the age and social contacts of the population. One of the main outcomes shows that, as expected, the spread of the disease is closely related to the mean number of contacts of individuals. The model permits to forecast thanks to an uncertainty quantification approach and in the short time horizon, the average number and the confidence bands of expected hospitalized classified by age and to test different options for an effective vaccination campaign with age-decreasing priority.
L. Pareschi, T. Trimborn, M. Zanella
International Journal for Uncertainty Quantification, to appear. Preprint arXiv, 2021.
In this paper, we extend a recently introduced multi-fidelity control variate for the uncertainty quantification of the Boltzmann equation to the case of kinetic models arising in the study of multiagent systems. For these phenomena, where the effect of uncertainties is particularly evident, several models have been developed whose equilibrium states are typically unknown. In particular, we aim to develop efficient numerical methods based on solving the kinetic equations in the phase space by Direct Simulation Monte Carlo (DSMC) coupled to a Monte Carlo sampling in the random space. To this end, exploiting the knowledge of the corresponding mean-field approximation we develop novel mean-field Control Variate (MFCV) methods that are able to strongly reduce the variance of the standard Monte Carlo sampling method in the random space. We verify these observations with several numerical examples based on classical models , including wealth exchanges and opinion formation model for collective phenomena.
G. Dimarco, A. Tosin, M. Zanella
Preprint arXiv, 2021.
In this paper, we derive second order hydrodynamic traffic models from kinetic-controlled equations for driver-assist vehicles. At the vehicle level we take into account two main control strategies synthesising the action of adaptive cruise controls and cooperative adaptive cruise controls. The resulting macroscopic dynamics fulfil the anisotropy condition introduced in the celebrated Aw-Rascle-Zhang model. Unlike other models based on heuristic arguments, our approach unveils the main physical aspects behind frequently used hydrodynamic traffic models and justifies the structure of the resulting macroscopic equations incorporating driver-assist vehicles. Numerical insights show that the presence of driver-assist vehicles produces an aggregate homogenisation of the mean flow speed, hich may also be steered towards a suitable desired speed in such a way that optimal flows and traffic stabilisation are reached
R. Borsche, A. Klar, M. Zanella
Preprint arXiv, 2020
We study the application of a recently introduced hierarchical description of traffic flow control by driver-assist vehicles to include lane changing dynamics. Lane-dependent feedback control strategies are implemented at the level of vehicles and the aggregate trends are studied by means of Boltzmann-type equations determining three different hydrodynamics based on the lane switching frequency. System of first order macroscopic equations describing the evolution of densities along the lanes are then consistently determined through a suitable closure strategy. Numerical examples are then presented to illustrate the features of the proposed hierarchical approach.
G. Dimarco, B. Perthame, G. Toscani, M. Zanella
Journal of Mathematical Biology, 83, 4, 2021. (Preprint arXiv)
We introduce a mathematical description of the impact of sociality in the spread of infectious diseases by integrating an epidemiological dynamics with a kinetic modeling of population-based contacts. The kinetic description leads to study the evolution over time of Boltzmann-type equa- tions describing the number densities of social contacts of susceptible, infected and recovered indi- viduals, whose proportions are driven by a classical SIR-type compartmental model in epidemiology. Explicit calculations show that the spread of the disease is closely related to moments of the con- tact distribution. Furthermore, the kinetic model allows to clarify how a selective control can be assumed to achieve a minimal lockdown strategy by only reducing individuals undergoing a very large number of daily contacts. We conduct numerical simulations which confirm the ability of the model to describe different phenomena characteristic of the rapid spread of an epidemic. Motiv- ated by the COVID-19 pandemic, a last part is dedicated to fit numerical solutions of the proposed model with infection data coming from different European countries.
L. Preziosi, G. Toscani, M. Zanella
Journal of Theoretical Biology, 514: 110579, 2021. (Preprint arXiv)
The mathematical modeling of tumor growth has a long history, and has been mathematically formulated in several different ways. Here we tackle the problem in the case of a continuous distribution using mathematical tools from statistical physics. To this extent, we introduce a novel kinetic model of growth which highlights the role of microscopic transitions in determining a variety of equilibrium distributions. At variance with other approaches, the mesoscopic description in terms of elementary interactions allows to design precise microscopic feedback control therapies, able to influence the natural tumor growth and to mitigate the risk factors involved in big sized tumors. We further show that under a suitable scaling both the free and controlled growth models correspond to Fokker–Planck type equations for the growth distribution with variable coefficients of diffusion and drift, whose steady solutions in the free case are given by a class of generalized Gamma densities which can be characterized by fat tails. In this scaling the feedback control produces an explicit modification of the drift operator, which is shown to strongly modify the emerging distribution for the tumor size. In particular, the size distributions in presence of therapies manifest slim tails in all growth models, which corresponds to a marked mitigation of the risk factors. Numerical results confirming the theoretical analysis are also presented.
G. Albi, L. Pareschi, M. Zanella
Mathematical Biosciences and Engineering, 18(6): 7161-7190, 2021. (Preprint arXiv)
After the introduction of drastic containment measures aimed at stopping the epidemic contagion from SARS-CoV2, many governments have adopted a strategy based on a periodic relaxation of such measures in the face of a severe economic crisis caused by lockdowns. Assessing the impact of such openings in relation to the risk of a resumption of the spread of the disease is an extremely difficult problem due to the many unknowns concerning the actual number of people infected, the actual reproduction number and infection fatality rate of the disease. In this work, starting from a compartmental model with a social structure and stochastic inputs, we derive models with multiple feedback controls depending on the social activities that allow to assess the impact of a selective relaxation of the containment measures in the presence of uncertain data. Specific contact patterns in the home, work, school and other locations have been considered. Results from different scenarios concerning the first wave of the epidemic in some major countries, including Germany, France, Italy, Spain, the United Kingdom and the United States, are presented and discussed.
G. Albi, L. Pareschi, M. Zanella
Journal of Mathematical Biology, 82, 63, 2021. (Preprint arXiv)
The adoption of containment measures to reduce the amplitude of the epidemic peak is a key aspect in tackling the rapid spread of an epidemic. Classical compartmental models must be modified and studied to correctly describe the effects of forced external actions to reduce the impact of the disease. In addition, data are often incomplete and heterogeneous, so a high degree of uncertainty must naturally be incorporated into the models. In this work we address both these aspects, through an optimal control formulation of the epidemiological model in presence of uncertain data. After the introduction of the optimal control problem, we formulate an instantaneous approximation of the control that allows us to derive new feedback controlled compartmental models capable of describing the epidemic peak reduction. The need for long-term interventions shows that alternative actions based on the social structure of the system can be as effective as the more expensive global strategy. The importance of the timing and intensity of interventions is particularly relevant in the case of uncertain parameters on the actual number of infected people. Simulations related to data from the recent COVID-19 outbreak in Italy are presented and discussed.