Many plasma physics problems in astrophysics, laboratory and space systems require the development of suitable relativistic plasma simulation methods. Collisionless particle-in-cell (PIC) methods are often used to this purpose. We investigate here the possibility of a relativistic extension of an existing implicit PIC code (CELESTE), based on a moment equation closure model (1). Four possible ways of formulating the relativistic moment expansion are analysed from a purely numerical point of view, studying truncation errors and stability properties. Based on the analysis a new method is implemented into the 1D pilot code, PARSEK (2). The classical equation of motion from CELESTE is modified in order to keep a similar implicit integration method. The field solver is also presented, taking into account the previous results of the classical version. \\ \\ (1) J.U. Brackbill, D.W. Forslund, J. Computat. Physics, 46, 271, 1982. (2) S. Markidis, G. Lapenta, W.B. VanderHeyden, Z. Budimlic, Concurrency Comput Practice Experience, 17, 821, 2005.
Tronci, C., Zuccaro, G., Relativistic Implicit Moment Method, Abstract de <<APS Division of Plasma Physics Meeting Abstracts, 2005>>, (Denver (Colorado - U.S.A.), 24-28 October 2005 ), The American Physical Society, Denver (Colorado - U.S.A.) 2005:47 69-69 [https://hdl.handle.net/10807/226827]
Relativistic Implicit Moment Method
Zuccaro, Gianluca
Secondo
2005
Abstract
Many plasma physics problems in astrophysics, laboratory and space systems require the development of suitable relativistic plasma simulation methods. Collisionless particle-in-cell (PIC) methods are often used to this purpose. We investigate here the possibility of a relativistic extension of an existing implicit PIC code (CELESTE), based on a moment equation closure model (1). Four possible ways of formulating the relativistic moment expansion are analysed from a purely numerical point of view, studying truncation errors and stability properties. Based on the analysis a new method is implemented into the 1D pilot code, PARSEK (2). The classical equation of motion from CELESTE is modified in order to keep a similar implicit integration method. The field solver is also presented, taking into account the previous results of the classical version. \\ \\ (1) J.U. Brackbill, D.W. Forslund, J. Computat. Physics, 46, 271, 1982. (2) S. Markidis, G. Lapenta, W.B. VanderHeyden, Z. Budimlic, Concurrency Comput Practice Experience, 17, 821, 2005.I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.