Numerical Solution of the Generalised Poisson Equation on Parallel Computers
|
|
- Hector Fleming
- 6 years ago
- Views:
Transcription
1 Numerical Solution of the Generalised Poisson Equation on Parallel Computers Diplomanden- und Dissertantenkolloquium Universität Wien Hannes Grimm-Strele Projektgruppe ACORE, Fakultät für Mathematik, Universität Wien April 15, 2010 Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
2 Introduction Subject of this talk Generalised Poisson Equation (GPE: Main questions considered: (κ u + cu = f on Ω, (1 κ(x κ 0 > 0, c(x 0 x Ω. (2 Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
3 Introduction Subject of this talk Generalised Poisson Equation (GPE: (κ u + cu = f on Ω, (1 κ(x κ 0 > 0, c(x 0 x Ω. (2 Main questions considered: Origin and practical relevance of GPE Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
4 Introduction Subject of this talk Generalised Poisson Equation (GPE: (κ u + cu = f on Ω, (1 κ(x κ 0 > 0, c(x 0 x Ω. (2 Main questions considered: Origin and practical relevance of GPE Importance of parallel computing Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
5 Introduction Subject of this talk Generalised Poisson Equation (GPE: (κ u + cu = f on Ω, (1 κ(x κ 0 > 0, c(x 0 x Ω. (2 Main questions considered: Origin and practical relevance of GPE Importance of parallel computing Approach to numerical solution on parallel computers Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
6 Heat conduction Physical Background Assumption: Particles move according to Brownian motion. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
7 Heat conduction Physical Background Assumption: Particles move according to Brownian motion. Fourier s law: Energy flux goes from regions of high temperature to regions of low temperature. J = κ T (3 T temperature κ heat conductivity of the material (> 0 e energy J energy flux Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
8 Heat conduction Physical Background Assumption: Particles move according to Brownian motion. Fourier s law: Energy flux goes from regions of high temperature to regions of low temperature. J = κ T (3 T temperature κ heat conductivity of the material (> 0 e energy J energy flux Conservation of energy e t = J Stationarity e t = 0 J = (κ T = 0 (4 Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
9 Euler Equations 1/2 Physical Background The Euler Equations govern the dynamics of flow without friction. Mass and momentum equations in two dimensions: ( ρu ρv t ( ρu 2 + ρuv x ( ρuv + ρuv 2 ρ t + (ρ u = 0, (5 ( px + = 0, (6 y p y ρ u = (u, v p mass density velocity field pressure Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
10 Euler Equations 1/2 Physical Background The Euler Equations govern the dynamics of flow without friction. Mass and momentum equations in two dimensions: ( ρu ρv t ( ρu 2 + ρuv x ( ρuv + ρuv 2 ρ t + (ρ u = 0, (5 ( px + = 0, (6 y p y ρ u = (u, v p mass density velocity field pressure Incompressibility assumption: Mass density does not change along a trajectory. Flow is incompressible (5 u = 0. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
11 Euler Equations 2/2 Physical Background Set (ρ u = (ρ u n t ( ( ρu 2 ρuv x ( ρuv + ρuv 2 y Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
12 Euler Equations 2/2 Physical Background Set (ρ u = (ρ u n t ( ( ρu 2 ρuv x ( ρuv + ρuv 2 y Euler forward step (ρ u t (ρ un+1 (ρ u t = p Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
13 Euler Equations 2/2 Physical Background Set (ρ u = (ρ u n t ( ( ρu 2 ρuv x ( ρuv + ρuv 2 y Euler forward step (ρ u t (ρ un+1 (ρ u t = p Divide by ρ n+1 = ρ un+1 u t = p ρ n+1 Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
14 Physical Background Euler Equations 2/2 Set (ρ u = (ρ u n t ( ( ρu 2 ρuv x ( ρuv + ρuv 2 y Euler forward step (ρ u t (ρ un+1 (ρ u t = p Divide by ρ n+1 = ρ un+1 u t = p ρ n+1 ( Take divergence u n+1 = u t p ρ n+1 Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
15 Physical Background Euler Equations 2/2 Set (ρ u = (ρ u n t ( ( ρu 2 ρuv x ( ρuv + ρuv 2 y Euler forward step (ρ u t (ρ un+1 (ρ u t = p Divide by ρ n+1 = ρ un+1 u t = p ρ n+1 ( Take divergence u n+1 = u t p ρ n+1 ( Incompressible flow p = u ρ n+1 t Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
16 Parallel Computing Why parallel computing is important Simulating 3s of the solar granulation with ANTARES requires approx CPUh and 1.1 TB memory. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
17 Parallel Computing Example: Deutsches Klimarechenzentrum SC1: Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
18 Parallel Computing Example: Deutsches Klimarechenzentrum SC1: Die globale Umwelt und insbesondere das Klima sind äußerst komplexe Systeme, deren Dynamik und künftige Entwicklung nur mit weitreichenden Untersuchungen und aufwendigen Modellrechnungen erfasst werden können. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
19 Parallel Computing Example: Deutsches Klimarechenzentrum SC1: Die globale Umwelt und insbesondere das Klima sind äußerst komplexe Systeme, deren Dynamik und künftige Entwicklung nur mit weitreichenden Untersuchungen und aufwendigen Modellrechnungen erfasst werden können. Vier Gründe, warum wir auf Höchstleistungsrechner angewiesen sind: Komplexität des Erdsystems Modellauflösung Ensemble -Rechnungen Integrationen über viele Jahrhunderte Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
20 Parallel Computing Example: Deutsches Klimarechenzentrum SC1: Die globale Umwelt und insbesondere das Klima sind äußerst komplexe Systeme, deren Dynamik und künftige Entwicklung nur mit weitreichenden Untersuchungen und aufwendigen Modellrechnungen erfasst werden können. Vier Gründe, warum wir auf Höchstleistungsrechner angewiesen sind: Komplexität des Erdsystems Modellauflösung Ensemble -Rechnungen Integrationen über viele Jahrhunderte HPC blizzard : 8448 cores, 158 TFlops, 20 TB memory Cp: VSC 3488 cores, 35.5 TFlops, 11.2 TB memory Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
21 Parallel Computing Example: Sauber Motorsport CFD (Computational Fluid Dynamics spielt insbesondere bei der Entwicklung von Front- und Heckflügeln sowie auch bei der Motor- und Bremskühlung eine wichtige Rolle. Seit dem Frühjahr 2008 ist Albert3 in Betrieb. Die jüngste Ausbaustufe verfügt über insgesamt 4224 Prozessorkerne. [...] Der Arbeitsspeicher wuchs auf 8448 GByte und die maximale Rechenleistung auf 50,7 TFlops. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
22 Why not just wait? Parallel Computing Moore sches Gesetz: Die Anzahl an Transistoren auf einem handelsüblichen Prozessor verdoppelt sich alle achtzehn Monate. Quelle: Der Spiegel 11/2005, S Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
23 Parallel Computing The concept of parallel programming Decomposition and distribution of data and work (e.g. by domain decomposition Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
24 Parallel Computing The concept of parallel programming Decomposition and distribution of data and work (e.g. by domain decomposition Synchronisation and communication Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
25 Parallel Computing The concept of parallel programming Decomposition and distribution of data and work (e.g. by domain decomposition Synchronisation and communication Criteria for good parallelisation: load balancing and optimal speedup (Amdahl s law Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
26 Parallel Computing The concept of parallel programming Decomposition and distribution of data and work (e.g. by domain decomposition Synchronisation and communication Criteria for good parallelisation: load balancing and optimal speedup (Amdahl s law Programming models: Message Passing Interface (MPI and OpenMP Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
27 Numerical Methods Finite Element Method (FEM Using FEM, GPE is transformed to the discrete system Au h = b, where u h is an approximate solution for u. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
28 Numerical Methods Finite Element Method (FEM Using FEM, GPE is transformed to the discrete system Au h = b, where u h is an approximate solution for u. The approximation error can be controlled by the choice of the finite-dimensional ansatz space. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
29 Numerical Methods Finite Element Method (FEM Using FEM, GPE is transformed to the discrete system Au h = b, where u h is an approximate solution for u. The approximation error can be controlled by the choice of the finite-dimensional ansatz space. Most times: space of linear splines on Ω. Figure: Basis function for the space of linear splines (2D. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
30 Numerical Methods Finite Element Method (FEM Using FEM, GPE is transformed to the discrete system Au h = b, where u h is an approximate solution for u. The approximation error can be controlled by the choice of the finite-dimensional ansatz space. Most times: space of linear splines on Ω. Figure: Basis function for the space of linear splines (2D. A is always symmetric and positive definite. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
31 Numerical Methods Preconditioned Conjugate Gradient Algorithm (PCG Iterative algorithm to invert linear systems Au h = b. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
32 Numerical Methods Preconditioned Conjugate Gradient Algorithm (PCG Iterative algorithm to invert linear systems Au h = b. Converges if A is symmetric and positive definite. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
33 Numerical Methods Preconditioned Conjugate Gradient Algorithm (PCG Iterative algorithm to invert linear systems Au h = b. Converges if A is symmetric and positive definite. Convergence speed depends on κ(a. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
34 Numerical Methods Preconditioned Conjugate Gradient Algorithm (PCG Iterative algorithm to invert linear systems Au h = b. Converges if A is symmetric and positive definite. Convergence speed depends on κ(a. κ(a can be lowered by preconditioning (e.g. Incomplete Cholesky Decomposition. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
35 Numerical Methods Schur Complement Method 1/2 If the matrix A is of the following form: ( A1,1 A A = 1,2 A 2,1 A 2,2 with certain submatrices A 1,1, A 1,2, A 2,1 and A 2,2 such that A 1,1 is invertible, one can define the Schur Complement S A by S A = A 2,2 A 2,1 A 1 1,1 A 1,2. (7 Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
36 Numerical Methods Schur Complement Method 1/2 If the matrix A is of the following form: ( A1,1 A A = 1,2 A 2,1 A 2,2 with certain submatrices A 1,1, A 1,2, A 2,1 and A 2,2 such that A 1,1 is invertible, one can define the Schur Complement S A by Lemma S A = A 2,2 A 2,1 A 1 1,1 A 1,2. (7 If A is symmetric and positive definite, then S A defined by (7 is also symmetric and positive definite. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
37 Numerical Methods Schur Complement Method 2/2 The system Au = b is equivalent to ( ( ( ( I 0 A1,1 A 1,2 u1 b1 A 2,1 A 1 =. 1,1 I 0 S A u 2 b 2 By the PCG algorithm these systems can be solved subsequently. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14
38 Numerical Methods Schur Complement Method 2/2 The system Au = b is equivalent to ( I 0 A 2,1 A 1 1,1 I ( A1,1 A 1,2 0 S A ( u1 u 2 = ( b1 By the PCG algorithm these systems can be solved subsequently. This can be done in parallel, if A 1,1 has block structure (can be achieved by reordering and renumbering of the nodes: A 1,1 = D D p. Then A 1,2 = (B 1,..., B p T, A 2,1 = (C 1,..., C p, A 2,2 = i E i. S A = p i=1 (E i C i D 1 i B i. Hannes Grimm-Strele (ACORE Parallel Solution of GPE April 15, / 14 b 2.
Contents. Preface... xi. Introduction...
Contents Preface... xi Introduction... xv Chapter 1. Computer Architectures... 1 1.1. Different types of parallelism... 1 1.1.1. Overlap, concurrency and parallelism... 1 1.1.2. Temporal and spatial parallelism
More informationCloud Computing mit mathematischen Anwendungen
Cloud Computing mit mathematischen Anwendungen Dr. habil. Marcel Kunze Engineering Mathematics and Computing Lab (EMCL) Institut für Angewandte und Numerische Mathematik IV Karlsruhe Institute of Technology
More informationMehrphasensimulationen zur Charakterisierung einer Aerosol-Teilprobenahme
Mehrphasensimulationen zur Charakterisierung einer Aerosol-Teilprobenahme G. Lindner 8.41, Modellierung und Simulation February 29, 2016 Review CFD Study for the Design of the Dilution Section The mixing
More informationMaster 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 informationOrdinals and Cardinals: Basic set-theoretic techniques in logic
Ordinals and Cardinals: Basic set-theoretic techniques in logic Benedikt Löwe Universiteit van Amsterdam Grzegorz Plebanek Uniwersytet Wroc lawski ESSLLI 2011, Ljubljana, Slovenia This course is a foundational
More informationHierarchical Parallel Solution of Stochastic Systems
Hierarchical Parallel Solution of Stochastic Systems Second M.I.T. Conference on Computational Fluid and Solid Mechanics Contents: Simple Model of Stochastic Flow Stochastic Galerkin Scheme Resulting Equations
More informationIncomplete Cholesky preconditioners that exploit the low-rank property
anapov@ulb.ac.be ; http://homepages.ulb.ac.be/ anapov/ 1 / 35 Incomplete Cholesky preconditioners that exploit the low-rank property (theory and practice) Artem Napov Service de Métrologie Nucléaire, Université
More informationNumerical modelling of phase change processes in clouds. Challenges and Approaches. Martin Reitzle Bernard Weigand
Institute of Aerospace Thermodynamics Numerical modelling of phase change processes in clouds Challenges and Approaches Martin Reitzle Bernard Weigand Introduction Institute of Aerospace Thermodynamics
More informationD-optimally Lack-of-Fit-Test-efficient Designs and Related Simple Designs
AUSTRIAN JOURNAL OF STATISTICS Volume 37 (2008), Number 3&4, 245 253 D-optimally Lack-of-Fit-Test-efficient Designs and Related Simple Designs Wolfgang Bischoff Catholic University ichstätt-ingolstadt,
More informationLocating RFID Tags. Matthias Benesch, Siegfried Bublitz Siemens AG. C-LAB Report. Vol. 8 (2009) No. 3
Locating RFID Tags Matthias Benesch, Siegfried Bublitz Siemens AG C-LAB Report Vol. 8 (009) No. 3 Cooperative Computing & Communication Laboratory ISSN 1619-7879 C-LAB ist eine Kooperation der Universität
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationAlgebra. Übungsblatt 10 (Lösungen)
Fakultät für Mathematik Sommersemester 2017 JProf. Dr. Christian Lehn Dr. Alberto Castaño Domínguez Algebra Übungsblatt 10 (Lösungen) Aufgabe 1. Es sei k ein Körper. Man zeige, dass es in k[x] unendlich
More informationInterpolation in h-version finite element spaces
Interpolation in h-version finite element spaces Thomas Apel Institut für Mathematik und Bauinformatik Fakultät für Bauingenieur- und Vermessungswesen Universität der Bundeswehr München Chemnitzer Seminar
More informationPDE 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 informationIterative Solvers in the Finite Element Solution of Transient Heat Conduction
Iterative Solvers in the Finite Element Solution of Transient Heat Conduction Mile R. Vuji~i} PhD student Steve G.R. Brown Senior Lecturer Materials Research Centre School of Engineering University of
More informationModule Module Component Aim LoI** Study programme Term expected calendar week/ time***
Call for Application Erasmus + program "Staff Training" Modules offered for participants at the Eberswalde University for Sustainable Development (EUSD) during winter term 2018* Module Module Component
More informationRobust 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 informationAlgebra. Übungsblatt 12 (Lösungen)
Fakultät für Mathematik Sommersemester 2017 JProf Dr Christian Lehn Dr Alberto Castaño Domínguez Algebra Übungsblatt 12 (Lösungen) Aufgabe 1 Berechnen Sie Minimalpolynome f i, i = 1,, 4, der folgenden
More informationOn 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 informationDiscrete 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 informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More informationNumerically Solving Partial Differential Equations
Numerically Solving Partial Differential Equations Michael Lavell Department of Applied Mathematics and Statistics Abstract The physics describing the fundamental principles of fluid dynamics can be written
More informationSS BMMM01 Basismodul Mathematics/Methods Block 1: Mathematics for Economists. Prüfer: Prof. Dr.
SS 2018 02.06.2018 1289BMMM01 Basismodul Mathematics/Methods Block 1: Mathematics for Economists Prüfer: Prof. Dr. Rainer Dyckerhoff Bitte füllen Sie die nachfolgende Zeile aus! Matrikelnummer (student
More informationJae Heon Yun and Yu Du Han
Bull. Korean Math. Soc. 39 (2002), No. 3, pp. 495 509 MODIFIED INCOMPLETE CHOLESKY FACTORIZATION PRECONDITIONERS FOR A SYMMETRIC POSITIVE DEFINITE MATRIX Jae Heon Yun and Yu Du Han Abstract. We propose
More informationwhen viewed from the top, the objects should move as if interacting gravitationally
2 Elastic Space 2 Elastic Space The dynamics and apparent interactions of massive balls rolling on a stretched horizontal membrane are often used to illustrate gravitation. Investigate the system further.
More informationJ.I. Aliaga 1 M. Bollhöfer 2 A.F. Martín 1 E.S. Quintana-Ortí 1. March, 2009
Parallel Preconditioning of Linear Systems based on ILUPACK for Multithreaded Architectures J.I. Aliaga M. Bollhöfer 2 A.F. Martín E.S. Quintana-Ortí Deparment of Computer Science and Engineering, Univ.
More informationAlgebra. Übungsblatt 8 (Lösungen) m = a i m i, m = i=1
Fakultät für Mathematik Sommersemester 2017 JProf. Dr. Christian Lehn Dr. Alberto Castaño Domínguez Algebra Übungsblatt 8 (Lösungen) Aufgabe 1. Es seien R ein Ring, m R ein maximales Ideal und M ein R-Modul.
More informationHalbeinfache Algebren
Prof. Dr. M. Rapoport Seminar (S2A1), WS 2016/17 Halbeinfache Algebren Ziel des Seminars ist es, halbeinfache Algebren über Körpern zu studieren. Beispiele sind Matrizenringe oder Schiefkörper wie die
More information1. Einleitung. 1.1 Organisatorisches. Ziel der Vorlesung: Einführung in die Methoden der Ökonometrie. Voraussetzungen: Deskriptive Statistik
1. Einleitung 1.1 Organisatorisches Ziel der Vorlesung: Einführung in die Methoden der Ökonometrie Voraussetzungen: Deskriptive Statistik Wahrscheinlichkeitsrechnung und schließende Statistik Fortgeschrittene
More information1. Positive and regular linear operators.
1. Positive and regular linear operators. Objekttyp: Chapter Zeitschrift: L'Enseignement Mathématique Band (Jahr): 19 (1973) Heft 3-4: L'ENSEIGNEMENT MATHÉMATIQUE PDF erstellt am: 1.0.018 Nutzungsbedingungen
More informationA 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 informationParallelization of Multilevel Preconditioners Constructed from Inverse-Based ILUs on Shared-Memory Multiprocessors
Parallelization of Multilevel Preconditioners Constructed from Inverse-Based ILUs on Shared-Memory Multiprocessors J.I. Aliaga 1 M. Bollhöfer 2 A.F. Martín 1 E.S. Quintana-Ortí 1 1 Deparment of Computer
More informationCyclic difference sets with parameters (511, 255, 127) BACHER, Roland
Article Cyclic difference sets with parameters (511, 255, 127) BACHER, Roland Reference BACHER, Roland. Cyclic difference sets with parameters (511, 255, 127). Enseignement mathématique, 1994, vol. 40,
More informationSPARSE SOLVERS POISSON EQUATION. Margreet Nool. November 9, 2015 FOR THE. CWI, Multiscale Dynamics
SPARSE SOLVERS FOR THE POISSON EQUATION Margreet Nool CWI, Multiscale Dynamics November 9, 2015 OUTLINE OF THIS TALK 1 FISHPACK, LAPACK, PARDISO 2 SYSTEM OVERVIEW OF CARTESIUS 3 POISSON EQUATION 4 SOLVERS
More informationOn deflation and singular symmetric positive semi-definite matrices
Journal of Computational and Applied Mathematics 206 (2007) 603 614 www.elsevier.com/locate/cam On deflation and singular symmetric positive semi-definite matrices J.M. Tang, C. Vuik Faculty of Electrical
More informationGrid Forming Control for Stable Power Systems with. up to 100 % Inverter Based Generation: A Paradigm. Scenario Using the IEEE 118-Bus System
Grid Forming Control for Stable Power Systems with up to 100 % Inverter Based Generation: Paradigm Scenario Using the IEEE 118-Bus System 17 Oct 2018 Mario Ndreko, Sven Rüberg, Wilhelm Winter Power System
More informationFinite Element Decompositions for Stable Time Integration of Flow Equations
MAX PLANCK INSTITUT August 13, 2015 Finite Element Decompositions for Stable Time Integration of Flow Equations Jan Heiland, Robert Altmann (TU Berlin) ICIAM 2015 Beijing DYNAMIK KOMPLEXER TECHNISCHER
More informationAn introduction to PDE-constrained optimization
An introduction to PDE-constrained optimization Wolfgang Bangerth Department of Mathematics Texas A&M University 1 Overview Why partial differential equations? Why optimization? Examples of PDE optimization
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationNumerical Solution of Partial Differential Equations governing compressible flows
Numerical Solution of Partial Differential Equations governing compressible flows Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
More informationOUTLINE 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 informationMultipole-Based Preconditioners for Sparse Linear Systems.
Multipole-Based Preconditioners for Sparse Linear Systems. Ananth Grama Purdue University. Supported by the National Science Foundation. Overview Summary of Contributions Generalized Stokes Problem Solenoidal
More informationThe 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 information4. Minkowski vertices and Chebyshev Systems
4. Minkowski vertices and Chebyshev Systems Objekttyp: Chapter Zeitschrift: L'Enseignement Mathématique Band (Jahr): 43 (1997) Heft 1-2: L'ENSEIGNEMENT MATHÉMATIQUE PDF erstellt am: 18.06.2018 Nutzungsbedingungen
More informationAlgebra. Übungsblatt 5 (Lösungen)
Fakultät für Mathematik Sommersemester 2017 JProf. Dr. Christian Lehn Dr. Alberto Castaño Domínguez Algebra Übungsblatt 5 (Lösungen) Aufgabe 1. Für eine Gruppe G und ein g G sei der Zentralisator von g
More informationCourse Okt 28, Statistische Mechanik plus. Course Hartmut Ruhl, LMU, Munich. Personen. Hinweis. Rationale
Okt 28, 2016 ASC, room A239, phone 089-21804210, email hartmut.ruhl@lmu.de. Patrick Böhl, ASC, room A205, phone 089-21804640, email patrick.boehl@physik.uni-muenchen.de. Das folgende Material ist nicht
More informationModelling and implementation of algorithms in applied mathematics using MPI
Modelling and implementation of algorithms in applied mathematics using MPI Lecture 3: Linear Systems: Simple Iterative Methods and their parallelization, Programming MPI G. Rapin Brazil March 2011 Outline
More informationDirichlet-Neumann and Neumann-Neumann Methods
Dirichlet-Neumann and Neumann-Neumann Methods Felix Kwok Hong Kong Baptist University Introductory Domain Decomposition Short Course DD25, Memorial University of Newfoundland July 22, 2018 Outline Methods
More informationScalable Non-blocking Preconditioned Conjugate Gradient Methods
Scalable Non-blocking Preconditioned Conjugate Gradient Methods Paul Eller and William Gropp University of Illinois at Urbana-Champaign Department of Computer Science Supercomputing 16 Paul Eller and William
More informationA Framework for Hybrid Parallel Flow Simulations with a Trillion Cells in Complex Geometries
A Framework for Hybrid Parallel Flow Simulations with a Trillion Cells in Complex Geometries SC13, November 21 st 2013 Christian Godenschwager, Florian Schornbaum, Martin Bauer, Harald Köstler, Ulrich
More informationSolving Large Nonlinear Sparse Systems
Solving Large Nonlinear Sparse Systems Fred W. Wubs and Jonas Thies Computational Mechanics & Numerical Mathematics University of Groningen, the Netherlands f.w.wubs@rug.nl Centre for Interdisciplinary
More informationThe Navier-Stokes Equations
s University of New Hampshire February 22, 202 and equations describe the non-relativistic time evolution of mass and momentum in fluid substances. mass density field: ρ = ρ(t, x, y, z) velocity field:
More informationNumerik 2 Motivation
Numerik 2 Motivation P. Bastian Universität Heidelberg Interdisziplinäres Zentrum für Wissenschaftliches Rechnen Im Neuenheimer Feld 368, D-69120 Heidelberg email: Peter.Bastian@iwr.uni-heidelberg.de April
More informationMixed Finite Elements Method
Mixed Finite Elements Method A. Ratnani 34, E. Sonnendrücker 34 3 Max-Planck Institut für Plasmaphysik, Garching, Germany 4 Technische Universität München, Garching, Germany Contents Introduction 2. Notations.....................................
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationThe notable configuration of inscribed equilateral triangles in a triangle
Elem. Math. 57 00) 6 1 001-6018/0/01006-6 $ 1.50+0.0/0 c irkhäuser Verlag, asel, 00 Elemente der Mathematik The notable configuration of inscribed equilateral triangles in a triangle las Herrera Gómez
More informationInternational Steam Tables - Properties of Water and Steam based on
International Steam Tables - Properties of Water and Steam based on the Industrial Formulation IAPWS-IF97 Tables, Algorithms, Diagrams, and CD-ROM Electronic Steam Tables - All of the equations of IAPWS-IF97
More informationIntroduction to PDEs and Numerical Methods Lecture 1: Introduction
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 1: Introduction Dr. Noemi Friedman, 28.10.2015. Basic information on the course Course
More informationThe Euler equations in fluid mechanics
The Euler equations in fluid mechanics Jordan Bell jordan.bell@gmail.com Department of Mathematics, niversity of Toronto April 14, 014 1 Continuity equation Let Ω be a domain in R n and let C (Ω R); perhaps
More informationDetailed 3D modelling of mass transfer processes in two phase flows with dynamic interfaces
Detailed 3D modelling of mass transfer processes in two phase flows with dynamic interfaces D. Darmana, N.G. Deen, J.A.M. Kuipers Fundamentals of Chemical Reaction Engineering, Faculty of Science and Technology,
More informationAdjoint-Fueled Advances in Error Estimation for Multiscale, Multiphysics Systems
Adjoint-Fueled Advances in Error Estimation for Multiscale, Multiphysics Systems Donald Estep University Interdisciplinary Research Scholar Department of Statistics Department of Mathematics Center for
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationOrganische Chemie IV: Organische Photochemie
rganische Chemie IV: rganische Photochemie Sommersemester 2006 Technische Universität München Klausur am 04.08.2006 ame; Vorname... Matrikel-r.... (Druckbuchstaben) geboren am... in...... (Eigenhändige
More informationParallelization of An Experimental Multiphase Flow Algorithm. Master Thesis Literature Review submitted to the Delft Institute of Applied Mathematics
Delft University of Technology Faculty of Electrical Engineering, Mathematics and Computer Science Delft Institute of Applied Mathematics Parallelization of An Experimental Multiphase Flow Algorithm Master
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationFEMLAB as a General Tool to Investigate the Basic Laws of Physics
FEMLAB as a General Tool to Investigate the Basic Laws of Physics Kurt Bräuer Universität Tübingen, Institut für Theoretische Physik, Auf der Morgenstelle 14, D-7076 Tübingen e-mail: kurt.braeuer@uni-tuebinen.de,
More informationOverview of Convection Heat Transfer
Overview of Convection Heat Transfer Lecture Notes for ME 448/548 Gerald Recktenwald March 16, 2006 1 Introduction CFD is a powerful tool for engineering design involving fluid flow and heat transfer.
More informationChapter 9 Implicit integration, incompressible flows
Chapter 9 Implicit integration, incompressible flows The methods we discussed so far work well for problems of hydrodynamics in which the flow speeds of interest are not orders of magnitude smaller than
More informationAnalysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems
Analysis of Hybrid Discontinuous Galerkin Methods for Incompressible Flow Problems Christian Waluga 1 advised by Prof. Herbert Egger 2 Prof. Wolfgang Dahmen 3 1 Aachen Institute for Advanced Study in Computational
More informationCourse Nov 03, Statistische Mechanik plus. Course Hartmut Ruhl, LMU, Munich. People involved. Hinweis.
Nov 03, 016 ASC, room A 38, phone 089-180410, email hartmut.ruhl@lmu.de Patrick Böhl, ASC, room A05, phone 089-1804640, email patrick.boehl@physik.uni-muenchen.de. Das folgende Material ist nicht prüfungsrelevant.
More informationStatic Program Analysis
Static Program Analysis Thomas Noll Software Modeling and Verification Group RWTH Aachen University https://moves.rwth-aachen.de/teaching/ws-1617/spa/ Software Architektur Praxis-Workshop Bringen Sie Informatik
More informationEF2000 Control Laws - Phase 4 Optimisation of Feedback Gains (Pitch)
Optimisation of Feedback Gains (Pitch) Some DA & IPA are flying with the Phase 4 Claws Pilot comments: (improved Handling Qualities/Agility (flying like a true fighter aircraft Optimisation program to
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More information3. Rational curves on quartics in $P^3$
3 Rational curves on quartics in $P^3$ Objekttyp Chapter Zeitschrift L'Enseignement Mathématique Band (ahr) 43 (997) Heft 3-4 L'ENSEIGNEMENT MATHÉMATIQUE PDF erstellt am 2904208 Nutzungsbedingungen Die
More informationFriction Temperature of POM PE Sliding Contacts
Friction Temperature of POM PE Sliding Contacts Ralf Bartsch, Jens Sumpf, André Bergmann Institut für Fördertechnik und Kunststoffe, Technische Universität Chemnitz, Reichenhainer Straße 70, 09126 Chemnitz,
More informationStatic Program Analysis
Static Program Analysis Thomas Noll Software Modeling and Verification Group RWTH Aachen University https://moves.rwth-aachen.de/teaching/ws-1617/spa/ Software Architektur Praxis-Workshop Bringen Sie Informatik
More information3-D FINITE ELEMENT MODELLING OF GRANULAR FLOW IN SILOS
3-D FINITE ELEMENT MODELLING OF GRANULAR FLOW IN SILOS Guenter A. Rombach and Frank Neumann 2 ABSTRACT A continuum approach based on the Finite Element Method is presented in this paper with which the
More informationSIMPLE Algorithm for Two-Dimensional Channel Flow. Fluid Flow and Heat Transfer
SIMPLE Algorithm for Two-Dimensional Channel Flow Fluid Flow and Heat Transfer by Professor Jung-Yang San Mechanical Engineering Department National Chung Hsing University Two-dimensional, transient, incompressible
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More information1. The Bergman, Kobayashi-Royden and Caratheodory metrics
1. The Bergman, Kobayashi-Royden and Caratheodory metrics Objekttyp: Chapter Zeitschrift: L'Enseignement Mathématique Band (Jahr): 33 (1987) Heft 1-2: L'ENSEIGNEMENT MATHÉMATIQUE PDF erstellt am: 18.06.2018
More informationNumerical 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 informationADDITIVE 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 informationThe Apollonius Circle as a Tucker Circle
Forum Geometricorum Volume 2 (2002) 175 182 FORUM GEOM ISSN 1534-1178 The Apollonius Circle as a Tucker Circle Darij Grinberg and Paul Yiu Abstract We give a simple construction of the circular hull of
More informationHydro-elastic Wagner impact using variational inequalities
Hydro-elastic Wagner impact using variational inequalities Thomas GAZZOLA, Alexander KOROBKIN, Šime MALENICA Introduction A simple model of water impact has been introduced by Wagner [6]. This model is
More informationPerformance Analysis of Parallel Alternating Directions Algorithm for Time Dependent Problems
Performance Analysis of Parallel Alternating Directions Algorithm for Time Dependent Problems Ivan Lirkov 1, Marcin Paprzycki 2, and Maria Ganzha 2 1 Institute of Information and Communication Technologies,
More informationObservation and Control for Operator Semigroups
Birkhäuser Advanced Texts Basler Lehrbücher Observation and Control for Operator Semigroups Bearbeitet von Marius Tucsnak, George Weiss Approx. 496 p. 2009. Buch. xi, 483 S. Hardcover ISBN 978 3 7643 8993
More informationFree energy concept Free energy approach LBM implementation Parameters
BINARY LIQUID MODEL A. Kuzmin J. Derksen Department of Chemical and Materials Engineering University of Alberta Canada August 22,2011 / LBM Workshop OUTLINE 1 FREE ENERGY CONCEPT 2 FREE ENERGY APPROACH
More informationMatroids on graphs. Brigitte Servatius Worcester Polytechnic Institute. First Prev Next Last Go Back Full Screen Close Quit
on K n on graphs Brigitte Servatius Worcester Polytechnic Institute Page 1 of 35 on K n Page 2 of 35 1. Whitney [9] defined a matroid M on a set E: M = (E, I) E is a finite set I is a collection of subsets
More informationNumerical Modeling of Inclined Negatively Buoyant Jets
Numerical Modeling of Inclined Negatively Buoyant Jets Presentation by: Hossein Kheirkhah Graduate Student in Civil Engineering Dep. of Civil Engineering University of Ottawa CANADA ICDEMOS April2014 Outline
More informationImplementation of Implicit Solution Techniques for Non-equilibrium Hypersonic Flows
Short Training Program Report Implementation of Implicit Solution Techniques for Non-equilibrium Hypersonic Flows Julian Koellermeier RWTH Aachen University Supervisor: Advisor: Prof. Thierry Magin von
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationApplications of Mathematical Economics
Applications of Mathematical Economics Michael Curran Trinity College Dublin Overview Introduction. Data Preparation Filters. Dynamic Stochastic General Equilibrium Models: Sunspots and Blanchard-Kahn
More informationMPI 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 information3. Clifford's Theorem The elementary proof
3. Clifford's Theorem The elementary proof Objekttyp: Chapter Zeitschrift: L'Enseignement Mathématique Band (Jahr): 30 (1984) Heft 1-2: L'ENSEIGNEMENT MATHÉMATIQUE PDF erstellt am: 03.12.2018 Nutzungsbedingungen
More informationA note on Steiner symmetrization of hyperbolic triangles
Elem. Math. 58 (2003) 21 25 0013-6018/03/010021-5 c Birkhäuser Verlag, Basel, 2003 Elemente der Mathematik A note on Steiner symmetrization of hyperbolic triangles François Guéritaud François Guéritaud
More information7.4 The Saddle Point Stokes Problem
346 CHAPTER 7. APPLIED FOURIER ANALYSIS 7.4 The Saddle Point Stokes Problem So far the matrix C has been diagonal no trouble to invert. This section jumps to a fluid flow problem that is still linear (simpler
More informationUniversität Dortmund UCHPC. Performance. Computing for Finite Element Simulations
technische universität dortmund Universität Dortmund fakultät für mathematik LS III (IAM) UCHPC UnConventional High Performance Computing for Finite Element Simulations S. Turek, Chr. Becker, S. Buijssen,
More informationCapacitance Matrix Method
Capacitance Matrix Method Marcus Sarkis New England Numerical Analysis Day at WPI, 2019 Thanks to Maksymilian Dryja, University of Warsaw 1 Outline 1 Timeline 2 Capacitance Matrix Method-CMM 3 Applications
More informationScenario Aggregation
Scenario Aggregation for Solvency Regulation Mathieu Cambou 1 (joint with Damir Filipović) Ecole Polytechnique Fédérale de Lausanne 1 SCOR ASTIN Colloquium The Hague, 22 May 2013 Outline Scenario Aggregation
More information