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.