Aligned Magnetic Field and Radiation Effects on Biomagnetic Fluid over an Unsteady Stretching Sheet with Various Slip Conditions

Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/). 1 Department of Mathematics, Comilla University, Comilla 3506, Bangladesh; thomasanikgomes@gmail.com (A.G.); jahangircu1994@gmail.com (J.A.); limonn@yahoo.com (G.M.) 2 Department of Science and Humanities, Military Institute of Science and Technology, Dhaka 1216, Bangladesh; tahmina2011@yahoo.com 3 Fluid Mechanics & Turbomachinery Laboratory, Department of Mechanical Engineering, University of the Peloponnese, 26334 Patras, Greece 4 Research Group of Fluid Flow Modeling and Simulation, Department of Applied Mathematics, University of Dhaka, Dhaka 1000, Bangladesh; ferdows@du.ac.bd * Correspondence: etzirtzilakis@uop.gr


Introduction
During the last few decades, due to its application in several areas in science and engineering, the study of flow and heat transfer over an unsteady stretching sheet has drawn significant attention to researchers.The study of rotating flow and heat transfer has received fervent interest in modern fluid dynamics research, with applications in geophysics, biomedical engineering, medical science, planetary science, thermal insulations, etc.
Biomagnetic Fluid Dynamics (BFD) is the study of the effects of an applied magnetic field on biological fluid flow [1].The most characteristic biomagnetic fluid is blood.Blood is a suspension of numerous cells such as red blood cells, white blood cells, and platelets in a liquid electrolyte solution called plasma.Plasma contains 7% of principal proteins and 90% of water, along with considerable concentration of ions.Blood as a whole is considered as a non-Newtonian fluid predominantly when the characteristic dimension of the flow is nearby the cell dimension.As far as the stretching sheet flows are concerned, Crane [2] computed an exact similarity solution for the boundary layer flow of a Newtonian fluid towards an elastic sheet.The sheet was stretched with the velocity proportional to the distance from the origin.Barozzi and Dumas [3] numerically studied the convective heat transfer in blood vessels of the circulatory system.They observed that the rheological behavior of blood does not significantly affect the heat transfer rate in small blood vessels.Pennes [4] studied the effects of blood perfusion and metabolic heat generation in living tissues using a simplified bio-heat transfer model.Although this model bears the potential to describe the effect of blood flow on tissue temperature, it has some considerable short comings.This is because uniform perfusion rate was assumed, and the direction of blood flow was not accounted for.Moreover, in his model, only the stream of venous blood as the fluid stream equilibrated with tissue was considered.
In recent years, the study of the magneto hydrodynamic (MHD) flow of blood through the arteries has gained considerable interest because of its important applications in physiology.Theoretical estimates of blood flow in arteries during the therapeutic procedure of electromagnetic hyperthermia used for cancer treatment were reported by Misra et al. [5].A few important discussions were also available in that paper.The effects of electromagnetic radiation/ultrasonic radiation on blood flow were studied by other investigations such as those of Inoue et al. [6], Nishimoto et al. [7], Bidin and Nazar [8], Irfan et al. [9], and Ishak [10].The effect of viscous dissipation and radiation on the unsteady flow of electrically conducting fluid passed over a stretching surface was considered by Brickman [11] and Chand et al. [12].The variable viscosity and thermal conductivity effects of combined heat and mass transfer in mixed convection over a UHF/UMF wedge in porous media were analyzed by Hassanien et al. [13] and Khan et al. [14].Moreover, the effects of variable viscosity and thermal conductivity on a thin film flow over a shrinking/stretching sheet were also studied.Pal and Mondal [15] studied the effects of temperature-dependent viscosity and variable thermal conductivity on an MHD, non-Darcy mixed convection diffusion of species over a stretching sheet.The effects of thermal radiation over a stretching sheet under several flow conditions have been also studied by several researchers [16][17][18][19][20].
Furthermore, a study of microploar fluid under the influence of a magnetic field through a stretched curved surface, using the Cattaneo-Christov heat model investigated by Khan et al. [21], had found that, in fluid velocity, the magnetic field parameter plays a significant role.The movement of peristaltic flow of a dusty fluid with elastic properties in a curved configuration was analyzed by Khan et al. [22].A fully developed model of non-Newtonian fluid through a 2D stretching sheet in the presence of Lorentz force and internal heat was presented by Vijaya et al. [23].The BFD model [1], which involves both ferrohydrodynamics (FHD) and magnetohydrodynamics (MHD) principles, was utilized for the study of the effect of thermal radiation through a two-dimensional unsteady stretching sheet by Alam et al. [24].Finally, the impact of MHD on non-Newtonian mass and heat transfer along a curved stretched sheet was numerically studied by Yasim et al. [25].
The aim of the present investigation is to study the flow and heat transfer in a stretching sheet with an angle α to the vertical plane in the presence of a non-uniform source/sink.The mathematical formulation of the effect of the magnetic field is that of BFD, involving both principles of ferrohydrodynamics (FHD) and MHD [1].The governing partial differential equations have been transformed by similarity transformations into a coupled system of nonlinear ordinary differential equations.The solution was attained by using a MATLAB package.The effects of various parameters on the momentum and heat transfer characteristics have been studied, and the numerical results are presented graphically for the various values of the parameters entering the problem into consideration.

Model Description
The unsteady two-dimensional BFD flow of a viscous incompressible fluid past a stretching sheet with an acute angle α to the vertical is considered.
Where u and v are the velocity components along X-direction and Y-direction, respectively.The X-axis is considered along the plate and Y-axis is taken normal to it (see Figure 1).Initially (t = 0), the sheet is stretched with velocity U w (x, t) = ax (1−ct) along the X-axis, whereas the origin is kept fixed in the fluid medium of ambient temperature T ∞ , and T w is the stretching surface temperature.A magnetic field of uniform strength B(t) is acting normal to the direction of the flow, with an acute angle ξ.The magnetic Reynold's number is assumed as much less than unity, and the flow is considered twodimensional, therefore the induced magnetic field can be neglected in comparison to the applied magnetic field.Moreover, the fluids exhibit polarization due to the incorporated principles of FHD, and the applied magnetic field is assumed to be strong enough to attain equilibrium magnetization.The velocity slip, thermal slip, viscous dissipation parameter, and ferromagnetic interaction parameter have been taken into account.The boundary layer equations of the fluid and energy equation for the problem can be written as [1,26,27].
tion, respectively.The X -axis is considered along the plate and Y -axis is taken normal to it (see Figure 1).Initially ( ) , the sheet is stretched with velocity along the X -axis, whereas the origin is kept fixed in the fluid medium of ambient temperature  T , and w T is the stretching surface temperature.A magnetic field of uniform strength ( ) is acting normal to the direction of the flow, with an acute angle .The magnetic Reynold's number is assumed as much less than unity, and the flow is considered two-dimensional, therefore the induced magnetic field can be neglected in comparison to the applied magnetic field.Moreover, the fluids exhibit polarization due to the incorporated principles of FHD, and the applied magnetic field is assumed to be strong enough to attain equilibrium magnetization.The velocity slip, thermal slip, viscous dissipation parameter, and ferromagnetic interaction parameter have been taken into account.The boundary layer equations of the fluid and energy equation for the problem can be written as [1,26,27].Conservation of mass: Conservation of mass: ∂u ∂x Momentum conservation: Energy (Heat) conservation: Here, ρ is the biomagnetic fluid density, , where B 0 is a constant representing the magnetic field strength at t = 0, k 1 (t) = k 2 (1 − ct) is the time dependent permeability parameter, k 2 is the constant permeability of the medium, g is the acceleration due to gravity, β * is the coefficient of thermal expansion, K is the thermal conductivity, C p is the specific heat at constant pressure, ϑ is the kinematic coefficient of viscosity, σ is the electrical conductivity, q r is the radiative heat flux, M 1 is the magnetization, H is the magnetic field of the fluid, and T is the temperature of the field.
The boundary conditions for the problem can be written as: [27,28] Here, U w (x, t) = ax (1−ct) is the stretching velocity, T w (x, t) = T ∞ + bx (1−ct) is the surface temperature.Where a, b, c are the constants such that a > 0, b > 0, c ≥ 0 and ct < 1.
In Equation ( 4), V w represents the blood velocity at the wall and is equal to injection/suction velocity given by As it is implied by Equation ( 5), the mass transfer at the sheet of the wall takes place with a velocity V w .For the case of injection, it is considered that V w > 0, whereas V w < 0 is considered for the case of suction.
is the thermal slip factor.The no-slip conditions hold when N = K = 0.
Using the Rosseland approximation [27], the radiation heat flux q r is simplified as where σ * and k * are the Stefan-Boltzman constant and the mean absorption co-efficient, respectively.
Considering that the temperature differences within the flow are such that the term T 4 may be expressed as a linear function of the temperature, T 4 is expanded in a Taylor series about T ∞ .By neglecting the higher order terms beyond the first degree in (T − T ∞ ), it is obtained that The non-uniform heat source/sink q m is defined as where A * and B * are the coefficient of a space-and temperature-dependent heat source/sink, respectively.The case of A * > 0, B * > 0 corresponds to internal heat generation, and that of A * < 0, B * < 0 corresponds to internal heat absorption.Physically, the role of a heat source in a fluid transport is to enhance its thermal conductivity, which consequently results in increased fluid temperature.On the other hand, heat sink decreases the thermal conductivity, which results in a decrease in the temperature of the fluid.By substituting ( 5) and ( 6) into (3), the energy equation is reduced to The term µ 0 M 1 ∂H ∂x in Equation (2) denotes the component of magnetic force per unit volume.This term is heavily dependent on the presence of a magnetic gradient, and, when the magnetic gradient is absent, this force vanishes.The heating due to adiabatic magnetization is represented by the second term, on the left hand side of the thermal energy Equation (8).The components H x and H y of the magnetic field H = H x , H y , which are due to a magnetic dipole, are given by [29,30] 2 where V = α 1 2π x x 2 +(y+d) 2 is a scalar potential of the magnetic dipole, γ = α 1 and α 1 is a dimensionless distance defined as Thus, the magnetic field strength intensity H = H is given by The corresponding gradients are given by The magnetization M 1 is generally determined by the fluid temperature provided that the applied magnetic field H is sufficiently strong enough to saturate the biomagnetic fluid.Anderson and Valnes [29] considered that the variation of magnetization M 1 with temperature T can be approximated by the linear equation where k is a constant.
To transform the momentum and energy equations, the following transformations are defined: Here, η is the similarity variable, ψ is the stream function, f and θ are dimensionless quantities.
The continuity Equation ( 1) is satisfied by the stream function ψ as u = ∂ψ ∂y and v = − ∂ψ ∂x Making use of Equation (9), Equations ( 2) and ( 8) can be written as The boundary conditions are transformed to: In Equation (12), S < 0 and S > 0 correspond to injection and suction, respectively.In the equation written above, primes denote derivatives with respect to η. S f = Nρ √ aϑ and S t = K a ϑ are the non-dimensional velocity slip factor and thermal slip factor, respectively.Furthermore, A = c a is the unsteadiness parameter, Pr = is the ferromagnetic interaction parameter, xϑ d is the dimensionless distance, and Re = xU w ϑ is the local Reynolds number.
The skin friction coefficient and the Nusselt number constitute important characteristics of the flow, defined as: where, the wall stress τ w and the heat transfer q w from the sheet are given by Using the similarity variables (9), it is obtained that:

Numerical Method
The numerical solution of the set of nonlinear ordinary differential Equations ( 10) and (11) subject to boundary conditions (12) is utilized by using bvp4c function technique MATLAB package.We consider f = y 1 , f = y 2 , f = y 3 , θ = y 4 , θ = y 5 .Then, the governing equations are transformed into the following system of first order ordinary differential equations: Along with the initial boundary conditions: Equations ( 13) and ( 14) are integrated numerically as an initial value problem to a given terminal point.All the calculations are made by using bvp4c function available in MATLAB software.

Parameter Estimation
In this study, the unsteady biomagnetic fluid flow along a two-dimensional stretching/shrinking sheet under the action of a magnetic field is investigated numerically.In order to achieve the numerical solution, it is necessary to determine some specific values for the dimensionless parameters, such as the Prandtl number, the unsteadiness parameter, the magnetic field parameter, the permeability parameter, the radiation parameter, the ferromagnetic interaction parameter, the Grashof number, the Eckert number, the suc-tion/injection parameter, the non-dimensional velocity-slip factor, the non-dimensional thermal slip factor, the acute angle of magnetic field, the inclination angle, the co-efficient of space and temperature.
Many researchers have, in scientific literatures, reported various values of the abovementioned dimensionless parameters.It is understood that, for human blood, the following data are considered: in [31][32][33], where, human body temperature is considered to be T w = 37 • C and the body Curie temperature is considered to be T ∞ = 41 • C. For this value of temperature, the dimensionless temperature is ε = 78. 5.

Results and Discussion
In order to assess the validity of the numerical results, the values of local Nusselt number −θ (0) have been compared with the existing works of Magyari and Keller [35], El-Aziz [36], Bidin and Nazar [8], and Anwar Ishak [10] by setting S = 0, 1 that the numerical scheme and the coding used give results in good agreement with the abovementioned, previously published studies.
A comparison of the local Nusselt number −θ (0) for various values R, M, Pr.Figures 2 and 3 show the velocity and temperature distributions with various values of the unsteadiness parameter A. From Figure 2, it is observed that the velocity profiles are decreased as the unsteadiness parameter is increasing.This is justified because the accompanying reduction in the thickness of the momentum in the boundary layer.Moreover, from Figure 3, it is obtained that the temperature profiles are decreased significantly as the unsteadiness parameter is increased.The fact is that, when the unsteadiness parameter is increased, less heat is transferred from the sheet to the fluid.2, it is observed that the velocity profiles are decreased as the unsteadiness parameter is increasing.This is justified because the accompanying reduction in the thickness of the momentum in the boundary layer.Moreover, from Figure 3, it is obtained that the temperature profiles are decreased significantly as the unsteadiness parameter is increased.The fact is that, when the unsteadiness parameter is increased, less heat is transferred from the sheet to the fluid.Figures 4 and 5 show the velocity and temperature profiles for various values of the radiation parameter R .From Figures 4 and 5, it is observed that an increment in radia- tion parameter R results in a decrement in the fluid velocity profile, whereas the tem- perature profile increases.The temperature profile is increased because the effect of the radiation parameter is to enhance heat transfer.The thermal boundary layer thickness is increased with the increment of the thermal radiation.f'() A=0.0,0.5,0.9    4 and 5, it is observed that an increment in radiation parameter R results in a decrement in the fluid velocity profile, whereas the temperature profile increases.The temperature profile is increased because the effect of the radiation parameter is to enhance heat transfer.The thermal boundary layer thickness is increased with the increment of the thermal radiation.
The velocity and temperature profiles with various values of Prandtl number Pr are shown in Figures 6 and 7.It is observed that an increment in Pr causes an increment in the velocity profile, whereas the temperature profile is decreased.This occurs because an increment in the Prandtl number means a decrement in the thermal diffusivity, and this phenomenon leads to the decreasing of energy ability that finally results in the reduction of the thermal boundary layer thickness.       of the thermal boundary layer thickness.
Figures 8 and 9 show the effect of the Grashof number Gr on the profiles of velocity and temperature.It was found that, with an increment in the Grashof number, which increases the velocity profile, the opposite is true for the temperature profile.This is due to the fact that an increase in the Grashof number means increment of the buoyancy forces which finally reduce the thermal boundary layer thickness.and temperature.It was found that, with an increment in the Grashof number, which increases the velocity profile, the opposite is true for the temperature profile.This is due to the fact that an increase in the Grashof number means increment of the buoyancy forces which finally reduce the thermal boundary layer thickness.Figures 8 and 9 show the effect of the Grashof number Gr on the profiles of velocity and temperature.It was found that, with an increment in the Grashof number, which increases the velocity profile, the opposite is true for the temperature profile.This is due to the fact that an increase in the Grashof number means increment of the buoyancy forces which finally reduce the thermal boundary layer thickness.
The velocity and temperature profiles for various values of the Eckert number Ec are shown in Figures 10 and 11.The relationship between the kinetic energy in the flow and the enthalpy is expressed by the Eckert number.It assimilates the conversion of kinetic energy into internal energy by the work done against the viscous fluid stresses.The positive Eckert number implies cooling of the sheet.Hence, greater viscous dissipative heat causes a rise in temperature as well as the velocity, both of which are evident in Figures 10 and 11    Figures 12 and 13 show the effect of the suction/injection parameter S on the ve- locity and temperature profiles.From the figures, it is observed that the momentum boundary layer thickness is decreased with increasing values of S .It is expected that the increment of the suction results in the decrement of the thickness of the hydrodynamic boundary layer.It also illustrates that an increment in S decreases the temperature pro- files in the flow region.This is due to the fact that, as the suction is increased, more warm fluid is taken away from the fluid region, causing a reduction in the thermal boundary layer thickness.Figures 12 and 13 show the effect of the suction/injection parameter S on the ve- locity and temperature profiles.From the figures, it is observed that the momentum boundary layer thickness is decreased with increasing values of S .It is expected that the increment of the suction results in the decrement of the thickness of the hydrodynamic boundary layer.It also illustrates that an increment in S decreases the temperature pro- files in the flow region.This is due to the fact that, as the suction is increased, more warm fluid is taken away from the fluid region, causing a reduction in the thermal boundary layer thickness.Figures 12 and 13 show the effect of the suction/injection parameter S on the velocity and temperature profiles.From the figures, it is observed that the momentum boundary layer thickness is decreased with increasing values of S. It is expected that the increment of the suction results in the decrement of the thickness of the hydrodynamic boundary layer.It also illustrates that an increment in S decreases the temperature profiles in the flow region.This is due to the fact that, as the suction is increased, more warm fluid is taken away from the fluid region, causing a reduction in the thermal boundary layer thickness.
Figures 14 and 15 show the effect of velocity slip parameter S f on the velocity and temperature profiles.From the Figure 14, it is observed that the presence of velocity slip within the boundary layer causes the velocity level along the sheet to decrease.This is happening because the quantity 1 − f (0) increases monotonically with S f .We also observe that the temperature profile decreases, as shown in Figure 15.14, it is observed that the presence of velocity slip within the boundary layer causes the velocity level along the sheet to decrease.This is happening because the quantity 1 −  ′ (0) increases monotonically with f S .We also ob- serve that the temperature profile decreases, as shown in Figure 15.
The effect of the thermal slip parameter t S on the velocity and temperature profiles is shown in Figures 16 and 17.In Figure 16, it can be observed that the presence of the thermal slip factor on the temperature profiles has a significant effect.It is clear that the temperature near the surface is decreased as the values of t S are increased.This is hap- pening because the increment in the thermal slip parameter results in the increment of the thermal coefficient, and the thermal diffusion towards the blood flow is reduced.The reverse behavior takes place for the velocity boundary layer, as shown in Figure 16.14, it is observed that the presence of velocity slip within the boundary layer causes the velocity level along the sheet to decrease.This is happening because the quantity 1 −  ′ (0) increases monotonically with f S .We also ob- serve that the temperature profile decreases, as shown in Figure 15.
The effect of the thermal slip parameter t S on the velocity and temperature profiles is shown in Figures 16 and 17.In Figure 16, it can be observed that the presence of the thermal slip factor on the temperature profiles has a significant effect.It is clear that the temperature near the surface is decreased as the values of t S are increased.This is hap- pening because the increment in the thermal slip parameter results in the increment of the thermal coefficient, and the thermal diffusion towards the blood flow is reduced.The reverse behavior takes place for the velocity boundary layer, as shown in Figure 16.The effect of the thermal slip parameter S t on the velocity and temperature profiles is shown in Figures 16 and 17.In Figure 16, it can be observed that the presence of the thermal slip factor on the temperature profiles has a significant effect.It is clear that the temperature near the surface is decreased as the values of S t are increased.This is happening because the increment in the thermal slip parameter results in the increment of the thermal coefficient, and the thermal diffusion towards the blood flow is reduced.The reverse behavior takes place for the velocity boundary layer, as shown in Figure 16.
Figures 18 and 19 show the effects of the permeability parameter k 3 on the velocity and temperature profiles.It is observed from Figure 18 that the presence of permeability parameter k 3 on the velocity profiles has a significant incremental effect.This is happening because the flow increases over the sheet as the permeability parameter is increased.The resistance on the flow above the sheet is decreased as the permeability of the sheet increases.From Figure 19, it can be noticed that the temperature profile declines when the permeability parameter k 3 enhances.temperature profiles.From Figure 20, we observe that the velocity profile decreases with an increment in the inclination parameter  .It seems that the angle of inclination de- creases the effect of the buoyancy force due to thermal diffusion by a factor of  cos .
Therefore, the driving force to the fluid decreases, and, as a result, the velocity is finally decreased.The reverse is happening in the temperature profile, which is shown in Figure 21.creases the effect of the buoyancy force due to thermal diffusion by a factor of  cos .
Therefore, the driving force to the fluid decreases, and, as a result, the velocity is finally decreased.The reverse is happening in the temperature profile, which is shown in Figure 21.Figures 20 and 21 show the effect of inclination parameter α on the velocity and temperature profiles.From Figure 20, we observe that the velocity profile decreases with an increment in the inclination parameter α.It seems that the angle of inclination decreases the effect of the buoyancy force due to thermal diffusion by a factor of cos α.Therefore, the driving force to the fluid decreases, and, as a result, the velocity is finally decreased.The reverse is happening in the temperature profile, which is shown in Figure 21.
Figures 22 and 23 show the effect of the inclination angle of the magnetic field ξ on the velocity and temperature profiles.It is noticed that the velocity profile is reduced and that the temperature profile is enhanced by the increment in the inclination angle.This may be due to the fact that a rise in the aligned angle makes the applied magnetic field stronger.
Figures 24 and 25 depict the effects of the magnetic parameter M on the velocity and temperature profiles.It is observed that the velocity profile is decreased as the magnetic parameter is increased.The increment of the magnetic parameter increases the introduced Lorentz force in the boundary layer, and, hence, the velocity profile in the boundary layer is decreased.An increment in the magnetic parameter would enhance the Lorentz force and, consequently, an augmentation of the Lorentz force opposes the flow, and the fluid motion is reduced.From Figure 25, it is noticed that the temperature profiles increase as the magnetic parameter increases.This indicates the fact that the introduction of the transverse magnetic field to an electrically conductive fluid gives rise to the Lorentz force.All these effects result in the increment of the temperature of the fluid.Figures 22 and 23 show the effect of the inclination angle of the magnetic field  on the velocity and temperature profiles.It is noticed that the velocity profile is reduced and that the temperature profile is enhanced by the increment in the inclination angle.This may be due to the fact that a rise in the aligned angle makes the applied magnetic field stronger.

OR PEER REVIEW 17 of 26
Figures 24 and 25 depict the effects of the magnetic parameter M on the velocity and temperature profiles.It is observed that the velocity profile is decreased as the magnetic parameter is increased.The increment of the magnetic parameter increases the introduced Lorentz force in the boundary layer, and, hence, the velocity profile in the boundary layer is decreased.An increment in the magnetic parameter would enhance the  Figures 22 and 23 show the effect of the inclination angle of the magnetic field  on the velocity and temperature profiles.It is noticed that the velocity profile is reduced and that the temperature profile is enhanced by the increment in the inclination angle.This may be due to the fact that a rise in the aligned angle makes the applied magnetic field stronger.
Figures 24 and 25 depict the effects of the magnetic parameter M on the velocity and temperature profiles.It is observed that the velocity profile is decreased as the magnetic parameter is increased.The increment of the magnetic parameter increases the introduced Lorentz force in the boundary layer, and, hence, the velocity profile in the boundary layer is decreased.An increment in the magnetic parameter would enhance the The velocity and temperatre profiles for various values of ferromagnetic interaction parameter β are shown in Figures 26 and 27.It is observed that the velocity of the fluid decreases with an increment of ferromagnetic number, whereas the temperature profile is increased in these cases.The region behind that ferromagnetic number is directly related to the celvin force, which is also known as the drug force.The results observed in Figures 24-27 are in accordance with those presented in [30,31].Figures 28-39 depict the skin friction coefficient and the rate of wall heat transfer with regard to the magnetic parameter for various values of the inclination angle of the sheet, angle of the magnetic field, radiation parameter, ferromagnetic parameter, unsteadiness parameter, and Eckert number.From the figures, it can be observed that skin friction decreases with increasing values of the inclination angle of the sheet and the acute angle of magnetic field, whereas the rate of the wall heat transfer is increased in these cases.Moreover, both skin friction and the rate of wall heat transfer are decreased with increasing values of the radiation parameter, ferromagnetic parameter, and unsteadiness parameter.Moreover, both skin friction and the rate of wall heat transfer are increased with increasing values of the Eckert number.
AppliedMath 2021, 1, FOR PEER REVIEW sheet, angle of the magnetic field, radiation parameter, ferromagnetic param iness parameter, and Eckert number.From the figures, it can be observed th decreases with increasing values of the inclination angle of the sheet and t of magnetic field, whereas the rate of the wall heat transfer is increased Moreover, both skin friction and the rate of wall heat transfer are decrease ing values of the radiation parameter, ferromagnetic parameter, and unstea eter.Moreover, both skin friction and the rate of wall heat transfer are increasing values of the Eckert number.AppliedMath 2021, 1, FOR PEER REVIEW sheet, angle of the magnetic field, radiation parameter, ferromagnetic param iness parameter, and Eckert number.From the figures, it can be observed th decreases with increasing values of the inclination angle of the sheet and t of magnetic field, whereas the rate of the wall heat transfer is increased Moreover, both skin friction and the rate of wall heat transfer are decrease ing values of the radiation parameter, ferromagnetic parameter, and unstea eter.Moreover, both skin friction and the rate of wall heat transfer are increasing values of the Eckert number.

Conclusions
The results concern the study of the effect of the magnetic field and r on the two-dimensional unsteady inclined stretching sheet with various s We observed that the suction/injection parameter has a good impact on th ber and skin friction coefficient.Such types of biological fluid flow problem ing in the biomedical and bioengineering sectors, especially in drug and cancer treatment, and MRIs.The important findings are given below: (1) Fluid velocity is reduced with ferromagnetic interaction parameter, ra eter, magnetic field parameter, inclination angle of the sheet, and acut netic field, whereas the temperature is increased in all cases;

Conclusions
The results concern the study of the effect of the magnetic field and radiation effects on the two-dimensional unsteady inclined stretching sheet with various slip conditions.We observed that the suction/injection parameter has a good impact on the Nusselt number and skin friction coefficient.Such types of biological fluid flow problems are interesting in the biomedical and bioengineering sectors, especially in drug and gene delivery, cancer treatment, and MRIs.The important findings are given below: (1) Fluid velocity is reduced with ferromagnetic interaction parameter, radiation parameter, magnetic field parameter, inclination angle of the sheet, and acute angle of magnetic field, whereas the temperature is increased in all cases; (2) Fluid velocity is enhanced with Prandtl number, Grashof number, permeability parameter, non-dimensional thermal slip factor, whereas temperature is decreased in all cases; (3) Both fluid velocity and temperature are decreased with unsteadiness parameter, suction/injection parameter, and non-dimensional velocity slip factor; (4) Both fluid velocity and temperature are increased with Eckert number, coefficient of space, and a temperature-dependent heat source/sink; (5) The effects of suction parameter on skin friction are enhanced with the increment of Eckert number, but decreased with an inclination angle of the sheet, acute angle of magnetic field, thermal slip factor, and unsteadiness parameter; (6) The effects of the injection parameter on the Nusselt number enhanced with the inclination angle of the sheet, acute angle of magnetic field, and Eckert number, but decreased with an increment in the radiation parameter and unsteadiness parameter; (7) In case of suction parameter, both skin friction and Nusselt number are decreased with increasing values of the ferromagnetic interaction parameter.

Figure 1 .
Figure 1.Physical sketch of the problem.

Figure 1 .
Figure 1.Physical sketch of the problem.

Figures 4
Figures4 and 5show the velocity and temperature profiles for various values of the radiation parameter R. From Figures4 and 5, it is observed that an increment in radiation

Figure 3 .
Figure 3. Temperature profile ( )   for various values of A .

Figure 5 . 5 
Figure 5. Temperature profile ( )   for various values of R .

Figure 4 .
Figure 4. Velocity profile f (η) for various values of R.

Figure 3 .
Figure 3. Temperature profile ( )   for various values of A .
Figures8 and 9show the effect of the Grashof number Gr on the profiles of velocity and temperature.It was found that, with an increment in the Grashof number, which increases the velocity profile, the opposite is true for the temperature profile.This is due to the fact that an increase in the Grashof number means increment of the buoyancy forces which finally reduce the thermal boundary layer thickness.The velocity and temperature profiles for various values of the Eckert number Ec are shown in Figures10 and 11.The relationship between the kinetic energy in the flow and the enthalpy is expressed by the Eckert number.It assimilates the conversion of kinetic energy into internal energy by the work done against the viscous fluid stresses.The positive Eckert number implies cooling of the sheet.Hence, greater viscous dissipative heat causes a rise in temperature as well as the velocity, both of which are evident inFigures 10 and 11.

Figure 9 .
Figure 9. Temperature profile ( )   for various values of Gr .

Figure 10 .
Figure 10.Velocity profile f (η) for various values of Ec.

Figure 11 .
Figure 11.Temperature profile ( )   for various values of Ec .

Figure 11 .
Figure 11.Temperature profile ( )   for various values of Ec .

Figure 13 .
Figure 13.Temperature profile ( )   for various values of S .

Figure 15 .
Figure 15.Temperature profile ( )   for various values of f S .

Figure 21 .
Figure 21.Temperature profile ( )   for various values of  .

Figure 21 .
Figure 21.Temperature profile ( )   for various values of  .

Figure 23 .
Figure 23.Temperature profile ( )   for various values of  .

Figure 23 .
Figure 23.Temperature profile ( )   for various values of  .

Figure 23 .
Figure 23.Temperature profile ( )   for various values of  .

Figure 27 .
Figure 27.Temperature profile ( )   for various values of  .

Figure 27 .
Figure 27.Temperature profile ( )   for various values of  .

Figure 27 .
Figure 27.Temperature profile ( )   for various values of  .

Figures 28 - 4 f 4 Figure 27 .
Figures 28-39 depict the skin friction coefficient and the rate of wall heat transfer with regard to the magnetic parameter for various values of the inclination angle of the

Figure 28 .
Figure 28.Skin friction coefficient f (0) with M for different values of α.

Figure 30 .
Figure 30.Skin friction coefficient f (0) with M for different values of ξ.

Figure 32 .
Figure 32.Skin friction coefficient f (0) with M for different values of R.

Figure 34 .
Figure 34.Skin friction coefficient f (0) with M for different values of β.

Figure 36 .
Figure 36.Skin friction coefficient f (0) with M for different values of A.

Figure 37 .
Figure 37. Local Nusselt number −θ (0) with M for different values of A.

Figure 38 .
Figure 38.Skin friction coefficient f (0) with M for different values of Ec.

Table 1 .
Comparisons with previous studies.Figures 2 and 3 show the velocity and temperature distributions with various values of the unsteadiness parameter A .From Figure

1 ,
[30,31]PEER REVIEW 18 of 26to the celvin force, which is also known as the drug force.The results observed in Figures 24-27 are in accordance with those presented in[30,31]