In this paper, we develop a gradient recovery based linear (GRBL) finite element method (FEM) and a Hessian recovery based linear (HRBL) FEM for second order elliptic equations in non-divergence form. The elliptic equation is casted into a symmetric non-divergence weak formulation, in which second order derivatives of the unknown function are involved. We use gradient and Hessian recovery operators to calculate the second order derivatives of linear finite element approximations. Although, thanks to low degrees of freedom (DOF) of linear elements, the implementation of the proposed schemes is easy and straightforward, the performances of the methods are competitive. The unique solvability and the $H^2$ seminorm error estimate of the GRBL scheme are rigorously proved. Optimal error estimates in both the $L^2$ norm and the $H^1$ seminorm have been proved when the coefficient is diagonal, which have been confirmed by numerical experiments. Superconvergence in errors has also been observed. Moreover, our methods can handle computational domains with curved boundaries without loss of accuracy from approximation of boundaries. Finally, the proposed numerical methods have been successfully applied to solve fully nonlinear Monge-Amp\`{e}re equations.
翻译:在本文中,我们开发了一种基于梯度回收的线性线性(GRBL)有限元素法(FEM)和一种基于黑森回收的线性线性(HRBL)线性线性(FEM)法(FEM),用于非diverence形式的第二顺序椭圆方程式。将椭圆方程式投入一种对称非diverence弱度配方,其中涉及未知函数的第二顺序衍生物。我们使用梯度和海森回收操作员计算线性元素近似线性线性元素的第二顺序衍生物。虽然由于线性元素的自由度较低,拟议的计划的实施既简单又直截了当,但方法的性能是竞争性的。此外,我们提出的方法可以成功地处理计算域,而不会导致MONUBL公式的精确度;最后,我们提出的方法可以顺利地使用数字曲线边界的精确度,最终可以解决MONBA的精确度。