Fourth-order compact difference and alternating direction implicit schemes for telegraph equations
Introduction
The telegraph equation is one of the important equations of mathematical physics with applications in many different fields such as transmission and propagation of electrical signals [1], [2], vibrational systems [3], random walk theory [4], and mechanical systems [5], etc. The heat diffusion and wave propagation equations are particular cases of the telegraph equation.
Recently, increasing attention has been paid to the development, analysis, and implementation of stable methods for the numerical solutions of second-order hyperbolic equations. There have been many numerical methods for hyperbolic equations, such as the finite difference, the finite element, and the collocation methods, etc. See [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16], [17], [18], [19], [20], [21] and literatures are therein. Mohanty et al. [17], [18] developed some alternating direction implicit schemes for the two- and three-dimensional linear hyperbolic equations. Most of these schemes are second-order accurate in both space and time. Some conditionally stable fourth-order compact difference schemes for the wave and telegraph equations were introduced in [8], [9], [13]. An advantage of compact difference schemes is that the schemes use narrower stencils, i.e., fewer neighboring nodes, and have less truncation error comparing to typical finite difference schemes. Various versions of compact schemes have been successfully developed and analyzed (cf. [22], [23]) and references are therein. The alternating direction implicit (ADI) methods were first introduced by Douglas, Peaceman, and Rachford [24], [25], [26], [27], which have an advantage of no increase in dimension of the coefficient matrices corresponding matrix equations. There are many extensions and a great variety of applications of ADI based on the finite difference or the finite element methods [28], [29], [30], [31], [32], [33].
In this paper, we propose two- and three-level compact difference and ADI compact difference schemes with high accuracy for numerical solutions of the one- and two-dimensional telegraph equations with Dirichlet boundary conditions. These compact schemes are unconditionally stable and have a truncation error of second- and fourth-order in time and space, respectively. The coefficient matrices of the matrix equations given by the schemes are symmetric and tridiagonal. The matrix equations can be effectively solved by using many linear solvers. Numerical experiments presented for the 1D and 2D telegraph equations show that the present schemes are of high accuracy. By comparing the numerical results obtained by the present methods with the ones of some other available methods, the present methods can be considered as practical and effective numerical techniques to solve telegraph equations.
In Section 2, we introduce the compact difference schemes for the 1D telegraph equation and consider the stability and error estimates of the schemes. In Section 3, we introduce the compact difference and ADI compact difference schemes for the 2D telegraph equation, and analyze the stabilities and error estimates of the schemes. Finally, we present the result of numerical experiments on test examples in Section 4.
Section snippets
One-dimensional telegraph equation
Let . Consider the 1D linear telegraph equation with the initial and boundary conditions where c, b, and p are constants with c, , and , and the functions f, φ, and ψ are assumed to be sufficiently smooth.
Two-dimensional telegraph equation
Consider the 2D linear telegraph equation with the initial and boundary conditions where c, b, p, and q are constants with c, , p, , and f, φ, ψ, and g are smooth functions, , and Γ is the boundary of Q.
The solution domain is discretized into grids described by the set of nodes, in which , , ,
Numerical experiments
In this section, we present the result of numerical experiments to test the accuracy of the present schemes. The accuracy and convergence rates are measured by using the discrete - and -norms defined in Sections 2.2 Stability and error estimates, 3.3 Stability and error estimates.
Conclusion and discussion
We proposed the two- and three-level compact difference and ADI compact difference schemes with high accuracy for the 1D and 2D telegraph equations with Dirichlet boundary conditions. The two-level schemes are unconditionally stable and the three-level compact difference and ADI compact difference schemes are unconditionally stable for parameters θ in the range of . The proposed schemes are second- and fourth-order accurate in time and space, respectively. For the telegraph equations
Acknowledgements
The authors wish to express their sincere gratitude to the referees for their valuable comments and suggestions on this paper.
References (34)
- et al.
Mixed problems for the time-dependent telegraph equation: Continuous numerical solutions with a priori error bounds
Math. Comput. Modelling
(1997) - et al.
Computing the variable coefficient telegraph equation using a discrete eigenfunction method
Comput. Math. Appl.
(2007) - et al.
A new fourth-order compact finite difference scheme for the two-dimensional second-order hyperbolic equation
J. Comp. Appl. Math.
(2009) - et al.
A numerical algorithm for the solution of telegraph equations
Appl. Math. Comput.
(2007) - et al.
Analytic-numerical solution with a priori error bounds for coupled time-dependent telegraph equations: Mixed problems
Math Comput. Modelling
(1999) An unconditionally stable difference scheme for the one-space-dimensional linear hyperbolic equation
Appl. Math. Lett.
(2004)- et al.
High order difference schemes for the system of two space second order nonlinear hyperbolic equations with variable coefficients
J. Comp. Appl. Math.
(1996) - et al.
On the use of high order difference methods for the system of one space second order non-linear hyperbolic equations with variable coefficients
J. Comp. Appl. Math.
(1996) - et al.
A fractional step method on a special mesh for the resolution of multidimensional evolutionary convection–diffusion problems
Appl. Numer. Math.
(1998) - et al.
High-order schemes for acoustic waveform simulation
Appl. Numer. Math.
(2007)
Numerical method for viscous and non-viscous wave equations
Appl. Numer. Math.
Fourier Analysis and Boundary Value Problems
Digital signal propagation in dispersive media
J. Appl. Phys.
Differential Equations Elementary and Boundary Value Problems
Singularly perturbed telegraph equations with applications in the random walk theory
J. Appl. Math. Stoch. Anal.
Equations of Mathematical Physics
Higher order compact implicit schemes for the wave equation
Math. Comp.
Cited by (35)
Stability and convergence of BDF2-ADI schemes with variable step sizes for parabolic equation
2024, Applied Numerical MathematicsAn efficient compact difference-proper orthogonal decomposition algorithm for fractional viscoelastic plate vibration model
2023, Computers and Mathematics with ApplicationsNumerical analysis for compact difference scheme of fractional viscoelastic beam vibration models
2022, Applied Mathematics and ComputationCitation Excerpt :It is of important significance to model its vibration from mathematics and numerics, and thus lower the rates of catastrophic failure, increase the safety and reduce the overall cost of system maintenance. In this section, we shall recall briefly some basic concepts and properties concerning the finite differences [8,16,33], the compact differences (see [7,11,12,14,20,24,30,33,41,42]) and the weighted Grünwald differences [25,34] for the temporal fractional operator, then propose a compact difference scheme for (2.4) and prove its solvability. The right inequalities of (4.6) and (4.7) can be proved respectively in an analogous manner, by Lemma 4.2.
A reduced-order extrapolated finite difference iterative scheme for uniform transmission line equation
2022, Applied Numerical MathematicsCitation Excerpt :But due to the complication of the uniform transmission line equation in the real-world applications, it has usually analytical solution so that we have to rely on numerical solutions. There have been many numerical methods for finding the uniform transmission line equation (i.e., telegraph equation), such as the spectral Garlerkin method presented by Hesameddini and Asadolahifard in [7], the time discontinuous space-time finite element (FE) method proposed by He and Li in [6], some unconditionally stable implicit FD schemes proposed by Mohanty [23], the numerical approximation of higher-order time-fractional uniform transmission line equation studied by Hashemiand and Baleanu in [5], the new homotopy perturbation method addressed by Biazar and Eslami in [2], the meshless collocation method with barycentric rational interpolation proposed by Ma et al. in [22], the higher-order shifted Gegenbauer pseudospectral method introduced by Elgindy in [3], the differential quadrature algorithm established by Jiwari et al. in [9], the cosine expansion-based differential quadrature method addressed by Verma and Jiwari in [27], and some fourth-order compact difference and alternating direction implicit schemes proposed by Xie et al. in [29]. However, these numerical methods mentioned above include lots of unknowns.
Unified compact ADI methods for solving nonlinear viscous and nonviscous wave equations
2018, Chinese Journal of Physics
- 1
The research of this author is partially supported by NNSF of China grants 10971204.
- 2
The research of this author was supported by Changwon National University in 2009.