Variational Asymptotic Preserving Scheme for the Vlasov-Poisson-Fokker-Planck System

by   José A. Carrillo, et al.

We design a variational asymptotic preserving scheme for the Vlasov-Poisson-Fokker-Planck system with the high field scaling, which describes the Brownian motion of a large system of particles in a surrounding bath. Our scheme builds on an implicit-explicit framework, wherein the stiff terms coming from the collision and field effects are solved implicitly while the convection terms are solved explicitly. To treat the implicit part, we propose a variational approach by viewing it as a Wasserstein gradient flow of the relative entropy, and solve it via a proximal quasi-Newton method. In so doing we get positivity and asymptotic preservation for free. The method is also massively parallelizable and thus suitable for high dimensional problems. We further show that the convergence of our implicit solver is uniform across different scales. A suite of numerical examples are presented at the end to validate the performance of the proposed scheme.



There are no comments yet.


page 1

page 2

page 3

page 4


An asymptotic-preserving IMEX method for nonlinear radiative transfer equation

We present an asymptotic preserving method for the radiative transfer eq...

Unconditional positivity-preserving and energy stable schemes for a reduced Poisson-Nernst-Planck system

The Poisson-Nernst-Planck (PNP) system is a widely accepted model for si...

A positivity-preserving and convergent numerical scheme for the binary fluid-surfactant system

In this paper, we develop a first order (in time) numerical scheme for t...

A Variational Integrator for the Discrete Element Method

A novel implicit integration scheme for the Discrete Element Method (DEM...

Drift-preserving numerical integrators for stochastic Poisson systems

We perform a numerical analysis of randomly perturbed Poisson systems. F...

An asymptotic preserving scheme for Lévy-Fokker-Planck equation with fractional diffusion limit

In this paper, we develop a numerical method for the Lévy-Fokker-Planck ...

The computational framework for continuum-kinematics-inspired peridynamics

Peridynamics (PD) is a non-local continuum formulation. The original ver...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.