Discrete Dynamics in Nature and Society

Volume 2011 (2011), Article ID 561015, 17 pages

http://dx.doi.org/10.1155/2011/561015

## An Operator-Difference Method for Telegraph Equations Arising in Transmission Lines

Department of Elementary Mathematics Education, Mevlana University, 42003 Konya, Turkey

Received 12 June 2011; Accepted 3 August 2011

Academic Editor: Hassan A. El-Morshedy

Copyright © 2011 Mehmet Emir Koksal. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

A second-order linear hyperbolic equation with time-derivative term subject to appropriate initial and Dirichlet boundary conditions is considered. Second-order unconditionally absolutely stable difference scheme in (Ashyralyev et al. 2011) generated by integer powers of space operator is modified for the equation. This difference scheme is unconditionally absolutely stable. Stability estimates for the solution of the difference scheme are presented. Various numerical examples are tested for showing the usefulness of the difference scheme. Numerical solutions of the examples are provided using modified unconditionally absolutely stable second-order operator-difference scheme. Finally, the obtained results are discussed by comparing with other existing numerical solutions. The modified difference scheme is applied to analyze a real engineering problem related with a lossy power transmission line.

#### 1. Introduction

Second-order linear hyperbolic partial differential equations with both constant and variable coefficients arise in many branches of science and engineering, for example, electromagnetic, electrodynamics, thermodynamics, hydrodynamics, elasticity, fluid dynamics, wave propagation, and materials science, see [1–13]. For example, they are used frequently for modelling power transmission lines [7–10, 13]. In numerical methods for solving these equations, the problem of stability has received a great deal of importance and attention (see [14–20]). Specially, a suitable model for analyzing the stability is provided by a proper unconditionally absolutely stable difference scheme with an unbounded operator.

In the literature, the work on new finite difference schemes has drawn important attention for numerical solutions of linear hyperbolic partial differential equations, (see [21–26] and the references given therein). However, these difference schemes are conditionally stable since the stability estimates of them are based on some restrictions on choice of the grid step sizes and with respect to the time and space variables, respectively.

Many scientists have investigated unconditionally stable difference schemes for linear hyperbolic differential equations. Such a difference scheme for approximately solving linear hyperbolic differential equations was studied for the first time in [14]. The first-order difference scheme was constructed for approximately solving the abstract Cauchy problem where is an unbounded selfadjoint positive linear operator with domain in an arbitrary Hilbert space . It is known (see [27, 28]) that various initial-boundary value problems for hyperbolic equations can be reduced to the initial value problem (1.1). Note that (1.1) is the well-known wave equation in the special case when is equal to Laplace operator . The stability estimates for the solution of this difference scheme and for the first- and second-order difference derivatives were established in Hilbert space.

Then, in the past decade, a huge variety of works on finite difference method for numerical solutions of linear hyperbolic partial differential equations were studied. For the problem (1.1) when , Ashyralyev and Sobolevskii [15] developed the first- and two different types of second-order difference schemes. The stability estimates for the solutions of these difference schemes and for the first- and second-order difference derivatives have been established. For the same problem, they developed also the high-order two-step difference methods generated by an exact difference scheme, and by the Taylor expansion on three points in [16]; here, the stability estimates for approximate solutions by these difference methods are also discussed. In [17, 18], two different types of second-order difference methods for the problem (1.1) were studied, and the stability estimates for the solutions of these difference schemes were established. But, the difference methods developed in these references are generated by square roots of . Further, an operator of which construction is not easy is required to realize these difference methods practically. Therefore, in spite of derived theoretical results, it is not very practical to apply these methods for solving an initial-value problem numerically. Another two different types of second-order difference methods generated by integer powers of for the problem (1.1) were developed and the stability estimates for the solutions of these difference methods were established in [19, 20]. We should mention that all difference methods in references [15–20] are unconditionally absolutely stable and applicable for multidimensional linear hyperbolic equations with space variable coefficient. Specially, the difference methods in [17–20] are also applicable for multidimensional linear hyperbolic equations with time and space variable coefficients.

For approximately solving the problem when and , a three-level implicit difference scheme of was discussed in [29] wheres and are real numbers. For , , and , (1.2) is called as a damped wave equation and as a telegraph equation, respectively, see [22]. In [30], also a three-level implicit difference scheme, whose order is the same as the previous literature for approximately solving the equation defined in the region , where with the same initial and boundary conditions of (1.2) at different points of the space variable, was developed. Gao and Chi. [31] developed two explicit difference schemes to solve the problem (1.2) numerically. The accuracy orders of these difference schemes are and , respectively. Recently, H.-W. Liu and L.-B. Liu. [32] developed a new implicit difference scheme based on quartic spline interpolation in space direction and finite discretization in time direction for solving the same problem. The accuracy order of this difference scheme is second order in time direction and fourth order in space direction. In [33], a difference scheme based on alternating direction implicit methods was presented for solving multidimensional telegraph equations with different coefficients numerically. The order of accuracy of the proposed scheme, in which two free parameters are introduced, is also . All the difference schemes in references [29–33] are unconditionally stable.

Note that Dehghan and Shokri [34] proposed a new numerical scheme based on radial-based function method (Kansa's method) for solving the equation with the same initial and Dirichlet boundary conditions in (1.2), where and are known constants.

Although there are other important unconditionally stable difference schemes in the literature [35, 36] to solve the equation in (1.2) numerically, it is not possible to consider all of them in detail in a single paper.

The only alternative is the application of stable numerical methods for solving the initial-boundary value problem (1.2) since the analytical solution cannot be determined for arbitrary . Being easy to implement and universally applicable, the finite difference method outstands among many of the available numerical methods. For solving a linear hyperbolic partial differential equation with time or space variable coefficient analytically, there is no specific method in the literature.

In the present paper, the second order in both time and space direction unconditionally absolutely stable operator-difference scheme in [20] generated by integer powers of is modified for approximately solving the equation Here, is an unbounded selfadjoint positive linear operator with domain in an arbitrary Hilbert space . The initial conditions of the above equation consists of Numerical solutions of the equations (1.2), (1.3), and (1.4) are computed with the modified difference scheme. The results are provided and compared with the difference schemes in [29–34]. The differences among all these difference schemes are illustrated with numerical examples by considering various initial-boundary value problems.

#### 2. Operator Difference Schemes-Stability Estimates

Using the finite difference formulas in the equation we obtain Further, we have Neglecting small terms , we obtain the following second-order difference scheme: for approximately solving (1.5).

Theorem 2.1. *For the solution of the difference scheme (2.5), the stability estimate
**
holds, where and does not depend on , and .*

Theorem 2.2. *For the solution of the difference scheme (2.5), the stability estimate
**
holds, where and does not depend on , and .*

*Proof. *Theorems of the existence and uniqueness of the classical solutions of the initial value problem
for the hyperbolic equation in a Hilbert space with the selfadjoint positive operator with a region of definition that does not depend on under some additional condition of differentiability on the nonlinear maps were established in [37].

The proofs of the above Theorems are based on the symmetry properties of the operators defined by formula (2.5) and the following theorem on the coercivity inequality for the solution of the elliptic difference problem in .

Theorem 2.3. *For the solutions of the elliptic difference problem
**
the following coercivity inequality holds (Sobolevskii) [38]*

#### 3. Numerical Examples and Discussion

We have not been able to obtain a sharp estimate for the constants figuring in the stability inequalities. So, in this section, the numerical results are presented by considering various initial-boundary value problems. The executions in all examples are carried out by MATLAB 7.01 and obtained by a PC Pentium (R) 2CPV, 2.00 6 Hz, 2.87 GB of RAM. The errors in the numerical solutions are computed by the root mean square (RMS) error [39] and by the absolute error, where and represent the exact and numerical solutions at final time, respectively.

*Example 1. *First, the following initial-boundary value problem [29]:
with the exact solution is considered. According to the problem (3.3), the differential operator can be defined follows:
Concerning the above differential operator, discrete operators can be defined as

For the approximate solution of the problem (3.3) and the ones appearing in the sequel, using the second order of approximation of second and fourth derivatives in the above discrete operators and applying the second-order difference scheme (2.5), the system of linear equations is obtained. Then, writing the system in the matrix form, a second-order difference equation with respect to with matrix coefficients is arrived. To solve this resulting difference equation, iterative method is applied (see the third Chapter of [40]).

The exact and numerical solution obtained for and by using the second-order difference scheme (2.5) with and is shown in Figures 1 and 2, respectively, as an example. The exact and numerical solutions of the next examples are shown in odd and even figure number, respectively.

The difference between Figures 1 and 2 is not clearly obvious though the largest step-size values which correspond to minimum grid numbers, , , are used to graph them. The errors should be computed for the accurate comparison of the numerical and exact solutions as well as for the comparison of the two different difference schemes. RMS errors at for different values of and are tabulated in Tables 1 and 2, and the results are compared with the solutions obtained by [29].

Though both difference schemes are of , considering the errors in the numerical results tabulated in Tables 1 and 2 for the different values of the parameters , and the step sizes ,, it can be seen easily in Table 1 that the present scheme is approximately 2 times better in the sense of accuracy than the scheme in [29] for , . For the values of parameters indicated in Table 2, the present scheme is a little bit better than the scheme in [29] for , but it yields approximately 1.5 times smaller error for .

Although both of the schemes compared in this example are of , the scheme in [29] requires adequate choice of its free parameters. Otherwise, the results come out to be worse than the results of this paper.

*Example 2. *When the following equation [31, 32, 34]:
with initial conditions
and boundary conditions
is considered, then the analytical solution of this problem is . We apply the same procedure (3.3) for the approximate solution of the above problem.

The exact solution and the numerical solution obtained for and with the largest step sizes and are shown in Figures 3 and 4.

The difference between the exact solution and numerical solutions for the given example is not obvious though the largest step-size values are used to graph them. For smaller values of and , there is apparently no difference between the graphs of the exact and numerical solutions without depending on the values of the parameters and . So again, the errors should be computed for the accurate comparison. The absolute errors involved by the scheme (2.5) and by the scheme in [31] at , are listed in Tables 3, 4 and 5 for different values of and .

All the results obtained by the scheme (2.5) are almost the same with those obtained by using the scheme in [31] except the results listed in the first rows of Tables 4 and 5. Although the order of accuracy of the scheme in [31] is much better than of the present paper, this is not apparent in the results of Tables 3, 4 and 5. In fact, its effect shows up when larger enough values (e.g., ) are used [31, 32].

The absolute errors given by the scheme (2.5) and by the scheme in [32] in numerical solutions of the same example are listed in Tables 6 and 7 for and , respectively, with the same step size for time variable.

For all values of , the errors at of the scheme in [32] are much smaller than those of the present scheme (2.5). This is naturally expected since the order of accuracy of the scheme in [32] is which is better than that of (2.5), namely, . This difference in the errors disappears at on behalf of the scheme (2.5). In other words, as increases to 2.0, the errors of the scheme (2.5) are decreasing quietly for all values of , but the errors of the scheme in [32] remain almost at the same level when compared with those at .

Now, consider the main equation in (1.2) in two cases , and , with the same exact solution and initial-boundary conditions of the problem (3.7). The errors are tabulated in Table 8 for and and compared with those found by using the radial basic function method in [34].

It is obvious that the scheme (2.5) is almost 5 and 10 times better in the sense of error size than the scheme in [34] for , and , respectively. However, the computation time of the solution in [34] is quite smaller than that of the present method.

*Example 3. *Third, consider the following equation [32, 33]:
over a region . The initial and boundary conditions are given by
The analytical solution given in [32, 33] is . The exact and numerical solutions obtained for and with and are shown in Figures 5 and 6.

The difference between the exact and numerical solutions for the given example is not obviously apparent though the largest step-size values are used again to graph them. Numerical solutions for different values of and when , for , and , obtained by using the scheme (2.5) are compared with those obtained by using the schemes in [32, 33]. The RMS errors in solutions are tabulated in Tables 9, 10 and 11 at .

The errors in Tables 9, 10 and 11 involved in the results obtained by using the schemes of [32, 33] are naturally much smaller than those in Table 11. The reason of this huge difference is based on that the order of accuracy of the schemes in [32, 33] is , whilst it is in the present paper.

*Example 4. *Finally, the following equation [30]:
with time and space variable coefficients that are defined in the region is solved numerically. The initial and boundary conditions of this problem are the same as the initial and boundary conditions of the first example because they have the same exact solution. RMS errors are tabulated at the points and for a fixed values , and the results are listed in Table 12.

Apparently, for the given equation with both time and space variable coefficients, the results of the scheme (2.5) are much more satisfactory than the results of the scheme in [30] for all values of shown in Table 12 although the order of accuracy for both schemes is . On the other hand, when applied to solve an equation with constant coefficients, both schemes result with the same order of errors.

#### 4. Application

As an example, consider the power transmission line having a length of with the resistance, inductance, conductance, and capacitance parameters ,, , and in [8].

The line is represented by the following telegraphic equation: When the line parameters are normalized with respect to the impedance frequency , and distance , then the above second-order partial differential equation becomes Due to the slight losses of the line the coefficients of and in the above equation are quite small. When the relaxed line is excited at one end (sending end) by a unit pulse of duration which corresponds to a normalized time of , the voltage at the receiving end is computed by solving (4.2) with the boundary and initial conditions The problem is solved by the difference scheme (2.5) for the intervals for time and for space with number of grids 1000 and 100, respectively. The plot of the computed receiving-end voltage is shown in Figure 7(a) including the exciting voltage as well. It is obvious that the input pulse is reached at the end of the line at , and it is doubled in magnitude by the reflection at the open end. The wave travels back to the sending end, it gets inverted due to the short-circuit property of the sending end, and reflects back to the receiving end after inversion. Since the travelling time along the line is 1, the wave reaches the receiving end at as the second time. Due to the very slight loss in the line parameters, the attenuation is not obvious in Figure 7(a). The same problem is solved by increasing the line losses ( and parameters) by times, though this is not a realistic case. The attenuation in the receiving-end voltage is now obvious as plotted in Figure 7(b).

#### 5. Concluding Remarks

In this work, the second-order unconditionally absolutely stable difference scheme in [20] generated by integer powers of is modified for solving the abstract Cauchy problem for the second-order linear hyperbolic differential equations containing the unbounded selfadjoint positive linear operators. The modified difference scheme is applied for solving various initial-boundary value problems and compared with other published papers in the literature. It is observed that this scheme is more accurate as compared to the other difference schemes whose order of accuracy is the same. Moreover, the present difference scheme (2.5) is compared with difference schemes whose order of accuracy is higher, and it is observed that the good results can be obtained using the difference scheme in the paper depending on the type of the problem and solution parameters of the problem. The difference scheme is applied to a linear hyperbolic equation with variable coefficients, and it is observed that it is much more satisfactory than the scheme obtained by radial basic function, having the same order of accuracy. The modified difference scheme (2.5) is shown to be well applied to analyze a real engineering problem involving a lossy power transmission line with constant coefficients. Being defined by an operator, it can be applied to multidimensional linear hyperbolic differential equations with time and space variables coefficients.

#### Acknowledgment

The author is very thankful to Professor Ashyralyev for his valuable comments.

#### References

- H. Lamb,
*Hydrodynamics*, Cambridge Mathematical Library, Cambridge University Press, Cambridge, UK, 6th edition, 1993. - J. Lighthill,
*Waves in Fluids*, Cambridge University Press, Cambridge, UK, 1978. - J. A. Hudson,
*The Excitation and Propagation of Elastic Waves*, Cambridge Monographs on Mechanics and Applied Mathematic, Cambridge University Press, Cambridge, UK, 1980. - D. S. Jones,
*Acoustic and Electromagnetic Waves*, Oxford Science Publications, The Clarendon Press/Oxford University Press, New York, NY, USA, 1986. - A. Taflove,
*Computational Electrodynamics: The Finite-Difference Time-Domain Method*, Artech House, Boston, Mass, USA, 1995. - S. Sieniutycz and R. S. Berry, “Variational theory for thermodynamics of thermal waves,”
*Physical Review E*, vol. 65, no. 4, Article ID 046132, 11 pages, 2002. View at Publisher · View at Google Scholar - M. S. Mamis and M. Koksal, “Numerical solutions of partial differential equations for transmission lines terminated by lumped components,” in
*Proceedings of the 5th International Colloquium on Numerical Analysis*, pp. 87–98, 1996. - M. S. Mamis and M. Koksal, “Transient analysis of nonuniform lossy transmission lines with frequency dependent parameters,”
*Electric Power Systems Research*, vol. 52, no. 3, pp. 223–228, 1999. View at Publisher · View at Google Scholar - M. S. Mamiş and M. Köksal, “Remark on the lumped parameter modeling of transmission lines,”
*Electric Machines and Power Systems*, vol. 28, no. 6, pp. 565–575, 2000. View at Google Scholar - S. Herdem and M. S. Mamiş, “Computation of corona effects in transmission lines using state-space techniques,”
*Computers and Electrical Engineering*, vol. 29, no. 5, pp. 603–611, 2003. View at Publisher · View at Google Scholar - I. Abu-Alshaikh and M. E. Koksal, “One-dimensional transient dynamic response in functionally gradient spherical multilayered media,”
*Proceedings of Dynamical Systems and Applications*, pp. 1–20, 2004. View at Google Scholar - I. Abu-Alshaikh, “One-dimensional wave propagation in functionally graded cylindrical layered media,” in
*Mathematical Methods in Engineering*, pp. 111–121, Springer, New York, NY, USA, 2007. View at Google Scholar - M. S. Mamiş, A. Kaygusuz, and M. Köksal, “State variable distributed-parameter representation of transmission line for transient simulations,”
*Turkish Journal of Electrical Engineering and Computer Sciences*, vol. 18, no. 1, pp. 31–42, 2010. View at Publisher · View at Google Scholar - P. E. Sobolevskiĭ and L. M. Čebotareva, “Approximate solution of the Cauchy problem for an abstract hyperbolic equation by the method of lines,”
*Izvestija Vysših Učebnyh Zavedeniĭ Matematika*, no. 5(180), pp. 103–116, 1977 (Russian). View at Google Scholar - A. Ashyralyev and P. E. Sobolevskii, “A note on the difference schemes for hyperbolic equations,”
*Abstract and Applied Analysis*, vol. 6, no. 2, pp. 63–70, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Ashyralyev and P. E. Sobolevskii, “Two new approaches for construction of the high order of accuracy difference schemes for hyperbolic differential equations,”
*Discrete Dynamics in Nature and Society*, vol. 2005, no. 2, pp. 183–213, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Ashyralyev and M. E. Koksal, “On the second order of accuracy difference scheme for hyperbolic equations in a Hilbert space,”
*Numerical Functional Analysis and Optimization*, vol. 26, no. 7-8, pp. 739–772, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Ashyralyev and M. E. Koksal, “Stability of a second order of accuracy difference scheme for hyperbolic equation in a Hilbert space,”
*Discrete Dynamics in Nature and Society*, vol. 2007, Article ID 57491, 25 pages, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Ashyralyev, M. E. Koksal, and R. P. Agarwal, “A difference scheme for Cauchy problem for the hyperbolic equation with self-adjoint operator,”
*Mathematical and Computer Modelling*, vol. 52, no. 1-2, pp. 409–424, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Ashyralyev, M. E. Koksal, and R. P. Agarwal, “An operator-difference scheme for abstract Cauchy problems,”
*Computers & Mathematics with Applications*, vol. 61, no. 7, pp. 1855–1872, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. Ciment and S. H. Leventhal, “A note on the operator compact implicit method for the wave equation,”
*Mathematics of Computation*, vol. 32, no. 141, pp. 143–147, 1978. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - E. H. Twizell, “An explicit difference method for the wave equation with extended stability range,”
*BIT*, vol. 19, no. 3, pp. 378–383, 1979. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S. I. Piskarëv, “Stability of difference schemes in Cauchy problems with almost periodic solutions,”
*Differentsial'nye Uravneniya*, vol. 20, no. 4, pp. 689–695, 1984. View at Google Scholar · View at Zentralblatt MATH - S. Piskarev,
*Principles of Discretization Methods III. Report*, Acoustic Institute, Academy of Science USSR, 1986. - S. Piskarev, “Approximation of holomorphic semigroups,”
*Tartu Riikliku Ülikooli Toimetised*, no. 492, pp. 3–14, 1979. View at Google Scholar · View at Zentralblatt MATH - R. K. Mohanty, M. K. Jain, and K. George, “On the use of high order difference methods for the system of one space second order nonlinear hyperbolic equations with variable coefficients,”
*Journal of Computational and Applied Mathematics*, vol. 72, no. 2, pp. 421–431, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - H. O. Fattorini,
*Second Order Linear Differential Equations in Banach Spaces*, vol. 108 of*North-Holland Mathematics Studies*, North-Holland, Amsterdam, The Netherlands, 1985. - S. G. Kreĭn,
*Lineikhye differentsialnye uravneniya v Banakhovom prostranstve*, Nauka, Moscow, Russia, 1967. - R. K. Mohanty, “An unconditionally stable difference scheme for the one-space-dimensional linear hyperbolic equation,”
*Applied Mathematics Letters*, vol. 17, no. 1, pp. 101–105, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - R. K. Mohanty, “An unconditionally stable finite difference formula for a linear second order one space dimensional hyperbolic equation with variable coefficients,”
*Applied Mathematics and Computation*, vol. 165, no. 1, pp. 229–236, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - F. Gao and C. Chi, “Unconditionally stable difference schemes for a one-space-dimensional linear hyperbolic equation,”
*Applied Mathematics and Computation*, vol. 187, no. 2, pp. 1272–1276, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - H.-W. Liu and L.-B. Liu, “An unconditionally stable spline difference scheme of $O({k}^{2}+{h}^{2})$ for solving the second-order 1D linear hyperbolic equation,”
*Mathematical and Computer Modelling*, vol. 49, no. 9-10, pp. 1985–1993, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - R. K. Mohanty, “New unconditionally stable difference schemes for the solution of multi-dimensional telegraphic equations,”
*International Journal of Computer Mathematics*, vol. 86, no. 12, pp. 2061–2071, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. Dehghan and A. Shokri, “A numerical method for solving the hyperbolic telegraph equation,”
*Numerical Methods for Partial Differential Equations*, vol. 24, no. 4, pp. 1080–1093, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - H.-W. Liu, L.-B. Liu, and Y. Chen, “A semi-discretization method based on quartic splines for solving one-space-dimensional hyperbolic equations,”
*Applied Mathematics and Computation*, vol. 210, no. 2, pp. 508–514, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - L.-B. Liu and H.-W. Liu, “Quartic spline methods for solving one-dimensional telegraphic equations,”
*Applied Mathematics and Computation*, vol. 216, no. 3, pp. 951–958, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - P. E. Sobolevskiĭ and V. A. Pogorelenko, “Hyperbolic equations in a Hilbert space,”
*Sibirskiĭ Matematičeskiĭ Žurnal*, vol. 8, pp. 123–145, 1967 (Russian). View at Google Scholar - P. E. Sobolevskii,
*Difference Methods for the Approximate Solution of Differential Equations*, Izdat. Voronezh. Gosud. Univ., Voronezh, Russia, 1975. - B. Bülbül and M. Sezer, “Taylor polynomial solution of hyperbolic type partial differential equations with constant coefficients,”
*International Journal of Computer Mathematics*, vol. 88, no. 3, pp. 533–544, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Ashyralyev and M. E. Koksal, “On the numerical solution of hyperbolic PDEs with variable space operator,”
*Numerical Methods for Partial Differential Equations*, vol. 25, no. 5, pp. 1086–1099, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH