On nonnegativity preservation in finite element methods for the heat equation with non-Dirichlet boundary conditions
Book chapter, 2018

By the maximum principle the solution of the homogeneous heat equation with homogeneous Dirichlet boundary conditions is nonnegative for positive time if the initial values are nonnegative. In recent work it has been shown that this does not hold for the standard spatially discrete and fully discrete piecewise linear finite element methods. However, for the corresponding semidiscrete and Backward Euler Lumped Mass methods, nonnegativity of initial data is preserved, provided the un- derlying triangulation is of Delaunay type. In this paper, we study the correspond- ing problems where the homogeneous Dirichlet boundary conditions are replaced by Neumann and Robin boundary conditions, and show similar results, sometimes requiring more refined technical arguments.

Neumann boundary condition

maximum principle,

Robin boundary condition

lumped mass method


Stig Larsson

Chalmers, Mathematical Sciences, Applied Mathematics and Statistics

Vidar Thomee

Chalmers, Mathematical Sciences, Applied Mathematics and Statistics

Contemporary Computational Mathematics - a celebration of the 80th birthday of Ian Sloan

978-3-030-10203-6 (ISBN)

Subject Categories

Computational Mathematics

Mathematical Analysis


Basic sciences



More information

Latest update