徐 瑞，田美美，徐衍聰

(杭州師范大學理學院，浙江 杭州 310036)

帶有交叉擴散項的反應擴散系統(tǒng)的 Turing 不穩(wěn)定性

徐 瑞，田美美，徐衍聰

(杭州師范大學理學院，浙江 杭州 310036)

通過研究含有 Lotka-Volterra 捕食和被捕食動力學行為的反應擴散方程的交叉擴散項來解釋斑圖的形成機制. 借此來說明能導致部分斑圖出現(xiàn)的交叉耗散項在失穩(wěn)機制中的重要性.在穩(wěn)定區(qū)域的附近采用弱非線性分析方法來研究斑圖的振幅， 得到其 Stuart-Landau 規(guī)范型振幅方程. 最后，當斑圖作為行波的波前侵入?yún)^(qū)域時，可以得到 Ginzburg-Landau 規(guī)范型振幅方程，這一方程常常用來研究波的形狀和波速.

交叉耗散；Turing 不穩(wěn)定；振幅方程；五次 Stuart-Landau 方程；Ginzburg-Landau 方程

1 預備知識

斑圖動力學是非線性科學領域內的一類重要的分支.作為一門應用科學，它的研究內容廣泛地涉及物理學、力學、化學、生物學、生態(tài)學等各方面. 斑圖(Pattern)是在空間和時間上具有某種規(guī)律性的非均勻宏觀結構.它普遍存在于自然界中，形形色色的斑圖結構構成了多姿多彩、千姿百態(tài)的世界. 因而越來越多的人將目光轉向生物學領域，并用一個反應擴散模型成功地解釋某些生物的體表所顯現(xiàn)的圖紋，如斑馬身上的斑紋形成機理.

文章的主要目的在于研究下列反應擴散系統(tǒng)的斑圖形成過程：

(1)

這里的U(z,τ) 和V(z,τ)，z∈Ω?Rn分別是捕食者和被捕食者的種群密度,R和K分別是捕食者的增長率和捕捉能力，M是被捕食者的死亡率，γ12和γ21分別是捕食和捕獲效率并且Γ 是空間區(qū)域的大小. 兩個種群的空間運動是由簡單線性耗散項和源于假設有高密度的捕食者區(qū)域的非線性項來描述的.通過無量綱化，把系統(tǒng) (1) 化為關于變量u和v的系統(tǒng)：

(2)

這里有

(3)

和

(4)

以下的討論限制在一維空間Ω=[0,2π]上并運用齊次Neumann邊值條件討論以上系統(tǒng).系統(tǒng)(2)是依賴于擴散項[1]和交叉耗散項來產(chǎn)生分離影響和空間生態(tài)位的一類模型[2-4].交叉擴散項是用來引入說明一類種群的密度梯度來減少另一類種群的變遷情形.這就是說，對于一大類捕食-被捕食或競爭動力學沒有自動催化項的系統(tǒng)來看，無論擴散率是多少，經(jīng)典的擴散項的擴散都是不能充分揭示斑圖的形成原因：在這些情況下，交叉擴散被用來解釋斑圖形成機理[3,5-6].值得一提的是，交叉擴散項的引入是可以嚴格地通過錨定到微觀世界的自相一致推導得到[7]，或從單一品種的變異和分裂推導[8-9].系統(tǒng)(2)中介紹該類型的交叉擴散方面已經(jīng)在文[3,10-11]出現(xiàn)，但是在其他方面比如趨化[12]、生態(tài)學[13]、社會系統(tǒng)[14]、電子電路[15]、半導體中的漂移擴散[16]、化學反應[17]、等離子體中的湍流[2]、粒狀材料[18]和細胞分裂的腫瘤生長[19]都已經(jīng)被廣泛地通過數(shù)值計算進行過研究[20-22].

第二部分，采用線性穩(wěn)定性分析研究共存奇點的動力學行為，表明交叉擴散是斑圖形成的重要原因之一.第三部分，通過弱非線性分析，得到三次Stuart-Landau規(guī)范型振幅方程(但是，在次臨界區(qū)域，必須得到五次Staurt-Landau規(guī)范型振幅方程來研究)，以此預測形狀和圖案的振幅.在次臨界區(qū)域，那里的線性分析沒有時間的振蕩規(guī)定而只有Turing模式，并且發(fā)現(xiàn)了數(shù)字并存與極限環(huán)圖靈模式.這些發(fā)現(xiàn)同樣證實了這一結果，即使是沒有Hopf分支或波不穩(wěn)定,Turing不穩(wěn)定性與交叉擴散的捕食模型也可以產(chǎn)生時空振蕩的解[23-24].在文[25-26]中也發(fā)現(xiàn)了振蕩圖靈斑圖的類似模型.最后，第四部分，在空間域很大的情況下考慮空間調制斑圖，即波前解的存在性，并且這些斑圖模式是有次序地侵入到整個區(qū)域.

2 Turing 不穩(wěn)定

在本節(jié)中, 將研究系統(tǒng) (2) 的 Turing 不穩(wěn)定性.易知該系統(tǒng)的共存非邊界平衡點是(u0,v0)≡(1,r-γ), 這一平衡點當且僅當r-γ>0 時才具有生物學意義. 當它存在時, 平衡點是穩(wěn)定的, 可以是一個有吸引力的節(jié)點(當 0γ+γ2/4). 顯然，此時該系統(tǒng)不會出現(xiàn) Hopf 分支.因此考慮在(u0,v0) 附近線性化該系統(tǒng):

(5)

這里有

(6)

給出特征值為波數(shù)的函數(shù)色散關系如下：

(7)

這里有

(8)

不失一般性，文中假設 det(D)>0.

靜態(tài)狀況下的線性不穩(wěn)定的空間擾動要求對于k≠0 有R(λ(k))>0.由于(u0,v0) 在動力學上是穩(wěn)定的，故有 tr(J)<0, 而且tr(D)>0.因此，尋找這些當h(k2)<0 情況下的k值.若式(8)有正實部的特征根，需保證其系數(shù)異號，而此時，唯一的可能性是當h(k2)<0 時q<0，即

(9)

因此，當線性擴散項起穩(wěn)定作用的時候，唯一的潛在不穩(wěn)定機制是交叉擴散項的存在.令d=d21-d12，在k=kc臨界條件是

(10)

(11)

對應的臨界波數(shù)是

(12)

圖1 當r>γ時的Turing不穩(wěn)定區(qū)域Fig.1 The region for Turing instability as r>γ

(13)

3 弱非線性分析

在本節(jié)中，將對系統(tǒng)(2)應用多尺度方法進行弱非線性分析.取ε2=(d-dc)/dc為無量綱化的距離閾值(threshold).為簡便起見，這部分用d1代替d12、d2代替d21.

首先，引入新的尺度變換：

(14)

因此偏導數(shù) ?t→?t+ε2?T和 ?x→?x+ε?X作為振蕩振幅的時間尺度和空間尺度來解耦，在這個階段，暫時忽略空間調制緩慢的可能性，這將在第五部分考慮.

分離系統(tǒng)的線性與非線性部分后，改寫帶有擾動w=(wu,wv) 的系統(tǒng)(2)為如下形式：

(15)

這里有一個依賴于分支參數(shù)d的線性算子Ld=ΓJ+Dd2和一個雙線性算子 Qk，這里的， Qk代表了系統(tǒng)的非線性部分，定義為:

其中x=(xu,xv) 和y=(yu,yv).最后，公式(15)的最后兩項是非線性擴散項.

按下列方式展開d1、d2和w:

因此，線性算子Ld和 Qk可以被展開為：

這里Ldc=ΓJ+Ddc2，且.將上述的展開式分別代入式(15)中并收集ε的各次項，則可以得到一系列關于wi的方程.

在O(ε) 次項中，得到線性問題Ldcw1=0，它的解滿足 Neumann 邊值條件，則有

由于Ldc不作用慢尺度T上，因此這里的A(T) 是在這個水平上的任意斑圖的振幅. 向量ρ=(ρu,ρv) 定義到一個常數(shù)，并且能夠按下列方式規(guī)范化：

(16)

在O(ε2) 次項中，得到下列系統(tǒng)：

(17)

通過 Fredholm 更替原理，方程(17)有解當且僅當=0，這里的·,·是在L2(0,2π/kc) 下的標量積，且，因此

(18)

因此 Fredholm 更替條件自動滿足.因此，對于線性部分，方程 (17) 得到一下滿足 Neumann 邊界條件的解:

這里的w2i(i=0,2) 是下列線性系統(tǒng)的解：

(19)

(20)

在O(ε3) 次項中得到：

(21)

這里有

因此，根據(jù)方程 (21) 的可解性條件，得到 Stuart-Landau 規(guī)范型振幅方程：

(22)

這里的系數(shù)σ和L如下：

(23)

如果d>dc，斑圖在此參數(shù)區(qū)域內可形成. 容易證明系數(shù)σ總是正的.另一方面， Landau 常數(shù)L可以是正的或者負的，這取決于系統(tǒng)參數(shù)的值.因此 Stuart-Landau 的動力學行為可以被分為兩種性質不同的情況： 超臨界情況(L>0)，次臨界情況(L<0). 而L作為所有參數(shù)的函數(shù).

3.1 超臨界狀態(tài)

如果方程 (22) 系數(shù)σ和L都是正的，則出現(xiàn)超臨界分支. 在這種情況下 Stuart-Landau 方程的靜態(tài)平衡點是A并且解的時間漸近行為由下列方程給出:

(24)

3.2 次臨界狀態(tài)

當L<0 時，系統(tǒng)出現(xiàn)次臨界分支，在這種情況下，弱非線性展開式必須推高到五階.因此，引入了多個時間尺度T和T1, 令

取時間導數(shù)解耦為 ?t→?t+ε2?T+ε4?T1+…， 并將d1、d2和w展開直到ε的五次項.

將展開式代入方程(15)，直到O(ε3)，將會得到和第三部分一樣的方程.在O(ε3) 項，對式(21)利用可解條件=0 同樣可以得到振幅方程(22)，盡管這里T的導數(shù)是偏導數(shù).如果滿足，則解是

這里w3i(i=1,2,3) 分別是下列線性系統(tǒng)的解：

在O(ε4) 次項：

(25)

這里的w4i(i=0,…,4) 可由下列線性系統(tǒng)解得：

在O(ε5) 次項可得

(26)

(27)

這里的系數(shù)分別是：

最后將式(22)加到式(27)，得到五次 Stuart-Landau 規(guī)范型振幅方程：

(28)

其中的系數(shù)分別是：

(29)

圖2 次臨界狀態(tài)下的分支圖Fig.2 Bifurcation diagram for subcritical state as <0

4 波 前

當物理區(qū)域很大時，斑圖會連續(xù)地形成并以一種波前的狀態(tài)侵入到整個區(qū)域. 為了描述這一現(xiàn)象，不得不考慮調制空間相關性的快慢程度，引入了弱非線性分析中定義的緩慢的空間尺度 (在式(14)中已有定義).隨后一個確定性的算子就出現(xiàn)在式(15)中.在O(ε) 中，得到線性算子Ldcw1=0，這里的Ldc=ΓJ+Dd?xx是線性算子的快部分算子.它的解是

這里的ρ由式(16)給出，A和前面的一樣仍然是任意的.

在O(ε2) 次項中，得到方程

并且可解性條件自動滿足，因此它的解是

這里的w2i(i=0,2)和wX1是下列線性系統(tǒng)的解：

在O(ε3) 次項中，可以得到：

利用可解條件得到關于振幅A(X,T) 的實 Ginzburg-Laudau 規(guī)范型振幅方程如下：

這里的σ和L由式(23)給出，耗散系數(shù)v的關系式如下：

這里的Ψ 由式(18)給出，(·,·) 是標準的標量積. 實 Ginzburg-Laudau 方程描述整個區(qū)域的入侵模式，并且它的解是一個斑圖的包絡解.

5 結 論

通過線性穩(wěn)定性分析和弱非線性分析來說明交叉擴散是斑圖形成的必要成分，進而說明即使是沒有 Hopf 分支不穩(wěn)定或Turing 不穩(wěn)定，該模型同樣可以產(chǎn)生時空振蕩的解.從本文看，引入對偶的交叉擴散項，對 Turing 分支沒有明顯的影響，但是會影響 Turing 斑圖發(fā)生的參數(shù)區(qū)域. 但實際上，加入對偶擴散項后，除了共存的3個奇點外，新的系統(tǒng)產(chǎn)生了至少2個新的奇點，這也說明對偶擴散項對系統(tǒng)的動力學行為的影響是顯而易見的， 限于本文的研究內容，我們暫時不予研究， 將在隨后的工作中刻畫這2個新奇點引發(fā)系統(tǒng)產(chǎn)生的動力學行為.

[1] MULONE G, RIONERO S, WANG W. The effect of density-dependent dispersal on the stability of populations[J]. Fuel & Energy Abstracts, 2011, 74(74):4831-4846.

[2] DEL-CASTILLO-NEGRETE D, CARRERAS B A, LYNCH V. Front propagation and segregation in a reaction-diffusion model with cross-diffusion[J]. Physica D Nonlinear Phenomena, 2002,168-169(2):45-60.

[3] PENG R, WANG M, YANG G. Stationary patterns of the Holling-Tanner prey-predator model with diffusion and cross-diffusion[J]. Applied Mathematics & Computation, 2008, 196(2):570-577.

[4] SHIGESADA N, KAWASAKI K, TERAMOTO E. Spatial segregation of interacting species[J]. Journal of Theoretical Biology, 1979, 79(1):83-99.

[5] GAMBINO G, LOMBARDO M C, SAMMARTINO M. Turing instability and traveling fronts for a nonlinear reaction-diffusion system with cross-diffusion[J]. Mathematics & Computers in Simulation, 2012, 82(6):1112-1132.

[6] GAMBINO G, LOMBARDO M C, SAMMARTINO M. Pattern formation driven by cross-diffusion in a 2D domain[J]. Nonlinear Analysis-Real World Applications, 2013, 14(3): 1755-1779.

[7] FANELLI D, CIANCI C, PATTI F D. Turing instabilities in reaction-diffusion systems with cross diffusion[J]. The European Physical Journal B, 2013, 86(4):1-8.

[8] CONFORTO F, DESVILLETTES L. Rigorous passage to the limit in a system of reaction-diffusion equations towards a system including cross diffusions[J]. Communications in Mathematical Sciences, 2014, 12(3):457-472.

[9] GALIANO G. On a cross-diffusion population model deduced from mutation and splitting of a single species[J]. Computers & Mathematics with Applications, 2012, 64(6):1927-1936.

[10] DUBEY B, DAS B, HUSSAIN J. A predator-prey interaction model with self and cross-diffusion[J]. Ecological Modelling, 2001, 141(1):67-76.

[11] KUTO K, YAMADA Y. Multiple coexistence states for a prey-predator system with cross-diffusion[J]. Journal of Differential Equations, 2004, 197(2):315-348.

[12] KELLER E F, SEGEL L A. Model for chemotaxis.[J]. Journal of Theoretical Biology, 1971, 30(2):225-234.

[13] GILAD E, VON H J, PROVENZALE A, et al. A mathematical model of plants as ecosystem engineers[J]. Journal of Theoretical Biology, 2007, 244(4):680-691.

[14] EPSTEIN J M. Nonlinear dynamics, mathematical biology and social science[M]. Reading:Addison-Wesley, 1997.

[15] BILOTTA E, PANTANO P, STRANGES F. A gallery of Chua attractors: Part II[J]. International Journal of Bifurcation & Chaos, 2007, 17(2):293-380.

[16] LI C, ANSGAR J. Analysis of a parabolic cross-diffusion semiconductor model with electron-hole scattering[J]. Communications in Partial Differential Equations, 2007, 32(1):127-148.

[17] GAMBINO G, LOMBARDO M C, SAMMARTINO M, et al. Turing pattern formation in the Brusselator system with nonlinear diffusion[J]. Physical Review E, 2013, 88(1):042925.

[18] ARANSON I S, TSIMRING L S. Continuum theory of partially fluidized granular flows[J]. Physical Review E, 2001, 65(1):061303.

[19] SHERRATT J A. Wavefront propagation in a competition equation with a new motility term modelling contact inhibition between cell populations[J]. Proceedings of the Royal Society A, 2000, 456:2365-2386.

[20] BARRETT J W, BLOWEY J F. Finite element approximation of a nonlinear cross-diffusion population model[J]. Numerische Mathematik, 2004, 98(2):195-221.

[21] GALIANO G, GARZN M L, JüNGEL A. Semi-discretization in time and numerical convergence of solutions of a nonlinear cross-diffusion population model[J]. Numerische Mathematik, 2003, 93(4):655-673.

[22] GAMBINO G, LOMBARDO M C, SAMMARTINO M. A velocity-diffusion method for a Lotka-Volterra system with nonlinear cross and self-diffusion[J]. Applied Numerical Mathematics, 2009, 59(5):1059-1074.

[23] MADZVAMUSE A, NDAKWO H S, BARREIRA R. Cross-diffusion-driven instability for reaction-diffusion systems: analysis and simulations[J]. Journal of Mathematical Biology, 2015, 70(4):709-743.

[24] LI A W, JIN Z, LI L, et al. Emergence of oscillatory Turing patterns induced by cross diffusion in a predator-prey system[J]. International Journal of Modern Physics B, 2012, 26(31):3847-3856.

[26] TULUMELLO E, LOMBARDO M C, SAMMARTINO M. Cross-diffusion driven instability in a predator-prey system with cross-diffusion[J]. Acta Applicandae Mathematicae, 2014, 132(1):621-633.

Turing Instability in a Predator-prey System with Cross-diffusion

XU Rui, TIAN Meimei, XU Yancong

(School of Science, Hangzhou Normal University, Hangzhou 310036, China)

The pattern formation system was explained by studying the cross-diffusion in a reaction-diffusion system with Lotka-Volterra predator-prey kinetics. It was indicated that the cross-diffusion term that leads to the emergence of spatial patterns is responsible of the destabilizing mechanism. The amplitude of patterns were studied by weakly nonlinear analysis near marginal stability, and the Stuart-Landau amplitude equation was obtained. Finally, when the patterns invaded the domain as a travelling wavefront, the Ginzburg-Landau amplitude equation which was able to describe the shape and speed of the wave was derived.

cross-diffusion; Turing instability; amplitude equation, Quintic Stuart-Landau equation; Ginzburg-Landau equation.

2017-03-20

國家自然科學基金項目 (11671114); 杭州師范大學科研基金項目 (HNUEYT) .

徐衍聰(1970—)，男，教授，博士，主要從事微分動力系統(tǒng)分支研究. E-mail：yancongx@163.com

10.3969/j.issn.1674-232X.2017.04.004

O175.29 MSC2010:35K57；37L15

A

1674-232X(2017)04-0360-08