Dynamic Dust Particle Confinement in Corona Discharge Plasma

Usually electrostatic fields created by electrodes of the trap are used for dust particle confinement in plasma devices. In this work the possibility of the dust particle confinement by electrodynamic fields is investigated for improved quadrupole traps. The behavior of dust particles is simulated by Brownian dynamics. Dust particle parameters and parameters of the traps needed for dust particle confinement have been obtained (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Particle‐in‐Cell Simulations for Ion Thrusters

The Particle‐in‐Cell (PIC) method was used to study two different ion thruster concepts: Hall Effect Thrusters (HETs) and High Efficiency Multistage Plasma Thrusters (HEMPs), in particular the plasma properties in the discharge chamber due to the different magnetic field configurations. Special attention was paid to the simulation of plasma particles fluxes on the thrusters inner surfaces. In both cases PIC proved itself as a powerful tool, delivering important insight into the basic physics of the different thruster concepts.The simulations demonstrated that the new HEMP thruster concept allows for a high thermal efficiency due to both minimal energy dissipation and high acceleration efficiency. In the HEMP thruster the plasma contact to the wall is limited only to very small areas of the magnetic field cusps, which results in much smaller ion flux to the thruster channel surface as compared to HET. (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effects of the System Temperature on Dust Cluster Formation in Plasma

Effects of the system temperature on dust aggregation in plasmas are investigated using two‐dimensional molecular dynamics simulations. It is shown that as the system temperature increases, the boundary of the clusters becomes gradually irregular (i.e., deviating from sphere‐like), and the cluster system gradually changes from solid to liquid and finally to gas state. The mean square displacement, mean nearest‐neighbor distance in the clusters, cluster size and coupling parameter of the system are obtained and the properties of the system structure and dynamics are investigated. The time τ needed for reaching equilibrium for different temperatures is obtained. It is shown that τ firstly decreases and then increases with the temperature, indicating that there is an optimum temperature allowing a dust aggregation to reach an equilibrium state most rapidly. The simulation results agree qualitatively with the experimental observations. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dust in Plasma II. Effects of Secondary Electrons: Ionization and Surface Emission

In this second paper, the effect of secondary electrons on the charge and potential of a dust particle immersed in plasma has been studied. The processes of electron‐induced ionization and those of photo‐electron and secondary electron emission from the particle surface as a function of primary electron temperature have been taken into account. Starting from temperatures as low as 6 eV in an Ar plasma, ionization produces an extra ion flux to the dust surface comparable to that of the ion charge exchange effect. For what concerns the surface emission, results show that a transition from negative to positive dust charge/potential takes place, and that the transition regime is characterized by a non‐monotonic behavior of the electric potential around the particle. In the case of photoelectric emission, the dust charge and potential are monotonic decreasing functions of the electron temperature, while in the case of emission induced by primary electrons a minimum charge/potential is reached before they grow towards positive values. In no case multiple dust charge states have been observed due to the presence of the potential well attached to the particle surface. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Ion‐Acoustic Waves Instability in a Three Components Magneto‐Plasma with Nonthermal Electrons

A theoretical investigation has been made on obliquely propagating ion‐acoustic (IA) solitary structures in a three components magneto‐plasma containing cold inertial ions, Boltzmann distributed positrons, and hot non‐thermal electrons. The Zakharov‐Kuznetsov equation has been derived by the reductive perturbation method, and its solitary wave solution has been analyzed. Multi‐dimensional instability has also studied by the small‐k (long wave‐length plane wave) perturbation expansion technique, which is found to exist in such a plasma. The effects of the external magnetic field, nonthermal electrons, obliqueness and temperature ratio have significantly modified the basic properties of small but finite‐amplitude IA solitary waves, such as amplitude, width, instability criterion and the growth rate. The present investigation contributes to the physics of the nonlinear IA waves in space and laboratory electron‐positron‐ion magneto‐plasmas in which wave damping produces an electron tail. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Surface‐Driven Asymmetry and Instability in the Acceleration Region of Hall Thruster

It has been shown experimentally that the channel wall material has a substantial effect on the behaviour of Hall discharges. For this reason, the radial profile inside the Hall thruster SPT‐100 is investigated in detail. This is done by a one‐dimensional fully kinetic self‐consistent Particle‐in‐Cell model between the two walls in the acceleration region of the channel. A detailed Monte Carlo probabilistic model for secondary electron emission is implemented as boundary module. Using the local field approximation, two different operative conditions (axial electric field E_z =100 V/cm and 300 V/cm) have been simulated. For high discharge voltage case, a strong radial asymmetry and a stream instability propagating all along the radial domain are detected, while in the low voltage case a stable classical situation is recovered. The critical parameters for triggering this unstable regime are the electron azimuthal drift energy and the induced secondary electron emission, while the saturation mechanism is the increasing of the temperature of the initially cold secondary‐electrons. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Self‐Force in 1D Electrostatic Particle‐in‐Cell Codes for NonEquidistant Grids

Effects of non‐equidistant grids on momentum conservation is studied for simple test cases of an electrostatic 1D PIC code. The aim is to reduce the errors in energy and momentum conservation. Assuming an exact Poisson solver only numerical errors for the particle mover are analysed. For the standard electric field calculation using a central‐difference scheme, artificial electric fields at the particle position are generated in the case when the particle is situated next to a cell size change. This is sufficient to destroy momentum conservation. A modified electric field calculation scheme is derived to reduce this error. Independent of the calculation scheme additional fake forces in a two‐particle system are found which result in an error in the total kinetic energy of the system. This contribution is shown to be negligible for many particle systems. To test the accuracy of the two electric field calculation schemes numerical tests are done to compare with an equidistant grid set‐up. All tests show an improved momentum conservation and total kinetic energy for the modified calculation scheme of the electric field. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Plasma Polarization Around Dust Particle in an External Electric Field

The effect of plasma polarization around a negatively charged dust particle is investigated with the help of Monte Carlo simulation of ion trajectories in the electric field of the dust particle and an external electric field. The induced dipole moment of such a system was estimated in a wide range of dust particle and plasma parameters. It is shown that the dipole moment is very large, and has a non‐monotonous dependence on the external electric field. For a small external electric field, it weakly depends on the charge of the dust particle. The dipole moment reduces with the decrease of ion mean free path (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

PIC Simulation of Current‐Driven Buneman Instability in Presence of Collisional and Thermal Effects

Collisional and thermal effects on the nonlinear evolution of Buneman instability in an unmagnetized current‐driven plasma are investigated by particle in cell simulations. These simulations show that, as the time passes, the electron distribution function profiles deviate from initial shifted Maxwellian distribution and the electron phase‐space holes can be formed. The electron distribution function profiles also indicate the counter‐streaming and plateau formation. Moreover, the contour of the electric field and profiles of the growth rate display the resonance condition for this instability. Finally, the time evolution of the electric field energy density and ion kinetic energy in the presence of collisional and thermal effects are presented. These evolutions illustrate that in the presence of collisions, the growth rate of the Buneman instability is smaller than the collisionless case. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Particle in Cell/Monte Carlo Collision Method for Simulation of RF Glow Discharges: Effect of Super Particle Weighting

A parallel Particle in Cell/Monte Carlo Collision (PIC/MCC) numerical code for glow discharge plasma simulations is developed and verified. This method is based on simultaneous solution of the Lorentz equations of motion of super particles, coupled with the Poisson's equation for electric field. Collisions between the particles are modelled by the Monte Carlo method. Proper choice of particle weighting is critically important in order to perform adequate and efficient PIC simulations of plasma. Herein, effects of particle weighting on the simulations of capacitive radio‐frequency argon plasma discharges are studied in details. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dust Self‐Organized Structures I. Role of Ion Drag and Ion Diffusion on Screened Grains

Self‐organized dust structures are investigated using a Gurevich‐Parker model for non‐linear dust screening. The non‐linear dust drag coefficients and non‐linear diffusion coefficients are calculated numerically as functions of nonlinear parameter for screening, dust density and ion flux drift velocity. Nonlinear ion dust drag inside the structures creates an electric field with potential well for ions at the structure center. The equilibrium dust structures confine both the dust grains and the plasma particles, have a finite size and have inside an enhanced dust and ion densities. The necessary conditions for existence of equilibrium dust structures are found. The equilibrium dust structures are determined by two global parameters related to the external plasma flux and to the power of ionization. The equilibrium exist only in a restricted phase space of these two parameters and depends on the the drag coefficient at the structure center. The equilibrium requirements are found using non‐linear drag coefficient calculated numerically. It is shown that this phase space area can be broad but it is systematically decreasing with an increase of the ionization rate. It is found that equilibrium exists for dust structures with large dust and ion density concentration at the center and that for these structures the ion diffusion is strongly suppressed by ion scattering on non‐linearly screened grains. The results of the theory can be used to interpret the recently observed compact dust structures in micro‐gravity experiments and can provide some recommendations for future micro‐gravity experiments in spherical chambers (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Controlled Levitation of Dust Particles in RF+DC Discharges

Experiments and particle‐based kinetic simulations were performed to obtain the equilibrium levitation height of dust particles in plane parallel electrode discharges in low pressure argon gas, established by combined RF and DC excitation. The computed values were compared to experimental data. The good overall agreement of the simulation results and the experimental data verifies our gas discharge, dust charging, as well as dust force balance models. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Investigations of Droplet‐Plasma Interaction using Multi‐Dimensional Coupled Model

Recently developed multi‐dimensional coupled fluid‐droplet model is used to investigate the behavior of complex interaction between the liquid precursor droplets and atmospheric pressure plasma (APP). The significance of this droplet‐plasma interaction is not well understood under diverse realm of working conditions in two‐phase flow. In this study, we explain the implication of vaporization of liquid droplets in APP which are subsequently responsible to control major characteristics of surface coating depositions. Coalescence of water droplets is more dominant than Hexamethyldisiloxane (HMDSO) droplets because of its sluggish rate of evaporation. A disparity in the performance of evaporation is identified in two independent mediums, such as gas mixture and discharge plasma using HMDSO precursor. The length of evaporation of droplets is amplified by an increment of gas flow rate indicating with a reduction in the gas temperature and electron mean energy. In particular, the spatio‐temporal density distributions of charged particles show a clear pattern in which the typical nitrogen impurity ions are primarily effective as compared to other helium ionic species along the pulse of droplets in APP. Finally, we contrast the behavior of discharge species in the pure helium and He‐N2 gas mixtures revealing the importance of stepwise and Penning ionization processes. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On‐line Determination of Particle Size Distributions in Flowing Dispersions

Emulsions are of great importance to industry. They are involved in many engineering operations, including chemical reactions, extraction, emulsification and suspension polymerization, etc. However, an important problem for these processes is how to control the size distribution of the dispersed phase. Indeed, off‐line analysis of the emulsion may generate uncertainties due to sampling and dilution of the product, which are likely to change the dispersion state and physico‐chemical properties. In this work, an on‐line optical method is proposed to characterize dispersed media in real flowing conditions. This method is based on the time‐analysis of back‐scattered light fluctuations. The present paper deals with the development of this method and its application to dispersions of alumina in water. The results obtained with the on‐line optical method are compared with those acquired by classical laser light scattering and microscopy.     

Elementary Many‐Particle Processes in Plasma Microfields

The effect of electric and magnetic plasma microfields on elementary many‐body processes in plasmas is considered. As detected first by Inglis and Teller in 1939, the electric microfield controls several elementary processes in plasmas as transitions, line shifts and line broadening. We concentrate here on the many‐particle processes ionization, recombination, and fusion and study a wide area of plasma parameters. In the first part the state of art of investigations on microfield distributions is reviewed in brief. In the second part, various types of ionization processes are discussed with respect to the influence of electric microfields. It is demonstrated that the processes of tunnel and rescattering ionization by laser fields as well as the process of electron collisional ionization may be strongly influenced by the electric microfields in the plasma. The third part is devoted to processes of microfield action on fusion processes and the effects on three‐body recombination are investigated. It is shown that there are regions of plasma densities and temperatures, where the rate of nuclear fusion is accelerated by the electric microfields. This effect may be relevant for nuclear processes in stars. Further, fusion processes in ion clusters are studied. Finally we study in this section three‐body recombination effects and show that an electric microfield influences the three‐body electron‐ion recombination via the highly excited states. In the fourth part, the distribution of the magnetic microfield is investigated for equilibrium, nonequilibrium, and non‐uniform magnetized plasmas. We show that the field distribution in a neutral point of a non‐relativistic ideal equilibrium plasma is similar to the Holtsmark distribution for the electrical microfield. Relaxation processes in nonequilibrium plasmas may lead to additional microfields. We show that in turbulent plasmas the broadening of radiative electron transitions in atoms and ions, without change of the principle quantum number, may be due to the Zeeman effect and may exceed Doppler and Stark broadening as well. Further it is shown that for optical radiation the effect of depolarization of a linearly polarized laser beams propagating through a magnetized plasma may be rather strong. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Particle‐in‐cell simulation of an optimized high‐efficiency multistage plasma thruster

Electric propulsion attracts increasing attention in contemporary space missions as an interesting alternative to chemical propulsion because of the high efficiency it offers. The High‐Efficiency Multistage Plasma thruster, a class of cusped field thruster, is able to operate at different anode voltages and operation points and thereby generate different levels of thrust in a stable and efficient way. Since experiments of such thrusters are inherently expensive, multi‐objective design optimization (MDO) is of great interest. Several optimized thruster designs have resulted from a MDO model based on a zero‐dimensional (0D) power balance model. However, the MDO solutions do not warrant self‐consistency due to their dependency on estimation from empirical modelling based on former experimental studies. In this study, one of the optimized thruster designs is investigated by means of particle‐in‐cell (PIC) analysis to examine the predicted performance characteristics with self‐consistent simulations. The 0D power balance model is used to develop additional diagnostics for the PIC simulations to improve the physics analysis. Using input parameters for the 0D power balance model from the PIC simulations allows further improvement for the design optimization.     

THz Radiation Generation via the Interaction of Ultra‐short Laser Pulses with the Molecular Hydrogen Plasma

In this paper, a two dimensional Particle In Cell‐Monte Carlo Collision simulation scheme is used to examine the THz generation via the interaction of high intensity ultra‐short laser pulses with an underdense molecular hydrogen plasma slab. The influences of plasma density, laser pulse duration and its intensity on the induced plasma current density and the subsequent effects on the generated THz signal characteristics are studied. It is observed that the induced current density in the plasma medium and THz spectral intensity are increased at the higher laser pulse intensities, laser pulse durations and plasma densities. Moreover, the generated THz electric field amplitude is reduced at the higher laser pulse durations. A wider frequency range for the generated THz signal is shown at the lower laser pulse durations and higher plasma densities. Additionally, it is found that the induced current density in hydrogen plasma medium is the dominant factor influencing the generation of THz pulse radiation. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Filamentation Instability in a Current‐Carrying Plasma in the Presence of Quantum Effects

The filamentation instability of a current‐carrying plasma under the diffusion condition is investigated taking into account the Bohm potential and the Fermi electron pressure. Using quantum hydrodynamic equations, the dispersion relation and growth rate of the instability is obtained. It is found that the filamentation instability, in the presence of quantum effects, depends on various characteristic parameters such as: electron Fermi velocity, plasma number density, ion thermal velocity and electron drift velocity. Moreover, the wavelength region in which the instability occurs is more restricted and the minimum size of filaments is larger, in comparison with the classical case. It is also found that the growth rate of the instability is smaller in the presence of quantum effects. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)