Parallel sparse direct solvers for Poisson's equation in streamer discharges

textabstractThe aim of this paper is to examine whether a hybrid approach of parallel computing, a combination of the message passing model (MPI) with the threads model (OpenMP) can deliver good performance in streamer discharge simulations. Since one of the bottlenecks of almost all streamer models is the solution of Poisson's equation, we focused on several direct solvers, which can solve large sparse systems in parallel. For this purpose, our basic thought was to concentrate on 'easy to get' performance improvements, or, without rewriting of the code. We have investigated in PARDISO, a shared memory solver, and CLUSTER SPARSE SOLVER and MUMPS, which both can apply hybrid parallelism; the latter two solvers can be called from a single core and do not require minor awareness of MPI. We show their performance for solving two- and three-dimensional Poisson's equations on the Dutch national supercomputer, called Cartesius. A runtime study of a code developed for streamer propagation nearby a dielectric rod is included. We discuss various issues that appear to be critical in a mixed MPI-OpenMP environment.

[1]  Konstantin Arturov,et al.  Asynchronous Approach to Memory Management in Sparse Multifrontal Methods on Multiprocessors , 2013 .

[2]  Patrick Amestoy,et al.  Hybrid scheduling for the parallel solution of linear systems , 2006, Parallel Comput..

[3]  Wolfgang Fichtner,et al.  Efficient Sparse LU Factorization with Left-Right Looking Strategy on Shared Memory Multiprocessors , 2000 .

[4]  Anna Dubinova,et al.  Modeling of streamer discharges near dielectrics , 2016 .

[5]  William Gropp,et al.  Mpi---the complete reference: volume 1 , 1998 .

[6]  C. Montijn.  Evolution of negative streamers in nitrogen : a numerical investigation on adaptive grids , 2005 .

[7]  Olaf Schenk,et al.  Fast Methods for Computing Selected Elements of the Green's Function in Massively Parallel Nanoelectronic Device Simulations , 2013, Euro-Par.

[8]  Olaf Schenk,et al.  Solving unsymmetric sparse systems of linear equations with PARDISO , 2004, Future Gener. Comput. Syst..

[9]  Chi-Fang Li,et al.  Joining particle and fluid aspects in streamer simulations , 2005 .

[10]  Olaf Schenk,et al.  Scalable parallel sparse LU factorization methods on shared memory multiprocessors , 2000 .

[11]  J. Verwer,et al.  Numerical solution of time-dependent advection-diffusion-reaction equations , 2003 .

[12]  Jean-Yves L'Excellent,et al.  Some Experiments and Issues to Exploit Multicore Parallelism in a Distributed-Memory Parallel Sparse Direct Solver , 2010 .

[13]  John K. Reid,et al.  The Multifrontal Solution of Indefinite Sparse Symmetric Linear , 1983, TOMS.

[14]  Jannis Teunissen,et al.  3D simulations and analysis of pulsed discharges , 2015 .

[15]  Jannis Teunissen,et al.  Afivo: A framework for quadtree/octree AMR with shared-memory parallelization and geometric multigrid methods , 2017, Comput. Phys. Commun..