Open source software package for simulation of viscous incompressible flows using the particle finite element method PFEM-2

A large number of complex real life fluid dynamics problems lead to the need of simulation of convection dominated flows (multiphase flows, high Reynolds number flows, fluid-structure interaction problems). As opposed to conventional approaches, specialized methods (such as PFEM-2) can be more effective for their solution. The PFEM-2 method is based on a splitting procedure which implies that convection is simulated by transport of Lagrangian particles with other processes simulated on a Eulerian mesh. Practical use of such promising contemporary methods is directly connected to the availability of their software implementations. In the case of PFEM-2 the only implementation is publicly available so far, which is developed within the Kratos multiphysics framework by the authors of this method. A new software package for simulation of viscous incompressible flows using PFEM-2 is developed by the author of this paper. The implementation described here gives more accurate results for commonly used test problems. This software package is based only on the open source software such as deal.II, p4est and Trilinos libraries, CMake package, UNV and VTK file formats. We develop parallel implementations of the algorithms using MPI and NVIDIA CUDA technologies with efficiency of parallelization studied for high performance computational systems of various architecture. The source code and test cases described in this paper are publicly available at the github repository.