TY - JOUR

T1 - Numerical simulation of three-dimensional blood flows using domain decomposition method on parallel computer

AU - Hwang, Feng Nan

AU - Wu, Chao Ying

AU - Cai, Xiao Chuan

PY - 2010/6

Y1 - 2010/6

N2 - A good numerical blood flow simulation tool based on patient-specific anatomy and physiological conditions can be clinically helpful for physicians or researchers to study vascular diseases, to enhance diagnoses, as well as to plan surgery procedures. Such a tool is computationally very expensive, and often requires the use of large scale supercomputers with many core processors. In this paper, we focus on developing parallel domain decomposition algorithms for solving nonlinear systems arising from the discretization of three-dimensional blood flow model equations with a stabilized finite element method for the spatial variables and an implicit backward Euler finite difference method for the temporal variable. More precisely speaking, at each time step, the resulting nonlinear system is solved by the Newton-Krylov-Schwarz algorithm. We implement the parallel fluid solver using PETSc and integrate it with other state-of-the-art software packages into a parallel blood flow simulation system, which includes Cubit, ParMETIS and ParaView for mesh generation, mesh partitioning, and visualization, respectively. We validate our parallel code and investigate the parallel performance of our algorithms for both a straight artery model and an end-to-side graft model.

AB - A good numerical blood flow simulation tool based on patient-specific anatomy and physiological conditions can be clinically helpful for physicians or researchers to study vascular diseases, to enhance diagnoses, as well as to plan surgery procedures. Such a tool is computationally very expensive, and often requires the use of large scale supercomputers with many core processors. In this paper, we focus on developing parallel domain decomposition algorithms for solving nonlinear systems arising from the discretization of three-dimensional blood flow model equations with a stabilized finite element method for the spatial variables and an implicit backward Euler finite difference method for the temporal variable. More precisely speaking, at each time step, the resulting nonlinear system is solved by the Newton-Krylov-Schwarz algorithm. We implement the parallel fluid solver using PETSc and integrate it with other state-of-the-art software packages into a parallel blood flow simulation system, which includes Cubit, ParMETIS and ParaView for mesh generation, mesh partitioning, and visualization, respectively. We validate our parallel code and investigate the parallel performance of our algorithms for both a straight artery model and an end-to-side graft model.

KW - Blood flow modeling

KW - Domain decomposition methods

KW - Newton-Krylov-Schwarz algorithm

KW - Parallel processing

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

M3 - 期刊論文

AN - SCOPUS:79955104635

SN - 0257-9731

VL - 31

SP - 199

EP - 208

JO - Journal of the Chinese Society of Mechanical Engineers, Transactions of the Chinese Institute of Engineers, Series C/Chung-Kuo Chi Hsueh Kung Ch'eng Hsuebo Pao

JF - Journal of the Chinese Society of Mechanical Engineers, Transactions of the Chinese Institute of Engineers, Series C/Chung-Kuo Chi Hsueh Kung Ch'eng Hsuebo Pao

IS - 3

ER -