Introduction to finite element exterior calculus
|
|
- Angelina Stone
- 6 years ago
- Views:
Transcription
1 Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway
2 Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl H(div, Ω) div L 2 (Ω) 0
3 Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl H(div, Ω) div L 2 (Ω) 0 and the corresponding elasticity complex T H 1 (Ω; R 3 ) ɛ H(J, Ω; S) J H(div, Ω; S) div L 2 (Ω; R 3 ) 0. Here J is the second order operator Jτ = curl(curl τ) T mapping symmetric matrix fields into symmetric matrix fields, and T is the six dimensional space of rigid motions.
4 Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl H(div, Ω) div L 2 (Ω) 0 and the corresponding elasticity complex T H 1 (Ω; R 3 ) ɛ H(J, Ω; S) J H(div, Ω; S) div L 2 (Ω; R 3 ) 0. Here J is the second order operator Jτ = curl(curl τ) T mapping symmetric matrix fields into symmetric matrix fields, and T is the six dimensional space of rigid motions. Reference: Arnold, Falk, Winther, Mixed finite elements for linear elasticity with weakly imposed symmetry, Math.Comp
5 Mapping of the domain Consider a map φ : R 3 R 3 mapping a domain Ω to Ω. Ω ϕ Ω
6 Mapping of the domain Consider a map φ : R 3 R 3 mapping a domain Ω to Ω. Ω ϕ Ω Desired commuting diagram: grad curl R H 1 (Ω ) H(curl, Ω ) H(div, Ω ) F 1 φ R H 1 (Ω) F c φ grad H(curl, Ω) F d φ curl H(div, Ω) div L 2 (Ω ) 0 F 0 φ div L 2 (Ω) 0
7 Piola transforms The proper Piola transforms are given by (Fφ 1 u)(x) = u(φ(x)) (F c φ u)(x) = (Dφ(x))T u(φ(x) (F d φ u)(x) = (det Dφ(x))(Dφ(x)) 1 u(φ(x)) (Fφ 0 u)(x) = (det Dφ(x))u(φ(x))
8 Alternative de Rham complex We replace the de Rham complex R H 1 (Ω) grad H(curl, Ω) curl H(div, Ω) div L 2 (Ω) 0 by an equivalent complex R HΛ 0 (Ω) d HΛ 1 (Ω) d HΛ 2 (Ω) d L 2 Λ 3 (Ω) 0 In fact, the second complex can be defined for any dimension n.
9 Exterior calculus and differential forms Instead of consider scalar fields and vector fields, i.e., functions of the form u : Ω R or u : Ω R 3 where Ω R 3, we consider instead functions of the form u : Ω Alt k R n Here Ω is allowed to be a subset of R n.
10 Exterior calculus and differential forms Instead of consider scalar fields and vector fields, i.e., functions of the form u : Ω R or u : Ω R 3 where Ω R 3, we consider instead functions of the form u : Ω Alt k R n Here Ω is allowed to be a subset of R n. For any vector space V, the space Alt k V is the space of alternating k-linear maps on V
11 The space Alt k V Assume dim V = n and that v 1, v 2,... v n is a basis. If ω Alt k V and j 1, j 2,..., j k are integers in {1, 2,..., n} then ω(v j1, v j2,..., v jk ) R.
12 The space Alt k V Assume dim V = n and that v 1, v 2,... v n is a basis. If ω Alt k V and j 1, j 2,..., j k are integers in {1, 2,..., n} then ω(v j1, v j2,..., v jk ) R. Furthermore, ω is linear in each argument, and if two neighboring vectors change position then the value of ω change sign. Therefore, it is enough to specify ω(v j1, v j2,..., v jk ) for an increasing sequence {j i }.
13 The space Alt k V Assume dim V = n and that v 1, v 2,... v n is a basis. If ω Alt k V and j 1, j 2,..., j k are integers in {1, 2,..., n} then ω(v j1, v j2,..., v jk ) R. Furthermore, ω is linear in each argument, and if two neighboring vectors change position then the value of ω change sign. Therefore, it is enough to specify ω(v j1, v j2,..., v jk ) for an increasing sequence {j i }. As a consequence Alt k V is a vector space with ( ) n dim Alt k V =. k
14 The space Alt k V Assume dim V = n and that v 1, v 2,... v n is a basis. If ω Alt k V and j 1, j 2,..., j k are integers in {1, 2,..., n} then ω(v j1, v j2,..., v jk ) R. Furthermore, ω is linear in each argument, and if two neighboring vectors change position then the value of ω change sign. Therefore, it is enough to specify ω(v j1, v j2,..., v jk ) for an increasing sequence {j i }. As a consequence Alt k V is a vector space with ( ) n dim Alt k V =. k Note that dim Alt k V = dim Alt n k V. The canonical map between these spaces is the Hodge star operator.
15 Interior and exterior products If ω Alt k V and v V then the interior product, or the contraction of ω by v, ω v Alt k 1 V is given by ω v(v 1, v 2,... v k 1 ) = ω(v, v 1, v 2,... v k 1 )
16 Interior and exterior products If ω Alt k V and v V then the interior product, or the contraction of ω by v, ω v Alt k 1 V is given by ω v(v 1, v 2,... v k 1 ) = ω(v, v 1, v 2,... v k 1 ) If ω Alt j V and η Alt k V then the exterior product, or wedge product, ω η Alt j+k is given by (ω η)(v 1,..., v j+k ) = σ (sign σ)ω(v σ(1),..., v σ(j) )η(v σ(j+1),..., v σ(j+k) ), v i V, where the sum is over all permutations σ of {1, 2,..., j + k} such that σ(1) <... < σ(j) and σ(j + 1) <... < σ(j + k).
17 Basis for Alt k R n We let {dx i } be the natural dual basis in Alt 1 R n given by dx i (e j ) = δ i,j
18 Basis for Alt k R n We let {dx i } be the natural dual basis in Alt 1 R n given by dx i (e j ) = δ i,j A basis for Alt 2 R n is given by dx i dx j, i < j,
19 Basis for Alt k R n We let {dx i } be the natural dual basis in Alt 1 R n given by dx i (e j ) = δ i,j A basis for Alt 2 R n is given by dx i dx j, i < j, while dx j1 dx j2 dx jk, where {j i } is increasing is a basis for Alt k R n.
20 Any ω Alt k R n can be represented uniquely on the form u σ dx σ(1) dx σ(2)... dx σ(k) σ where the sum is over all incresing maps σ mapping {1, 2,..., k} into {1, 2,..., n}, and u σ R.
21 Any ω Alt k R n can be represented uniquely on the form u σ dx σ(1) dx σ(2)... dx σ(k) σ where the sum is over all incresing maps σ mapping {1, 2,..., k} into {1, 2,..., n}, and u σ R. Furthermore, dx σ(1) dx σ(2)... dx σ(k) (v 1,... v k ) = det(dx σ(i) v j )
22 Correspondences Alt k and vector fields: Alt 0 R 3 = R Alt 1 R 3 = R 3 Alt 2 R 3 = R 3 c c u 1 dx 1 + u 2 dx 2 + u 3 dx 3 u Alt 3 R 3 = R c c dx1 dx 2 dx 3 u 3 dx 1 dx 2 u 2 dx 1 dx 3 + u 1 dx 2 dx 3 u
23 Correspondences Alt k and vector fields: Alt 0 R 3 = R Alt 1 R 3 = R 3 Alt 2 R 3 = R 3 c c u 1 dx 1 + u 2 dx 2 + u 3 dx 3 u Alt 3 R 3 = R c c dx1 dx 2 dx 3 Exterior product: u 3 dx 1 dx 2 u 2 dx 1 dx 3 + u 1 dx 2 dx 3 u : Alt 1 R 3 Alt 1 R 3 Alt 2 R 3 : R 3 R 3 R 3 : Alt 1 R 3 Alt 2 R 3 Alt 3 R 3 : R 3 R 3 R
24 Differential forms Differential forms are functions of the form ω : Ω Alt k Alt k R n.
25 Differential forms Differential forms are functions of the form ω : Ω Alt k Alt k R n. We let Λ k = Λ k (Ω) = C (Ω; Alt k ). For such smooth differentiable forms the exterior derivative, d : Λ k Λ k+1 is defined by k+1 dω x (v 1, v 2,..., v k+1 ) = ( 1) j+1 vj ω x (v 1,..., ˆv j,..., v k+1 ), j=1 ω Λ k, v 1,..., v k+1 R n
26 Differential forms Differential forms are functions of the form ω : Ω Alt k Alt k R n. We let Λ k = Λ k (Ω) = C (Ω; Alt k ). For such smooth differentiable forms the exterior derivative, d : Λ k Λ k+1 is defined by k+1 dω x (v 1, v 2,..., v k+1 ) = ( 1) j+1 vj ω x (v 1,..., ˆv j,..., v k+1 ), j=1 ω Λ k, v 1,..., v k+1 R n Using the fact that vi vj = vj vi it follows that d d = 0.
27 Check of d d = 0 k+1 dω x (v 1, v 2,..., v k+1 ) = ( 1) j+1 vj ω x (v 1,..., ˆv j,..., v k+1 ), j=1 ω Λ k, v 1,..., v k+1 R n Start with a zero form ω.
28 Check of d d = 0 k+1 dω x (v 1, v 2,..., v k+1 ) = ( 1) j+1 vj ω x (v 1,..., ˆv j,..., v k+1 ), j=1 ω Λ k, v 1,..., v k+1 R n Start with a zero form ω. dω x (v) = v ω x
29 Check of d d = 0 k+1 dω x (v 1, v 2,..., v k+1 ) = ( 1) j+1 vj ω x (v 1,..., ˆv j,..., v k+1 ), j=1 ω Λ k, v 1,..., v k+1 R n Start with a zero form ω. dω x (v) = v ω x (d d)ω(v 1, v 2 ) = v1 dω(v 2 ) v2 dω(v 1 ) = 0.
30 Representation by a basis If ω is given on the form ω x = σ u σ dx σ(1) dx σ(2)... dx σ(k), where each u σ is a scalar field on Ω then dω x = σ n j=1 u σ x j dx j dx σ(1) dx σ(2)... dx σ(k)
31 The de Rham complex and differential forms Using differential forms the de Rham complex (smooth version) can be written as R Λ 0 d Λ 1 d d Λ n 0.
32 Mapping of the domain Ω ϕ Ω
33 Mapping of the domain Ω ϕ Ω The corresponding pull backs, φ, which maps Λ k (Ω ) to Λ k (Ω), is given by (φ ω) x (v 1, v 2,..., v k ) = ω φ(x) (Dφ x (v 1 ), Dφ x (v 2 ),..., Dφ x (v k )), where Dφ x is the derivative of φ at x mapping T x Ω into T φ(x) Ω.
34 Mapping of the domain Ω ϕ Ω The corresponding pull backs, φ, which maps Λ k (Ω ) to Λ k (Ω), is given by (φ ω) x (v 1, v 2,..., v k ) = ω φ(x) (Dφ x (v 1 ), Dφ x (v 2 ),..., Dφ x (v k )), where Dφ x is the derivative of φ at x mapping T x Ω into T φ(x) Ω. This expression should be compared with the expressions we had for the Piola maps before.
35 Commuting diagram R Λ 0 (Ω ) R Λ 0 (Ω) d Λ 1 (Ω ) d d Λ n (Ω ) 0. φ φ φ d Λ 1 (Ω) d d Λ n (Ω) 0.
36 Commuting diagram R Λ 0 (Ω ) R Λ 0 (Ω) d Λ 1 (Ω ) d d Λ n (Ω ) 0. φ φ φ d Λ 1 (Ω) d d Λ n (Ω) 0. So the pullback commutes with the exterior derivative, i.e., φ (dω) = d(φ ω), ω Λ k (Ω ), and distributes with respect to the wedge product: φ (ω η) = φ ω φ η.
37 Commuting diagram R Λ 0 (Ω ) R Λ 0 (Ω) d Λ 1 (Ω ) d d Λ n (Ω ) 0. φ φ φ d Λ 1 (Ω) d d Λ n (Ω) 0. So the pullback commutes with the exterior derivative, i.e., φ (dω) = d(φ ω), ω Λ k (Ω ), and distributes with respect to the wedge product: φ (ω η) = φ ω φ η. Note this also gives us a way to define the exterior derivative on a manifold.
38 Commuting diagram R Λ 0 (Ω ) R Λ 0 (Ω) d Λ 1 (Ω ) d d Λ n (Ω ) 0. φ φ φ d Λ 1 (Ω) d d Λ n (Ω) 0. So the pullback commutes with the exterior derivative, i.e., φ (dω) = d(φ ω), ω Λ k (Ω ), and distributes with respect to the wedge product: φ (ω η) = φ ω φ η. Note this also gives us a way to define the exterior derivative on a manifold. Key identity: (ψ φ) = φ ψ
39 Exterior derivative and wedge product By combining the wedge product with the exterior derivative we obtain a Leibniz rule of the form d(ω η) = dω η + ( 1) j ω dη, ω Λ j, η Λ k.
40 Exterior derivative and wedge product By combining the wedge product with the exterior derivative we obtain a Leibniz rule of the form d(ω η) = dω η + ( 1) j ω dη, ω Λ j, η Λ k. The pullback of the inclusion Ω Ω is the trace map, Tr. The vanishing of Tr ω for ω Λ k (Ω) does not imply that ω x vanishes for x Ω, only that it vanishes when applied to k-tuples of tangent vectors to Ω.
41 Exterior derivative and wedge product By combining the wedge product with the exterior derivative we obtain a Leibniz rule of the form d(ω η) = dω η + ( 1) j ω dη, ω Λ j, η Λ k. The pullback of the inclusion Ω Ω is the trace map, Tr. The vanishing of Tr ω for ω Λ k (Ω) does not imply that ω x vanishes for x Ω, only that it vanishes when applied to k-tuples of tangent vectors to Ω. Since Tr is a pull back we will have Tr d = d Ω Tr
42 Integration of differential forms If ω Λ n (Ω), Ω R n, with the representation ω x = u x dx 1 dx 2... dx n, then we can define ω = u dx Ω Ω
43 Integration of differential forms If ω Λ n (Ω), Ω R n, with the representation ω x = u x dx 1 dx 2... dx n, then we can define ω = If φ : Ω Ω then we have φ ω = Ω Ω Ω u dx Ω ω
44 Integration of differential forms If ω Λ n (Ω), Ω R n, with the representation then we can define If φ : Ω Ω then we have ω x = u x dx 1 dx 2... dx n, Ω Ω ω = Ω φ ω = u dx Ω ω Note that this gives us a way to define the integral of an n-form on a manifold.
45 Stokes theorem With the three concepts exterior derivative, trace operator, and integral we can formulate Stokes theorem as dω = Tr ω, ω Λ n 1 (Ω) Ω Ω
46 Stokes theorem With the three concepts exterior derivative, trace operator, and integral we can formulate Stokes theorem as dω = Tr ω, ω Λ n 1 (Ω) Ω Ω Furthermore, by combining this with the Leibniz rule we obtain the following integration by parts identity dω η = ( 1) k 1 ω dη + Tr ω Tr η Ω where ω Λ k, η Λ n k 1. Ω Ω
47 L 2 differential forms If ω is a k form given by ω = σ u σ dx σ(1) dx σ(2)... dx σ(k) then ω is defined to be in L 2 Λ k (Ω) if each of the coefficients u σ are in L 2.
48 L 2 differential forms If ω is a k form given by ω = σ u σ dx σ(1) dx σ(2)... dx σ(k) then ω is defined to be in L 2 Λ k (Ω) if each of the coefficients u σ are in L 2. Furthernmore, we let HΛ k (Ω) = { ω L 2 Λ k (Ω) dω L 2 Λ k+1 (Ω) }
49 L 2 differential forms If ω is a k form given by ω = σ u σ dx σ(1) dx σ(2)... dx σ(k) then ω is defined to be in L 2 Λ k (Ω) if each of the coefficients u σ are in L 2. Furthernmore, we let HΛ k (Ω) = { ω L 2 Λ k (Ω) dω L 2 Λ k+1 (Ω) } We then ontain the L 2 de Rham complex: 0 HΛ 0 (Ω) d HΛ 1 (Ω) d d HΛ n (Ω) 0
50 Polynomial de Rham complex We define P r Λ k = {ω Λ k ω(v 1,... v k ) P r, v 1,... v k }
51 Polynomial de Rham complex We define P r Λ k = {ω Λ k ω(v 1,... v k ) P r, v 1,... v k } Alternatively, P r Λ k consists of all u σ dx σ(1) dx σ(2)... dx σ(k) σ where the coefficients u σ are in P r.
52 Polynomial de Rham complex We define P r Λ k = {ω Λ k ω(v 1,... v k ) P r, v 1,... v k } Alternatively, P r Λ k consists of all u σ dx σ(1) dx σ(2)... dx σ(k) σ where the coefficients u σ are in P r. The polynomial de Rham complex is exact. 0 P r Λ 0 d Pr 1 Λ 1 d d Pr n Λ n 0
53 The Kozul complex The Koszul differential κ of a k-form ω is the (k 1)-form given by (κω) x (v 1,..., v k 1 ) = ω x x(v 1,..., v k 1 ) = ω x ( x, v1,..., v k 1 ),
54 The Kozul complex The Koszul differential κ of a k-form ω is the (k 1)-form given by (κω) x (v 1,..., v k 1 ) = ω x x(v 1,..., v k 1 ) = ω x ( x, v1,..., v k 1 ), For each r, κ maps P r 1 Λ k to P r Λ k 1,
55 The Kozul complex The Koszul differential κ of a k-form ω is the (k 1)-form given by (κω) x (v 1,..., v k 1 ) = ω x x(v 1,..., v k 1 ) = ω x ( x, v1,..., v k 1 ), For each r, κ maps P r 1 Λ k to P r Λ k 1, and the Koszul complex 0 P r n Λ n κ Pr n+1 Λ n 1 κ κ Pr Λ 0 R 0, is exact.
56 The spaces P r Λ k The operators d and κ are related by the homotopy relation (dκ + κd)ω = (r + k)ω, ω H r Λ k, where H r denotes the homogeneous polynomials of degree r.
57 The spaces P r Λ k The operators d and κ are related by the homotopy relation (dκ + κd)ω = (r + k)ω, ω H r Λ k, where H r denotes the homogeneous polynomials of degree r. As a consequence we obtain the identity P r Λ k = P r 1 Λ k + κh r 1 Λ k+1 + dh r+1 Λ k 1
58 The spaces P r Λ k The operators d and κ are related by the homotopy relation (dκ + κd)ω = (r + k)ω, ω H r Λ k, where H r denotes the homogeneous polynomials of degree r. As a consequence we obtain the identity P r Λ k = P r 1 Λ k + κh r 1 Λ k+1 + dh r+1 Λ k 1 We then define the space P r Λ k by P r Λ k = P r 1 Λ k + κp r 1 Λ k+1.
59 The spaces P r Λ k The operators d and κ are related by the homotopy relation (dκ + κd)ω = (r + k)ω, ω H r Λ k, where H r denotes the homogeneous polynomials of degree r. As a consequence we obtain the identity P r Λ k = P r 1 Λ k + κh r 1 Λ k+1 + dh r+1 Λ k 1 We then define the space P r Λ k by P r Λ k = P r 1 Λ k + κp r 1 Λ k+1. We note that P r Λ 0 = P r Λ 0 and P r Λ n = P r 1 Λ n. Furthermore, 0 P r Λ 0 d P r Λ 1 d d P r Λ n 0 is an exact complex, and the space P r Λ k is affine invariant.
60 Characterization of affine invariant polynomial spaces Theorem Assume that the polynomial space QΛ k, satisfies P r 1 Λ k QΛ k P r Λ k. Then QΛ k = P r Λ k or QΛ k = P r 1 Λ k + dp r+1 Λ k 1
61 The four exact sequences ending with P r Λ 3 (T ) in 3D 0 P r+1 Λ 0 d P r+1 Λ 1 d P r+1 Λ 2 d Pr Λ P r+2 Λ 0 d Pr+1 Λ 1 d P r+1 Λ 2 d Pr Λ P r+2 Λ 0 d P r+2 Λ 1 d Pr+1 Λ 2 d Pr Λ P r+3 Λ 0 d Pr+2 Λ 1 d Pr+1 Λ 2 d Pr Λ 3 0
62 The four exact sequences ending with P r Λ 3 (T ) in 3D 0 P r+1 Λ 0 d P r+1 Λ 1 d P r+1 Λ 2 d Pr Λ P r+2 Λ 0 d Pr+1 Λ 1 d P r+1 Λ 2 d Pr Λ P r+2 Λ 0 d P r+2 Λ 1 d Pr+1 Λ 2 d Pr Λ P r+3 Λ 0 d Pr+2 Λ 1 d Pr+1 Λ 2 d Pr Λ 3 0 In general there are 2 n 1 exact sequences composed of the families of P and P spaces.
63 The four sequences ending with P 0 Λ 3 (T ) in 3D 0 grad curl div 0 0 grad curl div 0 0 grad curl div 0 0 grad curl div 0
64 Piecewise smooth differential forms It is a consequence of Stokes theorem that a piecewise smooth k form ω, with respect to a simplicial mesh T h of Ω, is in HΛ k (Ω) if and only if the trace of ω, Tr ω, is continuous on the interfaces. Here Tr ω is defined by restricting the spatial variable x to the interface, and by applying ω only to tangent vectors of the interface.
65 Degrees of freedom To obtain finite element differential forms not just pw polynomials we need degrees of freedom, i.e., a decomposition of the dual spaces (P r Λ k (T )) and (P r Λ k (T )) (with T a simplex), into subspaces associated to subsimplices f of T.
66 Degrees of freedom To obtain finite element differential forms not just pw polynomials we need degrees of freedom, i.e., a decomposition of the dual spaces (P r Λ k (T )) and (Pr Λ k (T )) (with T a simplex), into subspaces associated to subsimplices f of T. DOF for P r Λ k (T ): to a subsimplex f of dimension d we associate ω Tr f ω η, η P r+k d Λd k (f ) f
67 Degrees of freedom To obtain finite element differential forms not just pw polynomials we need degrees of freedom, i.e., a decomposition of the dual spaces (P r Λ k (T )) and (Pr Λ k (T )) (with T a simplex), into subspaces associated to subsimplices f of T. DOF for P r Λ k (T ): to a subsimplex f of dimension d we associate ω Tr f ω η, η P r+k d Λd k (f ) f DOF for Pr Λ k (T ): ω Tr f ω η, η P r+k d 1 Λ d k (f ) f
68 Degrees of freedom To obtain finite element differential forms not just pw polynomials we need degrees of freedom, i.e., a decomposition of the dual spaces (P r Λ k (T )) and (Pr Λ k (T )) (with T a simplex), into subspaces associated to subsimplices f of T. DOF for P r Λ k (T ): to a subsimplex f of dimension d we associate ω Tr f ω η, η P r+k d Λd k (f ) f DOF for Pr Λ k (T ): ω Tr f ω η, η P r+k d 1 Λ d k (f ) f Given a triangulation T, we can then define P r Λ k (T ), P r Λ k (T ). They are subspaces of HΛ k (Ω).
69 De Rham cohomology 0 Λ 0 d (Ω) Λ 1 d (Ω) Λ 2 (Ω) 0 0 C (Ω) grad C (Ω; R 2 ) rot C (Ω) 0
70 Cohomology The de Rham complex HΛ k 1 (Ω) is called exact if for all k, d k 1 HΛ k (Ω) d k HΛ k+1 (Ω) Z k := ker d k = range d k 1 =: B k.
71 Cohomology The de Rham complex HΛ k 1 (Ω) is called exact if for all k, d k 1 HΛ k (Ω) d k HΛ k+1 (Ω) Z k := ker d k = range d k 1 =: B k. In general, B k Z k and we assume throughout that the kth cohomology group Z k /B k is finite dimensional.
72 Cohomology The de Rham complex HΛ k 1 (Ω) is called exact if for all k, d k 1 HΛ k (Ω) d k HΛ k+1 (Ω) Z k := ker d k = range d k 1 =: B k. In general, B k Z k and we assume throughout that the kth cohomology group Z k /B k is finite dimensional. The space of harmonic k-forms, H k, consists of all q Z k such that q, µ = 0 µ B k.
73 Cohomology The de Rham complex HΛ k 1 (Ω) is called exact if for all k, d k 1 HΛ k (Ω) d k HΛ k+1 (Ω) Z k := ker d k = range d k 1 =: B k. In general, B k Z k and we assume throughout that the kth cohomology group Z k /B k is finite dimensional. The space of harmonic k-forms, H k, consists of all q Z k such that q, µ = 0 µ B k. This leads to the Hodge decomposition HΛ k (Ω) = Z k Z k = B k H k Z k. Note that H k = Z k /B k.
74 Hodge Laplace problem Λ k 1 (Ω) d k 1 Λ k (Ω) d k Λ k+1 (Ω)
75 Hodge Laplace problem Λ k 1 (Ω) d k 1 Λ k (Ω) d k Λ k+1 (Ω) Formally: Given f Λ k, find u Λ k such that (d k 1 δ k 1 + δ k d k )u = f. Here δ k is a formal adjoint of d k.
76 Hodge Laplace problem Λ k 1 (Ω) d k 1 Λ k (Ω) d k Λ k+1 (Ω) Formally: Given f Λ k, find u Λ k such that (d k 1 δ k 1 + δ k d k )u = f. Here δ k is a formal adjoint of d k. The following mixed formulation is always well-posed: Given f L 2 Λ k (Ω), find σ HΛ k 1, u HΛ k and p H k such that σ, τ dτ, u = 0 τ HΛ k 1 dσ, v + du, dv + p, v = f, v v HΛ k u, q = 0 q H k
77 Hodge Laplacian Well-posedness of the Hodge Laplace problem follows from the Hodge decomposition and Poincaré s inequality: ω L 2 c dω L 2, ω (Z k ). Special cases: Recall that if dim H k = 0: Find σ HΛ k 1, and u HΛ k σ, τ dτ, u = 0 τ HΛ k 1 dσ, v + du, dv = f, v v HΛ k
78 Hodge Laplacian Well-posedness of the Hodge Laplace problem follows from the Hodge decomposition and Poincaré s inequality: ω L 2 c dω L 2, ω (Z k ). Special cases: Recall that if dim H k = 0: Find σ HΛ k 1, and u HΛ k σ, τ dτ, u = 0 τ HΛ k 1 dσ, v + du, dv = f, v v HΛ k k = 0: ordinary Laplacian k = n: mixed Laplacian k = 1, n = 3: σ = div u, grad σ + curl curl u = f k = 2, n = 3: σ = curl u, curl σ grad div u = f,
79 Discretizations of the de Rham complex in three dimensions R H 1 (Ω) I 1 h grad H(curl, Ω) I c h curl H(div, Ω) I d h div L 2 (Ω) 0 I 0 h R H 1 h grad H h (curl) curl H h (div) div L 2 h 0.
80 Discretizations of the de Rham complex in three dimensions R H 1 (Ω) I 1 h grad H(curl, Ω) I c h curl H(div, Ω) I d h div L 2 (Ω) 0 I 0 h R H 1 h grad H h (curl) curl H h (div) div L 2 h 0. The operators I h are the canonical interpolation operators defined from the degrees of freedom. They commute with grad curl and div, but the are not bounded on the spaces H 1, H(curl) or H(div).
81 Discretizations of the de Rham complex in n dimensions R HΛ 0 (Ω) R d HΛ 1 (Ω) d d HΛ n (Ω) 0. I h I h I h Λ 0 h d Λ 1 h d d Λ n h 0.
82 Discretizations of the de Rham complex in n dimensions R HΛ 0 (Ω) R d HΛ 1 (Ω) d d HΛ n (Ω) 0. I h I h I h Λ 0 h d Λ 1 h d d Λ n h 0. The operators I h are the canonical interpolation operators defined from the degrees of freedom. They commute with d, but the are not bounded on the HΛ.
83 Discretization, Abstract setting Λ k 1 d k 1 Λ k Complex of Hilbert spaces with d k bounded and closed range.
84 Discretization, Abstract setting Λ k 1 d k 1 Λ k Complex of Hilbert spaces with d k bounded and closed range. Hodge decomposition and Poincaré inequality follow.
85 Discretization, Abstract setting Λ k 1 d k 1 Λ k Λ k 1 h d k 1 Λ k h Complex of Hilbert spaces with d k bounded and closed range. Hodge decomposition and Poincaré inequality follow. For discretization, construct a finite dimensional subcomplex.
86 Discretization, Abstract setting Λ k 1 d k 1 Λ k Λ k 1 h d k 1 Λ k h Complex of Hilbert spaces with d k bounded and closed range. Hodge decomposition and Poincaré inequality follow. For discretization, construct a finite dimensional subcomplex. Define H k h = (Bk h ) Z k h.
87 Discretization, Abstract setting Λ k 1 d k 1 Λ k Λ k 1 h d k 1 Λ k h Complex of Hilbert spaces with d k bounded and closed range. Hodge decomposition and Poincaré inequality follow. For discretization, construct a finite dimensional subcomplex. Define H k h = (Bk h ) Z k h. Discrete Hodge decomposition follows: Λ k h = Bk h Hk h (Zk h )
88 Discretization, Abstract setting Λ k 1 d k 1 Λ k Λ k 1 h d k 1 Λ k h Complex of Hilbert spaces with d k bounded and closed range. Hodge decomposition and Poincaré inequality follow. For discretization, construct a finite dimensional subcomplex. Define H k h = (Bk h ) Z k h. Discrete Hodge decomposition follows: Λ k h = Bk h Hk h (Zk h ) Galerkin s method: Λ k 1, Λ k, H k Λ k 1 h, Λ k h, Hk h When is it stable?
89 Bounded cochain projections Key property: Suppose that there exists a bounded cochain projection. Λ k 1 d k 1 Λ k π k 1 h Λ k 1 h π k h d k 1 Λ k h π k h uniformly bounded πh k ω ω 0. π k h a projection π k h d k 1 = d k 1 π k 1 h
90 Bounded cochain projections Key property: Suppose that there exists a bounded cochain projection. Λ k 1 d k 1 Λ k π k 1 h Λ k 1 h Theorem π k h d k 1 Λ k h π k h uniformly bounded πh k ω ω 0. π k h a projection π k h d k 1 = d k 1 π k 1 h If v π k h v < v v Hk, then the induced map on cohomology is an isomorphism. gap ( H k, H k h) sup v πh k v v H k, v =1 The discrete Poincaré inequality holds uniformly in h. Galerkin s method is stable and convergent.
91 Proof of discrete Poincaré inequality Theorem: There is a positive constant c, independent of h, such that ω c dω, ω Z k h. In fact, c c P c π, where c P is continuous Poincaré constant, and c π bounds π h.
92 Proof of discrete Poincaré inequality Theorem: There is a positive constant c, independent of h, such that ω c dω, ω Z k h. In fact, c c P c π, where c P is continuous Poincaré constant, and c π bounds π h. Proof: Given ω Z k h, define η Zk HΛ k (Ω) by dη = dω. By the Poincaré inequality, η c dω, so it is enough to show that ω c η. Now, ω π h η Λ k h and d(ω π hη) = 0, so ω π h η Z k h. Therefore ω 2 = ω, π h η + ω, ω π h η = ω, π h η ω π h η, whence ω π h η. The result follows from the uniform boundedness of π h.
93 Construction of bounded cochain projections The canonical projections, I h, determined by the degrees of freedom, commute with d. But they are not bounded on HΛ k.
94 Construction of bounded cochain projections The canonical projections, I h, determined by the degrees of freedom, commute with d. But they are not bounded on HΛ k. If we apply the three operations: extend (E) regularize (R) canonical projection (I h ) we get a map Qh k : HΛk (Ω) Λ k h which is bounded and commutes with d. But it is not a projection.
95 Construction of bounded cochain projections The canonical projections, I h, determined by the degrees of freedom, commute with d. But they are not bounded on HΛ k. If we apply the three operations: extend (E) regularize (R) canonical projection (I h ) we get a map Qh k : HΛk (Ω) Λ k h which is bounded and commutes with d. But it is not a projection. However the composition π k h = (Qk h Λ k h) 1 Q k h can be shown to be a bounded cochain projection.
96 Bounded projections, definitions We combine extension and smoothing into an operator R ɛ h : L2 Λ k (Ω) CΛ k (Ω) given by (Rh ɛ ω) x = ρ(y)(eω) x+ɛhy dy. B 1
97 Bounded projections, definitions We combine extension and smoothing into an operator R ɛ h : L2 Λ k (Ω) CΛ k (Ω) given by (Rh ɛ ω) x = ρ(y)(eω) x+ɛhy dy. B 1 Assume that the mesh is quasi uniform: Lemma For each ɛ > 0 sufficiently small there is a constant c(ɛ) such that for all h I h R ɛ h L(L 2 Λ k (Ω),L 2 Λ k (Ω)) c(ɛ). There is constant c, independent of ɛ and h such that I I h R ɛ h L(L 2 Λ k h,l2 Λ k h ) cɛ.
98 Scaling and compactness
99 Bounded cochain projections The operator π h = π k h defined by π h = (I h R ɛ h Λ k h) 1 I h R ɛ h, where ɛ is taken small, but not too small, has all the desired properties, i.e. π k h is bounded in L(L2 Λ k (Ω), L 2 Λ k (Ω)) and L(HΛ k (Ω), HΛ k (Ω)) π k h π k h commutes with the exterior derivative is a projection
100 Removing the assumption of a quasi uniform mesh If the family of meshes is shape regular, but not necessarily quasi uniform, we introduce a piecewise linear and uniformly Lipschitz continuous function g h (x) to represent the local mesh size.
101 Removing the assumption of a quasi uniform mesh If the family of meshes is shape regular, but not necessarily quasi uniform, we introduce a piecewise linear and uniformly Lipschitz continuous function g h (x) to represent the local mesh size. For each vertex x we let g h (x) be the average of the values of diam T for the simplices meeting x, and we define Φ ɛy h (x) = x + ɛg h(x)y.
102 Removing the assumption of a quasi uniform mesh If the family of meshes is shape regular, but not necessarily quasi uniform, we introduce a piecewise linear and uniformly Lipschitz continuous function g h (x) to represent the local mesh size. For each vertex x we let g h (x) be the average of the values of diam T for the simplices meeting x, and we define Φ ɛy h (x) = x + ɛg h(x)y. Then we replace the operator (Rh ɛ ω) x = ρ(y)(eω) x+ɛhy dy B 1 by (Rh ɛ ω) x = ρ(y)((φ ɛy h ) Eω) x dy. B 1
103 References The development of FEEC leans heavily on earlier results taken from Whitney, Bossavit, Raviart and Thomas, Nedelec, Hiptmair,... as well as on the theory of finite elements in general.
104 References The development of FEEC leans heavily on earlier results taken from Whitney, Bossavit, Raviart and Thomas, Nedelec, Hiptmair,... as well as on the theory of finite elements in general. The presentation here is mostly based on D.N. Arnold, R.S. Falk, R. Winther, Finite element exterior calculus, homological techniques, and applications, Acta Numerica D.N. Arnold, R.S. Falk and R. Winther Finite element exterior calculus: From Hodge theory to numerical stability Bulletin of the Amer. Math. Soc. 47 (2010),
Smoothed projections in finite element exterior calculus
Smooted projections in finite element exterior calculus Ragnar Winter CMA, University of Oslo Norway based on joint work wit: Douglas N. Arnold, Minnesota, Ricard S. Falk, Rutgers, and Snorre H. Cristiansen,
More informationMotivations. Outline. Finite element exterior calculus and the geometrical basis of numerical stability. References. Douglas N.
Outline Finite element exterior calculus and the geometrical basis of numerical stability Douglas N. Arnold joint with R. Falk and R. Winther School of Mathematics University of Minnesota April 009 Motivations
More informationPublished in Computer Methods in Applied Mechanics and Engineering 198 (2009)
Published in Computer Methods in Applied Mechanics and Engineering 198 (2009) 1660 1672 GEOMETRIC DECOMPOSITIONS AND LOCAL BASES FOR SPACES OF FINITE ELEMENT DIFFERENTIAL FORMS DOUGLAS N. ARNOLD, RICHARD
More informationFinite Element Exterior Calculus. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 7 October 2016
Finite Element Exterior Calculus Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 7 October 2016 The Hilbert complex framework Discretization of Hilbert complexes Finite element
More informationFinite element exterior calculus: A new approach to the stability of finite elements
Finite element exterior calculus: A new approach to the stability of finite elements Douglas N. Arnold joint with R. Falk and R. Winther Institute for Mathematics and its Applications University of Minnesota
More informationSpaces of Finite Element Differential Forms
Spaces of Finite Element Differential Forms Douglas N. Arnold Abstract We discuss the construction of finite element spaces of differential forms which satisfy the crucial assumptions of the finite element
More informationA posteriori error estimates in FEEC for the de Rham complex
A posteriori error estimates in FEEC for the de Rham complex Alan Demlow Texas A&M University joint work with Anil Hirani University of Illinois Urbana-Champaign Partially supported by NSF DMS-1016094
More informationA POSTERIORI ERROR ESTIMATES FOR FINITE ELEMENT EXTERIOR CALCULUS: THE DE RHAM COMPLEX
A POSTERIORI ERROR ESTIMATES FOR FINITE ELEMENT EXTERIOR CALCULUS: THE DE RHAM COMPLEX ALAN DEMLOW 1 AND ANIL N. HIRANI 2 Abstract. Finite element exterior calculus (FEEC) has been developed over the past
More informationAposteriorierrorestimatesinFEEC for the de Rham complex
AposteriorierrorestimatesinFEEC for the de Rham complex Alan Demlow Texas A&M University joint work with Anil Hirani University of Illinois Urbana-Champaign Partially supported by NSF DMS-1016094 and a
More informationFinite element exterior calculus
Finite element exterior calculus Douglas N. Arnold School of Mathematics, University of Minnesota NSF/CBMS Conference, ICERM, June 11 15, 2012 collaborators: R. Falk, R. Winther, G. Awanou, F. Bonizzoni,
More informationNew Model Stability Criteria for Mixed Finite Elements
New Model Stability Criteria for Mixed Finite Elements Andrew Gillette Department of Mathematics Institute of Computational Engineering and Sciences University of Texas at Austin http://www.math.utexas.edu/users/agillette
More informationFinite Element Exterior Calculus and Applications Part V
Finite Element Exterior Calculus and Applications Part V Douglas N. Arnold, University of Minnesota Peking University/BICMR August 15 18, 2015 De Rham complex 0 H 1 (Ω) grad H(curl, Ω) curl H(div, Ω) div
More informationarxiv: v2 [math.na] 8 Sep 2015
OMPLEXES OF DISRETE DISTRIBUTIONAL DIFFERENTIAL FORMS AND THEIR HOMOLOGY THEORY MARTIN WERNER LIHT arxiv:1410.6354v2 [math.na] 8 Sep 2015 Abstract. omplexes of discrete distributional differential forms
More informationAn abstract Hodge Dirac operator and its stable discretization
An abstract Hodge Dirac operator and its stable discretization Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at ICCA, Tartu,. Joint work with Ari Stern,
More informationMath 225B: Differential Geometry, Final
Math 225B: Differential Geometry, Final Ian Coley March 5, 204 Problem Spring 20,. Show that if X is a smooth vector field on a (smooth) manifold of dimension n and if X p is nonzero for some point of
More informationOutline. 1 Motivating examples. 2 Hilbert complexes and their discretization. 3 Finite element differential forms. 4 The elasticity complex 1 / 25
Outline 1 Motivating examples Hoge theory, Hilbert complexes, an finite element ifferential forms Douglas N. Arnol 2 Hilbert complexes an their iscretization 3 Finite element ifferential forms 4 The elasticity
More informationSurvey on exterior algebra and differential forms
Survey on exterior algebra and differential forms Daniel Grieser 16. Mai 2013 Inhaltsverzeichnis 1 Exterior algebra for a vector space 1 1.1 Alternating forms, wedge and interior product.....................
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationDivergence-free or curl-free finite elements for solving the curl-div system
Divergence-free or curl-free finite elements for solving the curl-div system Alberto Valli Dipartimento di Matematica, Università di Trento, Italy Joint papers with: Ana Alonso Rodríguez Dipartimento di
More informationFrom point cloud data to the continuum model of geometry
From point cloud data to the continuum model of geometry J. Harrison University of California, Berkeley July 22, 2007 Continuum model of geometry By the continuum model of geometry we refer to smooth manifolds
More informationDouble Complexes and Local Cochain Projections
Double Complexes and Local Cochain Projections Richard S. Falk, 1 Ragnar Winther 2 1 Department o Mathematics, Rutgers University, Piscataway, New Jersey 08854 2 Centre o Mathematics or Applications and
More informationSolving the curl-div system using divergence-free or curl-free finite elements
Solving the curl-div system using divergence-free or curl-free finite elements Alberto Valli Dipartimento di Matematica, Università di Trento, Italy or: Why I say to my students that divergence-free finite
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationMIXED FINITE ELEMENT APPROXIMATION OF THE VECTOR LAPLACIAN WITH DIRICHLET BOUNDARY CONDITIONS
MIXED FINITE ELEMENT APPROXIMATION OF THE VECTOR LAPLACIAN WITH DIRICHLET BOUNDARY CONDITIONS DOUGLAS N. ARNOLD, RICHARD S. FALK, AND JAY GOPALAKRISHNAN Abstract. We consider the finite element solution
More informationFinite Element Methods for Linear Elasticity
Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu
More informationA Generalization for Stable Mixed Finite Elements
A Generalization for Stable Mixed Finite Elements Andrew Gillette joint work with Chandrajit Bajaj Department of Mathematics Institute of Computational Engineering and Sciences University of Texas at Austin,
More informationEilenberg-Steenrod properties. (Hatcher, 2.1, 2.3, 3.1; Conlon, 2.6, 8.1, )
II.3 : Eilenberg-Steenrod properties (Hatcher, 2.1, 2.3, 3.1; Conlon, 2.6, 8.1, 8.3 8.5 Definition. Let U be an open subset of R n for some n. The de Rham cohomology groups (U are the cohomology groups
More informationMATH 4030 Differential Geometry Lecture Notes Part 4 last revised on December 4, Elementary tensor calculus
MATH 4030 Differential Geometry Lecture Notes Part 4 last revised on December 4, 205 Elementary tensor calculus We will study in this section some basic multilinear algebra and operations on tensors. Let
More informationHungry, Hungry Homology
September 27, 2017 Motiving Problem: Algebra Problem (Preliminary Version) Given two groups A, C, does there exist a group E so that A E and E /A = C? If such an group exists, we call E an extension of
More informationMath 868 Final Exam. Part 1. Complete 5 of the following 7 sentences to make a precise definition (5 points each). Y (φ t ) Y lim
SOLUTIONS Dec 13, 218 Math 868 Final Exam In this exam, all manifolds, maps, vector fields, etc. are smooth. Part 1. Complete 5 of the following 7 sentences to make a precise definition (5 points each).
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationFINITE ENERGY SOLUTIONS OF MIXED 3D DIV-CURL SYSTEMS
FINITE ENERGY SOLUTIONS OF MIXED 3D DIV-CURL SYSTEMS GILES AUCHMUTY AND JAMES C. ALEXANDER Abstract. This paper describes the existence and representation of certain finite energy (L 2 -) solutions of
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 informationNONCONFORMING MIXED ELEMENTS FOR ELASTICITY
Mathematical Models and Methods in Applied Sciences Vol. 13, No. 3 (2003) 295 307 c World Scientific Publishing Company NONCONFORMING MIXED ELEMENTS FOR ELASTICITY DOUGLAS N. ARNOLD Institute for Mathematics
More informationInequalities of Babuška-Aziz and Friedrichs-Velte for differential forms
Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms Martin Costabel Abstract. For sufficiently smooth bounded plane domains, the equivalence between the inequalities of Babuška Aziz
More informationMATH 263: PROBLEM SET 1: BUNDLES, SHEAVES AND HODGE THEORY
MATH 263: PROBLEM SET 1: BUNDLES, SHEAVES AND HODGE THEORY 0.1. Vector Bundles and Connection 1-forms. Let E X be a complex vector bundle of rank r over a smooth manifold. Recall the following abstract
More informationIs a dual mesh really necessary?
Is a dual mesh really necessary? Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at ICIAM, Vancouver, 2011. Joint work with Ari Stern, UCSD. 19 July 2011
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 informationDISCRETE EXTENSION OPERATORS FOR MIXED FINITE ELEMENT SPACES ON LOCALLY REFINED MESHES
DISCRETE EXTENSION OPERATORS FOR MIXED FINITE ELEMENT SPACES ON LOCALLY REFINED MESHES MAR AINSWORTH, JOHNNY GUZMÁN, AND FRANCISCO JAVIER SAYAS Abstract. The existence of uniformly bounded discrete extension
More informationTHE LOCAL THEORY OF ELLIPTIC OPERATORS AND THE HODGE THEOREM
THE LOCAL THEORY OF ELLIPTIC OPERATORS AND THE HODGE THEOREM BEN LOWE Abstract. In this paper, we develop the local theory of elliptic operators with a mind to proving the Hodge Decomposition Theorem.
More informationis the desired collar neighbourhood. Corollary Suppose M1 n, M 2 and f : N1 n 1 N2 n 1 is a diffeomorphism between some connected components
1. Collar neighbourhood theorem Definition 1.0.1. Let M n be a manifold with boundary. Let p M. A vector v T p M is called inward if for some local chart x: U V where U subsetm is open, V H n is open and
More informationSome Elliptic and Subelliptic Complexes from Geometry
Geometric and Nonlinear Partial Differential Equations, Xi An Jiaotong University p. 1/14 Some Elliptic and Subelliptic Complexes from Geometry Michael Eastwood [ based on joint work with Robert Bryant,
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 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 informationA posteriori error estimates for Maxwell Equations
www.oeaw.ac.at A posteriori error estimates for Maxwell Equations J. Schöberl RICAM-Report 2005-10 www.ricam.oeaw.ac.at A POSTERIORI ERROR ESTIMATES FOR MAXWELL EQUATIONS JOACHIM SCHÖBERL Abstract. Maxwell
More informationModeling, Structure and Discretization of Mixed-dimensional Partial Differential Equations
Modeling, Structure and Discretization of Mixed-dimensional Partial Differential Equations J. M. Nordbotten and W. M. Boon Department of Mathematics University of Bergen Norway Abstract Mixed-dimensional
More informationFoundations of finite element methods for wave equations of Maxwell type
Dept. of Math./CMA Univ. of Oslo Pure Mathematics No 5 ISSN 0806 2439 February 2009 Foundations of finite element methods for wave equations of Maxwell type Snorre H. Christiansen Abstract The first part
More informationSection 2. Basic formulas and identities in Riemannian geometry
Section 2. Basic formulas and identities in Riemannian geometry Weimin Sheng and 1. Bianchi identities The first and second Bianchi identities are R ijkl + R iklj + R iljk = 0 R ijkl,m + R ijlm,k + R ijmk,l
More informationSolutions of M3-4A16 Assessed Problems # 3 [#1] Exercises in exterior calculus operations
D. D. Holm Solutions to M3-4A16 Assessed Problems # 3 15 Dec 2010 1 Solutions of M3-4A16 Assessed Problems # 3 [#1] Exercises in exterior calculus operations Vector notation for differential basis elements:
More informationDifferential Forms JWR. May 2, 2005
Differential Forms JWR May 2, 2005 1 Algebra 1.1. For (real) vector spaces V and W denote by L(V, W ) the vector space of linear maps from V to W and by L k (V, W ) the space of k-multilinear maps from
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 informationDifferential Geometry: Discrete Exterior Calculus
Differential Geometry: Discrete Exterior Calculus [Discrete Exterior Calculus (Thesis). Hirani, 23] [Discrete Differential Forms for Computational Modeling. Desbrun et al., 25] [Build Your Own DEC at Home.
More informationk times l times n times
1. Tensors on vector spaces Let V be a finite dimensional vector space over R. Definition 1.1. A tensor of type (k, l) on V is a map T : V V R k times l times which is linear in every variable. Example
More informationNONLOCAL EXTERIOR CALCULUS ON RIEMANNIAN MANIFOLDS
The Pennsylvania State University The Graduate School Department of Mathematics NONLOCAL EXTERIOR CALCULUS ON RIEMANNIAN MANIFOLDS A Dissertation in Mathematics by Thinh Duc Le c 2013 Thinh Duc Le Submitted
More informationHigh Order Differential Form-Based Elements for the Computation of Electromagnetic Field
1472 IEEE TRANSACTIONS ON MAGNETICS, VOL 36, NO 4, JULY 2000 High Order Differential Form-Based Elements for the Computation of Electromagnetic Field Z Ren, Senior Member, IEEE, and N Ida, Senior Member,
More informationFinite Element Clifford Algebra: A New Toolkit for Evolution Problems
Finite Element Clifford Algebra: A New Toolkit for Evolution Problems Andrew Gillette joint work with Michael Holst Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/
More informationMULTILINEAR ALGEBRA MCKENZIE LAMB
MULTILINEAR ALGEBRA MCKENZIE LAMB 1. Introduction This project consists of a rambling introduction to some basic notions in multilinear algebra. The central purpose will be to show that the div, grad,
More informationSpace-Time Finite-Element Exterior Calculus and Variational Discretizations of Gauge Field Theories
21st International Symposium on Mathematical Theory of Networks and Systems July 7-11, 2014. Space-Time Finite-Element Exterior Calculus and Variational Discretizations of Gauge Field Theories Joe Salamon
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationDivergence-conforming multigrid methods for incompressible flow problems
Divergence-conforming multigrid methods for incompressible flow problems Guido Kanschat IWR, Universität Heidelberg Prague-Heidelberg-Workshop April 28th, 2015 G. Kanschat (IWR, Uni HD) Hdiv-DG Práha,
More informationProgress in Several Complex Variables KIAS 2018
for Progress in Several Complex Variables KIAS 2018 Outline 1 for 2 3 super-potentials for 4 real for Let X be a real manifold of dimension n. Let 0 p n and k R +. D c := { C k (differential) (n p)-forms
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 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 informationChap. 1. Some Differential Geometric Tools
Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationDIFFERENTIAL FORMS AND COHOMOLOGY
DIFFERENIAL FORMS AND COHOMOLOGY ONY PERKINS Goals 1. Differential forms We want to be able to integrate (holomorphic functions) on manifolds. Obtain a version of Stokes heorem - a generalization of the
More informationNumerical approximation of output functionals for Maxwell equations
Numerical approximation of output functionals for Maxwell equations Ferenc Izsák ELTE, Budapest University of Twente, Enschede 11 September 2004 MAXWELL EQUATIONS Assumption: electric field ( electromagnetic
More informationIn collaboration with Mathieu Desbrun, Anil N. Hirani, and Jerrold E. Marsden. Abstract
73 Chapter 3 Discrete Exterior Calculus In collaboration with Mathieu Desbrun, Anil N. Hirani, and Jerrold E. Marsden. Abstract We present a theory and applications of discrete exterior calculus on simplicial
More informationSelf-adjoint curl Operators
Eidgenössische Technische Hochschule Zürich Ecole polytechnique fédérale de Zurich Politecnico federale di Zurigo Swiss Federal Institute of Technology Zurich Self-adjoint curl Operators R. Hiptmair and
More informationCopyright by Andrew Kruse Gillette 2011
Copyright by Andrew Kruse Gillette 2011 The Dissertation Committee for Andrew Kruse Gillette certifies that this is the approved version of the following dissertation: Stability of Dual Discretization
More informationCoordinate free non abelian geometry I: the quantum case of simplicial manifolds.
Coordinate free non abelian geometry I: the quantum case of simplicial manifolds. Johan Noldus May 6, 07 Abstract We study the geometry of a simplicial complexes from an algebraic point of view and devise
More informationCohomology of Harmonic Forms on Riemannian Manifolds With Boundary
Cohomology of Harmonic Forms on Riemannian Manifolds With Boundary Sylvain Cappell, Dennis DeTurck, Herman Gluck, and Edward Y. Miller 1. Introduction To Julius Shaneson on the occasion of his 60th birthday
More informationGravitation: Tensor Calculus
An Introduction to General Relativity Center for Relativistic Astrophysics School of Physics Georgia Institute of Technology Notes based on textbook: Spacetime and Geometry by S.M. Carroll Spring 2013
More informationDISCRETE EXTERIOR CALCULUS
DISCRETE EXTERIOR CALCULUS MATHIEU DESBRUN, ANIL N. HIRANI, MELVIN LEOK, AND JERROLD E. MARSDEN Abstract. We present a theory and applications of discrete exterior calculus on simplicial complexes of arbitrary
More informationOrientation transport
Orientation transport Liviu I. Nicolaescu Dept. of Mathematics University of Notre Dame Notre Dame, IN 46556-4618 nicolaescu.1@nd.edu June 2004 1 S 1 -bundles over 3-manifolds: homological properties Let
More informationCHARACTERISTIC CLASSES
1 CHARACTERISTIC CLASSES Andrew Ranicki Index theory seminar 14th February, 2011 2 The Index Theorem identifies Introduction analytic index = topological index for a differential operator on a compact
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationA DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT
A DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT HUOYUAN DUAN, SHA LI, ROGER C. E. TAN, AND WEIYING ZHENG Abstract. To deal with the divergence-free
More informationW if p = 0; ; W ) if p 1. p times
Alternating and symmetric multilinear functions. Suppose and W are normed vector spaces. For each integer p we set {0} if p < 0; W if p = 0; ( ; W = L( }... {{... } ; W if p 1. p times We say µ p ( ; W
More informationHomology of Random Simplicial Complexes Based on Steiner Systems
Homology of Random Simplicial Complexes Based on Steiner Systems Amir Abu-Fraiha Under the Supervision of Prof. Roy Meshulam Technion Israel Institute of Technology M.Sc. Seminar February 2016 Outline
More informationLecture 4 - Lie Algebra Cohomology I
Lecture 4 - Lie Algebra Cohomology I January 25, 2013 Given a differentiable manifold M n and a k-form ω, recall Cartan s formula for the exterior derivative k dω(v 0,..., v k ) = ( 1) i x i (ω(x 0,...,
More informationLecture Notes: African Institute of Mathematics Senegal, January Topic Title: A short introduction to numerical methods for elliptic PDEs
Lecture Notes: African Institute of Mathematics Senegal, January 26 opic itle: A short introduction to numerical methods for elliptic PDEs Authors and Lecturers: Gerard Awanou (University of Illinois-Chicago)
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationPAijpam.eu NEW H 1 (Ω) CONFORMING FINITE ELEMENTS ON HEXAHEDRA
International Journal of Pure and Applied Mathematics Volume 109 No. 3 2016, 609-617 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: 10.12732/ijpam.v109i3.10
More informationDierential geometry for Physicists
Dierential geometry for Physicists (What we discussed in the course of lectures) Marián Fecko Comenius University, Bratislava Syllabus of lectures delivered at University of Regensburg in June and July
More informationLet V, W be two finite-dimensional vector spaces over R. We are going to define a new vector space V W with two properties:
5 Tensor products We have so far encountered vector fields and the derivatives of smooth functions as analytical objects on manifolds. These are examples of a general class of objects called tensors which
More informationCS 468 (Spring 2013) Discrete Differential Geometry
CS 468 (Spring 2013) Discrete Differential Geometry Lecture 13 13 May 2013 Tensors and Exterior Calculus Lecturer: Adrian Butscher Scribe: Cassidy Saenz 1 Vectors, Dual Vectors and Tensors 1.1 Inner Product
More informationMath 797W Homework 4
Math 797W Homework 4 Paul Hacking December 5, 2016 We work over an algebraically closed field k. (1) Let F be a sheaf of abelian groups on a topological space X, and p X a point. Recall the definition
More informationDr. Allen Back. Dec. 3, 2014
Dr. Allen Back Dec. 3, 2014 forms are sums of wedge products of the basis 1-forms dx, dy, and dz. They are kinds of tensors generalizing ordinary scalar functions and vector fields. They have a skew-symmetry
More informationA POSTERIORI ERROR ESTIMATES FOR MAXWELL EQUATIONS
A POSTERIORI ERROR ESTIMATES FOR MAXWELL EQUATIONS JOACHIM SCHÖBERL Abstract. Maxwell equations are posed as variational boundary value problems in the function space H(curl) and are discretized by Nédélec
More informationwith a given direct sum decomposition into even and odd pieces, and a map which is bilinear, satisfies the associative law for multiplication, and
Chapter 2 Rules of calculus. 2.1 Superalgebras. A (commutative associative) superalgebra is a vector space A = A even A odd with a given direct sum decomposition into even and odd pieces, and a map A A
More informationMultigrid Methods for Maxwell s Equations
Multigrid Methods for Maxwell s Equations Jintao Cui Institute for Mathematics and Its Applications University of Minnesota Outline Nonconforming Finite Element Methods for a Two Dimensional Curl-Curl
More informationOverlapping Schwarz Preconditioners for Spectral. Problem in H(curl)
Overlapping Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl) Technical Report TR2002-83 November 22, 2002 Department of Computer Science Courant Institute of Mathematical
More informationSheaf theory August 23, 2016
Sheaf theory August 23, 216 Chapter 1 The statement of de Rham s theorem Before doing anything fancy, let s start at the beginning. Let U R 3 be an open set. In calculus class, we learn about operations
More informationCosmology on Simplicial Complexes
Gravitation and Regge Calculus Astro Coee, Frankfurt, April 2015 Outline Gravitation and Regge Calculus 1 Gravitation and Regge Calculus Foundations of General Relativity Geometric Structure of Regge Calculus
More informationFinite Element Methods for Maxwell Equations
CHAPTER 8 Finite Element Methods for Maxwell Equations The Maxwell equations comprise four first-order partial differential equations linking the fundamental electromagnetic quantities, the electric field
More information4The Finite Element Method and Data Structures
4The Finite Element Method and Data Structures This chapter serves two purposes. The first is to point to some applied mathematics, in particular the finite element method and corresponding numerical linear
More informationOn the singular elements of a semisimple Lie algebra and the generalized Amitsur-Levitski Theorem
On the singular elements of a semisimple Lie algebra and the generalized Amitsur-Levitski Theorem Bertram Kostant, MIT Conference on Representations of Reductive Groups Salt Lake City, Utah July 10, 2013
More informationWITTEN HELLFER SJÖSTRAND THEORY. 1.DeRham Hodge Theorem. 2.WHS-Theorem. 3.Mathematics behind WHS-Theorem. 4.WHS-Theorem in the presence of.
WITTEN HELLFER SJÖSTRAND THEORY 1.DeRham Hodge Theorem 2.WHS-Theorem 3.Mathematics behind WHS-Theorem 4.WHS-Theorem in the presence of symmetry 1 2 Notations M closed smooth manifold, (V, V ) euclidean
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 information