The Vlasov-Poisson-Fokker-Planck system under the high field scaling describes the Brownian motion of a large system of particles in a surrounding bath where both collision and field effects (electrical or gravitational) are dominant. Numerically solving this system becomes challenging due to the stiff collision term and stiff nonlinear transport term with respect to the high field. We present a class of Asymptotic-Preserving scheme which is efficient in the high field regime, namely, large time steps and coarse meshes can be used, yet the high field limit is still captured. The idea is to combine the two stiff terms and treat them implicitly. Thanks to the linearity of the collision term, using the discretization described in [Jin S, Yan B. J. Comp. Phys., 2011, 230: 6420-6437] we only need to invert a symmetric matrix. This method can be easily extended to higher dimensions. The method is shown to be positive, stable, mass and asymptotic preserving. Numerical experiments validate its efficiency in both kinetic and high field regimes including mixing regimes.
Bibliographical noteFunding Information:
∗Received August 13, 2011. This research was partially supported by NSF grant No. DMS-0608720, and NSF FRG grant DMS-0757285. SJ was also supported by a Van Vleck Distinguished Research Prize and a Vilas Associate Award from University of Wisconsin-Madison.
- Asymptotic-preserving schemes
- High field limit
- Vlasov-Poisson-Fokker-Planck equation