**Academic Editor:**Muhammad Humayun, Huazhong University of Science and Technology, China.

**Checked for plagiarism:** Yes

**Review by:** Single-blind

Journal of Advances in Nanotechnology

ISSN: Coming Soon

Marangoni Effect in Second Grade Forced Convective Flow of Water Based Nanofluid

Advances in nanotechnology especially in nanofluids comprising of a typical base fluid saturated with nano-size metallic particles with enhanced thermophysical properties is one of the hot topics in industrial as well as engineering applications of applied mathematics. This article explores the impact of Lorentz forces and Maragoni effect on second grade nanofluid forced convective flow. A two phase model is chosen to validate the nanofluid. Series solutions are achieved through HAM after transformation of PDEs into ODEs. The Brownian motion effect, Thermophoresis and the Marangoni effect are the main influencing factors for present flow model. In addition, the influence of pertinent fluid parameters such as Schmidt number, Magnetic number and Prandtl number on the velocity, temperature and concentration profiles is discussed with the help of graphs. With an enhanced Marangoni factor the hydraulic boundary layer thickness shows enhancement.

The concept of nanofluid is very simple but it has broadly improved the efficiency of fluids especially from industrial point of view. The suspension of nanoparticles in typical base fluid results in a very highly conductive fluid that certainly reduces the human effort as well as the machine heating. These suspensions can meet the cooling requirement for any type of thermal system. The foundation of nanofluid was laid by Choi ^{1} and later on, the concept became famous in all the related fields of fluid dynamics due to its industrial applications. Ibanez et al. ^{2} discussed an analytic investigation for MHD nanoliquid flow with suction/injection and radiation effects. Gorji et al. ^{3} involved a variable magnetic field for investigation of a squeezed nanofluid flow considering the Brinkman model. Related articles can be seen in ^{4, }^{5, }^{6, }^{7, }^{8, }^{9, }^{10, }^{11} and references cited therein.

Surface tension as well as the gradients of temperature and concentration results in Marangoni convection appearing in fluid flow analysis. The topic is interesting due to its applications in crystal growth mechanism, beam melting and welding etc. Lin et al. ^{12} reported the Marangoni convection in nanoliquid flow via thermal gradient varying the magnetic effect. Surface tension is taken as a nonlinear temperature and concentration function. Aly and Ebaid ^{13} found exact solutions considering Marangoni effect in nanofluid flow using Laplace transformation method. Most recent articles can be seen in ^{14, }^{15, }^{16, }^{17} and reference cited therein.

Present inspiration is covered by the following novel aspects. Firstly, to model a second grade nanofluid forced convective flow due to the Lorentz forces instigated into the fluid by induction of a variable magnetic field. Secondly, to achieve series solutions through HAM. The PDEs are transformed into ODEs using usual similarity transformations. Finally, to discuss the results through graphs with sufficient concluding remarks.

We consider a second grade nanofluid forced convection due to the Lorentz forces instigated into the model by induction of an applied Magnetic field. The effect of magnetic field develops normal to the surface upon which the fluid is flowing. The Marangoni effect is utilized to apprehend the fluid flow in forward direction. The problem is considered in two dimensions such that the fluid flows along x-axis and the y-axis extends normal to the surface. There is no fluid motion along y-axis therefore, v=0 is taken. Along x-axis the fluid undergoes the Lorentz forces that are induced into the model due to the applied magnetic field neglecting the Hall effects. The interface temperature is taken as function of *x* as depicted in Figure 1. The interface condition and velocity components can be visualized in Figure 1.

Governing equations are therefore, as follows:

……. (2)

……. (3)

……. (4)

with following boundary conditions,

….. (5)

The surface tension σ, being a function of *T* and *C* can be defined as follows:

……… (6)

where

………. (7)

Here u , ν are the horizontal and vertical velocity components, respectively *v *_{fl}. is the kinematic viscosity, *ρ*_{fl} is the density of fluid B_{0} is the magnetic effect involved in the model for MHD, *σ* is the surface tension, *α*_{fl} is is thermal diffusivity of the fluid, the ratio between heat capacity of the nanoparticles (*ρ**c)*_{np} and heat capacity of base fluid (*ρ**c) *_{fl}. *D*_{Br}is Brownian diffusion, *D*_{Th} is Thermophoresis.

Define,

……… (8)

Using (8) in equations ^{1, }^{2, }^{3, }^{4}, we get

…………. (9)

……… (10)

………….. (11)

subject to the following boundary conditions:

……… (12)

where r=(C_{0} 𝜸_{C})/(T_{0} γ_{T} ), M=(L^{2} σ_{l} B_{0})/μ, Pr=ν_{fl} /α_{fl} is the Prandtl number N_{b}=(ρc)_{np} D_{B}_{𝜸} C_{0} x^{2}) / (ρc)_{lf} L^{2} a), N_{t}=(ρc)_{np} D_{Th} x^{2}) / (ρc)_{lf} L^{2} a), L= μν fl)/(σ_{0} T_{0} γ_{T}) are the Marangoni factor (r), modified Hartman number (M), Prandtl factor (Pr), Brownian motion parameter (Nb), Schmidt number (Sc), Thermophoresis parameter (Nt) and reference length (L), respectively. The second grade fluid parameter (α) is defined as α_{1}/(ρ_{fl} L^{2}).

Homotopy analysis method ^{18, }^{19, }^{20, }^{21, }^{22} for convergent series solutions is very convenient to obtain the approximated solutions for a given nonlinear system. The method is independent of small/large physical parameters. Thus it is an efficient method as compared to other conventional methods for solving nonlinear systems. Assuming the following initial guesses,

…. (13)

The auxiliary parameters can be defined as follows:

……. (14)

such that,

……. (15)

where *M*_{i} are constants for *i**=1-7*. Subsequently, the zeroth order problems of deformation are written as follows:

with following boundary conditions:

Therefore,

where p**∈** (0,1) is a typical embedding parameter and h ̂_{f }, h ̂_{θ}, h ̂_{ϕ} are so-called auxiliary parameters with N_{f}, N_{θ}, N_{ϕ} are the non-linear operators. For* p *=0,1 , we have:

……. (19)

The m^{th} problems of deformation are

where g_{m}=1 for m>1 otherwise 0.

Finally,

Thus,

…… (22)

are the general solutions where Mi are the arbitrary constants for i=1-7 and f_{m}^{*} (η) ,θ_{m}^{*} (η), ϕ_{m}^{*} (η) are special solutions.

The auxiliary parameters introduced in (series solution) for the velocity profile (*f*), temperature distribution (θ) and concentration distribution (ϕ) are termed as convergence control parameters. These parameters are significant to speed-up the convergence. Convergence intervals are sketched in Figure 2. The interval of interest for convergence of the aforementioned profiles is (-0.60, 0.30).

This section is concerned with the discussion on graphical results that are obtained through Mathematica based HAM code for specific values of pertinent fluid parameters including the second grade fluid parameter, the Magnetic number, the Marangoni ratio, the Prandtl and Schmidt numbers as well as the Brownian and Thermophoretic parameters. Figure 3. displays the effect of second grade fluid parameter on the velocity profile. The velocity shows decreasing behavior for incremental values of respective parameter. A certain increase in viscosity results in decrease of fluid motion. The velocity shows reduction with augmented values of Marangoni ratio as depicted in Figure 4. Similar is the case noticed for Magnetic number as displayed in Figure 5. A significant drop in velocity profile appears at first sight however, the variation slows down with the stronger effect of magnetic number. The temperature profile, as displayed in Figure 6, shows augmented behavior with augmented values of Marangoni ratio. However, after certain limitation, the effect can be seen in opposite nature. The Prandtl number shows enhancement in temperature profile due to an enhanced thermal diffiusivity as one can see in Figure 7. Figure 8 and Figure 9 are the display of variation in Temperature profile varying the values of Brownian motion parameter and Thermophoretic parameter, respectively. The rapid and in-predictive movement of nanoparticles as well as the more stronger Thermophoretic force results in incremental behavior of temperature profile. The effect of Marangoni ratio on concentration profile is displayed in Figure 10. The results persist with those appearing in Temperature profile for the same parameter. An opposite behavior is seen upon variation of Thermophoretic parameter and Schmidt number for concentration profile as displayed in Figure 11 and Figure 12, respectively.

**Figure 3.** Impact of second grade parameter on velocity profile

**Figure 4.** Impact of Marangoni ratio on velocity profile

**Figure 5.** Impact of magnetic number on velocity profile

**Figure 6.** Impact of Marangoni ratio on temperature profile

**Figure 7.** Impact of Prandtl number on temperature profile

**Figure 8.** Impact of Brownian motion on temperature profile

**Figure 9.** Impact of Thermophoresis on temperature profile

Figure 10. Impact of Marangoni ratio on concentration profile

**Figure 11.** Impact of Thermophoresis on concentration profile

**Figure 11.** Impact of Thermophoresis on concentration profile

**Figure 12.** Impact of Schmidt number on concentration profile

This study concludes with the impact of Marangoni effect and Lorentz force generated by MHD on second grade nanofluid flow. Flow model is formulated mathematically in PDEs which are transformed into ODEs using transformation and HAM is applied to get the convergent series solutions. We conclude that velocity profile reduces for stronger Marangoni effect in second grade nanofluid. The Brownian motion and Thermophoresis have significant impact on the flow profiles. Furthermore, the temperature and concentration profile shows augmented behavior with augmented Marangoni ratio. One can see that due to strong impact of Lorentz forces, a reduction in hydraulic boundary layer is noticed however, an opposite behavior is shown the other two profiles.