ISSDE: A Monte Carlo implicit simulation code based on Stratonovich SDE approach of Coulomb collision |
| |
| Affiliation: | School of Nuclear Science and Technology, University of Science and Technology of China, Hefei 230026, China |
| |
| Abstract: | A Monte Carlo implicit simulation program, Implicit Stratonovich Stochastic Differential Equations (ISSDE), is developed for solving stochastic differential equations (SDEs) that describe plasmas with Coulomb collision. The basic idea of the program is the stochastic equivalence between the Fokker-Planck equation and the Stratonovich SDEs. The splitting method is used to increase the numerical stability of the algorithm for dynamics of charged particles with Coulomb collision. The cases of Lorentzian plasma, Maxwellian plasma and arbitrary distribution function of background plasma have been considered. The adoption of the implicit midpoint method guarantees exactly the energy conservation for the diffusion term and thus improves the numerical stability compared with conventional Runge-Kutta methods. ISSDE is built with C++ and has standard interfaces and extensible modules. The slowing down processes of electron beams in unmagnetized plasma and relaxation process in magnetized plasma are studied using the ISSDE, which shows its correctness and reliability. |
| |
| Keywords: | Fokker-Planck equation Stratonovich SDE implicit slowing down process |
| 本文献已被 CNKI 等数据库收录！ |
| 点击此处可从《中国物理 B》浏览原始摘要信息 |
|
点击此处可从《中国物理 B》下载免费的PDF全文 |
|
