Two-Phase Flow of Blood with Magnetic Dusty Particles in Cylindrical Region:A Caputo Fabrizio Fractional Model

2021-12-16 06:37AneesImitazAaminaAaminaFarhadAliIlyasKhanandKottakkaranSooppyNisar
Computers Materials&Continua 2021年3期

Anees Imitaz,Aamina Aamina,Farhad Ali,Ilyas Khan and Kottakkaran Sooppy Nisar

1Department of Mathematics, City University of Science and Information Technology, Peshawar, 25000,Pakistan

2Computational Analysis Research Group, Ton Duc Thang University, Ho Chi Minh City, Vietnam

3Faculty of Mathematics and Statistics, Ton Duc Thang University, Ho Chi Minh City, Vietnam

4Department of Mathematics, College of Science Al-Zulfi, Majmaah University, Al-Majmaah, 11952,Saudi Arabia

5Department of Mathematics, College of Arts and Sciences,Prince Sattam bin Abdulaziz University, Wadi Aldawaser, 11991,Saudi Arabia

Abstract: The present study is focused on the unsteady two-phase flow of blood in a cylindrical region.Blood is taken as a counter-example of Brinkman type fluid containing magnetic (dust) particles.The oscillating pressure gradient has been considered because for blood flow it is necessary to investigate in the form of a diastolic and systolic pressure.The transverse magnetic field has been applied externally to the cylindrical tube to study its impact on both fluids as well as particles.The system of derived governing equations based on Navier Stoke’s,Maxwell and heat equations has been generalized using the well-known Caputo-Fabrizio (C-F) fractional derivative.The considered fractional model has been solved analytically using the joint Laplace and Hankel (L&H) transformations.The effect of various physical parameters such as fractional parameter,Gr, M and γ on blood and magnetic particles has been shown graphically using the Mathcad software.The fluid behaviour is thinner in fractional order as compared to the classical one.

Keywords: Two-phase blood flow; dusty fluid; Brinkman type model; magnetic dusty particles;heat transfer;C-F derivative

1 Introduction

Biomagnetic fluid dynamic (BFD) is a new area in fluid mechanics.It focuses on the usage of the magnetic particles as drug carriers in magnetic drug targeting, cancer tumor treatment and many more [1-3].The Biomagnetic fluid occurs in all living organisms and for its investigation, the BFD model was initially recommended by Haik et al.[4].Fluids that show non-linear relation between shear stress and strain are termed as non-Newtonian fluids e.g., blood.Blood is the only biological electrically conducting fluid and its mobility is influenced by an applied magnetic field.It contains plasma and red blood cells(RBC) that are oxides of iron and have hemoglobin fragments in high concentrations [5,6].Due to the oxygenated state, blood exhibits a magnetic nature [7,8].The non-Newtonian behaviour of blood due to the suspension of red blood cells in plasma and human thoracic aorta is analyzed by Caballero et al.[9].Tripathi et al.[10] have examined the Non-Newtonian blood in a channel and attained analytical solutions for the velocity, volumetric flow rate and wavelength.In the human left ventricle (LV), the significance of the non-Newtonian blood has examined by Doost et al.[11].Kumar et al.[12] evaluated the difference between Newtonian and non-Newtonian blood models and concluded that the non-Newtonian blood has more/less augmented wall shear stress as compared to the Newtonian blood.

Since blood is a biological fluid,biological heating is significant for metabolic heat generation[13].The phenomenon in biological fluids was first discussed by Bernard in 1876.Afterwards,bioheat transfer became a topic applied in the practice of biology in a wide variety of applications such as chemotherapy [14,15],human thermoregulation system [16] and others [17].Sharifi et al.[18] investigated the heat transfer applications in peripheral vascular disease using FHD principle through two inclined permanent magnets in a channel.Jimoh et al.[19] studied third-grade fluid in hematocrit with slip velocity.Dutta et al.[20]have developed an analytical study of heat propagation in biological tissues for constant and variable heat flux at the skin surface with hypothermia treatment.Fu et al.[21] reviewed the heat transfer modelling in thermoregulatory responses in the human body.Kengne et al.[22] discussed the bioheat transfer in the spherical biological tissues.Zhang et al.[23] discussed the heat transfer in LN2 cryoprobe systems and obtained effective results.David et al.[24] used the heat transfer in the warming of simulated blood by the generation of electronic components.Zainol et al.[25] investigated the heat transfer model for the prediction of human skin temperature using the bioheat equation.

The consideration of Two-phase flow is due to the presence of numerous interfaces separating two immiscible phases.The blood flow through a tiny tube at a very low shear is responsible for the twophase flow surrounded by a cell-depleted peripheral layer.Different types of particles have been considered as the second phase in blood flow, but the most recommended and suitable particles are magnetic particles.The magnetic particles in blood have a vital role in numerous medical applications[26,27].In drug delivery, a specific number of magnetic particles are used to transport the maximum number of a drug to the area of its choice.Due to the mentioned applications of magnetic particles,several researchers used the two-phase blood flow along with magnetic particles.Verma et al.[28]described a dual-phase blood flow model in thin pipes with the fundamental core of deferred erythrocytes and cell unrestricted film and found the results for the nonlinear problem numerically.The thermal and mass concentration effect of the multiphase blood model in a stenosed artery has been investigated by Tripathi et al.[29].An analytical approach has been used for the results to comprehend the comportment of blood flow rate, wall shear stress and flow resistivity.Arribas et al.[30] created a reliable two-phase RBC model for the blood vessel and calculated the viscosity, phase dispersals and volume fractions using the depletion theory.They have associated their results with numerical as well as experimental study and found extraordinary conclusions.A two-phase model of blood with mild stenosis magnetic field and thermal effects has been explored by Ponalagusamy et al.[31].They have concluded that thermal and shear stress slow down with increase in the levels of the plasma layer thickness and they are very effective for the diseased arterial treatment.Ali et al.[32] examined the two-phase dusty fluid with heat transfer in a fluctuating plate, and found that by enhancing the number of embedded particles, the dusty fluid velocity increases.

Due to multidimensional features,the non-integer order calculus is attracting the attention of scientists and researchers [33,34].Fractional calculus is an important and fruitful tool for describing many systems including memory effects.In the preceding few decades, fractional calculus is used for many purposes in various fields, such as electrochemistry, transportation of water in ground level, electromagnetism,elasticity, diffusion and in conduction of heat process [35].In 2015, Caputo et al.[36] worked together in the field of non-integer order calculus and presented a new expression for the non-integer order derivative with the non-singular kernel.So, keeping in sight the importance of CF operators, many researchers used the CF operator in their studies such as in physics, biological mathematics, and many more.Ali et al.[37]examined the magnetic flow of Walter’s-B fluid by using the CF non-singular operator.Salah Uddin et al.[38] investigated the CF model of blood flow with Ferro particles in cylindrical coordinates and their results were in agreement with the previously published works.Ali et al.[39,40] studied the fractionalized model of blood flow having magnetic particles in cylindrical coordinates.

There is no attempt found in the literature relevant to Caputo-Fabrizio fractional approach to find the closed-form solution for magnetite particles-based blood flow with thermal concentration.Hence, in the present article, the work of Saqib et al.[41] has been generalized by taking the flow of blood as a Brinkman type fluid with magnetic particles in cylinder.The governing equations for both fluid and particles are modelled and using the Caputo-Fabrizio fractional-order approach, the closed-form solutions have been obtained by using joint Laplace and Hankel transformations.The impact of different embedded parameters on blood and particles velocities have been examined through graphs.

2 Mathematical Modeling

The blood flow is considered in a vertical cylinder having a radius r0as represented in Fig.1.The magnetic particles are equally distributed throughout the blood flow.The cylinder has been considered along the z-axis and r1-axis is chosen perpendicular to it.The direction of the motion of the blood and magnetic particles are along the z-axis.The biological thermal effect has also been considered and the radiation has been neglected.The induced magnetic field due to a very slight Reynold number has been ignored [42].For a time t =0, the system is considered to be at rest with ambient temperature T∞.For t >0 the temperature rises to Tw.The force Femagis described by[43,44]

Figure 1: The geometry of the flow

The unsteady Brinkman-type blood flow in a cylinder is specified by:

the oscillating pressure gradient[45]is

where u(r1,t) is the blood velocity, up(r1,t) is the magnetic particles velocity.The termis the force between the fluid and particle due to relative motion and magnetic particles flow is conducted[46]:

The thermal equation is specified by:

subjected to the following IBCs

By incorporating the Non-dimensional variables

into Eqs.(2)-(5), then ignoring the* notation,we obtain:

For a generalized fractional model,the newly developed CF time-fractional derivative has been used to covert the linear model to the fractional model,therefore Eqs.(8)-(10) reduces to:

3 Solution of the Problem

For the solution of Eqs.(14)-(16)the non-dimensional IBC’s from Eq.(11)and the Laplace and Hankel transformations are utilized.

3.1 Energy Equation Solution

Applying the joint L&H transforms using Eqs.(11)-(14),we get

where

Applying inverse L&H transformations to Eq.(15),and by using Lorenzo and Hartley’s’and Robotnov and Hartley’s’ functions,respectively [43], yields:

3.1.1 Heat Transfer Rate(Nu)

The Non-dimensional Nusselt number is given by

3.2 Solution of the Blood and Particle Velocities

To obtain the solution for the blood velocity and Magnetic particles velocity, the Laplace and Hankel transforms have been applied on Eq.(15) using the corresponding transformed boundary conditions by lettingand we get:

Now for the blood velocity, Eq.(18) has been incorporated into Eq.(12) using the corresponding transformations and boundary condition(1, q)=0, which yields:

The simplified form of Eq.(19) is

After further simplification the Eq.(20) will take the following form as

where

In component form Eq.(21) is expressed as:

where

and

By applying inverse Laplace transform to Eq.(22),by using the Lorenzo and Hartley’s’respectively[46],

we get

where

Applying the finite Hankel transform of order zero to Eq.(23),we get

where

Now for the solution of magnetic particles velocity applying the inverse L&H transformations to Eq.(16),yields:

From Tab.1 it can be seen that by growing the fractional parameter an enhancement in (Nu) occurs for time.

Table 1: Time and α variation on Nusselt number

4 Graphical Results and Discussion

The considered work aims to study the generalized two-phase blood flow of Brinkman type fluid in a cylindrical tube.The analytical solutions have been attained for energy, velocity as well as for the magnetic particles contained in the blood.Various parameters have been discussed physically on velocities of the blood, particles and temperature.Fig.1 shows the physical model of the considered problem.Fig.2 shows the effect of fractional parameter on the temperature profile.It can be concluded from the figure that by using the fractional parameter, we obtained different temperature profiles by keeping the other parameters constant and this effect is called the memory effect, which is impossible in integer-order.In this graph, we obtained dual behaviour of temperature memory for different values of α and the same behaviour has been noticed as reported by Ali et al.[39].

Figure 2: Variation in Temperature for diverse values of α for Pr =22.64 at t =1

Fig.3 shows the effect of the fractional parameter α on velocity profiles.The corresponding results for regular blood and particles velocity are compared with the fractional order in a fixed time and with classical order and the fluid and particles memory has been discussed.

Figure 3: Blood and Particle velocities sketch of α at t= 2,Gr=3.2×102,Pr=22.64

Fig.4 shows the impact of the Magnetic parameter on blood velocity and particle velocity.From the figures, it has been concluded that for the higher values of M, the flow of blood, as well as magnetic particles, retards.It is physically true that by increasing M, the Lorentz forces increase, which produces resistive forces due to which the flow retards.This effect is very useful in various medical fields such as Magnetic drug targeting and separation techniques by the magnetic field for the cure of different types of diseases and to maintain the normal aspects of the human body.

Figure 4: Blood velocity sketch of M at t= 2, Gr=3.2×102,Pr =22.64

Fig.5 shows the impact of the Brinkman type parameter γ1on the velocities of blood and particles.It is observed that the velocity decreases with the increasing values of γ1.This is because γ1is the ratio between resistive forces and density.By increasing γ1,the opposite forces increase, which retards the fluid velocity and is directly related to the blood flow.These results are strongly in agreement with Saqib et al.[41].

Figure 5: Blood velocity graph of γ1 at t= 2, Gr=3.2×102,Pr =22.64

5 Conclusions

The Caputo-Fabrizio time-fractional derivative has been used.The effect of relative parameters has been shown graphically.Closed-form expressions have been obtained by using the Laplace transform and Hankel transform techniques.Based on the graphical study,it has been concluded that the velocity profile decreases in the response of an external applied magnetic field and Brinkman parameter.This phenomenon might play an important role in Magnetic wounds.Furthermore,by increasing the fractional parameter,the fluid memory becomes thicker.

Funding Statement:The author(s) received no specific funding for this study.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.