Parallel Implementation and Performance of an Immersed Boundary Method

An implementation for parallelization of an immersed boundary method is presented using the standard MPI communication. The immersed boundary method in the present work follows that in Zhang & Zheng [1]. The Poisson problem is solved with assistance of a scientific parallel computational library PETSc. A staggered grid mesh is employed to reduce the discretization error. The 2^nd order accuracy is obtained eventually. This method is validated by a case with a 2D uniform flow over a stationary or oscillating circular cylinder and a 3D uniform flow over a stationary sphere. Comparisons with previous studies are presented. Also, performances of parallelization, including parallelism and scalability, are tested.