The Scientific World Journal

Volume 2014 (2014), Article ID 362351, 11 pages

http://dx.doi.org/10.1155/2014/362351

Research Article

Finite Element Solution of Unsteady Mixed Convection Flow of Micropolar Fluid over a Porous Shrinking Sheet

Diksha Gupta, Lokendra Kumar, and Bani Singh

Department of Mathematics, Jaypee Institute of Information Technology, A-10, Sector 62, Noida, Uttar Pradesh 201307, India

Received 14 August 2013; Accepted 24 October 2013; Published 4 February 2014

Academic Editors: A. Al-Sarkhi, C. Bao, and A. Kosar

Copyright © 2014 Diksha Gupta et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The objective of this investigation is to analyze the effect of unsteadiness on the mixed convection boundary layer flow of micropolar fluid over a permeable shrinking sheet in the presence of viscous dissipation. At the sheet a variable distribution of suction is assumed. The unsteadiness in the flow and temperature fields is caused by the time dependence of the shrinking velocity and surface temperature. With the aid of similarity transformations, the governing partial differential equations are transformed into a set of nonlinear ordinary differential equations, which are solved numerically, using variational finite element method. The influence of important physical parameters, namely, suction parameter, unsteadiness parameter, buoyancy parameter and Eckert number on the velocity, microrotation, and temperature functions is investigated and analyzed with the help of their graphical representations. Additionally skin friction and the rate of heat transfer have also been computed. Under special conditions, an exact solution for the flow velocity is compared with the numerical results obtained by finite element method. An excellent agreement is observed for the two sets of solutions. Furthermore, to verify the convergence of numerical results, calculations are conducted with increasing number of elements.

1. Introduction

In the last few decades, the interest for non-Newtonian fluids has considerably increased due to their connection with applied sciences. The motion of these fluids plays essential role not only in theory but also in many industrial processes. Among the various non-Newtonian fluid models, the micropolar fluids have acquired the special attention in recent years due to their applications in polymeric fabrication, materials processing, and biotechnology. Flow and heat transfer behaviour of these fluids cannot be described by the classical theory of continuum mechanics. Eringen [1] has formulated the theory of micropolar fluids which describes the physics of such fluids. In the micropolar fluid theory, two new variables to the velocity are added which were not presented in the Navier-Stokes model. These variables are microrotations that represent spin and microinertia tensors which describe the distribution of atoms and molecules inside the microscopic fluid particles. This class of fluids represents, mathematically, many industrial important fluids such as paints, lubricants, polymers, human and animal blood, colloidal suspensions, and liquid crystals. The theory and applications of micropolar fluids can be found in the books by Eringen [2] and Bég et al. [3]. Later Eringen [4] extended the theory of micropolar fluids to thermo-microfluids. This theory takes into account thermal effects, that is, heat conduction, convection, and dissipation. These effects were not included in the classical field theories.

The boundary layer flow induced by stretching surface is of great practical interest because it occurs in a number of engineering processes. Boundary layer flow of a Newtonian fluid caused by a linearly stretching sheet was first examined by Crane [5]. He gave a similarity solution in closed analytical form for the steady two-dimensional problem. The flow and heat transfer of micropolar fluid past a continuously moving porous plate was analyzed by Takhar and Soundalgekar [6]. Hassanien and Gorla [7] conducted numerical computations to examine the heat transfer characteristics of micropolar fluid past a nonisothermal sheet with suction and blowing. Thereafter, various aspects of micropolar fluid flows from a stretching surface have been reported by El-Arabawy [8], Nazar et al. [9], Kumar [10], Ishak [11], and Rawat et al. [12].

The above-mentioned studies deal with a steady flow only. However, in certain cases, the flow and heat transfer can be unsteady due to a sudden stretching of the flat sheet or by a change of the temperature or heat flux of the sheet. Devi et al. [13] investigated the unsteady three-dimensional flow caused by a stretching flat surface. Using the finite difference scheme in combination with the quasi-linearization technique Rajeswari, and Nath [14] reported on the unsteady flow over a stretching surface in a rotating fluid. Exact similarity solution of the heat transfer in a liquid film on an unsteady stretching surface was obtained by Andersson et al. [15]. Abd El-Aziz [16] studied numerically the effect of radiation on the heat and fluid flow over an unsteady stretching sheet showing that the heat transfer rate is an increasing function of radiation and unsteadiness parameters. Lok et al. [17] obtained the numerical solution for the unsteady boundary layer flow of a micropolar fluid near the stagnation point of a plane surface. By taking into account strong and weak concentration of microelements, Hayat et al. [18] examined the effect of magnetic field on the time-dependent flow of micropolar fluid between stretching sheets. Bachok et al. [19] studied theoretically the unsteady boundary layer flow and heat transfer due to a stretching sheet, showing that the surface shear stress and the heat transfer rate are an increasing function of unsteadiness parameter.

Recently, the boundary layer flow of incompressible fluid over a shrinking sheet has attracted extensive attention due to its increasing applications in polymeric materials processing. This type of flow was first examined analytically by Wang [20]. Later, Miklavčič and Wang [21] proved the existence and uniqueness of the solution for the flow over a shrinking sheet. From the physical point of view, steady flow over a shrinking sheet is not possible since the generated vorticity is not confined within the boundary layer. To overcome this difficulty, the flow requires a certain amount of external opposite force at the sheet. This has been extensively discussed in the literature. Some important references are Miklavčič and Wang [21], Wang [22], and Bachok et al. [23]. Closed form exact solutions of MHD Newtonian flow over a shrinking sheet were derived by Fang and Zhang [24]. Using shooting method, Bhattacharyya and Layek [25] obtained the numerical solutions for the flow and heat transfer over a porous shrinking sheet in the presence of thermal radiation. The shrinking sheet problem was extended to micropolar fluids by Ishak et al. [26] who considered stagnation point flow. Yacob and Ishak [27] analyzed the flow and heat transfer over a shrinking sheet immersed in a micropolar fluid. Recently, numerical simulation of mixed convection flow of micropolar fluid over a shrinking sheet with thermal radiation was conducted by Gupta et al. [28].

Many researchers have investigated the unsteady flow over a shrinking sheet. Fang et al. [29] examined the unsteady viscous flow over a shrinking surface with mass suction and highlighted the deviation in flow behaviour for an unsteady shrinking sheet compared with an unsteady stretching sheet. Effect of radiation on the unsteady flow and heat transfer induced by a permeable shrinking sheet was reported by Ali et al. [30]. Bhattacharyya [31] examined the effects of radiation and heat source/sink on unsteady flow and heat transfer past a shrinking sheet with suction/injection. Numerical solution for the unsteady stagnation point flow and heat transfer over a stretching/shrinking sheet with prescribed heat flux was analyzed by Suali et al. [32]. They found that the skin friction and local Nusselt number increase with unsteadiness. Mahapatra and Nandy [33] studied the unsteady stagnation-point flow and heat transfer over a linearly shrinking sheet in the presence of velocity and thermal slips.

In many of the studies described above, buoyancy effect plays a significant role. The buoyancy force, developed from the temperature difference, induces a longitudinal pressure gradient which in turn modifies the flow field and the rate of heat transfer from the surface. A computational study of unsteady mixed convection flow in stagnation region adjacent to a vertical surface was conducted by Devi et al. [34]. Ishak et al. [35] presented the numerical solution for the unsteady mixed convection flow and heat transfer over a stretching vertical sheet. Both assisting and opposing flow cases were taken into consideration. Unsteady mixed convection boundary layer flow over a stretching vertical surface in the presence of velocity slip was presented by Mukopadhyay [36]. Sharma et al. [37] employed element-free Galerkin method to study the mixed convection flow and heat transfer of a viscous fluid over an unsteady stretching sheet in a porous medium. Using Keller-Box method, Vajravelu et al. [38] examined the effects of variable thermal conductivity, thermal radiation and the thermal buoyancy on the unsteady fluid flow and heat transfer at a porous stretching sheet.

The present paper investigates the unsteady mixed convection flow and heat transfer of an incompressible micropolar fluid over a vertical shrinking sheet with time-dependent suction at the sheet. Viscous dissipation effects are also included in the energy equation. The velocity and temperature of the sheet are assumed to vary with the horizontal coordinate and time . Similarity transformations are employed for the conversion of governing time-dependent boundary layer equations into ordinary differential equations. These equations are then solved numerically using finite element method. The influence of suction parameter, unsteadiness parameter, buoyancy parameter, and Eckert number has been depicted graphically. The skin friction and the rate of heat transfer have also been computed and tabulated for these parameters. The current study has applications in industrial polymeric materials processing and has not been considered so far to the knowledge of the authors.

2. Mathematical Model

Consider the unsteady, mixed convection, and boundary layer flow of an incompressible micropolar fluid past a permeable shrinking sheet. A schematic representation of the physical model and coordinate system is depicted in Figure 1. It is assumed that, for time , the fluid and heat flows are steady. The unsteady fluid and heat flow start at . The velocity of shrinking sheet is and the temperature of the sheet is , where , are constants (with , , wher ) and both have dimension , while is a constant and has dimensions (temperature/length). The coordinate system is such that -axis is taken along the shrinking sheet in a direction opposite to sheet motion and -axis is normal to it. Time-dependent suction is considered normal to the shrinking sheet. The physical properties of the fluid are assumed to be constant except density variation due to temperature difference which is used only to express the body force term as the buoyancy term. The effect of viscous dissipation is also included in the energy equation. Under these assumptions the governing boundary layer equations for unsteady flow over a shrinking sheet may be presented as follows.

 Continuity equation

 Momentum equation

 Angular momentum equation

 Energy equation

Figure 1: Physical model and coordinate system.

The associated boundary conditions are

It is assumed that is a variable distribution of suction through porous sheet and is given by (Bhattacharyya et al. [39]).

Introducing the similarity variable and the dimensionless functions , , and as follows:

where is a stream function defined as , which identically satisfies the continuity equation (1). On applying the transformations, (2)–(4) are reduced to

and the corresponding boundary conditions (5) now transform to

where prime denotes the differentiation with respect to only, (coupling constant parameter), (buoyancy parameter), (local Grashof number), (unsteadiness parameter), (Prandtl number), (Eckert number), and (suction parameter); corresponds to suction and corresponds to injection.

The quantities of physical interest, namely, the local skin friction coefficient and the rate of heat transfer, are, respectively, prescribed by

where the local wall shear stress and the heat transfer from the sheet are given by

Using the similarity transformations given in (6), we obtain the following:

where is the local Reynolds number.

3. Method of Solution

The set of differential equations given in (7) are non-linear and therefore, cannot be solved analytically. Thus for the solution of this problem, finite element method has been implemented. Comprehensive details of this method can be found in Reddy [40]. In order to apply finite element method first we assume

Using (12), equation (7) reduces to

and the corresponding boundary conditions now become

It has been observed that, for , there is no appreciable effect on the results. Therefore, for the computational purposes, can be fixed at 8.

3.1. Variational Formulation

The variational form associated with (12) and (13) over a typical two-noded linear element is given by

where , , , and are weight functions which may be viewed as the variation in , , , and , respectively.

3.2. Finite Element Formulation

The finite element model can be obtained from (15) by substituting finite element approximations of the form

with , where are the shape functions for a typical element and are taken as follows:

The finite element model of the equations thus formed is given by

where and () are the matrices of order , and respectively and are defined as follows:

where and are assumed to be known. After the assembly of element equations, a system of non-linear equations is obtained; therefore, an iterative scheme must be utilized to solve it. The system is linearized by incorporating the functions and , which are assumed to be known at lower iteration level and the computations for, , , and are then carried out for higher levels. This process is repeated until the desired accuracy of 0.00005 is attained. Convergence of the results with increasing number of elements is shown in Table 1. It is clear from the table that, for more than elements no significant variation in the values of , , and is observed. Thus the final results are reported for elements. This confirms the mesh-independence of the present computation.

Table 1: Convergence of results with the variation of number of elements ().

20
2.7194
0.0714
−0.1087
0.1857

40
2.6986
0.0524
−0.1198
0.1772

60
2.6955
0.0496
−0.1223
0.1756

80
2.6945
0.0487
−0.1231
0.1751

100
2.6941
0.0483
−0.1236
0.1748

120
2.6939
0.0481
−0.1238
0.1747

140
2.6937
0.0480
−0.1239
0.1746

160
2.6937
0.0479
−0.1240
0.1745

180
2.6936
0.0479
−0.1241
0.1745

For steady state , viscous fluid and in the absence of buoyancy force , the exact solution for , as obtained by Fang and Zhang [24] with , is given by , where .

The comparison of the flow velocity obtained by finite element method and the exact solution given by Fang and Zhang [24] is shown in Table 2. It is clear from the table that numerical results obtained are in complete agreement with the exact solution and thus confirm the validity and accuracy of the FEM.

Table 2: Comparison of the flow velocity obtained by analytical method [24] and FEM () in the special case.

[24]
FEM

0
−1
−1

1
−0.07295
−0.07290

2
−0.00532
−0.00531

3
−0.00039
−0.00039

4
−0.00003
−0.00003

5, 6, 7, 8
0.00000
0.00000

4. Results and Discussion

To study the behaviour of velocity, microrotation, and temperature functions, comprehensive numerical computations are carried out for various values of the parameters namely suction parameter , unsteadiness parameter, buoyancy parameter and Eckert number Ec. The other parameters such as coupling constant parameter and Prandtl number Pr are kept fixed at 2.0 and 0.733 respectively. The results obtained are presented through the graphs as shown in Figures 2–13. The skin friction and the rate of heat transfer have also been computed for these parameters and are tabulated in Tables 3 and 4.

Table 3: The skin friction coefficient for different values of , and .

3
2.49867
1
2.49867
0
1.31709
0.02
2.15971

3.5
2.58229
2
2.49630
1
1.74698
0.25
2.26044

4
2.70175
3
2.52259
2
2.12938
0.5
2.37579

4.5
2.85110
4
2.57288
3
2.49867
0.75
2.49867

5
3.02241
5
2.63686
4
2.87143
1
2.63105

Table 4: The local Nusselt number for different values of , and .

3
−0.24924
1
−0.24924
0
0.88895
0.02
1.99596

3.5
0.12914
2
0.14660
1
0.54444
0.25
1.39436

4
0.45990
3
0.46597
2
0.17095
0.5
0.63527

4.5
0.75067
4
0.71556
3
−0.24924
0.75
−0.24924

5
1.01279
5
0.91656
4
−0.73308
1
−1.28127

Figure 2: Velocity distribution for different .

Figure 3: Microrotation distribution for different .

Figure 4: Temperature distribution for different .

Figure 5: Velocity distribution for different .

Figure 6: Microrotation distribution for different .

Figure 7: Temperature distribution for different .

Figure 8: Velocity distribution for different .

Figure 9: Microrotation distribution for different .

Figure 10: Temperature distribution for different .

Figure 11: Velocity distribution for different .

Figure 12: Microrotation distribution for different .

Figure 13: Temperature distribution for different .

Figures 2–4 display the effect of suction parameter on the boundary layer flow induced by an unsteady shrinking sheet. Suction is the most suitable force to sustain the flow near the shrinking sheet by confining the generated vorticity inside the boundary layer. It is clear from Figure 2 that the effect of suction parameter on the velocity is negligible near the sheet whereas away from the sheet velocity decreases with increase in suction parameter. Negative values of velocity close to the sheet indicate the region of reverse flow. Figure 3 shows that as suction parameter increases the microrotation decreases near the sheet. After covering a small distance from the sheet, the opposite behaviour is observed. The negative values of microrotation show the reverse rotation of microparticles. Figure 4 illustrates the influence of suction parameter on the temperature distribution. It is evident from the figure that temperature of the fluid decreases with increase in suction parameter. As suction is applied, thermal boundary layer thickness decreases as a result of which temperature of the fluid in the boundary layer decreases. This observation is in agreement with Bhattacharyya and Layek [25]. Thus, suction parameter can be used effectively for controlling the flow and heat transfer characteristics.

The effect of unsteadiness parameter on the velocity, microrotation, and temperature functions is shown in Figures 5–7. Figure 5 shows that increasing the value of unsteadiness parameter tends to decrease the velocity in the boundary layer. Near the sheet flow velocity is negative whereas away from the sheet it becomes positive and finally satisfies the free stream boundary condition. Thus close to the sheet the flow is strongly reversed. Figure 6 demonstrates that, near the sheet the effect of unsteadiness on the microrotation is negligible, but away from the sheet it increases with increase in unsteadiness parameter. After covering a certain distance from the sheet, all profiles converge and finally satisfy the far field boundary condition. Figure 7 reveals that the temperature decreases with increase in unsteadiness parameter . Temperature at the sheet is invariant for higher values of unsteadiness parameter. These patterns are consistent with the findings of Ishak et al. [35]. Physically as the value of unsteadiness parameter increases the sheet loses more heat as a result of which temperature of the fluid decreases.

Figures 8–10 display the nature of velocity, microrotation and temperature distributions with buoyancy parameter . Physically positive values of are associated with cooling of the sheet (assisting flow), negative values of indicate heating of the sheet (opposing flow), and implies vanishing buoyancy effects. From Figure 8, it can be seen that velocity of the fluid increases with increase in buoyancy parameter. This is due to the reason that an increase in the value of buoyancy parameter leads to an increase in the temperature difference . This leads to an increase in the convection currents as a result of which fluid velocity increases and thus, boundary layer thickness decreases. This trend has also been identified by Shit and Haldar [41]. Figure 9 indicates that the microrotation decreases initially with increase in buoyancy parameter. For all profiles intersect and then increases with increases in buoyancy parameter. In all cases microrotation is negative which shows the reverse rotation of microelements. A crossing over point appears in the temperature profile as shown in Figure 10. This point is special where all temperature curves cross each other; that is, the temperature profile shows different behaviour before and after this point. It is observed that as the buoyancy parameter increases, fluid temperature increases up to this point and decreases after this point. Maximum temperature occurs near the sheet corresponding to . Temperature at the sheet is invariant for lower values of buoyancy parameter.

Figures 11–13 present the distributions of velocity, microrotation and temperature functions with Eckert number Ec. This parameter is called the fluid motion controlling parameter. From Figure 11 it is observed that the velocity increases with increase in Eckert number Ec. Boundary layer thickness will therefore decrease as Eckert number increases. Flow reversal arises near the sheet as testified by the negative values of velocity. Figure 12 demonstrates that in the vicinity of sheet microrotation decreases with increase in Eckert number. Minimum value of microrotation corresponds to and this value is approximately. After a small distance from the sheet all profiles converge and finally satisfy the far field boundary condition. Figure 13 reveals that the temperature as well as thermal boundary layer thickness increases with increase in Eckert number Ec. A temperature overshoot near the sheet has been observed for higher values of Eckert number. This is due to the fact that for higher values of Eckert number, there is a significant generation of heat due to viscous dissipation near the sheet, so that the temperature in the region close to the sheet exceeds the temperature of the wall .

Table 3 gives the skin friction for different values of , , and Ec. It is evident that the skin friction decreases slightly with a small increase in unsteadiness parameter while it increases with increase in suction parameter, buoyancy parameter, Eckert number and large values of unsteadiness parameters. Physically positive values of the skin friction show that the fluid exerts a drag force on the sheet. It has also been observed that the skin friction is lower in the absence of buoyancy parameter. Thus skin friction can be reduced effectively by assigning lower values to unsteadiness parameter and in the absence of buoyancy parameter. From Table 4 it may be noted that the rate of heat transfer decreases initially with increase in buoyancy parameter, Eckert number small values of suction, and unsteadiness parameters. It is also observed that heat transfer rate increases numerically with increase in suction, unsteadiness parameters, large values of buoyancy parameter, and Eckert number. Positive values of the heat transfer rate indicate that the heat is transferred from the surface of the sheet to the fluid and negative values mean the opposite. Thus effective cooling of the sheet can be achieved by the judicious selection of these parameters.

5. Conclusions

The present study has addressed theoretically and numerically the unsteady mixed convection flow and heat transfer of an incompressible micropolar fluid over a porous shrinking sheet in the presence of viscous dissipation. Using similarity transformations the governing time-dependent boundary layer equations are reduced to a set of nonlinear ordinary differential equations. A variational finite element method has been employed to solve the nondimensional momentum, angular momentum, and thermal boundary layer equations, subject to physically realistic boundary conditions. Under limiting cases the numerical results obtained for the flow velocity are compared very well with the exact solution available in the literature. Numerical computations have clearly demonstrated that the drag can be reduced effectively with the judicious selection of unsteadiness parameter, and buoyancy parameter. It has also been found that a fast rate of cooling can be achieved by implementing suction parameter, unsteadiness parameter, higher values of buoyancy parameter and Eckert number. The present study has neglected velocity and thermal slip effects at the sheet. These may be considered in the future. It is hoped that the results obtained from the present work may be useful for different model investigations.

Nomenclature

: | Specific heat at constant pressure
: | Skin friction coefficient
: | Dimensionless velocity
: | Dimensionless microrotation
: | Gravitational acceleration
: | Local Grashof number
: | Microinertia density
: | Coupling constant parameter
: | Microrotation component
: | Local Nusselt number
: | Prandtl number
: | Heat flux
: | Local Reynolds number
: | Constant characteristic of the fluid
: | Temperature of the fluid
: | Temperature of the ambient fluid
: | Time
: | Velocity in the -direction
: | Velocity of the sheet
: | Velocity in the -direction
: | Velocity at the wall
: | Weight functions
: | Distance along the surface
: | Distance normal to the surface.

Greek Symbols

: | Spin gradient viscosity
: | Similarity variables
: | The dynamic viscosity
: | Density of the fluid
: | Thermal conductivity
: | Kinematic viscosity
: | Dimensionless temperature
: | Suction parameter
: | Stream function
: | Buoyancy parameter
: | Unsteadiness parameter
: | Wall shear stress.

Subscripts

: | Surface condition
: | Conditions far away from the surface.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

A. C. Eringen, “Theory of micropolar fluids,” Journal of Mathematical Fluid Mechanics, vol. 16, pp. 1–18, 1966.

A. C. Eringen, Microcontinuum Field Theories II-Fluent Media, Springer, New York, NY, USA, 1st edition, 2001.

O. A. Bég, R. Bhargava, and M. M. Rashidi, Numerical Simulation in Micropolar Fluid Dynamics, Lambert Academic Publishing, 2011.

A. C. Eringen, “Theory of thermomicrofluids,” Journal of Mathematical Analysis and Applications, vol. 38, no. 2, pp. 480–496, 1972.

L. J. Crane, “Flow past a stretching plate,” Zeitschrift für angewandte Mathematik und Physik ZAMP, vol. 21, no. 4, pp. 645–647, 1970.

H. S. Takhar and V. M. Soundalgekar, “Flow and heat transfer of a micropolar fluid past a continuously moving porous plate,” International Journal of Engineering Science, vol. 23, no. 2, pp. 201–205, 1985.

I. A. Hassanien and R. S. R. Gorla, “Heat transfer to a micropolar fluid from a non-isothermal stretching sheet with suction and blowing,” Acta Mechanica, vol. 84, no. 1–4, pp. 191–199, 1990.

H. A. M. El-Arabawy, “Effect of suction/injection on the flow of a micropolar fluid past a continuously moving plate in the presence of radiation,” International Journal of Heat and Mass Transfer, vol. 46, no. 8, pp. 1471–1477, 2003.

R. Nazar, N. Amin, D. Filip, and I. Pop, “Stagnation point flow of a micropolar fluid towards a stretching sheet,” International Journal of Non-Linear Mechanics, vol. 39, no. 7, pp. 1227–1235, 2004.

L. Kumar, “Finite element analysis of combined heat and mass transfer in hydromagnetic micropolar flow along a stretching sheet,” Computational Materials Science, vol. 46, no. 4, pp. 841–848, 2009.

A. Ishak, “Thermal boundary layer flow over a stretching sheet in a micropolar fluid with radiation effect,” Meccanica, vol. 45, no. 3, pp. 367–373, 2010.

S. Rawat, R. Bhargava, S. Kapoor, and O. A. Beg, “Heat and mass transfer of a chemically reacting micropolar fluid over a linear stretching sheet in darcy forchheimer porous medium,” International Journal of Computational and Applied Mathematics, vol. 44, no. 6, pp. 40–51, 2012.

C. D. S. Devi, H. S. Takhar, and G. Nath, “Unsteady, three-dimensional, boundary-layer flow due to a stretching surface,” International Journal of Heat and Mass Transfer, vol. 29, no. 12, pp. 1996–1999, 1986.

V. Rajeswari and G. Nath, “Unsteady flow over a stretching surface in a rotating fluid,” International Journal of Engineering Science, vol. 30, no. 6, pp. 747–756, 1992.

H. I. Andersson, J. B. Aarseth, and B. S. Dandapat, “Heat transfer in a liquid film on an unsteady stretching surface,” International Journal of Heat and Mass Transfer, vol. 43, no. 1, pp. 69–74, 2000.

M. Abd El-Aziz, “Radiation effect on the flow and heat transfer over an unsteady stretching sheet,” International Communications in Heat and Mass Transfer, vol. 36, no. 5, pp. 521–524, 2009.

Y. Y. Lok, N. Amin, and I. Pop, “Unsteady boundary layer flow of a micropolar fluid near the rear stagnation point of a plane surface,” International Journal of Thermal Sciences, vol. 42, no. 11, pp. 995–1001, 2003.

T. Hayat, M. Nawaz, and S. Obaidat, “Axisymmetric magnetohydrodynamic flow of micropolar fluid between unsteady stretching surfaces,” Applied Mathematics and Mechanics, vol. 32, no. 3, pp. 361–374, 2011.

N. Bachok, A. Ishak, and R. Nazar, “Flow and heat transfer over an unsteady stretching sheet in a micropolar fluid,” Meccanica, vol. 46, no. 5, pp. 935–942, 2011.

C. Y. Wang, “Liquid film on an unsteady stretching sheet,” Quarterly of Applied Mathematics, vol. 48, no. 4, pp. 601–610, 1990.

M. Miklavčič and C. Y. Wang, “Viscous flow due to a shrinking sheet,” Quarterly of Applied Mathematics, vol. 64, no. 2, pp. 283–290, 2006.

C. Y. Wang, “Stagnation flow towards a shrinking sheet,” International Journal of Non-Linear Mechanics, vol. 43, no. 5, pp. 377–382, 2008.

N. Bachok, A. Ishak, and I. Pop, “Unsteady three-dimensional boundary layer flow due to a permeable shrinking sheet,” Applied Mathematics and Mechanics, vol. 31, no. 11, pp. 1421–1428, 2010.

T. Fang and J. Zhang, “Closed-form exact solutions of MHD viscous flow over a shrinking sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 7, pp. 2853–2857, 2009.

K. Bhattacharyya and G. C. Layek, “Effects of suction/blowing on steady boundary layer stagnation-point flow and heat transfer towards a shrinking sheet with thermal radiation,” International Journal of Heat and Mass Transfer, vol. 54, no. 1, pp. 302–307, 2011.

A. Ishak, Y. Y. Lok, and I. Pop, “Stagnation-point flow over a shrinking sheet in a micropolar fluid,” Chemical Engineering Communications, vol. 197, no. 11, pp. 1417–1427, 2010.

N. A. Yacob and A. Ishak, “Micropolar fluid flow over a shrinking sheet,” Meccanica, vol. 47, no. 2, pp. 293–299, 2012.

D. Gupta, L. Kumar, O. A. Bég, and B. Singh, “Finite element simulation of mixed convection flow of micropolar fluid over a shrinking sheet with thermal radiation,” Journal of Process Mechanical Engineering, 2013.

T. G. Fang, J. Zhang, and S. S. Yao, “Viscous flow over an unsteady shrinking sheet with mass transfer,” Chinese Physics Letters, vol. 26, no. 1, Article ID 014703, 2009.

F. M. Ali, R. Nazar, N. M. Arifin, and I. Pop, “Unsteady flow and heat transfer past an axisymmetric permeable shrinking sheet with radiation effect,” International Journal for Numerical Methods in Fluids, vol. 67, no. 10, pp. 1310–1320, 2011.

K. Bhattacharyya, “Effects of radiation and heat source/sink on unsteady MHD boundary layer flow and heat transfer over a shrinking sheet with suction/injection,” Frontiers of Chemical Engineering in China, vol. 5, no. 3, pp. 376–384, 2011.

M. Suali, N. M. A. N. Long, and A. Ishak, “Unsteady stagnation point flow and heat transfer over a stretching/shrinking sheet with prescribed surface heat flux,” Applied Mathematics and Computational Intelligence, vol. 1, no. 1, pp. 1–11, 2012.

T. R. Mahapatra and S. K. Nandy, “Slip effects on unsteady stagnation-point flow and heat transfer over a shrinking sheet,” Mecannica, vol. 48, no. 7, pp. 1599–1606, 2013.

C. D. S. Devi, H. S. Takhar, and G. Nath, “Unsteady mixed convection flow in stagnation region adjacent to a vertical surface,” Heat Mass Transfer, vol. 26, no. 2, pp. 71–79, 1991.

A. Ishak, R. Nazar, and I. Pop, “Boundary layer flow and heat transfer over an unsteady stretching vertical surface,” Meccanica, vol. 44, no. 4, pp. 369–375, 2009.

S. Mukhopadhyay, “Effects of slip on unsteady mixed convective flow and heat transfer past a stretching surface,” Chinese Physics Letters, vol. 27, no. 12, Article ID 124401, 2010.

R. Sharma, R. Bhargava, and I. V. Singh, “A numerical solution of MHD convection heat transfer over an unsteady stretching surface embedded in a porous medium using element free Galerkin method,” International Journal of Applied Mathematics and Mechanics, vol. 8, no. 10, pp. 83–103, 2012.

K. Vajravelu, K. V. Prasad, and C. O. Ng, “Unsteady convective boundary layer flow of a viscous fluid at a vertical surface with variable fluid properties,” Nonlinear Analysis: Real World Applications, vol. 14, no. 1, pp. 455–464, 2013.

K. Bhattacharyya, S. Mukhopadhyay, and G. C. Layek, “Unsteady MHD boundary layer flow with diffusion and first order chemical reaction over a permeable stretching sheet with suction or blowing,” Chemical Engineering Communications, vol. 200, no. 3, pp. 379–397, 2013.

J. N. Reddy, Introduction to the Finite Element Method, McGraw-Hill, 2005.

G. C. Shit and R. Haldar, “Effects of thermal radiation on MHD viscous fluid flow and heat transfer over nonlinear shrinking porous sheet,” Applied Mathematics and Mechanics, vol. 32, no. 6, pp. 677–688, 2011.