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

Multilevel Numerical Algorithms for Systems of Nonlinear Parabolic Partial Differential Equations

Aljohani, Mashael Atallah (2020) Multilevel Numerical Algorithms for Systems of Nonlinear Parabolic Partial Differential Equations. PhD thesis, University of Leeds.

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

Download (1584Kb) | Preview


This thesis is concerned with the development of efficient and reliable numerical algorithms for the solution of nonlinear systems of partial differential equations (PDEs) of elliptic and parabolic type. The main focus is on the implementation and performance of three different nonlinear multilevel algorithms, following discretisation of the PDEs: the Full Approximation Scheme (FAS), Newton-Multigrid (Newton-MG) and a Newton-Krylov solver with a novel pre- conditioner that we have developed based on the use of Algebraic Multigrid (AMG). In recent years these algorithms have been commonly used to solve nonlinear systems that arise from the discretisation of PDEs due to the fact that their execution time can scale linearly (or close to linearly) with the number of degrees of freedom used in the discretisation. We consider two mathematical models: a thin film flow and the Cahn-Hilliard-Hele-Shaw model. These mathematical models consist of nonlinear, time-dependent and coupled PDEs systems. Using a Finite Difference Method (FDM) in space and Backward Differentiation For- mulae (BDF) in time, we discrete the two models, to produce nonlinear algebraic systems. We are able to solve these nonlinear systems implicitly in computationally demanding 2D situa- tions. We present numerical results, for both steady-state and time-dependent problems, that demonstrate the optimality of the three numerical algorithms for the thin film flow model. We show optimality of the FAS and Newton-Krylov approaches for the time-dependent Cahn- Hilliard-Hele-Shaw (CHHS) problem. The main contribution is to address the question of which of these three nonlinear solvers is likely to be the best (i.e. computationally most effective) in practice. In order to asses this, we discuss the careful implementation and timing of these algorithms in order to permit a fair direct comparison of their computational cost. We then present extensive numerical results in order to make this comparison between these nonlinear multilevel methods. The conclusion emerging from this investigation is that it does not appear that there is a single superior approach, but rather that the best approach is problem dependent. Specifically, we find that our optimally preconditioned Newton-Krylov approach is best for the thin film flow model in the steady-state and time-dependent form, whilst the FAS solver appears best for the time-dependent CHHS model.

Item Type: Thesis (PhD)
Keywords: Thin film flow, nonlinear system, FAS, Newton-multigrid, Preconditioned Newton-Krylov.
Academic Units: The University of Leeds > Faculty of Engineering (Leeds) > School of Computing (Leeds)
Identification Number/EthosID: uk.bl.ethos.805362
Depositing User: Dr. Mashael Atallah Aljohani
Date Deposited: 14 May 2020 06:31
Last Modified: 11 Jun 2020 09:53
URI: http://etheses.whiterose.ac.uk/id/eprint/26747

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)