Numkin 2019 Talk

Solving the Vlasov equation with Julia

The Vlasov equation is of fundamental importance in plasma physics and particularly for simulations in Magnetic fusion and of Tokamak plasmas. I present here a way to use the Julia language to solve it numerically.

After a short introduction about the language, the first example showed the kinetic simulation of Vlasov-Ampere system. The next example proposes an implementation of the Particle In Cell method.

We are basing much of this effort on a previous implementation in the Fortran language. We have found that the translation into Julia is easy and it is interesting to look at what Julia has to offer without degrading performance.

Either use the link above to open the notebooks in mybinder.org or run them locally:

git clone https://github.com/JuliaVlasov/Numkin2019
cd Numkin2019
julia --project
julia> using Pkg
julia> Pkg.instantiate()
julia> using IJulia
julia> notebook(dir=joinpath(pwd(),"notebooks"))
[ Info: running ...

Slides

  1. Introduction
  2. Vlasov-Ampere with FFT
  3. Rotation on GPU
  4. Metaprogrammimg and Particle Group
  5. GEMPIC
  6. Conclusion

links