Development of the Backflow Model for Simplified Impurity Exhaust in Monte-Carlo Calculation

The Monte-Carlo (MC) approach has a lot of flexibility in impurity transport modeling in the SOL and divertor region. However, in the divertor plasma simulation with the noble impurity seeding, characteristic time of the impurity transport especially in the sub-divertor chamber is long because the MC calculation of the impurity gas transport can be finished only by exhaust. The impurity MC calculation for such long exhaust processes is difficult in a series of the iterative calculation of a suite of integrated divertor codes SONIC. In order to overcome such a problem, a backflow model has been developed. Amount of the backflow flux from the sub-divertor chamber to the divertor region is evaluated in advance, and then simulating impurity flux is injected from the exhaust slot to the divertor region like a backflow. By this model, the MC calculation time is reduced significantly and iterative calculation of SONIC becomes possible within a reasonable calculation time. As a demonstration, the SONIC code with the backflow model has been applied to investigation of power handling in JT-60SA divertor. The SONIC simulation showed that low divertor heat load (< 10 MW/m²) with the low SOL density (< 1.5 × 10¹⁹m^-3), which is required in the full non-inductive current drive scenario, was achieved by the Ar gas puffing of 0.86 Pa m³/s.