Effect of Activation Energy on Magnetized Couple Stress Fluid over an Inclined Stretching Permeable Cylinder

Document Type : Full Length Research Article

Authors

Department of Mathematics, University of Rajasthan, Jaipur, 302004, India

Abstract

Activation energy is of considerable significance in diverse applications such as chemical kinetics, catalyst development, enzymes, semiconductors, and systems sensitive to temperature, such as chemical reactors and engines. The objective of this research is to investigate the influence of activation energy on a magnetized couple stress fluid over an inclined stretching permeable cylinder in a non-Darcy porous medium. The effects of cross-diffusion and stratified mixed convection are also considered in fluid model. The boundary layer equations, which describe the flow, have been converted into dimensionless form through suitable transformable variables. Subsequently, these transformed equations are solved using fourth order Runge-Kutta mechanism along with the shooting technique. The outcomes comprise visual depictions and comprehensive explanations demonstrating the influence of relevant variables on thermal, concentration, and velocity fields. Observations reveal that the concentration profile is directly influenced by the Forchheimer number and activation energy parameter, whereas both temperature and concentration fields decrease with elevated thermal and solutal stratification parameters. Additionally, numerical outcomes for the skin-friction coefficient, Nusselt number, and Sherwood number are presented in tabular form.

Keywords

Main Subjects


1. Introduction

Stokes [1] introduced couple stresses in viscous fluids, envisioning size-dependent effects that are absent in the polar case as the primary characteristic of couple stresses. Examples of fluids exhibiting couple stresses include lubricants with trace amounts of polymer additives, synthetic blood fluids, microfibers, and electro-rheological fluids. These materials find diverse industrial applications, such as in the extrusion of polymer fluids, the solidification of animal blood, and the study of liquid crystals. Rani et al. [2] investigated the flow of couple stress fluid around an infinite vertical cylinder, while Gajjela and Garvandha [3] examined the impact of cross diffusion and chemical reactions on couple stress flow through a stretching cylinder with convective boundary conditions.

Ibrahim and Gadisa [4] studied the behavior of a couple stress fluid around a stretching cylinder under conditions of double stratified mixed convection. Palaiah et al. [5] examined the impact of thermal radiation on magnetized couple stress fluid flowing over a vertical cylinder. Awais and Salahuddin [6] determined that the concentration profile of a couple stress fluid passing a paraboloid is enhanced with the increase in activation energy.

Activation energy holds crucial significance in various applications, including chemical kinetics, reactor cooling, emulsification of oil and water, catalyst and enzyme design, semiconductor processes, and temperature-sensitive systems such as chemical reactors and engines. Coined by Arrhenius [7], activation energy represents the minimum amount of energy required for chemical reactants to undergo a chemical reaction. Bestman [8] delved into boundary layer flow accompanied by activation energy and chemical reactions. Abbas et al. [9] observed that the concentration profile of Casson fluid over a stretching sheet increases with the activation energy parameter. Waqas et al. [10] examined the impact of activation energy and magnetic fields on couple stress nanofluid. Ibrahim and Negera [11] investigated the Williamson nanofluid model with activation energy effects around a stretching cylinder. Mustafa et al. [12] analyzed buoyancy-driven MHD nanofluid considering the influence of chemical reaction and activation energy on a vertical surface. Hamid et al. [13] explored Williamson nanofluid under the influence of activation energy. Awad et al. [14] studied activation energy and chemical reactions in unsteady rotating fluid flow. Maleque [15] investigated the effects of activation energy and chemical reactions in the region of boundary layers.

Darcy's law asserts that velocity is directly proportional to the pressure gradient, applicable primarily in laminar flow conditions with low velocity and porosity. However, Darcy's law becomes inadequate when inertial effects become significant. To address high-velocity flow, a non-Darcy term known as the Forchheimer drag parameter is introduced. Saeed et al. [16] explored the Darcy-Forchheimer flow of MHD hybrid nanofluid with heat transfer around a porous cylinder. Nagaraja et al. [17] investigated non-Newtonian flow influenced by a non-Darcy medium around a circular cylinder. Saeed et al. [18] analyzed the behavior of a couple stress hybrid nanofluid in a non-Darcy porous medium with variable fluid properties. Ganesh et al. [19] examined Darcy-Forchheimer flow of MHD nanofluid with second-order slip conditions over a stretching sheet. Pal et al. [20] studied the impact of variable viscosity on the convective diffusion of species in a non-Darcy porous medium. Hayat et al. [21] examined Darcy-Forchheimer flow coupled with the Cattaneo-Christov model and variable thermal conductivity.

The application of heat transfer over a stretching cylinder has widespread applications in wire drawing, fire production, and hot rolling. Hayat et al. [22] analyzed the impact of thermal radiation and a non-uniform heat source on couple stress fluid in a thermally stratified medium over a stretching cylinder. Jain and Choudhary [23] studied the effects of suction and injection over a permeable horizontal cylinder in a porous medium. Ali A. et al. [24] investigated the effects of non-uniform heat flux and thermal radiation over a stretching cylinder with MHD hybrid nanofluid. Song et al. [25] analyzed Williamson nanofluid in the presence of a non-uniform heat source.

Stratification of the medium involves the deposition or formation of layers due to temperature or concentration differences in different fluids. This phenomenon is characteristic of all fluid bodies surrounded by heated sidewalls. Parmar and Jain [26] studied heat transfer in Casson fluid flow through a cylinder. Rehman et al. [27] investigated double stratification for Williamson fluid flow in the presence of chemical reactions and mixed convection. Rehman et al. [28] explored the effects of mixed convection and stratification in the stagnation region along an inclined cylinder on Williamson fluid flow. Ramzan et al. [29] analyzed the thermal and solutal stratification effects on Jeffrey fluid around an inclined stretching cylinder with heat generation/absorption and thermal radiation. Sohut et al. [30] investigated the effect of double stratification along a stretching cylinder in the presence of chemical reactions. Rehman et al. [31] studied Eyring-Powell fluid flow in a dual stratified mixed convective medium with the effect of heat generation/absorption. Abbasi et al. [32] studied Maxwell nanofluid flow with double stratified mixed convection.

After observing the above-cited literature,

  • It can be perceived that no research has been done to study the impact of activation energy on magnetized couple stress fluid over an inclined stretching permeable cylinder.
  • The objective of this research is to analyse the impact of activation energy on magnetized couple stress fluid with cross-diffusion and double stratification.
  • The boundary layer equations, which describe the flow, are converted into dimensionless form through suitable transformable variables. Subsequently, these transformed equations are solved using the fourth order Runge-Kutta mechanism along with shooting technique.
  • Graphical representations are used to analyse results related to different parameters for velocity, temperature and concentration profiles. Additionally, numerical outcomes for the skin-friction coefficient, Nusselt number and Sherwood number are presented in tabular form.

2. Mathematical Formulation

We have considered two-dimensional, magneto double-stratified couple stress fluid flows past an inclined permeable stretching cylinder in a non-Darcy porous medium. The impacts of Activation energy, nonlinear thermal radiation, cross-diffusion and of non-uniform heat source are also taken. The x-axis is taken along the axis of cylinder and r-axis is normal to it. The cylinder is longitudinally stretching with velocity U(x)=U0Lx, here U0 is reference velocity and L is characteristic length. Inclination of the cylinder is denoted by α. A uniform magnetic field B0 is also applied along r-direction. The temperature and concentration are maintained constant, symbolized by Tw and Cw respectively at surface and far from the boundary they are symbolized by T and C respectively.

Fig. 1. Coordinate system and schematic diagram of the problem

Under above suppositions the equations of continuity, momentum, energy and concentration that govern the flow are written as:

Continuity equation

(ru)x+(rv)r=0(1)

Momentum equation

uux+vur=υ(2ur2+1rur)-υ (4ur4+2r3ur3-1r22ur2+-1r3ur)+(gβ1(T-T)+gβ2(T-T)2)cos α-υk ′u-F0u2-𝜎B02𝜌u(2)

Energy equation

uTx+vTr=k1ρCp(1rTr+2Tr2)+DKTCsCp(2Cr2+1rCr)-1ρCp1r(rqr)r+q″1ρCp(3)

Concentration equation

uCx+vCr=D(1rCr+2Cr2)DKTT(2Tr2+Tr)-Kr2(C-C)(TT)eexp(-EaK2T)(4)

in which T is the fluid temperature, C is the fluid concentration, g is gravity, β1 is thermal expansion coefficient, β2 is the coefficient of volumetric expansion, k′ is the permeability of porous medium, σ is the electrical conductivity, qr=-16σ*T33k*Tr is the radiative heat flux, σ* is the Stefan-Boltzmann constant, k* is the mean absorption coefficient, q″ is the rate of internal heat generations/absorption coefficient mentioned as:

q‴=(kuwxν)[X*(Tw-T)+Y*(T-T)](5)

where X* and Y* are space/time dependent heat source/ sink parameter (the case X*>0 and Y*>0 holds for internal heat generation, while X*<0 and Y*<0 for internal heat absorption), k is the thermal conductivity of the fluid, Ea is activation energy, K2 is the Boltzmann constant.

The appropriate boundary conditions as per formulation are prescribed as:

At r=R,u=U(x)=U0Lx,v=v0,TT(x)=T0+bxL,X(x,r)=Cw(x)=C0+dxL.(6)as r,u=0,TT(x)=T0+cxL,CC(x)=C0+exL

here U0 is reference velocity, L is characteristic length, and b, c, d and e are positive constant.

The stream function ψ that satisfies the continuity equation is define as

u=1r(ψr)andv=-1r(ψx)(7)

The dimensionless similarity transformations are:

u=U0Lxf′(η),η=r2-R22R(U0υL)12v=Rr(U0υL)12f(η),ψ=(U0υx2L)12Rf(η),(8)θ(η)=T-TTw-T0,ϕ(η)=C-CCw-C0

Using these non-dimensional variables equations (1)-(4) are reduced into non-linear ordinary differential equations which are written as follows:

(1+2γη)f‴+f″+ff″-(1+F)f′2-KRe[(1+2γη)2f‴″+(1+2γη)f‴′+8γ2f‴]+(θ+γ1θ2)λ cos α-(d+M)f′=0(9)

Pr(f′θ-fθ′)=(1+43Rd)[θ′+(1+2γη)θ″]+PrDf[ϕ′+(1+2γη)ϕ″]+[X(1-δ1)f′+](10)

f′ϕ-fϕ′=1Sc(ϕ′+(1+2γη)ϕ″)+Sr(θ′+(1+2γη)θ″)-σϕ(1+ζθ)eexp(-E1+ζθ)(11)

Reduced boundary conditions are:

f(0)=-S,f′(0)=1,θ(0)=1-δ1,ϕ(0)=1-δ2f′=f″=f‴=θ=ϕ=0asη(12)

Where, F=xF0 Forchheimer number; K=υ′υR2 Couple stress parameter; Re=υ0R2υL Reynolds number; d=υLυ0k′ porosity parameter; γ=υLυ0R212 curvature parameter; γ1=β2(Tw-T0)β1 nonlinear temperature convection coefficient; λ=g(Tw-T0)β1L2U02x combined convection parameter; M=σB02LρU0 magnetic field parameter; Pr=μCpK Prandtl number; fw=ν0LU0ν suction/injection parameter Rd=4σsT3kk* thermal radiation parameter; Df=DKTCsCpυb Dufour number; X = space/time dependent heat source parameter; Y = space/time dependent heat sink parameter; δ1=cb thermal stratification parameter; Sc=υD Schmidt number; Sr=DKTυTbd Soret number; σ=Kr2dxL ratio rate parameter; ζ=Tw-T0T temperature difference parameter;

E=EaK2T activation energy parameter; and the powers denote differentiation with respect to η.

The skin-friction coefficient Cf, Nusselt number Nu and Sherwood number Sh that characterize this study are given by:

Skin-friction coefficient

Cf=2τwρUw2whereτwμ(ur)r=R(13)

Nusselt number

Nux=-xqwk(Tw-T0)whereqw=-k(Tr)r=R(14)

Sherwood number

Shx=-xjwk(Cw-C0)wherejw=-DB(Cr)r=R(15)

Using non-dimensional parameters equations (13-15) are written as:

CjfRex12=-f″(0),(16)

NuxRex-12=-θ′(0),(17)

ShxRex-12=-ϕ′(0)(18)

3. Solution of the Problem

The coupled ODEs (9)-(11) together with the transformed boundary conditions (12) are unravelled by fourth order Runge-Kutta mechanism together with Shooting technique, using a MATLAB boundary value problem solver.

In the first step, the non-linear ODEs (9-11) are altered to first order differential equations.

Let

f=f1,f′=f2,f″=f3,f‴=f4,f‴′=f5,θ=f6,θ′=f7,ϕ=f8,ϕ′=f9

For these assumptions equations (9)-(11) will take the following form:

f(v)=K-1Re-1(1+2γη)-2[(1+2γη)f4+2γf3+f1f3-(1+F)f22-KRe{(1+2γη)f5+2f4}+λcosα(f6+γ1f62)f2(M+d)](19)

θ″=(1+2γη)-1(1+43Rd)-1×[Pr(f2f6-f1f7)-(1+43Rd)f7-PrDf{2γf9+(1+2γη)ϕ″}-{X(1+δ1)f1+Yf6}](20)

ϕ″=(1+2γη)-1[Sc{f2f8-f1f9Sr(2γf7+(1+2γη)θ″)+σf8(1+ζf6)eexp(11+ζf6)}-2γf9](21)

The reduced conditions (12) are:

f1(0)=0,f2(0)=1,f6(0)=1-δ1,f8(0)=1-δ2andf2()=f3()=f4()=f6()=f8()0(22)

We guess the values of f(3), f(7) and f(9), these are not given at initial conditions. Using step size 0.01 the process is continued, until we get the minimum error term.

4. Validation of the Study

Table 1 illustrates a comparison between the current findings and prior published studies conducted by Gajjela and Garvandha [3] and Asad et al. [18]. The observed agreement between the present results and the earlier ones suggests a strong consistency. Furthermore, these findings serve to validate the accuracy of the current results.

γ -f’’(0) Pr -θ’(0)
Asad [18] Gajjela [3] Present Asad [18] Gajjela [3] Present
0.0 1.6114 1.61145 1.61157 00.7 0.1902 0.19025 0.1905
00.05 1.6214 1.62142 1.62176 1 0.2035 0.20356 0.2042
00.1 1.6307 1.63071 1.63066 11.3 0.2127 0.21274 0.2133
00.15 1.6392 1.63915 1.6412 11.5 0.2174 0.21744 0.2282
Table 1. Comparison for the f^'' (0) and θ^' (0)

5. Results and Discussion

To analyse the physical aspects of the considered problem, a comprehensive study have been done. The graphs for fluid velocity, concentration and temperature for the various values of involved parameters such as couple stress parameter, curvature parameter, Reynolds number, Forchheimer number, non-linear temperature convection coefficient, combined convection parameter, porosity parameter, thermal stratification parameter, Schmidt number and activation energy parameter are shown in Fig. 2 to 25. The numerical outcomes for local skin friction, Nusselt and Sherwood number are demonstrated graphically and in tabular form also (Table 1-2) with a good agreement with published literature. Fig. (2-4) sketched the impact of couple stress parameter on velocity, concentration and temperature fields respectively. Couple stress parameter K acts as a retardant agent because couple stress parameter increases the couple stress viscosity that causes denser fluid. It is perceived that the velocity profile fall off with increasing couple stress parameter K, temperature profile increases when K increases and concentration profile is increasing function of couple stress parameter.

Fig. 2. Velocity distribution for K

Fig. 3. Temperature distribution for K.

Fig. 4. Concentration distribution for K.

Figure (5-7) reveals the influence of curvature parameter. By raising the curvature parameter γ, the radius of cylinder reduces which cause to reduce contact area of fluid with cylinder. So the thickness of boundary layer reduces, thermal and concentration boundary layers become more thicker with curvature parameter γ. Hence as γ grow, the velocity field decreases near the surface of cylinder and enhances far away from surface and the temperature and concentration profiles increases. Figure (8-10) illustrates the impact of magnetic field on various fields. Magnetic field initiates a resistive force called Lorentz force. By increasing the value of magnetic field parameter M, the Lorentz force increases and the thickness of boundary layer reduces, therefore the velocity profiles fall off with an increase in the magnetic field. Thermal and concentration boundary layers become thicker with enhancing Lorentz force because it is a body retarding force that transverses the motion of the fluid and heat is evolved. It is shown in figure (9) and (10) that the profiles of temperature and concentrations are enhances with increasing magnetic field.

Fig. 5. Velocity distribution for γ.

Fig. 6. Temperature distribution for γ.

Fig. 7. Concentration distribution for γ.

Fig. 8. Velocity distribution for M.

Fig. 9. Temperature distribution for M.

Fig. 10. Concentration distribution for M.

Figures (11-13) illustrate the impact of inclination angle. In Figure (11), it is observed that as the angle of inclination (α) increases, the velocity field decreases due to a decrease in gravitational force. When the angle of inclination rises, the effective gravitational force in the direction of flow reduces, leading to a weakening of the force propelling the fluid. As a result, the velocity field, which denotes the speed and direction of fluid flow, experiences a decline. Furthermore, the increased angle of inclination can result in a greater expanse of the cylinder's surface being subjected to fluid flow. This expanded surface area fosters increased interaction between the fluid and the cylinder, consequently enhancing the efficiency of heat and mass transfer. Hence, the thermal and concentration field Intensifies by enhancing inclination angle, depicted in figure (12-13).

Fig. 11. Velocity distribution for α.

Fig. 12. Temperature distribution for α.

Fig. 13. Concentration distribution for α

Figure (14) anticipates effect of activation energy E on concentration profile. Concentration profile rises with increasing activation energy parameter because reaction rate constant decreases due to higher activation energy and low temperature, this increases the concentration of solute. By increasing the Schmidt number Sc concentration profile decreases, it is shown in figure (15). Schmidt number is stated as the ratio of momentum diffusivity to mass diffusivity. Hence mass diffusivity reduces on increasing Schmidt number, which decreases the concentration profile.

Fig. 14. Concentration distribution for E.

Fig. 15. Concentration distribution for Sc.

The effect of mixed convection parameter λ on velocity, thermal and concentration profiles is demonstrated in figures (16-18). As mixed convection grows, the buoyancy force dominates the viscous force that translates the flow from laminar to turbulent which hastens the flow of the fluid. Figure (16) display the mixed convection has the effect of enhancing the velocity profile. Thickness of thermal and concentration boundary layer reduces with enhancing buoyancy force. Hence the temperature and concentration profiles decreases with increasing mixed convection parameter λ.

Fig. 16. Velocity distribution for λ.

Fig. 17. Temperature distribution for λ.

Fig. 18. Concentration distribution for λ.

Figure (19-20) represents that the temperature and concentration profile decreases with thermal stratification δ1 and solutant stratification parameter δ2 because when stratification increase then the convective potential between the fluid and surface of cylinder declines.

Fig. 19. Temperature distribution for δ_1.

Fig. 20. Concentration distribution for δ_2.

In figure (21-23) it is observed that the velocity profile decreases, temperature and concentration increases with an increasing Forchheimer number F. In the realm of permeability, the Forchheimer number delineates the balance between inertial and viscous forces governing fluid flow within porous materials, offering insight into the complex dynamics of fluid behavior. A high Forchheimer number signifies the dominance of inertial effects. Higher value of inertial coefficient produces resistance in the flow path. Thus, the thickness of velocity boundary layer reduces, thermal and concentration boundary layers enhances with Darcy Forchheimer number.

Fig. 21. Velocity distribution for F.

Fig. 22. Temperature distribution for F.

Fig. 23. Concentration distribution for F.

Figure (24-26) depicts the effect of Soret and Dufour number on boundary profiles. The values of Soret and Dufour are taken in such way that their product is constant. When Sr is raised, there will be a higher thermal diffusion. It is observed in figure (24) that the velocity enhances for high Dufour number and temperature and concentration profiles (25-26) decreases for high dufour number.

Fig. 24. Velocity distribution for Sr and Du.

Fig. 25. Temperature distribution for Sr and Du.

Fig. 26. Concentration distribution for Sr and Du.

Fig. 27. Temperature distribution for Rd.

Figure (27) shows the effect of thermal radiation parameter on thermal profile. An increase in the radiation parameter Rd means a decrease in the Rosseland radiation absorptivity κ^* and this increases the divergence of radiative heat flux. Hence it is observed that the temperature increases as thermal radiation parameter upsurge. Figure (28) displays declined in temperature profile when the value of prandtl number (Pr) rise. Higher prandtl number has low thermal conductivity which results in heat conduction and thus the temperature profile decreases with higher Prandtl number.

Fig. 28. Temperature distribution for Pr

Numerical results of skin friction coefficient f^'' (0), Nusselt number θ^' (0) and Sherwood number ϕ'(0) for various physical parameters are presented in 2.

6. Conclusions

In this paper we studied the impact of activation energy on magnetized couple stress fluid over an inclined stretching permeable cylinder in a non-Darcy porous medium. The effects of cross-diffusion and stratified mixed convection are also considered in fluid model. The governing PDEs are remoulded into dimensionless nonlinear ODEs by using appropriate similarity transformations and then unravelled by fourth order Runge-Kutta mechanism along with shooting technique. The key results of this study are as follows :Upsurging values of couple stress parameter and magnetic parameters depreciate the velocity field, while improves the temperature and concentration field.

The curvature parameter has enhancing impact on velocity profile, while it reduces the temperature and concentration fields.

Activation energy has portionate impact on concentration field, although Schmidt number has an adverse effect.

The prevailing values of Convection parameter boost the velocity field, but have reverse effect on temperature and concentration fields.

Temperature and concentration fields depreciate with higher thermal and solutal stratification parameter respectively.

Velocity, temperature and concentration fields are strengthened by increasing the values of Forchheimer number.

Temperature field is an increasing function of thermal radiation parameter and decreasing function of Prandtl number.

C (x) Variable ambient concentration (Kmol/m3)
Cw (x) The prescribed surface concentration (Kmol/m3)
Df Dufour number
E Activation energy parameter
C0 Reference concentration (Kmol/m3)
F Forchheimer number
M Magnetic field parameter
T0 Reference temperature (K)
U0 Free stream velocity (ms-1)
Re Reynolds number
Pr Prandtl number
Rd Thermal radiation parameter
Sc Schmidt number
Sr Soret number
g Gravitational acceleration (ms-2)
k Thermal conductivity (Wm-1 K-1)
k' Permeability of porous medium (m2)
cp Specific heat at constant pressure (Jkg-1 K-1)
K Couple stress parameter
L Reference length (m)
B0 Uniform magnetic field (Wbm-2)
X* Space/time dependent heat source parameter
Y* Space/time dependent heat sink parameter
T (x) Variable ambient temperature (K)
Tw (x) The prescribed surface temperature (K)
U(x) Stretching velocity (ms-1)
d Porosity parameter
D Mass diffusivity (m2 s-1)
Ea Activation energy (kJmol-1)
kr Chemical reaction rate constant
Greek symbols
δ1 Thermal stratification parameter
α Inclination of the cylinder
γ Curvature parameter
ζ Temperature difference parameter
λ Combined convection parameter
σ Electrical conductivity (m-1 Ohm-1)
ν Kinematic viscosity (m2 s)
ν' Couple stress viscosity (m2 s)
γ1 Nonlinear temperature convection coefficient
Nomenclature

Funding Statement

This research has been conducted with financial assistance from UGC, India, in the form of a senior research fellowship under UGC-Ref. No.: 1294/ (CSIR-UGC NET JUNE 2019 granted to Suman Sharma.

Conflicts of Interest

The author declares that there is no conflict of interest regarding the publication of this article.

References

  1. Stokes, V. K., 1966. Couple stresses in fluids. Phys. Fluids, 9(9), pp. 1709–1715. doi: 10.1063/1.1761925..DOI
  2. Rani, H.P., G. J. Reddy, G. J., and Kim C. N., 2011. Numerical analysis of couple stress fluid past an infinite vertical cylinder. Eng. Appl. Comput. Fluid Mech., 5(2), pp. 159–169. doi: 10.1080/19942060.2011.11015360..DOI
  3. Gajjela, N. and Garvandha, M., 2020. The influence of magnetized couple stress heat, and mass transfer flow in a stretching cylinder with convective boundary condition, cross-diffusion, and chemical reaction. Therm. Sci. Eng. Prog., 18, pp. 100517. doi: 10.1016/j.tsep.2020.100517..DOI
  4. Ibrahim, W. and Gadisa, G., 2020. Double Stratified Mixed Convective Flow of Couple Stress Nanofluid past Inclined Stretching Cylinder Using Cattaneo-Christov Heat and Mass Flux Model. Advances in Mathematical Physics, 2020(1), p. 4890152. doi: 10.1155/2020/4890152.DOI
  5. Palaiah, S. S., Basha, H., and Reddy G. J., 2021. Magnetized couple stress fluid flow past a vertical cylinder under thermal radiation and viscous dissipation effects. Nonlinear Eng., 10(1), pp. 343–362. doi: 10.1515/nleng-2021-0027.DOI
  6. Awais, M. and Salahuddin, T., 2023. Natural convection with variable fluid properties of couple stress fluid with Cattaneo-Christov model and enthalpy process. Heliyon, 9(8), pp. 18546. doi: 10.1016/j.heliyon.2023.e18546.DOI
  7. Arrhenius, S., 1889. Über die Dissociationswärme und den Einfluss der Temperatur auf den Dissociationsgrad der Elektrolyte / About the heat of dissociation and the influence of temperature on the degree of dissociation of the electrolytes (in German). Zeitschrift für Phys. Chemie.
  8. Bestman, A. R., 1990. Natural convection boundary layer with suction and mass transfer in a porous medium. Int. J. Energy Res., 14(4), pp. 389–396. doi: 10.1002/er.4440140403.DOI
  9. Abbas, Z., Sheikh, M. and Motsa, S. S., 2016. Numerical solution of binary chemical reaction on stagnation point flow of Casson fluid over a stretching/shrinking sheet with thermal radiation. Energy, 95, pp. 12–20. doi: 10.1016/j.energy.2015.11.039.DOI
  10. Waqas, H., Wakif, A., Al-Mdallal, Q., Zaydan, M., Farooq, U. and Hussain, M., 2022. Significance of magnetic field and activation energy on the features of stratified mixed radiative-convective couple-stress nanofluid flows with motile microorganisms. Alexandria Eng. J., 61(2), pp. 1425–1436. doi: 10.1016/j.aej.2021.06.047.DOI
  11. Ibrahim, W., and Negera, M. 2020. The Investigation of MHD Williamson Nanofluid over Stretching Cylinder with the Effect of Activation Energy. Adv. Math. Phys. doi: 10.1155/2020/9523630.DOI
  12. Mustafa, M., Khan, J.A., Hayat, T., and Alsaedi, A., 2017. Buoyancy effects on the MHD nanofluid flow past a vertical surface with chemical reaction and activation energy. Int. J. Heat Mass Transf., 108, pp. 1340–1346. doi: 10.1016/j.ijheatmasstransfer.2017.01.029.DOI
  13. Hamid, A., Khan, M. and Khan, U., 2018. Thermal radiation effects on Williamson fluid flow due to an expanding/contracting cylinder with nanomaterials: Dual solutions. Phys. Lett. Sect. A Gen. At. Solid State Phys., 382(30), pp. 1982–1991. doi: 10.1016/j.physleta.2018.04.057.DOI
  14. Awad, F. G., Motsa, S. and Khumalo, M., 2014. Heat and mass transfer in unsteady rotating fluid flow with binary chemical reaction and activation energy. PLoS One, 9(9). doi: 10.1371/journal.pone.0107622.DOI
  15. Reddy, K. J. and Reddy, M. S., 2013. Effects Of Chemical Reaction And Heat Generation On MHD Boundary Layer Flow Of A Moving Vertical Plate With Suction And Dissipation. i-manager’s J. Futur. Eng. Technol., 8(4), pp. 30–40. doi: 10.26634/jfet.8.4.2359.DOI
  16. Saeed, A., Kumam, P., Gul, T., Alghamdi, W., Kumam, W. and Khan, A., 2021. Darcy–Forchheimer couple stress hybrid nanofluids flow with variable fluid properties. Sci. Rep., 11(1), pp. 1–13. doi: 10.1038/s41598-021-98891-z.DOI
  17. Nagaraja, L., Reddy, M. S. N., 2017. Heat Transfer in a Non-Newtonian Fluid Past a Horizontal Circular Cylinder in Non-Darcy Porous Medium with Suction\Injection Effects. Indian J. Sci. Technol., 10(32), pp. 1–7. doi: 10.17485/ijst/2017/v10i32/104601.DOI
  18. Hayat, T., Saeed, Y., Asad, S. and Alsaedi, A., 2016. Convective heat and mass transfer in flow by an inclined stretching cylinder. J. Mol. Liq., 220, pp. 573–580. doi: 10.1016/j.molliq.2016.03.047.DOI
  19. Ganesh, N. V., Abdul Hakeem, A. K. and Ganga, B., 2018. Darcy–Forchheimer flow of hydromagnetic nanofluid over a stretching/shrinking sheet in a thermally stratified porous medium with second order slip, viscous and Ohmic dissipations effects. Ain Shams Eng. J., 9(4), pp. 939–951. doi: 10.1016/j.asej.2016.04.019.DOI
  20. Pal, D. and Mondal, H., 2010. Effect of variable viscosity on MHD non-Darcy mixed convective heat transfer over a stretching sheet embedded in a porous medium with non-uniform heat source/sink, Commun. Nonlinear Sci. Numer. Simul., 15(6), pp. 1553–1564. doi: 10.1016/j.cnsns.2009.07.002.DOI
  21. Hayat T., Asad, S. and Alsaedi, A., 2017. Non-uniform heat source/sink and thermal radiation effects on the stretched flow of cylinder in a thermally stratified medium. J. Appl. Fluid Mech., 10(3), pp. 915–924. doi: 10.18869/acadpub.jafm.73.240.24008.DOI
  22. Hayat, T., Asad, S., Alsaedi, A. and Alsaadi, F.E., 2015. Radiative flow of jeffrey fluid through a convectively heated stretching cylinder. J. Mech., 31(11), pp. 69–78. doi: 10.1017/jmech.2014.49.DOI
  23. Jain, S., and Choudhary, R., 2017. Soret and Dufour effects on MHD fluid flow due to moving permeable cylinder with radiation. Glob. Stoch. Anal, 4, pp. 75–84.
  24. Khan, N. A., Khan, H. and Ali, S. A., 2016. Exact solutions for MHD flow of couple stress fluid with heat transfer. J. Egypt. Math. Soc., 24(1), pp. 125–129. doi: 10.1016/j.joems.2014.10.003.DOI
  25. Jain, S. and Parmar, A., 2018. Radiation Effect on MHD Williamson Fluid Flow over Stretching Cylinder Through Porous Medium with Heat Source. Lect. Notes Mech. Eng., 0, pp. 61–78. doi: 10.1007/978-981-10-5329-0_5.DOI
  26. Bilal, S., Rehman, K. U. and Malik, M. Y., 2017. Numerical investigation of thermally stratified Williamson fluid flow over a cylindrical surface via Keller box method. Results Phys., 7, pp. 690–696. doi: 10.1016/j.rinp.2017.01.032.DOI
  27. Rehman, K. U., Khan, A. A., Malik, M. Y. and Ali, U., 2017. Mutual effects of stratification and mixed convection on Williamson fluid flow under stagnation region towards an inclined cylindrical surface. MethodsX, 4, pp. 429–444. doi: 10.1016/j.mex.2017.10.007.DOI
  28. Ramzan, M., Bilal, M. and Chung, J. D., 2017. Effects of thermal and solutal stratification on jeffrey magneto-nanofluid along an inclined stretching cylinder with thermal radiation and heat generation/absorption. Int. J. Mech. Sci., 131–132, pp. 317–324. doi: 10.1016/j.ijmecsci.2017.07.012.DOI
  29. Mohd Sohut, N. F. H., Abd Aziz, A. S. and Ali, Z. M., 2017. Double stratification effects on boundary layer over a stretching cylinder with chemical reaction and heat generation. J. Phys. Conf. Ser., 890(1). doi: 10.1088/1742-6596/890/1/012019.DOI
  30. Rehman, A. and Nadeem, S., 2013. Heat Transfer Analysis of the Boundary Layer Flow over a Vertical Exponentially Stretching Cylinder. Glob. J. Sci. Front. Res. Math. Decis. Sci., 13(11)..
  31. Rashid, U., Liang, H. Ahmad, H., Abbas, M., Iqbal, A., and Hamed, S., 2021. Results in Physics Study of (Ag and TiO2)/ water nanoparticles shape effect on heat transfer and hybrid nanofluid flow toward stretching shrinking horizontal cylinder. Results Phys., 21. doi: 10.1016/j.rinp.2020.103812..DOI