MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 724385 10.1155/2013/724385 724385 Research Article Prandtl's Boundary Layer Equation for Two-Dimensional Flow: Exact Solutions via the Simplest Equation Method Aziz Taha 1 Fatima A. 1 Khalique C. M. 2 Mahomed F. M. 1 Jafari H. 1 Centre for Differential Equations Continuum Mechanics and Applications School of Computational and Applied Mathematics University of the Witwatersrand Wits 2050 South Africa wits.ac.za 2 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 nwu.ac.za 2013 10 4 2013 2013 12 03 2013 18 03 2013 2013 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.

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 , the sine-cosine method , the tanh method [3, 4], the extended tanh-function method , the Backlund transformation method , the transformed rational function method , the (G/G)-expansion method , the exponential function rational expansion method , and the Adomian's decomposition method .

Prandtl  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  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  generalized the Blasius boundary layer problem by considering the boundary layer flow over a wedge inclined at a certain angle. Sakiadis  initiated the study of the boundary layer flow over a continuously moving rigid surface with a uniform speed. Crane  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 γ=0. P. S. Gupta and A. S. Gupta  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  and later an analytic study was made by Bickley . Riley  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ář . Naz et al.  and Mason  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.  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 . 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 ψ(x,y) for an incompressible, steady two-dimensional flow with uniform or vanishing mainstream velocity is  (1)ψy2ψxy-ψx2ψy2-ν3ψy3=0.

Here (x,y) denote the Cartesian coordinates parallel and perpendicular to the boundary y=0 and ν is the kinematic viscosity. The velocity components u(x,y) and v(x,y), in the x and y directions, are related to stream function ψ(x,y) as (2)u(x,y)=ψy,v(x,y)=-ψx.

By the use of Lie group theoretic method of infinitesimal transformations , the general form of similarity solution for (1) is (3)ψ(x,y)=x1-γF(η),η=yxγ, 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 F(η) results, namely, in (4)νd3Fdη3+(1-γ)Fd2Fdη2+(2γ-1)(dFdη)2=0.

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: (5)E[F,dFdη,d2Fdη2,d3Fdη3,]=0, where F 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: (6)F(η)=i=0MAi(G(η))i, where G(η) 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 A0,A1,,AM 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 (7)dGdη=dG(η)+eG2(η), where d and e are constants independent of η. The solution of (7) is (8)G(η)=d[cosh[d(η+C)]+sinh[d(η+C)]1-ecosh[d(η+C)]-esinh[d(η+C)]], where C is a constant of integration.

For the Riccati equation (9)dGdη=dG(η)+eG2(η)+f, where d, e, and f are constants, we will use the solutions (10)G(η)=-d2e-θ2etanh[12θ(η+C)],G(η)=-d2e-θ2etanh(12θη)+sech(θη/2)Ccosh(θη/2)-(2e/θ)sinh(θη/2), where (11)θ2=d2-4ef>0, and C is a constant of integration.

Step 3.

One of the main steps in using the simplest equation method is to determine the positive number M in (6). The positive number M 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 G(η). Equating each coefficient of the polynomial to zero yields a set of algebraic equations for Ai, d, e, f.

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 M=1. Thus we search for a solution of (4) of the form (12)F(η)=A0+A1G(η), where G(η) satisfies the Bernoulli equation and A0 and A1 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 Gi to zero, we obtain an algebraic system of equations in terms of A0 and A1. Solving this system of algebraic equations, we obtain the values of the constants A0 and A1. Therefore the solution of Prandtl's boundary layer equation (4) with γ=2/3 is given by (13)F(η)=-3νd-6νed[cosh[d(η+C)]+sinh[d(η+C)]1-ecosh[d(η+C)]-esinh[d(η+C)]], and hence the corresponding stream function becomes (14)ψ(x,y)c=-3νdx1/3-6νedx1/3c×[-esinh[d(x(-2/3)y+C)])-1](cosh[d(x-2/3y+C)]+sinh[d(x-2/3y+C)])cccc×(1-ecosh[d(x-2/3y+C)]cccccci-esinh[d(x-2/3y+C)])-1].

Special Cases. By taking d=-1 and e=1 in the previous solution, we obtain a special solution given by (15)ψ(x,y)=3νx1/3coth[12(x-2/3y+C)].

Likewise, if we take d=-1 and e=-1, we deduce (16)ψ(x,y)=3νx1/3tanh[12(x-2/3y+C)].

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

The balancing procedure yields M=1. Thus the solution of (4) is written in the form (17)F(η)=A0+A1G(η).

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 A0 and A1. Solving this system, we obtain the solutions of Prandtl's boundary layer equation (4) for γ=2/3 as (18)F(η)=-3νd-6νe[-d2e-θ2etanh(12θ(η+C))],F(η)=-3νd-6νe[-d2e-θ2etanh(ηθ2)i+sech(θη/2)Ccosh(θη/2)-(2e/θ)sinh(θη/2)], and the solutions for corresponding stream functions are (19)ψ(x,y)=-3νdx1/3-6νex1/3[-d2e-θ2etanh(12θ(x-2/3y+C))],(20)ψ(x,y)=-3νdx1/3=-6νex1/3[sinh(θx-2/3y2))-1-d2e-θ2etanh(θx-2/3y2)=-6νex1/3b+(sech(θx-2/3y2))=-6νex1/3b×(Ccosh(θx-2/3y2)-2eθ=-6νex1/3bCco×sinh(θx-2/3y2))-1], where θ2=d2-4ef and C is a constant of integration.

By taking d=3, e=1, and f=1 in (19), we deduce a special solution of stream function ψ, given by (21)ψ(x,y)=-9νx1/3+3νx1/3[3+5tanh(52(x-2/3y+C))].

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.

Olver P. J. Applications of Lie Groups to Differential Equations 1993 107 2nd New York, NY, USA Springer Graduate Texts in Mathematics 10.1007/978-1-4612-4350-2 MR1240056 ZBL0823.58019 Wazwaz A.-M. The sine-cosine method for obtaining solutions with compact and noncompact structures Applied Mathematics and Computation 2004 159 2 559 576 10.1016/j.amc.2003.08.136 MR2096724 ZBL1061.35121 Malfliet W. Solitary wave solutions of nonlinear wave equations American Journal of Physics 1992 60 7 650 654 10.1119/1.17120 MR1174588 ZBL1219.35246 Wazwaz A.-M. The tanh method: solitons and periodic solutions for the Dodd-Bullough-Mikhailov and the Tzitzeica-Dodd-Bullough equations Chaos, Solitons & Fractals 2005 25 1 55 63 10.1016/j.chaos.2004.09.122 MR2123629 Fan E. Extended tanh-function method and its applications to nonlinear equations Physics Letters A 2000 277 4-5 212 218 10.1016/S0375-9601(00)00725-8 MR1827770 ZBL1167.35331 Miura M. R. Backlund Transformation 1978 Berlin, Germany Springer Ma W.-X. Lee J.-H. A transformed rational function method and exact solutions to the (1+3) dimensional Jimbo-Miwa equation Chaos, Solitons & Fractals 2009 42 3 1356 1363 10.1016/j.chaos.2009.03.043 MR2547030 ZBL1198.35231 Abazari R. Application of (G/G)-expansion method to travelling wave solutions of three nonlinear evolution equation Computers & Fluids 2010 39 10 1957 1963 10.1016/j.compfluid.2010.06.024 MR2720168 Xin H. The exponential function rational expansion method and exact solutions to nonlinear lattice equations system Applied Mathematics and Computation 2010 217 4 1561 1565 10.1016/j.amc.2009.06.056 MR2727398 ZBL1207.34097 Li J.-L. Adomian's decomposition method and homotopy perturbation method in solving nonlinear equations Journal of Computational and Applied Mathematics 2009 228 1 168 173 10.1016/j.cam.2008.09.007 MR2514274 ZBL1165.65345 Prandtl L. Uber Flussigkeitsbewegungen bei sehr kleiner Reibung Verhandlungen des III. Internationalen Mathematiker Kongresses 1904 Heidelberg, Germany 484 491 Blasius H. Grenzschichten in Flussigkeiten mit kleiner Reibung Zeitschrift für Mathematik und Physik 1908 56 1 37 Falkner V. M. Skan S. W. Some approximate solutions of the boundary layer equations Philosophical Magazine 1931 12 865 896 Sakiadis B. C. Boundary-layer behavior on continuous solid surface. I. Boundary-layer equations for two-dimensional and axisymmetric flow AIChE Journal 1961 7 26 28 Crane L. J. Flow past a stretching plate Zeitschrift für angewandte Mathematik und Physik 1970 21 4 645 647 2-s2.0-33646099212 10.1007/BF01587695 Gupta P. S. Gupta A. S. Heat and mass transfer on a stretching sheet with suction and blowing The Canadian Journal of Chemical Engineering 1977 55 744 746 Schlichting H. Laminare Strahlausbreitung Zeitschrift für Angewandte Mathematik und Mechanik 1933 13 260 263 Bickley W. G. The plane jet Philosophical Magazine 1937 23 727 731 Riley N. Asymptotic expansions in radial jets Journal of Mathematical Physics 1962 41 132 146 Kolář V. Similarity solution of axisymmetric non-Newtonian wall jets with swirl Nonlinear Analysis: Real World Applications 2011 12 6 3413 3420 10.1016/j.nonrwa.2011.06.004 MR2832981 ZBL1231.35170 Naz R. Mahomed F. M. Mason D. P. Symmetry solutions of a third-order ordinary differential equation which arises from Prandtl boundary layer equations Journal of Nonlinear Mathematical Physics 2008 15 1 179 191 10.2991/jnmp.2008.15.s1.16 MR2434716 Mason D. P. Group invariant solution and conservation law for a free laminar two-dimensional jet Journal of Nonlinear Mathematical Physics 2002 9 2 92 101 10.2991/jnmp.2002.9.s2.8 MR1923373 Naz R. Khan M. D. Naeem I. Nonclassical symmetry analysis of boundary layer equations Journal of Applied Mathematics 2012 2012 7 938604 MR2991602 10.1155/2012/938604 Kudryashov N. A. Simplest equation method to look for exact solutions of nonlinear differential equations Chaos, Solitons & Fractals 2005 24 5 1217 1231 10.1016/j.chaos.2004.09.109 MR2123270 ZBL1069.35018 Kudryashov N. A. Exact solitary waves of the Fisher equation Physics Letters A 2005 342 1-2 99 106 10.1016/j.physleta.2005.05.025 MR2152840 ZBL1222.35054 Jafari H. Kadkhoda N. Khalique C. M. Travelling wave solutions of nonlinear evolution equations using the simplest equation method Computers & Mathematics with Applications 2012 64 6 2084 2088 10.1016/j.camwa.2012.04.004 MR2960827 Taghizadeh N. Mirzazadeh M. The simplest equation method to study perturbed nonlinear Schrödinger's equation with Kerr law nonlinearity Communications in Nonlinear Science and Numerical Simulation 2012 17 4 1493 1499 10.1016/j.cnsns.2011.09.023 MR2855441 ZBL1245.35121 Vitanov N. K. Dimitrova Z. I. 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 2010 15 10 2836 2845 10.1016/j.cnsns.2009.11.029 MR2610628 ZBL1222.35201 Rosenhead L. Laminar Boundary Layers 1963 Oxford, UK Clarendon Press pp. 254–256