Unsteady MHD free convection boundary-layer flow of a nanofluid along a stretching sheet with thermal radiation and viscous dissipation effects

Background

Magnetohydrodynamics (MHD) boundary-layer flow of nanofluid and heat transfer over a linearly stretched surface have received a lot of attention in the field of several industrial, scientific, and engineering applications in recent years. Nanofluids have many applications in the industries since materials of nanometer size have unique chemical and physical properties. With regard to the sundry applications of nanofluids, the cooling applications of nanofluids include silicon mirror cooling, electronics cooling, vehicle cooling, transformer cooling, etc. This study is more important in industries such as hot rolling, melt spinning, extrusion, glass fiber production, wire drawing, and manufacture of plastic and rubber sheets, polymer sheet and filaments, etc.

Sakiadis [ ¹ ] was the first author to analyze the boundary-layer flow on a continuous surface. Crane [ ² ] obtained an exact solution of the boundary-layer flow of the Newtonian fluid caused by the stretching of an elastic sheet moving in its own plane linearly. Gorla et al. [ ³ , ⁴ ] solved the non-similar problem of free convective heat transfer from a vertical plate embedded in a saturated porous medium with an arbitrarily varying surface temperature. Cheng and Minkowycz [ ⁵ ] also studied free convection from a vertical flat plate with applications to heat transfer from a dick.

Dissipation is the process of converting mechanical energy of downward-flowing water into thermal and acoustical energy. Various devices are designed in streambeds to reduce the kinetic energy of flowing waters, reducing their erosive potential on banks and river bottoms. Vajravelu and Hadjinicalaou [ ⁶ ] analyzed the heat transfer characteristics over a stretching surface with viscous dissipation in the presence of internal heat generation or absorption.

Takhar et al. [ ⁷ ] studied the radiation effects on the MHD free convection flow of a gas past a semi-infinite vertical plate. Ghaly [ ⁸ ] considered the thermal radiation effect on a steady flow, whereas Rapits and Massalas [ ⁹ ] and El-Aziz [ ¹⁰ ] analyzed the unsteady case. Sattar and Alam [ ¹¹ ] presented unsteady free convection and mass transfer flow of a viscous, incompressible, and electrically conducting fluid past a moving infinite vertical porous plate with thermal diffusion effect. Na and Pop [ ¹² ] analyzed an unsteady flow due to a stretching sheet. In the case of unsteady boundary-layer flow, Singh et al. [ ¹³ ] investigated the thermal radiation and magnetic field effects on an unsteady stretching permeable sheet in the presence of free stream velocity. The study of convective instability and heat transfer characteristics of nanofluids was considered by Kim et al. [ ¹⁴ ]. Jang and Choi [ ¹⁵ ] obtained nanofluid thermal conductivity and the effect of various parameters on it.

The natural convective boundary-layer flows of a nanofluid past a vertical plate have been described by Kuznetsov and Nield [ ¹⁶ , ¹⁷ ]. In this model, the Brownian motion and thermophoresis are accounted with the simplest possible boundary conditions. They also studied the Cheng-Minkowycz problem for natural convective boundary-layer flow in a porous medium saturated by a nanofluid. Bachok et al. [ ¹⁸ ] have shown the steady boundary-layer flow of a nanofluid past a moving semi-infinite flat plate in a uniform free stream. It was assumed that the plate is moving in the same or opposite directions to the free stream to define the resulting system of non-linear ordinary differential equations.

Khan and Pop [ ¹⁹ , ²⁰ ] formulated the problem of laminar boundary-layer flow of a nanofluid past a stretching sheet. They also expressed free convection boundary-layer nanofluid flow past a horizontal flat plate. Hamad and Pop [ ²¹ ] discussed the boundary-layer flow near the stagnation-point flow on a permeable stretching sheet in a porous medium saturated with a nanofluid. Hamad et al. [ ²² ] investigated free convection flow of a nanofluid past a semi-infinite vertical flat plate with the influence of magnetic field. Very recently, Hady et al. [ ²³ ] investigated the effects of thermal radiation on the viscous flow of a nanofluid and heat transfer over a non-linearly stretching sheet.

However, the aim of the present work is to study the unsteady free convection boundary-layer nanofluid flows along a stretching surface with the influence of magnetic field and radiation effect. An explicit finite difference procedure [ ²⁴ ] has been taken to solve the obtained non-similar equations with stability and convergence analysis.

Methods

Presentation of the hypothesis

An unsteady two-dimensional MHD free convection laminar boundary-layer flow of a viscous incompressible and electrically conducting nanofluid along a vertical stretching sheet under the influence of thermal radiation and viscous dissipation is considered. The sketch of the physical configuration and coordinate system is shown in Figure ¹ . Introducing the Cartesian coordinate system, the x -axis is taken along the stretching sheet in the vertically upward direction, and the y -axis is taken as normal to the sheet. Two equal and opposite forces are introduced along the x -axis so that the sheet is stretched, keeping the origin fixed.

Figure 1

Instantaneously at time t  > 0, the temperature of the plate and the species concentration are raised to T _w (>  T _∞ ) and C _w (>  C _∞ ), respectively, which are thereafter maintained constant, where T _w and C _w are the temperature and species concentration at the wall, respectively, and T _∞ and C _∞ are the temperature and species concentration far away from the plate, respectively.

A strong magnetic field is applied in the y direction. The uniform magnetic field strength (magnetic induction) B ₀ can be taken as B  = (0,  B ₀ , 0). The Rosseland approximation is used to describe the radioactive heat flux q _r in the energy equation. Under the above assumptions and the usual boundary layer approximation, the MHD free convection unsteady nanofluid flow and heat and mass transfer with the radiation effect are governed by the following equations:

• The continuity equation

  ∂u∂x+∂v∂y=0

• The momentum equation

  ∂u∂t+u∂u∂x+v∂u∂y=υ∂2u∂y2+gβT-T∞-σB02uρ

• The energy equation

  ∂T∂t+u∂T∂x+v∂T∂y=kρcp∂2T∂y2-1ρcp∂qr∂y+υcp∂u∂y2+τDB∂T∂y·∂C∂y+DTT∞∂T∂y2

• The concentration equation

  ∂C∂t+u∂C∂x+v∂C∂y=DB∂2C∂y2+DTT∞∂2T∂y2

The initial and boundary conditions are

t=0,u=Dx,v=0,T=T∞,C=C∞everywheret≥0,u=0,v=0,T=T∞,C=C∞atx=0u=Dx,v=0,T=Tw,C=Cwaty=0u=0,v=0,T→T∞,C→C∞aty→∞,

where u and v are the velocity components in the x and y directions, respectively, v is the kinematic viscosity, k is the thermal conductivity, D _B is the Brownian diffusion coefficient, D _T is the thermophoresis diffusion coefficient, D (> 0) is the stretching constant, g is the acceleration due to gravity, ρ is the density of the fluid, and c _p is the specific heat at constant pressure.

The Rosseland approximation [ ²⁵ ] is expressed for radiative heat flux and leads to the form

qr=-4σ3κ*∂T4∂y,

where σ is the Stefan-Boltzmann constant and κ ^* is the mean absorption coefficient. The temperature difference within the flow is sufficiently small such that T ⁴ may be expressed as a linear function of the temperature, then Taylor's series for T ⁴ is about T _∞ after neglecting higher order terms:

T4=4T∞3-3T∞4.

Introducing the following non-dimensional variables, X=xU0υ,Y=yU0υ,U=uU0,V=vU0, τ=tU02υ,T―=T-T∞Tw-T∞,C―=C-C∞Cw-C∞.

Then, Equations ¹ to ⁵ become

∂U∂X+∂V∂Y=0

∂U∂τ+U∂U∂X+V∂U∂Y=∂2U∂Y2+GrT―-MU

∂T―∂τ+U∂T―∂X+V∂T―∂Y=1+RPr·∂2T―∂Y2+Ec∂U∂Y2+Nb∂T―∂Y·∂C―∂Y+Nt∂T―∂Y2

∂C―∂τ+U∂C―∂X+V∂C―∂Y=1Le∂2C―∂Y2+NtNb∂2T―∂Y2.

The non-dimensional boundary conditions are

τ≤0,U=0,V=0,T―=0,C―=0everywhereτ>0,U=0,V=0,T―=0,C―=0atX=0

U=1,V=0,T―=1,C―=1atY=0U=0,V=0,T―=0,C―=0asY→∞,

where the magnetic parameter M=σB02νρU02 , Grashof number Gr=gβTw-T∞νUo3 , radiation parameter R=16σT∞33kκ* , Prandtl number Pr=υα , Eckert number Ec=U02cpTw-T∞ , Lewis number Le=υDB , Brownian parameter Nb=τDBCw-C∞υ , and thermophoresis parameter Nt=DTT∞τνTw-T∞ .

Methods: numerical computation

In order to solve the non-similar unsteady coupled non-linear partial differential equations (Equations 8, 9, 10, and 11), the explicit finite difference method has been developed. For this, a rectangular region of the flow field is chosen, and the region is divided into a grid of lines parallel to X and Y axes, where the X -axis is taken along the plate and the Y -axis is normal to the plate.

Here, the plate of height X _max (=100) is considered, i.e., X varies from 0 to 100 and assumed Y _max (=25) as corresponding to Y  →  ∞ , i.e., Y varies from 0 to 25. There are m (=125) and n (=125) grid spacing in the X and Y directions, respectively, as shown in Figure ² . It is assumed that ΔX and ΔY are constant mesh sizes along the X and Y directions, respectively, and taken as follows: Δ X  = 0.8(0 ≤  X  ≤ 100) and Δ Y  = 0.2(0 ≤  Y  ≤ 25) with the smaller timestep Δ τ  = 0.005.

Figure 2

Let U ^′ , V ^′ , T―′ , and C―′ denote the values of U , V , T― , and T― at the end of a time step, respectively. Using the explicit finite difference approximation, the following appropriate sets of finite difference equations are obtained:

U′i,j-U′i-1,jΔX+Vi,j-Vi,j-1ΔY=0

U′i,j-Ui,jΔτ+Ui,jUi,j-Ui-1,jΔX+Vi,jUi,j+1-Ui,jΔY=Ui,j+1-2Ui,j+Ui,j-1ΔY2+GrT′-MUi,j

T―′i,j-T―i,jΔτ+Ui,jT―i,j-T―i-1,jΔX+Vi,jT―i,j+1-T―i,jΔY=1+RPrT―i,j+1-2T―i,j+T―i,j-1ΔY2+EcUi,j+1-Ui,jΔY2+NbT―i,j+1-T―i,jΔY·C―i,j+1-C―i,jΔY+NtT―i,j+1-T―i,jΔY2

C′―i,j-C―i,jΔτ+Ui,jC―i,j-C―i-1,jΔX+Vi,jC―i,j+1-C―i,jΔY=1LeC―i,j+1-2C―i,j+C―i,j-1ΔY2+NtNbT―i,j+1-2T―i,j+T―i,j-1ΔY2

with initial and boundary conditions

Ui,j0=0,Vi,j0=0,T―i,j0=0,C―i,j0=0

U0,jn=0,V0,jn=0,T―0,j0=0,C―0,jn=0

Ui,0n=1,Vi,0n=0,T―i,0n=1,C―i,0n=1

Ui,Ln=0,Vi,Ln=0,T―i,Ln=0,C―i,Ln=0,whereL→∞,

where the subscripts i and j designate the grid points with X and Y coordinates, respectively, and the superscript n represents a value of time, τ  =  n  ·  Δτ , where n  = 0, 1, 2, ….

Results and discussion

In order to investigate the problem under consideration, the results of numerical values of non-dimensional velocity, temperature, and species concentration within the boundary layer have been computed for different values of magnetic parameter M , radiation parameter R , Prandtl number P _r , Eckert number E _c , Lewis number L _e , Brownian motion parameter N _b , thermophoresis parameter N _t , and Grashof number G _r , respectively. To obtain the steady-state solutions of the computation, the calculations have been carried out up to non-dimensional time τ  = 5 to 80. The velocity, temperature, and concentration profiles do not show any change after non-dimensional time τ  = 50. Therefore, the solution for τ  ≥ 50 is the steady-state solution. The graphical representation of the problem has been shown in Figures ³ , ⁴ , ⁵ , ⁶ , ⁷ , ⁸ , ⁹ , and ¹⁰ .

Figure 3

Figure 4

Figure 5

Figure 6

Figure 7

Figure 8

Figure 9

Figure 10

In order to assess the accuracy of the numerical results, the present results (non-similar solution) are compared with the result obtained by Khan and Pop [ ¹⁹ ] (similar solution) and the values of magnetic parameter M , radiation parameter R , Eckert number E _c , and Grashof number G _r are considered zero (Table ¹ ). From the comparison, excellent agreement is observed.

In Figures ³ , ⁴ , ⁵ , ⁶ , ⁷ , ⁸ , ⁹ , and ¹⁰ , the dimensionless velocity, temperature, and concentration distributions are plotted against Y for the different non-dimensional time τ  = 5 to 50 and corresponding values of Grashof number G _r , magnetic parameter M , radiation parameter R , Eckert number E _c , Prandtl number P _r , Brownian motion parameter N _b , thermophoresis parameter N _t , and Lewis number L _e , respectively.

In Figures ³ and ⁴ , the velocity distribution is plotted respectively for different values of G _r and M . The non-dimensional time considered here is τ  = 5, 20 and 50 and displays the entire step with a different pattern. Here, it is observed that when the values of G _r increase, then the velocity profiles increase and when the values of M increase, then the velocity profiles decrease.

In Figures ⁵ , ⁶ , ⁷ , and ⁸ , the temperature distribution is plotted respectively for different values of R , E _c , P _r , and N _b . The non-dimensional time considered here is τ  = 5, 20 and 50 and displays the entire step with a different pattern. Here, it is observed that when the values of R increase, then the temperature profiles increases; when the values of E _c increase, then the temperature profiles also increase; when the values of P _r increase, then the temperature profiles decrease; and when the values of N _b increase, then the temperature profiles increase.

In Figures ⁹ and ¹⁰ , the concentration distribution is plotted respectively for different values of N _t and L _e . The non-dimensional time considered here is τ  = 5, 20 and 50 and displays the entire step with a different pattern. Here, it is observed that when the values of N _t increase, then the concentration profiles increase, but when the values of L _e increase, then the concentration profiles decrease.

From the present results and the result obtained by Khan and Pop [ ¹⁹ ], it was observed that the flow field shows the same trend with the variation of magnetic parameter M , radiation parameter R , Prandtl number P _r , Eckert number E _c , Lewis number L _e , Brownian motion parameter N _b , thermophoresis parameter N _t , and Grashof number G _r . However, the important part of this work is its comparison with the previous work, i.e., the present study is the unsteady case of Khan and Pop's [ ¹⁹ ] study when the values of magnetic parameter M , radiation parameter R , Eckert number E _c , and Grashof number G _r are considered zero.

Conclusions

An unsteady free convection boundary-layer flow of a nanofluid due to a stretching sheet is studied with the influence of magnetic field and thermal radiation. The explicit finite difference [ ²⁴ ] technique with stability and convergence analysis has been employed as a solution technique to complete the formulation of the unsteady model. For the unsteady case (time-dependent), the non-dimensional time considered here is τ  = 5, 20 and 50 and displays with the entire step and a different line pattern. The results are presented for the effect of various parameters. The velocity, temperature, and concentration effects on the sheet are studied and shown graphically.

From the present study, the concluding remarks have been taken as follows:

• Larger values of the Grashof number showed a significant effect on momentum boundary layer.

• The effect of the Brownian motion and thermophoresis stabilizes the boundary layer growth.

• The boundary layers are highly influenced by the Prandtl number.

• Using magnetic field, the flow characteristics could be controlled.

• The thermal boundary layer thickness increases as a result of increasing radiation.

• The presence of heavier species (large Lewis number) decreases the concentration in the boundary layer.

• The Eckert number has a significant effect on the boundary layer growth.

Authors’ information

AI and LEA are assistant professors of Mathematics. MSK and IK are researchers.