In April 2016 Manchester eScholar was replaced by the University of Manchester’s new Research Information Management System, Pure. In the autumn the University’s research outputs will be available to search and browse via a new Research Portal. Until then the University’s full publication record can be accessed via a temporary portal and the old eScholar content is available to search and browse via this archive.

Parallel Block Preconditioning for Multi-Physics Problems

Muddle, Richard Louden

[Thesis]. Manchester, UK: The University of Manchester; 2011.

Access to files

Abstract

In this thesis we study efficient parallel iterative solution algorithms for multi-physics problems. In particular, we consider fluid structure interaction (FSI) problems, a type of multi-physics problem in which a fluid and a deformable solid interact.All computations were performed in Oomph-Lib, a finite element library for the simulation of multi-physics problems. In Oomph-Lib, the constituent problems in a multi-physics problem are coupled monolithically, and the resulting system of non-linear equations solved with Newton's method. This requires the solution of sequences of large, sparse linear systems, for which optimal solvers are essential. The linear systems arising from the monolithic discretisation of multi-physics problems are natural candidates for solution with block-preconditioned Krylov subspace methods.We developed a generic framework for the implementation of block preconditioners within Oomph-Lib. Furthermore the framework is parallelised to facilitate the efficient solution of very large problems. This framework enables the reuse of all of Oomph-Lib's existing linear algebra infrastructure and preconditioners (including block preconditioners). We will demonstrate that a wide range of block preconditioners can be seamlessly implemented in this framework, leading to optimal iterative solvers with good parallel scaling.We concentrate on the development of an effective preconditioner for a FSI problem formulated in an arbitrary Lagrangian Eulerian (ALE) framework with pseudo-solid node updates (for the deforming fluid mesh). We begin by considering the pseudo-solid subsidiary problem; the deformation of a solid governed by equations of large displacement elasticity, subject to a prescribed boundary displacement imposed with Lagrange multiplier. We present a robust, optimal, augmented-Lagrangian type preconditioner for the resulting saddle-point linear system and prove analytically tight bounds for the spectrum of the preconditioned operator with respect to the discrete problem size.This pseudo-solid preconditioner is incorporated into a block preconditioner for the full FSI problem. One key feature of the FSI preconditioner is that existing optimal single physics preconditioners (such as the well known Navier-Stokes Least Squares Commutator preconditioner) can be employed to approximately solve the linear systems associated with the constituent sub-problems. We evaluate its performance on selected 2D and 3D problems. The preconditioner is optimal for most problems considered. In cases when sub-optimality is detected, we explain the reasons for such behavior and suggest potential improvements.

Bibliographic metadata

Type of resource:
Content type:
Form of thesis:
Type of submission:
Degree type:
Doctor of Philosophy
Degree programme:
PhD Computer Science
Publication date:
Location:
Manchester, UK
Total pages:
228
Abstract:
In this thesis we study efficient parallel iterative solution algorithms for multi-physics problems. In particular, we consider fluid structure interaction (FSI) problems, a type of multi-physics problem in which a fluid and a deformable solid interact.All computations were performed in Oomph-Lib, a finite element library for the simulation of multi-physics problems. In Oomph-Lib, the constituent problems in a multi-physics problem are coupled monolithically, and the resulting system of non-linear equations solved with Newton's method. This requires the solution of sequences of large, sparse linear systems, for which optimal solvers are essential. The linear systems arising from the monolithic discretisation of multi-physics problems are natural candidates for solution with block-preconditioned Krylov subspace methods.We developed a generic framework for the implementation of block preconditioners within Oomph-Lib. Furthermore the framework is parallelised to facilitate the efficient solution of very large problems. This framework enables the reuse of all of Oomph-Lib's existing linear algebra infrastructure and preconditioners (including block preconditioners). We will demonstrate that a wide range of block preconditioners can be seamlessly implemented in this framework, leading to optimal iterative solvers with good parallel scaling.We concentrate on the development of an effective preconditioner for a FSI problem formulated in an arbitrary Lagrangian Eulerian (ALE) framework with pseudo-solid node updates (for the deforming fluid mesh). We begin by considering the pseudo-solid subsidiary problem; the deformation of a solid governed by equations of large displacement elasticity, subject to a prescribed boundary displacement imposed with Lagrange multiplier. We present a robust, optimal, augmented-Lagrangian type preconditioner for the resulting saddle-point linear system and prove analytically tight bounds for the spectrum of the preconditioned operator with respect to the discrete problem size.This pseudo-solid preconditioner is incorporated into a block preconditioner for the full FSI problem. One key feature of the FSI preconditioner is that existing optimal single physics preconditioners (such as the well known Navier-Stokes Least Squares Commutator preconditioner) can be employed to approximately solve the linear systems associated with the constituent sub-problems. We evaluate its performance on selected 2D and 3D problems. The preconditioner is optimal for most problems considered. In cases when sub-optimality is detected, we explain the reasons for such behavior and suggest potential improvements.
Thesis main supervisor(s):
Thesis advisor(s):
Language:
en

Institutional metadata

University researcher(s):

Record metadata

Manchester eScholar ID:
uk-ac-man-scw:120215
Created by:
Muddle, Richard
Created:
18th March, 2011, 23:51:17
Last modified by:
Muddle, Richard
Last modified:
2nd November, 2011, 15:16:35

Can we help?

The library chat service will be available from 11am-3pm Monday to Friday (excluding Bank Holidays). You can also email your enquiry to us.