Figures 2 and 3 correspond to different simulation parameters than those stated. To reproduce Figs. 2 and 3 the rotor and ambient density are as stated, but the pressure should be set to 0.5, the ratio of specific heats should be set to 5/3, the rotor's toroidal velocity at a radius of 0.1 should be
A high-order finite-volume algorithm for Fokker–Planck collisions in magnetized plasmas
✍ Scribed by Z. Xiong; R.H. Cohen; T.D. Rognlien; X.Q. Xu
- Publisher
- Elsevier Science
- Year
- 2008
- Tongue
- English
- Weight
- 295 KB
- Volume
- 227
- Category
- Article
- ISSN
- 0021-9991
No coin nor oath required. For personal study only.
✦ Synopsis
A high-order finite-volume algorithm is developed for the Fokker-Planck Operator (FPO) describing Coulomb collisions in strongly magnetized plasmas. The algorithm uses a generic fourth-order reconstruction scheme on an unstructured grid in the velocity space spanned by parallel velocity and magnetic moment. By analytically mapping between different coordinates, it produces an accurate and density-conserving numerical FPO for an arbitrary choice of velocity space coordinates. A linearized FPO in constants-of-motion coordinates is implemented as an example of the present algorithm combined with a cut-cell merging procedure. Numerical tests include the thermalization of a test distribution with a background Maxwellian at a different temperature, and the return to isotropy for a distribution initialized with a velocity space loss-cone. Utilization of the method for a nonlinear FPO is straightforward but requires evaluation of the Trubnikov-Rosenbluth potentials.
📜 SIMILAR VOLUMES