1D fluid model of LHD divertor plasma and hydrogen recycling

One dimensional plasma and neutral models of the divertor plasma in Large Helical Device is presented. The plasma is described by stationary fluid equations for electron and ion. The atomic processes such as dissociation of hydrogen molecules released from the divertor plate, ionization of hydrogen atoms, charge exchange and recombination are included in equations of neutrals. This model is intended to be employed in an integrated simulation where an equilibrium of the upstream plasma and plasma-surface interactions at the divertor plate are solved in different numerical codes separately. From the computational point of view, the numerical code for the divertor plasma is developed for 1D flux tube where the boundary conditions of both ends are specified. The calculation time is less than one second and reasonably short to use in future integrated simulations. In the results, interactions between plasma and neutrals and dependence of the energy loss on the plasma density are studied. In low density case, the energy is lost through ionization and charge exchange but the total amount of the loss is small and the impurity loss is negligibly small. In high density, or high recycling case, the ionization loss and impurity cooling increases much larger than the charge exchange loss and causes the drop of the heat flux at the divertor plate. (author)