A versatile incompressible Navier-Stokes solver for blood flow application

M. Garbey, F. Pacull

Research output: Contribution to journalArticlepeer-review

9 Scopus citations


We present in this paper an integrated approach to compute quickly an incompressible Navier-Stokes (NS) flow in a section of a large blood vessel using medical imaging data. The goal is essentially to provide a first-order approximation of some main quantities of interest in cardiovascular disease: the shear stress and the pressure on the wall. The NS solver relies on the L2 penalty approach pioneered by Caltagirone and co-workers and combines nicely with a level set method based on the Mumford-Shah energy model. Simulations on stenosis cases based on angiogram are run in parallel with MatlabMPI on a shared-memory machine. While MatlabMPI communications are based on the load and save functions of Matlab and have high latency indeed, we show that our Aitken-Schwarz domain decomposition algorithm provides a good parallel efficiency and scalability of the NS code.

Original languageEnglish (US)
Pages (from-to)473-496
Number of pages24
JournalInternational Journal for Numerical Methods in Fluids
Issue number5
StatePublished - Jun 20 2007


  • Blood flow
  • Image analysis
  • Incompressible flow
  • Parallel algorithm

ASJC Scopus subject areas

  • Computational Mechanics
  • Mechanics of Materials
  • Mechanical Engineering
  • Computer Science Applications
  • Applied Mathematics


Dive into the research topics of 'A versatile incompressible Navier-Stokes solver for blood flow application'. Together they form a unique fingerprint.

Cite this