Multi-GPU Based Lattice Boltzmann Method for Hemodynamic Simulation in Patient-Specific Cerebral Aneurysm

Conducting lattice Boltzmann method on GPU has been proved to be an effective manner to gain a significant performance benefit, thus the GPU or multi-GPU based lattice Boltzmann method is considered as a promising and competent candidate in the study of large-scale complex fluid flows. In this work, a multi-GPU based lattice Boltzmann algorithm coupled with the sparse lattice representation and message passing interface is presented. Some numerical tests are also carried out, and the results show that a parallel efficiency close to 90% can be achieved on a single-node cluster equipped with four GPU cards. Then the proposed algorithm is adopted to study the hemodynamics of patient-specific cerebral aneurysm with stent implanted. It is found that the stent can apparently reduce the aneurysmal inflow and improve the hemodynamic environment. This work also shows that the lattice Boltzmann method running on the GPU platform is a powerful tool to study the fluid mechanism within the aneurysms and enable us to better understand the pathogenesis and treatment of cerebral aneurysms.

[1]  Manfred Krafczyk,et al.  TeraFLOP computing on a desktop PC with GPUs for 3D CFD , 2008 .

[2]  R. Benzi,et al.  The lattice Boltzmann equation: theory and applications , 1992 .

[3]  Thomas Zeiser,et al.  Performance evaluation of a parallel sparse lattice Boltzmann solver , 2008, J. Comput. Phys..

[4]  J. Boon.  The Lattice Boltzmann Equation for Fluid Dynamics and Beyond , 2003 .

[5]  R. Sclabassi,et al.  Saccular aneurysm formation in curved and bifurcating arteries. , 1999, AJNR. American journal of neuroradiology.

[6]  Shiyi Chen,et al.  LATTICE BOLTZMANN METHOD FOR FLUID FLOWS , 2001 .

[7]  Shi Bao-Chang,et al.  Simulating high Reynolds number flow in two-dimensional lid-driven cavity by multi-relaxation-time lattice Boltzmann method , 2006 .

[8]  Joon Sang Lee,et al.  The Effect of Stent Porosity and Strut Shape on Saccular Aneurysm and its Numerical Analysis with Lattice Boltzmann Method , 2010, Annals of Biomedical Engineering.

[9]  Dinan Wang,et al.  Non-Newtonian blood flow simulation in cerebral aneurysms , 2009, Comput. Math. Appl..

[10]  Bernard Tourancheau,et al.  Multi-GPU implementation of the lattice Boltzmann method , 2013, Comput. Math. Appl..

[11]  L. Luo,et al.  Theory of the lattice Boltzmann method: From the Boltzmann equation to the lattice Boltzmann equation , 1997 .

[12]  Zhenhua Chai,et al.  Non-Newtonian Effect on Hemodynamic Characteristics of Blood Flow in Stented Cerebral Aneurysm , 2013 .

[13]  Stefan Turek,et al.  Simulation of Intra-Aneurysmal Blood Flow by Different Numerical Methods , 2013, Comput. Math. Methods Medicine.

[14]  F Mut,et al.  Clinical application of image‐based CFD for cerebral aneurysms , 2011, International journal for numerical methods in biomedical engineering.

[15]  B. Shi,et al.  An extrapolation method for boundary conditions in lattice Boltzmann method , 2002 .

[16]  Massimo Bernaschi,et al.  A flexible high-performance Lattice Boltzmann GPU code for the simulations of fluid flows in complex geometries , 2010 .

[17]  Bernard Tourancheau,et al.  Global Memory Access Modelling for Efficient Implementation of the Lattice Boltzmann Method on Graphics Processing Units , 2010, VECPAR.

[18]  Ernst Rank,et al.  Parallelization Strategies and Efficiency of CFD Computations in Complex Geometries Using Lattice Boltzmann Methods on High-Performance Computers , 2002 .

[19]  L. Luo,et al.  Analytic solutions of simple flows and analysis of nonslip boundary conditions for the lattice Boltzmann BGK model , 1997 .

[20]  Cass T. Miller,et al.  A high-performance lattice Boltzmann implementation to model flow in porous media , 2003 .

[21]  Alejandro F Frangi,et al.  Reproducibility of haemodynamical simulations in a subject-specific stented aneurysm model--a report on the Virtual Intracranial Stenting Challenge 2007. , 2008, Journal of biomechanics.

[22]  Y. Qian,et al.  Lattice BGK Models for Navier-Stokes Equation , 1992 .

[23]  Ge Wei,et al.  Efficient parallel implementation of the lattice Boltzmann method on large clusters of graphic processing units , 2012 .

[24]  David A Steinman,et al.  Finite-element modeling of the hemodynamics of stented aneurysms. , 2004, Journal of biomechanical engineering.

[25]  Rafik Ouared,et al.  A lattice Boltzmann simulation of clotting in stented aneursysms and comparison with velocity or shear rate reductions , 2006, Math. Comput. Simul..

[26]  Miki Hirabayashi,et al.  Characterization of flow reduction properties in an aneurysm due to a stent. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  C. Putman,et al.  Characterization of cerebral aneurysms for assessing risk of rupture by using patient-specific computational hemodynamics models. , 2005, AJNR. American journal of neuroradiology.

[28]  G. Duckwiler,et al.  Lattice Boltzmann simulation of cerebral artery hemodynamics , 2009 .

[29]  C. Putman,et al.  Aneurysm Rupture Following Treatment with Flow-Diverting Stents: Computational Hemodynamics Analysis of Treatment , 2010, American Journal of Neuroradiology.

[30]  Bernard Tourancheau,et al.  The TheLMA project: Multi-GPU implementation of the lattice Boltzmann method , 2011, Int. J. High Perform. Comput. Appl..

[31]  Christian Obrecht,et al.  LBM based flow simulation using GPU computing processor , 2010, Comput. Math. Appl..

[32]  Massimo Bernaschi,et al.  Petaflop hydrokinetic simulations of complex flows on massive GPU clusters , 2013, Comput. Phys. Commun..

[33]  Rainald Löhner,et al.  Simulation of intracranial aneurysm stenting: Techniques and challenges , 2009 .