Achieving large-scale kinetic modeling is a crucial task for the development and optimization of modern plasma devices. With the trend of decreasing pressure in applications, such as plasma etching, kinetic simulations are necessary to self-consistently capture the particle dynamics. The standard, explicit, electrostatic, momentum-conserving particle-in-cell method suffers from restrictive stability constraints on spatial cell size and temporal time step, requiring resolution of the electron Debye length and electron plasma period, respectively. This results in a very high computational cost, making the technique prohibitive for large volume device modeling. We investigate the direct implicit algorithm and the explicit energy conserving algorithm as alternatives to the standard approach, both of which can reduce computational cost with a minimal (or controllable) impact on results. These algorithms are implemented into the well-tested EDIPIC-2D and LTP-PIC codes, and their performance is evaluated via 2D capacitively coupled plasma discharge simulations. The investigation reveals that both approaches enable the utilization of cell sizes larger than the Debye length, resulting in a reduced runtime, while incurring only minor inaccuracies in plasma parameters. The direct implicit method also allows for time steps larger than the electron plasma period; however, care must be taken to avoid numerical heating or cooling. It is demonstrated that by appropriately adjusting the ratio of cell size to time step, it is possible to mitigate this effect to an acceptable level.

We study collective processes for an electron beam propagating through a background plasma using simulations and analytical theory. A new regime where the instability of a Langmuir wave packet can grow locally much faster than ion frequency is clearly identified. The key feature of this new regime is an electron modulational instability that rapidly creates a local Langmuir wave packet, which in its turn produces local charge separation and strong ion density perturbations because of the action of the ponderomotive force, such that the beam-plasma wave interaction stops being resonant. Three evolution stages of the process and observed periodic burst features are discussed. Different physical regimes in the plasma and beam parameter space are demonstrated for the first time.

Electron-beam plasma interaction has long been a topic of great interest. Despite the success of the quasilinear and weak turbulence theories, their validities are limited by the requirements of a sufficiently dense mode spectrum and a small wave amplitude. In this paper, we extensively study the collective processes of a mono-energetic electron beam emitted from a thermionic cathode propagating through a cold plasma by performing high-resolution two-dimensional particle-in-cell simulations and using analytical theories. We confirm that, during the initial stage of two-stream instability between the beam and background cold electrons, it is saturated due to the well-known wave-trapping mechanism. Further evolution occurs due to strong wave-wave nonlinear processes. We show that the beam-plasma interaction can be classified into four different physical regimes in the parameter space for the plasma and beam parameters. The differences between the regimes are analyzed in detail. We identify a new regime in the strong Langmuir turbulence featured by what we call electron modulational instability (EMI) that could create a local Langmuir wave packet growing faster than the ion plasma frequency. Ions do not have time to respond to EMI in the initial growing stage. On a longer timescale, the action of the ponderomotive force produces very strong ion density perturbations, and eventually, the beam-plasma wave interaction stops being resonant due to the strong ion density perturbations. Consequently, in this EMI regime, electron beam-plasma interaction occurs in a repetitive (intermittent) process. The beam electrons are strongly scattered by waves, and the Langmuir wave spectrum is significantly broadened, which in turn gives rise to strong heating of bulk electrons. Associated energy transfer from the beam to the background plasma electrons has been studied. A resulting kappa (κ) distribution and a wave-energy spectrum $E^2\left(k\right)\sim k^{-5}$ are observed in the strong turbulent regime.

In order to study the properties of short carbon arcs, a self-consistent model was implemented into a CFD code ANSYS-CFX. The model treats the transport of heat and electric current in the plasma and electrodes in a coupled manner and accounts for gas convection in the chamber. Multiple surface processes at the electrodes are modeled, including the formation of space-charge limited sheaths, ablation and deposition of carbon, and emission and absorption of radiation and electrons. The simulations show that the arc is constricted near the cathode and anode front surfaces, leading to the formation of electrode spots. The cathode spot is a well-known phenomenon, and mechanisms of its formation were reported elsewhere. However, the anode spot formation mechanism discovered in this work was not previously reported. We conclude that the spot formation is not related to plasma instability, as commonly believed in the case of constricted discharge columns, but rather occurs due to the highly nonlinear nature of heat balance in the anode. We additionally demonstrate this property with a reduced anode heat transfer model. We also show that the spot size increases with the arc current. This anode spot behavior was also confirmed in our experiments. Due to the anode spot formation, a large gradient of carbon gas density occurs near the anode, which drives a portion of the ablated carbon back to the anode at its periphery. This can consequently reduce the total ablation rate. Simulation results also show that the arc can reach the local chemical equilibrium state in the column region, while the local thermal equilibrium state is not typically achieved for experimental conditions. It shows that it is important to account for different electron and gas temperatures in the modeling of short carbon arcs.

Direct-current ionization breakdown in argon at very high voltage on the low-pressure branch of the Paschen curve is studied by means of kinetic simulation. The employed particle-in-cell/Monte Carlo collision model tracks the motion of electrons, ions and fast neutral atoms, as well as their collisions with the background gas. Secondary electron emission due to the bombardment by ions and fast atoms, as well as backscattering of different species on the electrode surfaces are also consistently accounted for. To validate the model, we compare the Paschen curve predicted by kinetic simulations against published experimental results as well against the prediction of a reduced (fluid) model. The model is then applied to predict, for the first time, the low-pressure branch of the Paschen curve for argon gas in the 10–300 kiloVolt range of applied voltage. One of the findings is the existence of a turning point on the Paschen curve at V b r ≈ 20 k V and p d ≈ 0.08 T o r r ⋅ c m , where V b r is the breakdown voltage and the reduced pressure p d is the product of the gas pressure with the inter-electrode distance. According to our analysis, the presence of the turning point is due to the phenomena of electron and ion runaway in high electric field. We further investigate the role of the heavy species (i.e. ions and fast neutral atoms) in sustaining the Townsend discharge at high voltage and thus setting the breakdown threshold. Specifically, the role of fast atoms is significant over the entire low-pressure branch of the Paschen curve due to their contribution to both ionization and secondary electron emission. Also, the ions, as the precursor for the fast atoms, show an important influence on the latter process. The newly predicted low-pressure branch of the Paschen curve for argon at high voltage should provide an essential reference for the design of high-voltage devices, e.g. circuit breakers in electrical power transmission.

Azimuthal structures emerging in beam-generated partially magnetized plasmas are investigated using three-dimensional particle in cell-Monte Carlo collision simulations. Two distinct instability regimes are identified at low pressures. When the gas pressure is sufficiently high, quasineutrality is attained and 2D spiral-arm structures form as a result of the development of a lower-hybrid instability, resulting in enhanced cross-field transport. At lower pressures, quasineutrality is not achieved and a 3D helical-rotating plasma structure forms due to development of the diocotron instability. Analytical formulas are proposed for the critical threshold pressure between these regimes and for the rotation frequency of the helical structures. Preliminary experimental verification is provided.