Parsek - Parallel Code

Kinetic plasma simulation is the ultimate tool for plasma analysis. One of the prime tools for kinetic simulation is the particle in cell (PIC) method. The explicit or semi-implicit (i.e. implicit only on the fields) PIC method requires exceedingly small time steps and grid spacing, limited by the necessity to resolve the electron plasma frequency, the Debye length and the speed of light (for fully explicit schemes). A different approach is to consider fully implicit PIC methods where both particles and fields are discretized implicitly. This approach allows radically larger time steps and grid spacing, reducing the cost of a simulation by orders of magnitude while keeping the full kinetic treatment. In our previous work, simulations impossible for the explicit PIC method even on massively parallel computers have been made possible on a single processor machine using the implicit PIC code CELESTE3D. We propose here another quantum leap: PARSEK, a parallel cousin of CELESTE3D, based on the same approach but sporting a radically redesigned software architecture (object oriented, where CELESTE3D was structured and written in FORTRAN77/90) and fully parallelized using MPI for both particle and grid communication.

LAST PAPERS

Paolo Ricci, J.U. Brackbill, W. Daughton, Giovanni Lapenta A new role of the lower hybrid drift instability in magnetic reconnection, Phys. Plasmas, to appear

Giovanni Lapenta, D. Krauss-Varban, H. Karimabadi, J.D. Huba, L.I. Rudakov, Paolo Ricci Kinetic simulations of X-line propagation in 3D reconnection , Phys. Rev. Lett., submitted

Giovanni Lapenta, Philipp P. Kronberg Simulation of astrophysics jets: collimation and expansion into radio lobes, Astrophys. J, to appear

Giovanni Lapenta, D. Knoll Effect of a converging flow on the genesis of the slow solar wind, Astrophys. J, to appear

Enrico Camporeale, Giovanni Lapenta Model of bifurcated current sheets in the Earth's magnetotail: equilibrium and stability, J. Geophys. Res, to appear