Space weather can be caused by solar high-energy particles or by dynamic variations of the solar wind that can cause extended periods of major disturbances on ground and space, affecting technological systems (e.g., telecommunication and weather spacecraft at geostationary orbit, and ground-based power grids).

In Vlasiator, ions are represented as velocity distribution functions, while electrons are magnetohydrodynamic fluid, enabling a self-consistent global plasma simulation that can describe multi-temperature plasmas to resolve non-MHD processes that currently cannot be self-consistently described by the existing global space weather simulations. The novelty is that by modelling ions as velocity distribution functions the outcome will be numerically noiseless.

Due to the multi-dimensional approach at ion scales, Vlasiator's computational challenges are immense. We use advanced high performance computing techniques to allow massively parallel computations on tens of thousands of cores.

More details about the implementation can be found below, and in the presentations and publications, especially this paper.

Vlasiator simulates the dynamics of plasma using a hybrid-Vlasov model, where protons are described by their distribution function f(r,v,t) in ordinary (r) and velocity (v) space, and electrons are a charge-neutralising fluid. This approach neglects electron kinetic effects but retains ion kinetics. The time-evolution of f(r,v,t) is given by Vlasov's equation,

which is coupled self-consistently to Maxwell's equations giving the evolution of the electric and magnetic fields E and B. Maxwell's equations neglect the displacement current. The equations are closed by a generalised Ohm's law including the Hall term.

Vlasiator propagates the distribution function forward in time with a conservative fifth-order accurate Semi-Lagrangian algorithm. This algorithm allows using long time steps even in the presence of strong magnetic fields, as the propagation in velocity space is not limited by the Courant-Friedrichs-Levy (CFL) condition. The field solver is a second-order accurate divergence-free upwind-constrained transport method.

Vlasiator has a parallel Cartesian mesh in ordinary space. In each spatial cell there is a 3-dimensional sparse velocity grid, modelling the in the full 6-dimensional distribution function. Empty velocity space cells are neither stored nor propagated, which in a typical case reduces the total number of phase space cells by a factor of at least 100 while impacting mass conservation to a relative level of no more than 10−6.

The ordinary space grid is implemented using the open source DCCRG grid library developed by the group. It is parallelized using MPI-based domain decomposition and OpenMP-based threading is used to further parallelise the work done by each process. The Vlasov solver is vectorised using AVX intrinsics. The load is balanced with the Zoltan library using its recursive coordinate bisection partitioner. I/O is performed using our own parallel VLSV file format, which can be analyzed using VisIt or by using the python based Analysator package.