On nonnegativity preservation in finite element methods for the heat equation with non-Dirichlet boundary conditions
Kapitel i bok, 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.

maximum principle,

Robin boundary condition

lumped mass method

Neumann boundary condition


Stig Larsson

Chalmers, Matematiska vetenskaper, Tillämpad matematik och statistik

Vidar Thomee

Chalmers, Matematiska vetenskaper, Tillämpad matematik och statistik

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



Matematisk analys


Grundläggande vetenskaper