White Rose University Consortium logo
University of Leeds logo University of Sheffield logo York University logo

Efficient Iterative Solution Algorithms For Numerical Models of Multiphase Flow

Alrehaili, Ahlam Hamdan S (2018) Efficient Iterative Solution Algorithms For Numerical Models of Multiphase Flow. PhD thesis, University of Leeds.

Alrehaili-AHS-Computing-PhD-2018.pdf - Final eThesis - complete (pdf)
Available under License Creative Commons Attribution-Noncommercial-Share Alike 2.0 UK: England & Wales.

Download (8Mb) | Preview


This thesis is concerned with the development and application of optimally efficient numerical methods for the simulation of vascular tumour growth, based upon the multiphase fluid model introduced by Hubbard and Byrne [57]. This multiphase model involves the flow and interaction of four different, but coupled, phases which are each treated as incompressible fluids. Following a short review of models for tumour growth we describe in detail the model of Hubbard and Byrne [57], and introduce the discretization schemes used. This involves a finite volume scheme to approximate mass conservation and conforming finite element schemes to approximate momentum conservation and a reaction-diffusion equation for the background nutrient concentration. The momentum conservation system is represented as a Stokes-like flow of each phase, with source terms that reflect the phase interactions. It will be demonstrated that the solution of these coupled momentum equations, approximated using a Taylor-Hood finite element method in two dimensions, is the most computationally intensive component of the solution algorithm. The nonlinear system arising from the nutrient equation is the second most computationally expensive component. The solvers presented in this work for the discretized systems are based on preconditioned Krylov methods. Algebraic multigrid (AMG) preconditioner and a novel block preconditioner are used with Krylov methods for solving the linear systems arising from the nutrient equation at each Newton step and from the momentum equation, respectively. In each case these are shown to be very efficient algorithms: when the preconditioning strategies are applied to practical problems, the CPU time and memory are demonstrated to scale almost linearly with the problem size. Finally, the basic multiphase tumour model is extended to consider drug delivery and the inclusion of additional phases. To solve this extended model our preconditioning strategy is extended to cases with more than four phases. This is again demonstrated to perform optimally.

Item Type: Thesis (PhD)
Keywords: Multiphase system, Tumour model, Numerical simulations, Preconditioning.
Academic Units: The University of Leeds > Faculty of Engineering (Leeds) > School of Computing (Leeds)
Identification Number/EthosID: uk.bl.ethos.741195
Depositing User: Ahlam Alrehaili
Date Deposited: 14 May 2018 15:09
Last Modified: 25 Jul 2018 09:57
URI: http://etheses.whiterose.ac.uk/id/eprint/20292

You do not need to contact us to get a copy of this thesis. Please use the 'Download' link(s) above to get a copy.
You can contact us about this thesis. If you need to make a general enquiry, please see the Contact us page.

Actions (repository staff only: login required)