*Article* **Entropy Generation Analysis and Radiated Heat Transfer in MHD (Al2O3-Cu/Water) Hybrid Nanofluid Flow**

**Nabeela Parveen <sup>1</sup> , Muhammad Awais <sup>1</sup> , Saeed Ehsan Awan <sup>2</sup> , Wasim Ullah Khan 3,\*, Yigang He 3,\* and Muhammad Yousaf Malik <sup>4</sup>**


**Abstract:** This research concerns the heat transfer and entropy generation analysis in the MHD axisymmetric flow of Al2O<sup>3</sup> -Cu/H2O hybrid nanofluid. The magnetic induction effect is considered for large magnetic Reynolds number. The influences of thermal radiations, viscous dissipation and convective temperature conditions over flow are studied. The problem is modeled using boundary layer theory, Maxwell's equations and Fourier's conduction law along with defined physical factors. Similarity transformations are utilized for model simplification which is analytically solved with the homotopy analysis method. The h-curves up to 20th order for solutions establishes the stability and convergence of the adopted computational method. Rheological impacts of involved parameters on flow variables and entropy generation number are demonstrated via graphs and tables. The study reveals that entropy in system of hybrid nanofluid affected by magnetic induction declines for *β* while it enhances for *Bi*, *R* and λ. Moreover, heat transfer rate elevates for large *Bi* with convective conditions at surface.

**Keywords:** hybrid nanofluid; entropy generation; induced magnetic field; convective boundary conditions; thermal radiations; stretching disk

#### **1. Introduction**

MHD boundary layer flows of electrically conducting fluids and heat transfer caused by a stretching sheet have gained immense importance due to their ample applications and significant bearings on several engineering and technological processes. Major applications include heat exchangers, metals' spinning, and power generators, spinning of fiber, aerodynamic extrusion of plastic sheets, polymer industries, and condensation process of metallic sheets inside cooling glass. The quality of the resulting product greatly depends on the stretching process as well as its rate of cooling. Boundary-layer flow combined with heat transfer and followed by a stretching sheet was primarily proposed by Sakiadis [1]. Crane [2] also investigated flow caused by stretching sheets of plastic material. Boundary layer equations describing air motion due to a plate were solved analytically. Rates of flow and heat transfer were analyzed in terms of coefficients of friction and thermal conductivity, respectively. Since then, many attempts have been made in regard to their application in different areas. The analysis of viscous dissipation, thermal radiations and convective wall conditions in fluid flow has its importance in view of its broad coverage including chemical, mechanical, and aerospace engineering, paper production, continuous casting, stretching of plastic films, cooling of electronic chips and crystal growing etc. Khan et al. [3] explored analysis of heat and mass transfer in three-dimensional nanofluids flowing on a linear

**Citation:** Parveen, N.; Awais, M.; Awan, S.E.; Khan, W.U.; He, Y.; Malik, M.Y. Entropy Generation Analysis and Radiated Heat Transfer in MHD (Al2O3-Cu/Water) Hybrid Nanofluid Flow. *Micromachines* **2021**, *12*, 887. https://doi.org/10.3390/mi12080887

Academic Editors: Jin-yuan Qian and Lanju Mei

Received: 25 May 2021 Accepted: 19 July 2021 Published: 27 July 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

stretching sheet under convective wall conditions and thermal radiations. It was deduced that heat and mass transfer rates enhance with the stretching parameter. Gireesha et al. [4] investigated the influence of non-linear radiation on MHD boundary layer dynamics of a Jeffrey nanofluid past a non-linear porous stretching plate. Three-dimensional flow of fluid was considered. It was found that magnetic field diminishes the fluid velocity while it enhances temperature. Hayat et al. [5] analyzed the influence of magnetic induction on dynamics of second-grade nanofluid due to stretching sheet with convective wall conditions. The viscoelastic parameter was observed to upsurge the fluid parameter. Rafiq et al. [6] exposed the effects of non-linear thermal radiation towards boundary layer dusty fluid dynamics close to a rotating isothermally heated blunt-nosed object similar to a hemisphere. It was determined that skin friction coefficient shows an asymptotic trend while heat transfer coefficient increases significantly corresponding to large radiation parameter. Khan et al. [7] analytically studied a mixed convection hybrid nanofluid consisting of Al2O<sup>3</sup> and Ag nanoparticles affected by induced magnetic field in order to analyze entropy generation under viscous dissipation and heat generation effects. They observed that viscous dissipation dominantly increases flow and heat transfer rate due to the no-slip surface condition. Few other attempts in this regime are cited here [8–14].

Technological developments and increasing demand of optical and electronic equipment required an improved cooling performance of the products which is acquired by utilizing heat transfer fluids with the improvements in their thermo physical characteristics. In order to obtain modified results, there have been plenty of endeavors by researchers to synthesize these fluids for an efficient heat transfer rate using the composition of several fluids as well as the dispersion of metallic particles of different sizes and shapes etc. Recently, an upsurge of embedding thermal resistive and conductive nano-particles, initially introduced by Choi [15], has been implemented to enhance thermal characteristics of working fluids. Liquids containing suspended nanoparticles, named as nanofluids, were experimentally guaranteed to possess their higher thermal conductivities than that of base fluids and may enable the operation of cooling systems for practical use in many fields such as in micro-electromechanical systems, pharmaceutical procedures, heat transfer, hybridpowered engines and in field of nanotechnology. Dynamics of a magneto convective Casson nanofluid caused by Stefen blowing with bio active mixers was theoretically inspected by Puneeth et al. [16]. Awais et al. [17] theoretically exposed the rheological behavior of copperwater nanofluid peristaltic flowthrough generalized compliant walls in order to inspect influence of Hall and slip with temperature dependent viscosity. It was deduced that first and second order velocity slip parameters significantly increase flow velocity whereas rate of heat transfer is maximum in the vicinity of channel boundaries. An experimental study on the rheological characteristics of nanofluids manufactured by dispersing multi-walled carbon nanotubes in liquid paraffin was carried out by Liu et al. [18]. It was revealed that velocity components enlarge for velocity slip parameters while temperature-dependent viscosity has shown an impact of increasing temperature. In order to characterize the solar energy storage, improvement in thermal capacity of binary nitrate eutectic salt-silica nanofluid was studied by Hu et al. [19]. Relevant literature in the regime of nanofluids under several aspects can be found in [20–22]. Regardless of the noteworthy consequences of researchers' endeavors, authentic applications require transaction in dissimilar characteristics of nanofluids and therefore hybrid nanofluids were synthesized by embedding special nanoparticles in base fluid. Such fluids possess improved physical and chemical properties in a homogeneous phase. Waini et al. [23] inspected MHD flow dynamics and heat transport of a hybrid nanofluid over a porous stretching/shrinking wedge. A drop of the heat transport rate was determined with the rise in radiation parameter. The temporal stability analysis was presented to evaluate the dual solutions' stability, and it was revealed that one of the dual stables is reliable physically. Parveen et al. [24] utilized computational intelligence techniques in order to analyze heat transfer rate and pressure rise behavior in hybrid nanofluid dynamics influenced by magnetic induction effects past an endoscope. It was shown that coefficient of heat transfer accelerates toward Br and χ. Accuracy and stabil-

ity of experimental data were established by employing neural network algorithm and very reliable results were obtained. Radhika et al. [25] explained the effects of magnetic field and melting heat transfer in the dynamics of dusty fluid suspended with hybrid nanoparticles. It was revealed that thermal gradient enhances for high values of magnetic parameter and Prandtl number. Reddy et al. [26] carried out theoretical analysis for heat transfer in dusty fluid dynamics suspended with hybrid nanoparticles by taking the Cattaneo-Christov heat flux model. Khan et al. [27] studied sustainability based performance evaluation of hybrid nanofluid assisted machining. Their theoretical study showed that a very small portion of nanoparticles affect the cost of each industrial product and the study was in complete accordance with the industrial applications of nanofluids. Literature in this area is shown in the references [28–31]. Moreover, the shape factor can approximately portray the difference of shape among non-spherical and spherical nanoparticles. In general, the nanoparticles possess polyhedral shapes, and their surface is made up of various planes. Thermophysical properties of nanoparticles directly depend on shape of nanoparticles. Therefore, flow and heat transport rates are examined in terms of coefficient of skin friction and Nusselt number for the nanoparticles of platelets shape with *s* = 5.7 in this analysis.

Entropy generation analysis is one of the vital factors in fluid mechanics. Performance of thermal devices directly depends upon the available amount of work which degrades by flow irreversibility and causes more disorder. Therefore, the study of dynamics behind entropy production becomes necessary in order to optimize thermal efficiencies of devices. Dogonchi et al. [32] inspected entropy generation behavior in natural convective hybrid nanofluid rheology by considering the effects of applied magnetic field past a porous cavity with wavy walls embedded in three circular cylinders. Sahoo et al. [33] carried out the analysis of entropy optimization, with dissipative heat transfer, in mixed convective MHD Casson nanofluid dynamics under the influences of Hall current and thermal radiation. Results showed that entropy generation amplifies significantly for diffusive variable, Brinkman number, and concentration ratio parameter whereas Bejan number decreases for all these parameters. In this regard, some investigations on entropy generation analysis for different flows and geometries under various physical aspects are reviewed (see articles [34–36]). Moreover, use of an analytical technique for the solution of the mathematical model is aimed by using homotopy analysis method (HAM). HAM, intended by Shi Jun Liao in 1992, depends on the fundamental concept of topology and differential geometry, homotopy. Being an analytical technique, HAMs are able to solve algebraic, ordinary/partial differential and differential-integral, and linear/non-linear equations in terms of series sum. It provides a broader way for selection of its arguments like initial guess, linear operator and convergence control parameter, which can be highly effective to control convergence rate of the solutions. This characteristic of HAM makes it preferable toother analytical techniques.

The objectives of this study are to theoretically analyze entropy generation and rate of heat transfer in steady flow of (Al2O3-Cu/H2O) hybrid nanofluid induced due to radially stretching disk by imposing convective-type thermal conditions. Flow is axisymmetric in which all the flow variables are independent of angular dimension. The influence of induced magnetic field is taken into account. Flow is considered in the presence of viscous dissipation and thermal radiation effects. The complete mechanism of the present study is explored in terms of a workflow diagram in Figure 1.

**Figure 1.** Workflow chart.

## **2. Modeling and Problem**

Steady, boundary layer flow of viscous (Al2O3-Cu/H2O) hybrid nanofluid induced due to stretching disk in radial directions is assumed. The volume concentration of Al2O<sup>3</sup> and Cu nanomaterials is taken to be 0.05%. The stretching velocity of the disk surface is *Uw*(*r*) *= ar* where a represents positive constant. The surface of disk is convectively heated by fluid having temperature T and held in plane *z* = 0 while hybrid nanofluid is flowing in the region *z* > 0 as shown in Figure 2. Moreover, a scheme for manufacturing of nanofluid and hybrid nanofluid for the nanoparticles in the present study is presented in Figure 3.

**Figure 2.** Schematic representation.

**Figure 3.** Preparation scheme of nanofluid and hybrid nanofluids.

The constitutive governing model along with induced magnetic field, thermal radiations and viscous dissipation effects under boundary layer approximation is:

$$
\frac{\partial u}{\partial r} + \frac{u}{r} + \frac{\partial w}{\partial z} = 0,
\tag{1}
$$

$$\frac{\partial H\_1}{\partial r} + \frac{H\_1}{r} + \frac{\partial H\_3}{\partial z} = 0,\tag{2}$$

$$u\frac{\partial u}{\partial r} + w\frac{\partial u}{\partial z} - \frac{\hat{\mu}}{4\pi\rho\_{\text{f}}} \left(H\_{1}\frac{\partial H\_{1}}{\partial r} + H\_{3}\frac{\partial H\_{1}}{\partial z}\right) = -\frac{\hat{\mu}}{4\pi\rho\_{\text{f}}}H\_{\text{e}}\frac{dH\_{\text{e}}}{dr} + \left(\frac{\mu\_{\text{lnf}}}{\rho\_{\text{lnf}}}\right)\frac{\partial^{2}u}{\partial z^{2}},\tag{3}$$

$$
\mu \frac{\partial H\_1}{\partial r} + w \frac{\partial H\_1}{\partial z} - H\_1 \frac{\partial u}{\partial r} - H\_3 \frac{\partial u}{\partial z} = \mu\_\varepsilon \frac{\partial^2 H\_1}{\partial z^2} \,, \tag{4}
$$

$$u\frac{\partial T}{\partial r} + w\frac{\partial T}{\partial z} = a\_{\ln f}\frac{\partial^2 T}{\partial z^2} + \left(\frac{\mu\_{\ln f}}{(\rho c\_p)\_{\ln f}}\right)\left(\frac{\partial u}{\partial z}\right)^2 - \left(\frac{1}{(\rho c\_p)\_{\ln f}}\right)\frac{\partial q\_r}{\partial z},\tag{5}$$

Radial, axial and azimuthal components of the induced magnetic field vector are, *H*1, *H*<sup>2</sup> and *H*3, respectively. Corresponding boundary conditions are:

$$\begin{array}{ll} \mu = \mathcal{U}\_{\mathcal{U}}(r) = ar, & -\chi\_{\text{inf}}\frac{\partial T}{\partial \overline{z}} = h(T\_w - T), \; w = 0, \; \frac{\partial H\_1}{\partial \overline{z}} = 0, \\\ H\_2 = H\_3 = 0, \; at \; z = 0, \\\ u = w = 0, \; T \to T\_{\infty \prime} \; H\_1 = H\_\ell(r) = H\_0 r, \; as \; z \to \infty. \end{array} \tag{6}$$

In Equation (5), expression for *q<sup>r</sup>* by using Roseland approximation and Taylor series expansion of *T* <sup>4</sup> about *T*<sup>∞</sup> can be expressed as:

$$q\_r = -\frac{4\sigma}{3k} \frac{\partial T^4}{\partial z} = -\frac{16\sigma T\_\infty s^3}{3k} \frac{\partial T}{\partial z},\tag{7}$$

Now, using similarity transformation:

$$\begin{aligned} u(r,z) &= arf'(\eta), & w(r,z) &= -2\sqrt{\pi v\_f} f(\eta), & \eta &= \sqrt{\frac{a}{v\_f}} z, \\ H\_1 &= H\_0 r g'(\eta), & H\_3 &= -2\sqrt{\pi v\_f} g(\eta), & \theta(\eta) &= \frac{T - T\_{\infty}}{T\_w - T\_{\infty}}. \end{aligned} \tag{8}$$

For the above transformations, Equations (1) and (2) are identically satisfied while Equations (3)–(5) within boundary conditions (6) and Equation (7) gives:

$$f''' - \Phi\_1 \Phi\_2 \left\{ f'^2 - 2ff'' - \beta(g'^2 - 2\text{gg''} - 1) \right\} = 0,\tag{9}$$

$$
\lambda \mathbf{g}''' + 2f \mathbf{g}'' - 2f'' \mathbf{g} = \mathbf{0},
\tag{10}
$$

$$(\Phi\_4 + \frac{4}{3}R)\theta'' + 2\Phi\_3 Pr f \theta' + \frac{Pr \cdot Ec}{\Phi\_1} f''^2 = 0 \tag{11}$$

The transformed boundary conditions are:

$$f(\eta) = \mathbf{g}(\eta) = 0, \quad f'(\eta) = 1, \quad \mathbf{g}''(\eta) = 0, \; \theta'(\eta) = -\frac{\mathrm{Bi}}{\Phi\_4} (1 - \theta(0)), \; at \; \eta = 0,\tag{12}$$

$$f'(\eta) \to 0, \quad \mathbf{g}'(\eta) \to 1, \quad \theta(\eta) \to 0, \quad \text{as } \eta \to \infty.$$

where prime represents differentiation with *η*; respect to

Moreover, expressions for thermophysical properties are:

$$\begin{split} \Phi\_{1} &= (1-\varrho\_{1})^{2.5} (1-\varrho\_{2})^{2.5}, \quad \Phi\_{2} = (1-\varrho\_{2}) \Big[ (1-\varrho\_{1}) + \varrho\_{1} \Big(\frac{\rho\_{1}}{\rho\_{f}}\Big) \Big] + \varrho\_{2} \Big(\frac{\rho\_{2}}{\rho\_{f}}\Big), \\ \Phi\_{3} &= (1-\varrho\_{2}) \Big[ (1-\varrho\_{1}) + \varrho\_{1} \Big(\frac{\left(\rho c\_{\mathbb{P}}\right)\_{\mathbb{s\_{1}}}}{\left(\rho c\_{\mathbb{P}}\right)\_{f}}\Big) \Big] + \varrho\_{2} \frac{\left(\rho c\_{\mathbb{P}}\right)\_{\mathbb{s\_{2}}}}{\left(\rho c\_{\mathbb{P}}\right)\_{f}}, \\ \Phi\_{4} &= \frac{\kappa\_{2} + (s-1)\kappa\_{bf} - (s-1)\varrho\_{2} \left(\kappa\_{bf} - \kappa\_{2}\right)}{\kappa\_{2} + (s-1)\kappa\_{bf} + \varrho\_{2} \left(\kappa\_{bf} - \kappa\_{2}\right)} \frac{\kappa\_{\kappa\_{1}} + (s-1)\kappa\_{f} - (s-1)\varrho\_{1} \left(\kappa\_{f} - \kappa\_{\kappa\_{1}}\right)}{\kappa\_{\kappa\_{1}} + (s-1)\kappa\_{f} + \varrho\_{1} \left(\kappa\_{f} - \kappa\_{\kappa\_{1}}\right)}, \end{split} \tag{13}$$

Moreover, all the thermophysical characteristics of nanoparticles and base fluid are mentioned in Table 1 while expressions for dimensionless parameters are:

$$\begin{aligned} Pr &= \frac{\left(\mu c\_p\right)\_f}{\kappa\_f}, \quad \lambda = \frac{\mu\_c}{\upsilon\_f}, \quad \beta = \frac{\mu}{4\pi\rho\_f} \left(\frac{H\_0}{a}\right)^2, \\ Ec &= \frac{\left(\mu\_w\right)^2}{\left(c\_p\right)\_f \left(T\_w - T\_\infty\right)}, \; R = \frac{4\sigma T\_\infty^3}{k\kappa\_f}, \; Bi = \frac{\hbar}{\kappa\_f} \sqrt{\frac{\upsilon\_f}{a}} \end{aligned} \tag{14}$$


**Table 1.** Experimental values of various thermophysical properties for base fluid and nanoparticles [27].

Expressions for coefficient of skin friction *C<sup>f</sup>* and local Nusselt number *Nu* are:

$$\mathcal{C}\_f = \frac{\mathfrak{T}\_w}{\rho\_f \left(\mathcal{U}\_w\right)^2} \quad \text{,} \mathcal{N}u = -\frac{rq\_{\mathcal{U}}}{\kappa\_f \left(T\_w - T\_\infty\right)}.\tag{15}$$

In the above expressions, *τ<sup>w</sup>* and *q<sup>w</sup>* symbolize shear stress and heat flux for wall, respectively. The dimensionless form by substituting Equation (8) is:

$$\operatorname{Re}\_{\tau}^{\frac{1}{2}} \mathbb{C}\_{f} = \frac{-1}{\Phi\_{1}} f''(0), \quad \operatorname{Re}\_{\tau}^{-\frac{1}{2}} Nu = -\left(\Phi\_{4} + \frac{4}{3}\mathbb{R}\right) \theta'(0), \tag{16}$$

where, *Re<sup>r</sup>* = *Uwr υf* indicates local Reynolds number.

By adopting the second law of thermodynamics, volumetric entropy generation rate in existence of radiative and dissipative factors can be expressed as:

$$\dot{S}'''\_{\text{Gen}} = \left(\Phi\_4 + \frac{4}{3}R\right)\frac{\kappa\_f}{T\_\infty^2}\left(\frac{\partial T}{\partial y}\right)^2 + \frac{\mu\_f}{\Phi\_1 T\_\infty}\left(\frac{\partial u}{\partial y}\right)^2 + \frac{\mu\_f}{\Phi\_1 T\_\infty}\left(\frac{\mu}{4\pi\rho\_f}\right)\left(\frac{\partial H\_1}{\partial y}\right)^2. \tag{17}$$

The expression for characteristic entropy generation rate is:

$$\dot{\boldsymbol{S}}^{\prime\prime}\_{~0} = \frac{\kappa\_f}{T\_\infty^2 L^2} (T\_w - T\_\infty)^2. \tag{18}$$

Utilizing similarity transformation with Equation (18) in Equation (17), we have:

$$N\_{\rm S} = \frac{\dot{\mathcal{S}}'''\_{\rm Gen}}{\dot{\mathcal{S}}'''\_{0}} = \left(\Phi\_{4} + \frac{4}{3}\mathcal{R}\right) \text{Re}\_{L}\mathcal{\theta}^{\prime 2} + \frac{\varepsilon.Pr.\text{Ec.}\text{Re}\_{L}}{\Phi\_{1}}f''^{2} + \frac{\varepsilon.Pr.\text{Ec.}\text{Re}\_{L}\mathcal{\theta}}{\Phi\_{1}}g''^{2}.\tag{19}$$

where, *Ns*is non-dimensional form of entropy generation number, *Re<sup>L</sup>* <sup>=</sup> *aL*<sup>2</sup> *υf* and *ε* = *T*∞ *Tw*−*T*<sup>∞</sup> demonstrate local Reynolds number and temperature ratio parameter, respectively.

#### **3. Homotopy Analysis Method and Convergence of Solutions**

The dimensionless governing model mentioned in Equations (9)–(11) and subjected boundary conditions of Equation (12) for boundary layer flow of hybrid nanofluid over the stretching disk is analytically solved by employing the homotopy analysis method [7,37–39]. For flow variables, initial guesses are:

$$f\_0(\eta) = 1 - \exp(-\eta) \, \_\prime \mathcal{g}\_0(\eta) = \eta \, \_\prime \theta\_0(\eta) = \frac{Bi}{1 + Bi} \exp(-\eta) . \tag{20}$$

Expressions for linear operators are:

$$T\_1(\eta) = f'''(\eta) - f'(\eta), \ T\_2(\eta) = \mathbf{g'''}(\eta) + \mathbf{g''}(\eta), \ T\_3(\eta) = \theta''(\eta). \tag{21}$$

and

$$T\_1 \left( A\_1 + A\_2 e^{-\eta} - A\_3 e^{-\eta} \right) = T\_2 \left( A\_4 + A\_5 \eta + A\_6 e^{-\eta} \right) = T\_3 (A\_7 + A\_8 \eta) = 0. \tag{22}$$

where, *A1–A<sup>8</sup>* represents constants in general solutions.

#### *3.1. Convergence-Control Parameters*

Suppose an *h* is auxiliary parameter in the frame of HAM which directly affects the convergence of solutions. Let *ξ* ∈ [0, 1] be an embedding parameter, then the problem for zeroth order deformation is constructed as:

$$(1 - \xi)T\_1[f(\eta, \xi) - f\_0(\eta)] = \xi h \mathcal{R}\_1[f(\eta, \xi), g(\eta, \xi), \theta(\eta, \xi)],\tag{23}$$

$$(1 - \mathfrak{F})T\_2[\mathcal{g}(\eta, \mathfrak{f}) - \mathcal{g}\_0(\eta)] = \mathfrak{F}h\mathcal{R}\_2[f(\eta, \mathfrak{f}), \mathcal{g}(\eta, \mathfrak{f}), \theta(\eta, \mathfrak{f})],\tag{24}$$

$$(1 - \xi)T\_3[\theta(\eta, \xi) - \theta\_0(\eta)] = \xi h \mathcal{R}\_3[f(\eta, \xi), g(\eta, \xi), \theta(\eta, \xi)].\tag{25}$$

### Furthermore,

at *η* = *η*<sup>1</sup> = 0, *f*(*η*1, *ξ*) = *g*(*η*1, *ξ*) = 0, *f* 0 (*η*1, *ξ*) = 1, *g* <sup>00</sup>(*η*1, *ξ*) = 0, *θ* 0 (*η*1, *<sup>ξ</sup>*) <sup>=</sup> <sup>−</sup> *Bi A*4 (1 − *θ*(*η*1, *ξ*)). At *η* = *η*<sup>2</sup> → ∞, *f* 0 (*η*2, *ξ*) → 0, *g* 0 (*η*2, *ξ*) → 1, *θ*(*η*2, *ξ*) → 0. (26)

Using the aforementioned quantities, the *m*th order solution series is constructed as:

$$T\_1[f\_m(\eta\_\prime \xi) - \zeta\_m f\_{m-1}(\eta\_\prime \xi)] = h S\_m^1[(\eta\_\prime \xi)]\_\prime \tag{27}$$

$$\frac{1}{2}T\_2[\mathcal{g}\_m(\eta,\xi) - \zeta\_m \mathcal{g}\_{m-1}(\eta,\xi)] = hS\_m^2[(\eta,\xi)],\tag{28}$$

$$T\_3[\theta\_m(\eta\_\prime \xi) - \varrho\_m \theta\_{m-1}(\eta\_\prime \xi)] = hS\_m^3[(\eta\_\prime \xi)].\tag{29}$$

In the above equations,

$$
\varsigma\_m = \begin{array}{c} 0m \le 1 \\ 1m > 1 \end{array} \right]$$

Subjected boundary conditions are:

$$\begin{array}{c} \text{At } \eta = \eta\_1 = 0, \, f\_{\mathfrak{m}}(\eta\_1, \mathfrak{xi}) = \operatorname{g}\_{\mathfrak{m}}(\eta\_1, \mathfrak{xi}) = 0, \, f'\_{\mathfrak{m}}(\eta\_1, \mathfrak{z}) = 0, \, \operatorname{g}\_{\mathfrak{m}}'(\eta\_1, \mathfrak{z}) = 0, \\\ \theta\_{\mathfrak{m}}(\eta\_1, \mathfrak{z}) = 0, \\\ \text{At } \eta = \eta\_2 \to \infty, \, f'\_{\mathfrak{m}}(\eta\_2) \to 0, \, \operatorname{g}'\_{\mathfrak{m}}(\eta\_2) \to 0, \, \theta\_{\mathfrak{m}}(\eta\_2) \to 0. \end{array} \tag{30}$$

Then, we can write:

$$\begin{array}{l} f\_m(\eta, \mathfrak{F}) = F\_m(\eta, \mathfrak{F}) + A\_1 + A\_2 e^{-\eta} - A\_3 e^{-\eta}, \\ g\_m(\eta, \mathfrak{F}) = G\_m(\eta, \mathfrak{F}) + A\_4 + A\_5 \eta + A\_6 e^{-\eta}, \\ \Theta\_m(\eta, \mathfrak{F}) = \Theta\_m(\eta, \mathfrak{F}) + A\_7 + A\_8 \eta. \end{array} \tag{31}$$

#### *3.2. Convergence of Solutions*

For stability and convergence of analytical solutions, h-curves for, *f* 0 (0) *g* 0 (0) and *θ*(0) are prepared upto 20th order solutions by plotting the interval of convergence and are shown in Figures 4–6. Convergence control parameter *h* is able to control and adjust convergence region of HAM solutions. It is observed from convergence analysis that valid convergence region is *R<sup>h</sup>* = [−0.4, 0.2]. For *f* 0 (0), *R<sup>h</sup>* = [−0.35, 0.2] for *g* 0 (0), and *R<sup>h</sup>* = [−0.35, 0.15] for *θ*(0). Moreover, testing various values of *h* from corresponding regions, the error is minimum at *h* = −0.1 which guarantees that calculated solutions are very accurate with a negligible error and become more accurate for higher-order solutions.

**Figure 4.** *h*-curves for solution of *f* 0 (0).

**Figure 5.** *h*-curves for solution of *g* 0 (0).

**Figure 6.** *h*-curves for the solution of *θ*(0).

#### **4. Discussion of Results**

Analytical solutions for the Al2O3-Cu/water hybrid nanofluid flow obtained by HAM are discussed in this section. Figures 7 and 8 plot the variation in magnetic parameter *β* and the reciprocal of magnetic Prandtl number *λ* for the velocity profile. The graph indicates that rise in values of *β* enhances *f* 0 (*η*) due to more dominant induction effects than magnetic diffusion and the flow rate increases. A drop in *f* 0 (*η*) is noticed for *λ* because magnetic diffusivity rises with rise in *λ*. This effect causes enhancement of frictional force and the velocity boundary layer thickness reduces.

**Figure 7.** *β* verses velocity profile.

**Figure 8.** *λ* verses velocity profile.

The behavior of induced magnetic field profile against parameters *β* and *λ* is explored in Figures 9 and 10, respectively. Graphs demonstrate that with augmentation in *β*, magnetic effects become strong as due to fast advection process, therefore increasing flow rate amplifies magnetic induction profile. Consecutively, *g* 0 (*η*) is a decreasing function of *λ* which is mainly due to enhancing magnetic diffusivity with high values of *λ*. It is noteworthy that the influences of parameters on *f* 0 (*η*) and *g* 0 (*η*) are more prominent at the interface and very little variation is shown near stretching disk due to the no slip boundary condition.

**Figure 9.** *β* verses induced magnetic field profile.

**Figure 10.** *λ* verses induced magnetic field profile.

Figures 11–14 are plotted in order to express the variation in temperature against magnetic parameter, reciprocal of magnetic Prandtl number, radiation parameter and Biot number. Results reveal that *θ*(*η*) decreases with increment in *β* which is caused by increasing heat transfer rate near the disk in existence of induced magnetic field due to adding the flow mechanism as explored in Figure 11. The impact of radiation parameter in Figure 12 depicts the dual behavior of the *θ*(*η*) profile which is depressed close to the disk and elevated away from it showing dominant effects at free stream. Since the disk is convectively heated and these effects trim down away from its surface, temperature enhances at the disk surface due to rising intensity of convective heating in Figure 13.

**Figure 11.** *β* verses temperature profile.

**Figure 12.** *λ* verses temperature profile.

**Figure 13.** *Bi* verses temperature profile.

**Figure 14.** *R* verses temperature profile.

An opposite trend is observed in the free stream region due to the fact that additional heat released to the coolant at the surface. The impact of *λ* on temperature of hybrid nanofluid is exposed in Figure 14 which signifies that temperature rises against *λ*. The reason is high diffusivity and a small flow rate due to low magnetic induction for large magnitude of *λ*. Moreover, variation in *θ*(*η*) directly associated with values of parameters *Bi*, *R* and *Ec* which are taken to be small with *Pr* = 6.8 in this study.

The most persuasive part of this section is entropy generation analysis. The variation in entropy generation number Ns against emerging parameters is pointed out graphically in Figures 15–18. Impact of magnetic parameter in Figure 15 demonstrates that rise in *β* decreases entropy production as strong effects of magnetic field close to the disk reduce the frictional effects and rise in heat transfer rate. It is of true significance physically since the rate of fluid flow enhances instantly and heat transfer rate increases. The number slightly increases in the free stream region due to decreasing fluid flow. Figure 16 expresses an enhancement in Ns as values of reciprocal of magnetic Prandtl number rises. This is because of the fact that as *λ* enhances, reduction in viscosity occurs while magnetic diffusivity accelerates, which produce disorder in the system. In Figure 17, entropy generation number is plotted against *Bi* which illustrates that for somewhat large values of parameter *Bi*, Ns have high magnitude due to strong influences of thermal convection and maximum radial gradient. Additionally, thermal radiation effects on Ns in Figure 18 displays it as an increasing function of R due to increasing emitting radiations which boosts frictional irreversibility that encourage entropy generation. Furthermore, variations in entropy generation number are drawn for small values of radiation parameter *R*, Eckert number Ec and Biot number *Bi* which are directly related to entropy production. Also, extensive behavior for Ns against parameters is observed at the interface due to large velocity and temperature gradients caused by no slip at the surface and convective wall conditions but it is rarely affected by these parameters away from it.

**Figure 15.** *β* verses entropy generation number.

**Figure 16.** *λ* verses entropy generation number.

**Figure 17.** *Bi* verses entropy generation number.

**Figure 18.** *R* verses entropy generation number.

Moreover, empirical formulas for the thermophysical properties of nanofluid and hybrid nanofluid are displayed in Table 2.

**Table 2.** Formulas for thermophysical properties of nanofluid and hybrid nanofluid [40].


The impacts of involving parameters on skin friction coefficient and heat transfer rate through Tables 3 and 4 as well as bar graphs in Figures 19–25 are explained in this section. Observations conclude that the surface velocity gradient is depressed as values of magnetic parameter *β* rises but enhances with elevation in λ. Moreover, an expansion in Nusselt number is noticed for rising values of the reciprocal of magnetic Prandtl number, Biot number, radiation parameter, and Eckert number whereas it is a decreasing function of magnetic parameter.

**Table 3.** Behavior of skin friction coefficient for Al2O<sup>3</sup> -Cu/H2Oagainst several parameters.



**Table 4.** Behavior of Nusselt number for Al2O<sup>3</sup> -Cu/H2O against several parameters.

**Figure 19.** *β* verses *C<sup>f</sup>* .

**Figure 20.** λ verses Cf.

**Figure 21.** *β* verses *Nu*.

**Figure 22.** λ verses *Nu*.

**Figure 23.** *Ec* verses *Nu*.

**Figure 24.** *R* verses *Nu*.

#### **5. Concluding Remarks**

A theoretical discussion of MHD viscous flow of Al2O3-Cu/H2O hybrid nanofluid due to stretching of the disk is carried out. Some noteworthy influences of important emerging parameters on flow profiles and entropy generation are as follows:


**Author Contributions:** Conceptualization, N.P.; Data curation, N.P.; Formal analysis, M.A. and M.Y.M.; Investigation, S.E.A.; Methodology, W.U.K.; Project administration, Y.H.; Resources, W.U.K. and Y.H.; Software, S.E.A. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by theNational Natural Science Foundation of Chinaunder Grant No.51977153,51977161,51577046, State Key Program ofNational Natural Science Foundation of Chinaunder Grant No.51637004,National Key Research and Development Plan (China)"important scientific instruments and equipment development" Grant No.2016YFF010220,Equipment research project in advance (China)Grant No.41402040301, Wuhan Science and Technology Plan Project Grant No 20201G01.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study is available in article.

**Acknowledgments:** All authors sincerely thank the editor and the referees for the useful suggestions regarding improvement of our work. M. Y. Malik present his appreciations to the Deanship of Scientific Research at King Khalid University, Abha, 61413, Saudi Arabia for support through R.G.P-1-223-42.

**Conflicts of Interest:** Authors have no conflict of interest.

## **Nomenclature**


#### **References**


*Article*
