Research Article | Open Access
Computational Dynamics of Arterial Blood Flow in the Presence of Magnetic Field and Thermal Radiation Therapy
We conduct a numerical study to determine the influence of magnetic field and thermal radiation on both velocity and temperature distributions in a single blood vessel. The model here assumes that blood is a Newtonian incompressible conducting fluid with radially varying viscosity due to hematocrit variation. The transient equations of momentum and energy transport governing the flow in an axisymmetric configuration are solved numerically using a semi-implicit finite difference method. Results are presented graphically and discussed both qualitatively and quantitatively from the physiological point of view. The results of this work may enhance current understanding of the factors that determine the effects of hyperthermia treatment on tumor tissues.
Blood flow in a large blood vessel has a profound influence on the efficiency of thermal therapy treatment [1–3]. In pathological situations, thermal radiation therapy is one of the treatments employed by medical practitioners [4–6]. The procedure involves transmitting heat below the skin surface into tissues and muscles. Deep heat speeds up healing by increasing blood flow to the injury. Electromagnetic heat, such as shortwaves and microwaves, sends heat up to 2 inches into the tissue and muscles. It works best for injuries in joints, muscles, and tendons. Heat therapy may help reduce pain. Moreover, hyperthermia treatment has been demonstrated as effective during cancer therapy in recent years. Its objective is to raise the temperature of pathological tissues above cytotoxic temperatures (41–45°C) without overexposing healthy tissues [7, 8]. Temperature distribution within tissues primarily depends on tissue thermal conductivity, the heating source’s power deposition pattern characteristics, and heat transfer resulting from blood flow [9–11]. An important source of temperature nonuniformity is the presence of large vessels entering the heated volume and carrying blood at a lower systemic temperature (37°C). The design of delivered power devices and numerous theoretical, experimental, and clinical studies have demonstrated that large blood vessels may produce localized cooling regions within heated tissues during hyperthermia treatment [12, 13]. Kolios et al.  demonstrated the effect of large blood vessel in heated tissues and showed that the dissipation of heat from heated tissues was carried out by convection through blood flow and also by conduction process. A numerical study on the impact of large vessels on the temperature uniformity during hyperthermia treatment assuming steady-state condition was conducted by Creeze and Lagendijk . They reported that the presence of a large vessel may result in nonuniform temperature resulting in possible underdosage. Cho and Hyun  assumed a sinusoidal variation of the velocity at the pipe inlet in a numerical study of pulsatile flow and heat transfer characteristics within a pipe. Seo Young Kim et al.  analyzed numerically the heat transfer characteristics of fully developed pulsatile flow in a channel, assuming a sinusoidal variation of the velocity at the inlet of the channel. In recent time, several authors have also investigated the flow structure and heat transfer characteristics of Newtonian and non-Newtonian fluid models in order to understand the physiological flow systems [18–21]. However, in the large blood vessel it is a good approximation to consider whole blood as a Newtonian viscous fluid .
So far, the theoretical studies dealing with the combined effects of thermal radiation absorption and magnetic field on the transient flow of blood in a large artery with viscous dissipation have received very little attention. It was observed that the heart rate decreases by exposing biological systems to an external magnetic field. The electrocardiography (ECG) pattern taken in the presence of a magnetic field not only provides information on blood flow but also offers a new noninvasive method of studying the cardiac performance . Moreover, the red blood cell (RBC) is a major biomagnetic substance, and the blood flow may be influenced by the magnetic field. The magnetic properties of RBCs play important roles in the increase in blood viscosity during exposure to a static magnetic field. Motivated by the above studies, our aim is to investigate the combined effects of thermal radiation and magnetic field on transient flow of blood in a large blood vessel. In Section 2 the mathematical formulation of the problem is established. The numerical solutions to the problem are obtained in Section 3. Qualitative and quantitative discussions of the solutions are given in Section 4.
2. Formulation of the Problem
Consider unsteady, axisymmetric incompressible Newtonian blood flow through a rigid vessel of length and radius . The geometry of the problem and the coordinate system are shown in Figure 1.
For the development of mathematical model, the following assumptions are made.(i)The flow is incompressible so that the simple equation of state for a Boussinesq fluid holds in gravitational field acting in reverse direction to .(ii)Viscosity of blood varies radially with maximum magnitude along the arterial centerline due to the accumulation of red blood cells (RBC) in this region.(iii)The temperature difference of the blood and artery is high enough for radiative heat transfer to be valid.(iv)Magnetic field applied externally is of constant strength.(v)The electromagnetic force produced and the electrical conductivity are very small.(vi)The flow through the intimal and media layers and the deformation of blood vessels are neglected.
Under the above mentioned assumption, the one-dimensional equations of motion and energy balance for the unsteady and axisymmetric flow of blood through an artery provided are [22–25] where is the blood axial velocity, is the blood temperature, is the arterial wall temperature, is the blood reference temperature, is the gravitational acceleration, is the density, is the thermal expansion coefficient, is the intensity of magnetization, is the electric conductivity of the blood, is the radius of the normal blood vessel, is the blood thermal conductivity, and is the blood pressure. The most general differential approximation for the radiative flux as given by Cheng  is where is the absorption coefficient and is the Stefan-Boltzman constant. For the optically thin case that we adopt in the present study and an approximate form of the radiative flux equation could be taken as [1, 2, 26] Following Lih , the transverse variation in the blood viscosity is taken as where (a constant), and is the hematocrit function given as where is the maximum hematocrit at the center of the blood vessel and (≥2) is a parameter determining the exact shape of the profile. Equation (5) is valid only for dilute suspension of red cells which are assumed spherical [14, 15]. We introduce the following dimensionless quantities: Substituting (5) into (1), we obtain the following dimensionless governing equations: where is the radiation absorption parameter, is the Grashof number, is the temperature difference parameter, is the axial pressure gradient parameter, is the magnetic field parameter, is the Brinkman number, is the Reynolds number, is the Prandtl number, , are hematocrit constants, and is the hematocrit shape parameter. The volumetric flow of blood in the artery is given by The dimensionless shear stress () and the rate of heat transfer () at the arterial wall are given by In the following section, (7) are solved numerically and the skin-friction together with the arterial wall heat transfer rate is computed.
3. Numerical Solution
Our numerical algorithm is based on the semi-implicit finite difference scheme given in  for the isothermal viscoelastic case. As in [21, 28–31], we extend the algorithm to the temperature equation and take the implicit terms at the intermediate time level , where . The algorithm employed in  uses ; we will, however, follow the formulation in [20, 28–31] and thus take in this paper so that we can use larger time steps. The discretization of the governing equations is based on a linear Cartesian mesh and uniform grid on which finite differences are taken. We approximate both the second and first spatial derivatives with second-order central differences. The equations corresponding to the first and last grid point are modified to incorporate the boundary conditions. The semi-implicit scheme for the velocity component reads The equation for then becomes where . The solution procedure for thus reduces to inversion of tridiagonal matrices which is an advantage over a full implicit scheme. The semi-implicit integration scheme for the temperature equation is similar to that for the velocity component. Unmixed second partial derivatives of the temperature are treated implicitly: The equation for thus becomes where . The solution procedure again reduces to inversion of tridiagonal matrices. The schemes (11) and (13) were checked for consistency. For , these are first order accurate in time but second order in space. The schemes in  have which improves the accuracy in time to second order. We use here so that we are free to choose larger time steps and still converge to the steady solutions. The algorithm was also tested for both spatial and temporal convergence and shown to be independent of both mesh size and time step size.
4. Results and Discussion
Unless otherwise stated, we employ the parameter values: , , , , , , , , , , , , , and .
4.1. General Profiles of Flow Quantities
In the early stage of thermal therapy treatment (say ), it is interesting to note that the blood velocity profile in a large vessel converges to a standard parabolic profile much earlier than does the temperature profile as illustrated in Figure 2. This is because large vessels exhibit large temperature gradients than smaller vessels. As such, the dissipation of heat from heated tissues, which is carried out by convection and conduction, is greater for large vessels due to a strong cooling effect. Consequently, the cooling effect of the blood flow in large vessels may result in underdosage of select tumor parts, since part of the surrounding tumor may not reach the desired treatment temperature. However, at much later time of thermal dosage administration (say ), the blood temperature profile increases transversely towards the arterial centerline in the same manner as the velocity profile (see Figure 3). This enables the surrounding tumor parts to attain the desired treatment temperature.
We need to point out early on that the parameter will need to be carefully controlled. Physically, the thermal radiation therapy is administered on the patient for time . However, the therapy dose needs to be controlled in order to prevent development of hot spots (represented by blowup) in the flow system; see Figure 4 which plots the steady state maximum blood temperature against each value of . Occurrence of such blowup phenomena may severely damage normal body tissue.
Significantly large dosage of radiation therapy (i.e., large ) would lead to finite time temperature blowup since the terms associated with are strong heat sources as with, say, the Frank-Kamenetskii parameter in reacting flows . At our current material values, we notice that such blowup of solutions occurs below, say, . Our choices of in this paper will thus be informed by Figure 4.
4.2. Transient Solutions
As an immediate followup to the results of Figures 2 and 3, we display the transient solutions in Figure 5. As expected, Figure 5 shows that the blood velocity increases with time during the thermal therapy treatment.
As explained earlier, Figure 5 shows a progressive increase of arterial blood temperature with increasing time, until a steady state is reached which depicts an expected parabolic profile. The time evolution of the velocity and temperature with varying values of the radiation parameter along the arterial centerline is illustrated in Figure 6.
Since, as mentioned earlier, the radiation acts as a heat source within the blood, the arterial blood temperature (and hence also that along the centerline) should gradually increase with increasing radiation dosage; this is well illustrated in Figure 6. The resultant increase in temperature correspondingly increases the Boussinesq source terms in the momentum equation and hence also increases the blood velocity. This increase of velocity with radiation dosage is illustrated in Figure 6.
4.3. Velocity Dependence on Magnetic Field and Hematocrit Shape
Increasing the damping magnetic force predictably decreases the blood velocity as alluded to earlier and illustrated in Figure 7. A change in the hematocrit shape from the profile also decreases the velocity. This is expected since higher values of correspond to increased blood viscosity and hence reduced velocity.
4.4. Wall Shear Stress and Wall Heat Transfer
The wall shear stress dependence on hematocrit shape is illustrated in Figure 9 for varying values of the radiation parameter. Similarly, Figure 10 shows the wall shear stress dependence on the magnetic field for varying values of .
Since increased values of correspond to increased viscosities and thus reduced velocities, the magnitude of the slope of the velocity correspondingly decreases (with increasing ) and hence also the wall shear stress, Figure 9. Increases in the thermal radiation dosage would necessarily also increase the blood temperature as mentioned earlier. Thus, for fixed , increasing the radiation dosage () would indirectly increase the blood velocity and hence also the wall shear stress as illustrated in Figure 9. This is so since the increase in increases the blood temperature and hence also the source terms in the momentum equation.
The velocity profiles have been shown to decrease with increasing magnetic field force, . This in turn means that the magnitude of the slope of the velocity near the wall would also decrease with increasing , hence explaining the observed decrease in wall shear stress with increasing magnetic field force in Figure 10. At each value of the magnetic field, an increase in the Grashof number would directly increase the velocity via the increased Boussinesq source terms and hence also the magnitude of the slope of the velocity at the wall.
The variation of the wall heat transfer with hematocrit shape for varying radiation parameter values is shown in Figure 11 and is similarly explained as with Figure 9. In particular, the velocity (and hence also the magnitude of the velocity gradient at the wall) decreases with increasing . This in turn leads to a decrease in the viscous heating terms in the temperature equation leading to reduced blood temperatures. The decrease in the blood temperature (and hence also the magnitude of the temperature gradient at the wall) then leads to the corresponding decrease in wall heat transfer (with increasing ) observed in Figure 11. At fixed , increases in directly lead to corresponding increases in the blood temperature as already noted. The increase in blood temperature (and hence also the magnitude of the temperature gradient at the wall) then leads to the corresponding increase in wall heat transfer (with increasing ) observed in Figure 11.
The variation of the wall heat transfer with Brinkman number for varying values of the Grashof number is shown in Figure 12. Increases in the Brinkman number would correspondingly increase the magnitude of the viscous heating terms and thus also increase the arterial blood temperature. This leads to increases in the magnitudes of the temperature gradients at the wall and hence the observed increases in the wall heat transfer; see Figure 12. For each fixed Brinkman number, we again refer to the increase in the velocity with increasing Grashof number. This in turn similarly increases the magnitude of the viscous heating terms for each fixed value of the Brinkman number leading to the increased wall heat transfer shown in Figure 12.
We computationally investigate the influence of magnetic field force and thermal radiation on the dynamics of arterial blood with transversely varying viscosity and subjected to axisymmetric one-dimensional pressure driven flow under Boussinesq approximations. We observe that there is a transient increase in the blood velocity with an increase in the radiation absorption. Such increases in blood velocity and hence wall shear stress together with the possible development of hot spots in the flow system mean that the therapy dose needs to be carefully controlled. We also note that the blood velocity decreases with an increase in magnetic field intensity.
Conflict of Interests
The authors declare that there is no conflict of interests.
- A. Ogulu and A. R. Bestman, “Deep heat muscle treatment a mathematical model—I,” Acta Physica Hungarica, vol. 73, no. 1, pp. 3–16, 1993.
- A. Ogulu and A. R. Bestman, “Deep heat muscle treatment a mathematical model—II,” Acta Physica Hungarica, vol. 73, no. 1, pp. 17–27, 1993.
- T.-C. Shih, H.-L. Liu, and A. T.-L. Horng, “Cooling effect of thermally significant blood vessels in perfused tumor tissue during thermal therapy,” International Communications in Heat and Mass Transfer, vol. 33, no. 2, pp. 135–141, 2006.
- S. A. Sapareto and W. C. Dewey, “Thermal dose determination in cancer therapy,” International Journal of Radiation Oncology Biology Physics, vol. 10, no. 6, pp. 787–800, 1984.
- C. Damianou and K. Hynynen, “Focal spacing and near-field heating during pulsed high temperature ultrasound therapy,” Ultrasound in Medicine and Biology, vol. 19, no. 9, pp. 777–787, 1993.
- J. van der Zee, “Heating the patient: a promising approach?” Annals of Oncology, vol. 13, no. 8, pp. 1173–1184, 2002.
- W. Levin, M. D. Sherar, B. Cooper, R. P. Hill, J. W. Hunt, and F.-F. Liu, “Effect of vascular occlusion on tumour temperatures during superficial hyperthermia,” International Journal of Hyperthermia, vol. 10, no. 4, pp. 495–505, 1994.
- R. B. Roemer, “The local tissue cooling coefficient: a unified approach to thermal washout and steady-state perfusion calculations,” International Journal of Hyperthermia, vol. 6, no. 2, pp. 421–430, 1990.
- H. H. Pennes, “Analysis of tissue and arterial blood temperatures in the resting human forearm,” Journal of Applied Physiology, vol. 1, no. 2, pp. 93–122, 1948.
- O. I. Craciunescu and S. T. Clegg, “Pulsatile blood flow effects on temperature distribution and heat transfer in rigid vessels,” Journal of Biomechanical Engineering, vol. 123, no. 5, pp. 500–505, 2001.
- R. Siegel and M. Perlmutter, “Heat transfer for pulsating laminar duct flow,” Journal of Heat Transfer, vol. 8, pp. 111–116, 1962.
- J. C. Chato, “Heat transfer to blood vessels,” Journal of Biomechanical Engineering, vol. 102, no. 2, pp. 110–118, 1980.
- T. C. Shih, H. S. Kou, and W. L. Lin, “Effect of effective tissue conductivity on thermal dose distributions of living tissue with directional blood flow during thermal therapy,” International Communications in Heat and Mass Transfer, vol. 29, no. 1, pp. 115–126, 2002.
- M. C. Kolios, A. E. Worthington, D. W. Holdsworth, M. D. Sherar, and J. W. Hunt, “An investigation of the flow dependence of temperature gradients near large vessels during steady state and transient tissue heating,” Physics in Medicine and Biology, vol. 44, no. 6, pp. 1479–1497, 1999.
- J. Creeze and J. J. W. Lagendijk, “Temperature uniformity during hyperthermia: the impact of large vessels,” Physics in Medicine and Biology, vol. 37, no. 6, pp. 1321–1337, 1992.
- H. W. Cho and J. M. Hyun, “Numerical solutions of pulsating flow and heat transfer characteristics in a pipe,” International Journal of Heat and Fluid Flow, vol. 11, no. 4, pp. 321–330, 1990.
- S. Y. K. Seo Young Kim, B. H. K. Byung Ha Kang, and J. M. H. Jae Min Hyun, “Heat transfer in the thermally developing region of a pulsating channel flow,” International Journal of Heat and Mass Transfer, vol. 36, no. 17, pp. 4257–4266, 1993.
- J. Prakash, O. D. Makinde, and A. Ogulu, “Magnetic effect on oscillatory blood flow in a constricted tube,” Botswana Journal of Technology, vol. 13, no. 1, pp. 45–50, 2004.
- O. D. Makinde, “Analysis of non-newtonian reactive flow in a cylindrical pipe,” Journal of Applied Mechanics, Transactions ASME, vol. 76, no. 3, pp. 1–5, 2009.
- T. Chinyoka, Y. Y. Renardy, M. Renardy, and D. B. Khismatullin, “Two-dimensional study of drop deformation under simple shear for Oldroyd-B liquids,” Journal of Non-Newtonian Fluid Mechanics, vol. 130, no. 1, pp. 45–56, 2005.
- T. Chinyoka, “Computational dynamics of a thermally decomposable viscoelastic lubricant under shear,” Journal of Fluids Engineering, Transactions of the ASME, vol. 130, no. 12, Article ID 121201, 7 pages, 2008.
- T. J. Pedley, Fluid Dynamics of Large Blood Vessel, Cambridge Press, Cambridge, UK, 1980.
- L. Pauling and C. D. Coryell, “The magnetic properties and structure of hemoglobin, oxyhemoglobin and carbonmonoxyhemoglobin,” Proceedings of the National Academy of Sciences of the United States of America, vol. 22, pp. 210–216, 1936.
- A. R. Rao and K. S. Deshikachar, “Physiological-type flow in a circular pipe in the presence of a transverse magnetic field,” The Journal of the Indian Institute of Science, vol. 68, pp. 247–260, 1988.
- K. Haldar and S. N. Ghosh, “Effect of a magnetic field on blood flow through an indented tube in the presence of erythrocytes,” Indian Journal of Pure and Applied Mathematics, vol. 25, no. 3, pp. 345–352, 1994.
- P. Cheng, “Two dimensional radiating gas flow by a moment method,” The American Institute of Aeronautics and Astronautics, vol. 2, pp. 1662–1664, 1964.
- M. M. Lih, Transport Phenomenon in Medicine and Biology, John Wiley and Sons, New York, NY, USA, 1969.
- T. Chinyoka, “Poiseuille flow of reactive Phan-Thien-Tanner liquids in 1D channel flow,” Journal of Heat Transfer, vol. 132, no. 11, Article ID 111701, 2010.
- T. Chinyoka, “Suction-injection control of shear banding in non-isothermal and exothermic channel flow of Johnson-Segalman liquids,” Journal of Fluids Engineering, Transactions of the ASME, vol. 133, no. 7, Article ID 071205, 2011.
- I. E. Ireka and T. Chinyoka, “Non-isothermal flow of a Johnson-Segalman liquid in a lubricated pipe with wall slip,” Journal of Non-Newtonian Fluid Mechanics, vol. 192, pp. 20–28, 2013.
- T. Chinyoka, S. Goqo, and B. I. Olajuwon, “Computational analysis of gravity driven flow of a variable viscosity viscoelastic fluid down an inclined plane,” Computers & Fluids, vol. 84, no. 15, pp. 315–326, 2013.
Copyright © 2014 T. Chinyoka and O. D. Makinde. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.