Shape Functions and Basis Functions

In scientific computing the main aim of the tesselation is to provide a number of pairwise disjoint sets. Shape functions are transformed to single archetypical elements of the tesselation and basis functions for the function space are obtained. Each of these basis functions is non-zero on the given set (local support) whereas it is zero otherwise. The shape of these functions is typically the same for all sets of the tesselation.

Applying this method on all elements of the tesselation, a set of basis functions can be found based on the underlying tesselation of the simulation domain. In the following the construction of base functions is performed in detail for linearly bounded triangular elements. It is assumed that the domain is divided into subdomains of the same geometrical archetype, namely triangles. In this case it is possible to use standard functions which are based on the geometrical archetype of the elements. The single elements are transformed into the the standard element in order to evaluate the function value. For a triangle $ \mathcal{I}$ a sensible choice are local linear shape functions $ f^{\mathrm{l}}$ as basis functions is (Fig. 2.2):

$\displaystyle \mathcal{I}: 0 \leq x \leq 1, \; 0 \leq y \leq 1-x.$ (2.7)
$\displaystyle f^{\mathrm{l}}_1(x, y) = x$ (2.8)
$\displaystyle f^{\mathrm{l}}_2(x, y) = y$ (2.9)
$\displaystyle f^{\mathrm{l}}_3(x, y) = 1-x-y$ (2.10)

Figure: Shape of the base functions of $ \mathcal{F}$ .
\includegraphics[width=14cm]{DRAWINGS/spaceF.eps}

It can be easily seen, that each function yields unity in one corner point and zero in all other corner points. This implies that the weighting coefficients of the single functions have the value of the function when evaluated the respective corner point.

These functions are applied to the standard triangle and this triangle is back-transformed in its original position. Of course, such a transformation can be applied on many different archetypes of elements and also the local shape functions can be altered. From this point of view the shape of these functions $ f^{\mathrm{l}}_i$ can be chosen arbitrarily.

For the triangle example, the three shape functions are transformed to each triangular element. One obtains a function space with linear basis functions defined on local elements of the tesselation. However, it has to be mentioned that for such a function space even very essential properties such as continuity are not given.

In order to guarantee continuity it is necessary to use a subspace of the original function space defined by the local shape functions. The main aim is to search for a subspace in which the base function and therewith all functions are continuous. To form such a function space, some weighting coefficients have to be coupled in a predefined way. In this case it turns out that the collocation of functions and points as described above is convenient, because in the corner points, where several triangles intersect, the function values have to be identical in all triangles. For this reason the function space is restricted so that each shape function which is collocated with the same point has the same weighting factor. Of course, the coupling of the weighting factors restricts the possible functions of the function space and forms a subspace $ \mathcal{F}'$ of the initial space $ \mathcal{F}$ . The basis functions $ f'$ of the newly defined subspace $ \mathcal{F}'$ of the original base space can be written in the following manner, where the points of the intersection of the triangles are denoted as $ \mathcal{Q}_i$ .

$\displaystyle f'_i(\mathcal{P}) = \sum_{j : f_j(\mathcal{Q}_j) \neq 0} \frac{f_j(\mathcal{P})}{f_j(\mathcal{Q}_j)}$ (2.11)

If the values of the shape functions in the collocation points are identically unity, the definition of $ \mathcal{F}'$ can be simplified to (see Fig. 2. 3 )

$\displaystyle f'_i(\mathcal{P}) = \sum_{j : f_j(\mathcal{Q}_i) \neq 0} f_j(\mathcal{P})$ (2.12)

The resulting functions $ f'_i$ have a different local support than the base functions of $ \mathcal{F}$ . These functions are defined within a neighborhood of the respective intersection point as can be seen in Fig. 2.3

Figure: Local support and shape of the base functions of $ \mathcal{F}'$ .
\includegraphics[width=14cm]{DRAWINGS/spaceFprime.eps}

In some cases it is favorable to use more shape functions on a base element. An example could be the use of third order shape functions which form the function space $ \mathcal{G}$ . Such a function space can be based on $ 10$ basis functions, where three are collocated with the corner points of the triangle, $ 6$ are collocated with edge points, and one is collocated with a point the interior of the triangle.

It can be seen that the function space $ \mathcal{F}$ is a subspace of the function space $ \mathcal{G}$ . In analogy to the space $ \mathcal{F}$ , the property of continuity is not given in this case, because each function is based on a single triangle. In this case, however, the fitting of the functions is not restricted to a single point. While it is clear that continuity is given on the edge between two triangles, where the two end points of the edge are fit to the same value, this does not hold true for higher-order polynomial shape functions. For this reason, it is necessary to guarantee that continuity is preserved on the edge. Moreover, it can be stated that if more shape functions are available more restrictions on continuity as well as smoothness can be imposed by reducing the number of base functions.

Again, the collocation of functions and points on the triangle implies an easy fitting of neighboring triangles. In this case the collocation fitting is not only performed on the corner points of the triangle but also on two points of each edge. From this point of view it is clear that even though the use of collocation of points and shape functions is not coercive, other local shape functions imply further equations in the final equation system and, therefore, an unnecessary use of memory and computation time.

In analogy to the measures which have been taken to derive the global shape functions of the local linear shape functions the third order local shape functions yield the global shape functions $ g'_i$ of the function space $ \mathcal{G}'$ by

$\displaystyle g'_i(\mathcal{P}) = \sum_{j : g_j(\mathcal{Q}_i) \neq 0} g_j(\mathcal{P}) \; .$ (2.13)

It has to be mentioned that in this case the collocation points $ \mathcal{Q}_j$ are the corner points, points on the edges and points in the interior of the triangles. The functions which are collocated with the interior of the triangle have their local support only on the respective triangle. The functions which are collocated with an edge point of a triangle have their local support on both triangles which are supersets of the common edge. As already mentioned, the functions which are collocated with the corner points of the triangle have their local support on all triangles which cover the respective point.

It should be noted that the method which has been shown for triangles can also be used on many other archetypes such as tetrahedra, squares, prisms, and pyramides as long as a transformation into a standard element and local element functions are given.

Michael 2008-01-16