Abstract
The aim of the article is to present a mathematical model that predicts tissue temperature during thermal therapy (thermal ablation). The mathematical formulation of variable thermal conductivity of thermal wave bioheat transfer model for a generalized coordinate system due to external heat source during thermal ablation has been given. The governing bioheat transfer equation is simplified by Kirchhoff’s transformation. The numerical solution of the problem has been achieved using a finite difference scheme to discretize in space coordinate; the reduced system of second order ordinary differential equation is solved by the finite element Legendre wavelet Galerkin method (FELWGM). The obtained result from FELWGM is compared with exact analytical solution and shows good agreement. Parametric study is performed to evaluate the effect of variable thermal conductivity on tissue temperature distribution for different physical parameters and illustrated graphically. Particular cases are also deduced from present investigation. The effect of variability of thermal conductivity parameter, variability of time, generalized coordinate system, lagging time and external heat source coefficient on tissue temperature is discussed in detail. Thus, the consideration of the variable thermal conductivity is extremely beneficial for the clinical therapeutic application in the treatment of cancerous cells.
1. Introduction
For image-guided thermal treatments delivering of heat via focused ultrasound (FUS), radiofrequency (RF), microwave (MW), or laser, temperature induced changes in tissue properties are of relevance in relation to predicting tissue temperature profile, monitoring during treatment, and evaluation of treatment results. These thermal treatments cause reversible or irreversible changes in the properties of tissues. Therefore, for thermal therapies the adequate consideration of tissue properties and their temperature dependence is necessary to achieve accurate results.
Mathematical modelling of thermal therapies has been used extensively to predict and optimize clinical treatments and medical devices. Pennes (Citation1948) model is widely used to model such problems. Kumar, Kumar, and Rai (Citation2015a) studied the heat transfer process in biological tissue during thermal ablation in a generalized coordinate system. Non-linear dual-phase-lag model is utilized by Kumar, Kumar, and Rai (Citation2016) in living tissue during thermal ablation. Kumar and Rai (Citation2017) studied the time fractional DPL model in presence of temperature dependent metabolic and electromagnetic radiation heat source. Kashcooli, Salimpour, and Shirani (Citation2017) investigated the effects of blood vessels on temperature distribution in a skin tissue subjected to various thermal therapy conditions.
Recently, Kumar, Damor, and Shukla, (Citation2018) studied the heat transfer thermal damage in triple layer skin tissue using a fractional bioheat model. For therapeutic treatment of cancerous cells, Kumar, Singh, Sharma, and Rai (Citation2018) used non-linear dual-phase-lag model under Dirichlet boundary conditions. Tan, Zou, Dong, Ding, and Zhao (Citation2018) investigated the coupled effects of blood vessels and relaxation time on tissue temperature and thermal lesion region during HIFU hyperthermia.
Very few authors considered the temperature-dependent thermal parameters over wide temperature ranges in tissues. Valvano, Cochran, and Diller (Citation1985) measured thermal conductivity and diffusivity of biomaterials over a wide temperature range. Bhattacharya and Mahajan (Citation2003) presented the temperature-dependent thermal conductivity of sheep collagen and cow liver in vitro. Watanabe, Kobayashi, Hashiszume, and Fujie (Citation2009) investigated the effects of the temperature dependence of thermophysical properties of the liver on the temperature distribution during RFA. Trujillo and Berjano (Citation2013) reviewed the mathematical functions most commonly used for modelling the temperature dependence of electrical and thermal conductivities of biological tissue in RFA.
Guntur, Lee, Paeng, Coleman, and Choi (Citation2013) found that the thermal properties of ex vivo porcine liver tissue during RFA vary significantly with temperature. Rossmann and Haemmerich (Citation2014) reviewed temperature-dependent thermal and electrical properties of biological tissues in the supraphysiological temperature regime. Guntur and Choi (Citation2015) examined the effects of temperature-dependent thermal parameters on the temperature distribution in liver tissue exposed to a clinical high-intensity focused ultrasound (HIFU) field, and found that alteration of tissue thermal parameters significantly affects the prediction of temperature distributions and suggest that the dependence of thermal parameters should be considered in estimation of the thermal dose in tissues exposed to a clinical HIFU field.
Very useful work has been carried out by Ellahi et al. on thermal conductivity; notable among them are Bhatti, Zeeshan, and Ellahi (Citation2017), Bhatti, Zeeshan, Ellahi, and Shit (Citation2018), Bhatti, Zeeshan, Tripathi, and Ellahi (Citation2018), Ellahi, Alamri, Sultan, Basit, and Majeed, (Citation2018), Ellahi, Zeeshan, Shehzad, and Alamri (Citation2018), Hassan, Marin, Ellahi, and Alamri (Citation2018), Ijaz, Zeeshan, Bhatti, and Ellahi (Citation2018). Akbar, Raza, and Ellahi (Citation2016a) studied the effects of entropy and magnetic field in an endoscope for a thermal conductivity model. Akbar, Raza, and Ellahi (Citation2016b) studied the anti-bacterial applications for new thermal conductivity model in arteries with CNT suspended nanofluid.
The present study aims to investigate the effects of variable thermal conductivity during thermal ablation for a thermal wave bioheat transfer model with modified Gaussian external heat source in a generalized coordinate system. The governing bioheat transfer equation is simplified by Kirchhoff’s transformation and then finite difference scheme, and Legendre wavelet Galerkin approach is used to solve the problem. The problem is converted into a system of second order differential equation with initial conditions by using the finite difference scheme in space coordinate. Then, this reduced system is solved by the Wavelet Galerkin approach with Legendre wavelet as basis function, which reduces the problem into the Sylvester matrix equation. The solution of the Sylvester matrix equation gives the temperature distribution during thermal ablation with variable thermal conductivity and a significant effect of the variable conductivity has been observed on tissue temperature distribution.
2. Formulation of the problem
Consider a tissue of finite length (L), which is heated by modified Gaussian external heat source with variable thermal conductivity during thermal ablation in generalized coordinate system. Pennes (Citation1948) performed a series of experiments that measured temperatures on human forearms and derived a thermal energy conservation equation known as bioheat transfer equation as follows:
From the principle of conservation of energy to tissue control volume (V), we have (1) (1)
where Ugain, Ustorage, and Uloss are given below.
The general internal heat generation per unit volume due to metabolic heat (qm) and the interior heat caused by outer heat source (qs) is (2) (2)
The total rate of stored thermal energy over the control volume is (3) (3)
Heat loss to adjacent tissues by convection and conduction is given as (4) (4)
where (5) (5)
is the rate production through the control volume and (6) (6)
the total amount of thermal energy over the tissue control volume.
Using EquationEquations (2)–(6) in EquationEquation (1)(1) (1) , yields (7) (7)
Applying the divergence theorem to the surface integral and observing that EquationEquation (7)(7) (7) must hold for any arbitrary volume element, we obtain (8) (8)
Cattaneo (Citation1958), Vernotte (Citation1958) proposed the following modified heat flux model that is a linear extension of the Fourier’s law in order to take in the wave-like behaviour in heat conduction process, (9) (9)
where k(T) is the thermal conductivity with variable temperature.
Using EquationEquation (9)(9) (9) in EquationEquation (8)(8) (8) with temperature dependent specific heat yields (10) (10)
EquationEquation (10)(10) (10) in generalized coordinate system can be written as (11) (11)
where corresponds to Cartesian, axis-symmetric, and spherical symmetric coordinate system and (12) (12)
The term qb can be expressed as (13) (13)
The term qm can be approximated as a linear function of local tissue temperature taken in Gupta, Singh, and Rai (Citation2010), (14) (14)
Following (Kumar, Kumar, & Rai, Citation2015b) modified Gaussian heat source qs is taken as (15) (15)
To solve the problem, the following initial and boundary conditions are taken:
Initial conditions (16) (16)
Boundary conditions (17) (17)
Symmetric condition (18) (18)
Considering the mapping (Kirchhoff’s transformation), following Hetnarski (Citation1986) as (19) (19)
yields (20) (20)
An important factor to achieve realistic models is the use of mathematical functions to describe the temperature dependence of thermal properties of tissue. A commonly used approach for modelling the temperature dependence of thermal properties for temperature below is based on linear equations and employs constant temperature coefficient such as: (21) (21)
From the above, we have the following relations: (22) (22)
and (23) (23)
EquationEquations (11)(11) (11) and Equation(16–18) with the aid of EquationEquations (19–21) in linear form become (24) (24)
(25) (25) (26) (26)3. Solution of the problem
Introducing the non-dimensional variable and similarity criteria (27) (27)
With the help of EquationEquation (27)(27) (27) , EquationEquations (24–26) reduce to the following form (28) (28)
(29) (29) (30) (30) (31) (31)where .
3.1. Finite difference scheme for spatial-discretization
The space coordinate domain [0,1] divided into n + 1 sub intervals of equal step length (h) i.e. , where . The first and second order derivatives are approximated by using central difference formula i.e. (32) (32)
(33) (33)Using this spatial discretization scheme, EquationEquation (28)(28) (28) becomes (34) (34)
, with initial condition (35) (35)
boundary condition (36) (36)
The EquationEquation (34)(34) (34) is converted into a system of second order ordinary differential equation in vector matrix form as follows: (37) (37)
subjected to initial condition (38) (38)
where , I is the identity matrix of order n, A is a square matrix of order ; B, , and are the column matrix of which are given as follows
where
and
3.2. Legendre wavelet Galerkin approach
To solve the differential EquationEquation (37)(37) (37) , let us assume that (39) (39)
where
and (40) (40)
and CT is unknown matrix, and is a matrix defined as
and
The Legendre wavelet , k is any positive integer, defined on the interval [0,1) by Razzaghi, and Yousefi (Citation2001) (41) (41)
where is the well known Legendre polynomials of degree m which are given as follows:
and
Integrating EquationEquation (39)(39) (39) with respect to F0 from 0 to F0, we have (42) (42)
where P is operational matrix of integration (Razzaghi, & Yousefi, Citation2001) given as follows (43) (43)
where L and F are M × M matrices given by: (44) (44)
and (45) (45)
Again integrating EquationEquation (42)(42) (42) with respect to F0 from 0 to F0, (46) (46)
where .
Using EquationEquations (39)(39) (39) , Equation(42)(42) (42) , and Equation(46)(46) (46) in EquationEquation (37)(37) (37) , we have (47) (47)
where
In is the identity matrix of order n and d1 is matrix such that .
EquationEquation (47)(47) (47) is Sylvester matrix equation, after solving this equation, we obtain the element of matrix CT. Substituting the value of CT in EquationEquation (46)(46) (46) gives tissue temperature distribution . The tissue temperature in terms of T can be obtained from EquationEquation (23)(23) (23) .
Particular cases:
If we take in EquationEquation (47)(47) (47) , the problem is reduced to the non-Fourier bioheat transfer model during thermal ablation and the corresponding results are the same as obtained by Kumar, Kumar, and Rai (Citation2015a).
If we take in EquationEquation (47)(47) (47) , then we obtain results for the Pennes’ bioheat transfer model with variable thermal conductivity during thermal ablation.
If we take and in EquationEquation (47)(47) (47) , then we obtain results for the Pennes’ bioheat transfer model in the absence of variable thermal conductivity during thermal ablation.
4. Numerical results and discussion
Yadav, Kumar, and Rai (Citation2014) showed that the Legendre wavelet Galerkin method produces a higher degree of accuracy and minimizes error. That is why we used this approach for numerical procedure. This approach combines both the multi-resolution and multi-scale computational property of Legendre wavelet with element-wise analysis. Multi-resolution analysis of Legendre wavelet localizes small scale variation of solution and fast switching of functional bases. A brief algorithm of the numerical method is given in .
The problem is solved analytically by Laplace transform in the absence of modified Gaussian external heat source, which is given as follows:
By taking the Laplace transform of EquationEquation (28)(28) (28) and applying the initial conditions (29), the following ordinary differential equation is obtained (48) (48)
where (49) (49)
The exact solution of EquationEquation (48)(48) (48) , by using the Laplace transform of boundary conditions (30–31), is (50) (50)
The inverse Laplace transform of can be obtained from the following Bromwich contour integration (Arpaci, Citation1996) (51) (51)
Using the inversion theorem, the inverse Laplace transform of EquationEquation (50)(50) (50) is (52) (52)
By using Bromwich contour integration (Arpaci, Citation1996), temperature distribution in tissue is (53) (53)
where , m = 1, 2, . The approximate solution (FELWGM) has compared with analytical solution given by EquationEquation (53)(53) (53) . shows good agreement between the results obtained by the numerical and exact solution.
The value of thermo-physical parameters is similar to the values of parameters as in the literature (Kumar, Kumar, & Rai, Citation2015a).
MATLAB(R2016a) software has been used for all computational work. Thermal ablation refers to the destruction of tissue by extreme hyperthermia (elevated tissue temperatures) or hypothermia (depressed tissue temperatures). The temperature change is concentrated to a focal zone in and around the tumour. The high temperature occurs at the target region for four to six minutes during the process.
shows the effect of variable thermal conductivity parameter K1 on tissue temperature distribution. It is noted that as the value of K1 increases, the temperature increases and the temperature attains maxima at the heating position . Thus, to achieve more accuracy in clinical treatments variable thermal conductivity plays an important role.
shows the effect of time on tissue temperature along the space coordinate when . It is evident that temperature increases as the time increases from four minutes to six minutes during thermal ablation. Thus, it is ned that time duration affects the treatment process.
shows the effect of a generalized coordinate system on the treatment of tumour during thermal ablation for and . It is noted that for all three coordinate systems, tissue temperature remains unaffected at the heating position for both the cases, but for and , it affects the tissue temperature.
shows the comparison of temperature distribution for Pennes () and thermal wave model . It is observed that temperature profile for the thermal wave model is less than the temperature profile predicted by the Pennes’ model for both the cases of and . But, independently, both models predict higher temperature with the variable thermal conductivity parameter. Thus, lagging time is affecting the temperature distribution during the treatment.
shows the effects of external heat source on temperature distribution for and . It is clear that as the value of K1 and increases, temperature also increases. Thus, the consideration of external heat source is important during the treatment. depicts the absolute error of the solutions obtained by finite element Legendre wavelet Galerkin method (FELWGM) and analytical method. Temperature (T) for time (t) 0 to 500 s is evaluated by both the methods and results are given in .
Therefore, we can say that the consideration of variable thermal conductivity in bioheat transfer model is very beneficial for the clinical therapeutic application in the treatment of cancerous cells.
5. Conclusions
The problem of modelling the variable thermal conductivity of thermal wave bioheat transfer model in a generalized coordinate system due to external heat source during thermal ablation has been investigated. Kirchhoff’s transformation is used to formulate the governing bioheat transfer equation. The numerical solution of the problem is achieved by using a finite difference scheme and Legendre wavelet Galerkin approach. FELWGM converts the problem into a system of algebraic equation, so it reduces the computational complexity and avoids large storage requirements, as in the case of conventional discretize schemes. The results obtained from the FELWGM are compared with the analytical results and give good accuracy in the absence of external heat source. The significant effect of variable thermal conductivity parameters, variability of time, lagging times, and external heat source is observed on tissue temperature distribution. Therefore, the temperature dependent changes need to be considered to enhance the accuracy of predicting and monitoring pre-clinical and clinical treatments and consideration of this model is useful in image-guided thermal treatments (e.g. hyperthermia and thermal ablation) and other biomedical treatments. It will give new direction in clinical therapeutic application.
Nomenclature | ||
Γ | = | the number of classify coordinate |
κ | = | tissue thermal diffusivity, m2s−1 |
ωb | = | blood perfusion rate, s−1 |
ρ | = | tissue density, kgm−3 |
ρb | = | blood mass density, kgm−3 |
τq | = | heat flux relaxation time, s |
Θ | = | dimensionless tissue temperature, in EquationEquation (28)(28) (28) |
θ | = | dimensionless tissue temperature |
= | dimensionless arterial blood temperature | |
= | dimensionless wall temperature at boundary | |
a0 | = | scattering parameter, m−1 |
c | = | specific heat of tissue, jkg−1 |
cb | = | specific heat of blood, jkg−1k−1 |
= | dimensionless heat flux relaxation time | |
F0 | = | dimensionless time |
k | = | variable thermal conductivity of tissue, Wm−1K−1 |
k0 | = | constant thermal conductivity of tissue, Wm−1K−1 |
K1 | = | dimensionless variable thermal conductivity parameter |
k1 | = | variable thermal conductivity parameter, °C−1 |
L | = | tissue slab length, m |
l | = | Bromwich contour integration line |
P | = | transmitted power, W |
Pf | = | dimensionless blood perfusion coefficient |
= | dimensionless metabolic heat source coefficient | |
= | dimensionless external heat source coefficient | |
qb | = | heat source due to blood perfusion, Wm−3 |
qm | = | heat source due to metabolic heat generation, Wm−3 |
qs | = | external heat source, Wm−3 |
= | constant metabolic heat source | |
= | reference external heat source, Wm−3 | |
r | = | space coordinate, m |
r0 | = | location of tumor, m |
S | = | antenna power, m−1 |
s | = | Laplace domain parameter |
T | = | tissue temperature, °C |
t | = | time, s |
T0 | = | initial temperature of body, °C |
Ta | = | arterial blood temperature, °C |
Tw | = | wall temperature at the boundary, °C |
x | = | dimensionless space coordinate |
= | dimensionless location of tumor where heat is applied |
Acknowledgement
The third author is thankful to the Council of Scientific and Industrial Research (CSIR), India, for providing fellowship under the SRF(21/06/2015/(i)EU-V) scheme.
Disclosure statement
The authors have no conflict of interest.
References
- Akbar, N. S., Raza, M., & Ellahi, R. (2016a). Anti-bacterial application for new thermal conductivity model in arteries with CNT suspended nanofluid. Journal of Mechanics in Medicine and Biology, 16(5), 1650063. doi: 10.1142/S0219519416500639
- Akbar, N. S., Raza, M., & Ellahi, R. (2016b). Endoscopic effects with entropy generation analysis in peristalsis for the thermal conductivity of H2O+Cu nanofluid. Journal of Applied Mathematics and Mechanics, 9(4), 1721–1730.
- Arpaci, V. S. (1996). Conduction heat transfer. Reading, MA: Addison-Wesley Publishing Co.
- Bhattacharya, A., & Mahajan, R. L. (2003). Temperature dependence of thermal conductivity of biological tissues. Physiological Measurement, 24(3), 769–783. doi: 10.1088/0967-3334/24/3/312
- Bhatti, M. M., Zeeshan, A., & Ellahi, R. (2017). Electromagnetohydrodynamic (EMHD) peristaltic flow of solid particles in a third-grade fluid with heat transfer. Mechanics & Industry, 18(3), 314. doi: 10.1051/meca/2016061
- Bhatti, M. M., Zeeshan, A., Ellahi, R., & Shit, G. C. (2018). Mathematical modeling of heat and mass transfer effects on MHD peristaltic propulsion of two-phase flow through a Darcy–Brinkman–Forchheimer porous medium. Advanced Powder Technology, 29(5), 1189–1197. doi: 10.1016/j.apt.2018.02.010
- Bhatti, M. M., Zeeshan, A., Tripathi, D., & Ellahi, R. (2018). Thermally developed peristaltic propulsion of magnetic solid particles in biorheological fluids. Indian Journal of Physics, 92(4), 423–430. doi: 10.1007/s12648-017-1132-x
- Cattaneo, C. (1958). A form of heat conduction equation which eliminates the paradox of instantaneous propagation. Compte Rendus, 247, 431–433.
- Ellahi, R., Alamri, Sultan, Z., Basit, A., & Majeed, A. (2018). Effects of MHD and slip on heat transfer boundary layer flow over a moving plate based on specific entropy generation. Journal of Taibah University for Science, 12(4), 476–482. doi: 10.1080/16583655.2018.1483795
- Ellahi, R., Zeeshan, A., Shehzad, N., & Alamri, S. Z. (2018). Structural impact of kerosene-Al2O3 nanoliquid on MHD Poiseuille flow with variable thermal conductivity: Application of cooling process. Journal of Molecular Liquids, 264, 607–615. doi: 10.1016/j.molliq.2018.05.103
- Guntur, S. R., & Choi, M. J. (2015). Influence of temperature-dependent thermal parameters on temperature elevation of tissue exposed to high-intensity focused ultrasound: Numerical simulation. Ultrasound in Medicine and Biology, 41(3), 806–813. doi: 10.1016/j.ultrasmedbio.2014.10.008
- Guntur, S. R., Lee, K. I., Paeng, D. G., Coleman, A. J., & Choi, M. J. (2013). Temperature-dependent thermal properties of ex vivo liver undergoing thermal ablation. Ultrasound in Medicine and Biology, 39, 1–14.
- Gupta, P. K., Singh, J., & Rai, K. N. (2010). Numerical simulation for heat transfer in tissues during thermal therapy. Journal of Thermal Biology, 35(6), 295–301. doi: 10.1016/j.jtherbio.2010.06.007
- Hassan, M., Marin, M., Ellahi, R., & Alamri, S. Z. (2018). Exploration of convective heat transfer and flow characteristics synthesis by Cu–Ag water hybrid-nanofluids. Heat Transfer Research, 49(18), 1837–1848. doi: 10.1615/HeatTransRes.2018025569
- Hetnarski, R. B. (1986). Thermal Stresses, Vol. I. North-Holland, Amsterdam: Taylor and Francis Group.
- Ijaz, N., Zeeshan, A., Bhatti, M. M., & Ellahi, R. (2018). Analytical study on liquid-solid particles interaction in the presence of heat and mass transfer through a wavy channel. Journal of Molecular Liquids, 250, 80–87. doi: 10.1016/j.molliq.2017.11.123
- Kashcooli, M., Salimpour, M. R., & Shirani, E. (2017). Heat transfer analysis of skin during thermal therapy using thermal wave equation. Journal of Thermal Biology, 64, 7–18. doi: 10.1016/j.jtherbio.2016.12.007
- Kumar, D., & Rai, K. N. (2017). Numerical simulation of time fractional dual-phase-lag model of heat transfer within skin tissue during thermal therapy. Journal of Thermal Biology, 67, 49–58. doi: 10.1016/j.jtherbio.2017.05.001
- Kumar, D., Kumar, P., & Rai, K. N. (2015a). Numerical Study on non-Fourier bioheat transfer during thermal ablation. Procedia Engineering, 127, 1300–1307. doi: 10.1016/j.proeng.2015.11.487
- Kumar, D., Singh, S., Sharma, N., & Rai, K. N. (2018). Verified non-linear DPL model with experimental data for analyzing heat transfer in tissue during thermal therapy. International Journal of Thermal Sciences, 133, 320–329. doi: 10.1016/j.ijthermalsci.2018.07.031
- Kumar, P., Kumar, D., & Rai, K. N. (2015b). A numerical study on dual-phase-lag model of bio-heat transfer during hyperthermia treatment. Journal of Thermal Biology, 49–50, 98–105. doi: 10.1016/j.jtherbio.2015.02.008
- Kumar, P., Kumar, D., & Rai, K. N. (2016). Non-linear dual-phase-lag model for analyzing heat transfer phenomena in living tissues during thermal ablation. Journal of Thermal Biology, 20, 204–212. doi: 10.1016/j.jtherbio.2016.07.017
- Kumar, S., Damor, R. S., & Shukla, A. K. (2018). Numerical study on thermal therapy of triple layer skin tissue using fractional bioheat model. International Journal of Biomathematics, 11(04), 1850052. doi: 10.1142/S1793524518500523
- Pennes, H. H. (1948). Analysis of tissue and arterial blood temperature in the resting forearm. Journal of Applied Physiology, 1(2), 93–122. doi: 10.1152/jappl.1948.1.2.93
- Razzaghi, M., & Yousefi, S. (2001). The Legendre wavelets operational matrix of integration. International Journal of Systems Science, 32(4), 495–502. doi: 10.1080/002077201300080910
- Rossmann, C., & Haemmerich, D. (2014). Review of temperature dependence of thermal properties, dielectric properties, and perfusion of biological tissues at hyperthermic and ablation temperatures. Critical Reviews in Biomedical Engineering, 42, 467–492. doi: 10.1615/CritRevBiomedEng.2015012486
- Tan, Q., Zou, X., Dong, H., Ding, Y., & Zhao, X. (2018). Influence of blood vessels on temperature during high-intensity focused ultrasound hyperthermia based on the thermal wave model of bioheat transfer. Advances in Condensed Matter Physics, 2018, 5018460. doi: 10.1155/2018/5018460
- Trujillo, M., & Berjano, E. (2013). Review of the mathematical functions used to model the temperature dependence of electrical and thermal conductivities of biological tissues in radiofrequency ablation. International Journal of Hyperthermia, 29(6), 590–597. doi: 10.3109/02656736.2013.807438
- Valvano, J. W., Cochran, J. R., & Diller, K. R. (1985). Thermal conductivity and diffusivity of biomaterials measured with self-heated thermistors. International Journal of Thermophysics, 6(3), 301–311. doi: 10.1007/BF00522151
- Vernotte, P. (1958). Les paradoxes de la theorie continue de l’equation de la chaleur. Compte Rendus, 246, 3154–3155.
- Watanabe, H., Kobayashi, Y., Hashiszume, M., & Fujie, M. G. (2009). Modeling the temperature dependence of thermophysical properties: Study on the effect of temperature dependence for RFA, 31st Annual International Conference of the IEEE EMBS Minneapolis, Minnesota, USA, September 2–6.
- Yadav, S., Kumar, D., & Rai, K. N. (2014). Finite element Legendre wavelet Galerkin approach to inward solidification in simple body under most generalized boundary condition. Zeitschrift für Naturforschung, 69a, 501–510.