Dynamics controlled by magnetic fields:
parallel astrophysical computations
R. Keppens
FOM-Institute for Plasma Physics, The Netherlands
invited talk at Parallel CFD, 2000, Trondheim, Norway
I present the development history of and the obtained flexibility
within the Versatile Advection Code [VAC, see
http://www.phys.uu.nl/~toth], a software package for
solving sets of - near - conservation laws in any dimensionality.
Using a finite volume discretization on a static structured grid, the
versatility resides in the choice of applications (Euler, Navier-Stokes,
Magneto-Hydro-Dynamics), geometry, computer platform, shock-capturing
spatial and temporal discretizations. For distributed memory execution,
the VAC source code is automatically preprocessed to High Performance
Fortran which results in fully scalable applications.
In its most challenging configuration, the code is used for simulating
magnetized plasma dynamics in astrophysically relevant, three-dimensional
geometries. The magnetic field can severely complicate plasma flow behaviour
both in shock-dominated and otherwise (un)steady regimes. We present
examples of (i) `kinking' magnetic fluxtubes, related to emerging sunspots
on the solar surface; (ii) Kelvin-Helmholtz unstable magnetized jet flows,
present at all scales in our universe; and (iii) coronal mass
ejections in the continuously expanding solar magnetic environment.
All these simulations would benefit greatly from dynamically
controlled grid adaptivity, a research topic we are currently investigating.