Modelling pipeline for subject‐specific arterial blood flow—A review

SUMMARY In this paper, a robust and semi-automatic modelling pipeline for blood flow through subject-specific arterial geometries is presented. The framework developed consists of image segmentation, domain discretization (meshing) and fluid dynamics. All the three subtopics of the pipeline are explained using an example of flow through a severely stenosed human carotid artery. In the Introduction, the state-of-the-art of both image segmentation and meshing is presented in some detail, and wherever possible the advantages and disadvantages of the existing methods are analysed. Followed by this, the deformable model used for image segmentation is presented. This model is based upon a geometrical potential force (GPF), which is a function of the image. Both the GPF calculation and level set determination are explained. Following the image segmentation method, a semi-automatic meshing method used in the present study is explained in full detail. All the relevant techniques required to generate a valid domain discretization are presented. These techniques include generating a valid surface mesh, skeletonization, mesh cropping, boundary layer mesh construction and various mesh cosmetic methods that are essential for generating a high-quality domain discretization. After presenting the mesh generation procedure, how to generate flow boundary conditions for both the inlets and outlets of a geometry is explained in detail. This is followed by a brief note on the flow solver, before studying the blood flow through the carotid artery with a severe stenosis. Copyright © 2011 John Wiley & Sons, Ltd.

[1]  Perumal Nithiarasu,et al.  Computational flow studies in a subject‐specific human upper airway using a one‐equation turbulence model. Influence of the nasal cavity , 2011 .

[2]  D. Ku,et al.  Pulsatile flow in the human left coronary artery bifurcation: average conditions. , 1996, Journal of biomechanical engineering.

[3]  Qiang Du,et al.  Centroidal Voronoi Tessellations: Applications and Algorithms , 1999, SIAM Rev..

[4]  D. Giddens,et al.  Oscillatory shear stress stimulates endothelial production of O2- from p47phox-dependent NAD(P)H oxidases, leading to monocyte adhesion. , 2003, The Journal of biological chemistry.

[5]  C F Dewey,et al.  Vascular endothelial cells respond to spatial gradients in fluid shear stress by enhanced activation of transcription factors. , 1999, Arteriosclerosis, thrombosis, and vascular biology.

[6]  D C Barber,et al.  Characterisation of the haemodynamics of the superior mesenteric artery. , 2007, Journal of biomechanics.

[7]  P. Nithiarasu.  A SIMPLE LOCALLY CONSERVATIVE GALERKIN (LCG) FINITE-ELEMENT METHOD FOR TRANSIENT CONSERVATION EQUATIONS , 2004 .

[8]  R. van Loon,et al.  Towards computational modelling of aortic stenosis , 2010 .

[9]  Xiangmin Jiao,et al.  Fluid–structure interactions of the mitral valve and left heart: Comprehensive strategies, past, present and future , 2010, International journal for numerical methods in engineering.

[10]  Toshiaki Akita,et al.  Three-dimensional numerical simulation of blood flow in the aortic arch during cardiopulmonary bypass. , 2008, European journal of cardio-thoracic surgery : official journal of the European Association for Cardio-thoracic Surgery.

[11]  P. Nithiarasu.  An efficient artificial compressibility (AC) scheme based on the characteristic based split (CBS) method for incompressible flows , 2003 .

[12]  T. O'donnell,et al.  Pulsatile flow and atherosclerosis in the human carotid bifurcation: Positive correlation between plaque location and low and oscillating shear stress: Ku DN, Giddens DP, Zarins CK, et al. Arteriosclerosis 1985; 5: 293–302 , 1986 .

[13]  Baba C. Vemuri,et al.  Shape Modeling with Front Propagation: A Level Set Approach , 1995, IEEE Trans. Pattern Anal. Mach. Intell..

[14]  Christophe Geuzaine,et al.  Quality meshing based on STL triangulations for biomedical simulations , 2010 .

[15]  T. Ishikawa,et al.  A three‐dimensional particle simulation of the formation and collapse of a primary thrombus , 2010 .

[16]  P. Nithiarasu.  A Matrix Free Fractional Step Method for Static and Dynamic Incompressible Solid Mechanics , 2006 .

[17]  Yuri Bazilevs,et al.  A fully-coupled fluid-structure interaction simulation of cerebral aneurysms , 2010 .

[18]  Q. Long,et al.  Reconstruction of blood flow patterns in human arteries , 1999, Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine.

[19]  Takami Yamaguchi,et al.  Construction of the CFD model of the aortic arch based on MR images and simulation of the blood flow , 2001, Proceedings International Workshop on Medical Imaging and Augmented Reality.

[20]  Peter L. Choyke,et al.  Vessel surface reconstruction with a tubular deformable model , 2001, IEEE Transactions on Medical Imaging.

[21]  S. Osher,et al.  Algorithms Based on Hamilton-Jacobi Formulations , 1988 .

[22]  Xianghua Xie,et al.  Geometric Potential Force for the Deformable Model , 2009, BMVC.

[23]  Jia Lu,et al.  Fluid–structure interaction methods in biological flows with special emphasis on heart valve dynamics , 2010 .

[24]  S. Sherwin,et al.  The spectral/hp element modelling of steady flow in non‐planar double bends , 2008 .

[25]  L. Antiga,et al.  Computational geometry for patient-specific reconstruction and meshing of blood vessels from MR and CT angiography , 2003, IEEE Transactions on Medical Imaging.

[26]  Alejandro F. Frangi,et al.  COMPUTATIONAL ANALYSIS OF BLOOD FLOW DYNAMICS IN CEREBRAL ANEURYSMS FROM CTA AND 3 D ROTATIONAL ANGIOGRAPHY IMAGE DATA , 2003 .

[27]  M. Yin,et al.  Effects of flow vortex on a chorded mitral valve in the left ventricle , 2010 .

[28]  Xianghua Xie,et al.  Geometrically Induced Force Interaction for Three-Dimensional Deformable Models , 2011, IEEE Transactions on Image Processing.

[29]  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.

[30]  Ian Marshall,et al.  MRI measurement of time‐resolved wall shear stress vectors in a carotid bifurcation model, and comparison with CFD predictions , 2003, Journal of magnetic resonance imaging : JMRI.

[31]  Christophe Geuzaine,et al.  High‐quality surface remeshing using harmonic maps , 2010 .

[32]  Thomas J. R. Hughes,et al.  Patient-Specific Vascular NURBS Modeling for Isogeometric Analysis of Blood Flow , 2007, IMR.

[33]  R. Glowinski,et al.  A fictitious domain method for simulating viscous flow in a constricted elastic tube subject to a uniform external pressure , 2010 .

[34]  Tayfun E. Tezduyar,et al.  Space–time finite element computation of arterial fluid–structure interactions with patient‐specific data , 2010 .

[35]  David A. Steinman,et al.  Image-Based Computational Fluid Dynamics Modeling in Realistic Arterial Geometries , 2002, Annals of Biomedical Engineering.

[36]  Daniel P. Huttenlocher,et al.  Distance Transforms of Sampled Functions , 2012, Theory Comput..

[37]  Anders Heyden,et al.  A fast algorithm for level set-like active contours , 2003, Pattern Recognit. Lett..

[38]  Carl Ollivier-Gooch,et al.  Tetrahedral mesh improvement using swapping and smoothing , 1997 .

[39]  Tzung K Hsiai,et al.  Pulsatile Versus Oscillatory Shear Stress Regulates NADPH Oxidase Subunit Expression: Implication for Native LDL Oxidation , 2003, Circulation research.

[40]  R. Löhner,et al.  Fast numerical solutions of patient‐specific blood flows in 3D arterial systems , 2010, International journal for numerical methods in biomedical engineering.

[41]  J. Womersley.  Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known , 1955, The Journal of physiology.

[42]  Wolfgang A. Wall,et al.  A computational strategy for prestressing patient‐specific biomechanical problems under finite deformation , 2010 .

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

[44]  Jong Beom Ra,et al.  A locally adaptive region growing algorithm for vascular segmentation , 2003, Int. J. Imaging Syst. Technol..

[45]  Jan Vierendeels,et al.  Simulation of fluid–structure interaction with the interface artificial compressibility method , 2010 .

[46]  Adam Huang,et al.  On Concise 3-D Simple Point Characterizations: A Marching Cubes Paradigm , 2009, IEEE Transactions on Medical Imaging.

[47]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[48]  A. Moura,et al.  Topological flow structures and stir mixing for steady flow in a peripheral bypass graft with uncertainty , 2010 .

[49]  Perumal Nithiarasu,et al.  An artificial compressibility based characteristic based split (CBS) scheme for steady and unsteady turbulent incompressible flows , 2006 .

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

[51]  Onkar Sahni,et al.  Automated adaptive cardiovascular flow simulations , 2009, Engineering with Computers.

[52]  D. Holdsworth,et al.  Characterization of common carotid artery blood-flow waveforms in normal human subjects , 1999, Physiological measurement.

[53]  Nigel P. Weatherill,et al.  A stitching method for the generation of unstructured meshes for use with co-volume solution techniques , 2006 .

[54]  B. Rutt,et al.  Reconstruction of carotid bifurcation hemodynamics and wall thickness using computational fluid dynamics and MRI , 2002, Magnetic resonance in medicine.

[55]  David A. Steinman,et al.  An image-based modeling framework for patient-specific computational hemodynamics , 2008, Medical & Biological Engineering & Computing.

[56]  Alex M. Andrew,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science (2nd edition) , 2000 .

[57]  Azriel Rosenfeld,et al.  Winding and Euler numbers for 2D and 3D digital images , 1991, CVGIP Graph. Model. Image Process..

[58]  Pablo J. Blanco,et al.  Multidimensional modelling for the carotid artery blood flow , 2006 .

[59]  Einar Heiberg,et al.  Feasibility of Patient Specific Aortic Blood Flow CFD Simulation , 2006, MICCAI.

[60]  D. Reutens,et al.  Geometric Model Generation for CFD Simulation of Blood and Air Flows , 2007, 2007 1st International Conference on Bioinformatics and Biomedical Engineering.

[61]  A D Augst,et al.  Various issues relating to computational fluid dynamics simulations of carotid bifurcation flow based on models reconstructed from three-dimensional ultrasound images , 2003, Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine.

[62]  Spencer J. Sherwin,et al.  Automatic reconstruction of a patient-specific high-order surface representation and its application to mesh generation for CFD calculations , 2008, Medical & Biological Engineering & Computing.

[63]  David A. Steinman,et al.  Image-Based Modeling of Blood Flow and Vessel Wall Dynamics: Applications, Methods and Future Directions , 2010, Annals of Biomedical Engineering.

[64]  Erik J. Bekkers,et al.  Multiscale Vascular Surface Model Generation From Medical Imaging Data Using Hierarchical Features , 2008, IEEE Transactions on Medical Imaging.

[65]  Xianghua Xie,et al.  MAC: Magnetostatic Active Contour Model , 2008, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[66]  Peter L. Choyke,et al.  Patient-Specific Simulation of Carotid Artery Stenting Using Computational Fluid Dynamics , 2001, MICCAI.

[67]  Nigel P. Weatherill,et al.  Steady flow through a realistic human upper airway geometry , 2008 .

[68]  Ronald M. Summers,et al.  Tracheal and central bronchial aerodynamics using virtual bronchoscopy and computational fluid dynamics , 2004, IEEE Transactions on Medical Imaging.

[69]  Phillip Colella,et al.  Vessel segmentation and blood flow simulation using Level-Sets and Embedded Boundary methods , 2004, CARS.

[70]  R. Löhner.  Regridding Surface Triangulations , 1996 .

[71]  Yuri Bazilevs,et al.  High-Fidelity Tetrahedral Mesh Generation from Medical Imaging Data for Fluid-Structure Interaction Analysis of Cerebral Aneurysms , 2009 .

[72]  Perumal Nithiarasu,et al.  Application of a locally conservative Galerkin (LCG) method for modelling blood flow through a patient‐specific carotid bifurcation , 2010 .

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

[74]  Kazuhiro Nakahashi,et al.  Robust generation of high‐quality unstructured meshes on realistic biomedical geometry , 2006 .

[75]  Bernd Hamann,et al.  The asymptotic decider: resolving the ambiguity in marching cubes , 1991, Proceeding Visualization '91.

[76]  Mahmoud R. El-Sakka,et al.  Watershed segmentation for carotid artery ultrasound images , 2005, The 3rd ACS/IEEE International Conference onComputer Systems and Applications, 2005..

[77]  N. Weatherill,et al.  Efficient three‐dimensional Delaunay triangulation with automatic point creation and imposed boundary constraints , 1994 .

[78]  M. Kaazempur-Mofrad,et al.  Hemodynamics and wall mechanics in human carotid bifurcation and its consequences for atherogenesis: investigation of inter-individual variation , 2004, Biomechanics and modeling in mechanobiology.

[79]  K S Sakariassen,et al.  Shear-induced platelet activation and platelet microparticle formation at blood flow conditions as in arteries with a severe stenosis. , 1997, Arteriosclerosis, thrombosis, and vascular biology.

[80]  O. Zienkiewicz,et al.  A general algorithm for compressible and incompressible flow—Part I. the split, characteristic‐based scheme , 1995 .

[81]  Nigel P. Weatherill,et al.  Enhanced remeshing from STL files with applications to surface grid generation , 2006 .

[82]  Nigel P. Weatherill,et al.  EQSM: An efficient high quality surface grid generation method based on remeshing , 2006 .

[83]  E Chernyaev,et al.  Marching cubes 33 : construction of topologically correct isosurfaces , 1995 .

[84]  Marcel Breeuwer,et al.  Computational Mesh Generation for Vascular Structures with Deformable Surfaces , 2006, International Journal of Computer Assisted Radiology and Surgery.

[85]  L. Herrmann.  Laplacian-Isoparametric Grid Generation Scheme , 1976 .

[86]  W. M. Gao,et al.  Improved 3D Thinning Algorithms for Skeleton Extraction , 2009, 2009 Digital Image Computing: Techniques and Applications.

[87]  A. Fogelson,et al.  Computational model of whole blood exhibiting lateral platelet motion induced by red blood cells , 2010, International journal for numerical methods in biomedical engineering.

[88]  P. Nithiarasu,et al.  A 1D arterial blood flow model incorporating ventricular pressure, aortic valve and regional coronary flow using the locally conservative Galerkin (LCG) method , 2008 .

[89]  A. Hughes,et al.  Reproducibility Study of Magnetic Resonance Image-Based Computational Fluid Dynamics Prediction of Carotid Bifurcation Flow , 2003, Annals of Biomedical Engineering.

[90]  Charles A. Taylor,et al.  Patient-specific modeling of cardiovascular mechanics. , 2009, Annual review of biomedical engineering.

[91]  Wolfgang A. Wall,et al.  Coupling strategies for biomedical fluid–structure interaction problems , 2010 .

[92]  S. Alper,et al.  Hemodynamic shear stress and its role in atherosclerosis. , 1999, JAMA.

[93]  P. Nithiarasu.  An arbitrary Lagrangian Eulerian (ALE) formulation for free surface flows using the characteristic‐based split (CBS) scheme , 2005 .

[94]  Takashi Hanakawa,et al.  Blood vessel segmentation for head MRA using branch-based region growing , 2005, Systems and Computers in Japan.

[95]  R.W. Dutton,et al.  Improving geometric model construction for blood flow modeling , 1999, IEEE Engineering in Medicine and Biology Magazine.

[96]  Xiangmin Jiao,et al.  Variational generation of prismatic boundary‐layer meshes for biomedical computing , 2009, International journal for numerical methods in engineering.

[97]  van de Fn Frans Vosse,et al.  Computational modelling of endoleak after endovascular repair of abdominal aortic aneurysms , 2010 .

[98]  Ismail Celik,et al.  CFD mesh generation for biological flows: Geometry reconstruction using diagnostic images , 2009 .

[99]  S Giordana,et al.  Automated classification of peripheral distal by-pass geometries reconstructed from medical data. , 2005, Journal of biomechanics.

[100]  J. Peiro,et al.  On the segmentation of vascular geometries from medical images , 2010 .

[101]  Toshio Kobayashi,et al.  Influence of wall thickness on fluid–structure interaction computations of cerebral aneurysms , 2010 .

[102]  William H. Frey,et al.  Mesh relaxation: A new technique for improving triangulations , 1991 .

[103]  P. Nithiarasu,et al.  Laminar and turbulent flow calculations through a model human upper airway using unstructured meshes , 2006 .

[104]  E. V. Chernyaev.  Construction of Topologically Correct Isosurfaces Marching Cubes 33 : , 2009 .

[105]  O. C. Zienkiewicz,et al.  The Finite Element Method for Fluid Dynamics , 2005 .

[106]  F. N. van de Vosse,et al.  Special Issue: Fluid–structure interaction in biomedical applications , 2010 .

[107]  O Smedby,et al.  Do plaques grow upstream or downstream?: an angiographic study in the femoral artery. , 1997, Arteriosclerosis, thrombosis, and vascular biology.

[108]  Rainald Löhner,et al.  Special issue on patient‐specific computational modelling , 2010 .

[109]  M. Kaazempur-Mofrad,et al.  Characterization of the Atherosclerotic Carotid Bifurcation Using MRI, Finite Element Modeling, and Histology , 2004, Annals of Biomedical Engineering.

[110]  Dimitrios V Papavassiliou,et al.  Carotid geometry effects on blood flow and on risk for vascular disease. , 2008, Journal of biomechanics.

[111]  P. Nithiarasu,et al.  Patient‐specific blood flow simulation through an aneurysmal thoracic aorta with a folded proximal neck , 2011 .