SOLPS-ITER predictive simulations of the impact of ion-molecule elastic collisions on strongly detached MAST-U Super-X divertor conditions

The role of ion-molecule ( D+ − D2 ) elastic collisions in strongly detached divertor conditions has been studied in the MAST-U Super-X configuration using SOLPS-ITER. Two strongly detached steady state solutions were compared, one obtained through a main-ion fuelling scan and the other through a nitrogen seeding scan at fixed fuelling rate. A significant difference in the electron-ion recombination (EIR) levels was observed; significant EIR in strongly detached conditions in the fuelling scan and negligible EIR throughout the seeding scan. This is partly because the fuelling scan achieves electron temperatures ( Te ) as low as 0.2 eV near the divertor target, compared to 0.8 eV in the seeding scan (EIR increases strongly below Te ≈ 1 eV), and partly due to higher divertor plasma densities achieved in fuelling scan. Features of the strongly detached seeded cases, i.e. higher temperatures and negligible EIR, are recovered in the fuelling scan by turning off D+ − D2 elastic collisions. Analysis suggests that dissipation mechanisms like line radiation and charge exchange (important for detachment initiation) become weak when Te falls below 1 eV, and that D+ − D2 elastic collisions are necessary for further heat dissipation and access to strongly recombining conditions in the fuelling scan. In the seeding scan, heat dissipation through D+ − D2 elastic collisions is weak. This could be because our nitrogen seeding simulations do not include interactions between nitrogen ions and neutrals, and the strongly detached cases contain high levels of N+ in the divertor. As a result, the N+ acts like a reservoir of energy and momentum which appears to weaken the impact of D+ − D2 elastic collisions on the divertor plasma energy and momentum balance, making it more difficult to access recombining conditions. This suggests that some of the differences between seeding and fuelling scans could be because energy and momentum exchange between impurities and neutrals is not sufficiently captured in our simulations.