Isogeometric Analysis:
|
|
- Everett Benson
- 6 years ago
- Views:
Transcription
1 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
2 Outline of the talk 1 Introduction to Isogeometric Analysis 2 B-splines and NURBS 3 Approximation theory for NURBS 4 Conclusions
3 What is Isogeometric Analysis? Novel technique for the discretization of PDEs T. J. R. Hughes, J. A. Cottrell, and Y. Bazilevs, Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement, CMAME (2005).
4 What is Isogeometric Analysis? Novel technique for the discretization of PDEs Aim: T. J. R. Hughes, J. A. Cottrell, and Y. Bazilevs, Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement, CMAME (2005). Improve the connection between numerical simulation of physical phenomena and Computer Aided Design (CAD). Eliminate/reduce the approximation of the computational domain and the need of remeshing
5 What is Isogeometric Analysis? Novel technique for the discretization of PDEs Aim: Tools: T. J. R. Hughes, J. A. Cottrell, and Y. Bazilevs, Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement, CMAME (2005). Improve the connection between numerical simulation of physical phenomena and Computer Aided Design (CAD). Eliminate/reduce the approximation of the computational domain and the need of remeshing Use the geometry provided by CAD Represent unknown fields in the same basis functions used for the geometry
6 What is Isogeometric Analysis? Variational formulation of a BVP: Find u X such that B(u, v) = F(v), v Y.
7 What is Isogeometric Analysis? Variational formulation of a BVP: Find u X such that B(u, v) = F(v), v Y. X, Y reflexive Banach spaces B: X Y R continuous bilinear form F: Y R continuous linear functional
8 What is Isogeometric Analysis? Variational formulation of a BVP: Find u X such that B(u, v) = F(v), v Y. X, Y reflexive Banach spaces B: X Y R continuous bilinear form F: Y R continuous linear functional Galerkin approximation: Find u n X n such that B(u n, v n ) = F(v n ), v n Y n. X n X, Y n Y subspaces of dimension n
9 What is Isogeometric Analysis? Variational formulation of a BVP: Find u X such that B(u, v) = F(v), v Y. X, Y reflexive Banach spaces B: X Y R continuous bilinear form F: Y R continuous linear functional Galerkin approximation: Find u n X n such that B(u n, v n ) = F(v n ), v n Y n. X n X, Y n Y subspaces of dimension n FEM: piecewise polynomials Spectral methods: Orthogonal (global) polynomials IGA: B-splines, Non-Uniform Rational B-Splines (NURBS)
10 Why to use Isogeometric Analysis? In engineering problems geometry is usually defined by Computer Aided Design (CAD). CAD and FEM use different geometry descriptions
11 Why to use Isogeometric Analysis? In engineering problems geometry is usually defined by Computer Aided Design (CAD). CAD and FEM use different geometry descriptions CAD and IGA use the same geometry description
12 B-splines in one dimension n N number of basis functions p N 0 degree of polynomials Ξ := {0 = ξ 1 ξ 2 ξ n+p+1 = 1} knot vector {ζ 1,..., ζ m } mesh in [0, 1] r j number of repetitions of ζ j in Ξ k = {k 1,..., k m }, k j = p r j + 1
13 B-splines in one dimension n N number of basis functions p N 0 degree of polynomials Ξ := {0 = ξ 1 ξ 2 ξ n+p+1 = 1} knot vector {ζ 1,..., ζ m } mesh in [0, 1] r j number of repetitions of ζ j in Ξ k = {k 1,..., k m }, k j = p r j + 1 B i,p piecewise polynomial of degree p and continuous derivatives up to the order k j 1 at knot ζ j, with compact support in [ξ i, ξ i+p+1 ] S p k (Ξ) = span{b 1,p,..., B n,p }
14 B-splines in one dimension n N number of basis functions p N 0 degree of polynomials Ξ := {0 = ξ 1 ξ 2 ξ n+p+1 = 1} knot vector {ζ 1,..., ζ m } mesh in [0, 1] r j number of repetitions of ζ j in Ξ k = {k 1,..., k m }, k j = p r j + 1 B i,p piecewise polynomial of degree p and continuous derivatives up to the order k j 1 at knot ζ j, with compact support in [ξ i, ξ i+p+1 ] S p k (Ξ) = span{b 1,p,..., B n,p }
15 B-splines in higher dimensions B-splines in d dimensions B i1...i d (x 1,..., x d ) = B i1,p 1 (x 1 )... B id,p d (x d ). S p 1,...,p d k 1,...,k d = S p 1 k 1 S p d k d = span{b i1...i d } n 1,...,n d i=1,...,i d =1.
16 B-splines in higher dimensions B-splines in d dimensions B i1...i d (x 1,..., x d ) = B i1,p 1 (x 1 )... B id,p d (x d ). S p 1,...,p d k 1,...,k d = S p 1 k 1 S p d k d = span{b i1...i d } n 1,...,n d i=1,...,i d =1. B-spline curves in R d F(x) = n i=1 B i,p(x)c i, C i R d control points.
17 B-splines in higher dimensions B-splines in d dimensions B i1...i d (x 1,..., x d ) = B i1,p 1 (x 1 )... B id,p d (x d ). S p 1,...,p d k 1,...,k d = S p 1 k 1 S p d k d = span{b i1...i d } n 1,...,n d i=1,...,i d =1. B-spline curves in R d F(x) = n i=1 B i,p(x)c i, C i R d control points. B-spline surfaces in R d F(x, y) = n 1,n 2 i,j=1 B ij(x, y)c ij, C ij R d control points.
18 NURBS in Ω R d NURBS in R d are conic projections of B-splines in R d+1.
19 NURBS in Ω R d NURBS in R d are conic projections of B-splines in R d+1. Ω = (0, 1) d : w := n 1,...,n d i 1 =1,...,i d =1 w i 1...i d B i1...i d, w i1...i d 1, weighting function R i1...i d = w i 1...i d B i1...i d w NURBS basis functions in Ω
20 NURBS in Ω R d NURBS in R d are conic projections of B-splines in R d+1. Ω = (0, 1) d : w := n 1,...,n d i 1 =1,...,i d =1 w i 1...i d B i1...i d, w i1...i d 1, weighting function R i1...i d = w i 1...i d B i1...i d w NURBS basis functions in Ω Ω R d : F : Ω Ω, F = n1,...n d i 1 =1,i d =1 C i 1...i d R i1...i d geometrical map K = F (Q), Q cartesian mesh in Ω N i1...i d = R i1...i d F 1 NURBS basis functions in Ω N p 1,...,p d k 1,...,k d (K) = span{n i1...i d }
21 Properties of B-splines and NURBS They form a partition of unity (for open knot vectors). They are C k continuous, with 0 k p 1. The support of each basis function is compact. NURBS represent exactly a wide class of curves, e.g. conic sections. Three kinds of refinement can be performed: h-refinement = mesh refinement p-refinement = degree elevation k-refinement = regularity adjustment h-, p- and k-refinement can be performed without changing the geometry
22 Geometry description and refinement Coarsest mesh: geometry description Parametric domain, Ω = (0, 1) 2 Physical domain Ω F { R ij = w } { ijb ij N ij = w ( ) } wij B ij F 1 w
23 Geometry description and refinement Coarsest mesh: geometry description Parametric domain, Ω = (0, 1) 2 Physical domain Ω F { R ij = w } { ijb ij N ij = w ( ) } wij B ij F 1 w The geometrical map F and the weight w are fixed at the coarsest level of discretization!
24 Geometry description and refinement First refinement Parametric domain, Ω = (0, 1) 2 Physical domain Ω F { R ij = w } { ijb ij N ij = w ( ) } wij B ij F 1 w The geometrical map F and the weight w are fixed at the coarsest level of discretization!
25 Geometry description and refinement Second refinement... and so on Parametric domain, Ω = (0, 1) 2 Physical domain Ω F { R ij = w } { ijb ij N ij = w ( ) } wij B ij F 1 w The geometrical map F and the weight w are fixed at the coarsest level of discretization!
26 A priori error estimates Lemma (Cea) Suppose X n X, is a family of finite-dimensional subspaces of a HIlbert space X. Suppose B: X X R is a bounded, coercive bilinear form and F: X R is a continuous functional. Then the problem of finding u n X n such that B(u n, v n ) = F(v n ), v n X n, has a unique solution. If u X is the solution of B(u, v) = F(v), v X, then there exists a constant C independent of u, u n and n such that u u n X C inf w n X n u w n X.
27 Previous results: hp-estimates for FEM Theorem Let Ω R 2 be a polygon, T a parallelogram mesh in Ω with at most one hanging node per edge and let h denote its diameter. Then, for any 2 s p + 1 and any u H s (Ω), there exists Πu S p 0 (T) such that u Πu H 1 (Ω) Ch s 1 p (s 1) u H s (Ω), where C is a constant independent of h and p. C. Schwab, p- and hp- Finite Element Methods. Theory and applications in Solid and Fluid Mechanics, Oxford University Press, Oxford (1998).
28 Previous results: h-estimates for NURBS in 2D Theorem Let Ω = F ( Ω) R 2, K = F (Q) a mesh in Ω and let h denote its diameter. Then, for any 0 l s p + 1 and for any u H s (Ω), there exists Πu N p k (K) such that u Πu H l (Ω) C(w, F, p, k)h s l u H s (Ω), where C is a constant independent of h, but possibly depending on p and k. Y. Bazilevs, L. Beirão da Veiga, J. A. Cottrell, T. J. R. Hughes, and G. Sangalli, Isogeometric analysis: approximation, stability and error estimates for h-refined meshes, Math. Models Methods Appl. Sci., 16 (2006), pp
29 Our result: hpk-estimates for NURBS in 2D Theorem Let Ω = F ( Ω) R 2, K = F (Q) a mesh in Ω with diameter h. Then, if 2k p + 1, for any u H s (Ω), with 2k s p + 1, and for any 0 l s p + 1 there exists Πu N p k (K) such that u Πu H l (Ω) C(w, F )h s l (p k + 1) (s l) u H s (Ω), where C is a constant independent of h, p and k. L. Beirão da Veiga, A. Buffa, J. Rivas, G. Sangalli, Some estimates for h p k refinement in Isogeometric Analysis, to appear in Numer. Math.
30 Construction of the new projection operator (!1,1) I i T i tensor product F
31 Legendre polynomials Definition (Legendre polynomial of degree i) L i (x) = 1 i! 2 i d i dx i ( (x 2 1) i), i = 0, 1,...
32 Legendre polynomials Definition (Legendre polynomial of degree i) L i (x) = 1 i! 2 i d i dx i ( (x 2 1) i), i = 0, 1,... Definition (L 2 ( 1, 1)-orthogonal projection of order N N) π N ϕ(x) = N i=0 ˆϕ i L i (x), where ˆϕ i = 2i ϕ(x)l i (x) dx.
33 Legendre polynomials Definition (Legendre polynomial of degree i) L i (x) = 1 i! 2 i d i dx i ( (x 2 1) i), i = 0, 1,... Definition (L 2 ( 1, 1)-orthogonal projection of order N N) π N ϕ(x) = N i=0 ˆϕ i L i (x), where ˆϕ i = 2i ϕ(x)l i (x) dx. Definition (Primitives of Legendre polynomials) For n 0, Ψ i,n is the n-th primitive of L i, Ψ i,0 (x) = L i (x), Ψ i,n (x) = x 1 Ψ i,n 1 (ξ) dξ, n = 1, 2,...
34 Projection operator in ( 1, 1) Definition p, k nonnegative integers; S p (Λ) set of polynomials of degree p in Λ = ( 1, 1); π p,k : H k (Λ) S p (Λ) is defined as: ( π p,k u) (k) (x) =π p k u (k) (x), x Λ, ( π p,k u) (l) ( 1) =u (l) ( 1), l = 0, 1,..., k 1,
35 Projection operator in ( 1, 1) Definition p, k nonnegative integers; S p (Λ) set of polynomials of degree p in Λ = ( 1, 1); π p,k : H k (Λ) S p (Λ) is defined as: ( π p,k u) (k) (x) =π p k u (k) (x), x Λ, ( π p,k u) (l) ( 1) =u (l) ( 1), l = 0, 1,..., k 1, If u (k) (x) = i=0 α il i (x), then p k π p,k u(x) = α i Ψ i,k (x) + i=0 k 1 m=0 u (m) ( 1) (x + 1)m m!
36 Projection operator in ( 1, 1) Definition p, k nonnegative integers; S p (Λ) set of polynomials of degree p in Λ = ( 1, 1); π p,k : H k (Λ) S p (Λ) is defined as: ( π p,k u) (k) (x) =π p k u (k) (x), x Λ, ( π p,k u) (l) ( 1) =u (l) ( 1), l = 0, 1,..., k 1, If u (k) (x) = i=0 α il i (x), then p k π p,k u(x) = α i Ψ i,k (x) + i=0 k 1 m=0 u (m) (x + 1)m ( 1) m! If p 2k 1, ( π p,k u) (l) (1) = u (l) (1), l = 0, 1,..., k 1.
37 Spline approximation on the reference domain [0, 1] Definition p, k nonnegative integers, with 2k 1 p; {0 = ζ 1 < ζ 2 < < ζ m = 1}, I i = (ζ i, ζ i+1 ), 1 i m 1; T i : Λ I i linear mapping. The (local) projection operator π i p,k : Hk (I i ) S p is defined as: π i p,k u T i = ( π p,k (u T i ) ).
38 Spline approximation on the reference domain [0, 1] Definition p, k nonnegative integers, with 2k 1 p; {0 = ζ 1 < ζ 2 < < ζ m = 1}, I i = (ζ i, ζ i+1 ), 1 i m 1; T i : Λ I i linear mapping. The (local) projection operator π i p,k : Hk (I i ) S p is defined as: π i p,k u T i = ( π p,k (u T i ) ). The (global) projection operator π p,k : H k (0, 1) S p k is defined as: (π p,k u) Ii = π i p,k, i = 0,..., m 1.
39 Spline approximation on the reference domain [0, 1] Definition p, k nonnegative integers, with 2k 1 p; {0 = ζ 1 < ζ 2 < < ζ m = 1}, I i = (ζ i, ζ i+1 ), 1 i m 1; T i : Λ I i linear mapping. The (local) projection operator π i p,k : Hk (I i ) S p is defined as: π i p,k u T i = ( π p,k (u T i ) ). The (global) projection operator π p,k : H k (0, 1) S p k is defined as: (π p,k u) Ii = π i p,k, i = 0,..., m 1. Remark (π p,k u) (l) (ζ i ) = u (l) (ζ i ), 1 i m, 0 l k 1.
40 Error estimate for π p,k Theorem p, k nonnegative integers, 2k 1 p; {0 = ζ 1 < < ζ m = 1}, I i = (ζ i, ζ i+1 ), h i = ζ i+1 ζ i ; u (k) H s (I i ) for some 0 s κ = p k + 1. Then, for l = 0,..., k, u (l) (π p,k u) (l) 2 L 2 (I i ) ( hi ) 2(s+k l) (κ s)! (κ (k l))! 2 (κ + s)! (κ + (k l))! u(k) 2 H s (I i ).
41 Error estimate for π p,k Theorem p, k nonnegative integers, 2k 1 p; {0 = ζ 1 < < ζ m = 1}, I i = (ζ i, ζ i+1 ), h i = ζ i+1 ζ i ; u (k) H s (I i ) for some 0 s κ = p k + 1. Then, for l = 0,..., k, u (l) (π p,k u) (l) 2 L 2 (I i ) ( hi ) 2(s+k l) (κ s)! (κ (k l))! 2 (κ + s)! (κ + (k l))! u(k) 2 H s (I i ). Consequently, for u H σ (I i ), k σ p + 1, and l = 0,..., k, there exists a constant C independent of u, l, σ, h i, p and k, s.t. u π p,k u H l (I i ) Ch i σ l (p k + 1) (σ l) u H σ (I i ).
42 Error estimate for π p,k Theorem p, k nonnegative integers, 2k 1 p; {0 = ζ 1 < < ζ m = 1}, I i = (ζ i, ζ i+1 ), h i = ζ i+1 ζ i ; u (k) H s (I i ) for some 0 s κ = p k + 1. Then, for l = 0,..., k, u (l) (π p,k u) (l) 2 L 2 (I i ) ( hi ) 2(s+k l) (κ s)! (κ (k l))! 2 (κ + s)! (κ + (k l))! u(k) 2 H s (I i ). Consequently, for u H σ (0, 1), k σ p + 1, and l = 0,..., k, there exists a constant C independent of u, l, σ, h, p and k, s.t. u π p,k u H l (0,1) Ch σ l (p k + 1) (σ l) u H σ (0,1).
43 Spline approximation on the reference domain Ω = [0, 1] 2 Definition p = (p 1, p 2 ), k = (k 1, k 2 ), 2k d 1 p d, d = 1, 2; Q ij = I i J j = (ζ i,1, ζ i+1,1 ) (ζ j,2, ζ j+1,2 ) Q h, H k 1,k 2 (Q ij ) = H k 1(I i, H k 2(J j )) S p 1,p 2 (Q ij ) = {u: Q ij R : u(, y) S p 1(I i ), u(x, ) S p 2(J j )}. The (local) projection operator Π ij p,k : Hk 1,k 2 (Q ij ) S p 1,p 2 (Q ij ) is: Π ij p,k = πi p 1,k 1 π j p 2,k 2.
44 Spline approximation on the reference domain Ω = [0, 1] 2 Definition p = (p 1, p 2 ), k = (k 1, k 2 ), 2k d 1 p d, d = 1, 2; Q ij = I i J j = (ζ i,1, ζ i+1,1 ) (ζ j,2, ζ j+1,2 ) Q h, H k 1,k 2 (Q ij ) = H k 1(I i, H k 2(J j )) S p 1,p 2 (Q ij ) = {u: Q ij R : u(, y) S p 1(I i ), u(x, ) S p 2(J j )}. The (local) projection operator Π ij p,k : Hk 1,k 2 (Q ij ) S p 1,p 2 (Q ij ) is: Π ij p,k = πi p 1,k 1 π j p 2,k 2. The (global) projection operator Π p,k : H k 1,k 2 ( Ω) S p 1,p 2 k 1,k 2 (Q h ) is: (Π p,k v) Qij = (Π ij p,k v), Q ij Q h.
45 Error estimate for Π p,k Theorem p 1 = p 2 = p; k 1, k 2 be nonnegative integers, 2k d 1 p for d = 1, 2 k = min{k 1, k 2 }, k = max{k 1, k 2 }; Q ij = (ζ i,1, ζ i+1,1 ) (ζ j,2, ζ j+1,2 ), h ij = max{ζ i+1,1 ζ i,1, ζ j+1,2 ζ j,2 },, h = max h ij ; v H σ (Q ij ) with k 1 + k 2 σ p + 1. Then, for all integers 0 l k, there exists a positive constant C, independent of v, σ, l, h, p, k 1 and k 2, such that, v Π p,k v H l ( Ω) Chσ l (p k + 1) (σ l) v H σ ( Ω).
46 NURBS approximation in the physical domain Ω R 2 Definition p = (p 1, p 2 ), k = (k 1, k 2 ), 2k d 1 p d, d = 1, 2; K ij = F (Q ij ), Q ij Q h, The (local) projection operator for functions defined on K ij is: Π ij N : Hσ (K ij ) N p 1,p 2 k 1,k 2, σ k 1 + k 2 Π ij Πij p,k (w (v F )) N (v) F =. w
47 NURBS approximation in the physical domain Ω R 2 Definition p = (p 1, p 2 ), k = (k 1, k 2 ), 2k d 1 p d, d = 1, 2; K ij = F (Q ij ), Q ij Q h, The (local) projection operator for functions defined on K ij is: Π ij N : Hσ (K ij ) N p 1,p 2 k 1,k 2, σ k 1 + k 2 Π ij Πij p,k (w (v F )) N (v) F =. w The (global) projection operator for functions defined on Ω is Π N : H σ (Ω) N p 1,p 2 k 1,k 2, σ k 1 + k 2, Π N (v) Kij = Π ij N (v K ij ) K ij K h.
48 Error estimate for Π N Theorem p 1 = p 2 = p; k 1, k 2 nonnegative integers, 2k d 1 p for d = 1, 2 k = min{k 1, k 2 }, k = max{k 1, k 2 }; w and F fixed at the coarsest level of discretization; K K h, h K = diam K, h = max h K ; v H σ (K) with k 1 + k 2 σ p + 1; Then for l = 0,..., k, there exists a constant C = C(w, F ), independent of v, σ, l, h, p, k 1 and k 2 such that v Π N (v) H l (Ω) Ch σ l (p k + 1) (σ l) v H σ (Ω).
49 Concluding remarks We have constructed a new projection operator onto NURBS spaces in 2 dimensions and given error estimates in Sobolev norms which are explicit in the three discretization parameters: degree p, regularity k and mesh size h.
50 Concluding remarks We have constructed a new projection operator onto NURBS spaces in 2 dimensions and given error estimates in Sobolev norms which are explicit in the three discretization parameters: degree p, regularity k and mesh size h. A restriction on the regularity must be imposed, namely 2k 1 p.
51 Concluding remarks We have constructed a new projection operator onto NURBS spaces in 2 dimensions and given error estimates in Sobolev norms which are explicit in the three discretization parameters: degree p, regularity k and mesh size h. A restriction on the regularity must be imposed, namely 2k 1 p. The case of higher regularity, up to k = p remains open.
52 Concluding remarks We have constructed a new projection operator onto NURBS spaces in 2 dimensions and given error estimates in Sobolev norms which are explicit in the three discretization parameters: degree p, regularity k and mesh size h. A restriction on the regularity must be imposed, namely 2k 1 p. The case of higher regularity, up to k = p remains open. Thank you
An introduction to Isogeometric Analysis
An introduction to Isogeometric Analysis A. Buffa 1, G. Sangalli 1,2, R. Vázquez 1 1 Istituto di Matematica Applicata e Tecnologie Informatiche - Pavia Consiglio Nazionale delle Ricerche 2 Dipartimento
More informationCondition number estimates for matrices arising in the isogeometric discretizations
www.oeaw.ac.at Condition number estimates for matrices arising in the isogeometric discretizations K. Gahalaut, S. Tomar RICAM-Report -3 www.ricam.oeaw.ac.at Condition number estimates for matrices arising
More informationIsogeometric mortaring
Isogeometric mortaring E. Brivadis, A. Buffa, B. Wohlmuth, L. Wunderlich IMATI E. Magenes - Pavia Technical University of Munich A. Buffa (IMATI-CNR Italy) IGA mortaring 1 / 29 1 Introduction Splines Approximation
More informationBenchmarking high order finite element approximations for one-dimensional boundary layer problems
Benchmarking high order finite element approximations for one-dimensional boundary layer problems Marcello Malagù,, Elena Benvenuti, Angelo Simone Department of Engineering, University of Ferrara, Italy
More informationMultipatch Space-Time Isogeometric Analysis of Parabolic Diffusion Problems. Ulrich Langer, Martin Neumüller, Ioannis Toulopoulos. G+S Report No.
Multipatch Space-Time Isogeometric Analysis of Parabolic Diffusion Problems Ulrich Langer, Martin Neumüller, Ioannis Toulopoulos G+S Report No. 56 May 017 Multipatch Space-Time Isogeometric Analysis of
More informationIsogeometric modeling of Lorentz detuning in linear particle accelerator cavities
Isogeometric modeling of Lorentz detuning in linear particle accelerator cavities Mauro Bonafini 1,2, Marcella Bonazzoli 3, Elena Gaburro 1,2, Chiara Venturini 1,3, instructor: Carlo de Falco 4 1 Department
More informationPhD dissertation defense
Isogeometric mortar methods with applications in contact mechanics PhD dissertation defense Brivadis Ericka Supervisor: Annalisa Buffa Doctor of Philosophy in Computational Mechanics and Advanced Materials,
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 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 informationDual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations
Dual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations Christoph Hofer and Ulrich Langer Doctoral Program Computational Mathematics Numerical
More informationDomain Decomposition Preconditioners for Isogeometric Discretizations
Domain Decomposition Preconditioners for Isogeometric Discretizations Luca F Pavarino, Università di Milano, Italy Lorenzo Beirao da Veiga, Università di Milano, Italy Durkbin Cho, Dongguk University,
More informationThe space of C 1 -smooth isogeometric spline functions on trilinearly parameterized volumetric two-patch domains. Katharina Birner and Mario Kapl
The space of C 1 -smooth isogeometric spline functions on trilinearly parameterized volumetric two-patch domains Katharina Birner and Mario Kapl G+S Report No. 69 March 2018 The space of C 1 -smooth isogeometric
More informationNumerical 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 informationEFFICIENT MULTIGRID BASED SOLVERS FOR ISOGEOMETRIC ANALYSIS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK EFFICIENT MULTIGRID BASED SOLVERS FOR ISOGEOMETRIC
More informationSpace-time isogeometric analysis of parabolic evolution equations
www.oeaw.ac.at Space-time isogeometric analysis of parabolic evolution equations U. Langer, S.E. Moore, M. Neumüller RICAM-Report 2015-19 www.ricam.oeaw.ac.at SPACE-TIME ISOGEOMETRIC ANALYSIS OF PARABOLIC
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 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 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 informationNumerical 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 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 informationSerendipity Finite Element Methods for Cardiac Electrophysiology Models
Serendipity Finite Element Methods for Cardiac Electrophysiology Models Andrew Gillette Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/ Andrew Gillette -
More informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More informationISOGEOMETRIC METHODS IN STRUCTURAL DYNAMICS AND WAVE PROPAGATION
COMPDYN 29 ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, N.D. Lagaros, M. Fragiadakis (eds.) Rhodes, Greece, 22-24 June 29 ISOGEOMETRIC
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 informationOn a Discontinuous Galerkin Method for Surface PDEs
On a Discontinuous Galerkin Method for Surface PDEs Pravin Madhavan (joint work with Andreas Dedner and Bjo rn Stinner) Mathematics and Statistics Centre for Doctoral Training University of Warwick Applied
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 informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
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 informationThe Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods
The Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods by Hae-Soo Oh Department of Mathematics, University of North Carolina at Charlotte, Charlotte, NC 28223 June
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 informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationIsogeometric Analysis with Geometrically Continuous Functions on Two-Patch Geometries
Isogeometric Analysis with Geometrically Continuous Functions on Two-Patch Geometries Mario Kapl a Vito Vitrih b Bert Jüttler a Katharina Birner a a Institute of Applied Geometry Johannes Kepler University
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 informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationGAKUTO International Series
GAKUTO International Series Mathematical Sciences and Applications, Vol.28(2008) Proceedings of Fourth JSIAM-SIMMAI Seminar on Industrial and Applied Mathematics, pp.139-148 A COMPUTATIONAL APPROACH TO
More informationMultiscale methods for time-harmonic acoustic and elastic wave propagation
Multiscale methods for time-harmonic acoustic and elastic wave propagation Dietmar Gallistl (joint work with D. Brown and D. Peterseim) Institut für Angewandte und Numerische Mathematik Karlsruher Institut
More informationISOGEOMETRIC COLLOCATION METHOD FOR ELASTICITY PROBLEMS CONTAINING SINGULARITIES. Puja Rattan
ISOGEOMETRIC COLLOCATION METHOD FOR ELASTICITY PROBLEMS CONTAINING SINGULARITIES by Puja Rattan A dissertation proposal submitted to the faculty of The University of North Carolina at Charlotte in partial
More informationMatrix assembly by low rank tensor approximation
Matrix assembly by low rank tensor approximation Felix Scholz 13.02.2017 References Angelos Mantzaflaris, Bert Juettler, Boris Khoromskij, and Ulrich Langer. Matrix generation in isogeometric analysis
More informationMapping Techniques for Isogeometric Analysis of Elliptic Boundary Value Problems Containing Singularities
Mapping Techniques for Isogeometric Analysis of Elliptic Boundary Value Problems Containing Singularities Jae Woo Jeong Department of Mathematics, Miami University, Hamilton, OH 450 Hae-Soo Oh Department
More informationFunctional type error control for stabilised space-time IgA approximations to parabolic problems
www.oeaw.ac.at Functional type error control for stabilised space-time IgA approximations to parabolic problems U. Langer, S. Matculevich, S. Repin RICAM-Report 2017-14 www.ricam.oeaw.ac.at Functional
More informationIsogeometric Analysis for the Fractional Laplacian
Isogeometric Analysis for the Fractional Laplacian Kailai Xu & Eric Darve CME500 May 14, 2018 1 / 26 Outline Fractional Calculus Isogeometric Analysis Numerical Experiment Regularity for the Fractional
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 informationarxiv: v1 [math.na] 7 Jan 2019
Efficient p-multigrid Methods for Isogeometric Analysis R. Tielen a,, M. Möller a, D. Göddeke b, C. Vuik a arxiv:9.685v [math.na] 7 Jan 29 a Delft Institute of Applied Mathematics, Delft University of
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 informationA posteriori error estimation in the FEM
A posteriori error estimation in the FEM Plan 1. Introduction 2. Goal-oriented error estimates 3. Residual error estimates 3.1 Explicit 3.2 Subdomain error estimate 3.3 Self-equilibrated residuals 3.4
More informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More informationNotes for Math 663 Spring Alan Demlow
Notes for Math 663 Spring 2016 Alan Demlow 2 Chapter 1 Introduction and preliminaries 1.1 Introduction In this chapter we will introduce the main topics of the course. We first briefly define the finite
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 informationImplementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationA C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes. Marco Verani
A C 1 virtual element method for the Cahn-Hilliard equation with polygonal meshes Marco Verani MOX, Department of Mathematics, Politecnico di Milano Joint work with: P. F. Antonietti (MOX Politecnico di
More informationThe variational collocation method
The variational collocation method Hector Gomez a,, Laura De Lorenzis b a Departamento de Métodos Matemáticos, Universidade da Coruña, Campus de A Coruña, 15071, A Coruña, Spain. b Institut für Angewandte
More informationThe Generalized Empirical Interpolation Method: Analysis of the convergence and application to data assimilation coupled with simulation
The Generalized Empirical Interpolation Method: Analysis of the convergence and application to data assimilation coupled with simulation Y. Maday (LJLL, I.U.F., Brown Univ.) Olga Mula (CEA and LJLL) G.
More informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
More informationAn introduction to the mathematical theory of finite elements
Master in Seismic Engineering E.T.S.I. Industriales (U.P.M.) Discretization Methods in Engineering An introduction to the mathematical theory of finite elements Ignacio Romero ignacio.romero@upm.es October
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
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 informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
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 informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
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 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 informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
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 informationMath 5520 Homework 2 Solutions
Math 552 Homework 2 Solutions March, 26. Consider the function fx) = 2x ) 3 if x, 3x ) 2 if < x 2. Determine for which k there holds f H k, 2). Find D α f for α k. Solution. We show that k = 2. The formulas
More informationAn interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes
An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type interpolation
More informationApproximation of Geometric Data
Supervised by: Philipp Grohs, ETH Zürich August 19, 2013 Outline 1 Motivation Outline 1 Motivation 2 Outline 1 Motivation 2 3 Goal: Solving PDE s an optimization problems where we seek a function with
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationGeometric Decomposition of Serendipity Finite Element Spaces
Geometric Decomposition of Serendipity Finite Element Spaces Andrew Gillette Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/ Joint Mathematics Meetings 2013
More informationISOGEOMETRIC ANALYSIS OF FREE VIBRATION OF TRUSSES AND PLANE STRESS
XXXVI Ibero-La n American Congress on Computa onal Methods in Engineering Rio de Janeiro, 22-25 Nov ISOGEOMETRIC ANALYSIS OF FREE VIBRATION OF TRUSSES AND PLANE STRESS Mateus Rauen Roberto Dalledone Machado
More informationSparse Tensor Galerkin Discretizations for First Order Transport Problems
Sparse Tensor Galerkin Discretizations for First Order Transport Problems Ch. Schwab R. Hiptmair, E. Fonn, K. Grella, G. Widmer ETH Zürich, Seminar for Applied Mathematics IMA WS Novel Discretization Methods
More informationAn axisymmetric PIC code based on isogeometric analysis
Author manuscript, published in "CEMRACS 2010, Marseille : France (2010" DOI : 10.1051/proc/2011016 An axisymmetric PIC code based on isogeometric analysis A. Back A. Crestetto A. Ratnani E. Sonnendrücker
More informationProblem Set 6: Solutions Math 201A: Fall a n x n,
Problem Set 6: Solutions Math 201A: Fall 2016 Problem 1. Is (x n ) n=0 a Schauder basis of C([0, 1])? No. If f(x) = a n x n, n=0 where the series converges uniformly on [0, 1], then f has a power series
More informationOptimal Error Estimates for the hp Version Interior Penalty Discontinuous Galerkin Finite Element Method
Report no. 03/06 Optimal Error Estimates for the hp Version Interior Penalty Discontinuous Galerkin Finite Element Method Emmanuil H. Georgoulis and Endre Süli Oxford University Computing Laboratory Numerical
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
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 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 informationDiese Arbeit wurde vorgelegt am Lehrstuhl für computergestützte Analyse technischer Systeme
Diese Arbeit wurde vorgelegt am Lehrstuhl für computergestützte Analyse technischer Systeme Dreidimensionale Isogeometrische Elastische Festkörper für Fluid-Struktur Interaktionen Three-Dimensional Isogeometric
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 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 informationVirtual Element Methods for general second order elliptic problems
Virtual Element Methods for general second order elliptic problems Alessandro Russo Department of Mathematics and Applications University of Milano-Bicocca, Milan, Italy and IMATI-CNR, Pavia, Italy Workshop
More informationThe Proper Generalized Decomposition: A Functional Analysis Approach
The Proper Generalized Decomposition: A Functional Analysis Approach Méthodes de réduction de modèle dans le calcul scientifique Main Goal Given a functional equation Au = f It is possible to construct
More informationHigh order Galerkin appoximations for parametric second order elliptic partial differential equations
Research Collection Report High order Galerkin appoximations for parametric second order elliptic partial differential equations Author(s): Nistor, Victor; Schwab, Christoph Publication Date: 2012 Permanent
More information2.3 Variational form of boundary value problems
2.3. VARIATIONAL FORM OF BOUNDARY VALUE PROBLEMS 21 2.3 Variational form of boundary value problems Let X be a separable Hilbert space with an inner product (, ) and norm. We identify X with its dual X.
More informationMulti-level hp-fem and the Finite Cell Method for the Navier-Stokes equations using a Variational Multiscale Formulation
Department of Civil, Geo and Environmental Engineering Chair for Computation in Engineering Prof. Dr. rer. nat. Ernst Rank Multi-level hp-fem and the Finite Cell Method for the Navier-Stokes equations
More informationDIFFERENTIAL EQUATIONS
NUMERICAL PARTIAL DIFFERENTIAL EQUATIONS Lecture Notes, Winter 2013/14 Christian Clason July 15, 2014 Institute for Mathematics and Scientific Computing Karl-Franzens-Universität Graz CONTENTS I BACKGROUND
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationA Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form
A Primal-Dual Weak Galerkin Finite Element Method for Second Order Elliptic Equations in Non-Divergence Form Chunmei Wang Visiting Assistant Professor School of Mathematics Georgia Institute of Technology
More informationNUMERICAL PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL PARTIAL DIFFERENTIAL EQUATIONS Lecture Notes, Winter 2011/12 Christian Clason January 31, 2012 Institute for Mathematics and Scientific Computing Karl-Franzens-Universität Graz CONTENTS I BACKGROUND
More informationJOHANNES KEPLER UNIVERSITY LINZ. Institute of Computational Mathematics
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics Multipatch Space-Time Isogeometric Analysis of Parabolic Diffusion Problems Ulrich Langer Institute of Computational Mathematics,
More informationMimetic Finite Difference methods
Mimetic Finite Difference methods An introduction Andrea Cangiani Università di Roma La Sapienza Seminario di Modellistica Differenziale Numerica 2 dicembre 2008 Andrea Cangiani (IAC CNR) mimetic finite
More informationA BIVARIATE SPLINE METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS IN NON-DIVERGENCE FORM
A BIVARIAE SPLINE MEHOD FOR SECOND ORDER ELLIPIC EQUAIONS IN NON-DIVERGENCE FORM MING-JUN LAI AND CHUNMEI WANG Abstract. A bivariate spline method is developed to numerically solve second order elliptic
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 information