A hybrid fluid-kinetic model for hydrogenic atoms in the plasma edge of tokamaks based on a micro-macro decomposition of the kinetic equation

Abstract Monte Carlo (MC) simulations of the full kinetic equation for the neutral particles in the plasma edge become computationally costly for reactor-relevant regimes. To accelerate the simulations, we propose a hybrid fluid-kinetic approach that is based on a micro-macro decomposition of the kinetic equation. This leads to a macro/fluid model with kinetic corrections that follow from an MC simulation of the micro/kinetic part. We distinguish three hybrid models with different underlying fluid equations: (i) a pure pressure-diffusion equation with equal neutral and ion temperatures ( T n = T i ); (ii) a continuity and parallel momentum equation with pressure-diffusion transport retained in the directions perpendicular to the magnetic field lines, with T n = T i ; and (iii) the same model as (ii), but with a separate neutral energy equation ( T n ≠ T i ). To facilitate the future integration in more complete plasma edge codes, we neglect some kinetic correction terms. Hence, the hybrid model is not exactly equivalent to the full kinetic equation. We assess the hybrid performance on the basis of the reduction of the CPU time compared to an MC simulation of the full kinetic equation for the same statistical error on a certain plasma source. This is done for a high recycling slab case. Only the models with parallel momentum equation ((ii)-(iii)) are able to significantly reduce the CPU time. However, due to the incomplete kinetic corrections there is a remaining hybrid-kinetic discrepancy that mainly pops up in the ion energy source from the model with energy equation (iii).

[1]  G. Samaey,et al.  An analytical study of Monte Carlo source term estimators in plasma edge simulations of fusion reactors , 2017 .

[2]  Ronald L. Wasserstein,et al.  Monte Carlo: Concepts, Algorithms, and Applications , 1997 .

[3]  Shi Jin,et al.  Efficient Asymptotic-Preserving (AP) Schemes For Some Multiscale Kinetic Equations , 1999, SIAM J. Sci. Comput..

[4]  Lorenzo Pareschi,et al.  Implicit-Explicit Runge-Kutta Schemes for Hyperbolic Systems and Kinetic Equations in the Diffusion Limit , 2013, SIAM J. Sci. Comput..

[5]  M. Baelmans,et al.  Accuracy and convergence of coupled finite-volume/Monte Carlo codes for plasma edge simulations of nuclear fusion reactors , 2016, J. Comput. Phys..

[6]  Luc Mieussens,et al.  A New Asymptotic Preserving Scheme Based on Micro-Macro Formulation for Linear Kinetic Equations in the Diffusion Limit , 2008, SIAM J. Sci. Comput..

[7]  N. Crouseilles,et al.  Kinetic/fluid micro-macro numerical schemes for Vlasov-Poisson-BGK equation using particles , 2012 .

[8]  Shi Jin,et al.  Micro-macro decomposition based asymptotic-preserving numerical schemes and numerical moments conservation for collisional nonlinear kinetic equations , 2018, J. Comput. Phys..

[9]  Pierre Degond,et al.  A hybrid kinetic/fluid model for solving the gas dynamics Boltzmann-BGK equation , 2004 .

[10]  Luc Mieussens,et al.  Uniformly stable numerical schemes for the Boltzmann equation preserving the compressible Navier-Stokes asymptotics , 2008, J. Comput. Phys..

[11]  S. I. Braginskii.  Transport Processes in a Plasma , 1965 .

[12]  Martine Baelmans,et al.  Development and assessment of 2D fluid neutral models that include atomic databases and a microscopic reflection model , 2017 .

[13]  Giacomo Dimarco,et al.  Asymptotic-Preserving Monte Carlo Methods for Transport Equations in the Diffusive Limit , 2017, SIAM J. Sci. Comput..

[14]  Graham V. Candler,et al.  A hybrid continuum/particle approach for modeling subsonic, rarefied gas flows , 2004 .

[15]  Stéphane Cordier,et al.  An asymptotic preserving scheme for hydrodynamics radiative transfer models , 2007, Numerische Mathematik.

[16]  A. Klar.  An Asymptotic-Induced Scheme for Nonstationary Transport Equations in the Diffusive Limit , 1998 .

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

[18]  Jun S. Liu,et al.  Monte Carlo strategies in scientific computing , 2001 .

[19]  G. Samaey,et al.  Hybrid fluid-kinetic model for neutral particles in the plasma edge , 2019, Nuclear Materials and Energy.

[20]  Charles F. F. Karney,et al.  Modeling of neutral plasma in a divertor in the fluid‐kinetic transition , 1998 .

[21]  Nicolas Crouseilles,et al.  A particle micro-macro decomposition based numerical scheme for collisional kinetic equations in the diffusion scaling , 2017, 1701.05069.

[22]  Giacomo Dimarco,et al.  Hybrid Multiscale Methods II. Kinetic Equations , 2008, Multiscale Model. Simul..

[23]  D. Coster,et al.  Effect of conditions for gas recirculation on divertor operation in ITER , 2007 .

[24]  A. Wootton,et al.  Fluctuations and anomalous transport in tokamaks , 1990 .

[25]  Donald L. Iglehart,et al.  Importance sampling for stochastic simulations , 1989 .

[26]  N. Crouseilles,et al.  An asymptotic preserving scheme based on a micro-macro decomposition for collisional Vlasov equations: diffusion and high-field scaling limits. , 2011 .

[27]  H. Grad.  On the kinetic theory of rarefied gases , 1949 .

[28]  Wojciech Fundamenski.  Power Exhaust in Fusion Plasmas , 2009 .

[29]  P. Stangeby.  Plasma sheath transmission factors for tokamak edge plasmas , 1984 .

[30]  Michael Junk,et al.  Domain of Definition of Levermore's Five-Moment System , 1998 .

[31]  D. Reiter,et al.  The EIRENE and B2-EIRENE Codes , 2005 .

[32]  A. Bobylev,et al.  The Chapman-Enskog and Grad methods for solving the Boltzmann equation , 1982 .

[33]  S. Voskoboynikov,et al.  Simulation of tokamak edge plasma including self-consistent electric fields , 2001 .

[34]  Paul H Rutherford,et al.  Introduction to Plasma Physics , 1995 .

[35]  Dmitriy Y. Anistratov,et al.  Computational transport methodology based on decomposition of a problem domain into transport and diffusive subdomains , 2012, J. Comput. Phys..

[36]  G. Ciraolo,et al.  Two-phases hybrid model for neutrals , 2019, Nuclear Materials and Energy.