[论文解读] Effects of cross-diffusion on Turing patterns in a reaction-diffusion Schnakenberg model

G. Gambino, Simona Lupo+1 more|arXiv (Cornell University)|Jan 20, 2015

Nonlinear Dynamics and Pattern Formation参考文献 38被引用 11

一句话总结

本文研究了线性交叉扩散对双物种反应-扩散Schnakenberg模型中图灵斑图形成的影响。通过弱非线性多尺度分析，识别出在交叉扩散促进次临界分岔的参数区域——产生稳定且有限振幅的斑图，尤其当激活剂的交叉扩散（d_v）增大或抑制剂的交叉扩散（d_u）减小时更为显著。

ABSTRACT

In this paper the Turing pattern formation mechanism of a two component reaction-diffusion system modeling the Schnakenberg chemical reaction coupled to linear cross-diffusion terms is studied. The linear cross-diffusion terms favors the destabilization of the constant steady state and the mechanism of pattern formation with respect to the standard linear diffusion case, as shown in Madzvamuse et al. (J. Math. Biol. 2014). Since the subcritical Turing bifurcations of reaction-diffusion systems lead to spontaneous onset of robust, finite-amplitude localized patterns, here a detailed investigation of the Turing pattern forming region is performed to show how the diffusion coefficients for both species (the activator and the inhibitor) influence the occurrence of supercritical or subcritical bifurcations. The weakly nonlinear (WNL) multiple scales analysis is employed to derive the equations for the amplitude of the Turing patterns and to distinguish the supercritical and the subcritical pattern region, both in 1D and 2D domains. Numerical simulations are employed to confirm the WNL theoretical predictions through which a classification of the patterns (squares, rhombi, rectangle and hexagons) is obtained. In particular, due to the hysteretic nature of the subcritical bifurcation, we observe the phenomenon of pattern transition from rolls to hexagons, in agreement with the bifurcation diagram.

研究动机与目标

• 理解线性交叉扩散如何改变反应-扩散系统中均匀稳态的稳定性。
• 确定在存在交叉扩散时，次临界与超临界图灵分岔发生的条件。
• 对一维和二维区域中产生的空间斑图（如条纹、六边形、方形）进行分类。
• 研究交叉扩散系数d_u和d_v在稳定或不稳定化系统以形成有限振幅斑图中的作用。
• 通过数值模拟验证理论预测，并分析斑图转变（如条纹到六边形）过程。

提出的方法

• 采用弱非线性（WNL）多尺度分析，推导1D和2D区域中图灵斑图的振幅方程。
• 推导Stuart-Landau型振幅方程，根据非线性系数的符号区分超临界与次临界分岔。
• 利用Fredholm择一性条件消除摄动展开中的非共振项，确保解有界。
• 在时间和空间上进行双尺度渐近展开，引入慢时间尺度T_2和空间尺度，以捕捉斑图演化过程。
• 采用Schnakenberg动力学形式化控制方程：f(u,v) = a - u + u²v，g(u,v) = b - u²v，其中两个物种均包含线性交叉扩散项。
• 开展数值模拟以验证理论预测，并观察斑图转变，包括条纹到六边形的转变。

实验结果

研究问题

• RQ1 线性交叉扩散如何影响反应-扩散系统中的分岔类型（超临界与次临界）？
• RQ2 激活剂的交叉扩散系数d_v对次临界斑图区域的大小和稳定性有何影响？
• RQ3 抑制剂的交叉扩散系数d_u如何调节次临界不稳定性发生的时机？
• RQ4 在二维区域中，交叉扩散下会涌现出何种空间斑图（如六边形、方形、矩形）？
• RQ5 系统如何表现出斑图转变（如从条纹到六边形），其驱动力是什么？

主要发现

• 激活剂中的交叉扩散（d_v）增强了次临界分岔区域，促进产生稳定且有限振幅的斑图。
• 抑制剂中的交叉扩散（d_u）缩小了次临界区域，使通过此途径形成斑图的可能性降低。
• 弱非线性分析成功通过多尺度展开导出的Stuart-Landau方程预测了斑图的振幅与形状。
• 在二维区域中，系统表现出丰富的斑图类型，包括条纹、六边形、方形、菱形以及混合模态结构。
• 当最不稳定的模态具有双重重数且发生共振时，分岔总是次临界的，导致具有滞后特性的转变。
• 数值模拟证实了理论预测：当分岔参数变化时，系统确实发生从条纹到六边形的转变，与分岔图一致。