双层耦合非对称反应扩散系统中的振荡图灵斑图*

刘雅慧 董梦菲 刘富成2)† 田淼 王硕 范伟丽

1) (河北大学物理科学与技术学院, 保定 071002)

2) (河北大学生命科学与绿色发展研究院, 保定 071002)

1 引 言

自组织斑图是非线性动力学系统的一个典型特征, 普遍存在于自然界以及众多非线性系统中,例如物理、化学以及生物系统等[1-7].在过去的大半个世纪内, 人们对斑图形成机理的认识有了突飞猛进的发展.众所周知, 非线性系统中时空斑图的形成源于系统的失稳, 且总是伴随着相应的对称性破缺现象, 不同的对称性破缺(或者说分岔类型)对应着不同的时空斑图类型.其中最为常见的分岔类型有三种: 霍普夫分岔、图灵分岔以及短波分岔, 产生的斑图分别为时序振荡态、静态图灵斑图以及时空振荡斑图[1,4].

近年来, 时空振荡斑图由于其丰富的动力学行为而越来越受到人们的关注[2,8,9].产生时空振荡斑图的方式有很多, 除了上面提到的短波分岔外, 不同分岔引发的模式之间的相互作用也可以形成时间上和空间上均具有周期性的斑图结构.例如当系统处于图灵和霍普夫切空间点时, 就可以形成行波[10]、驻波[11]、振荡四边形[12]、振荡六边形、螺旋波等[13].

耦合反应扩散系统是目前研究不同模式之间相互作用的一种最为常用的方法[14-16].例如,2002年Yang等[17]在研究两个图灵模相互作用时发现, 当一个图灵长波模的次谐波处于霍普夫区域时, 通过次谐波图灵模与霍普夫模相互作用, 可以形成一种由三套子六边形相互嵌套的弛豫型振荡的振荡六边形斑图, 其中三个子六边形之间的相位相差 2 π/3.2003年, 该小组还详细研究了图灵和短波模之间的相互作用, 获得了一种眨眼六边形斑图, 同样, 该斑图由三套相位相差 2 π/3 的六边形相互嵌套而成, 但是其振荡属于正弦振荡.除了眨眼六边形斑图外, 还获得了各种相互嵌套的其他类型的时空振荡斑图, 例如螺旋波和同心圆波混合但是呈现六边形分布的时空斑图[18].Anguelov和Stoltz[19]对比了单层和双层耦合反应扩散系统中图灵-霍普夫相互作用的区别, 得到了反六边形振荡以及反条纹振荡斑图, 并给出了各种时空斑图的空间相图.最近, Li等[20]通过非线性耦合双层反应扩散模型, 获得了振荡四边形斑图.Pal等[21]研究发现在两个耦合反应系统中也存在着时空反共振现象.

然而, 这些研究大多集中于相同反应扩散模型之间的耦合, 对于不同模型之间耦合的研究还比较少.在实际的非线性系统中, 不同系统的局部动力学行为是不同的[22-24].例如, 在生物系统中, 热带鲶鱼的皮肤结构主要由表层表皮和底层真皮组成.由于不同皮肤层是由不同类型的细胞构成的, 即不同皮肤层的局部动力学行为是不同的[22].在复杂网络系统中, 系统是由具有不同功能的动态单元(例如大脑中的神经元或社会中的个体)间的耦合相互作用而成[23].介质阻挡放电系统也是一种典型的多层耦合斑图形成系统[24-26].它由放电层和电介质层组成, 放电产生的体电荷部分积累到介质表面形成表面电荷, 反过来表面电荷也会直接影响体放电, 两者相互耦合, 从而可以产生丰富多彩的时空斑图, 例如, 四边形斑图、六边形斑图、超点阵斑图、螺旋波斑图等[27-29].研究表明, 气体放电系统一种特殊的反应扩散系统, 可以唯象地用反应扩散模型来描述[25,26,30].为了进一步弄清介质阻挡放电系统中各类时空斑图形成的机理, 尤其是体电荷和表面电荷之间的耦合, 最近本小组线性耦合了两种不同的反应扩散模型, 获得了许多与实验相符合的模拟结果[31].除此以外, 还发现了一种新型的振荡六边形结构, 与以前发现的振荡六边形不同, 该振荡六边形是由图灵模与处于霍普夫区域的高阶模相互作用而成的, 但耦合机制还不甚清楚[31].为了进一步研究其动力学行为以及产生条件, 本文通过线性耦合两个不同的反应扩散模型, 系统研究了超临界图灵模、次临界图灵模、高阶图灵模以及霍普夫模之间的相互作用, 并分析了它们对振荡图灵斑图形成机理和时空动力学行为的影响.

2 物理模型

线性耦合两个不同的反应扩散模型, 构建一个非对称的双层耦合反应扩散系统模型.在无量纲的情况下, 该模型具有以下形式[17]:

其中, ( u1,v1) 和 ( u2,v2) 分别表示第一层和第二层的子反应扩散系统, 这里u和v分别为反应扩散系统中的活化子浓度和阻塞子浓度, Du和 Dv分别为它们的扩散系数, ∇2为拉普拉斯算符.α 为两个子系统活化子之间的耦合强度, 方程 f (u,v) 和g(u,v)为每层的子反应扩散系统的局部动力学方程, 不同的反应系统有着不同的具体表达式.本文选取Brusselator模型和Lengyel-Epstein模型来分别充当两个子系统.在无量纲的情况下, Brusselator模型的局部动力学方程具有下列形式[32]:

其中, 系数a和b是系统的控制参数.

Lengyel-Epstein模型在无量纲情况下形式如下[33]:

式中, 系数c和d是系统的控制参数.对于子系统1和子系统2, 其均匀定态解分别为(u10,v10)=通过系统的色散关系可以分析失稳模的性质以及系统的分岔类型.对于模型以及数值算法的详细描述可参见文献[31].

图1(a)为双层线性耦合反应扩散系统的色散关系曲线图, 而图1(b)和图1(c)分别给出了单层Brusselator模型和单层Lengyel-Epstein模型的色散关系图.可以看出, 双层耦合系统中存在着两个图灵模式: 一个为超临界图灵模, 其波数较小(长波模), 我们称其为 k1, 相应的本征值大小为 h1, 此模式是由第二层子系统Lengyel-Epstein模型所激发的失稳模; 另一个为次临界图灵模, 其波数较大(短波模), 我们称其为 k2, 相应的本征值大小为 h2,此模式是由第一层子系统Brusselator模型所激发的稳定模.对于图灵模, 其波数大小与系统变量的扩散系数成反比, 而本征值高度则与系统内阻塞子和活化子的扩散系数之比正相关, 所以通过调节模型中各变量的扩散系数, 可以改变两个图灵模的大小和本征值高度.

图1 双层线性耦合系统以及单层系统的色散关系 (a) 双层线性耦合系统; (b) 单层Brusselator模型; (c) 单层Lengyel-Epstein模型.参数取值为: ( a,b)=(3,9) ,(c,d)=(15,9) , D u1=2.2 , D v1=4.0 , D u2=21.9 ,Dv2=400,α=0.15Fig.1.Dispersion relationship of two-layered linear coupling system and single layer system: (a) Two-layered linear coupling system; (b) single layer Brusselator model; (c) single layer Lengyel-Epstein model.Parameter: ( a,b)=(3,9) ,(c,d)=(15,9), D u1=2.2 , D v1=4.0 , D u2=21.9 ,Dv2=400 , α =0.15.

3 模拟结果与讨论

采用欧拉向前差分的方法进行积分, 数值模拟在一个含有 1 28×128 个空间格点的二维平面上进行, 取时间积分步长 Δ t=0.01 个时间单位, 取空间积分步长 Δ x=Δy=1.0 个空间单位, 扩散项在数值计算中选择五点差分格式.边界条件选用周期性边界条件, 初始条件为均匀定态加上一个很小的随机扰动.扩散参数 Du和 Dv的乘积决定图灵模波数k的大小, 而 Du和 Dv的相对大小决定图灵模的本征值高度.本文先确定波数k的大小, 再通过改变 Du和 Dv的相对大小来调节相应的本征值高度.

3.1 同步振荡六边形斑图的基本性质

图2给出了双层非对称耦合反应扩散系统在图1(a)色散关系下的同步振荡六边形斑图.图2(a)为该同步振荡六边形斑图的振幅分布图, 显然, 振幅呈现六边形阵列排布, 并且幅值大小跟第二层活化子浓度 u2成反比, 也就是说图灵模 u2的静态效应会抑制 u1的振幅.该振幅六边形阵列的空间波数为 k1=0.2 , 由于六边形阵列中的每一个斑点的振荡都是同相位的, 所以图2(b)给出了六边形阵列中任意一个斑点处(图2(a)中A点)活化子 u1随时间t的变化图, 可以看出该同步振荡六边形斑图呈现正弦式的振荡, 振荡周期 T =2.37 , 相应的圆频率 ω =2.66.从单层Brusselator模型的色散关系图(图1(b))可知, 圆频率 ω =2.66 对应的波数为0.33, 此波数约为因此我们判断时间振荡来源于模式.在单层Brusselator模型中, 此模式是次临界模, 原本是稳定的, 但是由于受到另一层Lengyel-Epstein系统中 k1模式的驱动, 而被激发出来了.也就是说, 受到Lengyel-Epstein系统中的超临界图灵模 k1的激励, Brusselator系统中同时激发了图灵模 k1和处于霍普夫区域的高阶超临界图灵模此两种模式相互作用形成了同步振荡六边形斑图.

图2 同步振荡六边形斑图(参数与图1相同) (a) 振幅分布; (b) A点 u1 的时间变化关系图; (c) 半个振荡周期内的斑图演化过程Fig.2.Synchronous oscillatory hexagon pattern: (a) Amplitude distribution; (b) time variation of u1 at position A; (c) evolution of pattern in half an oscillating period.Parameters are the same as those in Fig.1.

图2 (c)给出了该同步振荡六边形斑图在半个振荡周期内的时空演化过程.在 t1=1.48 时刻, 六边形亮斑结构中心处浓度最低, 外围浓度最高, 再向外浓度降低, 从而形成了一种具有圆环结构的超六边形斑图.在 t2=1.78 时刻, 中心斑点处浓度最高, 从中心向外浓度逐渐降低, 形成了一种超点阵六边形斑图.随着时间演化, 中心斑点处浓度开始降低, 在 t5=2.38 时刻, 从中心斑点向外, 浓度先降低再升高, 呈现白眼六边形斑图.在 t6=2.75 时刻, 中心斑点处浓度降至最低, 远远低于最外环浓度, 此时形成了简单的蜂窝状六边形斑图.

3.2 霍普夫模式对振荡六边形斑图的影响

为了研究霍普夫模式在同步振荡六边形斑图形成过程中所起的作用, 固定其他参数不变, 通过改变控制参数b来调节霍普夫模本征值的大小.图3(a)为不同参数b下的色散关系曲线图及其局部放大图.控制参数b的改变, 除了可以调节霍普夫本征值的大小, 还对图灵模 k2和 h2有影响, 对k1和 h1影响很小.但由于 h2的值非常小, 也就是说图灵模 k2的作用很微弱, 因此这里主要考虑图灵模k1与霍普夫之间的相互作用, 后面将详细讨论图灵模 k2和 h2对斑图的影响.在 8.5＜b＜9.33 范围内,霍普夫模本征值的大小随着参数b的增大呈现下降趋势; 在 9.33＜b＜10 范围内, 霍普夫模本征值的大小随着参数b的增大而增大.图3(b)为不同参数b下获得的时空斑图及相应的活化子浓度u1随时间t的变化图.当 b =8.5 时, 系统距离霍普夫分岔点还比较远, 霍普夫模本征值比较小, 此时图灵模 k1及其高阶模均处于图灵区域内, 此时系统产生了静态蜂窝状六边形斑图, 如图3(b1)所示.这与介质阻挡放电系统中的静态蜂窝状六边形斑图一致[34,35].在 8.86＜b＜9.32 范围内, 随着参数b的增大, 霍普夫模本征值的大小呈现下降趋势,但是高阶图灵模进入到霍普夫区域内, 此时系统产生了振荡六边形斑图(图3(b2)和图3(b3)),这就说明该振荡六边形斑图是由超临界图灵模k1与处于霍普夫区域的高阶模相互作用而成的.当b ＜ 9.23时, 振荡六边形具有单一的振荡周期T0, 如图3(b2)所示.此斑图与介质阻挡放电实验中得到的点环斑图相符合[5,28].当b ＞ 9.23时, 振荡六边形经历倍周期分岔而进入双倍周期振荡, 即T=2T0, 如图3(b3)所示.在 9.32＜b＜9.33 这个很小的范围内, 系统处于近似5倍周期振荡, 但每次振荡的周期各不相同, 即T=T1+T2+T3+T4+T5≈5T0, 尽管振幅依然呈现六边形阵列排布, 但是振荡相位各不相同, 即系统产生的是非同步六边形斑图, 如图3(b4)所示.在参数9.33＜b＞9.7范围内, 霍普夫模本征值的大小随着参数b的增大而增大, 霍普夫模式开始被激发.当b=9.5时, 非同步斑图的六边形空间结构开始破缺, 其振荡周期呈现为近似三倍周期, 且有T=T1+2T2≈3T0, 如图3(b5)所示.当 b ＞9.7 时, 霍普夫模式引起的振荡起主导作用, 此时无论是空间上还是时间上, 都进入到混沌状态(图3(b6)).目前人们已经在介质阻挡放电系统中实验观测到了辉光模式下的倍周期分岔和混沌现象[36,37], 尽管斑图模式下还未见有实验报道, 但是在等离子体流体模拟中获得了相似的结果[38].

图3 参数 b对时空斑图的影响 (a) 不同参数b下的色散关系曲线; (b1)静态蜂窝状六边形斑图, b =8.5 ; (b2)单倍周期同步振荡六边形斑图, b =9 ; (b3) 2倍周期非同步振荡六边形斑图, b =9.25 ; (b4) 5倍周期非同步振荡六边形斑图, b =9.33 ; (b5) 3倍周期非同步振荡斑图, b =9.5 ; (b6)时空混沌,b=10Fig.3.Influence of parameter bon spatio-temporal patterns: (a) Dispersion curves under different parameters b; (b1) static honeycomb hexagon pattern at b =8.5 ; (b2) single-period synchronous oscillation hexagon pattern at b =9 ; (b3) non-synchronous oscillation hexagon pattern of 2 times the period at b =9.25 ; (b4) non-synchronous oscillation hexagon pattern of 5 times the period at b=9.33 ; (b5) non-synchronous oscillation pattern of 3 times the period at b =9.5 ; (b6) spatio-temporal chaos at b =10.

3.3 图灵模强度 h2 对振荡六边形斑图的影响

由前文可知, 振荡六边形斑图是由图灵模 k1与其高阶模相互作用而成.接下来研究Brusselator系统的图灵模强度 h2对振荡六边形斑图的影响.分两种情况来讨论: 一是两个图灵模之间满足空间共振关系; 二是不满足空间共振关系.首先讨论不满足空间共振的情况.以波数比为1∶5为例, 固定参数 Du2=21.9 , Dv2=400 , 使图灵模始终保持 k1=0.2 , h1=0.3 , 并且固定耦合强度α=0.15.通过调节参数 Du1和 Dv1的大小, 保证图灵模 k2=1 不变, 而只改变图灵模强度 h2的大小,探究图灵模强度 h2对振荡六边形斑图的影响.图4给出了不同 h2下的各种时空斑图.从图4(b)可以看出: 在 h2＜-4.8 范围内, 形成的是静态图灵斑图;在 - 4.8＜h2＜-0.74 范围内, 形成的是动态斑图;在 h2＞-0.74 范围内, 形成的是静态图灵斑图.当h2=-5.5 时, 本征值高度 h2很低, 系统中的次临界图灵模 k2不能被激发出来, 且此时高阶图灵模处于图灵区域内, 所以系统形成了静态蜂窝状六边形斑图, 如图4(b1)所示.增大本征值高度h2至—3.44时, 系统中的高阶图灵模处于霍普夫区域内, 且其本征值高于次临界图灵模的本征值h2, 此时系统呈现如图4(b2)所示的同步振荡六边形斑图, 同样证实了振荡六边形斑图是由图灵模k1与其高阶模相互作用而成的.继续增大本征值高度 h2, 当 h2=-0.56 时, 次临界图灵模的本征模 h2高于高阶模尽管高阶模也处于霍普夫区域, 但是系统优先选择 k2, 因此在两个图灵模k1和 k2的相互作用下形成了一种具有两种空间尺度的静态超六边形斑图, 该斑图中每个晶胞均由四个相邻的暗点排列而成, 如图4(b3)所示.该斑图与介质阻挡放电实验中观察到的超六边形斑图相类似[39].继续升高 h2至—0.21, 与前面所述相同, 系统仍然优先选择 k2, 此时系统呈现六边形和条纹混合的静态叠加斑图, 晶胞整体呈现六边形排列, 其内部有两种结构, 一种是六边形结构, 另一种是条纹结构, 如图4(b4)所示.在 0 ＜h2＜0.6 范围内, 此时图灵模 k2转变为失稳模, 因此Brusselator模型的斑图类型主要由图灵模 k2决定, 但是依然受到图灵模 k1的影响, 因此系统中呈现斑点与条纹的静态混合斑图, 如图4(b5)所示.继续增大本征值高度h2, 系统呈现单一波长的静态条纹斑图, 即系统完全由其自身的本征模 k2所决定, 如图4(b6)所示.

图4 非共振时不同本征值 h 2 下的时空斑图 (a) 不同本征值 h 2 下的色散关系曲线; (b1) 静态蜂窝状六边形斑图, h 2=-5.50 ,Du1=5.25 , D v1=6.5 ; (b2) 同步振荡六边形斑图, h 2=-3.44 , D u1=3.62 , D v1=5 ; (b3) 超六边形斑图, h 2=-0.56 ,Du1=2.2, D v1=4.6 ; (b4) 叠加斑图, h 2=-0.21 , D u1=2.1 , D v1=4.8 ; (b5)混合斑图, h 2=0.26 , D u1=2.03 ,Dv1=5.3 ; (b6) 条纹斑图, h 2=0.62 , D u1=1.98 ,Dv1=5.8Fig.4.Complex patterns under different eigenvalues h 2 at non-resonance: (a) Dispersion curves under different eigenvalues h 2 ;(b1) static honeycomb hexagon pattern, h 2=-5.50 , D u1=5.25 , D v1=6.5 ; (b2) synchronous oscillation hexagon pattern,h2=-3.44 , D u1=3.62 , D v1=5 ; (b3) super-hexagon pattern, h 2=-0.56 , D u1=2.2 , D v1=4.6 ; (b4) superposition pattern, h 2=-0.21 , D u1=2.1 , D v1=4.8 ; (b5) hybrid pattern, h 2=0.26 , D u1=2.03 , D v1=5.3 ; (b6) stripe pattern,h2=0.62, D u1=1.98 , D v1=5.8.

由此可知, 在非共振条件下, 要想形成振荡六边形斑图, 除了系统中的高阶图灵模处于霍普夫区域内外, 还需满足Brusselator模型中的次临界图灵模 k2对应的本征值高度低于处于霍普夫区域的高阶图灵模的本征值高度.此外, 在某些参数条件下, 处于霍普夫区域模式也可能被激发, 它与基模相互作用可以形成振荡四边形斑图, 但是该振荡四边形是不稳定的, 最终会演变为振荡六边形斑图.

接下来讨论图灵模满足空间共振的情况.固定参数 Du2=21.9 , Dv2=400 , 保证图灵模大小k1=0.2 及其相应的本征值高度 h1=0.3 不变, 并且固定耦合强度 α =0.15.通过调节参数 Du1和Dv1的大小, 探究图灵模强度 h2对振荡六边形斑图的影响.图5给出了不同 h2下的各种时空斑图.当h2＜-1.5 时, 本征值高度 h2比较低, 且此时高阶模处于图灵区域内, 因此在超临界图灵模k1的调制下, 系统形成了静态蜂窝状六边形斑图(图5(b1)).增大本征值高度 h2, 在-1.5＜h2＜0范围内, 两个子系统中均呈现六边形阵列排布, 对称性相同, 满足空间共振关系.此时虽然高阶图灵模被激发, 并高于次临界图灵模 k2, 但是由于空间共振作用, 系统依然优先选择 k1和 k2两个图灵模间的相互作用, 因此在系统中呈现静态黑眼斑图, 如图5(b2)所示.在 h2＞ 0范围内, 图灵模k2为失稳模, 此时系统完全由其自身的本征模 k2决定, 形成了如图5(b3)所示的具有单一波长的条纹斑图, 该条纹斑图与介质阻挡放电系统中观察到的条纹斑图相符合[25,28].

图5 共振时不同本征值 h 2 下的时空斑图 (a) 不同本征值 h 2 下的色散关系曲线; (b1)静态蜂窝状六边形斑图, h 2=-1.7 ,Du1=14.55, D v1=24 ; (b2) 静态黑眼斑图, h 2=-1.1 , D u1=13.6 , D v1=25 ; (b3) 条纹斑图, h 2=1.0 , D u1=12 ,Dv1=40Fig.5.Spatiotemporal patterns under different eigenvalues h 2 at resonance: (a) Dispersion curves under different eigenvalues h 2 ;(b1) static honeycomb hexagon pattern, h 2=-1.7 , D u1=14.55 , D v1=24 ; (b2) static black eye pattern, h 2=-1.1 ,Du1=13.6, D v1=25 ; (b3) stripe pattern, h 2=1.0 , D u1=12 , D v1=40.

3.4 耦合强度对振荡六边形斑图的影响

耦合强度也是影响不同模式之间相互作用的一个重要因素.图6(a)给出了不同耦合强度下的色散关系曲线图及其局部放大图.可以看出, 霍普夫模本征值的大小随着耦合强度的增大先减小后增大.图6(b)给出了不同耦合强度下获得的相应的时空斑图.当耦合强度 α =0.1 时, 与图3(b1)的情况类似, 系统距离霍普夫分岔点还比较远, 霍普夫模本征值比较小, 且高阶图灵模处于图灵区域内, 因此系统形成静态蜂窝状六边形斑图, 如图6(b1)所示.在 0.12＜α＜0.42 范围内, 霍普夫模本征值的大小随着耦合强度的增大呈现下降趋势, 但是高阶图灵模进入到霍普夫区域内, 系统产生了振荡六边形斑图(图6(b2)和图6(b3)).这再一次证实了振荡六边形斑图是由图灵模 k1与其高阶模相互作用而成.继续增大耦合强度,在 0.42＜ α ＜1 范围内, 系统中的 k1由原来的超临界图灵模变为了次临界图灵模, 系统呈现均匀态.因此系统要想产生振荡六边形斑图, k1必须是超临界图灵模, 且形成六边形阵列排布.

图6 不同耦合强度下的时空斑图(其他参数同图1) (a) 不同耦合强度下的色散关系曲线图; (b1)静态蜂窝状六边形斑图,α=0.1 ; (b2)同步振荡六边形斑图, α =0.17 ; (b3)非同步振荡六边形斑图, α =0.3 ; (b4)均匀态,α=0.6Fig.6.Spatio-temporal patterns under different coupling intensities at non-resonance: (a) Dispersion relationship curve diagram under different coupling intensities; (b1) static honeycomb hexagon pattern at α =0.1 ; (b2) synchronous oscillation hexagon pattern at α =0.17 ; (b3) non-synchronous oscillation hexagon pattern at α =0.3 ; (b4) uniform state at α =0.6.Other parameters are the same as those in Fig.1.

3.5 振荡黑眼斑图

接下来讨论两个图灵模 k1, k2以及霍普夫模共同作用下生成的振荡复杂斑图.图7给出了两个图灵模满足空间共振时获得的振荡黑眼斑图.色散关系如图7(a)所示.图7(b)为该振荡黑眼斑图的振幅分布图, 可以看出, 振幅呈现六边形阵列排布,且振幅中心斑点处的浓度远小于外围浓度.由于该振荡黑眼斑图为非同步振荡, 所以在图7(c)给出了3个不同斑点处(分别标记为A, B, C )活化子u1随时间t的变化图, 可以看出该振荡黑眼斑图呈现正弦型时间振荡, 3个位置处的振荡周期相同,均为 T =2.84 , 但是相位不相同, A, B, C相位依次延迟, 且3点的振荡幅度也不相同, A点的振荡幅度最大.图7(d)为该振荡黑眼斑图的时空演化过程, 显然, 在振荡过程中, 各个斑点为非同步振荡,但在整体上呈现六边形阵列排布.

图7 振荡黑眼斑图( ( a,b)=(3,10.5) , ( c,d)=(15,9) , D u1=15 , D v1=23 , D u2=21.9 , D v2=400 , α =0.45 ) (a) 耦合系统的色散关系曲线; (b) 振幅分布; (c) 3个位置处 u1 的时间变化关系图; (d) 斑图演化过程Fig.7.Oscillatory black-eye pattern ( ( a,b)=(3,10.5) , ( c,d)=(15,9) , D u1=15 , D v1=23 , D u2=21.9 , D v2=400 ,α=0.45 ): (a) Dispersion curve of coupled system; (b) amplitude distribution; (c) time variation of u1 at three positions; (d) evolution of pattern.

该振荡黑眼斑图的形成机理与图5(b2)的黑眼斑图有所不同.与图5(b2)相比, 在这种情况下,系统经历了霍普夫分岔, 激发了霍普夫模式, 因此该振荡黑眼斑图是在超临界图灵模 k1、次临界图灵模 k2以及霍普夫模的共同作用下形成的.

4 结 论

通过线性耦合Brusselator模型和Lengyel-Epstein模型, 数值研究了双层耦合非对称反应扩散系统中振荡图灵斑图的动力学行为.模拟结果表明, 时空振荡图灵斑图的形成和选择取决于霍普夫模、图灵模及其高阶模之间的相互作用.Lengyel-Epstein模型中激发的超临界图灵模 k1, 不仅在Brusselator模型中产生了图灵模 k1, 还激励了处于霍普夫区域的次临界模在这两个模式的相互作用下产生了同步振荡六边形斑图.此振荡六边形斑图形成的条件是Brusselator模型中的次临界图灵模 k2的本征值高度低于处于霍普夫区域的高阶图灵模的本征值高度, 且两个图灵模之间不存在空间共振关系.随着控制参数b的增加,该振荡六边形斑图首先经历倍周期分岔进入双倍振荡周期, 经历多倍振荡周期后, 在霍普夫模式的参与下, 最终进入时空混沌态.当系统的两个图灵模满足空间共振时, 系统优先选择空间共振模式,从而产生静态超点阵斑图.霍普夫模和图灵模共同作用下只能产生非同步振荡图灵斑图, 例如振荡黑眼斑图.此外, 耦合强度也是影响振荡图灵斑图的一个重要因素.耦合强度不仅可以影响斑图的振幅大小, 还影响空间结构的具体形式, 随着耦合强度逐渐增大, 系统依次经历了静态蜂窝状六边形斑图、同步振荡六边形斑图以及非同步振荡六边形斑图, 最终进入均匀状态.

本文的研究结果有助于人们进一步认知各种非线性系统中不同模式之间的相互作用及相应的斑图选择机理.