J M Donoso et al 2005 J. Phys. A: Math. Gen. 38 9145 doi:10.1088/0305-4470/38/41/021
J M Donoso1, J J Salgado2 and M Soler
Show affiliationsA recent method to obtain short-time propagators for finding path-integral solutions of Fokker–Planck equations is applied here to numerically solve the non-linear kinetic Fokker–Planck equation in plasma physics. Furthermore, we extend the use of this method to solve non-homogeneous equations. Cylindrical geometry in velocity space is used and two-species plasma is considered with no linearization of the exact conservative collisional operator. Numerical singularities in the diffusion tensor determinant are avoided by the splitting of the collisional operator into two parts, each one leading to different multiplicative integral operators which describe electron–electron and electron–ion interactions separately. The accurate advancing path-integral numerical formalism preserves conservative physical properties making this procedure a promising alternative to the classical linearized collisional operators used in kinetic theory. Here, we show the feasibility of the method by giving a new calculation of Spitzer's transport coefficients.
52.65.Ff Fokker-Planck and Vlasov equation
82C31 Stochastic methods (Fokker-Planck, Langevin, etc.) (See also 60H10)
Issue 41 (14 October 2005)
Received 5 June 2005, in final form 31 August 2005
Published 28 September 2005
J M Donoso et al 2005 J. Phys. A: Math. Gen. 38 9145
David S. Rupke et al. 2002 ApJ 570 588
D M Connor et al 2006 Phys. Med. Biol. 51 3283
Paola Porcari et al 2006 Phys. Med. Biol. 51 3141
César Esteban et al. 2009 ApJ 700 654
T P Flanagan 1957 J. Sci. Instrum. 34 450
S K Bhattacharya and S -I Chu 1985 J. Phys. B: At. Mol. Phys. 18 L275
Edson D Leonel and P V E McClintock 2006 J. Phys. A: Math. Gen. 39 11399
Pei-Qing Jin et al 2006 J. Phys. A: Math. Gen. 39 7115
O E Barndorff-Nielsen and R D Gill 2000 J. Phys. A: Math. Gen. 33 4481