arxiv: v1 [math.na] 27 Jan 2014
|
|
- Trevor Kennedy
- 5 years ago
- Views:
Transcription
1 L 2 -ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF BOUNDARY FLUXES MATS G. LARSON AND ANDRÉ MASSING arxiv: v1 [mat.na] 27 Jan 2014 Abstract. We prove quasi-optimal a priori error estimates for finite element approximations of boundary normal fluxes in te L 2 -norm. Our results are valid for a variety of different scemes for weakly enforcing Diriclet boundary conditions including Nitsce s metod, and Lagrange multiplier metods. Te proof is based on an error representation formula tat is derived by using a discrete dual problem wit L 2 -Diriclet boundary data and combines a weigted discrete stability estimate for te dual problem wit anisotropic interpolation estimates in te boundary zone. Key words. Boundary flux, L 2 -error estimates, discete dual problem, Nitsce s metod, Lagrange multipliers AMS subject classifications. 65N12, 65N15, 65N30 1. Introduction. Te normal flux at te boundary or on interior interfaces is in general of great interest in applications. Examples include surface stresses in mecanics, eat transfer troug interfaces, and transport of fluid in Darcy flow. Recently Melenk and Wolmut [15] as sown quasi-optimal order estimates for fluxes in a mortar setting were continuity and boundary conditions is enforced using a mortaring space of Lagrange multipliers. More precisely, tey sown tat te L 2 -norm of te error in te normal flux is of order ln for piecewise linear polynomials and of order k for piecewise polynomials of order k. In contrast only k 1/2 will be obtained if a trace inequality is used in combination wit standard convergence teory for saddle point problems, see [4]. In tis contribution we give an alternative proof of tis result and, focusing on te case k = 1, we also consider a wider variety of metods for weakly enforcing Dirclet boundary conditions, including Nitsce s metod and stable and stabilized Lagrange multipliers metods. Our proof is based on an error representation formula were te error in te normal flux is represented in terms of te interpolation error and te solution to a discrete dual problem wit L 2 -Diriclet boundary data. Key to te error estimate is a stability estimate for te discrete dual problem in terms of te L 2 -norm of te Diriclet data. In te continuous case suc an estimate is known, see Cabrowsky [8] and [9], and provides control of te gradient weigted wit te distance to te boundary as well as a max-norm control of te L 2 -norms of te solution on manifolds close to and parallel wit te boundary. We prove a corresponding stability estimate for our discrete dual problem. In contrast to te approac by Melenk and Wolmut [15], we avoid using a Besov space framework. Our error representation formula is related to te one derived in te Carey et al. [7], Giles et al. [12], Pelivanov et al. [17] were various estimates for functionals of te normal flux are derived and [10] were adaptive metods based on dual problems targeting te flux in a coupled problem are developed. Note owever tat in our setting were we seek an a priori estimate, we employ a discrete dual problem wile in te a posteriori setting, te corresponding continuous dual problem is used. Here we also establis te stability of te discrete dual problem using analytical tecniques wile in te duality based a posteriori error estimates, stability is often estimated using computational tecniques or a known analytical stability result. Te remainder of tis work is organized as follows. In Section 2 we introduce te model problem and its variational formulations we will consider trougout tis work. Corresponding finite element discretizations are presented in Section 3 togeter wit te definition of te discrete boundary fluxes. In Section 5 we prove stability bounds for te discrete dual problem and provide interpolation estimates of te solution close to te boundary. Combining tese results allows us to prove L 2 -error estimates for te boundary flux approximations in Section 6. In Section 7, we finally present numerical results illustrating te teoretical findings. Department of Matematics, Umeå University, SE Umeå, Sweden. mats.larson@mat.umu.se Simula Researc Laboratory, P.O. Box 134, 1325 Lysaker, Norway. massing@simula.no 1
2 2. Model Problem. Let Ω be a polygonal domain in R d, d = 2, 3 wit boundary Ω. We consider te elliptic model problem: find u : Ω R suc tat u = f in Ω 2.1 u = g on Ω 2.2 were f and g are given data. Ten te boundary flux σ n for te solution u is defined by σ n = n u 2.3 were n is te outwards pointing unit normal to Ω. In wat follows, we consider te standard Sobolev spaces H s U, s 0 on some domain U, endowed te te usual norms s,u and semi-norms s,u. More generally, te space W s,p U is defined as te Sobolev space consisting of all functions aving p-integrable derivates up to order s on U. As usual, H 0 Ω = L 2 Ω and H 1/2 Ω denotes te dual space of H 1/2 Ω. Moreover, for a function g H 1/2 Ω we introduce te notation H 1 g Ω = {v H 1 Ω : v Ω = g}. Te scalar product in H s U is written as, s,u and to simplify te notation, we generally omit te domain designation if U = Ω and te Sobolev index if s = 0 in bot norm and scalar product expressions. Using tis notation, a weak formulation of te elliptic boundary value problem is to seek u H 1 g Ω suc tat were au, v = lv v H 1 0 Ω 2.4 au, v = u, v 2.5 lv = f, v 2.6 Here, te boundary condition u Ω = g is already incorporated into te trial space H 1 g Ω. Alternatively, te boundary condition 2.2 can be enforced weakly by using a Lagrange multiplier approac. Introducing te bilinear form bµ, v = µ, v Ω 2.7 te resulting variational formulation is given by te saddle point problem: find u, λ H 1 Ω H 1/2 Ω suc tat au, v + bλ, v + bµ, u = lv + bµ, g v, µ H 1 Ω H 1/2 Ω 2.8 For brevity, we migt denote te left-and side by Au, λ; v, µ and te rigt-and side Lv, µ. It is well-known [1, 5, 18, 22], tat te saddle point problem 2.8 satisfies te Babuška-Brezzi condition, in particular bλ, v sup µ 1/2, Ω µ H 1/2 Ω 2.9 v H 1 Ω\{0} v 1,Ω Consequently, problem 2.8 possesses a unique solution u, λ, were u solves in a weak sense and te Lagrange multiplier λ represents te negative of te normal flux of u, i.e. λ = σ n. 3. Finite Element Discretizations of te Model Problem. In tis section, we introduce te finite element discretizations of problem we will consider trougout tis work. Te discretizations are defined on a quasi-uniform partition T of Ω into sape regular triangles in two or tetraedra in tree space dimensions wit mes parameter. For a given mes T, let te associated finite element space of piecewise linear continuous functions be denoted by V. We do not assume V H 1 g Ω and consequently, te discretizations to be considered will enforce te boundary condition 2.2 weakly. For eac discretization we will define a discrete counterpart Σ n of te boundary flux
3 3.1. Nitsce s Metod. Te Nitsce [16] finite element metod takes te form: find u V suc tat were te forms are defined by a u, v = l v v V 3.1 a u, v = au, v n u, v Ω n u, v Ω + β 1 u, v Ω 3.2 l v = lv g, n v Ω + β 1 g, v Ω 3.3 wit β being a positive parameter. Introducing te energy norm we recall tat te bilinear form a, is continuous v 2 = v 2 Ω + n v 2 Ω + 1 v 2 Ω 3.4 a u, v u v 3.5 and tat if te stabilization parameter β is large enoug, a coercivity condition is satisfied, yielding te standard error estimate v 2 a v, v v V 3.6 u u u Here and trougout, we use te notation a b for a Cb for some generic constant C wic vary wit te context but is always independent of te mes size. For proofs of 3.6 and 3.7, we refer to [14, 16]. To Nitsce s metod 3.1, we associate te discrete variational normal flux of te form Σ n, v Ω = u, v Ω u g, n v Ω f, v Ω v V 3.8 were Σ n is te so-called Nitsce flux Σ n = n u β 1 u g Lagrange Multiplier Metod. To formulate a finite element discretization of te saddle point problem 2.8, we assume tat a discrete function space Λ L 2 Ω H 1/2 Ω is given, and we equip V, Λ and te total approximation space V Λ wit te natural norms u 2 = u 2 Ω + 1/2 u 2 Ω 3.10 λ 2 = 1/2 λ 2 Ω 3.11 u, λ 2 = u 2 + λ respectively, see Pitkäranta [18]. Employing te discrete norms v and µ, it is well-known [18, 19] tat te approximation space Λ as to be designed carefully in order to satisfy te discrete equivalent of te inf-sup condition 2.9. Terefore, a stabilized Lagrange multiplier metod as been proposed by Barbosa and Huges [2, 3] were residual terms were added to circumvent te inf-sup condition Recently, a generalized approac based on projection stabilized Lagrange multipliers as been proposed by Burman [6]. To cover a broad range of stable and stabilized Lagrange multiplier metods, we assume tat te discrete saddle point problem is of te following form: find u, λ V Λ suc tat A u, λ ; v, µ = L v, µ v, µ V Λ
4 were A u, λ; v, µ = au, v + bλ, v + bµ, u c u, λ; v, µ 3.14 L v, µ = lv + g, µ Ω 3.15 Ten, te approximation of te normal flux 2.3 is naturally defined by te negative of te discrete Lagrange multiplier: Σ n = λ 3.16 In te variational form 3.14, te bilinear form c, represents a consistent, possibly vanising stabilization form suc te inf-sup condition olds, as well as te continuity condition and te error estimate A u, λ; v, µ sup u, λ 3.17 v,µ V Λ \{0,0} v, µ A u, λ; v, µ u, λ v, µ 3.18 u u, λ λ u 2,Ω + 3/2 λ 1, Ω 3.19 Well-known Lagrange multiplier discretizations wic are covered by tese assumptions are described and analyzed in [18, 19] and [2, 3, 22]. In [18, 19], Pitkäranta proved certain local stability conditions, rougly stating tat te pairing P 1 c T P 0 dc Γ H is stable, if te mes size H of a given discretization Γ H of te boundary Ω = Γ satisfies te condition ch for some c > 1. To avoid additional mesing of te boundary and to use te natural space Λ = {µ L 2 Ω µ P 0 F F T } defined on te trace mes T, a stabilized symmetric Lagrange Multiplier approac was proposed by Barbosa and Huges [2, 3]. Stenberg [22] simplified te approac by sowing tat te weak formulation 3.13 combined wit te stabilization form c u, λ; v, µ = αλ + n u, µ + n v Ω 3.20 satisfies te inf-sup condition 3.17, te continuity condition 3.18 and tus te error estimate 3.19 wen 0 < α < C I, wit C I being te constant in Finally, we would like to mention te general approac by Burman [6]. In tis metod, te stabilization operator is given by some symmetric form c λ, µ wic, rougly speaking, controls te distance between a given discretization space Λ and anoter discrete space L were V L presents an inf-sup stable pairing. Generally, te stabilization form is only required to be optimal weakly consistent and to not clutter te presentation, we skip te details for te trivial adaption of our approac to tis variant. 4. Error Representation Formulas. In tis section, we establis te error representation formulas for te discrete boundary fluxes. Te representation formula will later allow us to bound te L 2 -error of te boundary flux approximations in terms of interpolation errors and a stability estimate for te discrete solution to a suitable dual problem Nitsce s Metod. For given boundary data ψ L 2 Ω, we define te discrete dual problem for Nitsce s metod as follows: find φ V suc tat were a, is defined in 3.2 and a v, φ = m ψ, v v V 4.1 m ψ, v = β 1 ψ, v Ω ψ, n v Ω 4.2 4
5 Setting v = e = π u u we obtain a e, φ = m ψ, e 4.3 Using Galerkin ortogonality, we note tat te left and side can be written and for te rigt and side were te second term takes te form a e, φ = a π u u, φ 4.4 m ψ, e = m ψ, π u u + m ψ, u u 4.5 m ψ, u u = β 1 u u, ψ n u u, ψ Ω 4.6 = β 1 g u + n u n u, ψ Ω 4.7 = Σ n u σ n u, ψ Ω 4.8 Collecting tese identities, we arrive at te error representation formula σ n u Σ n u, ψ Ω = a u π u, φ m ψ, u π u 4.9 Tus we ave te following Lemma 4.1. Wit σ n u and Σ n u defined by 2.3 and 3.8 it olds σ n u Σ n u Ω sup ψ L 2 Ω\{0} 1 a u π u, φ + m ψ, u π u ψ Ω Lagrange Multiplier Metod. We consider te following discrete dual problem: find φ, θ V Λ suc tat were A, is defined as in 3.14 and A v, µ; φ, θ = m ψ, µ v, µ V Λ 4.11 m ψ, µ = ψ, µ Ω 4.12 Setting v, µ = π u u, π λ λ and using Galerkin ortogonality, we obtain m ψ, π λ λ = A π u u, π λ λ ; φ, θ = A π u u, π λ λ; φ, θ If we write λ λ, ψ = m ψ, λ π λ + m ψ, π λ λ, we arrive at an error representation form similar to 4.9: λ λ, ψ Ω = A π u u, π λ λ; φ, θ m ψ, λ π λ Consequently, te flux error σ n u Σ n u Ω = λ λ Ω can be estimated via following Lemma 4.2. It olds λ λ Ω sup ψ L 2 Ω\{0} 1 A π u u, π λ λ; φ, θ + m ψ, λ π λ ψ Ω
6 5. Stability Bounds for te Discrete Dual Problem. From te error representation formula stated in Lemma 4.1 and Lemma 4.2, we note tat in order to prove estimates for te flux in te L 2 -norm, we need to consider stability bounds in terms of te L 2 -norm of ψ. Cabrowski [9] proved suc estimates for te corresponding continuous problem: find φ : Ω R suc tat φ = 0 in Ω 5.1 u = ψ on Ω 5.2 wit ψ L 2 Ω. To state te basic energy type estimate, we sall introduce some notation tat will also be needed in our fortcoming developments. Let ρx = distx, Ω be te minimal distance between x Ω and Ω and px Ω be te point closest to x Ω. We note tat px = x + npxρx, were npx is te exterior unit normal to Ω at px, and tat tere is a constant δ 0 > 0, only dependent on te curvature of te boundary, suc tat for eac x Ω wit ρx δ 0 tere is a unique px Ω. Next, we define te sets Ω δ = {x Ω : ρx > δ} 5.3 were 0 δ δ 0, and we note tat te closest point mapping p : Ω δ Ω is a bijection wit inverse denoted by p 1 δ. Referring to [11, Lemma 14.16], we recall tat tat ρ Ck Ω δ for k 2 for δ 0 cosen small enoug. If we define a weigted norm v ρ,ω by v 2 ρ,ω = v 2 ρ dx ten Cabrowski [9] proved te following result for te continuous problem: if φ W 1,2 loc satisfies problem in te sense tat ten φ p 1 δ ψ Ω 0 wen δ 0 + φ 2 ρ,ω + φ 2 ρ,ω + sup φ 2 Ω δ ψ 2 Ω. 0 δ δ 0 We sall now prove a corresponding estimate for te discrete dual problems 4.1 and In order to formulate our results, we introduce te sifted weigt function and we let Ω ρ δ = max0, ρ δ, δ δ δ δ = C 5.5 wit te constant C > 0 cosen suc tat ρ δ = 0 on all elements wit a face on te boundary Ω, see Figure 5.2. Te existence of suc a constant C follows from te assumed quasi-uniformity of te mes. In te case were Ω is not a C 2 -domain but rater a convex polyedral domain described by faces {F i } N i=1, we define stripes S δf i = {x R : x 0 F i t s.t. x = x 0 + t n 0 t δ}, cf. Figure 5.1. Ten te analysis presenting in tis work carries over by considering eac stripe at a time and te fact tat locally only a finite number of stripes overlaps. We state now te main result of tis section. Proposition 5.1. Let φ V be te solution of te discrete dual problem 4.1. Ten φ satisfies te stability estimate φ 2 ρ δ,ω + φ 2 Ω + sup φ 2 Ω δ + φ 2 Ω ψ 2 Ω δ δ 0 Alternatively, if φ, θ V Λ is te solution of te discrete dual problem 4.11, ten φ 2 ρ δ,ω + φ 2 Ω + 2 θ 2 Ω + sup φ 2 Ω δ + φ 2 Ω ψ 2 Ω δ δ 0 Before we present te elaborated proof of Proposition 5.1 in Section 5.2, te next section collects useful inequalities and interpolation estimates wic will be used trougout te remaining work. 6
7 F 1 δ F 2 F 6 Ω F 3 F 5 F 4 Fig Polyedral domain Ω wit faces {F i } 6 i=1 wit corresponding stripes S δf i carted in different colors Interpolation Error Estimates. We recall te following trace inequality for v H 1 Ω: v T 1/2 T v T + 1/2 T v T T T 5.8 v T Ω 1/2 T v T + 1/2 T v T T T 5.9 See Hansbo and Hansbo [13] for a proof of 5.9. inverse estimates for v V : We will also need te following well-known v T 1 T v T T T /2 n v F v T T T 5.11 Let π : L 2 Ω V be te standard Scott Zang interpolation operator [21] and recall te interpolation error estimates v π v r,t s r v s,ωt 0 r s 2 T T 5.12 v π v r,f s r 1/2 v s,ωt 0 r s 2 F F 5.13 were ωt is te patc of neigbors of element T ; tat is, te domain consisting of all elements saring a vertex wit T. Recalling te definition 5.3 of Ω δ, we introduce te -band T Ωδ for a mes T by T Ωδ = {T T : T Ω δ } 5.14 Tis is illustrated in Figure 5.2. We note tat tanks to te quasi-uniformity T Ωδ d Ω δ d 1 wit d and d 1 denoting te volume and area of te corresponding sets. Te trace inequality 5.9 allows to generalize te interpolation estimate 5.13 to If we in addition assume tat v π v r,t Ωδ s r 1/2 v s,ωt 0 r s 2 T T 5.15 sup D s u Ωδ ωt δ δ 1 7
8 for some δ 1 suc tat ωt 0 δ δ 1 Ω δ, an order 1/2 can be recovered in estimate 5.13 and 5.15 by applying Hölder s inequality in normal direction to Ω δ : v π v r,t Ωδ s r sup 0 δ δ 0 D s u Ωδ ωt 0 r s 2 T T 5.17 We summarize our observations in te following global, anisotropic interpolation estimate: Proposition 5.2. Let u H 2 Ω and suppose tat sup 0 δ δ1 2 u Ωδ 1 for some δ 1 suc tat 0 δ δ 0 T Ωδ 0 δ δ 1 Ω δ. Ten te interpolation error satisfies sup u π u Ωδ + sup u π u Ωδ 2 sup 2 u Ωδ δ δ 0 0 δ δ 0 0 δ δ 1 Note tat te previous interpolation estimates olds if u W 2, Ω is te finite element solution of 2.4 wit strongly imposed boundary conditions, see [20]. Here owever, we only require tat, rougly speaking, 2 u L 2 on manifolds close and parallel to te boundary and 2 u L in normal direction as quantified by assumption Weigted Energy Stability. In tis section, we finally prove Proposition 5.1. Te main idea of te proof is to divide te domain into an interior region and a boundary layer of tickness O. Away from te boundary, a weigted stability estimate can be proven by testing te discrete dual problems wit a weigted test function. Tis function is cosen suc tat it is identically zero in a layer of elements next to te boundary and tus te boundary terms in te discrete bilinear forms vanis. Since te desired weigted test function does not reside in V we approximate it wit a Lagrange interpolant and estimate te reminder. Witin te boundary layer, an estimate for te discrete energy stability emanating from te coercivity of te finite element metod is establised. Tis stability scales wit since te boundary data only resides in L 2 but it olds all te way out to te boundary. More specifically, te following lemma olds: Lemma 5.3 Discrete Energy Stability. Let φ V be te solution of te discrete dual problem Ten for any κ 0 it olds φ 2 Ω + 2 n φ 2 Ω + κ φ 2 Ω + sup φ 2 Ω δ ψ 2 Ω δ δ Alternatively, assume φ, θ V Q is te solution of te discrete dual problem Ten for any κ 0 it olds φ 2 Ω + 2 θ 2 Ω + κ φ 2 Ω + sup φ 2 Ω δ C ψ 2 Ω δ δ Proof. We note tat te estimate φ 2 = φ 2 Ω + κ φ 2 Ω + 2 n φ 2 Ω + φ 2 Ω ψ 2 Ω 5.21 follows directly by setting v = φ in 5.26, using coercivity 3.1, and multiplying by. Furtermore, wit 0 δ δ 0 we find tat Tus for 0 δ δ, we ave φ 2 Ω δ φ 2 Ω + δ φ 2 Ω\Ω δ 5.22 sup φ 2 Ω δ φ 2 Ω + φ 2 Ω\Ω δ ψ 2 Ω δ δ Combining 5.21 and 5.23 we arrive at te desired estimate. Te second estimate 5.20 can be sown similarly. Setting v, µ = φ, θ in 4.11, using te inf-sup condition 3.17 and multiplying wit, we directly obtain φ, θ 2 ψ 2 Ω
9 In particular, we ave φ 2 Ω + 2 θ 2 Ω + κ φ 2 Ω + φ 2 Ω ψ 2 Ω 5.25 and using te estimate 5.22 once more, we arrive at te desired estimate. Proposition 5.4. If φ V satisfies or φ, θ V Λ satifies a v, φ + κv, φ = m ψ, v v V 5.26 A v, µ; φ, θ + κv, φ = m ψ, v, µ v, µ V Λ 5.27 wit a constant large enoug parameter κ > 0. estimate Ten, in bot cases, φ satisfies te stability φ 2 ρ δ,ω + φ 2 Ω + sup φ 2 Ω δ + φ 2 Ω ψ 2 Ω δ δ 0 Proof. First, we note tat discrete energy stability estimate provides sufficient control for δ. Let now δ be cosen suc tat 0 < δ < δ. Coosing te test function v = I ρ δ φ = ρ δ φ + I Iρ δ φ, δ δ δ were I is te Lagrange interpolant, in 5.26 we obtain te identity 0 = a φ, I ρ δ φ + κφ, I ρ δ φ Ω 5.30 = φ, I ρ δ φ Ω + κφ, I ρ δ φ Ω = φ, I Iρ δ φ + κφ, I Iρ δ φ Ω }{{} I + φ, ρ δ φ + κφ, ρ δ φ Ω }{{} II = I + II 5.31 Note tat, due to our coice of δ, I ρ δ φ = 0 on all elements wit a face on Ω and tus m ψ and te boundary terms in a, and vanis. Term I. We divide te set of elements in te mes T into tree disjoint subsets T 0 = {T T : ρ δ = 0 on T } T Ωδ = {T T : T suppρ δ } T Ωδ = T \ T 0 T Ωδ For eac element, term I can be estimated in te following way: T T 0 : Clearly φ, I Iρ δ φ K = 0. T T Ωδ : Using a standard interpolation error estimate for te Lagrange interpolant, we conclude tat φ, I Iρ δ φ T φ T ρ δ φ 2,T φ T ρ δ W 2, T φ H 1 T φ 2 T + φ 2 T T T Ωδ 5.32 T T Ωδ : In tis case ρ δ is discontinuous in T and to deal wit tis difficulty, we use Green s formula as follows φ, I Iρ δ φ T = n φ, I Iρ δ φ T n φ T I Iρ δ φ T 1/2 φ T I Iρ δ φ T 1/2 φ T ρ δ φ T ɛ 1 φ 2 T + ɛ ρ δ φ 2 T
10 for eac ɛ > 0. Here we used an inverse inequality and te interpolation estimate v I v v F F v F v F on eac of te faces F T of element T. Here, F is te tangent gradient F v = P F associated wit te face F and P F = I n F n F, were n F is te unit normal to F, te projection onto te tangent space of F. Now ρ δ φ T can be estimated by observing tat ρ δ L T since T T Ωδ. Using Hölder s inequality, we ave ρ δ φ 2 T ρ δ 2 L T φ 2 T + ρ δ 2 L T φ 2 T φ 2 T + φ 2 T 1 φ 2 T + φ 2 T + φ 2 T were we again used a trace inequality and an inverse estimate. Combining 5.33 and 5.34, we tus ave 1 φ 2 T + φ 2 T 5.34 φ, I Iρ δ φ T ɛ 1 φ 2 T + ɛ 1 φ 2 T T T Ωδ 5.35 for all 0 < ɛ 1. Summing over te elements and using 5.32 and 5.35, we obtain I φ 2 T + φ 2 T + ɛ 1 φ 2 T + ɛ 1 φ 2 T T T Ωδ ɛ 1 T T T T Ωδ φ 2 T + φ 2 T + ɛ T T Ωδ sup 0 d δ 0 φ 2 T Ω d ɛ 1 ψ 2 Ω + ɛ sup φ 2 Ω d 5.36 δ d δ 0 for all 0 < ɛ 1. Term II. An application of Green s formula gives te following identity II = φ, ρ δ φ Ωδ + κρ δ φ, φ Ωδ = ρ δ φ, φ Ωδ + φ φ, ρ δ Ωδ + κρ δ φ, φ Ωδ = ρ δ φ, φ Ωδ + κρ δ φ, φ Ωδ 1 2 φ2, ρ δ Ωδ φ2, n ρ δ Ωδ We tus obtain te estimate φ 2 ρ δ,ω + κ ρ δ φ 2 Ω φ 2 Ω δ ρ δ L Ω δ + φ 2 Ω δ ρ δ L Ω δ + I 5.37 φ 2 Ω δ + φ 2 Ω δ + ɛ 1 ψ 2 Ω + ɛ sup δ d δ 0 φ 2 Ω d 5.38 φ 2 Ω δ + ɛ 1 ψ 2 Ω + ɛ sup φ 2 Ω d δ d δ 0 for δ δ δ. Here we used te estimate 5.36 for Term I in 5.37 and te estimate 5.23 to bound φ 2 Ω δ for δ δ δ in Tus, letting δ δ we obtain φ 2 ρ δ,ω + κ φ 2 ρ δ,ω φ 2 Ω δ + ɛ 1 ψ 2 Ω + ɛ sup φ 2 Ω d 5.39 δ d δ 0 Using te fact n ρ c > 0 for δ 0 small enoug, we also obtain te bound sup φ 2 Ω δ φ 2 ρ δ δ,ω + κ φ 2 ρ δ,ω + φ 2 Ω δ + I δ δ0 φ 2 Ω δ + ɛ 1 ψ 2 Ω + ɛ sup φ 2 Ω d 5.40 δ d δ 0 10
11 were we used 5.36 and 5.39 in te second inequality. Coosing an appropriate ɛ and combining 5.39 and 5.40, we arrive at φ 2 ρ δ,ω + κ φ 2 ρ δ,ω + sup φ 2 Ω δ φ 2 Ω δ δ + ψ 2 Ω 5.41 δ δ0 To conclude te proof, we first note tat φ 2 Ω can be estimated by φ 2 Ω = φ 2 Ω\Ω d + φ 2 Ω d 5.42 d sup φ 2 Ω δ + d 1 φ 2 ρ,ω d δ d Applying te same argument for te domain Ω δ and te sifted distance function ρ δ, we note tat by coosing δ < d δ 0 small enoug and κ large enoug, te term φ 2 Ω δ can be absorbed in te left and side of Tus we finally ave te estimate φ 2 ρ δ,ω + κ φ 2 ρ δ,ω + sup φ 2 Ω δ ψ 2 Ω 5.44 δ δ δ0 Te proof now follows from combining 5.19 and We are now in te position to finalize te proof of of Proposition 5.1: Proof. Proposition 5.1 We decompose te solution φ to 4.1 into a sum φ = φ,0 + φ,1 were and Setting v = φ,0 in 5.46 we find tat a φ,1, v + κφ,1, v = m ψ, v v V 5.45 a φ,0, v = κφ,1, v v V 5.46 φ,0 2 κ φ,1 Ω φ,0 Ω ψ Ω φ, were we used Caucy-Scwarz, Poincaré, and Proposition 5.4. Tus Using tis estimate, we also obtain Collecting te estimates we conclude tat te estimate φ,0 2 Ω + 1 φ,0 2 Ω ψ 2 Ω 5.48 sup φ,0 2 Ω δ ψ 2 Ω δ δ 0 φ,0 2 Ω + sup φ,0 2 Ω δ ψ 2 Ω δ δ 0 olds. Observing tat tis estimate is stronger compared to te desired estimate and te triangle inequality, te estimate 5.50 for φ 0 and te estimate for φ,1 given by Proposition 5.4 we obtain te desired result. 6. L 2 Error Estimates for te Boundary Flux. Te previous results on te weigted stability estimate and te anisotropic interpolation error enable us to prove te main result of our work: Teorem 6.1. Let Σ n be te discrete boundary flux defined by eiter 3.9 or 3.16 and suppose u satisfies te assumption of Proposition 5.2. Ten te following error estimate olds σ n Σ n Ω ln 6.1 Proof. We start wit te proof of te estimate for te Nitsce flux 3.9. Recalling te estimate 4.10 in Lemma 4.1, we need to estimate I = a u π u, φ m ψ, u π u sup, II = sup 6.2 ψ L 2 Ω\{0} ψ Ω ψ L 2 Ω\{0} ψ Ω 11
12 δ Ω Ω δ Ω δ Ω δ T 0 δ Ω δ δ 0 Ω δ0 T Ωδ T Ωδ Fig Left Decomposition of te domain Ω into two boundary layers of widt and te far-field domain Ω δ0. Rigt Decomposition of te mes T wit respect to Ω δ consisting of a -band T Ωδ blue, an inner mes T Ωδ gray and a boundary zone T 0 wite. Estimate of I. We ave a u π u, φ = u π u, φ Ω n u π u, φ Ω 6.3 u π u, n φ Ω + β 1 u π u, φ Ω Te boundary terms may be directly estimated using a trace inequality followed by te interpolation error estimate 5.18 and te stability estimate 5.6. For instance, n u π u, φ Ω n u π u Ω φ Ω ψ Ω 6.4 u π u, n φ Ω 1 u π u Ω n φ Ω ψ Ω 6.5 and te oter terms may be estimated in te same way. To estimate te interior term we first split te integral as follows u π u, φ Ω = u π u, φ Ω\Ωδ u π u, φ Ωδ \Ω δ0 + u π u, φ Ωδ0 = I 1 + I 2 + I Term I 1. Using Caucy-Scwarz in te tangent direction and Hölders inequality in te normal direction we ave I = u π u, φ Ω\Ωδ 6.8 δ sup u π u Ωs φ Ωs ds s δ 0 1/2 δ 1/2 δ sup u π u Ωs ds φ 2 Ω s ds s δ 0 0 sup u π u Ωs δ 1/2 φ Ω\Ωδ s δ Since δ = C we may employ Proposition 5.1 as follows δ φ 2 Ω\Ω δ φ 2 Ω\Ω δ φ 2 Ω ψ 2 Ω
13 Using te interpolation error estimate 5.18 we get te estimate we get I ψ Ω 6.13 Term I 2. Proceeding in te same way as for Term I and observing tat ρ δ x = s δ, x Ω s 6.14 II = u π u, φ Ωδ \Ω δ 6.15 δ0 sup u π u Ωs φ Ωs ds 6.16 δ s δ 0 δ sup u π u Ωs 6.17 δ s δ 0 1/2 δ0 1/2 δ0 s δ 1 ds s δ φ 2 Ω s ds δ δ ln δ 1/2 φ ρδ,ω\ω δ 6.18 ln δ 1/2 ψ Ω 6.19 were we used te interpolation error estimate 5.18 and te stability estimate in Proposition 5.1. Term I 3. Using Caucy-Scwarz we obtain I = u π u, φ Ω\Ωδ 6.20 u π u Ω\Ωδ δ 1/2 φ ρ,ω\ωδ 6.21 wic can be directly estimated using standard interpolation error estimates and te stability bound. Estimate of II. Using Caucy-Scwarz and te interpolation estimate 5.18 we obtain II β 1 u π u Ω ψ Ω ψ Ω 6.22 wic concludes te proof. Following te same line of reasoning, we now state and prove te corresponding L 2 -error estimate wen te boundary flux is approximated by te Lagrange multiplier, cf Referring to te variational problem 3.13, te stabilization form is supposed te following localized version of te continuity condition 3.18 u Ωδ + 1/2 n u Ω + 1/2 u Ω + 1/2 λ Ω c u, λ; v, µ v Ωδ + 1/2 n v Ω + 1/2 v Ω + 1/2 µ Ω 6.23 Tis assumptions is trivially satisfies by te stabilization form 3.20 and merely quantifies tat te region of influence of te stabilization is located on or close to te boundary. Teorem 6.2. Let Σ n be te discrete boundary flux defined 3.16 and assume te u satisfies te assumption of Proposition 5.2 and tat λ H 1 Ω. Ten te following error estimate olds σ n Σ n Ω ln 6.24 Proof. Starting from te error representation formula 4.13, we need to estimate I = A π u u, π λ λ; φ, θ m ψ, λ π λ sup, II = sup, ψ L 2 Ω\0 ψ Ω ψ L 2 Ω\0 ψ Ω 13
14 By definition, A π u u, π λ λ; φ, θ = π u u, v Ω + π λ λ; φ Ω + θ, π u u Ω c π λ λ, π u u; φ, θ Since te estimate for first term as already been derived in te previous proof, it remains to bound te contribution from te boundary terms and te stabilization form. An application of te interpolation estimates and te discrete energy stability 5.20 yields π λ λ; φ Ω λ 1, Ω ψ Ω θ, π u u Ω θ Ω 1 π u u Ω ψ Ω Because of assumption 6.23, te contribution from te stabilization form can be estimated similarly. Finally, tanks to an interpolation estimate, term II trivially satisfies II λ Ω. 7. Numerical Results. We consider te elliptic model problem on te domain Ω = [0, 1] [0, 1] R 2. To examine te convergence rate of te normal flux approximations, we employ te metod of manufactured solution and coose ux, y = cos2πx cos2πy + sin2πx sin2πy as a reference solution, g = u Ω and f = u as te corresponding boundary data and source function, respectively. As discretization scemes, we pick Nitsce s metod 3.1 and a stabilized Lagrange multiplier metod 3.13 wit te stabilization form given by For te stabilization parameters we take α = β = 10. Te approximations for te boundary flux are ten computed on a sequences of uniform meses {T } wit mes sizes k for k = 0,..., 15. Te numerical results are depicted in Figure 7.1. In te pre-asymptotic regime ranging from 0.35 to 0.1, te convergence rate of bot metods deviates significantly from te optimal slope 1.0. Consequently, te fitted slopes indicate a sligtly sub-optimal convergence rate for te Nitsce flux, wile te convergence rate for Lagrange multiplier metod is iger ten te teoretical prediction. If we discard te pre-asymptotic regime as sown in te rigt plot of Figure 7.1, te approximation error σ n Σ Ω exibits optimal convergence rate for bot metods and corroborates te teoretical findings of our work. 14
15 σ n Σ n 0, Ω max A: 0.96 B: σ n Σ n 0, Ω A: 0.99 B: max Fig L 2 Ω convergence study for various flux computations. A Nitsce flux for CG1 elements. B Lagrange multiplier computed wit te stabilized metod by Barbosa and Huges based on a CG1 DG0 discretization. Te legend gives te fitted slope for eac approximation error. Left Approximation error for te entire mes sequence revealing different beavior in te pre-asymptotic regime. Rigt Asymptotic regime. Starting from 0.1 bot metods give optimal first order convergence. 15
16 Acknowledgments. Tis work is supported by a Center of Excellence grant from te Researc Council of Norway to te Center for Biomedical Computing at Simula Researc Laboratory. References. [1] I. Babuška. Te finite element metod wit Lagrangian multipliers. Num. Mat., 203: , June [2] H.J.C. Barbosa and T.J.R. Huges. Te finite element metod wit Lagrange multipliers on te boundary: circumventing te Babuška-Brezzi condition. Computer Metods in Applied Mecanics and Engineering, 851: , [3] H.J.C. Barbosa and T.J.R. Huges. Boundary Lagrange multipliers in finite element metods: error analysis in natural norms. Numer. Mat., 621:1 15, [4] J.H. Bramble. Te Lagrange multiplier metod for Diriclet s problem. Mat. Comput., :1 11, [5] F. Brezzi. On te Existence, Uniqueness and Approximation of Saddle-Point Problems Arising from Lagrangian Multipliers. RAIRO Anal. Numér., R 2: , [6] Erik Burman. Projection stabilization of Lagrange multipliers for te imposition of constraints on interfaces and boundaries. Numerical Metods for Partial Differential Equations, pages n/a n/a, [7] G.F. Carey, S.S. Cow, and M.K. Seager. Approximate boundary-flux calculations. Computer Metods in Applied Mecanics and Engineering, 502: , [8] J. Cabrowski. Note on te Diriclet problem wit L2-boundary data. Manuscripta Mat., 10840:91 108, [9] J. Cabrowski. Te Diriclet Problem wit L2-Boundary Data for Elliptic Linear Equations, volume 1482 of Lecture Notes in Matematics. Springer, [10] D. Estep, S. Tavener, and T. Wildey. A posteriori error estimation and adaptive mes refinement for a multiscale operator decomposition approac to fluid solid eat transfer. Journal of Computational Pysics, 22911: , [11] D. Gilbard and N.S. Trudinger. Elliptic partial differential equations of second order, volume 224 of Classics in Matematics. Springer, [12] M. Giles, M.G. Larson, M. Levenstam, and E. Süli. Adaptive error control for finite element approximations of te lift and drag coefficients in viscous flow. Tecnical report, Te Matematical Institute, University of Oxford, [13] A. Hansbo and P. Hansbo. An unfitted finite element metod, based on Nitsce s metod, for elliptic interface problems. Comput. Metods Appl. Mec. Engrg., : , [14] P. Hansbo. Nitsce s metod for interface problems in computational mecanics. GAMM- Mitt, 282: , [15] J.M. Melenk and B. Wolmut. Quasi-optimal approximation of surface based Lagrange multipliers in finite element metods. SIAM J. Numer. Anal., 504: , [16] J. Nitsce. Über ein Variationsprinzip zur Lösung von Diriclet-Problemen bei Verwendung von Teilräumen, die keinen Randbedingungen unterworfen sind. Abandlungen aus dem Matematiscen Seminar der Universität Hamburg, 361:9 15, July [17] A.I. Pelivanov, R.D. Lazarov, G.F. Carey, and S.S. Cow. Superconvergence analysis of approximate boundary-flux calculations. Numerisce Matematik, 631: , [18] J. Pitkäranta. Boundary subspaces for te finite element metod wit Lagrange multipliers. Numer. Mat., 28933: , [19] J. Pitkäranta. Local stability conditions for te Babuška metod of Lagrange multipliers. Mat. Comp., 35152: , [20] R. Rannacer and R. Scott. Some optimal error estimates for piecewise linear finite element approximations. Mat. Comp., 38158: , [21] R. Scott and S. Zang. Finite element interpolation of nonsmoot functions satisfying boundary conditions. Mat. Comp., 54190: , [22] R. Stenberg. On some tecniques for approximating boundary conditions in te finite element metod. J. Comput. Appl. Mat., 631: ,
Different Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More information2009 Elsevier Science. Reprinted with permission from Elsevier.
P. Hansbo and M. Juntunen. 009. Weakly imposed Diriclet boundary conditions for te Brinkman model of porous media flow. Applied Numerical Matematics, volume 59, number 6, pages 174 189. doi:10.1016/j.apnum.008.07.003.
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationHigh-Order Extended Finite Element Methods for Solving Interface Problems
Hig-Order Extended Finite Element Metods for Solving Interface Problems arxiv:1604.06171v1 [mat.na] 21 Apr 2016 Fei Wang Yuanming Xiao Jincao Xu ey words. Elliptic interface problems, unfitted mes, extended
More informationarxiv: v1 [math.na] 11 May 2018
Nitsce s metod for unilateral contact problems arxiv:1805.04283v1 [mat.na] 11 May 2018 Tom Gustafsson, Rolf Stenberg and Jua Videman May 14, 2018 Abstract We derive optimal a priori and a posteriori error
More informationVariational Localizations of the Dual Weighted Residual Estimator
Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)
More informationarxiv: v1 [math.na] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More informationarxiv: v1 [math.na] 19 Mar 2018
A primal discontinuous Galerkin metod wit static condensation on very general meses arxiv:1803.06846v1 [mat.na] 19 Mar 018 Alexei Lozinski Laboratoire de Matématiques de Besançon, CNRS UMR 663, Univ. Bourgogne
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 informationAN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS
AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous
More informationarxiv: v1 [math.na] 9 Sep 2015
arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on
More informationA SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALSIS AND MODELING Volume 3, Number 4, Pages 6 626 c 26 Institute for Scientific Computing and Information A SPLITTING LEAST-SQUARES MIED FINITE ELEMENT METHOD FOR
More informationarxiv: v1 [math.na] 20 Jul 2009
STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN DOUGLAS N. ARNOLD AND MARIE E. ROGNES arxiv:0907.3438v1 [mat.na] 20 Jul 2009 Abstract. Te stability properties of simple element coices for te mixed
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationA FINITE ELEMENT METHOD ON COMPOSITE GRIDS BASED ON NITSCHE S METHOD
ESAIM: MAN Vol. 37, N o 3, 3, pp. 495 54 DOI:.5/man:339 ESAIM: Matematical Modelling and Numerical Analysis A FINITE ELEMENT METHOD ON COMPOSITE GRIDS BASED ON NITSCHE S METHOD Anita Hansbo, Peter Hansbo
More informationMIXED AND STABILIZED FINITE ELEMENT METHODS FOR THE OBSTACLE PROBLEM
MIXED AND STABILIZED FINITE ELEMENT METHODS FOR THE OBSTACLE PROBLEM TOM GUSTAFSSON, ROLF STENBERG, AND JUHA VIDEMAN Abstract. We discretize te Lagrange multiplier formulation of te obstacle problem by
More informationKey words. Sixth order problem, higher order partial differential equations, biharmonic problem, mixed finite elements, error estimates.
A MIXED FINITE ELEMENT METHOD FOR A SIXTH ORDER ELLIPTIC PROBLEM JÉRÔME DRONIOU, MUHAMMAD ILYAS, BISHNU P. LAMICHHANE, AND GLEN E. WHEELER Abstract. We consider a saddle point formulation for a sixt order
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 informationA UNIFORM INF SUP CONDITION WITH APPLICATIONS TO PRECONDITIONING
A UNIFORM INF SUP CONDIION WIH APPLICAIONS O PRECONDIIONING KEN ANDRE MARDAL, JOACHIM SCHÖBERL, AND RAGNAR WINHER Abstract. A uniform inf sup condition related to a parameter dependent Stokes problem is
More informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationarxiv: v2 [math.na] 5 Jul 2017
Trace Finite Element Metods for PDEs on Surfaces Maxim A. Olsanskii and Arnold Reusken arxiv:1612.00054v2 [mat.na] 5 Jul 2017 Abstract In tis paper we consider a class of unfitted finite element metods
More informationAnalysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem
Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for
More informationA SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION
A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION GABRIEL R. BARRENECHEA, LEOPOLDO P. FRANCA 1 2, AND FRÉDÉRIC VALENTIN Abstract. Tis work introduces and analyzes novel stable
More informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationA posteriori error analysis of multiscale operator decomposition methods for multiphysics models
A posteriori error analysis of multiscale operator decomposition metods for multipysics models D. Estep, V. Carey 2, V. Ginting 3, S. Tavener 2, T. Wildey 4 5 Department of Matematics and Department of
More informationA symmetric mixed finite element method for nearly incompressible elasticity based on. biorthogonal system
A symmetric mixed finite element metod for nearly incompressible elasticity based on biortogonal systems Bisnu P. Lamicane and Ernst P. Stepan February 18, 2011 Abstract We present a symmetric version
More informationERROR ESTIMATES FOR THE DISCONTINUOUS GALERKIN METHODS FOR PARABOLIC EQUATIONS. 1. Introduction. We consider the parabolic PDE of the form,
ERROR ESTIMATES FOR THE DISCONTINUOUS GALERKIN METHODS FOR PARABOLIC EQUATIONS K. CHRYSAFINOS AND NOEL. J. WALKINGTON Abstract. We analyze te classical discontinuous Galerkin metod for a general parabolic
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
More informationError estimates for a semi-implicit fully discrete finite element scheme for the mean curvature flow of graphs
Interfaces and Free Boundaries 2, 2000 34 359 Error estimates for a semi-implicit fully discrete finite element sceme for te mean curvature flow of graps KLAUS DECKELNICK Scool of Matematical Sciences,
More informationAn approximation method using approximate approximations
Applicable Analysis: An International Journal Vol. 00, No. 00, September 2005, 1 13 An approximation metod using approximate approximations FRANK MÜLLER and WERNER VARNHORN, University of Kassel, Germany,
More informationA SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS
A SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS CONSTANTIN BACUTA AND KLAJDI QIRKO Abstract. We investigate new PDE discretization approaces for solving variational formulations wit different types
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationSome Error Estimates for the Finite Volume Element Method for a Parabolic Problem
Computational Metods in Applied Matematics Vol. 13 (213), No. 3, pp. 251 279 c 213 Institute of Matematics, NAS of Belarus Doi: 1.1515/cmam-212-6 Some Error Estimates for te Finite Volume Element Metod
More informationHigher order unfitted isoparametric space-time FEM on moving domains
Higer order unfitted isoparametric space-time FEM on moving domains Master s tesis by: Janosc Preuß Supervisor: Jun.-Prof. Dr. Cristop Lerenfeld Second Assessor: Prof. Dr. Gert Lube prepared at te Institute
More informationWeierstraß-Institut. im Forschungsverbund Berlin e.v. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stocastik im Forscungsverbund Berlin e.v. Preprint ISSN 0946 8633 Stability of infinite dimensional control problems wit pointwise state constraints Micael
More informationAPPROXIMATION BY QUADRILATERAL FINITE ELEMENTS
MATHEMATICS OF COMPUTATION Volume 71, Number 239, Pages 909 922 S 0025-5718(02)01439-4 Article electronically publised on Marc 22, 2002 APPROXIMATION BY QUADRILATERAL FINITE ELEMENTS DOUGLAS N. ARNOLD,
More informationA Weak Galerkin Method with an Over-Relaxed Stabilization for Low Regularity Elliptic Problems
J Sci Comput (07 7:95 8 DOI 0.007/s095-06-096-4 A Weak Galerkin Metod wit an Over-Relaxed Stabilization for Low Regularity Elliptic Problems Lunji Song, Kaifang Liu San Zao Received: April 06 / Revised:
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J NUMER ANAL Vol 4, No, pp 86 84 c 004 Society for Industrial and Applied Matematics LEAST-SQUARES METHODS FOR LINEAR ELASTICITY ZHIQIANG CAI AND GERHARD STARKE Abstract Tis paper develops least-squares
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More information1. Introduction. Consider a semilinear parabolic equation in the form
A POSTERIORI ERROR ESTIMATION FOR PARABOLIC PROBLEMS USING ELLIPTIC RECONSTRUCTIONS. I: BACKWARD-EULER AND CRANK-NICOLSON METHODS NATALIA KOPTEVA AND TORSTEN LINSS Abstract. A semilinear second-order parabolic
More informationarxiv: v1 [math.na] 6 Dec 2010
MULTILEVEL PRECONDITIONERS FOR DISCONTINUOUS GALERKIN APPROXIMATIONS OF ELLIPTIC PROBLEMS WITH JUMP COEFFICIENTS BLANCA AYUSO DE DIOS, MICHAEL HOLST, YUNRONG ZHU, AND LUDMIL ZIKATANOV arxiv:1012.1287v1
More informationEffect of Numerical Integration on Meshless Methods
Effect of Numerical Integration on Mesless Metods Ivo Babuška Uday Banerjee Jon E. Osborn Qingui Zang Abstract In tis paper, we present te effect of numerical integration on mesless metods wit sape functions
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationMixed Finite Element Methods for Incompressible Flow: Stationary Stokes Equations
Mixed Finite Element Metods for Incompressible Flow: Stationary Stoes Equations Ziqiang Cai, Carles Tong, 2 Panayot S. Vassilevsi, 2 Cunbo Wang Department of Matematics, Purdue University, West Lafayette,
More informationComputing eigenvalues and eigenfunctions of Schrödinger equations using a model reduction approach
Computing eigenvalues and eigenfunctions of Scrödinger equations using a model reduction approac Suangping Li 1, Ziwen Zang 2 1 Program in Applied and Computational Matematics, Princeton University, New
More informationA Finite Element Primer
A Finite Element Primer David J. Silvester Scool of Matematics, University of Mancester d.silvester@mancester.ac.uk. Version.3 updated 4 October Contents A Model Diffusion Problem.................... x.
More informationA MULTILEVEL PRECONDITIONER FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD
A MULTILEVEL PRECONDITIONER FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD KOLJA BRIX, MARTIN CAMPOS PINTO, AND WOLFGANG DAHMEN Abstract. In tis article we present a multilevel preconditioner for
More informationarxiv: v3 [math.na] 25 Jun 2017
ANALYSIS OF A HIGH ORDER UNFITTED FINITE ELEMENT METHOD FOR ELLIPTIC INTERFACE PROBLEMS CHRISTOPH LEHRENFELD AND ARNOLD REUSKEN arxiv:1602.02970v3 [mat.na] 25 Jun 2017 Abstract. In te context of unfitted
More informationWEIGHTED ERROR ESTIMATES OF THE CONTINUOUS INTERIOR PENALTY METHOD FOR SINGULARLY PERTURBED PROBLEMS
WEIGHTED ERROR ESTIMATES OF THE CONTINUOUS INTERIOR PENALTY METHOD FOR SINGULARLY PERTURBED PROBLEMS ERIK BURMAN, JOHNNY GUZMÁN, AND DMITRIY LEYKEKHMAN Abstract. In tis paper we analyze local properties
More informationComputers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations
Computers Matematics wit Applications 59 () 5 4 Contents lists available at ScienceDirect Computers Matematics wit Applications journal omepage: www.elsevier.com/locate/camwa A nonlinear weigted least-squares
More informationCONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL SINGULARLY PERTURBED DIFFERENTIAL EQUATIONS ON EQUIDISTRIBUTING MESHES
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Number, Pages 57 74 c 5 Institute for Scientific Computing and Information CONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL
More informationIsoparametric finite element approximation of curvature on hypersurfaces
Isoparametric finite element approximation of curvature on ypersurfaces C.-J. Heine Abstract Te discretisation of geometric problems often involves a point-wise approximation of curvature, altoug te discretised
More informationOn convergence of the immersed boundary method for elliptic interface problems
On convergence of te immersed boundary metod for elliptic interface problems Zilin Li January 26, 2012 Abstract Peskin s Immersed Boundary (IB) metod is one of te most popular numerical metods for many
More informationCELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION SCHEME WITHOUT FICTITIOUS DOMAINS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 7, Number 1, Pages 1 9 c 010 Institute for Scientific Computing and Information CELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More informationDownloaded 11/15/17 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 55, No. 6, pp. 2787 2810 c 2017 Society for Industrial and Applied Matematics EDGE ELEMENT METHOD FOR OPTIMAL CONTROL OF STATIONARY MAXWELL SYSTEM WITH GAUSS LAW IRWIN YOUSEPT
More informationChemnitz Scientific Computing Preprints
G. Of G. J. Rodin O. Steinbac M. Taus Coupling Metods for Interior Penalty Discontinuous Galerkin Finite Element Metods and Boundary Element Metods CSC/11-02 Cemnitz Scientific Computing Preprints Impressum:
More informationA trace finite element method for a class of coupled bulk-interface transport problems
Numerical Analysis and Scientific Computing Preprint Seria A trace finite element metod for a class of coupled bulk-interface transport problems S. Gross M.A. Olsanskii A. Reusken Preprint #28 Department
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationChapter 10. Function approximation Function approximation. The Lebesgue space L 2 (I)
Capter 1 Function approximation We ave studied metods for computing solutions to algebraic equations in te form of real numbers or finite dimensional vectors of real numbers. In contrast, solutions to
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationA New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions
Numerisce Matematik manuscript No. will be inserted by te editor A New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions Wang Ming 1, Zong-ci Si 2, Jincao Xu1,3 1 LMAM, Scool of Matematical
More informationA posteriori error analysis for time-dependent Ginzburg-Landau type equations
A posteriori error analysis for time-dependent Ginzburg-Landau type equations Sören Bartels Department of Matematics, University of Maryland, College Park, MD 74, USA Abstract. Tis work presents an a posteriori
More informationERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS
ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS YONGYONG CAI, AND JIE SHEN Abstract. We carry out in tis paper a rigorous error analysis
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationOSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix
Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra
More informationDiscontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System
Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics
More informationAPPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS. Introduction
Acta Mat. Univ. Comenianae Vol. LXVII, 1(1998), pp. 57 68 57 APPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS A. SCHMIDT Abstract. Te pase transition between solid and liquid in an
More informationA posteriori error estimates for non-linear parabolic equations
A posteriori error estimates for non-linear parabolic equations R Verfürt Fakultät für Matematik, Rur-Universität Bocum, D-4478 Bocum, Germany E-mail address: rv@num1rur-uni-bocumde Date: December 24 Summary:
More informationAdaptive Finite Element Method
39 Capter 3 Adaptive inite Element Metod 31 Introduction As already pointed out in capter 2, singularities occur in interface problems Wen discretizing te problem (221) wit inite Elements, te singularities
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationDownloaded 06/08/17 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 55, No. 3, pp. 357 386 c 7 Society for Industrial and Applied Matematics Downloaded 6/8/7 to 67.96.45.. Redistribution subject to SIAM license or copyrigt; see ttp://www.siam.org/journals/ojsa.pp
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationA Local Projection Stabilization/Continuous Galerkin Petrov Method for Incompressible Flow Problems
A Local Projection Stabilization/Continuous Galerkin Petrov Metod for Incompressible Flow Problems Naveed Amed, Volker Jon, Gunar Matties, Julia Novo Marc 22, 218 Abstract A local projection stabilization
More informationA Demonstration of the Advantage of Asymptotic Preserving Schemes over Standard Finite Volume Schemes
A Demonstration of te Advantage of Asymptotic Preserving Scemes over Standard Finite Volume Scemes Jocen Scütz Berict Nr. 366 Juni 213 Key words: conservation laws, asymptotic metods, finite volume metods,
More informationUniform estimate of the constant in the strengthened CBS inequality for anisotropic non-conforming FEM systems
Uniform estimate of te constant in te strengtened CBS inequality for anisotropic non-conforming FEM systems R. Blaeta S. Margenov M. Neytceva Version of November 0, 00 Abstract Preconditioners based on
More informationNumerical Analysis of the Double Porosity Consolidation Model
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real 21-25 septiembre 2009 (pp. 1 8) Numerical Analysis of te Double Porosity Consolidation Model N. Boal
More informationEXTENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS WITH APPLICATION TO AN AEROACOUSTIC PROBLEM
SIAM J. SCI. COMPUT. Vol. 26, No. 3, pp. 821 843 c 2005 Society for Industrial and Applied Matematics ETENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS
More informationA SHORT INTRODUCTION TO BANACH LATTICES AND
CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,
More informationDecay of solutions of wave equations with memory
Proceedings of te 14t International Conference on Computational and Matematical Metods in Science and Engineering, CMMSE 14 3 7July, 14. Decay of solutions of wave equations wit memory J. A. Ferreira 1,
More informationContinuity. Example 1
Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)
More informationPosteriori Analysis of a Finite Element Discretization for a Penalized Naghdi Shell
International Journal of Difference Equations ISSN 0973-6069, Volume 8, Number 1, pp. 111 124 (2013) ttp://campus.mst.edu/ijde Posteriori Analysis of a Finite Element Discretization for a Penalized Nagdi
More informationAnisotropic, adaptive finite elements for the computation of a solutal dendrite
Interfaces and Free Boundaries 5 (23), 13 127 Anisotropic, adaptive finite elements for te computation of a solutal dendrite E. BURMAN AND M. PICASSO Institut de Matématiques, École Polytecnique Fédérale
More informationParallel algorithm for convection diffusion system based on least squares procedure
Zang et al. SpringerPlus (26 5:69 DOI.86/s464-6-3333-8 RESEARCH Parallel algoritm for convection diffusion system based on least squares procedure Open Access Jiansong Zang *, Hui Guo 2, Hongfei Fu 3 and
More informationERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING
ERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING JINGYUE WANG AND BRADLEY J. LUCIER Abstract. We bound te difference between te solution to te continuous Rudin Oser Fatemi
More informationarxiv: v1 [math.na] 25 Jul 2014
A second order in time, uniquely solvable, unconditionally stable numerical sceme for Can-Hilliard-Navier-Stokes equation Daozi Han, Xiaoming Wang November 5, 016 arxiv:1407.7048v1 [mat.na] 5 Jul 014 Abstract
More informationAdvancements In Finite Element Methods For Newtonian And Non-Newtonian Flows
Clemson University TigerPrints All Dissertations Dissertations 8-3 Advancements In Finite Element Metods For Newtonian And Non-Newtonian Flows Keit Galvin Clemson University, kjgalvi@clemson.edu Follow
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationRobust approximation error estimates and multigrid solvers for isogeometric multi-patch discretizations
www.oeaw.ac.at Robust approximation error estimates and multigrid solvers for isogeometric multi-patc discretizations S. Takacs RICAM-Report 2017-32 www.ricam.oeaw.ac.at Robust approximation error estimates
More information