An Effective Finite Difference Method for Simulation of Bioheat Transfer in Irregular Tissues

A three-dimensional (3D) simulation of bioheat transfer is crucial to analyze the physiological processes and evaluate many therapeutic/diagnostic practices spanning from high to low temperature medicine. In this paper we develop an efficient numerical scheme for solving 3D transient bioheat transfer equations based on the alternating direction implicit finite-difference method (ADI-FDM). An algorithm is proposed to deal with the boundary condition for irregular domain which could capture accurately the complex boundary and reduce considerably the staircase effects. Furthermore, the local adaptive mesh technology is introduced to improve the computational accuracy for irregular boundary and the domains with large temperature gradient. The detailed modification to ADI-FDM is given to accommodate such special grid structure, in particular. Combination of adaptive-mesh technology and ADI-FDM could significantly improve the computational accuracy and decrease the computational cost. Extensive results of numerical experiments demonstrate that the algorithm developed in the current work is very effective to predict the temperature distribution during hyperthermia and cryosurgery. This work may play an important role in developing a computational planning tool for hyperthermia and cryosurgery in the near future. [DOI: 10.1115/1.4024064]

[1]  Frédéric Gibou,et al.  A Supra-Convergent Finite Difference Scheme for the Poisson and Heat Equations on Irregular Domains and Non-Graded Adaptive Cartesian Grids , 2007, J. Sci. Comput..

[2]  Jean-Louis Dillenseger,et al.  Fast FFT-based bioheat transfer equation computation , 2010, Comput. Biol. Medicine.

[3]  Jing Liu,et al.  Mathematical modeling of temperature mapping over skin surface and its implementation in thermal disease diagnostics , 2004, Comput. Biol. Medicine.

[4]  Zhong-Shan Deng,et al.  Modeling of multidimensional freezing problem during cryosurgery by the dual reciprocity boundary element method , 2004 .

[5]  C. Chan,et al.  Boundary element method analysis for the bioheat transfer equation. , 1992, Journal of biomechanical engineering.

[6]  M. Wernick,et al.  A physical model of multiple-image radiography , 2006, Physics in medicine and biology.

[7]  J. K. Chen,et al.  A dual reciprocity boundary element method for photothermal interactions in laser-induced thermotherapy , 2008 .

[8]  J. L. Arce-Diego,et al.  Efficient 3D numerical approach for temperature prediction in laser irradiated biological tissues , 2009, Comput. Biol. Medicine.

[9]  James C. Lin,et al.  Power density and temperature distributions produced by interstitial arrays of sleeved-slot antennas for hyperthermic cancer therapy , 2003 .

[10]  Youssef Belhamadia,et al.  Numerical prediction of freezing fronts in cryosurgery: Comparison with experimental results , 2005, Computer methods in biomechanics and biomedical engineering.

[11]  Peng Zeng,et al.  Parallel Algorithms for Freezing Problems during Cryosurgery , 2011 .

[12]  Zhong-Shan Deng,et al.  MONTE CARLO METHOD TO SOLVE MULTIDIMENSIONAL BIOHEAT TRANSFER PROBLEM , 2002 .

[13]  Luisa Consiglieri,et al.  Theoretical analysis of the heat convection coefficient in large vessels and the significance for thermal ablative therapies. , 2003, Physics in medicine and biology.

[14]  Daniele Bertaccini,et al.  Fast simulation of solid tumors thermal ablation treatments with a 3D reaction diffusion model , 2007, Comput. Biol. Medicine.

[15]  Jun Zhang,et al.  A numerical study of a 3D bioheat transfer problem with different spatial heating , 2005, Math. Comput. Simul..

[16]  Guirong Liu,et al.  An efficient algorithm for phase change problem in tumor treatment using αFEM , 2010 .

[17]  E Y K Ng,et al.  Boundary element method with bioheat equation for skin burn injury. , 2009, Burns : journal of the International Society for Burn Injuries.

[18]  John G. Webster,et al.  Finite-element analysis of hepatic multiple probe radio-frequency ablation , 2002, IEEE Transactions on Biomedical Engineering.

[19]  Win-Li Lin,et al.  Effects of pulsatile blood flow in large vessels on thermal dose distribution during thermal therapy. , 2007, Medical physics.

[20]  E Neufeld,et al.  Novel conformal technique to reduce staircasing artifacts at material boundaries for FDTD modeling of the bioheat equation , 2007, Physics in medicine and biology.

[21]  A distance-function-based Cartesian (DIFCA) grid method for irregular geometries , 2008 .

[22]  Haifeng Zhang,et al.  Lattice Boltzmann method for solving the bioheat equation , 2008, Physics in medicine and biology.

[23]  Ayda Saïdane,et al.  Thermal analysis of a three-dimensional breast model with embedded tumour using the transmission line matrix (TLM) method , 2011, Comput. Biol. Medicine.

[24]  Wen-Qiang Lu,et al.  Simulation of the thermal wave propagation in biological tissues by the dual reciprocity boundary element method , 1998 .

[25]  E. Y. K. Ng,et al.  Bioheat transfer in the human eye: a boundary element approach , 2007 .

[26]  Derek M. Causon,et al.  Developments in Cartesian cut cell methods , 2003, Math. Comput. Simul..

[27]  A Saidane,et al.  Transmission line matrix modelling of thermal injuries to skin. , 2008, Burns : journal of the International Society for Burn Injuries.

[28]  Kenji Shimada,et al.  An efficient numerical technique for bioheat simulations and its application to computerized cryosurgery planning , 2007, Comput. Methods Programs Biomed..

[29]  N Kuster,et al.  Effects of geometry discretization aspects on the numerical solution of the bioheat transfer equation with the FDTD technique , 2006, Physics in medicine and biology.

[30]  Jun Zhang,et al.  A two level finite difference scheme for one dimensional Pennes' bioheat equation , 2005, Appl. Math. Comput..

[31]  K. Muldrew,et al.  A finite element model for cryosurgery with coupled phase change and thermal stress aspects , 2008 .

[32]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

[33]  Jing Liu,et al.  Theoretical evaluation of burns to the human respiratory tract due to inhalation of hot gas in the early stage of fires. , 2006, Burns : journal of the International Society for Burn Injuries.

[34]  R B Roemer,et al.  Engineering aspects of hyperthermia therapy. , 1999, Annual review of biomedical engineering.

[35]  Mario Chiampi,et al.  A boundary element approach to relate surface fields with the specific absorption rate (SAR) induced in 3-D human phantoms , 2011 .

[36]  J. Douglas,et al.  A general formulation of alternating direction methods , 1964 .

[37]  J. Bischof,et al.  Quantitative measurement and prediction of biophysical response during freezing in tissues. , 2000, Annual review of biomedical engineering.

[38]  Ning Zhao,et al.  An RBF–MFS model for analysing thermal behaviour of skin tissues , 2010 .

[39]  Qing Hua Qin,et al.  FE approach with Green's function as internal trial function for simulating bioheat transfer in the human eye , 2010 .