arxiv: v1 [math.na] 5 Oct 2017
|
|
- Owen Aldous Greer
- 5 years ago
- Views:
Transcription
1 Lowest order Virtual Element approximation o magnetostatic problems L. Beirão da Veiga a,b, F. Brezzi b,, F. Dassi a, L.D. Marini c,b, A. Russo a,b a Dipartimento di Matematica e Applicazioni, Università di Milano Bicocca, Via Cozzi 53, I-20153, Milano, Italy b IMATI CNR, Via Ferrata 1, I avia, Italy c Dipartimento di Matematica, Università di avia, Via Ferrata 5, I avia, Italy Abstract arxiv: v1 [math.na] 5 Oct 2017 We give here a simpliied presentation o the lowest order Serendipity Virtual Element method, and show its use or the numerical solution o linear magneto-static problems in three dimensions. The method can be applied to very general decompositions o the computational domain (as is natural or Virtual Element Methods) and uses as unknowns the (constant) tangential component o the magnetic ield H on each edge, and the vertex values o the Lagrange multiplier p (used to enorce the solenoidality o the magnetic induction B = µh). In this respect the method can be seen as the natural generalization o the lowest order Edge Finite Element Method (the so-called irst kind Nédélec elements) to polyhedra o almost arbitrary shape, and as we show on some numerical examples it exhibits very good accuracy (or being a lowest order element) and excellent robustness with respect to distortions. Keywords: Finite Element Methods, Virtual Element Methods, Magnetostatic problems, Serendipity. AMS Subject Classiication: 65N30 1. Introduction In this paper we introduce a simpliied version o the Serendipity Virtual Element Methods (VEMs) presented in [9] and [10] and we show how they can be used or the numerical solution o linear magnetostatic problems in the so-called Kikuchi ormulation (see e.g. [21]). Serendipity VEMs are a recent variant o Virtual Element Methods that allow (as is the case o classical Serendipity Finite Elements (FEMs) on quadrilaterals and hexahedra) to eliminate a certain number o degrees o reedom (internal to aces and volumes) without compromising the order o accuracy. In the Virtual Element ramework they are particularly useul since the original ormulations o VEMs (as in [6] or [8]) oten use more degrees o reedom than their FEM counterpart (when it exists). The advantage o VEMs, when it comes to Serendipity variants, is that, contrary to FEMs, they do not use a reerence element: an inevitable sacriice, i you want to be able to deal with very general geometries. Corresponding author addresses: lourenco.beirao@unimib.it (L. Beirão da Veiga), brezzi@imati.cnr.it (F. Brezzi), ranco.dassi@unimib.it (F. Dassi), marini@imati.cnr.it (L.D. Marini), alessandro.russo@unimib.it (A. Russo) reprint submitted to Computer Methods in Applied Mechanics and Engineering November 11, 2018
2 Such a sacriice, that requires additional computations on the current element, has however the advantage o being much more robust with respect to distortions, whereas Serendipity FEMs can lose orders o accuracy already on innocent quadrilaterals that are not parallelograms (as is well known, and has been analyzed e.g. in [2], [3]). Here, as we said, we present a variant o the general theories o [9] and [10], that is specially designed or lowest order cases and comes out to be simpler, both or the theoretical presentation and the practical implementation. Then we apply it to the classical model magnetostatic problem, in a smooth-enough simply connected bounded domain in R 3 : given j H(div; ) (with divj = 0 in ), and µ = µ(x) µ 0 > 0, ind H H(curl; ) and B H(div; ) such that: (1.1) curlh = j and divb = 0, with B = µh in with the boundary conditions H n = 0 on. In particular we shall deal with the variational ormulation introduced in [21], that reads ind H H 0 (curl; ) and p H0 1 () such that: curlh curlv d + p µv d = j curlv d v H 0 (curl; ) (1.2) q µh d = 0 q H0 1 (). For many other dierent approaches to the same problem see e.g. [23], [18], [13] and the reerences therein. In our discretization, the scalar variable p (Lagrange multiplier or the condition div(µh) = 0) will be discretized using only vertex values as degrees o reedom, and the magnetic ield H will be discretized using only one degree o reedom (= constant tangential component) per edge. In its turn the current j (here a given quantity) will be discretized by its lowest order Face Virtual Element interpolant j I, individuated by its constant normal component on each ace. On tetrahedrons this would correspond to use a piecewise linear scalar or p, a lowest-order Nédélec o the irst kind or H, and a lowest order Raviart-Thomas or j: in a sense, nothing new. But already on prisms, pyramids, or hexahedra we start gaining, as we can allow more general geometries and more dramatic distortions, and there are no diiculties in using much more general polyhedrons. On polyhedrons the present approach could also be seen as being close to previous works on Mimetic Finite Dierences (the ancestor o Virtual Elements) like [16] or [22]. Here however the approach is more simple and direct, allowing a thorough analysis o convergence properties. Also the use o an explicit stabilizing term, reminiscent o Hybrid Discontinuous Galerkin methods (see e.g. [17] and the reerences therein) contributes, in our opinion, to the user-riendliness o the presentation. 2
3 A layout o the paper is as ollows. The next section will be dedicated to recall the basic notation o unctional spaces and dierential operators. Then in Section 3 we will introduce and discuss the two-dimensional VEMs (nodal and edge) that will be used on the aces o the three-dimensional decompositions. As usual, we will present irst the spaces on a single two dimensional element (local spaces). In Section 4 we will inally present our Simpliied Serendipity Spaces in three dimensions. We irst deal with a single element (polyhedron) and then discuss the spaces on a general decomposition. In Section 5 we will use these spaces to discretize the linear magneto-static problem, and briely discuss their convergence and the a-priori error analysis. Finally, in Section 6 we will present some numerical results. 2. Notation In any dimension, or an integer s 1 we will denote by s the space o polynomials o degree s. Following a common convention, 1 {0} and 0 R. Moreover, Π s,o will denote the L 2 (O)-orthogonal projection onto s (or ( s ) 2, or ( s ) 3 ). When no conusion can occur, this will be simply denoted by Π s Basic notation in 2 and 3 dimensions In two and three dimensions we will denote by x the indipendent variable. In two dimensions we will also use x = (x, y) or (more oten) x = (x 1, x 2 ) ollowing the circumstances. We will also use x := ( x 2, x 1 ) (2.1) In two dimensions, or a vector v and a scalar q we will write We observe that rotv := v 2 x v 1 y, q rot q := ( y, q x )T. (2.2) div(p 0 x) = rot(p 0 x ) = 2p 0, p 0 R. (2.3) The ollowing decompositions o polynomial vector spaces are well known and will be useul in what ollows. In two dimensions we have, or all s 0: ( s ) 2 = rot s+1 x s 1 and ( s ) 2 = grad s+1 x s 1. (2.4) In three dimensions the analogues o (2.4) (always or all s 0) are ( s ) 3 = curl(( s+1 ) 3 ) x s 1, and ( s ) 3 = grad( s+1 ) x ( s 1 ) 3. (2.5) We also note that by direct computation we have, similarly to (2.3): div(p 0 x) = 3p 0 p 0 0 and curl(p 0 x) = 2p 0 p 0 ( 0 ) 3. (2.6) 3
4 Finally, on a polyhedron we set x = x b with b the barycenter o. Analogously, or each ace, we set x = x b, with b = barycenter o. Note that clearly x d = 0, (2.7) as well as 2.2. olynomial spaces: Raviart-Thomas and Nédélec x d = 0. (2.8) We recall the deinition o the classical lowest order Raviart-Thomas local spaces in d space dimensions RT 0 = ( 0 ) d x 0, (2.9) and also the classical lowest order Nédélec irst-type local spaces in two and three space dimensions N 0 = ( 0 ) 2 x 0 or N 0 = ( 0 ) 3 x ( 0 ) 3. (2.10) In what ollows, when dealing with the aces o a polyhedron (or o a polyhedral decomposition) we shall use two-dimensional dierential operators that act on the restrictions to aces o scalar unctions that are deined on a three-dimensional domain. Similarly, or vector valued unctions we will use two-dimensional dierential operators that act on the restrictions to aces o the tangential components. In many cases, no conusion will be likely to occur; however, to stay on the sae side, we will oten use a superscript τ to denote the tangential components o a three-dimensional vector, and a subscript to indicate the twodimensional dierential operator. Hence, to ix ideas, i a ace has equation x 3 = 0 then x τ := (x 1, x 2 ) and, say, div v τ := v1 x 1 + v2 x Some Functional Spaces We recall some commonly used unctional spaces: H(div; ) = {v [L 2 ()] 3 with divv L 2 ()}, H 0 (div; ) = {ϕ H(div; ) with ϕ n = 0 on }, H(curl; ) = {v [L 2 ()] 3 with curlv [L 2 ()] 3 }, H 0 (curl; ) = {v H(curl; ) with H n = 0 on }, H 1 () = {q L 2 () with gradq (L 2 ()) 2 }, H 1 0 () = {q H 1 () with q = 0 on }. 3. Two-dimensional Simpliied Serendipity Spaces We begin with the deinition o the local spaces. Let denote a generic polyhedron, and let be a ace o such polyhedron. For the time being, we only assume that all aces are simply connected. 4
5 3.1. The local spaces on aces The ollowing spaces are a simpliied version o the local lowest-order nodal and edge Serendipity spaces on aces introduced in [7] The local nodal space on aces The present local nodal Virtual Element space on aces can be seen as an extension to polygons o the simplest space o piecewise linear unctions on triangles, as well as a simpliied version o the basic VEM nodal spaces. Indeed, in [9] (ormula (2.5) with k = 1, k = 0) the basic lowest order local nodal space on aces was introduced as { } Ṽ1 node () := q C 0 () : q e 1 (e) e, q 0 (). (3.1) It is clear that we always have 1 Ṽ node 1 (). The original degrees o reedom proposed in [6] or in [1] were The values o q at the vertices (that is : or each vertex ν, the value q(ν)); (3.2) 1 The mean value o q on (that is: or each ace, q d). (3.3) However, it was pointed out in [7] that the degrees o reedom (3.3) could be replaced by the integral q x d. (3.4) Indeed, the obvious identity 2 q d = q divx d = q x d + qx n ds (3.5) shows that given boundary values and mean value (3.3) on one can compute (3.4) and conversely, given the boundary values and (3.4) one can compute the mean value on. At this point, remembering (2.8), we easily see that q x d = 0 q 1, (3.6) so that the linear subspace o Ṽ node 1 () made o those q that veriy (3.6) will still contain 1. Such a subspace can obviously be written as { } V1 node () := q C 0 () : q e 1 (e) e, q 0 (), q x d = 0, (3.7) and will be our local scalar nodal VEM space, with the degrees o reedom (3.2) (i.e., the vertex values). Then rom the unisolvence o the degrees o reedom (3.2) and (3.3) in Ṽ node 1 () it ollows rather obviously that the degrees o reedom (3.2) are unisolvent in V node 1 (). Indeed: or a q V node 1 (), the degrees o 5
6 reedom (3.2) give the value o q along the whole. Then, ollowing (3.5), and using the inormation that q x d = 0 (that is included in the deinition (3.7) o q V1 node ()) we have 2 q d = q x n ds, (3.8) and then rom the boundary values o q we can easily compute the mean value o q on. Now we know both the degrees o reedom (3.2) and (3.3), and we are back on the track o the space Ṽ node 1 (). We point out that, i needed, out o the d.o.. in (3.2) we can compute the L 2 -projection Π 1 o q onto ( 1 ) 2. Indeed, rom the deinition o projection and (3.5) we deduce Π 1 q p 1 d := q p 1 d = q divp 1 d + q p 1 n ds, (3.9) and the last two terms are computable or every p 1 ( 1 ) The local edge space on aces The present local edge Virtual Element space on aces can be seen as an extension to polygons o the simplest space o lowest order Nédélec elements o the irst kind N 0 (see (2.10)), as well as a simpliied version o the basic VEM edge spaces in [7, 10]. Indeed, in [10] (ormula (6.1) with k = 0, k d = 0, k r = 0) the basic lowest order local edge space on aces was introduced as { } Ṽ edge 0 () := v (L 2 ()) 2 : divv 0 (), rotv 0 (), v t e 0 (e) e, (3.10) with the degrees o reedom given by the (constant) tangential components on each edge, plus the integral v x d. (3.11) The space given in (3.10) clearly contains all constant vector ields. However, or these unctions the degree o reedom (3.11) is identically zero (due to (2.8)). Hence we could consider the subspace { V edge 0 () := v (L 2 ()) 2 : divv 0 (), rotv 0 (), v t e 0 (e) e, } and v x d = 0, (3.12) that we will take, rom now on, as our edge VEM space, and have that it always contains the constant vector ields. In V edge 0 () we will thereore have the degrees o reedom: on each e, the moments v t e de. (3.13) Note that the number o degrees o reedom o V node 1 () (=number o vertexes) and o V edge 0 () (=number o edges) obviously coincide. 6 e
7 Remark 1. It is immediate to check that the space (3.12) contains, together with constant vector ields, also all vector ields o the orm v = p 0 + x p 0 with p 0 R 2 and p 0 R (that is, the lowest order Nédélec elements o the irst kind, that we recalled in (2.10)). We observe that the d.o.. (3.13) allow to compute, or each v V edge 0 (), the (constant) value o rotv on the ace by the usual Stokes theorem rotv = rotv d = v t ds, (3.14) as well as the L 2 -orthogonal projection Π 1 : V edge 0 () ( 1 ()) 2. Indeed, using (2.4) and integrating by parts we have Π 1 v p 1 d := v p 1 d = v (rot p 2 + x p 0 ) d = rot v p 2 d + v t p 2 ds + v x p 0 d = rot v p 2 d + v t p 2 ds, and all the terms in the right-hand side are computable. We close this section with a simple but important result. (3.15) roposition 3.1. It holds V node 1 () = {v V edge 0 () : rotv = 0}. (3.16) roo. We start by noting that, or any unction q V node 1 (), it holds q V edge 0 () and rot q = 0. Indeed, it is immediate to check that q satisies all the requirements in the deinition o V edge 0 () and, being a gradient, it also has vanishing rotor. Thereore V node 1 () {v V edge 0 () : rotv = 0}, that combined with dim ( V1 node () ) = dim ( V1 node () ) ( ) 1 = dim V edge 0 () 1 ( ) = dim {v V edge 0 () : rotv = 0} yields the proo. 4. Three-dimensional spaces Let be a polyhedron. For the time being, we just assume that and all o its aces are simply connected. Let N v be the number o vertices, N e the number o edges, and N the number o aces o. 7
8 4.1. The local spaces on polyhedrons For each ace we are going to use the spaces V node 1 () and V edge 0 () as deined in (3.7) and (3.12), respectively. Then we introduce their three-dimensional analogues The local nodal spaces The three dimensional local nodal space is deined as { } V1 node () := q C 0 () : q V1 node () ace, q = 0 in. (4.1) In V node 1 () the degrees o reedom are simply or each vertex ν, the nodal value q(ν). (4.2) We note that 1 () V node 1 (), since irst order polynomials clearly satisy all the conditions in (4.1) The local edge spaces In analogy with the two-dimensional case we start by recalling the three dimensional edge space deined in [10] as { Ṽ edge 0 () := v (L 2 ()) 3 : v τ V edge 0 () ace, v t e continuous at each edge e, divv = 0 in, curl(curlv) ( 0 ()) 3}, (4.3) with the degrees o reedom given by the values o the (constant) tangential components on each edge plus the integrals (curlv) (x p 0 ) d p 0 ( 0 ()) 3. (4.4) We observe, irst, that the constant vectors are contained in Ṽ edge 0 (), and or them the integral in (4.4) is always equal to zero. Hence, ollowing the path that is becoming usual here, we set { V edge 0 () := v (L 2 ()) 3 : v τ V edge 0 () ace, v t e continuous at each edge e, divv = 0 in, curl(curlv) ( 0 ()) 3, (curlv) (x p 0 ) d = 0 p 0 ( 0 ()) 3}, (4.5) and observe that all the constant vector ields are contained in V edge 0 (). Remark 2. It is easy to check that the space (4.5), together with constant vector ields, contains all vector ields o the orm v = p 0 + x q 0 with p 0 and q 0 in R 3 (that is, the lowest order Nédélec elements o the irst kind in three dimensions, as deined in (2.10)). 8
9 In V edge 0 () we have now the degrees o reedom, on each edge e the moments v t e de. (4.6) Out o the above degrees o reedom we can compute the (L 2 ()) 3 -orthogonal projection Π 0 rom V edge 0 () to ( 0 ()) 3. Indeed, by deinition o projection, (2.6), and an integration by parts we have: Π 0 v p 0 d := v p 0 d = v curl(x q 0 ) d (or q 0 = 1 2 p 0) = curlv (x q 0 ) d + (v n) (x q 0 ) ds ( ) = 0 + n (x q 0 ) v ds = ( τ n (x q 0 )) v τ d, that is computable as in (3.15). Hence, we can deine a scalar product [v, w] edge, := (Π 0 v, Π 0 w) 0, + h 2 [(v Π 0 v) t e ] [(w Π 0 w) t e ] de (4.8) e and we note that (assuming very mild mesh regularity conditions, or instance as in Section 5.1) we have e e (4.7) α (v, v) 0, [v, v] edge, α (v, v) 0, v V edge 0 () (4.9) or some constants α, α independent o h. We observe that [v, p 0 ] edge, = v p 0 d = (v, p 0 ) 0, v V edge 0 (), p 0 ( 0 ()) 3. (4.10) The local ace spaces In three dimensions we will also need a Virtual Element ace space. For it, we proceed as in the previous case. We start with the space deined in [10] { Ṽ0 ace ():= ψ (L 2 ()) 3 : ψ n 0 () ace, divψ 0 (), curlψ ( 0 ()) 3}, (4.11) where the degrees o reedom are given by the (constant) values o the normal components on the aces plus the value o the integrals: ψ (x p 0 ) d p 0 [ 0 ()] 3. (4.12) We note that the constant vector ields are inside this space, but also that the value o the integral in (4.12), or ψ constant, is always equal to zero, due to (2.7), so that we can deine { V0 ace ():= ψ (L 2 ()) 3 : ψ n 0 () ace, divψ 0 (), curlψ ( 0 ()) 3, and ψ (x p 0 ) d = 0 p 0 [ 0 ()] 3}. (4.13) 9
10 Remark 3. It is easy to check that the space (4.13), together with constant vector ields, contains all vector ields o the orm ( 0 ) 3 + x 0 (that is, the lowest order Raviart-Thomas space, as deined in (2.9)). For V0 ace () we have thereore the degrees o reedom or each ace the moments ψ n d. (4.14) Clearly, out o the degrees o reedom (4.14) we can easily compute the (constant) value o divψ divψ = 1 divψ d = 1 ψ n ds. (4.15) According to [10] we also have, now, that rom the above degrees o reedom we can compute the (L 2 ()) 3 - orthogonal projection rom V ace 0 () to ( 0 ()) 3 (and, actually, to ( 1 ()) 3 ). Indeed, using (2.5), an integration by parts, and (4.13), we have: Π 1 ψ p 1 d := ψ p 1 d = ψ ( p 2 + x p 0 ) d = divψ p 2 d + ψ n p 2 ds + ψ (x p 0 ) d = divψ p 2 d + ψ n p 2 ds, (4.16) where, using (4.15) and (4.14), all the terms in the right-hand side are computable. Hence, we can deine a scalar product [ψ, ϕ] ace, := (Π 0 ψ, Π 0 ϕ) 0, + h [(I Π 0 )ψ n][(i Π 0 )ϕ n] d, (4.17) and we again note that, assuming very mild mesh regularity conditions, or instance as in Section 5.1, α 1 (ψ, ψ) 0, [ψ, ψ] ace, α 2 (ψ, ψ) 0, ψ V ace 0 () (4.18) or positive constants α 1, α 2 independent o h. Moreover, we also have [ψ, p 0 ] ace, = ψ p 0 d = (ψ, p 0 ) 0, ψ V0 ace (), p 0 ( 0 ()) 3. (4.19) Remark 4. Using Remarks 3 and 2, by a simple dimensional count we see that the present local Nodal, Edge, and Face Virtual Element spaces coincide, whenever the element is a tetrahedron, with the classical 1, N 0 and RT 0 elements (respectively). With a minor additional eort we could see that the same is true when is a rectangular box. Note that the methods do not coincide, due to a dierent choice o the scalar products. However, the two types o scalar products are equivalent, coincide on constants, and in practice the results are not signiicantly dierent. Hence, the present setting might be considered as a sort o natural extension o the 1 -N 0 -RT 0 approach to (much) more general element geometries. 10
11 Exact sequence properties Here we present two important results. We irst have roposition 4.1. It holds V node 1 () = {v V edge 0 () : curlv = 0}. (4.20) roo. We irst point out that, with the above deinitions and using roposition 3.1, or every q V node 1 (), we easily have that the tangential gradient (applied ace by ace) will belong to V edge 0 (). Moreover, rom the deinition o V node 1 (), we immediately have that curl q = 0 and div q = q = 0. Thereore, or any q V node 1 (), it holds q {v V edge 0 () : curlv = 0} and thus V node 1 () {v V edge 0 () : curlv = 0}. Conversely, i v V edge 0 () with curlv = 0, then v = q or some q H 1 (). Since, or any ace, it holds rot (v ) = curlv n = 0, using roposition 3.1 yields that q restricted to the boundary o belongs to V node 1 (). Finally, since divv = 0 or all v V edge 0 (), we have q = 0. Thus q belongs to V node 1 (). roposition 4.2. It holds curl V edge 0 () := {ψ V ace 0 () : divψ = 0}. (4.21) roo. For every v V edge 0 () we have that ψ := curlv belongs to V ace 0 (). Indeed, on each ace we have that ψ n ( rot v ), rom (3.12) belongs to 0 () (as required in (4.13) ); moreover divψ = 0 (obviously) and curlψ ( 0 ()) 3 rom (4.5). Hence, curl V edge 0 () {ψ V ace 0 () : divψ = 0}. (4.22) We show the equality o the two sets by a dimensional count. It is immediate that the condition divψ = 0, or ψ V ace 0 (), is equivalent to and thus, recalling (4.14), ψ n d = 0 dim{ψ V ace 0 () : divψ = 0} = N 1. (4.23) By classical properties o linear operators, using (4.20) and inally the Euler ormula on polyhedrons, we get now dim(curl(v edge 0 ())) = dim(v edge 0 ()) dim{ψ V edge 0 () : curlψ = 0} and the result ollows rom (4.23) and (4.24). = N e dim( V node 1 ) = N e (N v 1) = N 1, (4.24) 11
12 4.2. The global spaces Let T h be a decomposition o the computational domain into polyhedrons. Again, or the time being, we assume just that all polyhedrons and all their aces are simply connected. More detailed assumptions will be presented in Section 5.1. Here we assume that the permeability µ is constant on each. (4.25) We deine the global spaces as ollows. with the obvious degrees o reedom { V1 node V1 node () := q H0 1 () such that q V1 node () T h }, (4.26) or each vertex ν: the nodal value q(ν). (4.27) For the global edge space we have { V edge 0 V edge 0 () := v H 0 (curl; ) such that v V edge 0 () T h }, (4.28) with the obvious degrees o reedom on each edge e: e v t e de. (4.29) Finally, or the ace space we have: { V0 ace V0 ace () := ψ H 0 (div; ) such that ψ V0 ace () T h }, (4.30) with the degrees o reedom It is important to point out that or each ace : ψ n d. (4.31) V node 1 V edge 0. (4.32) In particular, recalling the local results (4.20), it is easy to check that V node 1 {v V edge 0 such that curlv = 0}. (4.33) Similarly, we have curl V edge 0 V ace 0, (4.34) and it can be checked (recalling the local results (4.22)) that curl V edge 0 {ψ V ace 0 such that divψ = 0}. (4.35) 12
13 Remark 5. We point out that the inclusions (4.32), and (4.34) are (in a sense) also practical, and not only theoretical. By this, more speciically, we mean that: given the degrees o reedom o a q V node 1 we can compute the corresponding degrees o reedom o q in V edge 0 ; and given the degrees o reedom o a v V edge 0 we can compute the corresponding degrees o reedom o curl v in V ace 0. From (4.8) we can also deine a global scalar product: We note that, recalling (4.9), [v, w] edge := T h [v, w] edge,. (4.36) α (v, v) 0, [v, v] edge α (v, v) 0, v V edge 0. (4.37) It is also important to point out that, using (4.10) we have [v, p] edge = (v, p) 0, := v p d v V edge 0, p piecewise in ( 0 ) 3. (4.38) From (4.17) we can also deine a scalar product in V ace 0 in the obvious way [ψ, ϕ] ace := T h [ψ, ϕ] ace, (4.39) and we note that α 1 (ψ, ψ) 0, [ψ, ψ] ace α 2 (ψ, ψ) 0, ψ V ace 0. (4.40) It is also important to point out that, using (4.19) we have [ψ, p] ace = (ψ, p) 0, := ψ p d ψ V0 ace, p piecewise in ( 0 ) 3. (4.41) 5. Discretization o the Magneto-static roblem We are now ready to present the discretization o our Magneto-static roblem (1.2) that we recall here: given j H(div; ) (with divj = 0 in ), and µ = µ(x) µ 0 > 0, ind H H 0 (curl; ) and p H0 1 () such that: curlh curlv d + p µv d = j curlv d v H 0 (curl; ) q µh d = 0 q H0 1 (). It is easy to check, by the usual theory o mixed methods, that (5.1) has a unique solution (H, p). Then we check that H and µh give the solution o (1.1) and p = 0. Checking that p = 0 is immediate, just taking v = p in the irst equation. Once we know that p = 0, the irst equation gives curlh = j, and then the second equation gives divµh = (5.1)
14 Remark 6. We observe that an alternative variational ormulation could be given j H(div; ) (with divj = 0 in ), and µ = µ(x) µ 0 > 0, ind H H 0 (curl; ) and p H0 1 () such that: curlh curlv d + p µv d = j curlv d v H 0 (curl; ) q µh d p q d = 0 q H0 1 (), or other possible variants mimicking, one way or another, the Hodge-Laplacian approach (see [4]). observe that the discretization that we are going to introduce or (5.1) will apply to (5.2) as well. We irst construct the interpolant j I V0 ace o j that matches the degrees o reedom (4.31): or each ace : (j j I ) n d = 0. (5.3) From the d.o.. (5.3) and an integration by parts it ollows that div(j j I ) d = 0 T h. (5.4) Moreover, they also imply that j I H 0 (div; ) and that divj I = 0 in. Hence, according to (4.35), we have that j I will be the curl o some w V edge 0 : (5.2) w V edge 0 such that curlw = j I. (5.5) We Then we can introduce the discretization o (1.2): ind H h V edge 0 and p h V1 node such that: [curlh h, curlv] ace + [ p h, µv] edge = [j I, curlv] ace [ q, µh h ] edge = 0 q V node 1. v V edge 0 (5.6) We point out that both curlh h and curlv (as well as j I ) are ace Virtual Elements in V ace 0 () in each polyhedron, so that (taking also into account Remark 5) their ace scalar products are computable as in (4.39). Similarly, rom the degrees o reedom o a q V node 1 we can compute the degrees o reedom o q, as an element o V edge 0, so that the two edge-scalar products that appear in (5.6) are computable as in (4.36). roposition 5.1. roblem (5.6) has a unique solution (H h, p h ), and p h 0. roo. Taking v = p h (as we did or the continuous problem (5.1)) in the irst equation, and using (4.37) we easily obtain p h 0 or (5.6) as well. To prove uniqueness o H h, set j I = 0, and let H h be the solution o the homogeneous problem. From the irst equation we deduce that curl H h = 0. Hence, rom (4.33) we have H h = q h or some q h V node 1. The second equation and (4.37) give then H h = 0. 14
15 Once we know that p h = 0, the irst equation o (5.6) reads [curlh h, curlv] ace = [j I, curlv] ace v V edge 0, (5.7) that in view o (5.5) becomes [curlh h curlw, curlv] ace = 0 v V edge 0. (5.8) Using v = H h w and (4.40), this easily implies curlh h = curlw = j I. (5.9) 5.1. Error estimates For the theoretical derivations we consider the ollowing mesh assumptions, that are quite standard in the VEM literature. We assume the existence o a positive constant γ such that any polyhedron (o diameter h ) satisies the ollowing conditions: 1. is star shaped with respect to a ball o radius bigger than γh ; 2. any ace is star shaped with respect to a ball o radius bigger than γh, and every edge o has length bigger than γh. We note that condition 1 (and 2) implies that (and any ace o ) is simply connected. At the theoretical level, some o the above conditions could be avoided by using more technical arguments. At the practical level, as shown by the numerical tests o the Section 6, condition 2 is negligible since the method seems essentially impervious to degeneration o aces and edges. On the contrary, although the scheme is quite robust to distortion o the elements, condition 1 is more relevant since extremely anisotropic element shapes can lead to poor results. We also recall that we assumed µ to be piecewise constant (see (4.25)). Let us bound the error H H h. We start by deining the interpolant H I V edge 0 o H, deined through the degrees o reedom (4.29): on each edge e: (H H I ) t e ds = 0. (5.10) e roposition 5.2. With the choices (5.3) and (5.10) we have curlh I = j I. (5.11) roo. From (4.35) we know that curlh I V ace 0. To prove (5.11) we should just show that the ace degrees o reedom (4.31) o the dierence curlh I j I are zero, that is: : (curlh I j I ) n d = 0. (5.12) 15
16 Since j = curlh, rom the interpolation ormulas (5.3) we see that : (curlh j I ) n d = 0, so that we can replace j I with curlh in (5.12), that becomes : curl(h I H) n d = 0. (5.13) Observing that (5.10) implies that rot (H H I ) d = 0, and recalling that on each the normal component o curl(h I H) is equal to the rot o the tangential components (H I H), we deduce curl(h I H) n d rot (H I H) d = 0. Hence, (5.13) is satisied, and the proo is concluded. From (5.9) and (5.11) it ollows then curl(h I H h ) = 0 (5.14) and thereore, rom (4.33), H I H h = q h or some q h V node 1. (5.15) Now we deine an alternative (L 2 ()) 3 inner product and norm that take into account the (piecewise constant) value o the permeability µ. We set, or v in (L 2 ()) 3, v 2 0, := µ v 2 d. (5.16) When v = H in (5.16) we get that H 2 0, = B H d, showing the connection between the new norm and the energy. We now note that, using (4.37), we have α H I H h 2 0, [H I H h, µ(h I H h )] edge, (5.17) and also (rom (4.38) and (4.25)) [p 0, µv] edge =(p 0, µv) 0, v V edge 0 and p 0 piecewise constant vector. (5.18) 16
17 Then, starting with (5.17) we have: α H I H h 2 0, [H I H h, µ(h I H h )] edge = (using (5.15)) [H I H h, µ q h] edge = (using the second o (5.6)) [H I, µ q h] edge = (adding and subtracting Π 0 H) [H I Π 0 H, µ q h] edge + [Π 0 H, µ q h] edge = (using (5.18)) [H I Π 0 H, µ q h] edge + (Π 0 H, µ q h) 0, = (adding and subtracting H) [H I Π 0 H, µ q h] edge + (Π 0 H H, µ q h) 0, + (H, µ q h) 0, = (rom the second o (1.2)) [H I Π 0 H, µ q h] edge + (Π 0 H H, µ q h) 0, ) (using Cauchy-Schwarz and (4.37)) (α H I Π 0 H 0, + Π 0 H H 0, qh 0, ) (using again (5.15)) (α H I Π 0 H 0, + Π 0 H H 0, H I H h 0, that implies immediately (since α 1) H I H h 0, α α ( H I Π 0 H 0, + Π 0 H H 0, ). (5.19) We can summarize the result in the ollowing theorem. Theorem 5.3. roblem (5.6) has a unique solution and the ollowing estimate holds: ) H H h 0, C ( H H I 0, + H Π 0 H 0,, (5.20) with C a constant independent o the mesh size. Moreover, thanks to (5.9) we also have curl(h H h ) 0, = j j I 0,. (5.21) The above result can be combined with standard polynomial approximation estimates on star shaped polyhedra in order to estimate the terms involving the L 2 projection on polynomials. Moreover, approximation estimates or the VEM interpolants H I and j I can be derived (under the mesh assumptions at the beginning o this section) by an extension o the arguments in [12], [11], [15], [24]. We thereore obtain, provided that H and j are suiciently regular, H H h 0, Ch curl(h H h ) 0, Ch. 6. Numerical results In this section we provide some numerical results. We irst provide a test on a polyhedral domain with Dirichlet boundary data, then a test on a cylindrical domain with jumping coeicients and Neumann boundary data, and inally a benchmark example with a more complex geometry. 17
18 6.1. Test case 1: h-analysis with homogeneous Dirichlet boundary conditions In this subsection we set µ = 1 and take as exact solution o (1.1) sin(πy) sin(πz) 1 H(x, y, z) := sin(πz) sin(πx) π sin(πx) sin(πy) the vector ield. We consider as domain Ω the truncated octahedron [28] and three dierent discretizations, see Fig. 1: Structured, a mesh composed by structured cubes inside Ω and arbitrarily shaped elements close to the boundary; CVT, a Centroidal Voronoi Tessellation o Ω obtained via a standard Lloyd algorithm [19]; Random, a mesh obtained by the constrained Voronoi Tessellation o points randomly put inside Ω. The irst type o mesh can be generated by using a Voronoi generation algorithm with seeds disposed in a regular ormation, see [5]. This approach allows to approximate complex geometries and still inherits many among the computational advantages o regular cubical meshes. The last type o meshes is instead interesting in order to check the robustness o the method, since it exhibits edges and aces (see the details in Fig. 1) that can be very small with respect to the size o the parent element. To get such discretizations, we use the c++ library voro++ [25]. Structured CVT Figure 1: Three dierent discretizations o the truncated octahedron. We associate to each mesh a mesh-size h deined as h := N 1 X h, N i=1 18 Random
19 where N is the number o polyhedrons in the mesh and h is the diameter o the polyhedron. We compute the L 2 -error on H as H Π 0 H h 0, H 0,, where Π 0 H h is the piecewise constant projection o H h deined in (4.7). Fig. 2 (let) shows the convergence rate or each type o mesh. The slopes are coherent with the theory, see Equation (5.20). Moreover, rom this graph we can better appreciate the robustness o the VEM with respect to element distortions. Indeed, the convergence lines or the three meshes are very close to each other. l norm step Structured CVT Random e e e e e e e e e e e e-07 Figure 2: Test case 1: L 2 error or H (let), and l norm o do(p h ) (right). The results or the Lagrange multiplier p h are shown in Fig. 2 (right). As shown in Section 5, the exact solution p h o the discrete problem is identically zero. Clearly, the roundo errors generate, out o the computer, a p h that is almost identically zero. In some sense we could then take the value o the computed p h as a measure (or, better, a rough indicator) o the conditioning o the inal linear system. In particular we see that the Random meshes generate a worse conditioning Test case 2: discontinuous µ and homogeneous Neumann conditions In this subsection we consider an example taken rom [14]. The geometry consists in a cylindrical domain with two concentric connected components, S 1 and S 2, separated by a magnetic material, M (see a crosssection in Fig. 3). A current o the same intensity I =70 000A but opposite direction passes along S 1 and S 2. We assign the ollowing current intensity (with e z := (0, 0, 1) t ) I πa 2 e z i x S 1 j(x, y, z) := 0 i x M (6.1) I π(c 2 b 2 ) e z i x S 2. 19
20 a = 0.5 m b = 1.0 m c = 1.25 m Figure 3: Test case 2: meridian cross section o the computational domain. We apply homogeneous Neumann boundary conditions and we set µ S1 = µ S2 = 1.0, µ M = Then the exact solution is given by I 2πa 2 e θ i x S 1 I H(x, y, z) := 2πr 2 e θ [ I 2π(c 2 b 2 ) + 1 r 2 i x M ( I 2π + Ib 2 )] 2π(c 2 b 2 e θ i x S 2 ) (6.2) where e θ := ( y, x, 0) t, r = x 2 + y 2 and a = 0.5, b = 1, c = 1.25 (see Fig. 3). To deal with Neumann boundary conditions, problem (5.1) must be modiied, looking or H H(curl; ) and p H 1 (). The magnetic ield H will still be unique, while p will be determined only up to a constant. In the code the constant is ixed by requiring the average o the vertex values o p h to be zero. We build two meshes o the cylinder by extruding two planar two-dimensional meshes: one made o polygons [27] and one made o triangles [26]. We reer to the ormer as voro-extrusion and to the latter as tria-extrusion (or an example o both see Fig. 4). Then we make two sets o meshes with decreasing meshsize. In Fig. 5 we depict the convergence curves on the let, and on the right we provide the l norm o the vectors do(p h ). Both quantities behave as expected: indeed, we get a convergence rate equal to 1, or the L 2 error on the vector ield H h, while the values o l norm vanish up to machine algebra errors Test case 3: a cylindrical electromagnet In this subsection we consider a typical benchmark problem, see e.g. [13, 14, 20]. The geometry consists in a erromagnetic cylindrical core, C, surrounded by a toroidal coil with a rectangular cross section, T, with air, A, around these two structures. In Fig. 6 we show a meridian cross section o the domain where we speciy the dimensions o the cylindrical core, the toroidal coil and the bounding box o the domain. 20
21 Figure 4: Test case 2: two polyhedral decompositions o the cylinder: voro-extrusion (let), and tria-extrusion (right) l norm step voro-extrusion tria-extrusion e e e e e e e e Figure 5: Test case 2: L 2 error or H (let), and l norm o do(p h ) (right). 0.2 m 0.05 m 1. m 0.2 m 0.5 m 2. m 0.5 m 0.3 m 1. m Figure 6: Test case 3: meridian cross section o the computational domain. 21
22 One looks or the magnetic lux generated by a constant current passing along the toroidal coil. More speciically, we assume a constant current o 1 Ampère, i.e., 1 j(x, y, z) := A e θ e θ i x T 0 otherwise where e θ = ( y, x, 0) t, e θ denotes the norm o e θ and A is the area o the cross section o the coil. The relative magnetic permeability is taken as: µ T = 1.0 or the coil, µ C = in the erromagnetic core, and µ A = 1.0 or the air. Moreover, we suppose that the artiicial boundary is suiciently ar rom both the erromagnetic cylinder and the toroidal core so that we can apply the Neumann boundary conditions H n = 0. In Fig. 7 we show one o the meshes used in this example. Here too the meshes are constructed by extruding two-dimensional ones. Figs. 8 and 9 show a qualitative comparison between the solution provided Figure 7: Test case 3: On the right a clip o the 3d mesh obtained extruding the 2d mesh on the let. in [13] (let), and that obtained by the present method (right). More speciically, in Fig. 8, we show the modulus o the magnetic lux density, i.e. the modulus o B = µh, inside the cylindrical core. Then, in Fig. 9 we show B along some cross sections o the cylindrical core and the toroidal coil. In both cases the results o two methods show a comparable behavior. Finally, to have a more quantitative validation o this example, we compute the so-called magnetic energy W := B H d = µ H 2 d, D D or each sub-domain C, T and A (see also (5.16)). We consider a sequence o three nested meshes to veriy 22
23 Figure 8: Test case 3: B in a section o C: rom [13] (let) and that rom the present approach (right). Figure 9: Test case 3: B in a section o C and T : rom [13] (let), and rom the present approach (right). the convergence rate o these energies. We reer to these meshes as mesh1, mesh2 and mesh3, the irst mesh being the coarsest. 23
24 Since we do not have the exact solution o this problem, in order to assess the accuracy o our results or the magnetic energy W we proceed as in [14], i.e. we consider as exact values o the magnetic energies the values obtained by the sotware FLUX2D on a very ine mesh. The numerical solution o FLUX2D ollows a scalar potential ormulation and exploits the symmetry o the domain via a two-dimensional cylindric coordinate system. The data are collected in Table 1 and in parenthesis we show the relative error. The method behaves as expected. W in A C T FLUX2D 9.09e e e-08 mesh1 VEM 9.70e-07 (6.7%) 7.54e-10 (59.4%) 3.29e-08 (8.8%) mesh2 VEM 9.22e-07 (1.4%) 5.53e-10 (16.9%) 3.81e-08 (5.5%) mesh3 VEM 9.11e-07 (0.2%) 4.98e-10 (5.2%) 3.75e-08 (3.8%) Table 1: Test case 3: Behavior o W, in the three regions, computed with the present method Acknowledgments The irst and third authors were partially supported by the European Research Council through the H2020 Consolidator Grant (grant no ) CAVE Challenges and Advancements in Virtual Elements. This support is grateully acknowledged. Reerences [1] B. Ahmad, A. Alsaedi, F. Brezzi, L. D. Marini, and A. Russo, Equivalent projectors or virtual element methods, Comput. Math. Appl. 66 (2013), no. 3, [2] D.N. Arnold, D. Boi, and R.S. Falk, Approximation by quadrilateral inite elements, Math. Comp. 71 (2002), no. 239, [3], Quadrilateral H(div) inite elements, SIAM J. Numer. Anal. 42 (2005), no. 6, [4] D.N. Arnold, R.S. Falk, and R. Winther, Finite element exterior calculus, homological techniques, and applications, Acta Numer. 15 (2006), [5] L. Beirão da Veiga, F. Dassi, and A. Russo, High-order virtual element method on polyhedral meshes, Computers & Mathematics with Applications (2017). [6] L. Beirão da Veiga, F. Brezzi, A. Cangiani, G. Manzini, L. D. Marini, and A. Russo, Basic principles o virtual element methods, Math. Models Methods Appl. Sci. 23 (2013), no. 1, [7] L. Beirão da Veiga, F. Brezzi, F. Dassi, L. D. Marini, and A. Russo, Virtual element approximation o 2d magnetostatic problems, Comput. Methods Appl. Mech. Engrg. (2017), [8] L. Beirão da Veiga, F. Brezzi, L. D. Marini, and A. Russo, H(div) and H(curl)-conorming VEM, Numer. Math. 133 (2016), no. 2,
25 [9], Serendipity nodal VEM spaces, Comp. Fluids 141 (2016), [10], Serendipity ace and edge VEM spaces, Rend. Lincei Mat. Appl. 28 (2017), no. 1, [11] L. Beirão da Veiga, C. Lovadina, and A. Russo, Stability analysis or the virtual element method, Math. Models Methods Appl. Sci. (2017), /S X. [12] L. Beirão da Veiga, D. Mora, G. Rivera, and R. Rodríguez, A virtual element method or the acoustic vibration problem, Numerische Mathematik 136 (2017), [13] A. Bermúdez, D. Gómez, and. Salgado, Mathematical models and numerical simulation in electromagnetism, Unitext, vol. 74, Springer, [14] A. Bermúdez, R. Rodríguez, and. Salgado, A inite element method or the magnetostatic problem in terms o scalar potentials, SIAM Journal on Numerical Analysis 46 (2008), no. 3, [15] S. Brenner, Q. Guan, and L. Sung, Some estimates or virtual element methods, Computational Methods in Applied Mathematics, doi: /cmam , [16] F. Brezzi and A. Bua, Innovative mimetic discretizations or electromagnetic problems, J. Comput. Appl. Math 234 (2010), [17] B. Cockburn, D. Di ietro, and A. Alexandre Ern, Bridging the hybrid high-order and hybridizable discontinuous galerkin methods, ESAIM Math. Model. Numer. Anal. 50 (2016), [18] L. Demkowicz, J. Kurtz, D. ardo, M. aszenski, W. Rachowicz, and A. Zdunek, Computing with hp-adaptive inite elements. vol. 2. rontiers: Three dimensional elliptic and maxwell problems with applications, Applied Mathematics and Nonlinear Science, Chapman & Hall/CRC, Boca Raton, [19] Qiang Du, V. Faber, and M. Gunzburger, Centroidal voronoi tessellations: Applications and algorithms, SIAM Rev. 41 (1999), no. 4, [20] M. Fontana, Adaptive inite element computation o 3D magnetostatic problems in potential ormulation, Tech. Report reprint , University o Chalmers, Göteborg, [21] H. Kanayama, R. Motoyama, K. Endo, and F. Kikuchi, Three dimensional magnetostatic analysis using nedelec s elements, IEEE Transactions on Magnetics 26 (1990), [22] K. Lipnikov, G. Manzini, F. Brezzi, and A. Bua, The mimetic inite dierence method or the 3d magnetostatic ield problems on polyhedral meshes, J. Comput. hys. 230 (2011), [23]. Monk, Finite element methods or Maxwell s equations, Numerical Mathematics and Scientiic Computation, Oxord University ress, New York, [24] D. Mora, G. Rivera, and R. Rodríguez, A virtual element method or the Steklov eigenvalue problem, Math. Models Methods Appl. Sci. 25 (2015), no. 8, [25] C.H. Rycrot, Voro++: A three-dimensional voronoi cell library in c++, Chaos 19 (2009), no. 4, [26] J.R. Shewchuk, Triangle: Engineering a 2d quality mesh generator and delaunay triangulator, Applied computational geometry towards geometric engineering, Springer, 1996, pp [27] C. Talischi, G.H. aulino, A. ereira, and I.FM Menezes, olymesher: a general-purpose mesh generator or polygonal elements written in matlab, Structural and Multidisciplinary Optimization 45 (2012), no. 3, [28] R. Williams, The geometrical oundation o natural structure: A source book o design, Dover ublications,
H(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 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 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 informationarxiv: v1 [math.na] 3 Jun 2016
SRNDIPITY FAC AND DG VM SPACS L. BIRÃO DA VIGA[1], F. BRZZI [2], L.D. MARINI [3], A. RUSSO [1] arxiv:1606.01048v1 [math.na] 3 Jun 2016 Abstract. We extend the basic idea of Serendipity Virtual lements
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 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 informationVirtual Element Methods for general second order elliptic problems
Virtual Element Methods for general second order elliptic problems Alessandro Russo Department of Mathematics and Applications University of Milano-Bicocca, Milan, Italy and IMATI-CNR, Pavia, Italy Workshop
More informationThe Virtual Element Method: an introduction with focus on fluid flows
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.
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 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 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 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 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 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 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 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 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 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 informationWEAK GALERKIN FINITE ELEMENT METHODS ON POLYTOPAL MESHES
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 1, Pages 31 53 c 2015 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHODS ON POLYOPAL MESHES LIN
More 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 informationA Simple Explanation of the Sobolev Gradient Method
A Simple Explanation o the Sobolev Gradient Method R. J. Renka July 3, 2006 Abstract We have observed that the term Sobolev gradient is used more oten than it is understood. Also, the term is oten used
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 informationNonstandard finite element de Rham complexes on cubical meshes arxiv: v1 [math.na] 12 Apr 2018
Nonstandard inite element de Rham complexes on cubical meshes arxiv:184.439v1 [math.na] 12 Apr 218 Andrew Gillette aibo Hu Shuo Zhang Abstract We propose two general operations on inite element dierential
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 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 informationGeneral solution of the Inhomogeneous Div-Curl system and Consequences
General solution o the Inhomogeneous Div-Curl system and Consequences Briceyda Berenice Delgado López Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional March 07 2017 CINVESTAV
More information2.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics. References - geostatistics. References geostatistics (cntd.
.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics Spline interpolation was originally developed or image processing. In GIS, it is mainly used in visualization o spatial
More informationANDREA TOSELLI. Abstract. Two-level overlapping Schwarz methods are considered for nite element problems
OVERLAPPING SCHWARZ METHODS FOR MAXWELL'S EQUATIONS IN THREE DIMENSIONS ANDREA TOSELLI Abstract. Two-level overlapping Schwarz methods are considered for nite element problems of 3D Maxwell's equations.
More informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
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 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 informationCoupling of eddy-current and circuit problems
Coupling of eddy-current and circuit problems Ana Alonso Rodríguez*, ALBERTO VALLI*, Rafael Vázquez Hernández** * Department of Mathematics, University of Trento ** Department of Applied Mathematics, University
More informationLocal flux mimetic finite difference methods
Local flux mimetic finite difference methods Konstantin Lipnikov Mikhail Shashkov Ivan Yotov November 4, 2005 Abstract We develop a local flux mimetic finite difference method for second order elliptic
More informationNumerische Mathematik
Numer. Math. (2012) 122:61 99 DOI 10.1007/s00211-012-0456-x Numerische Mathematik C 0 elements for generalized indefinite Maxwell equations Huoyuan Duan Ping Lin Roger C. E. Tan Received: 31 July 2010
More informationFinite element approximation on quadrilateral meshes
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2001; 17:805 812 (DOI: 10.1002/cnm.450) Finite element approximation on quadrilateral meshes Douglas N. Arnold 1;, Daniele
More informationA priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces
A priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces A. Bespalov S. Nicaise Abstract The Galerkin boundary element discretisations of the
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 informationSTABILIZED HP-DGFEM FOR INCOMPRESSIBLE FLOW. Math. Models Methods Appl. Sci., Vol. 13, pp , 2003
Mathematical Models and Methods in Applied Sciences c World Scientiic Publishing Company STABILIZED HP-DGFEM FOR INCOMPRESSIBLE FLOW D. SCHÖTZAU1 C. SCHWAB 2 A. TOSELLI 2 1 Department o Mathematics, University
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 informationarxiv: v2 [math.na] 23 Apr 2016
Improved ZZ A Posteriori Error Estimators for Diffusion Problems: Conforming Linear Elements arxiv:508.009v2 [math.na] 23 Apr 206 Zhiqiang Cai Cuiyu He Shun Zhang May 2, 208 Abstract. In [8], we introduced
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 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 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 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 informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
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 informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationScattered Data Approximation of Noisy Data via Iterated Moving Least Squares
Scattered Data Approximation o Noisy Data via Iterated Moving Least Squares Gregory E. Fasshauer and Jack G. Zhang Abstract. In this paper we ocus on two methods or multivariate approximation problems
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 informationStrong Lyapunov Functions for Systems Satisfying the Conditions of La Salle
06 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 49, NO. 6, JUNE 004 Strong Lyapunov Functions or Systems Satisying the Conditions o La Salle Frédéric Mazenc and Dragan Ne sić Abstract We present a construction
More informationGlobal Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance
Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam) Abstract. A reormulation o the planetary
More informationFluctuationlessness Theorem and its Application to Boundary Value Problems of ODEs
Fluctuationlessness Theorem and its Application to Boundary Value Problems o ODEs NEJLA ALTAY İstanbul Technical University Inormatics Institute Maslak, 34469, İstanbul TÜRKİYE TURKEY) nejla@be.itu.edu.tr
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 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 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 information1 Relative degree and local normal forms
THE ZERO DYNAMICS OF A NONLINEAR SYSTEM 1 Relative degree and local normal orms The purpose o this Section is to show how single-input single-output nonlinear systems can be locally given, by means o a
More informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
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 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 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 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 information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationThe Relationship between Discrete Calculus Methods and other Mimetic Approaches
The Relationship between Discrete Calculus Methods and other Mimetic Approaches Blair Perot Theoretical and Computational Fluid Dynamics Laboratory Woudschoten Conerence Oct. 2011 October 7, 2011 Background
More informationA note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations
A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations Bernardo Cockburn Guido anschat Dominik Schötzau June 1, 2007 Journal of Scientific Computing, Vol. 31, 2007, pp.
More informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
More informationarxiv: v3 [math.na] 8 Sep 2015
A Recovery-Based A Posteriori Error Estimator for H(curl) Interface Problems arxiv:504.00898v3 [math.na] 8 Sep 205 Zhiqiang Cai Shuhao Cao Abstract This paper introduces a new recovery-based a posteriori
More informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More 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 informationarxiv: v1 [math.na] 15 Nov 2017
Noname manuscript No. (will be inserted by the editor) An HDG method with orthogonal projections in facet integrals Issei Oikawa arxiv:1711.05396v1 [math.na] 15 Nov 2017 Received: date / Accepted: date
More informationReview of Prerequisite Skills for Unit # 2 (Derivatives) U2L2: Sec.2.1 The Derivative Function
UL1: Review o Prerequisite Skills or Unit # (Derivatives) Working with the properties o exponents Simpliying radical expressions Finding the slopes o parallel and perpendicular lines Simpliying rational
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
More 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 informationHEXAHEDRAL H(DIV) AND H(CURL) FINITE ELEMENTS
HEXAHEDRAL H(DIV) AND H(CURL) FINITE ELEMENTS RICHARD S. FAL, PAOLO GATTO, AND PETER MON Abstract. We study the approximation properties of some finite element subspaces of H(div; Ω) and H(curl ; Ω) defined
More informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
More informationEXISTENCE OF SOLUTIONS TO SYSTEMS OF EQUATIONS MODELLING COMPRESSIBLE FLUID FLOW
Electronic Journal o Dierential Equations, Vol. 15 (15, No. 16, pp. 1 8. ISSN: 17-6691. URL: http://ejde.math.txstate.e or http://ejde.math.unt.e tp ejde.math.txstate.e EXISTENCE OF SOLUTIONS TO SYSTEMS
More informationINNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD
INNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD Abstract. Finite element methods for the Reissner Mindlin plate theory are discussed. Methods in which both the tranverse displacement and
More informationDiscontinuous Galerkin Methods: Theory, Computation and Applications
Discontinuous Galerkin Methods: Theory, Computation and Applications Paola. Antonietti MOX, Dipartimento di Matematica Politecnico di Milano.MO. X MODELLISTICA E CALCOLO SCIENTIICO. MODELING AND SCIENTIIC
More informationAn introduction to Isogeometric Analysis
An introduction to Isogeometric Analysis A. Buffa 1, G. Sangalli 1,2, R. Vázquez 1 1 Istituto di Matematica Applicata e Tecnologie Informatiche - Pavia Consiglio Nazionale delle Ricerche 2 Dipartimento
More 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 informationFind (u,p;λ), with u 0 and λ R, such that u + p = λu in Ω, (2.1) div u = 0 in Ω, u = 0 on Γ.
A POSTERIORI ESTIMATES FOR THE STOKES EIGENVALUE PROBLEM CARLO LOVADINA, MIKKO LYLY, AND ROLF STENBERG Abstract. We consider the Stokes eigenvalue problem. For the eigenvalues we derive both upper and
More informationIterative Methods for Stokes/Darcy Coupling
Iterative Methods or Stokes/Darcy Coupling Marco Discacciati Ecole Polytechnique Fédérale de Lausanne Institut d Analyse et Calcul Scientiique Chair o Modelling and Scientiic Computing marco.discacciati@epl.ch
More informationThe Clifford algebra and the Chevalley map - a computational approach (detailed version 1 ) Darij Grinberg Version 0.6 (3 June 2016). Not proofread!
The Cliord algebra and the Chevalley map - a computational approach detailed version 1 Darij Grinberg Version 0.6 3 June 2016. Not prooread! 1. Introduction: the Cliord algebra The theory o the Cliord
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 informationBANDELET IMAGE APPROXIMATION AND COMPRESSION
BANDELET IMAGE APPOXIMATION AND COMPESSION E. LE PENNEC AND S. MALLAT Abstract. Finding eicient geometric representations o images is a central issue to improve image compression and noise removal algorithms.
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 informationNumerical Solution of Ordinary Differential Equations in Fluctuationlessness Theorem Perspective
Numerical Solution o Ordinary Dierential Equations in Fluctuationlessness Theorem Perspective NEJLA ALTAY Bahçeşehir University Faculty o Arts and Sciences Beşiktaş, İstanbul TÜRKİYE TURKEY METİN DEMİRALP
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 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 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 informationMathematical Notation Math Calculus & Analytic Geometry III
Name : Mathematical Notation Math 221 - alculus & Analytic Geometry III Use Word or WordPerect to recreate the ollowing documents. Each article is worth 10 points and can e printed and given to the instructor
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
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 informationStandard Finite Elements and Weighted Regularization
Standard Finite Elements and Weighted Regularization A Rehabilitation Martin COSTABEL & Monique DAUGE Institut de Recherche MAthématique de Rennes http://www.maths.univ-rennes1.fr/~dauge Slides of the
More informationOn a Closed Formula for the Derivatives of e f(x) and Related Financial Applications
International Mathematical Forum, 4, 9, no. 9, 41-47 On a Closed Formula or the Derivatives o e x) and Related Financial Applications Konstantinos Draais 1 UCD CASL, University College Dublin, Ireland
More informationMixed Discontinuous Galerkin Methods for Darcy Flow
Journal of Scientific Computing, Volumes 22 and 23, June 2005 ( 2005) DOI: 10.1007/s10915-004-4150-8 Mixed Discontinuous Galerkin Methods for Darcy Flow F. Brezzi, 1,2 T. J. R. Hughes, 3 L. D. Marini,
More informationWEAK GALERKIN FINITE ELEMENT METHOD FOR SECOND ORDER PARABOLIC EQUATIONS
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 13, Number 4, Pages 525 544 c 216 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND ORDER PARABOLIC
More informationTelescoping Decomposition Method for Solving First Order Nonlinear Differential Equations
Telescoping Decomposition Method or Solving First Order Nonlinear Dierential Equations 1 Mohammed Al-Reai 2 Maysem Abu-Dalu 3 Ahmed Al-Rawashdeh Abstract The Telescoping Decomposition Method TDM is a new
More informationNON-AUTONOMOUS INHOMOGENEOUS BOUNDARY CAUCHY PROBLEMS AND RETARDED EQUATIONS. M. Filali and M. Moussi
Electronic Journal: Southwest Journal o Pure and Applied Mathematics Internet: http://rattler.cameron.edu/swjpam.html ISSN 83-464 Issue 2, December, 23, pp. 26 35. Submitted: December 24, 22. Published:
More informationMathematical Notation Math Calculus & Analytic Geometry III
Mathematical Notation Math 221 - alculus & Analytic Geometry III Use Word or WordPerect to recreate the ollowing documents. Each article is worth 10 points and should be emailed to the instructor at james@richland.edu.
More information