Heat transfer and the influence of electromagnetic field and thermal time during temperature distribution in human tissue
Khaled Lotfy1 , Alla El-Bary2
1Department of Mathematics, Faculty of Science, Zagazig University, P.O. Box44519, Zagazig, Egypt
1Department of Mathematics, College of Science, Taibah University, Madinah, Saudi Arabia
2Arab Academy for Science, Technology and Maritime Transport, P.O. Box 1029, Alexandria, Egypt
2National Committee for Mathematics, Academy of Scientific Research and Technology, Cairo, Egypt
Journal of Engineering and Thermal Sciences, (in Press).
Received 20 May 2021; received in revised form 3 August 2021; accepted 17 August 2021; published 12 October 2021
The Penne’s bio-heat transfer equation can be used with normal mode technique to describe the characterizing of the temperature fluctuation in muscles. The analytical solutions of heating pattern is obtained in a closed-form when the propagation of ultrasonic waves in tissue system are taken into consideration. The impact of external electromagnetic field is used to investigate the influence of temporal and spatial distributions of temperature. The numerical simulations during the 2D and 3D graphs can be obtained for human tissue in simplified geometry in the context of the derived method. The normal mode analysis is used as a mathematical technique to solve the bioheat transfer equation analytically with some conditions to get the complete solution of the main variables in this model.
- The Penne's bio-heat transfer equation and normal mode technique are used to describe the characterizing of the temperature fluctuation in muscles.
- The impact of external electromagnetic field is used to investigate the influence of temporal and spatial distributions of temperature
- The numerical simulations during the 2D and 3D graphs can be obtained for human tissue in simplified geometry in the context of the derived method.
Keywords: human tissue, Bioheat transfer equation, normal mode, analytical solution.
Increased in the present uses of technological devices and equipment such as computers, electrical appliances and home and led electric and magnetic fields emanating from these devices to the increasing environmental pollution mail which plays a major role in a rapid decrease in the intensity of the magnetic field of the Earth and human creatures adapted itself with the continuing decline in magnetic energy, but in turn lost a similar amount of capacity of vital functions within objects. The researchers proved that the decrease in the intensity of the Earth’s magnetic field associated damages arising from the impact of the electronic environment in which it operates to break down the cellular composition of the cells within the body, symptoms of feeling pain and roughness, arthritis, headaches and fatigue.
The magnetism is one of the fundamental forces, all of human civilization was born and lived her life under the magnetism arising from the Earth’s magnetic field, it is known that the space is filled with cosmic rays in the form of particles of nuclear high-power consists of the nuclei of atoms of elements light and heavy electrons moving at high speeds emitted from the sun, stars and galaxies across the universe. He explains that the cosmic rays are radiation serious are booked together in layers of the atmosphere under the affected spin in the belts , “Van Allen” radiation , which revolve in which ions high energy coming from space just between 4 thousand to 16 thousand kilometers from the Earth's surface and occur reservations for these rays thanks to the influence of the magnetic field of the planet, according to information monitored by the U.S. satellite “Explorer 1” in 1958, when he came with certain information about the this pouring barrage of deadly radiation surging in space cards awesome! If they hit us what carried us on this planet life. Were it not for the skies the atmosphere and the Earth’s magnetic perished all organisms on the planet.
On the other side, the magnetic positive impacts in our daily lives, and as the use of magnetic forces us due to the ancient civilizations.
Heat transfer in biological systems is relevant in many diagnostic and therapeutic applications that involve changes in temperature. For example, in hyperthermia the tissue temperature is elevated to 42–43°C using ultrasound by Seip and Ebbini . Stoll  studied the thermal properties of skin to understand conditions leading to thermal damage (burns) to skin. Stoll  investigated the contacting case of skin with a hot objects. Diederich et al , studied the effect of electromagnetic fields and ultrasound waves for skin heating and deeper tissue.
Carslaw and Jaeger  and Lienhard  predicted the heat transport by analytical technique with the numerical methodology. Riu  studies the rise of temperature in the context of the bioheat transfer method which can be solved analytically and numerically by finite element technique for simple geometries Bowman and Martin , . Davies  investigated a new models analytically during the temperature-dependent increases but this model is very difficult in perfusion, in this case the dependent of linear temperature can be used, on the other hand Martin  used a numerical simulations to describe the analytical model. Akrin  used the equation of bioheat transfer which it can be obtained in a wide range with many applications in bio-engineering to describe the blood heat transport during a perfuse tissues. Erdmann  investigated the finite element method to optimize the nonlinear form of the bioheat equation for optimizing regional hyper-thermia. The two-dimensional (2D) bio-thermal technique with the ultrasound applicators based on the bioheat equation can be solved by a finite difference equation, Yreus . The finite difference methods at the boundary element have been used to solve the bioheat equation -.
Othman and Lotfy  studied transient disturbance in a half-space under generalized magneto-thermoelasticity with moving internal heat source. Othman and Lotfy  studied the plane waves in generalized thermo-microstretch elastic half-space by using a general model of the equations of generalized thermo-microstretch for a homogeneous isotropic elastic half space. Othman and Lotfy  studied the generalized thermo-microstretch elastic medium with temperature dependent properties for different theories. Othman and Lotfy  studied the effect of magnetic field and inclined load in micropolar thermoelastic medium possessing cubic symmetry under three theories. The normal mode analysis was used to obtain the exact expression for the temperature distribution, thermal stresses, and the displacement components.
In this paper, the bioheat transfer equation can be solved by the normal mode technique analytically to obtain the complete solution for the basic quantities in biological muscles. In this problem, the soft tissue is used as a viscoelastic medium with different relaxation times in the bioheat transfer equation.
2. Bioheat transfer equation
During the biological muscles (see Fig. 1), the temperature evaluation can be obtained in the context of the Penne’s bioheat equation, which is:
where, , , and represent the temperature distribution, the density of human tissue, the heat capacity of human tissue, the blood flow diffusion respectively. refers to the blood heat capacity, represents the blood flow perfusion, expresses the density of human blood. On the other hand, , and express the arterial blood temperature, the absorbed power density and relaxation time respectively.
3. Magnetic field equations
We consider rectangular coordinate system having origin on the surface 0 and -axis pointing vertically into the muscles. A magnetic field with constant intensity acts parallel to the bounding plane (table as the direction of the -axis (Fig. 2)). Due to the application of initial magnetic field , there are results of an induced magnetic field and an induced electric field . The simplified linear equations of electrodynamics of slowly moving medium for a homogeneous,
Thermally and electrically conducting viscoelastic medium is:
where is the partied velocity of the muscles, and the small effect of temperature gradient on is ignored. The dynamic displacement vector is actually measured from a steady state deformed position and the deformation is supposed to be small.
Fig. 1. Geometry of the problem
Fig. 2. The direction of magnetic field
The components of the magnetic intensity vector in the viscoelastic (blood in muscles) medium are:
The electric intensity vector is normal to both the magnetic intensity and the displacement vectors. Thus, it has the components:
The current density vector be parallel to , thus:
If we restrict our analysis to plane parallel to -plane with displacement vector (distance between cells in tissue) . Body couples and heat sources can be written in visco-elasticity tissue by following the equations given by Minagawa et al. , Green Lindsay  and Othman and Baljeet  as:
Introducing potential functions defined by:
The field Eqs. (10)-(11), we can reduce to:
where , is coefficient of linear expansion, , and are representing effect of viscosity constants, is the time, is the initial uniform magnetic intensity vector, is the induced magnetic field vector, is the induced electric field vector, = is magnetic permeability, is the electric permeability, is a constant modulus of visco-elasticity, and .
The following dimensionless parameters were defined as:
where , , and is the tissue length.
Eqs. (13) and (14) take the following form (dropping the dashed for convenience):
The dimensional form of Eq. (1) can be obtained as:
4. The normal mode analysis
The normal mode technique is used to solve Eq. (1) as:
where, is the wave number in the -direction, is a complex time constant, , , and are the amplitude of the functions , , and respectively.
Substituting from Eqs. (21) into Eqs. (17), (18) and (20) we obtain:
where , , and .
Eliminating between Eqs. (22) and (24), we get the following fourth order ordinary differential equation satisfied by :
Eq. (25) can be factorized as:
where , are the roots of the following characteristic equation:
The solution of Eq. (28) is given by:
In a similar manner, we get:
where are parameters depending on and .
The solution of Eq. (23) has the form:
Using Eqs. (30) and (32), in order to obtain the amplitude of the displacement components and , which are bounded as , then Eqs. (33) and (34) become:
5. Boundary conditions
The plane boundary subjects to an instantaneous normal point force and the boundary surface is isothermal, the boundary conditions at the vertical plan and in the beginning of the operation of acting magnetic field at are (Othman and Lotfy 2011):
Using Eqs. (13), (16), on the non-dimensional boundary conditions and using Eqs. (37), (38), (39), we obtain the expressions of displacements and temperature distribution for the body under the effect radiation of magnetic field as follows:
where , .
Invoking the boundary conditions Eqs. (37-39) at the surface of the plate, we obtain a system of three equations. After applying the inverse of matrix method, we have the values of the four constants , . Hence, we obtain the expressions of displacements, and temperature distribution for the muscles:
where and .
6. Computational results
The properties of typical tissue and blood are used in these calculations as follow: 0.01, 1000 Kg/m3, 4200 J/Kg° C, 0.01, 0.5 W/mC, 0.5 Kg/m3 s, 0.2×10-15 cm3, 9.4×1011 dyne/cm2, 4.0×1011 dyne/cm2, , 0.1.
On the other hand, the parameter values of this problem are taken in dimensionless which can be given by = 1.0, 0.01, and the dimension of tissue is 3 cm. Fig. 3 displays the distribution of temperature at different three values of times against the distance. In this figure notice that the temperature distributions in tissue at the early stage of heating decreases with increasing the depth of tissue. The soft tissue can be considered as a viscous liquid for ultrasound propagation. The amount of heat which it is generated in the human tissue is more dependent on the tissue absorption coefficient. Acoustic absorption is due to the shear relaxation invariably cannot be accounted for by one relaxation time instead a distribution in continuous form of thermal memories must be used. In this case, this cause the relaxation process arises from the finite time taken for molecules to diffuse between adjacent shearing layers in the medium. Fig. 4 show the horizontal displacement distribution (distance between cells in the tissue) at different times as function of propagation distance. From this figure it can seen that the displacement distribution of tissue increases sharp in the start when the magnetic field acting on the tissue and hence smooth decreasing and converge to zero value (stable state) with the increase of the distance, after some time the displacement distribution more equilibrium this illustrated from the figure.
Fig. 3. Temperature distribution as function of distance for different times ( 0.5 W/m°C, 0.5 Kg/m3s)
Fig. 4. Horizontal displacement distribution u as function of distance for different times ( 0.5 W/m°C, 0.5 Kg/m3s)
Fig. 5. Vertical displacement distribution as function of distance for different times ( 0.5 W/m°C, 0.5 Kg/m3s)
Fig. 6. Temperature distribution as function of distance for different values of relaxation time ( 0.5 W/m°C, 0.5 Kg/m3s)
Fig. 5 show the vertical displacement distribution at various times which it taken as a function of distance propagation. From this figure the displacement distribution of tissue increases to arrive the maximum value (due to the magnetic field) in the start and hence smooth decreasing and converge to zero value (stable state) with the increase of the distance, after some time the displacement distribution more equilibrium. Fig. 6 displays the distribution of temperature at various dimensionless relaxation time as function of propagation distance. It can be seen that the temperature of tissue decreases with the increase of the thermal memories.
Fig. 7. Horizontal displacement distribution as function of distance for different values of relaxation time ( 0.5 W/m°C, 0.5 Kg/m3s)
Fig. 8. Vertical displacement distribution as function of distance for different values of relaxation time ( 0.5 W/m°C, 0.5 Kg/m3s)
Fig. 9. Effect of blood perfusion levels on temperature response as function of distance ( 0.5 W/m°C, 0.1)
Fig. 10. Effect of blood perfusion levels on Horizontal displacement response as function of distance ( 0.5 W/m°C, 0.1)
Fig. 7 shows horizontal displacement distribution at various dimensionless thermal memories which taken as a function of distance propagation. It can seen that the distance between the cell in tissue sharp increases in the first rang and then smooth decreases with the increase of the distance and converge to zero value when distance increases. Also, the amplitudes horizontal displacement decreases with the increase of the relaxation time. Fig. 8 shows vertical displacement distribution v at different non-dimensional relaxation time as function of propagation distance. From this figure it can seen that the distance between the cell in tissue smooth increases in the first rang and then smooth decreases also with the increase of the distance and converge to zero value when distance increases. Also, the amplitudes horizontal displacement decreases with the increase of the thermal memories. Fig. 9 displays the temperature of tissue responses corresponding to three various blood perfusion levels. For the case of 0.5 Kg/m3s, the tissue temperature appears much smaller than that of using 0.3 Kg/m3s at the same position. It can be seen that large blood perfusion tends to prevent the biological body from burn injury. Fig. 10 shows the horizontal displacement distribution at different three different blood perfusion levels. It clear from this figure the displacement between the cells in tissue increasing with blood perfusion increases. Fig. 11 shows the vertical displacement distribution v at different three different blood perfusion levels. It clears from this figure the displacement between the cells in this direction in tissue decreases with blood perfusion increases (it called spasm).
Fig. 11. Effect of blood perfusion levels on vertical displacement response function of distance ( 0.5W/m°C, 0.1)
Fig. 12 shows the distribution of temperature in the human tissue, and depicted the surface tissue temperature and displacement increases immediately after the exposure, while in the deeper tissues, the temperature decreases slightly until after propagation in the - plane. Fig. 13 depicts the tissue blood perfusion levels corresponding to the temperature responses and distance. The tissue temperature increases with the increase of the thermal conductivity of tissue It can be seen that large blood perfusion tends to prevent the biological body from burn injury.
Fig. 12. Temperature distribution in 3-D (0.5 W/m°C, 0.5 Kg/m3s, 0.1)
Fig. 13. Effect of temperature response on blood perfusion levels as function of distance ( 0.5 W/m°C, 0.1)
Effects of the individual’s physiological variables (such as blood perfusion, thermal conductivity and relaxation time) were investigated in detail. The present solution in this paper is very useful for a variety of bio-thermal studies. The curves in the context of the theories decrease or increase exponentially with increasing , this indicate that the thermoelastic waves are unattenuated and nondispersive, where purely thermoelastic waves undergo both attenuation and dispersion. The presence of magnetic field plays a significant role in all the physical quantities. Analytical solutions based upon normal mode analysis for themoelastic problem in solids have been developed and utilized. The value of all the physical quantities converges to zero with an increase in distance and all functions are continuous. The numerical results using normal mode analysis are in a good agreement with previously published data. Indeed, the most important factor is the extent to penetrate the magnetic field and the arrival of magnetic energy into the cells of the body in the areas sought to be accessed, and the duration of exposure to the magnetic field is very important, as the intensity of penetration of the field directly related to the mass of material magnetic-generating field and how to focus the magnetic field and treatment domain magnetic helps to get rid of pain in general also helps diabetics to lower blood sugar and cholesterol in the blood and reduce the swelling and works to calm the nervous system
- R. Seip and E. Ebbini, “Studies on the three-dimensional temperature response to heating fields using diagnostic ultrasound,” Transactions on Biomedical Engineering, Vol. 42, pp. 828–839, 1995. [Publisher]
- A. M. Stoll, “Thermal properties of human skin related to nondestructive measurement of epidermal thickness,” Journal of Investigative Dermatology, Vol. 69, No. 3, pp. 328–332, Sep. 1977, https://doi.org/10.1111/1523-1747.ep12507865 [Publisher]
- Stoll Am, Chianta Ma, and Piergallini Jr, “Thermal conduction effects in human skin,” Aviation, Space, and Environmental Medicine, Vol. 50, No. 8, pp. 778–787, Aug. 1979. [Search CrossRef]
- C. J. Diederich and K. Hynynen, “Ultrasound technology for hyperthermia,” Ultrasound in Medicine and Biology, Vol. 25, No. 6, pp. 871–887, Jul. 1999, https://doi.org/10.1016/s0301-5629(99)00048-4 [Publisher]
- D. R. Poirier and G. H. Geiger, “Conduction of heat in solids,” in Transport Phenomena in Materials Processing, Cham: Springer International Publishing, 2016, pp. 281–327, https://doi.org/10.1007/978-3-319-48090-9_9 [Publisher]
- J. Lienhard, A Heat Transfer Textbook. Englewood Cliffs: Prentice-Hall, 1987. [Search CrossRef]
- Pere J. Riu, Kenneth R. Foster, Dennis W. Blick, and Eleanor R. Adair, “A thermal model for human thresholds of microwave-evoked warmth sensations,” Bioelectromagnetics, Vol. 18, No. 8, pp. 578–583, Dec. 1998, https://doi.org/10.1002/(sici)1521-186x(1997)18:8 [Publisher]
- H. F. Bowman, E. G. Cravalho, and M. Woods, “Theory, measurement, and application of thermal properties of biomaterials,” Annual Review of Biophysics and Bioengineering, Vol. 4, No. 1, pp. 43–80, Jun. 1975, https://doi.org/10.1146/annurev.bb.04.060175.000355 [Publisher]
- G. Martin, H. Bowman, and W. Newman, “Basic element method for computing the temperature field during hyperthermia therapy planning,” Advanced Bio Heat Mass Transfer, Vol. 231, pp. 75–80, 1992. [Search CrossRef]
- C. R. Davies, G. M. Saidel, and H. Harasaki, “Sensitivity analysis of one-dimensional heat transfer in tissue with temperature-dependent perfusion,” Journal of Biomechanical Engineering, Vol. 119, No. 1, pp. 77–80, Feb. 1997, https://doi.org/10.1115/1.2796068 [Publisher]
- H. Arkin, L. X. Xu, and K. R. Holmes, “Recent developments in modeling heat transfer in blood perfused tissues,” IEEE Transactions on Biomedical Engineering, Vol. 41, No. 2, pp. 97–107, 1994, https://doi.org/10.1109/10.284920 [Publisher]
- B. Erdmann, J. Lang, and M. Seebass, “Optimization of temperature distributions for regional hyperthermia based on a nonlinear heat transfer modela,” Annals of the New York Academy of Sciences, Vol. 858, No. 1 BIOTRANSPORT, pp. 36–46, Sep. 1998, https://doi.org/10.1111/j.1749-6632.1998.tb10138.x [Publisher]
- P. D. Tyréus and C. J. Diederich, “Theoretical model of internally cooled interstitial ultrasound applicators for thermal therapy,” Physics in Medicine and Biology, Vol. 47, No. 7, pp. 1073–1089, Apr. 2002, https://doi.org/10.1088/0031-9155/47/7/306 [Publisher]
- K. R. Diller, “Development and solution of finite-difference equations for burn injury with spreadsheet software,” Journal of Burn Care and Rehabilitation, Vol. 20, No. 1 Pt 1, pp. 25–32, Jan. 1999, https://doi.org/10.1097/00004630-199901001-00005 [Publisher]
- K. R. Diller, “Modeling thermal skin burns on a personal computer,” Journal of Burn Care and Rehabilitation, Vol. 19, No. 5, pp. 420–429, Sep. 1998, https://doi.org/10.1097/00004630-199809000-00012 [Publisher]
- S. C. Jiang, N. Ma, H. J. Li, and X. X. Zhang, “Effects of thermal properties and geometrical dimensions on skin burn injuries,” Burns, Vol. 28, No. 8, pp. 713–717, Dec. 2002, https://doi.org/10.1016/s0305-4179(02)00104-3 [Publisher]
- C. L. Chan, “Boundary element method analysis for the bioheat transfer equation,” Journal of Biomechanical Engineering, Vol. 114, No. 3, pp. 358–365, Aug. 1992, https://doi.org/10.1115/1.2891396 [Publisher]
- B. Mochnacki and E. Majchrzak, “Sensitivity of the skin tissue on the activity of external heat sources,” Computer Modeling in Engineering and Sciences, Vol. 4, No. 3-4, pp. 431–438, Jun. 2003. [Search CrossRef]
- M. Othman and Kh. Lotfy, “Two-dimensional problem of generalized magneto-thermoelasticity under the effect of temperature dependent properties for different theories,” Multidiscipline Modeling in Materials and Structures, Vol. 5, pp. 235–242, 2009. [Publisher]
- M. Othman and Kh. Lotfy, “On the plane waves in generalized thermo – microstretch elastic half-space,” International Communication in Heat and Mass Transfer, Vol. 37, pp. 192–200, 2010. [Publisher]
- M. A. Othman, K. Lotfy, and R. M. Farouk, “Transient disturbance in a half-space under generalized magneto-thermoelasticity with internal heat source,” Acta Physica Polonica A, Vol. 116, No. 2, pp. 185–192, Aug. 2009, https://doi.org/10.12693/aphyspola.116.185 [Publisher]
- M. Othman and Kh. Lotfy, “Effect of magnetic field and inclined load in micropolar thermoelastic medium possessing cubic symmetry,” International Journal of Industrial Mathematics, Vol. 1, pp. 87–104, 2009. [Search CrossRef]
- K. Lotfy and M. I. A. Othman, “Effect of rotation on plane waves in generalized thermo-microstretch elastic solid with a relaxation time,” Meccanica, Vol. 47, No. 6, pp. 1467–1486, Aug. 2012, https://doi.org/10.1007/s11012-011-9529-7 [Publisher]