Temporal parallelization of edge plasma simulations using the parareal algorithm and the SOLPS code

Abstract It is shown that numerical modelling of edge plasma physics may be successfully parallelized in time. The parareal algorithm has been employed for this purpose and the SOLPS code package coupling the B2.5 finite-volume fluid plasma solver with the kinetic Monte-Carlo neutral code Eirene has been used as a test bed. The complex dynamics of the plasma and neutrals in the scrape-off layer (SOL) region makes this a unique application. It is demonstrated that a significant computational gain (more than an order of magnitude) may be obtained with this technique. The use of the IPS framework for event-based parareal implementation optimizes resource utilization and has been shown to significantly contribute to the computational gain.

[1]  A. Loarte,et al.  Status and physics basis of the ITER divertor , 2009 .

[2]  Ronald H. Cohen,et al.  Scrape‐Off Layer Turbulence Theory and Simulations , 1997 .

[3]  Raúl Sánchez,et al.  Parallelization in time of numerical simulations of fully-developed plasma turbulence using the parareal algorithm , 2010, J. Comput. Phys..

[4]  D. P. Coster.  Reduced Physics Models in SOLPS for Reactor Scoping Studies , 2016 .

[5]  O E Garcia,et al.  Bursting and large-scale intermittency in turbulent convection with differential rotation. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[6]  Xavier Bonnin,et al.  The two-mesh grid refinement method for the B2 code , 2002 .

[7]  Charbel Farhat,et al.  Time‐decomposed parallel time‐integrators: theory and feasibility studies for fluid, structure, and fluid–structure applications , 2003 .

[8]  Andrew J. Christlieb,et al.  A Parallel Space-Time Algorithm , 2012, SIAM J. Sci. Comput..

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

[10]  Y. Maday,et al.  A “Parareal” Time Discretization for Non-Linear PDE’s with Application to the Pricing of an American Put , 2002 .

[11]  A. Kukushkin,et al.  Divertor modelling and extrapolation to reactor conditions , 2002 .

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

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

[14]  Einar M. Rønquist,et al.  Stability of the Parareal Algorithm , 2005 .

[15]  F. Sardei,et al.  3D Edge Modeling and Island Divertor Physics , 2004 .

[16]  Rolf Krause,et al.  A space-time parallel solver for the three-dimensional heat equation , 2013, PARCO.

[17]  Mohamed Kamel Riahi,et al.  Parareal in time 3D numerical solver for the LWR Benchmark neutron diffusion transient model , 2014, J. Comput. Phys..

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

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

[20]  Eric Aubanel.  Scheduling of tasks in the parareal algorithm , 2011, Parallel Comput..

[21]  A. Kirschner,et al.  Simulation of the plasma-wall interaction in a tokamak with the Monte Carlo code ERO-TEXTOR , 2000 .

[22]  Xavier Bonnin,et al.  Electric Fields and Currents in an Island Divertor Configuration , 2001 .

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

[24]  Samantha S. Foley,et al.  A dependency-driven formulation of parareal: parallel-in-time solution of PDEs as a many-task application , 2011, MTAGS '11.

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

[26]  Shuichi Takamura,et al.  Chapter 4: Power and particle control , 2007 .

[27]  Y Maday,et al.  Parallel-in-time molecular-dynamics simulations. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[28]  Taina Kurki-Suonio,et al.  Paritcle simulation of the neoclassical plasmas , 2001 .

[29]  Raúl Sánchez,et al.  Mechanisms for the convergence of time-parallelized, parareal turbulent plasma simulations , 2012, J. Comput. Phys..

[30]  Ralf Schneider,et al.  B2-EIRENE code modelling of an island divertor , 1999 .

[31]  Yvon Maday,et al.  A Parareal in Time Semi-implicit Approximation of the Navier-Stokes Equations , 2005 .

[32]  T. W. Petrie,et al.  ELM energy scaling in DIII-D , 2002 .

[33]  Guido Ciraolo,et al.  Applications of SOLEDGE-2D code to complex SOL configurations and analysis of Mach probe measurements , 2011 .

[34]  Raúl Sánchez,et al.  Event-based parareal: A data-flow based implementation of parareal , 2012, J. Comput. Phys..

[35]  Martin J. Gander,et al.  Nonlinear Convergence Analysis for the Parareal Algorithm , 2008 .

[36]  P. C. Stangeby,et al.  Calculation of observable quantities using a divertor impurity interpretive code, DIVIMP , 1992 .

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

[38]  Martine Baelmans,et al.  Accuracy and Convergence of Coupled Finite-Volume / Monte-Carlo Codes for Plasma Edge Simulations , 2016 .

[39]  D. Coster,et al.  Benchmarking of a 1D scrape-off layer code SOLF1D with SOLPS and its use in modelling long-legged divertors , 2013, 1305.3144.

[40]  Michael L. Minion,et al.  TOWARD AN EFFICIENT PARALLEL IN TIME METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS , 2012 .

[41]  Guillaume Bal,et al.  PARALLELIZATION IN TIME OF ( STOCHASTIC ) ORDINARY DIFFERENTIAL EQUATIONS , 2006 .

[42]  Michael L. Minion,et al.  A HYBRID PARAREAL SPECTRAL DEFERRED CORRECTIONS METHOD , 2010 .

[43]  A. S. Kukushkin,et al.  Finalizing the ITER divertor design: The key role of SOLPS modeling , 2011 .

[44]  Magne S. Espedal,et al.  A Convergent Algorithm for Time Parallelization Applied to Reservoir Simulation , 2005 .

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