In this thesis we develop 2D parallel unstructured mesh adaptationmethods for the solution of partial differential equations (PDEs) bythe finite element method (FEM). Additionally, we develop a novelblock preconditioner for the iterative solution of the linear systemsarising from the finite element discretisation of the Föpplvon Kàrmànequations.Two of the problems arising in the numerical solution of PDEs by FEMare the memory constraints that limit the solution of large problems,and the inefficiency of solving the associated linear systems bydirect or iterative solvers.We initially focus on mesh adaptation, which is a memory demandingtask of the FEM. The size of the problem increases by adding moreelements and nodes to the mesh during mesh refinement. In problemsinvolving a large number of elements, the problem size is limited bythe memory available on a single processor.In order to be able to work with large problems, we use a domaindecomposition approach to distribute the problem over multipleprocessors. One of the main objectives of this thesis is thedevelopment of 2D parallel unstructured mesh adaptation methods forthe solution of PDEs by the FEM in a variety of problems; includingdomains with curved boundaries, holes and internal boundaries. Ournewly developed methods are implemented in the software libraryoomphlib, an opensource object oriented multiphysicssoftware library implementing the FEM. We validate and demonstratetheir utility in a set of increasingly complex problems ranging fromscalar PDEs to fully coupled multiphysics problems.Having implemented and validated the infrastructure which facilitatesthe finiteelementbased discretisation of PDEs in a distributedenvironment, we shift our focus to the second problem concerning thisthesis and one of the major challenges in the computational solutionof PDEs: the solution of the large linear systems arising from theirdiscretisation.For sufficiently large problems, the solution of their associatedlinear system by direct solvers becomes impossible or inefficient,typically because of memory and time constraints. We therefore focuson preconditioned Krylov subspace methods whose efficiency dependscrucially on the provision of a good preconditioner. Thesepreconditioners are invariably problem dependent. We illustrate theirapplication and development in the solution of two elasticity problemswhich give rise to relatively large problems.First we consider the solution of a linear elasticity problem andcompute the stress distribution near a crack tip where strong localmesh refinement is required. We then consider the deformation of thinplates which are described by the nonlinear Föpplvon Kàrmànequations. A key contribution of this work is the development of anovel block preconditioner for the iterative solution of theseequations, we present the development of the preconditioner anddemonstrate its practical performance.
Date of Award  31 Dec 2016 

Original language  English 

Awarding Institution   The University of Manchester


Supervisor  Matthias Heil (Supervisor) & Andrew Hazel (Supervisor) 

 oomphlib
 Elasticity
 Preconditioner
 Föpplvon Kàrmàn
 Parallel unstructured mesh adaptation
Parallel Unstructured Mesh Adaptation and Applications
Perez Sansalvador, J. (Author). 31 Dec 2016
Student thesis: Phd