# Order of accuracy of DGFEM or FEM

I know that it is possible to determine the theoretical order of accuracy (B) of numerical solutions in FVM (for instance for a steady problem in which only the central differencing scheme (CDS) is used to interpolate, B = 2 since CDS is 2nd order accurate). I know that with FVM the order of accuracy is obtained from the Taylor series.

I interested to know whether there is a method to obtain the convergence rate for FEM solutions.

I am also particularly interested to find out if there is a way of determining the order of accuracy for DGFEM. I have read in many papers that the convergence rate is normally p+1 where p is degree of the approximating polynomial. However I am uncertain as why this is not effected by the interpolation scheme used for the flux terms.

Thank you.

• Depends on what you mean by simple. In all cases, you make use of the approximation properties of polynomials to obtain a result of the form "if you use (local) polynomials of degree at most $p$ and if(!) your exact solution is $p$ times (weakly) differentiable, then the error is of order $p$". It's a matter of taste whether Taylor expansion (for FDM) is an easier way than the Bramble-Hilbert lemma (for FEM) to derive such an estimate. (It's true that DGFEM is a bit more complicated than FEM, though, but so is FVM compared to FDM). Dec 22, 2015 at 23:40
• Or to put it another way, it's unclear what you're asking here. Could you try narrowing things down to a concrete question? Dec 22, 2015 at 23:45
• Thank you. Please forget about which one is simpler, I guess that's a relative thing. Should not have said that ... I've changes the question slightly now. Dec 23, 2015 at 12:58
• It's still not clear what answer you expect besides "Yes, of course, look at any textbook on FEM (e.g., by Braess) or DGFEM (e.g., by Di Pietro and Ern)." And of course the flux terms are part of the DGFEM discretization, so whether you obtain the optimal order or not depends on whether they are chosen appropriately (or not). However, these are usually not derived via interpolation of an independent quantity but by suitable modifications of the bilinear form, into which the (locally) polynomial approximation to the solution is inserted. Dec 23, 2015 at 13:27

• To add to Reid's comment, order $p+1$ convergence is typically expected for the $L^2$ norm. In general, if you want convergence in a Sobolev norm, you lose one order for each degree of regularity (i.e. $H^1$ convergence is usually $p$, etc...) Dec 22, 2015 at 23:56