Enforcing the Discrete Maximum Principle for Linear Finite Element Solutions of Second-Order Elliptic Problems

Enforcing the Discrete Maximum Principle for Linear Finite Element Solutions of Second-Order Elliptic Problems

Year:    2008

Communications in Computational Physics, Vol. 3 (2008), Iss. 4 : pp. 852–877

Abstract

The maximum principle is a basic qualitative property of the solution of second-order elliptic boundary value problems. The preservation of the qualitative characteristics, such as the maximum principle, in discrete model is one of the key requirements. It is well known that standard linear finite element solution does not satisfy maximum principle on general triangular meshes in 2D. In this paper we consider how to enforce discrete maximum principle for linear finite element solutions for the linear second-order self-adjoint elliptic equation. First approach is based on repair technique, which is a posteriori correction of the discrete solution. Second method is based on constrained optimization. Numerical tests that include anisotropic cases demonstrate how our method works for problems for which the standard finite element methods produce numerical solutions that violate the discrete maximum principle.

You do not have full access to this article.

Already a Subscriber? Sign in as an individual or via your institution

Journal Article Details

Publisher Name:    Global Science Press

Language:    English

DOI:    https://doi.org/2008-CiCP-7878

Communications in Computational Physics, Vol. 3 (2008), Iss. 4 : pp. 852–877

Published online:    2008-01

AMS Subject Headings:    Global Science Press

Copyright:    COPYRIGHT: © Global Science Press

Pages:    26

Keywords: