# Hydromagnetic Flow of Prandtl Nanofluid Past Cylindrical Surface with Chemical Reaction and Convective Heat Transfer Aspects

Authors: Kottakkaran Sooppy Nisar; S. Bilal; Imtiaz Ali Shah; M. Awais; Khalil-Ur-Rehman; Ilyas Khan; Phatiphat Thonthong
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

Thermodynamical attributes of non-Newtonian fluids over stretched surfaces have gained pervasive essence due to extensive utilization in extruding plastic sheet procedures, liquid film condensation, glass blowing, paper production, biopolymer cylinder coatings, and so forth. So, currently communication is aimed to candidly explicate flow characteristic of Prandtl fluid generated by axial stretching of cylindrical surface. Mathematical modelling by using conservation laws of momentum, energy and concentration fields containing the aspects of magnetic field, convective heating, and chemical reaction are presented initially in the form of partial differential expressions. Later on, these attained PDEs are transmuted into nonlinear ordinary differential equations with implementation of similarity variables. Numerical approach renowned as shooting technique with improved coefficient of the Runge–Kutta (R–K) method by Cash and Karp is used to access accurate solution. Linear curved fitting analysis is also performed to analyze results. Influence of flow-controlling parameters on associated profiles is revealed through graphical visualization. Stream line plots representing flow behavior of Prandtl fluid versus different magnitudes of the curvature parameter are adorned. Variation in friction drag force at wall, heat flux, and concentration gradient are evaluated through numerical data and with interpolation of linear curved fittings. It is deduced from results that increasing curvature parameter momentum and temperature distributions enriches whereas skin-friction coefficient depicts reverse pattern. It is also inferred that temperature shows incrementing deviation in the absence of chemical reaction whereas concentration profiles exhibit reduction with consideration of influence of chemical reaction parameter. Magnetic field tends to reduce the velocity and create thinness of boundary layer thickness.


---

## Body

## 1. Introduction

Analysis of non-Newtonian fluids has attained superb attention of researchers due to numerous scientific implementations in multiple fields such as food mixing, multigrade oils, composition of materials, wire drawing, hot rolling, petroleum productions, metallurgical procedures, manufacturing of materials, preventive coating, lubricating products, polymerization processes, ink-jet printers, geophysical flows, liquid crystallizations, and several others. Taking into consideration the aforementioned extensive applicability in different technological and industrial frameworks, a lot of analysis has been conducted in this direction in recent years. However,, the complexity about the non-Newtonian model is its mathematical modelling on behalf of different rheological features. Hence, there is no single constitutive mathematical model which is present, which expresses the fluid characteristics of all non-Newtonian fluids single handedly. Consequently, non-Newtonian fluid models are categorized into time-dependent, time-independent, viscoelastic, and viscoinelastic fluids. Among these subclasses, the most practicing fluids are viscoinelastic fluids. Dunn [1] presented pioneering investigation on the properties of viscoinelastic fluids and depicted behavior of strain of such fluids against application of stress. In the literature, a lot of fluid models are presented to elaborate characteristics of viscoinelastic fluids, but it is found that Prandtl fluid is the most fittest to explicate features of viscoinelastic fluids. Peristaltic movement of Prandtl fluid in endoscope was investigated by Nadeem et al. [2] by calculating perturbation solution of modelled problem with implementation of analytical approach. Akbar [3] considered blood as viscoinelastic fluid and delineated flow features in stenosed tapered arteries by computing perturbation solution. Jothi et al. [4] discussed the physical aspects of magnetically effected Prandtl fluid flow in a symmetric channel. Nazari et al. [5] analyzed convective heat transfer of Prandtl nanofluid flow saturated in lid-driven enclosure. Hossainy et al. [6] discussed the influence of yielding stress on magneto hydrodynamic boundary layer fluid flow of 3D Casson nanofluid over stretched surface entrenched in a porous medium. Bilal et al. [7] adumbrated flow features of the Prandtl fluid model over stretching sheet along with effectiveness of double-diffusive heating.The mass transport phenomenon in many scientific disciplines involves molecular and convective transport of atoms and molecules. The driving force for producing such type of molecular diffusion is the concentration differences which can generate to application of chemical reaction within flow domain. The process of chemical reaction arises with addition of multiple dissimilar natured particles, and an outcome interaction of particles is produced. Mass transport procedures with effectiveness of chemical reactions have many practical applications, e.g., combustion system, metallurgy, distribution of fertilizers in agriculture field, nuclear reactor safety, and many more. In view of such extensive applications several interesting studies are conducted. For the sake of brevity, some of them are mentioned. Effect of chemical reactions on viscous fluid flow over horizontal plate was discussed by Anjalidevi and Kandasamy [8]. Zhang and Zheng [9] reported influence of chemical reactive species in thermosolutal Marangoni flow of electrified viscous fluid and predicted reduction in nanoparticle concentration with respect to uplift in chemical reaction magnitude. Chaudhary and Merkin [10] investigated thermophysical attributes of viscous fluid influenced by chemical reaction effects. Nandkeolyar et al. [11] described stagnant dissipative nanofluid in the attendance of chemically reaction and magnetic field. Hossainy and Eid [12] analyzed hydrothermal efficiency of non-Newtonian hybrid nanofluid in a heat-exchange channel by generation of chemical reactions in domain. Eid and Makinde [13] studied collaborative aspects of radiative heat energy and chemical reaction on electrically conducting nanofluid flow over a stretching sheet immersed in a permeable medium. The effects of slip velocity and chemical reaction generated in stagnant flow of nanofluid over a stretching sheet embedded in a porous medium were studied by Eid [14]. In precise, investigators are still working to explore hidden features of chemical reaction on Newtonian and non-Newtonian flows [7, 15, 16].In the today technological world, nanotechnology is considered one of momentous inventions which have brought advancement in industrial rebellion. Nanotechnology based on nanofluids is engineered colloidal suspension of nanometer-sized particles in host fluids. Before the discovery of nanoparticles, liquids with low thermophysical characteristics are capitalized in multiple procedures which produce deficiencies in output. So, Choi [17] proposed the idea of nanofluids by experimentally measuring their intrinsic features possessed by nanoparticles. Buongiorno [18] evaluated enhancement in thermal conductance of ordinary fluids with the induction of nanoparticles. Flow analysis of viscous fluid generated by an exponential stretching of surface by adding nanoparticles was probed by Nadeem et al. [19]. Reza et al. [20] interrogated influence of carbon nanotubes in peristaltic movement of water in a permeable channel under the impact of magnetic field. Peyghambarzadeh et al. [21] disclosed application of nanofluids in increasing automotive cooling. In recent years, extensive work on nanofluid flow analysis in numerical computational domains along with applications was manifest. So, for interest of readers and getting their motivation in this work few references are mentioned. Alaidrous and Eid [22] adumbrated nanofluid flow of non-Newtonian liquid under impedance of electromagnetic radiations and highly rated reactions. Eid and Mabood [23] irreversibility entropy generated aspects of carbon nanotube (CNT) placement in micropolar dusty fluid mixed with kerosene oil. The uplift in thermal characteristics of blood flow with induction of gold plated hybridized nanoparticles past a stretching sheet was analyzed by Eid [24]. Also, [25–30] presented the literature discusses by various researchers on nanoliquids.Analysis of fluid flows under impression of magnetic field has been demanding research area among researchers due to overwhelming implications in industry, engineering, medical technology, and so many. In recent years, magneto hydrodynamics of electrically conducting flows has been investigated extensively in multiple physical configurations. Likewise, Lahmar et al. [31] examined transportation of an unsteady ferro magnetic nanofluid between two parallel plates by considering thermal heat conduction as a function of temperature in the presence of magnetic source. Some recent investigations in this direction are enclosed in [32, 33].The abovementioned literature survey reflects prominence about thermophysical features of non-Newtonian fluid flow. However, it is seen that analysis of the Prandtl fluid model in the presence of chemical reaction and magnetic field over melted heating cylinder has not been scrutinized yet. So, current effort is made to fill this gap. In this regards, mathematical modelling of the non-Newtonian fluid rheological model with heat transfer aspects is attained in the form of partial differential equations by applying boundary layer approach. Afterwards, suitable scaling transformations are used for conversion into ordinary differential expressions. Improved numerical simulations by means of the Runge–Kutta method with improved coefficients presented by Cash and Karp are used. Computed results are presented through graphical and tabular manners.


## 2. Mathematical Formulation

Assume 2-D, in-compressed, time-independent, and electrified Prandtl nanoliquid persuaded by an axially stretchable cylinder. Letv and u be the velocity constituents along r and x directions, respectively, as shown in Figure 1. The cylinder is stretched axially with linear velocity. The magnetic field of strength B°is applied perpendicular to cylinder axis. Furthermore, the generated magnetic field is very small in comparison with applied one due to low magnetic Reynolds number assumption.Figure 1
Diagram of physical configuration of problem.Prandtl proposed the following rheological model [16]:(1)S=Aγ˙sin−11cγ˙A1,where fluid parameters are denoted with A and c, A1 is first Rivlin–Ericksen tensor, and γ˙ is defined as(2)γ˙=12trA12.By expandingsin−1⋅ function up to second-order approximation using Taylor’s series,(3)sin−1γ˙c=γ˙c+γ˙36c3,where1cγ˙<1.By means of usual boundary layered approximations, the constitutive expressions become(4)∂ru∂x+∂rv∂r=0,(5)u∂u∂x+v∂u∂r=Aρc∂2u∂r2+A2ρc3∂u∂r2∂2u∂r2+Arρc∂u∂r+A6rρc3∂u∂r3−σB02ρu,(6)u∂T∂x+v∂T∂r=α˙r∂∂rr∂T∂r+τDB∂C∂r∂T∂r+τDTT∞∂T∂r2+Acρcp∂u∂r2+A6c3ρcp∂u∂r3,(7)u∂C∂x+v∂C∂r=DBr∂∂rr∂C∂r+DTrT∞∂∂rr∂T∂r−KC−C∞.Accompanied to boundary conditions,(8)u=Uwx+L∂u∂r,v=0,kT∂T∂r=−hwTw−T,C=Cmx,atr=R,u⟶0,T⟶T∞x,C⟶C∞,asr⟶∞.In the above equationsσ is the electrical conductance, κ is the thermal conductivity, R is the radius of the cylinder, α˙=υ/ρcp denotes the thermal diffusivity, cp is the specific heat at constant pressure, DB is Brownian motion diffusion coefficient, DT is the thermophoresis diffusion coefficient, K is chemical reaction rate constant,Uwx=ax is linear stretching velocity, and wall heat flux coefficient is denoted with hw.The stream function which satisfied the continuity equation is mentioned as follows:(9)u=1r∂ψ∂r,v=−1r∂ψ∂x.Equations (4)–(7) can be transformed into nonlinearized ordinary differential equation by employing the following transmutations:(10)η=r2−R22RUwν1/2ψ=Uwνx1/2Rfηθη=T−T∞Tw−T∞φη=C−C∞Cw−C∞.After applying the abovementioned transformations, the governing equations, i.e., equations (4)–(7), are transformed into following form:(11)α1+2γηf‴+αβ1+2γη2f″2f‴+43γαβ1+2γηf″3+2γαf″−f′2+ff″−Mf2f′=0,1+2γηθ″+2γθ′+Prfθ′+Pr1+2γηNbθ′ϕ′+Ntθ′2=0,1+2γηϕ″+2γϕ′+LePrfϕ′−δcϕ+NtNb2rθ′+1+2γηθ″=0.The associated boundary conditions take the following form:(12)fη=0,f′η=1+λf″0,atη=0,f′η⟶0,atη⟶∞,θ′η=−k1−θη,atη=0,θη⟶0,asη⟶∞,ϕη=1,atη=0,ϕη=0,atη⟶∞,where αis the fluid parameter, β is the elastic parameter,γ is the curvature parameter, Nb is the Brownian motion parameter, Mf is the magnetic strength parameter, Pr is the Prandtl number, Re is the Reynolds number, Le is the Lewis number, Nt is the thermophoresis parameter, and δc is the chemical reactive parameter defined as follows:(13)α=Aμcβ=a3x22c2υγ=1RνaNb=τDBCw−C∞ρcpυδc=KLaMf=σxB02ρUwPr=να˙Re=UwxυLe=α˙DBNt=τDTTw−T∞T∞ρcpυ.The physical parameters of interest, i.e., coefficients of wall drag, wall heat flux, and wall mass flux, are defined as(14)Cf=τw1/2ρUw2,Nuw=xqwkT∞−Tm,Sh=xjwC∞−Cm,where τw is the surface drag force, qwis the surface heat flux, and jw is the surface mass flux given as follows:(15)τw=AC∂u∂r+A6C3−∂u∂r3r=R,qw=−k∂T∂rr=R,Jw=−D∂C∂r.Coefficients of wall drag, thermal flux, and mass flux are defined as(16)12CfRe−1/2=αf″η+αβ3f″η3,atη⟶0NuwRe−1/2=−θ′η,atη⟶0ShRe−1/2=−ϕ′η,atη⟶0.


## 3. Computational Algorithm

Shooting technique along with the Runge–Kutta–Fehlberg method (by using Cash and Carp coefficients) is adopted to solve flow-govern problem (Table1). For this process, flow-govern equations are transformed into(17)f″=f′2−ff″−2γαf″+Mf2f−4/3γαβ1+2γηf″3α1+2γη+αβ1+2γη2f″2,θ″=−Prfθ′+Pr1+2γηNbθ′ϕ′+Ntθ′2+2γθ′1+2γη,ϕ″=−2γϕ′+LePrfϕ′−δcϕ+Nt/Nb2rθ′+1+2γηθ″1+2γη.Table 1
Numerical values of the coefficient given by Cash and Karp.
IAiBijCiDi1––––––37/3782825/2764821/51/5––––0033/103/409/40–––250/62118575/4838443/53/10−9/406/5––125/59413525/5529651−11/545/2−70/2735/27–0277/1433667/81631/55296175/512575/1382444275/110592253/4096512/17711/4By letting a new set of variables,(18)e1=f,e2=f′,e3=f″,e4=θ,e5=θ′,e6=ϕ,e7=ϕ′.After using new variables, the governing equations take the following form:(19)e1′=e2,(20)e2′=e3,(21)e3′=e22−e1e3−2γαe3+Mf2e1−4/3γαβ1+2γηe33α1+2γη+αβ1+2γη2e32,(22)e4′=e5,(23)e5′=−Pre1e5+Pr1+2γηNbe5e7+Nte52+2γe31+2γη,(24)e6′=e7,(25)e7′=−2γe7+LePre1e7−δce6+Nt/Nb2γe5+1+2γηe5′1+2γη.along boundary conditions(26)e1η=0,e2η=1+λe3η,e4η=−1/k,e5η−1,e6η=1,atη⟶0,e2η⟶0,e4η⟶0,e6η⟶0,asη⟶∞.The Runge–Kutta–Fehlberg method is applicable on initial value problem, and in addition, Cash and Karp coefficients are used to raise the accuracy of solution. So, the given boundary conditions are transformed into initial conditions:(27)e10=0,e20=1+λe30,e30=ω1,unknown initial condition,e40=−1ke50−1,e50=ω2,unknown initial condition,e60=1,e70=ω3,unknown initial condition.Then, initial value problem is solved and error is computed by comparing given and calculated boundary values. If error is larger than tolerance, i.e.,10−6, then initial values are modified with Newton’s method and process is repeated unless differences are less than 10−6.


## 4. Graphs and Tables

Impacts and influences of different variables in absence ofδcare shown in Figures 2–14.Figure 2
Impacts ofα on f′ηin absence/presence of δc.Figure 3
Influence ofβ on f′ηin absence/presence of δc.Figure 4
γ influences on f′ηin absence/presence of δc.Figure 5
Mf impacts on f′ηin absence/presence of δc.Figure 6
γ impacts on θηin the absence/presence of δc.Figure 7
Pr impacts on θηin the absence/presence of δc.Figure 8
Nb influences on θηin the absence/presence of δc.Figure 9
Nt influences onθηin the absence/presence of δc.Figure 10
Nb impacts on ϕηin the absence/presence of δc.Figure 11
Nt influences onϕηin the absence/presence of δc.Figure 12
Le influences onϕηin the absence/presence of δc.Figure 13
Pr impacts on θηin the absence/presence of δc.Figure 14
(a) Stream lines against curvature parameter. (b) Stream lines against curvature parameter. (c) Stream lines against curvature parameter.
(a)(b)(c)


## 5. Straight Line Curve-Fitting Analysis

The behavior of wall drag, wall heat flux, and wall mass flux is analyzed via straight linearized curved fitting for Prandtl fluid parameter, elastic parameter, curvature parameter, velocity slip, magnetic field, Prandtl number, thermophoretic, Brownian motion, Biot number, Lewis number, and chemical reactive parameters. The Guass and Legendre method was adopted. The expressions for linear approximations are defined such as(28)na0+a1∑χi=∑ξi,a0∑χii+a1∑χi2=∑χiξi.We trace out linearized fitting for wall drag coefficient towards Prandtl fluid, elastic, curvature, velocity slip, and magnetic parameters, i.e.,α,β,γ,λ,andMf. Let χi=αi and ξi=αf″0+αβ/3f″03i, then we get(29)∑αi=0.7,∑αi2=0.25,∑αf″0+αβ3f″03i=−1.3085,∑αiαf″0+αβ3f″03i=−0.4826.For straight line approximation, we inserted above expressions in equation (23) and got(30)2a0+0.7a1=−1.3085,0.7a0+0.25a1=−0.4826.By evaluating systemized expression given by (24), we get(31)αf″0+αβ3f″03=a0+a1α,where a0=1.0694anda1=−4.9249.Linearized curve fitting for wall shear stress versus elastic parameterβ is found such as(32)∑βi=0.7,∑βi2=0.29,∑αf″0+αβ3f″03i=−1.1708,∑βiαf″0+αβ3f″03i=−0.4632,2b0+0.7b1=−1.1708,0.7b0+0.29b1=−0.4632.The linearized approximation for fluid elastic variable is specified by(33)αf″0+αβ3f″03=Pβ=b0+b1β,where b0=−0.1699andb1=−1.1871.The straight line-curved fitted model for the dragged wall number for curvature parameterγ is disclosed as follows:(34)∑γi=1.5,∑γi2=1.25,∑αf″0+αβ3f″03i=−1.5347,∑βiαf″0+αβ3f″03i=−1.3299,2c0+1.5c1=−1.5347,1.5c0+1.25c1=−1.3299.Then, the linear approximations for the curvature parameter are given by(35)αf″0+αβ3f″03=Pγ=c0+c1γ,where c0=−0.3059andc1=−1.4310.Straight line curve fitting for wall drag coefficient against velocity slip parameterλ is defined as(36)∑λi=0.4,∑λi2=0.1,∑αf″0+αβ3f″03i=−0.6685,∑λiαf″0+αβ3f″03i=−0.1191,2d0+0.4d1=−0.6685,0.4d0+0.1d1=−0.1191.Then, the straight line curves versus the curvature parameter are defined as follows:(37)αf″0+αβ3f″03=Pλ=d0+d1λ,where d0=−0.4802andd1=0.7300.Straight line curve fitting for skin-friction coefficient towards velocity slip parameterλ is computed as(38)∑Mfi=0.4,∑Mfi2=0.1,∑αf″0+αβ3f″03i=−0.9699,∑Mfiαf″0+αβ3f″03i=−0.2095,2e0+0.4e1=−0.9699,0.4e0+0.1e1=−0.2095.Then, the straight line approximations towards the curvature parameter are given by(39)αf″0+αβ3f″03=PMf=e0+e1Mf,where e0=−0.3297ande1=−0.7759.Now, to trace curve-fitted approximation for convectional wall transfer against the curvature parameter, Prandtl number, thermophoresis parameter, Brownian movement parameter, and Biot number are calculated as follows:(40)∑γi=1.5,∑γi2=1.25,∑−θ′0i=0.5973,∑γi−θ′0i=0.4764,2f0+1.5f1=0.5973,1.5f0+1.25f1=0.4764.After simplifying, we get(41)−θ′0=Pγ=f0+f1γ,where f0=0.1281andf1=0.2274.Wall heat flux straight line approximation towards the Prandtl number is assessed as(42)∑Pri=3,∑Pri2=5,∑−θ′0i=0.6737,∑Pri−θ′0i=2.8851,2g0+8g1=0.6737,8g0+34g1=0.4320.After computing by usual algebraic practice, we attain(43)−θ′0=PPr=g0+g1Pr,where g0=4.8624andg1=1.1314.Linear approximation for heat transferal rate towards the thermophoresis parameter is calculated by(44)∑Nti=0.1,∑Nti2=0.2,∑−θ′0i=0.4558,∑Nti−θ′0i=0.0669,2h0+0.3h1=0.4558,0.3h0+0.05h1=0.0669.After solving by common algebraic practice, we get(45)−θ′0=PNt=h0+h1Nt,where h0=4.8624andh1=1.1314.Straight line approximation by varying thermophoresis parameter for wall heat flux is found as(46)∑Nbi=0.1,∑Nbi2=0.2,∑−θ′0i=0.2417,∑Nbi−θ′0i=0.2212,2j0+0.3j1=0.4829,0.3j0+0.05j1=0.0644.After simplifying, we get(47)−θ′0=PNb=j0+j1Nb.where j0=0.4825andj1=−1.6070.Straight line approximations of wall heat flux against the Biot number are found as(48)∑ki=0.4,∑ki2=0.1,∑−θ′0i=0.2417,∑Nbi−θ′0i=0.2212,2l0+0.4l1=1.2262,0.4l0+0.1l1=0.2257.After computing by usual algebraic practice, we attain(49)−θ′0=Pk=l0+l1k,where l0=1.0694andl1=−4.9242.Wall mass flux straight line approximations versus Lewis number, Prandtl number, thermophoresis parameter, Brownian motion parameter, and chemical reaction parameter are examined as follows:(50)∑Lei=1.5,∑Lei2=1.25,2m0+1.5m1=1.3397,1.5m0+1.25m1=1.0665.After simplifying, we get(51)−ϕ′0=PLe=m0+m1k,where m0=0.2995andm1=0.4938.Straight line approximation for surface mass flux for the Prandtl number is calculated as(52)∑Pri=1.5,∑Pri2=1.25,∑−ϕ′0i=1.1471,∑Pri−ϕ′0i=0.8739,2n0+1.5n1=1.2262,1.5n0+0.1n1=0.8739.After computing by usual algebraic way, we attain(53)−ϕ′0=PPr=n0+n1Pr,where n0=0.4921andn1=1.0847.Linear approximation for mass flux towards thermophoresis variable is scrutinized as follows:(54)∑Nti=0.3,∑Nti2=0.05,∑−ϕ′0i=1.1587,∑Nti−ϕ′0i=0.1771,2q0+0.3q1=1.1587,0.3q0+0.05q1=0.1771.After computing by usual algebraic way, we attain(55)−ϕ′0=PNt=q0+q1Nt,where q0=0.4825andq1=0.6589.Surface mass flux straight line approximation for chemical reaction parameter is computed as(56)∑δci=0.3,∑δci2=0.05,∑−ϕ′0i=1.1505,∑δci−ϕ′0i=0.1754,2s0+0.3s1=1.1505,0.3s0+0.05s1=0.1754.After simplifying, we get(57)−ϕ′0=Pδc=s0+s1δc,where s0=0.4905ands1=0.5649.


## 6. Physical Outcomes

This section is fascinated to examine impact of embedded parameters on fluid concerning profiles. Moreover, physical quantities of interest, i.e., coefficients of wall friction, wall heat flux, and wall mass flux, are reported by means of both straight line curve fitting and tabular structure. In detail, impacts ofα and β on fluid velocity in the absence/presence of δc are presented in Figures 1–3, respectively. From graphical trends, it is seen that both α and β cause declination in momentum of fluid flow and associated boundary layer thickness. The reason behind this variation is that with increase in α and β, viscosity of fluid increases, and hence, movement of fluid decays. It is also divulged that flow regime depicts higher magnitude in the absence of chemical reaction as compared to the presence of chemical reaction. These trends in similar were reported by Eid et al. [14] and Nadeem and Lee [19]. Figure 4 reveals the variation of fluid velocity towards γ. The fluid velocity accelerated by increasing γ (see Figure 4) because larger values of γ correspond to reduce cylindrical radius which in turn reduces surface area; hence, fluid particles experienced less resistance due to which motion is accelerated. Figure 5 elucidates the impact of Mf on velocity. It is clear from these curves that the fluid velocity shows decline trend because an increase inMf leads to an increase in Lorentz force and it has tendency to resist the movement of fluid particles, so that velocity decreases. Figure 6 interprets the variations in temperature profile by varyingγ. This figure predicts that the fluid temperature increases by enlarging the curvature parameter. This is due to the fact that when γ enlarges, the surface heat flux increases due to which temperature drops adjacent to the cylindrical surface but temperature away from the cylindrical surface enhances. Figure 7 shows impact of Pr on fluid temperature. The graphical plotting depicts decrement in temperature distribution against Pr. It is also noticed that δc=0 (in absence of chemical reactions), and fluid temperature declines by enlarging the values of Pr, but for δc>0, the decline in fluid temperature is more rapid by enlarging the values of Pr. Figure 8 deliberates the fluctuations in fluid temperature by varying Nb in the absence/presence of the chemical reaction phenomenon. Since Brownian motion accelerates the particle movement, thus collision process is expedited and hence produces more thermal energy (which can be observed from the graph that fluid temperature rises for larger values of Nb in both cases). On the other hand, presence of chemical reactive species drops the fluid temperature. Figure 9 shows the temperature profile variations versus Nt in the absence/presence of chemical reaction. Temperature profile shows increase versus Nt both zero and nonzero δc. Figure 10 displays the variations in concentration profile versus combined effect of Nb and δc. This graph reveals that the larger values of Nblessen the nanoparticle concentration because enhancement in Nb corresponds to accelerate random motion of nanoparticles. However, in presence of destructive chemical reaction case, the concentration profile increases. Figure 11 reveals the effect of Nt in absence and presence of the chemical reaction parameter. It is clear from these curves that both parameters have the same impact on nanoparticle concentration. Figure 12 divulges the influence of Le on nanoparticle volume fraction in absence/presence of chemical reaction. This figure predicts that the concentration profile enlarges versus Le irrespective of absence/presence of chemical reaction. Figure 13 examines the physical significance of Pr on concentration profile in absence/presence of chemical reaction. It is evident that by increasing Pr and δc, nanoparticle volume fraction increases significantly. Figures 14(a)–14(c) are provided to inspect the flow pattern of Prandtl fluid over the stretching cylinder by sketching stream lines against curvature parameter γ. From stream line plots, it is seen that intensification in magnitude of lines is found with increase in γ because as curvature parameter γ increases, contact of fluid particles reduces, and as an outcome, less resistance is provided. So, rapid movement in fluid is observed versus γ. Furthermore, these curves show that γenhances the fluid momentum versus positive alteration in velocity. Figure 15 deliberates variation in skin-friction coefficient againstα,β,γ, λ,andMf. It is observed that with increase in α,β,andMf, wall shear stress increases because all of these mentioned parameters uplift the viscosity forces for reduction in velocity whereas reverse pattern is depicted in case of λandγ. The reason behind increase in skin friction is the uplift of viscosity force in flow domain against these mentioned parameters. Figure 16 displays fluctuations in wall heat flux by varying γ,Pr while Nt,Nb,andk. It is found that wall heat flux coefficient increases by enhancing γ,Pr while forNt,Nbandk reduce. The reason behind enhancement in temperature gradient at wall against γandPr is due to the fact that by increasing these parameters, resistance is provided by surface to molecules decrease, and hence, average kinetic energy of molecules increases which raise thermal flux. Figure 17 explicates variation in mass flux against Le,Pr,Nt,Nb,andδc. It is seen that all mentioned parameters diminish concentration distribution and create motion in fluid molecules, which raises the mass flux with in flow domain. Table 2 delineates the variations in wall friction coefficient versusα,β,γ,λ,and,Mfand increasing pattern of wall friction coefficient against α,β,γ while both λandMf decrease it. Table 3 is designed to explain the fluctuations in wall heat flux by varying governing parameters. Wall heat flux coefficient increases by enhancing γ,Pr while Nt,Nb,andk reduce it. Table 4 depicts the variations in the local Sherwood number viaLe,Pr,Nt,Nb,andδc, and we observed that with respect to defined parameters, enhancement in the local Sherwood number is manipulated. Tables 5 and 6 provide comparison with the previously published literature and given assurance about collected values in this work.Figure 15
Straight line curve fitting for wall shear stress.Figure 16
Straight line curve fitting for heat transfer rate.Figure 17
Straight line curve fitting for wall mass flux.Table 2
Numerical variation of wall friction coefficient for different parameters.
αβγλMfαf″0+αβ/3f″030.30.20.50.10.1−0.40730.4––––−0.90120.5––––−1.41210.30.2–––−0.4073–0.4–––−0.6481–0.5–––−0.7635–0.20.5––−0.4096––1––−1.1251––1.5––−1.7481––0.50.1–−0.4073–––0.3–−0.2612–––0.5–−0.1815–––0.10.1−0.4073––––0.3−0.5626––––0.7−0.2253Table 3
Numerical variation of HTR for different parameters.
γPrNtNbk−θ′00.530.10.10.50.24171––––0.35561.5––––0.45170.53–––0.2417–5–––0.4320–7–––0.7121–30.1––0.2417––0.2––0.2141––0.3––0.1782––0.10.1–0.2417–––0.2–0.2212–––0.3–0.2031–––0.10.10.7104––––0.30.5158––––0.70.1715Table 4
Numerical variation of MTR for different parameters.
LePrNtNbδc−ϕ′00.50.50.10.10.10.54631––––0.79341.5––––0.99460.50.5–––0.5463–1–––0.6008–1.5–––0.6458–0.50.1––0.5463––0.2––0.6124––0.3––0.7314––0.10.1–0.5463–––0.2–0.8201–––0.3–1.1241––––0.10.5463––––0.20.6042––––0.30.6513Table 5
Comparison of wall shear stress by varying magnetic field parameterMf while considering other parameters to zero.
MfAkbar at el. [32]Malik at el. [33]Present0.0−1−1−10.5−1.11803−1.11803−1.11801.0−1.41421−1.41419−1.41375.0−2.44949−2.44945−2.449510−3.31663−3.31657−3.3166100−10.04988−10.04981−10.0500500−22.38303−22.38294−22.38351000−31.63859−31.63851−31.6391Table 6
Comparison of results for the reduced Nusselt number−θ′0against Pr and considering all other parameters zero.
PrPresent resultsWang [30]Gorla and Sidawi [31]0.070.06630.06560.06560.200.16190.16190.16190.700.45390.45390.45392.000.91130.91140.91147.001.89541.89541.980520.003.35393.35393.353970.006.46216.46226.4622


## 7. Key Results

The current study investigates the magnetic field influences on partial slip flow of Prandtl nanofluid over the stretching cylinder. The fluid flow regime is manifested with the chemical reaction phenomenon. The modelled mathematical equations are reduced in terms of the system of ODEs, and then, the numerical method shooting is implemented to interpret physical aspects of involved parameters. In the presence or absence of the chemical reaction phenomenon, the main findings of present analysis are itemized as follows:(1)
Prandtl fluid velocity reflects an inciting trend towards both Prandtl fluid and elastic parameters(2)
Both fluid velocity and temperature are increasing function of curvature parameter(3)
Magnetic field parameter declines both fluid velocity and temperature(4)
Fluid temperature is increasing function of both Brownian motion and thermophoresis parameters(5)
Nanoparticle concentration shows diminishing variations towards Brownian motion parameter while the opposite trend is noticed versus thermophoresis parameter(6)
Nanoparticle concentration remarks inciting values towards Lewis and Prandtl numbers(7)
The skin-friction coefficient shows increasing nature forα,β,γ, but opposite behavior is observed for λandMf(8)
Wall heat flux coefficient increases againstγ,Pr while it decreases for Nt,Nb,andk(9)
The local Sherwood number shows an inciting value viaLe,Pr,Nt,Nb,andδc



---

*Source: 5162423-2021-01-19.xml*