RESUMEN
In this study, we have established a space-time fractional DPL bioheat transfer model in the presence of temperature-dependent metabolic and space-time dependent electromagnetic heat sources. Applying the Legendre wavelet collocation method, the fractional order partial differential equation is reduced into the system of algebraic equations, which has been solved using the Newton iteration method. The error bound as well as stability analysis and numerical scheme validation are provided. The time to achieve for the position of hyperthermia is discussed in three cases: the DPL model, the time-fractional DPL model, and the space-time-fractional DPL model. The effect of variability of time and space fractional derivative orders (α and ß), transmitted power (P) and lagging times on the temperature profile in biological tissue at a different time are discussed in detail. We conclude that a suitable value of α, ß, τT, τq, and P provides a desirable temperature at a particular time in thermal therapies. Such knowledge will be very useful in the clinical therapeutic application.
Asunto(s)
Temperatura Corporal , Hipertermia Inducida/métodos , Modelos Teóricos , Conductividad Térmica , Humanos , Neoplasias/terapia , Especificidad de Órganos , Análisis Espacio-TemporalRESUMEN
Hyperthermia is a process that uses heat from the spatial heat source to kill cancerous cells without damaging the surrounding healthy tissues. Efficacy of hyperthermia technique is related to achieve temperature at the infected cells during the treatment process. A mathematical model on heat transfer in multilayer tissues in finite domain is proposed to predict the control temperature profile at hyperthermia position. The treatment technique uses dual-phase-lag model of heat transfer in multilayer tissues with modified Gaussian distribution heat source subjected to the most generalized boundary condition and interface at the adjacent layers. The complete dual-phase-lag model of bioheat transfer is solved using finite element Legendre wavelet Galerkin approach. The present solution has been verified with exact solution in a specific case and provides a good accuracy. The effect of the variability of different parameters such as lagging times, external heat source, metabolic heat source and the most generalized boundary condition on temperature profile in multilayer tissues is analyzed and also discussed the effective approach of hyperthermia treatment. Furthermore, we studied the modified thermal damage model with regeneration of healthy tissues as well. For viewpoint of thermal damage, the least thermal damage has been observed in boundary condition of second kind. The article concludes with a discussion of better opportunities for future clinical application of hyperthermia treatment.
Asunto(s)
Calor/efectos adversos , Hipertermia Inducida/efectos adversos , Modelos Biológicos , Neoplasias/fisiopatología , Neoplasias/terapia , Conductividad Térmica , Análisis de Elementos Finitos , Humanos , Distribución Normal , Regeneración , Análisis de OndículasRESUMEN
This paper theoretically investigates the thermal behavior in a living biological tissue under various coordinate systems and different non-Fourier boundary conditions with the dual-phase-lag bioheat transfer model during thermal therapy. The properties of Legendre wavelets together with the finite difference scheme are used to find an approximate analytical solution of the present problem. It has been observed that surrounding healthy tissues are less affected in second and third kind of boundary condition when applied along with spherical symmetric coordinate system. Also greater temperature rise and fast achievement of peak hyperthermia temperature is achieved when second and third kind of boundary conditions are used in combination with Cartesian coordinate system. It is observed that due to the presence of blood perfusion and temperature dependent metabolic heat generation term, the dual-phase-lag bioheat transfer model reduces to Pennes bioheat transfer model only when τq=τT=0s, not for arbitrary τq=τT. Further, in case of dual-phase-lag bioheat transfer model wave-like or diffusion-like behavior will dominate depends whether the ratio τq/τT > 1 or τq/τT < 1. Effect of temperature dependent metabolic heat generation rate, thermal conductivity and blood perfusion rate on dimensionless temperature are discussed in details. The whole analysis is presented in dimensionless form.
Asunto(s)
Transferencia de Energía , Hipertermia Inducida , Modelos Biológicos , HumanosRESUMEN
In this article, a non-linear dual-phase-lag (DPL) bio-heat transfer model based on temperature dependent metabolic heat generation rate is derived to analyze the heat transfer phenomena in living tissues during thermal ablation treatment. The numerical solution of the present non-linear problem has been done by finite element Runge-Kutta (4,5) method which combines the essence of Runge-Kutta (4,5) method together with finite difference scheme. Our study demonstrates that at the thermal ablation position temperature predicted by non-linear and linear DPL models show significant differences. A comparison has been made among non-linear DPL, thermal wave and Pennes model and it has been found that non-linear DPL and thermal wave bio-heat model show almost same nature whereas non-linear Pennes model shows significantly different temperature profile at the initial stage of thermal ablation treatment. The effect of Fourier number and Vernotte number (relaxation Fourier number) on temperature profile in presence and absence of externally applied heat source has been studied in detail and it has been observed that the presence of externally applied heat source term highly affects the efficiency of thermal treatment method.
Asunto(s)
Simulación por Computador , Análisis de Elementos Finitos , Hipertermia Inducida , Modelos Biológicos , Neoplasias/terapia , Conductividad Térmica , Regulación de la Temperatura Corporal , Calor , Humanos , Hipertermia Inducida/métodosRESUMEN
In this paper, dual-phase-lag bioheat transfer model subjected to Fourier and non-Fourier boundary conditions for bi-layer tissues has been solved using finite element Legendre wavelet Galerkin method (FELWGM) during magnetic fluid hyperthermia. FELWGM localizes small scale variation of solution and fast switching of functional bases. It has been observed that moderate hyperthermia temperature range (41-46°C) can be better achieved in spherical symmetric coordinate system and treatment method will be independent of the Fourier and non-Fourier boundary conditions used. The effect of phase-lag times has been observed only in tumor region. FCC FePt magnetic nano-particle produces more effective treatment with respect to other magnetic nano-particles. The effect of variability of magnetic heat source parameters (magnetic induction, frequency, diameter of magnetic nano-particles, volume fractional of magnetic nano-particles and ligand layer thickness) has been investigated. The physical property of these parameters has been described in detail during magnetic fluid hyperthermia (MFH) treatment and also discussed the clinical application of MFH in Oncology.
Asunto(s)
Calor , Hipertermia Inducida , Membrana Dobles de Lípidos/química , Campos Magnéticos , Nanopartículas de Magnetita/química , Modelos QuímicosRESUMEN
The success of hyperthermia in the treatment of cancer depends on the precise prediction and control of temperature. It was absolutely a necessity for hyperthermia treatment planning to understand the temperature distribution within living biological tissues. In this paper, dual-phase-lag model of bio-heat transfer has been studied using Gaussian distribution source term under most generalized boundary condition during hyperthermia treatment. An approximate analytical solution of the present problem has been done by Finite element wavelet Galerkin method which uses Legendre wavelet as a basis function. Multi-resolution analysis of Legendre wavelet in the present case localizes small scale variations of solution and fast switching of functional bases. The whole analysis is presented in dimensionless form. The dual-phase-lag model of bio-heat transfer has compared with Pennes and Thermal wave model of bio-heat transfer and it has been found that large differences in the temperature at the hyperthermia position and time to achieve the hyperthermia temperature exist, when we increase the value of τT. Particular cases when surface subjected to boundary condition of 1st, 2nd and 3rd kind are discussed in detail. The use of dual-phase-lag model of bio-heat transfer and finite element wavelet Galerkin method as a solution method helps in precise prediction of temperature. Gaussian distribution source term helps in control of temperature during hyperthermia treatment. So, it makes this study more useful for clinical applications.