An entropy based thermalization scheme for hybrid simulations of Coulomb collisions

JOURNAL OF COMPUTATIONAL PHYSICS

卷:273

Article

Ricketson, L. F.¹  Rosin, M. S.¹  Caflisch, R. E.¹  Dimits, A. M.² 

[1] Univ Calif Los Angeles, Dept Math, Los Angeles, CA 90036 USA

[2] Lawrence Livermore Natl Lab, Livermore, CA 94511 USA

关键词: Plasma;    Coulomb collisions;    Particle collisions;    Monte Carlo;    Hybrid;    Entropy;   

DOI  :  10.1016/j.jcp.2014.04.059

来源: Elsevier

【 摘 要 】

We formulate and test a hybrid fluid-Monte Carlo scheme for the treatment of elastic collisions in gases and plasmas. While our primary focus and demonstrations of applicability are for moderately collisional plasmas, as described by the Landau-Fokker-Planck equation, the method is expected to be applicable also to collision processes described by the Boltzmann equation. This scheme is similar to the previously discussed velocity-based scheme (R. Caflisch et al. (2008) [7]) and the scattering-angle-based scheme (A.M. Dimits et al. (2010) [14])], but with a firmer theoretical basis and without the inherent limitation to the Landau-Fokker-Planck case. It gives a significant performance improvement (e.g., error for a given computational effort) over the velocity-based scheme. These features are achieved by assigning passive scalars to each simulated particle and tracking their evolution through collisions. The method permits a formal error analysis that agrees with numerical results. The tests performed are for the evolution of an anisotropic Maxwellian and a bump-on-tail distribution. (C) 2014 Elsevier Inc. All rights reserved.

【 授权许可】

Free