OEDGE modeling for the planned tungsten ring experiment on DIII-D

Abstract

The OEDGE code is used to model tungsten erosion and transport for experiments with toroidal rings of high-Z metal tiles in the DIII-D tokamak. Such modeling is needed for both experimental and diagnostic design to have estimates of the expected core and edge tungsten density and to understand the various factors contributing to the uncertainties in these calculations. OEDGE simulations are performed using the planned experimental magnetic geometries and plasma conditions typical of both L-mode and inter-ELM H-mode discharges in DIII-D. OEDGE plasma reconstruction based on specific representative discharges for similar geometries is used to determine the plasma conditions applied to tungsten plasma impurity simulations. A new model for tungsten erosion in OEDGE was developed which imports charge-state resolved carbon impurity fluxes and impact energies from a separate OEDGE run which models the carbon production, transport and deposition for the same plasma conditions as the tungsten simulations. These values are then used to calculate the gross tungsten physical sputtering due to carbon plasma impurities which is then added to any sputtering by deuterium ions; tungsten self-sputtering is also included. The code results are found to be dependent on the following factors: divertor geometry and closure, the choice of cross-field anomalous transport coefficients, divertor plasma conditions (affecting both tungsten source strength and transport), the choice of tungsten atomic physics data used in the model (in particular ionization rate for W-atoms), and the model of the carbon flux and energy used for calculating the tungsten source due to sputtering. Core tungsten density is found to be of order 1015m−3 (excluding effects of any core transport barrier and with significant variability depending on the other factors mentioned) with density decaying into the scrape off layer. For the typical core density in the plasma conditions examined of 2 to 4 ×1019m−3, this represents a concentration on the order of 5 × 10−5.