Abstract

In the present study, the effect of thermal stratification and heat generation in the boundary layer flow of the Eyring-Powell fluid over the stratified extending surface due to convection has been investigated. The governing equations of the flow are transformed from partial differential equations into a couple of nonlinear ordinary differential equations via similarity variables. The optimal homotopy asymptotic method (OHAM) is used to acquire the approximate analytical solution to the problems. Impacts of flow regulatory parameters on temperature, velocity, skin friction coefficient, and Nusselt number are examined. It is discovered that the fluid velocity augments with a greater value of material parameter , mixed convection parameter , and material fluid parameter . The result also revealed that with a higher value of the Prandtl number Pr and the stratified parameter , the temperature and the velocity profile decreases, but the opposite behavior is observed when the heat generation/absorption parameter increases. The results are compared with available literature and are in good harmony. The present study has substantial ramifications in industrial, engineering, and technological applications, for instance, in designing various chemical processing equipment, distribution of temperature and moisture over agricultural fields, groves of fruit trees, environmental pollution, geothermal reservoirs, thermal insulation, enhanced oil recovery, and underground energy transport.

1. Introduction

The study of the dynamics of non-Newtonian fluids has inconceivable importance on account of its numerous biological, industrial, and engineering applications. They have a vital role in material production, mixture, paper coating, bioengineering, oil repository designing, petroleum production, food engineering, chemical engineering, and the study of electroosmatic flow, blood flow, etc. [18]. Non-Newtonian fluid can be found in numerous materials for example in shampoos, paints, cosmetic productions, body fluids, grease, and animal blood to name a few. Several models including the Eyring-Powell model (one among the non-Newtonian fluid models) have been proposed by different scientists for the analysis of dynamics of non-Newtonian fluids. The Eyring-Powell model is derived from the kinetic molecular theory of liquids rather than empirical relation and properly shows the Newtonian behavior at low and high shear rates; owing to this, it gained the attention of many scholars. This class of fluid does not have constant viscosity with temperature alterations [913].

The Eyring-Powell fluid model has essential applications in designing various chemical processing equipment, distribution of temperature and moisture over agricultural fields, groves of fruit trees, environmental pollution, geothermal reservoirs, thermal insulation, enhanced oil recovery, and underground energy transport [11]. The Eyring-Powell model with heat transfer plays a crucial role in diverse industrial, geophysical, and natural processes. The understanding of the heat transfer effect in a variety of materials and boundary layer flow overstretching sheets is also imperative as heat has substantial applications in diverse power station engineering, biological phenomena, and industrial branches like paper manufacture, metal extrusion, and bubble absorption [1418].

Heat and mass transfer in the Eyring-Powell liquid model undertakes an essential role in the procedures which comprise the creation and spread of haze, plotting of mixture handling instrumentation, environmental pollution, drying of porous slides, raised oil recuperation, warm protection, and underground energy transport [1921]. Ishaq et al. [22] explored the 2D nanofluid film stream of the Eyring-Powell liquid with changeable heat conduction in the existence of MHD over the unsteady permeable extending sheet and asserted that the porosity parameter lessens the movement of the fluid.

Hayat et al. [23] conducted a study on MHD nonlinear extending flow of the Eyring-Powell nanofluid, where a zero nanoparticle mass flux stipulation was applied at the sheet. Their finding illustrates that with magnetic parameters, skin friction coefficient raises though the heat transmits rate at the surface declines. The heat and mass transfer analysis of non-Newtonian fluid overstretching sheets has been done by different scholars; Megahed et al. [24] conducted a study on heat transfer for MHD fluid as a result of the unsteady stretching plate with extended heat flux by varying viscosity and thermal conductivity with temperature; Reddy et al. [25] and Megahed and Gnaneswara [26] have also investigated the heat and mass transfer of different non-Newtonian flow across a stretched sheet.

Mixed convection flow plays a significant role, for example, in air limit layer flows, atomic reactors, solar gatherers, and in electronic hardware, just to name a few. Wubshet and Temesgen [19] used Cattaneo–Christov heat and mass flux model to inspect the mixed convection flow of the Eyring-Powell nanofluid through a porous medium on a linearly extending sheet in the company of permeability and Joule heating effects. They noticed that with the advanced value of the Eyring-Powell fluid parameter, the fluid velocities, local Nusselt number, and skin friction coefficient promotes. Malik et al. [27] have examined a mixed convection flow for the MHD Eyring-Powell fluid over the extending plate solving the problem numerically by the shooting method. Their outcome indicates swelling of the Eyring-Powell material parameter results in a meaningful decrease in both heat and mass transfer. Moreover, Rahimi et al. [28] presented a collocation method solution for an Eyring-Powell fluid flow over a prolonging plate. They found that the velocity enhances with mounting the Eyring-Powell fluid material parameter, but the opposite behavior is observed when the fluid material parameter is augmented.

Stratification plays a vigorous role in several industrial and natural phenomena. Owing to its wide uses in fluid mechanics, several investigators studied the stratification phenomenon. Zubair et al. [29] investigated the attributes of a chemical reaction and twofold stratification in the Eyring-Powell fluid flow persuaded by a nonlinear prolonging plate with changeable thickness. They confirmed the temperature decline with thermal stratified parameters; Salahuddin et al. [30] inspected the vital aspects of the stratification procedure with linear convection. They discovered that with an enhanced value of the thermal stratification parameter the velocity, the temperature profile reduces. Hayat et al. [31] examined the MHD effect for heat and mass transfer through the double stratified sheet. They reported that thermal and solutal stratification parameters have a conflicting effect on temperature distributions. Furthermore, Ogunseye et al. [32] examined the effect of thermal radiation on the flow and heat transfer of the Eyring-Powell nanofluid. The understanding of the heat generation/absorption effect of boundary layer flow is also crucial in different fields of study and applications. Owing to this, a lot of scholars [3335] conducted a study on the heat generation effect on boundary layer flow and heat transfer due to different non-Newtonian fluids. The study of the boundary layer flow of the Eyring-Powell fluid over a stretching surface plays a vital role in industrial applications. A stretching surface has abundant significance in numerous industrial applications including electronic chips, fiber yarn, polymer industries, and glass blowing.

Motivated by the above significant studies, we consider the problem of the Eyring-Powell fluid flow over a stratified stretching surface subject to mixed convection and heat generation/absorption effect. The problem has plenty of applications in various fields of science and engineering but it has not been analytically investigated so far. Moreover, the effect of heat generation and thermal stratification on the flow of the fluid under consideration is not investigated. Hence, the novelties of the current work are the investigation of the heat generation and thermal stratification effects in the boundary layer flow of the Eyring-Powell fluid over the stratified extending surface due to convection, the use of OHAM for solving the problem, and the analysis of effects of different parameters on velocity and temperature of the considered fluid. The effectiveness and reliability of the proposed method have been proved in many research articles, and it is used for solving both nonlinear and linear problems by several authors [34, 3638].

2. Mathematical Formulation

In the model, a two-dimensional incompressible steady boundary layer mixed convection flow of the Eyring-Powell fluid over the stratified extending surface has been pondered. The sheet is elongated at invariable speed along the -axis in a vertical direction and the fluid temperature at and far away from the plate surface is given, respectively, by . The influence of heat generation/absorption is considered in the energy equation. The geometric representation of the problem under consideration and the coordinate system are displayed in Figure 1.

Applying the boundary layer estimation the governing equation of the problem is condensed into the following form (see [21, 39, 40]):

Continuity equation:

Momentum equation: where and , and , respectively, represent the horizontal and vertical components of fluid velocity and the coordinate system, and are the dimensional Eyring-Powell fluid material parameters, and represent kinematic viscosity and gravitational acceleration, respectively.

Energy equation:

The corresponding boundary conditions of the governing flow problem are (see [16, 31]) as follows: where is a linear stretching rate, and are positive dimensional constants, and is the stretching sheet reference temperature.

To make the arithmetical analysis simple, we convert Equations (1)–(4) into ODEs, by introducing the stream function , such that where

The dimensionless similarity variable , stream function , and temperature are defined as follows (see [21, 39, 40]): where .

Now, Equation (1) is identically pleased with and of Equation (5). Furthermore, using Equation (5) and (8), Equations (2) and (3) become converted into the following ODEs: where represent a material parameter, fluid parameter, mixed convection parameter, Prandtl number, heat generation/absorption parameter, and thermal stratification parameter, respectively. Moreover, the prescribed side conditions Equation (4) become reduced to

The physical quantities of interest: local skin friction coefficient and the Nusselt number , are described as follows (see [30, 40]): where represents free stream speed, and given by (see [21]) denote skin shear stress on a surface and heat flux at the surface.

From Equations (5), (8), (13), and (15), we obtain where denotes the Reynolds number and is the distance from the reference point at which a boundary layer begins to form.

3. Solution of the Problem via the OHAM

Employing the basic principle of OHAM on Equations (10) and (11), we construct the homotopy equation: where and are expanded as follows: where are nonzero ancillary functions for and zero for . Using Equation (18) in Equation (17) and collecting like terms based on powers of , we get the next zero-, first-, second-, and third-order problems with their boundary conditions.

The zeroth-order problem:

Boundary conditions:

The solution for Equation (19) with BCs (20) is

The 1st-order problem:

Boundary conditions:

Solution for Equation (22) with BC Equation (23) is

The 2nd-order problem:

Boundary conditions:

The solution for the problem is

The 3rd-order problem:

Boundary conditions:

Solution for Equation (29) subjected to BCs Equation (30) is achieved in the same way as the zeroth-, first-, and second-order problems. However, there are very bulky expressions to be included herein this paper. Collecting the solution of Equations (21), (24), (28), and (29), we get the four terms solution given as

4. Result and Discussion

The OHAM is used to solve the transformed nonlinear ODEs describing the problem. The effect of a variety of parameters on skin friction coefficient , Nusselt number Nu, temperature, and velocity has been investigated. The obtained results are displayed with the help of graphical illustrations and tables. Table 1 displays the obtained value of Nu and for particular values of certain parameters and a comparison has been made with published results of others’ work. From the table, it is detected that increases with swelling values of the Eyring-Powell fluid material parameter and mixed convection parameter , but it falls off with an augmented value of thermal stratification parameter and material fluid parameter . It is also found that Nu increases with enhancing mixed convection and Eyring-Powell fluid parameter .

The impact of the thermal stratification parameter and the Prandtl number Pr on Nusselt number Nu and the skin friction coefficient for the fixed value of the other parameters is available in Table 2. The table illustrates is a decreasing function of Prandtl number Pr and the thermal stratification parameter . Table 2 also indicates that the local Nusselt number Nu enhances with the increasing value of Pr and .

The effect of a variety of parameters on the fluid velocity profile is displayed in Figures 27. In Figure 2, the impact of mounting the material parameter on dimensionless fluid velocity is displayed. The figure points out that the velocity increase with augmenting value of . Physically, has an inverse relation with the non-Newtonian fluid’s dynamic viscosity. Hence, with increased value of , the flow resistance reduces, and consequently, the fluid velocity intensifies. Figure 3 is sketched to illustrate the influence of mixed convection parameter on velocity profile. It depicts that an intensification in results in the rise in the velocity profile. In reality, while fluids are heated their density reduces and buoyancy increases, as a consequence their velocity boosts. Thus, by advancing value , one can enhance the velocity within a boundary layer when required.

The influence of thermal stratification parameter on the fluid velocity is illustrated in Figure 4. It is witnessed that with higher value of , the velocity profile declines. Actually, when increases, the convective potential between the heated surface and ambient temperature reduces. The effect of heat generation parameter on velocity is exhibited through Figure 5. As can be seen from the figure, the velocity profile enlarges with the value of . In reality, while the heat generation parameter rises, additional heat is produced in the fluid; consequently, the fluid turn out to be hotter and less dense, and as a result, its viscosity decreases. Figure 6 displays the effect of material fluid parameter on the velocity profile. It is confirmed that the velocity graph increases with the mounting value of . In actual fact, is inversely related to its viscosity. So, the higher the fluid parameter leads to the deterioration of the viscous forces; for this reason, the fluid flows rapidly. Figure 7 describes the impact of the Pr on the velocity profile. It is detected that an increment in Pr corresponds to a decrease in the velocity profile of the fluid. Physically, Pr is directly related to the viscosity of the fluid. Because of this fact, as the Pr increases, the fluid turns out to be denser, which causes a decrease in the velocity profile. This may be accredited to the fact that convection currents became diminutive at higher Pr and leads to a decline in velocity.

The effect of a variety of parametric quantities on the nondimensional temperature distribution of the fluid under consideration is displayed in Figures 813. The characteristic of stratification parameter on the temperature is depicted in Figure 8. Physically, the rise in declines the temperature distinction between the fluid on the surface and the surrounding fluid, this in turn trim down the temperature profile. Thus, the temperature distribution in the fluid drops with an augmentation in the value of . Figure 9 reports the effect of the Prandtl number Pr on the temperature profile. It is revealed that when the value of Pr is augmented, the temperature graph reduces. Mathematically, Pr is inversely related to thermal conductivity; hence, the fluid with greater Pr has a low energy diffusion rate. So, the raise in the Pr produces a sturdy reduction in the temperature, which upshots the thermal boundary layer thickness. Therefore, in heat transfer, Pr is used to rheostat the thicknesses of thermal and momentum boundary layers.

Figure 10 is plotted to scrutinize the impact of heat generation/absorption parameter on the fluid temperature. The figure confirms that when the value intensifies, the temperature profile also increases. Physically, higher generates extra heat in the fluid system, due to this the temperature profile enhanced. Thus, when the heat generation parameter upshots, considerable heat which can increase the temperature profile of the fluid will be produced. Therefore, increasing the heat generation/absorption parameter is used to enhance the temperature of the fluid, but sometimes, the augmented values of heat generation parameter results in overmounting in the temperature boundary layer and it should be controlled by using suitable mechanisms.

The characteristic of material fluid parameter on the fluid temperature is displayed in Figure 11. It is witnessed that an increment in the value of corresponds to a turn down in the temperature profile. This behavior is revealed as is inversely associated to the viscosity of the fluid. So, when increases, viscosity reduces which in turn increases the heat transfer rate. Consequently, the temperature of the fluid reduces. Figure 12 portrays the behavior of material parameter on temperature profile. As we can observe from the figure, an improvement of value of results in a decrease in the fluid temperature graph. Physically, when is increased, fluid resistance reduces and contributes to the reduction of temperature distribution. Figure 13 depicts the impact of varied convection parameter on temperature. It has been found that when the value of increases, the temperature profile declines. This happens attributable to the reduction of the fluid density, as hotter fluids are less dense. Thus, with enhanced value of heat transfer rate increases, thereby temperature profile decreases.

5. Conclusion

In this work, we analyze heat transfer and flow of the incompressible, steady Eyring-Powell fluid flow over a two-dimensional stratified prolonging surface with mixed convection. The nonlinear PDEs governing the flow problem and their respective prescribed side conditions have been converted into a couple of nonlinear ODEs via variable similarity transformation and then solved analytically by the OHAM. The quantities of attention in this work are velocity , temperature , skin friction coefficient , and Nusselt number Nu. Thus, the effect of varied parameters embedded in the model on , , , and Nu has been examined and described with the aid of tables and graphs. The main observations drawn from this finding are the following: (i)With an intensification in value of the thermal stratification parameter , both the fluid velocity and the temperature profile decreases. Hence, in various manufacturing processes can help to control heat and velocity of fluids(ii)As the value of the Eyring-Powell fluid material parameter upsurges, the velocity profile boosts, while the temperature profile declines(iii)With snowballing values of the mixed convection parameter and the material fluid parameter , the velocity profile enhances(iv)As the value of Prandtl number Pr enlarges, the velocity profile dwindles(v)The temperature distribution turn down with enhancing values of the Prandtl number Pr, the material fluid parameter , and the mixed convection parameter (vi)As the value of the heat generation/absorption parameter amplifies, both the temperature and the velocity profile increase(vii)As the values of the mixed convection and the Eyring-Powell fluid material parameters heighten, both the Nusselt number Nu and the skin friction coefficient augment(viii)With mounting values of the thermal stratification and the material fluid parameters, the skin friction coefficient falloffs(ix)The local Nusselt number Nu is an increasing function of Prandtl number Pr and the thermal stratification parameter , the mixed convection parameter , and the Eyring-Powell fluid material parameter

Abbreviations

:Linear stretching rate, (1/s)
:Dimensional constants, (K/m)
:Dimensional Eyring-Powell fluid material parameter
:The local skin friction coefficient
:Specific heat of the fluid at a constant temperature, (J/kg·K)
:The thermal stratification parameter
:The dimensionless stream function
:Gravitational acceleration, (m/s2)
k:The fluid thermal conductivity, (W/m·K)
:Nusselt number
Pr:Prandtl number
:Heat generation/absorption coefficient
:Radiation heat flux at the plate surface, (W/m2)
:The fluid temperature, (K)
:Stretching sheet reference temperature, (K)
:The fluid temperature at the plate surface, (K)
:The temperature far away from the plate surface, (K)
:Velocity component along the -axis, (m/s)
:The velocity of the stretching sheet, (m/s)
:Velocity component along the -axis, (m/s)
:Vertical and horizontal coordinate system
:Dimensionless similarity variable
:Kinematic viscosity of the fluid, (m2/s)
:The fluid density, (kg/m3)
:Dimensional Eyring-Powell fluid material parameter
:Dimensionless temperature
:Dimensionless stream function
:Material parameter
:Mixed convection parameter
:Heat generation/absorption parameter, (mK2s3/kg)
:Skin shear stress (N/m2)
:Fluid parameter
:Free stream condition
:Condition at the plate surface.

Data Availability

The data that assist the findings of this work are integrated in the paper.

Conflicts of Interest

The authors declare that they have no conflict of interest.