Introduction to finite elements/Choosing weight function

From Wikiversity
Jump to navigation Jump to search

Choice of weighting functions[edit]

As you have seen, we need to know the weighting functions (also called test functions) in order to define the weak (or variational) statement of the problem more precisely.

Consider the second model problem (3). If is a weighting function, then the weak form of equation (3) is

Let the possible weighting functions for this equation be such that they satisfy the boundary conditions, that is and . This set of weighting functions will be denoted by the symbol .

The exact weak form of the problem can then be stated as:


To get to an approximate weak form, we need to choose trial functions. Let these trial functions () reside in the function space . The approximate weak form of the problem is


One requirement that these trial functions () have to satisfy is that

The trial functions also have to satisfy the boundary conditions.

If you look at the weak form shown above, you will notice that the highest derivatives of are of the second order while there are no derivatives of . This leads to a lack of symmetry in the formulation that is usually avoided for numerical reasons.

To get a symmetric weak form, we integrate the first term inside the integral by parts to get

(Recall that the formula for integration by parts is .)

Applying the boundary conditions on , we get

Therefore the problem can be restated in symmetric form as :


Since the same order of derivatives appear in both the trial and the weighting functions, we can choose our functions from the same set , where . As you can see, the smoothness requirements on are further weakened. This boundary value problem is also called the variational boundary value problem.

The variational boundary value problem (18) can be (and often is) written using bilinear function notation as


where

Characteristics of [edit]

Let us now characterize , the class of admissible functions for the problem. We know that if and are irregular then their derivatives are even more irregular. Therefore, the most irregular term in the problem statement is

Now, both and are functions in the space . So we can choose to be equal to . Hence the minimum requirement for is

Functions like the above are called square integrable.

Thus, our weighting functions have to be square integrable. This is a stricter condition than just requiring to be piecewise continuous, and reduces the number of functions that you can choose from.

There are two other interesting properties of the space . These are :

  • It is a linear space which means that if and are two arbitrary functions in this space, then the following operations hold :
  • It is an infinite dimensional space. This means that an infinity of parameters are needed to uniquely determine a weighting function . Thus,
where, are linearly independent functions (also called basis functions) and are constants. If we choose only of these basis functions then this is a -dimensional (i.e., finite dimensional) subspace of and we denote this as .