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

Authors: Taha Aziz; A. Fatima; C. M. Khalique; F. M. Mahomed
Journal: Mathematical Problems in Engineering
 (2013)

---

## 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.


---

## Body

## 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 (G′/G)-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 γ=0. 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ψ(x,y) for an incompressible, steady two-dimensional flow with uniform or vanishing mainstream velocity is [29]
(1)∂ψ∂y∂2ψ∂x∂y-∂ψ∂x∂2ψ∂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 [1], 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=0M‍Ai(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 numberM 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 Step4 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 yieldsM=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 taked=-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 yieldsM=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 takingd=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))].


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

The balancing procedure yieldsM=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 taked=-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 yieldsM=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 takingd=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.



---

*Source: 724385-2013-04-10.xml*