LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
|
|
- Kristin Walker
- 5 years ago
- Views:
Transcription
1 LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1
2 Main Problem in Mathematical Modeling and Numerical Analysis One of the most challenging problems in mathematical modeling and numerical analysis, which also has great importance for most of applications in industry is RELIABLE VERIFICATION OF ACCURACY OF APPROXIMATE SOLUTIONS OBTAINED IN COMPUTER SIMULATIONS Mathematically, this task is related to the so-called a posteriori error estimates, giving computable bounds for errors of various types and detecting zones, where such errors are excessively high and some mesh-refinement algorithm should be used 2
3 SOURCES OF ERRORS IN MODELING AND COMPUTATIONS In this section we discuss sources of the main types of errors in modeling and computational phases 3
4 Modeling Error In modeling phase a physical process (or object) of interest is usually described by means of a certain mathematical model. Then, the so-called modeling error ε 1 = U u arises, where U is the value characterizing the process (object) and u is the corresponding value obtained from the mathematical model. This error exists due to, e.g., the fact that usually second-order phenomena are neglected, there is also an undeterminacy in the problem data, often dimensional reduction is used to simplify models, etc. The symbol denotes a convenient measure of the difference between U and u (e.g., the absolute value of the difference or a norm in a suitable functional space) 4
5 Discretization Error In most of cases, the exact solution u cannot be obtained in analytical (explicit) form due to a high complexity of models constructed for real-life problems. Indeed, adequate models for complicated processes normally involve several differential and integral equations, various algebraic relations, etc. The solution u of such complicated models is often understood in an abstract sense as an element of a certain functional space Qualitative properties of u can be studied by pure mathematical methods, but the quantitative analysis normally requires a replacement of the original problem by a sequence of simpler problems whose solutions can be found relatively easily. Let u h be solution of approximate problem obtained on mesh with characteristic size h. Then discretization error ε 2 = u u h appears 5
6 Computational Error In turn, those approximate problems are themselves solved approximately, using concrete computers and concrete software packages, so the third type error, called the computational error, ε 3 = u h u ε h, appears, where uε h is what we really obtain in computer simulations The error ε 3 includes roundoff errors, errors due to forcibly stopped iterative processes and errors caused by bugs in computer codes. It is clear that estimation of ε 3 is an extremely difficult task All said above is presented in the following diagram 6
7 U Physical process ε 1 Modeling error u Mathematical model Au = f ε 2 Discretization error u h Discrete model A h u h = f h ε 3 Computational error u ɛ h Numerical solution A ɛ h uɛ h = f h ɛ 7
8 Two Principal Relations Computations on the base of a reliable (certified) model. Here ε 1 is assumed to be small and computed u ɛ h gives a desired information on U U u ɛ h ε 1 + ε 2 + ε 3 (1) Verification of a mathematical model. Here physical data U and computed results u ɛ h are compared to judge on the quality of mathematical model ε 1 U u ɛ h + ε 2 + ε 3 (2) 8
9 Thus, two major problems of mathematical modeling: reliable computer simulation verification of mathematical models by comparing physical and mathematical experiments require efficient methods able to provide COMPUTABLE AND REALISTIC estimates of ε 2 + ε 3 9
10 No Error Control - No Reliability In order to have really reliable computations we have to do both things: obtain an aproximate solution and also explicitly control the error ε 2 + ε 3 In practical computations the second part of the total work is often ignored. Mostly it happenes because the problem is very complicated and computable error estimates are simply unknown. In other cases, it happens because the analysts believe that approximations obtained with help of sufficiently fine meshes and powerful computers have values ε 2 and ε 3 so small that they should not be taken into an account at all 10
11 Such an approach does not provide with reliable computer simulations. Approximation solutions are quite sensitive to mesh perturbation and restructuring mesh can lead to very different results exact solution mesh with 741 nodes error = mesh with 785 nodes error = mesh with 855 nodes error =
12 Mathematical Model: solution u and norm Approximate solutions without a real error control is just a part of the whole business. However, to develop tools for the error control we must first establish rules of the game : What do we mean as a solution of the mathematical model (e.g., of the given BVP)? How do we measure the difference between the exact solution of the given model and computed approximations? In order to answer the above two questions we shall introduce main definitions and denotations ( mathematical tools ) needed 12
13 MATHEMATICAL TOOLS: DEFINITIONS and TERMINOLOGY 13
14 Basic Facts from Functional Analysis Definition: A set V is called a linear vector space, if it has the following properties: If v, w V are arbitrary and α is a real (or complex) number then v + w and αv also belong to V ; moreover, these operations obey the usual rules of algebra, i.e.: v + w = w + v, v + (w + z) = (v + w) + z α(v + w) = αv + αw, α(βv) = (αβ)v, if v + w = v + z (α + β)v = αv + βv 1 v = v then w = z 14
15 Definition: A norm on linear vector space V is a real function denoted by V, which satisfies the following conditions: v + w V v V + w V (triangle inequality) αv V = α v V v V 0 if v 0 for any v, w V and any scalar α (by a scalar we mean real or complex number) It follows easily from above that a norm is a nonnegative function on V, i.e., v V 0 for any v V A linear vector space, which is equipped with a norm, is called a normed linear vector space 15
16 Definition: Banach space V is a normed linear vector space which is complete, that is, for each Cauchy sequence {v j } j=1 V there exists an element v V such that v j v V 0 as j Definition: We say that T is a linear operator from Banach space V to Banach space U if T (v + w) = T (v) + T (w) and T (αv) = αt (v) for any v, w V and any scalar α Definition: A linear operator T : V U is continuous if there exists a constant C > 0 such that T v U C v V v V 16
17 Let V be a linear vector space and let α C denotes the complex conjugate number to α Definition: A scalar (inner) product on V V is a complex function denoted by (, ) V which satisfies the following conditions: (v + y, w) V = (v, w) V + (y, w) V (αv, w) V = α(v, w) V (v, w) V = ((w, v) V ) C and (v, v) V 0 (v, v) V 0 if v 0 for any v, w, y V and any scalar α If V is a real space then (, ) V ( ) C can be omitted is a real function and the symbol 17
18 For the so-called induced norm v V = (v, v) V, v V we have the well-known Cauchy Schwarz inequality (v, w) V v V w V v, w V Moreover v V = sup w 0 (v, w) V w V A linear operator whose values are scalars is called a linear form (linear functional) A Banach space with a scalar product is called a Hilbert space 18
19 Riesz Theorem: Let V be a Hilbert space. Then for any linear continuous form F defined on V there exists exactly one element u V such that F (v) = (v, u) V v V Definition: A scalar mapping a(, ) defined on V V, where V is a linear vector space, is said to be a sesquilinear form, if for any fixed v V the mappings a(, v) and (a(v, )) C are linear. If moreover a(, ) is only real-valued then a(, ) is called a bilinear form. A sesquilinear form a(, ) is said to be continuous if there exists a constant C 1 > 0 such that a(v, w) C 1 v V w V v, w V 19
20 Lax Milgram Lemma: Let V be a Hilbert space and let a(, ) be a continuous sesquilinear form for which there exists a constant C 2 > 0 such that a(v, v) C 2 v 2 V v V (V -ellipticity condition) Then for any linear continuous form F defined on V there exists exactly one element u V such that a(v, u) = F (v) v V ( ) Actually, the Lax Milgram lemma says that a(, ) is a norm in V equivalent to the V -norm 20
21 A sesquilinear form a(, ) is said to be Hermitian if a(v, w) = (a(w, v)) C v, w V ( ) Theorem: Let the assumptions of Lax-Milgram lemma be satisfied and let a sesquilinear form a(, ) be Hermitian. Then the problem ( ) is equivalent to the problem: Find u V such that J(u) = inf v V where J is a quadratic functional given by J(v) ( ) J(v) = 1 a(v, v) Re(F (v)), v V ( ) 2 The functional ( ) is called the energy functional. Note that in the real case, the relation ( ) reduces to a(v, w) = a(w, v) and the corresponding bilinear form is then said to be symmetric 21
22 Functional Spaces L p (Ω) and C k (Ω) Let R d stand for the Euclidean space equipped with the norm d x = j=1 x 2 j 1/2, x = (x 1,..., x d ) T R d The Lebesgue space of real or complex-valued functions defined over an open set Ω R d, which are integrable with the power p [1, ), is denoted by L p (Ω) and equipped with the norm v 0,p,Ω = 1/p v p dx, v L p (Ω) Ω 22
23 When p = 2, we write shortly 0,Ω = 0,2,Ω and define the scalar (inner) product as follows (v, w) 0,Ω = vw C dx Recall the well-known Hölder inequality vw dx v 0,p,Ω w 0,q,Ω, v L p (Ω), w L q (Ω) Ω which holds for any p, q (1, ) satisfying the equality Ω 1 p + 1 q = 1 23
24 The Lebesgue space of measurable essentially bounded functions over Ω is denoted by L (Ω) and equipped with the norm v 0,,Ω = ess sup x Ω v(x) For any p, q [1, ], p q, the following algebraic imbedding L q (Ω) L p (Ω) holds. Moreover, we have also the topological imbedding; namely, there exists a constant C > 0 such that v 0,p,Ω C v 0,q,Ω v L q (Ω) We will not specify the term real or complex in further definitions of function spaces, since these definitions are essentially the same for both cases 24
25 By Ω we denote the closure of Ω and by Ω the boundary of Ω: Ω = Ω Ω, Ω = Ω (R d \ Ω) Recall that a domain Ω is an open and connected set in R d The symbol d {1, 2,...} is solely reserved for the dimension of Ω If Ω is a bounded domain then the space of continuous functions over Ω is denoted by C(Ω) and equipped with the norm Obviously we have v C(Ω) = max x Ω v(x) v C(Ω) = v 0,,Ω v C(Ω) 25
26 The symbol C k (Ω), k {0, 1,...}, stands for the space of functions whose classical derivatives up to order k belong to C(Ω) Moreover, we set C (Ω) = k=1 C k (Ω) i.e., C (Ω) is the space of infinitely differentiable functions over Ω Finally by C 0 (Ω) we denote the space of infinitely differentiable functions with a compact support in Ω, that is C 0 (Ω) = { v C (Ω) supp v Ω } where supp v = {x Ω v(x) 0} 26
27 Sobolev Spaces H k (Ω) and W k p (Ω) The (weak) solutions of many problems are looked for in Sobolev spaces. First we briefly recall some basic properties of the Sobolev space H k (Ω), k = {0, 1, 2,...}, i.e., classes of real or complex functions defined on a domain Ω R d whose generalized derivatives up to the k-th order belong to L 2 (Ω) We shall consider Sobolev spaces defined only on bounded domains with a Lipschitz continuous boundary which is a sufficiently wide class of domains for practical purpouses Definition: A bounded set Ω R d is said to have a Lipschitz continuous boundary if for any z Ω there exists a neighbourhood U = U(z) such that the set U Ω can be expressed, in some Cartesian coordinate system (x 1,..., x d ), by the inequality x d < F (x 1,..., x d 1 ), where F is a Lipschitz continuous function 27
28 Denote by L the set of all bounded domains with a Lipschitz continuous boundary If Ω L then the outward normal exists almost everywhere (a.e.) on Ω. This is the reason why we will confine ourselves to domains with a Lipschitz continuous boundary, since later we shall deal with, e.g., the normal derivative, normal component of flux, etc Another reason is that there are several definitions of the Sobolev spaces and in the case Ω / L the Sobolev spaces need not be uniquely defined Also, some important theorems formulated later on need not be valid when Ω / L Let us point out that there are other (nonequivalent) definitions of the Lipschitz continuity of Ω From now on, we always assume that Ω L 28
29 Definition: For any v C (Ω) and the index m = (m 1,..., m d ) we define the m-th classical derivative as follows D m v = m v x m 1 1 x m d d where m 1,..., m d are non-negative integers and m = m m d Definition: A function v L 2 (Ω) is said to have the m-th generalized derivative in L 2 (Ω) if there exists a function z L 2 (Ω) such that zw dx = ( 1) m vd m w dx w C0 (Ω) Ω Ω Then z is called the m-th generalized derivative of v and we set D m v = z. One may easily check that D m v is well-defined. Note that the generalized derivative is sometimes looked for in the larger space L 1,loc (Ω) 29
30 Now for k = 0, 1,..., the Sobolev spaces H k (Ω) are defined as H k (Ω) = { v L 2 (Ω) D m v L 2 (Ω), m k } It can be verified that H k (Ω), equipped with the scalar product (v, w) k,ω = D m v(d m w) C dx, v, w H k (Ω) m k Ω is a Hilbert space. Let us further introduce the induced norm and seminorm v k,ω = 1/2 D m v 2 dx, v H k (Ω) v k,ω = m k m =k Ω Ω D m v 2 dx 1/2, v H k (Ω) 30
31 The same symbols will be used also for vector functions, i.e., ( q ) 1/2 v k,ω = v l 2 k,ω for v = (v 1,..., v q ) T (H k (Ω)) q l=1 Moreover, the subscript Ω will be often omitted, i.e., Clearly, and (, ) k = (, ) k,ω, k = k,ω, k = k,ω v k 1 v k v H k (Ω), k = 1, 2,... L 2 (Ω) = H 0 (Ω) H 1 (Ω) H 2 (Ω). Note that each classical derivative is also the generalized derivative, and thus we have C k (Ω) H k (Ω), k = 0, 1,
32 Definition: A set Γ Ω is said to be a (relatively) open set in Ω if for any x Γ there exists an open ball B R d containing x such that B Ω Γ The Lebesgue space of square integrable function over an open set Γ Ω is denoted by L 2 (Γ) and equipped with the standard norm v 0,Γ = 1/2 v 2 ds, v L 2 (Γ) Γ Further we recall some important properties of Sobolev spaces 32
33 Trace Theorem: Let Ω L. Then there exists exactly one linear continuous operator γ : H 1 (Ω) L 2 ( Ω) such that γv = v v C (Ω) Ω The function γv for v H 1 (Ω) is called the trace of v and we denote it, for simplicity, by v Ω Trace Theorem, in fact, says that there exists a constant C > 0 such that v 0, Ω C v 1,Ω v H 1 (Ω) The trace theorem enables us to define the space H 1 0 (Ω) = { v H 1 (Ω) v = 0 on Ω } 33
34 Note that the spaces H 1 0 (Ω) and H 1 (Ω) can also be defined as the completion of C 0 (Ω) and C (Ω), respectively, under the norm 1,Ω, i.e., H 1 0 (Ω) = C 0 (Ω), H1 (Ω) = C (Ω) Denote by H 1/2 ( Ω) the space of traces of all functions from H 1 (Ω). Then we have the following density relation L 2 ( Ω) = H 1/2 ( Ω) where the closure is taken under the 0, Ω norm 34
35 Rellich Theorem: Let Ω L. Then the identity mapping from H 1 (Ω) to L 2 (Ω) is compact (i.e., any bounded sequence in H 1 (Ω) contains a subsequence converging in L 2 (Ω)) Green Theorem: Let Ω L. Then for each i {1,..., d} w i v dx + v i w dx = n i vw ds v, w H 1 (Ω) Ω Ω where n i are the components of the outward unit normal to Ω and i v = v x i Sobolev Imbedding Theorem: Let k be an integer such that 2k > d. Then H k (Ω) C(Ω) Ω and there exists a constant C > 0 such that v k C v C(Ω) v H k (Ω) 35
36 Generalized Poincaré s Inequality: Let Ω L and let ω be an open set either in Ω or in Ω 0, where Ω 0 Ω, Ω 0 L. Then v 1 C ( v ω where dω stands for dx or ds v dω 2 ) 1/2 v H 1 (Ω) Corollary: By the Cauchy Schwarz inequality 2 v dω meas ω v 2 dω v H 1 (Ω) ω Hence, we have ω v 1 C( v v 2 0,ω) 1/2 v H 1 (Ω) This inequality is usually called Friedrichs inequality when ω Ω. In particular, v 1 C v 1 v H 1 0 (Ω) 36
37 Finally in this section we briefly introduce the Sobolev spaces Wp k (Ω) of functions whose generalized derivatives up to order k {0, 1,...} belong to L p (Ω), p [1, ]. For p < they are equipped with the norm 1 v k,p,ω = p D m v p dx, v Wp k (Ω) m k Ω and seminorm v k,p,ω = m =k Ω D m v p dx 1 p, v W k p (Ω) 37
38 For p = v k,,ω = max m k Dm v 0,,Ω, We will again write k,p = k,p,ω, v k,,ω = max m =k Dm v 0,,Ω k,p = k,p,ω The Sobolev space W k p (Ω) is the completion of C (Ω) in the k,p -norm. Sobolev spaces are Banach spaces and for p = 2 they become Hilbert spaces, i.e., H k (Ω) = W k 2 (Ω) The Sobolev imbedding theorem for p [1, ] and k {1, 2,...} takes the form Wp k (Ω) C(Ω) if pk > d Note that the condition pk > d is not necessary but only sufficient. For instance, we have W1 2 (Ω) C(Ω) for d = 2 38
39 Classical and Variational Formulation of 2 nd order Elliptic Problem Using the model elliptic boundary value problem with mixed boundary conditions (which describes, for example, an electric or magnetic potential of a stationary linear problems) we shall demonstrate first the main idea of what we understand as the solution of BVPs, i.e., what our u is Classical Formulation: Find u C 2 (Ω) such that div (A u) = f in Ω (1) u = u 0 on Γ D (2) n T A u = g on Γ N (3) 39
40 Conditions: Assume that Ω L, n is the outward unit normal, A = (a ij ) (C 1 (Ω)) d d is a symmetric matrix which is uniformly positive definite, i.e., there exist constants C 1, C 2 > 0 such that C 1 ξ 2 ξ T A(x)ξ C 2 ξ 2 ξ R d x Ω f C(Ω), u 0 C 2 (Ω), g C(Γ N ), and Γ D, Γ N are relatively open sets with respect to the topology on Ω, Γ D Γ N = Ω and meas d 1 (Γ D Γ N ) = 0 (i.e., Γ D and Γ N are disjoint). In addition, we assume from that Γ D and Γ N have a finite number of components 40
41 Function u C 2 (Ω) satisfying (1) (3) is called the classical solution Note that the classical solution is sometimes looked for in the space C 1 (Ω) C 2 (Ω), i.e., u C 1 (Ω) and u Ω C 2 (Ω) The conditions (2) and (3) are called, respectively, Dirichlet and Neumann boundary conditions The problem (1) (3) is referred as the Dirichlet and Neumann problem, respectively, if Γ N = and Γ D = The matrix A describes physical properties of the medium Ω. If A is independent of x Ω we call the medium homogeneous (otherwise inhomogeneous) If a ii (x) = a 11 (x) for i = 2,..., d, and a ij (x) = 0 for i j, the medium is said to be isotropic (otherwise anisotropic) 41
42 If A is the identity matrix, (1) is called the Poisson equation If moreover f = 0, then (1) is called the Laplace equation and the associated operator u = div u = 11 u dd u is called the Laplace operator, where ii u = 2 u x 2 i In practical problems the coefficients a ij are often not continuous (they are, e.g., piecewise constant). Then the classical divergence operator in (1) cannot be employed. That is why we introduce a weak formulation of the problem (1) (3), which enables us to consider also nonsmooth functions a ij, or f, g 42
43 Weak Formulation: Let us rewrite the equation (1) as follows d i (a ij j u) = f i,j=1 and introduce the space of test functions V = { v H 1 } (Ω) v = 0 on Γ D ( ) (V ) Multiplying ( ) by an arbitrary test function v V and then integrating over Ω, we obtain i (a ij j u)v dx = fv dx Ω i,j Ω Applying now Green Theorem we find that i v a ij j u dx = f v dx + Ω i,j Ω Ω n i a ij j u v ds i,j 43
44 Since v = 0 on Γ D, we obtain A u v dx = Ω Ω fv dx + Γ N gv ds Defining now the bilinear and linear forms a(v, w) = A w v dx, v, w H 1 (Ω) (a) F (v) = Ω Ω fv dx + Γ N gv ds, v V (F ) we see that any classical solution u of (1) (3) (if it exists) satisfies the equation a(v, u) = F (v) v V From now on we shall assume that a ij L (Ω), f L 2 (Ω), g L 2 (Γ N ) and u 0 H 1 (Ω). Moreover let a ij (x) = a ji (x) and let the ellipticity condition on A holds only for a.e. x Ω. We see then that the integrals in (a) and (F ) are still well-defined 44
45 First, we recall that an integral form of a physical law is more natural than its differential form From formula (V) we begin to use implicitly the concept of completion. Instead of looking for u C 2 (Ω) satisfying the differential equation (1) with the boundary conditions (2) (3), we shall solve the integral form equation a(v, u) = F (v) v V Its solution will be searched in the Hilbert space H 1 (Ω), which is, of course, complete and bigger than C 2 (Ω). This allows us to employ some useful results from functional analysis (e.g., the Lax Milgram Lemma) and, e.g., to apply the finite element method to find approximations to u 45
46 Definition: Let Ω L and Γ D. A function u H 1 (Ω) is called a weak (or generalized) solution of the problem (1) (3) if u u 0 V and a(v, u) = F (v) v V For convenience weak solution is denoted by u as classical one Theorem: Let Ω L and Γ D. Then there exists exactly one weak solution u H 1 (Ω) of the problem (1) (3) Note that the classical solution u C 2 (Ω) of (1) (3) need not exist, in general. However, if it exists then, it is also the weak solution 46
47 Now, after we have defined what we mean as the solution of the model and discussed its existence and uniqueness, we can define how to measure the accuracy of computed approximations It is clear that the most preferable situation is when we construct approximations in the same class of functions to which the exact solution belongs, or, if it is not so, that we can easily map the computed approximations to this class. For example for the above considered problem such class of functions is H 1 (Ω), i.e., we shall measure the error in the corresponding norm 1,Ω (or in some norm equivalent to it) In most of cases the error estimation is not easy at all since the exact solution is normally unknown in analytical form. However, several different approaches aimed at the control of the accuracy of approximations have been developed. On the base of our model problem we shall survey the most important ones 47
48 A PRIORI ERROR ESTIMATES 48
49 Error Estimate via Energy Functional In this section we shall discuss so-called a priori error estimates, which are essentially based on the fact that approximations are exactly computed by the finite element method, i.e., the error ε 3 is neglected in this approach The weak solution u can be equivalently found as follows: Find a function u which minimizes the quadratic functional J(v) = 1 2 a(v, v) F (v) + a(v, u 0) ) over space V, and then set u = u 0 + u ( u 0 + V If w = u 0 + w is another function from u 0 + V then ( ) a(u w, u w) = 2 J(w ) J(u ) (EE) 49
50 Projection Estimate It is clear that we can rarely find the exact (weak) solution u of the problem (1) (3) in an analytical form. Usually we can only solve it numerically, e.g., by the finite element method For this purpose we introduce some finite-dimensional subspace V h of V and solve the problem: Find a function u h such that u h u 0 V h and a(v h, u h ) = F (v h ) v h V h 50
51 The approximate solution u h can be equivalently found as follows: Find a function u h which minimizes the quadratic functional J(v) = 1 2 a(v, v) F (v) + a(v, u 0) over space V h, and then set u h = u 0 + u h Thus, if v h = u 0 + v h is another function from u 0 + V h then Thus, we have J(u h ) J(v h ) a(u u h, u u h )=2(J(u h ) J(u )) 2(J(v h ) J(u ))=a(u v h, u v h ) From above we get the basic relation for the difference between the finite element approximation and the exact solution a(u u h, u u h )= inf vh u 0 +V h a(u v h, u v h ) 51
52 We can rewrite it also as follows a(u u h, u u h )= inf v h V h a(u v h, u v h ) The above presented projection estimate serves as a basis for deriving a priori error estimates since it can be obviously written in the following form, called the projection theorem or Céa Lemma u u h 1 C inf v h V h u v h 1 = C u π h u 1 where π h u is the projection of u onto V h The estimate above shows that the error is bounded from above by the distance between the exact solution and its orthogonal projection onto the finite-dimensional space V h 52
53 If the approximation spaces V h are constructed in such a way that for any v V inf vh V h v v h 1 0 as h 0 then it follows that the finite element approximations converge to the exact solution Estimates of the convergence rate require a detailed analysis of interpolation properties of the spaces V h. Obviously, piecewise polynomial functions form the most convenient class for constructing V h Since u is an element of some Sobolev space, the problem turns to the question: how accurately can the elements of a Sobolev space be interpolated by piecewise polynomial functions 53
54 This question is studied by the interpolation theory (see, e.g., the well-known monograph by Ciarlet) which together with projection type estimates produces a posteriori error estimates qualified in terms of a mesh size h For our model problem they have the following form provided that u u h 1 Ch k u k+1, k {1, 2,... } Exact solution is sufficiently smooth Approximation u h is the exact finite element solution Mesh elements do not degenerate during refinement process 54
55 Drawbacks of A Priori Error Estimation A priori error estimates cannot guarantee that the error monotonically decreases as h 0 The constants C s in those estimates are often not known at all or highly overestimated A priori error estimates have only theoretical meaning they show that the approximation method is correct in principle. However, in practice we are strongly interested in the error for the concrete approximation on the concrete mesh Because of these reasons, starting from 70-th a different approach for the error control has been developing. It is called as a posteriori error control for partial differential equations 55
Numerical 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 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 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 informationAppendix A Functional Analysis
Appendix A Functional Analysis A.1 Metric Spaces, Banach Spaces, and Hilbert Spaces Definition A.1. Metric space. Let X be a set. A map d : X X R is called metric on X if for all x,y,z X it is i) d(x,y)
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More 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 informationPDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM.
PDE & FEM TERMINOLOGY. BASIC PRINCIPLES OF FEM. Sergey Korotov Basque Center for Applied Mathematics / IKERBASQUE http://www.bcamath.org & http://www.ikerbasque.net 1 Introduction The analytical solution
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More 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 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 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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More 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 informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More 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 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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationPartial Differential Equations 2 Variational Methods
Partial Differential Equations 2 Variational Methods Martin Brokate Contents 1 Variational Methods: Some Basics 1 2 Sobolev Spaces: Definition 9 3 Elliptic Boundary Value Problems 2 4 Boundary Conditions,
More informationWeak Formulation of Elliptic BVP s
Weak Formulation of Elliptic BVP s There are a large number of problems of physical interest that can be formulated in the abstract setting in which the Lax-Milgram lemma is applied to an equation expressed
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
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 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 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 informationExercises - Chapter 1 - Chapter 2 (Correction)
Université de Nice Sophia-Antipolis Master MathMods - Finite Elements - 28/29 Exercises - Chapter 1 - Chapter 2 Correction) Exercise 1. a) Let I =], l[, l R. Show that Cl) >, u C Ī) Cl) u H 1 I), u DĪ).
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationAn introduction to some aspects of functional analysis
An introduction to some aspects of functional analysis Stephen Semmes Rice University Abstract These informal notes deal with some very basic objects in functional analysis, including norms and seminorms
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More 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 informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationLectures on. Sobolev Spaces. S. Kesavan The Institute of Mathematical Sciences, Chennai.
Lectures on Sobolev Spaces S. Kesavan The Institute of Mathematical Sciences, Chennai. e-mail: kesh@imsc.res.in 2 1 Distributions In this section we will, very briefly, recall concepts from the theory
More informationMath Tune-Up Louisiana State University August, Lectures on Partial Differential Equations and Hilbert Space
Math Tune-Up Louisiana State University August, 2008 Lectures on Partial Differential Equations and Hilbert Space 1. A linear partial differential equation of physics We begin by considering the simplest
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationMath 5520 Homework 2 Solutions
Math 552 Homework 2 Solutions March, 26. Consider the function fx) = 2x ) 3 if x, 3x ) 2 if < x 2. Determine for which k there holds f H k, 2). Find D α f for α k. Solution. We show that k = 2. The formulas
More informationA metric space X is a non-empty set endowed with a metric ρ (x, y):
Chapter 1 Preliminaries References: Troianiello, G.M., 1987, Elliptic differential equations and obstacle problems, Plenum Press, New York. Friedman, A., 1982, Variational principles and free-boundary
More informationA DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT
A DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT HUOYUAN DUAN, SHA LI, ROGER C. E. TAN, AND WEIYING ZHENG Abstract. To deal with the divergence-free
More informationLinear Analysis Lecture 5
Linear Analysis Lecture 5 Inner Products and V Let dim V < with inner product,. Choose a basis B and let v, w V have coordinates in F n given by x 1. x n and y 1. y n, respectively. Let A F n n be the
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationWeak Solutions, Elliptic Problems and Sobolev Spaces
3 Weak Solutions, Elliptic Problems and Sobolev Spaces 3.1 Introduction In Chapter 2 we discussed difference methods for the numerical treatment of partial differential equations. The basic idea of these
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationVector Spaces. Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms.
Vector Spaces Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms. For each two vectors a, b ν there exists a summation procedure: a +
More informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
More informationSOBOLEV SPACES AND ELLIPTIC EQUATIONS
SOBOLEV SPACES AND ELLIPTIC EQUATIONS LONG CHEN Sobolev spaces are fundamental in the study of partial differential equations and their numerical approximations. In this chapter, we shall give brief discussions
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationSobolev spaces. May 18
Sobolev spaces May 18 2015 1 Weak derivatives The purpose of these notes is to give a very basic introduction to Sobolev spaces. More extensive treatments can e.g. be found in the classical references
More informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationWeak Convergence Methods for Energy Minimization
Weak Convergence Methods for Energy Minimization Bo Li Department of Mathematics University of California, San Diego E-mail: bli@math.ucsd.edu June 3, 2007 Introduction This compact set of notes present
More informationS chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.
Sep. 1 9 Intuitively, the solution u to the Poisson equation S chauder Theory u = f 1 should have better regularity than the right hand side f. In particular one expects u to be twice more differentiable
More informationHilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space.
Hilbert Spaces Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Vector Space. Vector space, ν, over the field of complex numbers,
More information1. Subspaces A subset M of Hilbert space H is a subspace of it is closed under the operation of forming linear combinations;i.e.,
Abstract Hilbert Space Results We have learned a little about the Hilbert spaces L U and and we have at least defined H 1 U and the scale of Hilbert spaces H p U. Now we are going to develop additional
More informationUniversité de Metz. Master 2 Recherche de Mathématiques 2ème semestre. par Ralph Chill Laboratoire de Mathématiques et Applications de Metz
Université de Metz Master 2 Recherche de Mathématiques 2ème semestre Systèmes gradients par Ralph Chill Laboratoire de Mathématiques et Applications de Metz Année 26/7 1 Contents Chapter 1. Introduction
More informationIntroduction and Preliminaries
Chapter 1 Introduction and Preliminaries This chapter serves two purposes. The first purpose is to prepare the readers for the more systematic development in later chapters of methods of real analysis
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationEquations paraboliques: comportement qualitatif
Université de Metz Master 2 Recherche de Mathématiques 2ème semestre Equations paraboliques: comportement qualitatif par Ralph Chill Laboratoire de Mathématiques et Applications de Metz Année 25/6 1 Contents
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationSolutions of Selected Problems
1 Solutions of Selected Problems October 16, 2015 Chapter I 1.9 Consider the potential equation in the disk := {(x, y) R 2 ; x 2 +y 2 < 1}, with the boundary condition u(x) = g(x) r for x on the derivative
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 informationLectures on. Weak Solutions of Elliptic Boundary Value Problems
Lectures on Weak Solutions of Elliptic Boundary Value Problems S. Kesavan The Institute of Mathematical Sciences, CIT Campus, Taramani, Chennai - 600 113. e mail: kesh@imsc.res.in 2 1 Some abstract variational
More information1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:
Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not
More 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 informationINF-SUP CONDITION FOR OPERATOR EQUATIONS
INF-SUP CONDITION FOR OPERATOR EQUATIONS LONG CHEN We study the well-posedness of the operator equation (1) T u = f. where T is a linear and bounded operator between two linear vector spaces. We give equivalent
More informationFinite Element Methods for Fourth Order Variational Inequalities
Louisiana State University LSU Digital Commons LSU Doctoral Dissertations Graduate School 2013 Finite Element Methods for Fourth Order Variational Inequalities Yi Zhang Louisiana State University and Agricultural
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More informationTools from Lebesgue integration
Tools from Lebesgue integration E.P. van den Ban Fall 2005 Introduction In these notes we describe some of the basic tools from the theory of Lebesgue integration. Definitions and results will be given
More informationLecture 17. Higher boundary regularity. April 15 th, We extend our results to include the boundary. Let u C 2 (Ω) C 0 ( Ω) be a solution of
Lecture 7 April 5 th, 004 Higher boundary regularity We extend our results to include the boundary. Higher a priori regularity upto the boundary. Let u C () C 0 ( ) be a solution of Lu = f on, u = ϕ on.
More informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationOverview of normed linear spaces
20 Chapter 2 Overview of normed linear spaces Starting from this chapter, we begin examining linear spaces with at least one extra structure (topology or geometry). We assume linearity; this is a natural
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
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 informationSpectral Theory, with an Introduction to Operator Means. William L. Green
Spectral Theory, with an Introduction to Operator Means William L. Green January 30, 2008 Contents Introduction............................... 1 Hilbert Space.............................. 4 Linear Maps
More informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
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 information1 Continuity Classes C m (Ω)
0.1 Norms 0.1 Norms A norm on a linear space X is a function : X R with the properties: Positive Definite: x 0 x X (nonnegative) x = 0 x = 0 (strictly positive) λx = λ x x X, λ C(homogeneous) x + y x +
More informationDomain Perturbation for Linear and Semi-Linear Boundary Value Problems
CHAPTER 1 Domain Perturbation for Linear and Semi-Linear Boundary Value Problems Daniel Daners School of Mathematics and Statistics, The University of Sydney, NSW 2006, Australia E-mail: D.Daners@maths.usyd.edu.au
More informationA Note on the Variational Formulation of PDEs and Solution by Finite Elements
Advanced Studies in Theoretical Physics Vol. 12, 2018, no. 4, 173-179 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/astp.2018.8412 A Note on the Variational Formulation of PDEs and Solution by
More informationAnalysis in weighted spaces : preliminary version
Analysis in weighted spaces : preliminary version Frank Pacard To cite this version: Frank Pacard. Analysis in weighted spaces : preliminary version. 3rd cycle. Téhéran (Iran, 2006, pp.75.
More informationVectors. January 13, 2013
Vectors January 13, 2013 The simplest tensors are scalars, which are the measurable quantities of a theory, left invariant by symmetry transformations. By far the most common non-scalars are the vectors,
More informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
More informationNUMERICAL PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL PARTIAL DIFFERENTIAL EQUATIONS Lecture Notes, Winter 2011/12 Christian Clason January 31, 2012 Institute for Mathematics and Scientific Computing Karl-Franzens-Universität Graz CONTENTS I BACKGROUND
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationTraces and Duality Lemma
Traces and Duality Lemma Recall the duality lemma with H / ( ) := γ 0 (H ()) defined as the trace space of H () endowed with minimal extension norm; i.e., for w H / ( ) L ( ), w H / ( ) = min{ ŵ H () ŵ
More 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 informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationSobolev Spaces. Chapter Hölder spaces
Chapter 2 Sobolev Spaces Sobolev spaces turn out often to be the proper setting in which to apply ideas of functional analysis to get information concerning partial differential equations. Here, we collect
More informationB. Appendix B. Topological vector spaces
B.1 B. Appendix B. Topological vector spaces B.1. Fréchet spaces. In this appendix we go through the definition of Fréchet spaces and their inductive limits, such as they are used for definitions of function
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
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 informationON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS
Bendikov, A. and Saloff-Coste, L. Osaka J. Math. 4 (5), 677 7 ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS ALEXANDER BENDIKOV and LAURENT SALOFF-COSTE (Received March 4, 4)
More informationMATHS 730 FC Lecture Notes March 5, Introduction
1 INTRODUCTION MATHS 730 FC Lecture Notes March 5, 2014 1 Introduction Definition. If A, B are sets and there exists a bijection A B, they have the same cardinality, which we write as A, #A. If there exists
More informationPOTENTIAL THEORY OF QUASIMINIMIZERS
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 28, 2003, 459 490 POTENTIAL THEORY OF QUASIMINIMIZERS Juha Kinnunen and Olli Martio Institute of Mathematics, P.O. Box 1100 FI-02015 Helsinki University
More informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
More informationMATH 220: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY
MATH 22: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY When discussing separation of variables, we noted that at the last step we need to express the inhomogeneous initial or boundary data as
More informationANALYSIS QUALIFYING EXAM FALL 2017: SOLUTIONS. 1 cos(nx) lim. n 2 x 2. g n (x) = 1 cos(nx) n 2 x 2. x 2.
ANALYSIS QUALIFYING EXAM FALL 27: SOLUTIONS Problem. Determine, with justification, the it cos(nx) n 2 x 2 dx. Solution. For an integer n >, define g n : (, ) R by Also define g : (, ) R by g(x) = g n
More informationLecture Notes 1: Vector spaces
Optimization-based data analysis Fall 2017 Lecture Notes 1: Vector spaces In this chapter we review certain basic concepts of linear algebra, highlighting their application to signal processing. 1 Vector
More informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
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 information************************************* Applied Analysis I - (Advanced PDE I) (Math 940, Fall 2014) Baisheng Yan
************************************* Applied Analysis I - (Advanced PDE I) (Math 94, Fall 214) by Baisheng Yan Department of Mathematics Michigan State University yan@math.msu.edu Contents Chapter 1.
More information