We present an efficient and accurate immersed boundary (IB) finite element (FE) method for numerically solving the incompressible Navier–Stokes equations. Particular emphasis is given to internal flows with complex geometries (blood flow in vasculature system). IB methods are computationally costly flows, mainly due large percentage of grid points that lie outside domain. In this study, we appl...