A Numerical Computation of a System of Linear and Nonlinear Time Dependent Partial Differential Equations Using Reduced Differential Transform Method

This paper deals with an analytical solution of an initial value system of time dependent linear and nonlinear partial differential equations by implementing reduced differential transform (RDT) method. The effectiveness and the convergence of RDT method are tested by means of five test problems, which indicates the validity and great potential of the reduced differential transform method for solving system of partial differential equations.


Introduction
The reduced differential transform method has been successfully employed to solve various types of linear and nonlinear, homogeneous or nonhomogeneous, equations appearing in science and engineering.Partial differential equations have also been applied in modeling many physical engineering problems and differential equations in nonlinear dynamics [1][2][3].Most of the partial differential equations cannot be solved exactly, and so, developing schemes for getting accurate and efficient numerical solution differential equations have been an active research area.Burgers' equation [4], a system of nonlinear fractional differential equations [1], and nonlinear Klein-Gordon equation with a quadratic nonlinear term [2] have been solved using Adomian decomposition method.A system of nonlinear fractional partial differential equations has been solved using homotopy analysis method by Jafari and Seifi [3] and Bataineh et al. [5], using variational iteration method by Wazwaz [6].In [7], Wang and Cheng adopted variational method and finite element approach to solve damped nonlinear Klein-Gordon equations.
Keskin and Oturanc ¸ [23] have developed reduced differential transform method to solve partial differential equations of integer order [24] as well as fractional order.After Keskin and Oturanc ¸, RDT method has been implemented for the numerical computation of various physical models of engineering and sciences [25][26][27].

International Journal of Differential Equations
The main goal of this paper is to provide an analytical solution of initial value system of time dependent partial differential equations obtained by employing RDT method developed by Keskin and Oturanc ¸ [23].

Reduced Differential Transform Method
The basic properties of the fractional reduced differential transform method are described in this section.Let (, ) be a function of two variables such that (, ) = ()().By using the properties of the one-dimensional differential transform (DT) method (, ) can be written as where Ψ(, ) is referred to as the spectrum of (, ) and is defined by For more details on DT method, see [28] and the references therein.Denote the lowercase (, ) as the original function while its fractional reduced transformed function is denoted by the uppercase Ψ  ().Definition 1.If (, ) is analytic and continuously differentiable with respect to  and , then RDT of  is given by The reduced inverse differential transform of   () is defined as follows: Equations ( 3) and (4) together reduce to The basic properties of RDT method are found in [1,4] and can be deduced from ( 3) and ( 4), given in the following.

Some Basic Properties and Notation of RDT Method.
In this section, the properties of RDT method as in [23][24][25] have been revisited to complete our study.

Results and Discussion
In this section, we give five test problems of linear and nonlinear partial differential equations (PDEs) using reduced differential transform (RDT) method.
Example 2. Consider the initial value system of linear PDEs: On using RDT method (8) reduces to a set of recurrence relations as follows: On solving system (9), we get Using inverse RDT method (4), we get The same solution is obtained by using homotopy analysis method [3], variational iteration method [6], and homotopy perturbation method [21].The solution behavior of , V is depicted in Figure 1.
Example 3. Consider the following initial value system of nonlinear PDEs: On using RDT method (12) reduces to a set of recurrence relations as follows: On solving system (13), we have International Journal of Differential Equations  12) in (0, 1.5).
Example 4. Consider the following form of IVS of twodimensional coupled viscous Burgers' equation: On using RDT method (16) reduces to a set of recurrence relations as follows: On solving system (17), we get    16) in (−, ) at different time levels.
On using inverse RDT method (4), we get  (, ) =  − sin () , V (, ) =  − sin () .(19) This is the required exact solution of the initial value system of coupled viscous Burgers' equation (16).The same solution is obtained by homotopy perturbation method [21] and variational iteration method [22] for integer-order time derivatives.The physical behavior of , V is depicted in Figure 4.
Example 5. Consider the following form of initial value system of two-dimensional coupled viscous Burgers' equation: On using RDT method (20) reduces to a set of recurrence relations as follows: Solving the recurrence relation ( 21), we get . . .
Using inverse RDT method, the approximate solution is given as This is the required approximate solution of system (20).The approximate solution of system ( 20) is obtained by homotopy perturbation method [21] and variational iteration method [22] for integer-order time derivatives.
Example 6.Consider initial value system as follows: Setting  =  + V,   ( This is the desired approximate solution of the initial value system (24), which is the same as that obtained in [28] using DTM.

Concluding Remark
In this paper, reduced differential transform method has been implemented successfully to five test problems of the initial value systems of time dependent partial differential equations including two-dimensional coupled viscous Burgers' equations.The obtained results agreed well with homotopy perturbation method [21], homotopy analysis method [3], variational iteration method [6,22], and differential transform method [28].Easiness and effectiveness are the strength of RDT method.

Figure 4 :
Figure 4: The behavior of  and V of initial value system of coupled viscous Burgers' equation (16) in (−, ) at different time levels.