×

Close

- Finite Element Methods - FEM
- Note
**B. K birla - Rtu****46 Views**- Uploaded 11 months ago

Simulationsmethoden I WS 09 10 Lecture Notes Finite element methods applied to solve PDE Joan J. Cerdà ∗ December 14, 2009 ICP, Stuttgart Contents 1 In this lecture we will talk about 2 2 FDM vs FEM 2 3 Perspective: different ways of solving approximately a PDE. 2 4 Basic steps of any FEM intended to solve PDEs. 4 5 5 FEM in 1-D: heat equation for a cylindrical rod. 6 FEM in 2-D: the Poisson equation. 10 7 Sparse Matrixes (band matrixes) and FEM. 17 8 18 Other tricks for FEM and beyond. 9 Bibliography ∗ 19 jcerda put here an at symbol icp.uni-stuttgart.de 1

1 In this lecture we will talk about 1. FDM vs FEM. 2. Perspective: different ways of solving approximately a PDE. 3. Basic steps of any FEM intended to solve PDEs. 4. FEM in 1-D: heat equation for a cylindrical rod. 5. FEM in 2-D: the Poisson equation. 6. Sparse Matrixes (band matrixes) and FEM. 7. Other tricks for FEM and beyond. 8. Bibliography. 2 FDM vs FEM 1. With FEM you can handle better boundary problems with odd geometries. 2. With FEM it is easier to make more finer subdivisions of the space in regions where you need more accuracy. This is not so easy to do with FDM. A word of caution: FEM as FDM are suitable for linear PDE’s. If you have non-linear PDEs. You will have first to linearize it. 3 Perspective: different ways of solving approximately a PDE. I have a PDE with certain bc (boundary conditions) to be solved, which options do I have: 1. Analytical solution: the best, but not always available. 2. Approximate solutions. There is always an error. The difference between the exact value and the approximate value is called the residual, we will denote it by R. For instance, let’s suppose I have to solve the following ODE: dφˆ ˆ + φ(x) = 0 dx (1) if my approximate solution is φ(x) = 1 + a1 x + a2 x2 , where a1 and a2 must be determined to get the φ as close to the exact solution φˆ as I can. In that example the residual is R(x; a1 , a2 ) ≡ dφ + φ(x) = 1 + (1 + x) a1 + (2x + x2 ) a2 dx 2 (2)

There are several ways of solving approximately a PDE, the most usual are: 1. Ritz method (aka Goodman’s method): we get a1 and a2 via asking Z R(x; a1 , a2 )dx = 0 (3) Domain 2. Variational method (Rayleight-Ritz method): from the PDE we get its variational integral. The function φˆ that minimizes that variational integral is the solution to the problem. [FALTA] 3. Weighted residual methods: in this methods we assume the parameters a1 , a2 (, ..., an if we have n in a general case) to be determined by asking the residual to obey a set of equations: Z R(x; a1 , .., an ) ωi (x) dx = 0 (4) Domain where ωi (x) with i = 1, ..n are n arbitrary functions called weighting functions. Some usual choices of ωi (x) are a) Collocation method: we ask the residual R to be zero at n points {xi }i=n i=1 inside the domain. which is equivalent to say ωi (x) = δ(x − xi ). b) Sub-domain method: we split the domain in n regions and we ask the averaged R, < R >, to be zero inside each region. ∂R c) Least-Squares method: we set ωi (x) = ∂a i d) Galerkin method: if we use approximate solutions of the type φ(x) = 1 + a1 x + a2 x2 , this means that Ni (x) = x and N2 (x) = x2 are a set of basis vectors for all our possible solutions. We set ωi (x) = Ni (x),i.e., we have to solve n equations (in our example n = 2) like Z R(x; a1 , .., an ) Ni (x) dx = 0. (5) Domain 4. FDM: Finite Difference Methods. 5. FEM: Finite Element Methods. 6. FVM: Finite Volume Methods. A refined FDM popular in Computational fluid dynamics. 7. MOM:method of moments. You convert your differential equation into an integral equation. Used specially in electromagnetics. 3

4 Basic steps of any FEM intended to solve PDEs. In all FEM variants there are always the same sequence of steps to be taken 1. Discretize the continuum: divide the solution into smaller regions that we call elements. The elements contain inside a certain number of points we call nodes. There are lots of shapes the elements can have. From segments of lines, triangles, squares, etc, to curved elements. The one/ones you will use depends on the problem you want to solve. For instance, for a 1D problem, like a cylindrical rod with radial symmetry, the most simple is to take the elements as linear segments with two nodes per segment (see Fig.1) and discretize it is as shown in figure 2. For a 2D problem, the most simple can be using triangular elements (see figure 3). 2. Select the type of trial function to use, and in turn the shape functions: We select what kind of functions we will take to describe the variation of the function φ inside each element (the trial function). This is equivalent to say, that we select the basis set of functions that will describe our solution. One of the usual choices is to take a polynomial like for instance φ(x) = a0 + a1 x (known as linear element) or φ(x) = a0 + a1 x + a2 x2 (known as quadratic element). If we have n unknown coefficients a0 , a1 , ..., an−1 we will need the element to have n nodes to be able to determine them. Is it easy to show that the value of our trial function φ for a given position x inside an element can be written as a function of the values of φ at the N nodes the of element ([φ1 , φ2 , ..., φN ]) , i.e. φ(x) = [φ] · [N] = [φ1 , φ2 , ..., φN ] · [N1 (x), N2 (x), ..., NN (x)] (6) where the {Ni (x)}i=N i=1 are known as the Shape Functions. 3. The formulation: Given the PDE you want to solve, now you must find a system of algebraic equations for each element "e" such that by solving it you got the values of of φ at the position of nodes of the element "e" ([φ1 , φ2 , ..., φN ] ≡ [φ]e ), i.e., you must find for each element "e" the matrix [K]e and the vector [f ]e such that, [K]e · [φ]e = [f ]e . (7) This is one of the more tricky parts of the FEM. There are different ways of getting the matrix [K]e , we will see it later. 4. Assembling the equations for different elements: one has to assemble the equations for all elements. The "problem" is that usually contiguous elements have nodes in common (see for instance figure 2, and therefore two algebraic equations may refer to the same node, and that has to be taken into account. At the end, if we have a total of M effective nodes1 in the system, then we must build up a global matrix [K] of size M × M and a 1 For instance in figure 1 we have a total of 5 elements that have each one two nodes, nonetheless due to the nodes in common, the number of effective nodes we have (the number of unknowns) is just M = 6 and not 10. 4

## Leave your Comments