Nanomaterials considered a main factor in industry development. Nanofluids are an important branch of nanomaterials, which were firstly referred by Choi^[1] in $1995$. Nanofluids are identified as a base fluid contains suspended small particles (1$-$100) nm. Water, oil, and alcohols are commonly base fluids. The importance of nanofluids is due to their unusual thermophysical properties. Nanofluids exhibit high ability to conduct electricity and heat, so it plays a vital role in industry. Before long nanofluid components have expanded vital centralization of researchers inferable from their entrancing warm transport in a couple of calm disapproved of fields. There are many applications such as engine cooling, cooling of electronics, refrigeration, solar water heating, microprocessors, laser applications and super conducting magnets. Jang and Choi^[2] discussed the role of Brownian motion in the enhanced thermal conductivity of nanofluids. After that Brownian dynamics simulation to determine the effective thermal conductivity of nanofluids is reported by Bhattacharya $et$ $al.$^[3] Buongiorno^[4] presented complete model to analyze the aspects of thermophoresis and Brownian motion. Some continuous explores on nanofluid stream subject can be directed through the examinations.^[5-30]

Flow of liquid on stretching sheet is now massively acknowledge by the researchers because of their large engineering and industrial application in rubber sheets, manufacture of food, glass fiber, hot rolling, paper production and many others. However liquid flow due to non-linear stretching sheet is scare. Rahimi $et$ $al.$^[31] considered collocation method to explore the solutions of an Eyring-Powell fluid caused by linear stretching sheet. Combined properties of viscous dissipation and MHD on the micropolar nanofluid on stretching sheet have been examined by Hsiao.^[32] Zhang $et$ $al.$^[33] analyzed the unsteady flow of Oldroyd-B nanofluid because of stretching sheet. Hayat $et$ $al.$^[34] investigated stretched flow of Jeffery fluid. Recently researchers have investigated the flow of non-Newtonian and Newtonian liquids over non-linear stretching sheet. Seth and Mishra^[35] used the Navier's slip condition to study the transient flow of nanofluid past a non-linear stretching sheet. Hussain $et$ $al.$^[36] studied the characteristics of tangent hyperbolic fluid along non-linear stretching sheet by using convective boundary conditions. Hayat $et$ $al.$^[37] explored the non-linear stretched flow of second grade fluid. Nanofluid flow with variable thickness comprising electrical MHD in the non-linear stretched sheet is discussed by Daniel $et$ $al.$^[38]

Magnetohydrodynamic is the study of magnetic behavior in electrically conducting fluids. Plasmas, salt water and liquid metals are examples of such fluids. In Physics, Hannes Alfvén achieved Nobel Prize in $1970$ for his great work on MHD. Magnetohydrodynamic is important in astrophysics, space plasma physics, cancer tumor treatment, solar physics, blood pump machine, and laboratory plasma experiments. Ishak $et$ $al.$^[39] considered MHD flow past a radially shrinking or stretching disk. Huang and Liu^[40] attempted to combine fluid hammer effect with MHD effect. Hayat $et$ $al.$^[41] computed numerical results for MHD flow with Soret and Dufour effects. Hayat $et$ $al.$^[42] studied three-dimensional flow due to exponentially stretching surface in the existence of an applied magnetic field and Joule heating effects. Second grade nanofluid with MHD over a nonlinear stretching surface is studied by Hayat $et$ $al.$^[43] Tamoor $et$ $al.$^[44] discussed MHD Casson flow between stretching cylinder.

The prime purpose of present topic is to illustrate hydromagnetic 3D stream of Prandtl liquid^[45-48] inside seeing nanoparticles. Thermal and mass trade properties are portrayed through random spread and thermophoresis. Prandtl liquid is taken driving through uniform associated alluring field. Thermal convective condition^[49-50] and a condition related with zero nanoparticles change^[51-52] are completed at the farthest point. The obtained nonlinear differential systems are solved through optimal homotopic analysis method (OHAM).^[53-59] Effects of a couple of physical variables are inspected. In addition the coefficients of surface drag and warmth conversion standard are explored graphically.

We inspect steady hydromagnetic three-dimensional (3D) flow of Prandtl nanoliquid by a linear deformable surface. Thermal condition and as of late made necessity are requiring zero nanoparticles movement are constrained at the point of confinement. Brownian advancement and thermophoretic effects are investigated. The fluid is assumed to be conducted electrically with magnetic field $B_{0}$ applied parallel to $z$-direction. For very small Reynolds number\ current hall and magnetic field effects are ignored. Cartesian coordinate system is incorporated. The sheet is stretched along $x$- and $y$-directions at $z=0$ with velocities $U_{w}$ and $V_{w}$. Boundary layer expressions governing the flow of Prandtl nanofluid in the absence of viscous dissipation and thermal radiation are written as follows:^[16,48]

Here one has the following conditions:^[16,48]

Here $u$, $v$, and $w$ speak to the speeds in $x$-, $y$-, and $z$-bearings, $\mu$ the dynamic consistency, $\nu\left( =\mu /\rho _{f}\right)$ the kinematic thickness, $k$ the warm conductivity, $\rho _{f}$ the thickness, $A$ and $c$ the material constants of Prandtl fluid model,$\sigma $ the electrical conductivity, $\alpha ^{\star }=k/(\rho c)_{f}$ the warm diffusivity, $(\rho c)_{f}$ the warmth capability of the fluid, $(\rho c)_{p} $ the powerful warmth capability of nanoparticles, $T$ the temperature, $D_{B}$ the Brownian development, $C$ the focus, $D_{T}$ the thermophoretic dispersion and $T_{\infty }$ and $C_{\infty }$ the encompassing liquid temperature and fixation. Considering

Expression $(1)$ is naturally fulfilled and Eqs. (2)$-$(7) have the accompanying structures

Here $\beta _{1}$ remains for Prandtl fluid number, $\beta _{2}$ for flexible number, $Ha$ for Hartman parameter, $\alpha $ for ratio parameter, $%Pr$ for Prandtl number, $\gamma $ for Biot number, $N_{b}$ for Brownian development parameter, $N_{t}$ for thermophoresis number and $Sc$ for Schmidt number. These parameters are characterized by:

The physical quantities are given by

It is seen that mass motion spoken to by Sherwood number is presently indistinguishably evaporates and $Re_{x}$=$U_{w}x/\nu$ and $Re_{y}$=$V_{w}y/\nu$ delineate nearby Reynolds parameters. It is also noticed that the Prandtl fluid model reduces to viscous fluid case when $\beta_{1}=1$ and $\beta _{2}=0.$

It has been noted that Eqs. (9)--(12) along with boundary conditions (13) and (14) are four non-linear ordinary differential equations whose optimal series arrangements have been developed by employing OHAM. The initial deformations ($f_{0}$, $g_{0}$, $\theta _{0}$, $\phi _{0}$) and auxiliary linear operators ($\mathcal{L}_{f}$, $\mathcal{L}_{g}$, $\mathcal{L}_{\theta},\mathcal{L}_{\phi }$) are

The above linear operators obey

in which $H_{j}^{\star \ast }$ ($j$ = 1$-$10) illustrate the arbitrary constants. The zeroth and $m$-th order deformation problems can be easily defined in view of above operators via BVPH$2.0$ of software $Mathematica$.

We have unwound the power, essentialness and center verbalizations with the help of BVPh2.0. These verbalizations contain cloud factors $\hbar _{f}$, $\hbar _{g}$, $\hbar _{\theta }$, and $\hbar _{\phi }$. We can process the base estimation of these elements by taking total mix-up pretty much nothing. In the packaging of HAM, these elements expect a basic employment.That is the reason these variables insinuate as association control parameter, which shifts HAM from other illustrative conjecture systems. With a particular ultimate objective to diminish the CPU time, we have used typical waiting errors at the $m$-th order of theory which are described by

Here $\mathcal{N}_{f}$, $\mathcal{N}_{g}$, $\mathcal{N}_{\theta }$, and $\mathcal{N}_{\phi }$ denote the non-linear operators corresponding to Eqs. (9)--(12) respectively. Following Liao:^[53]

where $\varepsilon_{m}^{t}$ indicates add up to leftover squared blunder,$k=20$ and $\delta \zeta=0.5$. Ideal information for assistant parameters at second request of approximations is $\hbar_{f}=-1.591 89$, $\hbar_{g}=-3.056 54$,$\hbar_{\theta }=-1.365 14$, $\hbar _{\phi }=-1.186 26$ and $\varepsilon_{m}^{t}=9.39\times10^{-4}$. Figure 2 speaks to related aggregate remaining mistake plot. Table 1 illustrates normal square residual errors. It has been dissected that the normal averaged square errors decrease with higher request disfigurements.

This section researches impacts of two or three significant physical stream factors like Prandtl liquid parameter $\beta _{1}$, adaptable parameter $\beta _{2}$, Hartman number $Ha$, extent number $\alpha $, Biot parameter $\gamma $, Prandtl parameter $Pr$, Schmidt parameter $Sc$, Brownian improvement parameter $N_{b}$ and thermophoresis number $N_{t}$ on temperature $\theta (\zeta )$ and focus $\phi (\zeta )$. Figures 2 and 3 are constructed to present $\theta (\zeta )$ for different estimations of $\beta _{1}$ and $\beta _{2}$. It is noted from these figures that increase in $\beta _{1}$ and $\beta _{2}$ leads to decrease in temperature. Figure 4 displays the variations of Hartman number $Ha$ on temperature profile $\theta (\zeta )$. Lorentz force arises in $Ha$ that resists the fluid motion therefore temperature field $\theta (\zeta )$ enhances. Figure 5 demonstrates that an adjustment in extent number $\alpha $ prompts a poor temperature $\theta (\zeta )$ and less layer of warm. Impact of Biot number $\gamma $ on $\theta (\zeta )$ is depicted in Fig. 6. Increase in $\gamma $ causes a powerful convection that display an increment in $\theta (\zeta )$. Figure 7 shows that temperature diminish for greater values of Prandtl number. As greater $Pr$ corresponds to lower thermal diffusivity $\alpha $ which causes decrease in temperature. Figure 8 is constructed to study the influence of thermophoresis parameter on the temperature field. This figure illustrates that increase in thermophoresis $N_{t}$ parameter tends to higher temperature. This parameter is occurred due to nanomaterials. The existence of nanomaterials raised the thermal conductivity of nanoliquids. Nanofluid thermal conductivity is an increasing function of temperature. That is why enhancement in temperature is observed for greater estimation of $N_{t}$. Figures 9 and 10 elucidate that nanoparticles concentration is smaller for greater values of $\beta _{1}$ and $\beta _{2}$ (material parameters). Figures 11 and 12 are plotted to analyze the change in $\phi \left( \zeta \right) $ for larger Hartman number and extent parameter $\alpha $. We observed that increasing and decaying impacts occur for both dimensionless parameters on concentration profile. Figure 13 shows the consequences of Schmidt number on $\phi \left( \zeta \right) .$ Schmidt number relates to the mass diffusion of a system. As $Sc$ is increased mass diffusion decreases due to which concentration shows decreasing trend. Brownian parameter $N_{b}$ when increased causes a change in the Brownian motion of nanoparticles which reduces the distribution of concentration as depicted by Fig. 14. Increasing $N_{t}$ causes increase in thermal conductivity of the system which contributes in increase of concentration as seen in Fig. 15. Figure 16 presents impact of $Ha$ and $\beta _{1}$ on $C_{f}Re_{x}^{1/2}$. It has been seen that $C_{f}Re_{x}^{1/2}$ improves for $Ha$. Figure 17 demonstrates the effects of $\alpha $ and $\beta _{1}$ on $C_{f}Re_{x}^{1/2}$. Obviously $C_{f}Re_{x}^{1/2}$ demonstrates expanding conduct for $\alpha$ and $\beta _{1}$. Figure 18 demonstrates the impacts of $Ha$ and $\beta _{1}$ on $C_{g}Re_{y}^{1/2}$. An upgrade in $Ha$ indicates expanding conduct for $C_{g}Re_{y}^{1/2}$. Figure 19 demonstrates the impacts of $\alpha$ and $\beta_{1}$ on $C_{g}Re_{y}^{1/2}$. From this Figure it has been broke down that $C_{g}Re_{y}^{1/2}$ is a hoisting capacity of $\alpha$. Impact of $N_{b}$ and $N_{t}$ on $Nu_{x}Re_{x}^{1/2}$ are uncovered through Fig. 20. Here $Nu_{x}Re_{x}^{1/2}$ diminishes for $N_{t}$ while steady pattern is seen for $N_{b}$. Table 2 shows the comparison for different values of $\alpha $ with homotopy perturbation method (HPM) and exact solutions. Table 2 presents an excellent agreement of OHAM solutions with the existing homotopy perturbation method (HPM) and exact solutions in a limiting sense.

Here hydromagnetic 3D limit layer stream of Prandtl nanoliquid as a result of straightly deformable surface with convective surface condition is performed. Genuine consequences of the current analysis are sketched out as seeks after:

$\bullet$ Both temperature $\theta (\zeta )$ and fixation $\phi (\zeta )$ fields show decaying design for higher Prandtl liquid $\beta _{1}$ and adaptable $\beta _{2}$ parameters.

$\bullet$ An expansion in Hartman number $Ha$ demonstrates more grounded temperature $\theta (\zeta )$ and fixation $\phi (\zeta )$ fields.

$\bullet$ Higher proportion number $\alpha $ delineate lessening conduct for concentration $\phi (\zeta )$ and temperature $\theta (\zeta )$ fields.

$\bullet$ Higher Biot number $\gamma $ indicates more grounded temperature $\theta (\zeta )$ field.

$\bullet$ Similar behavior is observed for different values of $N_{t}$ on concentration $\phi (\zeta )$ and temperature $\theta (\zeta )$ fields.

$\bullet$ For higher estimations of Prandtl parameter $Pr$, temperature $\theta(\zeta )$ decreases.

$\bullet$ An increment in Schmidt number $Sc$ yields weaker Concentration $\phi(\zeta )$ field.

Concentration $\phi (\zeta )$ field exhibits decaying trend via Brownian advancement number $N_{b}$.