Numerical modeling of the thermal force in a plasma for test-ion transport simulation based on a Monte Carlo Binary Collision Model (II) - Thermal forces due to temperature gradients parallel and perpendicular to the magnetic field

Yuki Homma, Akiyoshi Hatayama

Research output: Contribution to journalArticle

9 Citations (Scopus)

Abstract

We have developed a numerical model of the thermal force for test-ion transport simulation in a magnetized background plasma, based on the Monte Carlo Binary Collision Model (BCM) [T. Takizuka, H. Abe, J. Comput. Phys. 25 (1977) 205]. The model is basically the same as presented in our previous paper [Y. Homma, A. Hatayama, J. Comput. Phys. 231 (2012) 3211-3227] for the case without magnetic field, but in the present paper, a more extended form of a distorted Maxwellian distribution is employed for the velocity distribution of background plasma ions to simulate the thermal force caused by parallel and perpendicular (with respect to the direction of magnetic field) temperature gradients. The model consists mainly of two steps: (i) choosing a background plasma ion velocity from a distorted Maxwellian distribution, and (ii) calculating a Coulomb collision between a test particle and the above chosen ion by using the BCM. In addition, equations of motion for charged test particle in the magnetic field are calculated by Buneman-Boris Algorithm.A series of test simulations has been done in a simple geometry with different temperature gradients and different strengths of magnetic field. Numerical results of the thermal force due to parallel and perpendicular temperature gradients agree well with the theoretical prediction for all test cases. Especially, it has been confirmed that the model reproduces the temperature screening effect of test particles (i.e. guiding center drift of test particle caused by the thermal force of perpendicular temperature gradient).

Original languageEnglish
Pages (from-to)206-223
Number of pages18
JournalJournal of Computational Physics
Volume250
DOIs
Publication statusPublished - 2013 Oct 1

Fingerprint

Thermal gradients
temperature gradients
Magnetic fields
Plasmas
collisions
Ions
magnetic fields
ions
simulation
Maxwell-Boltzmann density function
Velocity distribution
Coulomb collisions
Equations of motion
Numerical models
Screening
Hot Temperature
Geometry
equations of motion
screening
velocity distribution

Keywords

  • Distorted Maxwellian
  • Monte Carlo Binary Collision Model
  • Numerical model
  • Perpendicular temperature gradient
  • Thermal force

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy (miscellaneous)

Cite this

@article{c124a42a355f493690200c329d3c9169,
title = "Numerical modeling of the thermal force in a plasma for test-ion transport simulation based on a Monte Carlo Binary Collision Model (II) - Thermal forces due to temperature gradients parallel and perpendicular to the magnetic field",
abstract = "We have developed a numerical model of the thermal force for test-ion transport simulation in a magnetized background plasma, based on the Monte Carlo Binary Collision Model (BCM) [T. Takizuka, H. Abe, J. Comput. Phys. 25 (1977) 205]. The model is basically the same as presented in our previous paper [Y. Homma, A. Hatayama, J. Comput. Phys. 231 (2012) 3211-3227] for the case without magnetic field, but in the present paper, a more extended form of a distorted Maxwellian distribution is employed for the velocity distribution of background plasma ions to simulate the thermal force caused by parallel and perpendicular (with respect to the direction of magnetic field) temperature gradients. The model consists mainly of two steps: (i) choosing a background plasma ion velocity from a distorted Maxwellian distribution, and (ii) calculating a Coulomb collision between a test particle and the above chosen ion by using the BCM. In addition, equations of motion for charged test particle in the magnetic field are calculated by Buneman-Boris Algorithm.A series of test simulations has been done in a simple geometry with different temperature gradients and different strengths of magnetic field. Numerical results of the thermal force due to parallel and perpendicular temperature gradients agree well with the theoretical prediction for all test cases. Especially, it has been confirmed that the model reproduces the temperature screening effect of test particles (i.e. guiding center drift of test particle caused by the thermal force of perpendicular temperature gradient).",
keywords = "Distorted Maxwellian, Monte Carlo Binary Collision Model, Numerical model, Perpendicular temperature gradient, Thermal force",
author = "Yuki Homma and Akiyoshi Hatayama",
year = "2013",
month = "10",
day = "1",
doi = "10.1016/j.jcp.2013.04.039",
language = "English",
volume = "250",
pages = "206--223",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",

}

TY - JOUR

T1 - Numerical modeling of the thermal force in a plasma for test-ion transport simulation based on a Monte Carlo Binary Collision Model (II) - Thermal forces due to temperature gradients parallel and perpendicular to the magnetic field

AU - Homma, Yuki

AU - Hatayama, Akiyoshi

PY - 2013/10/1

Y1 - 2013/10/1

N2 - We have developed a numerical model of the thermal force for test-ion transport simulation in a magnetized background plasma, based on the Monte Carlo Binary Collision Model (BCM) [T. Takizuka, H. Abe, J. Comput. Phys. 25 (1977) 205]. The model is basically the same as presented in our previous paper [Y. Homma, A. Hatayama, J. Comput. Phys. 231 (2012) 3211-3227] for the case without magnetic field, but in the present paper, a more extended form of a distorted Maxwellian distribution is employed for the velocity distribution of background plasma ions to simulate the thermal force caused by parallel and perpendicular (with respect to the direction of magnetic field) temperature gradients. The model consists mainly of two steps: (i) choosing a background plasma ion velocity from a distorted Maxwellian distribution, and (ii) calculating a Coulomb collision between a test particle and the above chosen ion by using the BCM. In addition, equations of motion for charged test particle in the magnetic field are calculated by Buneman-Boris Algorithm.A series of test simulations has been done in a simple geometry with different temperature gradients and different strengths of magnetic field. Numerical results of the thermal force due to parallel and perpendicular temperature gradients agree well with the theoretical prediction for all test cases. Especially, it has been confirmed that the model reproduces the temperature screening effect of test particles (i.e. guiding center drift of test particle caused by the thermal force of perpendicular temperature gradient).

AB - We have developed a numerical model of the thermal force for test-ion transport simulation in a magnetized background plasma, based on the Monte Carlo Binary Collision Model (BCM) [T. Takizuka, H. Abe, J. Comput. Phys. 25 (1977) 205]. The model is basically the same as presented in our previous paper [Y. Homma, A. Hatayama, J. Comput. Phys. 231 (2012) 3211-3227] for the case without magnetic field, but in the present paper, a more extended form of a distorted Maxwellian distribution is employed for the velocity distribution of background plasma ions to simulate the thermal force caused by parallel and perpendicular (with respect to the direction of magnetic field) temperature gradients. The model consists mainly of two steps: (i) choosing a background plasma ion velocity from a distorted Maxwellian distribution, and (ii) calculating a Coulomb collision between a test particle and the above chosen ion by using the BCM. In addition, equations of motion for charged test particle in the magnetic field are calculated by Buneman-Boris Algorithm.A series of test simulations has been done in a simple geometry with different temperature gradients and different strengths of magnetic field. Numerical results of the thermal force due to parallel and perpendicular temperature gradients agree well with the theoretical prediction for all test cases. Especially, it has been confirmed that the model reproduces the temperature screening effect of test particles (i.e. guiding center drift of test particle caused by the thermal force of perpendicular temperature gradient).

KW - Distorted Maxwellian

KW - Monte Carlo Binary Collision Model

KW - Numerical model

KW - Perpendicular temperature gradient

KW - Thermal force

UR - http://www.scopus.com/inward/record.url?scp=84879144171&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84879144171&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2013.04.039

DO - 10.1016/j.jcp.2013.04.039

M3 - Article

AN - SCOPUS:84879144171

VL - 250

SP - 206

EP - 223

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -