Publication:
Aplicaciones del método LDG a sistemas químicos y biológicos con estructuras de Turing

Thumbnail Image
Authors
Henríquez-Rivera, Iván Y.
Embargoed Until
Advisor
Castillo, Paul E.
College
College of Arts and Sciences - Sciences
Department
Department of Mathematics
Degree Level
M.S.
Publisher
Date
2013-08
Abstract
En su trabajo pionero sobre morfogénesis [34] en 1952, A. Turing propuso un modelo teórico de dos especies químicas en el cual la evolución de un estado estacionario heterogéneo (estructura de Turing) desde un estado estacionario homogéneo es debido a la difusión. Desde la primera evidencia física de la existencia de tales estructuras de Turing el modelamiento teórico y diseño experimental de sistemas biológicos y químicos con estructuras de Turing han sido testigos de un progreso considerable. Típicamente tales modelos son matemáticamente formulados como un sistema de ecuaciones no lineales de reacción difusión ∂/∂tc(t,x) = (D ⊗ ∆)c(t,x) + r(c(t,x)), (2) donde c es un vector en R p que representa la concentración de p especies químicas, Dp es la matriz de coeficientes de difusión la cual generalmente se asume ser diagonal y constante; y rp es el operador de reacción que acopla el sistema. La naturaleza no lineal de las reacciones químicas hace difícil (de no ser imposible) expresar la solución del sistema con una formula exacta. Esto ha motivado el interés en el desarrollo de técnicas computacionales para la simulación numérica de tales sistemas de ecuaciones diferenciales en derivadas parciales. En este trabajo se realiza una comparación de dos versiones del método Local discontinuous Galerkin (LDG) para la discretización en espacio: la clásica y la de disipación mínima; la integración en tiempo es llevada a cabo mediante métodos de descomposición de operador o un método IMEX. Usando un método de Galerkin discontinuo para la desratización en espacio el sistema global no lineal es desacoplado en cada paso de tiempo en un conjunto de pequeños sistemas no lineales que pueden ser resueltos en paralelo. Un marco computacional fue implementado para discretizaciones espaciales de alto orden en 2D para sistemas no lineales de reacción difusión con un número arbitrario de especies. La comparación es llevada a cabo usando un sistema lineal de reacción difusión, el modelo de Schnakenberg y un modelo de pigmentación de pieles de animales.

In the pioneering work on morphogenesis [34] in 1952, A. Turing proposed a theoretical model consisting of two chemical species in which the evolution of a spatially stable heterogeneous steady state (Turing structures) from a homogeneous stationary state is driven by diffusion. Since the first physical evidence of the existence of such structures theoretical modeling and experimental design of chemical and biological systems revealing Turing structures has witnessed considerable progress. Typically such models are mathematically formulated as a system of non-linear reaction diffusion equations ∂/∂tc(t, x) = (D ⊗ ∆)c(t, x) + r(c(t, x)), (2) where c is a vector in Rp representing the concentration of p chemical species, Dp is a matrix of diffusion coefficients which is generally assumed diagonal and constant; and rp couples the system with p chemical reactions. The non-linearity nature of chemical kinetics makes it difficult (if not impossible) to express the solution of the system in closed form. This motivated an increase of interest in the development of computational techniques for the numerical simulation of such systems of partial differential equations. In this work a numerical comparison of two versions of the Local Discontinuous Galerkin (LDG) method for the spatial discretization: the classical and the minimal dissipation LDG; coupled with an operator splitting or an IMEX time stepping technique is performed. By using a discontinuous Galerkin discretization the global non-linear system of equations is decoupled at each time step into a set of small non-linear systems which can be solved in parallel. A computational framework was implemented for high order discretizations in 2D for non-linear reaction diffusion systems with an arbitrary number of species. The comparison is carried out using a linear reaction diffusion, the classical Schnakenberg model problem and a coat pigmentation model.
Keywords
Turing structures,
Theoretical model,
Diffusion equations
Cite
Henríquez-Rivera, I. Y. (2013). Aplicaciones del método LDG a sistemas químicos y biológicos con estructuras de Turing [Thesis]. Retrieved from https://hdl.handle.net/20.500.11801/225