Finite Element Method  Advanced Concepts  
In this page we will discuss several advanced concepts
that appear in the development of Finite Element Method. We start with
some basical concepts that we learned in our Calculus, such as contiuous
functions, derivative of a function, etc.
Continuous Function: f(x) is continuous at x_{0} if f(x) is defined at x_{0} and the limit of f(x) at x_{0} as x goes to x_{0} is f(x_{0}). Mathematically, we have
For example, the following functions are continuous functions in their domain (all x values, where the function is defined). All Polynomial Functions: f(x) = 2x + 3; g(x) = 5x^{6} + 4x^{5}  2x^{3} + 1 All Trignometric Functions: f(x) = sin(x); g(x) = cos(x); h(x) = tan(x) All Exponential Functions and Logarithmic Functions: f(x) = e^{x}; g(x) = 2^{x}; h(x) = ln(x) The following function is not continuous at x = 2 as shown in the following image, where the curve has as jump at x = 2:
Differentiable Function: f(x) is differentiable at x_{0} if f(x) is continuous at x_{0} and the following limit exits, which is called the derivative of f(x) at x_{0} and is denoted by f¢(x_{0}) or
Partial Derivatives are defined for multivariable functions. A multivariable function is a function that has two or more parameters. For example f(x, y) = x^{2} + y^{2} is a twovariable function, it represents a parabolic surface in three dimension space. When we take a derivative respect to one variable, we have a so called partial derivative. So for f(x, y) = x^{2} + y^{2} we can have partial derivative with respect to x or y variables. They are denoted as
Frequently, we also use f_{x}(x, y) and f_{y}(x, y) to denote the partial derivative against x and y variables. Example Let f(x, y) = x^{2} + 3xy + 2y^{3} + 5, find f_{x}(x, y) and f_{y}(x, y) Solution: To find f_{x}(x, y), we treat y as constant, just like the constant term 5 in f(x, y). So f_{x}(x, y) = 2x + 3y. Similarly, to find f_{y}(x, y), we treat x as constant. So f_{y}(x, y) = 3x + 6y. f(x) is said to be continuous over a closed interval [a, b] if f(x) is continuous at each point in the interval. f(x) is said to be differentiable over open an interval (a, b) if f(x) is differentiable at each point in the interval. C[a, b] is the set of all continuous over interval [a, b] C^{1}[a, b] is the set of all differential functions whose first derivative is continuous over interval [a, b]. C^{n} [a, b] is the set of all differential functions whose n^{th} derivative is continuous over interval [a, b], here n is an postive integer. C^{n}_{0} [a, b] is the set of all differential functions u whose n^{th} derivative is continuous over interval [a, b] and u(a) = u(b) = 0. Linear Space V: A linear space V is a set of elements with the addition of elements and multiplication by a scalar defined. And the addition and scalar multiplication satisfies a set of specifical properties. For example, C[a, b] , C^{1}[a, b] and C^{n}[a, b] are all linear spaces. Since from Algebra and Calculus, we know we can add two functions and mutiply a function by a real number. Also all polynomial functions forms a linear space with the usual function addition and scalar multiplication. Linearly Independent Elements: Let {v_{1}, v_{2}, v_{3} , ¼, v_{n}} be a set of elements of linear space V, we say v_{1}, v_{2}, v_{3} , ¼, v_{n} are linearly independent if
if and only if c_{1} = c_{2} = c_{3} = ¼ = c_{n} = 0. For example x^{1}, x^{2}, x^{3}, ¼, x^{n} is a set of linearly independent vectors in the linear space of all polynomials with degreen less or equal to n (denoted by P_{n} (x)). Basis of Linear Space: {v_{1}, v_{2}, v_{3} , ¼, v_{n}} Ì V is a basis of a linear space V if all element of V are of form
In the finite element method, we will construct piecewise polynomial functions as a basis for the approximation to the solutions of our problem. An Ordinary Differential Equation(ODE) is an equation that involves only one independent variable and an unknown one variable function and its derivatives. If an ODE contains upto k^{t}h order derivatives of the unknown function y(x), we say the ODE is of order k, or k^{t}h order ODE. Generally, for unknown function y(x), we can write the first order ODE equation in the form
and the second order ODE as
Example: Here are example of ODEs
Example: A boundary value problem of ODE used in beam stress analysis, here y(x) describes the deflection of a beam of length 1 with variable cross section represented by q(x). The deflection is due to the added stresses p(x) and f(x).
with boundary conditions y(0) = y(1) = 0 Functional: A functional I defined on a linear space V (Especially linear spaces of functions, like C[a, b] , C^{1}[a, b] , C^{n}[a, b] and C^{n}_{0}[a, b]) is a function whose argument is element of V and whose value is real number. We normally write as I : V ® R, where R is the set of all real numbers. Example: For the above second order ODE we can associate a functional
One can show that one way to solve this ODE is to find u_{0}(x) Î C^{1}_{0}[a, b], so that I reaches its minimum at u_{0}(x). That is I[u_{0}(x)] £ I[u(x)] for all u(x) in C^{1}_{0}[a, b]. We call u_{0}(x) a minimizer of I[u]. If the ODE has an unique solution, then the minimizer is also unique if it exists. The Finite Element Method helps us to solve the minizing
problem approximately. Instead of finding minimizer over the infinite
space C^{1}_{0}[a, b] , we search for the minimizer over
a finite dimensional subspace V, usually a linear space of all polynomials
with certain specific properties. This way, as we will see, will transform
a problem of solving differential equation into a problem of solving an
Algebraic system of linear equations as we have done in PreCalculus, Algebra,
and Linear Algebra.
