Magnetic Field, Variable Thermal Conductivity, Thermal Radiation, and Viscous Dissipation Effect on Heat and Momentum of Fractional Oldroyd-B Bio Nano-Fluid within a Channel

: This study deals with the analysis of the heat and velocity proﬁle of the fractional-order Oldroyd-B bio-nanoﬂuid within a bounded channel. The study has a wide range of scope in modern ﬁelds of basic science such as medicine, the food industry, electrical appliances, nuclear as well as industrial cooling systems, reducing pollutants, ﬂuids used in the brake systems of vehicles, etc. Oldroyd-B ﬂuid is taken as a bio-nanoﬂuid composed of base ﬂuid (blood) and copper as nanoparticles. Using the fractional-order Oldroyd-B parameter, the governing equation is generalized from an integer to a non-integer form. A strong approach, i.e., a ﬁnite difference scheme, is applied to discretize the model, because the fractional approach can well address the physical phenomena and memory effect of the ﬂow regime. Therefore, a Caputo fractional differentiation operator is used for the purpose. The transformations for the channel ﬂow are utilized to transﬁgure the fractional-order partial differential equations (PDEs) into non-dimension PDEs. The graphical outcomes for non-integer ordered Oldroyd-B bio-nanoﬂuid dynamics and temperature proﬁles are navigated using the numerical technique. These results are obtained under some very important physical conditions applied as a magnetic ﬁeld effect, variable thermal conductivity, permeable medium, and heat source/sink. The results show that the addition of (copper) nanoparticles to (blood) base ﬂuids enhances the thermal conductivity. For a comparative study, the obtained results are compared with the built-in results using the mathematical software MAPLE 2016.


Introduction
Most researchers are currently taking a keen interest in non-Newtonian fluids because of their vast and fruitful applications in emerging and modern technologies, especially in medical science.Some examples of non-Newtonian fluids include polymer solutions, honey, toothpaste, paints, blood, jams, and greases.There are many models to predict the features of these fluids.For non-Newtonian fluids, the relation between shear stress and shear rate is nonlinear.One most important class of rate-type fluids is Oldroyd-B fluid, which is the most appropriate fluid model for investigating the memory and elastic effects [1].It is also a well-known fact that nanofluids are the next exciting frontier of the modern era.Nanofluids are considered superior fluids because of their high thermal conductivity, viscosity, and high heat transfer rate.The most important domain of thermal engineering is heat transfer and a variety of information about how the body transfers heat has been explored successfully by many researchers.
The role of mathematical models is very important in heat transfer problems in medical science such as in cryosurgery, laser surgery, cancer therapy, and in physiotherapy.In cancer therapy, the cancerous tissues are subject to a thermal environment at 42 • C and must maintain surrounding normal tissues at a suitable temperature.A mathematical model of heat transfer for the blood flow is of significant importance [2].The mathematical model was developed to calculate heat loss for locally applied stimuli and to analyze the relation between skin blood flow and temperature [3].He et al. [4] investigated the relationship between blood circulation and peripheral temperature.
Nanofluids have numerous applications in medical science in the reduction of bleeding during surgery, the rapid targeted transport of drugs, tracing of diseased areas, and cancer tumor treatment.Saqib et al. [5] explored the application of magnetite dust particles in magnetic blood flow through a cylindrical tube.For cancer therapy, one of the best-known methods is to develop a magnetic field near the cancer tumor insight to capture copper or gold nanoparticles at the tumor site.Copper or gold particles act as heat sources under the effect of a magnetic field or laser beam [6].In hypothermia cases, the temperature of tumor cells has a key role.The cancerous tissues can be distracted when heated from 42 • C to 45 • C. Lin et al. [7] explain in detail how the time for cancer therapy can be reduced to half by improving the temperature to 1 • C. Misra et al. [8] have produced literature about damage to body tissues and damage in blood plasma when the temperature goes above 42 • C. When nanoparticles are injected into blood veins the are termed bio-magnetic nanofluid.By using nanofluid, the damaging effect on blood plasma can be minimized by improving the heat transfer rate.
To overcome the problem of entropy and heat management systems, the idea of the nanofluid was first given in 1995 by Choi et al. [9].Nanofluids consist of nano-sized particles of metals, metal oxides, carbides, silica, and carbon nano tubes ranging from 1 nm to 100 nm in base fluids such as water, alcohol, oil, or blood.Tiwari et al. [10] investigated the shape and size effect of nanoparticles on the flow of nanofluids in a square cavity.Sherm et al. [11] investigated the effects of thermal dispersion in a porous cavity filled with nanofluid by using the Tiwari and Das model.Wong et al. [12] gave a detailed description of the application of nanofluids in the medical field, especially in drug delivery.James et al. [13] elaborated the effects of thermo-physical properties on the thermal transport of bio-nanofluid.Elkhair et al. [14] investigated the effect of magnetic fields and nonlinear thermal radiations on the thermal transport of hybrid bio-nanofluid in detail.Eid et al. [15] discussed the NP-shape effect of gold nano-particles on the thermal transport of bio-nanofluid in the suction/blowing process.
The focus of modern researchers is on fluids, which are modeled fractionally.Currently, a new version of the order of the PDEs has been introduced for a deep and brief study of different physical problems such as chemical reactions, the fluid flow regime, conduction as well as convection, etc.The fractional (non-integer) order of PDEs gives the best illustration of physical phenomena, can control it, and provides the memory effect.The role of mathematical models in describing the different physical phenomena of daily life is more important.These models are composed of differential equations and some initial and boundary conditions.These differential equations are of different orders and may be integer orders or non-integer fractional derivatives.But non-integer fractional derivatives have a vibrant scope in the modern industrial era.Fractional models are close to the geometry of physical systems.In addition, the fractional parameter is part of the solution which elaborated the memory effect more briefly, helping to control the flow regime, consequently, fractional models are considered superior to ordinary models.They were used to establish the models of the physical systems and many engineering processes by fractional differential equations.Such models were found more accurate describe with the help of fractional differential equations.Many physical phenomena of thermal transport, engineering, circuits, biotechnology, and the signal process can be modified with the help of fractional derivatives.Viscoelastic damping in polymers can be well modeled by using a fractional derivative.The shape memory effect during deformation can be accurately described by using fractional calculus.Tanveer et al. [16] use a fractional approach for the thermal analysis of non-Newtonian bio-nanofluid with an application of carbon nanotubes between two parallel plates.Imran et al. [17] used the Prabhakar fractional approach to solve the Casson nanofluid model.Basit et al. [18] solved the Maxwell nanofluid by using a fractional approach with the application of different nanoparticles.The application of the fractional approach for different nanofluid models can be seen in [19][20][21].Imran et al. [22] explored the fractional derivative operator and its applications in bio-heat transfer problems.Anwar et al. [23] investigated the effect of nonlinear thermal radiation, Newtonian heating, and slip conditions on the thermal transport of fractional Oldroyd-B fluids.Wang et al. [24] investigated the comparative study of heat and mass transfer for Oldroyd-B bio-nanofluids in a porous medium; and that viscous stresses destroying the fluctuating velocity gradient is called viscous dissipation.The viscous dissipation effect is an irreversible process, which heated the fluid because the viscosity of fluid takes energy from the motion of fluid and transforms it into internal energy of the fluid.
After a keen study of the literature we detected research gaps, such as the lack of a tensor-based mathematical modeling of fractional-order Oldroyd-B bio-nanofluids within a bounded channel.There exists a deficiency in the elaboration of physical justification behind results obtained by the addition of nanoparticles of copper in base fluid blood.The viscous dissipation effect of the temperature profile has been less discussed.Additionally, numerical results were not obtained by using the strong technique of the finite difference method.The influence of magnetic forces on flow patterns and numerical predictions of the local Nusselt number and skin friction coefficient was not discussed in detail.Nonlinear thermal radiations also affect the temperature profile.The effect of nonlinear thermal radiations and variable thermal conductivity is not taken into account or not discussed in detail in the discussed literature.
To deliberate these research gaps, a mathematical model based on fractional partial differential equations was obtained using the Caputo derivative between two prepared parallel plates.The results were obtained by applying a strong magnetic field and adding a heat source.The variable thermal conduction was considered for the problem in a permeable medium.Blood was taken as a base fluid and copper particles as a nanoparticle.The governing equations and boundary conditions were transformed into dimensionless forms by using suitable transformations.Subsequently, using the strong technique of finite different method a discretization was made.The mathematical software MAPLE 2016 was used for developing and executing the code to obtain the graphical results.

Mathematical Formulation
Consider the unsteady natural convection flow of an incompressible fractional Oldroyd-B bio nano-fluid within a channel.A strong magnetic field B 0 is applied across the channel and ignores the induced magnetic field.One of the plates is fixed along the x − direction.Both plates are separated by the distance d between them.A magnetic field is applied along the y − direction that is normal to plates.At the time t = 0 the plates as well as bio-nanofluids are assumed to have the temperature θ d .Over time, the temperature rises to θ 0 .The geometry of the deliberated problem is represented in Figure 1.
The proposed model contains the following assumptions.

•
Neglecting the pressure gradient Applying magnetic field (neglecting induced magnetic field) Viscous dissipation is significant.
The equation of continuity for the foregoing assumptions may be seen in [25].
The stress tensor for the Oldroyd-B fluid is given in [26].The proposed model contains the following assumptions.

•
Neglecting the pressure gradient • Thermal conductivity is considered a variable

•
Applying magnetic field (neglecting induced magnetic field)

•
Viscous dissipation is significant.
The equation of continuity for the foregoing assumptions may be seen in [25].
Here T is the Cauchy stress tensor,  is the dynamic pressure, I is the identity tensor, and S is the extra stress tensor.The relation for the extra stress tensor is given in [26].
where  1 ,   , ,  1 are the time relaxation parameter, time retardation parameter, coefficient of dynamic viscosity, and first Reivlin-Ericksen tensor respectively.
The relation for the first Reivlin-Ericksen tensor is given in [26].
= (∇) + (∇)  . ( The Navier Stokes equation for fractional Oldroyd-B bio-nanofluids flowing through porous media is given in [27].Here T is the Cauchy stress tensor, p is the dynamic pressure, I is the identity tensor, and S is the extra stress tensor.The relation for the extra stress tensor is given in [26].
where λ 1 , λ r , µ, A 1 are the time relaxation parameter, time retardation parameter, coefficient of dynamic viscosity, and first Reivlin-Ericksen tensor respectively.The relation for the first Reivlin-Ericksen tensor is given in [26].
The Navier Stokes equation for fractional Oldroyd-B bio-nanofluids flowing through porous media is given in [27].
where ρ n f , V, T, J, B, G, β n f , θ, θ 0 , R x are the viscosity of the bio-nanofluid, the velocity of fluid, Cauchy stress tensor, current density, applied magnetic field, gravitational acceleration, thermal expansion coefficient, the temperature of bio-nanofluid, ambient temperature and Darcy's resistance, respectively.In Equation ( 5) Because the fluid is magneto-hydro-dynamics that is subjected to a strong magnetic field taking into account that B = B 0 + b 0 (induced magnetic field which is ignored) given in [27].
The constitutive relation for the Oldroyd-B fluid is given in [28], Relation for Darcy resistance for Oldroyd-B fluid is given in [29], ( Using Equations ( 10) and (11) in Equation ( 9), The energy equation for this deliberated problem is given in [30], where ρ, C p , Q, and K * n f (θ) represent the density, specific heat capacity, heat source/sink, and variable thermal conductivity.
The relation of the variable thermal conductivity is followed from [31], And by using the Rosland approximation for the thick bio-nanofluid the heat flux q r may be calculated in [32] as, Using Equations ( 14) and (15) in Equation (13), We applied both the sides to Equation ( 16) The appropriate initial and boundary conditions are The following transformations and dimensionless parameters are used for channel flow Thermo-physical properties of the nanofluid are given in [33], Using transformations and relations given in Equation (20) in Equation ( 12), ( 17)-( 19), the momentum and energy equations are given respectively as, Initial and boundary conditions in the dimensionless form are given as, Following are dimensionless parameters that appeared in Equations ( 21) and (22), where, M, Gr, E c , Pr, Nr, Q s and K 1 represent the magnetic field parameter, thermal Grashof number, Eckert number, Prandtl number, the heat generation/absorption parameter, thermal radiation parameter, and the porosity parameter.Equations ( 23) and ( 24) are modeled by using the Caputo derivative operator in appropriate initial conditions and boundary conditions.Coupled nonlinear equations are solved by a powerful numerical scheme which is known as finite difference method.The code is developed and executed successfully and graphs are plotted by using MAPLE 2016.Thermo physical values of nanoparticles (copper) and the base fluid (blood) are given in Table 1.
Table 1.Thermo-physical properties of blood and copper nanoparticles as given in [24].

Skin Friction and Nusselt Number
In an ordinary integer order system, the local skin friction and Nusselt number are defined in [34], By using Equation ( 10), the skin friction coefficient and local Nusselt number for fractional Oldroyd-B bio-nano fluid can be written as(detail can be seen in [35]), Non-dimensional forms of Equations ( 26) and ( 27) are

Numerical Procedure
Finite-difference schemed to handle the governing set of fractional-order fluid problems and heat transfers are given in this section.Discretization of the derivative of fractional order of u, u t and u yy are given in [36] as In the above expressions, d α l = −l 1−α + (1 + l) 1−α for l = 1, 2, 3, . . ., j.The rectilinear grid is subsequently considered to examine the solution of the governing set of fractionalorder fluid problems and heat transfer having grid spacing ∆y > 0, ∆t > 0 in the direction of space and time separately, where ∆t = T N , ∆x = L M for ∆y, ∆t from Z + .The inner points y i , t j in the discussed domain Ω = [0, T] × [0, L] are given as i∆y = y i and j∆t = t j .The discretization of the governing set of fractional-order fluid problems and heat transfer at y i , t j is given as Fractal Fract.2022, 6, 712 9 of 20 for j = 1, 2, 3, . . ., N − 1, i = 1, 2, 3, . . ., N − 1, with the following initial and boundary conditions, Test Problem.In order to check the efficiency of the proposed algorithm, suppose the following problem as Here, the initial and boundary conditions are specified below and an inhomogeneous term can be designated against the choice of fractional-order derivatives.
The exact solution of this considered problem is U(x, t) = x(x − t)t 2 .Several simulations have been accomplished to check the accuracy of the suggested method.Figure 2a,b is plotted between various values of N and the maximum absolute error (MAE), and computational order of convergence (COC) given below:

Numerical Analysis and Discussion
In this manuscript, the velocity profile (, ) and temperature profile (, ) of the fractional Oldroyd It is well known that the method is convergent beside each fractional-order derivative and its convergence order increases as α → 1 .Figure 2c,d comprises the L ∞ -norm among the consecutive solutions i.e., U j+1 − U j ∞ and |U i+1 − U i | ∞ when 0 ≤ i, j ≤ N, M = 500.It has been established that the suggested method is very effective, reliable, and accurate for the considered problem.It also validates that the solution is unchanging against the fractional order mesh parameters.

Numerical Analysis and Discussion
In this manuscript, the velocity profile u(y, t) and temperature profile θ(y, t) of the fractional Oldroyd-B bio-nanofluid are plotted against the different physical parameters Ha the square of the Hartmann number (magnetic field parameter), Nr thermal radiation parameter, Pr Prandtl number, Gr thermal Grashof number, φ volumetric fraction of nanoparticles (of Copper), α, β fractional-order (non-integer order) parameters and Q s heat absorption/source parameter in detail.The skin friction C f and local Nusselt number Nu are calculated against different physical parameters.The partial differential equations of fractional order (using the Caputo fractional order derivative) are discretized by the finite difference method.Finally, the simulation is formed by using mathematical software MAPLE 2016.
The fractional parameter plays an important role to control the thermal and momentum boundary layers.The fluid initially is at rest but over time there is an increase in the value of the velocity profile against increasing values of time; the relaxation parameter λ α 1 can be seen in Figure 3. in the value of the velocity profile against increasing values of time; the relaxation parameter  1  can be seen in Figure 3.The fluid velocity increases by increasing the value of (Gr) as thermal buoyancy forces dominated the viscous force.Advancement in the fractional parameter α boosted the velocity profile because the momentum boundary layer was reduced by increasing the value of the fractional parameter α. Results for Gr against the fractional parameter α are portrayed in Figure 4.It can be noticed that the velocity profile increased for increasing values of α; in addition, the same behavior of the velocity profile is seen for increasing values of Gr.
The fluid velocity increases by increasing the value of () as thermal buoyancy forces dominated the viscous force.Advancement in the fractional parameter α boosted the velocity profile because the momentum boundary layer was reduced by increasing the value of the fractional parameter α. Results for Gr against the fractional parameter α are portrayed in Figure 4.It can be noticed that the velocity profile increased for increasing values of ; in addition, the same behavior of the velocity profile is seen for increasing values of .With an enhancement in the strength of the magnetic field parameter (Ha), there is a decrease in fluid velocity.It can be justified as, by enhancing the magnetic field strength, the Lorentz force increased, which shows high resistance to fluid flow.Enhancing the With an enhancement in the strength of the magnetic field parameter (Ha), there is a decrease in fluid velocity.It can be justified as, by enhancing the magnetic field strength, the Lorentz force increased, which shows high resistance to fluid flow.Enhancing the magnetic field parameter amounts to a stronger Lorentz force and hence the rate of heat transfer in the boundary layer is higher.By enhancing the value of the fractional parameter α the momentum boundary layer decreased and became thinnest at α = 1 , consequently, the velocity increased for increasing values of the fractional parameter α.In Figure 5, the velocity profile is portrayed for Ha = 0 and for α = 0.4, 0.7, 1, and also for Ha = 2 against the fractional parameter α = 0.4, 0.7, 1 and Ha = 5 against α = 0.4, 0.7, 1 .In Figure 5, it can be seen that the fractional parameter enhanced results in enhancing the fluid velocity profile.But for increasing values of the Hartmann number the fluid velocity decreased.
Fractal Fract.2022, 6, 712 12 of 20 magnetic field parameter amounts to a stronger Lorentz force and hence the rate of heat transfer in the boundary layer is higher.By enhancing the value of the fractional parameter α the momentum boundary layer decreased and became thinnest at  = 1 , consequently, the velocity increased for increasing values of the fractional parameter α.In Figure 5, the velocity profile is portrayed for Ha = 0 and for  = 0.4, 0.7, 1, and also for Ha = 2 against the fractional parameter  = 0.4, 0.7, 1 and Ha = 5 against  = 0.4, 0.7, 1 .
In Figure 5, it can be seen that the fractional parameter enhanced results in enhancing the fluid velocity profile.But for increasing values of the Hartmann number the fluid velocity decreased.The effect of the volume fraction parameter () of nanoparticles on the fluid velocity profile is depicted in Figure 6.Advancement in control volume fraction parameter () reduced the fluid velocity because by increasing () the fluid becomes more viscous.Because the addition of nanoparticles increases the viscosity, consequently, the velocity profile decreases.But the opposite behavior of fluid velocity is observed for the fractional parameter .Enlarging the values of α boosted the velocity profile as can be seen in Figure 6.The effect of the volume fraction parameter of nanoparticles on the fluid velocity is reported in [37], and has a good agreement with present results.The effect of the volume fraction parameter (φ) of nanoparticles on the fluid velocity profile is depicted in Figure 6.Advancement in control volume fraction parameter (φ) reduced the fluid velocity because by increasing (φ) the fluid becomes more viscous.Because the addition of nanoparticles increases the viscosity, consequently, the velocity profile decreases.But the opposite behavior of fluid velocity is observed for the fractional parameter α.Enlarging the values of α boosted the velocity profile as can be seen in Figure 6.
The effect of the volume fraction parameter of nanoparticles on the fluid velocity is reported in [37], and has a good agreement with present results.
The effect of the volume fraction parameter () of nanoparticles on the fluid v ity profile is depicted in Figure 6.Advancement in control volume fraction param () reduced the fluid velocity because by increasing () the fluid becomes more cous.Because the addition of nanoparticles increases the viscosity, consequently velocity profile decreases.But the opposite behavior of fluid velocity is observed fo fractional parameter .Enlarging the values of α boosted the velocity profile as ca seen in Figure 6.The effect of the volume fraction parameter of nanoparticles on the id velocity is reported in [37], and has a good agreement with present results.The plots of Figure 7 express the results for the velocity profile u(u, t) of the fractional Oldroyd-B bio-nanofluid against the porosity parameter K 1 of the viscoelastic bio-nanofluid.The overall flow regime within a porous medium is well described by Darcy's law which states that the flow rate of fluids is proportional to the potential energy gradient, because it is the property of porous media that it contains the pores or so-called voids that gives resistance to flowing particles.These expected results are shown in the following Figure 7 an increasing values of porosity parameter K 1 decreases the velocity of the fractional Oldroyd-B bio-nanofluid against increasing values of fractional parameter α.
The plots of Figure 7 express the results for the velocity profile (, )of the f tional Oldroyd-B bio-nanofluid against the porosity parameter  1 of the viscoela bio-nanofluid.The overall flow regime within a porous medium is well described Darcy's law which states that the flow rate of fluids is proportional to the potential ergy gradient, because it is the property of porous media that it contains the pore so-called voids that gives resistance to flowing particles.These expected results shown in the following Figure 7 an increasing values of porosity parameter  1 decre the velocity of the fractional Oldroyd-B bio-nanofluid against increasing values of f tional parameter α.For the Oldroyd-B bio-nanofluid the velocity profile decreased with increasing ues of time; the retardation parameter   =  3 can be seen in Figure 8.For the Oldroyd-B bio-nanofluid the velocity profile decreased with increasing values of time; the retardation parameter λ r = λ 3 can be seen in Figure 8.For the Oldroyd-B bio-nanofluid the velocity profile decreased with increasing values of time; the retardation parameter   =  3 can be seen in Figure 8.By increasing the volume fraction parameter of nanoparticle () the heat transfer becomes better.Because nanoparticles improved, the heat transfer rate of the base fluid resulted in the temperature of the system decreasing because more heat transferred from the system and cooled the system more quickly, as can be seen in Figure 9.But the op- By increasing the volume fraction parameter of nanoparticle (φ) the heat transfer becomes better.Because nanoparticles improved, the heat transfer rate of the base fluid resulted in the temperature of the system decreasing because more heat transferred from the system and cooled the system more quickly, as can be seen in Figure 9.But the opposite behavior of fractional parameter α is seen in Figure 9.The thermal boundary layer becomes thin for increasing values of α and becomes thinnest for α = 1.As values of fractional parameter α increased the fluid temperature also increased, as can be observed in Figure 9.In the physical phenomena involving heat transfer, energy loss is the major issue.The idea of the addition of nanoparticles to the base fluid is to control the enthalpy of the system.These expected results are obtained, as depicted in the following Figure 9.The reason behind this is that Copper nanoparticles are good conductors of heat.The addition of Copper nanoparticles in flowing fluids (blood) reduced the energy loss that is the decreased temperature profile.
Fractal Fract.2022, 6, 712 14 of 20 posite behavior of fractional parameter α is seen in Figure 9.The thermal boundary layer becomes thin for increasing values of α and becomes thinnest for  = 1.As values of fractional parameter α increased the fluid temperature also increased, as can be observed in Figure 9.In the physical phenomena involving heat transfer, energy loss is the major issue.The idea of the addition of nanoparticles to the base fluid is to control the enthalpy of the system.These expected results are obtained, as depicted in the following Figure 9.
The reason behind this is that Copper nanoparticles are good conductors of heat.The addition of Copper nanoparticles in flowing fluids (blood) reduced the energy loss that is the decreased temperature profile.The relative effect of the kinetic energy of flow to the enthalpy difference across the thermal boundary layer is measured by an Eckert number (Ec).An increase in the value of (Ec) means increasing the viscous dissipation and consequently increasing the fluid temperature can be seen in Figure 10.Physically it can be justified that viscous dissipation is an irreversible conversion of kinetic energy into internal heating through an effort that opposes the viscous strain.It is a process where energy is utilized or lost without performing meaningful work.Similar behavior for fractional parameter  can be observed in Figure 10.Similar results are reported in [38].The relative effect of the kinetic energy of flow to the enthalpy difference across the thermal boundary layer is measured by an Eckert number (Ec).An increase in the value of (Ec) means increasing the viscous dissipation and consequently increasing the fluid temperature can be seen in Figure 10.Physically it can be justified that viscous dissipation is an irreversible conversion of kinetic energy into internal heating through an effort that opposes the viscous strain.It is a process where energy is utilized or lost without performing meaningful work.Similar behavior for fractional parameter α can be observed in Figure 10.Similar results are reported in [38].
The relative effect of the kinetic energy of flow to the enthalpy difference across thermal boundary layer is measured by an Eckert number (Ec).An increase in the v of (Ec) means increasing the viscous dissipation and consequently increasing the f temperature can be seen in Figure 10.Physically it can be justified that viscous diss tion is an irreversible conversion of kinetic energy into internal heating through an e that opposes the viscous strain.It is a process where energy is utilized or lost with performing meaningful work.Similar behavior for fractional parameter  can be served in Figure 10.Similar results are reported in [38].In the following Figure 12 graphical results are drawn for temperature profile θ(y, t) against physical parameter Q s .It is obvious that the high value of Q s , the heat source/sink term, increases the temperature, because the heat source rises to the overall heat transfer rate of the base fluid (blood).In the following Figure 12 graphical results are drawn for tempe (, ) against physical parameter   .It is obvious that the high value of source/sink term, increases the temperature, because the heat source rises heat transfer rate of the base fluid (blood).Figures 14 and 15 are plotted to compare the results between the proposed scheme and the Maple built-in command of the dimensionless velocity and temperature for various values of Nr and Qs number for α = β = 1.Results obtained by the proposed scheme and built-in commands are in good agreement as followed from the figure.It can be deduced that the numerical solutions obtained via the proposed method are in excellent agreement with the Maple built-in command results.This shows the correctness and effectiveness of the planned method.Figures 14 and 15 are plotted to compare the results between the proposed scheme and the Maple built-in command of the dimensionless velocity and temperature for various values of Nr and Qs number for  =  = 1.Results obtained by the proposed scheme and built-in commands are in good agreement as followed from the figure.It can be deduced that the numerical solutions obtained via the proposed method are in excellent agreement with the Maple built-in command results.This shows the correctness and effectiveness of the planned method.Tables 2 and 3 shows the numerical results for variations in the skin friction coefficient and local Nusselt number.Tables 2 and 3 shows the numerical results for variations in the skin friction coefficient and local Nusselt number.

Conclusions
In this article, a numerical procedure of the finite difference method was adopted.An investigation was performed for the fractional Oldroyd-B bio-nanofluid within a permeable medium.The bio-nanofluid was subjected to a strong magnetic field (neglecting the induced magnetic field).The Caputo fractional derivative was used for fractional modeling of partial differential equations for u(y, t), θ(y, t), co − efficient of skin friction C f and local Nusselt number Nu.The finite difference approach was applied to discretize the dimensional governing equations.Taking blood as the base fluid and gold as nanoparticles, a graphical illustration was made using the mathematical software MAPLE 2016.The key findings of this study are as follows: a The addition of a nanoparticle (copper) to the base fluid (blood) enhances the thermal conductivity of fractional Maxwell bio-nanofluid; this causes the temperature profile θ(y, t) to decrease by the increase in volumetric fractional nanoparticles (Copper).b The temperature profile θ(y, t) gains high value by increasing the values of the variable thermal conductivity constant .c The thermal radiation parameter Nr increases the heat transfer rate of the aforementioned physical problem.d The permeable medium decreases the velocity profile of the viscoelastic fractional bio-nanofluid.e There is an increase in the temperature profile observed by increasing the values of the viscous dissipation parameter Ec.
The solution obtained via the finite difference approach was excellent and in agreement with the built-in result of the test problem and the existing literature.Therefore, the finite difference scheme is an effective and consistent scheme to study such types of multifaceted geometry.This study of fluid dynamics gives a strong base to the biological study of blood flow within the blood arteries (cylindrical geometry).

Figure 1 .
Figure 1.Geometry of the problem.

FractalFigure 2 .
Figure 2. (a) Order of convergence; (b) Maximum absolute error; (c) L-norm for mesh size j; (d) L-norm for mesh size i.

9 Figure 2 .
Figure 2. (a) Order of convergence; (b) Maximum absolute error; (c) L-norm for mesh size j; (d) L-norm for mesh size i.

Figure 4
Figure 4  elaborated on the effect of the Grashof number  on the fluid velocity profile.The Grashof number is the ratio of thermal buoyancy forces to the viscous forces.The fluid velocity increases by increasing the value of () as thermal buoyancy forces dominated the viscous force.Advancement in the fractional parameter α boosted the velocity profile because the momentum boundary layer was reduced by increasing the value of the fractional parameter α. Results for Gr against the fractional parameter α are portrayed in Figure4.It can be noticed that the velocity profile increased for increasing values of ; in addition, the same behavior of the velocity profile is seen for increasing values of .

Figure 4
Figure 4  elaborated on the effect of the Grashof number Gr on the fluid velocity profile.The Grashof number is the ratio of thermal buoyancy forces to the viscous forces.The fluid velocity increases by increasing the value of (Gr) as thermal buoyancy forces dominated the viscous force.Advancement in the fractional parameter α boosted the velocity profile because the momentum boundary layer was reduced by increasing the value of the fractional parameter α. Results for Gr against the fractional parameter α are portrayed in Figure4.It can be noticed that the velocity profile increased for increasing values of α; in addition, the same behavior of the velocity profile is seen for increasing values of Gr.

Figure 11
Figure 11  depicted the effects of variable thermal conductivity on the fluid temperature profile.The fluid temperature increased for increasing values of variable thermal conductivity parameter ∈.Physically it can be justified because when values of ∈ are increased it means the thermal conductivity of fluid increases and consequently more heat is conducted from the surface; this raises the internal temperature of the fluid.Similar behavior of fractional parameter α can be observed in Figure11.Similar results are demonstrated in[31].

Figure 13 Figure 12 .
Figure 13 depicts the effects of thermal radiation parameter  on the capability of the fractional Oldroyd-B bio-nanofluid.Because su bio-nanofluid to thermal radiation give hype to Brownian motion (a rand particles within the fluid field) this increases kinetic energy of nanoparticl

Figure 13
Figure 13 depicts the effects of thermal radiation parameter Nr on the heat transfer capability of the fractional Oldroyd-B bio-nanofluid.Because subjecting the bio-nanofluid to thermal radiation give hype to Brownian motion (a random motion of particles within the fluid field) this increases kinetic energy of nanoparticles during motion; this enhances the thermal conductivity of the base fluid, and consequently the temperature profile increases for the increasing value of Nr.

Figures 14 Figure 13 .
Figures 14 and 15 are plotted to compare the results between the proposed and the Maple built-in command of the dimensionless velocity and temperature ious values of Nr and Qs number for  =  = 1.Results obtained by the p

Figure 14 .Figure 14 .Figure 15 .
Figure 14.Comparison of dimensionless velocity between Proposed Method and Maple Built-in command against (a)  and (b)   .

Figure 15 .
Figure 15.Comparison of dimensionless temperature between Proposed Method and Maple Built-in command against (a) Nr and (b) Q s .

Table 2 .
Variation in the skin friction coefficient as varying the physical parameters and .

Table 2 .
Variation in the skin friction coefficient as varying the physical parameters and α.

Table 3 .
Variation in the local Nusselt number as varying the physical parameters and α.