DISPATCH: A Numerical Simulation Framework for the Exa-scale Era. I. Fundamentals

We introduce a high-performance simulation framework that permits the semi-independent, task-based solution of sets of partial differential equations, typically manifesting as updates to a collection of `patches' in space-time. A hybrid MPI/OpenMP execution model is adopted, where work tasks are controlled by a rank-local `dispatcher' which selects, from a set of tasks generally much larger than the number of physical cores (or hardware threads), tasks that are ready for updating. The definition of a task can vary, for example, with some solving the equations of ideal magnetohydrodynamics (MHD), others non-ideal MHD, radiative transfer, or particle motion, and yet others applying particle-in-cell (PIC) methods. Tasks do not have to be grid-based, while tasks that are, may use either Cartesian or orthogonal curvilinear meshes. Patches may be stationary or moving. Mesh refinement can be static or dynamic. A feature of decisive importance for the overall performance of the framework is that time steps are determined and applied locally; this allows potentially large reductions in the total number of updates required in cases when the signal speed varies greatly across the computational domain, and therefore a corresponding reduction in computing time. Another feature is a load balancing algorithm that operates `locally' and aims to simultaneously minimise load and communication imbalance. The framework generally relies on already existing solvers, whose performance is augmented when run under the framework, due to more efficient cache usage, vectorisation, local time-stepping, plus near-linear and, in principle, unlimited OpenMP and MPI scaling.Comment: 17 pages, 8 figures. Accepted by MNRA

In-situ Particle Acceleration in Collisionless Shocks

The outflows from gamma ray bursts, active galactic nuclei and relativistic jets in general interact with the surrounding media through collisionless shocks. With three dimensional relativistic particle-in-cell simulations we investigate such shocks. The results from these experiments show that small--scale magnetic filaments with strengths of up to percents of equipartition are generated and that electrons are accelerated to power law distributions N(E)~E^{-p} in the vicinity of the filaments through a new acceleration mechanism. The acceleration is locally confined, instantaneous and differs from recursive acceleration processes such as Fermi acceleration. We find that the proposed acceleration mechanism competes with thermalization and becomes important at high Lorentz factors.Comment: 4 pages, 2 figures, submitted to Il nuovo cimento (4th Workshop Gamma-Ray Bursts in the Afterglow Era, Rome, 18-22 October 2004

Non-Fermi Power law Acceleration in Astrophysical Plasma Shocks

Collisionless plasma shock theory, which applies for example to the afterglow of gamma ray bursts, still contains key issues that are poorly understood. In this paper we study charged particle dynamics in a highly relativistic collisionless shock numerically using ~10^9 particles. We find a power law distribution of accelerated electrons, which upon detailed investigation turns out to originate from an acceleration mechanism that is decidedly different from Fermi acceleration. Electrons are accelerated by strong filamentation instabilities in the shocked interpenetrating plasmas and coincide spatially with the power law distributed current filamentary structures. These structures are an inevitable consequence of the now well established Weibel-like two-stream instability that operates in relativistic collisionless shocks. The electrons are accelerated and decelerated instantaneously and locally; a scenery that differs qualitatively from recursive acceleration mechanisms such as Fermi acceleration. The slopes of the electron distribution power laws are in concordance with the particle power law spectra inferred from observed afterglow synchrotron radiation in gamma ray bursts, and the mechanism can possibly explain more generally the origin of non-thermal radiation from shocked inter- and circum-stellar regions and from relativistic jets.Comment: 4 pages accepted for publication in ApJ Letters. High resolution figures are available online at http://www.astro.ku.dk/users/hededal/040855