GPU acceleration of Volumetric Lattice Boltzmann Method for patient-specific computational hemodynamics

Zhiqiang Wang, Ye Zhao, Alan P. Sawchuck, Michael Dalsing, Huidan Yu

Research output: Contribution to journalArticle

8 Citations (Scopus)

Abstract

Volumetric Lattice Boltzmann Method (VLBM) has been recently developed for solving complex flow with arbitrary curved boundaries. The VLBM regards fluid particles are uniformly distributed in cells and distinguishes fluid, solid, and boundary cells by introducing a volumetric parameter P(. x, t) defining the percentage of solid volume in each cell. The advantages of VLBM stem from the self-regulation of P(. x, t) in the volumetric lattice Boltzmann equation (VLBE) for particle collision and streaming with no spatial interpolation when dealing with an arbitrarily curved boundary with or without motion. First, the VLBE satisfies mass conservation strictly. Second, the implementation of VLBM is rather simple after the solid volume percentages are determined in boundary cells. And third, no-slip boundary condition is integrated in the streaming formulation thus significantly enhances the capability of parallelization. In this paper, we perform GPU (Graphics Processing Unit) parallelization for VLBM using a uniform computing scheme for both fluid and boundary cells. In contrast to the traditional LBM acceleration, the boundary conditions have to be imposed over boundary nodes, where branching operations are required to identify boundary nodes from others, the VLBM implementation does not need to distinguish fluid and boundary cells in the computation so that branching is minimized and the GPU kernel execution is accelerated. Furthermore, the algorithmic steps are optimized to improve coalesced access of GPU memory and avoid race condition. An application study is on a pulsatile blood flow in a patient-specific carotid artery segmented from an anonymous clinical CT image and more than 30 times speedup over the serial counterpart. Simulations of fluid dynamics and wall shear stress (WSS) are presented and known velocity skewness and WSS distributions are captured. The GPU accelerated VLBM is promising to perform patient-specific computational hemodynamics within clinical accepted time frame and is expected to reveal quantitative real-time blood flow in living human arteries to aid clinical assessment of cardiovascular diseases.

Original languageEnglish (US)
Pages (from-to)192-200
Number of pages9
JournalComputers and Fluids
Volume115
DOIs
StatePublished - Jul 2 2015

Fingerprint

Hemodynamics
Fluids
Boltzmann equation
Shear stress
Blood
Hazards and race conditions
Boundary conditions
Fluid dynamics
Stress concentration
Conservation
Interpolation
Data storage equipment
Graphics processing unit

Keywords

  • Graphics processor unit
  • Parallel computing
  • Patient-specific computational hemodynamics
  • Volumetric Lattice Boltzmann Method

ASJC Scopus subject areas

  • Computer Science(all)
  • Engineering(all)

Cite this

GPU acceleration of Volumetric Lattice Boltzmann Method for patient-specific computational hemodynamics. / Wang, Zhiqiang; Zhao, Ye; Sawchuck, Alan P.; Dalsing, Michael; Yu, Huidan.

In: Computers and Fluids, Vol. 115, 02.07.2015, p. 192-200.

Research output: Contribution to journalArticle

Wang, Zhiqiang ; Zhao, Ye ; Sawchuck, Alan P. ; Dalsing, Michael ; Yu, Huidan. / GPU acceleration of Volumetric Lattice Boltzmann Method for patient-specific computational hemodynamics. In: Computers and Fluids. 2015 ; Vol. 115. pp. 192-200.
@article{e12cd725d6894e7bbda6282f45a393ae,
title = "GPU acceleration of Volumetric Lattice Boltzmann Method for patient-specific computational hemodynamics",
abstract = "Volumetric Lattice Boltzmann Method (VLBM) has been recently developed for solving complex flow with arbitrary curved boundaries. The VLBM regards fluid particles are uniformly distributed in cells and distinguishes fluid, solid, and boundary cells by introducing a volumetric parameter P(. x, t) defining the percentage of solid volume in each cell. The advantages of VLBM stem from the self-regulation of P(. x, t) in the volumetric lattice Boltzmann equation (VLBE) for particle collision and streaming with no spatial interpolation when dealing with an arbitrarily curved boundary with or without motion. First, the VLBE satisfies mass conservation strictly. Second, the implementation of VLBM is rather simple after the solid volume percentages are determined in boundary cells. And third, no-slip boundary condition is integrated in the streaming formulation thus significantly enhances the capability of parallelization. In this paper, we perform GPU (Graphics Processing Unit) parallelization for VLBM using a uniform computing scheme for both fluid and boundary cells. In contrast to the traditional LBM acceleration, the boundary conditions have to be imposed over boundary nodes, where branching operations are required to identify boundary nodes from others, the VLBM implementation does not need to distinguish fluid and boundary cells in the computation so that branching is minimized and the GPU kernel execution is accelerated. Furthermore, the algorithmic steps are optimized to improve coalesced access of GPU memory and avoid race condition. An application study is on a pulsatile blood flow in a patient-specific carotid artery segmented from an anonymous clinical CT image and more than 30 times speedup over the serial counterpart. Simulations of fluid dynamics and wall shear stress (WSS) are presented and known velocity skewness and WSS distributions are captured. The GPU accelerated VLBM is promising to perform patient-specific computational hemodynamics within clinical accepted time frame and is expected to reveal quantitative real-time blood flow in living human arteries to aid clinical assessment of cardiovascular diseases.",
keywords = "Graphics processor unit, Parallel computing, Patient-specific computational hemodynamics, Volumetric Lattice Boltzmann Method",
author = "Zhiqiang Wang and Ye Zhao and Sawchuck, {Alan P.} and Michael Dalsing and Huidan Yu",
year = "2015",
month = "7",
day = "2",
doi = "10.1016/j.compfluid.2015.04.004",
language = "English (US)",
volume = "115",
pages = "192--200",
journal = "Computers and Fluids",
issn = "0045-7930",
publisher = "Elsevier Limited",

}

TY - JOUR

T1 - GPU acceleration of Volumetric Lattice Boltzmann Method for patient-specific computational hemodynamics

AU - Wang, Zhiqiang

AU - Zhao, Ye

AU - Sawchuck, Alan P.

AU - Dalsing, Michael

AU - Yu, Huidan

PY - 2015/7/2

Y1 - 2015/7/2

N2 - Volumetric Lattice Boltzmann Method (VLBM) has been recently developed for solving complex flow with arbitrary curved boundaries. The VLBM regards fluid particles are uniformly distributed in cells and distinguishes fluid, solid, and boundary cells by introducing a volumetric parameter P(. x, t) defining the percentage of solid volume in each cell. The advantages of VLBM stem from the self-regulation of P(. x, t) in the volumetric lattice Boltzmann equation (VLBE) for particle collision and streaming with no spatial interpolation when dealing with an arbitrarily curved boundary with or without motion. First, the VLBE satisfies mass conservation strictly. Second, the implementation of VLBM is rather simple after the solid volume percentages are determined in boundary cells. And third, no-slip boundary condition is integrated in the streaming formulation thus significantly enhances the capability of parallelization. In this paper, we perform GPU (Graphics Processing Unit) parallelization for VLBM using a uniform computing scheme for both fluid and boundary cells. In contrast to the traditional LBM acceleration, the boundary conditions have to be imposed over boundary nodes, where branching operations are required to identify boundary nodes from others, the VLBM implementation does not need to distinguish fluid and boundary cells in the computation so that branching is minimized and the GPU kernel execution is accelerated. Furthermore, the algorithmic steps are optimized to improve coalesced access of GPU memory and avoid race condition. An application study is on a pulsatile blood flow in a patient-specific carotid artery segmented from an anonymous clinical CT image and more than 30 times speedup over the serial counterpart. Simulations of fluid dynamics and wall shear stress (WSS) are presented and known velocity skewness and WSS distributions are captured. The GPU accelerated VLBM is promising to perform patient-specific computational hemodynamics within clinical accepted time frame and is expected to reveal quantitative real-time blood flow in living human arteries to aid clinical assessment of cardiovascular diseases.

AB - Volumetric Lattice Boltzmann Method (VLBM) has been recently developed for solving complex flow with arbitrary curved boundaries. The VLBM regards fluid particles are uniformly distributed in cells and distinguishes fluid, solid, and boundary cells by introducing a volumetric parameter P(. x, t) defining the percentage of solid volume in each cell. The advantages of VLBM stem from the self-regulation of P(. x, t) in the volumetric lattice Boltzmann equation (VLBE) for particle collision and streaming with no spatial interpolation when dealing with an arbitrarily curved boundary with or without motion. First, the VLBE satisfies mass conservation strictly. Second, the implementation of VLBM is rather simple after the solid volume percentages are determined in boundary cells. And third, no-slip boundary condition is integrated in the streaming formulation thus significantly enhances the capability of parallelization. In this paper, we perform GPU (Graphics Processing Unit) parallelization for VLBM using a uniform computing scheme for both fluid and boundary cells. In contrast to the traditional LBM acceleration, the boundary conditions have to be imposed over boundary nodes, where branching operations are required to identify boundary nodes from others, the VLBM implementation does not need to distinguish fluid and boundary cells in the computation so that branching is minimized and the GPU kernel execution is accelerated. Furthermore, the algorithmic steps are optimized to improve coalesced access of GPU memory and avoid race condition. An application study is on a pulsatile blood flow in a patient-specific carotid artery segmented from an anonymous clinical CT image and more than 30 times speedup over the serial counterpart. Simulations of fluid dynamics and wall shear stress (WSS) are presented and known velocity skewness and WSS distributions are captured. The GPU accelerated VLBM is promising to perform patient-specific computational hemodynamics within clinical accepted time frame and is expected to reveal quantitative real-time blood flow in living human arteries to aid clinical assessment of cardiovascular diseases.

KW - Graphics processor unit

KW - Parallel computing

KW - Patient-specific computational hemodynamics

KW - Volumetric Lattice Boltzmann Method

UR - http://www.scopus.com/inward/record.url?scp=84927780662&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84927780662&partnerID=8YFLogxK

U2 - 10.1016/j.compfluid.2015.04.004

DO - 10.1016/j.compfluid.2015.04.004

M3 - Article

AN - SCOPUS:84927780662

VL - 115

SP - 192

EP - 200

JO - Computers and Fluids

JF - Computers and Fluids

SN - 0045-7930

ER -