# Magnetohydrodynamic Blood Flow in a Cylindrical Tube with Magnetic Particles: A Time Fractional Model

Authors: Farhad Ali; Samina Majeed; Anees Imtiaz
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The present work theoretically investigates the natural convection blood flow as a Brinkman-type fluid with uniformly distributed magnetic particles in a circular cylindrical tube with the applied external magnetic field. The classical model for the blood flow is generalized by using the definition of Caputo time-fractional derivative. The exact solutions are obtained by using the Laplace and Henkel transforms. Unlike the classical model, the obtained general results are expressed in the form of “Lorenzo and Hartley’s” and “Robotnov and Hartley’s” functions. Graphs are plotted to show the effects of different parameters on the blood flow. Furthermore, the velocity and temperature distributions are discussed in terms of memory. The effect of fractional parameterα for a long and short time has also been observed. It is noticed that blood velocity can be controlled using the fractional parameter. It is also found that, for τ>0, fluid and particles motion increased, and reverse behavior is observed for τ<0. It has been noticed that increasing values of particle mass parameter Pm and magnetic parameter M slow down the motion of blood and magnetic particles. These results are helpful for effective drug delivery and regulating blood flow.


---

## Body

## 1. Introduction

Biomagnetic fluid dynamics (BFD) studies the magnetic effects on biofluids. BFD is a new area of attraction for researchers [1–3]. This new field has many applications in biomedical and engineering [4–7]. Blood is one of the essential biofluids, which is also biomagnetic [8–10]. Since many diseases, including cancer, there is a need to enhance the therapy methods which are less harmful and low cost. Therefore, scientists have made successful use of magnetic particles for this purpose, for example, hyperthermia where magnetic particles are transferred to the cancer cells and heated by the use of an external magnetic field. This helps to destroy the tumor cells without harming the normal tissues [11]. Magnetic particles are utilized in MRI, which is used for diagnostics of tumor cells [12]. Similarly, the use of magnetic particles is helpful in drug transport and bioseparation [13]. Due to wide applications of magnetic particles for the diagnostic and treatment of diseased cells, many researchers have developed mathematical models to observe the effects of magnetic particles in the blood. Tzirtzilakis [14] observed Newtonian fluid’s flow behavior in a rectangular duct under magnetic field application. Sharma et al. [9] studied the transverse magnetic field’s impact on blood flow containing magnetic particles in a cylindrical tube. Ghasemi et al. [15] analyzed the two-phase flow of blood with uniformly distributed nanoparticles through arteries in the existence of the magnetic field. The 2D flow of blood through a porous artery with stenosis is examined by Ellahi et al. [16]. Bhatti et al. [17] examined the MHD flow of blood containing microscopic size metallic particles in a hollow cylinder. Variation in blood flow by using radiation in a small magnetized vessel is discussed by Sinha and Shit [18]. Khan et al. [19] investigated the magnetic field and heat transfer effects in Casson fluid flow. Their observation reported the reduction of fluid flow due to the increase in magnetic parameter. Gul et al. [20] analytically and numerically analyzed Casson blood flow with carbon nanotubes considering Marangoni convection and magnetic effects.The joint applications of MHD and natural convection in non-Newtonian fluids have attracted many researchers due to vast applications in the real world, from industries to biological sciences [21, 22]. Fluid motion due to natural sources, i.e., buoyancy forces, is known as free convection or natural convection. Buoyancy forces arise because of the difference in densities. Rahbari et al. [23] investigated analytically and numerically the influence of MHD and natural convection fluid flow in a parallel plate channel. The influence of free convection heat transfer on the velocity of blood in a vessel is reported by Parkash and Ogulu [24]. Akbar [25] worked on the effects of mixed convection on blood flow in a stenosed artery.In the present work, blood is taken as a Brinkman-type fluid. Brinkman-type fluid is a special kind of viscous fluid derived from the flow past on a high permeable medium [26–28]. Brinkman gave the idea of calculation of viscous forces which are exerted by fluid [29]. Ali et al. [30], with the help of Laplace transform, obtained the exact solutions of the Brinkman-type fluid model. The exact solution for Brinkman-type fluid with the help of Fourier sine transform is obtained by Fetecau et al. [31]. Anwar et al. [32]investigated Brinkman-type fluid model to inspect the enhancement of heat transfer due to iron oxide nanoparticles. They considered water and kerosene oil as base fluid. It was observed that, in water, nanoparticles efficiently worked rather than kerosene oil. Saqib et al. [33] studied the blood as a Brinkman type with composite magnetic particles in the existence of a perpendicularly applied magnetic field. Kumar et al. [34] reported the impact of radiation on the natural convection flow of Brinkman-type nanofluid on a vertical plate under the application of a magnetic field.The subject of fractional calculus was introduced in 1965. It studies the derivative and integral with the order that is fractional [35]. Fractional calculus contains derivatives and integrals that are not local, making it much attractive for researchers. Fractional calculus briefly describes the memory effects [36, 37]. This memory helps to observe the past process. These properties make fractional calculus dominant over classical calculus. The research on fractional-order PDEs has shown that the results are comparatively more accurate and general than integer order [33]. Fractional-order derivative and integral has vast applications in science and technology and is a highly emerging area in calculus, having application in biochemistry and fluid mechanics [38–41]. Fractional calculus is used in bioengineering, mathematics, finance, and disease dynamics [42, 43]. Imtiaz et al. [44] proposed a fractional model for blood flow containing metallic particles in a cylindrical tube. Singha and Sengupta [40] used fractional derivative to observe fluid flow with dusty particles. The effect of MHD on blood flow passing through human artery is observed by Sud and Sekhon [45]. Ahmad et al. [46] examined Jeffery fluid flow over exponentially of an infinite length using Caputo time-fractional derivative and Caputo–Fabrizio fractional model. In this study, dual behavior for velocity is observed by the change in the memory parameter. A time-fractional model that deals with cancer therapy is developed by Asjad [47]. Butt et al. [48] used the fractional derivative to model MHD Jeffery fluid over a vertical plane with heat and mass transfer. It was found that temperature, velocity, and concentration can be enhanced by decreasing the value of the fraction parameter.In the view of all the literature mentioned above, no one has considered the role of natural convection and heat transfer in blood flow. The present study reports the exact solutions of the unsteady free convection blood flow as a Brinkman-type fluid using the transverse magnetic field. Furthermore, the magnetic particles with different velocities are evenly distributed in the blood, passing through a heated cylindrical tube. The fluid contains blood as a liquid phase and particles in a solid phase. Flow is considered a two-phase flow like a liquid, and solid phases are flowing with different velocities. The classical fractional derivative is replaced by Caputo time-fractional derivative to obtain the fractional model. The Laplace and Finite Hankel transformations are used to getting velocities and temperature distributions for blood flow and suspended particles.


## 2. Problem Formulation

Consider the blood flow as a Brinkman-type fluid in a horizontal heated cylinder with radiusr, as shown in Figure 1. The motion of the fluid is taken in the z-direction. It is also considered that blood contains evenly distributed magnetic particles. Flow is taken as a two-phase flow. Therefore, momentum equations for particles and blood are defined separately. The transverse magnetic field is applied along with heat transfer through natural convection. Smooth flow of fluid is considered, which implies that the Reynolds number is very low, and the induced magnetic field can be neglected [44]. At time t=0, fluid and cylinder both are at rest with ambient temperature T∞. At the time, t=0+, the fluid starts motion with velocity U0, and temperature varies from T∞ to Tw.Figure 1
The geometry of the problem.The developed model contains momentum equation for fluid flow, obtained by Navier–Stokes equations, Maxwell’s set of relations for the magnetic field, and Newton’s second law for the particle velocity.Maxwell set of equations are defined as [49, 50](1)∇⟶×E⟶=−∂B⟶1∂t,divB⟶1=0,∇⟶×B⟶1=μ0J⟶.According to Ohm’s law [51, 52],(2)J⟶=σ0E⟶+V⇀×B⟶,where σ0 is electrical conductivity and V⇀ is the velocity field. The electromagnetic force is defined as(3)J⟶×B⟶=−σB02U.The momentum equation for blood flow is defined as [33, 53](4)∂ur,t∂t+βur,t=−1ρ∂p∂z+ν∂2ur,t∂r2+1r∂ur,t∂r+KNρupr,t−ur,t+gβTT−T∞−1ρσB02ur,t.The oscillating pressure gradient [54] is defined as(5)−∂p∂z=λ0+λ1Cosωt.Incorporating equation (5) in (4), (6) becomes(6)∂ur,t∂t+βur,t=1ρλ0+λ1Cosωt+ν∂2ur,t∂r2+1r∂ur,t∂r+KNρupr,t−ur,t+gβTT−T∞−1ρσB02ur,t,where KN/ρupr,t−ur,t indicates the force between particles and fluid.The momentum equation for magnetic particles is defined as [55](7)m∂upr,t∂t=Kur,t−upr,t.The equation for energy is as follows(8)ρCp∂T∂t=k∂2T∂r2+1r∂T∂r,t>0,r∈0,r0,subjected to the following initial and boundary conditions:(9)ur,0=0,upr,0=0,Tr,0=T∞,ur0,t=HtU0,upr0,t=HtU01−e−K/mt,Tr0,t=Tw,∂u∂r|r=0=0,using the following dimensionless variables,(10)f=uu0,g=upu0,τ=tνr02,r1=rr0,A0=λ0r02μu0,A1=λ1r02μu0,θ=T−T∞Tw−T∞.The dimensionless form of the system with initial and boundary conditions are(11)∂fr1,τ∂τ=A0+A1Cosωt+∂2fr1,τ∂r12+1r1∂fr1,τ∂r1+Pcgr1,τ−fr1,τ−M0fr1,τ+Grθr1,τ,(12)∂gr1,τ∂τ=Pmfr1,τ−gr1,τ,(13)∂θr1,τ∂τ=1Pr∂2θr1,τ∂r12+1r1∂θr1,τ∂r1,τ>0,r1∈0,1,(14)fr1,0=0,gr1,0=0,θr1,0=0,f1,τ=1,g1,τ=1−e−bt,θ1,t=1,∂fr1,τ∂r1|r1=0=0,where(15)Pc=KNr02μ,M=σB0r02μ,β1=βr02ν,Pm=νKr02,M0=M−β1,Gr=gr02βTTw−T∞u0ν,Pr=μCpk,b=Kr02mν,where Pcis the dimensionless particle concentration parameter.Now, applying Caputo time derivative to equations (11)–(13), we obtain(16)Dtαfr1,τ=A0+A1Cosωt+∂2fr1,τ∂r12+1r1∂fr1,τ∂r1+Pcgr1,τ−fr1,τ−M0fr1,τ+Grθr1,τ,(17)Dtαgr1,τ=Pmfr1,τ−gr1,τ,(18)Dtαθr1,τ=1Pr∂2θr1,τ∂r12+1r1∂θr1,τ∂r1,where(19)Dtαfr,τ=1Γ1−α∫0tf′τt−ταdτ,0<α<1,which is defined as Caputo time-fractional derivative [56].


## 3. The Solution of the Problem

The Laplace and Henkel transform is employed to get the required solution for the blood and particle velocities. Both transforms are defined as [33](20)Lfr,tq=f¯q=∫0∞fr,te−qtdt,Hf¯r,qrn=f¯Hrn,q=∫01rf¯r,tJ0rrndr,where J0rn is Bessel’s function of the first kind of order zero and rn are positive roots of the equation J0x=0.


### 3.1. Temperature Field Calculation

Applying Laplace transform to equations (14) and (18), obtained equations are(21)qαθ¯r1,q=1Prd2θ¯r1,qdr12+1r1dθ¯r1,qdr1,(22)θ¯1,q=1q,where θ¯r1,q is Laplace transform of θr1,t, where q denotes transformation variable.By the using Finite Hankel transform to equation (21) and putting values from equation (22), we obtain(23)θ¯Hr1n,q=r1nJr1nPr1qqα+a1,where r1n2/Pr=a1.θ¯Hr1n,q shows Hankel transform of θ¯r,q. Equation (23) reduces to(24)θ¯Hr1n,q=J1r1nr1n1q−J1r1nr1nq−1−αqα+a1.Applying inverse Laplace transform by using Lorenzo and Hartley’sRβ0,υ−λ,τ functions [57] to equation (24),(25)θHr1n,τ=J1r1nr1n−J1r1nr1nRα,1−ατ,−a1,L−1s−υsβ0+λ=Rβ0,υ−λ,τ=∑n=0∞−λnτn+1β0−1−υΓn+1β0−υ.Applying the inverse Henkel transform to equation (25), we obtain(26)θr1,τ=1−2∑n=1∞J0r1r1nr1nJ1r1nRα,1−α−a1,τ.


### 3.2. The Calculation for Blood Flow

Applying Laplace transform to equations (16) and (17) and using equation (13), we obtain(27)qαf¯r1,q=A0q+A1qq2+ω2+∂2f¯r1,q∂r12+1r1∂f¯r1,q∂r1+Pcg¯r1,q−f¯r1,q−M0f¯r1,q+Grθ¯r1,q,(28)g¯r1,q=f¯r1,qPmqα+1.Applying the Hankel transform to equations (27) and (28) using equation (22), we obtain(29)qαf¯Hr1n,q=A0q+A1qq2+ω2J1r1nr1n+−r1n2f¯Hr1n,q+r1nJ1r1nf¯1,q+Pcg¯Hr1n,q−f¯Hr1n,q−M0f¯Hr1n,q+Grθ¯Hr1n,q,(30)g¯Hr1n,q=1Pm1qα+1/Pmf¯Hr1n,q,where(31)∫01∂2f¯r1,q∂r12+1r1∂f¯r1,q∂r1dr1=−r1n2⋅f¯Hr1n,q+r1nJ1r1nf¯1,q.Simplifying equation (29) and using equation (22) leads to(32)q2α+k0qα+k1qα+k2f¯Hr1n,q=F¯0nq+J1r1nq+Grθ¯r1n,q,where(33)F¯0nq=A0q+A1qq2+ω2J1r1nr1n,k1=r1n2+M0Pm,k0=Pmr1n2+1+PmPc+M0Pm,k2=1Pm.Substituting equations (25) into (32), then, after simplification, it takes the form(34)f¯Hr1n,q=r1nJ1r1nq+F¯0nqJ1r1nr1n+Gr1q−qα−1qα+a1J1r1nr1nqα+k2qα+k3qα+k4.Equation (34) can be written as(35)f¯Hr1n,q=J1r1nr1n1q−1q+H0q−1qα+k3+H1q−1qα+k4−H2A0q−1qα+k3+H3A0q−1qα+k4+H2A11qα+k3qq2+ω2+H3A11qα+k4qq2+ω2−GrH9q−1qα+k3+H10q−1qα+k4−H8q−1qα+a1J1r1nr1n,where(36)k3=k0+k0−4k12,k4=k0−k0−4k12,k5=k3+k4,k6=k3k4,k7=r1n2−k5,k8=r1n2k2−k6,H0=k8−k7k3−k32k3−k4,H1=−k8+k7k4+k42k3−k4,H2=k1−k2k3−k4,H3=k2−k4k3−k4,H4=a1a1−k3,H5=k3a1−k3,H6=a1a1−k4,H7=k4a1−k4,H8=H4+H6,H9=H2+H5,H10=H3+H7.Taking the Laplace inverse of equation (35), we obtain(37)fHr1n,τ=J1r1nr1n−J1r1nr1n1+N1Rα,−1−k3,τ+N2Rα,−1−k4,τ+N3cosωt∗Fα−k3,τ+N4cosωt∗Fα−k4,τ+N5Rα,−1−a1,τ,where Robotnov Hartley’s function [58] is defined as(38)L−11sβ0+λ=Fβ0−λ,τ=∑n=0∞−λnτn+1β0−1Γn+1β0,N1=H0−H2A0−GrH9,N2=H1−H3A0−GrH10,N3=−H2A1,N4=−H3A1,N5=GrH8.Applying inverse Hankel transform to equation (37) leads to(39)fr1,τ=1−2∑n=1∞J0r1r1nr1nJ1r1n1+N1Rα,−1−k3,τ+N2Rα,−1−k4,τ+N3cosωt∗Fα−k3,τ+N4cosωt∗Fα−k4,τ+N5Rα,−1−a1,τ.Taking inverse Laplace transform of equation (30) and by using Lorenzo–Hartley function,(40)gHr1n,τ=1PmFα−1Pm,τ∗fHr1n,τ.Taking inverse Finite Hankel transform of equation (40) reduces to(41)gr1,τ=1PmFα−1Pm,τ∗fr1,τ,


### 3.3. Limiting Cases

Case 1: takeα=1.Whenα⟶1 the Lorenzo and Hartley’s and Robotnov and Hartley’s function reduces to [58](42)R1,−1−λ,τ=L−1s−1s+λ=∑n=0∞−λnτn+1k!=1−eλτλ,R1,0=e−λτ,L−11s+λ=e−λτ=F1−λ,τ=∑n=0∞−λnτnk!.Equations (26), (39), and (41) reduce to(43)θr1,τ=1−2∑n=1∞J0r1r1nr1nJ1r1ne−α1τ,fr1,τ=1−2∑n=1∞J0r1r1nr1nJ1r1n1+N011−ek3τ+N021−ek4τ+N03e−k3τ+cosωt+w0sinωt+N04e−k4τ+cosωt+w1sinω+N051−ea1τ,gr1,τ=1Pme−1/Pmτ∗fr1,τ.Case 2: whenβ⟶0 in present general solutions, the problem reduces to the result obtained by Ali et al. [59], when 1/β⟶0:(44)N01=N1k3,N02=N2k4,N03=N3k3k32k32−w2,w0=wk32,N04=N4k4k42k42−w2,w1=wk42,N05=N5a1.


## 3.1. Temperature Field Calculation

Applying Laplace transform to equations (14) and (18), obtained equations are(21)qαθ¯r1,q=1Prd2θ¯r1,qdr12+1r1dθ¯r1,qdr1,(22)θ¯1,q=1q,where θ¯r1,q is Laplace transform of θr1,t, where q denotes transformation variable.By the using Finite Hankel transform to equation (21) and putting values from equation (22), we obtain(23)θ¯Hr1n,q=r1nJr1nPr1qqα+a1,where r1n2/Pr=a1.θ¯Hr1n,q shows Hankel transform of θ¯r,q. Equation (23) reduces to(24)θ¯Hr1n,q=J1r1nr1n1q−J1r1nr1nq−1−αqα+a1.Applying inverse Laplace transform by using Lorenzo and Hartley’sRβ0,υ−λ,τ functions [57] to equation (24),(25)θHr1n,τ=J1r1nr1n−J1r1nr1nRα,1−ατ,−a1,L−1s−υsβ0+λ=Rβ0,υ−λ,τ=∑n=0∞−λnτn+1β0−1−υΓn+1β0−υ.Applying the inverse Henkel transform to equation (25), we obtain(26)θr1,τ=1−2∑n=1∞J0r1r1nr1nJ1r1nRα,1−α−a1,τ.


## 3.2. The Calculation for Blood Flow

Applying Laplace transform to equations (16) and (17) and using equation (13), we obtain(27)qαf¯r1,q=A0q+A1qq2+ω2+∂2f¯r1,q∂r12+1r1∂f¯r1,q∂r1+Pcg¯r1,q−f¯r1,q−M0f¯r1,q+Grθ¯r1,q,(28)g¯r1,q=f¯r1,qPmqα+1.Applying the Hankel transform to equations (27) and (28) using equation (22), we obtain(29)qαf¯Hr1n,q=A0q+A1qq2+ω2J1r1nr1n+−r1n2f¯Hr1n,q+r1nJ1r1nf¯1,q+Pcg¯Hr1n,q−f¯Hr1n,q−M0f¯Hr1n,q+Grθ¯Hr1n,q,(30)g¯Hr1n,q=1Pm1qα+1/Pmf¯Hr1n,q,where(31)∫01∂2f¯r1,q∂r12+1r1∂f¯r1,q∂r1dr1=−r1n2⋅f¯Hr1n,q+r1nJ1r1nf¯1,q.Simplifying equation (29) and using equation (22) leads to(32)q2α+k0qα+k1qα+k2f¯Hr1n,q=F¯0nq+J1r1nq+Grθ¯r1n,q,where(33)F¯0nq=A0q+A1qq2+ω2J1r1nr1n,k1=r1n2+M0Pm,k0=Pmr1n2+1+PmPc+M0Pm,k2=1Pm.Substituting equations (25) into (32), then, after simplification, it takes the form(34)f¯Hr1n,q=r1nJ1r1nq+F¯0nqJ1r1nr1n+Gr1q−qα−1qα+a1J1r1nr1nqα+k2qα+k3qα+k4.Equation (34) can be written as(35)f¯Hr1n,q=J1r1nr1n1q−1q+H0q−1qα+k3+H1q−1qα+k4−H2A0q−1qα+k3+H3A0q−1qα+k4+H2A11qα+k3qq2+ω2+H3A11qα+k4qq2+ω2−GrH9q−1qα+k3+H10q−1qα+k4−H8q−1qα+a1J1r1nr1n,where(36)k3=k0+k0−4k12,k4=k0−k0−4k12,k5=k3+k4,k6=k3k4,k7=r1n2−k5,k8=r1n2k2−k6,H0=k8−k7k3−k32k3−k4,H1=−k8+k7k4+k42k3−k4,H2=k1−k2k3−k4,H3=k2−k4k3−k4,H4=a1a1−k3,H5=k3a1−k3,H6=a1a1−k4,H7=k4a1−k4,H8=H4+H6,H9=H2+H5,H10=H3+H7.Taking the Laplace inverse of equation (35), we obtain(37)fHr1n,τ=J1r1nr1n−J1r1nr1n1+N1Rα,−1−k3,τ+N2Rα,−1−k4,τ+N3cosωt∗Fα−k3,τ+N4cosωt∗Fα−k4,τ+N5Rα,−1−a1,τ,where Robotnov Hartley’s function [58] is defined as(38)L−11sβ0+λ=Fβ0−λ,τ=∑n=0∞−λnτn+1β0−1Γn+1β0,N1=H0−H2A0−GrH9,N2=H1−H3A0−GrH10,N3=−H2A1,N4=−H3A1,N5=GrH8.Applying inverse Hankel transform to equation (37) leads to(39)fr1,τ=1−2∑n=1∞J0r1r1nr1nJ1r1n1+N1Rα,−1−k3,τ+N2Rα,−1−k4,τ+N3cosωt∗Fα−k3,τ+N4cosωt∗Fα−k4,τ+N5Rα,−1−a1,τ.Taking inverse Laplace transform of equation (30) and by using Lorenzo–Hartley function,(40)gHr1n,τ=1PmFα−1Pm,τ∗fHr1n,τ.Taking inverse Finite Hankel transform of equation (40) reduces to(41)gr1,τ=1PmFα−1Pm,τ∗fr1,τ,


## 3.3. Limiting Cases

Case 1: takeα=1.Whenα⟶1 the Lorenzo and Hartley’s and Robotnov and Hartley’s function reduces to [58](42)R1,−1−λ,τ=L−1s−1s+λ=∑n=0∞−λnτn+1k!=1−eλτλ,R1,0=e−λτ,L−11s+λ=e−λτ=F1−λ,τ=∑n=0∞−λnτnk!.Equations (26), (39), and (41) reduce to(43)θr1,τ=1−2∑n=1∞J0r1r1nr1nJ1r1ne−α1τ,fr1,τ=1−2∑n=1∞J0r1r1nr1nJ1r1n1+N011−ek3τ+N021−ek4τ+N03e−k3τ+cosωt+w0sinωt+N04e−k4τ+cosωt+w1sinω+N051−ea1τ,gr1,τ=1Pme−1/Pmτ∗fr1,τ.Case 2: whenβ⟶0 in present general solutions, the problem reduces to the result obtained by Ali et al. [59], when 1/β⟶0:(44)N01=N1k3,N02=N2k4,N03=N3k3k32k32−w2,w0=wk32,N04=N4k4k42k42−w2,w1=wk42,N05=N5a1.


## 4. Graphical Results and Discussion

This section briefly highlights the influence of different parametersα, β, M,Pc,andPm on fluid temperature, velocities of blood, and magnetic particles by plotting various graphs.Figures2(a) and 2(b) show the impact of the fractional parameter α on temperature distribution. With the variation in α, different graphs for temperature distribution are developed. These graphs are drawn for fractional parameters keeping all the other physical parameters and the time parameter constant. This variation is termed as a memory effect that cannot be considered in the integer ordered models. The obtained curves will help the experimentalist to best fit the curve with the curve drawn for the experimental data. The behavior of blood’s velocity and magnetic particles has been observed in Figures 3 and 4 for different values of the fractional parameter. It is clear from these figures that the velocity shows opposite behavior for a long and a short time, which describes the effects of memory consideration by the fractional model as discussed in the description of Figures 2(a) and 2(b). These layers can help the researchers to investigate the fluid and the particles’ behavior at any point. Hence, fractional derivatives provide the best way to describe the velocity of blood at a fixed time to validate the experimental or real-time flow behavior.Figure 2
Temperature graph for different values ofα for long and short time when Pr=22.64.
(a)(b)Figure 3
Variation in blood and magnetic particles profile againstr1 due to α for long time: β=0.5,M=0.5,A0=0.5,A1=0.6,Pm=0.8,Pc=0.5,ωt=5π/8,Pr=22.64,and Gr=3.2×102.
(a)(b)Figure 4
Variation in blood and magnetic particles profile againstr1 due to α for short time: τ<1.β=0.5,M=0.5,A0=0.5,A1=0.6,Pm=0.8,Pc=0.5,ωt=5π/8,Pr=22.64, and Gr=3.2×102.
(a)(b)Figure5 portrays the influence of the Brinkman-type fluid parameter on blood and particle velocities. From the figure, it is noticed that when the Brinkman-type fluid parameter is increased, the velocities of both the fluid and particles decreased. Since the Brinkman parameter is the ratio of drag force and density of fluid [30], the increase in the Brinkman-type fluid parameter increases in friction/drag forces and reduces the blood and particles’ velocities. Figure 6 displays the impact of the magnetic parameter on blood and particle velocity. When the magnetic parameter is increased, velocities of both the fluid and particle show decreasing behavior. This happens due to the Lorentz force, which is opposite to the flow direction. When the values of the magnetic parameter are increased, the Lorentz force increases, which slows down both the velocities. This result indicates that the magnetic field helps control the blood flow during treatment without affecting the healthy tissues. Figures 7 and 8 display the behavior of fluid and particle velocity for various values Pc for a long and short time. Decreasing behavior is observed for a long and short time. The physics behind this is, due to an increase in particle concentration, the blood becomes denser and thicker, which slows down the motion of fluid and particles. Since the flow is fully developed, fluid gets its maximum velocity at the center, while near the wall, its velocity is relatively slow. The same behavior is observed by Saqib et al. [33]. Figure 9 depicts that rising values of magnetic particle mass parameter decreases the velocities of blood and particle. The reason for such behavior is that large size magnetic particles make the blood more viscous and hence retards the flow. Nandkeolyar and Das [60] noticed the similar behavior for fluid velocity by variation in the particle mass parameter. This shows that magnetic particles with small size are more effective for treatments as can easily move through blood and easily reach the target area, including liver and heart cells.Figure 5
Variation in blood and magnetic particles profile againstr1 due to variation in the Brinkman parameter: α=0.5,M=0.5,A0=0.5,A1=0.6,Pm=0.8,Pc=0.5,ωt=5π/8, Pr=22.64,and Gr=3.2×102.
(a)(b)Figure 6
Variation in blood and magnetic particles profile againstr1 due to variation in the magnetic parameter: α=0.5,β=0.5,A0=0.5,A1=0.6,Pm=0.8,Pc=0.5, ωt=5π/8,Pr=22.64,and Gr=3.2×102.
(a)(b)Figure 7
Variation in blood and magnetic particles profile againstr1 due to variation in Pc for long time: α=0.5,β=0.5,A0=0.5,A1=0.6,Pm=0.8,ωt=5π/8,Pr=22.64,and Gr=3.2×102.
(a)(b)Figure 8
Variation in blood and magnetic particles profile againstr1 due to variation in Pc for a short time when α=0.5,β=0.5,A0=0.5,A1=0.6,Pm=0.8, ωt=5π/8,Pr=22.64,and Gr=3.2×102.
(a)(b)Figure 9
Variation in blood and magnetic particles profile againstr1 due to variation in the particle concentration parameter: α=0.5,β=0.5,A0=0.5,A1=0.6,Pc=0.8, ωt=5π/8,Pr=22.64,and Gr=3.2×102.
(a)(b)


## 5. Conclusions

In the present paper, we have considered Brinkman-type fluid in a heated cylindrical tube with MHD effects. Flow equation is modeled by using Caputo time-fractional derivative. Laplace and Hankel’s transforms are employed to get the exact solutions to the problem. Furthermore, the impact of the different parameter on fluid temperature and velocity of blood and particles are discussed in detail. The present investigation plays a significant role in the study of blood flow in presence of magnetic field. Some of the useful findings are listed below:For long time, increase in the fractional parameter, a reduction in the fluid and particle velocity has been observed, while for short time, reverse behavior has been observed.A decrease in velocity of blood and magnetic particles has been observed for larger value of the magnetic number, which implies blood flow can be controlled by using magnetic field.Retardation in velocity distribution of fluid and particles have been noticed when the Brinkman parameter is increased.Increasing the values of concentration parameter causes a decrease in fluid and particles’ velocityA decrease in both the velocities of blood and particles has been observed for rising values of particle mass parameter. It is true because from the law of inertia if the applied forces are kept constant, then the increasing mass of particles slows down the motion of the subject fluid.



---

*Source: 6624912-2021-07-09.xml*