On positivity preservation in some finite element methods for the heat equation

Paper in proceeding, 2015

We consider the initial boundary value problem for the homogeneous heat equation, with homogeneous Dirichlet boundary conditions. By the maximum principle the solution is nonnegative for positive time if the initial data are nonnegative. We study to what extent this property carries over to some piecewise linear finite element discretizations, namely the Standard Galerkin method, the Lumped Mass method, and the Finite Volume Element method. We address both spatially semidiscrete and fully discrete methods.