An Efficient Meshless Method for Hyperbolic Telegraph Equations in(1+1)Dimensions

Fuzhang Wang,Enran Hou,Imtiaz Ahmad,Hijaz Ahmad and Yan Gu

1College of Mathematics and Statistics,Xuzhou University of Technology,Xuzhou,221018,China

2College of Mathematics,Huaibei Normal University,Huaibei,235000,China

3Department of Mathematics,University of Swabi,Swabi,25000,Pakistan

4Department of Basic Sciences,University of Engineering and Technology,Peshawar,25000,Pakistan

5College of Mathematics,Qingdao University,Qingdao,266071,China

ABSTRACT Numerical solutions of the second-order one-dimensional hyperbolic telegraph equations are presented using the radial basis functions.The purpose of this paper is to propose a simple novel direct meshless scheme for solving hyperbolic telegraph equations.This is fulfilled by considering time variable as normal space variable.Under this scheme,there is no need to remove time-dependent variable during the whole solution process.Since the numerical solution accuracy depends on the condition of coefficient matrix derived from the radial basis function method.We propose a simple shifted domain method,which can avoid the full-coefficient interpolation matrix easily.Numerical experiments performed with the proposed numerical scheme for several second-order hyperbolic telegraph equations are presented with some discussions.

KEYWORDS Radial basis functions; telegraph equation; shifted domain method; meshless method

1 Introduction

The telegraph equation,which has been used to describe phenomena in various fields,belongs to the hyperbolic partial differential equation scope.For example,the telegraph equation in (1+ 1)dimensions can model the vibrations of structures,the digital propagation and also has applications in the other fields [1–3].Several methods are used to get the analytical/exact solutions of the telegraph equations [4–6].However,it is almost impossible to get the analytical solutions for relatively complex problems.Thus,numerical approximations to the telegraph equation is a better choice.Some numerical methods have been developed and compared to deal with the hyperbolic telegraph equations [7–9].

Especially,there are several numerical methods concentrate on the second-order 1D linear hyperbolic telegraph equations.For example,Mohanty et al.[10] investigated an unconditionally stable schemes for solving the hyperbolic equations.Based on the finite difference approximation,Dehghan et al.[11–13] presented several methods to simulate the linear hyperbolic telegraph equations.Lakestani et al.[14] used interpolating scaling function technique to solve the 1D hyperbolic telegraph equation.The boundary integral equation accompanied with the dual reciprocity method is used to solve the hyperbolic telegraph equations by Dehghan et al.[15].Pekmen and Tezer-Sezgin et al.[16] and Jiwari et al.[17] applied the differential quadrature method for the approximate solution of hyperbolic telegraph equations in one-and two space-dimensions.Zerarka et al.[18] considered the 2D generalized differential quadrature method for solving the hyperbolic telegraph equation later.The homotopy analysis method [19] is used to obtain the approximate analytical solution solutions of the second-order 1D linear hyperbolic telegraph equations.A pseudospectral method is proposed by Elgindy [20] for the second-order 1D hyperbolic telegraph equations.The B-spline collocation method is improved by Mittal et al.[21] to get the numerical solutions of the second order 1D hyperbolic telegraph equations.These numerical techniques are based on two-level difference or integral approximations.

Based on the above-mentioned investigations,we propose a direct meshless scheme with onelevel approximation for the second-order 1D linear hyperbolic telegraph equations.This is fulfilled by considering time variable as normal space variable.There is no need to remove time-dependent variable during the whole solution process.Under this scheme,we can solve the hyperbolic telegraph equations in a direct way.The rest paper is organized as follows.The formulation of the direct radial basis function is briefly introduced with the methodology for the hyperbolic telegraph equations in Section 2.To cope with the full coefficient matrix derived from the radial basis function method,we propose a simple shifted domain method in Section 3.Section 4 presented some numerical examples to validate the applicability of the proposed direct meshless scheme.Finally,some conclusions are given in Section 5.

2 The Direct Radial Basis Function

The general mathematical formulation of second-order linear hyperbolic telegraph equation in(1 + 1)dimensions is

in terms with the initial condition

and boundary conditions

Here,the coefficientα,βandδare non-negative constant,f (x,t)is the non-homogeneous/source term,g1(x),g2(x),andh(t)are prescribed functions.We aim to seek for the solution of unknown functionu(x,t).

Almost all numerical techniques for Eqs.(1)–(3)are based on the two-level approximations,most of which are based on the finite difference approximations.Here,we propose a direct collocation scheme by using radial basis function (RBF)under Euclidean space.

2.1 Direct Radial Basis Function

As is known to all,the traditional RBF methods are mostly used to solve 2D or higher dimensional problems.However,there is only one space variablexfor Eqs.(1)–(3),we propose a simple direct radial basis function (DRBF)by combining the space variablexand time variabletas a point(x,t)for (1 + 1)dimensional problems.More specifically,the interval [0,L] is evenly divided into segments firstly 0=x0

Figure 1:Configuration of the space-time coordinate

Actually,there is another definition of notradialnonmetric space-time radial basis functions with non-geometrical relationship between the space and the time.More details can be found in [22–24].

2.2 Methodology for Hyperbolic Telegraph Equations

According to the definition of DRBF,the above-mentioned Eqs.(1)–(3)can be solved directly in a one level approximation.More specifically,the numerical solution of a functionu(x,t)can be obtained from the following approximation

We should seek for the unknown coefficientsλj(j=1,2,...,n).

The interpolation scheme upon which the Eqs.(1)–(3)collocation is based is as below:

with

Here,we useNIto denote domain point number,N1is the point number on boundaryt=0 andN2is the point number on boundaryx=0 andx= 1 with 2N1+N2=NBdenotes the boundary point number.Eqs.(6)–(9)have the matrix form as where A is aN×Nknown square matrix and b is aN×1 vectors.This can be directly solved by the backslash computation in MATLAB codes.

3 The Shifted Domain Method

It should be noted that for a relatively large physical domain (with largeLorT),more collocation numbers are needed to ensure accuracy.For the problems considered in this paper,collocation methods will lead to a full coefficient matrix of linear algebraic equations.This has effect on the numerical solution accuracy.We propose a simple shifted domain method (SDM),which can deal with this problem easily.

The procedure of the SDM is shown by the above-mentioned physical domain Ω1=[0,1]×[0,T),which is also considered as the standard scope.For a larger domain [0,1]×[0,2T)with a larger time 2T,we can first consider the half domain Ω1,the other half domain Ω2=[0,1]×[T,2T)is considered as the shifted domain of Ω1.We can get the numerical solutions on Ω1and Ω2,respectively.Fig.2 presents the configuration of the shifted domain in the horizontal direction,wheret=Tis considered as an artificial boundary.The corresponding solution in the shifted domain Ω2can be get from the following equations:

Figure 2:Configuration of the shifted domain in the horizontal direction

This procedure is same as presented in the above-mentioned physical domain Ω1=[0,1]×[0,T).Eqs.(12)–(15),which can be directly solved by the backslash computation in MATLAB codes,have the matrix form as

where ˙A is aN×Nknown square matrix and ˙b is aN×1 vectors.

For the other cases,the configuration of the shifted domain in the horizontal direction [0,2]×[0,T)and both directions [0,2]×[0,2T)(with four sub-domains Ω1,Ω2,Ω3,Ω4)are shown in Figs.3 and 4,respectively.

4 Numerical Experiments

In this section,three examples are considered to validate the DRBF.For fair comparison with the other numerical methods,we use the maximum absolute error (MAE),absolute error and root mean square error (RMSE).The RMSE is defined as [25,26]

whereuis the analytical solution at test points,j=1,2,...,Ntandisthe numerical solutions at the test points,j=1,2,...,Nt.Ntis the number of test points on the physical domain.Since the parameterεin the DRBF method is non-sensitive,we fix parameterεfor all the following cases.The optimal choice of DRBF parameter is similar with the other radial basis functions.For more details about this topic,one can be found in [27,28] and references therein.

Figure 3:Configuration of the shifted domain in the vertical direction

Figure 4:Configuration of the shifted domain in both directions

4.1 Example 1

In order to investigate the DRBF method with the shifted domain method,we consider the hyperbolic telegraph Eq.(1)with the initial conditions

and boundary condition

The corresponding coefficients areα=4,β=2 with analytical/exact solution

The source term is

In this example,the physical domain Ω=[0,π]×[0,1] is divided into three sub-domains,i.e.,with shifted domainsIt should be noted that the physical domain division is non-unique.For fixed finenessand time step,we compare the DRBF results with the other numerical methods,detailed results with root mean square errors are listed in Tab.1.We note that our time step,which leads to less computations,is larger than the one Δt=0.0001 in [21,29,30] and Δt=0.01 in [15].Meanwhile,the finenessis also larger than the reference cases.However,the root mean square errors of DRBF is smaller than the other methods,i.e.,the DRBF is more accurate.We note that the DRBF method without the shifted domain method for this case performs not well.

Table 1:The root mean square errors (RMSE)of numerical methods for Example 1

For different mesh sizesH= 1/h,Fig.5 illustrates the root mean square error curve for Example 1 at timeT=1.It is seen that the DRBF solutions consistently converge very quickly.The convergence rate is high before reaching the minimum relative error value.It should be noted that the CPU time is 0.27 for mesh sizeH=1/4 and 1.92 for mesh sizeH=1/30.

Figure 5:Root mean square error curve for example 1 at time T=1

4.2 Example 2

In order to see the performance of the DRBF with different coefficients,we consider the hyperbolic telegraph Eq.(1)with the corresponding initial conditions

and boundary condition

The corresponding analytical/exact solution

with source term

For fair comparison with the other methods in [31–33],we consider the maximum absolute errors (MAE)in this example.For space finenessand time step,numerical results of the DRBF are listed in Tab.2 with different coefficientsαandβ.From which we can find that the DRBF results are stable and accurate for different coefficients.Forα=10 andβ=5,the solution accuracy of the DRBF is similar with the other methods.For the largerα=20 andβ=10,the solution accuracy of the DRBF performs the best.

Table 2:The maximum absolute errors (MAE)of numerical methods for Example 2

In order to see the difference between the exact and DRBF approximate solutions,Fig.6 shows the configuration for the larger coefficientsα=20 andβ=10 at timeT=0.5.From which we can see that the numerical solution coincides with the exact solution very well.

Figure 6:Exact and DRBF approximate solutions for Example 2 at time T=0.5

5 Conclusions

A new direct meshless scheme is presented for the second-order hyperbolic telegraph equations in (1 + 1)dimensions.The present numerical procedure,in which the time variable is considered as normal space variable,is based on the time-dependent radial basis function.There is no need to remove time-dependent variable during the whole solution process.Besides,a simple shifted domain method is proposed to cope with the solution accuracy related to the ill-conditioned coefficient matrix.From the numerical results in Section 4,we find that the proposed meshless method is superior to the other numerical methods.Besides,the direct meshless method can be extended to solve nonlinear problems with Newton iterative method considered.The DRBF with the shifted domain method is promising in dealing with the other types of time-dependent problems,fractional problems [34–36] as well as developing a parallel algorithm for large-scale problems.

Availability of Data and Materials:The data and material used to support the findings of this study are available from the corresponding author upon request.

Funding Statement:The first author is supported by the Natural Science Foundation of Anhui Province (Project No.1908085QA09)and the University Natural Science Research Project of Anhui Province (Project Nos.KJ2019A0591 & KJ2020B06).

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.