Solving the Vlasov equation with Julia
Numerical Methods for the Kinetic Equations of Plasma Physics, Garching Oct 14-18, 2019.
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.
Slides
- Introduction
- Vlasov-Ampere with FFT
- Rotation on GPU
- Metaprogrammimg and Particle Group
- GEMPIC
- Conclusion
Jupyter notebooks
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 ...