Efficient Preconditioned Solution Methods For Elliptic Partial Differential Equations

This e-book presents several research areas of elliptical problems solved by differential equations. The mathematical models explained in this e-book have been contributed by experts in the field and can be applied to a wide range of real life examples. M

During the last decade essential progress has been achieved in the analysis and implementation of multilevel/rnultigrid and domain decomposition methods to explore a variety of real world applications. An important trend in mod ern numerical simulations is the quick improvement of computer technology that leads to the well known paradigm (see, e. g. , [78,179]): high-performance computers make it indispensable to use numerical methods of almost linear complexity in the problem size N, to maintain an adequate scaling between the computing time and improved computer facilities as N increases. In the h-version of the finite element method (FEM), the multigrid iteration real izes an O(N) solver for elliptic differential equations in a domain n c IRd d with N = O(h- ) , where h is the mesh parameter. In the boundary ele ment method (BEM) , the traditional panel clustering, fast multi-pole and wavelet based methods as well as the modern hierarchical matrix techniques are known to provide the data-sparse approximations to the arising fully populated stiffness matrices with almost linear cost O(Nr log?Nr), where 1 d Nr = O(h - ) is the number of degrees of freedom associated with the boundary. The aim of this book is to introduce a wider audience to the use of a new class of efficient numerical methods of almost linear complexity for solving elliptic partial differential equations (PDEs) based on their reduction to the interface.

Preconditioning and the Conjugate Gradient Method in the Context of Solving PDEs is about the interplay between modeling, analysis, discretization, matrix computation, and model reduction. The authors link PDE analysis, functional analysis, and calculus of variations with matrix iterative computation using Krylov subspace methods and address the challenges that arise during formulation of the mathematical model through to efficient numerical solution of the algebraic problem. The book?s central concept, preconditioning of the conjugate gradient method, is traditionally developed algebraically using the preconditioned finite-dimensional algebraic system. In this text, however, preconditioning is connected to the PDE analysis, and the infinite-dimensional formulation of the conjugate gradient method and its discretization and preconditioning are linked together. This text challenges commonly held views, addresses widespread misunderstandings, and formulates thought-provoking open questions for further research.

/homepage/sac/cam/na2000/index.html7-Volume Set now available at special set price ! Over the second half of the 20th century the subject area loosely referred to as numerical analysis of partial differential equations (PDEs) has undergone unprecedented development. At its practical end, the vigorous growth and steady diversification of the field were stimulated by the demand for accurate and reliable tools for computational modelling in physical sciences and engineering, and by the rapid development of computer hardware and architecture. At the more theoretical end, the analytical insight into the underlying stability and accuracy properties of computational algorithms for PDEs was deepened by building upon recent progress in mathematical analysis and in the theory of PDEs. To embark on a comprehensive review of the field of numerical analysis of partial differential equations within a single volume of this journal would have been an impossible task. Indeed, the 16 contributions included here, by some of the foremost world authorities in the subject, represent only a small sample of the major developments. We hope that these articles will, nevertheless, provide the reader with a stimulating glimpse into this diverse, exciting and important field. The opening paper by Thomée reviews the history of numerical analysis of PDEs, starting with the 1928 paper by Courant, Friedrichs and Lewy on the solution of problems of mathematical physics by means of finite differences. This excellent survey takes the reader through the development of finite differences for elliptic problems from the 1930s, and the intense study of finite differences for general initial value problems during the 1950s and 1960s. The formulation of the concept of stability is explored in the Lax equivalence theorem and the Kreiss matrix lemmas. Reference is made to the introduction of the finite element method by structural engineers, and a description is given of the subsequent development and mathematical analysis of the finite element method with piecewise polynomial approximating functions. The penultimate section of Thomée's survey deals with `other classes of approximation methods', and this covers methods such as collocation methods, spectral methods, finite volume methods and boundary integral methods. The final section is devoted to numerical linear algebra for elliptic problems. The next three papers, by Bialecki and Fairweather, Hesthaven and Gottlieb and Dahmen, describe, respectively, spline collocation methods, spectral methods and wavelet methods. The work by Bialecki and Fairweather is a comprehensive overview of orthogonal spline collocation from its first appearance to the latest mathematical developments and applications. The emphasis throughout is on problems in two space dimensions. The paper by Hesthaven and Gottlieb presents a review of Fourier and Chebyshev pseudospectral methods for the solution of hyperbolic PDEs. Particular emphasis is placed on the treatment of boundaries, stability of time discretisations, treatment of non-smooth solutions and multidomain techniques. The paper gives a clear view of the advances that have been made over the last decade in solving hyperbolic problems by means of spectral methods, but it shows that many critical issues remain open. The paper by Dahmen reviews the recent rapid growth in the use of wavelet methods for PDEs. The author focuses on the use of adaptivity, where significant successes have recently been achieved. He describes the potential weaknesses of wavelet methods as well as the perceived strengths, thus giving a balanced view that should encourage the study of wavelet methods.

One of the first things a student of partial differential equations learns is that it is impossible to solve elliptic equations by spatial marching. This new book describes how to do exactly that, providing a powerful tool for solving problems in fluid dynamics, heat transfer, electrostatics, and other fields characterized by discretized partial differential equations. Elliptic Marching Methods and Domain Decomposition demonstrates how to handle numerical instabilities (i.e., limitations on the size of the problem) that appear when one tries to solve these discretized equations with marching methods. The book also shows how marching methods can be superior to multigrid and pre-conditioned conjugate gradient (PCG) methods, particularly when used in the context of multiprocessor parallel computers. Techniques for using domain decomposition together with marching methods are detailed, clearly illustrating the benefits of these techniques for applications in engineering, applied mathematics, and the physical sciences.

Mathematics of Computing -- Parallelism.

Domain decomposition (DD) methods provide powerful tools for constructing parallel numerical solution algorithms for large scale systems of algebraic equations arising from the discretization of partial differential equations. These methods are well-established and belong to a fast developing area. In this volume, the reader will find a brief historical overview, the basic results of the general theory of domain and space decomposition methods as well as the description and analysis of practical DD algorithms for parallel computing. It is typical to find in this volume that most of the presented DD solvers belong to the family of fast algorithms, where each component is efficient with respect to the arithmetical work. Readers will discover new analysis results for both the well-known basic DD solvers and some DD methods recently devised by the authors, e.g., for elliptic problems with varying chaotically piecewise constant orthotropism without restrictions on the finite aspect ratios.The hp finite element discretizations, in particular, by spectral elements of elliptic equations are given significant attention in current research and applications. This volume is the first to feature all components of Dirichlet-Dirichlet-type DD solvers for hp discretizations devised as numerical procedures which result in DD solvers that are almost optimal with respect to the computational work. The most important DD solvers are presented in the matrix/vector form algorithms that are convenient for practical use.

Domain decomposition methods are divide and conquer computational methods for the parallel solution of partial differential equations of elliptic or parabolic type. The methodology includes iterative algorithms, and techniques for non-matching grid discretizations and heterogeneous approximations. This book serves as a matrix oriented introduction to domain decomposition methodology. A wide range of topics are discussed include hybrid formulations, Schwarz, and many more.