Entropy Generation and Heat Transfer Analysis in MHD Unsteady Rotating Flow for Aqueous Suspensions of Carbon Nanotubes with Nonlinear Thermal Radiation and Viscous Dissipation Effect

Abstract

The impact of nonlinear thermal radiations rotating with the augmentation of heat transfer flow of time-dependent single-walled carbon nanotubes is investigated. Nanofluid flow is induced by a shrinking sheet within the rotating system. The impact of viscous dissipation is taken into account. Nanofluid flow is assumed to be electrically conducting. Similarity transformations are applied to transform PDEs (partial differential equations) into ODEs (ordinary differential equations). Transformed equations are solved by the homotopy analysis method (HAM). The radiative source term is involved in the energy equation. For entropy generation, the second law of thermodynamics is applied. The Bejan number represents the current investigation of non-dimensional entropy generation due to heat transfer and fluid friction. The results obtained indicate that the thickness of the boundary layer decreases for greater values of the rotation parameter. Moreover, the unsteadiness parameter decreases the temperature profile and increases the velocity field. Skin friction and the Nusselt number are also physically and numerically analyzed.

1. Introduction

Applications of nanofluids in technology and science are increasing day by day, and they play an important role in various machinery and engineering applications such as detergents, microchip technology, transferences, micromechanical systems and biomedical applications. In light of these applications, researchers have applied modern techniques and have modified the base fluids by adding ultra-fine solid particles. A fast-growing field of research is micro channel cooling, floor heating and heat renewal systems in various industries, which have been flourishing in the current era. In comparison to base liquids, nanofluids have higher single-phase thermal conductivity and heat transfer coefficients.

Khan [1] scrutinized nanofluid flow for Buongiorno’s model with the transfer of heat and mass. Mahdy et al. [2] depicted Buongiorno’s model for an unsteady nanofluid flowing in a contracting cylinder in the presence of heat transfer. The flow of nanofluid through a vertical annular pipe was deliberated by Malvandi et al. [3]. Researchers have examined nanofluid flow over a stretching sheet [4,5,6]. Non-Newtonian MHD nanofluid flow through a pipe was depicted by Ellahi [7]. Jawad et al. [8] studied the nanofluid thin film flow of Sisko fluid. Nanofluid thin film flow through a stretching sheet with heat transfer was studied by Fakour et al. [9]. The flow of nanofluids with analytical techniques was investigated by Abolbashari et al. [10]. Choi [11] was the innovator who developed the word nanofluid. Nanofluids are a mixture of metallic, nanoscale, suspended particles with a base fluid. For mass, momentum and heat transport, non-homogeneous equilibrium models in nanofluids embody four equations with two components, as expressed by Buongiorno [12]. Some interesting results pertaining to the use of nanofluids can be found in [13,14,15,16]. Recently, insufficient expressive efforts for nonlinear thermal radiations have been observed in these investigations [17]. Kumar et al. [18] scrutinized the problem of entropy generation in rotating nanofluid flow. Nadeem et al. [19] examined the issue of flow of rotating fluid, including nanoparticles of titanium and copper oxide. Mabood et al. [20] explored the flow of rotating nanofluids with the effects of radiation, magnetism, the dissipation of viscosity and a heat source. Shah et al. [21,22] investigated the flow of nanofluids within a rotating system under the influences of hall current and thermal radiation. Further theoretical investigations of nanofluids using modern applications were performed by Sheikholeslami [23,24,25]. Gireesha et al. [26] explored a single-walled nanotube in an unsteady rotating flow weight transfer. Ishaq et al. [27] examined an unsteady nanofluid thin film flow with non-dimensional entropy generation through a stretched surface.

In the study of fluid flow involving nanoparticles, nanofluids have received more attention with the arrival of nanoscience. Nanofluids are 10⁹ nm sized materials such as nanotubes, nanofibers, droplets, nanoparticles, etc. The solid phase and liquid phase create a two-phase system. Moreover, stable fluids have good writing and spreading properties on hard surfaces. Using nanofluids, the thermal conductivity of fluids can be enhanced [28,29]. Sheikholeslami et al. [30,31,32] stressed the significance of nanofluids in nanotechnology. Yadav et al. [33,34] explored MHD nanofluid flow with dissimilar phenomena such as heat transfer enhancement, stability, instability, and linear and nonlinear properties in a nanofluid model.

Recently, it has been suggested that the analysis of thin film flow has pointedly contributed in different areas, such as industry, engineering and technology, etc. The study of nanofluids was improved due to their vast applications. The latest investigations of thin film flow using different models in different geometries can be seen in [35,36,37,38,39,40,41].

Entropy is the amount of unobtainable energy in a thermodynamically closed system. The total entropy remains constant in a steady-state system. Chemical reactions, thermal resistance, joule heating, diffusion and friction concerning fluid viscosity and hard surfaces within a system are irreversible processes. Entropy production minimization is necessary for the maximal usefulness of equipment [42,43,44,45,46,47,48,49]. The entropy generated by squeezing the nanofluid flow in three-dimensions between two parallel plates is considered here. Shah et al. [50,51] proposed the heat transfer model for nanofluids. Ahmad et al. [52] deliberated squeezing the time-dependent flow of viscous nanoparticles under five dissimilar shapes. Rehman et al. [53] deliberated the flow of rotating nanofluid with entropy generation included in the thermal slip. The impacts of dissimilar constraints on liquid flow and heat transfer in a rotating fluid over a stretched disk were studied recently [54]. Alnaqi et al. [55] explored nanofluid flow with the effects of a magnetic field on the convective heat transfer rate and entropy generation through an inclined square cavity equipped with a conductor fin. Moradikazerouni et al. [56] studied CPU heat sinks in computers using a structural stability method. Hajizadeh et al. [57] investigated the thermal conductivity enhancement of nano-antifreeze containing single-walled carbon nanotubes. Vo et al. [58] deliberated γ-AlOOH nano-fluid convection performance by using various shapes of nano-additives. Alsarraf et al. [59] explored nanofluid flow with different nanoparticle shapes in a mini-channel heat exchanger using a two-phase mixture model. Moradikazerouni et al. [60] investigated the effects of five different channel forms of a micro-channel heat sink in forced convection, with application to cooling a supercomputer circuit board.

The aim of the current research is to obtain an analytical solution using the homotopy analysis method (HAM) for an unsteady, MHD, and the incompressible rotating flow of carbon nano tubes nanofluid over a shrinking surface with nonlinear thermal radiation and viscous dissipation effect. Here, we consider three types of nanofluids: CuO-water, Ag-water and Au-water, where water is used as the base fluid. The impact of the first order chemical reaction is also deliberated. The problem is formulated, solved and the corresponding results are examined in detail. Finally, the impact of the physical parameters on temperature and concentration profiles are presented and analyzed.

2. Mathematical Formulation of the Problem

Single-walled carbon nanotube nanofluid unsteady laminar incompressible three-dimensional rotating flow is considered over a shrinking surface. The Cartesian coordinates are chosen in the x, y, and z dimensions. The nanofluid rotates with an angular velocity about the z-axis, which is denoted by $\overline{\mathsf{\Omega }}\left(t\right)$. The surface velocity is represented by $u_w\left(x,t\right)$ and given as $u_w\left(x,t\right)=\frac{bx}{\left(1-\delta t\right)}$ in the x direction, $v_w \left(x, t\right)$ in y direction and $w_w\left(x,t\right)$ in the z direction, and $w_w\left(x,t\right)$ is the wall mass flux velocity. The nanofluid flow is assumed to be thermally conductive. Radiative and viscous dissipation effects are taken into account.

Using all these assumptions, the governing equations are written as [9,25,44,45]:

$$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}+\frac{\partial w}{\partial z}=0$$

(1)

$$\frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}+w\frac{\partial u}{\partial z}+\frac{2\overline{\mathsf{\Omega }}v}{1-\delta t}=-\frac{1}{\rho }\frac{\partial p}{\partial x}+\frac{{\mu }_{nf}}{{\rho }_{nf}}\frac{{\partial }^{2}u}{\partial z^2}-\frac{\sigma {\beta }_0^2}{1-\delta t}u-\frac{{\mu }_{nf}}{k^{*}}\frac{u}{1-\delta t}$$

(2)

$$\frac{\partial v}{\partial t}+u\frac{\partial v}{\partial x}+v\frac{\partial v}{\partial y}+w\frac{\partial v}{\partial z}+\frac{2\overline{\mathsf{\Omega }}u}{1-\delta t}=-\frac{1}{\rho }\frac{\partial p}{\partial y}+\frac{{\mu }_{nf}}{{\rho }_{nf}}\frac{{\partial }^{2}v}{\partial z^2}-\frac{\sigma {\beta }_0^2}{1-\delta t}v-\frac{{\mu }_{nf}}{k^{*}}\frac{v}{1-\delta t}$$

(3)

$$\frac{\partial w}{\partial t}+u\frac{\partial w}{\partial x}+v\frac{\partial w}{\partial y}+w\frac{\partial w}{\partial z}=-\frac{1}{\rho }\frac{\partial p}{\partial z}+\frac{{\mu }_{nf}}{{\rho }_{nf}}\frac{{\partial }^{2}w}{\partial z^2}$$

(4)

$$\frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}+w\frac{\partial T}{\partial z}={\alpha }_{nf}\frac{{\partial }^{2}T}{\partial z^2}+\frac{1}{{\left(\rho c_p\right)}_{nf}}\frac{\partial q_r}{\partial z}+\frac{{\mu }_{nf}}{{\left(\rho c_p\right)}_{nf}}\left[{\left(\frac{\partial u}{\partial z}\right)}^2+{\left(\frac{\partial v}{\partial z}\right)}^2+{\left(\frac{\partial w}{\partial z}\right)}^2\right]$$

(5)

The boundary conditions are:

$$\begin{array}{l}u=u_w\left(x,t\right),v=0,w=0,T=T_w\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}at\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}z=0\\ u\to 0,v\to 0,w\to 0,T\to T_{\infty }\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}at\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}z\to \infty \end{array}$$

(6)

where $x,y$ and z are the directions of the velocity components; $\overline{\mathsf{\Omega }}$ denotes a constant angular velocity; ${\mu }_{nf}$ denotes the nanofluid dynamic viscosity; ${\rho }_{nf}$ denotes the nanofluid density ${\alpha }_{nf}$; $T$ represents the temperature of the nanofluid; $T_w$ and $T_{\infty }$ are the wall and the outside surface temperatures, respectively. The radiative heat flux in Equation (5) can be shown as [9,25]:

$$q_r=-\frac{4{\sigma }^{*}}{3{\left(\rho c_p\right)}_{nf}{k}^{*}}\frac{\partial T^4}{\partial z}=-\frac{16{\sigma }^{*}}{3k^{*}}{T}^3\frac{\partial T}{\partial z}$$

(7)

where the Stefan–Boltzmann constant and the mean absorption coefficient are denoted by ${\sigma }^{*}$ and $k^{*}$, respectively. By substituting Equation (7) into Equation (4), it can be written as [6,7,8,9,25]:

$$\begin{array}{l}\frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}+w\frac{\partial T}{\partial z}\\ ={\alpha }_{nf}\frac{{\partial }^{2}T}{\partial z^2}+\frac{16{\sigma }^{*}}{{\left(\rho c_p\right)}_{nf}{k}^{*}}\left[T^2\frac{{\partial }^{2}T}{\partial z^2}+3T^2{\left(\frac{\partial T}{\partial z}\right)}^2\right]+\frac{{\mu }_{nf}}{{\left(\rho c_p\right)}_{nf}}\left[{\left(\frac{\partial u}{\partial z}\right)}^2+{\left(\frac{\partial v}{\partial z}\right)}^2+{\left(\frac{\partial w}{\partial z}\right)}^2\right]\end{array}$$

(8)

Other parameters with nanoparticle volume fraction are mathematically presented as [6,7,22,23]:

$$\begin{array}{l}{\rho }_{nf}=\left(1-\varphi +\varphi \left(\frac{{({\rho }_s)}_{CNT}}{{\rho }_f}\right)\right),{\mu }_{nf}=\frac{{\mu }_f}{{\left(1-\varphi \right)}^{2.5}},{\alpha }_{nf}\\ =\left(1-\varphi \right){\left({\rho }_f\right)}_f+\varphi {({\rho }_s)}_{CNT},\\ {\left(\rho C_p\right)}_{nf}=\left(1-\varphi \right){\left(\rho C_p\right)}_f+\varphi {\left(\rho C_p\right)}_{CNT},\frac{k_{nf}}{k_f}\\ =\frac{1-\varphi +2\varphi \left(\frac{k_{CNT}}{k_{CNT}-k_f}\right)\ln \left(\frac{k_{CNT}+k_f}{2k_f}\right)}{1-\varphi +2\varphi \left(\frac{k_f}{k_{CNT}-k_f}\right)\ln \left(\frac{k_{CNT}+k_f}{2k_f}\right)}.\end{array}$$

(9)

In Equation (9), for the base fluid the volumetric heat capacity is denoted as ${\left(\rho C_p\right)}_f$ and for CNTs as ${\left(\rho C_p\right)}_{CNT}$. The thermal conductivity of the nanofluid, base fluid, and CNTs are denoted as $k_{nf}$, $k_f$ and $k_{CNT}$, respectively. The nanoparticle volume fraction is denoted by $\varphi$; the density viscosity of CNTs and the base fluid are represented by ${\rho }_{CNT}$ and ${\rho }_f$, respectively.

Similarity transformations [26] are introduced as:

$$\begin{array}{l}u=\frac{bx}{\left(1-\delta t\right)}{ƒ}^{\prime }\left(\mathsf{\eta }\right), v=\frac{bx}{\left(1-\delta t\right)}g\left(\mathsf{\eta }\right),\mathrm{w}=-\sqrt{\frac{b\nu }{\left(1-\alpha t\right)}}f\left(\mathsf{\eta }\right), \mathsf{\eta }=\sqrt{\frac{b}{\nu \left(1-\alpha t\right)}}z,\\ T=T_{\infty }\left(1+\left(1-\delta t\right)\theta \left(\mathsf{\eta }\right)\right).\end{array}$$

(10)

Using Equation (10) and Equations (2)–(6), we obtain:

$$\begin{array}{l}\frac{1}{{\left(1-\varphi \right)}^{2.5}\left(\left(1-\varphi \right)+\varphi \frac{{({\rho }_s)}_{CNT}}{{\rho }_f}\right)}{f}^{iv}-\frac{{\beta }^2}{\lambda }{f}^{″}-\frac{{\mu }_{nf}}{{\left(1-\varphi \right)}^{2.5}{k}^{*}b}{f}^{″}\\ -\left[\lambda \left(f^{\prime }+\frac{\eta }{2}{f}^{″}\right)+f^{\prime }{f}^{″}-f{f}^{″}-2\frac{\mathsf{\Omega }}{b}{g}^{\prime }\right]=0\end{array}$$

(11)

$$\begin{array}{l}\frac{1}{{\left(1-\varphi \right)}^{2.5}\left(\left(1-\varphi \right)+\varphi \frac{{({\rho }_s)}_{CNT}}{{\rho }_f}\right)}{g}^{‴}-\frac{{\beta }^2}{\lambda }{g}^{\prime }+\frac{{\mu }_{nf}}{{\left(1-\varphi \right)}^{2.5}{k}^{*}b}{g}^{\prime }\\ -\left[\lambda \left(g^{\prime }+\frac{\eta }{2}{g}^{″}\right)+g{f}^{″}-f{g}^{″}-2\frac{\mathsf{\Omega }}{b}{f}^{″}\right]=0\end{array}$$

(12)

$$\begin{array}{l}{\theta }^{″}+R\left[{\left(1+\left({\theta }_w-1\right)\theta \right)}^3{\theta }^{″}+3{\left(1+\left({\theta }_w-1\right)\theta \right)}^2\left({\theta }_w-1\right){{\theta }^{\prime }}^2\right]+\frac{Ec}{\mathrm{Pr}}\left({f^{″}}^2+{g^{\prime }}^2\right)\\ -\frac{1}{\mathrm{Pr}}\left[\lambda \frac{\eta }{2}{\theta }^{\prime }-f{\theta }^{\prime }\right]=0\end{array}$$

(13)

The transformed boundary conditions are written as:

$$\begin{array}{l}f\left(0\right)=0,f^{\prime }\left(0\right)=1,g\left(0\right)=0,\theta \left(0\right)=1\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}at\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\eta =0\\ f^{\prime }\left(\eta \right)\to 0,g\left(\eta \right)\to 0,f\left(\eta \right)\to 0,\theta \left(\eta \right)\to 0\hspace{0.17em}\hspace{0.17em}at\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\eta \to \infty .\end{array}$$

(14)

where $\mathsf{\Omega }=\frac{\omega }{b}$ is the rotation parameter, $\lambda =\frac{\delta }{b}$ is the unsteadiness parameter and $R=\frac{16{\sigma }^{*}{T}_{\infty }^3}{3k_{nf}{k}^{*}}$ is the radiation parameter. The Prandtl number is denoted by $\mathrm{Pr}=\frac{{\alpha }_{nf}}{{\nu }_{nf}}$, the Eckert number is denoted by $\mathrm{Ec}=\frac{u_w^2}{c_p\left(T-T_{\infty }\right)}$ and the temperature ratio parameter is denoted by ${\theta }_w=\frac{T_w}{T_{\alpha }}$.

2.1. Physical Quantities of Interest

Skin friction in the x and y directions is denoted as $C_{fx}$ and $C_{fy}$, respectively, and the Nusselt number is $N{u}_x$. These are defined as:

$$C_{fx}=\frac{{\tau }_{wx}}{{\rho }_f{\mathrm{u}}_w^2(x,t)},C_{fy}=\frac{{\tau }_{wy}}{{\rho }_f{\mathrm{u}}_w^2(x,t)},N{u}_x=\frac{x{q}_w}{\left(T_w-T_{\infty }\right)},$$

(15)

where ${\tau }_{wx}$ and ${\tau }_{wy}$ are the surface shear stress in the x and y directions, respectively; $q_w$ is the surface heat flux. These can be defined as:

$${\tau }_{wx}={\mu }_{nf}{\left(\frac{\partial u}{\partial z}\right)}_{z=0},\hspace{0.17em}{\tau }_{wy}={\mu }_{nf}{\left(\frac{\partial v}{\partial z}\right)}_{z=0},q_w=-k_{nf}\left(\frac{\partial T}{\partial z}\right)+{\left(q_r\right)}_{z=0}.$$

(16)

Using Equations (16) and (17), we obtain:

$$\sqrt{{\mathrm{Re}}_x}{C}_{fx}=\frac{1}{{\left(1-\varphi \right)}^{2.5}}{f}^{″}\left(0\right),\sqrt{{\mathrm{Re}}_x}{C}_{fy}=\frac{1}{{\left(1-\varphi \right)}^{2.5}}{g}^{\prime }\left(0\right),\hspace{0.17em}\frac{N{u}_x}{\sqrt{{\mathrm{Re}}_x}}=\frac{k_{nf}}{k_f}\left(-\left[1+R{\theta }_w^3\right]{\theta }^{\prime }\left(0\right)\right).$$

(17)

The local Reynolds number is denoted by ${\mathrm{Re}}_x=u_{w}x/v$.

2.2. Entropy Generation and Bejan Number

The dimensional local entropy rate per unit volume for a nanofluid is given by [42,43,44,45,46,47,48,49,50,51,52,53]:

$$S_{g,t}=\frac{K_{nf}}{T_{\alpha }^2}\left[{(\nabla T)}^2+\frac{16{\sigma }_e{T}_{\alpha }^3}{3{\beta }_R{K}_f}{(\nabla T)}^2\right]+\frac{{\mu }_{nf}}{T_{\alpha }}\mathsf{\Phi },$$

(18)

where $\nabla T=\frac{\partial T}{\partial x}+\frac{\partial T}{\partial y}+\frac{\partial T}{\partial z}$ and $\mathsf{\Phi }$ represent viscous dissipation.

In this instance,

$$S_{g,t}=\frac{K_{nf}}{T_{\alpha }^2}\left[1+\frac{16{\sigma }_e{T}_{\alpha }^3}{3{\beta }_R{K}_f}\right]{\left(\frac{\partial T}{\partial z}\right)}^2+\frac{{\mu }_{nf}}{T_{\alpha }}\left[{\left(\frac{\partial u}{\partial z}\right)}^2+{\left(\frac{\partial v}{\partial z}\right)}^2\right]$$

(19)

$$S_{g,t}=S_h+S_R+S_f.$$

(20)

We have three sub-generators of entropy, as deduced from Equation (20). The heat transfer dimensional entropy generator is represented by $S_h$. Due to thermal radiation, the dimensional entropy generator is represented by $S_R$ and the inter-friction of the fluid layers is represented by $S_f$. $S_{g,\mathrm{c}}$ is defined as:

$$S_{g,\mathrm{c}}=\frac{K_{nf}{(\Delta T)}^2}{L^2{T}_{\alpha }^2}.$$

(21)

Now, the non-dimensional $Ns$ (Nusselt number) is defined as:

$$Ns=\frac{S_{g,t}}{S_{g,c}}.$$

(22)

To evaluate the non-dimensional $Ns$ we use Equations (19) and (21) combined with Equations (10) and (22) to obtain:

$$Ns=\frac{S_{g,t}}{S_{g,c}}={\mathrm{Re}}_L\left[1+\frac{4}{3}R\right]{\left({\theta }^{\prime }\left(\eta \right)\right)}^2+{\mathrm{Re}}_L\frac{{\beta }_r}{{\left(1-\varphi \right)}^{2.5}\mathsf{\Omega }}\left[{\left(f^{″}\left(\eta \right)\right)}^2+{\left(g^{\prime }\left(\eta \right)\right)}^2\right].$$

(23)

Here $A=\frac{k_{nf}}{k_f},$ $Br$ and ${\mathrm{Re}}_L$ are the Brinkman and Reynolds numbers, respectively, and $\mathsf{\Omega }$ is the non-dimensional temperature, which can be shown as:

$$T_c=\frac{T_w-T_{\alpha }}{T_{\alpha }}$$

(24)

Equation (23) can be rewritten as:

$$Ns=N_h,+N_R+N_f$$

(25)

where the fluid friction, thermal radiation and heat transfer non-dimensional entropy generators are denoted by $N_h$, $N_R$ and $N_f$, respectively.

The mathematical description of the Bejan number ($Be$) is:

$$Be=\frac{N_h}{N_h+N_R+N_f}$$

(26)

$$Be=\frac{K_{nf}{\mathsf{\Omega }}^2}{L^2}{\mathrm{Re}}_l{{\theta }^{\prime }}^2/{\mathrm{Re}}_L\left[1+\frac{4}{3}R\right]{{\theta }^{\prime }}^2+\frac{{\beta }_r}{{\left(1-\varphi \right)}^{2.5}\mathsf{\Omega }}{\mathrm{Re}}_L\left[\left({f^{″}}^2\right)+\left({g^{\prime }}^2\right)\right].$$

(27)

From Equation (26), it is clear that the Bejan number is limited to the unit interval [0, 1].

3. Solution Procedure

The modeled Equations (11)–(13) with boundary conditions from Equation (14), together with the conditions from Equations (23) and (27), are solved with HAM. The homotopy analysis method is applied due to its outstanding results in boundary layer equations. Several researchers [46,47,48,49,50] have used HAM due to it fast convergence. The preliminary guesses are selected as follows:

$L_{\widehat{f}},\hspace{0.17em}{L}_{\widehat{g}}{ \mathrm{and} \mathrm{L}}_{\widehat{\theta }}$ are linear operators which are represented as.

$${\mathsf{\Upsilon }}_f(f)=\frac{d^{3}f}{d{f}^3},{\mathsf{\Upsilon }}_g(g)=\frac{d^{2}g}{d{g}^2},{\mathsf{\Upsilon }}_{\theta }(\theta )=\frac{d^{2}g}{d{g}^2}{\widehat{\theta }}^{″}$$

(28)

They have the following applicability:

$$\begin{array}{l}{L}_f({\mathsf{\Upsilon }}_1+{\mathsf{\Upsilon }}_2\eta +{\mathsf{\Upsilon }}_3{\eta }^2+{\mathsf{\Upsilon }}_4{\eta }^3)=0,L_g({\mathsf{\Upsilon }}_5+{\mathsf{\Upsilon }}_6\eta +{\mathsf{\Upsilon }}_7{\eta }^3)=0\\ L_{\theta }({\mathsf{\Upsilon }}_8+{\mathsf{\Upsilon }}_9\eta )=0\end{array}$$

(29)

where ${\mathsf{\Upsilon }}_k(k=1,2,3,\dots ,9)$ is constant.

4. Results and Discussion

In this section, the physical outcome of dissimilar parameters of the modeled problems and their effects on $f^{\prime }(\eta )$, $g(\eta )$ and $\theta \left(\eta \right)$ are discussed in detail. The effect of $\mathsf{\Omega }$, $\beta$, $\varphi$ and $\lambda$ on the velocity profile is shown in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8. The impact of $\mathsf{\Omega }$ on $f^{\prime }(\eta )$ and $g(\eta )$ is presented in Figure 1 and Figure 2. It can be seen that for larger values of $\mathsf{\Omega }$ the velocity profile ($f^{\prime }(\eta )$) is increased while $g(\eta )$ is decreased. Actually, increasing the rotation parameter enhances the kinetic energy, which consequently increases the velocity profile, whereas the transverse velocity ($g(\eta )$) is reduced with higher values of the rotation parameter. Figure 3 and Figure 4 represent the influence of $\varphi$ on $f^{\prime }\left(\eta \right)$ and $g\left(\eta \right)$. The higher values of $\varphi$ reduce the velocity profiles. This is because the increase in $\varphi$ further increases the density of the nanofluid, and as a result slows down the fluid velocity profile. Figure 5 and Figure 6 describe the effect of $\lambda$ on $f^{\prime }(\eta )$ and $g\left(\eta \right)$. It was perceived that increases in $\lambda$ reduce the velocity profile. It is also indicated from the figure that the velocity intensifies with increasing $\lambda$, whereas we observed the opposite influence of $\lambda$ on the fluid velocity inside the nanofluid and the thickness of the layer. Figure 7 and Figure 8 show the influence of $\beta$ on $f^{\prime }(\eta )$ and $g(\eta )$. With an increase in $\beta$ the velocity profile of the fluid film is decreased. It was also detected that an increase in $\beta$ results in a decrease in the fluid velocity of the nanofluid and the layer thickness. The purpose behind this influence of $\beta$ by the stimulation of a lingering body force, stated as the Lorentz force, is due to the existence of $\beta$ in an electrically conducting nanofluid layer. The action of this force is perpendicular to both fields. Since $\beta$ represents the ratio of the viscous force to the hydromagnetic body force, a larger value of $\beta$ specifies a higher hydromagnetic body force, due to which the fluid flow is reduced. The Lorentz force theory states that $\beta$ has a converse consequence on $f^{\prime }(\eta )$ and $g(\eta )$. Therefore, the greater values of $\beta$ reduce $f^{\prime }(\eta )$ and $g(\eta )$.

The influence of the physical parameters $\lambda$, $Ec$, $Rd$ and $\mathrm{Pr}$ on $\theta \left(\eta \right)$ is shown in Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13. Figure 9 presents the impact of $\lambda$ on the $\theta \left(\eta \right)$ profile. Figure 9 shows that a decrease in $\lambda$ reduces the boundary layer thickness. It can be seen that when unsteadiness in the stretching increases, the thin film fluid temperature and the free surface temperature are consequently reduced. Under consitions of stirring, it was revealed that greater values of $\lambda$ cause the fluid temperature to fall radically, while the thickness of the thermal boundary layer is increased. The influence of $Rd$ on $\theta \left(\eta \right)$ is shown in Figure 10. By increasing $Rd$, the temperature of the nanofluid boundary layer area increased. In fact, when $Rd$ is raised, then it is obvious that it increases $\theta \left(\eta \right)$ in the boundary layer area in the fluid layer. It is shown in Figure 11 that $\theta \left(\eta \right)$ increases with a rise in $Rd$. Thermal radiation has a dominating role in the comprehensive surface heat diffusion when the coefficient of convection heat transmission is small. Increasing $Rd$ then increases the temperature in the boundary layer area in the fluid layer. This increase leads to a drop in the rate of cooling for nanofluid flow. Therefore, the fluid $\theta \left(\eta \right)$ is increased. The graphical representation shows that $\theta \left(\eta \right)$ is increased when we increase the ratio strength and thermal radiation temperature. Thermal radiation has an important role in heat conduction when the coefficient of convection heat transmission is small. The impact of $\mathrm{Pr}$ on $\theta \left(\eta \right)$ given in Figure 11. It was observed that $\theta \left(\eta \right)$ decreases with larger values of $\mathrm{Pr}$, while it rises for smaller values. The variation of $\theta \left(\eta \right)$ with respect to the variation of $\mathrm{Pr}$ is illustrated and shows that $\mathrm{Pr}$ specifies the ratio of momentum diffusivity to thermal diffusivity. It can be concluded that $\theta \left(\eta \right)$ decreases with increasing $\mathrm{Pr}$. The nanofluids have a greater thermal diffusivity with a small $\mathrm{Pr}$, but this influence does not hold for larger values of $\mathrm{Pr}$; hence, $\theta \left(\eta \right)$ of a fluid displays a reducing behavior. Actually, the fluids having a smaller $\mathrm{Pr}$ have a greater thermal diffusivity, and this impact is the reverse for greater values of $\mathrm{Pr}$. Based on this, a very large value of $\mathrm{Pr}$ causes the thermal boundary layer to drop. Figure 13 shows that with increasing $Ec$, $\theta \left(\eta \right)$ is enlarged, which is supported by the physics. By increasing $Ec$, heat stored in the liquid is dissipated, causing the temperature to be enhanced. $\theta \left(\eta \right)$ is increased with greater values of $Ec$ and the thermal boundary layer thickness of the nanofluid becomes larger.

Now, we analyze the impact of the parameters that perform a role in entropy generation and the Bejan number (Equations (23) and (27)). The influence of $\mathrm{Re}$, $Br$, $Rd$ and $\mathsf{\Omega }$ on $Ns$ and $Be$ are examined and displayed in Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19 and Figure 20. Figure 14, Figure 15, Figure 16 and Figure 17 examine one of the significant features of this study, i.e., volumetric entropy generation for $Br$ and $\mathrm{Re}$. The influence of $Ns$ becomes increasingly important to all these parameters. Higher $Ns$ and $Be$ values are due to an increase of Br. The higher $Ns$ and $Be$ values are also generated by the role of $\mathrm{Re}$. $Ns$ and $Be$ strongly depend on $\mathrm{Re}$. We observed that increasing $\mathrm{Re}$ also increases $Ns$. As Re increases, hectic motion occurs, the fluid moves more vigorously and thus the impact of heat transfer and fluid friction on $Ns$ and $Be$ tends to increase entropy generation. Figure 18 shows that entropy generation is reduced with increased $Rd$. From Figure 20, it can be seen that $Be$ increases with an increase in $Rd$. From Figure 19 and Figure 20, it can be observed that $Be$ is reduced near the lower plate of the channel where $\mathsf{\Omega }$ is more intense; meanwhile, farther from the plate the drift is reversed due to further contribution from the irreversible heat transfer on $Ns$ and $Be$, which reduces the nearby upper plate of the channel with an increase in $\mathsf{\Omega }$.

Figure 21, Figure 22 and Figure 23 show the influence of $\mathsf{\Omega },\lambda$ and $\varphi$ on $C_{fx}$ and $C_{fy}$. From Figure 21 it can be observed that the unsteady parameter increases $C_{fx}$ and $C_{fy}$. However, this trend is reversed for greater values of $\mathsf{\Omega }$. It can be seen in Figure 22 that the skin friction coefficient reduces for increasing values of $\mathsf{\Omega }$. Figure 23 shows that for increasing values of $\varphi$, the skin friction coefficient increases. From the convergence of the series given in Equation (25), $f(\eta ), \mathrm{g}(\eta ),\theta (\eta )$ depends entirely upon the auxiliary parameters ${\hslash }_f, {\hslash }_g, {\hslash }_{\theta }$ and the so-called $\hslash$-curve. It is selected in such a way that it controls and converges on the series solution. The probable selection of $\hslash$ can be found by plotting $\hslash$-curves of $f^{″}(0),g^{\prime }(0),{\theta }^{\prime }(0)$ for the $20\mathrm{th}$ order approximated HAM solution, as shown in Figure 24 and Figure 25. The valid region of $\hslash$ is $-0.1<{\hslash }_f<0.3,-0.5<{\hslash }_g<0.1,-0.5<{\hslash }_{\theta }<0.1$.

5. Conclusions

The exploration of nanoparticles preparations has introduced more deliberation in mechanical and industrial engineering owing to their probable use for increasing the continuous phase fluid thermal performance of cooling devices. A significant source of renewable energy is thermal radiation, which can be beneficial to govern overall population levels. In the present work, the second law of thermodynamics is applied in terms of the impact of nanoparticles on non-dimensional entropy for rotating flow with suggested thermal radiation. Mathematical modeling is established by modeling five different types of nanoparticles with the purpose of achieving an appropriate mechanism to enhance the thermal conductivity of continuous phase fluid. The following conclusions can be made:

• The unsteadiness parameter decreases the temperature profile and increases the velocity field.
• The thermal boundary layer thickness is reduced for larger values of the rotation rate parameter.
• The heat transfer rate rises for greater values of $Rd$ and ${\theta }_w$.
• With increasing values of $\mathrm{Pr}$, the heat profile $\theta \left(\eta \right)$ reduces.
• The performance of $Be$ is examined for the optimal values of the parameters at which $Ns$ decreases.
• Entropy generation is increased with the increase of $\mathrm{Pr}$, $Ec$ and radiative heat flux.
• Velocity and temperature profiles decrease due to the increased unsteadiness parameter.
• Greater values of $\varphi$ increase the frictional force within the fluid motion.