`Journal of Applied MathematicsVolume 2012, Article ID 675287, 17 pageshttp://dx.doi.org/10.1155/2012/675287`
Research Article

## Lie Group Analysis of Unsteady Flow and Heat Transfer over a Porous Surface for a Viscous Fluid

1Department of Mechanical Engineering, Celal Bayar University, Muradiye, 45140 Manisa, Turkey
2Applied Mathematics and Computation Center, Celal Bayar University, Muradiye, 45140 Manisa, Turkey

Received 27 September 2012; Accepted 27 October 2012

Copyright © 2012 M. B. Akgül 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 problem of a two-dimensional, unsteady flow and a heat transfer of a viscous fluid past a surface in the presence of variable suction/injection is analyzed. The unsteadiness is due to the time dependent free stream flow. The governing equations are derived with the usual boundary layer approximation. Using Lie group theory, a group classification of the equations with respect to the variable free stream flow and suction/injection velocity is performed. Restrictions imposed by the boundary conditions on the symmetries are discussed. Adopting the obtained symmetry groups, governing partial differential equations are converted into ordinary differential equations and then solved numerically. Effects of the dimensionless problem parameters on the velocity and temperature profiles are outlined in the figures.

#### 1. Introduction

Boundary layer flow and heat transfer for an incompressible viscous fluid past a heated porous surface have attracted a great interest because of the importance in engineering applications such as cooling of nuclear reactors and turbine blades, food processing, crystal growth, flow control on airfoils, and electronic cooling. Blasius [1] first presented a similarity solution for velocity distribution within the boundary layer for a viscous incompressible fluid over a flat plate. Khaled and Vafai [2] studied hydromagnetic effects on dynamical and thermal boundary layer characteristics. This analysis is concerned with a family of squeezed flows. They displayed effects of the magnetic parameter, Prandtl number, suction/injection parameter, and squeezing parameter on velocity and temperature profiles graphically. Elbashbeshy and Bazid [3] studied the heat transfer over an unsteady stretching surface. They obtained a new similarity solution for the temperature field. They investigated the effect of the Prandtl number and the unsteadiness parameter which determine the velocity and temperature profiles and heat transfer coefficient. Mahmood et al. [4] studied flow and heat transfer over a permeable sensor surface placed in a squeezing channel. They displayed effects of the squeezing parameter on values of velocity and temperature and the variation of skin friction coefficient and heat transfer coefficient against the transpiration parameter with the graphics. Bataller [5] investigated the flow influenced by nonlinearly stretching of the sheet with heat transfer. They analyzed the effects of various physical parameters on heat transfer phenomena in a viscous flow over a nonlinearly stretching sheet. Mukhopadhyay [6] investigated the effects of thermal radiation on unsteady boundary layer mixed convection heat transfer problem from a vertical porous stretching surface embedded in a porous medium. Their study indicated that the flow and temperature field are significantly influenced by the unsteadiness parameter, buoyancy force, and suction parameter in both porous and nonporous media. Aziz [7] studied thermal boundary layer flow over a flat plate considering convective surface heat flux at the lower surface of the plate and determined the condition which the convection heat transfer coefficient must meet for a similarity solution to exist. Aziz [8] studied the boundary layer flow over a flat plate with slip flow and constant heat flux surface condition velocity and shear stress distributions were presented for a range of values of the parameter characterizing the slip flow. Ishak [9] considered the problem of steady laminar boundary layer flow and heat transfer over a stationary flat plate immersed in a uniform free stream with convective boundary condition. They found that suction increases the surface shear stress and as a consequence increases the heat transfer rate at the surface. Yao et al. [10] investigated heat transfer flow over a stretching/shrinking sheet with convective boundary conditions. They found that temperature at the surface of the plate increases with the increase of the surface convection parameter. Rahman [11] studied locally similar solutions for the hydromagnetic thermal boundary layer flow of a viscous incompressible electrically conducting fluid over a flat plate with a partial slip at the surface of boundary in the presence of the convective surface boundary condition numerically and displayed results for the velocity, temperature, and Prandtl number within the boundary layer delineating the effect of various parameters characterizing the flow.

In most of the published studies in the literature [1216], the partial differential equations were converted to ordinary differential equations via a similarity transformation found by ad hoc methods. The resulting ordinary differential equations were then solved by numerical methods. Lie group methods are performed to investigate similarity transformation of the boundary layer equations in only a few studies [1719]. However, a detailed investigation of the complete symmetries for unsteady boundary layer flow and heat transfer equations are lacking in the literature.

In this paper, Lie group analysis of unsteady flow and heat transfer over a porous surface in the presence of suction/injection velocity is analyzed. Group classification of the governing equations with respect to the variable free stream flow and suction/injection velocity is performed. The infinitesimal generators are obtained. The restrictions imposed by the boundary conditions on the generators are discussed. Using the resulting generators, the governing equations are transformed into ordinary differential equations and solved numerically. Numerical results are presented for the selected dimensionless problem parameters and discussed in detail. Present study reveals that injection decreases temperature gradient whereas suction increases the temperature gradient. The increase in the Prandtl number causes increase in the temperature gradient with increasing the slope.

#### 2. Mathematical Formulation

Consider the unsteady two-dimensional, incompressible laminar flow over a porous surface. Under the boundary layer approximation, the continuity, momentum, and energy equations can be written as where is the coordinate along the surface and is the coordinate vertical to . is time, and are the velocity components in and coordinates. is the free stream velocity outside the boundary layer, is the temperature, and are the kinematic viscosity and thermal diffusivity of the fluid, respectively.

The boundary conditions are where is the suction/injection velocity of the permeable surface. is the constant surface temperature and is the temperature of the ambient fluid.

The nondimensional form of (2.1) and boundary conditions are as follows: The nondimensional parameters are related to the dimensional ones through the following relations: where , , , and are the average free stream velocity, length of the horizontal surface, Reynolds number, and Prandtl number, respectively. is used for simplicity in the equations.

#### 3. Symmetry Analysis

The symmetry groups of (2.3)–(2.5) are calculated by using classical Lie group approach. Details of the theory can be found in Bluman and Kumei [20] and Stephani [21]. The one-parameter infinitesimals Lie group of transformations are defined as The prolonged infinitesimal generator which includes higher order derivatives is The invariance conditions for (2.3)–(2.5) are obtained by applying the generator to the equations: MathLie program [22] is used to calculate higher order infinitesimals and then equations are separated with respect to higher order variables. The initial restrictions on the infinitesimals are and the remaining determining equations are Equations (3.5) are an overdetermined partial differential system. Solving the system, the infinitesimals are finally obtained: is to be determined from the following equation: The generator is applied to the boundary conditions which should also be invariant under the transformation leading to restrictions in (3.6): or and functions should satisfy the following equations: Different functions can be obtained from (3.10) for and via setting the parameters as follows.

(I) Parameter a: selecting parameter without a loss of generality and all other parameters zero yields In this case, suction or injection velocity is an arbitrary function of time. Solving (3.11) for yields , where is an arbitrary function of time for this special choice.

(II) Parameter c: selecting parameter and all other parameters zero yields For this choice, and yield and from (3.10).

(III) Parameter b: selecting parameter and all other parameters zero yields In this case, suction or injection velocity and free stream velocity are arbitrary functions of time.

(IV) Parameter e: selecting parameter and all other parameters zero yields In this case, suction or injection velocity and free stream velocity are arbitrary functions of the spatial variable .

#### 4. Similarity Transformations and Reductions

Similarity transformations and reductions of the partial differential system to an ordinary differential system will be considered.

##### 4.1. A Similarity Transformation

In this section, we choose the first case (case I) for obtaining similarity transformations as follows: Similarity variables and functions are obtained from (4.1): The substitution of these variables and functions into the original equations and boundary conditions yields The partial differential system with three independent variables is reduced to a system with two independent variables.

##### 4.2. The Second Reduction

To reduce (4.3) and (4.4) into an ordinary differential system, one-parameter infinitesimal Lie group of transformations is defined as The prolonged infinitesimal generator is The invariance conditions for (4.3) are Substituting the higher order infinitesimals, separating the equations with respect to higher order variables yields from which the above equations are solved: The determining equation for which is selected to , outer velocity, is obtained as follows: The boundaries and boundary conditions should be invariant under the transformations also (Bluman and Kumei [20]): Two parameters are involved in Lie group Transformations.

and functions should satisfy the following equations: Different functions can be obtained from the above equations for and via setting the parameters. The general case in which both parameters are involved will be considered.

Solving (4.14) results The equations determining the similarity transformations are from which the similarity variables and functions are derived: Inserting these variables and functions into (4.3) and (4.4) yields the following ordinary differential system: Further, we scale variables and as The substitution of these variables into (4.18) yields The associated boundary conditions are where is the suction/injection parameter; and are the unsteadiness parameters.

#### 5. Result and Discussion

The nonlinear reduced ordinary differential equation system (4.20) with the boundary conditions in (4.21) is solved numerically by applying a finite difference code that implements the three-stage Lobatto IIIA formula that provides continuous solutions of fourth order accuracy in the problem domain. Description of the numerical method is given by Shampine et al. [23] and Kierzenka and Shampine [24]. Results are displayed graphically for various values of problem parameters to see developments in the velocity and temperature fields.

The variation of the dimensionless -component of velocity is shown for various values of unsteadiness parameter in the absence of suction/injection in Figure 1. It is clear that an increase in the causes a reduction in the velocity until a certain point. After this point, reverse effect in velocity is observed for increasing values of .

Figure 1: Variation of for different values of unsteadiness parameter (, ).

The effect of suction/injection parameter on the velocity profile is presented in Figure 2. It is noted that suction decreases the boundary layer thickness whereas injection increases.

Figure 2: Variation of for different values of suction/injection (, ).

Effects of Prandtl number on the temperature profiles are given in Figure 3. An increase in the Prandtl number decreases the temperature profiles. Effects of unsteadiness parameter on the temperature profiles are depicted in Figure 4. An increase in decreases the temperature profiles. Thermal boundary layer is narrower for higher unsteadiness parameter. In Figure 5 effects of suction-injection parameter on the temperature profiles are depicted. corresponds to no suction/injection case. Suction decreases the thickness of thermal boundary layer whereas injection increases the thickness of thermal boundary layer.

Figure 3: Variation of for different values of the Prandtl number (, ).
Figure 4: Variation of for different values of (, ).
Figure 5: Variation of for different values of suction/injection (, ).

In Figure 6, we displayed the velocity gradient against suction/injection parameter for the case of flow along a horizontal porous surface for different values of the unsteadiness parameter, . Injection decreases the velocity gradient whereas suction increases the velocity gradient. The reduction in causes an increase in velocity gradient with decreasing the slope.

Figure 6: Variation of the velocity gradient with ().

In Figure 7, we show the temperature gradient against suction/injection parameter for different -unsteadiness parameter. Injection decreases the temperature gradient whereas suction increases the temperature gradient. In Figure 8, we displayed the temperature gradient against unsteadiness parameter for different Prandtl numbers. The increase in the Prandtl number causes increase in temperature gradient.

Figure 7: Variation of the temperature gradient with ().
Figure 8: Variation of the temperature gradient with ().

#### 6. Conclusions

In this paper, we have studied the problem of two-dimensional, unsteady flow and heat transfer of a viscous fluid past a surface in the presence of variable suction/injection. The unsteadiness is due to time dependent free stream flow.

The governing equations are derived with the usual boundary layer approximation. Using Lie group theory, a group classification of the equations with respect to the variable free stream flow and suction/injection velocity is performed. Restrictions imposed by the boundary conditions on the symmetries are discussed. Adopting the obtained symmetry groups, governing partial differential equations are converted into ordinary differential equations and then solved numerically.

Effects of dimensionless problem parameters on the velocity and temperature profiles are outlined in the figures. From the present investigation the following results are obtained.(i)An increase in the unsteadiness parameter causes a reduction in the velocity.(ii)Suction decreases the boundary layer thickness whereas injection increases.(iii)An increase in the Prandtl number and unsteadiness parameter decreases the temperature profiles.(iv)Suction decreases the thickness of thermal boundary layer whereas injection increases the thickness of boundary layer. (v)Injection decreases the velocity gradient whereas suction increases. Moreover the reduction of causes an increase in the velocity gradient.(vi)Injection decreases the temperature gradient whereas suction increases the temperature gradient. The increase in the Prandtl number causes increase in the temperature gradient with increasing the slope.

#### References

1. H. Blasius, “Grenzschichten in Flussigkeiten mit kleiner Reibung,” Zeitschrift für Mathematik und Physik, vol. 56, no. 1, pp. 1–37, 1908.
2. A. R. A. Khaled and K. Vafai, “Heat transfer and hydromagnetic control of flow exit conditions inside oscillatory squeezed thin films,” Numerical Heat Transfer A, vol. 43, no. 3, pp. 239–258, 2003.
3. E. M. A. Elbashbeshy and M. A. A. Bazid, “Heat transfer over an unsteady stretching surface,” Heat and Mass Transfer, vol. 41, no. 1, pp. 1–4, 2004.
4. M. Mahmood, S. Asghar, and M. A. Hossain, “Squeezed flow and heat transfer over a porous surface for viscous fluid,” Heat and Mass Transfer, vol. 44, no. 2, pp. 165–173, 2007.
5. R. C. Bataller, “Similarity solutions for flow and heat transfer of a quiescent fluid over a nonlinearly stretching surface,” Journal of Material Processing Technology, vol. 203, no. 1–3, pp. 176–183, 2008.
6. S. Mukhopadhyay, “Effect of thermal radiation on unsteady mixed convection flow and heat transfer over a porous stretching surface in porous medium,” Heat and Mass Transfer, vol. 52, no. 13-14, pp. 3261–3265, 2009.
7. A. Aziz, “A similarity solution for laminar thermal boundary layer over a flat plate with a convective surface boundary condition,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1064–1068, 2009.
8. A. Aziz, “Hydrodynamic and thermal slip flow boundary layers over a flat plate with constant heat flux boundary condition,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 3, pp. 573–580, 2010.
9. A. Ishak, “Similarity solutions for flow and heat transfer over a permeable surface with convective boundary condition,” Applied Mathematics and Computation, vol. 217, no. 2, pp. 837–842, 2010.
10. S. Yao, T. Fang, and Y. Zhong, “Heat transfer of a generalized stretching/shrinking wall problem with convective boundary conditions,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 2, pp. 752–760, 2011.
11. M. M. Rahman, “Locally similar solutions for hydromagnetic and thermal slip flow boundary layers over a flat plate with variable fluid properties and convective surface boundary condition,” Meccanica, vol. 46, no. 5, pp. 1127–1143, 2011.
12. A. Ishak, R. Nazar, and I. Pop, “Heat transfer over an unsteady stretching permeable surface with prescribed wall temperature,” Nonlinear Analysis. Real World Applications, vol. 10, no. 5, pp. 2909–2913, 2009.
13. G. C. Layek, S. Mukhopadhyay, and S. A. Samad, “Heat and mass transfer analysis for boundary layer stagnation-point flow towards a heated porous stretching sheet with heat absorption/generation and suction/blowing,” International Communications in Heat and Mass Transfer, vol. 34, pp. 347–356, 2007.
14. R. Tsai, K. H. Huang, and J. S. Huang, “Flow and heat transfer over an unsteady stretching surface with non-uniform heat source,” International Communications in Heat and Mass Transfer, vol. 35, no. 10, pp. 1340–1343, 2008.
15. S. Sivasankaran, M. Bhuvaneswari, P. Kandaswamy, and E. K. Ramasami, “Lie group analysis of natural convection heat and mass transfer in an inclined surface,” Lithuanian Association of Nonlinear Analysis, vol. 11, no. 2, pp. 201–212, 2006.
16. S. Sivasankaran, M. Bhuvaneswari, P. Kandaswamy, and E. K. Ramasami, “Lie group analysis of natural convection heat and mass transfer in an inclined porous surface with heat generation,” International Journal of Applied Mathematics and Mechanics, vol. 2, no. 1, pp. 34–40, 2006.
17. M. Pakdemirli, “Similarity analysis of boundary layer equations of class of non-Newtonian fluids,” International Journal of Nonlinear Mechanics, vol. 29, no. 2, pp. 187–196, 1994.
18. M. Pakdemirli, M. Yürüsoy, and A. Küçükbursa, “Symmetry groups of boundary layer equations of class of non-Newtonian fluids,” International Journal of Nonlinear Mechanics,, vol. 31, no. 3, pp. 267–276, 1996.
19. M. Yürüsoy and M. Pakdemirli, “Symmetry reductions of unsteady three-dimensional boundary layers of some non-Newtonian fluids,” International Journal of Engineering Science, vol. 35, no. 8, pp. 731–740, 1997.
20. G. W. Bluman and S. Kumei, Symmetries and Differential Equations, Springer, New York, NY, USA, 1989.
21. H. Stephani, Differential Equations: Their Solution Using Symmetries, Cambridge University Press, New York, NY, USA, 1989.
22. Wolfram Mathematica, MathLie 3.0, Dr. Gerd Baumann, 1992.
23. L. F. Shampine, I. Gladwell, and S. Thompson, Solving ODEs with MATLAB, chapter 3, Cambridge University Press, 2003.
24. J. Kierzenka and L. F. Shampine, “A BVP solver based on residual control and the MATLAB PSE,” ACM Transactions on Mathematical Software, vol. 27, no. 3, pp. 299–316, 2001.