Numerical simulation for wire X-pinch plasma on 2D/3D geometry

Abstract

This paper investigates wire X-pinch (WXP) evolutions by the Eulerian resistive magneto-hydrodynamic code, STHENO, on 2D/3D geometry. A single-fluid two-temperature model is applied to pinch plasmas in local thermal equilibrium. The equation of state based on the Thomas–Fermi model is used to determine the ionization degree of the plasma. Electron internal energy is determined by the local density, temperature, and ionization potential with the average ion charge state. Lee–More–Desjarlais transport models are employed to obtain the thermal conductivities and resistivity for a non-ideal plasma. The radiation loss rate is calculated by the Bremsstrahlung and recombination emissivity within the ionization balance. The crossing point, which is the central part of the X-pinch, is assumed to be an axisymmetric configuration on a small computational domain in the RZ plane. The 2D simulation demonstrates that the micrometer size plasma column is elongated axially with the onset of the neck cascading structure. The radiation power is calculated and compared with the measured x-ray power from a modular X-pinch device (120-kA in 650-ns) at Seoul National University. The time evolution of the radiation power reproduces the trend of the measured x ray. 3D analyses are performed for the aluminum WXP configurations by varying wire numbers and cross-angles. The relation between the radiation performance and the numbers of wires reveals that the current density, rather than the line density, determines the central pinching condition. In addition, the multiple plasma instabilities (m = 0) near the central regions are found to degrade the radiation performance on the small cross-angle WXP.