# American Institute of Mathematical Sciences

## Turing instability and pattern formations for reaction-diffusion systems on 2D bounded domain

 School of Mathematics, Harbin Institute of Technology, Harbin 150001, China

* Corresponding author: Weihua Jiang

Received  December 2020 Revised  January 2021 Published  March 2021

Fund Project: The work was supported in part by the National Natural Science Foundation of China (No. 11871176, 11671110) and the Fundamental Research Funds for the Central Universities

In this article, Turing instability and the formations of spatial patterns for a general two-component reaction-diffusion system defined on 2D bounded domain, are investigated. By analyzing characteristic equation at positive constant steady states and further selecting diffusion rate $d$ and diffusion ratio $\varepsilon$ as bifurcation parameters, sufficient and necessary conditions for the occurrence of Turing instability are established, which is called the first Turing bifurcation curve. Furthermore, parameter regions in which single-mode Turing patterns arise and multiple-mode (or superposition) Turing patterns coexist when bifurcations parameters are chosen, are described. Especially, the boundary of parameter region for the emergence of single-mode Turing patterns, consists of the first and the second Turing bifurcation curves which are given in explicit formulas. Finally, by taking diffusive Schnakenberg system as an example, parameter regions for the emergence of various kinds of spatially inhomogeneous patterns with different spatial frequencies and superposition Turing patterns, are estimated theoretically and shown numerically.

Citation: Weihua Jiang, Xun Cao, Chuncheng Wang. Turing instability and pattern formations for reaction-diffusion systems on 2D bounded domain. Discrete & Continuous Dynamical Systems - B, doi: 10.3934/dcdsb.2021085
##### References:

show all references

##### References:
The graphs of functions $\varepsilon = \varepsilon_1$ and $\varepsilon = \varepsilon_*(\mathbf{k}(i),d),d\ge d_i,\; i = 1,2,3\cdots$ in $d \rm{-}\varepsilon$ plane
The first Turing bifurcation curve
$\varepsilon = \varepsilon_{**}(d)$ is the second Turing bifurcation curve. The grey area enclosed by the first and the second Turing bifurcation curves, represents $\mathfrak{D}^1$. The region below the second Turing bifurcation curve $\varepsilon = \varepsilon_{**}(d)$ is $\mathfrak{D}^2$, and the blue area denoted by $D_{\mathbf{k}(3), \mathbf{k}(7)}$ is one component of $\mathfrak{D}^2$
">Figure 4.  Turing patterns for system (18) with different values of $(d,\varepsilon)$ given in Table 2
(a), (b): For $\mathbf{k}(i) = (4,0)$ and $\mathbf{k}(j) = (0,4)$, there exists superposition pattern of (18) when $(d,\varepsilon) = (0.09,0.025)$; (c): The graph of function $z(x,y) = 0.9+1.8\cos(4\cdot\pi x)\cos(0\cdot\frac{\pi^2}{3} y)+\cos(0\cdot\pi x)\cos(4\cdot\frac{\pi^2}{3} y)$
(a), (b): Superposition pattern of (18) for $\mathbf{k}(i) = (3,0)$ and $\mathbf{k}(j) = (0,3)$; (c): The graph of function $z(x,y) = 1+1.8\cos(3\cdot\pi x)\cos(0\cdot\frac{\pi^2}{3} (\frac{\pi}{3}-y))+\cos(0\cdot\pi x)\cos(3\cdot\frac{\pi^2}{3} (\frac{\pi}{3}-y))$
(a), (b): Superposition of three kinds of spatial patterns with spatial wave numbers $(2,0)$, $(0,2)$ and $(2,2)$ for system (18); (c): The graph of function $z(x,y) = 0.9+\cos(2\cdot\pi x)\cos(0\cdot\frac{\pi^2}{3} y)+5\cos(0\cdot\pi x)\cos(2\cdot\frac{\pi^2}{3} y)+2\cos(2\cdot\pi x)\cos(2\cdot\frac{\pi^2}{3} y)$
These values of $\mathbf{k}(i),\; \mu_i/\pi^2$ and $d_{i,i+1}$ when parameters are chosen as (21), $i = 1,2,\ldots,28$
 $i$ $\mathbf{k}(i)$ $\mu_i/\pi^2$ $d_{i,i+1}$ $i$ $\mathbf{k}(i)$ $\mu_i/\pi^2$ $d_{i,i+1}$ 1 (1, 0) 1 0.2834 15 (4, 0) 16 0.0179 2 (0, 1) 1.0966 0.2034 16 (4, 1) 17.0966 0.0171 3 (1, 1) 2.0966 0.1064 17 (0, 4) 17.5460 0.0164 4 (2, 0) 4 0.0709 18 (1, 4) 18.5460 0.0160 5 (0, 2) 4.3865 0.0629 19 (3, 3) 18.8696 0.0151 6 (2, 1) 5.0966 0.0566 20 (4, 2) 20.3865 0.0142 7 (1, 2) 5.3865 0.0449 21 (2, 4) 21.5460 0.0128 8 (2, 2) 8.3865 0.0342 22 (5, 0) 25 0.0117 9 (3, 0) 9 0.0315 23 (4, 3) 25.8696 0.0114 10 (0, 3) 9.8696 0.0297 24 (5, 1) 26.0966 0.0113 11 (3, 1) 10.0966 0.0283 25 (3, 4) 26.5460 0.0110 12 (1, 3) 10.8696 0.0247 26 (0, 5) 27.4156 0.0106 13 (3, 2) 13.3865 0.0218 27 (1, 5) 28.4156 0.0103 14 (2, 3) 13.8696 0.0199 28 (5, 2) 29.3865 0.0098
 $i$ $\mathbf{k}(i)$ $\mu_i/\pi^2$ $d_{i,i+1}$ $i$ $\mathbf{k}(i)$ $\mu_i/\pi^2$ $d_{i,i+1}$ 1 (1, 0) 1 0.2834 15 (4, 0) 16 0.0179 2 (0, 1) 1.0966 0.2034 16 (4, 1) 17.0966 0.0171 3 (1, 1) 2.0966 0.1064 17 (0, 4) 17.5460 0.0164 4 (2, 0) 4 0.0709 18 (1, 4) 18.5460 0.0160 5 (0, 2) 4.3865 0.0629 19 (3, 3) 18.8696 0.0151 6 (2, 1) 5.0966 0.0566 20 (4, 2) 20.3865 0.0142 7 (1, 2) 5.3865 0.0449 21 (2, 4) 21.5460 0.0128 8 (2, 2) 8.3865 0.0342 22 (5, 0) 25 0.0117 9 (3, 0) 9 0.0315 23 (4, 3) 25.8696 0.0114 10 (0, 3) 9.8696 0.0297 24 (5, 1) 26.0966 0.0113 11 (3, 1) 10.0966 0.0283 25 (3, 4) 26.5460 0.0110 12 (1, 3) 10.8696 0.0247 26 (0, 5) 27.4156 0.0106 13 (3, 2) 13.3865 0.0218 27 (1, 5) 28.4156 0.0103 14 (2, 3) 13.8696 0.0199 28 (5, 2) 29.3865 0.0098
Parameter values of $(d,\varepsilon)$ in $\mathfrak{D}^1$ satisfying that (18) has $\mathbf{k}(i)-$mode Turing patterns
 $i$ $\mathbf{k}(i)$ $\varepsilon$ $(d_i^-(\varepsilon), d_i^+(\varepsilon))$ $d$ $(d,\varepsilon)\in$ Figure 3 (1, 1) 0.09 (0.08870, 0.2925) 0.1300 $D_{\mathbf{k}(3), \mathbf{k}(5)}$ 4(a) 6 (2, 1) 0.09 (0.0365, 0.12034) 0.0660 $D_{\mathbf{k}(4), \mathbf{k}(9)}$ 4(b) 8 (2, 2) 0.09 (0.02218, 0.0731) 0.0434 $D_{\mathbf{k}(5), \mathbf{k}(14)}$ 4(c) 9 (3, 0) 0.09 (0.0207, 0.0681) 0.0340 $D_{\mathbf{k}(8), \mathbf{k}(17)}$ 4(d) 12 (1, 3) 0.08 (0.0171, 0.0564) 0.0270 $D_{\mathbf{k}(8), \mathbf{k}(25)}$ 4(e) 13 (3, 2) 0.09 (0.0139, 0.0458) 0.0260 $D_{\mathbf{k}(8), \mathbf{k}(21)}$ 4(f)
 $i$ $\mathbf{k}(i)$ $\varepsilon$ $(d_i^-(\varepsilon), d_i^+(\varepsilon))$ $d$ $(d,\varepsilon)\in$ Figure 3 (1, 1) 0.09 (0.08870, 0.2925) 0.1300 $D_{\mathbf{k}(3), \mathbf{k}(5)}$ 4(a) 6 (2, 1) 0.09 (0.0365, 0.12034) 0.0660 $D_{\mathbf{k}(4), \mathbf{k}(9)}$ 4(b) 8 (2, 2) 0.09 (0.02218, 0.0731) 0.0434 $D_{\mathbf{k}(5), \mathbf{k}(14)}$ 4(c) 9 (3, 0) 0.09 (0.0207, 0.0681) 0.0340 $D_{\mathbf{k}(8), \mathbf{k}(17)}$ 4(d) 12 (1, 3) 0.08 (0.0171, 0.0564) 0.0270 $D_{\mathbf{k}(8), \mathbf{k}(25)}$ 4(e) 13 (3, 2) 0.09 (0.0139, 0.0458) 0.0260 $D_{\mathbf{k}(8), \mathbf{k}(21)}$ 4(f)
Parameter values of $(d,\varepsilon)$ in $\mathfrak{D}^2$ satisfying that (18) has superposition patterns
 $(d,\varepsilon)$ $\mathbf{k}(i)$ $\mathbf{k}(j)$ $\mathbf{k}(l)$ Figure (0.025, 0.09) $\in D_{\mathbf{k}(8), \mathbf{k}(21)}$ (4, 0) (0, 4) --- 5 (0.044, 0.09) $\in D_{\mathbf{k}(5), \mathbf{k}(14)}$ (3, 0) (0, 3) --- 6 (0.094, 0.09) $\in D_{\mathbf{k}(3), \mathbf{k}(7)}$ (2, 0) (0, 2) (2, 2) 7
 $(d,\varepsilon)$ $\mathbf{k}(i)$ $\mathbf{k}(j)$ $\mathbf{k}(l)$ Figure (0.025, 0.09) $\in D_{\mathbf{k}(8), \mathbf{k}(21)}$ (4, 0) (0, 4) --- 5 (0.044, 0.09) $\in D_{\mathbf{k}(5), \mathbf{k}(14)}$ (3, 0) (0, 3) --- 6 (0.094, 0.09) $\in D_{\mathbf{k}(3), \mathbf{k}(7)}$ (2, 0) (0, 2) (2, 2) 7
 [1] Yuncherl Choi, Taeyoung Ha, Jongmin Han, Sewoong Kim, Doo Seok Lee. Turing instability and dynamic phase transition for the Brusselator model with multiple critical eigenvalues. Discrete & Continuous Dynamical Systems, 2021  doi: 10.3934/dcds.2021035 [2] José Raúl Quintero, Juan Carlos Muñoz Grajales. On the existence and computation of periodic travelling waves for a 2D water wave model. Communications on Pure & Applied Analysis, 2018, 17 (2) : 557-578. doi: 10.3934/cpaa.2018030 [3] Liangliang Ma. Stability of hydrostatic equilibrium to the 2D fractional Boussinesq equations. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021068 [4] Peng Tong, Xiaogang Ma. Design of differentiated warranty coverage that considers usage rate and service option of consumers under 2D warranty policy. Journal of Industrial & Management Optimization, 2021, 17 (4) : 1577-1591. doi: 10.3934/jimo.2020035 [5] Xin-Guang Yang, Rong-Nian Wang, Xingjie Yan, Alain Miranville. Dynamics of the 2D Navier-Stokes equations with sublinear operators in Lipschitz-like domains. Discrete & Continuous Dynamical Systems, 2021, 41 (7) : 3343-3366. doi: 10.3934/dcds.2020408 [6] Qixiang Wen, Shenquan Liu, Bo Lu. Firing patterns and bifurcation analysis of neurons under electromagnetic induction. Electronic Research Archive, , () : -. doi: 10.3934/era.2021034 [7] Yuta Ishii, Kazuhiro Kurata. Existence of multi-peak solutions to the Schnakenberg model with heterogeneity on metric graphs. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2021035 [8] Izumi Takagi, Conghui Zhang. Existence and stability of patterns in a reaction-diffusion-ODE system with hysteresis in non-uniform media. Discrete & Continuous Dynamical Systems, 2021, 41 (7) : 3109-3140. doi: 10.3934/dcds.2020400 [9] Yu Jin, Xiao-Qiang Zhao. The spatial dynamics of a Zebra mussel model in river environments. Discrete & Continuous Dynamical Systems - B, 2021, 26 (4) : 1991-2010. doi: 10.3934/dcdsb.2020362 [10] Wenbin Yang, Yujing Gao, Xiaojuan Wang. Diffusion modeling of tumor-CD4$^+$-cytokine interactions with treatments: asymptotic behavior and stationary patterns. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021090 [11] Mingxin Wang, Qianying Zhang. Dynamics for the diffusive Leslie-Gower model with double free boundaries. Discrete & Continuous Dynamical Systems, 2018, 38 (5) : 2591-2607. doi: 10.3934/dcds.2018109 [12] Weiyi Zhang, Zuhan Liu, Ling Zhou. Dynamics of a nonlocal diffusive logistic model with free boundaries in time periodic environment. Discrete & Continuous Dynamical Systems - B, 2021, 26 (7) : 3767-3784. doi: 10.3934/dcdsb.2020256 [13] Teddy Pichard. A moment closure based on a projection on the boundary of the realizability domain: 1D case. Kinetic & Related Models, 2020, 13 (6) : 1243-1280. doi: 10.3934/krm.2020045 [14] Guodong Wang, Bijun Zuo. Energy equality for weak solutions to the 3D magnetohydrodynamic equations in a bounded domain. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021078 [15] Luigi C. Berselli, Jishan Fan. Logarithmic and improved regularity criteria for the 3D nematic liquid crystals models, Boussinesq system, and MHD equations in a bounded domain. Communications on Pure & Applied Analysis, 2015, 14 (2) : 637-655. doi: 10.3934/cpaa.2015.14.637 [16] Brian Ryals, Robert J. Sacker. Bifurcation in the almost periodic $2$D Ricker map. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021089 [17] Reza Mazrooei-Sebdani, Zahra Yousefi. The coupled 1:2 resonance in a symmetric case and parametric amplification model. Discrete & Continuous Dynamical Systems - B, 2021, 26 (7) : 3737-3765. doi: 10.3934/dcdsb.2020255 [18] Lin Yang, Yejuan Wang, Tomás Caraballo. Regularity of global attractors and exponential attractors for $2$D quasi-geostrophic equations with fractional dissipation. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021093 [19] Yangrong Li, Fengling Wang, Shuang Yang. Part-convergent cocycles and semi-convergent attractors of stochastic 2D-Ginzburg-Landau delay equations toward zero-memory. Discrete & Continuous Dynamical Systems - B, 2021, 26 (7) : 3643-3665. doi: 10.3934/dcdsb.2020250 [20] Shihu Li, Wei Liu, Yingchao Xie. Large deviations for stochastic 3D Leray-$\alpha$ model with fractional dissipation. Communications on Pure & Applied Analysis, 2019, 18 (5) : 2491-2509. doi: 10.3934/cpaa.2019113

2019 Impact Factor: 1.27