This preview shows page 1. Sign up to view the full content.
Unformatted text preview: on's rule weights. The diagonal
matrix D is sometimes called a \lumped" approximation of the consistent mass matrix
M. Both nite di erence and nite element solutions behave similarly for the present
problem and have the same order of accuracy at the nodes of a uniform mesh.
Example 1.3.1. Consider the nite element solution of
;u 00 +u=x which has the exact solution 0<x<1 u(0) = u(1) = 0 u(x) = x ; sinh x :
sinh 1 1.3. A Simple Finite Element Problem 17 Relative to the more general problem (1.3.1), this example has p = q = 1 and f (x) = x.
We solve it using the piecewiselinear nite element method developed in this section on
uniform meshes with spacing h = 1=N for N = 4 8 : : : 128. Before presenting results,
it is worthwhile mentioning that the load vector (1.3.15) is exact for this example. Even
though we replaced f (x) by its piecewise linear interpolant according to (1.3.14), this
introduced no error since f (x) is a linear function of x.
Letting e(x) = u(x) ; U (x) (1.3.21) denote the discretization error, in Table 1.3.1 we display the maximum error of the nite
element solution and of its r...
View
Full
Document
This document was uploaded on 03/16/2014 for the course CSCI 6860 at Rensselaer Polytechnic Institute.
 Spring '14
 JosephE.Flaherty

Click to edit the document details