Skip to main content
Article thumbnail
Location of Repository

Parallel unstructured solvers for linear partial differential equations

By Dulcenéia Becker

Abstract

This thesis presents the development of a parallel algorithm to solve symmetric systems of linear equations and the computational implementation of a parallel partial differential equations solver for unstructured meshes. The proposed method, called distributive conjugate gradient - DCG, is based on a single-level domain decomposition method and the conjugate gradient method to obtain a highly scalable parallel algorithm. An overview on methods for the discretization of domains and partial differential equations is given. The partition and refinement of meshes is discussed and the formulation of the weighted residual method for two- and three-dimensions presented. Some of the methods to solve systems of linear equations are introduced, highlighting the conjugate gradient method and domain decomposition methods. A parallel unstructured PDE solver is proposed and its actual implementation presented. Emphasis is given to the data partition adopted and the scheme used for communication among adjacent subdomains is explained. A series of experiments in processor scalability is also reported. The derivation and parallelization of DCG are presented and the method validated throughout numerical experiments. The method capabilities and limitations were investigated by the solution of the Poisson equation with various source terms. The experimental results obtained using the parallel solver developed as part of this work show that the algorithm presented is accurate and highly scalable, achieving roughly linear parallel speed-up in many of the cases tested

Publisher: Cranfield University
Year: 2006
OAI identifier: oai:dspace.lib.cranfield.ac.uk:1826/4140
Provided by: Cranfield CERES

Suggested articles

Citations

  1. An accurate, flexible, parallel Poisson solver with direct substructuring for complex geometries. doi
  2. (1995). An introduction to computational fluid dynamics: the finite volume method. Harlow: Longman Scientific and Technical,
  3. (1992). Bi-CGSTAB: a fast and smoothly converging variant of Bi-CG for the solution of nonsymmetric linear systems. doi
  4. (2003). Development of a Discrete Adaptive Gridless Method for the Solution of Elliptic Partial Differential Equations.
  5. (2004). Iterative methods for large linear systems, lectures notes.
  6. (1950). Iterative Methods for Solving Partial Difference Equations of Elliptic Type. doi
  7. (1971). Iterative Solution of Large Linear Systems. doi
  8. (1962). Matrix Iterative Analysis. doi
  9. Methods for PDE-constrained optimisation with equality geometric constrains. Personal communications - second review report.
  10. Parallel adaptive unstructured finite element schemes for 3D compressible and incompressible flows. doi

To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.