Novel approaches for the reliable and efficient numerical evaluation of Landau-type operators
For researchers simulating plasma physics or kinetic equations, this work provides more efficient numerical methods for Landau-type operators, though it is an incremental improvement over existing Fourier spectral methods.
This work introduces and compares novel numerical approaches for the reliable and efficient evaluation of Landau-type collision operators, focusing on three-dimensional Coulomb interactions. The methods achieve significant computational savings by transferring a large portion of work to precomputations independent of the density function.
Numerical approximations of Landau-type operators represent fundamental components of time integration methods for demanding problems such as inhomogeneous Vlasov-Landau-type equations. Substantial computational issues arise from the treatment of the physically most relevant three-dimensional case with Coulomb-type interaction. This work is concerned with the introduction and numerical comparison of novel approaches for the reliable and efficient evaluation of Landau-type collision operators, where the focus is on the treatment of integral operators involving general singular kernels. In the spirit of collocation, common tools are the identification of fundamental integrals, series expansions of the integral kernel and the density function on the main part of the velocity domain, and interpolation as well as quadrature approximation nearby the singularity of the kernel. Focusing on the favourable choice of the Fourier spectral method, their practical implementation uses the reduction to basic integrals, fast Fourier techniques, and summations along certain directions. Moreover, an important observation is that a significant percentage of the overall computational effort can be transferred to precomputations which are independent of the density function. For the purpose of exposition and numerical validation, the cases of constant, regular, and singular integral kernels are distinguished, and the procedure is adapted accordingly to the increasing complexity of the problem.