Thesis

Mixed formulations for the convection-diffusion equation

Creator
Rights statement
Awarding institution
  • University of Strathclyde
Date of award
  • 2020
Thesis identifier
  • T15495
Person Identifier (Local)
  • 201286015
Qualification Level
Qualification Name
Department, School or Faculty
Abstract
  • This thesis explores the numerical stability of the stationary Convection-Diffusion-Reaction (CDR) equation in mixed form, where the second-order equation is expressed as two first-order equations using a second variable relating to a derivative of the primary variable. This first-order system uses either a total or diffusive flux formulation. Westart by numerically testing the unstabilised Douglas and Roberts classical discretisation of the mixed CDR equation using Raviart-Thomas elements. The results indicate that,as expected, for both total and diffusive flux, the stability of the formulation degrades dramatically as diffusion decreases.Next, we investigate stabilised formulations that are designed to improve the ability of the discrete problem to cope with problems containing layers. We test the Masud and Kwack method that uses Lagrangian elements but whose analysis has not been developed.We then significantly modify the formulation to allow us to prove existence of a solution and facilitate the analysis. Our new method, which uses total flux, is then tested for convergence with standard tests and found to converge satisfactorily over a range of values of diffusion.Another family of first-order methods called First-Order System of Least-Squares (FOSLS/LSFEM) is also investigated in relation to solving the CDR equation. These symmetric,elliptic methods do not require stabilisation but also do not cope well with sharp layers and small diffusion. Modifications have been proposed and this study includes aversion of Chen et al. which uses diffusive flux, imposing boundary conditions weakly in a weighted formulation.We test our new method against all the aforementioned methods, but we find that other methods do not cope well with layers in standard tests. Our method compares favourably with the standard Streamline-Upwind-Petrov-Galerkin method (SUPG/SDFEM), but overall is not a significant improvement. With further fine-tuning, our method could improve but it has more computational overhead than SUPG.
Advisor / supervisor
  • Barrenechea, Gabriel
Resource Type
DOI
Date Created
  • 2020
Former identifier
  • 9912789192302996

Relations

Items