Numerical solution of biological reaction diffusion models on fixed domains by the finite element method


  • Juan Carlos Vanegas Universidad Nacional de Colombia
  • Libardo Andrés González Universidad Nacional de Colombia
  • Diego Alexander Garzón Universidad Nacional de Colombia


Reaction diffusion models, pattern formation, finite element method, biomathematics


Several biological phenomena have been described using mathematical models based on reaction diffusion equations. The solution of this type of equations gives rise to formation of spatial-temporal patterns, in agreement with the biological reality of the simulated phenomenon. This article describes the numerical implementation of a set of three well-known reaction diffusion models: the morphogenesis Schnakenberg model, and the GiererMeinhardt and Thomas reaction kinetics models. The aim is to analyze the set of parameters associated with the spatial-temporal pattern formation. The numerical implementation was performed using the finite element method in one dimensional and two dimensional domains. It was concluded that spatialtemporal pattern formation in reaction diffusion models depends on the constant parameters of the model, the initial conditions and the implementation technique. The analysis of these dependences is useful in the formulation and validation of new mathematical models describing biological phenomena.

= 73 veces | PDF (ESPAÑOL (ESPAÑA))
= 61 veces|


Download data is not yet available.

Author Biographies

Juan Carlos Vanegas, Universidad Nacional de Colombia

Grupo de Modelamiento Matemático y Métodos Numéricos GNUM. Facultad de Ingeniería

Libardo Andrés González, Universidad Nacional de Colombia

Grupo de Modelamiento Matemático y Métodos Numéricos GNUM. Facultad de Ingeniería

Diego Alexander Garzón, Universidad Nacional de Colombia

Grupo de Modelamiento Matemático y Métodos Numéricos GNUM. Facultad de Ingeniería


P. K. Maini. “Mathematical models in morphogenesis”. Mathematics Inspired by Biology. Ed. Springer Berlin- Heidelberg. 1999. pp. 151-189.

D. A. Garzón. Simulación de procesos de reacción difusión: aplicación a la morfogénesis del tejido óseo. Tesis de Doctorado. Zaragoza. España: Centro Politécnico Superior de la Universidad de Zaragoza. 2007. pp. 38-70.

D. Benson, P. Maini, J. A. Sherratt. “Unravelling the Turing bifurcation using spatially varying diffusion coefficient”. J. Math. Biol. Vol. 37. 1998. pp. 381- 417.

P. K. Maini, K. Painter, H. Chau. “Spatial pattern formation in chemical and biological systems”. J. Chem. Soc. Faraday Trans. Vol. 93. 1997. pp. 3601- 3610.

A. Turing. “The chemical basis of morphogenesis”. Philos. Trans. Roy. Soc. Vol. 237. 1952. pp. 37-72.

M. Dolnik, A. M. Zhabotinsky, A. B. Rovinsky, I. R. Epstein. “Spatio-temporal patterns in a reaction diffusion system with wave instability”. Chemical Engineering Science. Vol. 55. 2000. pp. 223-231.

J. D. Murray. Mathematical Biology II. Spatial models and biomedical applications. 3a. ed. Ed. Springer-Verlag. New York. 2003. pp. 75-97.

K. Page, P. K. Maini, N. Monk. “Patter formation in spatially heterogeneous Turing reaction diffusion models”. Physics D. Vol. 181. 2003. pp. 80-101

A. Madzvamuse, A. J. Wathen, P. K. Maini. “A moving grid finite element method applied to a model biological patter generador”. Journal of Computational Physics. Vol. 190. 2003. pp. 478-500.

T. Miura, K. Shiota. “TGFβ-2 acts as an activator molecule in reaction-diffusion model and is involved in cell sorting phenomenon in mouse limb micromass culture”. Dev. Dyn. Vol. 217. 2000. pp. 241-249.

S. Kondo, R. Asai. “A reaction diffusion wave on the skin of the marine angelfish Pomacanthus”. Nature. Vol. 376. 1995. pp. 765-768.

H. Jung, P. Francis-West, R. Widelitz, T. Jiang, S. Ting- Berreth, C. Tickle, L. Wolpert, C. Chuong. “Local inhibitory action of BMPs and their relationships with activators in feather formation: implications for periodic patterning”. Dev. Biol. Vol. 196. 1998. pp. 11-23.

T. Miura, K. Morriss, P. K. Maini. “Mixed-mode pattern in doublefoot mutant mouse limb: Turing reaction diffusion model on a growing domain during limb development”. J. Theor Biol. Vol. 240. 2006. pp. 562-573.

J. D. Murray. Mathematical Biology I. An introduction.3a ed. Ed. Springer-Verlag. New York. 2002. pp. 405-509.

K. Painter. Chemotaxis as a mechanism for morphogenesis. PhD Thesis. United Kingdom. Oxford University. 1997. pp. 3-25.

J. C. Vanegas A., N. S. Landinez, D. A. Garzón-Alvarado. “Análisis de la inestabilidad de Turing en modelos biológicos”. Revista DYNA, Universidad Nacional de Colombia - Medellín. Aprobado para publicación. 2009.

J. D. Murray, G. F. Ester. “Cell traction models for generation pattern and form in morphogenesis”. J. Math. Biology. Vol. 19. 1984. pp. 265-279.

A. Madzvamuse. A numerical approach to the study of spatial pattern formation. PhD Thesis. Oxford, UK: Computing Laboratory. University of Oxford. 2000. pp. 10-40.

E. Crampin. Reaction diffusion patterns on growing domains. PhD Thesis. Oxford, UK: Magdalen College. University of Oxford. 2000. pp.1-34.

E. Oñate, J. Miquel, F. Zárate. “Stabilized solution of the multidimensional advection–diffusion–absorption equation using linear finite elements”. Computers and Fluids. Vol. 36. 2007. pp. 92–112.

A. Madzvamuse. “Time-stepping schemes for moving grid finite elements applied to reaction diffusion systems on fixed and growing domains”. J. Comput Phys. Vol. 214. 2006. pp. 239-263.



How to Cite

Vanegas, J. C. ., González, L. A., & Garzón, D. A. (2013). Numerical solution of biological reaction diffusion models on fixed domains by the finite element method. Revista Facultad De Ingeniería Universidad De Antioquia, (48), 65–75. Retrieved from