Abstract:We extend the finite element method introduced by Lakkis and Pryer [2011] to approximate the solution of second order elliptic problems in nonvariational form to incorporate the discontinuous Galerkin (DG) framework. This is done by viewing the NVFEM as a mixed method whereby the finite element Hessian is an auxiliary variable in the formulation. Representing the finite element Hessian in a discontinuous setting yields a linear system of the same size and having the same sparsity pattern of the compact DG methods for variational elliptic problems. Furthermore, the system matrix is very easy to assemble, Thus this approach greatly reduces the computational complexity of the discretisation compared to the continuous approach.
We conduct a stability and consistency analysis making use of the unified framework set out in Arnold et. al. [2001]. We also give an apriori analysis of the method. The analysis applies to any consistent representation of the finite element Hessian, thus is applicable to the previous works making use of continuous Galerkin approximations.
What problem does this paper attempt to address?
This paper attempts to solve the numerical solution of non - variational problems, especially the discretization method of second - order elliptic partial differential equations (PDEs) in non - variational form. Specifically, the main objective of the paper is to apply the Discontinuous Galerkin (DG) method to non - variational problems to improve computational efficiency and numerical stability.
### Main Problem Background
1. **Non - variational Problems**:
Non - variational partial differential equations are those that cannot be derived from variational principles. For example, equations of the form
\[
-A : D^2 u = f
\]
where \( A \) is a matrix, \( D^2 u \) is the Hessian matrix of \( u \), and \( f \) is the given source term. In some cases, standard finite element methods cannot be directly applied to such problems because they do not have a corresponding variational structure.
2. **Limitations of Traditional Methods**:
For non - variational problems, traditional continuous finite element methods may encounter difficulties, especially when dealing with coefficient matrices \( A \) with singular derivatives, which may lead to numerical instability or high computational complexity.
### The Method Proposed in the Paper
The paper proposes a new method, that is, to solve non - variational problems within the Discontinuous Galerkin framework. The specific steps are as follows:
1. **Introducing the Discrete Hessian**:
Discretize the Hessian matrix \( D^2 u \) in the Discontinuous Galerkin framework and denote it as \( H[u] \). This allows the problem to be transformed into a mixed - type problem, where the Hessian matrix serves as an auxiliary variable.
2. **Constructing the DG Scheme**:
Utilize a unified DG framework (such as the framework proposed by Arnold et al. [2001]), and construct a DG discretization scheme suitable for non - variational problems. This scheme not only maintains the same sparsity and matrix size as the continuous finite element method but also greatly simplifies the system assembly process.
3. **Stability and Consistency Analysis**:
Through strict mathematical analysis, it is proved that the proposed DG method is stable and consistent under specific conditions. In particular, for coercive non - variational problems, variational techniques can be used for analysis, and the optimal convergence rate is proved.
### Main Contributions
- **Improving Computational Efficiency**: Compared with the continuous finite element method, the DG method has higher computational efficiency when dealing with non - variational problems, especially in large - scale parallel computing environments.
- **Expanding the Scope of Application**: This method is not only applicable to non - variational problems but can also be extended to the study of more general fully nonlinear PDEs.
- **Theoretical Guarantee**: Provides detailed stability and convergence analysis, providing a solid theoretical foundation for practical applications.
In conclusion, this paper provides an efficient and stable numerical solution scheme for non - variational problems by introducing the Discontinuous Galerkin method, and broadens the numerical research field of non - variational PDEs.