A stochastic differential equation code for multidimensional Fokker-Planck type problems

Abstract We present a newly developed numerical code that integrates Fokker–Planck type transport equations in four to six spatial dimensions (configuration plus momentum space) and time by means of stochastic differential equations. In contrast to other, similar approaches our code is not restricted to any special configuration or application, but is designed very generally with a modular structure and, moreover, allows for Cartesian, cylindrical or spherical coordinates. Depending on the physical application the code can integrate the equations forward or backward in time. We exemplify the mathematical ideas the method is based upon and describe the numerical realisation and implementation in detail. The code is validated for both cases against an established finite-differences explicit numerical code for a scenario that includes particle sources as well as a linear loss term. Finally we discuss the new possibilities opened up with respect to general applications and newly developed hardware.

[1]  M. Pohl,et al.  Cosmic-Ray Propagation Properties for an Origin in Supernova Remnants , 2004, astro-ph/0410352.

[2]  D. Thomson.  Criteria for the selection of stochastic models of particle trajectories in turbulent flows , 1987, Journal of Fluid Mechanics.

[3]  A. Gourlay,et al.  On the Structure of Alternating Direction Implicit (A.D.I.) and Locally One Dimensional (L.O.D.) Difference Methods , 1972 .

[4]  Ming Zhang,et al.  PROPAGATION OF SOLAR ENERGETIC PARTICLES IN THREE-DIMENSIONAL INTERPLANETARY MAGNETIC FIELDS: IN VIEW OF CHARACTERISTICS OF SOURCES , 2009 .

[5]  Takuji Nishimura,et al.  Mersenne twister: a 623-dimensionally equidistributed uniform pseudo-random number generator , 1998, TOMC.

[6]  J. Birn,et al.  Quantitative study of the nonlinear formation and acceleration of plasmoids in the Earth's magnetotail , 1990 .

[7]  M. Kac.  On distributions of certain Wiener functionals , 1949 .

[8]  M. Zhang,et al.  ENERGY SPECTRUM OF ENERGETIC PARTICLES ACCELERATED BY SHOCK WAVES: FROM FOCUSED TRANSPORT TO DIFFUSIVE ACCELERATION , 2011 .

[9]  Eugene N. Parker,et al.  THE PASSAGE OF ENERGETIC CHARGED PARTICLES THROUGH INTERPLANETARY SPACE , 1965 .

[10]  N. Pogorelov,et al.  FOUR-DIMENSIONAL TRANSPORT OF GALACTIC COSMIC RAYS IN THE OUTER HELIOSPHERE AND HELIOSHEATH , 2009 .

[11]  B. Heber,et al.  On the definition and calculation of a generalised McIlwain parameter , 2010 .

[12]  David Potter.  Computational physics , 1973 .

[13]  J. Clem,et al.  A general time-dependent stochastic method for solving Parker's transport equation in spherical coordinates , 2010 .

[14]  C. Gardiner.  Handbook of Stochastic Methods , 1983 .

[15]  John Schoenmakers,et al.  Transition density estimation for stochastic differential equations via forward-reverse representations , 2004 .

[16]  Reinhard Schlickeiser,et al.  Cosmic Ray Astrophysics , 2002 .

[17]  Ming Zhang.  A Markov Stochastic Process Theory of Cosmic-Ray Modulation , 1999 .

[18]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .

[19]  R. D. Strauss,et al.  MODELING THE MODULATION OF GALACTIC AND JOVIAN ELECTRONS BY STOCHASTIC PROCESSES , 2011 .

[20]  M. Potgieter.  The Modulation of Galactic Cosmic Rays in the Heliosphere: Theory and Models , 1998 .

[21]  A. Etheridge.  A course in financial calculus , 2002 .

[22]  M. E. Muller,et al.  A Note on the Generation of Random Normal Deviates , 1958 .

[23]  T. Krüger,et al.  A Fisk-Parker Hybrid Heliospheric Magnetic Field with a Solar-Cycle Dependence , 2008 .

[24]  E. Renshaw,et al.  STOCHASTIC DIFFERENTIAL EQUATIONS , 1974 .