Technische Universität Darmstadt Fachbereich Mathematik   Alf Gerisch

in: M. Deville and R. Owens (Eds.). Proceedings of the 16th IMACS World Congress 2000, Lausanne, Switzerland. ISBN 3-9522075-1-9.

Abstract: Angiogenesis -- the process by which new blood vessels grow into a tissue from surrounding parent vessels - is an important process in many areas of medicine. Here we consider the numerical simulation of a PDE model of tumor-induced angiogenesis. It contains convection (migration), diffusion and reaction terms. Despite the restriction to one specific model, the observations should also be relevant for the solution of similar problems.
Our general approach is the method of lines and we use a positivity preserving spatial discretization resulting in a large and in general stiff ODE system. For the solution of this system we consider splitting methods (approximate matrix factorization, Strang-type and source splitting) and Krylov-W-methods. The aim is to reduce the complexity of the implicit relations in the solution process.
Advantages and disadvantages of the different approaches are discussed. We compare the methods with respect to efficiency and accuracy of the solution. A Strang-type splitting method combined with approximate matrix factorization is found to be most efficient in the low to modest accuracy range and this range is of interest for the model.

MSC(1991):

65M06 Finite difference methods

65M20 Method of lines

65Y20 Complexity and performance of numerical algorithms, see also 68Q25

65Y05 Parallel computation

Keywords: tumor angiogenesis, chemotaxis, time-dependent convection-reaction-diffusion system, numerical methods, stiff ODE, splitting, approximate factorization, linearly-implicit, parallel method

Notes:

Last Modification: 06.09.2010 17:36 Author: Alf Gerisch