Date of Award

8-2011

Degree Type

Dissertation

Degree Name

Doctor of Philosophy

Major

Mathematics

Major Professor

Ohannes A. Karakashian

Committee Members

Steven M. Wise, Xiaobing H. Feng, Gao Yanfei

Abstract

This PhD dissertation concentrates on the development and application of adaptive Discontinuous Galerkin Finite Element (DG-FE) methods for the numerical solution of a Cahn-Hilliard-type diffuse interface model for biological growth. Models of this type have become popular for studying cancerous tumor progression in vivo. The work in this dissertation advances the state-of-the-art in the following ways: To our knowledge the work here contains the first primitive-variable, completely discontinuous numerical implementations of a 2D scheme for the Cahn-Hilliard equation as well as a di ffuse interface model of cancer growth. We provide numerical evidence that the schemes above are convergent, with the optimal order. The efficiency of the numerical algorithms depends largely on the implementation of fast solvers for the systems of equations resulting from the DG-FE discretizations. We have developed such capabilities based on multigrid and sparse direct solver techniques. We demonstrate proof-of-concept regarding the implementation of a practical spatially adaptive meshing algorithm for the numerical schemes just mentioned and the effective use of a very simple, but powerful, marking strategy based on an inverse estimate. We demonstrate proof-of-concept for a novel simplified diffuse interface model of tumor growth. This model is essentially the Cahn-Hilliard equation with an added source term that is specialized for the context of cancerous tumor progression. We devise and analyze a mixed DG-FE scheme of convex splitting (CS) type for the Cahn-Hilliard equation in any space dimension. Specifically, we prove that our scheme is unconditionally energy stable and unconditionally uniquely solvable. Likewise, we devise and analyze a CS, mixed DG-FE scheme for our diff use interface cancer model. This scheme is energy stable for any (positive) time step size and for any (positive) space step size that is sufficiently small.

Files over 3MB may be slow to open. For best results, right-click and select "save as..."

Share

COinS