The Virtual Element Method: an introduction with focus on fluid flows
|
|
- Nathaniel Gardner
- 5 years ago
- Views:
Transcription
1 The Virtual Element Method: an introduction with focus on fluid flows L. Beirão da Veiga Dipartimento di Matematica e Applicazioni Università di Milano-Bicocca VEM people (or group representatives): P. Antonietti, B. Ayuso, S. Berrone, S. Bordas, F. Brezzi, A. Cangiani, R. Falk, G. Gatica, C. Lovadina, G. Manzini, D. Marini, D. Mora, G. Paulino, P. Pietra, R. Rodriguez, S. Scacchi, A. Russo, G. Vacca, M. Verani, P. Wriggers,... L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
2 The Virtual Element Method The Virtual Element Method (VEM) is a generalization of the Finite Element Method. VEM allow to use very general polygonal and polyhedral meshes, also for high polynomial degrees and without overwhelming complications. The flexibility of VEM is not limited to the mesh: an example will be shown later. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
3 Why polygons/polyhedrons? The interest (and use in commercial codes ) for polytopes is growing. Immediate combination of tets and hexahedrons Easier/better meshing of domain and data features (cracks, sliding meshes, inclusions...) Automatic use of hanging nodes Adaptivity: more efficient mesh refinement/coarsening Generate meshes with more local rotational simmetries Robustness to distortion... for example CD-ADAPCO and ANSYS. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
4 Virtual elements: current range of problems Poisson problem (initial paper) Diffusion-convection-reaction in primal and mixed form Stokes problem Plate bending problems Linear elasticity, nonlinear elasticity and inelasticity Discrete fracture networks Cahn-Hilliard equation Contact problems Eigenvalue problems Topology optimization Time dependent problems... L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
5 The model problem We consider the Poisson problem in two dimensions u = f in Ω, u = 0 on Ω, where Ω R 2 is a polygonal domain; the loading f is assumed in L 2 (Ω). Variational formulation: where a(v, w) = find u V := H0 1 (Ω) such that a(u, v) = f v dx v V, Ω Ω v w dx, v, w V. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
6 A Virtual Element Method [BdV,Brezzi,Cangiani,Manzini,Marini,Russo, M3AS 2013] We will build a discrete problem in following form { find uh V h such that where a h (u h, v h ) =< f h, v h > v h V h, V h V is a finite dimensional space; a h (, ) : V h V h R is a discrete bilinear form approximating the continuous form a(, ); < f h, v h > is a right hand side term approximating the load L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
7 A Virtual Element Method [BdV,Brezzi,Cangiani,Manzini,Marini,Russo, M3AS 2013] We will build a discrete problem in following form { find uh V h such that where a h (u h, v h ) =< f h, v h > v h V h, V h V is a finite dimensional space; a h (, ) : V h V h R is a discrete bilinear form approximating the continuous form a(, ); < f h, v h > is a right hand side term approximating the load Let m 1 be a fixed integer index. Such index will represent the degree of accuracy of the method. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
8 The local spaces V h E Let Ω h be a simple polygonal mesh on Ω. This can be any decomposition of Ω in non overlapping polygons E with straight faces. The space V h will be defined element-wise, by introducing local spaces V h E ; the associated local degrees of freedom. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
9 The local spaces V h E Let Ω h be a simple polygonal mesh on Ω. This can be any decomposition of Ω in non overlapping polygons E with straight faces. The space V h will be defined element-wise, by introducing local spaces V h E ; the associated local degrees of freedom. For all E Ω h : V h E = { v H 1 (E) : v P m 2 (E), v e P m (e) e E }. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
10 The local spaces V h E Let Ω h be a simple polygonal mesh on Ω. This can be any decomposition of Ω in non overlapping polygons E with straight faces. The space V h will be defined element-wise, by introducing local spaces V h E ; the associated local degrees of freedom. For all E Ω h : V h E = { v H 1 (E) : v P m 2 (E), v e P m (e) e E }. the functions v V h E are continuous (and known) on E; the functions v V h E are unknown inside the element E! it holds P m (E) V h E L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
11 Depiction of the degrees of freedom for V h E m=1 m=2 m=3 Green dots stand for vertex pointwise values Red squares represent edge pointwise values (m 2 per edge) Blue squares represent internal (volume) moments v h p m 2 p m 2 P m 2 (E). E L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
12 The space V h The global space V h is built by assembling the local spaces V h E as usual: V h = {v H 1 0 (Ω) : v E V h E E Ω h } The total d.o.f.s are one per internal vertex, m 1 per internal edge and m(m 1)/2 per element. The choice of degrees of freedom guarantees the global continuity of the functions in V h. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
13 The bilinear form a h (, ) The bilinear form a h (, ) is built element by element a h (v h, w h ) = E Ω h a E h (v h, w h ) v h, w h V h, where a E h (, ) : V h E V h E R are symmetric bilinear forms that mimic a E h (, ) a(, ) E by satisfying stability and consistency conditions. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
14 A projection operator We introduce the following energy projector Π : V h E P m (E) defined by, for all v h V h E, { a E (Πv h, p) = a E (v h, p) p P m (E)/R P 0 (Πv h ) = P 0 (v h ) (projection on constants). L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
15 A projection operator We introduce the following energy projector Π : V h E P m (E) defined by, for all v h V h E, { a E (Πv h, p) = a E (v h, p) p P m (E)/R P 0 (Πv h ) = P 0 (v h ) (projection on constants). An integration by parts gives a E (p, v h ) = p v h dx E = ( p) v h dx + E E ( p n E ) v h ds. for all p P m (E), v h V h E. Thus Π is COMPUTABLE. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
16 The local bilinear forms It is immediate to check that v h, w h V h E a E (v h, w h ) = a E (Πv h, Πw h ) + a E ((I Π)v h, (I Π)w h ). L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
17 The local bilinear forms It is immediate to check that v h, w h V h E a E (v h, w h ) = a E (Πv h, Πw h ) + a E ((I Π)v h, (I Π)w h ). We introduce the discrete bilinear form ah E (v h, w h ) = a E (Πv h, Πw h ) + s E ((I Π)v h, (I Π)w h ) where the symmetric bilinear form s E : V h E V h E R is positive definite; scales (in h) like the original bilinear form s E (v h, v h ) a E (v h, v h ) v h V h E, Πv h = 0. NOTE: the above discrete form is exact for polynomials in P m. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
18 The discrete load term We consider m 2 first. Let, for all E Ω h, the approximated load f h E be the L 2 -projection of f E on P m 2 (E). Then (f h, v h ) h := E Ω h E f h v h dx that is computable due to the internal dofs of V h E. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
19 The discrete load term We consider m 2 first. Let, for all E Ω h, the approximated load f h E be the L 2 -projection of f E on P m 2 (E). Then (f h, v h ) h := E Ω h E f h v h dx that is computable due to the internal dofs of V h E. In the case m = 1 a simple integration rule based on the vertex values of the polygon can be used, for instance (f h, v h ) h := ( ) 1 f dx v h (ν) dx. N E E Ω h E ν E Note: for m = 2 better choices can be made [Ahmad, Alsaedi, Brezzi, Marini, Russo, CMA 2013], [BdV, Brezzi, Marini, SINUM 2013]. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
20 A convergence result Let the sequence {Ω h } h satisfy the following mesh assumptions: each element E in Ω h is star-shaped with respect to a ball of uniform radius ( or suitable union of ) ; for each element E in Ω h, the lenght of all edges is comparable with its diameter h E ( not needed by paying (1 + log(he /h e )) ). Then the following holds [... volley team..., M3AS 2013]. Theorem Let the stability and consistency assumptions hold. Then, if f H s (Ω h ) and u H s+1 (Ω h ), we have u u h H 1 (Ω) C ( hs u H s+1 (Ω h ) + f ) H s (Ω h ) for 0 s m and with C independent of h. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
21 k-plain div(k u)+αu = f exact solution: ue(x, y) =y x +log y 3 + x +1 xy xy 2 + x 2 y + x 2 + x 3 +sin(5x) sin(7y) 1 diffusion: K(x, y) =1 zero-order term: α(x, y) =1 right-hand-side: f(x, y) = log y 3 + x +1 1 y 5 x + (y 3 + x +1) y 4 (y 3 + x +1) 2 xy xy2 + x 2 y + x 2 + x 3 6 y +75sin(5x) sin(7y) 3 y 3 + x +1 L2 norm of the exact solution: ue 0,Ω = H1 seminorm of the exact solution: ue 1,Ω = Generated by genue on 09-Nov :08:59
22 0 h=hmean slope: Π uh u e 0,Ω = 4.011, Π u h u e 0,Ω = h 5 h 4 h 3 h 2 0 Π uh u e 0,Ω Π u h u e 0,Ω 0 h=hmean slope: Π uh u e 0,Ω = 3.010, Π u h u e 0,Ω = h 5 h 4 h 3 h 2 0 Π uh u e 0,Ω Π u h u e 0,Ω 10 6 h h L 2 error k =3 H 1 error triangle polygons, h max = , h mean = triangle polygons, h max = , h mean = triangle polygons, h max = , h mean = triangle polygons, h max = , h mean = triangle polygons, h max = , h mean =
23 0 h=hmean slope: Π uh u = 5.009, Π u u = e 0,Ω h e 0,Ω h=hmean slope: Π uh u e 0,Ω = 4.003, Π u h u e 0,Ω = Π 0 u h u e 0,Ω Π 0 u h u e 0,Ω 10 5 h 5 Π u h u e 0,Ω 10 3 h 5 Π u h u e 0,Ω 10 6 h h 4 h 3 h h 2 h h 1 h L 2 error k =4 H 1 error triangle polygons, h max = , h mean = triangle polygons, h max = , h mean = triangle polygons, h max = , h mean = triangle polygons, h max = , h mean = triangle polygons, h max = , h mean =
24 0 h=hmean slope: Π uh u e 0,Ω = 5.928, Π u h u e 0,Ω = h=hmean slope: Π uh u e 0,Ω = 5.020, Π u h u e 0,Ω = Π 0 u h u e 0,Ω Π 0 u h u e 0,Ω 10 6 Π u h u e 0,Ω 10 5 h 5 Π u h u e 0,Ω 10 7 h h h 4 h h 3 h h h h L 2 error k =5 H 1 error triangle polygons, h max = , h mean = triangle polygons, h max = , h mean = triangle polygons, h max = , h mean = triangle polygons, h max = , h mean = triangle polygons, h max = , h mean =
25 Convection-diffusion: an example We consider the problem on Ω = [0, 1] 2 ε u + b u = 0 in Ω, u(0, y) = (1 cos (2πy))/2 y [0, 1], u(1, y) = 0 y [0, 1], homogeneous Neumann BCs in remaining boundary, where ε = 10 2 and b = ( 1, 1/3). The data on the left edge will be transported (with some amplitude loss) and generate a layer near the right edge. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
26 Convection-diffusion: an example Triangular mesh with VEM: joined L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
27 Convection-diffusion: an example Solution plot (polynomial degree one): VEM solution k= L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
28 The three dimensional case In principle the 3D case is analogous; the degrees of freedom are one point value per vertex (m 1) point values per edge M0 f, Mf 1,.., Mf m 1 moments per face M0 E, ME 1,.., ME m 2 moments per element. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
29 The three dimensional case In principle the 3D case is analogous; the degrees of freedom are one point value per vertex (m 1) point values per edge M0 f, Mf 1,.., Mf m 1 moments per face M0 E, ME 1,.., ME m 2 moments per element. A more efficient formulation with less degrees of freedom per face M0 f, Mf 1,.., Mf m 2 moments per face can be built using the L 2 projector in [Ahmed, Alsaedi, Brezzi, Marini, Russo, CMA 2013]. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
30 The Stokes problem Find (u, p) such that ν u p = f in Ω, div u = 0 in Ω, u = 0 on Ω, VARIATIONAL FORM: Find (u, p) V Q, such that ν u : v + p div v = f v for all v V, Ω Ω Ω q div u = 0 for all q Q, Ω where the velocity and pressure spaces V = [H 1 0 (Ω)]2, Q = L 2 0 (Ω). L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
31 A first VEM approximation where Find (u h, p h ) V h Q h, such that a h (u h, v h ) + p h div v h = (f, v) h for all v h V h, Ω q h div u h = 0 for all q h Q h, Ω V h = [V h ] 2 (of degree m), Q h = { q L 2 0 (Ω) : q } E P m 1 E Ω h, a h (v h, w h ) ν v h : w h (vector version of the previous one!) Ω NOTE: The second form is computable (integrate by parts); the inf-sup condition holds! (m 2) L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
32 A new virtual velocity space (m 2) Let W h = { v [H0 1 } (Ω)]2 : v E W h E E Ω h where W h E := { v [H 1 (E)] 2 : v e [P m (e)] 2 e E, { ν v s [Pm 2 (E)] 2, s L 2 0 (E), div v P m 1 (E). The functions in W h E are solution of a local Stokes problem with polynomial data; in particular, the divergence of functions in W h E are polynomials of degree m 1 (known explicitly); the dimension of W h E is the same as that of V h E. (!!) L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
33 Degrees of freedom We take as degrees of freedom for W h E : vertex pointwise values; (m-2) pointwise values for every edge e E; L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
34 Degrees of freedom We take as degrees of freedom for W h E : vertex pointwise values; (m-2) pointwise values for every edge e E; divergence moments (div v h ) p m 1 p m 1 P m 1 (E)/R; E L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
35 Degrees of freedom We take as degrees of freedom for W h E : vertex pointwise values; (m-2) pointwise values for every edge e E; divergence moments (div v h ) p m 1 p m 1 P m 1 (E)/R; E orthogonal moments (x = [ x 2, x 1 ] ) v h x p m 3 p m 3 P m 3 (E). E L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
36 New VEM approximation Let [BdV, Lovadina, Vacca, ArXiv and submitted] Find (u h, p h ) W h Q h, such that a h (u h, v h ) + p h div v h = (f, v) h for all v h W h, Ω q h div u h = 0 for all q h Q h, Ω By using that [P m 2 ] 2 = P m 1 x P m 3, it can be shown that we can still compute v h : p m for all p m [P m (E)] 2, E and thus build the form a h (, ) (loading term analogous). L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
37 Advantages The second row of the discrete problem implies p m 1 div u h = 0 E Ω h, p m 1 P m 1 (E). E Moreover we remind that THEREFORE: (divu h ) E P m 1 (E) E Ω h. the solution u h is pointwise divergence free; the second row of the problem implies, in particular, that all the divergence moments degrees of freedom of u h vanish! L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
38 Reduced problem Let the reduced problem Find (u h, p h ) Ŵh Q h, such that a h (u h, v h ) + p h div v h = (f, v) h for all v h Ŵh, Ω q h div u h = 0 for all q h Q h, where Ω Ŵ h W h with all the divergence moment DoFs set to zero; Q h are piecewise constant functions. NOTE: the solution u h is the same as that of the original problem but many degrees of freedom (both for velocities and pressures) are saved. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
39 A convergence result (m 2) Let the sequence {Ω h } h satisfy the following mesh assumptions: each element E in Ω h is star-shaped with respect to a ball of uniform radius; for each element E in Ω h, the lenght of all edges is comparable with its diameter h E. Then the following holds [LBdV, Lovadina, Vacca, ArXiv and submitted]. Theorem If the right hand side makes sense, we have u u h H 1 (Ω) + p p h L 2 (Ω) C ( hs u H s+1 (Ω h ) + f ) H s (Ω h ) + p H s (Ω h ) for 0 s m and with C independent of h. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
40 Numerical tests TEST: Ω = [0, 1] 2 and loading in accordance with ( ) 1 u(x, y) = 2 cos2 (x) cos(y) sin(y) 1 p(x, y) = sin(x) sin(y). 2 cos2 (y) cos(x) sin(x) We consider the error quantities δ(u) := ui u h 1,h u I and δ(p) := pi p h L 2 1,h p I. L 2 and the mesh families L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
41 Numerical tests: Voronoi meshes L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
42 Numerical tests: triangular meshes L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
43 Numerical tests: quad meshes L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
44 Conclusions The Virtual Element Method is a novel generalization of FEM (2013) that allows for general polytopal meshes The flexibility of the method is reflected also in other aspects, such as the divergence-free scheme shown Virtual elements need (more development and) to approach and be tested on more complex problems (Navier-Stokes,...) and benchmarks We believe the idea is valuable and opens up various interesting research directions both for theoretical and applicative people. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
45 Construction of the stiffness matrix We introduce the following energy projector defined by, for all v h V h E, Π : V h E P m (E) { a E (Πv h, p) = a E (v h, p) p P m (E)/R P 0 (Πv h ) = P 0 (v h ) The operator P 0 is a projection on constants, that for m 2 is simply the average, and for m = 1 the vertex value average. NOTE: due to the consistency assumption, the projection operator above is computable (more later). L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
46 Construction of the stiffness matrix It is immediate to check that v h, w h V h E a E (v h, w h ) = a E (Πv h, Πw h ) + a E ((I Π)v h, (I Π)w h ). Then the bilinear form a E h (v h, w h ) = a E (Πv h, Πw h ) + s E ((I Π)v h, (I Π)w h ) is consistent and stable, provided the positive bilinear form s E : V h E V h E R scales like the original bilinear form a. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
47 Construction of the stiffness matrix It is immediate to check that v h, w h V h E a E (v h, w h ) = a E (Πv h, Πw h ) + a E ((I Π)v h, (I Π)w h ). Then the bilinear form a E h (v h, w h ) = a E (Πv h, Πw h ) + s E ((I Π)v h, (I Π)w h ) is consistent and stable, provided the positive bilinear form s E : V h E V h E R scales like the original bilinear form a. Let now the local stiffness matrix M E M E ij := a E h (ϕ i, ϕ j ) i, j = 1, 2,..., N with {ϕ i } the canonical basis associated to the degrees of freedom. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
48 Construction of the stiffness matrix We introduce also {m α } n α=1 a basis for the polynomial space P m = span{m α } n α=1. Since P m V h E, the matrix dof 1 (m 1 ) dof 1 (m 2 )... dof 1 (m n ) dof 2 (m 1 ) dof 2 (m 2 )... dof 2 (m n ) D = dof N (m 1 ) dof N (m 2 )... dof N (m n ) expresses the {m α } in terms of the V h E basis. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
49 Construction of the stiffness matrix We have also a matrix P 0 ϕ 1... P 0 ϕ N a E (ϕ 1, m 2 )... a E (ϕ N, m 2 ) B =..... a E (ϕ 1, m n )... a E (ϕ N, m n ) expressing (in terms of the bases) the right hand side in the definition { a E (Πv h, p) = a E (v h, p) p P m (E)/R P 0 (Πv h ) = P 0 (v h ) NOTE: such matrix is computable! (integration by parts and d.o.f.s definition) L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
50 Construction of the stiffness matrix Let P 0 m 1 ( P 0 m 2... P 0 m n 0 m2, m 2 )0,P... ( ) m2, m n G = BD =..... (. 0 mn, m 2 )0,P... ( ) mn, m n 0,P 0,P. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
51 Construction of the stiffness matrix Let P 0 m 1 ( P 0 m 2... P 0 m n 0 m2, m 2 )0,P... ( ) m2, m n G = BD =..... (. 0 mn, m 2 )0,P... ( ) mn, m n Compute the matrices corresponding to the projection operator: Π = (G) 1 B, Π = DΠ. 0,P 0,P. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
52 Construction of the stiffness matrix Let P 0 m 1 ( P 0 m 2... P 0 m n 0 m2, m 2 )0,P... ( ) m2, m n G = BD =..... (. 0 mn, m 2 )0,P... ( ) mn, m n Compute the matrices corresponding to the projection operator: Π = (G) 1 B, Π = DΠ. Finally compute the local stiffness matrix M E = Π T G Π + (I Π) T (I Π). 0,P 0,P. L. Beirão da Veiga (Univ. of Milan-Bicocca) The Virtual Element Method PARIS / 31
arxiv: 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 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 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 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 informationDivergence free Virtual Elements for the Stokes problem on polygonal meshes
Divergence free Virtual Elements for the Stokes problem on polygonal meshes L. Beirão da Veiga, C. Lovadina and G. Vacca Abstract arxiv:1510.01655v1 [math.na] 6 Oct 2015 In the present paper we develop
More informationH(div) and H(curl)-conforming Virtual Element Methods
Noname manuscript No. (will be inserted by the editor) H(div) and H(curl)-conforming Virtual Element Methods L. Beirão da Veiga F. Brezzi L. D. Marini A. Russo the date of receipt and acceptance should
More informationVEM s for the numerical solution of PDE s. F. Brezzi
VEM s for the numerical solution of PDE s F. Brezzi IUSS-Istituto Universitario di Studi Superiori, Pavia, Italy IMATI-C.N.R., Pavia, Italy From Joint works with L. Beirão da Veiga, L.D. Marini, and A.
More informationarxiv: v1 [math.na] 11 Sep 2017
Dual virtual element method in presence of an inclusion Alessio Fumagalli arxiv:1709.03519v1 [math.na] 11 Sep 2017 Abstract We consider a Darcy problem for saturated porous media written in dual formulation
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 informationVirtual Elements for plates and shells
Virtual Elements for plates and shells Claudia Chinosi Dipartimento di Scienze e Innovazione Tecnologica, Università del Piemonte Orientale, Alessandria, Italia E-mail: claudia.chinosi@unipmn.it eywords:
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 informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationFinite Element Multigrid Framework for Mimetic Finite Difference Discretizations
Finite Element Multigrid Framework for Mimetic Finite ifference iscretizations Xiaozhe Hu Tufts University Polytopal Element Methods in Mathematics and Engineering, October 26-28, 2015 Joint work with:
More informationTHE NONCONFORMING VIRTUAL ELEMENT METHOD
THE NONCONFORMING VIRTUAL ELEMENT METHOD BLANCA AYUSO DE DIOS, ONSTANTIN LIPNIOV, AND GIANMARCO MANZINI Abstract. We introduce the nonconforming Virtual Element Method (VEM) for the approximation of second
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 informationFebruary 18, :41 WSPC/INSTRUCTION FILE Brezzi-Corrections VIRTUAL ELEMENT METHOD FOR GENERAL SECOND ORDER ELLIPTIC PROBLEMS ON POLYGONAL MESHES
February 8, 206 7:4 WSPC/INSTRUCTION FIL Brezzi-Corrections Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company VIRTUAL LMNT MTHOD FOR GNRAL SCOND ORDR LLIPTIC PROBLMS
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 informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for Kirchhoff plate elements Jarkko Niiranen Laboratory of Structural
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES On the approximation of a virtual coarse space for domain decomposition methods in two dimensions Juan Gabriel Calvo Preprint No. 31-2017 PRAHA 2017
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
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 informationBACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigid-body spring deformation expressed in terms of field variables assembly of rigid-bodies connected by spring Distinct Element Method (DEM) simple
More informationarxiv: v1 [math.na] 8 Jan 2019
arallel solvers for virtual element discretizations of elliptic equations in mixed form F. Dassi a, S. Scacchi b a Dipartimento di Matematica, Università degli Studi di Milano-Bicocca, Via Cozzi, 2133
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More informationLibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu and Roy H. Stogner roystgnr@cfdlab.ae.utexas.edu Univ. of Texas at Austin September 9, 2008 1 Introduction 2 Weighted Residuals
More informationarxiv: v1 [math.na] 11 Feb 2015
A C 1 VIRTUAL LMNT MTHOD FOR TH CAHN-HILLIARD QUATION WITH POLYGONAL MSHS PAOLA F. ANTONITTI, LOURNCO BIRÃO DA VIGA, SIMON SCACCHI, AND MARCO VRANI arxiv:1502.03259v1 [math.na] 11 Feb 2015 Abstract. In
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 informationA Virtual Element Method for the Transmission Eigenvalue Problem
A Virtual Element Method for the Transmission Eigenvalue Problem David Mora a,b, Iván Velásquez b,c a GIMNAP, Departamento de Matemática, Universidad del Bío-Bío, Casilla 5-C, Concepción, Chile. arxiv:183.1979v1
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES A virtual overlapping Schwarz method for scalar elliptic problems in two dimensions Juan Gabriel Calvo Preprint No. 25-2017 PRAHA 2017 A VIRTUAL
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
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 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 informationDiscontinuous Petrov-Galerkin Methods
Discontinuous Petrov-Galerkin Methods Friederike Hellwig 1st CENTRAL School on Analysis and Numerics for Partial Differential Equations, November 12, 2015 Motivation discontinuous Petrov-Galerkin (dpg)
More informationICES REPORT A Generalized Mimetic Finite Difference Method and Two-Point Flux Schemes over Voronoi Diagrams
ICS RPORT 15-17 July 2015 A Generalized Mimetic Finite Difference Method and Two-Point Flux Schemes over Voronoi Diagrams by Omar Al-Hinai, Mary F. Wheeler, Ivan Yotov The Institute for Computational ngineering
More informationLibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu Univ. of Texas at Austin January 12, 2007 1 Introduction 2 Weighted Residuals 3 Poisson Equation 4 Other Examples 5 Essential
More informationFinite element approximation of eigenvalue problems
Acta NumeTica (2010), pp. 1-120 doi: 10.1017/80962492910000012 Cambridge University Press, 2010 Finite element approximation of eigenvalue problems Daniele Boffi Dipartimento di Matematica, Universita
More informationarxiv: v1 [math.na] 8 Feb 2018
arxiv:180.094v1 [math.na] 8 Feb 018 The nonconforming virtual element method for eigenvalue problems F. Gardini a, G. Manzini b, and G. Vacca c a Dipartimento di Matematica F. Casorati, Università di Pavia,
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 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 informationarxiv: v1 [math.na] 3 Jun 2018
SUPG stabilization for the nonconforming virtual element method for advection-diffusion-reaction equations S Berrone a, A Borio a, G Manzini b,c arxiv:186879v1 [mathna] 3 Jun 18 Abstract a Dipartimento
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 informationComputer simulation of multiscale problems
Progress in the SSF project CutFEM, Geometry, and Optimal design Computer simulation of multiscale problems Axel Målqvist and Daniel Elfverson University of Gothenburg and Uppsala University Umeå 2015-05-20
More informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationA P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS
A P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS SHANGYOU ZHANG DEDICATED TO PROFESSOR PETER MONK ON THE OCCASION OF HIS 6TH BIRTHDAY Abstract. On triangular grids, the continuous
More informationTMA4220: Programming project - part 1
TMA4220: Programming project - part 1 TMA4220 - Numerical solution of partial differential equations using the finite element method The programming project will be split into two parts. This is the first
More informationarxiv: v1 [math.na] 21 Dec 2018
The p- and hp-versions of the virtual element method for elliptic eigenvalue problems O. Čertík, F. Gardini, G. Manzini, L. Mascotto, G. Vacca arxiv:8.090v [math.na] Dec 08 Abstract We discuss the p- and
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 informationSolutions of Selected Problems
1 Solutions of Selected Problems October 16, 2015 Chapter I 1.9 Consider the potential equation in the disk := {(x, y) R 2 ; x 2 +y 2 < 1}, with the boundary condition u(x) = g(x) r for x on the derivative
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 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 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 informationMATH 676. Finite element methods in scientific computing
MATH 676 Finite element methods in scientific computing Wolfgang Bangerth, Texas A&M University Lecture 33.25: Which element to use Part 2: Saddle point problems Consider the stationary Stokes equations:
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 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 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 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 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 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 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 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 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 informationA $C^1$ virtual element method for the Cahn-Hilliard equation with polygonal meshes
MOX-Report No. 12/2015 A $C^1$ virtual element method for the Cahn-Hilliard equation with polygonal meshes Antonietti, P. F.; Beirao da Veiga, L.; Scacchi, S.; Verani, M. MOX, Dipartimento di Matematica
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 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 informationComputational Methods for Engineering Applications II. Homework Problem Sheet 4
S. Mishra K. O. Lye L. Scarabosio Autumn Term 2015 Computational Methods for Engineering Applications II ETH Zürich D-MATH Homework Problem Sheet 4 This assignment sheet contains some tasks marked as Core
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 informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
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 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 informationLinear algebra II Homework #1 due Thursday, Feb A =
Homework #1 due Thursday, Feb. 1 1. Find the eigenvalues and the eigenvectors of the matrix [ ] 3 2 A =. 1 6 2. Find the eigenvalues and the eigenvectors of the matrix 3 2 2 A = 2 3 2. 2 2 1 3. The following
More informationH(div) and H(curl)-conforming VEM
H(div) and H(curl)-conforming VM L. Beirão da Veiga, F. Brezzi, L. D. Marini, and A. Russo arxiv:1407.6822v1 [math.na] 25 Jul 2014 Abstract In the present paper we construct Virtual lement Spaces that
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationBASIC PRINCIPLES OF VIRTUAL ELEMENT METHODS
Mathematical Models and Methods in Applied Sciences Vol. 23, No. 1 (2013) 199 214 c World Scientific Publishing Company DOI: 10.1142/S0218202512500492 BASIC PRINCIPLES OF VIRTUAL ELEMENT METHODS L. BEIRÃO
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
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 informationSome New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
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 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 informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationPart D: Frames and Plates
Part D: Frames and Plates Plane Frames and Thin Plates A Beam with General Boundary Conditions The Stiffness Method Thin Plates Initial Imperfections The Ritz and Finite Element Approaches A Beam with
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
MATHEMATICS OF COMPUTATION Volume 75, Number 256, October 2006, Pages 1659 1674 S 0025-57180601872-2 Article electronically published on June 26, 2006 ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED
More informationInstitute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I
Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix
More informationAnalysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems
Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems Christian Waluga 1 advised by Prof. Herbert Egger 2 Prof. Wolfgang Dahmen 3 1 Aachen Institute for Advanced Study in Computational
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 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 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 informationSome remarks on the stability coefficients and bubble stabilization of FEM on anisotropic meshes
Some remarks on the stability coefficients and bubble stabilization of FEM on anisotropic meshes Stefano Micheletti Simona Perotto Marco Picasso Abstract In this paper we re-address the anisotropic recipe
More informationA non-standard Finite Element Method based on boundary integral operators
A non-standard Finite Element Method based on boundary integral operators Clemens Hofreither Ulrich Langer Clemens Pechstein June 30, 2010 supported by Outline 1 Method description Motivation Variational
More informationHigher-Order Compact Finite Element Method
Higher-Order Compact Finite Element Method Major Qualifying Project Advisor: Professor Marcus Sarkis Amorn Chokchaisiripakdee Worcester Polytechnic Institute Abstract The Finite Element Method (FEM) is
More informationInstitute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I
Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix
More informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
More informationQUADRILATERAL H(DIV) FINITE ELEMENTS
QUADRILATERAL H(DIV) FINITE ELEMENTS DOUGLAS N. ARNOLD, DANIELE BOFFI, AND RICHARD S. FALK Abstract. We consider the approximation properties of quadrilateral finite element spaces of vector fields defined
More informationA Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators
A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationRobust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations
Robust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with
More informationThe Mortar Boundary Element Method
The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract
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 informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
More information