Many high power electronic devices operate in a regime where the current they draw is limited by the self-fields of the particles. This space-charge-limited current poses particular challenges for numerical modeling where common techniques like over-emission or Gauss Law are computationally inefficient or produce nonphysical effects. In this paper we show an algorithm using the value of the electric field in front of the surface instead of attempting to zero the field at the surface, making the algorithm particularly well suited to both electromagnetic and parallel implementations of the PIC algorithm. We show how the algorithm is self-consistent within the framework of finite difference (for both electrostatics and electromagnetics). We show several 1D and 2D benchmarks against both theory and previous computational results. Finally we show application in 3D to high power microwave generation in a 13 GHz magnetically insulated line oscillator.
This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. These results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.
A massively parallel simulation code, called textit{dHybrid}, has been developed to perform global scale studies of space plasma interactions. This code is based on an explicit hybrid model; the numerical stability and parallel scalability of the code are studied. A stabilization method for the explicit algorithm, for regions of near zero density, is proposed. Three-dimensional hybrid simulations of the interaction of the solar wind with unmagnetized artificial objects are presented, with a focus on the expansion of a plasma cloud into the solar wind, which creates a diamagnetic cavity and drives the Interplanetary Magnetic Field out of the expansion region. The dynamics of this system can provide insights into other similar scenarios, such as the interaction of the solar wind with unmagnetized planets.
Upon inclusion of collisions, the speed-limited particle-in-cell (SLPIC) simulation method successfully computed the Paschen curve for argon. The simulations modelled an electron cascade across an argon-filled capacitor, including electron-neutral ionization, electron-neutral elastic collisions, electron-neutral excitation, and ion-induced secondary-electron emission. In electrical breakdown, the timescale difference between ion and electron motion makes traditional particle-in-cell (PIC) methods computationally slow. To decrease this timescale difference and speed up computation, we used SLPIC, a time-domain algorithm that limits the speed of the fastest electrons in the simulation. The SLPIC algorithm facilitates a straightforward, fully-kinetic treatment of dynamics, secondary emission, and collisions. SLPIC was as accurate as PIC, but ran up to 200 times faster. SLPIC accurately computed the Paschen curve for argon over three orders of magnitude in pressure.
Experimental observations have long-established that there exists a smooth roll-off or knee transition between the temperature-limited (TL) and full-space-charge-limited (FSCL) emission regions of the emission current density-temperature J-T (Miram) curve, or the emission current density-voltage J-V curve for a thermionic emission cathode. In this paper, we demonstrate that this experimentally observed smooth transition does not require frequently used a priori assumptions of a continuous distribution of work functions on the cathode surface. Instead, we find the smooth transition arises as a natural consequence of the physics of nonuniform thermionic emission from a spatially heterogeneous cathode surface. We obtain this smooth transition for both J-T and J-V curves using a predictive nonuniform thermionic emission model that includes 3-D space charge, patch fields (electrostatic potential nonuniformity on the cathode surface based on local work function values), and Schottky barrier lowering physics and illustrate that a smooth knee can arise from a thermionic cathode surface with as few as two discrete work function values. Importantly, we find that the inclusion of patch field effects is crucial for obtaining accurate J-T and J-V curves, and the further inclusion of Schottky barrier lowering is needed for accurate J-V curves. This finding, and the emission model provided in this paper have important implications for modeling electron emission from realistic, heterogeneous surfaces. Such modeling is important for improved understanding of the interplay of emission physics, cathode materials engineering, and design of numerous devices employing electron emission cathodes.
Ethyl-hexyl substituted polyfluorene (PF) with its high level of molecular disorder can be described very well by one-carrier space-charge-limited conduction for a discrete set of trap levels with energy $sim$ 0.5 eV above the valence band edge. Sweeping the bias above the trap-filling limit in the as-is polymer generates a new set of exponential traps, which is clearly seen in the density of states calculations. The trapped charges in the new set of traps have very long lifetimes and can be detrapped by photoexcitation. Thermal cycling the PF film to a crystalline phase prevents creation of additional traps at higher voltages.