We deal with the numerical solution of the nonstationary heat conduction equation with mixed Dirichlet/Neumann boundary conditions. The backward Euler method is employed for the time discretization and the interior penalty discontinuous Galerkin method for the space discretization.
Assuming shape regularity, local quasi-uniformity, and transition conditions, we derive both a posteriori upper and lower error bounds. The analysis is based on the Helmholtz decomposition, the averaging interpolation operator, and on the use of cut-off functions.
Numerical experiments are presented.