2,599 research outputs found
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
- …