Maximum principle for the finite element solution of time dependent anisotropic diffusion problems

Abstract

Preservation of the maximum principle is studied for the combination of the linear finite element method in space and the $\theta$-method in time for solving time dependent anisotropic diffusion problems. It is shown that the numerical solution satisfies a discrete maximum principle when all element angles of the mesh measured in the metric specified by the inverse of the diffusion matrix are non-obtuse and the time step size is bounded below and above by bounds proportional essentially to the square of the maximal element diameter. The lower bound requirement can be removed when a lumped mass matrix is used. In two dimensions, the mesh and time step conditions can be replaced by weaker Delaunay-type conditions. Numerical results are presented to verify the theoretical findings.

25 pages, 7 figures, 4 tables

References in corpus (2)

• Preserving Monotonicity in Anisotropic Diffusion
• Maximum principle in linear finite element approximations of anisotropic diffusion-convection-reaction problems

Cited by in corpus (8)

• The cutoff method for the numerical computation of nonnegative solutions of parabolic PDEs with application to anisotropic diffusion and lubrication-type equations
• A Study on Moving Mesh Finite Element Solution of the Porous Medium Equation
• Lumped finite element method for reaction-diffusion systems on compact surfaces
• An Efficient High-order Numerical Solver for Diffusion Equations with Strong Anisotropy
• Monotone finite difference schemes for anisotropic diffusion problems via nonnegative directional splittings
• A study on nonnegativity preservation in finite element approximation of Nagumo-type nonlinear differential equations
• Anisotropic Mesh Adaptation for Finite Element Solution of Anisotropic Porous Medium Equation
• Anisotropic mesh adaptation for 3D anisotropic diffusion problems with application to fractured reservoir simulation