Variational Formulations
|
|
- Aubrey Lyons
- 6 years ago
- Views:
Transcription
1 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 provide a rigorous understanding of how to solve (1.4) using the finite element method. 2.1 Computational domains We shall consider so-called Lipschitz domains as computational domains. An open, bounded, connected set R d, d = 1, 2, 3, is called Lipschitz domain, if its boundary := \ is a Lipschitz boundary, which is roughly speaking the case if there exists a finite covering of with open subsets of R d (rectangles), and is locally (in each open rectangle) the graph of a Lipschitz continuous function. This particularly implies that the boundary of a Lipschitz domain is of finite length (excluding fractal boundaries), and slightly smoother than continuous, i. e. excluding slit and cusp domains. Figure 2.1: Domains that are not Lipschitz: slit domain (left), cusp domain (middle) and domain whose boundary is the limit of a Koch curve and has infinite length (right). In almost all numerical computations only a special type of Lipschitz domains are relevant, namely domains that can be described in the widely used CAD (computer-aided design) software packages. This set of domains is usually identical to the domains that we shall define in the following. Definition 2.1. In the case d = 2 a connected domain is called a curvilinear Lipschitz polygon, if is a Lipschitz domain, and there is a finite number of open subsets Γ k, k = 1,..., P, P N (where open refers to the d 1-dimensional manifold of R d ) such that := Γ 1 Γ P, Γ k Γ l = if k l, and for each k {1,..., P } there is a C 1 -diffeomorphism ( invertible mapping of smooth manifolds) Φ k : [0, 1] Γ k. The boundary segments are called edges, their endpoints are the vertices of. 11
2 A tangential direction can be defined for all points of an edge, which allows us to define angles of the vertices, see Figure 2.2. Analogously, a normal direction can be defined, which implies that the exterior unit vectorfield n : R d is well-defined almost everywhere on. The mappings Φ k can be regarded as smooth parametrizations of the edges. Figure 2.2: Curvilinear Lipschitz polygon with added angles at vertices. An analogous notion exists in three dimensions: curved Lipschitz polyhedron. Definition 2.2. A subset R d is called a computational domain if it is bounded and its boundary is of class C k, k 1, or if it is a bounded connected interval for d = 1, a curvilinear Lipschitz polygon for d = 2, a curved Lipschitz polyhedron for d = Linear differential operators Let α N d 0 be a multi-index, i. e., a vector of d non-negative integers α = (α 1,..., α d ) N d 0. We set α := α α d and call α α := α1 d x α1 1 x α d n the partial derivative of order α. Remember that for sufficiently smooth functions all partial derivatives commute. Provided that the derivatives exist, the gradient of a function f : R d R is the column vector ( f f(x) := (x),..., f ) (x), x. x 1 x d The divergence of a vector field f = (f 1,..., f d ) : R d R d is the function f(x) := d k=1 f k x k (x), x. The differential operator := is known as Laplacian. 2.3 Integration by parts Below we assume that R d, d = 1, 2, 3, is bounded and an interval for d = 1, a Lipschitz polygon for d = 2, and a Lipschitz polyhedron for d = 3. Throughout this section we adopt the notation n = (n 1,..., n d ) T for the exterior unit normal vectorfield that is defined almost everywhere on Γ :=, see Section
3 Theorem 2.3 (Gauß theorem). If f (C 1 ()) d (C 0 ()) d, then f dx = f, n ds. By the product rule (u f) = u f + u, f for u C 1 (), f (C 1 ()) d, we deduce the first Green s formula f, u + f u dx = f, n uds Γ Γ (FGF) for all f (C 1 ()) d (C 0 ()) d and all u C 1 () C 0 (). Plugging in the special f = fe k, k = 1,..., d, where e k is the k-th unit vector, we get f u + f u dx = fu (n) k ds (2.1) x k x k for f, u C 1 () C 0 (). We may also plug f = v into (FGF), which yields v, u + v u dx = v, n u ds (2.2) for all v C 2 () C 1 (), u C 1 () C 0 (). 2.4 Distributional derivatives Definition 2.4 (Test function space). For a non-empty open set R d we denote C 0 () := {v C () : supp v := {x : v(x) 0} } the space of test functions, which is the space of functions C () with compact support. Lemma 2.5. Two integrable functions f and g defined in the bounded set are almost everywhere equal if and only if fv dx = gv dx for all test functions v C 0 (). We consider in the following functions to be equivalent if they are equal almost everywhere. Definition 2.6. Let u L 2 () and α N d 0. A function w L 2 () is called the weak derivative or distributional derivative α u (of order α ) of u, if wv dx = ( 1) α u α v dx v C0 (). Based on this definition, all linear differential operators introduced in Section 2.2 can be given a weak/distributional interpretation. For example, the weak gradient u of a function u L 2 () will be a vectorfield w (L 2 ()) d with w, v dx = u v dx v (C0 ()) d. (2.3) This can be directly concluded from (FGF). The same is true of the weak divergence w L 2 () of a vectorfield u (L 2 ()) d w v dx = u, v dx v C0 (). (2.4) Theorem 2.7. If u C m (), then all weak derivatives of order m agree in L 2 () with the corresponding classical derivatives. Proof. Clear by a straightforward application of (2.1). Hence, without changing notations, all derivatives will be understood as weak derivatives in the sequel. Straight from the definition we also infer that all linear differential operators in weak sense commute. Γ Γ 13
4 Extension of integration by parts formula For Lipschitz domains the Gauß theorem (Theorem 2.3) and hence, also the first Green s formula (FGF) can be extended to integrable functions with weak derivatives. Weak derivative in subdomains The weak derivative in and a subdomain 1 coincide in 1 since the corresponding spaces of test functions satisfy C 0 ( 1 ) C 0 (). Continuity over interfaces of subdomains For the partition = 1 2, 1 2 =, where both subdomains are supposed to have a Lipschitz boundary, a function u, which possesses a weak gradient in 1 and 2, possesses a weak gradient in if and only if [u] Σ = 0 on the common interface Σ := 1 2. Proof. The only candidate for weak gradient in can only be the piecewise weak gradient in 1 and 2. Let n Σ denote the normal vector on Σ pointing outwards 1. We test with v (C0 ()) d and so u v dx = u v dx 1 2 = u v dx + (u 1 v n Σ u 2 v n Σ )ds(x) 1 2 Σ = u v dx, if and only if Σ [u] Σ vds(x) = 0 v C 0 (Σ), which is by Lemma 2.5 equivalent to [u] Σ = 0 almost everywhere. Remark 2.8. The fact that a function has classical derivatives in subdomains and even the fact that the classical derivatives can be extended to a function in, which is in L 2 (), does not guarantee that it has a weak derivative in the whole domain. For example a function which is piecewise C ( l ), l = 1, 2, but discontinuous over Σ := 1 2, does not have a weak gradient. 2.5 Variational formulations of elliptic boundary value problems Dirichlet boundary conditions Interpreting the divergence in (1.4a) as weak derivative we can equivalently write a u, v + c uv dx = fv dx (2.5) for test functions v C 0 (). Additionally the solution u has to fulfill the Dirichlet boundary conditions (1.4b), which all trial functions functions to try if they solve (2.5) have to fulfill. Boundary conditions which are present in constraint to the space of trial or test functions are called essential Neumann boundary conditions We enrich the space of test functions and choose v C () C(). Using (FGF) we get a u, v + c uv dx a u, n v ds = fv dx (2.6) for all v C () C(). We can incorporate the Neumann boundary condition (1.4c) in the variational formulation. For pure Neumann boundary conditions, i. e., Γ N =, this reads a u, v + c uv dx = fv dx + h v ds (2.7) 14
5 for all v C () C(). Boundary conditions which are present in the variational formulation are called natural. To see that (2.7) is equivalent to (1.4a) and (1.4c) we choose first test functions v C0 (). Then, the term on disappears, and with the definition of the weak divergence (2.4) we get (1.4a). Now, let v C () C 0 ( ). Integrating by parts and using the fact that (1.4a) holds we get a u, n v ds = hvds which is by Lemma 2.5 equivalent to (1.4c) Robin boundary conditions For the case of pure Robin boundary conditions, where Neumann boundary conditions are included with β = 0, inserting (1.4d) into (2.6) leads to a u, v + c uv dx + βu v ds = fv dx + h v ds. (2.8) Also Robin boundary conditions are natural. 2.6 Linear and bilinear forms Let V, W be real/complex vector spaces. Linear operator is a mapping of the form T : V W, that satisfies T(λv + µw) = λ T v + µ T w for all v, w V, λ, µ R/C. Linear form is a linear operator of the form l : V R/C. Bilinear form is a mapping of the form b : V V R/C, where v b(v, w) for every w and w b(v, w) for every v are both linear forms on V, e. g. a u, v + c uv dx = fv dx } {{ } }{{ } b(u,v) l(v) Antilinear operator is a mapping of the form T : V W that satisfies T(λv + µw) = λ T v + µ T w for all v, w V, λ, µ R/C. Antilinear form is an antilinear operator of the from l : V R/C. Sesquilinear form is a mapping of the form b : V V R/C, where v b(v, w) is a linear form for every w and w b(v, w) is an antilinear form for every v, e. g. a u, v + c uv dx = fv dx } {{ } }{{ } b(u,v) l(v) Linear variational problem Given some linear form l seek u V such that b(u, v) = l(v) v V. (LVP) Symmetric bi-/sesquilinear form if b(v, w) = b(w, v) for all v, w V. Positive definite bi-/sesquilinear form if for all v V we have b(v, v) > 0 v 0. 15
6 V -continuous bi-/sesquilinear form satisfies b(u, v) b u V v V for all u, v V with continuity constant b := sup u,v V \{0} u V v V b(u,v). (This property is also called boundedness, i. e., b is said to be bounded in V with constant b.) V -elliptic bi-/sesquilinear form satisfies b(v, v) γ v 2 V for all v V with ellipticity constant γ. (This property is also called coercivity.) A symmetric positive definite bi-/sesquilinear form b is an inner product that induces a norm through The fundamental Cauchy-Schwarz inequality v b := b(v, v) 1 2 v V. b(v, w) v b w b v, w V ensures that b is continuous with continuity constant 1 with respect to its induced norm. Moreover, we have Pythagoras theorem a(v, w) = 0 v 2 a + w 2 a = v + w 2 a. In the context of elliptic partial differential equations a norm that can be derived from a V -elliptic bi-/sesquilinear form is often dubbed energy norm, denoted by e. 2.7 Sobolev spaces Definition 2.9. A normed vector space V is complete, if every Cauchy sequence {v k } k V has a limit v in V. A complete normed vector space is called a Banach space. Definition A Hilbert space is a Banach space whose norm is induced by an inner product. In Section 2.5 we learned that the formal variational problem associated with the pure homogeneous Neumann problem for (1.4) is: seek u : R such that a u, v + c uv dx = fv dx (2.9) for all test functions v. Ideal space of u and v in (2.9) is a Hilbert space H whose inner product coincides with the bi-/sesquilinear form, and consequently, which is equipped with the energy norm v 2 e := a v 2 + c v 2 dx and the space is defined as H := {v : R : weak gradient v exists and v e < }. The existence of a weak gradient u is built in the space, so we restrict the search only to those functions. The same space applies for test functions, that are now more functions than C () (equivalence due to density by Meyers-Serrin theorem, see Theorem 2.15). Definition 2.11 (Sobolev space H 1 ()). The Sobolev space H 1 () is the space of square integrable functions R with square integrable weak gradients with norm H 1 () := {v L 2 () : (the weak gradient) v exists and v L 2 ()} v 2 H 1 () := v 2 L 2 () + D2 v 2 L 2 () where D = diam() is introduced to match units, which is often omitted in the definition. 16
7 Note, that the H 1 ()-norm and the energy norm are equivalent if the latter is based on an H 1 ()- elliptic and H 1 ()-continuous bilinear form, i. e., there exist two constants C 1, C 2 > 0 such that C 1 v e v H 1 () C 2 v e v H 1 (). Definition For m N 0 and R d we define the Sobolev space of order m as H m () := {v L 2 () : α v L 2 (), α m}, equipped with the norm v Hm () := ( α m α v 2 L 2 () ) 1 2. Notation For all m N 0 and R d denotes a semi-norm on H m (). v H m () := ( α =m α v 2 L 2 () ) 1 2 Theorem The Sobolev spaces H m (), m N 0, are Hilbert spaces with the inner product (u, v) Hm () := ( α u, α v) L2 () u, v H m (). α m Theorem 2.15 (Meyers-Serrin theorem). The space C () is a dense subspace of H m () for all m N 0. Definition H 1 0 () is defined as the completion of C 0 () with respect to the H 1 ()-norm. 2.8 Theory of variational formulations Ellipticity and the lemma of Lax-Milgram Definition 2.17 (Ellipticity of bilinear forms). A bilinear form b(, ) is called V -elliptic, if there exists a γ > 0 such that for all v V b(v, v) γ v 2 V. Definition 2.18 (Ellipticity of sesquilinear forms). A sesquilinear form b(, ) is called V -elliptic, if there exists a γ > 0 and a θ C with θ = 1, such that for all v V Re(θ b(v, v) γ v 2 V. we rotate the sesquilinear form by θ to have positive real part Definition The dual V of a normed vector space V is the normed vector space of continuous linear forms on V. The dual space is equipped with the (operator) norm f V = sup 0 v V f(v) v V. Example Linear forms fvdx with f L2 () are in the dual of H 1 () since by the Cauchy- Schwarz inequality sup fvdx f L sup 2 () v L 2 () f L 0 v H 1 () v H1 () 0 v H 1 () v 2 (). H1 () Lemma 2.21 (Lax-Milgram). Let V be a reflexive Banach space. Let the bilinear form b : V V R or sesquilinear form b : V V C be bounded and V -elliptic. Then, the variational problem (LVP) has for any f V a unique solution u V with u V 1 γ f V. 17
8 Proof. The V -ellipticity of the sesquilinear form b implies for u V γ u 2 V Re(θb(u, u)) θb(u, u) = b(u, u) = f(u) f V u V, and so the stability estimate and the uniqueness (f = 0 u = 0) follow. The existence follows by the Riesz representation theorem where the bi-/sesquilinear form is taken as inner product. Theorem 2.22 (Riesz representation theorem). Let H be a Hilbert space and ϕ H an arbitrary linear form on H. Then there exists a unique element u H such that Moreover, ϕ H = u H The inf-sup conditions ϕ(v) = (u, v) H for all v H. A more general result on the existence and uniqueness of weak solutions is provided by the inf-sup conditions. Theorem Let U, V be a reflexive Banach spaces with norms U and V. Furthermore let the bi-/sesquilinear form b be bounded in U, V. Then the following statements are equivalent: (i) For all f V the linear variational problem (LVP) has a unique solution u f U that satisfies with γ > 0 independent of f. u f U 1 γ f V, (ii) The bilinear form b satisfies the inf-sup conditions γ > 0 : inf sup w U\{0} v V \{0} b(w, v) v V w U γ, (IS1) v V \ {0} : sup b(w, v) > 0. (IS2) w U\{0} 2.9 Well-posedness of elliptic boundary value problems Boundary value problems with Dirichlet data For the variational formulation of the PDE (1.4) the Dirichlet boundary conditions are essential and are considered in the space of trial functions H 1 Γ D,g() := {v H 1 () : v = g on Γ D } H 1 (), which is equipped with the (usual) H 1 ()-norm. Note that HΓ 1 D,g() is not a vector space, i. e., if g 0, u HΓ 1 D,g () and λ 1 then λu H1 Γ D,g (). For the test functions we use the associated homogeneous space H 1 Γ D,0 (). If Γ D = we can write shortly H 1 g () and H 1 0 (). Variational formulation for pure Dirichlet boundary conditions: Seek u Hg 1 () such that a(x) u(x) v(x) + cu(x)v(x)dx = f(x)v(x)dx v H0 1 (). It is equivalent to write u = u 0 + u g with u g Hg 1 () some extension of g into and the new unknown (offset) u 0 H0 1 (): Seek u 0 H0 1 () such that b(u 0, v) = a u 0 v + cu 0 vdx = fv a u g v + cu g vdx = l(v) v H0 1 (). For well-posedness it is enough to prove H 1 0 ()-ellipticity b(v, v) γ v H1 () v H 1 0 (), which is a weaker statement than H 1 ()-ellipticity. 18
9 Case c c 0 > 0 Since by assumption a a 0 > 0, we have b(v, v) min{a 0, c 0 } v H 1 () v H 1 (), from which H 1 0 ()-ellipticity directly follows. Case c = 0 Now we only have b(v, v) a 0 v L2 () v H 1 (). However, we shall apply the following result to obtain H 1 Γ D,0 ()-ellipticity. Lemma 2.24 (Friedrichs inequality). Let be an open bounded domain and Γ D > 0. Then, for all u H 1 Γ D,0 () there exists a constant C F (, Γ D ) such that Now we can conclude for all v H 1 Γ D,0 (). u 2 L 2 () C F (, Γ D ) diam() 2 u 2 L 2 (). b(v, v) a 0 v 2 L 2 () a 0 2 v 2 L 2 () + a 0 2C F diam() 2 v 2 L 2 () a { } 0 2 min 1 1, C F diam() 2 v 2 H 1 () For which functions g can we find an extension u g H 1 ()? Any restriction trace of a function v H 1 () to Γ D can be extended to v H 1 (). We can choose u g to be any of the functions v for which v = g on Γ D. For Γ D = the trace space of H 1 () is H 1 /2 ( ) L 2 ( ) with the norm v H 1 /2 ( ) = inf{ w H 1 () : w C (), w Γ = v}. The index 1 /2 corresponds to the definition of Sobolev spaces with rational indices (as interpolation spaces), where a direct definition of the norm exists (Sobolev-Slobodeckij norm). For all functions in g H 1 /2 ( ) an extension u g H 1 () exists, but not for any g L 2 ( ) (e. g., discontinuous functions). Theorem The trace operator R 1 : H 1 () H 1 /2 (Γ) is continuous and surjective (it is onto) and has a bounded right inverse (extension operator) E 1 : H 1 /2 (Γ) H 1 (), i. e., R 1 E 1 = Id Boundary value problems with Neumann data For the variational formulation of the boundary value problem (1.4) the Neumann boundary conditions are natural and we use in case of pure Neumann boundary conditions H 1 () as space of trial and test functions: Seek u H 1 () such that b(u, v) = a(x) u(x) v(x) + c(x)u(x)v(x)dx = f(x)v(x)dx + h(x)v(x)ds(x) = l(v) v H 1 (). (2.10) Analogously to the discussion of boundary value problems with Dirichlet data in Section we have H 1 ()-ellipticity if a a 0 > 0 and c c 0 > 0. 19
10 Case c 0 Let us test (2.10) with v = 1. Then b(u, 1) = 0 for all u H 1 (), and hence, l(1) = 0 must hold as well. This leads us to the compatibility condition for the data of the problem f dx + h ds = 0. (2.11) The solution is not unique. Since b(1, v) = 0 for any v H 1 () we can add any constant. The constant can be fixed by demanding a vanishing mean value of u, i. e., we switch to the subspace H () 1 := {v H 1 () : v dx = 0}. The variational formulation is: Seek u H 1 () such that b(u, v) = l(v) v H 1 (). The bilinear form b is H 1 ()-elliptic. For the proof we apply Lemma 2.26 (Poincaré s inequality). Let be a open bounded domain. Then, for all u H 1 () there exists a constant C P () such that where u = 1 u(x)dx is the mean of u in. u u 2 L 2 () C P () diam() 2 u 2 L 2 (), Case c c 0 > 0 in some subdomain c There is no need of a compatibility condition (2.11) since there exist functions u H 1 () such that b(u, 1) = c c(x)u(x)dx 0 The constant is also fixed by the term c c(x)u(x)v(x)dx as we can apply Lemma Let f : H 1 () R be a continuous quadratic functional which fixes constants and for which we have f(λv) = λ 2 f(v) for λ R, i. e., any constant function vanishes if f(1) = 0. Then, it exists a constant C = C(f, ) > 0 such that for all u H 1 () we have ( ) u 2 L 2 () C u 2 L 2 () + f(u). Since c c 0 > 0 in c, the (quadratic) functional c c(x) u(x) 2 dx fixes constants. Then, we have H 1 ()-ellipticity, since { b(v, v) a 0 v 2 L 2 () + c(x) v(x) 2 dx min 1, a } ( ) 0 2 v 2 L c 2 2 () + c(x) v(x) 2 dx c { min 1, a } ( 0 v 2 L 2 2 () + 1 C v 2 L ()) γ v 2 2 H 1 () for all v H 1 () where γ = min { 1, 1 C, a0 2, a0 2C } Boundary value problems with Robin data Variational formulation for pure Robin boundary conditions: Seek u H 1 () such that b(u, v) = a(x) u(x) v(x) + c(x)u(x)v(x)dx + = f(x)v(x)dx β(x)u(x)v(x)ds(x) h(x)v(x)ds(x) = l(v) v H 1 (). Analogously to above, we can deduce H 1 ()-ellipticity if a a 0 > 0 and c c 0 > 0 (like for Dirichlet boundary conditions) or if a a 0 > 0, c 0 and β β 0 > 0 (like for Neumann boundary conditions with c 0 in c ). 20
11 2.10 Discrete variational formulations A first step towards finding a practical algorithm for the approximate solution of (LVP) is to convert it into a discrete variational problem. The attribute discrete means that the solution can be characterized by a finite number of real (or complex) numbers. Galerkin discretization of (LVP) We replace V, e. g., H 1 () or H 1 0 (), in (LVP) by finite dimensional subspaces equipped with the same norm. The most general approach relies on two subspaces of V, i. e., If W n V : trial space, dim W n = N N, and V n V : test space, dim V n = N N. W n = V n we speak of a classical Galerkin discretization, if furthermore the bilinear form provides an inner product, we call it Ritz-Galerkin discretization, W n V n the discretization is called Petrov-Galerkin discretization. The subscript n labels discrete entities, but it can also deal as index referring to sequences of trial and test spaces in a convergence study. Discrete variational problem corresponding to (LVP): Seek u n W n such that b(u n, v n ) = l(v n ) v n V n. (DVP) Galerkin orthogonality of the discretization error e n := u u n to the test space V n We can insert v n V n into (DVP) and (LVP), the difference gives b(u u n, v n ) = 0 v n V n. The error is the smallest possible when measured in the energy norm (b-optimality). u e n V 0 u n V n Figure 2.3: b-orthogonality of the error e n = u u n with respect to V n, if b is an inner product. Well-posedness of the discretized problem For a classical Galerkin discretization with V n V V n -ellipticity follows directly from the V - ellipticity with the same constant γ and hence, by the Lax-Milgram lemma, existence and uniqueness follow. The error can be estimated by Lemma 2.28 (Céa s lemma). Let the bi-/sesquilinear form b be V -elliptic with constant γ and V - continuous with constant b. Then, the solutions u of (LVP) and u n of (DVP) satisfy the quasioptimality estimate u u n V b γ inf u v n V. v n V n 21
12 Proof. Using V -ellipticity, Galerkin orthogonality and continuity of b we have where v n V n was arbitrarily chosen. γ u u n 2 V b(u u n, u u n ) = b(u u n, u v n + v n u n ) b(u u n, u v n ) + b(u u n, v n u n ) }{{} =0 by Galerkin orthog. b u u n V u v n V For the energy norm we have γ = b = 1 and thus what is known as b-optimality. u u h b inf v n V n u v n b For general bi-/sesquilinear forms we need besides (IS1), (IS2) that the discrete inf-sup-condition is fulfilled γ n > 0 : inf sup w n W n v n V n\{0} b(w n, v n ) v n V w n V γ n. (DIS) One cannot conclude (DIS) from (IS1) and (IS2) because the supremum is taken over a much smaller set. If the constant γ n is uniformly bounded away from zero an asymptotic quasi-optimality of the Galerkin solution is guaranteed ( u u n V 1 + b ) V V R inf u w n γ n w V, n W n which is automatically satisfied in case of ellipticity, see Céa s lemma (Lemma 2.28) Algebraic setting (DVP) can be solved on a computer by introducing ordered bases B V := {p 1 n,..., p N n } of V n, B W := {q 1 n,..., q N n } of W n, N := dim V n = dim W n. Linear system of equations Bu = f (LSE) with B := ( b(q k n, p j n) ) N j,k=1 RN,N, f := ( l(p k n) ) N k=1 RN. Unknown are the coefficients u (coefficient vector) in the basis representation. Equivalence: (LSE) has a unique solution if and only if (DVP) has a unique solution. j is the row index, whereas k is the column index: element b(q k n, p j n) of the matrix B is highlighted: j k 22
13 (Continuous) variational problem (LVP): b(u, v) = l(v) Choose finite dimensional trial and test spaces W n and V n. Discrete variational problem (DVP): b(u n, v n ) = l(v) Choose bases of W n and V n Algebraic problem (LSE): Bu = f Figure 2.4: Overview of stages involved in the Galerkin discretization of an abstract variational problem. 23
14 24
Finite 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 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 informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More 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 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 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 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 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 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 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 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 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 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 informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
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 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 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 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 information2.3 Variational form of boundary value problems
2.3. VARIATIONAL FORM OF BOUNDARY VALUE PROBLEMS 21 2.3 Variational form of boundary value problems Let X be a separable Hilbert space with an inner product (, ) and norm. We identify X with its dual X.
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More 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 informationFriedrich symmetric systems
viii CHAPTER 8 Friedrich symmetric systems In this chapter, we describe a theory due to Friedrich [13] for positive symmetric systems, which gives the existence and uniqueness of weak solutions of boundary
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 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 informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
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 informationt y n (s) ds. t y(s) ds, x(t) = x(0) +
1 Appendix Definition (Closed Linear Operator) (1) The graph G(T ) of a linear operator T on the domain D(T ) X into Y is the set (x, T x) : x D(T )} in the product space X Y. Then T is closed if its graph
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 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 information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationCSE386C METHODS OF APPLIED MATHEMATICS Fall 2014, Final Exam, 9:00-noon, Tue, Dec 16, ACES 6.304
CSE386C METHODS OF APPLIED MATHEMATICS Fall 2014, Final Exam, 9:00-noon, Tue, Dec 16, ACES 6.304 1. Let A be a closed operator defined on a dense subspace D(A) of a Hilbert space U. (a) Define the adjoint
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 informationThe following definition is fundamental.
1. Some Basics from Linear Algebra With these notes, I will try and clarify certain topics that I only quickly mention in class. First and foremost, I will assume that you are familiar with many basic
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 informationVariational Principles for Equilibrium Physical Systems
Variational Principles for Equilibrium Physical Systems 1. Variational Principles One way of deriving the governing equations for a physical system is the express the relevant conservation statements and
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 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 informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
More informationDIFFERENTIAL EQUATIONS
NUMERICAL PARTIAL DIFFERENTIAL EQUATIONS Lecture Notes, Winter 2013/14 Christian Clason July 15, 2014 Institute for Mathematics and Scientific Computing Karl-Franzens-Universität Graz CONTENTS I BACKGROUND
More 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 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 informationInequalities of Babuška-Aziz and Friedrichs-Velte for differential forms
Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms Martin Costabel Abstract. For sufficiently smooth bounded plane domains, the equivalence between the inequalities of Babuška Aziz
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 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 informationEigenvalues and Eigenfunctions of the Laplacian
The Waterloo Mathematics Review 23 Eigenvalues and Eigenfunctions of the Laplacian Mihai Nica University of Waterloo mcnica@uwaterloo.ca Abstract: The problem of determining the eigenvalues and eigenvectors
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationApproximation in Banach Spaces by Galerkin Methods
2 Approximation in Banach Spaces by Galerkin Methods In this chapter, we consider an abstract linear problem which serves as a generic model for engineering applications. Our first goal is to specify the
More 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 informationc 2004 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 64, No. 3, pp. 902 931 c 2004 Society for Industrial and Applied Mathematics COMPLETE ELECTRODE MODEL OF ELECTRICAL IMPEDANCE TOMOGRAPHY: APPROXIMATION PROPERTIES AND CHARACTERIZATION
More informationThe Dirichlet-to-Neumann operator
Lecture 8 The Dirichlet-to-Neumann operator The Dirichlet-to-Neumann operator plays an important role in the theory of inverse problems. In fact, from measurements of electrical currents at the surface
More informationTechnische Universität Graz
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 Technische
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 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 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 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 information1. Bounded linear maps. A linear map T : E F of real Banach
DIFFERENTIABLE MAPS 1. Bounded linear maps. A linear map T : E F of real Banach spaces E, F is bounded if M > 0 so that for all v E: T v M v. If v r T v C for some positive constants r, C, then T is bounded:
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 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. Preliminaries. The purpose of this chapter is to provide some basic background information. Linear Space. Hilbert Space.
Chapter 1 Preliminaries The purpose of this chapter is to provide some basic background information. Linear Space Hilbert Space Basic Principles 1 2 Preliminaries Linear Space The notion of linear space
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 informationProblem Set 6: Solutions Math 201A: Fall a n x n,
Problem Set 6: Solutions Math 201A: Fall 2016 Problem 1. Is (x n ) n=0 a Schauder basis of C([0, 1])? No. If f(x) = a n x n, n=0 where the series converges uniformly on [0, 1], then f has a power series
More informationANALYSIS HW 5 CLAY SHONKWILER
ANALYSIS HW 5 CLAY SHONKWILER Let X be a normed linear space and Y a linear subspace. The set of all continuous linear functionals on X that are zero on Y is called the annihilator of Y and denoted by
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 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 informationA WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS
A WEAK GALERKIN MIXED FINITE ELEMENT METHOD FOR BIHARMONIC EQUATIONS LIN MU, JUNPING WANG, YANQIU WANG, AND XIU YE Abstract. This article introduces and analyzes a weak Galerkin mixed finite element method
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationCategories and Quantum Informatics: Hilbert spaces
Categories and Quantum Informatics: Hilbert spaces Chris Heunen Spring 2018 We introduce our main example category Hilb by recalling in some detail the mathematical formalism that underlies quantum theory:
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 informationSobolevology. 1. Definitions and Notation. When α = 1 this seminorm is the same as the Lipschitz constant of the function f. 2.
Sobolevology 1. Definitions and Notation 1.1. The domain. is an open subset of R n. 1.2. Hölder seminorm. For α (, 1] the Hölder seminorm of exponent α of a function is given by f(x) f(y) [f] α = sup x
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 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 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 REVIEW OF OPTIMIZATION
1 OPTIMAL DESIGN OF STRUCTURES (MAP 562) G. ALLAIRE December 17th, 2014 Department of Applied Mathematics, Ecole Polytechnique CHAPTER III A REVIEW OF OPTIMIZATION 2 DEFINITIONS Let V be a Banach space,
More 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 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 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 informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
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 informationand finally, any second order divergence form elliptic operator
Supporting Information: Mathematical proofs Preliminaries Let be an arbitrary bounded open set in R n and let L be any elliptic differential operator associated to a symmetric positive bilinear form B
More informationChapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method
Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 On the divergence constraint in mixed finite element methods for incompressible
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 information16 1 Basic Facts from Functional Analysis and Banach Lattices
16 1 Basic Facts from Functional Analysis and Banach Lattices 1.2.3 Banach Steinhaus Theorem Another fundamental theorem of functional analysis is the Banach Steinhaus theorem, or the Uniform Boundedness
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 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 informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More information1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.
1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the
More informationPartial Differential Equations and Sobolev Spaces MAT-INF4300 autumn Snorre H. Christiansen November 10, 2016
Partial Differential Equations and Sobolev Spaces MAT-INF4300 autumn 2016 Snorre H. Christiansen November 10, 2016 1 2 Contents 1 Introduction 5 2 Notations for differential operators 7 3 Generalities
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 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 informationRegularity Theory a Fourth Order PDE with Delta Right Hand Side
Regularity Theory a Fourth Order PDE with Delta Right Hand Side Graham Hobbs Applied PDEs Seminar, 29th October 2013 Contents Problem and Weak Formulation Example - The Biharmonic Problem Regularity Theory
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 informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationPROBLEMS. (b) (Polarization Identity) Show that in any inner product space
1 Professor Carl Cowen Math 54600 Fall 09 PROBLEMS 1. (Geometry in Inner Product Spaces) (a) (Parallelogram Law) Show that in any inner product space x + y 2 + x y 2 = 2( x 2 + y 2 ). (b) (Polarization
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 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 information