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.

Robin boundary condition

maximum principle,

lumped mass method

Neumann boundary condition

Author

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

Subject Categories

Computational Mathematics

Mathematical Analysis

Roots

Basic sciences

DOI

10.1007/978-3-319-72456-0_35

More information

Latest update

12/28/2018