DOCSLIB.ORG

Arxiv:1605.07198V1 [Math.NA] 23 May 2016 Htaerbs in Robust Are That H Xeso Prtri Iia Otedrcdlafnto.W Note Trace We a Function

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Arxiv:1605.07198V1 [Math.NA] 23 May 2016 Htaerbs in Robust Are That H Xeso Prtri Iia Otedrcdlafnto.W Note Trace We a Function EFFICIENT PRECONDITIONERS FOR SADDLE POINT SYSTEMS WITH TRACE CONSTRAINTS COUPLING 2D AND 1D DOMAINS ∗ MIROSLAV KUCHTA †, MAGNE NORDAAS ‡, JORIS C. G. VERSCHAEVE †, MIKAEL MORTENSEN †‡, AND KENT-ANDRE MARDAL †‡ Abstract. We study preconditioners for a model problem describing the coupling of two ellip- tic subproblems posed over domains with different topological dimension by a parameter dependent constraint. A pair of parameter robust and efficient preconditioners is proposed and analyzed. Ro- bustness and efficiency of the preconditioners is demonstrated by numerical experiments. Key words. preconditioning, saddle-point problem, Lagrange multipliers AMS subject classifications. 65F08 1. Introduction. This paper is concerned with preconditioning of multiphysics problems where two subproblems of different dimensionality are coupled. We assume that Γ is a sub-manifold contained within Ω Rn and consider the following problem: ∈ ∆u + ǫδ p = f in Ω, (1.1a) − Γ ∆v p = g on Γ, (1.1b) − − ǫu v =0 onΓ, (1.1c) − where δΓ is a function with properties similar to the Dirac delta function as will be discussed later. To allow for a unique solution (u,v,p) the system must be equipped with suitable boundary conditions and we shall here, for simplicity, consider homoge- neous Dirichlet boundary conditions for u and v on ∂Ω and ∂Γ respectively. We note that the unknowns u, v are here the primary variables, while the unknown p should be interpreted as a Lagrange multiplier associated with the constraint (1.1c). The two elliptic equations that are stated on two different domains, Ω and Γ, are coupled and therefore the restriction of u to Γ and the extension of p to Ω are crucial. When the codimension of Γ is one, the restriction operator is a trace operator and the extension operator is similar to the Dirac delta function. We note that ǫ (0, 1) and that the typical scenario will be that ǫ 1. We will therefore focus on methods∈ that are robust in ǫ. ≪ The problem (1.1a)–(1.1c) is relevant to biomedical applications [18, 15, 2, 17] arXiv:1605.07198v1 [math.NA] 23 May 2016 where it models the coupling of the porous media flow inside tissue to the vascular bed through Starlings law. Further, problems involving coupling of the finite element method and the boundary element method, e.g. [24, 26], are of the form (1.1). The system is also relevant for domain decomposition methods based on Lagrange multi- pliers [32]. Finally, in solid mechanics, the problem of plates reinforced with ribs, cf. for example [44, ch. 9.11], can be recast into a related fourth order problem. We also †Department of Mathematics, Division of Mechanics, University of Oslo {mirok, joris, mikaem}@math.uio.no ‡Center for Biomedical Computing, Simula Research Laboratory, {magneano, kent-and}@simula.no ∗The work of Magne Nordaas, Mikael Mortensen and Kent-Andre Mardal has been supported by a Center of Excellence grant from the Research Council of Norway to the Center for Biomedical Computing at Simula Research Laboratory. 1 2 Preconditioning for trace constrained systems note that the techniques developed here to address the constraint (1.1c) are applica- ble in preconditioning fluid-structure interaction problems involving interactions with thin structures, e.g. filaments [22]. One way of deriving equations (1.1) is to consider the following minimization problem ( u)2 2uf dx ∇ − ZΩ min (1.2) 2 → ( v) 2vg ds ∇ − ZΓ subject to the constraint ǫu v =0 onΓ. (1.3a) − Using the method of Lagrange multipliers, the constrained minimization problem will be re-cast as a saddle-point problem. The saddle-point problem is then analyzed in terms of the Brezzi conditions [13] and efficient solution algorithms are obtained using operator preconditioning [35]. A main challenge is the fact that the constraint (1.3a) necessitates the use of trace operators which leads to operators in fractional Sobolev spaces on Γ. An outline of the paper is as follows: Section 2 presents the necessary notation and mathematical framework needed for the analysis. Then the mathemathical analysis as well as the numerical experiments of two different preconditioners are presented in §3 and §4, respectively. Section 5 discusses computational efficiency of both methods. 2. Preliminaries. Let X be a Hilbert space of functions defined on a domain D 2 and let X denote its norm. The L inner product on a domain D is denoted ( , )D k·k · · or D , while , D denotes the corresponding duality pairing between a Hilbert · h· ·i m m space X and its dual space X∗. We will use H = H (D) to denote the Sobolev spaceR of functions on D with m derivatives in L2 = L2(D). The corresponding norm m m is denoted m,D. In general, we will use H0 to denote the closure in H of the space of smoothk·k functions with compact support in D and seminorm is denoted as m,D. | · | The space of bounded linear operators mapping elements of X to Y is denoted (X, Y ) and if Y = X we simply write (X) instead of (X,X). If X and Y are HilbertL spaces, both continuously containedL in some largerL Hilbert space, then the intersection X Y and the sum X + Y are both Hilbert spaces with norms given by ∩ 2 2 2 2 2 2 x X Y = x X + x Y and z X+Y = inf ( x X + y Y ). k k ∩ k k k k k k x X,y Y k k k k ∈z=x+∈y In the following Ω Rn is an open connected domain with Lipschitz boundary ⊂ ∂Ω. The trace operator T is defined by Tu = u Γ for u C(Ω) and Γ a Lipschitz submanifold of codimension one in Ω. The trace| operator∈ extends to bounded and 1 1 surjective linear operator T : H (Ω) H 2 (Γ), see e.g. [1, ch. 7]. The fractional 1 → Sobolev space H 2 (Γ) can be equipped with the norm 2 2 2 u(x) u(y) u 1 = u 2 + | − | dxdy. (2.1) H 2 (Γ) L (Γ) n+1 k k k k Γ Γ x y Z × | − | 1 1 2 However, the trace is not surjective as an operator from H0 (Ω) into H (Γ), in 1 particular the constant function 1 H 2 (Γ) is not in the image of the trace operator. ∈ 3 1 1 1 2 2 2 Note that H0(Γ) does not characterize the trace space, since H0(Γ) = H (Γ), see [30, 1 2 ch. 2, thm. 11.1]. Instead, the trace space can be identified as H00(Γ), defined as 1 1 the subspace of H 2 (Γ) for which extension by zero into H 2 (Γ)˜ is continuous, for some suitable extension domain Γ˜ extending Γ (e.g. Γ=Γ˜ ∂Ω). To be precise, the space 1 2 ∪ H00(Γ) can be characterized with the norm u(x) x Γ 1 1 u 2 = u˜ H 2 (Γ)˜ , u˜(x)= ∈ (2.2) k kH00(Γ) k k 0 x / Γ. ( ∈ 1 2 ˜ The space H00(Γ) does not depend on the extension domain Γ, since the norms induced by different choices of Γ˜ will be equivalent. The above norms (2.1)–(2.2) for the fractional spaces are impractical from an implementation point of view, and we will therefore consider the alternative construc- 1 tion following [30, ch. 2.1] and [16]. For u, v H0 (Γ), set Lu(v) = (u, v)Γ. Then 1 ∈ Lu is a bounded linear functional on H0 (Γ) and in accordance with the Riesz-Fr´echet theorem there is an operator S H1(Γ) such that ∈ L 0 1 (Su,w) 1 = Lu(w) = (u, w) , u,w H (Γ). (2.3) H0 (Ω) Γ ∈ 0 The operator S is self-adjoint, positive definite, injective and compact. Therefore the spectrum of S consists of a nonincreasing sequence of positive eigenvalues λk ∞ { }k=1 such that 0 < λk λk and λk 0, see e.g. [48, ch. X.5, thm. 2]. The eigenvectors +1 ≤ → φk ∞ of S satisfy the generalized eigenvalue problem { }k=1 1 Aφk = λk− Mφk where operators A, M are such that Au, v Γ = ( u, v)Γ and Mu,v Γ = (u, v)Γ. h i ∇1 ∇ h i The set of eigenvectors φk k∞=1 forms a basis of H0 (Γ) orthogonal with respect the 1 { } 2 inner product of H0 (Γ) and orthonormal with respect to the inner product on L (Γ). Then for u = ckφk span φk ∞ and s [ 1, 1], we set k ∈ { }k=1 ∈ − P 2 s u H = c λ− (2.4) k k s k k k sX 2 and define Hs to be the closure of span φk k∞=1 in the above norm. Then H0 = L (Γ) 1 1 { } 2 1 and H1 = H0 (Γ), with equality of norms. Moreover, we have H 2 = H00(Γ) with 1 2 1 equivalence of norms. This essentially follows from the fact that H 2 and H00(Γ) are closely related interpolation spaces, see [16, thm. 3.4]. Note that we also have 1 1 1 1 2 1 2 H 1 = (H0 (Γ))∗ = H− (Γ) and H 2 = (H00(Γ))∗ = H− (Γ). − As the preceeding paragraph suggests− we shall use normal font to denote linear operators, e.g. A. To signify that the particular operator acts on a vector space with multiple components we employ calligraphic font, e.g. Vectors and matrices are denoted by the sans serif font, e.g., A and x. In case theA matrix has a block structure it is typeset with the blackboard bold font, e.g.
Load more

Recommended publications
  • L P and Sobolev Spaces
    NOTES ON Lp AND SOBOLEV SPACES STEVE SHKOLLER 1. Lp spaces 1.1. Definitions and basic properties. Definition 1.1. Let 0 < p < 1 and let (X; M; µ) denote a measure space. If f : X ! R is a measurable function, then we define 1 Z p p kfkLp(X) := jfj dx and kfkL1(X) := ess supx2X jf(x)j : X Note that kfkLp(X) may take the value 1. Definition 1.2. The space Lp(X) is the set p L (X) = ff : X ! R j kfkLp(X) < 1g : The space Lp(X) satisfies the following vector space properties: (1) For each α 2 R, if f 2 Lp(X) then αf 2 Lp(X); (2) If f; g 2 Lp(X), then jf + gjp ≤ 2p−1(jfjp + jgjp) ; so that f + g 2 Lp(X). (3) The triangle inequality is valid if p ≥ 1. The most interesting cases are p = 1; 2; 1, while all of the Lp arise often in nonlinear estimates. Definition 1.3. The space lp, called \little Lp", will be useful when we introduce Sobolev spaces on the torus and the Fourier series. For 1 ≤ p < 1, we set ( 1 ) p 1 X p l = fxngn=1 j jxnj < 1 : n=1 1.2. Basic inequalities. Lemma 1.4. For λ 2 (0; 1), xλ ≤ (1 − λ) + λx. Proof. Set f(x) = (1 − λ) + λx − xλ; hence, f 0(x) = λ − λxλ−1 = 0 if and only if λ(1 − xλ−1) = 0 so that x = 1 is the critical point of f. In particular, the minimum occurs at x = 1 with value f(1) = 0 ≤ (1 − λ) + λx − xλ : Lemma 1.5.
    [Show full text]
  • Degenerate Parabolic Initial-Boundary Value Problems*
    JOURNAL OF DIFFERENTIALEQUATIONS 31,296-312 (1979) Degenerate Parabolic Initial-Boundary Value Problems* R. E. SHOWALTER Department of Mathematics, RLM 8.100, The University of Texas, Austin, Texas 78712 Received September 22, 1977 l. INTRODUCTION We consider a class of implicit linear evolution equations of the form d d-t J[u(t) + 5flu(t) = f(t), t > O, (1.1) in Hilbert space and their realizations in function spaces as initial-boundary value problems for partial differential equations which may contain degenerate or singular coefficients. The Cauchy problem consists of solving (1.1) subject to the initial condition Jdu(0) = h. We are concerned with the case where the solution is given by an analytic semigroup; it is this sense in which the Canchy problem is parabolic. Sufficient conditions for this to be the case are given in Theorem 1; this is a refinement of previously known results [15] to the linear problem and it extends the related work [13] to the (possibly) degenerate situation under consideration. Specifically, we do not assume ~ is invertible, but only that it is symmetric and non-negative. Our primary motivation for considering the Cauchy problem for (1.1) is to show that certain classes of mixed initial-boundary value problems for partial differential equations are well-posed. Theorem 2 shows that if the operators ~/ and 5fl have additional structure which is typical of those operators arising from (possibly degenerate) parabolic problems then the evolution equation (1.1) is equivalent to a partial differential equation d d~ (Mu(t)) -}- Lu(t) -= F(t) (1.2) (obtained by restricting (1.1) to test functions) and a complementary boundary condition d d-~ (8,~u(t)) + 8zu(t) = g(t) (1.3) in an appropriate space of boundary values.
    [Show full text]
  • Math 246B - Partial Differential Equations
    Math 246B - Partial Differential Equations Viktor Grigoryan version 0.1 - 03/16/2011 Contents Chapter 1: Sobolev spaces 2 1.1 Hs spaces via the Fourier transform . 2 1.2 Weak derivatives . 3 1.3 Sobolev spaces W k;p .................................... 3 1.4 Smooth approximations of Sobolev functions . 4 1.5 Extensions and traces of Sobolev functions . 5 1.6 Sobolev embeddings and compactness results . 6 1.7 Difference quotients . 7 Chapter 2: Solvability of elliptic PDEs 9 2.1 Weak formulation . 9 2.2 Existence of weak solutions of the Dirichlet problem . 10 2.3 General linear elliptic PDEs . 12 2.4 Lax-Milgram theorem, solvability of general elliptic PDEs . 14 2.5 Fredholm operators on Hilbert spaces . 16 2.6 The Fredholm alternative for elliptic equations . 18 2.7 The spectrum of a self-adjoint elliptic operator . 19 Chapter 3: Elliptic regularity theory 21 3.1 Interior regularity . 21 3.2 Boundary regularity . 25 Chapter 4: Variational methods 26 4.1 The Derivative of a functional . 26 4.2 Solvability for the Dirichlet Laplacian . 27 4.3 Constrained optimization and application to eigenvalues . 29 1 1. Sobolev spaces In this chapter we define the Sobolev spaces Hs and W k;p and give their main properties that will be used in subsequent chapters without proof. The proofs of these properties can be found in Evans's\PDE". 1.1 Hs spaces via the Fourier transform Below all the derivatives are understood to be in the distributional sense. Definition 1.1. Let k be a non-negative integer. The Sobolev space HkpRnq is defined as k n 2 n α 2 H pR q tf P L pR q : B f P L for all |α| ¤ ku: k n 2 k p 2 n Theorem 1.2.
    [Show full text]
  • Homogenization and Shape Differentiation of Quasilinear Elliptic Equations Homogeneización Y Diferenciación De Formas De Ecuaciones Elípticas Cuasilineales
    Homogenization and Shape Differentiation of Quasilinear Elliptic Equations Homogeneización y Diferenciación de Formas de Ecuaciones Elípticas Cuasilineales David Gómez Castro Director: Prof. Jesús Ildefonso Díaz Díaz Dpto. de Matemática Aplicada & Instituto de Matemática Interdisciplinar Facultad de Matemáticas Universidad Complutense de Madrid arXiv:1712.10074v1 [math.AP] 28 Dec 2017 Esta tesis se presenta dentro del Programa de Doctorado en Ingeniería Matemática, Estadística e Investigación Operativa Diciembre 2017 Quisiera dedicar esta tesis a mi abuelo, Ángel Castro, con el que tantos momentos compartí de pequeño. Acknowledgements / Agradecimientos The next paragraphs will oscillate between Spanish and English, for the convenience of the different people involved. Antes de empezar, me gustaría expresar mi agradecimiento a todas aquellas personas sin las cuales la realización de esta tesis no hubiese sido posible o, de haberlo sido, hubiese resultado en una experiencia traumática e indeseable. En particular, me gustaría señalar a algunas en concreto. Lo primero agradecer a mis padres (y a mi familia en general) su infinito y constante apoyo, sin el cual, sin duda, no estaría donde estoy. A Ildefonso Díaz, mi director, quién hace años me llamó a su despacho (entonces en el I.M.I.) y me dijo que, en lugar de la idea sobre aerodinámica que se me había ocurrido para el T.F.G., quizás unas cuestiones en las que él trabajaba sobre “Ingeniería Química” podían tener más proyección para nuestro trabajo conjunto. Nadie puede discutirle hoy cuánta razón tenía. Ya desde aquel primer momento he seguido su consejo, y creo que me ha ido bien. Además de su guía temática (y bibliográfica), le quedaré por siempre agradecido por haberme presentado a tantas (y tan célebres) personas, que han enriquecido mis conocimientos, algunas de las cuales se han convertido en colaboradores, y a las que me refiero a más adelante.
    [Show full text]
  • A Some Tools from Mathematical Analysis
    A Some Tools from Mathematical Analysis In this appendix we give a short survey of results, essentially without proofs, from mathematical, particularly functional, analysis which are needed in our treatment of partial differential equations. We begin in Sect. A.1 with a simple account of abstract linear spaces with emphasis on Hilbert space, including the Riesz representation theorem and its generalization to bilinear forms of Lax and Milgram. We continue in Sect. A.2 with function spaces, where after k a discussion of the spaces C , integrability, and the Lp-spaces, we turn to L2-based Sobolev spaces, with the trace theorem and Poincar´e’s inequality. The final Sect. A.3 is concerned with the Fourier transform. A.1 Abstract Linear Spaces Let V be a linear space (or vector space) with real scalars, i.e., a set such that if u, v ∈ V and λ, µ ∈ R, then λu + µv ∈ V .Alinear functional (or linear form) L on V is a function L : V → R such that L(λu + µv)=λL(u)+µL(v), ∀u, v ∈ V, λ,µ ∈ R. A bilinear form a(·, ·)onV is a function a : V × V → R, which is linear in each argument separately, i.e., such that, for all u, v, w ∈ V and λ, µ ∈ R, a(λu + µv, w)=λa(u, w)+µa(v,w), a(w, λu + µv)=λa(w, u)+µa(w, v). The bilinear form a(·, ·) is said to be symmetric if a(w, v)=a(v,w), ∀v,w ∈ V, and positive definite if a(v,v) > 0, ∀v ∈ V, v =0 .
    [Show full text]
  • Recent Progress in Elliptic Equations and Systems of Arbitrary Order With
    Recent progress in elliptic equations and systems of arbitrary order with rough coefficients in Lipschitz domains Vladimir Maz’ya Department of Mathematical Sciences University of Liverpool, M&O Building Liverpool, L69 3BX, UK Department of Mathematics Link¨oping University SE-58183 Link¨oping, Sweden Tatyana Shaposhnikova Department of Mathematics Link¨oping University SE-58183 Link¨oping, Sweden Abstract. This is a survey of results mostly relating elliptic equations and systems of arbitrary even order with rough coefficients in Lipschitz graph domains. Asymp- totic properties of solutions at a point of a Lipschitz boundary are also discussed. 2010 MSC. Primary: 35G15, 35J55; Secondary: 35J67, 35E05 Keywords: higher order elliptic equations, higher order elliptic systems, Besov spaces, mean oscillations, BMO, VMO, Lipschitz domains, Green’s function, asymp- totic behaviour of solutions Introduction The fundamental role in the theory of linear elliptic equations and systems is played by results on regularity of solutions up to the boundary of a domain. The following classical example serves as an illustration. arXiv:1010.0480v1 [math.AP] 4 Oct 2010 Consider the Dirichlet problem ∆u = 0 in Ω, ( Tr u = f on ∂Ω, where Ω is a bounded domain with smooth boundary and Tr u stands for the boundary value (trace) of u. Let u L1(Ω), 1 <p< , that is ∈ p ∞ u pdx < |∇ | ∞ ZΩ 1−1/p and let f belong to the Besov space Bp (∂Ω) with the seminorm f(x) f(y) p 1/p | −n+p−2| dσxdσy . ∂Ω ∂Ω x y Z Z | − | 1 It is well-known that 1 1−1/p Tr Lp(Ω) = Bp (∂Ω).
    [Show full text]
  • Acoustics, Mechanics, Mathematical Flnalqsis
    Acoustics, Mechanics, and the Related Topics of Mathematical flnalqsis This page intentionally left blank Proceedings of the International Conference to Celebrate Robert P. Gilbert's 70th Birthday Acoustics, Mechanics, and the Related Topics of Mathema tical Analqsis CAES du CNRS, Frejus, France 18 - 22 June 2002 Editor Armand Wirgin Laboratoire de Mecanique et dAcoustique Marseille, France orld Scientific Jersey London Singapore Hong Kong Published by World Scientific Publishing Co. Re. Ltd. 5 Toh Tuck Link, Singapore 596224 USA ofice: Suite 202,1060 Main Street, River Edge, NJ 07661 UK ofice: 57 Shelton Street, Covent Garden, London WCZH 9HE British Library Cataloguing-in-PublicationData A catalogue record for this book is available from the British Library. ACOUSTICS, MECHANCIS, AND THE. RELATED TOPICS OF MATHEMATICAL ANALYSIS Proceedings of the International Conference to Celebrate Robert P. Gilbert’s 70th Birthday Copyright 0 2002 by World Scientific Publishing Co. Pte. Ltd. All rights reserved. This book, or parts thereof, may not be reproduced in any form or by any means, eZectronic or mechanical, including photocopying, recording or any information storage and retrieval system now known or to be invented, without written permission from the Publisher. For photocopying of material in this volume, please pay a copying fee through the Copyright Clearance Center, Inc., 222 Rosewood Drive, Danvers, MA 01923, USA. In this case permission to photocopy is not required from the publisher. ISBN 981-238-264-X This book is printed on acid-free paper. Printed in Singapore by Mainland Press PREFACE The international conference Acoustics, Mechanics, and the Related Topics of Mathematical Analysis (AMRTMA) was held on 18-22 June, 2002 at the Villa Clythia, CAES du CNRS, in Frbjus, France.
    [Show full text]
  • Integral Methods in Science and Engineering, Volume 1 Theoretical Techniques
    Christian Constanda Matteo Dalla Riva Pier Domenico Lamberti Paolo Musolino Editors Integral Methods in Science and Engineering, Volume 1 Theoretical Techniques Christian Constanda • Matteo Dalla Riva Pier Domenico Lamberti • Paolo Musolino Editors Integral Methods in Science and Engineering, Volume 1 Theoretical Techniques Editors Christian Constanda Matteo Dalla Riva Department of Mathematics Department of Mathematics The University of Tulsa The University of Tulsa Tulsa, OK, USA Tulsa, OK, USA Pier Domenico Lamberti Paolo Musolino Department of Mathematics Systems Analysis, Prognosis and Control University of Padova Fraunhofer Institute for Industrial Math Padova, Italy Kaiserslautern, Germany ISBN 978-3-319-59383-8 ISBN 978-3-319-59384-5 (eBook) DOI 10.1007/978-3-319-59384-5 Library of Congress Control Number: 2017948080 © Springer International Publishing AG 2017 This work is subject to copyright. All rights are reserved by the Publisher, whether the whole or part of the material is concerned, specifically the rights of translation, reprinting, reuse of illustrations, recitation, broadcasting, reproduction on microfilms or in any other physical way, and transmission or information storage and retrieval, electronic adaptation, computer software, or by similar or dissimilar methodology now known or hereafter developed. The use of general descriptive names, registered names, trademarks, service marks, etc. in this publication does not imply, even in the absence of a specific statement, that such names are exempt from the relevant protective laws and regulations and therefore free for general use. The publisher, the authors and the editors are safe to assume that the advice and information in this book are believed to be true and accurate at the date of publication.
    [Show full text]
  • Computation and Visualization of Geometric Partial Differential Equations
    UNIVERSITY OF CALIFORNIA, SAN DIEGO Computation and Visualization of Geometric Partial Differential Equations A dissertation submitted in partial satisfaction of the requirements for the degree Doctor of Philosophy in Mathematics by Christopher L. Tiee Committee in charge: Professor Michael Holst, Chair Professor Bennett Chow Professor Ken Intriligator Professor Xanthippi Markenscoff Professor Jeff Rabin 2015 Copyright Christopher L. Tiee, 2015 All rights reserved. The dissertation of Christopher L. Tiee is approved, and it is acceptable in quality and form for publication on microfilm and electronically: Chair University of California, San Diego 2015 iii DEDICATION To my grandfathers, Henry Hung-yeh Tiee, Ph. D. and Jack Fulbeck, Ph. D. for their inspiration. iv EPIGRAPH L’étude approfondie de la nature est la source la plus féconde des découvertes mathématiques. [Profound study of nature is the most fertile source of mathematical discoveries.] —Joseph Fourier, Theorie Analytique de la Chaleur v TABLE OF CONTENTS Signature Page........................................ iii Dedication.......................................... iv Epigraph...........................................v Table of Contents...................................... vi List of Figures........................................ ix List of Tables......................................... xi List of Supplementary Files................................ xii Acknowledgements..................................... xiii Vita.............................................. xv Abstract
    [Show full text]
  • Numerical Methods for Partial Differential
    Numerical Methods for Partial Differential Equations Joachim Sch¨oberl October 14, 2007 2 Contents 1 Introduction 5 1.1 Classification of PDEs . 5 1.2 Weak formulation of the Poisson Equation . 6 1.3 The Finite Element Method . 8 2 The abstract theory 11 2.1 Basic properties . 11 2.2 Projection onto subspaces . 14 2.3 Riesz Representation Theorem . 15 2.4 Symmetric variational problems . 16 2.5 Coercive variational problems . 17 2.5.1 Approximation of coercive variational problems . 20 2.6 Inf-sup stable variational problems . 21 2.6.1 Approximation of inf-sup stable variational problems . 23 3 Sobolev Spaces 25 3.1 Generalized derivatives . 25 3.2 Sobolev spaces . 27 3.3 Trace theorems and their applications . 28 3.3.1 The trace space H1/2 ......................... 34 3.4 Equivalent norms on H1 and on sub-spaces . 37 4 The weak formulation of the Poisson equation 41 4.1 Shift theorems . 42 5 Finite Element Method 45 5.1 Finite element system assembling . 47 5.2 Finite element error analysis . 49 5.3 A posteriori error estimates . 55 5.4 Non-conforming Finite Element Methods . 63 6 Linear Equation Solvers 69 6.1 Direct linear equation solvers . 70 3 4 CONTENTS 6.2 Iterative equation solvers . 73 6.3 Preconditioning . 82 7 Mixed Methods 95 7.1 Weak formulation of Dirichlet boundary conditions . 95 7.2 A Mixed method for the flux . 96 7.3 Abstract theory . 97 7.4 Analysis of the model problems . 100 7.5 Approximation of mixed systems . 107 8 Applications 111 8.1 The Navier Stokes equation .
    [Show full text]
  • Sobolev Spaces on Subsets of Rn (Via E.G
    Sobolev spaces on non-Lipschitz subsets of Rn with application to boundary integral equations on fractal screens S. N. Chandler-Wilde, D. P. Hewett and A. Moiola Abstract. We study properties of the classical fractional Sobolev spaces n on non-Lipschitz subsets of R . We investigate the extent to which the properties of these spaces, and the relations between them, that hold in the well-studied case of a Lipschitz open set, generalise to non-Lipschitz cases. Our motivation is to develop the functional analytic framework in which to formulate and analyse integral equations on non-Lipschitz sets. In particular we consider an application to boundary integral equations for wave scattering by planar screens that are non-Lipschitz, including cases where the screen is fractal or has fractal boundary. 1. Introduction In this paper we present a self-contained study of Hilbert{Sobolev spaces de- n arXiv:1607.01994v3 [math.FA] 9 Jan 2017 fined on arbitrary open and closed sets of R , aimed at applied and numerical analysts interested in linear elliptic problems on rough domains, in particu- lar in boundary integral equation (BIE) reformulations. Our focus is on the ◦ s s s s s Sobolev spaces H (Ω), H0 (Ω), He (Ω), H (Ω), and HF , all described below, where Ω (respectively F ) is an arbitrary open (respectively closed) subset of Rn. Our goal is to investigate properties of these spaces (in particular, to provide natural unitary realisations for their dual spaces), and to clarify the nature of the relationships between them. Our motivation for writing this paper is recent and current work by two of the authors [8, 10{12] on problems of acoustic scattering by planar screens with rough (e.g.
    [Show full text]
  • Regularizing Effect of Homogeneous Evolution Equations with Perturbation
    REGULARIZING EFFECT OF HOMOGENEOUS EVOLUTION EQUATIONS WITH PERTURBATION DANIEL HAUER CONTENTS 1. Introduction and main results 1 2. Preliminaries 6 3. Homogeneous accretive operators 16 4. Homogeneous completely accretive operators 23 4.1. General framework 24 4.2. Regularizing effect of the associated semigroup 28 5. Application 31 5.1. An elliptic-parabolic boundary-value problem 31 5.2. Framework 31 5.3. Construction of the Dirichlet-to-Neumann operator 33 du 5.4. Global regularity estimates on dt 36 References 37 ABSTRACT. Since the pioneering works by Aronson & Bénilan [C. R. Acad. Sci. Paris Sér., 1979], and Bénilan & Crandall [Johns Hopkins Univ. Press, 1981] it is well-known that first-order evolution problems governed by a non- linear but homogeneous operator admit the smoothing effect that every corre- sponding mild solution is Lipschitz continuous at every positive time. More- over, if the underlying Banach space has the Radon-Nikodým property, then these mild solution is a.e. differentiable, and the time-derivative satisfies glo- bal and point-wise bounds. In this paper, we show that these results remain true if the homogeneous operator is perturbed by a Lipschitz continuous mapping. More precisely, we establish global L1 Aronson-Bénilan type estimates and point-wise Aronson- Bénilan type estimates. We apply our theory to derive global Lq-L∞-estimates on the time-derivative of the perturbed diffusion problem governed by the Dirichlet-to-Neumann operator associated with the p-Laplace-Beltrami oper- arXiv:2004.00483v4 [math.AP] 14 Apr 2021 ator and lower-order terms on a compact Riemannian manifold with a Lips- chitz boundary.
    [Show full text]