Research article

Evaluation of a radial basis function node refinement algorithm applied to bioheat transfer modeling

  • Received: 30 October 2020 Accepted: 17 December 2020 Published: 23 December 2020
  • Many bioheat transfer problems involve linear/non-linear equations with non-linear or time-dependent boundary conditions. For heat transfer problems, the presence of time and space-dependent functions under Neumann and Mixed type boundary conditions characterize trivial applications in bioengineering, such as thermotherapies, laser surgeries, and burn studies. This greatly increases the complexity of the numerical solution in several problems, requiring fast and accurate numerical solutions. This paper has a main objective evaluate an adaptive mesh refinement radial basis function method strategy for the classical Penne's bioheat transfer modeling. Our numerical results had errors of ~0.1% compared to analytical solutions. Thus, the proposed methodology is accurate and has a low computational cost. For step function heating, two RBF shape parameters were applied, again achieving excellent results. The distributions of the nodes in the solution domain show that the primary source of error in the numerical solutions came from the boundary conditions. This finding should arouse the interest of engineers and scientists in the development of new strategies for problems involving boundary conditions with periodic functions.

    Citation: Rafael Pinheiro Amantéa, Robspierre de Carvalho, Luiz Otávio Barbosa. Evaluation of a radial basis function node refinement algorithm applied to bioheat transfer modeling[J]. AIMS Bioengineering, 2021, 8(1): 36-51. doi: 10.3934/bioeng.2021005

    Related Papers:

  • Many bioheat transfer problems involve linear/non-linear equations with non-linear or time-dependent boundary conditions. For heat transfer problems, the presence of time and space-dependent functions under Neumann and Mixed type boundary conditions characterize trivial applications in bioengineering, such as thermotherapies, laser surgeries, and burn studies. This greatly increases the complexity of the numerical solution in several problems, requiring fast and accurate numerical solutions. This paper has a main objective evaluate an adaptive mesh refinement radial basis function method strategy for the classical Penne's bioheat transfer modeling. Our numerical results had errors of ~0.1% compared to analytical solutions. Thus, the proposed methodology is accurate and has a low computational cost. For step function heating, two RBF shape parameters were applied, again achieving excellent results. The distributions of the nodes in the solution domain show that the primary source of error in the numerical solutions came from the boundary conditions. This finding should arouse the interest of engineers and scientists in the development of new strategies for problems involving boundary conditions with periodic functions.

    Nomenclature

    Symbol Description Unit
    θc Coarse node parameter -
    θr Refine node parameter -
    ωb Blood perfusion m3/s/m3
    Ø Radial basis function -
    Δt Numerical time step s
    c Specific heat of tissue J/(kg·K)
    cb Specific heat of blood J/(kg·K)
    h0 Heat convection coeficient W/(m2·K)
    k Thermal conductivity of tissue W/(m·K)
    L Distance between skin surface and body core m
    Qm Metabolic rate of tissue W/m3
    Qr Spatial heating W/m3
    r Euclidean distance -
    t Time s
    T Tissue temperature °C
    T0 Steady state temperature °C
    Ta Arterial temperature °C
    Tc Blood temperature °C
    Tf Fluid temperature °C
    x Spatial coordinate m
    α Thermal diffusivity m2/s
    ϵ Multiquadric shape parameter -
    λ Radial basis function interpolator -
    ρ Density of tissue Kg/m3
    ρb Density of blood Kg/m3
    加载中


    [1] Pennes HH (1948) Analysis of tissue and arterial blood temperature in the resting human forearm. J Appl Phisiol 1: 93-122.
    [2] Dillenseger JL, Esneault S (2010) Fast FFT-based bioheat transfer equation computation. Comput Biol Med 40: 119-123.
    [3] Zhao JJ, Zhang J, Kang N, et al. (2005) A two level finite difference scheme for one dimensional Pennes' bioheat equation. Appl Math Comput 171: 320-331.
    [4] Huang HW, Horng TL (2015) Bioheat transfer and thermal heating for tumor treatment. Heat Transfer and Fluid Flow in Biological Processes Amsterdam: Elsevier, 1-42.
    [5] Bedin L, Bazán FSV (2014) On the 2D bioheat equation with convective boundary conditions and its numerical realization via a highly accurate approach. Appl Math Comput 236: 422-436.
    [6] Karaa S, Zhang J, Yang F (2005) A numerical study of a 3D bioheat transfer problem with different spatial heating. Math Comput Simulat 68: 375-388.
    [7] Askarizadeh H, Ahmadikia H (2015) Analytical study on the transient heating of a two-dimensional skin tissue using parabolic and hyperbolic bioheat transfer equations. Appl Math Model 39: 3704-3720.
    [8] Deng ZS, Liu J (2002) Analytical study on bioheat transfer problems with spatial or transient heating on skin surface or inside biological bodies. J Biomech Eng 124: 638-649.
    [9] Dutta J, Kundu B (2017) A revised approach for an exact analytical solution for thermal response in biological tissues significant in therapeutic treatments. J Therm Biol 66: 33-48.
    [10] Bojdi ZK, Hemmat AA (2017) Wavelet collocation methods for solving the Pennes bioheat transfer equation. Optik 130: 345-355.
    [11] Kansa EJ (1990) Multiquadrics—A scattered data approximation scheme with applications to computational fluid-dynamics—II solutions to parabolic, hyperbolic and elliptic partial differential equations. Comput Math Appl 19: 147-161.
    [12] Beatson RK, Levesley J, Mouat CT (2011) Better bases for radial basis function interpolation problems. J Comput Appl Math 236: 434-446.
    [13] Zerroukat M, Power H, Chen CS (1998) A numerical method for heat transfer problems using collocation and radial basis functions. Int J Numer Meth Eng 42: 1263-1278.
    [14] Cao L, Qin QH, Zhao N (2010) An RBF-MFS model for analysing thermal behaviour of skin tissues. Int J Heat Mass Tran 53: 1298-1307.
    [15] Jamil M, Ng EYK (2013) Evaluation of meshless radial basis collocation method (RBCM) for heterogeneous conduction and simulation of temperature inside the biological tissues. Int J Therm Sci 68: 42-52.
    [16] Hon YC, Mao XZ (1998) An efficient numerical scheme for Burgers' equation. Appl Math Comput 95: 37-50.
    [17] Sarra SA (2005) Adaptive radial basis function methods for time dependent partial differential equations. Appl Numer Math 54: 79-94.
    [18] Driscoll TA, Heryudono ARH (2007) Adaptive residual subsampling methods for radial basis function interpolation and collocation problems. Comput Math Appl 53: 927-939.
    [19] Zhang J, Chauhan S (2019) Real-time computation of bio-heat transfer in the fast explicit dynamics finite element algorithm (FED-FEM) framework. Numer Heat Transfer, Part B 75: 217-238.
    [20] Zhang J, Chauhan S (2019) Neural network methodology for real-time modelling of bio-heat transfer during thermo-therapeutic applications. Artif Intell Med 101: 101728.
    [21] Fahmy MA (2019) Boundary element modeling and simulation of biothermomechanical behavior in anisotropic laser-induced tissue hyperthermia. Eng Anal Bound Elem 101: 156-164.
    [22] Fahmy MA (2020) A new convolution variational boundary element technique for design sensitivity analysis and topology optimization of anisotropic thermo-poroelastic structures. Arab J Basic Appl Scis 27: 1-12.
    [23] Verma R, Kumar S (2020) Computational study on constant and sinusoidal heating of skin tissue using radial basis functions. Comput Biol Med 121: 103808.
    [24] Hoffman JD, Hoffman JD, Frankel S (2001)  Numerical Methods for Engineers and Scientists CRC Press.
    [25] Kengne E, Mellal I, Hamouda MB, et al. (2014) A mathematical model to solve bio-heat transfer problems through a bio-heat transfer equation with quadratic temperature-dependent blood perfusion under a constant spatial heating on skin surface. J Biomed Sci Eng 7: 721-730.
    [26] Zhang ZW, Wang H, Qin QH (2015) Meshless method with operator splitting technique for transient nonlinear bioheat transfer in two-dimensional skin tissues. Int J Mol Sci 16: 2001-2019.
    [27] Woolfe G, MacDonald AD (1944) The evaluation of the analgesic action of pethidine hydrochloride (Demerol). J Pharmacol Exp Ther 80: 300-307.
    [28] Gholami M, Saboory E, Mehraban S, et al. (2015) Time dependent antinociceptive effects of morphine and tramadol in the hot plate test: Using different methods of drug administration in female rats. Iran J Pharm Res 14: 303-311.
    [29] Da Silva S, França AS, Pinatti M (2011) One-dimensional simulation of heat transfer in the canine knee joint during therapeutic heating and cooling. Braz J Biom Eng 27: 163-174.
    [30] Tita B, Abdel-Haq H, Vitalone A, et al. (2001) Analgesic properties of Epilobium angustifolium, evaluated by the hot plate test and the writhing test. Farmaco 56: 341-343.
    [31] Narasimhan A, Jha KK (2012) Bio-heat transfer simulation of retinal laser irradiation. Int J Numer Meth Biomed Eng 28: 547-559.
    [32] Ezzat MA, El-Bary AA, Al-Sowayan NS (2016) Tissue responses to fractional transient heating with sinusoidal heat flux condition on skin surface. Anim Sci J 87: 1304-1311.
    [33] Shih TC, Kou HS, Liauh CT, et al. (2002) Thermal models of bioheat transfer equations in living tissue and thermal dose equivalence due to hyperthermia. Biomed Eng: Appl, Basis, Commun 14: 86-96.
    [34] Zhang J, Chauhan S (2020) Fast computation of soft tissue thermal response under deformation based on fast explicit dynamics finite element algorithm for surgical simulation. Comput Meth Prog Bio 187: 105244.
    [35] Zhang J, Lay RJ, Roberts SK, et al. (2020) Towards real-time finite-strain anisotropic thermo-visco-elastodynamic analysis of soft tissues for thermal ablative therapy. Comput Meth Prog Bio 198: 105789.
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2201) PDF downloads(139) Cited by(0)

Article outline

Figures and Tables

Figures(10)  /  Tables(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog