Predictions of the Huber–Chankin (HC) scaling for the upstream impurity fraction were verified in a series of EDGE2D-EIRENE ('code') runs for highly radiating plasmas with nitrogen injection. The main quantity extracted from the code was poloidally averaged, from X-point to X-point, separatrix impurity fraction cZ in the main scrape-off layer (SOL). Variation of the main working gas (H, D and T) revealed a qualitative agreement between the model and code results owing to the very large isotope difference in the predicted cZ values caused primarily by the inverse isotope mass dependence of the H-mode power threshold assumed in the HC model and implemented in the code. At the same time, the variation of the toroidal field and safety factor in deuterium cases yielded no correlation between the model predictions and code results. The code showed much higher local impurity fractions (fZ) in the divertor compared to the main SOL, as well as large case-to-case variations in the divertor to the main SOL ratio of impurity fractions. The analysis of code results has wide-ranging consequences not only for the HC model, but also for other similar 1D models which use simple geometry ignoring strong neutral recycling in the divertor/ Different topology makes plasma parameters in the divertor and main SOL very different, resulting in different impurity charge state composition. Missing mechanisms in 1D codes (e.g. friction and thermo-forces exerted on impurity ions by main working gas ions) lead to impurity density redistribution. Neglecting all above factors, 1D models assume a constant impurity fraction along field lines.
Evaluation of SPARC divertor conditions in H-mode operation using SOLPS-ITER