Acessibilidade / Reportar erro

Differential transformation method for solving one-space-dimensional telegraph equation

Abstract

In this research, the Differential Transformation Method (DTM) has been utilized to solve the hyperbolic Telegraph equation. This method can be used to obtain the exact solutions of this equation. In the end, some numerical tests are presented to demonstrate the effectiveness and efficiency of the proposed method. Mathematical subject classification: 35Lxx, 35Qxx.

differential transformation method; spectral method; telegraph equation; Taylor expansion


Differential transformation method for solving one-space-dimensional telegraph equation

B. Soltanalizadeh

Young Researchers Club, Sarab Branch, Islamic Azad University, Sarab, Iran. E-mail: babak.soltanalizadeh@gmail.com

ABSTRACT

In this research, the Differential Transformation Method (DTM) has been utilized to solve the hyperbolic Telegraph equation. This method can be used to obtain the exact solutions of this equation. In the end, some numerical tests are presented to demonstrate the effectiveness and efficiency of the proposed method.

Mathematical subject classification: 35Lxx, 35Qxx.

Key words: differential transformation method, spectral method, telegraph equation, Taylor expansion.

1 Introduction

In this paper we focus on the following one-space-dimensional Telegraphequation:

where α, β ∈ R, ψ : R × R R are known and u : R × RR is the unknown function. This equation with initial and boundary conditions is considered by some authors. Numerical solution of Telegraph equation with variable coefficient, was discussed in [16]. Mohanty and his coworkers [17, 18], also developed new three-level alternative direction implicit schemes for the two and three-space-dimensional linear hyperbolic equations that these are unconditionally stable. Saadatmandi and Dehghan [15], applied a shifted Chebyshev Tau method. Gao and Chi [22] solved Eq. (1) by a semi-discretization method. They transformed Eq. (1) into a system consisting of ordinary differential equations with respect to time t and then found an exact solution containing an infinite matrix series. Authors of [22], presented two unconditionally stable numerical schemes based on C3 quartic splines with accuracy orders of O(k5 + h4) and O(k7 + h4).

The concept of DTM was first proposed by Zhou [1], who solved linear and nonlinear problems in electrical circuit problems. Chen and Ho [2] developed this method for partial differential equations. Ayaz [3] applied it to the system of differential equations. During recent years, this method has been used for solving various types of equations by many authors. For example, this method has been used for differential algebraic equations [8], partial differential equations [5, 6, 7], fractional differential equations [10, 11], Volterra integral equations [24] and Difference equations [9]. Shahmorad et al. [19, 20] developed DTM to fractional-order integro-differential equations with nonlocal boundary conditions and a class of two-dimensional Volterra integral equations. Borhanifar and Abazari applied this method for Burgers and Schrödinger equations [4, 21, 23]. In [14], this method has been utilized for the Kuramoto-Sivashinsky equation with an initial condition. There exist similar problem. For example, authors of [12, 13, 14] presented several matrix formulation method for solving some equation with a boundary integral condition.

The aim of this paper is to extend the differential transformation method to solve the hyperbolic Telegraph equation. The method can be used to evaluate the approximating solution by the finite Taylor series and by an iteration procedure described by the transformed equations obtained from the original equation using the operations of differential transformation.

2 The definitions and operations of DT

2.1 The one-dimensional differential transform

The basic definitions and operations of one-dimensional DTM are introducedin [1, 2, 3] as follows:

Definition 2.1. If u(t) is analytic in the time domain T then

For t = ti, φ(t, k) = φ(ti, k), where k belong to the non-negative integer, denoted as the K domain. Therefore, Eq. (2) can be rewritten as

where Ui(k) is called the spectrum of u(t) at t = ti in the K domain.

Definition 2.2.If u(t) is analytic, then it can be shown as

Equation (4) is known as the inverse transformation of U(k). If U(k) isdefined as

then the function u(t) can be described as

where M(k) ≠ 0, q(t) ≠ 0. M(k) is called the weighting factor and q(t) is regarded as a kernel corresponding to u(t). If M(k) = 1 and q(t) = 1 then Eqs. (4) and (6) are equivalent. In this paper, the transformation with M(k) = 1/k! and q(t) = 1 is applied. Thus from Eq. (7), we have

Using the differential transform, a differential equation in the domain of interest can be transformed to be an algebraic equation in the K domain and u(t) can be obtained by finite-term Taylor series plus a remainder, as

In order to speed up the convergence rate and improve the accuracy of calculation, the entire domain of t needs to be split into sub-domains [9, 10].

2.2 The two-dimensional differential transform

Consider a function of two variables w(x, t): R × R R, and supposethat it can be represented as a product of two single-variable functions, i.e., w(x, t) = u(x)v(t). Based on the properties of one-dimensional differential transform, function w(x, t) can be represented as

where W(i, j) is called the spectrum of w(x, t). Now we introduce the basic definitions and operations of two-dimensional DT as follows [10].

Definition 2.3. If w(x, t) is analytic and continuously differentiable with respect to time t in the domain of interest, then

where the spectrum function W(k, h) is the transformed function, which is also called the T-function. Let w(x, t) as the original function while the uppercase W(k, h) stands for the transformed function. Now we define the differential inverse transform of W(k, h) as following:

Using Eq. (10) in (11), we have

where x0 = 0 and t0 = 0.

Now from the above definitions and Eqs. (11) and (12), we can obtain someof the fundamental mathematical operations performed by two-dimensionaldifferential transform in Table 1.

3 Application of the DTM

In this section, we apply the DTM to solve the presented Telegraph equation. Consider the equation (1),

with the initial conditions

and

and boundary conditions

and

Let U(k, h) as the differential transform of u(x, t). Applying Table 1, Eq. (2) and Definition 2.3 when x0 = t0 = 0, we get the differential transform version of Eq. (13) as following:

By the first initial condition we get

which implies

and from the second initial condition and Table 1, we have

which implies

Similarly, from the boundary conditions, we have

Therefore, the values of U(k, 0), U(k, 1), and for U(0, h) can be obtained from Eqs. (20), (22) and (23). By using Eqs. (18) and (24), the remainder values of U can be found as follows:

Example 3.1. Consider the Eqs. (13)-(15) with the

Applying Eqs. (20), (22) and (23) in initial and boundary conditions of this problem, we get

and

From Eqs. (24) and (25), we have

By continuing this process, we obtain

which implies that

which is the Taylor expansion of the

u(x, t) = x × exp(–t) 0 < x < L,

which is the exact solution of the Example 3.1. The computational results of Example 3.1 are presented in Table 3, and the plot of corresponding exact and approximate functions are shown in Figs. 1 and 2.



Example 3.2 ([11]) Consider the Eqs. (13)-(15) with the

From the above initial and boundary conditions and Eqs. (20), (22) and (23), we obtain the corresponding spectra as follows:

and

and

From Eqs. (24) and (25), we have

Thus we get

which implies that

which is the Taylor expansion of the

u(x, t) = exp(–t) sin(x),

which is the exact solution of the Example 3.2. The computational results of Example 4.2 are presented in Table 5, and the plot of corresponding exact and approximate functions are shown in Figs. 3 and 4.



4 Conclusions

In this paper, we solved the Telegraph equation with initial conditions by theDifferential Transformation method. By using this method, Numerical/analytical results obtained by a simple iterative process. The numerical results prove that this method is a powerful techniques for this case of problems. Consequently, it is seen that this method can be an alternative way for the solution of partial differential equations that have no analytic solutions.

Acknowledgments. The author would like to thank the reviewers for carefully reading the paper and for their constructive comments and suggestions that have improved the paper. The author is deeply grateful to the Young Researchers Club, Sarab Branch, Islamic Azad University, for the financial supports.

Received: 01/XI/10.

Accepted: 02/IV/11.

#CAM-281/10.

  • [1] J.K. Zhou, Differential Transformation and Its Applications for Electrical Circuits. Huazhong University Press, Wuhan, China, 1986 (in Chinese).
  • [2] C.K. Chen and S.H. Ho, Solving partial differential equations by two dimensional differential transform method. Appl. Math. Comput., 106 (1999), 171-179.
  • [3] F. Ayaz, Solutions of the systems of differential equations by differential transform method. Appl. Math. Comput., 147 (2004), 547-567.
  • [4] R. Abazari and A. Borhanifar, Numerical study of the solution of the Burgers and coupled Burgers equations by a differential transformation method. Comput. Math. Appl., 59 (2010), 2711-2722.
  • [5] C.K. Chen, Solving partial differantial equations by two dimensional differential transformation method. Appl. Math. Comput., 106 (1999), 171-179.
  • [6] M.J. Jang and C.K. Chen, Two-dimensional differential transformation method for partial differantial equations. Appl. Math. Comput., 121 (2001), 261-270.
  • [7] F. Kangalgil and F. Ayaz, Solitary wave solutions for the KDV and mKDV equations by differential transformation method. Choas Solitons Fractals, 41 (2009), 464-472.
  • [8] J.D. Cole, On a quasilinear parabolic equation occurring in aerodynamics. Quart. Appl. Math., 9 (1951), 225-236.
  • [9] A. Arikoglu and I. Ozkol, Solution of difference equations by using differential transformation method. Appl. Math. Comput., 174 (2006), 1216-1228.
  • [10] S. Momani, Z. Odibat and I. Hashim, Algorithms for nonlinear fractional partial differantial equations: A selection of numerical methods. Topol. Method Nonlinear Anal., 31 (2008), 211-226.
  • [11] A. Arikoglu and I. Ozkol, Solution of fractional differential equations by using differential transformation method. Chaos Solitons Fractals, 34 (2007), 1473-1481.
  • [12] H.R. Ghehsareh, B. Soltanalizadeh and S. Abbasbandy, A matrix formulation to the wave equation with nonlocal boundary condition. Inter. J. Comput. Math., 88 (2011), 1681-1696.
  • [13] B. Soltanalizadeh, Numerical analysis of the one-dimensional Heat equation subject to a boundary integral specification. Optics Communications, 284 (2011), 2109-2112.
  • [14] B. Soltanalizadeh and M. Zarebnia, Numerical analysis of the linear and nonlinear Kuramoto-Sivashinsky equation by using Differential Transformation method. Inter. J. Appl. Math. Mechanics, 7(12) (2011), 63-72.
  • [15] A. Saadatmandi and M. Dehghan, Numerical solution of hyperbolic Telegraph equation using the chebyshev Tau method. Numerical Methods for Partial Differential Equations. doi 10.1002/num.
  • [16] R. Aloy, M.C. Casaban, L.A. Caudillo-Mata and L. Jodar, Computing the variable coefficient telegraph equation using a discrete eigenfunctions method. Comput. Math. Appl., 54 (2007), 448-458.
  • [17] R.K. Mohanty and M.K. Jain, An unconditionally stable alternating direction implicit scheme for the two space dimensional linear hyperbolic equation. Numer Methods Partial Differential Equations, 17 (2001), 684-688.
  • [18] R.K. Mohanty, M.K. Jain and U. Arora, An unconditionally stable ADI method for the linear hyperbolic equation in three space dimensions. Int. J. Comput. Math., 79 (2002), 133-142.
  • [19] A. Tari, M.Y. Rahimi, S. Shahmoradb and F. Talati, Solving a class of two-dimensional linear and nonlinear Volterra integral equations by the differential transform method. J. Comput. Appl. Math., 228 (2009), 70-76.
  • [20] D. Nazari and S. Shahmorad, Application of the fractional differential transform method to fractional-order integro-differential equations with nonlocal boundary conditions. J. Comput. Appl. Math., 234 (2010), 883-891.
  • [21] A. Borhanifar and R. Abazari, Exact solutions for non-linear Schrödinger equations by differential transformation method. J. Appl. Math. Comput., 35 (2011), 37-51.
  • [22] F. Gao and C.M. Chi, Unconditionally stable difference schemes for a one-space-dimensional linear hyperbolic equation. Appl. Math. Comput., 187 (2007), 1272-1276.
  • [23] A. Borhanifar and R. Abazari, Numerical study of nonlinear Schrödinger and coupled Schrödinger equations by differential transformation method. Optics Communications, 283 (2010), 2026-2031.
  • [24] Z.M. Odibat, Differential transformation method for solving Volterra integral equations with separable kernels. Math. Comput. Modeling, 48 (2008), 1144-1149.

Publication Dates

  • Publication in this collection
    06 Jan 2012
  • Date of issue
    2011

History

  • Accepted
    02 Apr 2011
  • Received
    01 Nov 2010
Sociedade Brasileira de Matemática Aplicada e Computacional Sociedade Brasileira de Matemática Aplicada e Computacional - SBMAC, Rua Maestro João Seppe, nº. 900 , 16º. andar - Sala 163, 13561-120 São Carlos - SP Brasil, Tel./Fax: 55 16 3412-9752 - São Carlos - SP - Brazil
E-mail: sbmac@sbmac.org.br