@phdthesis{Ikkurthi2009,
  author    = {Venkata Ramana Ikkurthi},
  title     = {Kinetic Simulations of Dusty Plasmas},
  journal    = {Kinetische Simulationen Staubiger Plasmen},
  url       = {https://nbn-resolving.org/urn:nbn:de:gbv:9-000731-1},
  year      = {2009},
  abstract  = {This thesis constitutes a computational study of charge and ion drag force on micron-sized dust particles immersed in rf discharges. Knowledge of dust parameters like dust charge, floating potential, shielding and ion drag force is very crucial for explaining complex laboratory dusty plasma phenomena, such as void formation in microgravity experiments and wakefield formation in the sheaths. Existing theoretical models assume standard distribution functions for plasma species and are applicable over a limited range of flow velocities and collisionality. Kinetic simulations are suitable tools for studying dust charging and drag force computation. The main aim of this thesis is to perform three dimensional simulations using a Particle-Particle-Particle-Mesh (\$P^3M\$) model to understand how the dust parameters vary for different positions of dust in rf discharges and how these parameters on a dust evolve in the presence of neighboring dust particles. At first, rf discharges in argon have been modelled using a three-dimensional PIC-MCC code for the discharge conditions relevant to the dusty plasma experiments. All necessary elastic and inelastic collisions have been considered. The plasma background is found collisional, charge-exchange collisions between ions and neutrals being dominant. Electron and ion distributions are non-Maxwellian. The dominant heating mechanism is Ohmic. Then, simulations have been done to compute the dust parameters for various sizes of dust located at different positions in the rf discharges. Dust charge and floating potential in the presheath are slightly larger than the values in the bulk due to the higher electron flux to the dust particle in the presheath. From presheath to the sheath the charge and floating potential values decrease due to the decrease of the electron current to the dust. A linear dependence of dust potential on dust size has been found, which results in a nonlinear dependence of the dust charge with the dust size when the particle is assumed to be a spherical capacitor. This has been verified by independently counting the charges collected by the dust. \%where indeed it has been noted that the dust charge \%scales nonlinearly with the dust size. The computed dust parameters are also compared with theoretical models. Simulated dust floating potentials are comparable to values obtained from Allen-Boyd-Reynolds (ABR) and Khrapak models, but much smaller than the values obtained from Orbit Motion Limited (OML) model. The dust potential distribution behaves Debye-H\\"{u}ckel-like. The shielding lengths are in between ion and electron Debye lengths. \% indicating shielding by both ions and electrons. Further, the orbital drag force is typically larger than the collection drag force. The total drag force for the collisional case is larger than for the collisionless case and it scales nonlinearly with the dust size. The collection drag values and size-scaling agrees with Zobnin's model. The charging and drag force computation is then extended to two and multiple static dust particles in the rf discharge to study the influence of neighboring dust particles on the dust parameters. Initially, the dust parameters on two dust particles are computed for various interparticle separation distances and for dust particles placed at different locations in the rf discharge. It is observed that for dust separations larger than the shielding length the dust parameters for the two dust particles match with the single dust particle values. As the dust separation is equal to or less than the shielding length the ion drag force increases due to the buildup of a parallel drag force component. However, the main dust properties like charge, potential, vertical component of ion drag are not affected considerably. This is attributed to the smaller collection impact parameter values compared to the dust separation. \%This is because the \%collection impact parameter values in the sheath and the presheath are smaller \%than the smallest dust separation and in case of the dust in the bulk, the \%collection impact parameter is comparable with the dust separation. Then the dust charges on multiple dust particles located at different positions in the discharge and arranged along the discharge axis are also computed. It is found that the charges of the multiple dust particles in the bulk or presheath do not differ much from the single particle values at that location. But the dust charges of multiple dust particles located in the sheath drastically differ from the single dust parameter values. Due to ion focusing from dust particles in the upper layers, the ion current increases to dust particles in the lower layers resulting in smaller charge values. This is as well the case where dust particles are vertically aligned as in the standard experiments of dusty plasmas. In conclusion, this work used a fully kinetic (PIC and MD or \$P^3M\$) model to study the physics of dust charging in rf plasmas. Our simulations revealed that the dust parameters vary considerably from the bulk to the sheath. The CX collisions increase flux to the dust thereby affecting the dust parameters and their scaling with dust size. Also, a dust particle affects the charging dynamics of its neighbor only when their separation is within the shielding length. In the plasma sheath, ion focussing can cause great reduction in dust charges.},
  language  = {en}
}