A comparison of source term estimators in coupled finite-volume/Monte-Carlo methods with applications to plasma edge simulations in nuclear fusion

In many applications, such as plasma edge simulation of a nuclear fusion reactor, a coupled PDE/kinetic description is required, which is usually solved with a coupled finite-volume/Monte-Carlo method. Different procedures have been proposed to estimate the source terms in the finite volume part that appear from the Monte Carlo part of the simulation. In this paper, we present a systematic comparison of the variance and computational cost of a coherent set of such estimation procedures. We compare the different estimation procedures for mass in a simplified forward-backward scattering model problem, where an analytical comparison is possible, and for mass and momentum in a model problem with realistic scattering. Our results reveal a non-trivial dependence of the optimal choice of estimator on the model parameters and show that different estimation procedures prevail for different quantities of interest.

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

[2]  I. Lux.  Systematic study of some standard variance reduction techniques , 1978 .

[3]  I. Lux.  Monte Carlo Particle Transport Methods: Neutron and Photon Calculations , 1991 .

[4]  Thomas E. Booth,et al.  Verification of the history-score moment equations for weight-window variance reduction , 2011 .

[5]  P. K. Sarkar,et al.  Prediction of statistical error and optimization of biased Monte Carlo transport calculations. [Integral equations] , 1979 .

[6]  A. Kukushkin,et al.  Physics of ultimate detachment of a tokamak divertor plasma , 2017, Journal of Plasma Physics.

[7]  F. Sardei,et al.  Transport in island divertors: physics, 3D modelling and comparison to first experiments on W7-AS* , 2002 .

[8]  Variance Reduction Under Exponential and Scattering Angle Biasing: An Analytic Approach , 1986 .

[9]  G. Rybicki.  Radiative transfer , 2019, Climate Change and Terrestrial Ecosystem Modeling.

[10]  D. Harting,et al.  Preliminary Comparison of the Conventional and Quasi‐Snowflake Divertor Configurations with the 2D Code EDGE2D/EIRENE in the FAST Tokamak , 2014 .

[11]  H. Kahn.  APPLICATIONS OF MONTE CARLO , 1954 .

[12]  D. B. MacMillan.  Comparison of Statistical Estimators for Neutron Monte Carlo Calculations , 1966 .

[13]  D. N. Hill,et al.  UEDGE and DEGAS modeling of the DIII-D scrape-off layer plasma. Revision 1 , 1994 .

[14]  P. Stangeby.  The Plasma Boundary of Magnetic Fusion Devices , 2000 .

[15]  Daren P. Stotler,et al.  Neutral Gas Transport Modeling with DEGAS 2 , 1994 .

[16]  E.D. Denman,et al.  An introduction to invariant imbedding , 1977, Proceedings of the IEEE.

[17]  Jerome Spanier,et al.  An Analytic Approach to Variance Reduction , 1970 .

[18]  Giovanni Samaey,et al.  Simulating individual-based models of bacterial chemotaxis with asymptotic variance reduction , 2011, 1111.5321.

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

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

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