Two-dimensional self-consistent numerical simulation of the whole discharge region in an atmospheric argon arc

A 2D self-consistent numerical model of the whole argon-arc discharge region that includes electrodes is developed in this work to facilitate analysis of the physical processes occurring in atmospheric arc plasma. The 2D arc column model contains the ionization and thermal non-equilibrium, which is coupled with a 1D electrode sheath model. The influence of plasma-species diffusion near the electrode region is investigated based on Maxwell–Stefan equations and the generalized Ohm's law. The numerical results of argon free-burning arcs at atmospheric pressure are then investigated. The simulation shows that the plasma is obviously in the state of thermal and ionization equilibrium in the arc core region, while it deviates from thermal and ionization equilibrium in the arc fringe region. The actual electron density decreases rapidly in the near-anode and near-cathode regions due to non-equilibrium ionization, resulting in a large electron number gradient in these regions. The results indicate that electron diffusion has an important role in the near-cathode and near-anode regions. When the anode arc root gradually contracts, it is easy to obtain a positive voltage drop of the anode sheath (I = 50 A), while it remains difficult to acquire a positive anode sheath voltage drop (I = 150 A). The current–voltage characteristics predicted by our model are found to be identical to the experimental values.