In this article, the three-dimensional Magnetohydrodynamics flow of a nanofluid over a horizontal non-linearly stretching sheet in bilateral directions under boundary layer approximation is addressed. A two-phase model has been used for the nanofluid. The influences of thermophoresis, Brownian motion and thermal radiation on heat and mass transfers are considered. Two different cases for the heat and mass transfers are studied. In the first case, uniform wall temperature and zero nanoparticles flux due to thermophoresis are considered. In the second case, prescribed heat and mass fluxes at the boundary are considered. By using the appropriate transformations, a system of non-linear partial differential equations along with the boundary conditions is transformed into coupled non-linear ordinary differential equations. Numerical solutions of the self-similar equations are obtained using a Runge–Kutta method with a shooting technique. Our results for special cases are compared with the available results in the literature, and the results are found to be in good agreement. It is observed that the pertaining parameters have significant effects on the characteristics of flow, heat and mass transfer. The results are presented and discussed in detail through illustrations.
- Brownian motion
- thermal radiation
- power-law index
Boundary layer flow over a stretching surface has various trade and technical applications, such as metal and polymer extrusion, paper production, and many others.
Also, it has important applications in engineering and manufacturing processes. Such processes are wire and fibre covering, foodstuff dispensation, heat-treated materials travelling between a feed roll and a wind-up roll or materials manufactured by extrusion, glass fibre production, cooling of metallic sheets or electronic chips, crystal growing, drawing of plastic sheets and so on. The quality of the final product with the desired features in these processes depends on the degree of freezing in the process and the process of stretching.
Flow past a linearly stretching sheet was first studied analytically by Crane (1970). Gupta and Gupta (1977) mentioned that in all realistic situations, the stretching mechanism is not linear. Later on, several investigators, Hamad and Pop (2010), Turkyilmazoglu (2013), Pal and Mandal (2015), Waqas et al. (2016), studied such type of flow problems either with heat or with mass transfer effects.
The main characteristic of nanofluid (a combination of nanoparticles having size 1–50 nm and liquid) is to improve the base fluid's thermal conductivity. Usually, ethylene glycol, oil, water, and so on are taken as base fluids that have limited heat transfer properties due to their low thermal conductivity. By adding nanoparticles into the base fluid, Masuda et al. (1993) first detected the change in viscosity and thermal conductivity of such fluids. The abnormal increase in the thermal conductivity occurs by the dispersion of nanoparticle, which was first observed by Choi and Eastman (1995). In case of nanofluids, Buongiorno (2006) studied the convective heat transfer and indicated that the Brownian motion, thermophoresis as the most important means for the unpredicted heat transfer augmentation. Later on, different aspects of nanofluid flow were addressed by several researchers. Considering heat and mass transfer, Das (2012) pointed out the slip effects in nanofluid flow. By the finite element method, Rana and Bhargava (2012) studied the nanofluid flow towards a non-linearly stretched surface. Nadeem and Lee (2012) studied the steady nanofluid flow over an exponentially stretched sheet. Considering the convective boundary condition (CBC), Mustafa et al. (2013) discussed nanofluid flow over an exponentially stretched surface. Using the Jeffrey's model for non-Newtonian fluid, Nadeem et al. (2014) investigated nanofluid flow over a stretching surface to obtain a numerical solution. Magnetohydrodynamics (MHD) plays a great role in several areas particularly in geophysics, medicine, engineering and many others [Hayat et al. (2016)]. Magneto nanofluids have significant applications in drug delivery, cancer therapy, enhancement of magnetic resonance and magnetic cell separation (Hayat et al. 2016). In 2013, Khan et al. analysed the effects of a magnetic field, applied externally on nanofluid flow due to a stretching sheet considering a non-similar solution of the problem. Influences of the magnetic field on Cu-water flow were analysed by Sheikholeslami et al. (2014). Using the Buongiorno's model, in 2016, Sheikholeslami and his co-authors described the unsteady flow of nanofluid in the presence of a magnetic field. Ma et al. (2019) studied the consequences of the magnetic field on forced convection Ag-MgO nanofluid flow in a channel with heat transfer. Many other research papers on nanofluid with magnetic field effects can be found in the studies conducted by Ma et al. (2019) and Bhatti et al. (2018).
Thermal radiation affects the heat transfer characteristics significantly particularly in control of heat transfer, space technology and high-temperature processes (Das et al. 2014). Hayat et al. (2016) reported the combined effects due to Joule heating and solar emission on MHD thixotropic nanofluid flow. Radiative flows in companies with the external magnetic fields have wider applications in nuclear engineering, power technology, astrophysics, power generation and so on (Hayat et al. 2016).
The investigations described above are two-dimensional. Three-dimensional problem being more realistic and has already drawn the attention of several researchers. Wang (1984) first addressed the three-dimensional flow due to a bidirectional stretching sheet. Ariel (2007) obtained the solution for a three-dimensional flow caused by a stretching surface using a perturbation technique. Liu and Andersson (2008) analysed the heat transport of a 3-D flow due to a stretching surface. Ahamad and Nazar (2011) discussed the MHD 3-D flow of a viscoelastic fluid due to a stretching surface. In 2013, Hayat and his co-authors considered 3-D Oldroyd-B fluid flow past a stretching surface with a CBC. Considering the exponential temperature distribution, 3-D flow and heat transfer of a viscous fluid were investigated by Liu et al. (2013). Nadeem et al. (2013) analysed the 3-D Casson fluid flow due to a linearly stretching sheet in the presence of a magnetic field. Khan et al. (2014) obtained analytical and numerical solutions for a three-dimensional flow caused by a non-linearly stretching sheet. Khan et al. (2015) investigated three-dimensional flow of a nanofluid past a non-linearly stretching sheet and highlighted its applications to the solar energy industry. Existence of two solutions was reported by Raju et al. (2016) for flow due to a non-linearly stretching sheet. Considering thermal radiation and heat generation/absorption, Shehzad et al. (2016) analysed 3-D MHD flow of a Jeffrey nanofluid over a bidirectional stretching sheet.
The combined effects of MHD and thermal radiation on three-dimensional non-linear convective Maxwell nanofluid flow owing to a stretching sheet were reported by Hayat et al. (2017). The literature review reveals the fact that there is not so much work has been done on three-dimensional nanofluid flow caused by a non-linearly stretching sheet.
Motivated by this, in this article, we investigate the radiative heat transport of the MHD flow of a nanofluid past a stretching sheet with power-law velocity in two directions. The mathematical representations of the nanofluid model proposed by Buongiorno (2006) reflecting the contribution of nanoparticle volume fraction were used. This model also includes the effects of thermophoresis and the Brownian motion. The main objective here is to analyse the heat transport features of the flow of a nanofluid owing to the collective influences of magnetic field and radiation. Two different cases are considered: In Case I, uniform surface temperature is considered. Due to thermophoresis, a flux of nanoparticles volume fraction at the sheet is assumed to be zero. This assumption controls the volume fraction of nanoparticles submissively rather than actively (Mansur and Ishak 2013; Nield and Kuznetsov 2014); on the other hand, in Case II, heat and mass fluxes are considered at the sheet. Suitable similarity transformations are considered to transform the highly non-linear partial differential equations into ordinary ones. Numerical results are obtained by a Runge–Kutta method with a shooting technique. For a clear understanding of the flow and heat transport characteristics, the obtained results are analysed in detail with their physical interpretation through Tables and Figures. It is observed that the physical parameters affect the heat and mass transport significantly.
Let us consider a 3-D steady boundary layer nanofluid flow past a non-linearly stretching surface. In a usual notation, the governing equations of the problem can be written as
Using the Rosseland approximation, the expression for radiative heat flux can be written as
The boundary conditions can be written as
There are different types of heating process that specify the temperature distribution on the surface and ambient fluid. These are constant or prescribed wall temperature (CWT or PWT); constant or prescribed surface heat flux (CHF or PHF), CBC and Newtonian heating. In this problem, we have considered two cases of temperature distributions: CWT and CHF separately. Moreover, different methods are available for the active or passive control of nanoparticles, namely, constant or prescribed concentration at the wall (CHF or PWC); zero mass or normal flux of nanoparticles at the surface (ZMF); constant or prescribed mass flux (CMF or PMF). In this problem, two different boundary conditions, ZMF and CMF for mass transfer, are considered.
Therefore, the conditions at the boundary for temperature, volume fraction for nanoparticles are
Let us define the similarity variable and similarity transformations as
Using these transformations in Eqs (1)–(5) and also using the relation given by Eq. (6), we get the following self-similar equations
The boundary condition becomes
The differential Eqs (11)–(14) are coupled and non-linear. Equations (11)–(14), together with the conditions (15a–15d), form a two-point boundary value problem. By converting it into an initial value problem, those equations are numerically solved using the Runge–Kutta method and the shooting technique by choosing an appropriate value for
To test the validity of our numerical results, an assessment is made with the available results from the open literature for some special cases and presented in Table 1. Table 1 shows a complete agreement of our results with those of Khan et al. (2015) and Raju et al. (2016) in the absence of the magnetic field and thermal radiation.
Comparisons for the values of
Grid independence test has been performed and presented through Table 2, which shows that the values of
Comparisons of the values of
To obtain in-depth features, the numerical results are presented graphically in Figures 2–12.
Figures 2(a–f) depict the effects of magnetic field on velocity components along with the directions of
From Figures 3(a–f), it is clear that the fluid velocity in the direction of
The effect of
From Figures 5(a–d), it is seen that the temperature [Figures 5(a,b)] and nanoparticle volume fraction [Figures 5(c,d)] are increasing functions of thermophoresis parameter N
The influences of radiation on temperature and nanoparticle volume fraction are portrayed in Figures 6(a–d). It is observed that the temperature [Figures 6(a,b)] and the nanoparticle volume fraction [Figures 6(c,d)] both increase for a rise in R. As
The effects of the Prandtl number
Almost negligible effect in the temperature is seen due to changes in the Brownian motion parameter
To provide a clear understanding of the flow field, the streamlines are presented for several values of
Figures 11(a,b) indicate that the coefficient of the skin friction (local) in the directions of
The effects of N
Table 3 shows that −
Values of −
The current study provides solutions (numerical) for a steady MHD three-dimensional nanofluid flow past a bi-directional non-linearly stretching surface. Two types of conditions at the boundary for temperature and volume fraction for nanoparticles are considered. From the current investigation, the following observations are derived:
A rise in the magnetic strength causes a reduction in the resistance and the rate of heat transport. Both the thermal and the concentration boundary layer thicknesses increase with an increase in the magnetic parameter in addition to the radiation parameter. The volume fraction of the nanoparticles is a decreasing function of the Brownian motion parameter. In the presence of thermophoresis, both temperature and volume fraction of the nanoparticles increase significantly in the boundary layer.
A rise in the magnetic strength causes a reduction in the resistance and the rate of heat transport.
Both the thermal and the concentration boundary layer thicknesses increase with an increase in the magnetic parameter in addition to the radiation parameter.
The volume fraction of the nanoparticles is a decreasing function of the Brownian motion parameter.
In the presence of thermophoresis, both temperature and volume fraction of the nanoparticles increase significantly in the boundary layer.
Comparisons for the values of f″(0) and g″(0) in the absence of magnetic field
Comparisons of the values of f″(0) and g″(0) in the absence of the magnetic field for different grid sizes
Values of −θ′(0) and −ϕ′(0) for different values of thermophoresis parameter Nt and Prandtl number Pr with n = 4.0, λ = 0.9, Nb = 0.7, Pr = 2.0, R = 0.5, Sc = 15.0