Notes on Numerical Fluid Mechanics

Size: px
Start display at page:

Download "Notes on Numerical Fluid Mechanics"

Transcription

1 Notes on Numerical Fluid Mechanics Volume 31 Parallel Algorithm,s for Partial Differential Equations Edited by Wolfgang Hackbusch al V1eweg

2 Parallel Algorithms for Partial Differential Equations Proceedings of the Sixth GAMM-Seminar, Kiel, January 19-21,1990 Edited by Wolfgang Hackbusch vieweg UNIVERSITATSBIBLSOTHEK HANNOVER TECHN1SCHE INFORMATIONSBIBLIOTHEK

3 CIP-Titelaufnahme der Deutschen Bibliothek Parallel algorithms for partial differential equations: Kiel, January 19-21, 1990 / ed. by Wolfgang Hackbusch. - Braunschweig; Wiesbaden: Vieweg, 1991 (Proceedings of the... GAMM seminar; 6) (Notes on numerical fluid mechanics; Vol. 31) ISBN NE: Hackbusch, Wolfgang \Hrsg.} ; Gesellschaft für Angewandte Mathematik und Mechanik: Proceedings of the... ; 2. GT Manuscripts should have well over 100 pages. As they, ill be rep.rodu ed pho ome b.anicall they should be typed with utmost care on special stationary which will be supplied on reque In print, the size will be reduced linearly to approximat ely 75 per cen. F d~ should be lettered accordingly so as to produce letters not smaller than mm prin The is valid for handwritten formulae. Manuscripts (in English} or p roposals hould be n general editor Prof. Dr. E.H. Hirschel, Herzog-Heinrich-Weg 6, D-8011 Zorned Vieweg is a subsidiary company of the Bertelsman n Publi hing Group ln emationa All rights reserved Friedr. Vieweg & Sohn Verlagsgesellschaft mb H, Braun chwei 1991 No part of this publication may be reproduced, ored in a cetrie" system or transmitted, mechanical, photocop in or o ~ ;i ou prior p ermission of the copyright h older. Produced by W. Langelüddecke, Braunschweig Printed in the Federal Republic of Germany ISSN ISBN

4 FOREWORD The GAMM Committee for "Efficient Numerical Methods for Partial Differential Equations" organises workshops on subjects concerning the algorithmical treatment of partial differential equations. The topics are discretisation methods like the finite element and the boundary element method f or various types of applications in structural and fluid mechanics. Particular attention is devoted to the advanced solution techniques like e.g. multi-grid methods. The series of such workshops was continued in 1990, January 19-21, with the 6th Kiel-Seminar on the special topic "Parallel Algorithms for PDEs" at the Christian-Albrechts-University of Kiel. The seminar was attended by 90 scientists from 10 countries and 22 lectures were given. In many lectures applications of the domain decomposition are considered. Combinations with multi-grid methods, conjugate gradient methods, adaptive refinement, and boundary collocation methods are mentioned. But also the orginal domain decomposition algorithm is described. Another decomposition is the problem decomposition, which offers a new kind of parallelism. A parallel variant of the conjugate gradient method is described. Concerning multi-grid methods quite new features are introduced: correction schemes replacing the coarse-grid correction and so-called sparse grids. The latter co~struc~ion leads to completely different discrete problems, wh1ch give accurate results although the dimension of the system is strongly reduced. All topics mentioned above are contained in the 19 contributions to this volume. The editor thanks the Land Schleswig-Holstein for its support. He also likes to express his gratitude to all persons involved in the organisation of the seminar. Kiel, July 1990 W. Hackbusch

5 CONTENTS E. ALLGOWER, K. BOHMER, M. ZHEN: On a problem decomposition for semilinear nearly symmetric elliptic problems 1 P.BASTIAN, J. BURMEISTER, G. HORTON: Implementation of a parallel multi-grid method for parabolic partial differential equations 18 P.E. BJORSTAD, R. MOE, M. SKOGEN: Parallel decomposition and iterative refinement algorithms 28 Th. DREYER: Nonlinear ilu-decomposition as smoother in multigrid methods 47 M. DRYJA, O.B. WIDLUND: Multilevel additive methods for elliptic finite element problems 58 R.E. EWING, R.D.LAZAROV, P.S. VASSILEVSKI: Finite difference schemes on grids with local refinement in time and in space for parabolic problems. Optimal order two-grid iterative methods 70 M. GRIEBEL: A parallelizable and vectorizable multi-level algorithm on sparse grids 94 G. HAASE, U. LANGER: On the use of multigrid preconditioners in the domain decomposition method 101 W. HACKBUSCH: A parallel variant of the conjugate gradient method Ill G. HORTON: Parallel solution of the Navier-Stokes equations in complex geometries 120 B. KISS, G. MOLNARKA: A preconditioned domain decomposition algorithm for the solution of the elliptic Neumann problem. 129 J.A. KOLODZIEJ, G. MUSIELAK: Domain decomposition in boundary collocation method 140 A. KRECHEL, H.-J. PLUM: A parallel block preconditioned conjugate gradient method 157 R.W. LELAND, J.S. ROLLETT: A comparison of parallel extrapolation methods 177 V. MEHRMANN: Divide and conquer methods for tridiagonal linear systems 188

6 H. SBOSNY, K. WITSCH: Parallel multigrid using domain decomposition St. VANDEWALLE, R. PIESSENS: A parallel and vectorizable algorithm for solving parabolic partial differential equations G. WITTUM: An ILU-based smoothing correction scheme Ch. ZENGER: Sparse Grids 241 List of participants 252

Modeling and Computation in Environmental Sciences. Edited by Rainer Helmig Willi Jager Wolfgang Kinzelbach Peter Knabner and Gabriel Wittum

Modeling and Computation in Environmental Sciences. Edited by Rainer Helmig Willi Jager Wolfgang Kinzelbach Peter Knabner and Gabriel Wittum Modeling and Computation in Environmental Sciences Edited by Rainer Helmig Willi Jager Wolfgang Kinzelbach Peter Knabner and Gabriel Wittum Notes on Numerical Fluid Mechanics (NNFM) Volume 59 Series Editors:

More information

Wilfred W. J. Hulsbergen. Coniectures in Arithmetic Algebraic Geometry

Wilfred W. J. Hulsbergen. Coniectures in Arithmetic Algebraic Geometry Wilfred W. J. Hulsbergen Coniectures in Arithmetic Algebraic Geometry As~ tvbthematics Edited by Klas Diederich Vol. E1: Vol. E2: VoI.E3: Vol. E4: Vol. E5: VoI.E6: Vol. E7: Vol. E8: Vol. E 9: Vol. E 10:

More information

Jean-Pierre Serre. Lectures on the. Mordell-Weil Theorem

Jean-Pierre Serre. Lectures on the. Mordell-Weil Theorem Jean-Pierre Serre Lectures on the Mordell-Weil Theorem Aspeds cl Mathematics Aspekte der Mathematik Editor: Klas Diederich All volumes of the series are listed on pages 219-220 Jean-Pierre Serre Lectures

More information

Raymond Gerard Hidetoshi Tahara. Singular Nonlinear Partial Differential Equations

Raymond Gerard Hidetoshi Tahara. Singular Nonlinear Partial Differential Equations Raymond Gerard Hidetoshi Tahara Singular Nonlinear Partial Differential Equations Aspects Mathematics Edited by Klas Diederich Vol. E 3: Vol. E 5: Vol. E 6: Vol. E 7: Vol. E 9: Vol. E 10: G. Hector/U.

More information

Preface to the Second Edition. Preface to the First Edition

Preface to the Second Edition. Preface to the First Edition n page v Preface to the Second Edition Preface to the First Edition xiii xvii 1 Background in Linear Algebra 1 1.1 Matrices................................. 1 1.2 Square Matrices and Eigenvalues....................

More information

Notes on Numerical Fluid Mechanics

Notes on Numerical Fluid Mechanics Notes on Numerical Fluid Mechanics Josef Ballmann Richard Eppler Wolfgang Hackbusch Hrsg. Panel Methods in Fluid Mechanics with Emphasis on Aerodynamics Proceedings of the Third GAMM- Seminar Kiel, January

More information

C. Vuik 1 R. Nabben 2 J.M. Tang 1

C. Vuik 1 R. Nabben 2 J.M. Tang 1 Deflation acceleration of block ILU preconditioned Krylov methods C. Vuik 1 R. Nabben 2 J.M. Tang 1 1 Delft University of Technology J.M. Burgerscentrum 2 Technical University Berlin Institut für Mathematik

More information

Nonlinear Hyperbolic Equations - Theory, Computation Methods, and Applications

Nonlinear Hyperbolic Equations - Theory, Computation Methods, and Applications Josef Ballmann Rolf Jeltsch (Eds.) Nonlinear Hyperbolic Equations - Theory, Computation Methods, and Applications Proceedings of the Second International Conference on Nonlinear Hyperbolic Problems, Aachen,

More information

On the choice of abstract projection vectors for second level preconditioners

On the choice of abstract projection vectors for second level preconditioners On the choice of abstract projection vectors for second level preconditioners C. Vuik 1, J.M. Tang 1, and R. Nabben 2 1 Delft University of Technology 2 Technische Universität Berlin Institut für Mathematik

More information

Newton s Method and Efficient, Robust Variants

Newton s Method and Efficient, Robust Variants Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science

More information

Multigrid and Iterative Strategies for Optimal Control Problems

Multigrid and Iterative Strategies for Optimal Control Problems Multigrid and Iterative Strategies for Optimal Control Problems John Pearson 1, Stefan Takacs 1 1 Mathematical Institute, 24 29 St. Giles, Oxford, OX1 3LB e-mail: john.pearson@worc.ox.ac.uk, takacs@maths.ox.ac.uk

More information

D-,+O K.1K. b9 V 378

D-,+O K.1K. b9 V 378 CLSSIICATON ALD -A207 000 REPORT DOCUMENTATION PAGE ors Appove1 la. REPORT SECURITY CLASSIFICATION 1lb. RESTRICTIVE MARKINGS 2a. SECURITY CLASSIFICATION AUT _"ALECTE 3. DISTRIBUTION /AVAILABILITY OF REPOR

More information

A MULTILEVEL SUCCESSIVE ITERATION METHOD FOR NONLINEAR ELLIPTIC PROBLEMS

A MULTILEVEL SUCCESSIVE ITERATION METHOD FOR NONLINEAR ELLIPTIC PROBLEMS MATHEMATICS OF COMPUTATION Volume 73, Number 246, Pages 525 539 S 0025-5718(03)01566-7 Article electronically published on July 14, 2003 A MULTILEVEL SUCCESSIVE ITERATION METHOD FOR NONLINEAR ELLIPTIC

More information

Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations

Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg

More information

Book of Abstracts - Extract th Annual Meeting. March 23-27, 2015 Lecce, Italy. jahrestagung.gamm-ev.de

Book of Abstracts - Extract th Annual Meeting. March 23-27, 2015 Lecce, Italy. jahrestagung.gamm-ev.de GESELLSCHAFT für ANGEWANDTE MATHEMATIK und MECHANIK e.v. INTERNATIONAL ASSOCIATION of APPLIED MATHEMATICS and MECHANICS 86 th Annual Meeting of the International Association of Applied Mathematics and

More information

Numerical Approximation Methods for Elliptic Boundary Value Problems

Numerical Approximation Methods for Elliptic Boundary Value Problems Numerical Approximation Methods for Elliptic Boundary Value Problems Olaf Steinbach Numerical Approximation Methods for Elliptic Boundary Value Problems Finite and Boundary Elements Olaf Steinbach Institute

More information

Numerical Solution I

Numerical Solution I Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions

More information

Scientific Computing on Supercomputers III

Scientific Computing on Supercomputers III Scientific Computing on Supercomputers III Scientific Computing on Supercomputers III Edited by J ozef T. Devreese and Piet E. Van Camp Universiteit Antwerpen Antwerpen, Belgium SPRINGER SCIENCE+BUSINESS

More information

Robust Parallel Newton{Multilevel Methods. Bodo Heise. Institute for Mathematics, Johannes Kepler University Linz, A Linz, Austria.

Robust Parallel Newton{Multilevel Methods. Bodo Heise. Institute for Mathematics, Johannes Kepler University Linz, A Linz, Austria. Robust Parallel Newton{Multilevel Methods Bodo Heise Institute for Mathematics, Johannes Kepler University Linz, A - 4040 Linz, Austria. e-mail: heise@miraculix.numa.uni-linz.ac.at Michael Jung Faculty

More information

2 CAI, KEYES AND MARCINKOWSKI proportional to the relative nonlinearity of the function; i.e., as the relative nonlinearity increases the domain of co

2 CAI, KEYES AND MARCINKOWSKI proportional to the relative nonlinearity of the function; i.e., as the relative nonlinearity increases the domain of co INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2002; 00:1 6 [Version: 2000/07/27 v1.0] Nonlinear Additive Schwarz Preconditioners and Application in Computational Fluid

More information

Curriculum Vitae. Yvan Notay. (version December 2018)

Curriculum Vitae. Yvan Notay. (version December 2018) Curriculum Vitae Yvan Notay (version December 2018) Personal Place / date of birth: Anderlecht, Belgium / October 1, 1964 Citizenship: Belgian Work address Yvan Notay Université Libre de Bruxelles (U.L.B.)

More information

FINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET OPTIMAL CONTROL PROBLEMS

FINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET OPTIMAL CONTROL PROBLEMS Numerical Functional Analysis and Optimization, 28(7 8):957 973, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0163-0563 print/1532-2467 online DOI: 10.1080/01630560701493305 FINITE ELEMENT APPROXIMATION

More information

A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations

A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for

More information

1. Fast Iterative Solvers of SLE

1. Fast Iterative Solvers of SLE 1. Fast Iterative Solvers of crucial drawback of solvers discussed so far: they become slower if we discretize more accurate! now: look for possible remedies relaxation: explicit application of the multigrid

More information

Introduction to Numerical Analysis

Introduction to Numerical Analysis J. Stoer R. Bulirsch Introduction to Numerical Analysis Translated by R. Bartels, W. Gautschi, and C. Witzgall Springer Science+Business Media, LLC J. Stoer R. Bulirsch Institut fiir Angewandte Mathematik

More information

AMS526: Numerical Analysis I (Numerical Linear Algebra)

AMS526: Numerical Analysis I (Numerical Linear Algebra) AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 24: Preconditioning and Multigrid Solver Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 5 Preconditioning Motivation:

More information

Christian Okonek Michael Schneider Heinz SRindler. ector undies on omplex. rojective S aces

Christian Okonek Michael Schneider Heinz SRindler. ector undies on omplex. rojective S aces Christian Okonek Michael Schneider Heinz SRindler ector undies on omplex rojective S aces Progress in Mathe~natics Vol. 1: H. Gross, Quadratic Forms in Infinite-Dimensional Vector Spaces. XXII, 4!9 pages,!979

More information

Editors-in-Chief Anne Boutet de Monvel, Université Paris VII Denis Diderot, France Gerald Kaiser, Center for Signals and Waves, Austin, TX, USA

Editors-in-Chief Anne Boutet de Monvel, Université Paris VII Denis Diderot, France Gerald Kaiser, Center for Signals and Waves, Austin, TX, USA Progress in Mathematical Physics Volume 45 Editors-in-Chief Anne Boutet de Monvel, Université Paris VII Denis Diderot, France Gerald Kaiser, Center for Signals and Waves, Austin, TX, USA Editorial Board

More information

2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9

2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9 Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers

More information

Rudolf Taschner. The Continuum

Rudolf Taschner. The Continuum Rudolf Taschner The Continuum Rudolf Taschner The Continuum A Constructive Approach to Basic Concepts of Real Analysis ai vleweg Bibliographic information published by Die Deutsche Bibliothek Die Deutsche

More information

Partial Differential Equations with Numerical Methods

Partial Differential Equations with Numerical Methods Stig Larsson Vidar Thomée Partial Differential Equations with Numerical Methods May 2, 2003 Springer Berlin Heidelberg New York Barcelona Hong Kong London Milan Paris Tokyo Preface Our purpose in this

More information

Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction

Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de

More information

, /'/ NI'9 ; / Proceedings of the Third International. High Order Methods. Conference on Spectral and. Houston, Texas

, /'/ NI'9 ; / Proceedings of the Third International. High Order Methods. Conference on Spectral and. Houston, Texas ; /, /'/ NI'9 - Edited by Andrew V. Ilin and L. Ridgway Scott Proceedings of the Third International Conference on Spectral and High Order Methods Houston, Texas 5-9 June 1995 Houston Journal of Mathematics

More information

OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU

OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative

More information

The solution of the discretized incompressible Navier-Stokes equations with iterative methods

The solution of the discretized incompressible Navier-Stokes equations with iterative methods The solution of the discretized incompressible Navier-Stokes equations with iterative methods Report 93-54 C. Vuik Technische Universiteit Delft Delft University of Technology Faculteit der Technische

More information

- Part 4 - Multicore and Manycore Technology: Chances and Challenges. Vincent Heuveline

- Part 4 - Multicore and Manycore Technology: Chances and Challenges. Vincent Heuveline - Part 4 - Multicore and Manycore Technology: Chances and Challenges Vincent Heuveline 1 Numerical Simulation of Tropical Cyclones Goal oriented adaptivity for tropical cyclones ~10⁴km ~1500km ~100km 2

More information

The Conjugate Gradient Method

The Conjugate Gradient Method The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large

More information

Torge Geodesy. Unauthenticated Download Date 1/9/18 5:16 AM

Torge Geodesy. Unauthenticated Download Date 1/9/18 5:16 AM Torge Geodesy Wolfgang Torge Geodesy Second Edition W DE G Walter de Gruyter Berlin New York 1991 Author Wolfgang Torge, Univ. Prof. Dr.-Ing. Institut für Erdmessung Universität Hannover Nienburger Strasse

More information

OUTLINE 1. Introduction 1.1 Notation 1.2 Special matrices 2. Gaussian Elimination 2.1 Vector and matrix norms 2.2 Finite precision arithmetic 2.3 Fact

OUTLINE 1. Introduction 1.1 Notation 1.2 Special matrices 2. Gaussian Elimination 2.1 Vector and matrix norms 2.2 Finite precision arithmetic 2.3 Fact Computational Linear Algebra Course: (MATH: 6800, CSCI: 6800) Semester: Fall 1998 Instructors: { Joseph E. Flaherty, aherje@cs.rpi.edu { Franklin T. Luk, luk@cs.rpi.edu { Wesley Turner, turnerw@cs.rpi.edu

More information

Comparison of Fixed Point Methods and Krylov Subspace Methods Solving Convection-Diffusion Equations

Comparison of Fixed Point Methods and Krylov Subspace Methods Solving Convection-Diffusion Equations American Journal of Computational Mathematics, 5, 5, 3-6 Published Online June 5 in SciRes. http://www.scirp.org/journal/ajcm http://dx.doi.org/.436/ajcm.5.5 Comparison of Fixed Point Methods and Krylov

More information

/00 $ $.25 per page

/00 $ $.25 per page Contemporary Mathematics Volume 00, 0000 Domain Decomposition For Linear And Nonlinear Elliptic Problems Via Function Or Space Decomposition UE-CHENG TAI Abstract. In this article, we use a function decomposition

More information

A Space-Time Multigrid Solver Methodology for the Optimal Control of Time-Dependent Fluid Flow

A Space-Time Multigrid Solver Methodology for the Optimal Control of Time-Dependent Fluid Flow A Space-Time Multigrid Solver Methodology for the Optimal Control of Time-Dependent Fluid Flow Michael Köster, Michael Hinze, Stefan Turek Michael Köster Institute for Applied Mathematics TU Dortmund Trier,

More information

Parallelizing large scale time domain electromagnetic inverse problem

Parallelizing large scale time domain electromagnetic inverse problem Parallelizing large scale time domain electromagnetic inverse problems Eldad Haber with: D. Oldenburg & R. Shekhtman + Emory University, Atlanta, GA + The University of British Columbia, Vancouver, BC,

More information

Aspects of Multigrid

Aspects of Multigrid Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)

More information

Rakhim Aitbayev, Ph.D.

Rakhim Aitbayev, Ph.D. Curriculum Vitae 1 Rakhim Aitbayev, Ph.D. Department of Mathematics New Mexico Institute of Mining and Technology Socorro, New Mexico 87801 phone: (575) 835-5463 e-mail: aitbayev@nmt.edu webpage: www.nmt.edu/

More information

On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities

On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)

More information

Multigrid absolute value preconditioning

Multigrid absolute value preconditioning Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical

More information

Introduction to Computational Fluid Dynamics

Introduction to Computational Fluid Dynamics AML2506 Biomechanics and Flow Simulation Day Introduction to Computational Fluid Dynamics Session Speaker Dr. M. D. Deshpande M.S. Ramaiah School of Advanced Studies - Bangalore 1 Session Objectives At

More information

Transactions on Modelling and Simulation vol 18, 1997 WIT Press, ISSN X

Transactions on Modelling and Simulation vol 18, 1997 WIT Press,   ISSN X Application of the panel clustering method to the three-dimensional elastostatic problem K. Hayami* and S. A. Sauter** ^Department of Mathematical Engineering and Information Physics, School ofengineering,

More information

Initial Boundary Value Problems in Mathematical Physics

Initial Boundary Value Problems in Mathematical Physics Initial Boundary Value Problems in Mathematical Physics Initial Boundary Value Problems in Mathematical Physics Rolf leis University of Bonn Federal Republic of Germany Springer Fachmedien Wiesbaden GmbH

More information

Partial Differential Equations

Partial Differential Equations Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This

More information

DELFT UNIVERSITY OF TECHNOLOGY

DELFT UNIVERSITY OF TECHNOLOGY DELFT UNIVERSITY OF TECHNOLOGY REPORT 12-10 Fast linear solver for pressure computation in layered domains P. van Slingerland and C. Vuik ISSN 1389-6520 Reports of the Department of Applied Mathematical

More information

Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses

Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,

More information

PDE Solvers for Fluid Flow

PDE Solvers for Fluid Flow PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic

More information

Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions

Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz

More information

ADDITIVE SCHWARZ FOR SCHUR COMPLEMENT 305 the parallel implementation of both preconditioners on distributed memory platforms, and compare their perfo

ADDITIVE SCHWARZ FOR SCHUR COMPLEMENT 305 the parallel implementation of both preconditioners on distributed memory platforms, and compare their perfo 35 Additive Schwarz for the Schur Complement Method Luiz M. Carvalho and Luc Giraud 1 Introduction Domain decomposition methods for solving elliptic boundary problems have been receiving increasing attention

More information

Optimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36

Optimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36 Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,

More information

HIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

HIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume LIII, Number 1, March 2008 HIGHER-ORDER LINEARLY IMPLICIT ONE-STEP METHODS FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IOAN TELEAGA AND JENS

More information

IMPLEMENTATION OF A PARALLEL AMG SOLVER

IMPLEMENTATION OF A PARALLEL AMG SOLVER IMPLEMENTATION OF A PARALLEL AMG SOLVER Tony Saad May 2005 http://tsaad.utsi.edu - tsaad@utsi.edu PLAN INTRODUCTION 2 min. MULTIGRID METHODS.. 3 min. PARALLEL IMPLEMENTATION PARTITIONING. 1 min. RENUMBERING...

More information

From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes

From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,

More information

Nonlinear Parabolic and Elliptic Equations

Nonlinear Parabolic and Elliptic Equations Nonlinear Parabolic and Elliptic Equations Nonlinear Parabolic and Elliptic Equations c. V. Pao North Carolina State University Raleigh, North Carolina Plenum Press New York and London Library of Congress

More information

European Consortium for Mathematics in Industry E. Eich-Soellner and C. FUhrer Numerical Methods in Multibody Dynamics

European Consortium for Mathematics in Industry E. Eich-Soellner and C. FUhrer Numerical Methods in Multibody Dynamics European Consortium for Mathematics in Industry E. Eich-Soellner and C. FUhrer Numerical Methods in Multibody Dynamics European Consortium for Mathematics in Industry Edited by Leif Arkeryd, Goteborg Heinz

More information

Lecture Notes in Mathematics

Lecture Notes in Mathematics Lecture Notes in Mathematics A collection of informal reports and seminars Edited by A. Dold, Heidelberg and B. Eckmann, Zorich 309 David H. Sattinger University of Minnesota, Minneapolis, MN/USA Topics

More information

Robust Preconditioned Conjugate Gradient for the GPU and Parallel Implementations

Robust Preconditioned Conjugate Gradient for the GPU and Parallel Implementations Robust Preconditioned Conjugate Gradient for the GPU and Parallel Implementations Rohit Gupta, Martin van Gijzen, Kees Vuik GPU Technology Conference 2012, San Jose CA. GPU Technology Conference 2012,

More information

Short title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic

Short title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic Short title: Total FETI Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ-70833 Ostrava, Czech Republic mail: zdenek.dostal@vsb.cz fax +420 596 919 597 phone

More information

Master Thesis Literature Study Presentation

Master Thesis Literature Study Presentation Master Thesis Literature Study Presentation Delft University of Technology The Faculty of Electrical Engineering, Mathematics and Computer Science January 29, 2010 Plaxis Introduction Plaxis Finite Element

More information

Adaptive sparse grids

Adaptive sparse grids ANZIAM J. 44 (E) ppc335 C353, 2003 C335 Adaptive sparse grids M. Hegland (Received 1 June 2001) Abstract Sparse grids, as studied by Zenger and Griebel in the last 10 years have been very successful in

More information

MPI parallel implementation of CBF preconditioning for 3D elasticity problems 1

MPI parallel implementation of CBF preconditioning for 3D elasticity problems 1 Mathematics and Computers in Simulation 50 (1999) 247±254 MPI parallel implementation of CBF preconditioning for 3D elasticity problems 1 Ivan Lirkov *, Svetozar Margenov Central Laboratory for Parallel

More information

Computational Graph Theory

Computational Graph Theory Computing Supplementum 7 G. Tinhofer, E. Mayr, H. Noltemeier, M. M. Syslo (eds.) in cooperation with R. Albrecht Computational Graph Theory Springer-Verlag Wien New York Prof. Dr. G. Tinhofer Institut

More information

A Recursive Trust-Region Method for Non-Convex Constrained Minimization

A Recursive Trust-Region Method for Non-Convex Constrained Minimization A Recursive Trust-Region Method for Non-Convex Constrained Minimization Christian Groß 1 and Rolf Krause 1 Institute for Numerical Simulation, University of Bonn. {gross,krause}@ins.uni-bonn.de 1 Introduction

More information

Optimal Control of Partial Differential Equations I+II

Optimal Control of Partial Differential Equations I+II About the lecture: Optimal Control of Partial Differential Equations I+II Prof. Dr. H. J. Pesch Winter semester 2011/12, summer semester 2012 (seminar), winter semester 2012/13 University of Bayreuth Preliminary

More information

Dieter Kind Hermann Karner. High -Voltage Insulation Technology

Dieter Kind Hermann Karner. High -Voltage Insulation Technology Dieter Kind Hermann Karner High -Voltage Insulation Technology Dieter Kind Hermann Kamer High-Voltage Insulation Thchnology Thxtbook for Electrical Engineers lfanslated from the German by Y Narayana Rao

More information

APPLIED NUMERICAL LINEAR ALGEBRA

APPLIED NUMERICAL LINEAR ALGEBRA APPLIED NUMERICAL LINEAR ALGEBRA James W. Demmel University of California Berkeley, California Society for Industrial and Applied Mathematics Philadelphia Contents Preface 1 Introduction 1 1.1 Basic Notation

More information

Algebraic Multigrid as Solvers and as Preconditioner

Algebraic Multigrid as Solvers and as Preconditioner Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven

More information

A Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers

A Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned

More information

SOLVING ELLIPTIC PDES

SOLVING ELLIPTIC PDES university-logo SOLVING ELLIPTIC PDES School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 POISSON S EQUATION Equation and Boundary Conditions Solving the Model Problem 3 THE LINEAR ALGEBRA PROBLEM

More information

Auxiliary space multigrid method for elliptic problems with highly varying coefficients

Auxiliary space multigrid method for elliptic problems with highly varying coefficients Auxiliary space multigrid method for elliptic problems with highly varying coefficients Johannes Kraus 1 and Maria Lymbery 2 1 Introduction The robust preconditioning of linear systems of algebraic equations

More information

M.A. Botchev. September 5, 2014

M.A. Botchev. September 5, 2014 Rome-Moscow school of Matrix Methods and Applied Linear Algebra 2014 A short introduction to Krylov subspaces for linear systems, matrix functions and inexact Newton methods. Plan and exercises. M.A. Botchev

More information

The All-floating BETI Method: Numerical Results

The All-floating BETI Method: Numerical Results The All-floating BETI Method: Numerical Results Günther Of Institute of Computational Mathematics, Graz University of Technology, Steyrergasse 30, A-8010 Graz, Austria, of@tugraz.at Summary. The all-floating

More information

Numerical Methods for the Solution of Ill-Posed Problems

Numerical Methods for the Solution of Ill-Posed Problems Numerical Methods for the Solution of Ill-Posed Problems Mathematics and Its Applications Managing Editor: M.HAZEWINKEL Centre for Mathematics and Computer Science, Amsterdam, The Netherlands Volume 328

More information

SOLVING MESH EIGENPROBLEMS WITH MULTIGRID EFFICIENCY

SOLVING MESH EIGENPROBLEMS WITH MULTIGRID EFFICIENCY SOLVING MESH EIGENPROBLEMS WITH MULTIGRID EFFICIENCY KLAUS NEYMEYR ABSTRACT. Multigrid techniques can successfully be applied to mesh eigenvalue problems for elliptic differential operators. They allow

More information

Curriculum Vitae. Gender: Male Date of birth: 8th July Place of birth: Novi Sad, municipality of Novi Sad, Serbia Present Citizenship: Serbian

Curriculum Vitae. Gender: Male Date of birth: 8th July Place of birth: Novi Sad, municipality of Novi Sad, Serbia Present Citizenship: Serbian Curriculum Vitae Damir Ma darević PhD in Mechanical Engineering Address Work Trg Dositeja Obradovića 6. 21000 Novi Sad Republic of Serbia Contact Phone (work): +381 21 485 22 43 E-mail: damirm@uns.ac.rs

More information

Goal-oriented error control of the iterative solution of finite element equations

Goal-oriented error control of the iterative solution of finite element equations J. Numer. Math., Vol. 0, No. 0, pp. 1 31 (2009) DOI 10.1515/ JNUM.2009.000 c de Gruyter 2009 Prepared using jnm.sty [Version: 20.02.2007 v1.3] Goal-oriented error control of the iterative solution of finite

More information

The Mathematics of Computerized Tomography

The Mathematics of Computerized Tomography The Mathematics of Computerized Tomography The Mathematics of Computerized Tomography F. Natterer University of Münster Federal Republic of Germany B. G. TEUBNER Stuttgart @) JOHN WILEY & SONS Chichester.

More information

An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems

An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,

More information

NUMERICAL METHODS FOR ENGINEERING APPLICATION

NUMERICAL METHODS FOR ENGINEERING APPLICATION NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto

More information

Equations Introduction

Equations Introduction Chapter 9 90 Introduction Partial Differential Equations The numerical treatment of partial differential equations is, by itself, a vast subject Partial differential equations are at the heart of many,

More information

Seismicity in Mines. Edited by Slawomir J. Gibowicz. Birkhauser Verlag Basel Boston. Berlin

Seismicity in Mines. Edited by Slawomir J. Gibowicz. Birkhauser Verlag Basel Boston. Berlin Seismicity in Mines Edited by Slawomir J. Gibowicz 1989 Birkhauser Verlag Basel Boston. Berlin Reprint from Pure and Applied Geophysics (PAGEOPH), Volume 129 (1989), No. 3/4 Editor's address: S~awomir

More information

Linear and Non-Linear Preconditioning

Linear and Non-Linear Preconditioning and Non- martin.gander@unige.ch University of Geneva June 2015 Invents an Iterative Method in a Letter (1823), in a letter to Gerling: in order to compute a least squares solution based on angle measurements

More information

Research Article Evaluation of the Capability of the Multigrid Method in Speeding Up the Convergence of Iterative Methods

Research Article Evaluation of the Capability of the Multigrid Method in Speeding Up the Convergence of Iterative Methods International Scholarly Research Network ISRN Computational Mathematics Volume 212, Article ID 172687, 5 pages doi:1.542/212/172687 Research Article Evaluation of the Capability of the Multigrid Method

More information

System Reduction for Nanoscale IC Design (SyreNe)

System Reduction for Nanoscale IC Design (SyreNe) System Reduction for Nanoscale IC Design (SyreNe) Peter Benner February 26, 2009 1 Introduction Since 1993, the German Federal Ministry of Education and Research (BMBF Bundesministerium füa Bildung und

More information

Dune. Patrick Leidenberger. 13. May Distributed and Unified Numerics Environment. 13. May 2009

Dune. Patrick Leidenberger. 13. May Distributed and Unified Numerics Environment. 13. May 2009 13. May 2009 leidenberger@ifh.ee.ethz.ch 1 / 25 Dune Distributed and Unified Numerics Environment Patrick Leidenberger 13. May 2009 13. May 2009 leidenberger@ifh.ee.ethz.ch 2 / 25 Table of contents 1 Introduction

More information

Josef Ballmann Rolf Jeltsch (Eds.) Nonlinear Hyperbolic Equations - Theory. Computation Methods and Applications

Josef Ballmann Rolf Jeltsch (Eds.) Nonlinear Hyperbolic Equations - Theory. Computation Methods and Applications Josef Ballmann Rolf Jeltsch (Eds.) Nonlinear Hyperbolic Equations - Theory. Computation Methods and Applications Notes on Numerical Fluid Mechanics Volume 24 Series Editors: Ernst Heinrich Hirschel, MUnchen

More information

A PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY MATCHED LAYER

A PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY MATCHED LAYER European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 A PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY

More information

A Multiscale Mortar Method And Two-Stage Preconditioner For Multiphase Flow Using A Global Jacobian Approach

A Multiscale Mortar Method And Two-Stage Preconditioner For Multiphase Flow Using A Global Jacobian Approach SPE-172990-MS A Multiscale Mortar Method And Two-Stage Preconditioner For Multiphase Flow Using A Global Jacobian Approach Benjamin Ganis, Kundan Kumar, and Gergina Pencheva, SPE; Mary F. Wheeler, The

More information

From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes

From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element

More information

A multilevel algorithm for inverse problems with elliptic PDE constraints

A multilevel algorithm for inverse problems with elliptic PDE constraints A multilevel algorithm for inverse problems with elliptic PDE constraints George Biros and Günay Doǧan School of Engineering and Applied Science, University of Pennsylvania, Philadelphia, PA 19104 E-mail:

More information

Research Article A Rapid Numerical Algorithm to Compute Matrix Inversion

Research Article A Rapid Numerical Algorithm to Compute Matrix Inversion International Mathematics and Mathematical Sciences Volume 12, Article ID 134653, 11 pages doi:.1155/12/134653 Research Article A Rapid Numerical Algorithm to Compute Matrix Inversion F. Soleymani Department

More information

Fast algorithms for the inverse medium problem. George Biros University of Pennsylvania

Fast algorithms for the inverse medium problem. George Biros University of Pennsylvania Fast algorithms for the inverse medium problem George Biros University of Pennsylvania Acknowledgments S. Adavani, H. Sundar, S. Rahul (grad students) C. Davatzikos, D. Shen, H. Litt (heart project) Akcelic,

More information

Solving Sparse Linear Systems: Iterative methods

Solving Sparse Linear Systems: Iterative methods Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccs Lecture Notes for Unit VII Sparse Matrix Computations Part 2: Iterative Methods Dianne P. O Leary c 2008,2010

More information