Comparison of continuous and discontinuous finite element methods for parabolic differential equations employing implicit time integration

This article compares the computational cost, stability, and accuracy of continuous and discontinuous Galerkin Finite Element Methods (GFEM) for various parabolic differential equations including the advection–diffusion equation, viscous Burgers’ equation, and Turing pattern formation equation system. The results show that, for implicit time integration, the continuous GFEM is typically 5–20 times less computationally expensive than the discontinuous GFEM using the same finite element mesh and element order. However, the discontinuous GFEM is significantly more stable than the continuous GFEM for advection dominated problems and is able to obtain accurate approximate solutions for cases where the classic, un-stabilized continuous GFEM fails.