Modeling the heating of biological tissue based on the hyperbolic heat transfer equation

Therapeutic processes for eradicating cancerous or benign tumours by laser beams using the excitonic approach of peptide groups

The aim of the present study was to develop a protocol for the treatment of cancerous or benign tumours making use of laser rays, also demonstrating that the destruction process remains exclusively confined in the defective organ. Thermal effects of lasers on biological tissue have been elucidated using vibrational excitations approach of peptide groups (PGs). It was proposed a Hamiltonian which integrate excitations induced by laser pulses and it was shown that the system is governed by a nonlinear equation with strong nonlinearity. It was also exactly described what happens in polypeptide chain once the unwanted organ is irradiated by the Neodymium-doped yttrium aluminium garnet, chosen as incident laser. It was shown that, the advent of incident laser beams contributes to a sudden reinforcement of the vibrational excitations of PGs frequencies and amplitudes. It was also demonstrated that the heating process leads to transverse and longitudinal deformation of the polypeptide chain and these sudden changes lead to the denaturation and subsequently to the destruction of the bulky organ. The drawn curves make it possible to estimate the spatial expansion of the denaturation, in order to effectively control the spread of the heat. Laser irradiation leads to a drastic increase in the vibration amplitudes of the PGs and subsequently results in the destruction of the undesirable tissue. An appropriate choice of the laser can make it possible to circumscribe the destruction only in the defective zone and to protect healthy cells.

The convective instability of a Maxwell-Cattaneo fluid in the presence of a vertical magnetic field

We study the instability of a Bénard layer subject to a vertical uniform magnetic field, in which the fluid obeys the Maxwell-Cattaneo (MC) heat flux-temperature relation. We extend the work of Bissell (Proc. R. Soc. A 472, 20160649 (doi:10.1098/rspa.2016.0649)) to non-zero values of the magnetic Prandtl number p m . With non-zero p m , the order of the dispersion relation is increased, leading to considerably richer behaviour. An asymptotic analysis at large values of the Chandrasekhar number Q confirms that the MC effect becomes important when C Q 1/2 is O(1), where C is the MC number. In this regime, we derive a scaled system that is independent of Q. When CQ 1/2 is large, the results are consistent with those derived from the governing equations in the limit of Prandtl number p → ∞ with p m finite; here we identify a new mode of instability, which is due neither to inertial nor induction effects. In the large p m regime, we show how a transition can occur between oscillatory modes of different horizontal scale. For Q ≫ 1 and small values of p, we show that the critical Rayleigh number is non-monotonic in p provided that C > 1/6. While the analysis of this paper is performed for stress-free boundaries, it can be shown that other types of mechanical boundary conditions give the same leading-order results.

Thermal field and tissue damage analysis of moving laser in cancer thermal therapy

In this paper, a closed-form analytical solution of hyperbolic Pennes bioheat equation is obtained for spatial evolution of temperature distributions during moving laser thermotherapy of the skin and kidney tissues. The three-dimensional cubic homogeneous perfused biological tissue is adopted as a media and the Gaussian distributed function in surface and exponentially distributed in depth is used for modeling of laser moving heat source. The solution procedure is Eigen value method which leads to a closed form solution. The effect of moving velocity, perfusion rate, laser intensity, absorption and scattering coefficients, and thermal relaxation time on temperature profiles and tissue thermal damage are investigated. Results are illustrated that the moving velocity and the perfusion rate of the tissues are the main important parameters in produced temperatures under moving heat source. The higher perfusion rate of kidney compared with skin may lead to lower induced temperature amplitude in moving path of laser due to the convective role of the perfusion term. Furthermore, the analytical solution can be a powerful tool for analysis and optimization of practical treatment in the clinical setting and laser procedure therapeutic applications and can be used for verification of other numerical heating models.

Two-dimensional closed-form model for temperature in living tissues for hyperthermia treatments

This research article determines an exact analytical expression for 2-D thermal field in single layer living tissues under a therapeutic condition by means of Fourier and non-Fourier heat transfer approaches. An actual spatially dependent initial condition has been adopted to analyze the heat propagation in tissues. The exact analytical determination for this actual initial condition for temperature may be difficult. However, in this study, an approximate analytical method has newly been established for an appropriate initial condition. With this initial expression, an exact temperature distribution for 2-D heat conduction in plane co-ordinates has been investigated for the predefined therapeutic boundary condition to have knowledge for practical aspects of the thermal therapy. Laplace Transform Method (LTM) in conjunction with the Inversion Theorem is used for the analytical solution treatment. We have utilized both Pennes' bioheat equation (PBHE) and thermal wave model of bioheat equation (TWMBHE) for the analysis. The influence of thermo-biological behavior on 2-D heat conduction in tissues has been studied with the variation of several dependable parameters in relation to the Hyperthermia treatment protocol in a moderate temperature range (42-45°C). The result in the present study has been evidenced for the biological heat transfer for the enforcement of different circumstances and also has been validated with the published value where the maximum temperature deviation of 2.6% has been recorded. We conclude that the temperature curve for TWMBHE model shows a higher waveform nature for low thermal relaxation time and this wavy nature gradually diminishes with an increase in relaxation time. The maximum peak temperature attains 46.3°C for the relaxation time = 2s and with the increase in the relaxation time the peak temperature gradually falls. The impact of blood perfusion rate on the relaxation time has also been established in this paper.

Nonlinear analysis of a non-Fourier heat conduction problem in a living tissue heated by laser source

The effect of laser, as a heat source, on a one-dimensional finite living tissue was studied in this paper. The dual phase lagging (DPL) non-Fourier heat conduction model was used for thermal analysis. The thermal conductivity was assumed temperature-dependent, resulting in a nonlinear equation. The obtained equations were solved using the approximate-analytical Adomian decomposition method (ADM). It was concluded that the nonlinear analysis was important in non-Fourier heat conduction problems. Moreover, a good agreement between the present nonlinear model and experimental result was obtained.

Convective instabilities of Maxwell–Cattaneo fluids

Motivated by the need to understand better the dynamics of non-Fourier fluids, we examine the linear and weakly nonlinear stabilities of a horizontal layer of fluid obeying the Maxwell–Cattaneo relationship of heat flux and temperature using three different forms of the time derivative of the heat flux. Linear stability mode regime diagrams in the parameter plane have been established and used to summarize the linear instabilities. The energy balance of the system is used to identify the mechanism by which the Maxwell–Cattaneo effect (i) introduces overstability, (ii) leads to preferred stationary modes with the critical Rayleigh and wavelengths either both increasing or both decreasing, (iii) gives rise to instabilities in a layer heated from above, and (iv) enhances heat transfer. A formal weakly nonlinear analysis leads to evolution equations for the amplitudes of linear instability modes. It is shown that the amplitude of the stationary mode obeys an equation of the Landau–Stuart type. The two equally excitable overstable modes obey two equations of the same type coupled by an interaction term. The evolution of the different amplitudes leads to supercritical stability, supercritical instability or subcritical instability depending on the model and parameter values. The results are presented in regime diagrams.

Thermal convection in a magnetized conducting fluid with the Cattaneo-Christov heat-flow model

By substituting the Cattaneo-Christov heat-flow model for the more usual parabolic Fourier law, we consider the impact of hyperbolic heat-flow effects on thermal convection in the classic problem of a magnetized conducting fluid layer heated from below. For stationary convection, the system is equivalent to that studied by Chandrasekhar (Hydrodynamic and Hydromagnetic Stability, 1961), and with free boundary conditions we recover the classical critical Rayleigh number [Formula: see text] which exhibits inhibition of convection by the field according to [Formula: see text] as [Formula: see text], where Q is the Chandrasekhar number. However, for oscillatory convection we find that the critical Rayleigh number [Formula: see text] is given by a more complicated function of the thermal Prandtl number [Formula: see text], magnetic Prandtl number [Formula: see text] and Cattaneo number C. To elucidate features of this dependence, we neglect [Formula: see text] (in which case overstability would be classically forbidden), and thereby obtain an expression for the Rayleigh number that is far less strongly inhibited by the field, with limiting behaviour [Formula: see text], as [Formula: see text]. One consequence of this weaker dependence is that onset of instability occurs as overstability provided C exceeds a threshold value CT(Q); indeed, crucially we show that when Q is large, [Formula: see text], meaning that oscillatory modes are preferred even when C itself is small. Similar behaviour is demonstrated in the case of fixed boundaries by means of a novel numerical solution.

On oscillatory convection with the Cattaneo-Christov hyperbolic heat-flow model

Adoption of the hyperbolic Cattaneo-Christov heat-flow model in place of the more usual parabolic Fourier law is shown to raise the possibility of oscillatory convection in the classic Bénard problem of a Boussinesq fluid heated from below. By comparing the critical Rayleigh numbers for stationary and oscillatory convection, Rc and RS respectively, oscillatory convection is found to represent the preferred form of instability whenever the Cattaneo number C exceeds a threshold value CT≥8/27π2≈0.03. In the case of free boundaries, analytical approaches permit direct treatment of the role played by the Prandtl number [Formula: see text], which-in contrast to the classical stationary scenario-can impact on oscillatory modes significantly owing to the non-zero frequency of convection. Numerical investigation indicates that the behaviour found analytically for free boundaries applies in a qualitatively similar fashion for fixed boundaries, while the threshold Cattaneo number CT is computed as a function of [Formula: see text] for both boundary regimes.

Measurement of the anisotropic thermal conductivity of the porcine cornea

Accurate thermal models for the cornea of the eye support the development of thermal techniques for reshaping the cornea and other scientific purposes. Heat transfer in the cornea must be quantified accurately so that a thermal treatment does not destroy the endothelial layer, which cannot regenerate, and yet is responsible for maintaining corneal transparency. We developed a custom apparatus to measure the thermal conductivity of ex vivo porcine corneas perpendicular to the surface and applied a commercial apparatus to measure thermal conductivity parallel to the surface. We found that corneal thermal conductivity is 14% anisotropic at the normal state of corneal hydration. Small numbers of ex vivo feline and human corneas had a thermal conductivity perpendicular to the surface that was indistinguishable from the porcine corneas. Aqueous humor from ex vivo porcine, feline, and human eyes had a thermal conductivity nearly equal to that of water. Including the anisotropy of corneal thermal conductivity will improve the predictive power of thermal models of the eye.

Tipping points in Cattaneo–Christov thermohaline convection

A model is proposed for thermohaline convection when the heat flux satisfies a relaxation time law rather than the classical Fourier one. Here, we adopt the recent law due to Christov. That a Cattaneo-like law would be relevant in thermohaline convection in star evolution was proposed in 1995 by Herrera and Falcón. They do not develop a detailed model which we do here. We find that with the Cattaneo–Christov law, there is a transition curve (depending on the salt concentration) such that for Cattaneo numbers greater than this transition, the nature of convection changes dramatically.

Analytical validation of COMSOL Multiphysics for theoretical models of Radiofrequency ablation including the Hyperbolic Bioheat transfer equation

In this paper we outline our main findings about the differences between the use of the Bioheat Equation and the Hyperbolic Bioheat Equation in theoretical models for Radiofrequency (RF) ablation. At the moment, we have been working on the analytical approach to solve both equations, but more recently, we have considered numerical models based on the Finite Element Method (FEM). As a first step to use FEM, we conducted a comparative study between the temperature profiles obtained from the analytical solutions and those obtained from FEM. Regarding the differences between both methods, we obtain agreement in less than 5% of relative differences. Then FEM is a good alternative to model heating of biological tissues using BE and HBE in, for example, more complex and realistic geometries.