26.1 Definition and examples
|
|
- Gerald Atkinson
- 6 years ago
- Views:
Transcription
1 Part V, Chapter 26 Quadratures Evaluating integrals over cells and faces are frequent tasks when implementing a solution method based on finite elements. Such integrals are often approximately evaluated by so-called quadratures. In this chapter, we define the notion of quadratures and review one- and multi-dimensional quadratures frequently used in finite element codes. Then, as an example of application, we analyze the impact of the quadrature error in the finite element approximation of a model second-order elliptic problem (see Chapter 25). 26. Definition and examples Let φ be a smooth function and let T h be a mesh of a domain D in R d that covers D exactly. Suppose we want to evaluate the integral Dφ(x)dx. Since φ(x)dx = φ(x) dx, D K T h this problem reduces to evaluating the integral over each mesh cell. An effective way of doing this approximately is by means of quadratures. Definition 26. (Quadrature nodes and weights). Let K be a compact, connected, Lipschitz subset of R d with non-empty interior. Let l q be an integer. A quadrature in K with l q nodes is specified through a set of l q points {ξ l } l {:lq} in K, called Gaussnodes and a set of l q real numbers {ω l } l {:lq}, called quadrature weights. The quadrature then consists of the approximation φ(x)dx ω l φ(ξ l ). (26.) K K l {:l q} The largest integer k such that (26.) is an equality for any polynomial in P k,d is called the quadrature order and is denoted k q.
2 374 Chapter 26. Quadratures Given a quadrature on the reference element K and a mesh T h, a quadrature in any cell K T h can be generated by using the geometric map T K : K K. Let J K the Jacobian matrix of T K. Proposition 26.2 (Quadrature generation). Consider a quadrature in K with nodes { ξ l } l {:lq} and weights { ω l } l {:lq}. Then, setting ξ lk := T K ( ξ l ) and ω lk := ω l det ( J K ( ξ l ) ), (26.2) for all l {:l q }, generates a quadrature on K. If the quadrature on K is of order k q and the geometric map T K is affine, then the quadrature on K is also of order k q. Proof. Since T K is a C -diffeomorphism, the change of variables x = T K ( x) yields K φ(x)dx = φ ( T K K ( x) ) det ( J K ( x) ) d x, and we can apply the quadrature over K to the right-hand side. The statement on the quadrature order is immediate to verify; indeed, whenever T K is affine, J K is constant and φ T K is in P k,d if φ P k,d. Remark 26.3 (Surface quadrature). When generating a surface quadrature from a quadrature on a reference surface, Lemma 9.4 must be used to account for the transformation of the surface measure; see Exercise A large amount of literature is devoted to quadratures; see Abramowitz and Stegun [, Chap. 25], Hammer and Stroud [287], Stroud [467], Davis and Rabinowitz [90], Brass and Petras [96]. we refer the reader to 6..2 for onedimensional quadratures. Example 26.4 (Cuboids). Quadratures on cuboids can be deduced from one-dimensional quadratures by taking the Gauss nodes in tensor-product form. Note though that tensor-product formulas are not optimal in the sense of using the fewest function evaluations for a given order. Although no general formula for non-tensor-product quadratures for the cube is known, many quasi-optimal quadratures are available in the literature, see e.g., Cools [7, 2.3.] and Cools and Rabinowitz [72, 4.]. Example 26.5 (Quadratures on simplices). Table 26. lists some quadratures on triangles in two dimensions. In this table, we call multiplicity the number of permutations that must be performed on the barycentric coordinates to obtain the list of all the Gauss nodes of the quadrature. For instance, the first-order formula in the second line has three Gauss nodes; the corresponding barycentric coordinates and weights are {,0,0; 3 S}, {0,,0; 3 S}, {0,0,; 3 S}, where S denotes the surface of the triangle. Table 26.2 lists some quadratures on a tetrahedron. As in two dimensions, the multiplicity is the number of permutations to perform on the barycentric coordinates to obtain all the Gauss
3 Part V. Elliptic PDEs 375 nodes of the quadrature. For instance, the third-order formula has five Gauss nodes which are the node ( 4, 4, 4, 4 ) with the weight 4 5V and the four nodes ( 6, 6, 6, 2 ), ( 6, 6, 2, 6 ), ( 6, 2, 6, 6 ), ( 2, 6, 6, 6 ) with the weight 9 20V; here, V denotesthe volumeofthetetrahedron.aformulavalidonanysimplexin R d is K λα0 0...λα d d dx = K α 0!...α d!d! (α α d +d)!, where {λ 0,...,λ d } are the barycentric coordinates in K and α 0,...,α d are natural numbers. This formula is useful to numerically verify the order of quadratures. k q l q Barycentric coord. Multiplicity Weights ω l (,, 3 3 3) S 3 S 3 (,0,0) (,, ) 3 3 ( 2 3,,0) ( 3 4,, 3 3 3) 9 S 6 (,, ) 3 48 ( 3 7,, ) 20 (,,0) (,0,0) 3 S 5 S (a i,a i, 2a i) for i =,2 3 ω i for i =,2 5 7 a = ω = S a 2 = ω 2 = S (,, 3 3 3) (a i,a i, 2a i) for i =,2 3 a = a 2 = (a i,a i, 2a i) for i =, S 55 5 S S 200 a = S a 2 = S (a,b, a b) 6 a = S b = Table 26.. Nodes and weights for quadratures on a triangle of area S.
4 376 Chapter 26. Quadratures k q l q Barycentric coord. Multiplicity Weights ω l (,,, ) V 4 (,0,0,0) (a,a,a, 3a) a = (,,0,0) V 4 V 5 V (,0,0,0) 4 V 20 (,,, ) 4 V 5 (,,, ) 4 V 20 (,,, ) (a i,a i,a i, 3a i) for i =,2 4 a = a 2 = (a,a, a, a) a = V V V V Table Nodes and weights for quadratures on a tetrahedron of volume V Discrete problem with quadratures We now study how quadratures impact the well-posedness of a discrete model problem and the convergence properties of its finite element approximation. This is a crucial question since quadratures cannot be avoided in practice. For simplicity, we focus on the purely diffusive, homogeneous Dirichlet problem. As in 25.,this problem is approximatedusing H -conformingfinite elements.the discretespacev h isbuilt usinganaffinemeshfromashape-regular mesh sequence and a reference finite element ( K, P, Σ) of degree k. The noveltyisthat we nowtakeinto accountthe useofquadraturesto evaluatethe bilinear form a(v,w) = D v wdx and the linear form l(w) = D fwdx Quadrature error For all K T h, we generate a quadrature in each mesh cell K T h with nodes {ξ lk } l {:lq} and weights {ω lk } l {:lq} from a reference quadrature of order k q. For any function φ that is smooth enough to have point values, we define the quadrature error as follows: E K (φ) := K φ(x)dx l q l= ω lk φ(ξ lk ). (26.3)
5 Part V. Elliptic PDEs 377 Lemma 26.6 (Quadrature error with polynomial factor). Let m, n N, r [, ] and assume that k q n+m. Assume that the mesh sequence (T h ) h>0 is affine. Then, there is c, uniform, such that E K (φp) ch m K φ W m, (K) p L (K), (26.4) for all φ W m, (K), all p such that p T K P P n,d, and all K T h. Proof. Let K be the reference element. Let p such that p := p T K P n,d. Let us also denote φ = φ T K. Since T h is affine, after making a change of variable we obtain E K (φp) = det(j K ) Ê( φ p) (with obvious notation for Ê). The linear form φ Ê( φ p) is bounded on W m, ( K) since Ê( φ p) c φ C0 ( K) p C 0 ( K) c φ W m, ( K) p L ( K), where we have used (.8) and norm equivalence in P for p (since dim( P) < )). Moreover, Ê( φ p) = Ê(( φ ĝ) p) for all ĝ P m,d, since p P n,d and k q n+m. Owing to the Bramble Hilbert Lemma 0.8, we infer that Ê( φ p) c inf ĝ P n,d φ ĝ W m, ( K) p L ( K). c φ W m, ( K) p L ( K). We conclude by using Lemma Discrete problem and well-posedness For all K T h, we assume that the diffusion coefficients ij and the source term f are continuous in K, and we define a Q (v h,w h ) = ω lk v h (ξ lk ) (ξ lk ) w h (ξ lk ), l Q (w h ) = K T h l {:l q} K T h l {:l q} ω lk f(ξ lk )w h (ξ lk ), for all (v h,w h ) V h V h. Observe that it is not possible in general to extend a Q and l Q to H0 (D), since functions in H 0 (D) are not necessarily continuous. The discrete problem with quadratures is formulated as follows: { Find uh V h such that (26.5) a Q (u h,w h ) = l Q (w h ), w h V h. Lemma 26.7 (Well-posedness). Assume that P Pl,d for some integer l. Assume that k q 2l 2 and that W, (T h ;R d d ). Then, the discrete problem (26.5) is well-posed for h small enough.
6 378 Chapter 26. Quadratures Proof. We need to prove that there is α Q > 0 such that α Q v h H (D) sup wh V h a Q(v h,w h ) w h L 2 (D) for all v h V h. Since a Q (v h,w h ) a(v h,w h ) (a a Q )(v h,w h ) sup sup sup, w h V h w h L 2 (D) w h V h w h L 2 (D) w h V h w h L 2 (D) and since α v h L 2 (D) sup wh V h a(v h,w h ) w h L 2 (D) with α > 0 owing to the coercivity of a on V h, the assertion follows if we prove that (a a Q )(v h,w h ) ch v h L 2 (D) w h L 2 (D). We observe that (a a Q )(v h,w h ) = K T h E K ( v h w h ). For all i,j {:d}, we use Lemma 26.6 with φ = ij, p = i v h j w h, m =, and n = 2l 2 (so that n + m k q ). Note that indeed p T K P 2l 2,d, because the mesh is affine, i v h = d i = JT K,ii ( i v h) T K, i v h P l,d and a similar argument holds for j w h. Since i v h j w h L (K) i v h L 2 (K) j w h L 2 (K), we infer that E K ( v h w h ) ch K W, (K;R d d ) v h L 2 (K) w h L 2 (K). The desired bound follows by a discrete Cauchy Schwarz inequality (the constant c depends on W, (T h ;R )). d d Error analysis Theorem 26.8 (Error estimate). Assume that P k,d P P l,d for some integers l k. Assume that k q l + k 2. Assume that W k, (T h ;R d d ) and f W k, (T h ). Assume that (26.5) is well-posed and that u H k+ (D). Then, u u h H (D) ch k ( u H k+ (D)(+C Q ( ))+C Q (f)), (26.6) with C Q ( ) = W k, (T h ;R d d ), C Q (f) = max(c P,D l D f W k, (T h ), f W k, (T h )), and C P,D from the Poincaré inequality (2.0). Proof. We use Strang s First Lemma with the quasi-interpolation operator I g,av h0 introduced in Using the H -seminorm, (2.4) becomes (u u h ) L2 (D) C (u I g,av h0 u) L 2 (D) + g,av (a a Q )(Ih0 sup u,w h) (l l Q )(w h ), α Q w h V h w h L 2 (D) where α Q results from the well-posedness of the discrete problem. The supremum over w h V h represents the quadrature error that we now bound by estimating each of the two terms separately.
7 Part V. Elliptic PDEs 379 () Bound on (a a Q ). Letting v h = I g,av h0 u, we use Lemma 26.6 with φ = ij i v h,p = j w h,m = k,andn = l (sothat n+m = l+k 2 k q ) for all i,j {:d} to infer that E K ( ij i v h j w h ) ch k K ij i v h W k, (K) j w h L (K). Combining the Leibniz product rule with the inverse inequality (5.2) (with p = and r = 2) leads to This leads to ij i v h W k, (K) c k ij W m, (K) i v h W k m, (K) m=0 c K 2 ij W k, (K) i v h H k (K). E K ( v h w h ) ch k K W k, (K;R d d ) v h Hk (K) w h L 2 (K). As a result, (a a Q )(Ih0 av sup u,w h) w h V h w h L 2 (D) ch k (max K T h W k, (K;R d d )) u H k+ (D), where we have used the estimate I g,av h0 u H k+ (D) c u H k+ (D) which follows from Theorem 8.2. (2) Bound on (f f Q ). We have (l Q l)(w h ) = K T h E K (fw h ). We cannot apply Lemma 26.6 with φ = f, p = w h, m = k, and n = l since n+m = l+k may be larger than k q. Fortunately, we can use instead Exercise 26.2 with µ = k and ν = l (which satisfy ν +µ 2 = l+k 2 k q ) to infer that E K (fw h ) ch k K( f W k, (K) w h L (K) + f W k, (K) w h L (K)). We conclude using inverse inequalities for w h, a discrete Cauchy Schwarz inequality and the global Poincaré inequality. Remark 26.9 (Literature). The above analysis is inspired from Ciarlet [52, 4.], Ciarlet and Raviart [54], Dautray and Lions [88, XII.5]. It is possible to refine the analysis by assuming f W k,q (T h ) with q > d k and q 2, and that a surface quadrature of order k+l is used to approximate boundary integrals resulting from Neumann conditions; see [88, XII.5] Implementation This section addresses practical aspects of quadrature implementation in view of matrix assembling.
8 380 Chapter 26. Quadratures Nodesand weights. Letm {:N c }andletk m bethecorrespondingmesh cell. The nodes and weights of the quadrature are defined in Proposition Recall from Definition?? that the geometric map T Km is built from reference shape functions { ψ n } n {:ngeo}. This leads us to define the double-entry array psi(:n geo,:l q ) such that psi(n,l) = ψ n ( ξ l ). The k-th Cartesian component of the Gauss node ξ lkm = T Km ( ξ l ) is given by n geo (ξ lkm ) k = coord(k, j(n, m)) psi(n, l). n= We also need the the triple-entry arraydpsi dhatk(:d,:n geo,:l q ) providing the derivatives of the geometric shape functions at the Gauss nodes, i.e., dpsi dhatk(k,n,l) = ψ n x k ( ξ l ). Then,theentriesoftheJacobianmatrixJ Km at ξ l canbecomputedasfollows: ( ) n geo J Km ( ξ l ) = coord(k,j(n,m))dpsi dhatk(k 2,n,l), k,k 2 n= for all k,k 2 {:d}. Since the determinant of J Km ( ξ l ) is always multiplied by the weight ω l in the quadratures, it may be useful to store this product once and for all in a double-entry array of weights weight K({:l q },{:N c }): weight K(l,m) = ω l det ( J Km ( ξ l ) ). Notice that when the mesh is affine, the partial derivatives of the shape functions ψ n areconstant in K, so that the size of the arraydpsi dhatkis reduced to d n geo ; moreover, memory space is saved by storing separately the reference quadrature weights ω l and the determinants det(j Km ). Shape functions. Let { θ n } n N be the referenceshape functions. Since these functions and their derivatives generally need to be evaluated many times at the Gauss nodes in K, it can be useful to compute these values once and for all and store them in the double-entry array theta({:n sh },{:l q }) such that theta(n,l) = θ n ( ξ l ), and in the triple-entry array dtheta dhatk({:d},{:n sh },{:l q }) such that dtheta dhatk(k,n,l) = θ n x k ( ξ l ).
9 Part V. Elliptic PDEs 38 Let us assume for simplicity that the linear bijective map used to generate the local shape functions is the pullback by the geometric map. As a result, the values of the local shape functions at the Gauss nodes in the mesh cell K m are given by θ n (ξ lkm ) = θ n ( ξ l ) = theta(n,l), for all n {:n sh }, all l {:l q }, and all m {:N c } (notice that the value of θ n (ξ lkm ) is independent of K m ). Let us now consider the first-order derivatives of the local shape functions at the Gauss nodes. Using the chain rule yields, for all k {:d}, θ n x k (ξ lkm ) = = d k 2= d k 2= θ n ( ξ l ) (T K m ) k2 (ξ lkm ) x k2 x k θ ( ) n ( ξ l ) J K x m ( ξ l ). k2 k 2,k One can evaluate these derivatives once and for all and store them in an array dtheta dk({:d},{:n sh },{:l q },{:N c }) such that dtheta dk(k,n,l,m) = θ n x k (ξ lkm ). Notice that the size of this array, d n sh l q N c, can be extremely large. If the mesh is affine, an alternative strategy is possible once observing that the Jacobian matrix J Km and its inverse do not depend on the Gauss nodes. In this case, one may store the inverse of the Jacobian matrix in a triple-entry array inv jac K({:d},{:d},{:N c }) such that ( [JKm ] ) inv jac K(k,k 2,m) = k,k 2. Each time the quantity θn x k (ξ lkm ) is needed, the following operations must be performed: θ n x k (ξ lkm ) = d dtheta dhatk(k 2,n,l)inv jac K(k 2,k,m). k 2= The array inv jac K has d d N c entries; this number is much smaller than the number of entries of dtheta dk if d n sh l q. In this situation, storing inv jac K will save memory space at the prize of some additional computations. Assembling. Let us first consider the assembling of the stiffness matrix. We consider a slightly more general problem than in 26.2, namely a diffusionadvection-reaction problem for which
10 382 Chapter 26. Quadratures a Q (v h,w h ) = ω lkm A(ξ lkm,v h Km,w h Km ), m {:N c} l {:l q} with A(ξ,ϕ,ψ) = k,k 2 {:d} +ψ(ξ) ϕ x k (ξ) k k 2 (ξ) ψ x k2 (ξ) k {:d} β k (ξ) ϕ x k (ξ)+ϕ(ξ)µ(ξ)ψ(ξ). A general assembling procedure for a stiffness matrix (stored in the CSR format) has been outlined in Algorithm Our goal is now to detail the evaluation of the array tmp used in this algorithm by means of quadratures. We assume for simplicity that the coefficients ( k k 2 ) k,k 2 {:d}, (β k ) k {:d}, and µ are known analytically; see Exercise 26.5 for discrete data. The assembling procedure is shown in Algorithm 26.. Notice that we preliminarily evaluate and store the coordinates of the Gauss nodes ξ lkm since we need to evaluate the values of the coefficients at these nodes. Algorithm 26.. Assembling of A for analytic data. A = 0 for m {:N c} do for l {:l q} do; tmp = 0 for k {:d} do n geo xi l(k) = coord(k,j(n,m)) psi(n,l) n= for ni {:n sh } do for nj {:n sh } do d x = dtheta dk(k,nj,l,m) k k 2 (xi l) dtheta dk(k 2,ni,l,m) k,k 2 = x 2 = theta(ni,l) d β k (xi l) dtheta dk(k,nj,l,m) k = x 3 = theta(ni,l) µ(xi l) theta(nj,l) tmp(ni,nj) = tmp(ni,nj)+[x +x 2 +x 3] weight K(l,m) Accumulate tmp in A as in Algorithm 23.2 The assembling of the right-hand side vector can be performed similarly. Let us consider that l Q (w h ) = m {:N c} l {:l ω q} lk m F(ξ lkm,w h Km )
11 Part V. Elliptic PDEs 383 with F h (ξ,ψ) = f(ξ)ψ(ξ). The assembling procedure is presented in Algorithm Algorithm Assembling of B for analytic data. B = 0 for m {:N c} do for l {:l q} do; tmp = 0 for k {:d} do n geo xi l(k ) = coord(k,j(n,m))psi(n,l) n= for ni {:n sh } do tmp(ni) = tmp(ni)+f(xi l) theta(ni,l) weight K(l,m) for ni {:n sh } do; i = gdof(ni,m) B(i) = B(i)+tmp(ni) Exercises Exercise 26. (Quadrature error). Consider a quadrature oforderk q. Let p [, ] and let m N be such that k q + m > d p. Prove that there is c, uniform with respect to h, such that E K (φ) ch m+d( p ) K φ W m,p (K), for all φ W m,p (K) and all K T h. (Hint: use Lemma??.) Prove that D φ(x)dx K T h l {:l q} ω lk φ(ξ lk ) ch m φ W m,p (D) for all φ W m,p (D) and all h > 0. (Hint: use a discrete Hölder s inequality.) Exercise 26.2 (Quadrature error with polynomial). Let ν, µ and assume that k q ν +µ 2. The goal is to prove that E K (φp) ch µ K ( φ W µ, p L (K) + φ W µ, p L (K)), (26.7) for all φ W µ, (K), all p P ν,d, all K T h, and all h > 0. (Note that φ W µ, p L (K) is not needed if k q ν +µ ; see Lemma 26.6,) (i) Prove that E K (φp K ) ch µ K φ W µ, p L (K), where p K is the meanvalue of p over K.
12 384 Chapter 26. Quadratures (ii) Prove (26.7) for µ 2. (Hint: Apply Lemma 26.6 with m = µ.) (iii) Prove (26.7) for µ =. (Hint: prove first that Ê( φ( p p K)) c φ( p p K) W, ( K) ; note that p K = p K.) Exercise 26.3 (Asssembling). Let D = (0,) 2. Consider the problem u+u = in D and u D = 0. Approximate its solution with P H - conforming finite elements on the two meshes shown on the right. 0 0 Evaluate the discrete solution in both cases. (Hint: there is only one degree of freedom in both cases; see Exercise 25.5 for computing the gradient part of the stiffness coefficient and use a quadrature from Table 26. for the zeroorder term.) For a fine mesh composed of 800 elements, u h ( 2, 2 ) Comment. Exercise 26.4 (Surface quadrature). Let F be a face of a three-dimensional element. Let F R 2 be the reference face and let T F : F F be a geometric transformation for F. Let t (ŝ),t 2 (ŝ) be the two column vectors of the Jacobian matrix of T F (ŝ), say J F (ŝ) = [t (ŝ),t 2 (ŝ)]. (i) Compute the metric tensor F := J T F J F in terms of the dot products t i t j (ŝ), i,j 2. (ii) Show that ds = t (ŝ) t 2 (ŝ) l 2 (R 3 )dŝ (Hint: use Lagrange s identity). (iii) Let {(ŝ l,ŵ l } l {:l q } be a quadrature on F. What is the corresponding quadrature on F? Exercise 26.5 (Assembling with discrete data). Assume that the coefficients ( k k 2 ) k,k 2 {:d}, (β k ) k {:d}, and µ are in the discrete space V h. Let dif, beta, and mu be the corresponding coordinate vectors. Adapt Algorithm 26. to this situation. (Hint: µ(ξ lkm ) = n {:n sh } mu(gdof(n,m)) theta(n,l).) Exercise 26.6 (Assembling of RHS). Write the assembling algorithm for the right-hand side vector in the case where F h (ξ,w h ) = f(ξ)w h (ξ) + d k = β k (ξ) w h x k (ξ) with analytically-known data.
Numerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More 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 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 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 informationChapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method
Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More 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 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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationPDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM.
PDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM. Sergey Korotov Basque Center for Applied Mathematics / IKERBASQUE http://www.bcamath.org & http://www.ikerbasque.net 1 Introduction The analytical solution
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More 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 informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationFinite Element Interpolation
Finite Element Interpolation This chapter introduces the concept of finite elements along with the corresponding interpolation techniques. As an introductory example, we study how to interpolate functions
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
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 informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationFinite Element Method for Ordinary Differential Equations
52 Chapter 4 Finite Element Method for Ordinary Differential Equations In this chapter we consider some simple examples of the finite element method for the approximate solution of ordinary differential
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationChapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction
Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction The a posteriori error estimation of finite element approximations of elliptic boundary value problems has reached
More 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 informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationPart VIII, Chapter 39. Fluctuation-based stabilization Model problem
Part VIII, Capter 39 Fluctuation-based stabilization Tis capter presents a unified analysis of recent stabilization tecniques for te standard Galerkin approximation of first-order PDEs using H 1 - conforming
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationLECTURE 5: THE METHOD OF STATIONARY PHASE
LECTURE 5: THE METHOD OF STATIONARY PHASE Some notions.. A crash course on Fourier transform For j =,, n, j = x j. D j = i j. For any multi-index α = (α,, α n ) N n. α = α + + α n. α! = α! α n!. x α =
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationA multipoint flux mixed finite element method on hexahedra
A multipoint flux mixed finite element method on hexahedra Ross Ingram Mary F. Wheeler Ivan Yotov Abstract We develop a mixed finite element method for elliptic problems on hexahedral grids that reduces
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More information1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:
Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not
More informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More informationThomas Apel 1, Ariel L. Lombardi 2 and Max Winkler 1
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANISOTROPIC MESH REFINEMENT IN POLYHEDRAL DOMAINS: ERROR ESTIMATES WITH DATA IN
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationControllability of the linear 1D wave equation with inner moving for
Controllability of the linear D wave equation with inner moving forces ARNAUD MÜNCH Université Blaise Pascal - Clermont-Ferrand - France Toulouse, May 7, 4 joint work with CARLOS CASTRO (Madrid) and NICOLAE
More informationTheory of PDE Homework 2
Theory of PDE Homework 2 Adrienne Sands April 18, 2017 In the following exercises we assume the coefficients of the various PDE are smooth and satisfy the uniform ellipticity condition. R n is always an
More informationPart III. 10 Topological Space Basics. Topological Spaces
Part III 10 Topological Space Basics Topological Spaces Using the metric space results above as motivation we will axiomatize the notion of being an open set to more general settings. Definition 10.1.
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More 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 informationDIFFERENTIAL EQUATIONS
NUMERICAL PARTIAL DIFFERENTIAL EQUATIONS Lecture Notes, Winter 2013/14 Christian Clason July 15, 2014 Institute for Mathematics and Scientific Computing Karl-Franzens-Universität Graz CONTENTS I BACKGROUND
More informationNUMERICAL PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL PARTIAL DIFFERENTIAL EQUATIONS Lecture Notes, Winter 2011/12 Christian Clason January 31, 2012 Institute for Mathematics and Scientific Computing Karl-Franzens-Universität Graz CONTENTS I BACKGROUND
More informationAn interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes
An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type interpolation
More informationStatistics 612: L p spaces, metrics on spaces of probabilites, and connections to estimation
Statistics 62: L p spaces, metrics on spaces of probabilites, and connections to estimation Moulinath Banerjee December 6, 2006 L p spaces and Hilbert spaces We first formally define L p spaces. Consider
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationJUHA KINNUNEN. Harmonic Analysis
JUHA KINNUNEN Harmonic Analysis Department of Mathematics and Systems Analysis, Aalto University 27 Contents Calderón-Zygmund decomposition. Dyadic subcubes of a cube.........................2 Dyadic cubes
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationApproximation in Banach Spaces by Galerkin Methods
2 Approximation in Banach Spaces by Galerkin Methods In this chapter, we consider an abstract linear problem which serves as a generic model for engineering applications. Our first goal is to specify the
More informationJuan Vicente Gutiérrez Santacreu Rafael Rodríguez Galván. Departamento de Matemática Aplicada I Universidad de Sevilla
Doc-Course: Partial Differential Equations: Analysis, Numerics and Control Research Unit 3: Numerical Methods for PDEs Part I: Finite Element Method: Elliptic and Parabolic Equations Juan Vicente Gutiérrez
More informationWeak Formulation of Elliptic BVP s
Weak Formulation of Elliptic BVP s There are a large number of problems of physical interest that can be formulated in the abstract setting in which the Lax-Milgram lemma is applied to an equation expressed
More informationCONVERGENCE ANALYSIS OF A BALANCING DOMAIN DECOMPOSITION METHOD FOR SOLVING A CLASS OF INDEFINITE LINEAR SYSTEMS
CONVERGENCE ANALYSIS OF A BALANCING DOMAIN DECOMPOSITION METHOD FOR SOLVING A CLASS OF INDEFINITE LINEAR SYSTEMS JING LI AND XUEMIN TU Abstract A variant of balancing domain decomposition method by constraints
More informationBoundary Value Problems and Iterative Methods for Linear Systems
Boundary Value Problems and Iterative Methods for Linear Systems 1. Equilibrium Problems 1.1. Abstract setting We want to find a displacement u V. Here V is a complete vector space with a norm v V. In
More information1 Directional Derivatives and Differentiability
Wednesday, January 18, 2012 1 Directional Derivatives and Differentiability Let E R N, let f : E R and let x 0 E. Given a direction v R N, let L be the line through x 0 in the direction v, that is, L :=
More informationFinite Element Spectral Approximation with Numerical Integration for the Biharmonic Eigenvalue Problem
City University of New York (CUNY) CUNY Academic Works Publications and Research Kingsborough Community College 2014 Finite Element Spectral Approximation with Numerical Integration for the Biharmonic
More informationA WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS
A WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS LIN MU, JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This article introduces and analyzes a weak Galerkin mixed finite element method
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 informationSPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS
SPECTRAL THEOREM FOR COMPACT SELF-ADJOINT OPERATORS G. RAMESH Contents Introduction 1 1. Bounded Operators 1 1.3. Examples 3 2. Compact Operators 5 2.1. Properties 6 3. The Spectral Theorem 9 3.3. Self-adjoint
More informationA REVIEW OF OPTIMIZATION
1 OPTIMAL DESIGN OF STRUCTURES (MAP 562) G. ALLAIRE December 17th, 2014 Department of Applied Mathematics, Ecole Polytechnique CHAPTER III A REVIEW OF OPTIMIZATION 2 DEFINITIONS Let V be a Banach space,
More informationRepresentations of moderate growth Paul Garrett 1. Constructing norms on groups
(December 31, 2004) Representations of moderate growth Paul Garrett Representations of reductive real Lie groups on Banach spaces, and on the smooth vectors in Banach space representations,
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 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 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 informationFINITE ELEMENT METHODS
FINITE ELEMENT METHODS Lecture notes arxiv:1709.08618v1 [math.na] 25 Sep 2017 Christian Clason September 25, 2017 christian.clason@uni-due.de https://udue.de/clason CONTENTS I BACKGROUND 1 overview of
More informationINTEGRATION ON MANIFOLDS and GAUSS-GREEN THEOREM
INTEGRATION ON MANIFOLS and GAUSS-GREEN THEOREM 1. Schwarz s paradox. Recall that for curves one defines length via polygonal approximation by line segments: a continuous curve γ : [a, b] R n is rectifiable
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 informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationLecture Notes on Metric Spaces
Lecture Notes on Metric Spaces Math 117: Summer 2007 John Douglas Moore Our goal of these notes is to explain a few facts regarding metric spaces not included in the first few chapters of the text [1],
More informationOne-dimensional and nonlinear problems
Solving PDE s with FEniCS One-dimensional and nonlinear problems L. Ridgway Scott The Institute for Biophysical Dynamics, The Computation Institute, and the Departments of Computer Science and Mathematics,
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 informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationThe mortar element method for quasilinear elliptic boundary value problems
The mortar element method for quasilinear elliptic boundary value problems Leszek Marcinkowski 1 Abstract We consider a discretization of quasilinear elliptic boundary value problems by the mortar version
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationLecture 4 Lebesgue spaces and inequalities
Lecture 4: Lebesgue spaces and inequalities 1 of 10 Course: Theory of Probability I Term: Fall 2013 Instructor: Gordan Zitkovic Lecture 4 Lebesgue spaces and inequalities Lebesgue spaces We have seen how
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
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 informationTraces and Duality Lemma
Traces and Duality Lemma Recall the duality lemma with H / ( ) := γ 0 (H ()) defined as the trace space of H () endowed with minimal extension norm; i.e., for w H / ( ) L ( ), w H / ( ) = min{ ŵ H () ŵ
More informationICES REPORT Direct Serendipity Finite Elements on Convex Quadrilaterals
ICES REPORT 17-8 October 017 Direct Serendipity Finite Elements on Convex Quadrilaterals by Todd Arbogast and Zhen Tao The Institute for Computational Engineering and Sciences The University of Texas at
More informationNORMS ON SPACE OF MATRICES
NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system
More informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
More informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
More informationEfficient hp-finite elements
Efficient hp-finite elements Ammon Washburn July 31, 2015 Abstract Ways to make an hp-finite element method efficient are presented. Standard FEMs and hp-fems with their various strengths and weaknesses
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationConvergence analysis of a balancing domain decomposition method for solving a class of indefinite linear systems
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2000; 00:1 6 [Version: 2002/09/18 v1.02] Convergence analysis of a balancing domain decomposition method for solving a class of indefinite
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More informationj=1 [We will show that the triangle inequality holds for each p-norm in Chapter 3 Section 6.] The 1-norm is A F = tr(a H A).
Math 344 Lecture #19 3.5 Normed Linear Spaces Definition 3.5.1. A seminorm on a vector space V over F is a map : V R that for all x, y V and for all α F satisfies (i) x 0 (positivity), (ii) αx = α x (scale
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationA non-standard Finite Element Method based on boundary integral operators
A non-standard Finite Element Method based on boundary integral operators Clemens Hofreither Ulrich Langer Clemens Pechstein June 30, 2010 supported by Outline 1 Method description Motivation Variational
More information