Abstract
The objective of this paper is to present a fast parallel incompressible Navier Stokes (NS) simulator that has enough versatility to apply to various endovascular situations with ease in medical image integration. Our long term goal is to bring NS simulation to surgeons as an additional routine image modality that can complement existing medical imaging channels. Given a medical image describing the geometry of an artery, we are able to simulate the blood flow through an artery close to clinical conditions. To achieve this goal, our interface relies on the combination of the L2 penalty method to deal with complex geometry for solving the NS equations. We present here a high performance portable domain decomposition solver based on a tuning of the fastest linear solver depending on the subdomain and the processor architecture. This paper provides a detailed proof of concept in the two dimensional case. The three dimensional results will be reported in a companion paper.