skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A fast parallel Poisson solver on irregular domains applied to beam dynamics simulations

Abstract

We discuss the scalable parallel solution of the Poisson equation within a Particle-In-Cell (PIC) code for the simulation of electron beams in particle accelerators of irregular shape. The problem is discretized by Finite Differences. Depending on the treatment of the Dirichlet boundary the resulting system of equations is symmetric or 'mildly' nonsymmetric positive definite. In all cases, the system is solved by the preconditioned conjugate gradient algorithm with smoothed aggregation (SA) based algebraic multigrid (AMG) preconditioning. We investigate variants of the implementation of SA-AMG that lead to considerable improvements in the execution times. We demonstrate good scalability of the solver on distributed memory parallel processor with up to 2048 processors. We also compare our iterative solver with an FFT-based solver that is more commonly used for applications in beam dynamics.

Authors:
 [1]
  1. Paul Scherrer Institut, CH-5234 Villigen (Switzerland)
Publication Date:
OSTI Identifier:
21333914
Resource Type:
Journal Article
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 229; Journal Issue: 12; Other Information: DOI: 10.1016/j.jcp.2010.02.022; PII: S0021-9991(10)00102-6; Copyright (c) 2010 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA); Journal ID: ISSN 0021-9991
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; AGGLOMERATION; ALGORITHMS; BEAM DYNAMICS; COMPARATIVE EVALUATIONS; DIRICHLET PROBLEM; ELECTRON BEAMS; ITERATIVE METHODS; MATHEMATICAL SOLUTIONS; POISSON EQUATION; SIMULATION; SPACE CHARGE

Citation Formats

Adelmann, A., Arbenz, P., Ineichen, Y, and ETH Zuerich, Chair of Computational Science, Universitaetsstrasse 6, CH-8092 Zuerich. A fast parallel Poisson solver on irregular domains applied to beam dynamics simulations. United States: N. p., 2010. Web. doi:10.1016/j.jcp.2010.02.022.
Adelmann, A., Arbenz, P., Ineichen, Y, & ETH Zuerich, Chair of Computational Science, Universitaetsstrasse 6, CH-8092 Zuerich. A fast parallel Poisson solver on irregular domains applied to beam dynamics simulations. United States. doi:10.1016/j.jcp.2010.02.022.
Adelmann, A., Arbenz, P., Ineichen, Y, and ETH Zuerich, Chair of Computational Science, Universitaetsstrasse 6, CH-8092 Zuerich. Sun . "A fast parallel Poisson solver on irregular domains applied to beam dynamics simulations". United States. doi:10.1016/j.jcp.2010.02.022.
@article{osti_21333914,
title = {A fast parallel Poisson solver on irregular domains applied to beam dynamics simulations},
author = {Adelmann, A. and Arbenz, P. and Ineichen, Y and ETH Zuerich, Chair of Computational Science, Universitaetsstrasse 6, CH-8092 Zuerich},
abstractNote = {We discuss the scalable parallel solution of the Poisson equation within a Particle-In-Cell (PIC) code for the simulation of electron beams in particle accelerators of irregular shape. The problem is discretized by Finite Differences. Depending on the treatment of the Dirichlet boundary the resulting system of equations is symmetric or 'mildly' nonsymmetric positive definite. In all cases, the system is solved by the preconditioned conjugate gradient algorithm with smoothed aggregation (SA) based algebraic multigrid (AMG) preconditioning. We investigate variants of the implementation of SA-AMG that lead to considerable improvements in the execution times. We demonstrate good scalability of the solver on distributed memory parallel processor with up to 2048 processors. We also compare our iterative solver with an FFT-based solver that is more commonly used for applications in beam dynamics.},
doi = {10.1016/j.jcp.2010.02.022},
journal = {Journal of Computational Physics},
issn = {0021-9991},
number = 12,
volume = 229,
place = {United States},
year = {2010},
month = {6}
}