Special function form solutions of multi-parameter generalized Mittag-Leffler kernel based bio-heat fractional order model subject to thermal memory shocks

The primary objective of this research is to develop a mathematical model, analyze the dynamic occurrence of thermal shock and exploration of how thermal memory with moving line impact of heat transfer within biological tissues. An extended version of the Pennes equation as its foundational framework, a new fractional modelling approach called the Prabhakar fractional operator to investigate and a novel time-fractional interpretation of Fourier's law that incorporates its historical behaviour. This fractional operator has multi parameter generalized Mittag-Leffler kernel. The fractional formulation of heat flow, achieved through a generalized fractional operator with a non-singular type kernel, enables the representation of the finite propagation speed of heat waves. Furthermore, the dynamics of thermal source continually generates a linear thermal shock at predefined locations within the tissue. Introduced the appropriate set of variables to transform the governing equations into dimensionless form. Laplace transform (LT) is operated on the fractional system of equations and results are presented in series form and also expressed the solution in the form of special functions. The article derives analytical solutions for the heat transfer phenomena of both the generalized model, in the Laplace domain, and the ordinary model in the real domain, employing Laplace inverse transformation. The pertinent parameter's influence, such as α, β, γ, a0, b0, to gain insights into the impact of the thermal memory parameter on heat transfer, is brought under consideration to reveal the interesting results with graphical representations of the findings.

Numerical simulation of burn injuries with temperature-dependent thermal conductivity and metabolism under different surface heat sources

In the present paper, the phenomena of heat transport inside human forearm tissue are studied through a one-dimensional nonlinear bioheat transfer model under the influence of various boundary and interface conditions. In this study, we considered temperature-dependent thermal conductivity and metabolic heat to predict temperature distribution inside the forearm tissue. We have studied the temperature distribution inside inner tissue and bone because it has been found that burn injuries are mostly affected by layer thickness. The temperature distribution inside human forearm tissue is analyzed using the finite difference and bvp4c numerical techniques. To examine the accuracy of present numerical code, we compare the obtained numerical result with the exact analytical result in a specific case and find an excellent agreement with the exact results. We also validated our present numerical code with a hybrid scheme based on Runge-Kutta (4,5) and finite difference technique and found it in good compliance. From the obtained results, we observed that the homogeneous heat flux has a greater impact on the temperature at the outer surface of the skin, but the sinusoidal heat flux has a greater impact on the temperature of the subcutaneous layer and inner tissue. It is found that there is no burn injury in the first type of heat source (Tw=44°C), but it may occur in the second and third types of heat sources. It has been observed that by raising the blood perfusion rate and reducing the values of reference metabolic heat, coefficient of thermal conductivity, and heat fluxes, we can manage and reduce burn injuries and achieve hyperthermia temperature.

Analytical and numerical analysis of the dual-pulse lag heat transfer in a three-dimensional tissue subjected to a moving multi-point laser beam

An extensive algorithm based on both analytical and numerical solution methodologies is proposed to obtain transient temperature distributions in a three-dimensional living tissue subjected to a moving single-point and multi-point laser beam by considering metabolic heat generation and blood perfusion rate. Here, the dual-phase lag/Pennes equation is analytically solved by using the method of Fourier series and the Laplace transform. The ability to model single-point or multi-point laser beams as an arbitrary function of place and time is a significant advantage of the proposed analytical approach, which can be used to solve similar heat transfer problems in other living tissues. Besides, the related heat conduction problem is numerically solved based on the finite element method. The effects of laser beam transitional speed, laser power, and the number of laser points on the temperature distribution within the skin tissue are investigated. Moreover, the temperature distribution predicted by the dual-phase lag model is compared with that of the Pennes model under different working conditions. For the studied cases, it is observed that the maximum tissue temperature decreased about 63% by an increase of 6mm/s in the speed of the laser beam. An increase in the laser power from 0.8W/cm³ to 1.2W/cm³ results in a 28 °C increase in the maximum temperature of the skin tissue. It is observed that the maximum temperature predicted by the dual-phase lag model is always lower than that of the Pennes model and the temperature variations over time are sharper, while their results are entirely consistent over the simulation time. The obtained numerical results indicated that the dual-phase lag model is preferred in heating processes occurring at short intervals. Among the investigated parameters, the laser beam speed has the most considerable effect on the difference between the results of the Pennes and the dual-phase lag models.

Effects of non-Fourier bioheat transfer on bone drilling temperature in orthopedic surgery: Theoretical and in vitro experimental investigation

The mechanical drilling process is a typical step in treating bone fractures to fix broken parts with screws and plates. Drilling generates a significant amount of heat and elevates the temperature of the bone, which can cause thermal osteonecrosis and damage to the surrounding bone tissue and nerves. Thermal inertia between heat flux and temperature gradient in nonhomogeneous interior structural medium-like biological tissues is arguable. Therefore, this paper proposes an analytical model of heat propagation in bone drilling for orthopedic surgery based on the hyperbolic Pennes bioheat transfer equation (HPBTE). Drilling experiments in bovine cortical bone samples were also carried out using an infrared thermography approach to confirm the proposed analytical model. Around the drilled hole surface, thermal necrosis is spread out from 1 to 10 mm. Increased feed rate reduces necrosis penetration distance and increases intense bone necrosis. The HPBTE includes thermal relaxation time effect and internal convective function of tissue perfusion rate. As these factors are not considered in the parabolic heat transfer equation (PHTE), the results show that the HPBTE is more accurate in predicting temperature and thermal osteonecrosis than the PHTE. As a result, proposed analytical model is a handy tool for calculating temperature to avoid thermal damage while improving process efficiency. Furthermore, it has the capability of controlling the manual or robotic drilling procedure for minimally invasive operations.

Learning the Hydrophobic, Hydrophilic, and Aromatic Character of Amino Acids from Thermal Relaxation and Interfacial Thermal Conductance

In this study, the thermal relaxation of the 20 naturally occurring amino acids in water and in the protein lysozyme is investigated using transient nonequilibrium molecular dynamics simulations. By modeling the thermal relaxation process, the relaxation times of the amino acids in water occurs over a time scale covering 2-5 ps. For the hydrophobic amino acids, the relaxation time is controlled by the size of the hydrocarbon side chain, while for hydrophilic amino acids, the number of hydrogen bonds does not significantly affect the time scales of the heat dissipation. Our results show that the interfacial thermal conductance at the amino acid-water interface is in the range of 40-80 MW m^-2 K^-1. Hydrophobic and aromatic amino acids tend to have a lower interfacial thermal conductance. Notably, we show that amino acids can be correlated with their thermal relaxation times and molar masses, into simply connected phases with the same hydrophilicity, hydrophobicity, and aromaticity. The thermal relaxation slows down by a factor of up to five in the protein relative to that in water. In the case of the hydrophobic amino acids in the protein lysozyme, the slow down in the thermal relaxation relative to that in water appears to be controlled primarily by the size of the side chain.

Non-Fourier bioheat model for bone grinding with application to skull base neurosurgery

Bone grinding is used to remove the skull bone and access tumors through the nasal passage during cranial base neurosurgery. The generated heat of the spherical diamond tool propagates and could damage the nerves or coagulate the arteries blood. Little is known about the non-Fourier behavior of heat propagation during bone grinding. Therefore, this study develops an analytical model of the hyperbolic Pennes bioheat transfer equation (HPBTE) to calculate the three-dimensional temperature and necrosis in the grinding region. In vitro experimental investigations were carried out, and the contact zone temperature was measured using an infrared thermography system to validate the proposed thermal model. The results demonstrate that the HPBTE provides more reliable temperature evaluation and thermal damage than Fourier or parabolic heat transfer equation (PHTE). Due to the low thermal diffusivity of the bone, the lower grinding feed rate leads to higher temperature amplitude and a smaller radius of the affected zone in the surface and depth of the bone. Also, the intensity of bone necrosis decreases with the increase of the feed rate, and the shape of the damage zone becomes stretched. This analytical model can assess the potential risk of the surgery before clinical trials. Also, it could be used for comparing the different operating conditions to minimize bone necrosis and improve the control process in neurosurgeries.

Analysis of hyperbolic Pennes bioheat equation in perfused homogeneous biological tissue subject to the instantaneous moving heat source

The one-dimensional hyperbolic Pennes bioheat equation under instantaneous moving heat source is solved analytically based on the Eigenvalue method. Comparison with results of in vivo experiments performed earlier by other authors shows the excellent prediction of the presented closed-form solution. We present three examples for calculating the Arrhenius equation to predict the tissue thermal damage analysis with our solution, i.e., characteristics of skin, liver, and kidney are modeled by using their thermophysical properties. Furthermore, the effects of moving velocity and perfusion rate on temperature profiles and thermal tissue damage are investigated. Results illustrate that the perfusion rate plays the cooling role in the heating source moving path. Also, increasing the moving velocity leads to a decrease in absorbed heat and temperature profiles. The closed-form analytical solution could be applied to verify the numerical heating model and optimize surgery planning parameters.

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.

Numerical study on the effects of blood perfusion and body metabolism on the temperature profile of human forearm in hyperthermia conditions

The development of mathematical models for describing the thermal behavior of living tissues under normal or hyperthermia conditions is of increasing importance. In this research, a 3D forearm model based on anthropometric measurement of 25 samples in Tehran, Iran was developed. The tissue temperature distribution is obtained via the Finite Volume Method (FVM) by considering the appropriate boundary conditions, blood perfusion, body metabolism, and the application of hyperthermia conditions on the tissue. The Pennes Bioheat Transfer Equation (PBHTE) is considered in this regard. Also, various thermophysical properties are assumed for the model in order to clarify the effects of such parameters on the tissue temperature distribution. The results of this study indicate that it is possible to provide the desired conditions for many therapeutic processes by controlling the parameters such as blood perfusion, body metabolism and the type of external heat source applied on the tissue. Generally, by decreasing the body metabolism, increasing the blood perfusion rate in tissue and applying a fluctuating heat flux, instead of uniform heat flux on the surface of the forearm skin, it is possible to provide the hyperthermia conditions without causing damages such as burn injuries to the other parts of the tissue. By using the results of this study, the appropriate conditions of hyperthermia can be obtained.