Analysis and numerical simulation of tumor growth models

  1. Acosta Soba, Daniel
Supervised by:
  1. José Rafael Rodríguez Galván Director
  2. Jin Wang Co-director
  3. Francisco Manuel Guillén González Co-director

Defence university: Universidad de Cádiz

Fecha de defensa: 08 March 2024

  1. Kristoffer van der Zee Chair
  2. Concepción García Vázquez Secretary
  3. Christopher Cox Committee member

Type: Thesis

Teseo: 834736 DIALNET lock_openTESEO editor


In this dissertation we focus on the numerical analysis of tumor growth models. Due to the difficulty of developing physically meaningful approximations of such models, we divide the main problem into more simple pieces of work that are addressed in the different chapters. First, in Chapter 2 we present a new upwind discontinuous Galerkin (DG) scheme for the convective Cahn-Hilliard model with degenerate mobility which preserves the pointwise bounds and prevents non-physical spurious oscillations. These ideas are based on a well-suited piecewise constant approximation of convection equations. The proposed numerical scheme is contrasted with other approaches in several numerical experiments. Afterwards, in Chapter 3, we extend the previous ideas to a mass-conservative, positive and energy-dissipative approximation of the Keller Segel model for chemotaxis. Then we carry out several numerical tests in regimes of chemotactic collapse. These ideas are used later in Chapter 4 to develop a well-suited approximation of two different models related to chemotaxis: a generalization of the classical Keller-Segel model and a model of the neuroblast migration process to the olfactory bulb in rodents brains. Now we propose and study a phase-field tumor growth model in Chapter 5. Then, we develop an upwind DG scheme preserving the mass conservation, pointwise bounds and energy stability of the continuous model and we show both the good properties of the approximation and the qualitative behavior of the model in several numerical tests. Next, in Chapter 6, we present two new coupled and decoupled approximations of a Cahn Hilliard Navier Stokes model with variable densities and degenerate mobility that preserve the physical properties of the model. Both approaches are compared in different computational tests including benchmark problems. Consequently, we propose, in Chapter 7, an extension of the previous tumor model including the effects of the surrounding fluid by means of a Cahn Hilliard Darcy model for which obtaining a physically meaningful approximation seems rather plausible using the previous ideas. Finally, this and other future lines of research are described, along with the conclusions and the scientific production of the dissertation, in Chapter 8.