A meshless point collocation treatment of transient bioheat problems

A meshless numerical method is proposed for the solution of the transient bioheat equation in two and three dimensions. The Pennes bioheat equation is extended in order to incorporate water evaporation, tissue damage, and temperature-dependent tissue properties during tumor ablation. The conductivity of the tissue is not assumed constant but is treated as a local function to simulate local variability due to the existence of usually unclear interfacing of healthy and pathological segments. In this way, one avoids the need for accurate identification of the boundaries between pathological and healthy regions, which is a typical problem in medical practice, and sidesteps, evidently, the corresponding mathematical treatment of such boundaries, which is usually a tedious procedure with some inevitable degree of approximation. The numerical results of the new method for test applications of the bioheat transfer equation are validated against analytical predictions and predictions of other numerical methods. 3D simulations are presented that involve the modeling of tumor ablation and account for metabolic heat generation, blood perfusion, and heat ablation using realistic values for the various parameters. An evaluation of the effective medium approximation to homogenize conductivity fields for use with the bioheat equation is also provided.

[1]  Boris Gramatikov.  "Handbook of biomedical optics", edited by David A. Boas, Constantinos Pitris, and Nimmi Ramanujam , 2012, BioMedical Engineering OnLine.

[2]  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.

[3]  Shadi Mahjoob,et al.  Analytical characterization of heat transport through biological media incorporating hyperthermia treatment , 2009 .

[4]  Do Wan Kim,et al.  Meshfree method for the non-stationary incompressible Navier-Stokes equations , 2005 .

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

[6]  Diego Rossinelli,et al.  Vortex methods for incompressible flow simulations on the GPU , 2008, The Visual Computer.

[7]  J. Zee.  Heating the patient : a promising approach ? , 2002 .

[8]  Q. Nan,et al.  Thermal characteristics of microwave ablation in the vicinity of an arterial bifurcation , 2006, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

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

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

[11]  Piet Hut,et al.  A hierarchical O(N log N) force-calculation algorithm , 1986, Nature.

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

[13]  P. Lancaster,et al.  Surfaces generated by moving least squares methods , 1981 .

[14]  Mark A. Christon,et al.  Multi-scale meshfree parallel computations for viscous, compressible flows , 2000 .

[15]  Q.,et al.  APPLICATION OF HYBRID-TREFFTZ ELEMENT APPROACH TO TRANSIENT HEAT CONDUCTION ANALYSIS , 2003 .

[16]  Guirong Liu.  Mesh Free Methods: Moving Beyond the Finite Element Method , 2002 .

[17]  Qing Hua Qin,et al.  A meshless model for transient heat conduction in functionally graded materials , 2006 .

[18]  W. Dewey,et al.  Arrhenius relationships from the molecule and cell to the clinic , 2009, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[19]  H. H. Pennes.  Analysis of tissue and arterial blood temperatures in the resting human forearm. , 1948, Journal of applied physiology.

[20]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[21]  Jing Liu,et al.  Boundary information based diagnostics on the thermal states of biological bodies , 2000 .

[22]  B. Hooper.  Optical-thermal response of laser-irradiated tissue , 1996 .

[23]  George Nikiforidis,et al.  Truly meshless localized type techniques for the steady-state heat conduction problems for isotropic and functionally graded materials , 2011 .

[24]  H. Arkin,et al.  Recent developments in modeling heat transfer in blood perfused tissues , 1994, IEEE Transactions on Biomedical Engineering.

[25]  Chris J. Diederich,et al.  Hyperthermia classic commentary: ‘Arrhenius relationships from the molecule and cell to the clinic’ by William Dewey, Int. J. Hyperthermia, 10:457–483, 1994 , 2009, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[26]  Ivo F. Sbalzarini,et al.  Fast neighbor lists for adaptive-resolution particle simulations , 2012, Comput. Phys. Commun..

[27]  P R Stauffer,et al.  Phantom and animal tissues for modelling the electrical properties of human liver , 2003, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[28]  Mehdi Dehghan,et al.  A spectral element method for solving the Pennes bioheat transfer equation by using triangular and quadrilateral elements , 2012 .

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

[30]  Edward Ng,et al.  Evaluation of meshless radial basis collocation method (RBCM) for heterogeneous conduction and simulation of temperature inside the biological tissues , 2013 .

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

[32]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[33]  Nail A. Gumerov,et al.  Fast parallel Particle-To-Grid interpolation for plasma PIC simulations on the GPU , 2008, J. Parallel Distributed Comput..

[34]  Heung-Jin Chung,et al.  Adaptive crack propagation analysis with the element‐free Galerkin method , 2003 .

[35]  K. Delman,et al.  The role of hyperthermia in optimizing tumor response to regional therapy , 2008, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

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

[37]  Punit Prakash,et al.  Theoretical Modeling for Hepatic Microwave Ablation , 2010, The open biomedical engineering journal.

[38]  J W Hunt,et al.  Differential thermal sensitivity of tumour and normal tissue microvascular response during hyperthermia. , 1992, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[39]  G. Bourantas,et al.  Numerical simulation of the unsteady non-linear heat transfer problems. Application on nanosecond laser annealing of Si , 2012 .

[40]  Alain J. Kassab,et al.  A meshless method for conjugate heat transfer problems , 2005 .

[41]  John A. Pearce,et al.  Relationship between Arrhenius models of thermal damage and the CEM 43 thermal dose , 2009, BiOS.

[42]  X. X. Zhang,et al.  Effects of dynamic changes of tissue properties during laser-induced interstitial thermotherapy (LITT) , 2005, Lasers in Medical Science.

[43]  William Gropp,et al.  A Parallel Version of the Fast Multipole Method-Invited Talk , 1987, PPSC.

[44]  Gang Zhao,et al.  Effect of blood flow and metabolism on multidimensional heat transfer during cryosurgery. , 2007, Medical engineering & physics.

[45]  Dieter Haemmerich,et al.  Effect of variable heat transfer coefficient on tissue temperature next to a large vessel during radiofrequency tumor ablation , 2008, BioMedical Engineering OnLine.

[46]  Muhammad Jamil,et al.  To optimize the efficacy of bioheat transfer in capacitive hyperthermia : a physical perspective , 2013 .

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

[48]  Hui-Ping Wang,et al.  A Lagrangian reproducing kernel particle method for metal forming analysis , 1998 .

[49]  J C Bischof,et al.  Investigation of the thermal and tissue injury behaviour in microwave thermal therapy using a porcine kidney model , 2004, International journal of hyperthermia : the official journal of European Society for Hyperthermic Oncology, North American Hyperthermia Group.

[50]  D. A. G. Bruggeman.  Berechnung verschiedener physikalischer Konstanten von heterogenen Substanzen. I. Dielektrizitätskonstanten und Leitfähigkeiten der Mischkörper aus isotropen Substanzen , 1935 .

[51]  Dirk Diehl,et al.  Evaluation of MR-Induced Hot Spots for Different Temporal SAR Modes Using a Time-Dependent Finite Difference Method With Explicit Temperature Gradient Treatment , 2007, IEEE Transactions on Biomedical Engineering.

[52]  Muhammad Jamil,et al.  Ranking of parameters in bioheat transfer using Taguchi analysis , 2013 .

[53]  Jan J W Lagendijk,et al.  Prostate perfusion in patients with locally advanced prostate carcinoma treated with different hyperthermia techniques. , 2002, The Journal of urology.

[54]  John G. Webster,et al.  Measurement and Analysis of Tissue Temperature During Microwave Liver Ablation , 2007, IEEE Transactions on Biomedical Engineering.

[55]  Muhammad Jamil,et al.  Statistical modeling of electrode based thermal therapy with Taguchi based multiple regression , 2013 .

[56]  G. R. Johnson,et al.  A generalized particle algorithm for high velocity impact computations , 2000 .

[57]  Ping Yuan,et al.  Temperature response in biological tissue by alternating heating and cooling modalities with sinusoidal temperature oscillation on the skin , 2008 .

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