Abstract
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.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.