Numerical values of nanoparticles and water.

## Abstract

A computational simulation for two-dimensional steady flow of modified nanofluid over an exponential stretching surface in a porous medium with magnet hydrodynamics and variable viscosity is presented in this study. Modified nanofluids are generalization of both hybrid nanofluids and simple nanofluids. Here, we consider three nanoparticles which drastically enhance the thermal conductivity of nanofluid. The viscous model associated with variable viscosity and MHD flow is employed. Well-known similarity transformations are utilized to convert the partial differential equations to system of ordinary differential equations. These converted equations are solved by utilizing the numerical technique Runge–Kutta-Fehlberg method. The impacts of variable viscosity, porosity parameter, Nusselt number, thermal and velocity slip, skin friction coefficient, solid nanoparticle, and magnetic field are observed. The computational results accomplished in the present investigation are validated and felt to be a good agreement with decayed results. It is highlighted that modified nanofluid model enhances the heat transfer rate much higher than the case of hybrid nanofluid and simple nanofluid model.

### Keywords

- variable viscosity
- exponential stretching
- modified nanofluid
- MHD
- porous medium
- shooting method

## 1. Introduction

Porous medium is one of the most useful studies due to its applications in the industry and medical sciences. In the medical sciences, it is used in the transport process in the human lungs and kidneys, gall bladder in the presence of stone, clogging in arteries, and also little blood vessels which cannot be opposed. There are several examples of the naturally porous medium such as limestone, wood, seepage of water in river beds, etc. Many researchers are interested to discuss the porous medium due to scientific and technically importance such as earth’s science and metallurgy. Such kinds of the flow are analyzed at low Reynolds number in the presence of porous space theoretically. Few researchers were analyzed analytically and experimentally on the porous medium with respect to different aspects (see [1, 2, 3]). Recently, the Carreau fluid flow over porous medium in the presence of pressure-dependent viscosity has been discussed by Malik et al. [4]. Some significant results are analyzed on the porous medium for Newtonian fluids and non-Newtonian fluids with respect to different aspects (see [5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]) (Figure 1).

In the depth study, flow phenomenon focusing on the variable viscosity and exponentially stretching surface is an important rule in the study of fluid mechanics and has attracted the investigators after its valuable applications in the industry as well as flows detected over the tip of submarine and aircrafts. Numerous methods have been established in recent past years to enhance the fluid thermal conductivity which is suspended with micro−/nano-sized particle mix with base fluid. The nanoparticle possesses chemical and physical properties uniquely because it has been used widely in nanotechnology. The nano-sized particle which is suspended with fluid is called nanofluid. Many investigators investigate about the enhancement of thermal conductivity [16, 17, 18, 19, 20] by using the nano-sized particles.

Several experiments have been done in two types of the particles suspended in the base fluid, namely, “hybrid nanofluid.” Basically, such type of fluids is enhances thermal conductivity which was proven through experimental research. Suresh et al. [21, 22] were the first to discuss the idea of hybrid nanofluid through their experimental and numerical results. According to their views, the hybrid nanofluid boosts the heat transfer rate at the surface as compared to nanofluid and simple fluid. These results open a new horizon to the researchers to do a work in the field of hybrid nanofluid. Baghbanzadeh et al. [23] also discussed about the mixture of multiwall/spherical silica nanotube hybrid nanostructures and analysis of thermal conductivity of associated nanofluid. The analysis of

The physical characteristics of hybrid nanofluid and nanofluid are usually considered constant. It is prominent that the significant physical characteristics of nanofluid and hybrid nanofluid can vary with temperature. For lubricating fluids, heat generated by the internal friction and the corresponding rise in temperature affects the viscosity of the fluid and so the fluid viscosity can no longer be assumed constant. The increase of temperature leads to a local increase in the transport phenomena by reducing the viscosity across the momentum boundary layer and so the heat transfer rate at the wall is also affected. The impact of thermal radiation and dependent viscosity of fluid on free convective and heat transfer past a porous stretching surface were discussed by Mukhopadhyay and Layek [29]. They gain some significant results for the variable viscosity on the temperature profile and velocity profile. The velocity profile increases and temperature profile decreases for large values of the variable viscosity parameter. The existing literature survey on the variable fluid characteristics and hybrid nanofluid [30, 31, 32, 33] reveals that the work is not carried out for hybrid nanofluid over an exponentially stretching surface.

The investigation about the stretching surface has attracted the interest of scientists because of its several applications in the fields of engineering including glass blowing, cooling of microelectronics, quenching in metal foundries, wire drawing, polymer extrusion, rapid spray, etc. Crane [34] discussed about the theoretical boundary layer flow over stretching surface. Various researchers analyzed the exponentially stretching surface [35, 36, 37, 38], major applications in the industry and technology.

Here, we study the temperature-dependent viscosity effects on the modified nanofluid flow over an exponentially stretching porous medium in the presence of MHD and Biot number. It is highlighted here that the idea of modified nanofluid has been proposed by us from whom the hybrid nanofluid and simple nanofluid cases can be recovered as a special case. The temperature depends on the Biot number, nanoparticle, and variable viscosity. The system of the flow is illustrated in the form of partial differential equations (PDEs). The system of PDEs is converted into the form of ordinary differential equations (ODEs) by utilizing acceptable similarity transformations. These nonlinear ODEs are solved “numerically” through MATLAB built-in technique. The outcomes are represented through table and graphs.

## 2. Flow formulation

Investigation of steady laminar flow of two-dimensional electrically conducting modified nanofluid over exponentially stretching surface in the presence of variable viscosity has been taken into consideration which is revealed in Figure 1.

The fluid flows in the *x*-direction and is maintains at a constant wall temperature

Thermophysical properties | Ni | |||
---|---|---|---|---|

444 | ||||

8900 | ||||

40 | 90.7 |

Under these assumptions with the usual boundary layer approximation, the governing differential equations of mass, momentum, and energy for the problem under consideration are defined as follows:

The appropriated boundary conditions are stated as

where

Properties | Nanofluid | Hybrid nano-fluid |
---|---|---|

Density | ||

Heat capacity | ||

Viscosity | ||

Thermal conductivity | where |

Properties | Modified nanofluid |
---|---|

Density | |

Heat capacity | |

Viscosity | |

Thermal conductivity |

An extraordinary type of physical characteristics is acquainted in the present examination to investigate the boundary layer equations for modified nanofluid. Modified nanofluid is deliberated through taking the combination of *Ni*,

i.e.,

The mathematical model over exponentially stretching surface is chosen to allow the coupled non-linear partial differential equations are converted into coupled non-linear ordinary differential equations by using the suitable similarity transformation which is given above. Where

with boundary conditions

## 3. Numerical solution method

Boundary layer heat transfer and modified nanofluid flow of an exponentially stretching surface with (

subject to the boundary conditions

For brevity, the points of interest of the solution strategy are not performed here. The heat transfer and modified nanofluid are effected by dependent viscosity parameter and MHD; the fundamental focus of our investigation is to bring out the impacts of these parameters by the numerical solution. It is worth specifying that we have utilized the information displayed in Tables 1–3 for the thermophysical properties of the fluid, nanofluid, hybrid nanofluid, modified nanofluid, and nanoparticles. Three types of the nanoparticles are used, namely, *Ni*. The Nussle number and skin friction coefficient are the most important features of this study. For practical purposes, the functions

and skin friction coefficient

Here, the local Reynolds number is

## 4. Numerical results

The impact of dependent viscosity parameter *f*″(0) and θ′(0) reveals that the same behavior to be noted for large values of

−10 | −1.98532 | −1.7779 | −1.90976 | −1.98809 |

−5 | −2.07138 | −1.77154 | −1.99315 | −1.98175 |

−1 | −2.65068 | −1.72622 | −2.55803 | −1.93639 |

−0.1 | −5.21601 | −1.49115 | −5.08404 | −1.69636 |

1 | −0.768213 | −1.85951 | −0.731972 | −2.06964 |

5 | −1.70487 | −1.79801 | −1.63866 | −2.00813 |

10 | −1.80225 | −1.79113 | −1.7327 | −2.00127 |

The computational results are shown in Table 5. The velocity of the flow decreases due to increase in the solid nanoparticle of

0.0 | −1.50718 | −1.86066 | −1.44127 | −2.07585 | |||

0.5 | −1.6881 | −1.83996 | −1.61326 | −2.05562 | |||

1.0 | −1.8533 | −1.82077 | −1.77026 | −2.03689 | |||

0.5 | 0.0 | −1.58519 | −2.79411 | −1.52063 | −3.06462 | ||

0.2 | −1.6881 | −1.83996 | −1.61326 | −2.05562 | |||

0.4 | −1.73712 | −1.37364 | −1.65867 | −1.54835 | |||

0.2 | 0.0 | −1.50718 | −1.86066 | −1.44127 | −2.07585 | ||

0.5 | −1.6881 | −1.83996 | −1.61326 | −2.05562 | |||

1.0 | −1.8533 | −1.82077 | −1.77026 | −2.03689 | |||

0.5 | 0.005 | −1.6158 | −1.99761 | −1.5041 | −2.23451 | ||

0.04 | −1.6881 | −1.83996 | −1.61326 | −2.05562 | |||

0.08 | −1.74614 | −1.67941 | −1.70227 | −1.87409 |

## 5. Graphical results

The temperature profile shows the variation of solid nanoparticle in Figure 3. The nanoparticle dissipates energy in the form of heat. So, the mixture of more nanoparticles may exert more energy which increases the thickness of the boundary layer and temperature.

Figure 4 reveals the impacts of solid particle on velocity profiles. The velocity profile gets decelerated due to increase in solid nanoparticle for modified nanofluid. This phenomenon exists due to more collision with suspended nanoparticles.

Figure 5 reveals the effects of magnetic field on the velocity profile. Being there, the transverse magnetic field creates Lorentz force which arises from the attraction of electric field and magnetic field during the motion of an electrically conducting fluid. The velocity profile decreases for the positive values of magnetic field parameter. Because the resisting force increases and consequently velocity declines in the

Figure 6 reveals the variation of dimensionless quantity of Biot number on the temperature profile. The relative transport of internal and external resistances is called the Biot number. The thermal boundary layer increases with increasing in the biot number.

Figure 7 shows the impact of the porosity parameter on the velocity profile. It is noted that velocity profiles decreases for the higher values of the porosity parameter. The boundary layer thickness decreases for large values of porousity parameter.

## 6. Closing remarks

In this paper, the impacts of dependent viscosity parameter, magnetic field, and solid nanoparticle flow and the heat transfer of modified nanofluid flow at the exponential stretching surface have been analyzed numerically. The governing coupled partial differential equations are converted into ordinary coupled differential equations which are solved numerically by bvp4c method. The parametric analysis is executed to investigate the impacts of the governing physical parameters (magnetic field, variable viscosity (for both cases

## Nomenclature

Pr | Prandtl number |

Φ1 | nanoparticle volume fraction of Al2O3 |

Φ2 | nanoparticle volume fraction of Cu |

Φ3 | nanoparticle volume fraction of Ni |

Bi | Biot number |

θ | temperature profile |

R | permeability |

f | velocity profile along x-direction |

ρ | density |

f | fluid |

Tw | wall temperature |

T∞ | ambient temperature |

νf | fluid kinematic viscosity |

νnf | nanofluid kinematic viscosity |

νhnf | hybrid nanofluid kinematic viscosity |

νmnf | modified nanofluid kinematic viscosity |

ρCphnf | heat capacity of hybrid nanofluid |

ρCpmnf | heat capacity of modified nanofluid |

κf | thermal conductivity of fluid |

κnf | thermal conductivity of nanofluid |

κhnf | thermal conductivity of hybrid nanofluid |

κmnf | thermal conductivity of modified nanofluid |

μhnf | viscosity of hybrid nanofluid |

μmnf | viscosity of modified nanofluid |

μnf | viscosity of nanofluid |

ρCpnf | heat capacity of nanofluid |

αhnf | thermal diffusivity of hybrid nanofluid |

αmnf | thermal diffusivity of modified nanofluid |

αnf | thermal diffusivity of nanofluid |

U,V | velocity components |

X,Y | direction components |

θe | variable viscosity parameter |

γ | porosity parameter |