ISSDE:A Monte Carlo implicit simulation code based on Stratonovich SDE approach of Coulomb collision  

作　　者：Yifeng Zheng Jianyuan Xiao Yanpeng Wang Jiangshan Zheng Ge Zhuang 郑艺峰;肖建元;王彦鹏;郑江山;庄革(School of Nuclear Science and Technology,University of Science and Technology of China,Hefei 230026,China)

机构地区：[1]School of Nuclear Science and Technology,University of Science and Technology of China,Hefei 230026,China

出　　处：《Chinese Physics B》2021年第9期297-312,共16页中国物理B（英文版）

基　　金：Project supported by the National MCF Energy R&D Program of China(Grant No.2018YFE0304100);the National Key Research and Development Program of China(Grant Nos.2016YFA0400600,2016YFA0400601,and 2016YFA0400602);the National Natural Science Foundation of China(Grant Nos.NSFC-11805273 and NSFC-11905220).

摘　　要：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.

关 键 词：Fokker-Planck equation Stratonovich SDE IMPLICIT slowing down process 

分 类 号：O211.63[理学—概率论与数理统计]