Collisional transport

Predicting flows of particles, heat, and electric current in magnetized plasmas

In magnetized plasmas, charged particles have complicated trajectories, drifting across the magnetic field lines due to nonuniformity in B (figure 2 here) . Particles also experience collisions, causing them to change from one trajectory to another. At the same time, sources and sinks of particles and heat mean that there are gradients of density and temperature. The interplay of geometrically complicated trajectories with collisions and density/temperature gradients gives rise to complicated non-Maxwellian distribution functions, as shown in figure 1. These phenomena result in flows of particles and heat, with two important consequences. First, these flows set the minimum possible rate at which heat leaks out of the confinement region; any turbulence would result in even higher energy loss. Second, differences between ion and electron flows represent an electric current, and this self-generated current in the plasma modifies the magnetic field. Together, these phenomena are called “collisional transport” or “neoclassical transport.”

Figure 1. The distribution function in a magnetically confined plasma has complicated non-Maxwellian structure with internal boundary layers. The distribution function on the right was computed with the SFINCS code.

Computing these neoclassical effects is numerically challenging. One reason is the high dimensionality of the problem: not only real space but also velocity space must be discretized. Furthermore, high resolution in velocity space is required, since internal boundary layers in the distribution form at the boundaries between different classes of trajectories, as shown in figure 1. To solve this challenging problem, our group has developed the SFINCS code. SFINCS uses modern parallelized numerical libraries (PETSc and MUMPS) to compute the distribution function in stellarator and tokamak plasmas. The code also computes moments of the distribution function to find the fluxes of particles and heat, and the self-generated electric current in the plasma.

SFINCS is being applied to several stellarator experiments, including W7-X. Some papers describing SFINCS simulations of W7-X can be found here, here, and here. One example of experimental validation using SFINCS is shown in figure 2, which shows good agreement between the electric field computed by the code and the electric field measured in W7-X.

Figure 2. Comparison of the radial electric field measured in the W7-X stellarator (two diagnostics shown in red and green) with the prediction of the SFINCS code developed by our group. Figure from Pablant et al (2018).

Besides modeling of existing experiments, our group is pursuing many other applications of the SFINCS code. Graduate student Elizabeth Paul has used SFINCS to predict the torque in ITER, a large tokamak being constructed, caused by imperfections in its rotation symmetry. The collisional flux of impurities is outward for perfect quasisymmetry but inward if quasisymmetry is strongly broken, and PhD student Mike Martin is using SFINCS to study how accurately quasisymmetry must be achieved to avoid impurity accumulation. We are also using SFINCS to predict the self-generated plasma current (“bootstrap current”) as part of the design of future stellarator experiments, including WISTELL.