VlaPlEx (Vlasov Plasma Expansion) code aims to provide a simple way to solve the collisionless plasma expansion into a vacuum.

It solves the Vlasov equation for ions in a 1D spherical domain.

Electrons are treated as a fluid either isothermal (Boltzmann) or polytropic.

Self-consistent electric field is solved using a Newton-Raphson method to solve the Poisson equation.

A 2D table linking T_e, n_e and the average Z must be provided. Examples can be found in 'data/TNZ'.

Input files for different cases are in the 'input' folder.

Files with the boundary conditions are provided in 'data/boundary'.

Input files must specify which TNZ table and boundary file are using.

To run the code execute

./vlaplex <path-to-input-file>