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

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. AMS subject classifications: 35J25, 65N99

[1]  Randolph E. Bank,et al.  PLTMG - a software package for solving elliptic partial differential equations: users' guide 8.0 , 1998, Software, environments, tools.

[2]  P. G. Ciarlet,et al.  Basic error estimates for elliptic problems , 1991 .

[3]  Sergey Korotov,et al.  Discrete maximum principles for finite element solutions of nonlinear elliptic problems with mixed boundary conditions , 2005, Numerische Mathematik.

[4]  Louis J. Durlofsky,et al.  Unstructured grid optimization for improved monotonicity of discrete solutions of elliptic equations with highly anisotropic coefficients , 2006, J. Comput. Phys..

[5]  M. Shashkov,et al.  The Numerical Solution of Diffusion Problems in Strongly Heterogeneous Non-isotropic Materials , 1997 .

[6]  M. Powell.  On the quadratic programming algorithm of Goldfarb and Idnani , 1985 .

[7]  Todd F. Dupont,et al.  Failure of the discrete maximum principle for an elliptic finite element problem , 2004, Math. Comput..

[8]  Donald Goldfarb,et al.  A numerically stable dual method for solving strictly convex quadratic programs , 1983, Math. Program..

[9]  Michal Křížek,et al.  Finite element approximation of a nonlinear heat conduction problem in anisotropic media , 1998 .

[10]  Sergey Korotov,et al.  Discrete maximum principles for finite element solutions of some mixed nonlinear elliptic problems using quadratures , 2006 .

[11]  Robert Eymard,et al.  A mixed finite volume scheme for anisotropic diffusion problems on any grid , 2006, Numerische Mathematik.

[12]  P. G. Ciarlet,et al.  Maximum principle and uniform convergence for the finite element method , 1973 .

[13]  Tomás Vejchodský,et al.  Discrete maximum principle for a 1D problem with piecewise-constant coefficients solved by hp-FEM , 2007, J. Num. Math..

[14]  Mikhail Shashkov,et al.  The repair paradigm and application to conservation laws , 2004 .

[15]  M. Shashkov.  Conservative Finite-Difference Methods on General Grids , 1996 .

[16]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[17]  Raphaël Loubère,et al.  The repair paradigm: New algorithms and applications to compressible flow , 2006 .

[18]  K. Schittkowski,et al.  QL : A Fortran Code for Convex Quadratic Programming-User ’ s Guide , 2006 .

[19]  Ivar Aavatsmark,et al.  Monotonicity of control volume methods , 2007, Numerische Mathematik.

[20]  Alexandre Ern,et al.  Discrete maximum principle for Galerkin approximations of the Laplace operator on arbitrary meshes , 2004 .

[21]  J. Strikwerda.  Finite Difference Schemes and Partial Differential Equations , 1989 .

[22]  J. Erhel,et al.  The maximum principle violations of the mixed‐hybrid finite‐element method applied to diffusion equations , 2002 .

[23]  Daniil Svyatskiy,et al.  Monotone finite volume schemes for diffusion equations on unstructured triangular and shape-regular polygonal meshes , 2007, J. Comput. Phys..

[24]  Michal Křížek,et al.  On the maximum and comparison principles for a steady‐state nonlinear heat conduction problem , 2003 .

[25]  Kazufumi Ito,et al.  Three-dimensional elliptic solvers for interface problems and applications , 2003 .