ISSDE: First-principles Implicit Simulations Based on Stratonovich SDE Approach of Coulumb Collision

doi:10.19596/j.cnki.1001-246x.7843

Abstract

Abstract: A first-principles implicit simulation program, Implicit Stratonovich Stochastic Differential Equations (ISSDE), is constructed for stochastic differential equations which describes plasmas with Coulomb collision. Basic idea of the program is stochastic equivalence between Fokker-Planck equation and Stratonovich SDE. Implicit discrete method guarantees numerical stability and conservation of kinetic energy. ISSDE is built with C++ language, and is designed to possess standard interfaces and extendible modules. Slowing-down processes of electron beams in both unmagnetized and magnetized plasmas are studied, which shows correctness of ISSDE. It provides a powerful tool for collisional plasmas studies.

Key words: FP equation, Stratonovich SDE, implicit, slowing-down

CLC Number:

ZHENG Yifeng, WANG Yulei, LIU Jian, QIN Hong. ISSDE: First-principles Implicit Simulations Based on Stratonovich SDE Approach of Coulumb Collision[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2019, 36(3): 265-279.

References

[1] PERONA A, BORGOGNO D, ERIKSSON L G. A test electron model for the study of three dimensional magnetic reconnection effects[J]. Computer Physics Communications, 2014, 185(1):86-95.
[2] HAMILTON B, FLETCHER L, MCCLEMENTS K, et al. Electron acceleration at reconnectiong X-points in solar flares[J]. The Astrophyscal Journal, 2005, 625(1):496-505.
[3] PFEFFERLE D, COPPER W A, GRAVES J P, et al. VENUS-LEVIS and its spline-Fourier interpolation of 3D toroidal magnetic field representation for guiding-centre and full-orbit simulations of charged energetic particles[J]. Computer Physics Communications, 2014, 185(12):3127-3140.
[4] DALENA S, CHUYCHAI P, MACE R L. Streamline generation code for particle dynamics description in numerical models of turbulence[J]. Computer Physics Communications, 2012, 183(9):1974-1985.
[5] TAUTZ R C. A new simulation code for particle diffusion in anisotropic, large-scale and turbulent magnetic fields[J]. Computer Physics Communications, 2010, 181(1):71-77.
[6] QIN H, GUAN X. Variational symplectic integrator for long-time simulations of the guiding-center motion of charged particles in general magnetic fields[J]. Physcal Review Letter, 2008, 100(3):035006.
[7] HE Y, SUN Y, LIU J, et al. Volume-preserving algorithms for charged particle dynamics[J]. Journal of Computational Physics, 2015, 281:135.
[8] HE Y, SUN Y, LIU J, et al. Higher order volume-preserving schemes for charged particle dynamics[J]. Journal of Computational Physics, 2016, 305:172.
[9] ZHANG R, LIU J, QIN H, et al. Application of lie algebra in constructing volume-preserving algorithms for charged particles dynamics[J]. Communications in Computational Physics, 2016, 19(5):1397.
[10] ZHANG H, WU S Z, ZHOU C T, et al. Numerical method of relativistic Fokker-Planck equation for energy deposition of fast electron[J]. Chinese Journal of Computational Physics. 2017,34(5):555-562.
[11] XU H, ZHUO H B, YANG X H, et al. Hybrid particle-in-cell/fluid model for hot electron transport in dense plasmas[J]. Chinese Journal of Computational Physics, 2017, 34(5):505-525.
[12] ZHENG P W, HE L H, HUANG Q H. Electron cyclotron wave heating and current drive simulation based on three-dimensional Fokker-Planck equation[J]. Chinese Journal of Computational Physics, 2017,34(4):387-393.
[13] KERBEL G, MCCOY M. Kinetic-theory and simulation of multispecies plasmas in Tokamaks excited with electromagnetic-waves in the ion-cyclotron range of frequencies[J]. Physics of Fluids, 1985, 28(12):3629-3653.
[14] EPPERIEIN E M, RICKARD G J, BELL A R, A code for the solution of the Vlasov-Fokker-Planck equation in 1-D or 2-D[J]. Computer Physics Communications, 1988, 52(1):7-13.
[15] MOREL J E, Fokker-Planck calculations using standard discrete ordinates transport codes[J]. Nuclear Science and Engineering, 1981, 79(4):340-356.
[16] CADJAN M G, IVANOV M F. Langevin approach to plasma kinetics with Coulomb collisions[J]. Journal of Plasma Physics, 1997, 61(1):89-106.
[17] SHERLOCK M. A Monte-Carlo method for Coulomb collisions in hybrid plasma models[J]. Journal of Computational Physics, 2008, 227:2286-2292.
[18] ROSIN M S, RICKETSON L F, DIMITS A M, et al. Multilevel Monte Carlo simulation of Coulomb collisions[J]. Journal of Computation Physics, 2014274:140-157.
[19] ZHANG Guannan, DIEGO del-Castillo-Negrete. A backward Monte-Carlo method for time-dependent runaway electron simulations[J]. Physics of Plasma, 2017, 24:092511.
[20] STREVENS D E. A new approach to charged particle slowing down and dispersion[J]. LLAL-TR-686862, 2016.
[21] RUDAKOV L I. Collective slowing down of an intense beam of relativistic electrons in a dense plasma target[J]. Soviet Physics JETP, 1971, 32(6):1134-1140.
[22] HUSSEINY A A, FORSEN H K. Scattering probability for fast test particles in a plasma[J]. Physics Review A, 1970, 2(5):2019-2024.
[23] DAR A, GRUNZWEIG-GENSSAR J, PERES A, et al. Slowing down of ions by ultrahigh-density electron plasma[J]. Physical Review Letters, 1974, 32(23):1299-1301.
[24] KLOEDEN P E, PLATEN E. Numerical solution of stochastic differential equations[M]//Applications of Mathematics. Berlin, New York:Springer-Verlag, 1992:101-103.
[25] ÖTTINGER H C. Stochastic processes in polymeric fluids[M]. Berlin:Springer, 1995:112.
[26] ROSENBLUTH M N, MACDONALD W M, JUDD D L. Fokker-Planck equation for an inverse-square force[J]. The Physical Review, 1957, 107(1).
[27] COHEN R S, SPITZER L, Jr ROUTLY P M. The electrical conductivity of an ionized gas[J]. The Physical Review, 1950, 80:230.
[28] BEN-ISRAEL A. A Newton-Raphson method for the solution of systems of equations[J]. Journal of Mathematical Analysis and Applications, 1966, 15(2):243-252.
[29] ORTEGA J M, RHEINBOLDT W C. Iterative solution of nonlinear equations in several variables[M]. New York and London:Academic Press, 1970:572.
[30] BJÖRCK A. Numerical methods for least squares problems[M]. North-Holland, Amsterdam:SIAM, Philadelphia, 1996, 343-345.
[31] Anonym.Lua.[EB/OL].http://www.lua.org/.
[32] Anonym.HDf5[EB/OL]. http://hdfgroup.org/HDF5/.
[33] MONTGOMERY D C, TIDMAN D A. Plasma kinetic theory[M]. New York, San Francisco, Toronto, London:McGRAW-HILL, 1964:26-37.
[34] GOEBEL D M, KATZ I. Fundamentals of electric propulsion:Ion and Hall thrusters[M]. Jet Propulsion Laboratory, California Institute of Technology, 2008:479-482.
[35] KRAUS M. Variational integrators in plasma physics[J]. arXiv preprint arXiv:1307.5665, 2013.
[36] XIAO J, LIU J, QIN H, et al. A variational multi-symplectic particle-in-cell algorithm with smoothing functions for the Vlasov-Maxwell system[J]. Physcal of Plasmas, 2013, 20(10):102517.
[37] QIN H, LIU J, XIAO J, et al. Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov-Maxwell equations[J]. Nuclear Fusion, 2015, 56(1):014001.
[38] XIAO J, QIN H, LIU J, et al. Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems[J]. Physcal of Plasmas, 2015, 22(11):112504.