Electromagnetic gyrokinetic PIC simulation with an adjustable control variates method

In the last decade, it became clear that electromagnetic (gyro)kinetic particle-in-cell (PIC) simulations are very demanding in respect to numerical methods and the number of markers used. The Monte Carlo discretization of the gyrokinetic equations leads to a severe signal-to-noise problem: the statistical representation of the physically irrelevant but numerically dominant adiabatic current causing a high statistical noise level. The corresponding inaccuracy problem is very pronounced at high plasma @b and/or small perpendicular wave numbers k"@?. We derive several numerical schemes to overcome the problem using an adjustable control variates method which adapts to the dominant adiabatic part of the gyro-center distribution function. We have found that the inaccuracy problem is also present in the quasi-neutrality equation as a consequence of the p"@?-formulation [T.S. Hahm, W.W. Lee, A. Brizard, Nonlinear gyrokinetic theory for finite-beta plasmas, Phys. Fluids 31 (1988) 1940]. For slab simulations in the magnetohydrodynamic (MHD) limit k"@?->0, the number of markers can be reduced by more than four orders of magnitude compared to a conventional @df scheme. The derived schemes represent first steps on a road to fully adaptive control variates method which can significantly reduce the inherent statistical noise of PIC codes.

[1]  A. Aydemir.  A unified Monte Carlo interpretation of particle simulations and applications to non-neutral plasmas , 1994 .

[2]  R. Hatzky,et al.  Energy conservation in a nonlinear gyrokinetic particle-in-cell code for ion-temperature-gradient-driven modes in θ-pinch geometry , 2002 .

[3]  R. Hatzky,et al.  A revised δf algorithm for nonlinear PIC simulation , 2003 .

[4]  S. Parker,et al.  A fully nonlinear characteristic method for gyrokinetic simulation , 1993 .

[5]  Roman Hatzky,et al.  Conventional δf-particle simulations of electromagnetic perturbations with finite elements , 2004 .

[6]  Laurent Villard,et al.  Finite element approach to global gyrokinetic Particle-In-Cell simulations using magnetic coordinates , 1998 .

[7]  C. R. Deboor,et al.  A practical guide to splines , 1978 .

[8]  Roman Hatzky,et al.  Particle simulations with a generalized gyrokinetic solver , 2005 .

[9]  W. W. Lee,et al.  The split-weight particle simulation scheme for plasmas , 2000 .

[10]  A. Friedman.  Techniques for Robust Nonlinear Delta-F Simulations of Beams , 2001 .

[11]  W. Lee,et al.  Gyrokinetic Particle Simulation Model , 1987 .

[12]  Scott E. Parker,et al.  A δ f particle method for gyrokinetic simulations with kinetic electrons and electromagnetic perturbations , 2003 .

[13]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[14]  Alain J. Brizard,et al.  Nonlinear gyrokinetic theory for finite‐beta plasmas , 1988 .

[15]  John Van Wicheren Reynders.  Gyrokinetic simulation of finite-beta plasmas on parallel architecture , 1993 .

[16]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[17]  K. Höllig.  Finite element methods with B-splines , 1987 .

[18]  Laurent Villard,et al.  A drift-kinetic Semi-Lagrangian 4D code for ion turbulence simulation , 2006, J. Comput. Phys..

[19]  W. Lee,et al.  Partially linearized algorithms in gyrokinetic particle simulation , 1993 .

[20]  J. L. V. Lewandowski,et al.  Shear-Alfven Waves in Gyrokinetic Plasmas , 2000 .

[21]  C. D. Beidler,et al.  Physics and Engineering Studies for Wendelstein 7-X , 1991 .

[22]  Alain J. Brizard,et al.  Nonlinear gyrokinetic Maxwell-Vlasov equations using magnetic co-ordinates , 1989, Journal of Plasma Physics.

[23]  Roman Hatzky,et al.  Gyrokinetic simulations with a particle discretization of the field equations , 2005 .

[24]  W. A. Cooper,et al.  Global linear gyrokinetic simulations in quasi-symmetric configurations , 2001 .

[25]  J. Craggs.  Applied Mathematical Sciences , 1973 .

[26]  C. Birdsall,et al.  Plasma Physics via Computer Simulation , 2018 .