# American Institute of Mathematical Sciences

September  2018, 23(7): 2775-2801. doi: 10.3934/dcdsb.2018163

## Domain-growth-induced patterning for reaction-diffusion systems with linear cross-diffusion

 1 School of Mathematical and Physical Sciences, Department of Mathematics, University of Sussex, Pevensey Ⅲ, 5C15, Falmer, Brigton, BN1 9QH, England, UK 2 Escola Superior de Tecnologia do Barreiro/IPS, Rua Américo da Silva Marinho-Lavradio, 2839-001 Barreiro, Portugal

Received  March 2017 Revised  March 2018 Published  June 2018

In this article we present, for the first time, domain-growth induced pattern formation for reaction-diffusion systems with linear cross-diffusion on evolving domains and surfaces. Our major contribution is that by selecting parameter values from spaces induced by domain and surface evolution, patterns emerge only when domain growth is present. Such patterns do not exist in the absence of domain and surface evolution. In order to compute these domain-induced parameter spaces, linear stability theory is employed to establish the necessary conditions for domain-growth induced cross-diffusion-driven instability for reaction-diffusion systems with linear cross-diffusion. Model reaction-kinetic parameter values are then identified from parameter spaces induced by domain-growth only; these exist outside the classical standard Turing space on stationary domains and surfaces. To exhibit these patterns we employ the finite element method for solving reaction-diffusion systems with cross-diffusion on continuously evolving domains and surfaces.

Citation: Anotida Madzvamuse, Raquel Barreira. Domain-growth-induced patterning for reaction-diffusion systems with linear cross-diffusion. Discrete & Continuous Dynamical Systems - B, 2018, 23 (7) : 2775-2801. doi: 10.3934/dcdsb.2018163
##### References:

show all references

##### References:
Phase-diagrams corresponding to the non-autonomous system of ordinary differential equations (3.1)-(3.2) with exponential growth. A stable limit cycle or spiral point exists depending on the choice of the parameter values $a$ and $b$. (a) $a = 0.1$, $b = 0.75$, (b) $a = 0.15$, $b = 0.6$, (c) $a = 0.15$, $b = 0.5$, (d) $a = 0.1$, $b = 0.1$, (e) $a = 0.1$, $b = 0.75$, and (f) $a = 0.12$, $b = 0.5$. Other parameter values are fixed as follows: $r = 0.01$, $\gamma = 200$, $m = 2$, $\kappa = 2$ and $A = 1$
(a) Snap-shots of the evolving parameter space for exponential growth rate with diffusion coefficient $d = 10,$ linear cross-diffusion coefficients $d_u = d_v = 0$, with $\gamma = 200$ and growth rate $r = 0.01$. We select the model kinetic parameter values $a = 0.1$ and $b = 0.75$ from the green parameter space in (a). (b)- (e) Finite element numerical simulations exhibiting the formation of spatial structure corresponding to the chemical specie $u$ during growth development. In the absence of domain growth, these patterns are non-existent
(a) Snap-shots of the evolving parameter space for exponential growth rate with diffusion coefficient $d = 10,$ linear cross-diffusion coefficients $d_u = 1, d_v = 0$, with $\gamma = 200$ and growth rate $r = 0.01$. We select the model kinetic parameter values $a = 0.1$ and $b = 0.75$ from the green parameter space in (a). (b)- (h) Finite element numerical simulations exhibiting the formation of spatial structure corresponding to the chemical specie $u$ during growth development. Note that when the domain is not sufficiently large enough, no patterns are observed
(a) Snap-shots of the evolving parameter space for exponential growth rate with diffusion coefficient $d = 10,$ linear cross-diffusion coefficients $d_u = 0, d_v = 1$, with $\gamma = 200$ and growth rate $r = 0.01$. We select the model kinetic parameter values $a = 0.15$ and $b = 0.6$ from the green parameter space in (a). (b)- (e) Finite element numerical simulations exhibiting the formation of spatial structure corresponding to the chemical specie $u$ during growth development
(a) Snap-shots of the evolving parameter space for exponential growth rate with diffusion coefficient $d = 10,$ linear cross-diffusion coefficients $d_u = 1, d_v = 1$, with $\gamma = 200$ and growth rate $r = 0.01$. We select the model kinetic parameter values $a = 0.15$ and $b = 0.5$ from the green parameter space in (a). (b)- (f) Finite element numerical simulations exhibiting the formation of spatial structure corresponding to the chemical specie $u$ during growth development
(a) Snap-shots of the evolving parameter space for exponential growth rate with diffusion coefficient $d = 1,$ linear cross-diffusion coefficients $d_u = 1, d_v = 0.5$, with $\gamma = 200$ and growth rate $r = 0.01$. We select the model kinetic parameter values $a = 0.1$ and $b = 0.1$ from the green parameter space in (a). (b)- (e) Finite element numerical simulations exhibiting the formation of spatial structure corresponding to the chemical specie $u$ during growth development
(a) Snap-shots of the evolving parameter space for exponential growth rate with diffusion coefficient $d = 1,$ linear cross-diffusion coefficients $d_u = 1, d_v = 0.5$, with $\gamma = 200$ and growth rate $r = 0.01$. We select the model kinetic parameter values $a = 0.1$ and $b = 0.75$ from the green parameter space in (a). (b) - (f) Finite element numerical simulations exhibiting the formation of spatial structure corresponding to the chemical specie $u$ during growth development
(a) Snap-shots of the evolving parameter space for exponential growth rate with diffusion coefficient $d = 1,$ linear cross-diffusion coefficients $d_u = -0.8, d_v = 1$, with $\gamma = 200$ and growth rate $r = 0.01$. We select the model kinetic parameter values $a = 0.15$ and $b = 0.25$ from the green parameter space in (a). (b) - (f) Finite element numerical simulations exhibiting the formation of spatial structure corresponding to the chemical specie $u$ during growth development
(a) Snap-shots of the evolving parameter space for exponential growth rate with diffusion coefficient $d = 1,$ linear cross-diffusion coefficients $d_u = 1, d_v = 0.8$, with $\gamma = 200$ and growth rate $r = 0.01$. We select the model kinetic parameter values $a = 0.12$ and $b = 0.5$ from parameter space in $(a)$. (b) - (g) Finite element numerical simulations exhibiting the formation of spatial structure corresponding to the chemical specie $u$ during growth development
(a)-(d) Snap-shots of the finite element numerical simulations on the evolving unit sphere under exponential growth rate with diffusion coefficient $d = 1,$ linear cross-diffusion coefficients $d_u = -0.8, d_v = 1$, with $\gamma = 200$ and growth rate $r = 0.01$. We select the model kinetic parameter values $a = 0.15$ and $b = 0.2$ from the domain-induced parameter space shown in Figure 11(a)
(a)-(d) Snap-shots of the finite element numerical simulations on the evolving saddle like surface under exponential growth rate with diffusion coefficient $d = 1,$ linear cross-diffusion coefficients $d_u = 1, d_v = 0.5$, with $\gamma = 200$ and growth rate $r = 0.01$. We select the model kinetic parameter values $a = 0.1$ and $b = 0.1$ from the domain-induced parameter space shown in Figure 9(a)
Snap-shots of continuously evolving parameter spaces for an exponential evolution of the domain with diffusion coefficient $d = 10$ with varying linear cross-diffusion coefficients $d_u$ and $d_v$ and the growth rate $r$. Row-wise: (a)-(c) $d_u = d_v = 0$, (d)-(f) $d_u = 0$ and $d_v = 1$, (g)-(i) $d_u = 1$ and $d_v = 0$, (j)-(l) $d_u = d_v = 1$. Column-wise: First column: $r = 0.01$, second column, $r = 0.04$ and third column: $r = 0.08$. All plots are exhibited at times $t = 0$ (no growth), $t = 0.1$ and $t = 0.3$
Snap-shots of continuously evolving parameter spaces for an exponential evolution of the domain with diffusion coefficient $d = 1$ with variable linear cross-diffusion coefficients and growth rates: Row-wise: (a)-(c) $d_u = 1$, $d_v = 0.8$, (d)-(f) $d_u = 0.8$, $d_v = 1$, and (g)-(i) $d_u = -0.8$, $d_v = 1$. Column-wise: First column: $r = 0.01$, second column, $r = 0.04$ and third column: $r = 0.08$. All plots are exhibited at times $t = 0$ (no growth), $t = 0.1$ and $t = 0.3$. Such spaces do not exist if linear cross-diffusion is absent (with or without domain growth)
Super imposed parameter spaces with no growth, linear, logistic and exponential growth profiles of the domain with diffusion coefficient $d = 1$ and linear cross-diffusion coefficients $d_u = 1$ and $d_v = 0.8$. The growth rate $r$ is varied row-wise accordingly: (a)-(b) $r = 0.01$, (c)-(d) $r = 0.04$, and (e)-(f) $r = 0.08$. Snap-shots of the continuously evolving parameter spaces at times $t = 0$ and $t = 0.1$ (first column) and $t = 0$ and $t = 0.3$ (second column), respectively. Without linear cross-diffusion, these spaces do not exist
Table illustrating the function $h(t)$ for linear, exponential and logistic growth functions on evolving domains. Similar functions can be obtained for evolving surfaces. $\kappa$ is the carrying capacity (final domain size) corresponding to the logistic growth function. In the above $h (t)$ as defined by (2.6) (or (2.7)) on planar domains (or on surfaces)
 Type of growth Growth Function $\rho (t)$ $h(t)= \nabla \cdot \mathit{\boldsymbol{v}}$ $q(t)=e^{-\int_{t_0}^t h(\tau) d\tau }$ Linear $\rho(t)=rt+1$ $h(t)=\frac{m\, r}{rt+1}$ $q(t)= \left(\frac{1}{rt+1}\right)^m$ Exponential $\rho(t)=e^{rt}$ $h(t)=m\, r$ $q(t)=e^{-mrt}$ Logistic $\rho (t) = \frac{\kappa Ae^{\kappa rt}}{1+Ae^{\kappa rt}}, \, A = \frac{1}{\kappa-1}$ $h(t)= m\, r\, \Big(\kappa-\rho (t)\Big)$ $q(t) = \left(\frac{e^{-\kappa r t}+A}{1+A}\right)^m$
 Type of growth Growth Function $\rho (t)$ $h(t)= \nabla \cdot \mathit{\boldsymbol{v}}$ $q(t)=e^{-\int_{t_0}^t h(\tau) d\tau }$ Linear $\rho(t)=rt+1$ $h(t)=\frac{m\, r}{rt+1}$ $q(t)= \left(\frac{1}{rt+1}\right)^m$ Exponential $\rho(t)=e^{rt}$ $h(t)=m\, r$ $q(t)=e^{-mrt}$ Logistic $\rho (t) = \frac{\kappa Ae^{\kappa rt}}{1+Ae^{\kappa rt}}, \, A = \frac{1}{\kappa-1}$ $h(t)= m\, r\, \Big(\kappa-\rho (t)\Big)$ $q(t) = \left(\frac{e^{-\kappa r t}+A}{1+A}\right)^m$
 [1] Anotida Madzvamuse, Hussaini Ndakwo, Raquel Barreira. Stability analysis of reaction-diffusion models on evolving domains: The effects of cross-diffusion. Discrete & Continuous Dynamical Systems - A, 2016, 36 (4) : 2133-2170. doi: 10.3934/dcds.2016.36.2133 [2] Yuan Lou, Wei-Ming Ni, Shoji Yotsutani. Pattern formation in a cross-diffusion system. Discrete & Continuous Dynamical Systems - A, 2015, 35 (4) : 1589-1607. doi: 10.3934/dcds.2015.35.1589 [3] Hideki Murakawa. A relation between cross-diffusion and reaction-diffusion. Discrete & Continuous Dynamical Systems - S, 2012, 5 (1) : 147-158. doi: 10.3934/dcdss.2012.5.147 [4] Georg Hetzer, Anotida Madzvamuse, Wenxian Shen. Characterization of turing diffusion-driven instability on evolving domains. Discrete & Continuous Dynamical Systems - A, 2012, 32 (11) : 3975-4000. doi: 10.3934/dcds.2012.32.3975 [5] Joseph G. Yan, Dong-Ming Hwang. Pattern formation in reaction-diffusion systems with $D_2$-symmetric kinetics. Discrete & Continuous Dynamical Systems - A, 1996, 2 (2) : 255-270. doi: 10.3934/dcds.1996.2.255 [6] Yi Li, Chunshan Zhao. Global existence of solutions to a cross-diffusion system in higher dimensional domains. Discrete & Continuous Dynamical Systems - A, 2005, 12 (2) : 185-192. doi: 10.3934/dcds.2005.12.185 [7] Martino Prizzi. A remark on reaction-diffusion equations in unbounded domains. Discrete & Continuous Dynamical Systems - A, 2003, 9 (2) : 281-286. doi: 10.3934/dcds.2003.9.281 [8] Hongyan Zhang, Siyu Liu, Yue Zhang. Dynamics and spatiotemporal pattern formations of a homogeneous reaction-diffusion Thomas model. Discrete & Continuous Dynamical Systems - S, 2017, 10 (5) : 1149-1164. doi: 10.3934/dcdss.2017062 [9] Robert Stephen Cantrell, Xinru Cao, King-Yeung Lam, Tian Xiang. A PDE model of intraguild predation with cross-diffusion. Discrete & Continuous Dynamical Systems - B, 2017, 22 (10) : 3653-3661. doi: 10.3934/dcdsb.2017145 [10] Yuan Lou, Wei-Ming Ni, Yaping Wu. On the global existence of a cross-diffusion system. Discrete & Continuous Dynamical Systems - A, 1998, 4 (2) : 193-203. doi: 10.3934/dcds.1998.4.193 [11] Lianzhang Bao, Wenjie Gao. Finite traveling wave solutions in a degenerate cross-diffusion model for bacterial colony with volume filling. Discrete & Continuous Dynamical Systems - B, 2017, 22 (7) : 2813-2829. doi: 10.3934/dcdsb.2017152 [12] Peng Feng, Zhengfang Zhou. Finite traveling wave solutions in a degenerate cross-diffusion model for bacterial colony. Communications on Pure & Applied Analysis, 2007, 6 (4) : 1145-1165. doi: 10.3934/cpaa.2007.6.1145 [13] Lukas F. Lang, Otmar Scherzer. Optical flow on evolving sphere-like surfaces. Inverse Problems & Imaging, 2017, 11 (2) : 305-338. doi: 10.3934/ipi.2017015 [14] Peter E. Kloeden, Meihua Yang. Forward attracting sets of reaction-diffusion equations on variable domains. Discrete & Continuous Dynamical Systems - B, 2019, 24 (3) : 1259-1271. doi: 10.3934/dcdsb.2019015 [15] Yuan Lou, Wei-Ming Ni, Shoji Yotsutani. On a limiting system in the Lotka--Volterra competition with cross-diffusion. Discrete & Continuous Dynamical Systems - A, 2004, 10 (1&2) : 435-458. doi: 10.3934/dcds.2004.10.435 [16] Kousuke Kuto, Yoshio Yamada. Coexistence states for a prey-predator model with cross-diffusion. Conference Publications, 2005, 2005 (Special) : 536-545. doi: 10.3934/proc.2005.2005.536 [17] Kousuke Kuto, Yoshio Yamada. On limit systems for some population models with cross-diffusion. Discrete & Continuous Dynamical Systems - B, 2012, 17 (8) : 2745-2769. doi: 10.3934/dcdsb.2012.17.2745 [18] Daniel Ryan, Robert Stephen Cantrell. Avoidance behavior in intraguild predation communities: A cross-diffusion model. Discrete & Continuous Dynamical Systems - A, 2015, 35 (4) : 1641-1663. doi: 10.3934/dcds.2015.35.1641 [19] F. Berezovskaya, Erika Camacho, Stephen Wirkus, Georgy Karev. "Traveling wave'' solutions of Fitzhugh model with cross-diffusion. Mathematical Biosciences & Engineering, 2008, 5 (2) : 239-260. doi: 10.3934/mbe.2008.5.239 [20] Ciprian G. Gal, Mahamadi Warma. Reaction-diffusion equations with fractional diffusion on non-smooth domains with various boundary conditions. Discrete & Continuous Dynamical Systems - A, 2016, 36 (3) : 1279-1319. doi: 10.3934/dcds.2016.36.1279

2018 Impact Factor: 1.008