Efficient algorithm for simulating a multidimensional glow discharge

A normal glow discharge in nitrogen between two infinite plane-parallel electrodes is simulated numerically in three dimensions in the drift-diffusion model. An algorithm is proposed that is based on solving the model equations by a hybrid of the sweep and relaxation methods and that combines a satisfactorily fast rate of convergence with the efficiency of individual iterations. The algorithm developed was used to carry out a three-dimensional simulation of a discharge in nitrogen at the pressures p = 5 and 1 Torr, the source voltage and ballast resistance being 2 kV and 300 kΩ, respectively. It was found that, at the pressure p = 1 Torr, the current density and charged particle densities change to a torioidal distribution in the anode region.