Fast elliptic solvers in cylindrical coordinates and the Coulomb collision operator

Abstract

In this paper, we describe a new class of fast solvers for separable elliptic partial differential equations in cylindrical coordinates $(r,\theta ,z)$ with free-space radiation conditions. By combining integral equation methods in the radial variable $r$ with Fourier methods in $\theta$ and $z$, we show that high-order accuracy can be achieved in both the governing potential and its derivatives. A weak singularity arises in the Fourier transform with respect to $z$ that is handled with special purpose quadratures. We show how these solvers can be applied to the evaluation of the Coulomb collision operator in kinetic models of ionized gases.

20 pages, 5 figures

References in corpus (3)

• Fokker-Planck and Quasilinear Codes
• Linearized model Fokker-Planck collision operators for gyrokinetic simulations. I. Theory
• Differential Form of the Collision Integral for a Relativistic Plasma

Cited by in corpus (8)

• New velocity-space discretization for continuum kinetic calculations and Fokker-Planck collisions
• NORSE: A solver for the relativistic non-linear Fokker-Planck equation for electrons in a homogeneous plasma
• Landau Collision Integral Solver with Adaptive Mesh Refinement on Emerging Architectures
• An entropy based thermalization scheme for hybrid simulations of Coulomb collisions
• Differential formulation of the gyrokinetic Landau operator
• The Gaussian Radial Basis Function Method for Plasma Kinetic Theory
• A Spectral Transform Method for Singular Sturm-Liouville Problems with Applications to Energy Diffusion in Plasma Physics
• Kinetic Simulation of Collisional Magnetized Plasmas with Semi-Implicit Time Integration