Calculation of selfconsistent electrostatic field in kinetic simulation of dusty plasma

The correct account of electrons and ions behavior near the dust grain is extremely important in kinetic simulation of dusty plasma. It is crucial both for distribution function calculation and for computation of the force experienced by the grain. The accuracy of description of particle motion depends on the precision of field calculation. Classical methods (MD, PIC, P$^3$M) of field calculation face difficulties if applied to dusty plasma. The modified PDPM method which is both accurate and economical is developed. The main idea is to use different spatial and time steps near the grain and far from it.