Parallel Image-Based Hemodynamic Simulator

With cost effective distributed memory computer systems reaching high performances, it may become feasible in the near future to provide routinely reliable blood flow simulations during angiographic procedures to enhance standard medical imaging techniques. The long term goal of our work is to produce a fast parallel image base Navier-Stokes flow solver for these kinds of clinical procedures. To achieve such performance, our method relies on the proper combination of three techniques that are the L2 penalty method to deal with complex geometry, a level set method to extract the geometry of the arteries, and a fast scalable domain decomposition solver for the parallel processing of the discrete equations. We developed a web application with a database back-end to store both the original medical images and the results of these simulations. This will assist the endovascular surgeon by providing quickly blood flow indicators such as pressure or shear stresses at the wall of the artery, in addition to existing image analysis technique aimed at improving therapy.

[1]  Thomas J. R. Hughes,et al.  Finite Element Modeling of Three-Dimensional Pulsatile Flow in the Abdominal Aorta: Relevance to Atherosclerosis , 2004, Annals of Biomedical Engineering.

[2]  Wei Shyy,et al.  Moving boundaries in micro-scale biofluid dynamics , 2001 .

[3]  W. Shyy,et al.  Regular Article: An Accurate Cartesian Grid Method for Viscous Incompressible Flows with Complex Immersed Boundaries , 1999 .

[4]  Michael M. Resch,et al.  Efficient metacomputing of elliptic linear and non-linear problems , 2003, J. Parallel Distributed Comput..

[5]  Marc Garbey,et al.  A versatile incompressible Navier–Stokes solver for blood flow application , 2007 .

[6]  S. Osher,et al.  Geometric Level Set Methods in Imaging, Vision, and Graphics , 2011, Springer New York.

[7]  Tony F. Chan,et al.  Active contours without edges , 2001, IEEE Trans. Image Process..

[8]  Carlos Castro,et al.  Proceedings of the ACM Symposium on Applied Computing , 2003 .

[9]  David A. Steinman,et al.  Path-Dependent Hemodynamics of the Stenosed Carotid Bifurcation , 2003, Annals of Biomedical Engineering.

[10]  Marc Garbey,et al.  Acceleration of the Schwarz Method for Elliptic Problems , 2005, SIAM J. Sci. Comput..

[11]  Marc Garbey,et al.  On some Aitken-like acceleration of the Schwarz method , 2002 .

[12]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[13]  John W. Ward,et al.  Silent victories: the history and practice of public health in twentieth-century America. , 2007 .

[14]  Philippe Angot,et al.  A penalization method to take into account obstacles in incompressible viscous flows , 1999, Numerische Mathematik.

[15]  Kai Schneider,et al.  Numerical simulation of the transient flow behaviour in tube bundles using a volume penalization method , 2005 .

[16]  Marc Garbey,et al.  Fast elliptic solver for incompressible Navier Stokes Flow and Heat Transfer problems on the grid , 2005 .

[17]  Andrew Simpson,et al.  An Approach to the Storage of DICOM files for Grid−Enabled Medical Imaging Databases , 2004 .

[18]  Marc Garbey,et al.  WIMI - web interface to medical information , 2007, SEDE.

[19]  Eric Mestreau,et al.  On The Modeling Of Boundary Conditions For Embedded Schemes , 2003 .

[20]  Andrew C. Simpson,et al.  A relational approach to the capture of DICOM files for Grid-enabled medical imaging databases , 2004, SAC '04.

[21]  Isaac N. Bankman,et al.  Handbook of medical imaging , 2000 .