Figures
Abstract
This paper investigates the effect of thermal radiation on unsteady convection flow and heat transfer over a vertical permeable stretching surface in porous medium, where the effects of temperature dependent viscosity and thermal conductivity are also considered. By using a similarity transformation, the governing time-dependent boundary layer equations for momentum and thermal energy are first transformed into coupled, non-linear ordinary differential equations with variable coefficients. Numerical solutions to these equations subject to appropriate boundary conditions are obtained by the numerical shooting technique with fourth-fifth order Runge-Kutta scheme. Numerical results show that as viscosity variation parameter increases both the absolute value of the surface friction coefficient and the absolute value of the surface temperature gradient increase whereas the temperature decreases slightly. With the increase of viscosity variation parameter, the velocity decreases near the sheet surface but increases far away from the surface of the sheet in the boundary layer. The increase in permeability parameter leads to the decrease in both the temperature and the absolute value of the surface friction coefficient, and the increase in both the velocity and the absolute value of the surface temperature gradient.
Citation: Cai W, Su N, Liu X (2014) Unsteady Convection Flow and Heat Transfer over a Vertical Stretching Surface. PLoS ONE 9(9): e107229. https://doi.org/10.1371/journal.pone.0107229
Editor: Zhonghao Rao, China University of Mining and Technology, China
Received: July 2, 2014; Accepted: August 6, 2014; Published: September 29, 2014
Copyright: © 2014 Cai et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper.
Funding: This work is supported by National Natural Science Foundation of China(11271218). Cai obtained financial support from China Scholarship Council through a joint research program between Tsinghua University and Iowa State University (2013-2014). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Convection and heat transfer in porous medium appear in many disciplines, such as thermal and insulation engineering, geophysics and chemistry. In the past few decades, the study in this area has attracted extensive attention of many researchers. Raptis [1] discussed the influence of the radiation on free convection flow through bounded porous medium. Afify [2] analyzed non-Darcy free convection flow through a non-isothermal impermeable vertical plate embedded in a thermally stratified porous medium. Hayat et al. [3]. considered magnetohydrodynamic stagnation-point flow pass a stretching vertical plate with thermal radiation in a porous medium and analyzed the existence of solution by using homotopy analysis method. Convective heat transfer of incompressible viscous fluid over a porous wedge embedded in porous medium was investigated by Anbuchezhian et al. [4]. Mukhopadhyay and Layek [5] discussed the effect of variable viscosity on the boundary layer flow and heat transfer of a fluid through a porous medium. Moreover, Some researchers considered the influence of radiation on heat and mass transfer of steady MHD flow (see [6]–[11]).
The above-mentioned studies [1]–[11] concerned about steady fluid flow. However, the flow and heat transfer is at unsteady conditions in many practical situations, for example, because of a sudden stretching of the plate or due to temperature change of the plate. When the surface is extended suddenly with a certain speed, the flow in the viscous boundary layer near the plate is slowly developed and evolved into a fully developed steady flow after a period of time. So it is necessary to consider physical quantity related to time in mathematical modeling. Ishak et al. [12] studied unsteady laminar boundary layer flow over stretching permeable surface. Further, Tsai et al. [13] and Abd El-Aziz [14] analyzed the impact of radiation on the unsteady flow and heat transfer over stretching surface. In fact, the combined effects of radiation and magnetic field on unsteady viscous incompressible fluid flow have been considered by some scholars (see [15]–[17]).
In addition, when the unsteady stretching surface is located in porous medium, the impact of different factors on the heat transfer is discussed in some recent works. For instance, in the presence of a magnetic field, the viscous incompressible conductive fluid flow along a semi-infinite vertically porous moving plate is researched by Kim [18]. Israel-Cookey et al. [19] analyzed the influence of radiation and viscous dissipation on unsteady MHD free convection flow. For unsteady MHD flow pass a porous vertical plate immersed in porous medium, Samad and Mansur-Rahman [20] considered the combined impact of radiation and magnetic field on MHD free convection flow. Elbashbeshy et al. [21] studied unsteady convective flow over porous stretching surface in the porous medium in the presence of heat source or sick. On this basis, Pal and Hiremath [22] discussed the impact of the dissipation on the heat transfer of fluid flow.
Similar to steady condition, the process of unsteady fluid flow not only includes heat transfer but also mass transfer. When there exists the heat source or sink, Ibrahim et al. [23] analyzed the combined impact of chemical reaction and radiation on heat and mass transfer of MHD flow. Some researchers considered the influence of magnetic field, the radiation, chemical reaction and their combined effects on heat and mass transfer (see [24]–[27]).
According to physical properties of most of realistic fluids, the viscosity and the thermal conductivity are usually related to the temperature and may vary dramatically with temperature. Thus, it would be more reasonable to take into account the effects of temperature dependent viscosity and thermal conductivity in mathematical modeling so as to accurately predict the flow behaviour. In view of this, Seddeek [28] studied the effects of radiation and variable viscosity on the magnetic fluid flow. Solving controlling partial differential equations by using the finite difference method, he further analyzed heat transfer characteristics of fluid flow in the boundary layer. Mukhopadhyay [29] considered the impact of the radiation on unsteady mixed convection boundary layer flow, the effects of the various physical parameters on the velocity and temperature are investigated by using the numerical analysis. By considering assisting and opposing buoyant flow situations, Vajravelu et al. [30] discussed the influence of variable thermal conductivity and radiation on unsteady convective and heat transfer.
To the best of the author's knowledge, no attempt has been made to analyze the combined effects of both temperature dependent viscosity and thermal conductivity on unsteady convection and heat transfer over a vertical permeable stretching surface in porous medium in the presence of radiation. Hence, the aim of the present investigation is design a suitable physical model to describe unsteady two-dimensional incompressible viscous fluid flow past a vertical permeable sheet in a porous medium and solve such an issue numerically using Runge-Kutta fourth-fifth order method with secant shooting technique, which is important from both theoretical and practical point of view because of its wide application to polymer technology and metallurgy. And finally the impact of various physical parameters (such as viscosity variation parameter, the variable thermal conductivity parameter, unsteady parameter, permeability parameter, the suction or injection parameter, convection parameter, thermal radiation parameter, Prandtl number) on the velocity and temperature profiles is displayed in the form of tables and graphs. It is hoped that the results obtained from the present investigation will provide useful information for application and also serve as an effective complement to the previous studies.
The Construction and Analysis of the Model
Consider convection and heat transfer of the unsteady two-dimensional incompressible viscous fluid flow along a vertical permeable sheet in a porous medium. The origin of the Cartesian coordinates is the slot position, which is shown in Fig.1.
The and
axis are taken in the direction along the sheet and perpendicular to it, respectively. Assuming flow region is restricted in the plane
. It is assumed that the magnetic Reynolds number of the flow is extremely small so that the induced magnetic field is negligible, which is a valid assumption on a laboratory scale(see Muhaimin et al. [31]). Simultaneously, suppose that no external electric field and magnetic field are applied, Edge effect and hall effect are also negligible. All dependent variables will be independent of the y-direction (Gorla and Sidawi [32]). The sheet moves with the velocity
in its own plane, where
and
are constants with dimension reciprocal time and satisfying that
and
. The flow is generated by the movement of the sheet; surface temperature of the sheet is
where
is the temperature of the fluid outside the boundary layer;
is constant;
is the kinematic viscosity of the ambient fluid, here
is the dynamic viscosity of ambient fluid and
is the density of the fluid.
is the velocity of suction if
and injection if
. These particular forms of both
and
are chosen in order to obtain self-similar solutions of present problem.
is the gravitational acceleration. The ambient fluid is stationary, that is
. In addition, the temperature-dependent dynamic viscosity
is assumed to linearly vary with the temperature in the form
(1)where
is a small viscosity variation parameter;
;
is the temperature of the fluid inside the thermal boundary layer. Also, it is assumed that the temperature-dependent thermal conductivity
vary as a linear function of temperature given by the following (Vajravelu et al. [30])
(2)where
is thermal conductivity of the ambient fluid;
is a small variable thermal conductivity parameter;
.
Because of the behavior of the boundary layer, the temperature gradient along direction is much larger than that along
direction, therefore this paper only consider the velocity component of the thermal buoyancy which is vertical to the surface of sheet. Under Boussinesq's approximation and previous assumptions, the conservation equations governing this convective flow and heat transfer in boundary layer are
The appropriate boundary conditions to this physical problem are
(7)where
and
are the corresponding velocity components in
and
directions.
is the coefficient of thermal expansion.
is the specific heat at constant pressure.
is the permeability of the porous medium.
is the initial permeability.
is the radiative heat flux. The second term on RHS of the momentum equation (4) denotes the thermal buoyancy effect. Also the first and second terms on the RHS of energy equation (5) represent heat conduction and the effect of thermal radiation, respectively.
Using Rosseland approximation for radiative heat flux term (Raptis [1]), this paper takes
(8)where
and
are the Stefan-Boltzman constant and the mean absorption coefficient, respectively. As described in the literature of Raptis ([1]), suppose that the temperature differences within the flow are sufficiently small, so
may be expressed as a linear function of temperature. This is accomplished by expanding
in a Taylor series about
and neglecting higher-order terms, thus
In view of (9), equation (8) is now reduced to
The similarity transformations for Eqs.(4)–(5) are as follows:
(11)where the stream function
is defined by
and
thus, the continuity equation (3) is satisfied automatically. By calculating, it is easy to obtain
. In the above and behind equations,
denotes the differentiation with respect to
only. Here
and
are the dimensionless velocity and temperature. Substituting (10)–(11) into (4)–(7), a set of ordinary differential equations with variable coefficients is given by
In the above equations, variable physical parameters are defined as
where
is the local Grashof number,
is the local Reynolds number,
is unsteady parameter,
is the permeability parameter,
represents the convection parameter,
is the Prandtl number,
is thermal radiation parameter. All the symbols are defined in the Nomenclature.
The corresponding boundary conditions (6)–(7) for the velocity and temperature fields are transformed to
(14)where
is suction or injection parameter which is used to control the strength and direction of normal flow at the boundary (Vajravelu et al. [30]).
Two physical quantities of interest in this problem are the surface friction coefficient and the local Nusselt number
, which are used to describe the characteristics of surface shear stress and the rate of heat transfer at the surface, respectively. The definitions are as follows
where
and
are the shear stress along the surface and heat transfer from the sheet, respectively. Here
is viscosity coefficient of the ambient fluid and
is the thermal conductivity of the fluid. The results of the calculation can be expressed as
Numerical Procedure
The boundary value problem of higher-order ordinary differential equations (12)–(14) will be transformed to the following initial value problem of first-order ordinary differential equations and will be solved numerically by using an efficient Runge-Kutta fourth-fifth order method with secant shooting technique under Matlab. Now new variables are defined by the equations
Thus, the above two coupled higher order differential equations (12)–(13) may be reduced to five equivalent first-order differential equations as follows
The boundary conditions (14) are converted to the following initial value conditions
(17)where
and
are the suitable initial guess values for
and
, respectively. After choosing
and
, the initial value problems (16)–(17) are solved repeatedly by using Runge-Kutta fourth-fifth order numerical method with secant shooting technique until the boundary conditions
and
are satisfied. In the process of numerical solving (Ibrahim et al. [33]), the values of
and
are improved and the position of the edge of the boundary layer
depending on nine unknown parameters had to be regulated to reach the accuracy. The step size of
satisfies the convergence to the fifth decimal precision in almost all cases, which is sufficient for convergence of numerical solution.
In the next section, by using the above method, the detailed results of numerical simulation will be given to show the impact of various physical parameters on the absolute value of surface friction coefficient, the absolute value of the surface temperature gradient, the velocity and temperature profiles, respectively.
Results and Discussion
In order to validate the method of this paper, ignoring the effects of and
, the model in this paper is the same as that in [30]. Vajravelu et al. [30] used a second order finite difference scheme known as the Keller-Box method to solve nonlinear ordinary differential equations subject to appropriate boundary conditions. Unlike their approach, this paper transforms the boundary value problem of higher-order ordinary differential equations to initial value problem of first-order ordinary differential equations and then uses Runge-Kutta fourth-fifth order numerical method with secant shooting technique to solve these equations numerically. The numerical simulation is employed to investigate the influence of critical parameters on both the absolute values of the surface friction coefficient
and the absolute values of the surface temperature gradient
, and the existing results are in excellent agreement. These favorable comparisons give confidence in the numerical method employed and the numerical results to be presented subsequently. Moreover, this paper analyzes the impact of
and
on
and
. The obtained results are shown in Table 1.
More detailed results are shown in Tables 2–9, which reflects the effect of each parameter on both and
. And also a representative set of numerical results is shown graphically in Figs.2–9 in order to illustrate the effects of various physical parameters on the velocity and temperature profiles.
Fig.2 depicts the effects of various unsteady parameter on velocity field and temperature field. As can be seen from these profiles, as the distance
increases from 0, the velocity and the temperature monotonically decreasing tends to 0. As
increases, the velocity
decreases near the sheet surface, but increases far away from the surface of the sheet in the boundary layer. The temperature and the temperature boundary layer thickness decrease with the increase of
, and simultaneously the temperature gradient increases. So the heat transfer rate of surface increases.
Fig.3 exhibits the velocity and the temperature profiles for different values of the convection parameter . When
is very small, the velocity of the fluid decreases monotonically to 0. However, when the parameter
increases to a certain value, the velocity first monotonically increases to a peak and then monotonically decreasingly approaches to 0. The temperature decreases monotonically to 0 as the distance
increases from 0. The increase in the values of
has the tendency to decrease the temperature and thermal boundary layer thickness but increase the temperature gradient.
Fig.4 represents the variation of both the velocity and temperature profiles in response to a change in the thermal conductivity parameter . The graphs depict that both the velocity and the temperature decreasingly tends to 0 with the increase of values
. The velocity profiles increase slightly with an increase in
. The increase of thermal conductivity parameter values leads to the increase of the temperature and the thermal boundary layer thickness, but simultaneously the decrease of the temperature gradient. In general, it can be seen from Fig.4 that the impact of the thermal conductivity parameter on the temperature field is more noticeable than that on the velocity field. Meanwhile, this phenomenon can be roughly observed from Eq.(13). Actually,
is the diffusion coefficient of the temperature and has a direct impact on the temperature. However, the impact of
on the velocity is achieved through the coupling of the various terms, hence the effect may be weakened.
Figs.5–7 show the velocity and the temperature profiles with respect to the thermal radiation parameter , Prandtl number
and the suction (or injection) parameter
. As shown in these graphs, both the velocity and the temperature decrease monotonically to 0 with the increase of
. As
increases, both the velocity and the temperature increase, and both the velocity boundary layer and the temperature boundary layer become thicker, but the velocity gradient and the temperature gradient decrease. The varies of both
and
also affect the velocity, the velocity gradient, the velocity boundary layer thickness, the temperature, the temperature gradient and temperature boundary layer thickness, which is opposite to the effect of
on the corresponding physical quantities.
Fig.8 reveals the impact of viscosity variation parameter on the velocity and the temperature profiles. As can be seen from these profiles, as the distance
increases from 0, the velocity and the temperature monotonically decreasingly tend to 0. As
increases, the velocity
decreases near the sheet surface, but increases far away from the surface of the sheet in the boundary layer. The temperature profiles decrease slightly with an increase in
and simultaneously the temperature boundary layer thickness become thinner. From Fig.8, it is noticed that the effect of
on the velocity is more pronounced than that on the temperature, which is different from thermal conductivity parameter. Similarly, this phenomenon can be roughly observed from Eq.(12). Furthermore, by comparing Fig.2 and Fig.8, it draws a conclusion that the effect of
on the velocity and the temperature profiles is similar to that of
whereas is not as significant as that of
.
Fig.9 demonstrates that the influence of the permeability parameter on the velocity and the temperature profiles. The results indicate that both the velocity and the temperature decreasingly tend to 0. As
increases, the velocity increases and the velocity boundary layer thickens but the velocity gradient decreases. The impact of
on the temperature, the temperature gradient and the thermal boundary layer thickness is opposite to that of
on the corresponding velocity.
Conclusions
This paper investigates convection flow and heat transfer of an incompressible viscous fluid along a vertical permeable sheet through a porous medium. The resulting ordinary differential equations are solved numerically by using fourth-fifth order Runge-Kutta scheme with secant shooting method. The numerical results are presented for the major parameters including unsteady parameter , convection parameter
, variable thermal conductivity parameter
, thermal radiation parameter
, Prandtl number
, the suction or injection parameter
, viscosity variation parameter
and permeability parameter
. It finally analyzes the impact of the various parameters on flow and heat transfer characteristics.
Compared with the previous literature, this paper not only considers temperature-dependent viscosity and thermal conductivity, but also studies the effects of both the permeability of porous medium and radiation in mathematical modeling, which is helpful to accurately predict the flow behavior. Some novel numerical results are obtained as follows:
- A different numerical method is used to investigate the impact of the physical parameters (
,
,
,
,
and
) on the velocity and the temperature. The obtained results are similar to those in [30] where viscosity variation parameter and permeability parameter are not considered, which further verifies the accuracy of the obtained results.
- Numerical results show that the velocity and temperature monotonously decrease to 0 with the increase of
from the boundary, which means the velocity of the fluid far away from the boundary is close to the velocity of the ambient fluid which is assumed to be zero in this paper. This is quite consistent with the flow of actual fluid.
- With the increase of viscosity variation parameter
, both the absolute value of the surface friction coefficient
and the absolute value of the surface temperature gradient
increase, where the increase of
is very slow. Moreover,
and
reach the minimum in the absence of
. As
increases, the velocity decreases near the sheet surface, whereas increases far away from the surface of the sheet. During this process, the velocity gradient decreases but the velocity boundary layer thickness increases. The temperature profiles decrease slightly with an increase in
and simultaneously the temperature boundary layer becomes thinner. The influence of
on the velocity is more pronounced than that on the temperature.
- An increment in the permeability parameter
leads to the decrease of
whereas the increase of
. The increase in the values of
results in the increase of both the velocity and the velocity boundary layer thickness, which further leads to the decrease of the velocity gradient. But the effect of
on the temperature, the temperature gradient and thermal boundary layer thickness is opposite to that of
on the corresponding velocity.
Acknowledgments
The authors would like to thank Prof. Hailiang Liu and Ph.D Xijun He for their stimulating discussions. The authors are also very grateful to the reviewers and academic editor (Prof. Zhonghao Rao) for their careful reading, constructive comments and valuable suggestions, which leads to a definite improvement of the original manuscript.
Author Contributions
Conceived and designed the experiments: WLC NS XDL. Performed the experiments: WLC. Analyzed the data: WLC NS XDL. Contributed to the writing of the manuscript: WLC NS XDL.
References
- 1. Raptis A (1998) Radiation and free convection flow through a porous medium. International Communications in Heat and Mass Transfer 25: 289–295.
- 2. Afify AA (2007) Effects of variable viscosity on non-Darcy MHD free convection along a non-isothermal vertical surface in a thermally stratified porous medium. Applied Mathematical Modelling 31: 1621–1634.
- 3. Hayat T, Abbas Z, Pop I, Asghar S (2010) Effects of radiation and magnetic field on the mixed convection stagnation-point flow over a vertical stretching sheet in a porous medium. International Journal of Heat and Mass Transfer 53: 466–474.
- 4. Anbuchezhian N, Srinivasan K, Kandasamy R (2012) Magnetohydrodynamic laminar boundary-layer flow over a porous wedge with thermophoresis particle deposition in the presence of variable stream conditions. Meccanica 47: 515–528.
- 5. Mukhopadhyay S, Layek GC (2012) Effects of variable fluid viscosity on flow past a heated stretching sheet embedded in a porous medium in presence of heat source/sink. Meccanica 47: 863–876.
- 6. Seddeek MA (2000) The effect of variable viscosity on hydromagnetic flow and heat transfer past a continuously moving porous boundary with radiation. International Communications in Heat and Mass Transfer 27: 1037–1046.
- 7. Pop SR, Grosan T, Pop I (2004) Radiation effects on the flow near the stagnation point of a stretching sheet. Technische Mechanik 25: 100–106.
- 8. Abbas Z, Hayat T (2008) Radiation effects on MHD flow in a porous space. International Journal of Heat and Mass Transfer 51: 1024–1033.
- 9. Srinivas S, Muthuraj R (2010) Effects of thermal radiation and space porosity on MHD mixed convection flow in a vertical channel using homotopy analysis method. Communications in Nonlinear Science and Numerical Simulation 15: 2098–2108.
- 10. Pal D, Talukdar B (2010) Buoyancy and chemical reaction effects on MHD mixed convection heat and mass transfer in a porous medium with thermal radiation and Ohmic heating. Communications in Nonlinear Science and Numerical Simulation 15: 2878–2893.
- 11. Pal D, Chatterjee S (2010) Heat and mass transfer in MHD non-Darcian flow of a micropolar fluid over a stretching sheet embedded in a porous media with non-uniform heat source and thermal radiation. Communications in Nonlinear Science and Numerical Simulation 15: 1843–1857.
- 12. Ishak A, Nazar R, Pop I (2009) Heat transfer over an unsteady stretching permeable surface with prescribed wall temperature. Nonlinear Analysis: Real World Applications 10: 2909–2913.
- 13. Tsai R, Huang KH, Huang JS (2008) Flow and heat transfer over an unsteady stretching surface with non-uniform heat source. International Communications in Heat and Mass Transfer 35: 1340–1343.
- 14. Abd El-Aziz M (2009) Radiation effect on the flow and heat transfer over an unsteady stretching sheet. International Communications in Heat and Mass Transfer 36: 521–524.
- 15. Abd El-Naby MA, Elbarbary EM, AbdElazem NY (2004) Finite difference solution of radiation effects on MHD unsteady free-convection flow over vertical porous plate. Applied Mathematics and Computation 151: 327–346.
- 16. Jordán JZ (2007) Network simulation method applied to radiation and viscous dissipation effects on MHD unsteady free convection over vertical porous plate. Applied Mathematical Modelling 31: 2019–2033.
- 17. Singh P, Jangid A, Tomer NS, Sinha D (2010) Effects of thermal radiation and magnetic field on unsteady stretching permeable sheet in presence of free stream velocity. International Journal of Information and Mathematical Sciences 6: 63–69.
- 18. Kim YJ (2000) Unsteady MHD convective heat transfer past a semi-infinite vertical porous moving plate with variable suction. International Journal of Engineering Science 38: 833–845.
- 19. Israel-Cookey C, Ogulu A, Omubo-Pepple VB (2003) Influence of viscous dissipation and radiation on unsteady MHD free-convection flow past an infinite heated vertical plate in a porous medium with time-dependent suction. International Journal of Heat and Mass Transfer 46: 2305–2311.
- 20. Samad MA, Mansur-Rahman M (2006) Thermal radiation interaction with unsteady MHD flow past a vertical porous plate immersed in a porous medium. Journal of Naval Architecture and Marine Engineering 3: 7–14.
- 21. Elbashbeshy E, Yassmin DM, Dalia AA (2010) Heat transfer over an unsteady porous stretching surface embedded in a porous medium with variable heat flux in the presence of heat source or sink. African Journal of Mathematics and Computer Science Research 3: 68–73.
- 22. Pal D, Hiremath PS (2010) Computational modeling of heat transfer over an unsteady stretching surface embedded in a porous medium. Meccanica 45: 415–424.
- 23. Ibrahim FS, Elaiw AM, Bakr AA (2008) Effect of the chemical reaction and radiation absorption on the unsteady MHD free convection flow past a semi infinite vertical permeable moving plate with heat source and suction. Communications in Nonlinear Science and Numerical Simulation 13: 1056–1066.
- 24.
Shateyi S, Motsa SS (2010) Thermal radiation effects on heat and mass transfer over an unsteady stretching surface. Mathematical Problems in Engineering 2009.
- 25. Hayat T, Qasim M, Abbas Z (2010) Homotopy solution for the unsteady three-dimensional MHD flow and mass transfer in a porous space. Communications in Nonlinear Science and Numerical Simulation 15: 2375–2387.
- 26. Hayat T, Mustafa M, Asghar S (2010) Unsteady flow with heat and mass transfer of a third grade fluid over a stretching surface in the presence of chemical reaction. Nonlinear Analysis: Real World Applications 11: 3186–3197.
- 27. Kesavaiah DC, Satyanarayana PV, Venkataramana S (2011) Effects of the chemical reaction and radiation absorption on an unsteady MHD convective heat and mass transfer flow past a semi-infinite vertical permeable moving plate embedded in a porous medium with heat source and suction. International Journal of Applied Mathematics and Mechanics 7: 52–69.
- 28. Seddeek MA (2002) Effects of radiation and variable viscosity on a MHD free convection flow past a semi-infinite flat plate with an aligned magnetic field in the case of unsteady flow. International Journal of Heat and Mass Transfer 45: 931–935.
- 29. Mukhopadhyay S (2009) Effect of thermal radiation on unsteady mixed convection flow and heat transfer over a porous stretching surface in porous medium. International Journal of Heat and Mass Transfer 52: 3261–3265.
- 30. Vajravelu K, Prasad KV, Ng C (2013) Unsteady convective boundary layer flow of a viscous fluid at a vertical surface with variable fluid properties. Nonlinear Analysis: Real World Applications 14: 455–464.
- 31. Muhaimin I, Kandasamy R, Hashim I, Khamis AB (2009) On the effect of chemical reaction, heat and mass transfer on nonlinear MHD boundary layer past a porous shrinking sheet with suction. Theoretical and Applied Mechanics 36: 101–117.
- 32. Gorla RSR, Sidawi I (1994) Free convection on a vertical stretching surface with suction and blowing. Applied Scientific Research 52: 247–257.
- 33. Ibrahim W, Shankar B, Nandeppanavar MM (2013) MHD stagnation point flow and heat transfer due to nanofluid towards a stretching sheet. International Journal of Heat and Mass Transfer 56: 1–9.