Since 18 of December 2019 uses Nucleus credentials. Visit our help pages for information on how to Register and Sign-in using Nucleus.

Nov 4 – 7, 2019
IAEA Headquarters, Vienna, Austria
Europe/Vienna timezone
Meeting Material is now available and accessible from the left-menu

Status of Divertor/SOL modelling in PROCESS

Nov 6, 2019, 4:30 PM
Board Room C (C Building - 4th Floor) (IAEA Headquarters, Vienna, Austria)

Board Room C (C Building - 4th Floor)

IAEA Headquarters, Vienna, Austria

Poster Scrape-off Layer and Divertor Physics Poster Session III


Dr Stuart Muldrew (CCFE, UKAEA)


Systems codes, such as PROCESS, model all systems of a power plant to investigate large numbers of design points. They are used for scoping studies and to identify areas of feasible design points.

Multi-dimensional modelling of the plasma Scrape-Off-Layer (SOL), divertor and seeded impurities is too computationally intensive to be incorporated directly into a systems code. Divertor protection parameters such as $P_{\rm sep}/R_{\rm 0}$ and $P_{\rm sep}B_{\rm T}/q_{\rm 95}AR_{\rm 0}$ have been used as a constrait for capturing the divertor problem in previous studies instead. A 1-D SOL/divertor model has been implemented in PROCESS to try and produce more accurate information regarding the divertor conditions. The aim of the 1-D model is to determine if the divertor is detached, whether the power crossing the separatrix is consistent with required conditions at the target, and to model the loss of power and momentum along the 1-D flux tube.

The following physical processes are included: convected heat flux; thermal conduction; momentum conservation; radiation by deuterium, tritium and impurities; charge exchange; electron impact ionisation; and surface recombination. Pumping is not included – all particles striking the target are recycled. The strong shearing of the flux tube near the X-point is not taken into account.

As the seeded impurity concentration is increased a discontinuous transition is observed between an attached state, where the plasma temperature at the target is $50\, \rm{eV}$, and a state where the temperature at the target hits the lower bound of the simulation, $1.1\, \rm{eV}$. We interpret this as a detached state, within the limitations of the model.

The 1-D model has been compared to 2-D models (e.g. the Japanese code SONIC) for European DEMO-like machines. However, a large database of DEMO-like runs using the detailed codes is not readily available, so benchmarking the 1-D model against detailed codes is an ongoing process.

PROCESS now also allows for a double-null divertor configuration as an alternative to the single-null considered standard for conventional aspect ratio tokamaks, but to achieve worthwhile power sharing between the upper and lower divertors would seem to require a high degree of control of the plasma position.

This work has been carried out within the framework of the EUROfusion Consortium and has received funding from the Euratom research and training programme 2014-2018 and 2019-2020 under grant agreement No 633053 and from the RCUK [grant number EP/P012450/1]. The views and opinions expressed herein do not necessarily reflect those of the European Commission.

Country or International Organization United Kingdom

Primary author

Dr Stuart Muldrew (CCFE, UKAEA)


Mr Michael Kovari (CCFE, UKAEA) Mr James Morris (CCFE, UKAEA)

Presentation materials