An adaptive, implicit, conservative, 1D-2V multi-species Vlasov-Fokker-Planck multi-scale solver in planar geometry
Abstract
We consider a 1D-2V Vlasov-Fokker-Planck multi-species ionic description coupled to fluid electrons. We address temporal stiffness with implicit time stepping, suitably preconditioned. To address temperature disparity in time and space, we extend the conservative adaptive velocity-space discretization scheme proposed in [Taitano et al., J. Comput. Phys., 318, 391-420, (2016)] to a spatially inhomogeneous system. In this approach, we normalize the velocity-space coordinate to a temporally and spatially varying local characteristic speed per species. We explicitly consider the resulting inertial terms in the Vlasov equation, and derive a discrete formulation that conserves mass, momentum, and energy up to a prescribed nonlinear tolerance upon convergence. Our conservation strategy employs nonlinear constraints to enforce these properties discretely for both the Vlasov operator and the Fokker-Planck collision operator. Numerical examples of varying degrees of complexity, including shock-wave propagation, demonstrate the favorable efficiency and accuracy properties of the scheme.
- Publication:
-
Journal of Computational Physics
- Pub Date:
- July 2018
- DOI:
- 10.1016/j.jcp.2018.03.007
- arXiv:
- arXiv:1709.04502
- Bibcode:
- 2018JCoPh.365..173T
- Keywords:
-
- Conservative discretization;
- Thermal velocity based adaptive grid;
- 1D2V;
- Fokker-Planck;
- Rosenbluth potentials;
- Physics - Plasma Physics
- E-Print:
- 38 pages, 18 figures, pre-print submitted to J. Comp. Physics