A parallel 3D Poisson solver for space charge simulation in cylindrical coordinates

This paper presents the development of a parallel three-dimensional Poisson solver in cylindrical coordinate system for the electrostatic potential of a charged particle beam in a circular tube. The Poisson solver uses Fourier expansions in the longitudinal and azimuthal directions, and Spectral Ele...

Full description

Saved in:
Bibliographic Details
Published in:Computer physics communications 2008-02, Vol.178 (4), p.290-300
Main Authors: Xu, J., Ostroumov, P.N., Nolen, J.
Format: Article
Language:English
Subjects:
ACCELERATORS
APERTURES
BEAM DYNAMICS
BENCHMARKS
BOUNDARY CONDITIONS
CARTESIAN COORDINATES
CHARGED PARTICLES
Domain decomposition
ELECTROSTATICS
Fourier expansion
LINEAR ACCELERATORS
PARTICLE ACCELERATORS
Poisson equation
PROTONS
SIMULATION
SPACE CHARGE
Spectral element method
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:This paper presents the development of a parallel three-dimensional Poisson solver in cylindrical coordinate system for the electrostatic potential of a charged particle beam in a circular tube. The Poisson solver uses Fourier expansions in the longitudinal and azimuthal directions, and Spectral Element discretization in the radial direction. A Dirichlet boundary condition is used on the cylinder wall, a natural boundary condition is used on the cylinder axis and a Dirichlet or periodic boundary condition is used in the longitudinal direction. A parallel 2D domain decomposition was implemented in the ( r , θ ) plane. This solver was incorporated into the parallel code PTRACK for beam dynamics simulations. Detailed benchmark results for the parallel solver and a beam dynamics simulation in a high-intensity proton LINAC are presented. When the transverse beam size is small relative to the aperture of the accelerator line, using the Poisson solver in a Cartesian coordinate system and a Cylindrical coordinate system produced similar results. When the transverse beam size is large or beam center located off-axis, the result from Poisson solver in Cartesian coordinate system is not accurate because different boundary condition used. While using the new solver, we can apply circular boundary condition easily and accurately for beam dynamic simulations in accelerator devices.
ISSN:0010-4655
1879-2944
DOI:10.1016/j.cpc.2007.09.008