Next Article in Journal
Fractional-Order Model-Free Predictive Control for Voltage Source Inverters
Next Article in Special Issue
Fractional Gradient Optimizers for PyTorch: Enhancing GAN and BERT
Previous Article in Journal
Topological Subordination in Quantum Mechanics
Previous Article in Special Issue
Combined Liouville–Caputo Fractional Differential Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of Fractional Derivatives on Thermo-Mechanical Interaction in Biological Tissues during Hyperthermia Treatment Using Eigenvalues Approach

1
Mathematics Department, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2
Mathematics Department, Faculty of Science, Sohag University, Sohag 82524, Egypt
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(6), 432; https://doi.org/10.3390/fractalfract7060432
Submission received: 3 May 2023 / Revised: 20 May 2023 / Accepted: 24 May 2023 / Published: 26 May 2023
(This article belongs to the Special Issue Advances in Fractional Order Derivatives and Their Applications)

Abstract

:
This article studies the effects of fractional time derivatives on thermo-mechanical interaction in living tissue during hyperthermia treatment by using the eigenvalues approach. A comprehensive understanding of the heat transfer mechanism and the related thermo-mechanical interactions with the patient’s living tissues is crucial for the effective implementation of thermal treatment procedures. The surface of living tissues is traction-free and is exposed to a pulse boundary heat flux that decays exponentially. The Laplace transforms and their associated techniques are applied to the generalized bio-thermo-elastic model, and analytical procedures are then implemented. The eigenvalue approach is utilized to obtain the solution of governing equations. Graphical representations are given for the temperature, the displacement, and the thermal stress results. Afterward, a parametric study was carried out to determine the best method for selecting crucial design parameters that can improve the precision of hyperthermia therapies.

1. Introduction

The temperature behavior of living tissue is still not fully understood due to the difficulty in accurately measuring it in vivo. This is because necropsy can alter tissue temperature characteristics and examining tissue outside of the body lacks perfusion effects. Various in vivo techniques exist to measure thermal behavior, but they yield varying results, and precise measurement of tissue temperature in vivo remains elusive. For effective treatment, it is critical to control the body’s heat transmission pathways. Thermal therapies aim to freeze or heat tumors without damaging surrounding healthy tissue. If doctors could predict how tissue would react thermally, they could plan therapy dosages and durations before surgery. Tissue and blood perfusion’s temperature response is associated with a range of diseases and injuries, including diabetes, skin grafts, and frostbites. The severity of these diseases is influenced by the extent to which blood can reach a certain location. If the thermal characteristics of damaged tissue can be precisely monitored before problems arise, appropriate and efficient therapy can be provided immediately. Heat transport analysis in living tissue is challenging due to its diverse internal structure, which includes perfusion via capillary tubes, convection between blood and tissues, thermal conduction between solid tissues and blood arteries, and heating generation through metabolism, evaporation, and other factors. Pennes’ [1] bioheat model, which is based on Fourier’s law of thermal conduction, is used to represent heat transmission in living tissues. Various thermotherapy procedures, such as laser tissue welding [2], laser operations [3], and hyperthermia [4], are frequently utilized in modern medicine. Skin tissue temperature distributions are influenced by complex factors such as metabolic heating generation and blood circulation, so researchers have expanded on basic relationships by incorporating a wide range of phenomenological mechanisms, such as metabolic heating production, thermal conduction, blood perfusions, radiation, and phase change. Biological tissues undergo diverse stages of change, which can take many different forms. Modified versions of Pennes’ bioheat models, which employ various numerical methods, are available in the relevant literature. Several techniques have been used to study heat transmission in living tissues, including the homotopy perturbation technique [5], Legendre wavelets Galerkin approaches [6,7], and the finite element approach [8]. When a person’s body temperature was unusually low, Esneault and Dillenseger [9] used finite difference approaches to investigate the time-dependent increase in temperature. During thermal therapy, Ghanmi and Abbas [10] performed an analytical investigation of the fractional time derivative in skin tissue. Marin et al. [11] utilized finite element analysis to study the nonlinear bio-heat model in skin tissue caused by external heating sources, while Hobiny and Abbas [12] explored the analytical solution for the fractional bioheat model in tissues with a spherical shape. On the other hand, Keangin and Phadungsak [13] conducted an analysis on the heat transport in porous liver during microwave ablation, specifically focusing on local thermal non-equilibrium. In their research, Keangin et al. [14] investigated the analysis of heat transfer in a deformed liver cancer model treated with a microwave coaxial antenna. Andreozzi et al. [15] studied the effects of a pulsating heat source on interstitial fluid transport in tumor tissues in which the effects of modulating-heat strategies to influence interstitial fluid transport in tissues were analyzed.
Different methods can be used to solve the time-dependent heat transfer equation and model infinite thermal propagations based on classical Fourier heat conduction. Fractional computation has recently proved to be a successful method for modifying many physical models. Fractional derivatives have received significant attention, and various definitions and methods have been developed. The use of fractional time derivatives has enabled the successful modification of many physical models’ processes. Ezzat and colleagues introduced a novel bio-heat model based on the fractional heat conduction formulation, as described in their publications [16,17], while the investigation of transient heating in skin tissues caused by time-dependent thermal therapy, utilizing a heat transport law that incorporates memory, was carried out by Mondal and colleagues in reference [18]. Researchers have conducted several studies on the use of thermal transfer on living tissue [19] to improve treatment methods, develop more accurate temperature prediction technologies, and ultimately, find a cure for cancer. Over the past forty years, a wide range of researchers from various fields, including high-energy particle accelerators, continuum mechanics, acoustics, nuclear engineering, and aeronautics have expressed significant interest in generalized thermo-elastic models from both a technical and mathematical standpoint. The concepts of heat transfer and elasticity are linked in this model. Lord and Shulman [20] developed the thermo-elasticity hypothesis with multiple generalizations. Mondal et al. [21] employed the memory-dependent derivatives on a sliding interval within the framework of the Lord–Shulman model to investigate the heat transfer equation for this problem. Diaz et al. [22] used the finite element method to obtain the solutions of thermo-diffusion types present in living tissue to develop thermal damage. Zhu et al. [22] used the diffusion theory to consider rate process models for the results of thermal damages and the sedimentation of lighting energy in the tissue. When studying the actual phenomena of thermal transfer in finite media, the nonlinear and linear models of heating transfer are extended, and many authors have sought numerical or analytical solutions to the problems posed by these models [23,24,25,26,27,28,29,30,31,32,33]. Despite the growing popularity of laser, microwave, and other forms of thermal therapy in dermatology, thermo-mechanical interactions are seldom considered in current research, despite being central to the thermo-mechanically linked nature of these therapies. Generalized thermo-elastic models govern tissue thermo-elastic behaviors, which include the G-N model, the G-NII model, the DPL, the fractional model, and Li et al. [34,35,36], who further investigated the effects of heat-induced mechanical responses in skin tissues under temperature-dependent properties.
This research aims to create an analytical approach to examine the thermo-mechanical interactions with fractional time derivatives in living tissue that experiences instantaneous heating and has varying thermal and mechanical properties. A generalized thermo-elastic model is constructed that takes into account the tissue structure and variable thermal and mechanical characteristics within the bioheat transfer equation framework. The effect of fractional parameter in temperature, displacement, and thermal stress variations are shown in the graphics.

2. Materials and Methods

The field of bio-thermo-elasticity merges the principles of elasticity and bioheat conduction. The basic equations under fractional time derivatives in the living tissues can be given by [20,37,38]:
λ + μ u j , i j + μ u i , j j γ T , i = ρ 2 u i t 2 ,
k 2 T = 1 + τ o β Γ β + 1 β t β ρ c e T t + ω b ρ b c b T T b + γ T o 2 u t x Q m , 0 < β 1 ,
σ i j = μ u i , j + u j , i + λ u k , k γ T T o δ i j ,
where   t is the time,   τ o is the thermal relaxation time, T b is the blood temperature, c e refers to the specific heat at constant strain, ρ is the tissue mass density, γ = 3 λ + 2 μ α t , α t refers to the linear thermal expansion coefficient, λ ,   μ refer to the Lame’s constants, T is the tissue temperature, k is the tissue thermal conductivity, ω b is the blood perfusion rate, Q m is the metabolic heat generation in skin tissues, e i j are the strain components, σ i j are the stress components, u i are the displacement components, δ i j is the Kronecker symbol and c b is the blood specific heat. The definition of the fractional order derivative is as follows:
β h r , t t β = h r , t h r , 0 , β 0 , I β 1 h r , t t , 0 < β < 1 , h r , t t , β = 1 , ,
I ν h r , t = 0 t t s ν Γ ν h r , s d s ,   ν > 0 ,
lim ν 1 ν h r , t t ν = h r , t t ,
Equation (4) illustrates how the range of local thermal conduction can be characterized by two types of conductivity: standard thermal conduction and heat ballistic conduction. The fractional parameter is β , where 0 < β 1 is used to define these conductivities. For normal conductivity, β = 1 , while for low conductivity 0 < β < 1 . Here, we assume that the surface and the bottom boundary of a limited domain of tissues with thickness L . As a result, the displacement components can be expressed as follows:
u x = u x , t , u y = 0 , u z = 0 .
So that the model has the following form
λ + 2 μ 2 u x 2 γ T x = ρ 2 u t 2 ,
k 2 T x 2 = 1 + τ o β Γ β + 1 β t β ρ c T t + ω b ρ b c b T T b + γ T o 2 u t x Q m ,
σ x x = λ + 2 μ u x γ T T o .
To obtain a solution for the equations, it is necessary to have two sets of initial and boundary conditions that align with the physical model description:
T x , 0 = T b , T x , 0 t = 0 , u x , 0 = 0 , u x , 0 t = 0 ,
σ x x 0 , t = 0 ,   σ x x L , t = 0 , k T x , t x x = 0 = q o t 2 e t t p 16 t p 2 , k T x , t x x = L = 0 ,
where   t p   points to the characteristic time of pulsing heat flux and q o is constant. Now, the dimensionless quantities may be utilized for ease by employing:
T = T T o T o , t , τ o , t p = η c 2 t , τ o , t p ,   x , u = η c x , u ,
σ x x =   σ x x λ + 2 μ ,   q o = q o η c T o K ,   Q m = Q m T o K η 2 c 2 ,
where η = ρ c e K and c 2 = λ + 2 μ ρ .
The primary equations can be expressed in a non-dimensional form by removing the dashes and introducing appropriate parameters, as shown in Equation (13):
2 u x 2 ε 1 T x = 2 u t 2 ,
2 T x 2 = 1 + τ o β Γ β + 1 β t β T t + ε 2 T + ε 3 2 u t x Q m ,
σ x x = u x ε 1 T ,
u x , 0 = 0 , u x , 0 t = 0 , T x , 0 = 0 , T x , 0 t = 0 ,
σ x x 0 , t = 0 , σ x x L , t = 0 , T x , t x x = 0 = q o t 2 e t t p 16 t p 2 , T x , t x x = L = 0 ,
where ε 1 = T o γ e λ e + 2 μ e , ε 2 = ρ b c b ω b K η 2 c 2 , ε 3 = γ e ρ c e .

3. Analytical Solutions in the Transform Domain

Equations (14)–(18) can be transformed using Laplace transforms, as in
f ¯ x , s = L f x , t = 0 f x , t e s t d t .
Therefore, the following equations can be obtained
d 2 u ¯ d x 2 = s 2 u ¯ + ε 1 d T ¯ d x ,
d 2 T ¯ d x 2 = 1 + s β τ o β Γ β + 1 s + ε 2 T ¯ + s ε 3 d u ¯ d x Q m s ,
σ ¯ x x = d u ¯ d x ε 1 T ¯ ,
σ ¯ x x 0 , s = 0 ,   σ ¯ x x L , s = 0 , d T ¯ x , s d x x = 0 = q o t p 8 s t p + 1 3 ,     d T ¯ x , s d x x = L = 0 .
Equations (20) and (21) can be used to represent the vector–matrix differential equation as follows:
d M d x = B M f ,
where M = u ¯ T ¯ d u ¯ d x d T ¯ d x , B = 0 0 1 0 0 0 0 1 s 2 0 0 ε 1 0 1 + s β τ o β Γ β + 1 s + ε 2 1 + s β τ o β Γ β + 1 s ε 3 0 , f = 0 0 0 Q m s .
To solve Equation (24) using the eigenvalue techniques described in [39,40,41,42,43,44,45,46,47,48], Matrix B ’s characteristic equation is expressed as:
α 4 s 2 + 1 + s β τ o β Γ β + 1 s + ε 2 + 1 + s β τ o β Γ β + 1 s ε 3 ε 1 α 2 + 1 + s β τ o β Γ β + 1 s + ε 2 s 2 = 0 ,
Relation (25) has four roots, which are the eigenvalues of matrix B , which are defined by   ± α 1 ,   ± α 2 . The general solutions to the nonhomogeneous system (24) can be obtained by adding the complementary solution of the corresponding homogeneous system to a particular solution of the nonhomogeneous system. To find the particular solution to the nonhomogeneous Equation (24), we need to consider that the inhomogeneous terms in (24) contain functions of the Laplace parameter s. As a result, the particular solution can be expressed as:
M x , s = A 1 X 1 e α 1 x + A 2 X 2 e α 1 x + A 3 X 3 e α 2 x + A 4 X 4 e α 2 x + 0 R 0 0 .
where R = Q m s s + ε 2 1 + s β τ o β Γ β + 1 . Hence, in the Laplace domain, the general solutions of displacement, temperature, and stress can be given by:
u ¯ x , s = A 1 U 1 e α 1 x + A 2 U 2 e α 1 x + A 3 U 3 e α 2 x + A 4 U 4 e α 2 x .
T ¯ x , s = A 1 T 1 e α 1 x + A 2 T 2 e α 1 x + A 3 T 3 e α 2 x + A 4 T 4 e α 2 x + R .
σ ¯ x x = A 1 α 1 U 1 ε 1 T 1 e α 1 x + A 2 α 1 U 2 ε 1 T 2 e α 1 x + A 3 ( α 2 U 3 ε 1 T 3 ) e α 2 x + A 4 α 2 U 4 ε 1 T 4 e α 2 x + R ,
where U i , T i refer to the eigenvectors of displacement and temperature, respectively. The problem boundary conditions can be utilized to determine the values of A 1 ,   A 2 , A 3 , and A 4 . To obtain the final solutions for the displacement, temperature, and stress distributions, the Stehfest approach [49] can be used, which has an inverse function f x , t defined by a specific formulation:
f x , t = l n 2 t j = 1 M V j f ¯ x , j l n 2 t ,
where V j can be given by
V j = 1 M 2 + 1 k = i + 1 2 m i n i , M 2 k M 2 + 1 2 k ! M 2 k ! k ! i k ! 2 k 1 !

4. Numerical Outcomes and Discussion

To demonstrate the theoretical results discussed in the above sections, the numerical values of the physical parameters are presented. The material parameters for living tissues at the reference temperature, used in the following calculation, are denoted as follows [35]:
λ = 8.27 × 10 8 N m 2 , α t = 1 × 10 4 k 1 , T o = 310 k ,   t p = 0.15 ,
μ = 3.446 × 10 7 N   m 2 , c b = 3770   J k g 1 k 1 ,   ρ b = 1060 k g m 3 ,
τ o = 0.05 ,   K = 0.235   W m 1 k 1 ,   Q m = 1.19 × 10 3 W m 3 ,   ρ = 1190 k g m 3 , c e = 3600   J k g 1 k 1
The numerical values of the calculated physical quantities under the fractional biothermo-elastic model, considering one thermal relaxation time and using the previous parameters, are displayed in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12. Numerical calculations have been performed at time t = 0.2 to determine the temperature variations, the variation of displacement, and stress variation along the distance x under different values of the studied parameters as in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12. Figure 1, Figure 4, Figure 7 and Figure 10 display the temperature variation along the distance x . It can be observed from the figures that the temperature initially peaks at the tissue surface x = 0 due to the exponentially decaying pulse boundary heat flux. As the distance x increases, the temperature steadily decreases. Figure 2, Figure 5, Figure 8 and Figure 11 display the variations of displacement along the distance x . Observing the figures, it can be noted that the displacement reaches its highest negative values at the tissue surface x = 0 . Subsequently, it progressively increases towards peak values near the surface before decreasing back to zero. Figure 3, Figure 6, Figure 9 and Figure 12 show the variation of stress σ x x along the distance x . It can be observed that the stress σ x x starts from zero and ends at zero to comply with the boundary condition of the problem. Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12 can be classified into four distinct groups.
In the first group, Figure 1, Figure 2 and Figure 3 show the variations of temperature, displacement, and stress under various models. In Figure 1, Figure 2 and Figure 3, the solid line (—) refers to the Pennes model (Pennes model) without thermal relaxation time ( τ o = 0 ) and without fractional time derivative ( β = 1 ), the dashed line (---) points to the single-phase lag model (SPL model) with one thermal relaxation time ( τ o = 0.02 ) and without fractional time derivative ( β = 1 ), while the dotted line (…) refers to the single-phase lag model under fractional time derivative (FSPL model) with one thermal relaxation time ( τ o = 0.02 ) and with fractional time derivative ( β = 0.5 ). The significant effects on the quantities under consideration are evident due to the various models.
In the second group, Figure 4, Figure 5 and Figure 6 display the variation of temperature, displacement, and stress under different values of the fractional parameter ( β = 1 , 0.5 , 0.1 ) when ( τ o = 0.02 ). It is evident that the fractional time derivatives are responsible for the evident significant impacts on the quantities under consideration. A decrease in fractional time derivatives weakens the effect of thermo-mechanical propagation, as evidenced by a decrease in the maximum amplitude of temperature, displacement, and stress.
In the third group, Figure 7, Figure 8 and Figure 9 show the variation of temperature, displacement, and stress under various values of thermal relaxation time ( τ o = 0 ,   0.02 ,   0.2 ) under fractional time derivative when ( β = 0.5 ). It is evident that the thermal relaxation time has significant effects on the quantities under consideration. The maximum amplitude of the temperature, displacement, and stress decreases with the increase the thermal relaxation time, which means that the thermal relaxation time is apt to weaken the effect of thermo-mechanical propagation.
In the fourth group, Figure 10, Figure 11 and Figure 12 show the variations of temperature, displacement, and stress under various of the characteristic time of pulsing heat flux ( t p = 0.1 ,   0.15 ,   0.2 ) with one thermal relaxation time ( τ o = 0.02 ) under fractional time derivative when ( β = 0.5 ). It is observed that the characteristic time of pulsing heat flux has a significant effect on the quantities under consideration. An increase in the characteristic time of pulsing heat flux weakens the effect of thermo-mechanical propagation, as evidenced by a decrease in the maximum amplitude of temperature, displacement, and stress.

5. Conclusions

The aim of this research paper was to examine how biological tissues react to a sudden pulsing heat flux load by employing generalized thermo-elasticity under the fractional time derivatives framework, which incorporates one thermal relaxation time.
  • This study specifically focused on the impacts of fractional parameter, thermal relaxation time, and the pulsing heat flux characteristic time on bio-thermo-elastic behaviors.
  • A comparative analysis was conducted between the fractional single-phase lag model (FSPL model) and previous single-phase lag (SPL model) and Pennes (Pennes model) models.
  • The findings of this study, which presented a modified thermo-elasticity approach, offered a fresh perspective on the propagation of thermal waves, representing the first attempts in this area.
  • These results significantly contribute to enhancing our understanding of thermo-elastic behavior in living tissue.

Author Contributions

Conceptualization, A.H. and I.A.; methodology: A.H. and I.A.; validation: A.H. and I.A.; formal analysis: A.H. and I.A.; investigation: A.H. and I.A.; resources: I.A.; data curation: A.H. and I.A.; writing—original draft preparation: A.H. and I.A.; writing—review and editing: A.H.: visualization: I.A.; supervision: A.H. and I.A.; project administration: I.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research work was funded by Institutional Fund Projects under grant no. (IFPIP: 93-130-1443). The authors gratefully acknowledge the technical and financial support provided by the Ministry of Education and King Abdulaziz University, DSR, Jeddah, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pennes, H.H. Analysis of tissue and arterial blood temperatures in the resting human forearm. J. Appl. Physiol. 1948, 1, 93–122. [Google Scholar] [CrossRef] [PubMed]
  2. Gabay, I.; Abergel, A.; Vasilyev, T.; Rabi, Y.; Fliss, D.M.; Katzir, A. Temperature-controlled two-wavelength laser soldering of tissues. Lasers Surg. Med. 2011, 43, 907–913. [Google Scholar] [CrossRef] [PubMed]
  3. Zhou, J.; Chen, J.; Zhang, Y. Dual-phase lag effects on thermal damage to biological tissues caused by laser irradiation. Comput. Biol. Med. 2009, 39, 286–293. [Google Scholar] [CrossRef] [PubMed]
  4. Mahjoob, S.; Vafai, K. Analytical characterization of heat transport through biological media incorporating hyperthermia treatment. Int. J. Heat Mass Transf. 2009, 52, 1608–1618. [Google Scholar] [CrossRef]
  5. Gupta, P.K.; Singh, J.; Rai, K. Numerical simulation for heat transfer in tissues during thermal therapy. J. Therm. Biol. 2010, 35, 295–301. [Google Scholar] [CrossRef]
  6. Kumar, P.; Kumar, D.; Rai, K. A numerical study on dual-phase-lag model of bio-heat transfer during hyperthermia treatment. J. Therm. Biol. 2015, 49, 98–105. [Google Scholar] [CrossRef]
  7. Yadav, S.; Kumar, D.; Rai, K.N. Finite Element Legendre Wavelet Galerkin Approch to Inward Solidification in Simple Body under Most Generalized Boundary Condition. Z. Nat. A 2014, 69, 501–510. [Google Scholar] [CrossRef]
  8. Gupta, P.K.; Singh, J.; Rai, K.; Rai, S. Solution of the heat transfer problem in tissues during hyperthermia by finite difference–decomposition method. Appl. Math. Comput. 2013, 219, 6882–6892. [Google Scholar] [CrossRef]
  9. Dillenseger, J.-L.; Esneault, S. Fast FFT-based bioheat transfer equation computation. Comput. Biol. Med. 2010, 40, 119–123. [Google Scholar] [CrossRef]
  10. Ghanmi, A.; Abbas, I.A. An analytical study on the fractional transient heating within the skin tissue during the thermal therapy. J. Therm. Biol. 2019, 82, 229–233. [Google Scholar] [CrossRef]
  11. Marin, M.; Hobiny, A.; Abbas, I. Finite element analysis of nonlinear bioheat model in skin tissue due to external thermal sources. Mathematics 2021, 9, 1459. [Google Scholar] [CrossRef]
  12. Hobiny, A.; Abbas, I. Analytical solutions of fractional bioheat model in a spherical tissue. Mech. Based Des. Struct. Mach. 2021, 49, 430–439. [Google Scholar] [CrossRef]
  13. Keangin, P.; Rattanadecho, P. Analysis of heat transport on local thermal non-equilibrium in porous liver during microwave ablation. Int. J. Heat Mass Transf. 2013, 67, 46–60. [Google Scholar] [CrossRef]
  14. Keangin, P.; Wessapan, T.; Rattanadecho, P. Analysis of heat transfer in deformed liver cancer modeling treated using a microwave coaxial antenna. Appl. Therm. Eng. 2011, 31, 3243–3254. [Google Scholar] [CrossRef]
  15. Andreozzi, A.; Iasiello, M.; Netti, P. Effects of pulsating heat source on interstitial fluid transport in tumour tissues. J. R. Soc. Interface 2020, 17, 20200612. [Google Scholar] [CrossRef]
  16. Ezzat, M.A.; AlSowayan, N.S.; Al-Muhiameed, Z.I.; Ezzat, S.M. Fractional modelling of Pennes’ bioheat transfer equation. Heat Mass Transf. 2014, 50, 907–914. [Google Scholar] [CrossRef]
  17. Ezzat, M.A.; El-bary, A.A.; Al-sowayan, N.S. Tissue responses to fractional transient heating with sinusoidal heat flux condition on skin surface. Anim. Sci. J. 2016, 87, 1304–1311. [Google Scholar] [CrossRef]
  18. Mondal, S.; Sur, A.; Kanoria, M. Transient heating within skin tissue due to time-dependent thermal therapy in the context of memory dependent heat transport law. Mech. Based Des. Struct. Mach. 2019, 49, 271–285. [Google Scholar] [CrossRef]
  19. Andreozzi, A.; Brunese, L.; Iasiello, M.; Tucci, C.; Vanoli, G.P. Modeling Heat Transfer in Tumors: A Review of Thermal Therapies. Ann. Biomed. Eng. 2019, 47, 676–693. [Google Scholar] [CrossRef]
  20. Lord, H.W.; Shulman, Y. A generalized dynamical theory of thermoelasticity. J. Mech. Phys. Solids 1967, 15, 299–309. [Google Scholar] [CrossRef]
  21. Sur, A.; Mondal, S.; Kanoria, M. Influence of Moving Heat Source on Skin Tissue in the Context of Two-Temperature Caputo–Fabrizio Heat Transport Law. J. Multiscale Model. 2019, 11, 2050002. [Google Scholar] [CrossRef]
  22. Díaz, S.H.; Nelson, J.S.; Wong, B.J. Rate process analysis of thermal damage in cartilage. Phys. Med. Biol. 2002, 48, 19. [Google Scholar] [CrossRef] [PubMed]
  23. Ghazanfarian, J.; Saghatchi, R.; Patil, D. Implementation of Smoothed-Particle Hydrodynamics for non-linear Pennes’ bioheat transfer equation. Appl. Math. Comput. 2015, 259, 21–31. [Google Scholar] [CrossRef]
  24. Li, L.; Liang, M.; Yu, B.; Yang, S. Analysis of thermal conductivity in living biological tissue with vascular network and convection. Int. J. Therm. Sci. 2014, 86, 219–226. [Google Scholar] [CrossRef]
  25. Selvi, C.; Srinivas, A.; Sreenadh, S. Peristaltic transport of a power-law fluid in an elastic tube. J. Taibah Univ. Sci. 2018, 12, 687–698. [Google Scholar] [CrossRef]
  26. Khan, A.A.; Bukhari, S.R.; Marin, M.; Ellahi, R. Effects of chemical reaction on third-grade MHD fluid flow under the influence of heat and mass transfer with variable reactive index. Heat Transf. Res. 2019, 50, 1061–1080. [Google Scholar] [CrossRef]
  27. Lata, P.; Himanshi. Orthotropic magneto-thermoelastic solid with higher order dual-phase-lag model in frequency domain. Struct. Eng. Mech. 2021, 77, 315–327. [Google Scholar] [CrossRef]
  28. Abbas, I.A. A GN model for thermoelastic interaction in an unbounded fiber-reinforced anisotropic medium with a circular hole. Appl. Math. Lett. 2013, 26, 232–239. [Google Scholar] [CrossRef]
  29. Abbas, I.A. Analytical solution for a free vibration of a thermoelastic hollow sphere. Mech. Based Des. Struct. Mach. 2015, 43, 265–276. [Google Scholar] [CrossRef]
  30. Hobiny, A.; Abbas, I. A GN model on photothermal interactions in a two-dimensions semiconductor half space. Results Phys. 2019, 15, 102588. [Google Scholar] [CrossRef]
  31. Marin, M.; Ellahi, R.; Vlase, S.; Bhatti, M. On the decay of exponential type for the solutions in a dipolar elastic body. J. Taibah Univ. Sci. 2020, 14, 534–540. [Google Scholar] [CrossRef]
  32. Abo-Dahab, S.M.; Abouelregal, A.E.; Marin, M. Generalized thermoelastic functionally graded on a thin slim strip non-Gaussian laser beam. Symmetry 2020, 12, 1094. [Google Scholar] [CrossRef]
  33. Alzahrani, F.; Hobiny, A.; Abbas, I.; Marin, M. An Eigenvalues Approach for a Two-Dimensional Porous Medium Based upon Weak, Normal and Strong Thermal Conductivities. Symmetry 2020, 12, 848. [Google Scholar] [CrossRef]
  34. Li, X.; Li, C.; Xue, Z.; Tian, X. Analytical study of transient thermo-mechanical responses of dual-layer skin tissue with variable thermal material properties. Int. J. Therm. Sci. 2018, 124, 459–466. [Google Scholar] [CrossRef]
  35. Li, X.; Xue, Z.; Tian, X. A modified fractional order generalized bio-thermoelastic theory with temperature-dependent thermal material properties. Int. J. Therm. Sci. 2018, 132, 249–256. [Google Scholar] [CrossRef]
  36. Li, X.; Li, C.; Xue, Z.; Tian, X. Investigation of transient thermo-mechanical responses on the triple-layered skin tissue with temperature dependent blood perfusion rate. Int. J. Therm. Sci. 2019, 139, 339–349. [Google Scholar] [CrossRef]
  37. Xu, F.; Seffen, K.; Lu, T. Non-Fourier analysis of skin biothermomechanics. Int. J. Heat Mass Transf. 2008, 51, 2237–2259. [Google Scholar] [CrossRef]
  38. Ahmadikia, H.; Fazlali, R.; Moradi, A. Analytical solution of the parabolic and hyperbolic heat transfer equations with constant and transient heat flux conditions on skin tissue. Int. Commun. Heat Mass Transf. 2012, 39, 121–130. [Google Scholar] [CrossRef]
  39. Abbas, I.A. Eigenvalue approach on fractional order theory of thermoelastic diffusion problem for an infinite elastic medium with a spherical cavity. Appl. Math. Model. 2015, 39, 6196–6206. [Google Scholar] [CrossRef]
  40. Othman, M.I.A.; Abbas, I.A. Eigenvalue approach for generalized thermoelastic porous medium under the effect of thermal loading due to a laser pulse in DPL model. Indian J. Phys. 2019, 93, 1567–1578. [Google Scholar] [CrossRef]
  41. Kumar, R.; Miglani, A.; Rani, R. Eigenvalue formulation to micropolar porous thermoelastic circular plate using dual phase lag model. Multidiscip. Model. Mater. Struct. 2017, 13, 347–362. [Google Scholar] [CrossRef]
  42. Kumar, R.; Miglani, A.; Rani, R. Analysis of micropolar porous thermoelastic circular plate by eigenvalue approach. Arch. Mech. 2016, 68, 423–439. [Google Scholar]
  43. Gupta, N.D.; Das, N.C. Eigenvalue approach to fractional order generalized thermoelasticity with line heat source in an infinite medium. J. Therm. Stress. 2016, 39, 977–990. [Google Scholar] [CrossRef]
  44. Santra, S.; Lahiri, A.; Das, N.C. Eigenvalue Approach on Thermoelastic Interactions in an Infinite Elastic Solid with Voids. J. Therm. Stress. 2014, 37, 440–454. [Google Scholar] [CrossRef]
  45. Baksi, A.; Roy, B.K.; Bera, R.K. Eigenvalue approach to study the effect of rotation and relaxation time in generalized magneto-thermo-viscoelastic medium in one dimension. Math. Comput. Model. 2006, 44, 1069–1079. [Google Scholar] [CrossRef]
  46. Das, N.C.; Lahiri, A.; Giri, R.R. Eigenvalue approach to generalized thermoelasticity. Indian J. Pure Appl. Math. 1997, 28, 1573–1594. [Google Scholar]
  47. Abbas, I.A.; Abdalla, A.-E.-N.N.; Alzahrani, F.S.; Spagnuolo, M. Wave propagation in a generalized thermoelastic plate using eigenvalue approach. J. Therm. Stress. 2016, 39, 1367–1377. [Google Scholar] [CrossRef]
  48. Abbas, I.; Hobiny, A.; Marin, M. Photo-thermal interactions in a semi-conductor material with cylindrical cavities and variable thermal conductivity. J. Taibah Univ. Sci. 2020, 14, 1369–1376. [Google Scholar] [CrossRef]
  49. Stehfest, H. Algorithm 368: Numerical inversion of Laplace transforms [D5]. Commun. ACM 1970, 13, 47–49. [Google Scholar] [CrossRef]
Figure 1. The temperature variations T via x under three different models.
Figure 1. The temperature variations T via x under three different models.
Fractalfract 07 00432 g001
Figure 2. The variations of displacement u via x under three different models.
Figure 2. The variations of displacement u via x under three different models.
Fractalfract 07 00432 g002
Figure 3. The variations of stress σ x x via x under three different models.
Figure 3. The variations of stress σ x x via x under three different models.
Fractalfract 07 00432 g003
Figure 4. The temperature variations T via x under varying fractional parameter β .
Figure 4. The temperature variations T via x under varying fractional parameter β .
Fractalfract 07 00432 g004
Figure 5. The variations of displacement u via x under varying fractional parameter β .
Figure 5. The variations of displacement u via x under varying fractional parameter β .
Fractalfract 07 00432 g005
Figure 6. The variations of stress σ x x via x under varying fractional parameter β .
Figure 6. The variations of stress σ x x via x under varying fractional parameter β .
Fractalfract 07 00432 g006
Figure 7. The temperature variations T via x under varying thermal relaxation time τ o .
Figure 7. The temperature variations T via x under varying thermal relaxation time τ o .
Fractalfract 07 00432 g007
Figure 8. The variations of displacement u via x under varying thermal relaxation time τ o .
Figure 8. The variations of displacement u via x under varying thermal relaxation time τ o .
Fractalfract 07 00432 g008
Figure 9. The variations of stress σ x x via x under varying thermal relaxation time τ o .
Figure 9. The variations of stress σ x x via x under varying thermal relaxation time τ o .
Fractalfract 07 00432 g009
Figure 10. The temperature variations T via x under varying pulsing heat flux characteristic time t p .
Figure 10. The temperature variations T via x under varying pulsing heat flux characteristic time t p .
Fractalfract 07 00432 g010
Figure 11. The variations of displacement u via x under varying pulsing heat flux characteristic time t p .
Figure 11. The variations of displacement u via x under varying pulsing heat flux characteristic time t p .
Fractalfract 07 00432 g011
Figure 12. The variations of stress σ x x via x under varying pulsing heat flux characteristic time t p .
Figure 12. The variations of stress σ x x via x under varying pulsing heat flux characteristic time t p .
Fractalfract 07 00432 g012
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hobiny, A.; Abbas, I. The Effect of Fractional Derivatives on Thermo-Mechanical Interaction in Biological Tissues during Hyperthermia Treatment Using Eigenvalues Approach. Fractal Fract. 2023, 7, 432. https://doi.org/10.3390/fractalfract7060432

AMA Style

Hobiny A, Abbas I. The Effect of Fractional Derivatives on Thermo-Mechanical Interaction in Biological Tissues during Hyperthermia Treatment Using Eigenvalues Approach. Fractal and Fractional. 2023; 7(6):432. https://doi.org/10.3390/fractalfract7060432

Chicago/Turabian Style

Hobiny, Aatef, and Ibrahim Abbas. 2023. "The Effect of Fractional Derivatives on Thermo-Mechanical Interaction in Biological Tissues during Hyperthermia Treatment Using Eigenvalues Approach" Fractal and Fractional 7, no. 6: 432. https://doi.org/10.3390/fractalfract7060432

APA Style

Hobiny, A., & Abbas, I. (2023). The Effect of Fractional Derivatives on Thermo-Mechanical Interaction in Biological Tissues during Hyperthermia Treatment Using Eigenvalues Approach. Fractal and Fractional, 7(6), 432. https://doi.org/10.3390/fractalfract7060432

Article Metrics

Back to TopTop