Technische Universität Graz
|
|
- Gregory Nichols
- 6 years ago
- Views:
Transcription
1 Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2
2
3 Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2
4 Technische Universität Graz Institut für Numerische Mathematik Steyrergasse 30 A 8010 Graz WWW: c Alle Rechte vorbehalten. Nachdruck nur mit Genehmigung des Autors.
5 Stability of the Laplace single layer boundary integral operator in Sobolev spaces Olaf Steinbach Institut für Numerische Mathematik, TU Graz, Steyrergasse 30, 8010 Graz, Austria 1 Introduction As a model problem we consider the Dirichlet boundary value problem for the Laplace equation, u(x) = 0 for x, u(x) = g(x) for x =, (1.1) where R 3 is a bounded Lipschitz polyhedron. Using an indirect approach, the solution of (1.1) can be described as single layer potential u(x) = (Ṽw)(x) := U (x,y)w(y)ds y for x, (1.2) where U (x,y) = 1 1 4π x y is the fundamental solution of the Laplacian. It is well known, see, e.g. [2], that Ṽ : H 1/2 () H 1 (). The unknown density w H 1/2 () is then found by applying the interior Dirichlet trace operator γ int 0 : H 1 () H 1/2 () to (1.2) which results in the boundary integral equation (Vw)(x) := U (x,y)w(y)ds y = g(x) for x, (1.3) and which is equivalent to a Galerkin Bubnov formulation: Find w H 1/2 () such that Vw,v = g,v for all v H 1/2 (). (1.4) Since the single layer boundary integral operator V = γ int 0 Ṽ : H 1/2 () H 1/2 () is elliptic, [6], Vw,w c V 1 w 2 H 1/2 () for all w H 1/2 (), (1.5) 1
6 unique solvability of the variational formulation (1.4) follows. Moreover we can deduce a stability and error analysis of related boundary element discretization schemes, see, e.g., [7]. Error estimates then rely on the regularity of w = V 1 g, i.e. on the regularity of the given Dirichlet datum g, and on the mapping properties of the single layer boundary integral operator V. In the case of a Lipschitz domain we have that V : H 1/2+s () H 1/2+s () is bijective for all s [ 1, 1 ], see [2, 8], while in the case of a polyhedral bounded domain 2 2 this remains true for s < s 0 where s 0 > 1 is determined by the related interior and 2 exterior angles in corners and at edges, see, e.g., [5], and [4] for the two dimensional case. The error estimate for the Galerkin solution of the Galerkin Bubnov variational formulation (1.4) is given, due to Cea s lemma, in the energy norm in H 1/2 (). Hence, to derive error estimates in stronger norms, e.g. in L 2 (), we have to use an inverse inequality for the used boundary element space and where we have to assume a globally quasi uniform boundary element mesh, see, e.g., [7] for a more detailed discussion. In fact, this excludes non uniform and adaptive meshes as often used in practice. Instead of the Galerkin Bubnov variational formulation(1.4) we will consider a Galerkin Petrov variational formulation which allows the use of different trial and test spaces, both in the continuous and discrete setting. In this case, the ellipticity estimate (1.5) has to be replaced by an appropriate stability condition, also known as inf sup condition. While the analysis of the Galerkin Bubnov formulation (1.4) relies on a related domain variational formulation in H 1 (), our analysis is based on using a Galerkin Petrov domain variational formulation for which we have to introduce suitable Sobolev spaces. With this we can not only conclude known mapping properties of the single layer boundary integral operator, but we can establish a new stability condition which ensures unique solvability of the Galerkin Petrov variational formulation. In this note we will not consider a stability and error analysis of related Galerkin Petrov boundary element methods which will be a topic of further research. In fact, such an approach can also be used for Galerkin Petrov variational formulations in weaker Sobolev spaces, e.g., when the given Dirichlet data have reduced regularity, for example if we have g L 2 () only, see, e.g., [1]. However, the main focus of future work will be on the extension of this concept to the mathematical and numerical analysis of boundary integral equation and boundary element methods for time dependent problems such as the heat equation, see, e.g., [3] for related Galerkin Bubnov formulations. 2 Strong domain variational formulation For the Dirichlet boundary value problem (1.1) we consider, instead of a standard domain variationalformulationinh 1 ()whichisbasedongreen sfirstformula, agalerkin Petrov variational formulation. For this we introduce { } H () := v H 1 () : v L 2 () H 1 (), with the norm v 2 H () = v 2 L 2 () + v 2 L 2 () + v 2 L 2 (). 2
7 Then we have to find u H () satisfying u(x) = g(x) for x such that [ u(x)]v(x)dx = 0 for all v L 2 (), (2.1) where we have to assume that the given Dirichlet datum g is in the Dirichlet trace space γ int 0 H () H 1/2 (). In particular, let u g H () be a bounded and norm preserving extension of g γ int 0 H () with g γ int 0 H () = min v H ():v =g v H () = u g H (). (2.2) It remains to find u 0 X S := H () H0 1 () such that a S (u 0,v) := [ u 0 (x)]v(x)dx = [ u g (x)]v(x)dx for all v Y S := L 2 (). (2.3) Related to the trial and test spaces we introduce the associated norms u XS := [ u 2 L 2 () + u 2 L 2 ()] 1/2, v YS := v L 2 (). Lemma 2.1 The bilinear form of the variational problem (2.3), is bounded, i.e. and satisfies the stability condition a S (u,v) u XS v YS for all u X S, v Y S, c S u XS sup 0 v Y S a S (u,v) v YS for all u X S, c S = λ min () 1+λ min () where λ min () is the minimal eigenvalue of the Dirichlet eigenvalue problem u(x) = λu(x) for x, u(x) = 0 for x. Proof. The boundedness of the bilinear form a S (, ) is a direct consequence of the Cauchy Schwarz inequality, a S (u,v) = [ u(x)]v(x) dx u L 2 () v L 2 () u XS v YS. To prove the stability condition we consider u X S and choose v = u u Y S. By using the minimal Dirichlet eigenvalue for the Laplacian in, v 2 λ min () = min L2 (), 0 v H0 1() v 2 L 2 () 3
8 and Hölders inequality we have v YS = u u L 2 () u L 2 () + u L 2 () 1 λmin () u L 2 () + u L 2 () ( 1 ) 1/2( ) 1/2 λ min () +1 u 2 1+λ min () L 2 () + u 2 L 2 () = u XS. λ min () Then, a S (u,v) = a S (u,u u) = [ u(x)][u(x) u(x)]dx = [ u(x)]u(x) dx + [ u(x)] 2 dx = u(x) 2 dx+ [ u(x)] 2 dx = u 2 L 2 () + u 2 L 2 () = u 2 X S λ min () 1+λ min () u X S v YS implies the stability condition as claimed. As a consequence of Lemma 2.1 we conclude unique solvability of the variational problem (2.3) to obtain u = u 0 + u g H (). In particular, when chosing in (2.1) v = u L 2 (), this gives u 2 L 2 () = [ u(x)] 2 dx = 0. (2.4) For the solution u H () H 1 () of the variational formulation (2.1) we note that the interior Neumann trace γ int 1 u(x) := lim x x n x x u( x) = n x u(x) for x (2.5) is well defined, at least we have γ int 1 u H 1/2 () due to duality arguments and the use of Green s first formula. To do a more detailed analysis, for u H () we define ψ = u L 2 () and we consider the Dirichlet boundary value problem φ(x) = ψ(x) for x, φ(x) = 0 for x. In the case of a domain with a sufficient smooth boundary or in the case of a convex polyhedron we find φ H 2 (), and therefore H () = H 2 () follows. However, this is not true when the domain is polyhedral bounded, but non convex. In this case, H () includes harmonic functions which are not in H 2 () but in H s (), s < s i, for some s i > 3, see, for example, [5, Corollary 2.6.7]. In any case, the Neumann trace operator 2 4
9 γ int 1 : H () γ int 1 H () is well defined, implying the Neumann trace space γ int 1 H (), and satisfying γ int 1 v γ int 1 H () c N v H () for all v H (). (2.6) Lemma 2.2 Let u H () be the unique solution of the variational formulation (2.1). Then, c i γ int 1 u 2 γ1 inth () u 2 L 2 (). (2.7) Proof. Let u = u 0 +u g H () be the unique solution of (2.3). We then define ũ(x) = u(x) γ, γ = 1 g(x)ds x, ũ(x)ds x = 0, and where ũ is the unique solution of the Dirichlet boundary value problem ũ(x) = 0 for x, ũ(x) = g(x) γ for x. Obviously, ũ H (), and (2.6) together with (2.4) then implies 1 c N γ int 1 ũ 2 γ int 1 H () ũ 2 H () = ũ 2 L 2 () + ũ 2 L 2 () + ũ 2 L 2 () = ũ 2 L 2 () + ũ 2 L 2 () = ũ 2 H 1 (). Since an equivalent norm in H 1 () is given by [ ] 2 v 2 H 1() := v(x)ds x + v 2 L 2 (), we immediately conclude Inserting ũ = u γ concludes the proof. 1 c N γ int 1 ũ 2 γ 1 H () c ũ 2 H 1 () = c ũ 2 L 2 (). In addition to the interior Dirichlet boundary value problem (1.1) we also consider the exterior Dirichlet problem u(x) = 0 for x c := R n \, u(x) = g(x) for x, (2.8) where in addition we have to impose a suitable radiation condition, u(x) = O(1/ x ) as x. (2.9) When introducing the bounded domain r := B r \ with B r := {x R 3 : x < r}, and chosing r > 0 such that B r, we can proceed as in the case of the interior Dirichlet boundary value problem (1.1), when considering the limit r and the radiation condition (2.9). As in (2.7) we may define the exterior Neumann trace of the solution u satisfying c e γ ext 1 u 2 γ1 exth ( c ) u 2 L 2 ( c ). (2.10) Note that H ( c ) H s ( c ), s < s e, for some s e > 3 which may differ from s 2 i. 5
10 3 Ultra week domain variational formulation To derive mapping properties of the single layer boundary integral operator V we may also consider the ultra week domain variational formulation, see, e.g., [1]. Multiplying the partial differential equation in (1.1) with a test function v H () H0(), 1 integrating over, and applying integration by parts twice, this gives 0 = [ u(x)]v(x) dx = u(x) v(x)dx = u(x)[ v(x)]dx+ u(x) n x v(x)ds x. When inserting the Dirichlet boundary condition, this results in the Galerkin Petrov variational formulation to find u X U = L 2 () such that a U (u,v) := u(x)[ v(x)]dx = g(x) v(x)ds x (3.1) n x is satisfied for all v Y U = H () H0 1 (). In this case we have to assume that the given Dirichlet datum g is in the dual of the interior Neumann trace space, i.e. g [γ int 1 H ()]. We obviously have X U = Y S and Y U = X S, respectively, with the associated norms [ 1/2. u XU = u L 2 (), v YU = v 2 L 2 () + v 2 L ()] 2 Similar as in Lemma 2.1 we can prove boundedness, and the stability condition a U (u,v) u XU v YU for all u X U, v Y U, c S u XU sup 0 v Y U a U (u,v) v YU for all u X U. As a consequence, we conclude unique solvability of the variational problem (3.1) to obtain u X U = L 2 (). 4 Single layer potential Wenowconsiderthesinglelayerpotential(1.2),u(x) = (Ṽw)(x),x R3 \. Whendefining g = γ int 0 Ṽw, we observe that u is a solution of the Dirichlet boundary value problem (1.1) being also the unique solution of the strong Galerkin Petrov formulation (2.1). To ensure u H (), we chose ψ [H ()] and we consider Ṽw,ψ H () [H ()] = ψ(x) U (x,y)w(y)ds y dx = w(y) U (x,y)ψ(x)dxds y = ϕ,w 6
11 where the duality pairing has to be specified. Using the Newton potential ϕ(y) = (N 0 ψ)(y) = U (x,y)ψ(x)dx for y and the Dirichlet datum φ(y) = ϕ(y) = (N 0 ψ)(y) for y, we note that ϕ X U is the solution of the Dirichlet boundary value problem y ϕ(y) = ψ(y) for y, ϕ(y) = φ(y) for y. In fact, ϕ X U solves, due to ψ [H ()] Y U = [H () H0()] 1, the ultra week variational formulation ϕ(x) v(x) dx = φ(x) v(x)ds x ψ(x)v(x)dx for all v Y U. n x This variational formulation implies [ ] φ = ϕ γ int [ ] 1 Y U = γ int 1 [H0 1 () H ()] from which we further conclude as well as and finally, w γ int 1 [H 1 0 () H ()] = γ int 1 H () Ṽ : γ int 1 H () H (), V : γ int 1 H () γ int 0 H () (4.1) follows. In particular, when H () H s ()is satisfied for 3 2 < s < s i, we have V : H s 3 2 () H s 1 2 (). Since we can do the same considerations subject to the exterior problem, we finally conclude V : H s 1 () H s () for all s (1,min{s i 1/2,s e 1/2}). (4.2) Note that in the case of a polygonal bounded domain R 2 this result was already given in [4]. Due to the mapping properties (4.2) and when assuming g H s () for some s (1,min{s i,s e }) we may consider the Galerkin Petrov variational formulation to find w H s 1 () such that Vw,v = g,v for all v H s (). (4.3) To prove unique solvability of the Galerkin Petrov formulation (4.3) we need to establish an appropriate stability condition. 7
12 Theorem 4.1 Let w H s 1 () be given for some s (1,min{s i,s e }). Then there holds the stability condition Vw,v c V w H s 1 () (4.4) v H s () sup 0 v H s () with a positive constant c V > 0 independent of w. Proof. In fact, we follow the standard approach to prove the ellipticity estimate (1.5), see, e.g., [6, 7]. Since u = Ṽw is harmonic in, Green s first formula implies u(x)u(x)ds x = u(x) 2 dx. (4.5) n x With the jump relations for the interior Dirichlet and Neumann trace operators for the single layer potential we have γ int 0 u(x) = (Vw)(x), γ int 1 u(x) = 1 2 w(x)+(k w)(x) for x, where in addition to the single layer boundary integral operator V we used the adjoint double layer boundary integral operator, (K w)(x) = U (x,y)w(y)ds y, x. n x Hence, (4.5) gives ( 1 2 I +K )w,vw = u 2 L 2 (). When doing the same considerations subject to the exterior Dirichlet boundary value problem, this gives ( 1 2 I K )w,vw = u 2 L 2 ( c ). Note that in R 3 the single layer potential u = Now we conclude, by using (2.7) and (2.10), Ṽw satisfies the radiation condition (2.9). w,vw = u 2 L 2 () + u 2 L 2 ( c ) c i γ int 1 u 2 H s 1 () +c e γ ext 1 u 2 H s 1 (). On the other hand we have ( ) w 2 H s 1 () = γint 1 u γ ext 1 u 2 H s 1 () 2 γ int 1 u 2 H s 1 () + γext 1 u 2 H s 1 () and therefore we obtain w,vw c S w 2 H s 1 () c V w H s 1 () w H s () due to H s 1 () H s () for s > 1, which finally gives (4.4). 8
13 From the stability condition(4.4) we can conclude unique solvability of the Galerkin Petrov formulation (4.3). By using Vw,v = w,vv for w H s 1 (),v H s () we can define the single layer boundary integral operator V : H s () H 1 s (), which satisfies the following stability condition. Lemma 4.2 Let v H s be given for some s (1,min{s i,s e }). Then there holds the stability condition Vv,w c V v H s () (4.6) w H s 1 () sup 0 w H s 1 () with the positive constant c V > 0 as used in (4.4). Proof. For g H s () we find, by solving (4.3), w H s 1 (), and the stability condition (4.4) gives c V w H s 1 () sup 0 v H s () Vw,v v H s () = sup 0 v H s () g,v v H s () g H s (). Using duality we then conclude the stability estimate as claimed. v,g v H s () = sup 0 g H s () g H s () 1 c V sup 0 w H s 1 () v,vw w H s 1 () Using the indirect single layer potential u(x) = (Ṽv)(x), x we candescribe thesolution of the Dirichlet boundary value problem (1.1) with a given Dirichlet datum g H 1 s (), i.e. v H s () is the unique solution of the Galerkin Petrov formulation Vv,w = g,w for all w H s 1 (). (4.7) Due to s > 1 it is possible to consider g L 2 () H 1 s () within the variational formulation (4.7) which can be seen as the boundary integral equation counter part of the ultra week finite element formulation [1]. Remark 4.1 In the two dimensional case R 2 we need to have w H s 1 () with the constraint w,1 = 0 to satisfy the radiation condition (2.9) for the single layer potential u(x) = (Ṽw)(x) as x. To ensure solvability of the Galerkin Petrov formulation (4.3), g has to satisfy the solvability condition g,w eq = 0 with the natural density w eq = V 1 1. So solvability of the Dirichlet boundary value problem (1.1) for general g can always be guaranteed when considering an appropriate additive splitting of g(x) = γ g + g(x), γ g = g,w eq 1,weq, where we have to assume diam < 1 to ensure 1,w eq > 0. All other results then remain true when considering appropriate factor spaces. 9
14 References [1] Apel, T., Nicaise, S., Pfefferer, J.: A dual singular complement method for the numerical solution of the Poisson equation with L 2 boundary data in non-convex domains. Num. Methods Part. Diff. Eq., published electronically (2016). [2] Costabel, M.: Boundary integral operators on Lipschitz domains: Elementary results. SIAM J. Math. Anal., 19, (1988). [3] Costabel, M.: Boundary integral operators for the heat equation. Integral Eqns. Operator Th., 13, (1990). [4] Costabel, M., Stephan, E. P.: Boundary integral equations for mixed boundary value problems in polygonal domains and Galerkin approximations. In Mathematical Models and Methods in Mechanics, Banach Centre Publ., vol. 15, PWN, Warschau (1985), pp [5] Grisvard, P.: Singularities in boundary value problems, Research Notes in Applied Mathematics, vol. 22, Springer, New York (1992). [6] Hsiao, G. C., Wendland, W. L.: A finite element method for some integral equations of the first kind. J. Math. Anal. Appl., 19, (1977). [7] Steinbach, O.: Numerical Approximation Methods for Elliptic Boundary Value Problems. Finite and Boundary Elements, Springer, New York (2008). [8] Verchota, G.: Layer potentials and regularity for the Dirichlet problem for Laplace s equation in Lipschitz domains, J. Funct. Anal., 59, (1984). 10
Technische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationSpace time finite and boundary element methods
Space time finite and boundary element methods Olaf Steinbach Institut für Numerische Mathematik, TU Graz http://www.numerik.math.tu-graz.ac.at based on joint work with M. Neumüller, H. Yang, M. Fleischhacker,
More informationIntegral Representation Formula, Boundary Integral Operators and Calderón projection
Integral Representation Formula, Boundary Integral Operators and Calderón projection Seminar BEM on Wave Scattering Franziska Weber ETH Zürich October 22, 2010 Outline Integral Representation Formula Newton
More informationTechnische Universität Graz
Technische Universität Graz Error Estimates for Neumann Boundary Control Problems with Energy Regularization T. Apel, O. Steinbach, M. Winkler Berichte aus dem Institut für Numerische Mathematik Bericht
More informationTechnische Universität Graz
Technische Universität Graz Schur complement preconditioners for the biharmonic Dirichlet boundary value problem L. John, O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2013/4
More informationTechnische Universität Graz
Technische Universität Graz Robust boundary element domain decomposition solvers in acoustics O. Steinbach, M. Windisch Berichte aus dem Institut für Numerische Mathematik Bericht 2009/9 Technische Universität
More informationTechnische Universität Graz
Technische Universität Graz Modified combined field integral equations for electromagnetic scattering O. Steinbach, M. Windisch Berichte aus dem Institut für Numerische Mathematik Bericht 2007/6 Technische
More informationTechnische Universität Graz
Technische Universität Graz A non-symmetric coupling of the finite volume method and the boundary element method C. Erath, G. Of, F. J. Sayas Berichte aus dem Institut für Numerische Mathematik Bericht
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 informationDomain decomposition methods via boundary integral equations
Domain decomposition methods via boundary integral equations G. C. Hsiao a O. Steinbach b W. L. Wendland b a Department of Mathematical Sciences, University of Delaware, Newark, Delaware 19716, USA. E
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 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 informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
More informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
More 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 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 informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
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 informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
More informationThe Helmholtz Equation
The Helmholtz Equation Seminar BEM on Wave Scattering Rene Rühr ETH Zürich October 28, 2010 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness
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 informationThe Mortar Boundary Element Method
The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract
More informationA Space-Time Boundary Element Method for the Wave Equation
W I S S E N T E C H N I K L E I D E N S C H A F T A Space-Time Boundary Element Method for the Wave Equation Marco Zank and Olaf Steinbach Institut für Numerische Mathematik Space-Time Methods for PDEs,
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
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 informationFact Sheet Functional Analysis
Fact Sheet Functional Analysis Literature: Hackbusch, W.: Theorie und Numerik elliptischer Differentialgleichungen. Teubner, 986. Knabner, P., Angermann, L.: Numerik partieller Differentialgleichungen.
More informationSobolev Spaces 27 PART II. Review of Sobolev Spaces
Sobolev Spaces 27 PART II Review of Sobolev Spaces Sobolev Spaces 28 SOBOLEV SPACES WEAK DERIVATIVES I Given R d, define a multi index α as an ordered collection of integers α = (α 1,...,α d ), such that
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationSTEKLOV REPRESENTATIONS OF GREEN S FUNCTIONS FOR LAPLACIAN BOUNDARY VALUE PROBLEMS. 1. Introduction
STEKLOV REPRESENTATIONS OF GREEN S FUNCTIONS FOR LAPLACIAN BOUNDARY VALUE PROBLEMS. GILES AUCHMUTY Abstract. This paper describes different representations for solution operators of Laplacian boundary
More information4 Sobolev spaces, trace theorem and normal derivative
4 Sobolev spaces, trace theorem and normal derivative Throughout, n will be a sufficiently smooth, bounded domain. We use the standard Sobolev spaces H 0 ( n ) := L 2 ( n ), H 0 () := L 2 (), H k ( n ),
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationDIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION
Meshless Methods in Science and Engineering - An International Conference Porto, 22 DIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION Robert Schaback Institut für Numerische und Angewandte Mathematik (NAM)
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 informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
More informationTechnische Universität Graz
Technische Universität Graz On the Numerical Solution of Minimal Energy Problems G. Of, W. L. Wendland, N. Zorii Berichte aus dem Institut für Numerische Mathematik Bericht 2009/6 Technische Universität
More informationTechnische Universität Graz
Technische Universität Graz An energy space finite element approach for elliptic Dirichlet boundary control problems G. Of, T. X. Phan, O. Steinbach Berichte aus dem Institut für Numerische Mathematik
More informationTechnische Universität Graz
echnische Universität Graz Coercive space time finite element methods for initial boundary value problems O. Steinbach, M. Zank Berichte aus dem Institut für Angewandte Mathematik Bericht 18/7 echnische
More informationweak Galerkin, finite element methods, interior estimates, second-order elliptic
INERIOR ENERGY ERROR ESIMAES FOR HE WEAK GALERKIN FINIE ELEMEN MEHOD HENGGUANG LI, LIN MU, AND XIU YE Abstract Consider the Poisson equation in a polytopal domain Ω R d (d = 2, 3) as the model problem
More informationu = f in Ω, u = q on Γ. (1.2)
ERROR ANALYSIS FOR A FINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET BOUNDARY CONTROL PROBLEMS S. MAY, R. RANNACHER, AND B. VEXLER Abstract. We consider the Galerkin finite element approximation of
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationTraces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains
Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis
More informationFractional Operators with Inhomogeneous Boundary Conditions: Analysis, Control, and Discretization
Fractional Operators with Inhomogeneous Boundary Conditions: Analysis, Control, and Discretization Harbir Antil Johannes Pfefferer Sergejs Rogovs arxiv:1703.05256v1 [math.na] 15 Mar 2017 March 16, 2017
More informationTechnische Universität Graz
Technische Universität Graz Analysis of a kinematic dynamo model with FEM BEM coupling W. Lemster, G. Lube, G. Of, O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 0/3 Technische
More informationFinite Element Error Estimates in Non-Energy Norms for the Two-Dimensional Scalar Signorini Problem
Journal manuscript No. (will be inserted by the editor Finite Element Error Estimates in Non-Energy Norms for the Two-Dimensional Scalar Signorini Problem Constantin Christof Christof Haubner Received:
More informationMaximum-norm stability of the finite element Ritz projection with mixed boundary conditions
Noname manuscript No. (will be inserted by the editor) Maximum-norm stability of the finite element Ritz projection with mixed boundary conditions Dmitriy Leykekhman Buyang Li Received: date / Accepted:
More informationApproximation of the biharmonic problem using P1 finite elements
Approximation of the biharmonic problem using P1 finite elements Robert Eymard, Raphaèle Herbin, Mohamed Rhoudaf To cite this version: Robert Eymard, Raphaèle Herbin, Mohamed Rhoudaf. Approximation of
More informationThe All-floating BETI Method: Numerical Results
The All-floating BETI Method: Numerical Results Günther Of Institute of Computational Mathematics, Graz University of Technology, Steyrergasse 30, A-8010 Graz, Austria, of@tugraz.at Summary. The all-floating
More informationFractional Operators with Inhomogeneous Boundary Conditions: Analysis, Control, and Discretization
Fractional Operators with Inhomogeneous Boundary Conditions: Analysis, Control, and Discretization Harbir Antil Johannes Pfefferer Sergejs Rogovs arxiv:1703.05256v2 [math.na] 11 Sep 2017 September 12,
More informationEXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM. Saeid Shokooh and Ghasem A. Afrouzi. 1. Introduction
MATEMATIČKI VESNIK MATEMATIQKI VESNIK 69 4 (217 271 28 December 217 research paper originalni nauqni rad EXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM Saeid Shokooh and Ghasem A.
More informationNONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian UNDER NONHOMOGENEOUS NEUMANN BOUNDARY CONDITION
Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 210, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu NONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian
More informationWhen is the error in the h BEM for solving the Helmholtz equation bounded independently of k?
BIT manuscript No. (will be inserted by the editor) When is the error in the h BEM for solving the Helmholtz equation bounded independently of k? I. G. Graham M. Löhndorf J. M. Melenk E. A. Spence Received:
More informationDivergence Boundary Conditions for Vector Helmholtz Equations with Divergence Constraints
Divergence Boundary Conditions for Vector Helmholtz Equations with Divergence Constraints Urve Kangro Department of Mathematical Sciences Carnegie Mellon University Pittsburgh, PA 15213 Roy Nicolaides
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 informationMINIMAL GRAPHS PART I: EXISTENCE OF LIPSCHITZ WEAK SOLUTIONS TO THE DIRICHLET PROBLEM WITH C 2 BOUNDARY DATA
MINIMAL GRAPHS PART I: EXISTENCE OF LIPSCHITZ WEAK SOLUTIONS TO THE DIRICHLET PROBLEM WITH C 2 BOUNDARY DATA SPENCER HUGHES In these notes we prove that for any given smooth function on the boundary of
More informationA Remark on the Regularity of Solutions of Maxwell s Equations on Lipschitz Domains
A Remark on the Regularity of Solutions of Maxwell s Equations on Lipschitz Domains Martin Costabel Abstract Let u be a vector field on a bounded Lipschitz domain in R 3, and let u together with its divergence
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 informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More informationFINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET OPTIMAL CONTROL PROBLEMS
Numerical Functional Analysis and Optimization, 28(7 8):957 973, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0163-0563 print/1532-2467 online DOI: 10.1080/01630560701493305 FINITE ELEMENT APPROXIMATION
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationSobolev spaces, Trace theorems and Green s functions.
Sobolev spaces, Trace theorems and Green s functions. Boundary Element Methods for Waves Scattering Numerical Analysis Seminar. Orane Jecker October 21, 2010 Plan Introduction 1 Useful definitions 2 Distributions
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 informationOn second order sufficient optimality conditions for quasilinear elliptic boundary control problems
On second order sufficient optimality conditions for quasilinear elliptic boundary control problems Vili Dhamo Technische Universität Berlin Joint work with Eduardo Casas Workshop on PDE Constrained Optimization
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 informationAn Adaptive Space-Time Boundary Element Method for the Wave Equation
W I S S E N T E C H N I K L E I D E N S C H A F T An Adaptive Space-Time Boundary Element Method for the Wave Equation Marco Zank and Olaf Steinbach Institut für Numerische Mathematik AANMPDE(JS)-9-16,
More informationChapter 7: Bounded Operators in Hilbert Spaces
Chapter 7: Bounded Operators in Hilbert Spaces I-Liang Chern Department of Applied Mathematics National Chiao Tung University and Department of Mathematics National Taiwan University Fall, 2013 1 / 84
More informationFunctional Analysis. Martin Brokate. 1 Normed Spaces 2. 2 Hilbert Spaces The Principle of Uniform Boundedness 32
Functional Analysis Martin Brokate Contents 1 Normed Spaces 2 2 Hilbert Spaces 2 3 The Principle of Uniform Boundedness 32 4 Extension, Reflexivity, Separation 37 5 Compact subsets of C and L p 46 6 Weak
More informationThe continuity method
The continuity method The method of continuity is used in conjunction with a priori estimates to prove the existence of suitably regular solutions to elliptic partial differential equations. One crucial
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 informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
More informationNONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS
Fixed Point Theory, Volume 9, No. 1, 28, 3-16 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.html NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS GIOVANNI ANELLO Department of Mathematics University
More informationTechnische Universität Graz
echnische Universität Graz A Stabilized Space ime Finite Element Method for the Wave Equation Olaf Steinbach and Marco Zank Berichte aus dem Institut für Angewandte Mathematik Bericht 18/5 echnische Universität
More informationAnisotropic mesh refinement in polyhedral domains: error estimates with data in L 2 (Ω)
Anisotropic mesh refinement in polyhedral domains: error estimates with data in L 2 (Ω) Thomas Apel Ariel L. Lombardi Max Winkler February 6, 2014 arxiv:1303.2960v1 [math.na] 12 Mar 2013 Abstract. The
More informationA priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces
A priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces A. Bespalov S. Nicaise Abstract The Galerkin boundary element discretisations of the
More informationSolving Negative Order Equations by the Multigrid Method Via Variable Substitution
J Sci Comput (214) 59:371 385 DOI 1.17/s1915-13-9762-4 Solving Negative Order Equations by the Multigrid Method Via Variable Substitution George C. Hsiao Liwei Xu Shangyou Zhang Received: 26 December 212
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationAsymptotic behavior of infinity harmonic functions near an isolated singularity
Asymptotic behavior of infinity harmonic functions near an isolated singularity Ovidiu Savin, Changyou Wang, Yifeng Yu Abstract In this paper, we prove if n 2 x 0 is an isolated singularity of a nonegative
More informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
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 informationJOHANNES KEPLER UNIVERSITY LINZ. Shape-Explicit Constants for Some Boundary Integral Operators
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics Shape-Explicit Constants for Some Boundary Integral Operators Clemens Pechstein Institute of Computational Mathematics, Johannes Kepler
More informationA RELATIONSHIP BETWEEN THE DIRICHLET AND REGULARITY PROBLEMS FOR ELLIPTIC EQUATIONS. Zhongwei Shen
A RELATIONSHIP BETWEEN THE DIRICHLET AND REGULARITY PROBLEMS FOR ELLIPTIC EQUATIONS Zhongwei Shen Abstract. Let L = diva be a real, symmetric second order elliptic operator with bounded measurable coefficients.
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition Sukjung Hwang CMAC, Yonsei University Collaboration with M. Dindos and M. Mitrea The 1st Meeting of
More informationBesov regularity of solutions of the p-laplace equation
Besov regularity of solutions of the p-laplace equation Benjamin Scharf Technische Universität München, Department of Mathematics, Applied Numerical Analysis benjamin.scharf@ma.tum.de joint work with Lars
More informationMixed exterior Laplace s problem
Mixed exterior Laplace s problem Chérif Amrouche, Florian Bonzom Laboratoire de mathématiques appliquées, CNRS UMR 5142, Université de Pau et des Pays de l Adour, IPRA, Avenue de l Université, 64000 Pau
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 informationWhen all else fails, integrate by parts an overview of new and old variational formulations for linear elliptic PDEs
When all else fails, integrate by parts an overview of new and old variational formulations for linear elliptic PDEs E. A. Spence December 31, 2014 Abstract We give an overview of variational formulations
More informationPARTIAL DIFFERENTIAL EQUATIONS. Lecturer: D.M.A. Stuart MT 2007
PARTIAL DIFFERENTIAL EQUATIONS Lecturer: D.M.A. Stuart MT 2007 In addition to the sets of lecture notes written by previous lecturers ([1, 2]) the books [4, 7] are very good for the PDE topics in the course.
More informationLocal pointwise a posteriori gradient error bounds for the Stokes equations. Stig Larsson. Heraklion, September 19, 2011 Joint work with A.
Local pointwise a posteriori gradient error bounds for the Stokes equations Stig Larsson Department of Mathematical Sciences Chalmers University of Technology and University of Gothenburg Heraklion, September
More informationAn Introduction to Variational Inequalities
An Introduction to Variational Inequalities Stefan Rosenberger Supervisor: Prof. DI. Dr. techn. Karl Kunisch Institut für Mathematik und wissenschaftliches Rechnen Universität Graz January 26, 2012 Stefan
More informationSEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT
Electronic Journal of Differential Equations, Vol. 2012 (2012), No. 139, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SEMILINEAR ELLIPTIC
More informationA Posteriori Error Bounds for Meshless Methods
A Posteriori Error Bounds for Meshless Methods Abstract R. Schaback, Göttingen 1 We show how to provide safe a posteriori error bounds for numerical solutions of well-posed operator equations using kernel
More informationUniqueness of Positive Solutions for a Class of p-laplacian Systems with Multiple Parameters
Int. Journal of Math. Analysis, Vol. 2, 2008, no. 2, 005-03 Uniqueness of Positive Solutions for a Class of p-laplacian Systems with Multiple Parameters G. A. Afrouzi and E. Graily Department of Mathematics,
More informationCorner Singularities
Corner Singularities Martin Costabel IRMAR, Université de Rennes 1 Analysis and Numerics of Acoustic and Electromagnetic Problems Linz, 10 15 October 2016 Martin Costabel (Rennes) Corner Singularities
More informationAn ill-posed boundary value problem for the Helmholtz equation on Lipschitz domains
c de Gruyter 2009 J. Inv. Ill-Posed Problems 17 (2009), 703 711 DOI 10.1515 / JIIP.2009.041 An ill-posed boundary value problem for the Helmholtz equation on Lipschitz domains W. Arendt and T. Regińska
More informationOPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS
PORTUGALIAE MATHEMATICA Vol. 59 Fasc. 2 2002 Nova Série OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS J. Saint Jean Paulin and H. Zoubairi Abstract: We study a problem of
More information