DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Prediction of divertor heat flux width for ITER pre-fusion power operation using BOUT++ transport code

Abstract

Prediction of divertor heat flux width is performed for the first and the second pre-fusion power operation (PFPO) phases specified in the new ITER research plan using BOUT++ transport code (Li et al 2018 Comput. Phys. Commun.228 69–82). Here the initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients in transport code are calculated by inverting the plasma profiles inside the separatrix and are assumed to be constants in the scrape-off-layer. An anomalous thermal diffusivity scan is performed with E × B and magnetic drifts. The results in two scenarios identifying two distinct regimes: a drift-dominant regime when diffusivity is smaller than the respective critical thermal diffusivity χc and a turbulence-dominant regime when diffusivity is larger than it. The Goldston heuristic drift model and the ITPA multi-machine experimental scaling yield a lower limit of the width λq. From transport simulations, we obtain the critical thermal diffusivity χc = 0.5 m2 s–1 for the PFPO-1 scenario with toroidal magnetic field B = 1.77 T and plasma current Ip = 5 MA, and χc = 0.3 m2 s–1 for the PFPO-2 scenario with toroidal magnetic field B = 2.65 T and plasma current Ip = 7.5 MA. Separatrix temperature and collisionality also have a significant impact on the heat flux width in the drift-dominant regime. The investigation clearly yields a scaling for critical thermal diffusivity $${\chi }_{\text{c}}\propto {A}^{1/2}/(Z{\left(1+Z\right)}^{\frac{1}{2}}{B}_{\text{p}}^{2})$$ using ITER scenarios with fixed safety factor q95, major radius R, aspect ratio R/a, and the separatrix temperature Tsep, and establishes the connection with CFETR and C-Mod discharges. This scaling implies that for a given tokamak device with q95, R, R/a, and Tsep fixed, a reduction of poloidal magnetic field by a factor of 3 leads to a 9 times higher critical value of thermal diffusivity χc, possibly yielding a transition from turbulence- to drift-dominant regime.

Authors:
ORCiD logo [1]; ORCiD logo [2]; ORCiD logo [3]; ORCiD logo [2];  [4]
  1. Dalian Univ. of Technology (China); Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  2. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  3. Oak Ridge Associated Univ., Oak Ridge, TN (United States)
  4. Dalian Univ. of Technology (China)
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA); Chinese Scholarship Council; National Key Research and Development Program of China
OSTI Identifier:
1872278
Report Number(s):
LLNL-JRNL-815224
Journal ID: ISSN 0029-5515; 1024415; TRN: US2306856
Grant/Contract Number:  
AC52-07NA27344; 201906060127; 2018YFE0303102
Resource Type:
Accepted Manuscript
Journal Name:
Nuclear Fusion
Additional Journal Information:
Journal Volume: 62; Journal Issue: 5; Journal ID: ISSN 0029-5515
Publisher:
IOP Science
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; ITER; heat flux width; transport; drift

Citation Formats

He, Xiaoxue X., Xu, Xueqiao Q., Li, Zeyu Y., Zhu, Ben, and Liu, Y. Prediction of divertor heat flux width for ITER pre-fusion power operation using BOUT++ transport code. United States: N. p., 2022. Web. doi:10.1088/1741-4326/ac3e80.
He, Xiaoxue X., Xu, Xueqiao Q., Li, Zeyu Y., Zhu, Ben, & Liu, Y. Prediction of divertor heat flux width for ITER pre-fusion power operation using BOUT++ transport code. United States. https://doi.org/10.1088/1741-4326/ac3e80
He, Xiaoxue X., Xu, Xueqiao Q., Li, Zeyu Y., Zhu, Ben, and Liu, Y. Wed . "Prediction of divertor heat flux width for ITER pre-fusion power operation using BOUT++ transport code". United States. https://doi.org/10.1088/1741-4326/ac3e80. https://www.osti.gov/servlets/purl/1872278.
@article{osti_1872278,
title = {Prediction of divertor heat flux width for ITER pre-fusion power operation using BOUT++ transport code},
author = {He, Xiaoxue X. and Xu, Xueqiao Q. and Li, Zeyu Y. and Zhu, Ben and Liu, Y.},
abstractNote = {Prediction of divertor heat flux width is performed for the first and the second pre-fusion power operation (PFPO) phases specified in the new ITER research plan using BOUT++ transport code (Li et al 2018 Comput. Phys. Commun.228 69–82). Here the initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients in transport code are calculated by inverting the plasma profiles inside the separatrix and are assumed to be constants in the scrape-off-layer. An anomalous thermal diffusivity scan is performed with E × B and magnetic drifts. The results in two scenarios identifying two distinct regimes: a drift-dominant regime when diffusivity is smaller than the respective critical thermal diffusivity χc and a turbulence-dominant regime when diffusivity is larger than it. The Goldston heuristic drift model and the ITPA multi-machine experimental scaling yield a lower limit of the width λq. From transport simulations, we obtain the critical thermal diffusivity χc = 0.5 m2 s–1 for the PFPO-1 scenario with toroidal magnetic field B = 1.77 T and plasma current Ip = 5 MA, and χc = 0.3 m2 s–1 for the PFPO-2 scenario with toroidal magnetic field B = 2.65 T and plasma current Ip = 7.5 MA. Separatrix temperature and collisionality also have a significant impact on the heat flux width in the drift-dominant regime. The investigation clearly yields a scaling for critical thermal diffusivity ${\chi }_{\text{c}}\propto {A}^{1/2}/(Z{\left(1+Z\right)}^{\frac{1}{2}}{B}_{\text{p}}^{2})$ using ITER scenarios with fixed safety factor q95, major radius R, aspect ratio R/a, and the separatrix temperature Tsep, and establishes the connection with CFETR and C-Mod discharges. This scaling implies that for a given tokamak device with q95, R, R/a, and Tsep fixed, a reduction of poloidal magnetic field by a factor of 3 leads to a 9 times higher critical value of thermal diffusivity χc, possibly yielding a transition from turbulence- to drift-dominant regime.},
doi = {10.1088/1741-4326/ac3e80},
journal = {Nuclear Fusion},
number = 5,
volume = 62,
place = {United States},
year = {Wed Mar 16 00:00:00 EDT 2022},
month = {Wed Mar 16 00:00:00 EDT 2022}
}

Works referenced in this record:

Simulation of divertor heat flux width on EAST by BOUT++ transport code
journal, July 2020


Turbulent transport regimes and the scrape-off layer heat flux width
journal, April 2015

  • Myra, J. R.; D'Ippolito, D. A.; Russell, D. A.
  • Physics of Plasmas, Vol. 22, Issue 4
  • DOI: 10.1063/1.4919255

ITER divertor plasma response to time-dependent impurity injection
journal, August 2017


Tungsten dust impact on ITER-like plasma edge
journal, January 2015

  • Smirnov, R. D.; Krasheninnikov, S. I.; Pigarov, A. Yu.
  • Physics of Plasmas, Vol. 22, Issue 1
  • DOI: 10.1063/1.4905704

Edge-localized-modes in tokamaks
journal, September 2014


Simulations of divertor heat flux width using transport code with cross-field drifts under the BOUT++ framework
journal, January 2020

  • Li, N. M.; Xu, X. Q.; Hughes, J. W.
  • AIP Advances, Vol. 10, Issue 1
  • DOI: 10.1063/1.5126884

Simulations of tokamak boundary plasma turbulence transport in setting the divertor heat flux width
journal, October 2019


Progress toward ITER’s First Plasma
journal, June 2019


Calculation of two-dimension radial electric field in boundary plasmas by using BOUT++
journal, July 2018


Edge localized modes (ELMs)
journal, February 1996


BOUT++: A framework for parallel plasma fluid simulations
journal, September 2009


Investigation of key parameters for the development of reliable ITER baseline operation scenarios using CORSICA
journal, March 2018


Regime of Improved Confinement and High Beta in Neutral-Beam-Heated Divertor Discharges of the ASDEX Tokamak
journal, November 1982


Modelling one-third field operation in the ITER pre-fusion power operation phase
journal, September 2019


Prediction of divertor heat flux width for ITER using BOUT++ transport and turbulence module
journal, February 2019


Toroidal modeling of resonant magnetic perturbations in preparation for the initial phase of ITER operation
journal, November 2019


Gyrokinetic projection of the divertor heat-flux width from present tokamaks to ITER
journal, August 2017


Divertor heat flux simulations in ELMy H-mode discharges of EAST
journal, August 2017


Drift reduced Landau fluid model for magnetized plasma turbulence simulations in BOUT++ framework
journal, October 2021


Heuristic drift-based model of the power scrape-off width in low-gas-puff H-mode tokamaks
journal, December 2011


Six-field two-fluid simulations of peeling–ballooning modes using BOUT++
journal, May 2013


Progress towards modeling tokamak boundary plasma turbulence and understanding its role in setting divertor heat flux widths
journal, May 2018

  • Chen, B.; Xu, X. Q.; Xia, T. Y.
  • Physics of Plasmas, Vol. 25, Issue 5
  • DOI: 10.1063/1.5016582

Consequences of a reduction of the upstream power SOL width in ITER
journal, July 2013


SOLPS-ITER modelling of ITER edge plasma with drifts and currents
journal, March 2020


CORSICA modelling of ITER hybrid operation scenarios
journal, August 2016


Scaling of the tokamak near the scrape-off layer H-mode power width and implications for ITER
journal, August 2013


Impact of plasma density/collisionality on divertor heat flux width
journal, December 2020


Presentation of the New SOLPS-ITER Code Package for Tokamak Plasma Edge Modelling
journal, January 2016

  • Bonnin, Xavier; Dekeyser, Wouter; Pitts, Richard
  • Plasma and Fusion Research, Vol. 11, Issue 0
  • DOI: 10.1585/pfr.11.1403102

Scaling of L-mode heat flux for ITER and COMPASS-U divertors, based on five tokamaks
journal, May 2020


Scrape-off layer plasmas for ITER with 2nd X-point and convective transport effects
journal, June 2007


Low-to-high confinement transition simulations in divertor geometry
journal, May 2000

  • Xu, X. Q.; Cohen, R. H.; Rognlien, T. D.
  • Physics of Plasmas, Vol. 7, Issue 5
  • DOI: 10.1063/1.874044

Experimental study of heating scheme effect on the inner divertor power footprint widths in EAST lower single null discharges
journal, February 2018


Study of plasma current effect on divertor power footprint widths through experiments and modeling in EAST L-mode plasmas
journal, April 2017

  • Deng, G. Z.; Liu, X. J.; Wang, L.
  • Physics of Plasmas, Vol. 24, Issue 4
  • DOI: 10.1063/1.4979670

Assessment of ITER divertor performance during early operation phases
journal, December 2020


Principal physics developments evaluated in the ITER design review
journal, May 2009


Development of ITER non-activation phase operation scenarios
journal, June 2017


Analysis of pedestal plasma transport
journal, May 2010


Edge-localized modes - physics and theory
journal, May 1998


2D simulations of transport dynamics during tokamak fuelling by supersonic molecular beam injection
journal, March 2014


Impurity seeding in ITER DT plasmas in a carbon-free environment
journal, August 2015


Finalizing the ITER divertor design: The key role of SOLPS modeling
journal, December 2011


MHD stability of the ITER pedestal and SOL plasma and its influence on the heat flux width
journal, August 2015


Modeling of Tungsten Dust Transport in ITER with Multi-Physics Code DUSTT/UEDGE: Modeling of Tungsten Dust Transport in ITER with Multi-Physics Code DUSTT/UEDGE
journal, June 2014

  • Pigarov, A. Yu.; Smirnov, R. D.; Krasheninnikov, S. I.
  • Contributions to Plasma Physics, Vol. 54, Issue 4-6
  • DOI: 10.1002/ctpp.201410088

Chapter 4: Power and particle control
journal, June 2007