Formation of Self-Organized Anode Patterns in Arc Discharge Simulations

Pattern formation and self-organization are phenomena commonly observed experimentally in diverse types of plasma systems, including atmospheric-pressure electric arc discharges. However, numerical simulations reproducing anode pattern formation in arc discharges have proven exceedingly elusive. Time-dependent three-dimensional thermodynamic nonequilibrium simulations reveal the spontaneous formation of self-organized patterns of anode attachment spots in the free-burning arc, a canonical thermal plasma flow established by a constant DC current between an axi-symmetric electrodes configuration in the absence of external forcing. The number of spots, their size, and distribution within the pattern depend on the applied total current and on the resolution of the spatial discretization, whereas the main properties of the plasma flow, such as maximum temperatures, velocity, and voltage drop, depend only on the former. The sensibility of the solution to the spatial discretization stresses the computational requirements for comprehensive arc discharge simulations. The obtained anode patterns qualitatively agree with experimental observations and confirm that the spots originate at the fringes of the arc - anode attachment. The results imply that heavy-species - electron energy equilibration, in addition to thermal instability, has a dominant role in the formation of anode spots in arc discharges.

💡 Research Summary

This paper tackles a long‑standing challenge in plasma physics: reproducing the spontaneous formation of self‑organized anode attachment spots in a free‑burning electric arc using numerical simulations. The authors employ a fully three‑dimensional, time‑dependent, thermodynamic nonequilibrium model that treats electrons and heavy species (argon atoms) as separate energy carriers. The governing equations include continuity, momentum, separate electron and heavy‑species energy balances, and Maxwell’s equations, all solved without any imposed external forcing. The computational domain represents a canonical axial‑symmetric electrode configuration with a constant DC current flowing between a cathode and an anode.

A systematic parametric study is performed by varying the total current (30 A, 60 A, 90 A) and the spatial discretization (grid spacings of 0.5 mm, 0.25 mm, and 0.125 mm). The results reveal a pronounced sensitivity of the anode spot pattern to the mesh resolution: coarse grids capture only a few large spots, whereas fine grids resolve dozens of smaller, more uniformly distributed spots. As the current increases, the number of spots grows while their individual size shrinks, indicating that higher currents promote more frequent local overheating events. Conversely, macroscopic plasma properties—maximum temperature, average axial velocity, and overall voltage drop—depend almost exclusively on the applied current and remain essentially unchanged across different mesh resolutions. This decoupling demonstrates that the overall energy balance of the arc is governed by the global current, while the fine‑scale organization of attachment spots is dictated by local electron‑heavy‑species energy exchange dynamics.

The simulations show that spots originate at the fringe of the arc where it meets the anode surface. In these fringe regions, the electron–heavy‑species equilibration time becomes locally short, leading to rapid temperature spikes. The resulting localized increase in electrical conductivity focuses current density and further amplifies heating, creating a positive feedback loop that stabilizes the spot. This mechanism integrates two previously proposed ideas—thermal instability of the plasma column and surface charge instability of the electrode—into a unified picture where heavy‑species–electron energy equilibration plays a dominant role.

Qualitative comparison with experimental observations confirms that the simulated spot numbers, sizes, and spatial distributions closely resemble those seen in laboratory arcs. The agreement validates the model’s capability to capture real‑world physics despite the extreme computational demands. The authors also discuss the practical implications of their findings: accurate prediction of anode spot patterns is essential for the design of high‑current electrodes, arc welding processes, and plasma‑based material processing, where localized overheating can lead to electrode erosion or process instability.

Finally, the paper highlights the computational cost associated with high‑resolution three‑dimensional nonequilibrium simulations. The authors suggest that future work should explore adaptive mesh refinement, reduced‑order modeling, or hybrid kinetic‑fluid approaches to make such detailed simulations more tractable for engineering applications. In summary, this study provides the first comprehensive 3‑D time‑dependent simulation that reproduces self‑organized anode spot patterns in a free‑burning arc, elucidates the underlying physics, and bridges the gap between experimental observations and numerical modeling.