Dense plasma focus (DPF) devices are compact pulsed plasma systems capable of producing intense bursts of neutrons, X-rays and energetic particles, and their performance is strongly influenced by the geometric configuration of the electrodes. In this work, a two-dimensional relaxation magnetohydrodynamic (MHD) model is employed to systematically investigate the effects of key geometric parameters on the discharge dynamics and neutron production characteristics of a deuterium-filled DPF device, with particular emphasis on identifying optimal parameter ranges and elucidating the underlying physical mechanisms.
The numerical model self-consistently describes the macroscopic evolution of the plasma sheath driven by electromagnetic forces after breakdown, including the axial acceleration along the anode, the subsequent radial implosion toward the axis, and the formation of a transient high-density plasma column. The ionization degree is incorporated through a Saha equilibrium correction to ensure thermodynamic consistency of the plasma parameters, while the neutron yield is evaluated by considering only the thermonuclear D–D fusion reactions. Kinetic effects such as beam–target interactions and non-Maxwellian ion populations are not explicitly included, and their possible influence is discussed in the context of the model applicability.
For the optimal geometric configuration, the simulated discharge current, neutron production rate and plasma thermodynamic quantities exhibit clear temporal correlations. The evolution of the discharge current shows a rapid rise followed by a peak near the end of the radial compression stage, accompanied by a pronounced burst of neutron production concentrated within a narrow time window. This indicates that neutron generation is closely associated with the formation of a dense and hot plasma region during the late implosion stage. The neutron production rate spans several orders of magnitude over the discharge duration, highlighting the highly transient nature of the fusion process.
To assess the impact of ionization modeling, the temporal evolutions of the ionization degree and characteristic plasma temperature with and without the Saha correction are compared. The results show that, during the partially ionized phase, the Saha-corrected temperature is systematically higher than that obtained under the assumption of full ionization. This behavior arises from the redistribution of internal energy due to ionization equilibrium effects, which modify the effective particle number density. As the discharge progresses and the ionization degree approaches unity, the two temperature definitions gradually converge, indicating that ionization effects become less significant in the highly ionized regime. The characteristic temperature, defined as the average temperature at the plasma sheath front, exhibits a distinct peak concurrent with the neutron burst, reflecting strong energy localization during the compression process.
The simulated two-dimensional ion density distributions further reveal the staged evolution of the plasma, from axial motion of the current sheath to radial convergence and the formation of a dense axial structure. The spatial and temporal features obtained in the present simulations are in good qualitative agreement with typical DPF numerical and experimental observations reported in the literature.
Overall, although the present MHD model does not explicitly account for kinetic beam–target mechanisms, the results demonstrate that geometric parameters exert a decisive influence on the macroscopic compression dynamics and thermonuclear neutron production trends. The conclusions of this study are therefore expected to provide meaningful guidance for electrode design and parameter optimization of DPF devices, particularly with respect to improving plasma compression efficiency and fusion performance under conditions where thermal processes play a significant role.