Marangoni Boundary Layer Flow and Heat Transfer of Graphene–Water Nanofluid with Particle Shape Effects

Abstract

Graphene nanofluids have attracted the attention of many researchers because of a variety of remarkable properties such as extraordinary electronic transport properties, high thermal conductivity, and large specific surface areas. This paper investigates the shape effects of nanoparticles on the Marangoni boundary layer of graphene–water nanofluid flow and heat transfer over a porous medium under the influences of the suction parameter. The graphene–water nanofluid flow was contained with various shapes of nanoparticles, namely sphere, column, platelet, and lamina. The problem is modeled in form of partial differential equations (PDES) with boundary conditions. The governing transport equations are converted into dimensionless form with the help of some suitable nondimensional variables. The solution of the problem was found numerically using the NDSolve technique of Mathematica 10.3 software. In addition, the numerical solutions were also compared with analytical results. The homotopy analysis method (HAM) is used to calculate the analytical results. The results show that lamina-shaped nanoparticles have better performance on temperature distribution while sphere-shaped nanoparticles are more efficient for heat transfer than other shapes of nanoparticles.

1. Introduction

Due to the abovementioned outstanding thermophysical properties of graphene, it has become a solid candidate in the field of fluid heat transfer. The work of Novoselov et al. [1] in 2004 is pioneering in its discussion of graphene. Graphene-based nanofluids have high heat transfer and thermal conductivity characteristics as compared to other carbon materials [2]. Graphene is mostly used as a nanofluid. A water–graphene nanofluid is preferable because it is a good coolant for real-time applications of heat transfer [3]. Many researchers have worked on graphene nanofluid flow due to its outstanding thermophysical properties. Alay [4] examined the exact solution of water–graphene nanofluid over a shrinking/stretching sheet. Upadhya et al. [5] numerically discussed water–graphene nanofluid flow over a stretching cylinder and also performed a comparison with silver–water nanofluid. Upadhya et al. [6] studied theoretically the unsteady flow of Carreau fluid in a suspension of graphene and dust nanoparticles with Cattaneo–Christov heat flux. Noor et al. [7] discussed an analytical solution of graphene nanofluid and heat transfer in a thin film. Shit and Mukherjee [8] examined graphene–polydimethylsiloxane nanofluid between two plates with the thermal radiation effect. Arif et al. [9] studied the analytical results of Maxwell fluid flow by using molybdenum disulfide and graphene nanoparticles in engine oil. Mankiandan and Baskar [10] discussed the thermal conductivity of water–graphene nanofluid using factorial design experiments.

Very challenging scientific phenomena and problems are displayed by nonlinear ordinary differential equations (ODEs) or partial differential equations (PDEs). For instance, boundary layer fluid flow can be mentioned as an example. The research on several techniques applied for solving nonlinear differential equations (DEs) is a significant topic for analyzing mathematical problems related to engineering. There are various techniques for solving nonlinear equations that range from completely numerical to completely analytical ones. Besides the many advantages of applying numerical techniques, closed results seem more appealing because they explain physical insights through the physical problem. By using an analytical method, parametric studies become more convenient [11]. Some of these techniques include the variation iteration method [12], Galerkin method [13], optimal homotopy asymptotic method [14], homotopy analysis method (HAM) [15], differential transformation method [16], homotopy perturbation method [17], and Adomian decomposition method [18].

Marangoni convection flow appears with the change of surface tension along liquid–gas or liquid–liquid interfaces, which was first observed by Napolitan [19,20]. Moreover, it is also significant in daily life uses such as biomedical and industrial fields, coating flow technology, film drainage in emulsions, the drying of semi-conductor wafers in microelectronics, foams, microfluidics, and surfactant replacement therapy for neonatal infants [19]. Most significantly, Marangoni convection has been demonstrated using the boiling test in microgravity; likewise, it is also significant in the absence of gravity [21,22]. Pop et al. [23] examined Marangoni convection phenomena. Many researchers have discussed Marangoni convection by considering different geometry shapes [23,24,25,26,27].

Arafune and Hirata [28] discussed an analytical solution for the velocity profile of Marangoni convection by considering constant velocity at any point on the free surface. Various researchers studied the different physical effects of Marangoni convection boundary layer flow by taking the thermophysical properties model without any effect of the Prandtl number on base fluid and nanoparticles [29,30,31,32]. Marangoni convection is a significant part of the determination of fluid movement and material processing, especially low gravity hydrodynamics, to analyze mass and heat transfer interaction [33]. Many researchers have worked on heat transfer in Marangoni convection flow. Marangoni convection showed a dominant role played in heat transfer and flow properties of fluid due to the change in liquid surface tension in a temperature gradient field regardless of system size. Marangoni convection and heat transfer in non-Newtonian power-law fluid while imposing linear surface temperature have been explored [34]. The Prandtl number effect was studied by Christopher and Wang [25], and they ascertained the analytical solution of the temperature profile by considering Marangoni convection over a flat surface. Lin and Zheng [35] studied the heat transfer in Cu–water Marangoni boundary layer nanofluid flow over a disk.

Inspired by the above uses of graphene nanofluids, a graphene-water nanofluid was considered in this study. As far as we are aware, studies on the shape effects of graphene nanoparticles have never been done. The shape of nanoparticles is very significant in the enhancement of heat transfer. It is necessary to find the heat transfer under exact nanoparticle shapes [36]. The main objective of this study is to find the shape effects of graphene nanoparticles in Marangoni boundary layer flow and heat transfer. The effect of the suction parameter under the influences of thermal radiation and a magnetic field was also considered. Finally, graphs of velocity, temperature, and the Nusselt number were plotted and all their aspects are debated.

2. Description of the Problem

We consider a laminar, incompressible, steady, and two-dimensional Marangoni boundary layer flow of graphene–water nanofluid over a porous surface with the effects of the suction parameter. It is also considered that nanoparticles and base fluid water are in equilibrium and, also, that no slip occurs between them. It is assumed that flow takes place at y ≥ 0. Furthermore, it is also assumed that surface tension is linear with temperature. A uniform B magnetic field is applied in the normal direction to the surface. The partial governing equations of the problem are modeled as [30]. The thermophysical properties of graphene–water nanofluid are presented in

Table 1. Thermophysical properties of the water and graphene nanoparticles [4].

Physical Properties	Graphene	Pure Water
Cp (J/kg K)	2100	4179
ρ (kg/m3)	2250	997.1
k (W/m K)	2500	0.613
$\sigma$ (Ω m)−1	1 × 107	0.05

. The values of empirical shape factor and sphericity for different shapes of nanoparticles are presented in

Table 2. The values of shape parameters as [37,38].

Shape parameter	Column	Sphere	Platelet	Lamina
$\mathrm{\varphi }$	0.4710	1	0.526	0.1857
m	6.3698	3	5.7	16.1576

.

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

(1)

$$\mathrm{u}\frac{\partial \mathrm{u}}{\partial \mathrm{x}}+\mathrm{v}\frac{\partial \mathrm{u}}{\partial \mathrm{y}}=\frac{{\mathsf{\mu }}_{{\mathrm{n}}_{\mathrm{f}}}}{{\rho }_{{\mathrm{n}}_{\mathrm{f}}}}\frac{{\partial }^2\mathrm{u}}{\partial {\mathrm{y}}^2}-\frac{{\sigma }_{{\mathrm{n}}_{\mathrm{f}}}{\mathrm{B}}^2{u}^2}{{\rho }_{n_f}}-\frac{{\mathsf{\mu }}_{{\mathrm{n}}_{\mathrm{f}}}}{{\sigma }_{{\mathrm{n}}_{\mathrm{f}}}K}\mathrm{u}$$

(2)

$$u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\frac{k_{n_f}}{{\left(\rho Cp\right)}_{nf}}\frac{{\partial }^{2}T}{\partial y^2}-\frac{1}{{\left(\rho Cp\right)}_{nf}}\frac{\partial {\mathrm{q}}_{\mathrm{r}}}{\partial \mathrm{y}}$$

(3)

where $\mathrm{u}$ and $\mathrm{v}$ are components of velocity, along the x axis and y axis, respectively, and $T,{\mathsf{\mu }}_{{\mathrm{n}}_{\mathrm{f}}}{\rho }_{{\mathrm{n}}_{\mathrm{f}}},{\left(\rho Cp\right)}_{nf},{\mathsf{\alpha }}_{{\mathrm{n}}_{\mathrm{f}}}$ are temperature, effective viscosity, effective density, heat capacity, and thermal diffusivity of nanofluid, respectively. The quantities like viscosity, effective density, heat capacitance, and electrically conductivity of nanofluid are the following as in [4].

$$\{\begin{array}{c}{\rho }_{n_f}={\rho }_f\left(1-\varphi \right)+{\rho }_s\varphi ,{\mu }_{n_f}=\frac{{\mu }_f}{{(1-\varphi )}^{2\cdot 5}}, {\left(\rho Cp\right)}_{nf}={\left(\rho Cp\right)}_f\left(1-\varphi \right)+{\left(\rho Cp\right)}_s\varphi \\ {\sigma }_{{\mathrm{n}}_{\mathrm{f}}}=1+\frac{3\varphi \left(\frac{{\sigma }_s}{{\sigma }_f}-1\right)}{\left(\frac{{\sigma }_s}{{\sigma }_f}+2\right)-\varphi \left(\frac{{\sigma }_s}{{\sigma }_f}-1\right)}\end{array}$$

(4)

where $f$,$nf,$ and $s$ are denote the thermophysical properties of the base fluid, nanofluid, and nanoparticle solid, respectively. Furthermore, the Hamilton–Crosser model is used for mixtures of thermal conductivity and is effective for both nonspherical and spherical shapes of dispersed nanoparticles. According to Hamilton–Crosser model, the mixture of thermal conductivity can be calculated from the following formula

$${\mathrm{k}}_{{\mathrm{n}}_{\mathrm{f}}}={\mathrm{k}}_{\mathrm{f}} \left[\frac{{\mathrm{k}}_{\mathrm{s}}+{\mathrm{k}}_{\mathrm{f}}\left(m-1\right)-\left({\mathrm{k}}_{\mathrm{f}}-{\mathrm{k}}_{\mathrm{s}}\right)\left(m-1\right)\mathsf{\varphi }}{{\mathrm{k}}_{\mathrm{s}}+{\mathrm{k}}_{\mathrm{f}}\left(m-1\right)+\left({\mathrm{k}}_{\mathrm{f}}-{\mathrm{k}}_s\right)\mathsf{\varphi }}\right]$$

(5)

where $m,\mathsf{\varphi }$ are the shape factor and sphericity of solid nanoparticle.

We have used the Rosseland approximation for radiation. The radiative heat flux in simplified form is expressed as

$${\mathrm{q}}_{\mathrm{r}}=-\frac{{\mathsf{\sigma }}^{\ast }}{3k^{\ast }}\cdot \frac{\partial {\mathrm{T}}^4}{\partial \mathrm{y}}$$

(6)

where ${\mathsf{\sigma }}^{\ast }$ and $k^{\ast }$ are the Stefan Boltzmann constant and coefficient of mean absorption, respectively. Further expanding ${\mathrm{T}}^4$ in a Taylor series about ${\mathrm{T}}_{\infty }$ and neglecting terms of higher order, we get

$${\mathrm{T}}^4\approx 4{\mathrm{T}}_{\infty }^3\mathrm{T}-3{\mathrm{T}}_{\infty }^4.$$

(7)

Equation (3) reduces to

$$u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\frac{k_{n_f}}{{\left(\rho Cp\right)}_{nf}}\frac{{\partial }^{2}T}{\partial y^2}+\frac{16{\sigma }^{\ast }}{3{\left(\rho Cp\right)}_{nf}{K}^{\ast }}\frac{{\partial }^{2}T}{\partial y^2}\frac{1}{{\left(\rho Cp\right)}_{nf}}\frac{\partial {\mathrm{q}}_{\mathrm{r}}}{\partial \mathrm{y}}.$$

(8)

The boundary conditions subjected to the problems are

$${\mathsf{\mu }}_{{\mathrm{n}}_{\mathrm{f}}}\frac{\partial \mathrm{u}}{\partial \mathrm{y}}=\frac{\partial \sigma }{\partial \mathrm{T}}\frac{\partial T}{\partial \sigma },\mathrm{v}=0,\mathrm{T}=a{x}^2+{\mathrm{T}}_{\infty }\mathrm{at}y=0$$

$$\mathrm{u}=0, \mathrm{T}={\mathrm{T}}_{\infty }\mathrm{at}y\to \infty$$

(9)

where $a=\frac{\Delta T}{L}$. $\Delta T$ and $L$ are constant characteristic temperature and length of surface, respectively.

The stream function $\mathsf{\psi }$, similarity variable $\mathsf{\eta },$ dimensionless temperature $\theta \left(\mathsf{\eta }\right)$ and dimensionless velocity $f\left(\mathsf{\eta }\right)$ are defined as:

$$\{\begin{array}{c}\mathsf{\psi }\left(\mathsf{\eta }\right)={\mathsf{\xi }}_{2}xf\left(\mathsf{\eta }\right),\\ \theta \left(\mathsf{\eta }\right)=\frac{\mathrm{T}-{\mathrm{T}}_{\infty }}{a{x}^2},\\ \eta ={\mathsf{\xi }}_{1}y,\end{array}$$

(10)

where

$${\mathsf{\xi }}_1={\left(\frac{{\sigma }_{0\gamma \alpha {\rho }_f}}{{\mu }_f{}^2}\right)}^{\frac{1}{3}},{\mathsf{\xi }}_2={\left(\frac{{\sigma }_{0\gamma \alpha {\mu }_f}}{{\rho }_f{}^2}\right)}^{\frac{1}{3}},\mathrm{u}=\frac{\partial \mathsf{\psi }}{\partial \mathrm{y}},v=-\frac{\partial \mathsf{\psi }}{\partial \mathrm{x}}$$

where $T$ and $T_{\infty }$ are the reference temperature and nanofluid temperature far away the wall and $\gamma$ is the surface tension defined as

$$\mathsf{\gamma }=\left[1-\left(\mathrm{T}-{\mathrm{T}}_{\infty }\right)\overline{\mathsf{\gamma }}\right]{\mathsf{\gamma }}_0$$

(11)

where ${\mathsf{\gamma }}_0$ and $\overline{\mathsf{\gamma }}$ are surface tension and surface tension change with temperature, respectively.

Substituting Equations (4), (5), and (10) into Equations (1), (2), and (8), the following nonlinear ODES are obtained

$$A_1{f}^{‴}+ A_2\left({f^{\prime }}^2-f^{\prime }{f}^{″}\right)-f^{\prime }\left( A_{3}M+\frac{A_1}{\mathrm{Kp}}\right)=0$$

$$\frac{1}{A_5}\left(\frac{1}{Pr}{A}_4+Rd\right){\theta }^{″}+f{\theta }^{\prime }-2\theta f^{\prime }=0.$$

(12)

The boundary value equation becomes

$$f\left(0\right)=A, f^{″}\left(0\right)=-2{(1-\varphi )}^{2\cdot 5},f^{\prime }(\infty )=0, \theta \left(0\right)=1, \theta (\infty )=0$$

(13)

where $\mathrm{A}$ is constant mass transfer parameter, $A>0$ is for suction parameter $r, M=\frac{H_{0{\sigma }_f}^2}{{\rho }_f}\frac{1}{{\mathsf{\xi }}_1{\mathsf{\xi }}_2}$ denotes the magnetic field, $\mathrm{Kp}=\frac{{\mu }_f}{{\rho }_{f}k}\frac{1}{{\mathsf{\xi }}_1{\mathsf{\xi }}_2}$ denotes the permeability parameter, $Rd=\frac{16{\mathsf{\sigma }}^{\ast }{\mathrm{T}}_{\infty }^3}{3{\mathrm{k}}_{{\mathrm{n}}_{\mathrm{f}}}{k}^{\ast }}$ denotes the radiation parameters, and $\mathrm{Pr}=\frac{{\mathsf{\nu }}_{\mathrm{f}}}{{\mathsf{\alpha }}_{\mathrm{f}}}$ is the prandtl number.

$$\{\begin{array}{l} A_1=\frac{1}{{(1-\varphi )}^{2\cdot 5}} ,A_2=\frac{{\rho }_s}{{\rho }_f}\varphi +\left(1-\varphi \right), A_3=\frac{3\varphi \left(\frac{{\sigma }_s}{{\sigma }_f}-1\right)}{(2+\frac{{\sigma }_s}{{\sigma }_f})-\left(\frac{{\sigma }_s}{{\sigma }_f}-1\right)\varphi }+1,\\ A_4=\frac{\left[k_s+\left(m-1\right)k_f\right]-\left(m-1\right)\varphi \left(k_f-k_s\right)}{\left[k_S+\left(m-1\right)k_f\right]+\varphi \left(k_f-k_s\right)},A_5=\left(1-\varphi \right)+\varphi \frac{{\left(\rho Cp\right)}_s}{{\left(\rho Cp\right)}_f}.\end{array}$$

(14)

The quantity particle of interest Nusselt number $\mathrm{Nu}$ is given as

$$\mathrm{Nu}= \frac{x{q}_w\left(x\right)}{k_f\left[T- T_{\infty }\right]}$$

(15)

where $q_{\mathrm{w}}\left(x\right)$ is the wall heat flux given by

$$q_w\left(x\right)=-k_{n_f}{\left(\frac{\partial T}{\partial z}\right)}_{y=o}$$

(16)

Using Equation (6) in Equation (13), we get

$$\mathrm{Nu}=-(A_4+\mathrm{Rd}\ast \mathrm{Pr}){\mathsf{\xi }}_{1}x{\theta }^{\prime }\left(0\right)$$

(17)

3. Solution Methodology

The initial guesses are defined as

$${ℒ}_{\mathrm{f}}={\mathrm{f}}^{‴}-{\mathrm{f}}^{\prime },{ℒ}_{\mathsf{\theta }}={\mathsf{\theta }}^{″}-\mathsf{\theta }$$

(18)

which satisfied the following properties

$${ℒ}_{\mathrm{f}}\left[{\mathrm{Z}}_1+{\mathrm{Z}}_2{\mathrm{e}}^{\mathsf{\eta }}+{\mathrm{Z}}_3{\mathrm{e}}^{-\mathsf{\eta }}\right]=0,{ℒ}_{\mathsf{\theta }}\left[{\mathrm{Z}}_4{\mathrm{e}}^{\mathsf{\eta }}+{\mathrm{Z}}_5{\mathrm{e}}^{-\mathsf{\eta }}\right]=0$$

(19)

where $\mathrm{Z}$₁, $\mathrm{Z}$₂, $\mathrm{Z}$₃, $\mathrm{Z}$₄, and $\mathrm{Z}$₅ are arbitrary constants.

Let q ∈ [$0$,1] represent an embedding parameter and $\hslash$₁ ≠ 0 and $\hslash$₂ ≠0 are the auxiliary linear operators. The generalized homotopic equations corresponding to Equations (12) and (13)

$$\{\begin{array}{c}\left(1-q\right){ℒ}_{\mathrm{f}}\left[\widehat{\mathrm{f}}\left(\eta ,q\right)-{\mathrm{f}}_0\left(\mathsf{\eta }\right)\left]=q{\hslash }_1{\mathcal{N}}_{\mathrm{f}}\right[\widehat{\mathrm{f}}\left(\mathsf{\eta },\mathrm{q}\right),{\widehat{\mathrm{\theta }}}_0\left(\mathsf{\eta },\mathrm{q}\right)\right]\\ \widehat{\mathrm{f}}\left(0,\mathrm{q}\right)= f_w, {\widehat{\mathrm{f}}}^{″}\left(0,\mathrm{q}\right)=-2{(1-\varphi )}^{2\cdot 5}, \widehat{\mathrm{f}}\left(\infty ,\mathrm{q}\right)=0\end{array}$$

(20)

$$\{\begin{array}{c}\left(1-\mathrm{q}\right){ℒ}_{\mathsf{\theta }}\left[\widehat{\mathsf{\theta }}\left(\mathsf{\eta },\mathrm{q}\right)-{\mathsf{\theta }}_0\left(\mathsf{\eta }\right)\right]=q{\hslash }_2{\mathcal{N}}_{\mathsf{\theta }}\left[ {\widehat{\mathrm{\theta }}}_0\left(\mathsf{\eta },\mathrm{q}\right),\widehat{\mathrm{f}}\left(\mathsf{\eta },\mathrm{q}\right)\right]\\ \widehat{\mathrm{\theta }}\left(0,\mathrm{q}\right)=1, \widehat{\mathrm{\theta }}\left(\infty ,\mathrm{q}\right)=0\end{array}$$

(21)

and

$${\mathrm{N}}_{\mathrm{f}}\left[ \widehat{\mathrm{f}}\left(\mathsf{\eta },\mathrm{q}\right),\widehat{\mathrm{\theta }}\left(\mathsf{\eta },\mathrm{q}\right)\right]=A_1\frac{{\partial }^3 \widehat{\mathrm{f}}\left(\mathsf{\eta },\mathrm{q}\right)}{\partial {\mathsf{\eta }}^3}- A_2\widehat{\mathrm{f}}\left(\mathsf{\eta },\mathrm{q}\right)\frac{{\partial }^2 \widehat{\mathrm{f}}\left(\mathsf{\eta },\mathrm{q}\right)}{\partial {\mathsf{\eta }}^2}+ A_2\frac{{\partial }^2 \widehat{\mathrm{f}}\left(\mathsf{\eta },\mathrm{q}\right)}{\partial {\mathsf{\eta }}^2}\frac{{\partial }^2 \widehat{\mathrm{f}}\left(\mathsf{\eta },\mathrm{q}\right)}{\partial {\mathsf{\eta }}^2}-\left( A_{3}M+\frac{A_1}{\mathrm{Kp}}\right)\frac{\partial \widehat{\mathrm{f}}\left(\mathsf{\eta },\mathrm{q}\right)}{\partial \mathsf{\eta }}=0$$

(22)

$${\mathrm{N}}_{\mathsf{\theta }}\left[ \widehat{\mathrm{\theta }}\left(\mathsf{\eta },\mathrm{q}\right), \widehat{\mathrm{f}}\left(\mathsf{\eta },\mathrm{q}\right)\right]=\frac{1}{A_5}\left[\frac{1}{Pr}{A}_4+Rd\right] \frac{{\partial }^2\widehat{\mathrm{\theta }}\left(\mathsf{\eta },\mathrm{q}\right)}{\partial {\mathsf{\eta }}^2}+\widehat{\mathrm{f}}\left(\mathsf{\eta },\mathrm{q}\right)\frac{\partial \widehat{\mathrm{\theta }}\left(\mathsf{\eta },\mathrm{q}\right)}{\partial \mathsf{\eta }}-2\frac{\partial \widehat{\mathrm{f}}\left(\mathsf{\eta },\mathrm{q}\right)}{\partial \mathsf{\eta }}\widehat{\mathrm{\theta }}\left(\mathsf{\eta },\mathrm{q}\right)=0.$$

(23)

Due to Taylor’s series, with respect to q, we have

$$\mathrm{f}\left(\mathsf{\eta },\mathrm{q}\right)={\mathrm{f}}_0\left(\mathsf{\eta }\right)+{\displaystyle \sum }_{\mathrm{m}=1}^{\infty }{\mathrm{f}}_{\mathrm{m}}(\mathsf{\eta }){\mathrm{q}}^{\mathrm{m}},\mathsf{\theta }\left(\mathsf{\eta },\mathrm{q}\right)={\mathsf{\theta }}_0\left(\mathsf{\eta }\right)+{\displaystyle \sum }_{\mathrm{m}=1}^{\infty }{\mathsf{\theta }}_{\mathrm{m}}(\mathsf{\eta }){\mathrm{q}}^{\mathrm{m}}$$

(24)

where

$${\mathrm{f}}_{\mathrm{m}}\left(\mathsf{\eta }\right)={\frac{1}{\mathrm{m}!}\frac{{\partial }^{\mathrm{m}}\mathrm{f}\left(\mathsf{\eta },\mathrm{q}\right)}{\partial {\mathsf{\eta }}^{\mathrm{m}}}|}_{\mathrm{q}=0},{\mathsf{\theta }}_{\mathrm{m}}\left(\mathsf{\eta }\right)={\frac{1}{\mathrm{m}!}\frac{{\partial }^{\mathrm{m}}\mathsf{\theta }\left(\mathsf{\eta },\mathrm{q}\right)}{\partial {\mathsf{\eta }}^{\mathrm{m}}}|}_{\mathrm{q}=0}$$

(25)

And, thus, higher-order deformation problems are

$$\{\begin{array}{c}{ℒ}_{\mathrm{f}} \left[{\mathrm{f}}_{\mathrm{m}}\left(\mathsf{\eta }\right)-{\mathsf{\chi }}_{\mathrm{m}}{\mathrm{f}}_{\mathrm{m}-1}\left(\mathsf{\eta }\right)\right]={\hslash }_{\mathrm{f}}{ℛ}_{\mathrm{f}}^{\mathrm{m}}\left({\mathrm{f}}_{\mathrm{m}-1}\left(\mathsf{\eta }\right),{\mathsf{\theta }}_{\mathrm{m}-1}\left(\mathsf{\eta }\right)\right)\\ {\mathrm{f}}_{\mathrm{m}}\left(0\right)=0,{\mathrm{f}}_{\mathrm{m}}^{’’}\left(0\right)=0,{\mathrm{f}}_{\mathrm{m}}^{\text{’}}(\infty )=0\end{array}$$

(26)

$$\{\begin{array}{c}ℒ{\mathrm{L}}_{\mathrm{f}} \left[{\mathsf{\theta }}_{\mathrm{m}}\left(\mathsf{\eta }\right)-{\mathsf{\chi }}_{\mathrm{m}}{\mathsf{\theta }}_{\mathrm{m}-1}\left(\mathsf{\eta }\right)\right]={\hslash }_{\mathsf{\theta }}{ℛ}_{\mathsf{\theta }}^{\mathrm{m}}\left({\mathsf{\theta }}_{\mathrm{m}-1}\left(\mathsf{\eta }\right),{\mathrm{f}}_{\mathrm{m}-1}\left(\mathsf{\eta }\right)\right)\\ {\mathsf{\theta }}_{\mathrm{m}}\left(0\right)=0,{\mathsf{\theta }}_{\mathrm{m}}(\infty )=0.\end{array}$$

(27)

Here,

$${\mathsf{\chi }}_{\mathrm{m}}=\{\begin{array}{c}0 \mathrm{when} \mathsf{m} \le 1, \\ 1, \mathsf{m}>1\end{array}$$

(28)

where

$${ℛ}_{\mathrm{f}}^{\mathrm{m}}{\mathrm{f}}_{\mathrm{m}}\left(\mathsf{\eta }\right)=A_1{\mathrm{f}}_{\mathrm{m}-1}^{‴}\left(\mathsf{\eta }\right)- A_2{\displaystyle \sum }_{\mathrm{z}=0}^{\mathrm{m}-1}{\mathrm{f}}_{\mathrm{z}}^{\prime }{\mathrm{f}}_{\mathrm{m}-1-\mathrm{z}}^{″}+ A_2{\displaystyle \sum }_{\mathrm{z}=0}^{\mathrm{m}-1}{\mathrm{f}}_{\mathrm{z}}^{\prime }{\mathrm{f}}_{\mathrm{m}-1-\mathrm{z}}^{\prime }-\left(A_{3}M+\frac{A_1}{\mathrm{Kp}} \right){\mathrm{f}}_{\mathrm{m}-1}^{\prime }\left(\mathsf{\eta }\right)$$

(29)

$${ℛ}_{\mathrm{f}}^{\mathrm{m}}{\mathsf{\theta }}_{\mathrm{m}}\left(\mathsf{\eta }\right)=\frac{1}{A_5}\left[\frac{1}{Pr}{A}_4+Rd\right]{\mathsf{\theta }}_{\mathrm{m}-1}^{″}+{\displaystyle \sum }_{\mathrm{z}=0}^{\mathrm{m}-1}{\mathrm{f}}_{\mathrm{z}}{\mathsf{\theta }}_{\mathrm{m}-1-\mathrm{z}}^{\prime }-2{\displaystyle \sum }_{\mathrm{z}=0}^{\mathrm{m}-1}{\theta }_{\mathrm{z}}{\mathrm{f}}_{\mathrm{m}-1-\mathrm{z}}^{\prime }.$$

(30)

4. Comparison of Results and Discussion

The numerical and analytical solutions of shape effects of nanoparticles in Marangoni convection of graphene–water nanofluid under the influence of suction parameter are discussed here. The significant theme of the present section is to check the dimensionless physical quantities such as velocity and temperature of fluid for various values of pertinent parameters, namely suction parameter, permeability parameter, magmatic field, solid volume fraction, and radiation parameter.

Figure 1 presents the schematic model of the problem. Figure 2 reveals the effect of nanoparticles and also shows that the lamina-shaped nanoparticles have better performance in temperature distribution as compared to other shapes of nanoparticles in nanofluid. The reason is that the lamina-shaped nanoparticles have the lowest viscosity and thermal conductivity as compared to the other nanoparticle shapes, while an opposite performance is noted in sphere-shaped nanoparticles. The impacts of M on velocity and temperature of nanofluids are presented in Figure 3. It is clear from Figure 3 that the velocity of nanoparticles is inversely proportional but that the temperature is directly proportional to M. This is due to Lorentz force which is induced by a magnetic field that opposes the motion of nanofluid and also Lorentz force heating in the energy equation, which assists as an increment in heat to the boundary layer. The lamina nanoparticles act as though they have a dominant role in temperature distribution with the impact of M. The influences of Kp on velocity and temperature profiles are shown in Figure 4. This figure reveals that Kp has direct and inverse relations with velocity and temperature profiles, respectively. The reason is that with the enhancement of Kp, there is an increment in porous medium permeability which decreases the resistance in the nanofluid. As compared to other shapes of nanoparticles, lamina-shaped nanoparticles have a dominant role in temperature distribution under the impact of Kp. The variation effects of φ on velocity and temperature profiles are displayed in Figure 5. The velocity of nanofluid decreases, but the temperature profile intensifies with an increase in φ. The reason is that because of the nanoparticles, the thermal conductivity of the nanofluid increases which, in turn, boosts the thermal boundary layer thickness. The performance of lamina-shaped nanoparticles is dominant in temperature profiles with a variation in φ. Figure 6 describes the influences of A on velocity and temperature distribution and also shows that the velocity and temperature distribution of nanofluid decrease with intensifying of the suction parameter; physically, this happens because of the decreasing wall shear stress. The fact is that the heated nanofluid is pushed toward the wall where the buoyancy forces can increase the viscosity, so it decreases the wall shear stress. Under the influences of A, lamina-shaped nanoparticles have excellent performance in temperature profiles. Figure 7 portrays the potential of Rd on the temperature distribution and demonstrates that Rd has a direct relation with the temperature profile, physically, in the existence of thermal radiation significance increment in radiative heat, which assists in increasing the thermal stat of nanofluid, initiating its temperature to increase. It is also noted form Figure 7 that lamina-shaped nanoparticles show an effective role in temperature distribution. The performance of nanoparticles in the form of heat transfer is presented in Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13. From Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13, it is observed that the heat transfer rate of sphere shape nanoparticles in graphene–water is much better than for other nanoparticles. From Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13, it is noted that sphere > platelet > column > lamina. The accuracy of the numerical results are checked by comparison with the analytical results which are obtained by the homotopy analysis method. In order to provide a comparison of the solution in a graphical form, Figure 14, Figure 15, Figure 16 and Figure 17 were plotted. The validation of results for each nanoparticle shape is presented in

Table 3. Comparison results for sphere-shaped nanoparticles.

$\mathit{\eta }$	HAM Results	Numerical Results	Absolute Error
0	1.000000	1.000000	0.000000
1	0.873762	0.868390	0.005372
2	0.762481	0.745049	0.017432
3	0.655491	0.629456	0.026035
4	0.552687	0.521125	0.031562
5	0.453964	0.419600	0.034365
6	0.358922	0.324453	0.034469
7	0.267190	0.235283	0.031907
8	0.177983	0.511715	0.026268
9	0.088937	0.073397	0.015540
10	0.000045	−2.867620 × 10−16	0.000048

,

Table 4. Comparison results for platelet-shaped nanoparticles.

$\mathit{\eta }$	HAM Results	Numerical Results	Absolute Error
0	1.000000	1.000000	2.22050 × 10−16
1	0.879780	0.875368	0.004413
2	0.77126	0.756973	0.014289
3	0.665913	0.644504	0.021409
4	0.563643	0.537663	0.025970
5	0.464400	0.436170	0.028230
6	0.367944	0.339756	0.028187
7	0.274100	0.248168	0.025932
8	0.182831	0.161164	0.021668
9	0.091886	0.078512	0.013374
10	0.000045	−0.000012	0.0000057

,

Table 5. Comparison results for column-shaped nanoparticles.

$\mathit{\eta }$	HAM Results	Numerical Results	Absolute Error
0	1.00000	1.00000	1.110220 × 10−16
1	0.880803	0.876647	0.0004157
2	0.772721	0.759171	0.013550
3	0.667628	0.647293	0.020335
4	0.565433	0.540745	0.024688
5	0.466112	0.439275	0.026849
6	0.369449	0.342639	0.026810
7	0.275239	0.250608	0.024631
8	0.183575	0.612962	0.020612
9	0.92388	0.079491	0.012897
10	0.000045	−0.000015	0.0000060

and

Table 6. Comparison results for lamina-shaped nanoparticles.

$\mathit{\eta }$	HAM Results	Numerical Results	Absolute Error
0	1.00000	1.00000	1.110220 × 10−16
1	0.886622	0.886708	0.000085
2	0.779944	0.776596	0.003348
3	0.675410	0.669576	0.005834
4	0.573134	0.565560	0.007574
5	0.473222	0.464464	0.008757
6	0.375596	0.366207	0.009272
7	0.279981	0.270708	0.009272
8	0.186039	0.177891	0.008148
9	0.923294	0.087678	0.005616
10	0.000045	−0.000012	0.0000058

. Figure 18 displays the error decay for the 10th-order approximation.

5. Concluding Remarks

The shape effects of four different graphene nanoparticles (sphere, column, platelet, and lamina) on Marangoni boundary layer nanofluid flow and heat transfer with the effects of suction parameter were computed and discussed. On the basis of illustrative discussion, the following deductions are obtained:

• The accuracy of the numerical results has been verified with analytical results.
• The nanoparticles of lamina shape play a significant role in the disturbance of temperature profile.
• The nanoparticles of sphere shape play the smallest role in the disturbance of temperature profile.
• The nanoparticles of sphere shape play a remarkable role in the rate of heat transfer.
• The nanoparticles of lamina shape play the smallest role in the heat transfer rate.