Development of parallel implementation for the Navier-Stokes equation in doubly connected areas using the fictitious domain method

Abstract

This paper presents a numerical realization of the Navier-Stokes equations in irregular domains using the fictitious domain method with a continuation along with the lowest coefficient. To solve numerous connected issues in irregular regions, the fictitious domain method is broadly used. The advantage of the fictitious domain method is that the problem is solved not in the original complex domain, but in a few other, easier domains. Using the method, computation is done easily for a sufficiently wide class of problems with arbitrary computational domains. The problem is solved using two methods. The primary method is based on the development of a distinct issue in variables of the stream function and the vortex velocity using the pressure uniqueness condition. The second method is to understand the expressed issue by the fictitious domain method with a continuation by lower coefficients. Using the fictitious domain method, a computational algorithm is constructed based on the explicit finite difference schemes. The finite difference scheme is stable and has high computational accuracy and it gives the possibility to parallelize. Temperature distributions and stream functions are presented as numerical results. A parallel algorithm has been developed using Open Multi-Processing (hereinafter OpenMP) and Message Passing Interface (hereinafter MPI) technologies. Within the parallel approach, we used OpenMP technology for parallel calculation of vorticity and stream work, and for calculating temperature we applied MPI technology. The performance analysis on our parallel code shows favorable strong and weak scalability. The test results show that the code running in the parallel approach gives the expected results by comparing our results with those obtained while running the same simulation on the central processing unit (CPU)