Search

We provide the 21N-moment multi-temperature collision coefficients for the Boltzmann collision operator using the Sonine-Hermite polynomial ansatz in the style of Zhdanov et al. First, we outline the general derivation method. Then, we provide the collision coefficients in the most general form in terms of the Chapman-Cowling integrals for any potential of interaction for the 21N-moment approximation. Then, we provide the the collision coefficients for the specific case of the approximated Coulomb potential cross-section with Debye cutoff. These coefficients help in order to find easy implementation in current SOL/edge fluid packages which currently implement the 21N-moment single-temperature coefficients. This provides the completion of a missing link in the current literature., Comment: 19 pages, 0 figures. Accepted by Contributions to Plasma Physics (Dec 2021). arXiv admin note: substantial text overlap with arXiv:2103.02455

The derivation of the multi-temperature generalized Zhdanov closure is provided starting from the most general form of the left hand side of the moment averaged kinetic equation with the Sonine-Hermite polynomial ansatz for an arbitrary number of moments. The process of arriving at the reduced higher-order moment equations, with its assumptions and approximations, is explicitly outlined. The generalized multi-species, multi-temperature coefficients from the authors' previous article are used to compute values of higher order moments such as heat flux in terms of the lower order moments. Transport coefficients and the friction and thermal forces for magnetic confinement fusion relevant cases with the generalized coefficients are compared to the scheme with the single-temperature coefficients previously provided by Zhdanov et al. It is found that the 21N-moment multi-temperature coefficients are adequate for most cases relevant to fusion. Furthermore, the 21N-moment scheme is also tested against the trace approximation to determine the range of validity of the trace approximation with respect to fusion relevant plasmas. Possible refinements to the closure scheme are illustrated as well, in order to account for quantities which might be significant in certain schemes such as the drift approximation., Comment: Submitted to Plasma Physics and Controlled Fusion, July 15, 2021

Grad's method is used on the linearized Boltzmann collision operator to derive the most general expressions for the collision coefficients for a multi-component, multi-temperature plasma up to rank-2. In doing so, the collision coefficients then get expressed as series sum of pure coefficients of temperature and mass ratios multiplied by the cross-section dependent Chapman-Cowling integrals. These collisional coefficients are compared to previously obtained coefficients by Zhdanov et al [Zhdanov V.M., Transport processes in multi-component plasma, Taylor and Francis (2002)] for 13N-moment multi-temperature scheme. First, the differences in coefficients are compared directly, and then the differences in first approximation to viscosity and friction force are compared. For the 13N-moment multi-temperature coefficients, it is found that they behave reasonably similarly for small temperature differences, but display substantial differences in the coefficients when the temperature differences are high, both for the coefficients and for viscosity and friction force values. Furthermore, the obtained coefficients are compared to the 21N-moment single-temperature approximation provided by Zhdanov et al, and it is seen that the differences are higher than the 13N-moment multi-temperature coefficients, and have substantial differences even in the vicinity of equal temperatures, especially for the viscosity and friction force calculations., Comment: Corrected some typos from the previous version. Accepted by Plasma Physics and Controlled Fusion (April 9,2021)

In order to improve the presently used ad hoc flux limiter treatment of parallel heat flux transport in edge plasma fluid codes we consider here a generalized version of the Fourier law implementing a non-local kernel for the heat flux computation. The Bohm boundary condition at the wall is recovered introducing a volumetric loss term representing the contribution of suprathermal particles to the energy out flux. As expected, this contribution is negligible in the strongly collisional regime while it becomes more and more dominant for marginally and low collisional regimes. In the second part of the paper, we consider a kinetic approach where collisions are considered using the Multi-Particle-Collision (MPC) algorithm. Kinetic simulation results at medium and low collisionality are also reported., Comment: 13 Pages, 6 figures, Accepted for publication in Contributions to Plasma Physics