Self-consistent treatment of the sheath boundary conditions by introducing anisotropic ion temperatures and virtual divertor model

One-dimensional SOL-divertor plasma fluid simulation code which considers anisotropy of ion temperature has been developed so as to deal with sheath theory self-consistently. In our fluid modeling, explicit use of boundary condition for Mach number M at divertor plate, e.g., M = 1 , becomes unnecessary. In order to deal with the Bohm condition and the sheath heat transmission factors at divertor plate self-consistently, we introduced a virtual divertor (VD) model which sets an artificial region beyond divertor plates and artificial sinks for particle, momentum and energy there to model the effects of the sheath region in front of the divertor plate. Validity of our fluid model with VD model is confirmed by showing that simulation results agree well with those from a kinetic code regarding the Bohm condition, ion temperature anisotropy and supersonic flow. We also show that the strength of artificial sinks in VD region does not affect profiles in plasma region at least in the steady state and that sheath heat transmission factors can be adjusted to theoretical values by VD model. Validity of viscous flux is also investigated.

[1]  Kazuo Hoshino,et al.  Simulation study of power load with impurity seeding in advanced divertor “short super-X divertor” for a tokamak reactor , 2015 .

[2]  Bruce D. Scott,et al.  Free-energy conservation in local gyrofluid models , 2005 .

[3]  T. Takizuka,et al.  Kinetic modelling of impurity transport in detached plasma for integrated divertor simulation with SONIC (SOLDOR/NEUT2D/IMPMC/EDDY) , 2009 .

[4]  Boris B. Kadomtsev,et al.  Reviews of Plasma Physics , 2012 .

[5]  Chia-Jung Hsu.  Numerical Heat Transfer and Fluid Flow , 1981 .

[6]  B. Braams,et al.  Simulation of the Edge Plasma in Tokamaks , 2004 .

[7]  T. Takizuka,et al.  Particle simulation of detached plasma in the presence of diffusive particle loss and radiative energy loss , 2001 .

[8]  T. Takizuka,et al.  SOL–divertor plasma simulations introducing anisotropic temperature with virtual divertor model , 2015 .

[9]  Philippe Angot,et al.  A penalization method to take into account obstacles in incompressible viscous flows , 1999, Numerische Mathematik.

[10]  William H. Press,et al.  Book-Review - Numerical Recipes in Pascal - the Art of Scientific Computing , 1989 .

[11]  H. Berk,et al.  An alternative derivation of the parallel ion viscosity , 1982 .

[12]  J. Milovich,et al.  A fully implicit, time dependent 2-D fluid code for modeling tokamak edge plasmas , 1992 .

[13]  T. Takizuka,et al.  Kinetic Particle Simulation Study of Parallel Heat Transport in Scrape-off Layer Plasmas over a Wide Range of Collisionalities , 2010 .

[14]  T. Takizuka,et al.  Particle simulation of divertor plasma , 1984 .

[15]  P. Ghendrih,et al.  Penalization technique to model wall-component impact on heat and mass transport in the tokamak edge , 2013 .

[16]  T. Takizuka,et al.  Improvement of the detachment modelling in the SONIC simulation , 2011 .

[17]  W. Fundamenski.  Parallel heat flux limits in the tokamak scrape-off layer , 2005 .

[18]  Gregory W. Hammett,et al.  Toroidal gyrofluid equations for simulations of tokamak turbulence , 1996 .

[19]  T. Takizuka,et al.  Particle Simulation Study of the Effect of Radial Electric Field on Scrape‐off Layer Plasma and Sheath Formation , 2000 .

[20]  D. Coster,et al.  Assessment of edge modeling in support of ITER , 2011 .

[21]  Guido Ciraolo,et al.  Penalization modeling of a limiter in the Tokamak edge plasma , 2010, J. Comput. Phys..

[22]  Martine Baelmans,et al.  The new SOLPS-ITER code package , 2015 .

[23]  T. Takizuka,et al.  Simulation Study of an Extended Divertor Leg for Heat Control in the SlimCS DEMO Reactor , 2012 .

[24]  D. Coster,et al.  Comparison of 2D models for the plasma edge with experimental measurements and assessment of deficiencies , 2009 .

[25]  N. Ohyabu,et al.  One-dimensional simulation on stability of detached plasma in a tokamak divertor , 1999 .

[26]  F. Najmabadi,et al.  Generalized fluid equations for parallel transport in collisional to weakly collisional plasmas , 1986 .

[27]  B. Braams,et al.  Plasma Edge Physics with B2‐Eirene , 2006 .

[28]  F. Imbeaux,et al.  Numerical study of transition to supersonic flows in the edge plasma , 2014 .

[29]  S. I. Braginskii.  Reviews of Plasma Physics , 1965 .

[30]  Philippe Ghendrih,et al.  The Plasma Boundary of Magnetic Fusion Devices , 2001 .

[31]  T. Takizuka,et al.  Development of Integrated SOL/Divertor Code and Simulation Study in JAEA , 2006 .

[32]  A. Taroni,et al.  Models and Numerics in the Multi-Fluid 2-D Edge Plasma Code EDGE2D/U , 1994 .

[33]  J. Contributors,et al.  Influence of cross-field drifts and chemical sputtering on simulations of divertor particle and heat loads in ohmic and L-mode plasmas in DIII-D, AUG, and JET using UEDGE , 2010 .

[34]  Y. Ogawa,et al.  “Multi-layer” one-dimensional model for stability analysis on partially detached divertor plasmas , 2011 .

[35]  W. Press.  Numerical recipes in Fortran 77 : the art of scientific computing : volume 1 of fortran numerical recipes , 1996 .

[36]  J. Bucalossi,et al.  Investigation of Edge and SOL Particle Flux Patterns in High Density Regimes using SOLEDGE2D‐EIRENE Code , 2014 .

[37]  D. P. Coster,et al.  Detachment physics in SOLPS simulations , 2011 .