# Entropy Analysis on a Three-Dimensional Wavy Flow of Eyring–Powell Nanofluid: A Comparative Study

Authors: Arshad Riaz; Ahmed Zeeshan; M. M. Bhatti
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The thermal management of a system needs an accurate and efficient measurement of exergy. For optimal performance, entropy should be minimized. This study explores the enhancement of the thermal exchange and entropy in the stream of Eyring–Powell fluid comprising nanoparticles saturating the vertical oriented dual cylindrical domain with uniform thermal conductivity and viscous dissipation effects. A symmetrical sine wave over the walls is used to induce the flow. The mathematical treatment for the conservation laws are described by a set of PDEs, which are, later on, converted to ordinary differential equations by homotopy deformations and then evaluated on the Mathematica software tool. The expression of the pressure rise term has been handled numerically by using numerical integration by Mathematica through the algorithm of the Newton–Cotes formula. The impact of the various factors on velocity, heat, entropy profile, and the Bejan number are elaborated pictorially and tabularly. The entropy generation is enhanced with the variation of viscous dissipation but reduced in the case of the concentration parameter, but viscous dissipation reveals opposite findings for the Newtonian fluid. From the abovementioned detailed discussion, it can be concluded that Eyring–Powell shows the difference in behavior in the entropy generation and in the presence of nanoparticles due to the significant dissipation effects, and also, it travels faster than the viscous fluid. A comparison between the Eyring-Powell and Newtonian fluid are also made for each pertinent parameter through special cases. This study may be applicable for cancer therapy in biomedicine by nanofluid characteristics in various drugs considered as a non-Newtonian fluid.


---

## Body

## 1. Introduction

Nonlinear fluid models are the center of consideration of several theoretical and experimental evaluations due to their massive applicability in the engineering sciences, biomechanics, and mechanical manufacturing such as the petroleum industry, flow of blood in the body, and transport of sewage. Working in the domain of nonlinear fluid flows creates challenges to the mathematician and simulation engineers with its diversity and complexity. As the complexity of such fluids offers no unique constitutive equation which encounters all the properties, for such fluids, several non-Newtonian models are presented in [1–4]. The polymeric liquids reveal diverse rheological attributes apparent in the flow, such as polymer accumulation, mechanical mortification, and solvent composition. [5]. Kwack [6] explained that all these polymer properties influence the elasticity of the fluids. Mathematical models such as the Carreau fluid model [7], Ellis fluid model [8], and Powell–Eyring fluid model [9] effectively discuss elastic properties of the fluid. Out of these models, the Powell-Eyring model stands out because of its fluid elasticity property and ease for its applicability to both experimental and theoretical studies. These advantages lead many scientists to unveil these hidden properties of such fluids in different flow situations. Patel and Timol [10] found a numerical solution to debate on the characteristics of the Powell-Eyring type model. They used the so-called MSABC scheme to find an asymptotically convergent solution. Hayat et al. [11] found a series solution using the so-called HAM to define the momentum and heat transfer about the stagnation point for the Powell-Eyring fluid. They also analyzed melting heat transfer effects. Jalil et al. [12] discovered self-similar solutions of parallel free stream flow of this fluid. The flow is considered on a moving plate. Islam et al. [13] considered the flow of the Powell-Eyring fluid in slider bearings. They employed HPM to develop the analytic series solution of the current flow. Hina et al. [14] elaborated the study of Powell-Eyring fluids between curved compliant walls of the channel. The effects on temperature variations are also discussed. Hayat et al. [15] evaluated MHD Powell-Eyring liquid flow on a radially stretching cylinder. They exclusively focused on the effects of Newtonian heating.In recent times, the study of flows driven by peristaltic waves has become more and more popular in physiology and biomedical industry. The idea of peristalsis was driven from the food transport duct in the human body. Other than food transport from the mouth to stomach, this phenomenon can be seen in the spermatozoa movement in the male genital region and ovum in the females and also in urine travelling and chyme movement through the gastrointestinal tract, etc. In industry, machine-like rollers and finger pumps and blood filtration devices in dialysis employ the mechanism of peristalsis. Shapiro [16] and Latham [17] presented the idea of peristalsis theoretically and verified experimentally. Ever since, numerous studies, both theoretically and experimentally, have been published in the domain of nonlinear fluids [18–20]. The peristaltic mechanism of certain elastic fluids is a matter of great interest. Recent studies are indicative of the importance of these in industrial applications. Ellahi et al. [21] used 3D analysis of a peristaltic stream of Carreau fluid in the magnetic environment. Hayat et al. [22] worked with the Carreau-Yasuda fluid with nanosized particle theoretically. They considered that flow is driven due to the sinusoidal motion of the wall. Also, they examined both convective and no-mass flux on the wall. Prakash and Tripathi [23] examined the effects of the electric double layer on Williamson ionic liquid with homogeneously distributed particles of nanosize. They considered a tapered channel with peristaltic waves on the wall. Effects of thermal radiation are also discovered. Zeeshan et al. [24] focused on the flow of the bio-Jeffrey fluid and the sustainability character in the human body. Tanveer et al. [25] checked the flow through a curved channel. They revealed mixed convection impacts and elastic properties of the walls on pumping transport of Eyring-Powell liquid with nanomaterial. Bhatti et al. [26] presented the effects of the Darcy porous medium of a two-phase nanofluid in a channel.The applications such as paper drying, thermal coating, hot rolling, and glass fiber stretching divert one’s focus towards the heat transfer of non-Newtonian fluids. Motivated by this, it stimulates an extensive literature in this field [27–30] which includes analysis in flexible channels. Presence of nanoparticles shows a positive trend in enhancing the heat transfer of fluid. The theme of nanofluids is the concept of Choi [31]. The concept was to add nanosized particles in a fluid which has low thermal conductivity. Hence, fluid behavior and metallic conductivity both combine to obtain an enhanced heat conductor which can be transported like a fluid. It was one of the promising aspects of such fluids. Many mathematical models using continuum formulation [32] or two-phase suspension models [33] are developed over the years now. One such model due to Buongiorno [34] was developed in many ways. Some researchers [35–39] revealed the basics of nanofluids’ contribution in different geometries with the peristaltic wave.It is essential to generate the systems or the machines which perform and help in efficient energy transport. Such a machine is broadly utilized in power plants, manufacturing plants, and transportation. The warm productivity of the framework relies upon the material used for thermal vitality. It is fundamental to understand the factor which diminishes warm proficiency. As indicated by the laws of thermodynamics, vitality of the framework stays to save, yet can be changed over into different structures for the utility. All the more regularly, we state all the energies of the framework are spent in accomplishing work or to expand the temperature of the body. Numerous frictional powers emerging from the attractive field, permeable spaces, and so on are explained in the improvement of the temperature of the framework. The warm messiness in the framework is consistently on the ascent. Realizing the elements engaged with entropy rise is consistently imperative to get ideal warm effectiveness of the framework. Ascend being used of nanofluids is one approach to decrease the loss of vitality. Entropy generation and exergy are one of the most dynamic examination fields in thermodynamics. Most recently, Naz et al. [40] optimized entropy generated in the flow of pseudoplastic fluid. The fluid also comprises of microorganisms. Abbas et al. [41] analyzed the entropy generated theoretically in the stream of magnetized nanofluid. They use statistical techniques to get results. Abbas et al. [42] also analyzed entropy for peristaltic flow of nanofluids. Ellahi et al. [43] extended the work with a porous medium. Hayat et al. [44] use modified Darcy’s porous medium and the effects of the endoscope.The limited literature on entropy generation of nanoelastic fluid, useful characteristics of the Eyring-Powell model, and the applicability in the industrial usage of such flows is the main motivation for the current analysis. With the cited literature in mind, the existing study designed to investigate the distribution of nanoparticles in a Powell-Eyring viscoelastic fluid having enhanced elastic properties and viscous dissipation effects through a cylindrical structure with the peristaltic wave is focused. Energy loss in the transport is also an area of interest. The observations for the Newtonian fluid are also taken into consideration as a limiting case of the readings. The obtained PDEs are simplified with the help of laminar movement of the liquid and small wave amplitude approximations and then solved analytically. The results are elaborated comprehensively and displayed in the graphical way. Conclusions which highlight the key findings are placed in the end.


## 2. Mathematical Treatment

We take the peristaltic flow of an Eyring-Powell liquid situated in an annular part of two eccentric cylinders by introducing nanoparticles. We also consider the term of viscous dissipation and entropy generation occurring during the process due to the irreversibilities of viscous factors, thermal analysis, and mass exchanging phenomenon of nanoparticles. The geometry is assumed to be having flexible outer surfaces. The flow is taken inside the annular region of two eccentric annuli with eccentricityϵ. Moreover, the inner cylinder is considered as rigid with radius δ, and the outer cylinder with radius ais producing a sinusoidal wave propagating along the axial root. The walls of the inner and outer annuli have been described in [30]. The conservation laws for an incompressible and unsteady flow of non-Newtonian nanofluid are described in [30, 35]. The temperatures and nanoparticle concentrations at the inner cylinder are fixed at the amount of Toand Co, while those of the outer cylinder are maintained at T1and C1≃accordingly. The mathematical model of the stress matrix representing the Eyring-Powell fluid is described as follows [15]:(1)SLM=μ∂Vp∂xi+1βArc sinhc−1∂Vp∂xi,(2)Arc sinhc−1∂Vp∂xi≈c−1∂Vp∂xi−16c−1∂Vp∂xi3+Oc−1∂Vp∂xi5,where μrepresents the shear viscosity coefficient and βand care the material constants. Now, we introduce local wave frame coordinates r,θ,z running with the same wave velocity relative to the fixed frame x1,x2,x3 through the following transformations:(3)z=x3−c1t,r=x1,θ=x2,w=W¯−c1,p=P,u=U¯,SLM=slm.We use the following dimensionless parameters in the problem:(4)p′=a2μc1λp,w′=wc1,u′=λac1u,V′=Vc1,z′=zλ,r′=ra,θ′=θ,t′=cλt,ϕ=ba,ϵ′=ϵa,Re=ρc1aμ,δ′=δa,θ¯=T−ToT1−To,δo=aλ,r1′=r1,r2′=r2a,σ=C−CoC1−Co,Pr=μρα,Sc=μρDB,αf=kρcf,τ=ρcpρcf,Br=ρfgαa2μcT1−To,Gr=ρfgαa2μcC1−Co,Nb=τDBαfC1−Co,Nt=τDTToαfT1−To,S̖=aμc1S,Gc=μc2kT1−To,M=1βμc,K=Mc126c2a2,where V,ϕ,Re,δo,Pr,Nb,Nt,Gc,Gr, and Brevaluate the velocity of the inner rod, amplitude fraction, Reynolds number factor, nondimensional wave number, the Prandtl number, the Brownian motion parameter, the thermophoresis constant parameter, Brinkman number, localized temperature Grashof number, and localized nanoparticles Grashof number, separately. After employing the abovementioned nondimensional parameters and wave frame, we approach to the following dimensionless form of the equations of motion defined in [30] under the physiological limitations of long wavelength λ⟶∞ and low Reynolds number Re⟶0:(5)∂u′∂r′+∂w′∂z′+u′r′=0,(6)0=−∂p′∂z′+1r′∂r′Srz′∂r′+1r′∂Sθz′∂θ′+Grθ¯+Brσ,(7)0=∂2θ′¯∂r′2+1r′2∂2θ′¯∂θ′2+1r′∂θ′¯∂r′+Nb∂σ∂r′∂θ′¯∂r′+1r′2∂σ∂θ′∂θ′¯∂θ′+Nt∂θ′¯∂r′2+1r′2∂θ′¯∂θ′2+GcSrz′∂w′∂r′+1r′Szθ′∂w′∂θ′,(8)0=∂2σ∂r′2+1r′2∂2σ∂θ2+1r′∂σ∂r′+NtNb∂2θ′¯∂r′2+1r′2∂2θ′¯∂θ′2+1r′∂θ′¯∂r′.The needed nondimensional components of stress tensorSLMcalculated from relation (1) can be collected as follows:(9)Srz′=∂w′∂r′+M∂w′∂r′−K∂w′∂r′3,Szθ′=1r′∂w′∂θ′+Mr′∂w′∂θ′−K1r′∂w′∂θ′3,Szθ′=0.Now, we useSrz',Sθz', and Szθ' in (6) and (7) and also omit the prime symbols to have the resulting structure as follows:(10)11+Mp′z=∂2w∂r2+1r∂w∂r+1r2∂2w∂θ2−3K1+M∂w∂r2∂2w∂r2−Kr1+M∂w∂r3−3Kr41+M∂w∂θ2∂2w∂θ2+11+MGrθ¯+Brσ,(11)0=∂2θ¯∂r2+1r∂θ¯∂r+1r2∂2θ¯∂θ2+Nb∂σ∂r∂θ¯∂r+1r2∂σ∂θ∂θ¯∂θ+Nt∂θ¯∂r2+1r2∂θ¯∂θ2+Gc∂w∂r2+M∂w∂r2−K∂w∂r4,(12)0=∂2σ∂r2+1r∂σ∂r+1r2∂2σ∂θ2+NtNb∂2θ¯∂r2+1r∂θ¯∂r+1r2∂2θ¯∂θ2.The relevant no-slip walls’ conditions in a nondimensional format are described as follows [30]:(13)w=Vwhenr=r1=δ+εcosθ,0whenr=r2=1+ϕcos2πz,θ¯,σ=0whenr=r1=δ+εcosθ,1whenr=r2=1+ϕcos2πz.


## 3. Solution Procedure

The solutions of the overhead nonlinear partial equations, (10), (11), and (12), are achieved analytically. The deformation expressions for the considering problems are evaluated as follows [30, 37]:(14)1−q″ℒw˜−ℒw˜o+q″ℒw˜+1r2∂2w˜∂θ2−3K1+M∂w˜∂r2∂2w˜∂r2−Kr1+M∂w˜∂r3−3Kr41+M∂w˜∂θ2∂2w˜∂θ2+11+MGrΘ+BrΩ−p′z=0,(15)1−qℒΘ−ℒθ˜o+q″ℒΘ+1r2∂2Θ∂θ2+Nb∂Θ∂r∂Ω∂r+1r2∂Θ∂θ∂Ω∂θ+Nt∂Θ∂r2+1r2∂Θ∂θ2+Gc∂w˜∂r2+M∂w˜∂r2−K∂w˜∂r4=0,(16)1−qℒΩ−σ˜o+q″ℒΩ+1r2∂2Ω∂θ2+NtNb∂2Θ∂r2+1r2∂2Θ∂θ2+1r∂Θ∂r=0.Let us choose the linear operator asℒ=∂/∂r2+1/r∂/∂r. We observe the following initial choices for w,θ¯,andσ:(17)w˜o=Vlogr−logr2logr1−logr2−1,θ¯o=σ¯o=logr1−logrlogr1−logr2−1.Now, using the perturbation technique on the parameterq″, we have(18)w˜r,θ,z,q″=wo+q″w1+…,(19)Θ˜r,θ,z,q″=θ¯o+q″θ¯1+…,(20)Ω˜r,θ,z,q″=σo+q″σ1+….After substituting the abovementioned equations into equations (14)–(16), the solution of velocity w1 contains a very large output, so it is not displayed here, and the solutions of θ¯1 and σ1 are given in Appendix. Now, for (q″⟶1), we approach the final solutions. So, from equations (18)–(20), we get(21)wr,θ,z=wo+w1,θ¯or,θ,z=θo+θ1,σr,θ,z=σo+σ1.The instant value of volume rate of flowQ″¯ is taken by(22)Q″¯=2π∫r1r2rwdr.The periodic mean volumetric rate of streamQ is found mathematically as follows [30]:(23)Q″¯z,t=Qπ−ϕ22+2ϕcos2πz−t+ϕ2cos22πz.One can get the measure of pressure ratep′z by decoding (22) and (23). The dimensionless pressure rise _Δpis defined as follows [30]:(24)Δp=∫01p′zdz.The integration is performed numerically by using Newton–Cotes formulas through the NIntegrate tool on Mathematica 7.


## 4. Entropy Measures

The expression of entropy generation under the irreversibility of heat exchange, nanoparticles, and viscous effects of the nanofluid is defined as follows [30]:(25)Sgen′=KTo2∂T∂r2+1r2∂T∂θ2+∂T∂z2+DBCo∂C∂r2+1r2∂C∂θ2+∂C∂z2+DBTo∂C∂r∂T∂r+1r2∂C∂θ∂T∂θ+∂C∂z∂T∂z+1ToSrr∂u∂r+1rSrθ∂u∂θ+Srz∂u∂z+∂w∂r+1rSzθ∂w∂θ+Szz∂w∂z+urSθθ.The dimensionless factors used in these equations are considered as follows:(26)Ns=Sgen′SG,SG=KT1−To2a2To2,Γ=DBToC1−CoKT1−To,Λ=T1−ToTo,Ω=C1−CoCo.Now, we substitute the abovementioned manipulated factors in equation (25); we have(27)Ns=∂θ¯∂r′2+1r′2∂θ¯∂θ′2+δO2∂θ¯∂z′2+ΓΛΩ∂σ∂r′2+1r′2∂σ∂θ′2+δO2∂σ∂z′2+Γ∂σ∂r′∂θ¯∂r′+1r′2∂σ∂θ′∂θ¯∂θ′+δO2∂σ∂z′∂θ¯∂z′+GcΩδoSrr′∂u′∂r′+δo1r′Srθ′∂u′∂θ′+Srz′δO2∂u′∂z′+∂w′∂r′+1r′Szθ′∂w′∂θ′+δoSzz′∂w′∂z′+δou′r′Sθθ′.Incorporating the lubrication strategy in the abovementioned equation, it will take the form(28)Ns=∂θ¯∂r′2+1r′2∂θ¯∂θ′2+ΓΛΩ∂σ∂r′2+1r′2∂σ∂θ′2+Γ∂σ∂r′∂θ¯∂r′+1r′2∂σ∂θ′∂θ¯∂θ′+GcΩSrz′∂w′∂r′+1r′Szθ′∂w′∂θ′.Now, using the values ofSrz′ and Szθ′ in the abovementioned equation,(29)Ns=∂θ¯∂r′2+1r′2∂θ¯∂θ′2+ΓΛΩ∂σ∂r′2+1r′2∂σ∂θ′2+Γ∂σ∂r′∂θ¯∂r′+1r′2∂σ∂θ′∂θ¯∂θ′+GcΩ∂w′∂r′∂w′∂z′+M∂w′∂r′∂w′∂z′−K∂w′∂r′3∂w′∂z′.By ignoring the prime symbols of the abovementioned equation, we obtain(30)Ns=∂θ¯∂r2+1r2∂θ¯∂θ2+ΓΛΩ∂σ∂r2+1r2∂σ∂θ2+Γ∂σ∂r∂θ¯∂r+1r2∂σ∂θ∂θ¯∂θ+GcΩ∂w∂r∂w∂z+M∂w∂r∂w∂z−K∂w∂r3∂w∂z.Moreover, the Bejan number (Be) is defined as [30](31)Be=∂θ¯∂r2+1r2∂θ¯∂θ2∂θ¯∂r2+1r2∂θ¯∂θ2+ΓΛΩ∂σ∂r2+1r2∂σ∂θ2+Γ∂σ∂r∂θ¯∂r+1r2∂σ∂θ∂θ¯∂θ+GcΩ∂w∂r∂w∂z+M∂w∂r∂w∂z−K∂w∂r3∂w∂z−1.


## 5. Discussion on Results

This segment is assigned to giving the graphical appearance dedicating the key features of the study by drawing figures against the most pertinent parameters by taking the solution series up to the first two orders of the summation. These results also reflect the physical logic through a mathematical model of the real phenomenon. It also helps the engineers and scientist to work on such problems experimentally. The numerical data of velocity, temperature, entropy, and Bejan number profiles have been presented in Tables1–4, respectively. From Tables 1 and 2, one can clarify that this study reproduces the results of Nadeem et al. [45] which notify the validity of the current analysis. These tables have been prepared also to estimate the variations of physical properties and imagine a difference of the behavior of Newtonian and non-Newtonian fluids under similar circumstances.Table 1
Comparative data of velocityw.
rNadeem et al. [45]Current study whenK=M=0Eyring-Powell fluid whenK=M=0.20.10.936210.936211.073050.21.202801.202801.065570.31.338091.338091.237870.41.369551.369551.323230.51.313541.313541.308760.61.178961.178961.200330.70.970920.970921.004030.80.692510.692510.724430.90.345690.345690.364791.0−0.06832−0.06832−0.072581.1−0.54873−0.54873−0.58602Table 2
Comparative data of temperature profileθ.
rNadeem et al. [45]Current study whenK=M=0Eyring-Powell fluid whenK=M=0.20.1−2.65812−2.658120.290810.20.007170.007170.003230.30.940410.940410.590330.41.331871.331870.957610.51.484321.484321.155080.61.512571.512571.244050.71.469841.469841.263060.81.384111.384111.235710.91.271341.271341.176831.01.141251.141251.096101.11.000001.000001.00000Table 3
Comparative data of entropy generation profileNs.
rNewtonian fluid whenK=M=0Eyring-Powell fluid whenK=M=0.20.150.193862.96150.210.551812.11260.34.346144.791250.42.335292.512820.51.446901.531050.60.979601.023530.70.704570.728680.80.529480.542880.90.411340.418571.00.328000.331491.10.267080.26822Table 4
Comparative data of the Bejan number profileBe.
rNewtonian fluid whenK=M=0Eyring-Powell fluid whenK=M=0.20.10.183890.159350.20.331300.292250.30.463490.415170.40.583630.531090.50.695740.643700.60.803770.756910.70.911420.874780.81.022161.001750.91.139511.143041.01.267281.305301.11.409941.49768From here, we can clearly find out the viscous effects on the different viscosity model which are also manipulated latterly through graphics. The graphs for velocity of the nanofluid, temperature profile, entropy generation, and Bejan number have also been invoked in this study. The data for obtained quantities have been collected from the achieved analytical solutions. The main focus is kept on analyzing the difference of behavior shown by NewtonianM=K=0 and the Eyring-Powell fluid model M=K=0.2 under the variation of different flow factors. The velocity radial component w can be understood diagrammatically from Figures 1–3. Figures 4–6 exhibit the deportment of temperature distribution θ. The entropy generation analysis can be perceived from graphs of NS which are placed in Figures 7–10, and the Bejan number variation has been offered through Figures 11–14.Figure 1
Structure of velocity factorw with Gr for fixed θ=0.8,ϕ=0.1, Br=1,δ=0.1,Q=0.1,K=2,M=5.Figure 2
Structure of velocity factorw with Br for fixed θ=0.8, ϕ=0.1,δ=0.1,Gr=1,Q=0.1,K=2,M=5.Figure 3
Structure of velocity factorw with ϵ for fixed θ=1,ϕ=0.1,δ=0.1,Gr=1,Br=1,Q=0.25,K=0.1,M=0.1.Figure 4
Alteration of heat factorθ¯ for Nb when θ=0.1Br=0.5,δ=0.1,ϕ=0.5,V=0.3,K=1,M=1,Gc=1,Nt=0.2,ɛ=0.1.Figure 5
Alteration of heat factorθ¯ for Nt when θ=0.1, ϕ=0.5, Nb=0.2,δ=0.1,Br=0.5,V=0.3,K=1,M=1,Gc=1,ɛ=0.1.Figure 6
Alteration of heat factorθ¯ for Gc when θ=0.1, ϕ=0.5, Nt=0.2,δ=0.1,Br=0.5,V=0.3,K=1,M=1,Nb=1,ɛ=0.1.Figure 7
ForNs with fixed Gc where δ=0.01,θ=0.1,ϕ=0.8,ɛ=0.01,Γ=0.4,Nb=0.9,Nt=0.5,Br=1,Gr=1,V=0.3,K=2,M=1,Q=1,Λ=0.1,Ω=0.3.Figure 8
ForNs with fixed Ω where δ=1.5,θ=0.1,ϕ=2,ɛ=0.01,Γ=0.4,Nb=0.9,Nt=0.5,Br=1,Gr=5,V=0.3,K=1,M=1,Q=1,Λ=1,Gc=0.01.Figure 9
ForNs with fixed Γ where δ=0.01,θ=0.3,ϕ=0.8,ϵ=0.01,Gc=0.01,Nb=0.9,Nt=0.5,Br=1,Gr=3,V=0.3,K=1,M=1,Q=1,Λ=0.4,Ω=0.3.Figure 10
ForNs with fixed Λ where δ=0.01,θ=0.1,ϕ=0.8,ϵ=0.1,Gc=0.01,Nb=0.9,Nt=0.5,Br=1,Gr=3,V=0.1,K=1,M=1,Q=1,Γ=0.4,Ω=0.3.Figure 11
ForBe, Be with fixed Gc where δ=0.1,θ=0.1,ϕ=0.5,ϵ=0.01,Λ=1,Nb=0.9,Nt=0.5,Br=1,Gr=5,V=0.3,K=1,M=3,Q=1,Γ=0.4,Ω=0.3.Figure 12
ForBe, Be with fixed Ω where δ=0.1,θ=0.1,ϕ=0.5,ɛ=0.2,Λ=1,Nb=0.9,Nt=0.5,Br=1,Gr=5,V=0.3,K=1,M=3,Q=1,Γ=0.4,Gc=0.01.Figure 13
ForBe with fixed Γ where δ=0.1,θ=0.1,ϕ=0.5,ɛ=0.2,Λ=0.1,Nb=0.9,Nt=0.5,Br=5,Gr=1,V=0.2,K=1,M=1,Q=1,Ω=0.3,Gc=0.01.Figure 14
ForBe with fixed Λ where δ=0.1,θ=0.1,ϕ=0.5,ɛ=0.2,Γ=0.4,Nb=0.9,Nt=0.5,Br=1,Gr=5,V=0.3,K=1,M=1,Q=1,Ω=0.3,Gc=0.01.


### 5.1. Velocity Profile

It is observed from Figure1 that the Eyring-Powell fluid travels faster as compared to the Newtonian fluid under the variation of the Grashof numberGr. It is deeply observed that the rise in values of Grashof number Gr shows that the peak value tends to shift toward the upper cylinder. It is for the verity that uplift in the values of Gr buoyancy impacts raising the fluid, which happens to generate free convection. Normally, the density tends to rise with the reduction in temperature, but in case of the Newtonian fluid, the uniform viscosity measures the inverse behavior against the same factor. With the large value of Br, the energy dissipates due to thermal enhancement and, hence, rises the flow near the outer cylinder which are observed in Figure 2, and the same observation is reported for viscous and nonlinear fluids. In Figure 3, variation of velocity is observed for the change in eccentricity of the cylinders. It is observed that velocity declines when the eccentricity enhanced, whereas it shows the opposite behavior for the Newtonian fluid case. It is found due to the fact that Eyring-Powell fluid is higher in its viscosity with the detraction of the inner cylinder from the central axis which was not in the case of a linear model.


### 5.2. Temperature Dispersion Profile

The variation in temperature distribution is observed from Figures4–6. The parameters of interest are the Brownian motion factor Nb, thermophoresis coefficient Nt, and nanoparticles Grashof factor Gc. Examining the results, we observe that a more parabolic nature is exhibited in case of temperature profile than velocity curves. With the rise in the value of the Brownian motion parameter Nb, Brownian motion increases that leads to a rise in temperature of the fluid; it is validated from Figure 4. In Figure 5, the variation of the thermophoresis parameter Nt is observed. It is clear from here that temperature increases in the flow domain since the rise in thermal diffusion. The Brinkman number Gc shows a reduction in temperature in the flow domain due to the rise in energy dissipation due to viscosity diffusion which is contradictory to the Newtonian fluid case where an opposite behavior is provoked (see Figure 6). In the abovementioned three discussed results, it is very obvious that temperature variation for the viscous and nonlinear model of fluid shows a similar situation against Nb and Nt, but there are inverse calculations for the two different phases of fluid just because of the viscosity alteration aspect.


### 5.3. Entropy Measurements

The factor which is the most important one for thermal analysis of liquid flow is the dissipation analysis of energy which is defined by the entropy generation parameterNs. The variation in entropy generation is graphically displayed for various values GC,Ω,Γ, and Λ in Figures 7–10, respectively. Entropy generation decreases for larger values of the Brinkman number GC since temperature profile shows a rise for Gc (Figure 7). From Figure 8, the findings made are described as with the rise in the concentration parameter, entropy decreases. An opposite behavior of the temperature parameter Γ is observed from Figure 9. From Figure 10, it is reflected that the ratio of the concentration parameter Λ is responsible for higher entropy generation. From Figures 7–10, the interesting aspect which is pointed out is that only Gcputs opposite impacts on viscous fluid but other parameters make no prominent change for the Newtonian fluid, and an important notable thing is that entropy peaks are measured in the neighborhood of the inner cylinder which becomes flat at the surface of the outer one.


### 5.4. Bejan Number Curves

A graphical representation ofBeis displayed in Figures 11–14. It is defined as the entropy generation due to heat loss divided by total entropy. It is observed from Figure 11 that, with a rise in the Brinkman number Gc, the Bejan Number Be increases due to the increase in thermal entropy for the rise in temperature which comes from high values of the Brinkman number. It indicates that the contribution of Gc in entropy increases with r for higher values of Gc. Figure 12 shows a decrease in the entropy ratio for large values of Ω, which means that an increase in the concentration parameter produces less entropy due to thermal transfer and more entropy due to collective aspects generating entropy, whereas the temperature parameter Γ and temperature to concentration parameter Λ show an opposite behavior than the Brinkman number as seen from Figures 13 and 14, respectively. On the parallel side, these plots also investigate that the Bejan number curves are reflecting similar results for both the viscous and nonlinear fluids. The only thing that is prominent from these graphs is that the thermal entropy is maximum near the outer cylinder.


## 5.1. Velocity Profile

It is observed from Figure1 that the Eyring-Powell fluid travels faster as compared to the Newtonian fluid under the variation of the Grashof numberGr. It is deeply observed that the rise in values of Grashof number Gr shows that the peak value tends to shift toward the upper cylinder. It is for the verity that uplift in the values of Gr buoyancy impacts raising the fluid, which happens to generate free convection. Normally, the density tends to rise with the reduction in temperature, but in case of the Newtonian fluid, the uniform viscosity measures the inverse behavior against the same factor. With the large value of Br, the energy dissipates due to thermal enhancement and, hence, rises the flow near the outer cylinder which are observed in Figure 2, and the same observation is reported for viscous and nonlinear fluids. In Figure 3, variation of velocity is observed for the change in eccentricity of the cylinders. It is observed that velocity declines when the eccentricity enhanced, whereas it shows the opposite behavior for the Newtonian fluid case. It is found due to the fact that Eyring-Powell fluid is higher in its viscosity with the detraction of the inner cylinder from the central axis which was not in the case of a linear model.


## 5.2. Temperature Dispersion Profile

The variation in temperature distribution is observed from Figures4–6. The parameters of interest are the Brownian motion factor Nb, thermophoresis coefficient Nt, and nanoparticles Grashof factor Gc. Examining the results, we observe that a more parabolic nature is exhibited in case of temperature profile than velocity curves. With the rise in the value of the Brownian motion parameter Nb, Brownian motion increases that leads to a rise in temperature of the fluid; it is validated from Figure 4. In Figure 5, the variation of the thermophoresis parameter Nt is observed. It is clear from here that temperature increases in the flow domain since the rise in thermal diffusion. The Brinkman number Gc shows a reduction in temperature in the flow domain due to the rise in energy dissipation due to viscosity diffusion which is contradictory to the Newtonian fluid case where an opposite behavior is provoked (see Figure 6). In the abovementioned three discussed results, it is very obvious that temperature variation for the viscous and nonlinear model of fluid shows a similar situation against Nb and Nt, but there are inverse calculations for the two different phases of fluid just because of the viscosity alteration aspect.


## 5.3. Entropy Measurements

The factor which is the most important one for thermal analysis of liquid flow is the dissipation analysis of energy which is defined by the entropy generation parameterNs. The variation in entropy generation is graphically displayed for various values GC,Ω,Γ, and Λ in Figures 7–10, respectively. Entropy generation decreases for larger values of the Brinkman number GC since temperature profile shows a rise for Gc (Figure 7). From Figure 8, the findings made are described as with the rise in the concentration parameter, entropy decreases. An opposite behavior of the temperature parameter Γ is observed from Figure 9. From Figure 10, it is reflected that the ratio of the concentration parameter Λ is responsible for higher entropy generation. From Figures 7–10, the interesting aspect which is pointed out is that only Gcputs opposite impacts on viscous fluid but other parameters make no prominent change for the Newtonian fluid, and an important notable thing is that entropy peaks are measured in the neighborhood of the inner cylinder which becomes flat at the surface of the outer one.


## 5.4. Bejan Number Curves

A graphical representation ofBeis displayed in Figures 11–14. It is defined as the entropy generation due to heat loss divided by total entropy. It is observed from Figure 11 that, with a rise in the Brinkman number Gc, the Bejan Number Be increases due to the increase in thermal entropy for the rise in temperature which comes from high values of the Brinkman number. It indicates that the contribution of Gc in entropy increases with r for higher values of Gc. Figure 12 shows a decrease in the entropy ratio for large values of Ω, which means that an increase in the concentration parameter produces less entropy due to thermal transfer and more entropy due to collective aspects generating entropy, whereas the temperature parameter Γ and temperature to concentration parameter Λ show an opposite behavior than the Brinkman number as seen from Figures 13 and 14, respectively. On the parallel side, these plots also investigate that the Bejan number curves are reflecting similar results for both the viscous and nonlinear fluids. The only thing that is prominent from these graphs is that the thermal entropy is maximum near the outer cylinder.


## 6. Concluding Remarks

In the current paper, we have elaborated entropy measures occurring during the biological movement of a viscoelastic-type fluid in between an annular region of two cylinders having distinct centers. The nanoparticles phenomenon is also taken into consideration. The equations have been assembled in the light of the lubrication strategy and steady format by transformation used in a newly arranged wave frame of reference. In terms of solution, expressions have been analytically gathered with the help ofHPM and numerical integration is operated to the expression of pressure rise on mathematical computing the Mathematica tool. The expressions of the entropy generation factor and Bejan number are incorporated by obtained data of viscous dissipation, heat, and energy terms. Comparison of two different models of the fluid is kept under key observation in all figures. The final headings extracted from the whole analysis are reflected in the following:(i)
It is found that the velocity of the substance is getting lower by increasing eccentricity measures of the cylinders and the local temperature Grashof number but higher in front of the local nanoparticles Grashof number. It is also seen that the Newtonian fluid shows opposite readings against the eccentricity and local nanoparticles Grashof number but the same for the temperature Grashof number(ii)
Temperature of the fluid is enhanced for both of the factors, the Brownian motion parameter and thermophoresis parameter, but suppressed against the Brinkman number which is a contradictory result in the Newtonian fluid(iii)
It is noted that entropy generation gets lowered with the increase in the concentration parameter, but it is enhanced under the variation of the Brinkman number, temperature variance representative, and the factor showing the ratio of the thermal factor and concentration factors, but the difference of variation for the viscous fluid is only reported in the case of the Brinkman number(iv)
Bejan number curves are showing totally inverse characteristics than the entropy generation curves(v)
It can be resulted from abovementioned detailed discussion that the Eyring-Powell fluid shows the difference of behavior mostly for viscous dissipation effects in entropy generation and in the presence of nanoparticles, and it travels faster than the viscous fluid(vi)
This study is validated by the result of Nadeem et al. [45] and also compares the characteristics of both the viscous and Eyring-Powell fluid model through special cases K=M=0



---

*Source: 6672158-2021-03-05.xml*