Study of inertial electrostatic confinement fusion using a finite-volume scheme for the one-dimensional Vlasov...

While the majority of fusion energy research is focused on magnetic confinement, there have been several alternative confinement methods aimed at the development of smaller and less expensive reactors. A number of these alternative reactors are based on a spherically convergent beam of recirculating ions and include designs such as inertial electrostatic confinement (IEC), multigrid IEC, and the periodically oscillating plasma sphere concept. Here, a fully time-dependent GPU-based Vlasov solver was developed in order to study these spherically convergent devices. This code solves the Vlasov equation for a spherically symmetric system using a finite-volume method with a modified flux to account for electrode transparency. The solver accounts for secondary electron emission, interactions between the charged particles, and collisional effects such as ionization and charge exchange. This code was used to investigate a system similar to the ion-injected device described by Hirsch (see [R. L. Hirsch, J. Appl. Phys. 38, 4522 (1967)10.1063/1.1709162]), who had reported a neutron production rate for deuterium-deuterium reactions in the range of 10^{6}to10^{7} neutrons per second, which was attributed to the formation of a virtual electrode structure near the center of the chamber. Attempts to reproduce this experiment [B. J. Egle, Ph.D. thesis, 2010] yielded similar fusion rates, though the majority of the reactions were found not to occur near the center of the chamber. The results of this Vlasov solver, considering only beam-beam and beam-background fusion reactions, show that beam-background reactions would be dominant in such an ion-injected device. This result is consistent with work by Baxter and Stuart, who proposed a simplified steady-state Boltzmann model. However, the result of both models are inconsistent with the experimental results, which indicate a higher neutron production rate, and an inverse pressure scaling trend. It is shown that the higher experimental rates may be explained by beam-target fusion between the ion beam and deuterium embedded on the inner surface of the cathode.