Novel 3D GPU based numerical parallel diffusion algorithms in cylindrical coordinates for health care simulation

Modeling diffusion processes, such as drug deliver, bio-heat transfer, and the concentration change of cytokine for computational biology research, requires intensive computing resources as one must employ sequential numerical algorithms to obtain accurate numerical solutions, especially for real-time in vivo 3D simulation. Thus, it is necessary to develop a new numerical algorithm compatible with state-of-the-art computing hardware. The purpose of this article is to integrate the graphics processing unit (GPU) technology with the locally-one-dimension (LOD) numerical method for solving partial differential equations, and to develop a novel 3D numerical parallel diffusion algorithm (GNPD) in cylindrical coordinates based on GPU technology, which can be used in the neuromuscular junction research.