A High order Conservative Flux Optimization Finite Element Method for Diffusion Equations

TL;DR

This paper introduces a high order conservative flux optimization finite element method for elliptic diffusion equations, offering improved accuracy and efficiency through a constrained minimization approach with flux conservation.

Contribution

It develops a novel high order CFO finite element scheme that reduces degrees of freedom and enhances flux conservation compared to traditional methods.

Findings

Optimal error estimates for flux and primary variables.

Super-closeness of numerical and finite element solutions.

Numerical results demonstrate advantages in convergence and handling heterogeneous media.

Abstract

This article presents a high order conservative flux optimization (CFO) finite element method for the elliptic diffusion equations. The numerical scheme is based on the classical Galerkin finite element method enhanced by a flux approximation on the boundary of a prescribed set of arbitrary control volumes (either the finite element partition itself or dual voronoi mesh, etc). The numerical approximations can be characterized as the solution of a constrained-minimization problem with constraints given by the flux conservation equations on each control volume. The discrete linear system is a typical saddle-point problem, but with less number of degrees of freedom than the standard mixed finite element method, particularly for elements of high order. Moreover, the numerical solution of the proposed scheme is of super-closeness with the finite element solution. Error estimates of optimal order are established for the numerical flux as well as the primary variable approximations. We present several numerical studies in order to verify convergence of the CFO schemes. A simplified two-phase flow in highly heterogeneous porous media model problem will also be presented. The numerical results show obvious advantages of applying high order CFO schemes.