Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
|
|
- Albert Stone
- 6 years ago
- Views:
Transcription
1 Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman
2 Contents of the course Fundamentals of functional analysis Abstract formulation FEM Spatial (meshing) and functional discretization (the basis functions) Convergence, regularity Variational crimes Implementation Mixed formulations (e.g. Stokes) Stabilisation for flow problems Error indicators/estimation Adaptivity FEM and its convergence Dr. Noemi Friedman PDE2 Seite 2
3 Abstract formulation, examples FEM and piecewise polynomials Non-degenerate triangulation, refinement of triangulation Convergence using piecewise linear functions Convergence using piecewise higher order elements Implementation, sparsity of the stiffnes matrix Variational crime: numerical integration Variational crime: curved boundaries FEM and its convergence Dr. Noemi Friedman PDE2 Seite 3
4 Recap: boundedness of the error of Galerkin method a u u h, v = 0 v V h If a, is an inner product, it means that the approximation is an orthogonal projection to the subspace in the energy norm: error = u x u h x E u x v x E v x V h According to Céa s theorem (see prove at the lecture note), even without a, being symmetric, the error of the approximation of Galerkin will be allways bounded: u u h M δ u v v V h Where M and δ are constants from the conditions of boundedness and V-ellipticity of the bilinear term a, : a u, v M u v a u, u δ u 2 FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 4
5 Recap: what do we have to solve? [Chapter 4] N i=1 c i Ψ i (x) Ψ j x dω = Ω f x Ψ j x dω Ω Kc = f K ij f j FEM: Galerkin method where Ψ j are piecewise polynomials Main goals when implementing: efficient calculation of K efficient calculation of f solve Kc=f efficiently true solution is approximated well (error is small enough) piecewise polynomials FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 5
6 Piecewise polynomials and the FEM [Chapter 4.1] piecewise polynomials: a function that is defined by a polynomial on each subdomain mesh: the collection of subdomains nodal basis: 1 i = j ψ i x j = 0 i j ψ i x j, y j = 1 i = j 0 i j ψ 2 ψ 3 ψ 4 ψ 5 u x u h x = 4 i=1 u i ψ i x u h (x) P h 1 continious piecewise linear functions for Poisson equation we have to show that P h 1 H 1 when it s not satisfied it is not a conforming element method FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 6
7 Recap: Why linear piecewise polynomials defined on a [Chapter 4.1.1] triangulation are in the subspace Suppose, that uniqueness and existence can be shown in H 1 Show, that the space of all continious piecewise linear functions defined on the triangulation T h is a subspace of H 1. The continious piecewise linear functions can be writen in the general form: u = a i + b i x + c i y (x, y) T i a j + b j x + c j y = a i + b i x + c i y x, y e = T i T j, We have to show that all us are in H 1 : Ω u 2 < Ω = i T i Ω u x 2 + u y 2 < FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 7
8 Recap: Why linear piecewise polynomials defined on a triangulation are in the subspace 1) First let s show that u 2 < Ω u 2 = Ω i T i a i + b i x + c i y 2 < i T i max a i + b i x + c i y 2 T i < 2.)Show that Ω u 2 < Ω = i T i There is no derivatives in a strong sense, but let s see whether these give weak derivatives: g 1 : = u x = b i (x, y) int(t i ) g 2 : = u y = c i (x, y) int(t i ) FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 8
9 Recap: Why linear piecewise polynomials defined on a triangulation are in the subspace are weak derivatives if for any v C 0 Ω u v x = Ω vg 1 and Ω u v y = Ω vg 2 holds. We prove only this, because the approach of the derivation for g 2 is same. Let s first look at the left hand side of the equation: Ω u v x = i T i u v x = i T i uv n 1 + T i u x v uv n 1 = 0 on the edges of the boundary, and uv n 1 = T i T i On the common edges of neighboring triangles T i and T j T j uv n 1 i T i uv n 1 = 0 T i u x v = T i b i v FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 9
10 Recap: Why linear piecewise polynomials defined on a triangulation are in the subspace The left hand side of the equation: Ω u v x = vg 1 Ω Ω u v x = i T i b i v The left hand side of the equation: Ω vg 1 = i T i b i v FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 10
11 Piecewise polynomials and the FEM [Chapter 4.1] Derivatives in the classical sense: weak derivatives of v (see proof in Gockenbach Chapter 4.1) v x, y L 2 v x, y L 2 P h 1 H 1 Ψ 1, Ψ 19 0 K 1,19 = Ψ 1 x Ψ 19 x dω = 0 Ω FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 11
12 Sparsity of the stiffness matrix [Chapter 4.1] K ij = Ψ i x Ψ j x dω 0 Ω if nodes i and j are adjacent Example: K 13,j 0 if j = 7,8,12,13,14,18,19 Similarly: K i,13 0 if i = 7,8,12,13,14,18,19 max 7 nonzero elements/row and /column FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 12
13 Quadratic piecewise polynomials Ψ i x, y = Ψ 5 Ψ 18 [Chapter 4.2] 6 nodes needed FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 13
14 Quadratic piecewise polynomials [Chapter 4.2] nz = = 16% nz = = 5% FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 14
15 Cubic piecewise polynomials [Chapter 4.3] FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 15
16 Higher order piecewise polynomials Requirements for polynomial of degree d in 2D (two variables) number of nodes per edge (to guarantee continuity of the ansatz function): d + 1 (d 1 in between the vertices) 3d on the edges Number of parameters needed to define 2D polynomials x 3 x 1 x 2 xy y 2 y x 2 y xy 2 y 3 Co Lin x 4 x 3 y x 2 y 2 xy 3 y 4 d = 4 d = 6 d = 5 FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 16
17 Higher order piecewise polynomials Let s check setup for d = 4 in 2D (two variables) number of nodes per edge (to guarantee continuity of the ansatz function): d + 1 = 5 (d 1=3 in between the vertices) 3d = 12 on the edges Number of parameters needed to define 2D polynomials x 3 x 1 x 2 xy y 2 y x 2 y xy 2 y 3 Co Lin (d + 1)(d + 2) 2 = 15 Number of nodes in the middle x = 3 x 3 y x 2 y 2 xy 3 y 4 d = 4 FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 17
18 Higher order piecewise polynomials The a general piecewise polynomial takes the form u = (i) a 1 + (i) (i) + a 2 x + a3 y + (i) + a 4 x 2 (i) (i) + a 5 xy + a6 y 2 + (i) + a 7 x 3 (i) + a 8 x 2 (i) y + a 9 xy 2 (i) + a 10y 3 + (i) + a 11 x 4 (i) + a 12 x 3 (i) y + a 13 x 2 y 2 (i) + a 14 xy 3 + a 14 (i) y 4 (x, y) T i 12 Let s suppose I know the solution at the nodes x j, y j j=1 : u j = u(x j, y j ) Then the nodes should define uniquely the plane, that is, the values for a k (i) 12 j=1 d = 4 FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 18
19 Higher order piecewise polynomials The a general piecewise polynomial takes the form u 1 u 2 1 x 1 y 1 2 x 1 x 1 y 1 2 y 1 1 x 2 y 2 2 x 2 x 2 y 2 2 y 2 3 x 1 4 y 1 3 x 2 4 y 2 a 1 (i) a 2 (i) (x j, y j ) T i = u x 12 y 12 (i) a 12 W W is not singular d = 4 FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 19
20 Different discretisations of the functional space Shape functions in C 0 and in C 1 on 1D and or on 2D (quadratic or triangular) elements, the Lagrange polynomials and the Hermite polynomials, number of basis functions, conforming elements (see lecture, or more in [1], [2], [3]) Picture source: [1]: Brenner&Scott: The Mathematical Theory of FEM, Chapter 3 Construction of the finite element space [2]: Zienkiewicz&Taylor: The Finite Elment Method, Chapter 8. Standard and hierarchical element shape functions: some general families of C_0 continuity [3]: Logg&Mardal&Wells: Automated Solution of Differential Equations by the FEM- The Fenics Book, Chapter 3: Common and unusual finite elements FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 20
21 1D Example with linear nodal basis p(x) p x = ax l Strong form: l/5 l/5 l/5 l/5 l/5 u 0 = u l = 0 Weak form: Discretisation of the weak form: l EA du l dψ dx dx dx = p x ψ x dx 0 0 u x 4 i=1 u i ψ i (x) Not efficient to calculate all the elements of the stiffness matrix one by one! 4 i=1 u j EA l ψ i (x) x ψ j (x) dx = x l p(x)ψ j (x)dx Calculate element stiffness matrices and assemble K ij f j FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 21
22 1D Example with linear nodal basis l p(x) 1 l/5 2 l/5 3 l/5 4 l/5 5 l/5 6 instead: Compute stiffness matrix elementwisely and then assemble Global stiffness matrix ψ 2 K 4 e 1,1 = EA K 4 e 1,2 = EA K 4 e 2,1 = EA K 4 e 2,2 = EA ψ 3 ψ 4 (x) Ω 4 x ψ 4 (x) Ω 4 x ψ 5 (x) Ω 4 x ψ 5 (x) Ω 4 x ψ 4 (x) dx x ψ 5 (x) dx x ψ 4 (x) dx x ψ 5 (x) dx x ψ 4 ψ K e 2 1,1 K e K e 2 (1,2) 1 (2,2) K e 3 1,1 K e K e K e 2 (2,1) 2 (2,2) 3 (1,2) K e K e 4 (1,1) 4 (1,2) K e 3 (2,1) K e 3 (2,2) K 4 e (2,1) K 4 e (2,2) K 5 e (1,1) u 1 u 2 u 3 u 4 u 5 = 0 f 2 f 3 f 4 f 5 K e 4 = 4 K 4 e (1,1) K 4 e (1,2) 6 1 u K 4 e (2,1) K 4 e (2,2) K u FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 22 f
23 1D Example with linear nodal basis l p(x) instead: Compute stiffness matrix elementwisely and then assemble 1 l/5 2 l/5 3 l/5 4 l/5 5 l/ ψ 2 ψ 3 ψ 4 ψ u 1 = 0 f 4 e 1 = Ω 4 p(x)ψ 4 (x)dx 2 K e 2 1,1 K e K e 2 (1,2) 1 (2,2) u 2 f e 1 2 f e 2 1 f 4 e 2 = Ω 4 p(x)ψ 5 (x)dx 3 K e 3 1,1 K e 2 (2,1) K e 2 (2,2) K 3 e (1,2) u 3 f e 2 2 f e K e K e 4 (1,1) 4 (1,2) K e 3 (2,1) K e 3 (2,2) u 4 f 3 e 2 f 4 e 1 f e 4 = 4 f 4 e K 4 e (2,1) K 4 e (2,2) K 5 e (1,1) 1 u 5 u 6 f e 4 2 f e f 4 e 2 K u FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 23 f
24 Non-degenerate triangulation [Chapter 5.1.2] Diameter of a set: Diameter of a triangle: D T : length of longest side d T : largest circle contained in T d T D T : measures how skinny the triangle is Other definitions Τ h :triangulation (set of triangles), with h: maximal diameter of any triangle in Τ h (the length of longest side) Nondegenerate triangulation: for all the triangles in the triangulation. FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 24
25 Convergence using piecewise polynomials error = u u h E u v E v V h u u h M δ u v v V h Let s compare the best approximation u h with the proximodel with piecewise linear functions: u I n x = j=1 u Ij Ψ j x u I V h u u h E u u I E u u h M δ u u I If I can bound the expression in the r.h.s, I also bound the errors. FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 25
26 Convergence using piecewise linear functions [Chapter 5.1] Theorem: {T h } : non-degenerate family of triangulations of a polygonal domain Ω R 2 u H 2 u I : piecewise linear approximation There exists a constant C depending on Ω and the value ρ (see definition of nondegenerate triangulation) such that u u I u u I L2 Ch 2 u H2 H1 Ch u H2 where: (seminorm) FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 26
27 Convergence using piecewise linear functions [Chapter 5.1] Source: Gockenbach: Understanding and Implementing FEM FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 27
28 Convergence using piecewise higher-order polynomials Theorem: {T h } : non-degenerate family of triangulations of a polygonal domain Ω R 2 u H p+1 u I,p : piecewise d-order approximation There exists a constant C depending on Ω and the value ρ such that [Chapter 5.2] u u I u u I L2 Ch d+1 u Hd+1 H1 Ch d u Hd+1 where: See proof in [3]: Brenner&Scott: The Mathematical Theory of FEM, Chapter 4.4 FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 28
29 Convergence using piecewise higher-order polynomials d = 2 [Chapter 5.2] d = 4 Source: Gockenbach: Understanding and Implementing FEM FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 29
30 Convergence using piecewise higher-order polynomials u u I L 2 Ch d+1 u H d+1 u u I H 1 Ch d u H d+1 u u h E u u I E? convert to homogeneous problem: u u h H 1 M δ u u I H 1 M δ Chd u H d+1 = O(h d ) G:known function,g = g on Γ D u:new function that we look for u = G + u Ω κ(x) u x v x dω = Ω f x v x dω + Γ N hv x dγ Ω κ(x) G x v x dω from natural/neumann BC from essential/dirichlet BC FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 30
31 Convergence using piecewise higher-order polynomials First let s suppose homogenous Dirichlet condition: κ x u x v x dω = Ω Ω f x v x dω + Γ N hv x dγ u u h E u u I E If 0 < k 0 κ k 1 u u I E = a u u I, u u I M u u I H 1 C h d u H d+1 One can also show for inhomogeneous boundary conditions: u u h E 2C h d u H d+1 [Chapter 5.3] FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 31
32 Variational crime: curved boundary [Chapter 5.5] Source: Gockenbach: Understanding and Implementing FEM Variational crimes Céa s lemma and the error estimators may not be valid anymore but: additional errors can be also estimated (Strang) FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 32
33 Variational crime: numerical integration [Chapter 5.5] p(x) p x = ax Strong form: l u 0 = u l = 0 l/5 l/5 l/5 l/5 l/5 Weak form: Discretisation of the weak form: 4 l EA du l dψ dx dx dx = p x ψ x dx 0 0 ψ x u x u i ψ i (x) i=1 K ij f j 4 i=1 u j EA l ψ i (x) x ψ j (x) dx = x l p(x)ψ j (x)dx For more complicated trial functions, p(x) and nonconstant EA(x) difficult to calculate Example: f j = l p x ψ j x dx Numerical integration (example: Gauß-quadrature) n k=1 ω k p x k ψ j x k FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 33
34 Local/ coordinate system, isoparametric mapping 1D x Shape functions: coordinate transformation using the ansatzfunctions Transformation from local to global coordinates: x ξ = x i N 1 ξ + x i+1 N 2 ξ = N 1 ξ Stiffness matrix with isoparametric elements: K 4 e k, l = EA K 4 e k, l = EA ξ = [0,1] N 1 ξ = 1 ξ N 2 ξ = ξ ψ i (x) Ω 4 x 0 1 N k (ξ) ξ local coordinate global coordinate x i x isoparametric mapping 2 ψ j (x) dx = EA x dx dξ 1 Nl (ξ) ξ N 2 ξ ξ = 0 ξ = 1 x i x i+1 N k (ξ) Ω 4 ξ dx dξ 1 dx(ξ) dξ functions of lower order: subparametric functions of higher order: superparametric ξ N l (ξ) ξ x ξ x dx i, j [4,5] 1 k, l [1,2] FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 34 dx dξ 1 dξ dx dξ dx dξ = x dn 1 ξ i dξ dx dξ = dn 1 ξ dξ + x i+1 dn 2 ξ dξ dn 2 ξ dξ x i x i+1
35 Local/ coordinate system, isoparametric mapping 2D triangular elements Basis functions: 1 Transformation from local to global coordinates: x glob ξ, η y glob ξ, η Stiffness matrix: = N 1 ξ, η N 2 ξ, η N 1 ξ, η N 2 ξ, η N 3 ξ, η N 3 ξ, η x 1 y 1 x 2 y 2 x 3 y 3 i, j [1,2,3] FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 35
36 Local/ coordinate system, isoparametric mapping 2D triangular elements Stiffness matrix: i, j [1,2,3] Stiffness matrix with local coordinates: N j N i 1 1 η K ij = J T ξ J N T ξ 0 0 j N i η η where: J dξdη substitution rule determinant should not be negative or zero! i, j [1,2,3] 4 Ni ξ, η 4 Ni ξ, η J = i=1 ξ x i i=1 η x i i=1 4 Ni ξ, η ξ y i i=1 4 Ni ξ, η η y i FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 36
37 Local/ coordinate system, isoparametric mapping 2D triangular elements, example Transformation from local to global coordinates (isoparametric mapping): x glob ξ, η y glob ξ, η x ξ, η y ξ, η = = N 1 ξ, η N 2 ξ, η N 1 ξ, η 1 ξ η ξ 1 ξ η ξ N 2 ξ, η η η N 3 ξ, η N 3 ξ, η FEM and its Convergence Dr. Noemi Friedman PDE 2 Seite 37 x 1 y 1 x 2 y 2 x 3 y 3
38 Local/ coordinate system, isoparametric mapping 2D triangular elements, example Stiffness matrix with local coordinates: K ij = η J T 0 N j ξ N j η J T N i ξ N i η J dξdη i, j [1,2,3] FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 38
39 Condition number of the stiffness matrix Condition number What happens with the roundoff errors in K = LU = K + δk K Ku = f K + δk u = f + δf u u u λ max λ min δf f λ max λ min = κ(k) Condition number of K with nodal bases with 2D triangular mesh: For the Poisson equation N i=1 c i Ψ i (x) Ψ j x dω = Ω f x Ψ j x dω Ω turns il-conditioned for refined mesh!!! K ij FEM and its convergence Dr. Noemi Friedman PDE 2 Seite 39
Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationIntroduction to PDEs and Numerical Methods Tutorial 10. Finite Element Analysis
Patzhater für Bid, Bid auf Titefoie hinter das Logo einsetzen Introduction to PDEs and Numerica Methods Tutoria. Finite Eement Anaysis Dr. Noemi Friedman, 3..2 FROM STRONG FORM TO WEAK FORM inhomogeneous
More informationIntroduction to PDEs and Numerical Methods Tutorial 1: Overview of essential linear algebra and analysis
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Tutorial 1: Overview of essential linear algebra and analysis Dr. Noemi Friedman, 25.10.201.
More informationNumerical methods for PDEs FEM - abstract formulation, the Galerkin method
Patzhater für Bid, Bid auf Titefoie hinter das Logo einsetzen Numerica methods for PDEs FEM - abstract formuation, the Gaerkin method Dr. Noemi Friedman Contents of the course Fundamentas of functiona
More informationNumerical methods for PDEs FEM implementation: element stiffness matrix, isoparametric mapping, assembling global stiffness matrix
Platzhaltr für Bild, Bild auf Titlfoli hintr das Logo instzn Numrical mthods for PDEs FEM implmntation: lmnt stiffnss matrix, isoparamtric mapping, assmbling global stiffnss matrix Dr. Nomi Fridman Contnts
More informationIntroduction to PDEs and Numerical Methods Tutorial 5. Finite difference methods equilibrium equation and iterative solvers
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Tutorial 5. Finite difference methods equilibrium equation and iterative solvers Dr. Noemi
More informationIntroduction to PDEs and Numerical Methods Lecture 7. Solving linear systems
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 7. Solving linear systems Dr. Noemi Friedman, 09.2.205. Reminder: Instationary heat
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationResources. Introduction to the finite element method. History. Topics
Resources Introduction to the finite element method M. M. Sussman sussmanm@math.pitt.edu Office Hours: 11:1AM-12:1PM, Thack 622 May 12 June 19, 214 Strang, G., Fix, G., An Analysis of the Finite Element
More informationPlatzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen PARAMetric UNCertainties, Budapest STOCHASTIC PROCESSES AND FIELDS Noémi Friedman Institut für Wissenschaftliches Rechnen, wire@tu-bs.de
More informationPlatzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 6: Numerical solution of the heat equation with FD method: method of lines, Euler
More informationAlgorithms for Scientific Computing
Algorithms for Scientific Computing Finite Element Methods Michael Bader Technical University of Munich Summer 2016 Part I Looking Back: Discrete Models for Heat Transfer and the Poisson Equation Modelling
More information1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:
Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationIntroduction to PDEs and Numerical Methods Tutorial 11. 2D elliptic equations
Platzalter für Bild, Bild auf Titelfolie inter das Logo einsetzen Introduction to PDEs and Numerical Metods Tutorial 11. 2D elliptic equations Dr. Noemi Friedman, 3. 1. 215. Overview Introduction (classification
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationIntroduction to PDEs and Numerical Methods Lecture 1: Introduction
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 1: Introduction Dr. Noemi Friedman, 28.10.2015. Basic information on the course Course
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationIntroduction to PDEs and Numerical Methods Tutorial 4. Finite difference methods stability, concsistency, convergence
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Tutorial 4. Finite difference methods stability, concsistency, convergence Dr. Noemi Friedman,
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationIsogeometric Analysis:
Isogeometric Analysis: some approximation estimates for NURBS L. Beirao da Veiga, A. Buffa, Judith Rivas, G. Sangalli Euskadi-Kyushu 2011 Workshop on Applied Mathematics BCAM, March t0th, 2011 Outline
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationMath 660-Lecture 15: Finite element spaces (I)
Math 660-Lecture 15: Finite element spaces (I) (Chapter 3, 4.2, 4.3) Before we introduce the concrete spaces, let s first of all introduce the following important lemma. Theorem 1. Let V h consists of
More informationR T (u H )v + (2.1) J S (u H )v v V, T (2.2) (2.3) H S J S (u H ) 2 L 2 (S). S T
2 R.H. NOCHETTO 2. Lecture 2. Adaptivity I: Design and Convergence of AFEM tarting with a conforming mesh T H, the adaptive procedure AFEM consists of loops of the form OLVE ETIMATE MARK REFINE to produce
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationWEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
More informationConvergence Order Studies for Elliptic Test Problems with COMSOL Multiphysics
Convergence Order Studies for Elliptic Test Problems with COMSOL Multiphysics Shiming Yang and Matthias K. Gobbert Abstract. The convergence order of finite elements is related to the polynomial order
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationA brief introduction to finite element methods
CHAPTER A brief introduction to finite element methods 1. Two-point boundary value problem and the variational formulation 1.1. The model problem. Consider the two-point boundary value problem: Given a
More informationFinite Element Method for Ordinary Differential Equations
52 Chapter 4 Finite Element Method for Ordinary Differential Equations In this chapter we consider some simple examples of the finite element method for the approximate solution of ordinary differential
More informationChapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method
Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.
More informationBuilding a finite element program in MATLAB Linear elements in 1d and 2d
Building a finite element program in MATLAB Linear elements in 1d and 2d D. Peschka TU Berlin Supplemental material for the course Numerische Mathematik 2 für Ingenieure at the Technical University Berlin,
More informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationBoundary Value Problems and Iterative Methods for Linear Systems
Boundary Value Problems and Iterative Methods for Linear Systems 1. Equilibrium Problems 1.1. Abstract setting We want to find a displacement u V. Here V is a complete vector space with a norm v V. In
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More information2 Elliptic Differential Equations
2 Elliptic Differential Equations 2.1 Classical solutions As far as existence and uniqueness results for classical solutions are concerned, we restrict ourselves to linear elliptic second order elliptic
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationHierarchical Parallel Solution of Stochastic Systems
Hierarchical Parallel Solution of Stochastic Systems Second M.I.T. Conference on Computational Fluid and Solid Mechanics Contents: Simple Model of Stochastic Flow Stochastic Galerkin Scheme Resulting Equations
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationHIGH DEGREE IMMERSED FINITE ELEMENT SPACES BY A LEAST SQUARES METHOD
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 14, Number 4-5, Pages 604 626 c 2017 Institute for Scientific Computing and Information HIGH DEGREE IMMERSED FINITE ELEMENT SPACES BY A LEAST
More informationLecture 3: Function Spaces I Finite Elements Modeling. Bruno Lévy
Lecture 3: Function Spaces I Finite Elements Modeling Bruno Lévy Overview 1. Motivations 2. Function Spaces 3. Discretizing a PDE 4. Example: Discretizing the Laplacian 5. Eigenfunctions Spectral Mesh
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationFinite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems
Finite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems Introduction Till now we dealt only with finite elements having straight edges.
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationMesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain
Mesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain Stephen Edward Moore Johann Radon Institute for Computational and Applied Mathematics Austrian Academy of Sciences,
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017
ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip
More informationChapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction
Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction The a posteriori error estimation of finite element approximations of elliptic boundary value problems has reached
More informationNew class of finite element methods: weak Galerkin methods
New class of finite element methods: weak Galerkin methods Xiu Ye University of Arkansas at Little Rock Second order elliptic equation Consider second order elliptic problem: a u = f, in Ω (1) u = 0, on
More informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationIntroduction to PDEs: Assignment 9: Finite elements in 1D
Institute of Scientific Computing Technical University Braunschweig Dr. Elmar Zander Dr. Noemi Friedman Winter Semester 213/14 Assignment 9 Due date: 6.2.214 Introduction to PDEs: Assignment 9: Finite
More informationFinite Elements. Colin Cotter. January 18, Colin Cotter FEM
Finite Elements January 18, 2019 The finite element Given a triangulation T of a domain Ω, finite element spaces are defined according to 1. the form the functions take (usually polynomial) when restricted
More informationWeighted Residual Methods
Weighted Residual Methods Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents Problem definition. oundary-value Problem..................
More informationNumerical Methods for Partial Differential Equations
Numerical Methods for Partial Differential Equations Eric de Sturler University of Illinois at Urbana-Champaign Read section 8. to see where equations of type (au x ) x = f show up and their (exact) solution
More informationOptimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms
Optimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms Marcus Sarkis Worcester Polytechnic Inst., Mass. and IMPA, Rio de Janeiro and Daniel
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
More informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationFEM Convergence for PDEs with Point Sources in 2-D and 3-D
FEM Convergence for PDEs with Point Sources in 2-D and 3-D Kourosh M. Kalayeh 1, Jonathan S. Graf 2 Matthias K. Gobbert 2 1 Department of Mechanical Engineering 2 Department of Mathematics and Statistics
More informationOrthogonality of hat functions in Sobolev spaces
1 Orthogonality of hat functions in Sobolev spaces Ulrich Reif Technische Universität Darmstadt A Strobl, September 18, 27 2 3 Outline: Recap: quasi interpolation Recap: orthogonality of uniform B-splines
More informationTHE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER QUADRILATERAL MESHES. Zhong-Ci Shi and Xuejun Xu.
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 9, Number, January 2008 pp. 63 82 THE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER
More informationPDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM.
PDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM. Sergey Korotov Basque Center for Applied Mathematics / IKERBASQUE http://www.bcamath.org & http://www.ikerbasque.net 1 Introduction The analytical solution
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationMatrix construction: Singular integral contributions
Matrix construction: Singular integral contributions Seminar Boundary Element Methods for Wave Scattering Sophie Haug ETH Zurich November 2010 Outline 1 General concepts in singular integral computation
More informationWave Theory II (6) Finite Element Method
Wave Theory II (6) Finite Method Jun-ichi Takada (takada@ide.titech.ac.jp) In this lecture, the finite element method (FEM) is described. The Helmholtz equation and the boundary condition are transformed
More informationCIVL4332 L1 Introduction to Finite Element Method
CIVL L Introduction to Finite Element Method CIVL L Introduction to Finite Element Method by Joe Gattas, Faris Albermani Introduction The FEM is a numerical technique for solving physical problems such
More informationNumerical Simulation of Potential Flow using the Finite Element Method
Numerical Simulation of Potential Flow using the Finite Element Method Cynthia Vanessa Flores (CSUN), Eduardo Rogério Fávaro(Ibilce-Unesp) Gilberto Pereira Sassi (ICMC-USP) Leonardo Koller Sacht(UFSC)
More information