Abstract
In this paper, we describe a new class of fast solvers for separable elliptic partial differential equations in cylindrical coordinates (r, θ, z) with free-space radiation conditions. By combining integral equation methods in the radial variable r with Fourier methods in θ 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.
Original language | English (US) |
---|---|
Pages (from-to) | 7840-7852 |
Number of pages | 13 |
Journal | Journal of Computational Physics |
Volume | 230 |
Issue number | 21 |
DOIs | |
State | Published - Sep 1 2011 |
Keywords
- Biharmonic equation
- Collision operator
- Fast solvers
- Plasma physics
- Poisson equation
ASJC Scopus subject areas
- Numerical Analysis
- Modeling and Simulation
- Physics and Astronomy (miscellaneous)
- General Physics and Astronomy
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics