From arteriographies to computational flow in saccular aneurisms: the INRIA experience

Saccular aneurisms illustrate usefulness and possible techniques of image-based modeling of flow in diseased vessels. Aneurism flow is investigated in order to estimate the rupture risk, assuming that the pressure is the major factor and that high-pressure zones are correlated to within-wall strong-stress concentrations. Computational flow is also aimed at providing additional arguments for the treatment strategy. Angiographies of aneurismal vessels of large and medium size are processed to provide three-dimensional reconstruction of the vessel region of interest. Different reconstruction techniques are used for a side and a terminal aneurisms. Reconstruction techniques may lead to different geometries especially with poor input data. The associated facetisation is improved to get a computation-adapted surface triangulation, after a treatment of vessel ends and mesh adaptation. Once the volumic mesh is obtained, the pulsatile flow of an incompressible Newtonian blood is computed using in vivo non-invasive flowmetry and the finite element method. High pressure zones are observed in the aneurism cavity. The pressure magnitude in the aneurism, the location and the size of high pressure zones depend mainly on the aneurism implantation on the vessel wall and its orientation with respect to the blood flux in the upstream vessel. The stronger the blood impacts on the aneurismal wall the higher the pressure. The state of the aneurism neck, where a high-pressure zone can occur, and the location of the aneurism, with an easy access or not, give arguments for the choice between coiling and surgical clipping. Mesh size and 3D reconstruction procedure affect the numerical results. Helpful qualitative data are provided rather than accurate quantitative results in the context of multimodeling.

[1]  H. Borouchaki,et al.  Geometric surface mesh optimization , 1998 .

[2]  Frédéric Hecht,et al.  NSP1B3 : un logiciel pour resoudre les equations de Navier Stokes incompressible 3D , 1991 .

[3]  O. Pironneau.  On the transport-diffusion algorithm and its applications to the Navier-Stokes equations , 1982 .

[4]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[5]  Nicholas Ayache,et al.  Model Based Multiscale Detection and Reconstruction of 3D Vessels , 1998 .

[6]  Katsushi Ikeuchi,et al.  Shape representation and image segmentation using deformable surfaces , 1992, Image Vis. Comput..

[7]  A Fenster,et al.  Accuracy and variability assessment of a semiautomatic technique for segmentation of the carotid arteries from three-dimensional ultrasound images. , 2000, Medical physics.

[8]  Jean-Daniel Boissonnat,et al.  Shape reconstruction from planar cross sections , 1988, Comput. Vis. Graph. Image Process..

[9]  R. Glowinski,et al.  Numerical Methods for Nonlinear Variational Problems , 1985 .

[10]  J. W. Eastwood,et al.  Springer series in computational physics Editors: H. Cabannes, M. Holt, H.B. Keller, J. Killeen and S.A. Orszag , 1984 .

[11]  P. George.  Improvements on Delaunay-based three-dimensional automatic mesh generator , 1997 .

[12]  Jean-Daniel Boissonnat,et al.  Smooth surface reconstruction via natural neighbour interpolation of distance functions , 2000, SCG '00.

[13]  Grégoire Malandain,et al.  Modélisation numérique de l'écoulement dans les anévrismes cérébraux , 1998 .

[14]  J S Milner,et al.  Rapid three-dimensional segmentation of the carotid bifurcation from serial MR images. , 2000, Journal of biomechanical engineering.

[15]  Hanif M. Ladak,et al.  Rapid 3D Segmentation of the Carotid Bifurcation from Serial MR Images , 2000 .

[16]  Jean-Daniel Boissonnat,et al.  Three-dimensional reconstruction of complex shapes based on the Delaunay triangulation , 1993, Electronic Imaging.

[17]  Pascal Frey,et al.  YAMS A fully Automatic Adaptive Isotropic Surface Remeshing Procedure , 2001 .

[18]  William E. Lorensen,et al.  Marching cubes: A high resolution 3D surface construction algorithm , 1987, SIGGRAPH.

[19]  B. Rutt,et al.  Hemodynamics of human carotid artery bifurcations: computational studies with models reconstructed from magnetic resonance imaging of normal subjects. , 1998, Journal of vascular surgery.

[20]  M. Fortin,et al.  A stable finite element for the stokes equations , 1984 .

[21]  Jean-Frédéric Gerbeau,et al.  A Quasi-Newton Algorithm Based on a Reduced Model for Fluid-Structure Interaction Problems in Blood Flows , 2003 .

[22]  Mariette Yvinec,et al.  Algorithmic geometry , 1998 .

[23]  Rainald Löhner,et al.  From Medical Images to CFD Meshes , 1999, IMR.

[24]  Rainald Löhner,et al.  From medical images to anatomically accurate finite element grids , 2001 .

[25]  R. I. Issa,et al.  A pulsatile developing flow in a bend , 1992 .

[26]  James A. Sethian,et al.  Image Processing: Flows under Min/Max Curvature and Mean Curvature , 1996, CVGIP Graph. Model. Image Process..

[27]  Xavier Descombes,et al.  Estimating Gaussian Markov random field parameters in a nonstationary framework: application to remote sensing imaging , 1999, IEEE Trans. Image Process..

[28]  Mark de Berg,et al.  Computational geometry: algorithms and applications , 1997 .

[29]  James S. Duncan,et al.  Medical Image Analysis , 1999, IEEE Pulse.

[30]  Houman Borouchaki,et al.  Surface meshing using a geometric error estimate , 2003 .

[31]  Gabriel Taubin,et al.  Curve and surface smoothing without shrinkage , 1995, Proceedings of IEEE International Conference on Computer Vision.