1,028 research outputs found
Efficient numerical calculation of drift and diffusion coefficients in the diffusion approximation of kinetic equations
In this paper we study the diffusion approximation of a swarming model given
by a system of interacting Langevin equations with nonlinear friction. The
diffusion approximation requires the calculation of the drift and diffusion
coefficients that are given as averages of solutions to appropriate Poisson
equations. We present a new numerical method for computing these coefficients
that is based on the calculation of the eigenvalues and eigenfunctions of a
Schr\"odinger operator. These theoretical results are supported by numerical
simulations showcasing the efficiency of the method
Particle based gPC methods for mean-field models of swarming with uncertainty
In this work we focus on the construction of numerical schemes for the
approximation of stochastic mean--field equations which preserve the
nonnegativity of the solution. The method here developed makes use of a
mean-field Monte Carlo method in the physical variables combined with a
generalized Polynomial Chaos (gPC) expansion in the random space. In contrast
to a direct application of stochastic-Galerkin methods, which are highly
accurate but lead to the loss of positivity, the proposed schemes are capable
to achieve high accuracy in the random space without loosing nonnegativity of
the solution. Several applications of the schemes to mean-field models of
collective behavior are reported.Comment: Communications in Computational Physics, to appea
Uncertainty quantification for kinetic models in socio-economic and life sciences
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.Comment: To appear in "Uncertainty Quantification for Hyperbolic and Kinetic
Equations
Structure preserving stochastic Galerkin methods for Fokker-Planck equations with background interactions
This paper is devoted to the construction of structure preserving stochastic
Galerkin schemes for Fokker-Planck type equations with uncertainties and
interacting with an external distribution, that we refer to as a background
distribution. The proposed methods are capable to preserve physical properties
in the approximation of statistical moments of the problem like nonnegativity,
entropy dissipation and asymptotic behaviour of the expected solution. The
introduced methods are second order accurate in the transient regimes and high
order for large times. We present applications of the developed schemes to the
case of fixed and dynamic background distribution for models of collective
behaviour
Local sensitivity analysis for the Cucker-Smale model with random inputs
We present pathwise flocking dynamics and local sensitivity analysis for the
Cucker-Smale(C-S) model with random communications and initial data. For the
deterministic communications, it is well known that the C-S model can model
emergent local and global flocking dynamics depending on initial data and
integrability of communication function. However, the communication mechanism
between agents are not a priori clear and needs to be figured out from observed
phenomena and data. Thus, uncertainty in communication is an intrinsic
component in the flocking modeling of the C-S model. In this paper, we provide
a class of admissible random uncertainties which allows us to perform the local
sensitivity analysis for flocking and establish stability to the random C-S
model with uncertain communication.Comment: 32 page
Efficient numerical calculation of drift and diffusion coefficients in the diffusion approximation of kinetic equations
In this paper we study the diffusion approximation of a swarming model given by a system of interacting Langevin equations with nonlinear friction. The diffusion approximation requires the calculation of the drift and diffusion coefficients that are given as averages of solutions to appropriate Poisson equations. We present a new numerical method for computing these coefficients that is based on the calculation of the eigenvalues and eigenfunctions of a Schr¨odinger operator. These theoretical results are supported by numerical simulations showcasing the efficiency of the method
- …