Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 724385, 5 pages

http://dx.doi.org/10.1155/2013/724385

Research Article

Prandtl's Boundary Layer Equation for Two-Dimensional Flow: Exact Solutions via the Simplest Equation Method

Taha Aziz,¹ A. Fatima,¹ C. M. Khalique,² and F. M. Mahomed¹

¹Centre for Differential Equations, Continuum Mechanics and Applications, School of Computational and Applied Mathematics, University of the Witwatersrand, Wits 2050, South Africa

²Department of Mathematical Science, International Institute of Symmetry Analysis and Mathematical Modeling, North-West University, Mafikeng Campus, Private Bag X 2046, Mmabatho 2735, South Africa

Received 12 March 2013; Accepted 18 March 2013

Academic Editor: H. Jafari

Copyright © 2013 Taha Aziz 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 simplest equation method is employed to construct some new exact closed-form solutions of the general Prandtl's boundary layer equation for two-dimensional flow with vanishing or uniform mainstream velocity. We obtain solutions for the case when the simplest equation is the Bernoulli equation or the Riccati equation. Prandtl's boundary layer equation arises in the study of various physical models of fluid dynamics. Thus finding the exact solutions of this equation is of great importance and interest.

1. Introduction

Many scientific and engineering problems and phenomena are modeled by nonlinear differential equations. Therefore, the study of nonlinear differential equations has been an active area of research from the past few decades. Considerable attention has been devoted to the construction of exact solutions of nonlinear equations because of their important role in the study of nonlinear physical models. For nonlinear differential equations, we do not have the freedom to compute exact (closed-form) solutions and for analytical work we have to rely on some approximate analytical or numerical techniques which may be helpful for us to understand the complex physical phenomena involved. The exact solutions of the nonlinear differential equations are of great interest and physically more important. These exact solutions, if reported, facilitate the verification of complex numerical codes and are also helpful in a stability analysis for solving special nonlinear problems. In recent years, much attention has been devoted to the development of several powerful and useful methods for finding exact analytical solutions of nonlinear differential equations. These methods include the powerful Lie group method [1], the sine-cosine method [2], the tanh method [3, 4], the extended tanh-function method [5], the Backlund transformation method [6], the transformed rational function method [7], the -expansion method [8], the exponential function rational expansion method [9], and the Adomian's decomposition method [10].

Prandtl [11] introduced boundary layer theory in 1904 to understand the flow behavior of a viscous fluid near a solid boundary. Prandtl gave the concept of a boundary layer in large Reynolds number flows and derived the boundary layer equations by simplifying the Navier-Stokes equations to yield approximate solutions. Prandtl's boundary layer equations arise in various physical models of fluid mechanics. The equations of the boundary layer theory have been the subject of considerable interest, since they represent an important simplification of the original Navier-Stokes equations. These equations arise in the study of steady flows produced by wall jets, free jets, and liquid jets, the flow past a stretching plate/surface, flow induced due to a shrinking sheet, and so on. These boundary layer equations are usually solved subject to certain boundary conditions depending upon the specific physical model considered. Blasius [12] solved the Prandtl's boundary layer equations for a flat moving plate problem and gave a power series solution of the problem. Falkner and Skan [13] generalized the Blasius boundary layer problem by considering the boundary layer flow over a wedge inclined at a certain angle. Sakiadis [14] initiated the study of the boundary layer flow over a continuously moving rigid surface with a uniform speed. Crane [15] was the first one who studied the boundary layer flow due to a stretching surface and developed the exact solutions of boundary layer equations with parameter . P. S. Gupta and A. S. Gupta [16] extended the Crane's work and for the first time introduced the concept of heat transfer with the stretching sheet boundary layer flow. The numerical solution for a free two-dimensional jet was obtained by Schlichting [17] and later an analytic study was made by Bickley [18]. Riley [19] derived the solution for a radial liquid jet. Recently, the similarity solution of axisymmetric non-Newtonian wall jet with swirl effects was investigated by Kolář [20]. Naz et al. [21] and Mason [22] have investigated the general boundary layer equations for two-dimensional and radial flows by using the classical Lie group approach and very recently Naz et al. [23] have provided the similarity solutions of the Prandtl's boundary layer equations by implementing the nonclassical/conditional symmetry method.

The simplest equation method is a powerful mathematical tool for finding exact solutions of nonlinear ordinary differential equations. It has been developed by Kudryashov [24, 25] and used successfully by many researchers for finding exact solutions of nonlinear ordinary differential equations [26–28]. The purpose of the present work is to find the exact closed-form solutions of Prandtl's boundary layer equation for two-dimensional flow with constant or uniform main stream velocity by the use of simplest equation method.

Prandtl's boundary layer equation for the stream function for an incompressible, steady two-dimensional flow with uniform or vanishing mainstream velocity is [29]

Here denote the Cartesian coordinates parallel and perpendicular to the boundary and is the kinematic viscosity. The velocity components and , in the and directions, are related to stream function as

By the use of Lie group theoretic method of infinitesimal transformations [1], the general form of similarity solution for (1) is

where is the constant determined from the further conditions and is the similarity variable. By the substitution of (3) into (1), the third-order nonlinear ordinary differential equation in results, namely, in

Equation (4) gives the general form of Prandtl's boundary layer equation for two-dimensional flow of a viscous incompressible fluid. The boundary layer equation (4) is usually solved subject to certain boundary conditions depending upon the particular physical model considered. Here, we present the exact closed-form solutions of (4) using the simplest equation method. We organize the paper as follows. In Section 2, we describe briefly the simplest equation method. In Section 3, we apply this method to solve nonlinear Prandtl's boundary layer equation for two-dimensional flow. Finally, some closing remarks are presented in Section 4.

2. A Description of the Simplest Equation Method

Here we present a brief description of the simplest equation method for solving nonlinear ordinary differential equations.

Step 1. We first consider a general form of a nonlinear ordinary differential equation:

where is the dependent variable and is the independent variable.

Step 2. The basic idea of the simplest equation method consists in expanding the solutions of the previous ordinary differential equation in a finite series:

where is a solution of some ordinary differential equations. These ordinary differential equations are called the simplest equations. The main property of the simplest equation is that we know the general solution of the simplest equation or we at least know the exact analytical solutions of the simplest equation. The parameters are to be determined from the further conditions.

In this paper we use the Bernoulli and Riccati equations as the simplest equations. These equations are well-known nonlinear ODEs whose solutions can be expressed in terms of elementary functions.

For the Bernoulli equation

where and are constants independent of . The solution of (7) is

where is a constant of integration.

For the Riccati equation

where , , and are constants, we will use the solutions

where

and is a constant of integration.

Step 3. One of the main steps in using the simplest equation method is to determine the positive number in (6). The positive number can be determined by considering the homogeneous balance between the highest order derivatives and nonlinear terms appearing in (5).

Step 4. By the substitution of (6) into (5) and with (7) or (9), the left hand side of (5) is converted into a polynomial in . Equating each coefficient of the polynomial to zero yields a set of algebraic equations for , , , .

Step 5. By solving the algebraic equations obtained in Step 4 and substituting the results into (6), we obtain the exact solutions of ODE (5).

3. Application of the Simplest Equation Method

In this section, we employ the simplest equation method and obtain exact closed-form solutions of Prandtl's boundary layer equation (4).

3.1. Solutions of Boundary Layer Equation Using the Equation of Bernoulli as the Simplest Equation

The balancing procedure yields . Thus we search for a solution of (4) of the form

where satisfies the Bernoulli equation and and are the parameters to be determined.

By the substitution of (12) into (4) and making use of the Bernoulli equation (7) and then equating all coefficients of the functions to zero, we obtain an algebraic system of equations in terms of and . Solving this system of algebraic equations, we obtain the values of the constants and . Therefore the solution of Prandtl's boundary layer equation (4) with is given by

and hence the corresponding stream function becomes

Special Cases. By taking and in the previous solution, we obtain a special solution given by

Likewise, if we take and , we deduce

3.2. Solutions of Boundary Layer Equation Using the Equation of Riccati as the Simplest Equation

The balancing procedure yields . Thus the solution of (4) is written in the form

By the insertion of (17) into (4) and making use of the Riccati equation (9) and proceeding as above, we obtain algebraic system of equations in terms of and . Solving this system, we obtain the solutions of Prandtl's boundary layer equation (4) for as

and the solutions for corresponding stream functions are

where and is a constant of integration.

By taking , , and in (19), we deduce a special solution of stream function , given by

4. Concluding Remarks

In this study, we have utilized the method of simplest equation for obtaining exact closed-form solutions of the well-known Prandtl's boundary layer equation for two-dimensional flow with uniform mainstream velocity. As the simplest equations, we have used the Bernoulli and Riccati equations. Prandtl's boundary layer equations arise in various physical models of fluid dynamics and thus the exact solutions obtained may be very useful and significant for the explanation of some practical physical models dealing with Prandtl's boundary layer theory. We have also verified that the solutions obtained here are indeed the solutions of Prandtl's boundary layer equation.

Acknowledgments

T. Aziz and A. Fatima would like to thank the School of Computational and Applied Mathematics and the Financial Aid and Scholarship Office, University of the Witwatersrand, for financial support and research grant.

References

P. J. Olver, Applications of Lie Groups to Differential Equations, vol. 107 of Graduate Texts in Mathematics, Springer, New York, NY, USA, 2nd edition, 1993.

A.-M. Wazwaz, “The sine-cosine method for obtaining solutions with compact and noncompact structures,” Applied Mathematics and Computation, vol. 159, no. 2, pp. 559–576, 2004.

W. Malfliet, “Solitary wave solutions of nonlinear wave equations,” American Journal of Physics, vol. 60, no. 7, pp. 650–654, 1992.

A.-M. Wazwaz, “The tanh method: solitons and periodic solutions for the Dodd-Bullough-Mikhailov and the Tzitzeica-Dodd-Bullough equations,” Chaos, Solitons & Fractals, vol. 25, no. 1, pp. 55–63, 2005.

E. Fan, “Extended tanh-function method and its applications to nonlinear equations,” Physics Letters A, vol. 277, no. 4-5, pp. 212–218, 2000.

M. R. Miura, Backlund Transformation, Springer, Berlin, Germany, 1978.

W.-X. Ma and J.-H. Lee, “A transformed rational function method and exact solutions to the () dimensional Jimbo-Miwa equation,” Chaos, Solitons & Fractals, vol. 42, no. 3, pp. 1356–1363, 2009.

R. Abazari, “Application of ()-expansion method to travelling wave solutions of three nonlinear evolution equation,” Computers & Fluids, vol. 39, no. 10, pp. 1957–1963, 2010.

H. Xin, “The exponential function rational expansion method and exact solutions to nonlinear lattice equations system,” Applied Mathematics and Computation, vol. 217, no. 4, pp. 1561–1565, 2010.

J.-L. Li, “Adomian's decomposition method and homotopy perturbation method in solving nonlinear equations,” Journal of Computational and Applied Mathematics, vol. 228, no. 1, pp. 168–173, 2009.

L. Prandtl, “Uber Flussigkeitsbewegungen bei sehr kleiner Reibung,” in Verhandlungen des III. Internationalen Mathematiker Kongresses, pp. 484–491, Heidelberg, Germany, 1904.

H. Blasius, “Grenzschichten in Flussigkeiten mit kleiner Reibung,” Zeitschrift für Mathematik und Physik, vol. 56, pp. 1–37, 1908.

V. M. Falkner and S. W. Skan, “Some approximate solutions of the boundary layer equations,” Philosophical Magazine, vol. 12, pp. 865–896, 1931.

B. C. Sakiadis, “Boundary-layer behavior on continuous solid surface. I. Boundary-layer equations for two-dimensional and axisymmetric flow,” AIChE Journal, vol. 7, pp. 26–28, 1961.

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

P. S. Gupta and A. S. Gupta, “Heat and mass transfer on a stretching sheet with suction and blowing,” The Canadian Journal of Chemical Engineering, vol. 55, pp. 744–746, 1977.

H. Schlichting, “Laminare Strahlausbreitung,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 13, pp. 260–263, 1933.

W. G. Bickley, “The plane jet,” Philosophical Magazine, vol. 23, pp. 727–731, 1937.

N. Riley, “Asymptotic expansions in radial jets,” Journal of Mathematical Physics, vol. 41, pp. 132–146, 1962.

V. Kolář, “Similarity solution of axisymmetric non-Newtonian wall jets with swirl,” Nonlinear Analysis: Real World Applications, vol. 12, no. 6, pp. 3413–3420, 2011.

R. Naz, F. M. Mahomed, and D. P. Mason, “Symmetry solutions of a third-order ordinary differential equation which arises from Prandtl boundary layer equations,” Journal of Nonlinear Mathematical Physics, vol. 15, no. 1, pp. 179–191, 2008.

D. P. Mason, “Group invariant solution and conservation law for a free laminar two-dimensional jet,” Journal of Nonlinear Mathematical Physics, vol. 9, no. 2, pp. 92–101, 2002.

R. Naz, M. D. Khan, and I. Naeem, “Nonclassical symmetry analysis of boundary layer equations,” Journal of Applied Mathematics, vol. 2012, Article ID 938604, 7 pages, 2012.

N. A. Kudryashov, “Simplest equation method to look for exact solutions of nonlinear differential equations,” Chaos, Solitons & Fractals, vol. 24, no. 5, pp. 1217–1231, 2005.

N. A. Kudryashov, “Exact solitary waves of the Fisher equation,” Physics Letters A, vol. 342, no. 1-2, pp. 99–106, 2005.

H. Jafari, N. Kadkhoda, and C. M. Khalique, “Travelling wave solutions of nonlinear evolution equations using the simplest equation method,” Computers & Mathematics with Applications, vol. 64, no. 6, pp. 2084–2088, 2012.

N. Taghizadeh and M. Mirzazadeh, “The simplest equation method to study perturbed nonlinear Schrödinger's equation with Kerr law nonlinearity,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 4, pp. 1493–1499, 2012.

N. K. Vitanov and Z. I. Dimitrova, “Application of the method of simplest equation for obtaining exact traveling-wave solutions for two classes of model PDEs from ecology and population dynamics,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 10, pp. 2836–2845, 2010.

L. Rosenhead, Laminar Boundary Layers, pp. 254–256, Clarendon Press, Oxford, UK, 1963.