This course consists of three consequitively blocks. The first and second block are intended as stepping stone to the second and third block, respectively.
The course consists of two hours of lectures and four hours of computer lab sessions each week. We will use the Julia programing language julia and Jupiter notebooks to develop the assignments. We will use the gmsh mesh generation software. You are kindly requested to install Julia, the Jupiter notebook system and the Julia interface to gmsh on your laptop. You are most welcome to give a look at the references below.
Each block will be concluded with a graded homework assignment. The three homework assignments have to be defended during the oral exam.
In the first block we discuss the finite difference method to solve the Poisson equation in one and two spatial dimensions. We restrict ourselves to the interval and the square discretized by a uniform mesh. The Poisson equation models a diffusion process and typically appears in the computation of electrostatic and magneto-static fields. Our motivation is to introduce initial problem formulations and simplified numerical solution methods.
- first lab session: getting started
- second lab session: 1D Finite Difference Method
- third lab session: 2D Finite Difference Method
- wiki on the Laplace equation
- wiki on the Poisson equation: including an example on electrostatics
- wiki on boundary value problem
- wiki on electrostatics
- wiki on magnetostatics
- wiki on discrete Poisson matrix
- wiki on finite difference method
- Electromagnetic Fields and Energy by Hauss, Melcher, Zahn and Silva
In the second block we discuss the variational formulation of the Poisson equation in one spatial dimension. This variational formulation allows to subsequently solve the problem by a finite element method using non-uniform meshes. We show that the versatility of the finite element method in solving problems that are more complex than in the first block.
- fourth lab session: 1D Finite Element Method
- See Third Block
In the third and last block we discuss the finite element solution of the Poisson equation in two spatial dimension. We discuss the mesh generation using triangle and the construction of the discrete problem using a loop over all elements. We illustrate the method in the computation of magnetostatic fields in transformers and electrical machines.
- fifth lab session: using GMSH
- sixth lab session: unit square
- seventh lab session: more complex 2D geometries
- Introduction to Numerical Methods for Variational Problems by Hans Petter Langtangen and Kent-Andre Mardal. The book is freely available;
- Wolfgang Bangerth's video lectures;
- wiki Finite Element Method: Section 3 for the weak form and Section 4 for the finite element discretization;
- Comsol Multiphysics Finite Element Method: more information and illustrations;
- Comsol Multiphysics Brief Introduction to the Weak Form: good introduction to a theoretical concept that provides a basis for the finite element method;
- Ferrite Introduction to FEM
- Sphinx Finite Element Method: reference for implementation;