A revised model for Jeffrey nanofluid subject to convective condition and heat generation/absorption.

Here magnetohydrodynamic (MHD) boundary layer flow of Jeffrey nanofluid by a nonlinear stretching surface is addressed. Heat generation/absorption and convective surface condition effects are considered. Novel features of Brownian motion and thermophoresis are present. A non-uniform applied magnetic field is employed. Boundary layer and small magnetic Reynolds number assumptions are employed in the formulation. A newly developed condition with zero nanoparticles mass flux is imposed. The resulting nonlinear systems are solved. Convergence domains are explicitly identified. Graphs are analyzed for the outcome of sundry variables. Further local Nusselt number is computed and discussed. It is observed that the effects of Hartman number on the temperature and concentration distributions are qualitatively similar. Both temperature and concentration distributions are enhanced for larger Hartman number.

1. Introduction

The mixture of ultrafine nanometer sized particles and a convectional heat transfer base fluid is known as nanofluid. These nanometer sized particles have different physical and chemical characteristics. Such particles have thermal conductivities remarkably higher than base liquids. The prime use of nanofluids is for thermal conductivity improvement. Nanofluids are significant in various applications including paper and printing, paints and coatings, power generation, drug delivery, cancer therapy, ceramics and food products etc. Further magneto nanofluids are quite prevalent in MHD pumps and accelerators, hyperthermia, cancer tumor treatment, sink float separation, wound treatment and several others. Choi and Eastman [1] proposed the word nanofluid. They concluded that insertion of metallic nanoparticles in the ordinary fluids can dramatically enhance the thermal conductivities and improve the heat transfer performance of these fluids. A model for convective transport in nanofluids was presented by Buongiorno [2]. He pointed out that Brownian diffusion and thermophoresis are the most important slip mechanisms. Boundary layer flow of nanofluid induced by a linear stretching surface was discussed by Khan and Pop [3]. Turkyilmazoglu and Pop [4] examined unsteady natural convection flow of nanofluids by a vertical flat plate with radiation. Double stratification effect in boundary-layer flow of nanofluid by a vertical plate is reported by Ibrahim and Makinde [5]. Further relevant studies involving nanofluids can be seen through the investigations [6-25] and various studies therein.

The study of boundary layer flow caused by a stretchable surface is relevant in numerous industrial and engineering utilizations. Such applications include drawing of copper wires, condensation process, die forging and extrusion of polymer in melt spinning, polymer extrusion, continuous stretching of plastic films, metal extrusion, paper production and fiber production etc. It is noted that stretching of surface is not linear in all the cases. The stretching surface may be nonlinear. Gupta and Gupta [26] declared that the stretching of surface is not linear in plastic process. Vajravelu [27] addressed two-dimensional flow of viscous fluid over a nonlinear stretching surface. Here power law surface velocity distribution u(x) = cx is considered. Cortell [28] analyzed heat transfer in the flow past a nonlinear stretching surface. Here two different thermal boundary conditions on the surface namely constant surface temperature and prescribed surface temperature are employed. The boundary layer flow of viscous fluid induced by a nonlinear stretching surface with thermal radiation and viscous dissipation effects is addressed by Cortell [29]. Hydromagnetic flow generated by a nonlinear stretching surface through modified Adomian decomposition method and Pade approximation technique is demonstrated by Hayat et al. [30]. Rana and Bhargava [31] studied flow of nanofluid over a nonlinear stretching surface with heat transfer. Mukhopadhyay [32] addressed the flow and heat transfer characteristics in the flow of nanofluid over a permeable nonlinear stretching surface with partial slip condition. Mustafa et al. [33] explored axisymmetric flow of nanofluid over a nonlinear stretching surface. Magnetohydrodynamic flow of water-based nanofluid bounded by a nonlinear stretching surface with viscous dissipation is analyzed by Mabood et al. [34]. Magnetohydrodynamic flow of second grade nanofluid over a nonlinear stretching surface is reported by Hayat et al. [35].

Recently the non-Newtonian fluids have gained much attention due to their extensive industrial and engineering applications. These applications involve bioengineering and polymeric liquids, plastics manufacturing, food processing, petroleum production, annealing and thinning of copper wires, drawing of stretching sheet through quiescent fluid, aerodynamic extrusion of plastic films etc. The Navier-Stokes equations are not appropriate to characterize the flow of non-Newtonian fluids. A single relation is not sufficient to predict the characteristics of all the non-Newtonian materials. Therefore different types of relations are given in the literature. The fluid model under discussion is called Jeffrey material [36-41]. This model is linear viscoelastic fluid which exhibits the effects of ratio of relaxation to retardation times and retardation time. The Jeffrey fluid is a relatively simpler linear model considering time derivatives while in non-Newtonian fluid mechanics convective derivatives are assumed. Further the analysis of liquid-liquid two-phase flows are widely encountered in several industrial processes such as spray processes, lubrication, natural gas networks, nuclear reactor cooling etc. Thus Gao et al. [42] provided a multivariate weighted complex network analysis to characterize the nonlinear dynamic behavior in two-phase flow. Gao et al. [43] also addressed the multi-frequency complex network to uncover oil-water flows. Slug to churn flow transition with multivariate pseudo Wigner distribution and multivariate multiscale entropy is reported by Gao et al. [44]. Recently Gao et al. [45] provided a four-sector conductance method to explore the low-velocity oil-water two-phase flows.

Present communication explores magnetohydrodynamic (MHD) boundary-layer flow of Jeffrey nanofluid over a nonlinear stretching surface. Jeffrey fluid is assumed to be electrically conducting. We considered the simultaneous effects of heat and mass transfer in the presence of Brownian motion, thermophoresis and heat generation/absorption. Thermal convective [46, 47] and zero nanoparticles mass flux [48, 49] conditions are imposed at the stretching surface. These conditions are studied rarely and more realistic physically. To the best of our knowledge, no such consideration for the flow of Jeffrey nanofluid is made yet. Small magnetic Reynolds number and boundary layer are used in mathematical modelling. The governing nonlinear ordinary differential equations are solved by homotopy analysis method (HAM) [50-60]. This technique for the solutions development has advantages through three directions i.e., (i) It is independent of small/large physical parameters. (ii) It provides a simple way to ensure the convergence of series solutions. (iii) It provides freedom to choose the base functions and related auxiliary linear operators. Temperature and concentration profiles are examined via plots. The local Nusselt number is computed numerically and analyzed.

2. Statement

Two-dimensional (2D) flow of Jeffrey nanofluid induced by a surface stretching with nonlinear velocity is considered. Non-uniform magnetic field of strength B0 acts in the y− direction. Small magnetic Reynolds number justifies the absence of induced magnetic field. Non-uniform heat generation/absorption effect is considered. Brownian motion and thermophoresis are present. The x− and y− axes are along and transverse to the surface respectively. The stretching velocity is u(x) = ax (a, n > 0). The surface temperature is regulated by a convective heating process which is described by heat transfer coefficient h and temperature of hot fluid T under the surface. Resulting boundary layer problems are

Note that u and v depict the flow velocities in the horizontal and vertical directions respectively while v (= μ / ρ), μ and ρ show the kinematic viscosity, dynamic viscosity and density of base liquid respectively. The ratio of relaxation to retardation times and the retardation time are represented by λ1 and λ2. Here σ represents the electrical conductivity, the non-uniform magnetic field, T the temperature, α = k / (ρc), k, (ρc) and (ρc) the thermal diffusivity, thermal conductivity, heat capacity of liquid and effective heat capacity of nanoparticles respectively, Q(x) = Q0x the non-uniform heat generation/absorption coefficient, D the Brownian diffusivity, C the concentration, D the thermophoretic diffusion coefficient, a the positive constant and T∞ and C∞ the ambient fluid temperature and concentration respectively. Putting Eq (1) is trivially satisfied while Eqs (2)–(6) are reduced to

Here β1 denotes local Deborah number, Ha Hartman number, γ Biot number, S1 heat generation/absorption parameter, Pr Prandtl number, N Brownian motion parameter, N thermophoresis parameter and Le Lewis number. The definitions of these parameters are

Local Nusselt number is given by

The non-dimensional local Sherwood number is identically zero and Re = ux/v represents the local Reynolds number.

3. Solutions by HAM

The appropriate initial approximations and auxiliary linear operators for approximate series solutions by homotopy analysis method (HAM) are subject to where (j = 1–7) denote the arbitrary constants. Deformation problems at zeroth-order are

Here þ ∈[0,1] represents the embedding parameter, ℏ, ℏ and ℏ the non-zero auxiliary parameters and N, N and Nϕ the nonlinear operators. For þ = 0 and þ = 1 we have

When þ changes from 0 to 1 then , and vary from primary approximations f0(ζ), θ0(ζ) and ϕ0(ζ) to the desired solutions f(ζ), θ(ζ) and ϕ(ζ). The following expressions by Taylor's series expansion can be written as

The convergence of above series expressions strongly depends upon ℏ ℏ and ℏ. The values of ℏ ℏ and ℏ are chosen so that Eqs (28)–(30) converge at þ = 1 then

The th-order deformation problems can be expressed as follows:

General expressions of through special solutions are presented by the following expressions: in which the constants (j = 1–7) through the boundary condition (37) are given by

4. Convergence analysis

The expressions (31)–(33) contain ℏ ℏ and ℏ. Obviously the convergence is accelerated by the auxiliary parameters ℏ ℏ and ℏ for the series solutions. For appropriate values of ℏ ℏ and ℏ, the ℏ− curves at 15th order of approximations are sketched. It is apparent from Fig 1 that the admissible ranges of ℏ ℏ and ℏ are −1.35 ≤ ℏ ≤ −0.15, −1.50 ≤ ℏ ≤ −0.15 and −1.60 ≤ ℏ ≤ −0.15 respectively. The residual errors for velocity, temperature and concentration distributions are calculated through the following expressions:

Fig 1

The ℏ − curves for f, θ and ϕ.

To get the suitable range for ℏ, the ℏ− curves for the residual errors of velocity, temperature and concentration distributions are plotted in the Figs 2–4. It is observed that the correct results up to fifth decimal place are obtained for values of ℏ from this range. Table 1 presents that the 24th order of deformations is enough for the convergent series solutions of velocity, temperature and concentration distributions.

Fig 2

ℏ − curve for the residual error .

Fig 4

ℏ − curve for the residual error .

Table 1

Homotopic solutions convergence when n = 1.5, β1 = Ha = 0.2, N = 0.5, λ1 = γ = N = 0.3, Pr = 1.2, S1 = 0.1, and Le = 1.0.

Order of approximations	−f′′(0)	−θ′(0)	−ϕ′(0)
1	1.05680	0.21543	0.12926
8	1.04981	0.19237	0.11542
15	1.04981	0.19085	0.11451
24	1.04981	0.19065	0.11439
30	1.04981	0.19065	0.11439
40	1.04981	0.19065	0.11439
50	1.04981	0.19065	0.11439

5. Discussion

This portion organized the impacts of local Deborah number β1, Hartman number Ha, Brownian motion parameter N, ratio of relaxation to retardation times λ1, Biot number γ, thermophoresis parameter N, Prandtl number Pr, heat generation/absorption parameter S1 and Lewis number Le on the temperature θ(ζ) and concentration ϕ(ζ). Fig 5 illustrates that how local Deborah number β1 affects the temperature distribution θ(ζ). It is analyzed that temperature θ(ζ) and related thermal layer thickness are decreased for larger local Deborah number β1. Physically there exists a direct relationship between local Deborah number β1 and retardation time. Hence by increasing local Deborah number β1, the retardation time is also enhanced. Such enhancement in retardation time corresponds to lower temperature distribution θ(ζ) and thinner thermal layer thickness. Influence of λ1 on temperature distribution θ(ζ) is shown in Fig 6. For larger λ1, the relaxation time increases and retardation time decays. Thus temperature distribution θ(ζ) and thermal layer thickness are increased. Fig 7 presents variation in temperature distribution θ(ζ) for Hartman number Ha. An increase in Hartman number corresponds to more temperature θ(ζ) and thermal layer thickness. As expected the magnetic field introduces the retarding body force that acts transverse to the direction of an applied magnetic field. It retards the fluid motion and as a result the temperature distribution θ(ζ) enhances. This body force is known as Lorentz force. Fig 8 presents the impact of Biot number γ on temperature θ(ζ). Stronger convection is caused by increasing Biot number γ. Therefore the temperature θ(ζ) and thermal layer thickness are enhanced. Variation in temperature θ(ζ) due to heat generation/absorption parameter S1 is shown in Fig 9. Here S1 > 0 represents heat generation and S1 < 0 yields heat absorption. Temperature profile and related thermal layer thickness have increasing behavior for heat generation but it is not the case for heat absorption. Fig 10 demonstrates the variation of temperature θ(ζ) for Prandtl number Pr. It is observed that temperature θ(ζ) and thermal layer thickness are decreasing functions of Pr. Physically Prandtl number Pr is an integral part of thermal diffusivity. Thermal diffusivity is responsible for lower temperature θ(ζ) and thermal layer thickness. Higher values of Prandtl number yields weaker thermal diffusivity which corresponds to lower temperature and less thickness of thermal layer. Fig 11 is drawn for impact of thermophoresis parameter N on temperature θ(ζ). Larger thermophoresis parameter N lead to higher temperature and more thermal layer thickness. Actually an enhancement in N yields a stronger thermophoretic force which allows deeper migration of nanoparticles in the fluid. Far away from the surface there is higher temperature field and more thickness of thermal layer. Fig 12 is sketched to examine concentration field ϕ(ζ) for local Deborah number β1. Here concentration field is weaker for larger values of β1. Concentration field ϕ(ζ) enhances when λ1 increases (see Fig 13). Fig 14 shows impact of Hartman number Ha on concentration ϕ(ζ). The concentration ϕ(ζ) and associated layer thickness are enhanced for larger Hartman number. From Fig 15 we observed that an increase in Biot number γ yields an enhancement in concentration profile ϕ(ζ) and its related boundary layer thickness. Larger Lewis number Le indicate decay in the concentration field ϕ(ζ) (see Fig 16). Physically Lewis number is based on Brownian diffusivity. An increase in Lewis number Le yields weaker Brownian diffusivity. Such weaker Brownian diffusivity corresponds to lower concentration field ϕ(ζ). Fig 17 addresses variation of Prandtl number Pr on concentration ϕ(ζ). The concentration ϕ(ζ) and associated thickness of boundary layer are decreased for higher Prandtl number Pr. From Fig 18 it is clearly examined that a weaker concentration profile ϕ(ζ) is generated by higher Brownian motion parameter N. Fig 19 shows that the larger thermophoresis parameter N yields a higher concentration profile ϕ(ζ). Table 2 is calculated for numerical computations of local Nusselt number via β1, λ1, Ha, γ, S1, N, N Le and Pr when n = 1.5. Here we noticed that the local Nusselt number has higher values for larger Prandtl number Pr while opposite trend is noticed for Lewis number Le. It is also observed that λ1, S1 and Ha yield lower local Nusselt number. The local Deborah number β1 shows opposite behavior for local Nusselt number when compared with aforementioned parameters.

Fig 5

Plots of temperature θ(ζ) for local Deborah number β1.

Fig 6

Plots of temperature profile θ(ζ) for ratio of relaxation to retardation time λ1.

Fig 7

Plots of temperature profile θ(ζ) for Hartman number Ha.

Fig 8

Plots of temperature profile θ(ζ) for Biot number γ.

Fig 9

Plots of temperature profile θ(ζ) for heat generation/absorption parameter S1.

Fig 10

Plots of temperature profile θ(ζ) for Prandtk number Pr.

Fig 11

Plots of temperature profile θ(ζ) for thermophoresis parameter N

Fig 12

Plots of concentration profile ϕ(ζ) for local Deborah number β1.

Fig 13

Plots of concentration profile ϕ(ζ) for ratio of relaxation to retardation time λ1.

Fig 14

Plots of concentration profile ϕ(ζ) for Hartman number Ha.

Fig 15

Plots of concentration profile ϕ(ζ) for Boit number γ.

Fig 16

Plots of concentration profile ϕ(ζ) for Lewis number Le.

Fig 17

Plots of concentration profile ϕ(ζ) for Prandtl number Pr.

Fig 18

Plots of concentration profile ϕ(ζ) for Brownian motion parameter N.

Fig 19

Plots of concentration profile ϕ(ζ) for thermophoresis motion parameter N.

Table 2

Numerical calculations of local Nusselt number for different values of β1, λ1, Ha, γ, S1, N, N Le and Pr when n = 1.5.

β1	λ1	Ha	γ	S1	Nt	Nb	Le	Pr	Rex−1/2Nu_x
0.0	0.3	0.2	0.3	0.1	0.3	0.5	1.0	1.2	0.2060
0.3									0.2158
0.6									0.2210
0.2	0.0	0.2	0.3	0.1	0.3	0.5	1.0	1.2	0.2187
	0.5								0.2098
	1.0								0.2004
0.2	0.3	0.2	0.3	0.1	0.3	0.5	1.0	1.2	0.2133
		0.5							0.2088
		0.8							0.2003
0.2	0.3	0.2	0.1	0.1	0.3	0.5	1.0	1.2	0.0941
			0.6						0.3107
			1.2						0.4007
0.2	0.3	0.2	0.3	0.0	0.3	0.5	1.0	1.2	0.2282
				0.1					0.2133
				0.2					0.1839
0.2	0.3	0.2	0.3	0.1	0.0	0.5	1.0	1.2	0.2151
					0.5				0.2119
					1.0				0.2082
0.2	0.3	0.2	0.3	0.1	0.3	0.5		1.2	0.2133
						1.0	1.0		0.2133
						1.5			0.2133
0.2	0.3	0.2	0.3	0.1	0.3	0.5	0.5	1.2	0.2140
							1.0		0.2133
							1.5		0.2128
0.2	0.3	0.2	0.3	0.1	0.3	0.5	1.0	1.0	0.2017
								1.5	0.2260
								2.0	0.2403

6. Conclusions

Magnetohydrodynamic (MHD) flow of Jeffrey nanofluid bounded by a nonlinear stretching surface with heat generation/absorption is investigated. The observations are summarized in the following points.

An increase in local Deborah number β1 depicts a decreasing behavior for temperature θ(ζ) and concentration ϕ(ζ) profiles.

Both temperature θ(ζ) and concentration ϕ(ζ) profiles are enhanced when ratio of relaxation to retardation times λ1 is increased.

An increase in Hartman number Ha shows higher temperature θ(ζ) and concentration ϕ(ζ) profiles.

Biot number γ has similar effects for temperature θ(ζ) and concentration ϕ(ζ) profiles.

Prandtl number Pr indicates qualitatively similar behavior for both temperature θ(ζ) and concentration ϕ(ζ) profiles.

Temperature profile θ(ζ) and associated thermal layer thickness are increasing functions of heat generation/absorption parameter S1.

Concentration profile ϕ(ζ) decays for larger Brownian motion parameter N.

Increasing behavior is noted for temperature θ(ζ) and concentration ϕ(ζ) profiles for larger thermophoresis parameter N.

Local Nusselt number reduces for larger N but it remains constant for N.