On the monotonicity and discrete maximum principle of the finite difference implementation of \(C^0\)-\(Q^2\) finite element method (Q777500)

The authors solve a Dirichlet boundary value problem attached to a linear second order elliptic equation in divergence form, on a two dimensional rectangular domain. They are mainly concerned with the monotonicity and discrete maximum principle of the simplest finite difference implementation of the continuous finite element method with tensor product of quadratic polynomial basis. The scheme is fourth order accurate. The FEM discretization is based on triangular meshes under reasonable mesh constraints and the integrals in the weak formulation are computed using a tensor product of 3-point Gauss-Lobatto quadrature rule. Some numerical tests are carried out.