Exact Solution for a Class of Stiff Systems by Differential Transform Method

In this paper, Differential Transform Method (DTM) is proposed for the closed form solution of linear and non-linear stiff systems. First, we apply DTM to find the series solution which can be easily converted into exact solution. The method is described and illustrated with different examples and figures are plotted accordingly. The obtained result confirm that DTM is very easy, effective and convenient.

• Differential Transform Method; Stiff System; Inverse Differential Transform

1. Introduction

Consider the stiff initial value problem [1]:

on the finite interval, where

and are continuous.

The initial value problem of stiff differential equations occurs in almost every field of science [1-5], particularly, in the fields of:

1) Chemical Reactions: A famous chemical reaction is the Oregenator reaction between HBrO₂, Br⁻, and Ce (IV) described by Field and Noyes in 1984.

2) Reaction-diffusion systems: Problems in which the diffusion is modeled via the Laplace operator may become stiff as they are discretized in space by finite differences or finite elements. Well-known example of such systems which appear so often in mathematical biology.

Several further occurrences of stiffness can be found in electrical circuits, mechanics, meteorology, oceanography and vibrations.

Definition 1: If the solution of the system contains components which change at significantly different rates for given changes in the independent variable, then system is said to be stiff [2,3].

Stiff differential equations are characterized as those whose exact solution has a term of the form, where

is a large positive constant. The key features of stiff equations are that the derivative terms may increase rapidly as t increases [1].

In the last three decades numerous works have been focusing on the development of more advanced and efficient methods for stiff problems [1,2]. The situation becomes more complicated when stiffness coupled with nonlinearity. Carroll presents an exponential fitted scheme for solving stiff systems of initial value problems [3]. The numerical solution of linear and nonlinear system of stiff system can be found in [4-6].

Differential Transform Method (DTM) is a semi numerical method which gives series solution. But sometimes the series solution can be easily converted into closed form solution. This Method was first introduced by Zhou [7] who solved linear and non-linear problems in electrical circuit problems. After that the method become very popular to solved different kind of problems. The year wise detail of such problems is as follows.

In [8], two dimensional DTM is used to solve partial differential equations (PDEs). In [9], one dimensional DTM is applied to solve linear and nonlinear initial value problems. In [10], one and two dimensional DTM has been applied to solve eigenvalue problems and PDEs. In [11], two dimensional DTM is applied to solve initial value problems for PDEs. In [12], DTM is used to solve parameters identification problems. In [13], DTM was applied to transient advective-dispersive transport equation. In [14] used one dimensional DTM to solve differential-algebraic equations. In [15], three dimensional DTM is used to solve linear and nonlinear PDEs. In [16], generalized DTM is used to solve PDEs. In [17,18], DTM is used to solve boundary value problems for integro-differential equations. In [19], different approximate methods has been used for initial value problems. In [20], two dimensional DTM is utilized to solve linear and nonlinear PDEs. In [21], DTM is used to solve free vibrations equations of beam on elastic soil. In [22], DTM is applied to solve difference equations. In [23], DTM is applied to solve differential-difference equations. In [24, 25], Generalized DTM is used to solve multi order fractional differential equations and linear PDEs of fractional order respectively. In [26], DTM is used to solve fourth order boundary value problem. In [27], two dimensional DTM is used to solve non-linear oscillatory system. In [28], Idrees et al. has applied DTM for the exact solution of Goursat Problems, in [29] the authors have utilized an efficient method for stiff system, whilst in [30] the authors presented the numerical solution of the stiff system.

In this paper, we solve the linear and non-linear stiff system via DTM. In Section 2, we give some basic properties of one-dimensional DTM. In Section 3, we have applied the method to linear and non-linear stiff systems.

2. One-Dimensional Differential Transform

In this section, we first give some basic properties of one-dimensional differential transform method. Differential transform of a function is defined as follows:

where is the original function and is the transformed function for. The differential inverse transform of is defined as

From Equations (2) and (3) we get

which implies that the concept of DTM is derived from Taylor series expansion, but the method does not evaluate the derivative symbolically. However, relative derivatives are calculated by an iterative procedure which is described by the transformed equations of the original functions.

In this work, we use the lower case letters to represent the original functions and upper case letters to represent the transformed functions. In actual applications, the function is expressed by a finite series Equation (5) can be written as

Here is represented the convergence of natural frequency.

From Equations (2) and (3) we obtain Table 1 of the fundamental operations of one-dimensional differential transform method is given by.

3. Application to Stiff System

In this section, we apply DTM to both linear and nonlinear stiff systems.

Problem 1: Consider the linear stiff system:

with initial value,

This system has eigen values of large modulus lying closed to the imaginary axis.

By applying Differential Transformation, we have

The initial conditions of Differential Transformation are given by:

.

For, the series coefficients for and can be obtained as

Table 1. The fundamental operations of one-dimensional DTM.

We used MATHEMATICA to calculate the unknown coefficients and.

Using the inverse Transform, we get

Equations (10) and (11) can be written in the exponential form are given by

and

Thus we get the exact solution by differential transform method.

Problem 2: Consider the non-linear system in the form of initial value problems [1] is given by:

Applying Differential Transform, we have

For the series coefficients for and can be obtained as

We used MATHEMATICA to calculate the unknown coefficients and.

Using the inverse Transform:

This can be written as folows

and

Thus we get the exact solution by differential transform method.

Problem 3: Consider the system of initial value problems [29]:

with initial conditions, i = 1, 2, 3, 4.

Applying Differential Transform, we have

;

;

;

.

With the transformed initial conditions are,.

For the series coefficients for, , and can be obtained as:

Using the inverse Transform:

We obtain

, , ,.

In this section, we have presented three different linear and nonlinear stiff systems via Differential Transform Method and the series solution of Equations (13) and (22) have shown in Figures 1 and 2 respectively.

4. Conclusion

In this work, DTM has been applied to the exact solution of linear and non-linear stiff system. DTM is a semi numerical method and we obtained a closed form solution such as [11] and expressed the series form solutions in graphical form for the first two examples. It has been

observed that DTM is simpler, effective and reliable.

References

1. G. Hojjati, M. Y. Rahimi Arabili and S. M. Hosseini, “An Adaptive Method for Numerical Solution of Stiff System of Ordinary Differential Equations,” Mathematics and Computers in Simulation, Vol. 66, No. 1, 2004, pp. 33-41. doi:10.1016/j.matcom.2004.02.019

2. E. Harier and G. Wanner, “Solving Ordinary Differential Equations II, Stiff and Differential-Algebraic Problems,” Springer-Verlag, New York, 1996. doi:10.1007/978-3-642-05221-7

3. L. Lopidus and W. E. Schiesser, “Numerical Methods for Differential Systems,” Academic Press, New York, 1976.

4. J. Carroll, “A Metrically Exponentially Fitted Scheme for the Numerical Solution of Stiff Initial Value Problems,” Computers & Mathematics with Applications, Vol. 26, 1993, pp. 57-64. doi:10.1016/0898-1221(93)90034-S

5. C. H. Hsiao, “Haar Wavelet Approach to Linear Stiff Systems,” Mathematics and Computer in Simulation, Vol. 64, No. 1, 2004, pp. 561-567. doi:10.1016/j.matcom.2003.11.011

6. C. H. Hsiao and W. J. Wang, “Haar Wavelet Approach to Non-Linear Stiff Systems,” Mathematics and Computer in Simulation, Vol. 57, No. 6, 2001, pp. 347-353. doi:10.1016/S0378-4754(01)00275-0

7. J. K. Zhou, “Differential Transformation and Its Applications for Electrical Circuits,” Huazhong University Press, Wuhan, 1986 (in Chinese).

8. C. K. Chen and S. H. Ho, “Solving Partial Differential Equations by Two-Dimensional Differential Transform Method,” Applied Mathematics and Computation, Vol. 106, No. 2-3, 1999, pp. 171-179. doi:10.1016/S0096-3003(98)10115-7

9. M. J. Jang, C. L. Chen and Y. C. Liu, “Two-Dimensional Differential Transform for Partial Differential Equations,” Applied Mathematics and Computation, Vol. 121, No. 2-3, 2001, pp. 261-270. doi:10.1016/S0096-3003(99)00293-3

10. I. H. Abdel-Halim Hassan, “Different Applications for the Differential Transformation in the Differential Equations,” Applied Mathematics and Computation, Vol. 129, No. 2-3, 2002, pp. 183-201. doi:10.1016/S0096-3003(01)00037-6

11. F. Ayaz, “On the Two-Dimensional Differential Transform Method,” Applied Mathematics and Computation, Vol. 143, No. 2-3, 2003, pp. 361-374.

12. M. J. Jang, J. S. Wang and Y. C. Liu, “Applying Differential Transformation Method to Parameter Identification Problems,” Applied Mathematics and Computation, Vol. 139, No. 2, 2003, pp. 491-502. doi:10.1016/S0096-3003(02)00211-4

13. C. K. Chen and S. P. Ju, “Application of Differential Transformation Method to Transient Adjective-Disper-sive Transport Equations,” Applied Mathematics and Computation, Vol. 155, No. 1, 2004, pp. 25-38. doi:10.1016/S0096-3003(03)00755-0

14. F. Ayaz, “Application of Differential Transform Method to Differential-Algebraic Equations,” Applied Mathematics and Computation, Vol. 152, No. 3, 2004, pp. 649-657. doi:10.1016/S0096-3003(03)00581-2

15. F. Ayaz, “Solution of the System of Differential Equations by Differential Transform Method,” Applied Mathematics and Computation, Vol. 147, No. 2, 2004, pp. 547-567. doi:10.1016/S0096-3003(02)00794-4

16. A. Kurnaz, G. Oturanc and M. E. Kiris, “N-Dimensional Differential Transform Method for Solving PDEs,” International Journal of Computer Mathematics, Vol. 82, No. 3, 2005, pp. 369-380. doi:10.1080/0020716042000301725

17. A. Arikoglu and I. Ozkol, “Solution of Boundary Value Problems for Integro Differential Equation by Using Differential Transform Method,” Applied Mathematics and Computation, Vol. 168, No. 2, 2005, pp. 1145-1158. doi:10.1016/j.amc.2004.10.009

18. A. Arikoglu and I. Ozkol, “Solution of Fractional Differential Equation by Using Differential Transform,” Chaos Solitons & Fractials, Vol. 34, No. 5, 2006, pp. 1473-1381. doi:10.1016/j.chaos.2006.09.004

19. I. H. Abdel-Halim Hassan, “Comparison of Differential Transformation Technique with Adomian Decomposition Method for Linear and Nonlinear Initial Value Problems,” Communications in Nonlinear Science and Numerical Simulation, Vol. 36, No. 1, 2008, pp. 53-65.

20. N. Bildik, A. Konuralp, F. O. Bek and S. Kucukarslan, “Solution of Different Type of the Partial Differential Equations by Differential Transform Method and Adomian Decomposition Method,” Applied Mathematics and Computation, Vol. 172, No. 1, 2006, pp. 551-567. doi:10.1016/j.amc.2005.02.037

21. S. Catal, “Solution of Free Vibration Equations of Beam on Elastic Soil by Using Differential Transform Method,” Applied Mathematical Modelling, Vol. 32, No. 9, 2008, pp. 1744-1757. doi:10.1016/j.apm.2007.06.010

22. A. Arikoglu and I. Ozkol, “Solution of Difference Equation by Using Differential Transform Method,” Applied Mathematics and Computation, Vol. 174, No. 2, 2006, pp. 1216-1228. doi:10.1016/j.amc.2005.06.013

23. A. Arikoglu and I. Ozkol, “Solution of Differential-Difference Equation by Using Differential Transform Method,” Applied Mathematics and Computation, Vol. 181, No. 1, 2006, pp. 153-162. doi:10.1016/j.amc.2006.01.022

24. A. Arikoglu and I. Ozkol, “Comments on Application of Generalized Differential Transform Method to Multi-Order Fractional Differential Equation,” Communication in Nonlinear Science and Numerical Simulation, Vol. 13, No. 8, 2008, pp. 1337-1340.

25. Z. Odibat and S. Momani, “A Generalized Differential Transform Method for Linear Partial Differential Equations of Fractional Order,” Applied Mathematics Letters, Vol. 21, No. 2, 2007, pp. 194-199.

26. V. S. Erturk and S. Momani, “Comparing Numerical Methods for Solving Fourth-Order Boundary Value Problems,” Applied Mathematics and Computation, Vol. 188, No. 2, 2007, pp. 1963-1968. doi:10.1016/j.amc.2006.11.075

27. M. El-Shahed, “Application of Differential Transform Method to Non-Linear Oscillatory Systems,” Communications in Nonlinear Science and Numerical Simulation, Vol. 13, No. 8, 2007, pp. 1714-1720.

28. M. Idrees, S. Islam, R. A. Shah and M. Zeb, “Exact Solution of Goursat Problems Using Differential Transform Method,” Journal of Advanced Research in Scientific Computing, Vol. 3, No. 3, 2011, pp. 1-13.

29. G. A. F. Ismail and I. H. Ibrahim, “New Efficient Second Derivative Multistep Methods for Stiff Systems,” Applied Mathematical Modelling, Vol. 23, No. 4, 1999, pp. 279-288. doi:10.1016/S0307-904X(98)10086-0

30. N. Guzel and M. Bayram, “On the Numerical Solution of Stiff Systems,” Applied Mathematics and Computation, Vol. 170, No. 1, 2005, pp. 230-236. doi:10.1016/j.amc.2004.11.035