A Comparison between the Reduced Differential Transform Method and Perturbation-Iteration Algorithm for Solving Two-Dimensional Unsteady Incompressible Navier-Stokes Equations

In this work, approximate analytical solutions to the lid-driven square cavity flow problem, which satisfied two-dimensional unsteady incompressible Navier-Stokes equations, are presented using the kinetically reduced local Navier-Stokes equations. Reduced differential transform method and perturbation-iteration algorithm are applied to solve this problem. The convergence analysis was discussed for both methods. The numerical results of both methods are given at some Reynolds numbers and low Mach numbers, and compared with results of earlier studies in the review of the literatures. These two methods are easy and fast to implement, and the results are close to each other and other numerical results, so it can be said that these methods are useful in finding approximate analytical solutions to the unsteady incompressible flow problems at low Mach numbers.


Introduction
Fluid flow is one of the most important engineering phenomena that have received widespread attention in theoretical and practical scientific research.Many of these studies focus on simulated mathematical models which represent these phenomena.Therefore, the equations of Navier-Stokes, which are the basic model for describing the movement of fluid, have received considerable attention from researchers to find their analytical and numerical solutions.
In this work, unsteady viscous incompressible flows characterized by two-dimensional Navier-Stokes equations are studied.The non-dimensional momentum and continuity equations have the following form where t is the physical time, ( ) , , u x y t and ( ) , , v x y t are the fluid velocity components, ( ) , , p x y t is the pressure, and Re is the Reynolds number.Since, the Navier-Stokes equations are nonlinear partial differential equations and there is no explicit equation for calculating pressure, these equations are difficult to solve, so many studies have suggested the alternative thermodynamic description of incompressible fluid flows.One of these alternative formulas is the kinetically reduced local Navier-Stokes (KRLNS) equations [1] [2] [3] [4] [5] which is obtained by replacing the pressure by where Ma is the Mach number and ( ) , , g x y t is the grand potential.The time scale in INS equations is related to that of KRLNS equations; ( ) Then, the system of equations of KRLNS has the following form The KRlNS equations suggested in [1] of the reduced equations for the grand potential and the fluid momentum were derived from the compressible Navier-Stokes equations in order to present the thermodynamic description of incompressible fluid flows at low Mach numbers.The two-dimensional KRLNS system is simplified and compared with a Chorin's artificial compressibility method for steady state computation of flow in two-dimensional lid-driven cavity and Taylor-Green vortex flow in [2].In [3], KRLNS equations were applied to two-dimensional simulation of doubly periodic shear layers and decaying ho-is the first.In Section 2 and 3, we describe the reduced differential transform method and perturbation-iteration algorithm, and applied them to KRLNS equations.We derived the condition of convergence for both methods (Section 4).
We then present the approximate analytical solutions for two-dimensional lid-driven cavity flow, which are obtained by applying differential transform method and perturbation-iteration algorithm (Section 5).Next, we introduce the numerical results and compare these results with other works (Section 6).The last Section summarizes the major findings of this study.

Reduced Differential Transform Method (RDTM)
The RDTM was first introduced by Keskin [13].It is an iterative procedure based on the use of the Taylor series solution of differential equations.It has been successfully applied to solve various nonlinear partial differential equations [13]- [27].Since it does not require any parameter, discretization, linearization or small perturbations, thus it reduces the size of computations and can be easily used.The RDTM was used for solving the generalized Korteweg-de Vries equation [14], the fractional Benney-Lin equation [15], the Wu-Zhang equation [16], the equal width wave equation and the inviscid Burgers equation [17], the Sine-Gordon equation [18], the Burgers and Huxley equations [19], the time-fractional telegraph equation [20], the generalized Drinfeld-Sokolov equations and Kaup-Kupershmidt equation [21], the Zakharov-Kuznetsov equations [22], the heat-like equations [23], the coupled Ramani equations [25], two integral members of nonlinear Kadomtsev-Petviashvili hierarchy equations [26], and the second order hyperbolic telegraph equation [27].Few studies have been applied RDTM to solve the Navier-Stocks equations, which is one of the reasons for choosing it as a method for solving the lid-driven cavity flow.

X
x y = be a vector, if ( ) , u X t is analytic function and continuously differentiable with re- spect to time t and space in the domain of interest.Then, let is the t-dimensional spectrum function of ( ) , u X t which is the transformed function.The reduced differential inverse transform of ( ) from Equation (2.1) and Equation (2.2), we can conclude that ( ) ( ) 3) The fundamental mathematical operations performed by RDTM are readily obtained and listed in Table 1.
In order to apply this method with KRLNS equations to find approximate analytical solutions for INS equations, we suppose that ( ) , where

( )
, u X t and ( ) , v X t are the fluid velocity components in the x and y directions, and ( ) , g X t respectively.Then, we have , , ( ) Then the exact solution is obtained as follows: where This approach is referred to by (KRDTM) in this paper.

Perturbation-Iteration Algorithm (PIA)
Perturbation methods are important analytical methods which have been used to construct approximate analytical solutions of algebraic equations, differential equations, and integro-differential equations.The main limitation of using the perturbation methods is to install a small auxiliary parameter in the equation.
For this reason, the solutions of these methods are restricted by validity range of physical parameters, so many of perturbation techniques have been suggested by several authors.PIA is one of the techniques which was proposed by Pakdemirli and Boyac in [28], and used a combination of perturbation expansions and Taylor series expansions to construct an iteration scheme for using to generate root finding algorithms.It is applied by many authors to get the approximate analytical solution for differential equations.In [29], PIA was applied to obtain the solution of Bratu-type equations.In [30], PIA was utilized to find the solution first order differential equations.This algorithm was tested on three nonlinear heat equations in [31].Moreover, PIA was generalized to an arbitrary number of first-order coupled equations in [32].It was applied to Fredholm and Volterra integral equations in [33].Also, in [34], PIA was proposed for solving the Riccati differential equation.It was developed in [35] to obtain the solutions of Lotka-Volterra differential equations.In [36], some types of fractional differential equation systems were solved by using this method.PIA with Laplace transform method was combined in [37] to solve Newell-Whitehead-Segel equations.
In [38], PIA is used for solving the fractional Zakharov-Kuznetsov equation and compared with the residual power series method.By reviewing the previous literature, we have not found any research that has used this method to find a solution to the two-dimensional lid-driven cavity flow problem and which is an important reason to use this method to solve this problem.
In general, PIA is obtained by taking different numbers of terms in the perturbation expansions and different order of correction terms in the Taylor series expansions.Therefore, the perturbation-iteration algorithm is called PIA(m,n) where the m is the number of the correction terms in the perturbation expansion and n is the highest order derivative term in the Taylor series such that m should always be less than or equal to n.
To obtain approximate analytical solutions for two-dimensions Navier-Stokes equations, PIA (1, 1) will be applied to KRLNS equations and which will be referred to this article by (KPIA).Firstly, we write Equation (1.5) as follows: ( ) where  is a small perturbation parameter.Secondly, we define the following perturbation expansions with only one correction term: ( ) where n represents the n_th iteration and c u , c v and c g are the correction terms in the perturbation expansion.Thirdly, by replacing (3.2) into (3.1) and writing in the Taylor series expansion for first order derivative terms about , Finally, by substituting the above derivative in the formulas (3.3) and setting we obtain the following iteration equation formulas: .
The calculations start with initial condition ( ) ) 2) to obtain 1 u , 1 v and 1 g which are the solutions at the first iteration.So we can get ( ) iteration solutions by repeating this process and using the previous solution n as an initial guess.

Analysis of Convergence
We now study the convergence analysis of the approximate analytical solutions which are computed from the application KRDTM and KPIA.
Let us consider the Hilbert space We consider the KRINS equation in the following form which is equivalent to the following formula where  is the linear partial derivative with respect to τ ,  is a nonlinear operator,  is a linear operator, and  is a general nonlinear operator involving both linear and nonlinear terms.
Case 1: According to KRDTM, formula (4.1) can be written in the following form and the solutions where ( ) , , It is noted that the solutions by KRDTM is equivalent to determining the sequence ( ) .

S U
Case 2: To study the convergence of KPIA, we write the approximate solutions in different form.To do this, we define ( ) So the solutions, which are resulted from KPIA have the form such that ( ) γ γ γ γ = + + and ( ) is a Cauchy sequence in the Hilbert space 3 H .For this reason, we suppose that ( ) ( ) , .
Then, by using the triangle inequality, we find that ( ) ( ) ( ) , , , , , , ) , , =     be a nonlinear operator satisfies Lipschitz condition from a Hilbert space H 3 into H 3 and ( ) , , , X X τ τ u u , then we have Therefore, from the Banach fixed-point theorem, there is a unique solution of the problem (4.1).Now we have to prove that { } 0 n ∞ S converges to ( ) then we can say that the series approximate solutions { } 0 n ∞

S
converges to the exact solution ( )

The Two-Dimensional Lid-Driven Cavity Flow
In , where ( ) ( ) such that the stream function ψ and vorticity ω are defined as ( ) ( ) Then, by applying KRDTM with the initial conditions of this problem, we obtained the iterative solutions like the form (2.5), such that ( )  To make a decision on the convergence of the KRDTM, we computed k γ as: ( ) ( ) , and 10 Re = then Thus, the iterative solutions (3.2) for this problem, which are obtained by using KPIA, have the following form ) )

Results and Discussion
In this section, we introduce the numerical computations of velocity components u, v, vorticity function w and stream function ψ In Table 2 and Table 3, we reviewed the calculated values of u velocity along the vertical line and v velocity along the horizontal line through the geometric center of the square cavity by using KRDTM and KPIA at 0. Ref. [10] KRDTM KPIA  with the numerical results which have been evaluated by using the finite volume method and introduced by [10].
By comparison, we note that our solutions are remarkably good, and these results represent solutions for the second iteration step.The L ∞ -errors for stream function ψ and vorticity ω are given in Table 5 for Reynolds numbers 10,100 Re = and 1000 at 0.001 Ma = , are compared with the calculated errors by the rational fourth-order compact finite difference method in [43].We note that the calculated errors are small for all values of the Reynolds number and are not affected by the number of grid points.Also, the approximate solutions of velocity in the two directions at cavity center, which are obtained by KRDTM and KPIA, are shown in Figure 1

Conclusion
In this paper, we applied the reduced differential transform method and the perturbation-iteration algorithm on the kinetically reduced local Navier-Stokes  successful techniques and yield remarkably good results to solve unsteady viscous incompressible flow problems at low Mach numbers.Therefore, the application of KRDTM and KPIA could be expanded to include various and multi-dimensions of flow problems.In addition, these methods can be combined with other methods to increase the accuracy of solutions.
where these values are used as estimate values for ( ) 3.4), and then substitute the results of Equation (3.4) into Equation (3.
for both cases.The sufficient condition for convergence of the series solution { } 0 Cauchy sequence in the Hilbert space H 3 , thus, the series solution { } 0 The boundary conditions for the velocities u and v in this problem are of Dirichlet type, which are equal to zero everywhere except along the top surface where To obtain the approximate analytical solutions of the unsteady lid-driven cavity flow problem, we consider the analytical solutions to this problem, which are given in (5.1) as initial conditions for u, v and p.

2 0, 1 .
, which have been obtained by the application of KRDTM and KPIA.All calculations are run by Maple 2017 software with used various values of Reynolds numbers and Mach numbers in the domain [ ] . We observe that the results of KRDTM and KPIA methods at 0.01 Ma = are better than the results at other Mach numbers.Thus, the numerical results of both methods are good and close to each other at low values of Mach numbers.

Table 2 .
The approximate solutions by KRDTM for u and v at

Table 3 .
The approximate solution by KPIA for u and v at

Table 4 .
Comparison between the approximate solutions at

Table 5 .
Comparisons of the L ∞ -errors at