arxiv: v1 [math.na] 19 Nov 2016

Transcription

1 A new numerical algorithm for the solutions of hperbolic partial differential equations in ( + )-dimensional space Brajesh Kumar Singh arxiv:6697v [mathna] 9 Nov 6 Department of Applied Mathematics, School for Phsical Sciences, Babasaheb Bhimrao Ambedkar Universit Lucknow-6 5 (UP) INDIA bksingh584@gmailcom Abstract This paper deals with a construction of new algorithm: the modified trigonometric cubic B-Spline differential quadrature (MTB-DQM) for space discretization together with a time integration algorithm for numerical computation of the hperbolic equations Speciall, MTB-DQM has been implemented for the initial value sstem of the telegraph equations together with both Dirichlet and Neumann tpe boundar conditions The MTB-DQM is a DQM based on modified trigonometric cubic B-splines as new base functions The problem has been reduced into an amenable sstem of ordinar differential equations adopting MTB-DQM The resulting sstem of ordinar differential equations is solved using time integration algorithms Further, the stabilit of MTB-DQM is studied b computing the eigenvalues of the coefficients matrices for various grid points, which confirmed the stabilit of MTB-DQM for the telegraphic equations The accurac of the method has been illustrated in terms of the various discrete error norms for si test problems of the telegraph equation A comparison of computed numerical solutions with that obtained b the other methods has been carried out for various time levels considering various space sizes Kewords: Differential quadrature method, hperbolic telegraph equation, modified trigonometric cubic B- splines, MTB-DQM, Thomas algorithm Introduction The hperbolic partial differential equations have a great attention due to remarkable applications in fields of applied science and engineering, for instance, these equations model fundamental equations in atomic phsics [] and are ver useful in understanding various phsical phenomena in applied sciences and engineering It models the vibrations of structures (eg buildings, machines and beams) Consider the initial value linear sstem of second-order hperbolic telegraph equation in ( + ) dimension together with the following boundar conditions: u(,,t) t +α u(,,t) +β u(,,t) = u(,,t) t u(,,) = φ(,), u t (,,) = ψ(,), (,) Ω,t > + u(,,t) +f(,,t), () a) Dirichlet boundar condition: u(,,t) = φ (,t),u(,,t) = φ (,t),u(,,t) = φ 3 (,t),u(,,t) = φ 4 (,t),(,) Ω,t > () or b) Neumann boundar conditions: u (,,t) = ψ (,t), u (,,t) = ψ (,t), u (,,t) = ψ 3(,t), u (,,t) = ψ 4(,t),(,) Ω,t > (3) where Ω denotes the boundar of the computational domain Ω = [,] [,] R and α >,β are arbitrar constants, and ψ,φ,ψ i,φ i (i =,,3,4) are known smooth functions Eq () with β = is a damped wave equation while for β > it is known as hperbolic telegraph equation This equation is more convenient than

2 ordinar diffusion equation in modeling reaction diffusion for such branches of sciences [], and mostl used in wave propagation of electric signals in a cable transmission line [3] In the recent ears, various tpe of numerical techniques have been developed for solving partial differential equations [7,8,9,3,3,3,33,34], among others, initial value problems of telegraph equation, in one dimension, have been solved b Talor matri method [4], dual reciprocit boundar integral method [], unconditionall stable finite difference scheme [5], implicit difference scheme [8], variational iteration method [9], modified B- spline collocation method [6], Chebshev tau method [7], interpolating scaling function method [], cubic B- spline collocation method [] whereas two dimensional initial value problem of the telegraph equations have been solved b Talor matri method b Bülbül and Sezer [] which converts the telegraph equation into the matri equation, Two meshless methods- namel meshless local weak-strong (MLWS) and meshless local Petrov-Galerkin (MLPG) method b Dehghan and Ghesmati [], higher order implicit collocation method[3], A polnomial based differential quadrature method[4], modified cubic B-spline differential quadrature method[35], an unconditionall stable alternating direction implicit scheme [8], A hbrid method due to Dehghan and Salehi [5], compact finite difference scheme b Ding and Zhang [6] with accurac of order four in both space and time Two dimensional linear hperbolic telegraph equation with variable coefficients has been solved b Dehghan and Shorki [7] The differential quadrature method(dqm) has been developed b Bellman et al[9,] to obtain the numerical solutions of partial differential equations (PDEs), have received a great attention among the researchers After seminal work of Bellman et al[9,], and Quan and Chang [,], the DQMs have been emploed with various tpes of base functions, among others, cubic B-spline DQM [3,4], modified cubic B-spline differential quadrature method (MCB-DQM) [33,35], DQM based on fourier epansion and Harmonic function [36,37,38], sinc DQM [39], generalized DQM [4], polnomial based DQM [43,4], quartic B-spline based DQM [44], Quartic and quintic B-spline methods [45], eponential cubic B-spline based DQM [46], etended cubic B-spline DQM [47] In this paper our aim is to develop a new method so called modified trigonometric cubic-b-spline differential quadrature method (MTB-DQM) for hperbolic partial differential equations Speciall, MTB-DQM is emploed for the numerical computation of two dimensional second order linear hperbolic telegraph equation with both Dirichlet boundar conditions and Neumann boundar conditions The MTB-DQM is a new differential quadrature method based on modified trigonometric cubic-b-splines as base functions The MTB-DQM converts the initialboundar value sstem of the telegraph equation into an initial value sstem of ODEs, in time The resulting sstem of ODEs can be solved b using various time integration algorithm, among them, we prefer SSP-RK54 scheme (and SSP-RK43 scheme) [48,49] due to their reduce storage space, which results in less accumulation errors The accurac and adaptabilit of the method is illustrated b si test problems of second order linear hperbolic telegraphic equations in dimension ( + ) The rest of the paper is organized into five more sections, which follow this introduction Specificall, Section deals with the description of MTB-DQM Section 3 is devoted to the procedure for the implementation of MTB- DQM for the initial value problem () with the boundar conditions () and (3) The stabilit analsis of the MTB-DQM is studied in Section 4 Section 5 is concerned with the main aim, the numerical stud of si test problems, to establish the accurac of the proposed method in terms of the relative error norm (R e ), L and L error norms, considering various grid points Finall, Section 6 concludes the paper with reference to critical analsis and research perspectives Description of MTB-DQM The differential quadrature method is an approimation to derivatives of a function is the weighted sum of the functional values at certain nodes [] The weighting coefficients of the derivatives is depend onl on grids [4] This is the reason for taking the partitions P[Ω] of the problem domain Ω = {(,) R :, }

3 3 distributed uniforml as follows: P[Ω] = {( i, j ) Ω : h = i+ i,h = j+ j,i,j }, where = {,,,N }, = {,,,N }, and h = N and h = N are the discretization steps in both and directions, respectivel That is, a uniform partition in each, -direction with the grid points: = < < < i < < N < N = ; = < < < j < < N < N = For generic grid point ( i, j ), set u ij u ij (t) u( i, j,t), i,j Then r-th order derivative of u(,,t), for r {,}, with respect to, at ( i, j ) for i,j is approimated as follows: r ) u r = ij N l= a (r) il u lj, i ; r ) u r = ij N l= b (r) jl u il, j, (4) where the time dependent unknown quantities a (r) il, b(r) jl are the weighting functions of the rth-order derivatives, to be determine b adopting various tpe of base functions Trigonometric cubic B-spline function T i = T i () at node i in direction [5,6] read as: p T i = 3 ( i ), [ i, i+ ) p( i ){p( i )q( i+ )+p( i+ )q( i+3 )}+p ( i+ )q( i+4 ), [ i+, i+ ) ω q( i+4 ){p( i+ q( i+3 )+p( i+ )q( i+4 )}+p( i )q (5) ( i+3 ), [ i+, i+3 ) q 3 ( i+4 ), [ i+3, i+4 ) where p( i ) = sin( i ); q( i ) = sin( i ), ω = sin ( h ) sin(h )sin ( 3h ) If we take, Ti = ζ i+, then the set {ζ,,ζ N,ζ N+} forms a basis over the interval [a,b] Setting sin ( h ) a = sin(h )sin ( 3h ); a = +cos(h ) ; a 3 3 = 4sin ( 3h ); a 4 = 3+9cos(h ) 6sin ( h )( ( cos h ) ( +cos 3h )); a 5 = 3cos ( h ) sin ( h ) (+4cos(h )) The values of ζ i and its first and second derivatives in the grid point j, denoted b ζ ij := ζ i ( j ), ζ ij := ζ i ( j) and ζ ij := ζ i ( j), respectivel, read: a, if i j = { a ζ ij = a, if i j = ± ; ζ ij = ±a3, if i j = ± 5, if i j = ; ζ, otherwise ij = a 4, if i j = ± (6), otherwise otherwise The modified trigonometric cubic B-spline function are obtained b modifing the eponential cubic B-spline function (5) as follows [33]: ψ () = ζ ()+ζ () ψ () = ζ () ζ () ψ j () = ζ j (), for j = 3,4,,N (7) ψ N () = ζ N () ζ N+() ψ N () = ζ N ()+ζ N+()

4 4 Now, the set {ψ,ψ,,ψ N } be a basis over [a,b] The procedure for defining the modified trigonometric cubic B-splines in direction is followed analogousl The evaluation of the weighting coefficients a (r) ij and b (r) ij (r =,) In order to evaluate the weighting coefficients a () il of first order partial derivative in Eq (4), the modified trigonometric cubic B-spline ψ p (), p are used in DQM as base functions Setting ψ pi := ψ p ( i) and ψ pl := ψ p ( l ) The approimate values of the first-order derivative is obtained using TBM-DQM as follows: ψ pi = l a () il ψ pl, p,i (8) Setting Ψ = [ψ pl ], A = [a () il ] (the unknown weighting coefficient matri), and Ψ = [ψ pl ], then Eq (8) can be re-written as the following set of sstem of linear equations: where the coefficient matri Ψ of order N can be obtained from (6) and (7): Ψ = a +a a a a a a a and in particular the columns of the matri Ψ read: a 3 a 3 Ψ [] =,Ψ [] = a 3 a 3 ΨA T = Ψ, (9),,Ψ [N ] = a a a a a a a +a a 3 a 3, and Ψ [N ] = NoticethatcoefficientmatriΨ inequation(9)isdiagonalldominantforani,thecoefficientsa () i,a() i,,a() in have been determined b adopting Thomas Algorithm to solve sstem (9) The weighting coefficients b () il can be determined in similar manner, considering the grids in direction The weighting coefficients a (r) il and b (r) il, for r, can also be computed using the weighting functions in quadrature formula for second derivative on the given basis But, in the present paper, we prefer the following recursive formulae [4]: a 3 a 3

5 5 a (r) ij = r a (r) ii = b (r) ij = r ( ( a () ij a(r ) ii N i=,i j b () ij b(r ) ii ) a(r ) ij,i j : i,j, i j a (r) ij,i = j : i,j ) b(r ) ij,i j : i,j i j () b (r) ii = N i=,i j b (r) ij,i = j : i,j 3 Implementation of MTB-DQM for the telegraph equation Keeping all above in mind, the initial value sstem () of second order telegraph equation under the transformation u t = v (and so, u tt = v t ) transformed in to the initial value coupled sstem of first order differential equations as follows: u(,,t) t v(,,t) t = v(,,t) = αv(,,t) β u(,,t)+ u(,,t) u(,,) = φ(,), v(,,) = ψ(,), (,) Ω + u(,,t) +f(,,t),(,) Ω,t >, Setting f( i, j,t) = f ij The MTB-DQM transforms equation () to u ij = v ij t v ij = a () t il u lj + b () jl u il αv ij β u ij +f ij, l l u ij (t = ) = φ ij, v ij (t = ) = ψ ij, i, j Net, further simplification is not required in case of Dirichlet boundar conditions In this case the solution on the boundar can be read directl from the conditions () as: } u j = φ ( j,t) = φ (j), u Nj = φ ( j,t) = φ (j), j, t (3) u i = φ 3 ( i,t) = φ 3 (i), u in = φ 4 ( i,t) = φ 4 (i), i, But, for Neumann or mied tpe boundar conditions, further simplification is required In this case, the solutions at the boundar are obtained b using MTB-DQM on the boundar, which ields a sstem of linear equations On solving resulting sstem of the linear equation, the desired solution at the boundaries, is obtained The procedure for implementation of Neumann boundar conditions can be seen in [35], and is given below to complete the problem Now, the functions ψ,ψ from the boundar conditions (3) for r = ields the linear sstem in u j,u Nj, in compact matri form as [ ] a () a () [ ] [ ] N uj S (j) =, (4) S (j) a () N a() N N u Nj () ()

6 6 Equation (4) ields the solutions at the boundar of -ais u j = S (j)a () N N S (j)a () N a () a() N N a () N a() N, u Nj = S (j)a () S (j)a () N a () a() N N a () N a() N, j (5) where S (j) = ψ (j) N l u lj and S (j) = ψ (j) N N u l lj Similarl, the linear sstem obtained for ψ 3,ψ 4 from the boundar condition (3) ields the following l= a () u i = S 3(i)b () N N S 4 (i)b () N b () b() N N b () N b() N, u in = l= a () S 4 (i)b () S 3(i)b () N b () b() N N b () N b() N, i, (6) where S 3 (i) = ψ 3 (i) N l= b () l u il and S 4 (i) = ψ 4 (i) N l= b () N u l il Finall, on using boundar values u j,u Nj,u i and u in obtained from either boundar conditions (Dirichlet boundar conditions () or Neumann boundar conditions (3)), Eq () can be rewritten as follows: u ij = v ij t N where i N, j N and N a () il u lj + v ij = t l= l= u ij (t = ) = φ ij, v ij (t = ) = ψ ij, b () jl u il αv ij β u ij +K ij, K ij = f ij +a () i u j +a () in u Nj +b () j u i +b () jn u in (8) A lot time integration schemes have been proposed for the numerical computation of initial value sstem of differential equations (7), among others, the SSP-RK scheme allows low storage and large region of absolute propert [48,49] In what follows, we adopt SSP-RK43 scheme (and SSP-RK54 scheme) to integrate sstem from t m to t m +dt as given below, for getting the stable solution of hperbolic differential equations: SSP-RK43 Scheme: SSP-RK54 Scheme: k () = u m + t L(um ); k () = k () + t L(k() ); k (3) = 3 um + k() 3 + t 6 L(k() ) u m+ = k (3) + t L(k(3) ), k () = u m tl(u m ) k () = u m k () tl(k () ) k (3) = u m k () tl(k () ) k (4) = u m k (3) tl(k (3) ) u m+ = k () k (3) tl(k (3) ) k (4) tl(k (4) ) L() be linear differential operator as defined in (9) (7)

7 7 4 Stabilit analsis The sstem (7), in compact form, can be rewritten as: du = AU +G, dt (9) U(t = ) = U where [ ] [ ] [ ] [ ] O I O u φ ) A =, G =, U =, and U B αi K v = ψ ) O and O are null matrices; ) I is the identit matri of order (N )(N ); 3) U = (u,v) T the vector solution at the grid points: u = (u,u 3,,u (N ),u 3,u 33,,u 3(N ),,u (N ),,u (N )(N )) v = (v,v 3,,v (N ),v 3,v 33,,v 3(N ),,v (N ),,v (N )(N )) 4) K = (K,K 3,,K (N ),K 3,,K 3(N ),K (N ),K (N )(N ), where K ij, for i,j is calculated from Eq (8) 5) B = β I+B +B, where B and B are the following matrices (of order (N )(N )) of the weighting coefficients a () ij and b () ij : B = a () I a () 3 I a () (N I ) a () 3 I a () 33 I a () 3(N ) I a () (N I ) a () (N I )3 a () (N I )(N ) M = M O O O M O, B = O O M where identit matri, I, and null matri, O, both are of order (N ) and b () b () 3 b () (N ) b () 3 b () 33 b () 3(N ) b () (N ) b() (N )3 b() (N )(N ) The stabilit of MTB-DQM for the telegraph equation () depends on the stabilit of the sstem of ODEs defined in (9) and the proposed MTB-DQM for temporal discretization ma not converge to the eact solution whenever the sstem of ODEs (9) is unstable [5] The stabilit of (9) depends on the eigenvalues of the matri A [5] In fact, the stabilit region is the set S = {z C : R(z),z = λ A t}, where R() is the stabilit function and λ A is the eigenvalue of the matri A Speciall, the stabilit region for SSP-RK43 scheme [5, pp7] and SSP-RK54 scheme [5, Fig5], which shows that the sufficient condition for the stabilit of the sstem (9) is that for some t, λ A t S for each eigenvalue λ A the matri A, and so, the real part of each eigenvalue is necessaril less than or equal to zero Let λ A be an eigenvalue of A associated with the eigenvector (X,X ) T, where both X,X are vector of order (N )(N ), then from Eq (9), we get [ ] [ ][ ] [ ] X O I X X A = = λ X B αi X A, () X ()

8 8 The above equation ields the following two equations IX = λ A X, () and BX αx = λ A X (3) Eq () and Eq (3) ields BX = λ A (λ A + α)x, that is the eigenvalue λ B of the matri B is given b λ B = λ A (λ A +α), where the matri B is given b B = β I +B +B, (4) For different grid sizes: N N =,,3 3,4 4, the eigenvalues λ and λ have been depicted in Fig Eq (4) and Fig confirms that for different values of the grid sizes the computed eigenvalues λ B = λ β of B are real negative numbers, that is Re(λ B ) and Im(λ B ) =, (5) where Re(z) and Im(z) denote the real and the imaginar part of z, respectivel Setting λ A = +ι, then λ B = λ A (λ A +α) = +α+ι(+α) (6) Eq (5) and Eq (6) implies together that { +α < (+α) = (7) Eq(7) have the following possible solutions ) If, then = α, ) If =, then (+α) < α As α >, and so, in each case is negative Thus, for a given h one can find a sufficient small value of t so that for each eigenvalue λ A of matri A, tλ A belongs to the stabilit region S, and so, the MTB-DQM produces stable solutions for the sstem of second order hperbolic telegraph equations in ( + ) dimension 5 Numerical eperiments and discussion This section deals with the main goal of the paper is the computation of numerical solution of the telegraph equation obtained b new algorithm: MTB-DQM, with space step size h = h = h The accurac and the efficienc of the method are measured for si test problems in terms of the discrete error norms: L, L and relative error (R e ) [] which are defined as: L = ( h N k= e k) /, L = ma k N e k, and R e = ( N k= e k N k= u k where e k = u k u k ; u k and u k denote eact solution and numeric solution at node k, respectivel ) /

9 Problem Consider the telegraph equation () in the region Ω with α =,β =,f(,,t) = ep(+ t), andφ(,) = ep(+),ψ(,) = ep(+)inω,andthemiedboundarconditionsφ (,t) = ep( t),φ (,t) = ep(+ t) for and ψ 3 (,t) = ep( t),φ 4 (,t) = ep(+ t) for The eact solution [] is given b u(,,t) = ep(+ t) (8) In Table 6, the error norms: L,L,R e have been compared with the errors obtained due to MCB-DQM [35] taking the parameters h =, t = and h = 5, t = The surface plots of the MTB-DQM solutions and the eact solutions at different time levels t =,,4 have been depicted in Fig It is evident that the proposed solutions are more accurate in comparison to the results b Mittal and Bhatia [35] Problem Consider the telegraph equation () with α =,β =, f(,,t) = π ep( t)sinπ sinπ in region Ω,t > together with φ(,) = sinπsinπ,ψ(,) = sinπsinπ in Ω, and the mied boundar conditions ψ (,t) = πep( t)sin(π),φ (,t) = in, and φ 3 (,t) =,ψ 4 (,t) = πep( t)sin(π) in The eact solution as given in [] is u(,,t) = ep( t)sin(π)sin(π) The solutions are computed in terms of L,L error norms, for h =, t = and h = 5, t =, and are reported in Table 7 The surface plots of MTB-DQM solutions and eact solutions at t = 5,, are depicted in Fig 3 The findings from Table 7 concludes that the proposed MTB-DQM solutions are more accurate as compared to that obtained in [35] Problem 3 Consider the telegraph equation () with the eact solution u(,,t) = log(+++t) for α = β = as given in [], from which one can etract the functions φ and f directl The mied boundar conditions are taken as φ (,t) = log(++t),ψ (,t) = ++t for and ψ 3(,t) = ++t,φ 4(,t) = log(++t) for The error norms: L,L and R e has been computed with the parameters t =,h = 5 in the region Ω b adopting MTB-DQM (SSP-RK43), and are compared with the error norms due to MCB-DQM (SSP-RK43) [35] and Dehghan and Ghesmati [], in Table 8 It is seen from Table 8 that the computed MTB-DQM results are more accurate in comparison to the results obtained in [35,] The surface plots of the MTB-DQM solutions and the eact solutions at different time levels t =,,3 are depicted in Fig 4 Problem 4 Considerthetelegraphequation()intheregionΩ withα = β =,f(,,t) = (cost sint)sinsin, φ(,) = sinsin;ψ(,) =, and the Dirichlet boundar conditions: φ (,t) =,φ (,t) = costsin()sin,, φ 3 (,t) =,φ 4 (,t) = costsinsin(),, 9 (9) The eact solution [] is u(,,t) = costsinsin (3) We have compared the computed error norms: R e, L,L with the recent results of Mittal and Bhatia [35] at different time levels t, and reported them in Table for the parameters t =,h = while for t =,h = 5 in Table The findings shows that the proposed solution are much better than that in [35], and are in ecellent agreement with the eact solutions Problem 5 Consider the telegraph equation () with f(,,t) = ( α + β )ep( t)sinhsinh, φ(,) = sinhsinh,ψ(,) = sinhsinh, in Ω;φ (,t) = ;φ (,t) = ep( t)sinh()sinh for and ψ 3 (,t) = ;ψ 4 (,t) = ep( t)sinhsinh() for

10 The eact solution [4] is u(,,t) = ep( t)sinhsinh (3) Theerrornormsarecomputed forα =,β = 5andα =,β = withparameters t =,, h =,5 The computed L,L error norms for different time levels t are compared with that obtained b [35] in Table 3 taking t =,h = In Table 4, the computed error norms for t = and h = 5 are compared with that b Mittal and Bhatia [35] and Jiwari et al [4] The findings from the above tables confirms that the MTB-DQM results are more accurate in compression to that obtained b MCB-DQM [35] and PDQM [4] Problem 6 Considerthetelegraphequation()intheregionΩ withf(,,t) = ( 3cost+αsint+β cost)sinhsinh, and φ(,) = sinhsinh,ψ(,) = in Ω, and φ (,t) = o,φ (,t) = costsinh()sinh for, and φ 3 (,t) =,ψ 4 (,t) = costsinhsinh() for The eact solution [4] is u(,,t) = costsinhsinh (3) The error norms: L,L,R e at different time levels are computed for the parameters α =,β = 5 and α = 5,β = 5 taking t =,h = 5, and are reported in Table 5 The findings from Table 5 confirm that our results are comparabl more accurate than that obtained b MCB-DQM [35] and PDQM [4] 6 Conclusion In this paper, a new differential quadrature method based on modified trigonometric cubic B-splines as ne base functions has been developed, and so, call it modified trigonometric cubic B-spline differential quadrature method (MTB-DQM) for solving initial value sstem of linear hperbolic partial differential equations in (+ ) dimension Speciall, we implemented this method along with SSP-RK43 (and SSP-RK54) scheme for time integration for solving the second order hperbolic telegraph equation in ( + ) dimension subject to initial conditions, and each tpe of boundar conditions: Dirichlet, Neumann, mied boundar conditions The matri stabilit analsis method has been adopted to discus the stabilit of the MTB-DQM for the linear initial value sstem of second order hperbolic telegraph equations For different values of grid sizes: N N =,,3 3,4 4, we have computed the eigenvalues λ,λ to compute the eigenvalues λ B of the matri B for arbitrar β, which are found are real and negative Using this propert of λ B, we found that each eigenvalue λ A of the matri A has negative real part, which confirms that stabilit of the sstem of second order hperbolic telegraph equations in ( + ) dimension The findings from Section 5, shows that the MTB-DQM solutions are found comparativel more accurate than that obtained in [35,4,] Acknowledgement To be inserted References M Lakestani, BN Sara, Numerical solution of telegraph equation using interpolating scaling functions, Comput Math Appl 6 (7) () M Dehghan, A Ghesmati, Solution of the second-order one-dimensional hperbolic telegraph equation b using the dual reciprocit boundar integral equation (DRBIE) method, Eng Anal Bound Elem 34 () () D M Pozar, Microwave Engineering, Addison-Wesle, 99 4 B Bülbül, M Sezer, Talor polnomial solution of hperbolic tpe partial differential equations with constant coefficients, Int J Comput Math 88 (3) ()

11 5 F Gao, C Chi, Unconditionall stable difference schemes for a one-space-dimensional linear hperbolic equation, Appl Math Comput 87 () (7) RC Mittal, R Bhatia, Numerical solution of second order one dimensional hperbolic telegraph equation b cubic B-spline collocation method, Appl Math Comput (3) A Saadatmandi, M Dehghan, Numerical solution of hperbolic telegraph equation using the Chebshev tau method, Numer Methods Partial Differ Equ 6 () () RK Mohant, New unconditionall stable difference schemes for the solution of multi dimensional telegraphic equations, Int J Comput Math 86 () (9) M Dehghan, SA Yousefi, A Lotfi, The use of He s variational iteration method for solving the telegraph and fractional telegraph equations, Int J Numer Methods Biomed Eng 7 () () 9 3 S Sharifi, J Rashidinia, Numerical solution of hperbolic telegraph equation b cubic B-spline collocation method, Applied Mathematics and Computation 8 (6) 8-38 B Bülbül, M Sezer, A Talor matri method for the solution of a two-dimensional linear hperbolic equation, Appl Math Lett 4 () () 76-7 M Dehghan, A Ghesmati, Combination of meshless local weak and strong (MLWS) forms to solve the two dimensional hperbolic telegraph equation, Eng Anal Bound Elem 34 (4) () M Dehghan, A Mohebbi, High order implicit collocation method for the solution of two-dimensional linear hperbolic equation, Numer Methods Partial Differ Equ 5 () (9) R Jiwari, S Pandit, RC Mittal, A differential quadrature algorithm to solve the two dimensional linear hperbolic telegraph equation with Dirichlet and Neumann boundar conditions, Appl Math Comput 8 (3) () M Dehghan, R Salehi, A method based on meshless approach for the numerical solution of the two-space dimensional hperbolic telegraph equation, Math Methods Appl Sci 35 () () H Ding, Y Zhang, A new fourth-order compact finite difference scheme for the two-dimensional second-order hperbolic equation, J Comput Appl Math 3 () (9) M Dehghan, A Shokri, A meshless method for numerical solution of a linear hperbolic equation with variable coefficients in two space dimensions, Numer Methods Partial Differ Equ 5 () (9) RK Mohant, MK Jain, An unconditionall stable alternating direction implicit scheme for the two space dimensional linear hperbolic equation, Numer Methods Partial Differ Equ 7 (6) () R Bellman, B Kashef, ES Lee, R Vasudevan, Solving hard problems b eas methods: differential and integral quadrature, Comput Math Appl () (975) R Bellman, BG Kashef, J Casti, Differential quadrature: a technique for the rapid solution of nonlinear partial differential equations, J Comput Phs (97) 4-5 JR Quan, CT Chang, New insights in solving distributed sstem equations b the quadrature method-ii, Comput Chem Eng 3 (9) (989) 7-4 JR Quan, CT Chang, New insights in solving distributed sstem equations b the quadrature method-i, Comput Chem Eng 3 (7) (989) A Korkmaz, I Dag, Cubic B-spline differential quadrature methods and stabilit for Burgers equation, Eng Comput Int J Comput Aided Eng Software 3 (3) (3) A Korkmaz, I Dag, Numerical simulations of boundar - forced RLW equation with cubic B-spline-based differential quadrature methods, Arab J Sci Eng 38 (3) G Arora, V Joshi, A computational approach for solution of one dimensional parabolic partial differential equation with application in biological processes, Ain Shams Eng J (6), 6 Muhammad Abbas, Ahmad Abd Majid, Ahmad Izani Md Ismail, Abdur Rashid, The application of cubic trigonometric B-spline to the numerical solution of the hperbolic problems, Applied Mathematics and Computation 39 (4) Ozlem Erso and Idris Dag, Numerical solutions of the reaction diffusion sstem b using eponential cubic B-spline collocation algorithms, Open Phs 3(5) G Arora, RC Mittal and B K Singh, Numerical Solution of BBM-Burger Equation with Quartic B-spline collocation method, Journal of Engineering Science and Technolog, Special Issue, /4, V K Srivastava and B K Singh, A Robust finite difference scheme for the numerical solutions of two dimensional time-dependent coupled nonlinear Burgers equations, International Journal of Applied Mathematics and Mechanics (7) (4) B K Singh and G Arora, A numerical scheme to solve Fisher-tpe reaction-diffusion equations, Nonlinear Studies/MESA- MATH- EMATICS IN ENGINEERING, SCIENCE AND AEROSPACE 5() (4) B K Singh, G Arora, M K Singh, A numerical scheme for the generalized Burgers-Hule equation, Journal of the Egptian Mathematical Societ (6) 3 B K Singh and Carlo Bianca, A new numerical approach for the solutions of partial differential equations in three-dimensional space, Appl Math Inf Sci, No 5, - (6) 33 G Arora, B K Singh, Numerical solution of Burgers equation with modified cubic B-spline differential quadrature method Applied Math Comput 4 (3) 66-77

12 34 BK Singh, P Kumar, A novel approach for numerical computation of Burgers equation (+) and (+) dimension, Aleandria Eng J (6) 35 RC Mittal, Rachna Bhatia, A numerical stud of two dimensional hperbolic telegraph equation b modified B-spline differential quadrature method, Applied Mathematics and Computation 44 (4) C Shu, YT Chew, Fourier epansion - based differential quadrature and its application to Helmholtz eigenvalue problems, Commun Numer Methods Eng 3 (8) (997) C Shu, H Xue, Eplicit computation of weighting coefficients in the harmonic differential quadrature J Sound Vib 4(3) (997) AG Striz, X Wang, CW Bert, Harmonic differential quadrature method and applications to analsis of structural components, Acta Mech (-)(995) A Korkmaz, I Dag, Shock wave simulations using sinc differential quadrature method, Eng Comput Int J Comput Aided Eng Software 8 (6) () C Shu, B E Richards, Application of generalized differential quadrature to solve two dimensional incompressible navier-stokes equations, Int J Numer MethFluids 5 (99) A Korkmaz, AM Akso, I Dag, Quartic B-spline differential quadrature method, Int J Nonlinear Sci (4) () A Korkmaz, I Dag, Polnomial based differential quadrature method for numerical solution of nonlinear Burgers equation, J Franklin Inst 348 () () A Korkmaz Numerical algorithms for solutions of Kortewegde Vries equation Numerical methods for partial differential equations, 6(6) () A Bashan, S B G Karakoc, T Geikli, Approimation of the KdVB equation b the quintic B-spline differential quadrature method Kuwait J Sci 4 () (5) A Korkmaz, I Dag, Quartic and quintic Bspline methods for advection diffusion equation Applied Mathematics and Computation 74 (6) A Korkmaz, and H K Akmaz, Numerical Simulations for Transport of Conservative Pollutants Selcuk Journal of Applied Mathematics 6() (5) 47 A Korkmaz and HK Akmaz, Etended B-spline Differential Quadrature Method for Nonlinear Viscous Burgers Equation, Proceedings of International Conference on Mathematics and Mathematics Education, pp 33-33, Elaziǧ, Turke -4 Ma, 6 48 S Gottlieb, D I Ketcheson, C W Shu, High Order Strong Stabilit Preserving Time Discretizations, J Sci Comput 38 (9) J R Spiteri, S J Ruuth, A new class of optimal high-order strong stabilit-preserving time-stepping schemes, SIAM J NumerAnalsis 4 () () MK Jain, Numerical Solution of Differential Equations, nd ed, Wile, New York, NY, Ethan J Kubatko, Benjamin A Yeager, David I Ketcheson, Optimal strong-stabilit-preserving Runge-Kutta time discretizations for discontinuous Galerkin methods, ssp stabilitpdf

13 3 7 List of Tables and Figures Table Comparison of the error norms for Problem 4 for h =, t = MCB-DQM (SSP-RK43) [35] MTB-DQM SSP-RK54 SSP-RK43 t L L R e L L R e L L R e 997E-4 746E E E E E-5 486E E-6 649E-5 96E-3 876E-3 859E E E E-5 35E E-6 93E E-4 688E-4 747E E E-6 36E E E-6 78E E-3 994E-3 767E-3 443E E E E-6 494E-6 785E E-4 878E-3 357E-3 383E E-6 59E E E-6 85E E-4 558E-3 874E-3 385E E E-5 576E E-6 54E-5 Table Comparison of the error norms for Problem 4 for h = 5, t = MCB-DQM (SSP-RK43) [35] MTB-DQM SSP-RK54 SSP-RK43 t L L R e L L R e L L R e 9887E E E E E-7 399E E E-7 444E-6 48E-4 396E-4 5E E-7 676E E E E-7 38E E-5 993E-5 378E-4 789E-7 39E E-6 993E-7 99E-6 336E E-4 335E-4 54E E E E-6 887E-7 348E-7 788E E E-4 389E-4 546E-7 889E-7 486E-6 73E-7 39E E E-5 35E-4 45E-4 586E E-7 5E E-7 98E E-6 Table 3 Comparison of the error norms for Problem 5 with α =,β = 5 and t = and h = t MTB-DQM MCB-DQM [35] SSP-RK54 SSP-RK43 SSP-RK43 L L R e L L R e L L R e 5 69E-6 9E-5 68E E E-6 696E E-4 339E-3 89E-3 536E-6 748E E-5 363E-6 494E-6 953E-5 654E-4 597E-3 348E-3 37E-6 858E E-5 35E E-6 95E-5 467E E E E-7 6E-6 389E E E-7 83E E-5 79E E E-7 44E E-5 679E-8 933E-8 5E-5 68E-5 378E E-3

14 4 Table 4 Comparison of the error norms for Problem 5 with δt =,α =,β =,5 and h = 5 t MTB-DQM MCB-DQM [35] PDQM [4] SSP-RK54 SSP-RK43 SSP-RK43 β = 5 L L R e L L R e L L R e R e 5 835E-7 36E-6 385E-6 697E-7 5E-6 75E-6 69E E-4 88E-4 85E E E E E E E-6 593E-5 338E-4 366E-4 85E-4 35E-7 3E-7 437E-6 84E E E E-5 38E E E E-8 945E-7 47E E-8 587E E-6 994E E-5 34E-4 656E E-8 63E-8 46E E E-8 385E-6 75E-6 544E E-4 977E-4 β = 5 79E-7 98E E E E-6 369E E E E-4 98E-4 878E-7 67E-6 53E E-7 946E-7 43E E E E E E E-7 4E-5 356E E E-6 554E E E E E E E-5 733E-7 5E E E E-5 454E-4 44E E E E E E-8 395E-5 56E E E-4 883E-4 Table 5 Comparison of the error norms for Problem 6 with α =,5,β = 5 and h = 5, t = t MTB-DQM MCB-DQM [35] PDQM [4] SSP-RK43 SSP-RK54 SSP-RK43 α = L L R e L L R e L L R e R e 5 97E E-6 496E-6 867E E E-6 696E E E-5 49E E-7 533E-6 399E-6 549E E-6 4E-5 774E E E-5 463E-5 896E E E E-6 794E-6 795E E-4 598E-4 447E E E E-6 384E-6 44E-6 66E-6 336E E E-5 66E E E E E E-6 478E E E E E E E E-6 848E E-6 389E E-6 4E-4 473E-4 37E E-7 933E-7 35E E-6 464E E-6 4E-4 4E-4 654E-5 49E-5 α = E E-6 634E-6 3E-6 384E E-6 988E E-4 735E-5 874E-5 593E-7 79E E E E-6 46E E E-4 354E-5 764E-5 35E E-7 59E-6 47E E-6 37E-5 79E-4 557E E E E-7 567E-6 663E-6 67E E-6 394E-6 746E E E E E E E E E-6 648E-5 84E E E-5 84E E-7 54E-6 447E E-6 53E-6 8E-5 376E E-4 647E-5 796E-7 67E-6 439E E E-6 956E-6 69E E-4 577E-5 44E-5

15 5 Table 6 Comparison of the error norms for Problem with α =,β = for different values of t,h t MTB-DQM: t = h = MCB-DQM: t = h = [35] SSP-RK54 SSP-RK43 SSP-RK43 L L R e L L R e L L R e 3979E-4 676E-4 386E E E E-4 444E- 9996E- 89E E-5 8E-4 5E-4 475E-5 33E-4 574E E-3 397E E E E E E-5 88E E-4 38E-3 78E E E E E E E E-4 E-4 68E E E-7 446E-6 575E E-7 65E-6 63E E-5 35E-5 485E E-8 746E-8 64E E E-8 978E-4 448E E E-4 h = 5, t = 5 789E-4 67E E-5 998E E E E-3 959E-3 845E-5 533E-4 845E-4 884E-5 89E-4 863E-4 97E-5 335E E-3 96E-4 444E-5 373E-5 389E-5 54E E-5 399E-5 858E-4 36E E E-5 57E E-5 55E-5 8E-5 69E-5 38E E E-5 5 9E-6 399E-6 477E-5 364E-6 43E E E E E E E E-5 9E E E E-6 83E-6 484E-5 78E-8 778E-8 85E-5 98E-8 79E-8 97E E E-6 8E-5 Table 7 Comparison of the error norms for Problem with α = β = for different values of t,h t MTB-DQM: t = h = MCB-DQM: t = h = [35] SSP-RK54 SSP-RK43 SSP-RK43 L L L L L L 57445E-4 779E E E-4 644E-3 366E-3 74E-4 454E-4 853E E E E E-5 56E-5 87E-5 567E E-4 479E E E E E-6 48E-4 3E E-6 636E-6 83E-7 39E E-5 66E E E E E E-6 483E-6 h = 5, t = 5 369E E E E E-4 959E-4 55E-5 744E E E-5 335E E-4 554E-5 383E-5 49E E-5 858E-5 36E E-7 53E E-7 493E-6 38E E E E-7 45E-7 63E E E E E-7 34E-7 697E E-7 83E-7 783E E-9 564E E E E-8 Table 8 Comparison of the error norms for Problem 3 with α = β = taking h = 5, t = t MTB-DQM (SSP-RK43) MCB-DQM (SSP-RK43)[35] Dehghan and Ghesmati [] L L R e L L R e R e :MLWS R e :MLPG E E-5 544E-5 69E E-3 88E E-5 999E E-5 833E E-5 593E-3 338E-3 366E-3 998E-5 798E-5 979E-5 49E-5 4E E-4 38E E-4 875E E E-5 6E-5 734E-6 994E E-4 34E-4 993E-4 48E E-5 969E E-6 75E-4 544E E-5 473E-3 69E E-6 5E-5 493E-6 7E-4 348E-4 845E-5 73E E E E-6 94E-6 685E-6 449E-6 539E E-5 437E-4 964E-5

16 6 8 6 N N = N N = N N = 3 3 N N = Im(λ) Re(λ ) 8 6 N N = N N = N N = 3 3 N N = Im(λ) Re(λ ) Fig Eigenvalues λ and λ of B and B resp, for different grid sizes N N =,,3 3,4 4

17 7 Num Sol n (t= ) Num Sol n (t=) Num Sol n (t=4) Fig Plots of MTB-DQM solutions at different time levels for Problem Num Sol n (t=) Num Sol n (t=) Num Sol n (t=5) Fig3 Plots of MTB-DQM solutions at different time levels for Problem 4 8 Num Sol n (t=) Num Sol n (t=3) Num Sol n (t=5) Fig4 Plots of MTB-DQM solutions at different time levels for Problem 3