High frequency 2D acoustic scattering in distributed memory systems: An OpenMPI-Scalapack implementation.

EDMUNDO LAVIA

Congreso; Escuela Giambiagi; 2017

Facultad de Ciencias Exactas y Naturales. Universidad de Buenos Aires

When the problem of a plane acoustic wave incident on a 2D-obstacle immersed in a homogeneous fluid is considered, a solution can be obtained from an integral collocation method based on Green identities. The numerical fulfillment of that method involves the solution of a matrix system whose size depends quadratically on the wave number k = 2 pi/lambda, where lambda is the wave length. In the high-frequency regime, ka >> 1, being a a characteristic length of the obstacle that produces the acoustic scattering, the size of the associated matrix exceeds the capabilities of current desktop PC. This work presents the results of implementing a numerical solution of 2Dscattering under Neumann boundary conditions (sound hard case) in a computing cluster, using MPI and the ScaLAPACK library which allows for solving the matrix system in a distributed way. This means distributing the allocation in memory of the matrix elements among all the processors in the cluster, no matter in which host they really are located.