Development of an efficient parallel Poisson equation solver for the simulation of protoplanetary disk evolution

Abstract

A Poisson equation solver is proposed. The main idea of the solver is taking the features of the problem under study into account. Since the problem is nonstationary, it is possible to solve the Poison equation quickly by means of simple iteration methods. The consideration of the axial symmetry of the disk allows one to construct a domain decomposition with a very small amount of interprocessor communications. In addition, the infinitely thin disk approximation allows one to speed up the computations when using the iteration methods. The paper was prepared on the basis of the authors' report at the International Conference on Parallel Computing Technologies (PaVT-2009; http://agora.guru.ru/pavt).