Course detail
Numerical Methods III
FSI-SN3-A Acad. year: 2023/2024 Winter semester
The course gives an introduction to the finite element method as a general computational method for solving differential equations approximately. Throughout the course we discuss both the mathematical foundations of the finite element method and the implementation of the involved algorithms.
The focus is on underlying mathematical principles, such as variational formulations of differential equations, Galerkin finite element method and its error analysis. Various types of finite elements are introduced.
Language of instruction
English
Number of ECTS credits
4
Supervisor
Department
Entry knowledge
Differential and integral calculus for multivariable functions. Fundamentals of functional analysis. Partial differential equations. Numerical methods, especially interpolation, quadrature and solution of systems of ODE. Programming in MATLAB.
Rules for evaluation and completion of the course
Graded course-unit credit is awarded on the following conditions: Active participation in practicals and elaboration of assignments. Participation in the lessons may be reflected in the final mark.
If we measure the classification success in percentage points, then the grades are: A (excellent): 100--90, B (very good): 89--80, C (good): 79--70, D (satisfactory): 69--60, E (sufficient): 59--50, F (failed): 49--0.
Attendance at lectures is recommended, attendance at seminars is required. Absence from lessons may be compensated by the agreement with the teacher supervising the seminars.
Aims
The aim of the course is to acquaint students with the mathematical principles of the finite element method and an understanding of algorithmization and standard programming techniques used in its implementation.
In the course Numerical Methods III, students will be made familiar with the finite element method and its mathematical foundations and use this knowledge in several individual projects.
The study programmes with the given course
Programme N-MAI-A: Mathematical Engineering, Master's, compulsory
Type of course unit
Lecture
26 hours, optionally
Teacher / Lecturer
Syllabus
The first four lectures will be devoted to the explanation of the algorithm for solution of the model problem of type "stationary heat conduction" in a plane polygonal domain using linear triangular finite elements. This enables students from the very beginning of practicals to start experimenting with code programming. Only the following lectures will concentrate on the mathematical theory of finite elements.
1. Classical and variational formulation, triangulation, piecewise linear functions.
2. Discrete variational formulation, elementary matrices and vectors.
3. Elementary matrices and vectors – continuation.
4. Assembly of global system of equations, its solution, postprocessing.
5. Selected pieces of knowledge of functional analysis. The space W^k_2.
6. Traces of functions from the space W^k_2. Friedrich's and Poincare's inequality.
7. Bramble-Hilbert's lemma. Sobolev's imbedding theorem.
8. Formal equivalence of the elliptic boundary value problem and the related variational problem.
9. Finite element spaces of Lagrange's type. Definition of approximate solution. Existence and uniqueness theorem.
10. Transformation of a general triangle onto the reference triangle. Relations between norms on the general triangle and on the reference triangle.
11. Interpolation theorem.
12. Numerical integration.
13. Adaptivity in FEM.
Computer-assisted exercise
13 hours, compulsory
Teacher / Lecturer
Syllabus
Practicals will take place in a computer lab with the support of the MATLAB and Visual Studio. The algorithm for the elliptic problem will be explained during the first four lessons. The algorithms for the parabolic, hyperbolic and eigenvalue problems will be explained in brief on practicals. It is supposed that students will work individually with lecture notes (containing detailed descriptions of algorithms). Students are also expected to create and debug individually their own MATLAB programs.
1-2. Programming tools, first introduction.
3-4. Further details, preparation for writing of the program for solution of an elliptic problem (stationary heat conduction).
5-6. Developing the program for an elliptic problem. Explanation of the algorithm for the solution of the parabolic problem (nonstationary heat conduction).
7-8. Developing the program for a parabolic problem. Explanation of the algorithm for the solution of the hyperbolic problem (membrane vibrations).
9-10. Developing the program for a hyperbolic problem. Explanation of the algorithm for the solution of the eigenvalue problem.
11-12. Developing the program for an eigenvalue problem.
13. Teacher's reserve.