Multi-scale simulation of plasma-object interaction

We present first simulation results from 3D3V hybrid kinetic method. The combination of discrete and continuous representation of distribution function allows controlling numerical noise and diffusion. Adaptive 3D grid automatically follows regions of high gradient and discontinuities such as shocks, radiation fronts, double-layers and sheaths. Our semi-analytical approach [1] to fractional sub-steps exceeds implicit schemes in accuracy and numerical stability. [1] O.Batishchev, Semi-Analytical Adaptive Vlasov -- Fokker-Planck -- Boltzmann Methods, pp.237-315, in book (Ed. M.Shoucri) Eulerian Codes for the Numerical Solution of the Kinetic Equations of Plasmas, Nova Science, 2010.