# The Investigation of MHD Williamson Nanofluid over Stretching Cylinder with the Effect of Activation Energy

Authors: Wubshet Ibrahim; Mekonnen Negera
Journal: Advances in Mathematical Physics
 (2020)

---

## Abstract

In this paper, we discussed the effect of activation energy on mixed convective heat and mass transfer of Williamson nanofluid with heat generation or absorption over a stretching cylinder. Dimensionless ordinary differential equations are obtained from the modeled PDEs by using appropriate transformations. Numerical results of the skin friction coefficient, Nusselt number, and Sherwood number for different parameters are computed. The effects of the physical parameter on temperature, velocity, and concentration have been discussed in detail. From the result, it is found that the dimensionless velocity decreases whereas temperature and concentration increase when the porous parameter is enhanced. The present result has been compared with published paper and found good agreement.


---

## Body

## 1. Introduction

The heat and mass transfer of boundary layer flow of non-Newtonian fluids over a stretching cylinder is of great interest for scientists, engineers, and researchers due to its wide applications. Some examples of its applications are the extrusion process, extraction of metals, annealing, thinning of copper wire, and pipe industry. Nowadays, in most industries, the significance of non-Newtonian fluids governs the Newtonian fluids. The rheological properties of non-Newtonian fluids cannot be illuminated by the classical Naiver-Stokes equations. Also, non-Newtonian fluid characteristics cannot be modeled by a single model. To overcome this difficulty, ample models have come into being. The rheological models that were projected were Williamson, Cross, Ellis, power law, Carreau fluid model, etc. Typical of a non-Newtonian fluid model with shear retreating property is the Williamson fluid model and was first expected by Williamson [1]. Some recent investigations concerning the flow and heat and mass transfer of MHD flow of Williamson fluid can be mentioned through investigations [2–6].Heat and mass transfer in non-Newtonian fluid flow through porous medium in engineering has extensive application, such as ventilation procedure, oil production, solar collection, cooling of nuclear reactors, and electronic cooling. Accordingly, Maripala and Kishan [7] investigated the effect of thermal radiation and chemical reaction on time-dependent MHD flow and heat transfer of nanofluid over a permeable shrinking sheet. Their result indicates that with an increase of the suction parameter, the temperature profiles decrease whereas the concentration profiles upsurge. Moreover, Kairi and RamReddy [8], Vijaya et al. [9], Ambreen et al. [10], Aurangzaib et al. [11], Bal Reddy et al. [12], Al-Mamun et al. [13], and Ahmed et al. [14] studied mixed convective heat and mass transfer of MHD nanofluids embedded in porous medium over stretching surfaces. From their result, it can be seen that temperature distribution and thermal boundary layer reduce with an increase of dimensionless thermal free convection parameter, dimensionless mass free convection parameter, and Prandtl number.The miracle of heat generation/absorption in working fluid plays an important role in several engineering and industries. The influence of heat generation may increase the temperature distribution in moving fluids and consequently impacts the heat transfer rate. Some of its application in many industries is the geothermal system, thermal absorption, removal of heat from nuclear wreckage, food storage, microelectronics manufacturing, and cooling of nuclear reactors. Consequently, heat generation or absorption in MHD flow of Casson fluid over a stretching wedge with viscous dissipation and Newtonian heating was investigated by Ullah et al. [15]. Moreover, Khan et al. [16] and Padmavathi and Anitha Kumari [17] examined the effect of heat generation or absorption on mixed heat and mass transfer nanofluids embedded over stretching surfaces. Some researchers have studied the impact of heat source/sink on MHD non-Newtonian nanofluids over stretching surfaces [18–23].Activation energy is the smallest amount of energy needed by chemical reactants to endure a chemical reaction. The influence of activation energy on convective heat and mass transfer in the region of boundary layers was initially inspected by Bestman [24]. Later, many researchers have studied the impact of activation energy on heat and mass transfer of boundary layer flow of the fluids. Among those researchers, Awad et al. [25], Dhlamini et al. [26], Anuradha and Sasikala [27], and Hamid et al. [28] scrutinized the influence of activation energy on heat and mass transfer in unsteady fluid flow under different geometry. Moreover, Huang [29] and Mustafa [30] studied the effect of activation energy on MHD boundary layer flow of nanofluids past a vertical surface and permeable horizontal cylinder. Further investigations on the effect of activation energy on non-Newtonian fluid under different surfaces are stated in [31–34]. From their plots, it can be seen that as the values of the activation energy parameter upsurges, the concentration nanoparticles increase.In view of all the above-revealed studies, it is decided that the effect of activation energy on mixed convective heat and mass transfer of Williamson nanofluid over a stretching cylinder embedded in a porous medium with heat generation/absorption is not examined yet. Thus, to fill this gap, we aim to explore the effect of activation energy on mixed convective heat and mass transfer of Williamson nanofluid over a stretching cylinder embedded in a porous medium with heat generation/absorption. For information, our work is innovative in terms of the proposed fluids and incorporated parameters in the boundary layer flow. A numerical solution is obtained by using the Runge-Kutta method in conjugation with the shooting technique. Numerical results of the Nusselt number, skin friction coefficient, and Sherwood number for different values of the physical parameter are computed in tables. The effects of different physical parameters on dimensionless velocity, temperature, and concentration are presented in graphs.


## 2. Mathematical Formulation

Two-dimensional time-independent Williamson nanofluid over a stretching cylinder in the presence of a boundary layer slip has been considered. A porous cylinder is chosen as shown in Figure1. The magnetic field of constant strength B0 is applied normal to the flow along the r-axis in a radial direction with the assumption of a small Reynolds number so that the induced magnetic field is neglected. The fluid is infinite in magnitude of positive x-direction. The mixed convective mass and heat transfer phenomenon are included in the existence of heat generation/absorption. Tw and Cw are surface temperature and concentration, respectively. The ambient fluid temperature and concentration are denoted by T∞ and C∞, respectively. Further, u and v are the velocity components in the direction of x and r, respectively. All the physical properties of the fluid are considered constant.Figure 1
Physical model of the problem.Using the assumptions given above, the governing equations of continuity, momentum, energy, and mass are specified as follows:(1)∂∂xru+∂∂yrv=0,(2)u∂u∂x+v∂u∂r=ν1r∂u∂r+∂2u∂r2+νΓ2r∂u∂r2+2Γ∂u∂r∂2u∂r2−νk′u−σ1B02ρu+gβtT−T∞+gβcC−C∞=0,(3)u∂T∂x+v∂T∂r=kρcp1r∂T∂r+∂2T∂r2+τDB∂C∂r∂T∂r+DTT∞∂T∂r2+QρcpT−T∞,(4)u∂C∂x+v∂C∂r=DB1r∂C∂r+∂2C∂r2+DTT∞1r∂T∂r+∂2T∂r2−kr2C−C∞TT∞eexpEak1T.The boundary conditions are(5)u=u0xL+b0ν∂u∂r,v=−vw,T=Tw,C=Cwatr=R,u⟶u∞,T⟶T∞,C⟶C∞,asr⟶∞,where Tw=T∞+T0x/LβandCw=C∞+C0x/Lς are surface temperature and surface concentration, respectively; u0,T0,andC0 are reference velocity, temperature, and concentration, respectively; β is temperature exponent; and ς is concentration exponent signifying the change of amount of solute in x-direction.The partial differential equations given in equations (2)–(4) are converted to a system of ordinary differential equations by inspiring the following similarity transformations:
(6)η=r2−R22Ru0xν1/2,ψRfη,θη=T−T∞Tw−T∞,ϕη=C−C∞Cw−C∞,u=1r∂ψ∂r,v=−1r∂ψ∂x.Thus, the system of ordinary differential equations is(7)1+2ηAf′′′+2Af′′+321+2ηA1/2λf′′2+λ1+2ηA3/2f′′f′′′−f′2−d+Mf′+pθ+qϕ=0,(8)1+2ηAθ′′+PrNb1+2ηAθ′ϕ′+2Aθ′+Prf′θ′+PrNt1+2ηAθ′2+PrSθ−Prβf′θ=0,(9)1+2ηAϕ′′+Scfϕ′+2Aϕ′+2NtNbAθ′+NtNb1+2ηAθ′′−Scσ1+ζθeexp−E1+ζθ−Scςf′ϕ=0.With the boundary conditions(10)f0=δ,(11)f′0=1+γf′′0,(12)θ0=1,(13)ϕ0=1,(14)f′∞=0,(15)θ∞=0,(16)ϕ∞=0,where M=σ1B02L/ρu01/2 is the magnetic parameter, A=νL/u0R2 is the curvature parameter, γ=B0νu0/L is the slip parameter, λ=Γ2u03/νL is the Weissenberg number (Williamson parameter), p=gβtTw−T∞/Rexν is the local temperature buoyancy parameter, q=gβcCw−C∞/Rexν is the concentration buoyancy parameter, d=νL/u0k is the porosity parameter, Nb=τDBCw−C∞/v is the Brownian motion parameter, Nb=τDTTw−T∞/νT∞ is the thermophoresis parameter, S=QL/ρcu0 is the heat source/sink parameter, Pr=μcp/k is the Prandtl number, Sc=ν/DB is the Schmidt number, σ=kr2/c is the ratio rate parameter, E=Ea/k1T∞ is the activation energy parameter, and ζ=Tw−T∞/T∞ is the temperature difference parameter.Local skin friction coefficientCf, local Nusselt number Nux, and local Sherwood numberShx are defined as follows:
(17)Cf=2τwρuw2,Nux=xpwkTw−T∞,Shx=xqwDBCw−C∞,where
(18)τw=μ∂u∂r+Γ2∂u∂r2r=R,pw=−k∂T∂rr=R,qw=−DB∂C∂rr=R.The nondimensional forms of the local Nusselt number, local skin friction coefficient, and local Sherwood number are(19)CfRex2=f′′0+λ2f′′20,NuxRex−1/2=−θ0,ShxRex−1/2=−ϕ0,where Rex=u02x/v is the local Reynolds number.


## 3. Solution Methodology

The coupled nonlinear ordinary differential equations (7)–(9) subjected to the boundary conditions (10) and (14) are solved numerically using the Runge-Kutta method with shooting technique. For the minor change of initial guesses u1,u2,andu3, the Newton-Raphson method is applied subject to the toleranceε=10−5. On the origin of a number of computational experiments, there is no important change in the results afterη=5, so we are considering 0,5 as the domain of the problem in terms of0,∞. To solve this problem by using this method, equations (7)–(9) are converted to a system of first order ordinary differential equations. The system of first order ordinary differential equations is defined:
(20)f=y1,f′=y2,f″=y3,f‴=y′3,θ=y4,θ′=y5,θ″=y′5,ϕ=y6,ϕ′=y7,ϕ″=y7′.Thus, the systems of first order concurrent ODEs are(21)y′1=y2,y′2=y3,y′3=−2Ay3−3/21+2ηA1/2Aλy32−y1y3+y22+d+My2−py4−qy61+2ηA+λ1+2ηA3/2y3,y′4=y5,y′5=−PrNb1+2ηAy5y7−2Ay5−Pry1y5−PrNt1+2ηAy52+PrSy4+Prβy2y41+2ηA,y′6=y7,y′7=−Scy1y7+PrNt1+2ηAy5y7+2Ay7+4Nt/NbAy5−2Nt/NbPry1y5+2PrNt2/Nb1+2ηAy52+2Nt/NbPrSy4−2Nt/NbPrβy2y4+Scσ1+ζy4eexp−E/1+ζy4y6+Scςy2y61+2ηA.Here, the prime represents the derivative with respect toη, and the transformed boundary conditions are
(22)y10=δ,y20=1+γu1,y30=u1,y40=1,y50=u2,y6=1,y7=u3.


## 4. Result and Discussion

Figures2(a)–2(c) are plotted to show the distribution of velocity, temperature, and concentration for different values of the suction δ parameter, respectively. From these figures, it can be seen that when the values of the suction parameter increases, the velocity, temperature, and concentration boundary layer thickness are reduced. This is due to the fact that suction or blowing is the method of controlling the boundary layer. The effect of suction consists in the removal of decelerated fluid particles from the boundary layer before they are given a chance to cause separation. Figures 3(a) and 3(b) reveal the profiles of velocity and temperature for different values of injection −δ. The graph shows that with an increase of the injection parameter, velocity and temperature upsurge.Figure 2
The graph of velocity, temperature, and concentration for different values of suction (δ).
(a)(b)(c)Figure 3
The graph of velocity and temperature for different values of injection (−δ).
(a)(b)For different values of slip parameterγ, the variation of velocity and temperature is plotted in Figures 4(a) and 4(b). Figure 4(a) shows that when the value of the slip parameter raises, the velocity profiles reduce. This is because of the fact that when the slip parameter increases, slip velocity increases, and consequently, fluid velocity declines because under the slip condition, the pulling of the stretching wall can only be partly conveyed to the fluid. Figure 4(b) displayed that the temperature profiles raise with an increase of the slip parameter. Figures 5(a), 5(b), and 5(c) are plotted to show the velocity, temperature, and concentration profiles, respectively, for various values of the magnetic field. It is noticed in Figure 5(a).Figure 4
The graph of velocity and temperature for different values of slip parameterγ.
(a)(b)Figure 5
The graph of velocity, temperature, and concentration for different values of magnetic field (M).
(a)(b)(c)The graph of velocity, temperature, and concentration is plotted in Figures6(a), 6(b), and 6(c), respectively, for different values of the porosity parameter. The porous medium causes higher restriction to the fluid flow, which in turn decelerates its motion. Consequently, from Figure 6(a), we observed that with an increase of the porosity parameter, the velocity profiles decrease. But from Figures 6(b) and 6(c), it can be seen that with an increase of the permeability parameter, both temperature and concentration profiles are increased. The distribution of velocity and temperature for different values of curvature parameter (A) is plotted in Figures 7(a) and 7(b), respectively. Figure 7(a) displays that the velocity profiles are increased with an increase of the curvature parameter whereas Figure 7(b) speculates that when the values of the curvature parameter upsurge, the temperature boundary layer thickness declines. This is due to the fact that when the curvature parameter increases, the radius of the cylinder reduces which deduces that the contact area of the cylinder with the fluid reduces. Hence, less resistance is offered by the surface to the fluid motion, and as a result, the fluid velocity improved.Figure 6
The graph of velocity, temperature, and concentration for different values of porous parameter (d).
(a)(b)(c)Figure 7
The graph of velocity and temperature for different values of curvature parameter (A).
(a)(b)For different values of temperature buoyancy parameter (p), the graph of velocity, temperature, and concentration is presented in Figures 8(a)–8(c). From these figures, we observed that velocity profiles increase with an increase of the temperature buoyancy parameter, but temperature and concentration field decrease with an increase of the temperature buoyancy parameter. Figures 9(a)–9(c) show the graph of temperature for different values of heat source/sink parameter (S) and temperature exponent β. It is seen that the temperature field increases with an increase of the heat source, but an opposite trend is happening for large values of heat sink and temperature exponent β. Figures 10(a) and 10(b) display the graph of velocity, temperature, and concentration, respectively, for various values of concentration buoyancy parameter (q). The figures illustrate that the velocity field increases whereas temperature and concentration fields decline with an increase of the concentration buoyancy parameter.Figure 8
The graph of velocity, temperature, and concentration for different values of temperature buoyancy parameter (p).
(a)(b)(c)Figure 9
The graph of temperature for different values of heat source/sink parameter (S) and temperature exponent β.
(a)(b)(c)Figure 10
The graph of velocity, temperature, and concentration for different values of concentration buoyancy parameter (q).
(a)(b)(c)The effect of activation energy (E) and reaction rate parameter σ on the concentration field is plotted in Figures 11(a) and 11(b). From these figures, it is seen that when activation energy upsurges, the concentration field rises, but it diminishes with an increase of the reaction rate parameter. This is due to the fact that large activation energy and low temperature cause a smaller reaction rate constant and thus slow down the chemical reaction. So, the concentration of the solute increases. Figures 12(a), 12(b), and 12(c) execute the impact of the Williamson parameter on velocity, temperature, and concentration fields, respectively. Figure 12(a) tells us that with an increase of the Williamson parameter, the velocity profiles decrease. As we know that the Williamson parameter is the ratio of relaxation time to retardation time, hence, as the values of the retardation time decrease, the Williamson parameter increases which displays that there is a reduction in the velocity field, and in addition, boundary layer thickness decreases. On behalf of Figures 12(b) and 12(c), it is noted that the temperature and concentration fields upsurge with an enhancement of the Williamson parameter.Figure 11
The graph of concentration for different values of activation energy (E) and reaction rate parameter σ.
(a)(b)Figure 12
The graph of velocity, temperature, and concentration for different values of Williamson parameterλ.
(a)(b)(c)The variation of concentration profiles for various values of fitted rate constant (e) and concentration exponent ς is shown in Figures 13(a) and 13(b). It is noticed that the distribution of concentration profiles is decreased with an increase in both fitted rate constant and concentration exponent parameters. This is because when the values of the fitted rate constant increased, the factor σ1+ζθeexp−E/1+ζθ is augmented. This finally favors the destructive chemical reaction which causes concentration gradient rise. So, the diminution in the concentration field is accompanied by a higher concentration gradient at the wall. The impact of Brownian motion (Nb), thermophoresis parameter (Nt), and Prandtl number (Pr) on temperature is plotted in Figures 14(a)–14(c). Figure 14(a) indicates that when the Brownian motion Nbparameter increases, the movement of nanoparticles from the hot surface to the cold surface happens and ambient fluid occurred. Due to this, the temperature and thermal boundary layer thickness rise. From Figure 14(b), it is observed that the temperature field grows with an increase of the thermophoresis parameter. A phenomenon in which small particles are pulled away from the hot surface to the cold one is called thermophoresis. So, when the surface is heated, the large number of nanoparticles is moved away which raises the temperature of the fluid. Therefore, the temperature of fluid increases. Figure 14(c) displays the effect of the Prandtl number (Pr) on temperature. From this, it can be seen that when the values of Pr rise, the temperature profiles decline. This is because fluid with higher Pr has relatively low thermal conductivity, which results in heat conduction and thereby thermal boundary layer thickness and temperature drop. The effect of Brownian motion (Nb), thermophoresis parameterNt, Schmidt number Sc, and temperature difference parameterζ on concentration profiles is plotted in Figures 15(a)–15(d). Figure 15(a) indicates that when Brownian motionNb upsurges, volume fraction of nanoparticles within the boundary layer upsurges. It is interesting to note that Brownian motion of the nanoparticles at the molecular and nanoscale levels is a key mechanism in governing their thermal behavior. In the nanofluid system, due to the size of the nanoparticles, Brownian motion affects the heat transfer properties. As the particle size scale approaches to the nanometer scale, Brownian motion on the surrounding liquids plays an important role in heat transfer. Figure 15(b) shows that with an increase of the thermophoresis parameter, concentration profiles rise. From Figure 15(c), it is observed that the concentration profiles decrease with an increase of the Schmidt number. This is due to the fact that weaker molecular diffusivity appeared for a larger Schmidt number and stronger molecular diffusivity corresponds to a lower Schmidt number. Figure 15(d) forecasts that with an increase in the temperature difference parameter, the concentration profiles decline.Figure 13
The graph of concentration for different values of fitted rate constant (e) and concentration exponent ς.
(a)(b)Figure 14
The graph of temperature for different values of Brownian motion (Nb), thermophoresis parameter (Nt), and Prandtl number (Pr).
(a)(b)(c)Figure 15
The graph of concentration for different values of Brownian motion (Nb), thermophoresis parameter (Nt), Schmidt numberSc, and temperature difference parameter (ζ); the velocity declines with an increase of the magnetic field parameter. From (b, c), it can be seen that the temperature and concentration upsurge with an increase of the magnetic field parameter (M). The physical meaning of this behavior is an upsurge in the values of the magnetic field parameter develops the force opposite to the flow, which is known as Lorentz force. This body force reduced the boundary layer flow and congeals the velocity boundary layers. Furthermore, this force creates a resistance force that opposes the fluid motion. Therefore, for the higher magnetic field, heat is generated leading to a boost in the temperature and thermal boundary layer thickness.
(a)(b)(c)(d)Table1 shows the comparison of the skin friction coefficient with the available published results of Malik et al. [2] and found that it is an excellent agreement. Moreover, the comparisons of the values of the Nusselt number and Sherwood number with the other published result of Hayat et al. [33] are shown in Table 2. Here, from the table, it can be seen that there is an excellent agreement.Table 1
Comparison of values of skin friction coefficientf″0 for various values of λ with earlier scholar results.
ΛMalik et al. [2]f″0Present resultsf″00.1−0.9776−0.97790.2−0.9118−0.91280.3−0.8413−0.8417Table 2
Comparison of values of the Nusselt number−θ′0 and Sherwood number −ϕ′0 for various values of Pr and E with earlier scholar results.
PrEHayat et al. [33]−θ′0Present resultHayat et al. [33]−ϕ′0Present result110.519030.5191320.851090.8512151.515531.5155871.845401.84547102.260332.2603910.950640.9506920.753330.7533840.580150.5802060.535280.5353080.523680.52371Straightforwardly, the numerical values of the skin friction coefficient, Nusselt number, and Sherwood number appear in Table3 for different values of suction parameter δ, permeability parameter (d), curvature parameter (A), temperature buoyancy parameter (p), concentration buoyancy parameter (q), and slip parameter γ. It is noticed that the magnitude of the skin friction coefficient increases with an increase of suction parameter δ, permeability parameter (d), and curvature parameter (A) whereas with an upsurge of temperature buoyancy parameter (p), concentration buoyancy parameter (q), and slip parameter γ, the skin friction coefficient magnitude was reduced, but it remains constant with an increase of activation energy (E) and reaction rate parameter σ. Moreover, from this table, it can be seen that when the values of suction parameter δ, temperature buoyancy parameter (p), concentration buoyancy parameter (q), activation energy (E), and curvature parameter (A) increase, the magnitude of the Nusselt number increases whereas with an increase of permeability parameter (d), reaction rate parameter σ, and slip parameter γ, there is a magnitude of Nusselt number diminution. Furthermore, from this table, it can be seen that the Sherwood number declines when the values of permeability parameter (d), activation energy (E), and slip parameterγ rise, but an opposite trend is happening on the Sherwood number as the values of suction parameter δ, temperature buoyancy parameter (p), concentration buoyancy parameter (q), and reaction rate parameter σ.Table 3
Numerical results of the skin friction coefficient, Nusselt number, and Sherwood number for different values ofδ,d,p,q,γ,E,R,andσ when Nt=0.1,Nb=0.1,M=0.3,λ=0.1,e=0.1,ς=1.0,ζ=0.1,Sc=3,Pr=2,β=1,andS=0.1.
δdpQEAσγ−f″0−θ′0−ϕ′00.00.00.10.11.00.00.10.01.59441.38861.88930.11.71161.48232.05400.21.84821.58082.22770.32.01431.68292.40950.01.66841.60152.25300.21.84821.58052.22770.42.05001.56022.20290.62.30451.54022.17840.11.76281.58952.23870.31.61211.60602.25880.51.48021.62092.27711.01.20261.65342.34670.11.69921.59142.24610.31.58261.60662.26000.51.47701.61712.27311.51.24691.64072.30241.01.84821.55091.19092.01.84821.55341.17814.01.84821.55381.17118.01.84821.55391.16990.01.54911.52732.19910.151.84821.58052.22770.21.99301.59672.23520.252.20871.61172.24130.11.84821.55271.19090.31.84821.55041.23150.71.84821.54631.30801.51.84821.53931.44570.01.84821.58052.22770.21.15061.44132.03110.30.99581.39401.96430.40.88321.35441.9083


## 5. Conclusion

This article deals with numerical solution of MHD non-Newtonian nanofluid over a stretching cylinder embedded in a porous medium with the effect of activation energy and heat generation/absorption. The impact of various dimensionless physical parameters on velocity, temperature, and concentration profiles as well as skin coefficient frictionCf, heat transfer rate Nux, and mass transfer rate parameter Shx is discussed in graphs and tables. The summarized conclusion of numerical study is as follows:
(i)
The suctionδ parameter interaction reduces the velocity, temperature, and concentration profiles of the flow(ii)
Curvature parameter (A) improves velocity profiles, but it reduces temperature profiles(iii)
With an increase of porous parameter (d), the profiles of velocity decline whereas temperature and concentration profiles are improved(iv)
Raise in temperature and concentration buoyancy parameter enhances velocity whereas temperature and concentration are reduced(v)
Increase in activation energy (E) enhances concentration profiles(vi)
With an increase of reaction rate parameterσ, the magnitude of the skin friction coefficient remains constant but the magnitude of Nusselt number diminution and Sherwood number rises(vii)
Sherwood number declines when the value of slip parameterγ rises.



---

*Source: 9523630-2020-06-25.xml*