An assessment of the mathematical model for estimating of entropy optimized viscous fluid flow towards a rotating cone surface.

Entropy optimization in convective viscous fluids flow due to a rotating cone is explored. Heat expression with heat source/sink and dissipation is considered. Irreversibility with binary chemical reaction is also deliberated. Nonlinear system is reduced to ODEs by suitable variables. Newton built in shooting procedure is adopted for numerical solution. Salient features velocity filed, Bejan number, entropy rate, concentration and temperature are deliberated. Numerical outcomes for velocity gradient and mass and heat transfer rates are displayed through tables. Assessments between the current and previous published outcomes are in an excellent agreement. It is noted that velocity and temperature show contrasting behavior for larger variable viscosity parameter. Entropy rate and Bejan number have reverse effect against viscosity variable. For rising values of thermal conductivity variable both Bejan number and entropy optimization have similar effect.

Introduction

Influence of variable viscosity (temperature dependent viscosity) for flow of fluids is more realistic. Augmentation in temperature leads to decay of viscosity of liquids while gases viscosity enhances. In oiling liquids the enhancement in heat creates inner resistance which distresses the fluid viscosity, and therefore viscosity of liquid does not remain constant. Thus it is described to scrutinize the impact of different temperature variable viscosity. Mukhopadhyay and Layek[1] studied the radiative convective flow by a porous stretchable surface with temperature dependant viscosity. Impact of variable viscosity in an unsteady magnetohydrodynamic convection flow is investigated by Seddeek[2]. Salient features of variable properties for thin film flow is explored by Khan et al.[3]. Hayat et al.[4] studied unsteady convective viscous liquids flow. Effect of heat flux on unsteady magnetohydrodynamic viscous liquids flow over a rotating disk is discoursed by Turkyilmazoglu[5]. Hayat et al.[6] scrutinized the behavior of chemical reaction in Jeffrey liquid flow with variable thermal conductivity. Some relevant attempts about variable properties made in Refs.[7-10].

The ability of noteworthy improvement apparatus such as spinning cone columns, centrifugal disc atomizers, fluid degausser, rotating packed-bed reactors and centrifugal film evaporators etc. depends upon the nature of motion of liquid and pressure distributions. Rotating cone has utilizations in engineering field, advanced nanotechnology and industrial sites including nuclear reactor, liquid film evaporators and cooling system etc. Shevchuk[11] successfully presented the novel numerical and analytical simulations for the various rotating flows like system rotation, swirl flows associated with the swirl generators and surface curvature in bends as well as turns. The impact of centrifugal and Coriolis forces on the distinct flow pattern due to rotating flows was also successfully presented in this scientific continuation. The work of Shevchuk[12] visualized the impact of wall temperature in order to inspect the heat transfer characteristics in the laminar flow confined by rotating disk. The analytical solutions for the formulated rotating disk problems were also successfully addressed. In interesting another continuation, Shevchuk[13] modeled the turbulent flow problem in presence of heat transfer phenomenon due to rotating disk. The applications of heat and mass transfer pattern in rotating flow of cone and plate devices has been pointed out by Shevchuk[14]. Turkyilmazoglu[15] presented the analytical solutions for a rotating cone problem for viscous fluid. In another continuation, Turkyilmazoglu[16] inspected the heat transfer pattern in viscous fluid confined by a rotating cone. Behaviors of variable properties on mixed convection viscous liquid flow with dissipation over a rotating cone are deliberated by Malik et al.[17]. Turkyilmazoglu[18] analyzed the fluctuation in heat transfer mechanism for viscous fluid flow configured by rotating disk in with porous space. Impact of variable viscosity in magnetohydrodynamic flow of Carreau nanofluid by a rotating cone is illustrated by Ghadikolaei et al.[19]. Sulochana et al.[20] studied radiative magnetohydrodynamic flow of laminar liquid with Soret effect over a rotating cone. Salient behaviors of thermal flux in unsteady MHD convective flow due to a rotating cone are presented by Osalusi et al.[21]. Turkyilmazoglu[22] addressed the radially impacted flow of viscous fluid accounted by rotating disk. Asghar et al.[23] used Lie group approach to simulate the solution for a rotating flow problem in presence of heat transfer. Turkyilmazoglu[24] visualized the flow pattern of triggered fluid due to rotating stretchable disk. The fluid flow due to stationary and moving rotating cone subject to the magnetic force impact has been depicted by Turkyilmazoglu[25].

With excellent thermal effectiveness and multidisciplinary applications, the study of nanoparticles becomes the dynamic objective of scientists. The valuable importance of nano-materials in distinct processes includes solar systems, technological processes, engineering devices, nuclear reactors, cooling phenomenon etc. With less than 100 nm size and structure, the nanoparticles are famous due to extra-ordinary thermal performances in contrast to base liquids. In modern medical sciences, the nanoparticles are used to demolish the precarious cancerous tissues. Choi[26] presents the novel investigation on nanofluids and examined the extra-ordinary thermal activities of such materials. Later on, many investigations are claimed in the literature to analyze the thermal assessment of nano-materials. For example, Chu[27] explained the thermal aspects of third grade nanofluid with significances of activation energy and microorganisms. Majeed et al.[28] inspected the improvement in thermal properties of conventional base fluids with interaction of magnetic nano-fluid subject to the dipole effects. Hassan et al.[29] visualized the shape factor in ferrofluid with dynamic of oscillating magnetic force. The thermal inspection in Maxwell nanofluid with external impact of heat generation was directed numerically by Majeed et al.[30]. Khan[31] discussed the entropy optimized flow of hybrid nanofluid over a stretched surface of rotating disk. The enhanced features of metallic nanoparticles subject to the magnetic dipole phenomenon were addressed by Majeed et al.[32].

In microscopic level the entropy rate is caused due to heat transfer, molecular vibration, dissipation, spin movement, molecular friction, kinetic energy Joule heating etc. and heat loss occurs. For improvement the productivity of numerous thermal schemes, it is necessary to optimize the irreversibility. Thermodynamic second law redirects more significant behaviors in comparison to thermodynamic first law. Thermodynamics second law gives the entropy optimization and scientific tools for decrease of confrontation. It helps us to develop the ability of various engineering improvements. These processes encompass heat conduction and furthermore to calculate the entropy generation rate. Primary attention of entropy generation problems is done by Bejan[33]. Zhou et al.[34] discussed irreversibility analysis about convective flow of nanoliquids in a cavity. Salient characteristics of thermophoretic and Brownian diffusion in flow of Prandtl-Eyring liquids with entropy optimization are exemplified by Khan et al.[35]. Irreversibility analysis in magnetohydrodynamic flow of Carreau nanofluids through Buongiorno nanofluid model is validated by Khan et al.[36]. Jiang and Zhou[37] studied viscous nanoliquid flow with irreversibility. Some advancement about irreversibility analysis is given in Refs.[38-45].

The above presented research work, it is observed that no determination has been completed to investigate the irreversibility consideration for convective viscous fluid flow over a rotating cone. Therefore intension in this paper is to scrutinize the irreversibility for mixed convection reactive flow of viscous fluid by a rotating cone. Heat transfer is demonstrated with heat generation/absorption and dissipation. Furthermore a physical characteristic of entropy is considered. Nonlinear governing system is altered to ODEs. The given system is than tackled through NDSolve procedure. Prominent characteristics of different engineering variable on velocity field, entropy rate, Bejan number, concentration and temperature are realistically examined. The computational outcomes of surface drag force, heat transfer rate and gradient of concentration are scrutinized via different remarkable parameters.

Formulation

We examine mixed convective flow of incompressible laminar fluid over a rotating cone. Angular velocity is denoted by . Energy expression with heat source/sink and dissipation is considered. Innovative behaviors regarding entropy optimization is accounted. First order chemical reaction is deliberated. The resistive force arises owing to variation in concentration and temperature in the liquid and flow is axi-symmetric. The acceleration associated with gravitational force are assumed along the downward direction. Figure 1 describes the physical model[9,10].

Figure 1

Sketch of problem[9,10].

The related expressions are[15,16]:withhere viscosity and conductivity are employed in the forms[44]here denotes the density, the constant viscosity, and the velocity components, the semi-vertical angle, the coefficient of solutal expansion, the variable viscosity parameter, the thermal coefficient expansion, , the temperature, the constant thermal conductivity, the specific heat, the wall temperature, the variable thermal conductivity parameter, the ambient temperature, the heat generation/absorption coefficient, the concentration, the dimensionless angular velocity, the ambient concentration the mass diffusivity, the wall concentration and the chemical reaction rate.

Lettingone haswhere shows the mixed convection parameter, the Reynold number, the Grashoff number, the buoyancy ratio variable, the Eckert number, the Prandtl number, the heat generation variable, the chemical reaction variable and the Schmidt number.

Entropy modeling

Mathematically entropy optimization is given by[41-43]:while after utilization of Eq. (11) yields[41-43]:

Bejan number is given as[41-43]:orin which signifies the entropy rate, the Brinkman number, the concentration ratio parameter, the temperature difference variable, dimensionless parameter and the diffusion variable.

Physical quantities

Velocity gradient

Surface drag forces are given aswith and as shear stresses are given by Finally we can write

Nusselt number

It is expressed aswith heat flux represented bynow

Mass transfer rate

Sherwood number iswith as mass flux through following expression Finally we have

Validation of results

Tables 1 and 2 are provided to authenticate the precision of current outcome with aforementioned published outcomes in literature. These tables deliberated the evaluation of velocity gradient and Nusselt number versus increasing values of with those of Saleem and Nadeem[34] and Chamka et al.[35]. These outcomes are established in good agreement.

Table 1

Comparison of surface drag force with Saleem and Nadeem[44] and Chamka et al.[45].

Pr	λ	Saleem and Nadeem[44]		Chamka et al.[45]		Recent results
		C_fxRex1/2	12C_fyRex1/2	C_fxRex1/2	12C_fyRex1/2	C_fxRex1/2	12C_fyRex1/2
0.7	0.0	1.0255	0.6154	1.0255	0.6158	1.0255	0.6156
	1.0	2.2010	0.8493	2.2012	0.8496	2.2010	0.8494
	10.0	8.5042	1.3992	8.5041	1.3995	8.5043	1.3992
10.0	0.0	1.0255	0.6158	1.0256	0.6158	1.0256	0.6158
	1.0	1.5630	0.6835	1.5636	0.6837	1.5631	0.6835
	10.0	5.0820	0.9845	5.0821	0.9840	5.0822	0.9842

Table 2

Comparison of Nusselt number with Saleem and Nadeem[44] and Chamka et al.[45].

Pr	λ	Saleem and Nadeem[44]	Chamka et al.[45]	Recent results
0.7	0.0	0.4299	0.4299	0.4298
	1.0	0.6121	0.6120	0.6122
	10.0	1.3992	1.0097	1.3993
10.0	0.0	1.4111	1.4110	1.4119
	1.0	1.5661	1.5662	1.5664
	10.0	2.3581	2.3580	2.3583

Physical description

Noticeable performances of various sundry variables about entropy rate, temperature, velocity field, Bejan number and concentration and are deliberated through graphs. Velocity gradient and Nusselt and Sherwood numbers are numerically computed against various parameters. The analysis is performed for flow parameters with specified numerical values range like and

Velocity

Salient effects of , and on (tangential velocity) and (azimuthal velocity) are examined in Figs. 2, 3, 4, 5, 6 and 7. Figure 2 depicts characteristics of tangential velocity for viscosity parameter . For increasing values of an enhancement occurs in . Characteristic of on is exposed in Fig. 3. Clearly is a decaying function of viscosity parameter . In fact increments in leads to reduction in temperature difference (convective potential) between ambient fluid heated surface and as a result azimuthal velocity decays. Figures 4 and 5 scrutinize the behaviors of on (tangential velocity) and (azimuthal velocity). One can find that and have reverse effects via larger . In fact augmentation in makes the fluid viscous and consequently decreases. Characteristics of on and are demonstrated in Figs. 6 and 7. These figures demonstrates that higher estimation of improves the tangential velocity , while reverse effect holds for azimuthal velocity .

Figure 2

against A.

Figure 3

against A.

Figure 4

against N.

Figure 5

against N.

Figure 6

against .

Figure 7

against .

Temperature

Figures 8, 9, 10, 11 and 12 have been displayed to explore behavior of pertinent variables like , , , and on . Figure 8 studied effect of viscosity variable on . Clearly temperature is a decreasing function of . Outcome of (Br) on temperature is sketched in Fig. 9. Here the increasing values of corresponds to an augmentation in . For larger Brinkman number the slower heat transmission is produced by viscous force and therefore boosts up. Figure 10 interprets the behaviors of on temperature. We noted that temperature improves through . Variation of on is interpreted in Fig. 11. Temperature against rises. Figure 12 is devoted to see the outcome of on . Clearly larger the thermal layer reduces which improves and heat transfer rate improves. Therefore decays.

Figure 8

against A.

Figure 9

against Br.

Figure 10

against .

Figure 11

against .

Figure 12

against .

Concentration

Impact of on is plotted in Fig. 13. Through Schmidth number, the concentration decays. Figure 14 is depicts the characteristics of on concentration . Clearly is diminished for higher estimation of . The fluid acts thick for higher and so reduction in occurs.

Figure 13

against .

Figure 14

against .

Entropy and Bejan number

Figures 15, 16, 17, 18, 19, 20, 21 and 22 are devoted to scrutinize the behaviors of various interesting parameter like viscosity parameter , thermal conductivity parameter , diffusion parameter and Brinkman number on and . Figures 15 and 16 are depicted to explore the effect of on and . Here and have opposite impact for increasing values of . Variation of on and is shown in Figs. 17 and 18. Clearly increasing values of give rise to both the and . Figures 19 and 20 are devoted to see the behavior of on and . Clearly for larger both and have increasing behaviors. Figures 21 and 22 display impact of and for Brinkman number . Larger Brinkman number rises the entropy generation. Figure 22 shows that for rising values of the decays.

Figure 15

against A.

Figure 16

against A.

Figure 17

against .

Figure 18

against .

Figure 19

against L.

Figure 20

against L.

Figure 21

against Br.

Figure 22

against Br.

Analysis for engineering quantities

Here impacts of various influencing variables on gradient of velocity ( and ) along azimuthal and tangential direction respectively, mass transfer rate and gradient of temperature are discussed in Tables 3, 4, and 5.

Table 3

Computational outcomes of .

λ	A	Surface drag force
		C_fx	C_fy
1	0.2	1.1345	0.46536
3		2.02356	0.76543
5		3.0145	1.45362
2	0.2	0.89654	0.80983
	0.4	0.75643	0.69954
	0.6	0.65874	0.56432

Table 4

Computational outcomes of .

Br	Pr	A	δ	Nu_x
0.0	0.7	2.0	0.2	0.6126
0.5				0.5325
1.0				0.4765
1.0	1.0	2.0	0.2	1.6875
	3.0			1.7894
	5.0			1.9283
1.0	0.7	1.0	0.2	0.7865
		2.0		0.6923
		3.0		0.6198
1.0	0.7	2.0	0.2	0.7967
			0.4	0.6987
			0.6	0.6089

Table 5

Numerical value of .

N	γ	Sc	Sh_x
0.0	0.4	0.1	0.56796
0.5			0.60289
1.0			0.64156
0.5	0.1	0.1	0.45342
	0.3		0.49786
	0.5		0.53675
0.5	0.4	0.2	0.5745
		0.5	0.6745
		0.8	0.7981

Velocity gradient

The numerical results of via various interesting parameters like viscosity parameter and mixed convection parameter are analyzed in Table 3. Clearly one can find that an increment occurs in via increasing values of . From this table it is noted that for larger estimation of viscosity variable the are decreased.

Temperature gradient

Influences of different sundry variables like , , and on is scrutinized in Table 4. Nusselt number in enhanced for larger and . Further is decreased for higher viscosity parameter and thermal conductivity parameter .

Sherwood number

The computational outcomes of via various flow variables are studied in Table 5. Here has similar characteristics for larger and . We noticed that rises via .

Conclusions

The applications of entropy generation phenomenon in the convective transport of viscous nanofluid due to rotating cone have been addressed in presence of viscous dissipation and heat generation. The analysis is performed in presence of variable thermal conductivity and fluid viscosity. The key observations are given below.

The tangential velocity and azimuthal velocity have contradictory behavior for mixed convection parameters.

The applications of viscosity parameter show increasing effects on tangential velocity.

The tangential velocity boosts up via buoyancy ratio variable.

The nanofluid temperature is enhanced for larger heat generation variable it decreased for viscosity parameter.

The nanofluid concentration is decreased for higher values of chemical reaction variable and Schmidt number.

The entropy rate and Bejan number are enhanced for diffusion variable.

The entropy rate upsurges versus Brinkman number.

The entropy rate and Bejan number have reverse effects for viscosity parameter.

The wall shear force increase via higher mixed convection parameter.

The surface drag force is diminished against viscosity parameter as it is reversely related to the magnitude of drag force per unit area.

The Nusselt number is increased for larger Prandtl number.

Gradient of temperature versus Brinkman number decreases.