The hypercircle method for the error estimation of FEM solutions




This project describes how to obtain explicit a priori error estimation for FEM solution by using the hypercircle method.



  • DEMO

Explicit a priori error estimation for FEM solutions

Xuefeng LIU (Last updated: 2020/05/13)

Problem description

We consider the homogeneous boundary value problem of the Poisson equation over a bounded 2D domain.

$$ -\Delta u = f \quad \mbox{ in } \Omega; \quad u=0 \mbox{ on } \partial\Omega. $$

Let $u_h$ be the linear conforming FEM solution to the above problem. We provide an a priori error estimation to $u_h$ with an explicit value of $C_h$ such that $$ |\nabla(u-u_h)| \le C_h |f| \quad \forall f \in L^2(\Omega)\:. $$

The hypercircle method is adopted here. See detail in https://doi.org/10.1137/120878446 .


  • create_mat.ipynb: The python code to create the matrices.
  • a_priori_error_estimation.ipynb: Calculate the quantity $\kappa_h$ and the explicit value (upper bound) of $C_h$. Language: Octave.


  • Python3 + FEniCS for create_matrix.ipynb .
  • Matlab or Octave for a_priori_error_estimation.ipynb .

Setting of computing

In create_matatrix.ipynb, set the mesh size and degree of finite element method.

Here is a sample setting. It defines an uniform mesh with subdivision number 16 for the unit square domain. The FEM is selected as the linear conforming one.

N = 16
mesh = UnitSquareMesh(N,N) 
degree = 1

About the directory

Folders or files beginning with a dot are not displayed by default.

Virtual Machine Setting


You are starting the virtual machine as a visitor to current project. As a visitor, you can change files in the booted virtual machine, but the changed files will be aborted when the server is shut down.

(Please login first to start the virtual machine.)

About Machine Type

The machine with type as "n1-standard-1" has 1 CPU Core and 4GB memory. The Google app compute engine provides a detailed guide of the machine type. For more detailed information, please refer to More detail.
If you need a high-spec machine type, please contact the site manager.