On parallel PIC versatility and the structure of parallel PIC approaches

Abstract

We discuss insights gained from an effort to develop a versatile parallel PIC code (PICARD – PIC with arbitrary rebalancing and decompositions). Because of the highly dynamic nature of PIC, it is sometimes necessary to tailor decomposition strategies on a problem-by-problem basis instead of solely on the basis of the algorithm. The dynamics of a problem may dictate the particular strategy to be used. Thus, finding a general operational method for parallel PIC applications is important. We find approaches to partitioning the mesh and particle domains which provide some trade-off opportunities, but the relationship enforced between the partition layouts in addition to the way one chooses to manage the problem's dynamics defines categories of strategies which are independent of the partitioning technique used. The overall strategy enforced is as important as the individual techniques applied to partition the domains. Three strategy groups arise from the analysis and form the basis of widely disparate parallel PIC operational methods. These distinct solutions are obtained from rather simple parametric choices. We introduce the concept of rubber boundaries as a means to bring a compromise between what are commonly referred to as Lagrangian and Eulerian particle mappings. We show several parallel PIC strategies that are as yet unexplored. © 1997 John Wiley & Sons, Ltd.