Numerical Analysis of an H 1-Galerkin Mixed Finite Element Method for Time Fractional Telegraph Equation

We discuss and analyze an H 1-Galerkin mixed finite element (H 1-GMFE) method to look for the numerical solution of time fractional telegraph equation. We introduce an auxiliary variable to reduce the original equation into lower-order coupled equations and then formulate an H 1-GMFE scheme with two important variables. We discretize the Caputo time fractional derivatives using the finite difference methods and approximate the spatial direction by applying the H 1-GMFE method. Based on the discussion on the theoretical error analysis in L 2-norm for the scalar unknown and its gradient in one dimensional case, we obtain the optimal order of convergence in space-time direction. Further, we also derive the optimal error results for the scalar unknown in H 1-norm. Moreover, we derive and analyze the stability of H 1-GMFE scheme and give the results of a priori error estimates in two- or three-dimensional cases. In order to verify our theoretical analysis, we give some results of numerical calculation by using the Matlab procedure.

Recently, some people have made use of the method to obtain the numerical solution for some partial differential equations since Pani (in 1998) [33] proposed an 1 -GMFE method. This method includes some advantages, such as avoiding the LBB consistency condition, allowing different polynomial degrees of the finite element spaces, and obtaining the optimal a priori estimates in both 1 and 2 -norms. In [34], Pani and Fairweather discussed some detailed a priori error results of two numerical schemes based on the 1 -GMFE method for linear parabolic integrodifferential equations. In [35], Pani et al. gave some error analysis based on 1 -GMFE scheme for the partial differential equation of hyperbolic type. At the same time, a modified 1 -GMFE procedure was also proposed and analyzed for the case of two or three dimensions. Guo and Chen [36,37] obtained some theoretical error analysis and numerical results of 1 -GMFE method for RLW equation and Sobolev equation. In 2013, Liu et al. [38] introduced another auxiliary variable, which is different from the one in [36], and then proposed and studied an explicit multistep 1 -GMFE scheme for a RLW equation. Liu and Li [39] and Zhou [40] gave some different discussions on 1 -GMFE method for pseudohyperbolic equations (heat transport equation), respectively. In [41], Che et al. studied the 1 -GMFE method for a nonlinear integrodifferential equation. Recently, Shi et al. [42,43] proposed some nonconforming mixed scheme based on the 1 -GMFE method.
Based on the above review on 1 -GMFE method, we easily see that the method was studied based on the integer-order partial differential equations. However, the theoretical results of 1 -GMFE method for fractional telegraph equation have not been presented and analyzed.
In this paper, our aim is to give some detailed a priori error analysis and numerical results on the 1 -GMFE method for time fractional telegraph equation. We apply the difference schemes to approximate the time fractional derivatives and use the 1 -GMFE method to discretize spatial direction. We obtain some optimal a priori error results of one dimension for the scalar unknown in 2 and 1 -norms. Moreover, we also get an a priori error result of the optimal 2 -norm for the auxiliary variable. At the same time, we use the 1 -GMFE method to deal with the cases in several dimensions and analyze the stable results for the 1 -GMFE scheme. We calculate some numerical results to verify the theoretical analysis of 1 -GMFE method for fractional telegraph equation. The layout of the paper is as follows. In Section 2, we formulate an 1 -GMFE scheme for time fractional telegraph equation (5). In Section 3, we introduce some lemmas of two important projections and two difference approximations for time fractional derivatives and then analyze some a priori error results. In Section 4, some theoretical results are given in the cases of two and three dimensions. In Section 5, we choose a numerical example to verify the theoretical analysis of our method. In Section 6, we give some remarks and extensions about the 1 -GMFE method for fractional PDEs.

An 1 -GMFE Scheme in One Space Dimension
In this section, we first consider the 1 -GMFE method for the following time fractional telegraph equation in 1D case: with boundary condition and initial conditions where Ω = [ , ] ⊂ . In order to get the 1 -GMFE formulation, we first introduce an auxiliary variable = ( , )/ and split (5) into the following first-order system by Multiplying (8) by − / , ∈ 1 , integrating with respect to space from to , and using an integration by parts with ( , )/ = ( , )/ = 0 and 2 ( , )/ 2 = 2 ( , )/ 2 = 0, we easily get Multiply (9) by V/ , V ∈ 1 0 , and integrate with respect to space from to to obtain For formulating finite element scheme, we now choose the finite element spaces ℎ ⊂ 1 0 and ℎ ⊂ 1 , which satisfy The Scientific World Journal 3 the following approximation properties: for 1 ≤ ≤ ∞ and , positive integers [33], Based on the chosen finite element spaces, the semidiscrete 1 -GMFE scheme is described by

Two Projection Lemmas.
For a priori error estimates for fully discrete scheme, we introduce two projection operators [33,44] in Lemmas 1 and 2.

Lemma 1.
One defines a Ritz projection ℎ ∈ ℎ for the variable by Then the following estimates hold, for = 0, 1: Lemma 2. Further, one also defines an elliptic projection ℎ ∈ ℎ of as the solution of where B( , ) = ( , ) + ( , ). Here > 0 is chosen to satisfy Then the following estimates are found: for = 0, 1,

Approximation of Time-Fractional Derivative.
For formulating fully discrete scheme, let 0 = 0 < 1 < 2 < ⋅ ⋅ ⋅ < = be a given partition of the time interval [0, ] with step length Δ = / and nodes = Δ , for some positive integer . For a smooth function on [0, ], define = ( ). In the following analysis, for deriving the convenience of theoretical process, we now denote Now, we will introduce two lemmas on the approximations of time fractional derivatives.

Error Estimates for Fully Discrete Scheme. Based on the approximation formulas
The Scientific World Journal Now, we formulate a fully discrete procedure: find Making a combination of (24)-(25) with two projections (14) and (16), we get the following error equations: where In the following discussion, we will derive the proof for the fully discrete a priori error estimates.
Then there exists a positive constant ( , , , ) free of space-time discrete parameters ℎ and Δ such that, for and for → 1 − ℎ ≤ ( , , ) (Δ −1 ℎ min{ +1, +1} + Δ ) Proof. For the need of error analysis, we first consider the 2 -norm ‖ ‖ and the 1 -norm ‖ ‖ 1 . Taking V ℎ = in (26) and using Poincaré inequality based on the 1 0 -space and Cauchy-Schwarz inequality, we easily get In the next analysis, we will give the estimates of ‖ ‖ and ‖ ‖ in 2 -norm. Noting that The Scientific World Journal 5 then (27) may be rewritten as We take ℎ = in (35) and multiply by to arrive at By the simple calculation, we get the following equality: By applying the similar process of calculation to (38), we use Cauchy-Schwarz inequality to get Noting (33), we use Cauchy-Schwarz inequality and Young inequality to have Substitute (38), (39), and (40) into (37) and use Cauchy-Schwarz inequality to arrive at The Scientific World Journal By an application of Young inequality, we get Noting that 2 − 2 +1 > 0, we arrive at At the same time, noting that − +1 > 0, we use the similar method to have By a combination of (43) and (44) with (18) and noting that The Scientific World Journal 7 Substitute (45) into (42) and note that 0 In the following discussion, we apply mathematical induction to obtain the error result Take = 1 in (46) and note that 2 0 = 1; it is easy to find that the following inequality holds: Assuming that the inequalities ‖ ‖ ≤ ( ( , , , )/ 2 −1 )(ℎ +1 +ℎ +1 +Δ 2 +Δ 1+2 ) hold for = 1, 2, . . . , −1, we now prove that the inequality (47) holds. Noting that 1/ < 1/ +1 , = 1, 2, and combining (43) and (44) with (46), we have Based on the process of mathematical induction, we claim that (47) holds. For the need of the next proof, we have to estimate the term 1−2 / 2 −1 . We now use Taylor formula to arrive at The Scientific World Journal where 0 < < 1 and 1/2 < < 1. Noting that ( Δ ) 2 −1 ≤ 2 −1 and (50), we have Substitute (51) into (33) and use (18) to get By combining (51) and (18) with triangle inequality, the 2norm estimate ‖ − ℎ ‖ is got. Similarly, the estimate ‖ − ℎ ‖ in the 2 -norm and the estimate ‖ − ℎ ‖ 1 in the 1norm also are obtained by a combination of (52) and (15) with triangle inequality. Now we mainly analyze the case → 1. We can find that the error inequality (29) in this case has no meaning since the coefficient ( , , , )/(2 − 2 ) → ∞ as → 1. So, we have to look for another error estimate's process. Noting the fact that Δ ≤ ( = 1, 2, . . . , ), we can obtain the following error inequality: Now we can use induction to prove the inequality (53). The detailed proof is similar to the above process of analysis, so we do not give the detailed proof again. Making a combination of (53) and (18) with triangle inequality, we can get the estimate (31). A similar discussion for (32) can also be made.

An 1 -GMFE Scheme for Several Spaces Variables
In this section, we consider (1) with two and three spaces variables. Let L 2 (Ω) = ( 2 (Ω)) , ( = 2 or 3) with inner product and norm Further, let Taking q = ∇ , we use a similar process to the system (10) The corresponding time semidiscrete system is defined by where R 0 = (Δ ), which can be estimated by a similar process to 0 .
The Scientific World Journal 9 In order to get fully discrete mixed finite element scheme, we now choose the finite element spaces ℎ ⊂ 1 0 and ℎ ⊂ W, which satisfy the following approximation properties: for 1 ≤ ≤ ∞ and , positive integers [33], The fully discrete 1 -GMFE scheme is to find

The Analysis of Stability.
In the following discussion, we will give the stability for the system (60) and (61). First, we need to obtain an important lemma on two initial value conditions.

Lemma 6.
With q 0 = ∇ 0 and q (0) = ∇ 1 , the following inequality holds: Proof. By a similar discussion as in [32], we have Based on the given initial value conditions and (63), we arrive at Theorem 7. The following stable inequality for the system (60) and (61) holds: Proof. In (60), we take V ℎ = ∇ ℎ and use Cauchy-Schwarz inequality and Poincaré inequality to arrive at In (61), we choose w ℎ = q ℎ and use Cauchy-Schwarz inequality, Young inequality, and (66) to get By the simple simplification for (67), we easily get For the case = 1 in (68), we multiply (36) to get easily Noting that 0 = 1 = 2 0 and Lemma 6, we use Cauchy-Schwarz inequality to get By the simple calculation for (70), we arrive at 10 The Scientific World Journal For the case ≥ 1, using a similar process to the proof of Theorem 5 based on the mathematical induction, we can get By a combination of (72) and (66), we get the conclusion of Theorem 7.

A Priori Error Results.
For deriving the a priori error analysis, we define the Ritz projectioñℎ ∈ ℎ by Further, letq ℎ ∈ ℎ be the standard finite element interpolant of q. Let = q −q ℎ , = −̃ℎ; see [33,44]; we obtain ≤ ℎ +1− ‖ ‖ +1 , = 0,1, Now, we can get the following theorem of a priori error estimates based on the above contents. there exists a positive constant ( , q, , ) free of space-time meshes ℎ and Δ such that, for 1/2 < < 1, and for → 1 Proof. We can use a similar proof as in Theorem 5 to get the conclusion of Theorem 8, so we do not discuss that again.
In Table 1, we calculate the numerical results of a priori error results and orders of convergence in 2 -norm for both and with a fixed time step length Δ = 1/4000 and the changed spatial meshes ℎ 1 = 1/20, ℎ 2 = 1/40, and ℎ 3 = 1/80. From the calculated data, we can clearly see that the orders of convergence in 2 -norm keep unchanged with the changed = 0.6, 0.7, 0.8, 0.9, which confirm the optimal second-order convergence results of 1 -Galerkin mixed finite element method. Similarly, we obtain optimal first-order rates of convergence of 1 -norm for and in Table 2.
In Table 3, we calculate the numerical results of a priori error results and orders of convergence in 2 -norm for both and with different space-time step length ℎ = 5Δ = 1/ ( = 20, 40, 80). From the calculated data, we can find that the rates of convergence, which are higher than the results (Δ 3−2 + ℎ 2 ) of theory, gradually decrease with the increased (which is taken from 0.6 to 0.9 with interval 0.1). In Table 4, some first-order convergence results for both and in 1 -norm, which are unchanged with the changed = 0.6, 0.7, 0.8, 0.9, are given. In view of the above analysis on the numerical results, we now announce that the time fractional telegraph equation can be well solved by the 1 -GMFE method.

Some Concluding Remarks and Extensions
As far as we know, more and more people have proposed and analyzed a lot of numerical methods for fractional partial differential equations. However, the discussions on mixed finite element methods for solving fractional partial differential equations are fairly limited. The a priori error analysis of 1 -GMFE method for time fractional telegraph equation, especially, has not been made and discussed. In this paper, we give the detailed proof 's process of the error analysis on 1 -GMFE method for time fractional telegraph equation. Further, we provide a numerical procedure to verify the theoretical results of the studied method.

12
The Scientific World Journal  In the near future, our aim is to study an 1 -Galerkin moving mixed finite element method, which is based on a combination of 1 -GMFE methods and moving finite element methods.