# Passive and Active Control on 3D Convective Flow of Viscoelastic Nanofluid With Heat Generation and Convective Heating

^{1}Department of Mathematics, Dr. N.G.P. Arts and Science College, Coimbatore, India^{2}Department of Mathematics, Kongunadu Polytechnic College, Dindigul, India

The present article addresses the impact of passive control (PC) and active control (AC) on 3D flow of a viscoelastic nanofluid upon a stretching plate including heat generation and convective heating. The system of appearing non-linear PDE's are converted into a couple of ODE's by using suitable similarity transformations. Convergent series solutions are derived using the homotopy analysis method (HAM). Graphical results of velocity, nanoparticle volume fraction and temperature of different pertinent physical parameters with notable discussions are mentioned along with their physical significance.

## 1. Introduction

Fluids are much used to transfer heat in heat transfer equipment in industrial and engineering processes. Such processes are die casting, catalysis, distillation, synthesis in petrochemical industry, gas processing, hot-mix paving in concrete industry and steam generators in industrial laundry. In the above applications, thermal conductivity plays a vital role in heat transfer equipment. Conventional heat transfer of ordinary fluids, like water and oil, have low thermal conductivity, and poor heat transfer characters. The nanofluid is an advanced type of fluid incorporating nanometer-sized particles. This fluid is more stable, has high thermal conductivity, high mobility, a larger heat transfer surface between particles and fluids, low pumping power compared with ordinary fluids, low particle clogging, and low volume concentrations. Ahmed et al. (2014) addressed the uncertainties of dynamic viscosity and thermal conductivity of nanofluid flow in a permeable stretching tube under the influence of a heat sink/source. They proved that the skin friction coefficient of Ag–water nanofluid is smaller compared to the TiO_{2} nanofluid. The dual solutions of nanofluid flow past a moving semi-infinite flat plate are derived by Bachok et al. (2010), who found that a smaller heat transfer rate occurs in higher values of the Prandtl number. Khan and Pop (2010) explored the impact of Brownian motion and thermophoresis on a nanofluid flow in a flat surface. It is seen that Brownian motion suppresses the reduced Nusselt number and strengthens the reduced Sherwood number. The impact of Brownian motion and thermophoresis of nanofluid fluid flow between two parallel plates was recently investigated by Derakhshan et al. (2019). They proved that the Nusselt number is suppressed with enhanced thermophoretic and Brownian motion parameters. A few important works on this direction are (Hassani et al., 2011; Makinde and Aziz, 2011; Kasmani et al., 2015, 2016, 2017).

In recent years, many researchers have concerned themselves with the study of non-Newtonian materials because they have more usages in food processing, petroleum production, drawing of plastic films, metal spinning, etc. Viscoelastic fluid is one of the non-Newtonian fluids and contains viscous and elastic behaviors. Cortell (2006) discussed the MHD flow of a viscoelastic fluid past a stretching sheet with suction. He found that the viscoelastic parameter leads to an improvement of the heat transfer gradient. MHD viscoelastic boundary layer flow over a flat sheet with non-uniform heat source and radiation was analytically and numerically solved by (Abel and Mahesha, 2008). Eswaramoorthi et al. (2016) explored the effect of thermal radiation of a viscoelastic fluid over a stretching surface. They found that the velocity boundary layer thickness decays with a more enhanced viscoelastic parameter. All the above studies reported only viscoelastic fluid without nanofluid. However, in recent years, many researchers dealt with viscoelastic nanofluids because they have more industrial applications. Shit et al. (2016) examined the viscoelastic nanofluid over a stretching surface. Their result showed that the fluid temperature decays in the presence of the viscoelasticity of the nanofluid. Very recently, Hayat et al. (2018) studied the three-dimensional flow of viscoelastic nanofluid over a stretching surface with heat generation/absorption. They implemented the Cattaneo-Christov double-diffusion theory for their study. They observed that the skin friction coefficient enhances upon increasing the ratio parameter. Some useful investigations of this direction are (Ramzan and Yousaf, 2015; Ramzan et al., 2016; Seth et al., 2016). In some practical applications, such as cooling of nuclear reactors, underground disposal of radioactive waste material, and storage of foodstuffs, the heat absorption/generation effects are essential. MHD boundary layer flow of a viscoelastic nanofluid with source/sink was investigated by Goyal and Bhargava (2014) and they found that the thermal boundary layer thickens with large values for the heat source/sink parameter. Some recent discussions on heat absorption/generation can be seen in Azim et al. (2010), Eswaramoorthi et al. (2017), and Halim et al. (2017).

In nanofluids, the base fluid does not obey Newtonian fluid properties, so it becomes more justified to imagine them as viscoelastic fluids, as for example with Ethylene glycol-*ZnO*, Ethylene glycol-*CuO* and Ethylene glycol-*Al*_{2}*O*_{3}. This type of work can be explored on two- and three-dimensional viscoelastic fluids, but less research has been done for viscoelastic nanofluids, especially three-dimensional viscoelastic nanofluids. To overcome this, we investigated the impact of passive and active controls on 3D viscoelastic nanofluid flow upon a stretching plate with heat generation and convective heating.

## 2. Mathematical Formulation

Let us consider the time-independent 3D convective flow of a heat generating viscoelastic nanofluid upon a stretched plate. Assume that the nanofluid has a single phase as well as uniform shape and size. The thermophoresis and Brownian motion effects are considered in this model. The bottom side of the plate is heated with hot fluid of temperature *T*_{c} and makes a heat transfer coefficient *h*_{c} (see Eswaramoorthi et al., 2015, 2016). Under the above assumptions, the governing equations of the present model are given as (see Ramzan and Yousaf, 2015):

where, *u, v, w*, ν, *k*_{1}, *D*_{B}, *C, D*_{T}, *T, T*_{∞}, α_{m}, τ, *Q*, *ρ*, and *c*_{p} are the *x*− direction velocity, *y*− direction velocity, *z*− direction velocity, kinematic viscosity, material parameter of fluid, Brownian motion coefficient, fluid concentration, thermophoretic diffusion coefficient, temperature of the fluid, free stream temperature, thermal diffusivity, ratio of effective heat capacity of the nanoparticle material to heat capacity of the fluid, heat generation/absorption coefficient, density of the fluid and specific heat, respectively.

The corresponding boundary conditions are,

Now, we introduce the following dimensionless similarity variables,

Substituting Equation (7) in Equations (2-5), we have

Boundary conditions (6) in terms of *f*, *g*, ϕ and θ become:

where $K=\frac{{k}_{1}a}{\nu}$ is the viscoelastic parameter, $Pr=\frac{\nu}{{\alpha}_{m}}$ is the Prandtl number, $Nb=\frac{\tau {D}_{B}({C}_{w}-{C}_{\infty})}{\nu}$ is the Brownian motion parameter, $Nt=\frac{\tau {D}_{T}({T}_{c}-{T}_{\infty})}{\nu {T}_{\infty}}$ is the thermophoresis parameter $Hg=\frac{Q}{a\rho {c}_{p}}$ is the heat generation/absorption parameter, $Le=\frac{{\alpha}_{m}}{De}$ is the Lewis number, $c=\frac{b}{a}$ is the stretching ratio and $Bi=\frac{{h}_{c}\sqrt{\frac{v}{a}}}{k}$ is the Biot number.

The skin friction coefficients and Nusselt number are defined as follows:

Then, the dimensionless form of the skin friction coefficients (*C*_{fx} &* C*_{fy}) and Nusselt number (*Nu*) are defined as:

## 3. HAM Solutions

We define the initial assumptions of HAM as *f*^{0}(*η*) = 1 − *Exp*(−*η*), *g*^{0}(*η*) = *c*(1 − *Exp*(−*η*)), ${\varphi}^{0}(\eta )=-\frac{NtBiExp(-\eta )}{Nb(1+Bi)}$ (for PC), *ϕ*^{0}(*η*) = *Exp*(−*η*) (for AC) and ${\theta}_{0}(\eta )=\frac{BiExp(-\eta )}{1+Bi}$ and the linear operators are ${L}_{f}=\frac{{d}^{3}f}{d{\eta}^{3}}-\frac{df}{d\eta}$, ${L}_{g}=\frac{{d}^{3}g}{d{\eta}^{3}}-\frac{dg}{d\eta}$, ${L}_{\varphi}=\frac{{d}^{2}\varphi}{d{\eta}^{2}}-\varphi $ and ${L}_{\theta}=\frac{{d}^{2}\theta}{d{\eta}^{2}}-\theta $ with ${L}^{f}\left[{E}_{1}+{E}_{2}Exp(\eta )+{E}_{3}Exp(-\eta )\right]$, ${L}^{g}\left[{E}_{4}+{E}_{5}Exp(\eta )+{E}_{6}Exp(-\eta )\right]$, ${L}^{\varphi}\left[{E}_{7}Exp(\eta )+{E}_{8}Exp(-\eta )\right]$ and ${L}^{\theta}\left[{E}_{9}Exp(\eta )+{E}_{10}Exp(-\eta )\right]$, where *E*_{j}(*j* = 1 − 10) are constants. After simplifying the m^{th} order HAM equations, we get the following:

where *f*^{m}^{+}(*η*), *g*^{m}^{+}(*η*), *ϕ*^{m}^{+}(*η*) and *θ*^{m}^{+}(*η*) are the particular solutions and these solutions have the axillary parameters *h*_{f}, *h*_{g}, *h*_{ϕ} and *h*_{θ}. These parameters are important in the convergence of the HAM solutions. The range values are −2.0 ≤ *h*_{f}, *h*_{g} ≤ −0.2, −1.7 ≤ *h*_{ϕ} ≤ 0.1, (for PC), −1.4 ≤ *h*_{ϕ} ≤ −0.2 (for AC) and −1.5 ≤ *h*_{θ} ≤ 0.1, see Figure 1. We found that the *h* value for the whole region of *η* is -1.

**Figure 1**. *h* curves of *f*″(0) &* g*″(0) **(A)** and θ′(0) &* ϕ*′(0) for PC **(B)** and AC **(C)** with *K* = 0.1, *c* = 0.5, *Nb* = 0.5, *Nt* = 0.2, *Hg* = −0.3 and *Bi* = 0.5.

Table 1 represents the different order of −*f*″(0), −*g*″(0), *ϕ*′(0) and −*θ*′(0). It is clear the that 10^{th} order is adequate for both velocity profiles and the 25^{th} order is needed for both nanoparticle volume fraction and temperature profiles. Table 2 presents the comparison of −*f*″(0) and −*g*″(0) for various values of *c* (Qayyum et al., 2014) and found that our results are in excellent agreement.

**Table 2**. Comparison of −*f*″(0) and −*g*″(0) for different values of *c* with (Qayyum et al., 2014).

## 4. Correlation Analysis

For the design of the thermal system and to analyze their performance, the correlation equations are essential. Linear regression analysis is used to derive the correlation equations from the obtained numerical values. The correlation equations of the skin friction coefficients and Nusselt number for both controls are given by:

This equation is valid for 0 ≤ *K* ≤ 0.3, 0 ≤ *c* ≤ 1, −0.3 ≤ *Hg* ≤ 0.3, 0.2 ≤ *Nt* ≤ 1.2 and 0.4 ≤ *Nb* ≤ 1.2 with a maximum error of 0.05.

## 5. Results and Discussion

The impacts of pertinent parameters on *x*− direction velocity profile (*f*′(*η*)), *y*− direction velocity profile (*g*′(*η*)), nanoparticle volume fraction (*ϕ*(*η*)), temperature profile (*θ*(*η*)) and local Nusselt number ($Nu/\sqrt{Re}$) with fixed value of Prandtl number (*Pr* = 1.2) and Lewis number (*Le* = 1) are delineated in this section. Figures 2A,B portray the impact of viscoelastic parameter (*K*) on (*f*′(*η*)) and (*g*′(*η*)). In general, viscoelasticity generates tensile stress that opposes the fluid motion, and hence both velocities reduce when the values of *K* are enhanced. The consequences of thermophoresis parameter (*Nt*) and Biot number (*Bi*) on *ϕ*(*η*) for both controls are plotted in Figures 3A–D, and it was found that ϕ(η) is an enhancing function of both *Nt* and *Bi* for both controls. Physically, *Nt* creates a thermal conductivity of the nanoparticles. Larger thermal conductivity reveals more concentration between the nanoparticles. So, increasing the values of *Nt* would enrich the nanoparticle volume fraction boundary thickness. In addition, nanoparticle volume fraction is an increasing function of the Biot number. The nanoparticle volume fraction depends on the fluid temperature. A higher Biot number causes higher fluid temperature and this causes the nanoparticle volume fraction to be enhanced. It can be seen in Figures 4A,B that the nanoparticle volume fraction increased along with the rising heat absorption/generation parameter (*Hg*) and it dropped when the values of the Brownian motion parameter (*Nb*) increased. This is an important result for thermal engineering applications. Physically, larger values of *Nb* give more kinetic energy, and this energy increases the collision between the nanoparticles. This causes the reduction of nanoparticle volume fraction layer. The impacts of *Hg* and *Bi* on *θ*(*η*) profile for both controls are illustrated in Figures 5A–D. The positive values of *Hg* generate heat energy inside the boundary and this causes fluid temperature to increase. On the contrary, the negative values of *Hg* absorb heat energy inside the boundary and this causes a reduction of the fluid temperature. Physically, the large Biot number increases the heat transfer coefficient, and this creates a higher surface temperature. So, a larger Biot number escalates the thickness of the thermal boundary layer. This result is used to enhance heat transfer characteristics in the petroleum industry, thermal energy storage and heat exchangers. The changes of the local Nusselt number for various combinations of *Hg, K* and *Nt* are illustrated in Figures 6A–D for both controls. These figures shows that the heat transfer gradient is reduced when the values of *Hg, K* and *Nt* are high.

**Figure 2**. *x*− direction velocity **(A)** and *y*− direction velocity **(A)** profiles for different values of viscoelastic parameter with *c* = 0.5.

**Figure 3**. Nanoparticle volume fraction profile of different values of *Nt* **(A,B)** and *Bi* **(C,D)** for both passive and active cases with *K* = 0.1, *c* = 0.5, *Nb* = 0.5 and *Hg* = −0.3.

**Figure 4**. Nanoparticle volume fraction profile of different values of *Nb* **(A)** and *Hg* **(B)** for passive case with *K* = 0.1, *c* = 0.5, *Nt* = 0.2 and *Bi* = 0.5.

**Figure 5**. Temperature profile for different values of *Hg* **(A,B)** and *Bi* **(C,D)** for both passive and active cases with *K* = 0.1, *c* = 0.5, *Nb* = 0.5 and *Nt* = 0.2.

**Figure 6**. Local Nusselt number for different combination of *Hg, K* **(A,B)** and *Nt* **(C,D)** with *K* = 0.1, *c* = 0.5 and *Bi* = 0.5.

## 6. Conclusions

In this work, we address the impact of passive control and active control on 3D viscoelastic nanofluid flow upon a stretching membrane with heat absorption/generation and convective heating. The system of appearing non-linear PDEs is converted into a couple of ODE's by using suitable similarity transformations. Convergent series solutions are obtained using the homotopy analysis method (HAM). The following findings are observed in our study.

• Both *x*− direction and *y*− direction velocities are weakened with increasing *K* values.

• The nanoparticle volume fraction is strengthened when the values of *Nt* and *Bi* are high for both controls.

• The fluid temperature improves when the values of *Hg* and *Bi* increase for both controls.

• The heat transfer gradient drops when the values of *Hg, K* and *Nt* rise for both controls.

## Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

## Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## References

Abel, M. S., and Mahesha, N. (2008). Heat transfer in MHD viscoelastic fluid flow over a stretching sheet with variable thermal conductivity, non-uniform heat source and radiation. *Appl. Math. Model*. 32, 1965–1983. doi: 10.1016/j.apm.2007.06.038

Ahmed, S. E., Hussein, A. K., Mohammed, H. A., and Sivasankaran, S. (2014). Boundary layer flow and heat transfer due to permeable stretching tube in the presence of heat source/sink utilizing nanofluids. *Appl. Math. Comp.* 238, 149–162. doi: 10.1016/j.amc.2014.03.106

Azim, M. A., Mamun, A. A., and Rahman, M. M. (2010). Viscous Joule heating MHD-conjugate heat transfer for a vertical flat plate in the presence of heat generation. *Int. Commu. Heat Mass Transf.* 37, 666–674. doi: 10.1016/j.icheatmasstransfer.2010.02.002

Bachok, N., Ishak, A., and Pop, I. (2010). Boundary-layer flow of nanofluids over a moving surface in a flowing fluid. *Int. J. Ther. Sci.* 49, 1663–1668. doi: 10.1016/j.ijthermalsci.2010.01.026

Cortell, R. (2006). Effects of viscous dissipation and work done by deformation on the MHD flow and heat transfer of a viscoelastic fluid over a stretching sheet. *Physics Letters A*. 357, 298–305. doi: 10.1016/j.physleta.2006.04.051

Derakhshan, R., Ahmadreza Shojaei, Hosseinzadeh, K. H., Nimafar, M., and Ganji, D. D. (2019). Hydrothermal analysis of magneto hydrodynamic nanofluid flow between two parallel by AGM. *Case Studies Ther. Engin*. 14, 100439. doi: 10.1016/j.csite.2019.100439

Eswaramoorthi, S., Bhuvaneswari, M., Sivasankaran, S., and Rajan, S. (2015). Effect of radiation on MHD convective flow and heat transfer of a viscoelastic fluid over a stretching surface. *Procedia Eng.* 127, 916–923. doi: 10.1016/j.proeng.2015.11.364

Eswaramoorthi, S., Bhuvaneswari, M., Sivasankaran, S., and Rajan, S. (2016). Soret and Dufour effects on viscoelastic boundary layer flow, heat and mass transfer in a stretching surface with convective boundary condition in the presence of radiation and chemical reaction. *Scientia Iranica Transactions B: Mech. Engin*. 23, 2575–2586. doi: 10.24200/sci.2016.3967

Eswaramoorthi, S., Bhuvaneswari, M., Sivasankaran, S., and Rajan, S. (2017). Effect of partial slip and chemical reaction on convection of a viscoelastic fluid over a stretching surface with Cattaneo-Christov heat flux model. *IOP Conf. Ser. Mater. Sci. Eng.* 263, 1–8. doi: 10.1088/1757-899X/263/6/062009

Goyal, M., and Bhargava, R. (2014). Boundary layer flow and heat transfer of viscoelastic nanofluids past a stretching sheet with partial slip conditions. *Appl. Nanosci.* 4, 761–767. doi: 10.1007/s13204-013-0254-5

Halim, N. A., Sivasankaran, S., and Noor, N. F. M. (2017). Active and passive controls of the Williamson stagnation nanofluid flow over a stretching/shrinking surface. *Neural Comput. Appli.* 28, 1023–1033. doi: 10.1007/s00521-016-2380-y

Hassani, M., Mohammad Tabar, M., Nemati, H., Domairry, G., and Noori, F. (2011). An analytical solution for boundary layer flow of a nanofluid past a stretching sheet. *Int. J. Ther. Sci.* 50, 2256–2263. doi: 10.1016/j.ijthermalsci.2011.05.015

Hayat, T., Qayyum, S., Shehzad, S. A., and Alsaedi, A. (2018). Cattaneo-Christov double-diffusion theory for three-dimensional flow of viscoelastic nanofluid with the effect of heat generation/absorption. *Res. Phys*. 8, 489–495. doi: 10.1016/j.rinp.2017.12.060

Kasmani, R. M., Sivasankaran, S., Bhuvaneswari, M., and Hussein, A. K. (2017). Analytical and numerical study on convection of nanofluid past a moving wedge with Soret and Dufour effects. *Int. J. Numer. Methods Heat Fluid Flow*. 27, 2333–2354. doi: 10.1108/HFF-07-2016-0277

Kasmani, R. M., Sivasankaran, S., Bhuvaneswari, M., and Siri, Z. (2015). Effect of thermal radiation and suction on convective heat transfer of nanofluid along a wedge in the presence of heat generation/absorption. *AIP Conf. Proc.* 1682, 020008-1–020008-7. doi: 10.1063/1.4932417

Kasmani, R. M., Sivasankaran, S., Bhuvaneswari, M., and Siri, Z. (2016). Effect of chemical reaction on convective heat transfer of boundary layer flow in nanofluid over a wedge with heat generation/absorption and suction. *J. Appl. Fluid Mech.* 3, 379–388. doi: 10.18869/acadpub.jafm.68.224.24151

Khan, W. A., and Pop, I. (2010). Boundary-layer flow of a nanofluid past a stretching sheet. *Int. J. Heat Mass Transf.* 53, 2477–2483. doi: 10.1016/j.ijheatmasstransfer.2010.01.032

Makinde, O. D., and Aziz, A. (2011). Boundary layer flow of a nanofluid past a stretching sheet with a convective boundary condition. *Int. J. Ther. Sci.* 50, 1326–1332. doi: 10.1016/j.ijthermalsci.2011.02.019

Qayyum, A., Hayat, T., Alhuthali, M. S., and Malaikah, H. M. (2014). Newtonian heating effects in three-dimensional flow of viscoelastic fluid. *IChin. Phys. B*. 23, 054703. doi: 10.1088/1674-1056/23/5/054703

Ramzan, M., and Yousaf, F. (2015). Boundary layer flow of three-dimensional viscoelastic nanofluid past a bi-directional stretching sheet with Newtonian heating. *AIP Advances* 5, 057132-1-057132-15. doi: 10.1063/1.4921312

Ramzan, M., Yousaf, F., Farooq, M., and Chung, J. D. (2016). Mixed convective viscoelastic nanofluid flow past a porous media with Soret-Dufour effects. *Commu. Theoret. Phys.* 66, 133–142. doi: 10.1088/0253-6102/66/1/133.

Seth, G. S., Mishra, M. K., and Chamkha, A. J. (2016). Hydromagnetic convective flow of viscoelastic nanofluid with convective boundary condition over an inclined stretching sheet. *J. Nanofluids*. 5, 511–521. doi: 10.1166/jon.2016.1249

Keywords: viscoelastic nanofluid, heat generation, passive/active control, convective heating, convection, stretching plate, HAM

Citation: Eswaramoorthi S and Bhuvaneswari M (2019) Passive and Active Control on 3D Convective Flow of Viscoelastic Nanofluid With Heat Generation and Convective Heating. *Front. Mech. Eng.* 5:36. doi: 10.3389/fmech.2019.00036

Received: 04 October 2018; Accepted: 31 May 2019;

Published: 20 June 2019.

Edited by:

Dipankar Chatterjee, Central Mechanical Engineering Research Institute (CSIR), IndiaCopyright © 2019 Eswaramoorthi and Bhuvaneswari. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: M. Bhuvaneswari, msubhuvana@yahoo.com