On Local and Global Bisection-Type Mesh Refinements in C Programming Language

The finite element method (FEM) is a numerical method of resolution of many problems modeled in terms of partial differential equations. It is a powerful and widely used method in science and engineering. Its simulation requires in a first phase the construction of a mesh of the computational domain. In our paper, we present an efficient approach for local and global mesh refinements of triangular and quadrilateral two-dimensional meshes in C programming language. The proposed refinement algorithms are based on a bisection-type method which produces nested refinements of the triangulation. The algorithms are short, easy to understand and modify, moreover they can be easily integrated in a FEM C program.