Two New Error Estimates of a Fully Discrete Primal-Dual Mixed Finite Element Scheme for Parabolic Equations in Any Space Dimension

We consider a fully discrete implicit scheme in which the discretization in space is performed using the PDMFEM (Primal-Dual Mixed Finite Element Method) defined in Ciarlet (The finite element method for elliptic problems. Classics in applied mathematics, 40. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, pp 415–416, 2002) and Quarteroni and Valli (Numerical approximation of partial differential equations. Springer, Berlin, p 230, 2008) for the heat equation as a model of parabolic equations in any space dimension. The time discretization is performed using a uniform mesh. The scheme is formulated in a general setting in which the pairs of finite element spaces ( V h div , W h ) ⊂ H div ( Ω ) × L 2 ( Ω ) are arbitrary but satisfy the i n f - s u p hypothesis and another known condition. The discrete unknowns of the considered scheme are the set of couples ρ h n : = ( p h n , u h n ) ∈ V h div × W h which are expected to approximate ρ ( t n ) : = ( ∇ u ( t n ) , u ( t n ) ) where t n n are the mesh points of the time discretization and u is the exact solution. We justify rigorously the existence and uniqueness of the discrete solution. We then prove the new convergence results which state the optimal estimate for the error ∇ u ( t n ) - p h n n (resp. u t ( t n ) - ∂ 1 u h n n , where ∂ 1 denotes a discrete time derivative) in L ∞ ( H div ( Ω ) ) (resp. L ∞ ( L 2 ( Ω ) ) ), under assumption that the exact solution is smooth. These results are obtained thanks to some new well developed discrete a priori estimates. We applied these results to the particular case when the spaces ( V h div , W h ) are defined using the well-known Raviart–Thomas spaces of of order l ≥ 0 and we justified that the error ρ h n - ρ ( t n ) n is of order k + h l + 1 in the discrete norms of L ∞ ( H div ( Ω ) ) × W 1 , ∞ ( L 2 ( Ω ) ) under assumption that the solution u is satisfying u ∈ C 3 [ 0 , T ] ; H l + 3 ( Ω ) . The convergence results obtained in this work improve the known existing ones for PDMFEM for parabolic equations which state the convergence towards ( - ∇ u ( t ) , u ( t ) ) in the discrete norms of L ∞ ( L 2 ( Ω ) d ) × L ∞ ( L 2 ( Ω ) ) , see Johnson and Thomee (Error estimates for some mixed finite element methods for parabolic type problems. RAIRO Anal. Numér. 15(1):41–78, 1981).