Vlasov Simulation of Beams and HALO

This paper is devoted to the development and optimization of an adaptive grid Vlasov solver on a reduced axisymmetric model where the canonical angular momentum invariant is assumed to be 0. Whereas in PIC simulations most computational effort is put where the particle density is highest, the number of grid points in our adaptive Vlasov solver is higher in zones where the phase-space distribution function f varies most. Therefore, we expect that they might be more efficient for the investigation of non linear halo formation.