DOCSLIB.ORG

Computation and Visualization of Geometric Partial Differential Equations

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read 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 of the Dissertation................................ xvi Introduction.........................................1 0.1 The Main Problem: its Motivation and Antecedents......2 0.2 Part-by-Part Summary.......................7 I Background 10 Chapter 1 Boundary Value Problems........................ 11 1.1 Differential Forms.......................... 13 1.2 Integration of Differential Forms and Hodge Duality..... 19 1.3 Sobolev Spaces of Differential Forms............... 30 1.4 The Extended Trace Theorem................... 36 1.5 Boundary Value Problems with the Hodge Laplacian..... 42 1.6 The Hilbert Space Setting for Elliptic Problems......... 56 1.6.1 Recasting in terms of Sobolev Spaces.......... 57 1.6.2 The General Elliptic Problem............... 63 1.7 The Theory of Weak Solutions................... 69 1.7.1 The Lax-Milgram Theorem................ 71 1.7.2 Basic Existence Theorems................. 72 1.8 Hilbert Complexes......................... 76 1.9 Evolutionary Partial Differential Equations........... 87 1.9.1 Motivation: The Heat Equation............. 88 1.9.2 Bochner Spaces....................... 90 vi Chapter 2 Numerical Methods............................ 96 2.1 The Finite Element Method.................... 97 2.1.1 The Rayleigh-Ritz Method................ 98 2.1.2 The Galërkin Method................... 102 2.2 Details of the Finite Element Method.............. 103 2.2.1 The Basis........................... 104 2.2.2 Shape Functions...................... 106 2.2.3 Computation of the Stiffness Matrix........... 109 2.3 Adding Time Dependence..................... 110 2.4 Numerical Methods for Evolutionary Equations........ 112 2.4.1 Euler Methods....................... 113 2.4.2 Other Methods....................... 117 2.5 Error Estimates for the Finite Element Method......... 120 2.6 Discretization of Differential Forms............... 126 2.6.1 Approximation in Hilbert Complexes.......... 128 2.6.2 Approximation with Variational Crimes........ 130 2.6.3 Polynomial Spaces and Error Estimates for Differen- tial Forms.......................... 133 Chapter 3 Some Methods for Nonlinear Equations................ 141 3.1 Overview............................... 142 3.2 Linearizing the Equation...................... 154 3.3 Adding Time Dependence..................... 156 3.4 Newton’s Method.......................... 157 3.4.1 Kantorovitch’s Theorem.................. 160 3.4.2 Globalizing Newton’s Method.............. 161 II Applications to Evolution Problems 164 Chapter 4 Approximation of Parabolic Equations in Hilbert Complexes... 165 4.0 Abstract................................ 165 4.1 Introduction............................. 166 4.2 The Finite Element Exterior Calculus.............. 172 4.2.1 Hilbert Complexes..................... 172 4.2.2 Approximation of Hilbert Complexes.......... 177 4.2.3 Extension of Elliptic Error Estimates for a Nonzero Harmonic Part....................... 182 4.3 Abstract Evolution Problems................... 194 4.3.1 Overview of Bochner Spaces and Abstract Evolution Problems.......................... 195 4.3.2 Recasting the Problem as an Abstract Evolution Equa- tion.............................. 198 4.4 A Priori Error Estimates for the Abstract Parabolic Problem. 201 vii 4.5 Parabolic Equations on Compact Riemannian Manifolds... 213 4.6 Numerical Experiments and Implementation Notes...... 228 4.7 Conclusion and Future Directions................ 230 4.8 Acknowledgements......................... 232 Chapter 5 Finite Element Methods for Ricci Flow on Surfaces......... 233 5.1 Introduction............................. 233 5.2 Notation and Conventions..................... 234 5.3 The Ricci Flow on Surfaces..................... 238 5.4 Weak Form of the Equation.................... 242 5.5 Numerical Method......................... 245 5.6 A Numerical Experiment...................... 249 5.7 Conclusion and Future Work................... 252 5.8 Acknowledgements......................... 256 III Appendices 257 Appendix A Canonical Geometries.......................... 258 A.1 Introduction to Spectral Geometry................ 258 A.2 Solving Poisson’s Equation..................... 261 A.3 Finding Dirichlet Green’s Functions............... 263 A.4 The Dirichlet Problem....................... 265 A.5 The Neumann Problem...................... 268 Appendix B Examples of Green’s Functions and Robin Masses.......... 272 B.1 In One Dimension.......................... 272 B.2 Two-Dimensional Examples.................... 280 B.3 Two-Dimensional Example: The Hyperbolic Disk....... 292 B.4 Derivations for Neumann Boundary Conditions........ 298 B.5 The Finite Cylinder......................... 303 B.6 Domains with Holes in the Plane and the Bergman Metric.. 308 B.7 Conclusion and Future Work................... 321 Bibliography......................................... 322 Index............................................. 330 viii LIST OF FIGURES Figure 1.1: Vectors forming a parallelepiped..................... 15 Figure 1.2: A nonorientable manifold: the Möbius strip and transition charts; the left and right edge are identified in opposite directions as indi- cated by the black arrows. The interior of the charts are indicated with the respective colored arrows and dashed curve boundaries.. 22 Figure 1.3: Demonstration of the cone condition and its violation: (1.3a): The cone condition. Note that the nontrivial cone fits in the corners (and of course, everywhere else) nicely, although it occasionally requires a rigid motion. (1.3b):...................... 31 Figure 1.4: A 1-form ! (thin black level sets) whose hodge dual ?! (gray field lines) has vanishing trace on the boundary @U. This says the field lines of ?! are tangent to @U........................ 51 Figure 1.5: A form and pseudoform in R2 dual to each other, with the two kinds of boundary conditions in the annulus A {a r b}. (1.5a): dθ, a Æ Ç Ç harmonic form whose Hodge dual has vanishing trace on @A. (“dθ” actually is a form determined by overlaps, θ ( ¼,¼) and...... 54 2 ¡ Figure 1.6: Example of harmonic form on closed manifold (here, a torus).... 55 Figure 1.7: Two generators for the harmonic forms for H˚ 1(A) and H1(A), where A is the annulus {a r b} R2, reflecting the different kinds of Ç Ç ⊆ boundary conditions. Note how different they are, but at the same time, how they are dual in some sense, one having level sets that. 81 1 Figure 2.1: Example tent function constructed for the node 2 ; where the nodes in the mesh are are k , k 0,...,6..................... 104 4 Æ Figure 2.2: The heat equation on a piecewise linear approximation of a sphere (3545 triangles). The solution is graphed in the normal direction of the sphere. The spatial discretization uses a surface finite element method detailed in Chapter4 (based on [28]), and implemented. 117 Figure 2.3: The wave equation on a piecewise linear approximation of a sphere (3545 triangles). The solution is graphed in the normal direction of the sphere. The spatial discretization uses a surface finite element method detailed in Chapter4 (based on [28]), and implemented. 121 Figure 3.1: Graphical illustration of for Newton’s Method on a function f (the graph y f (x) is in blue). At each x on the x-axis, draw a vertical Æ i line (dashed red in the above) to the point (xi , f (xi )). From that point, draw a tangent line (in red). Then xi 1 is the intersection.. 158 Å ix Figure 4.1: A curve M with a triangulation (blue polygonal curve Mh) within a tubular neighborhood U of M. Some normal vectors º are drawn, in red; the distance function ± is measured along this normal. The intersection x of the normal with Mh defines a mapping a from x . 216 Figure 4.2: Approximation of a quarter unit circle (black) with a segment (blue) and quadratic Lagrange interpolation for the normal projection (red). Even though the underlying triangulation is the same (and thus also the mesh size), notice how much better the quadratic.. 217 Figure 4.3: Hodge heat equation for k 2 in a square modeled as a 100
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]
  • 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.
    [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]
  • 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]