AFEM01
AFEM01
Justify your answer and eventually compute the weak derivative. Exercise 2. Build the stiness matrix and the load vector for the linear nite element discretization of the Sturm-Liouville problem pu with p(x) = 1 + x2 , q(x) = 6, f (x) = 8 + 16x corresponding to a partition T consisting of 3 intervals of equal length and compute the nite element solution uT . Determine the values uT (0.25) and uT (0.375). Exercise 3. Decide which of the following functions are contained in H 1 () u(x, y) = 1 + x2 y 3 + x3 y, w(x, y) = 9x2 y xy 4 v(x, y) =
4
+ qu = f in (0, 1),
u(0) = u(1) = 0
Here, denotes the interior of a circle with radius 1 centred at the origin. Justify your answer and eventually compute the rst order weak derivatives. Exercise 4. Give the variational formulation of the dierential equation 2u 2u 2u 2 +2 3 2 + 2(x2 + y 2 )u = f in , x xy y u = 0 on . Exercise 5. Describe the basic steps of the nite element discretization of an elliptic dierential equation. Exercise 6. How does one handle Neumann boundary conditions in the nite element discretization of elliptic dierential equations?
1
Exercise 7. Build the system of linear equations corresponding to a linear nite element discretization of the dierential equation u = xy in u = 0 on in the unit square = (x, y) : 0 < x, y < 1 and compute the discrete solution uT . Here, the partition T is a Courant triangulation consisting of right angled isosceles triangles with 1 short sides of length h = 3 . Compute the values uT (0.25, 0.25) and uT (0.375, 0.5). Exercise 8. Write a program that solves the linear nite element discretization of the dierential equation u = f in u = 0 on in the unit square = (x, y) : 0 < x, y < 1 using a Courant triangulation consisting of right angled isosceles tri1 angles with short side of length h = n , n 2. Integrals should be evaluated using a suitable quadrature formula. Exercise 9. Compute the nodal shape functions z,1 of the triangle K with vertices (1, 1), (5, 2) and (3, 4) (1) using a transformation to the reference triangle, (2) using the geometrical data of K. Exercise 10. Compute the nodal shape functions z,1 of the parallelogram K with vertices (1, 1), (5, 2), (7, 5) und (3, 4) (1) using a transformation to the unit square, (2) using the geometrical data of K.