劉若琪 賈萌萌 范偉麗 賀亞峰 劉富成?

1) (河北大學(xué)物理科學(xué)與技術(shù)學(xué)院,保定 071002)

2) (河北大學(xué)環(huán)境工程研究所,保定 071002)

擴(kuò)散在圖靈斑圖形成和演化過程中起到了至為關(guān)鍵的作用.本文采用雙變量反應(yīng)擴(kuò)散模型,數(shù)值研究了異質(zhì)環(huán)境下各向異性擴(kuò)散對(duì)圖靈斑圖的影響.模擬結(jié)果表明: 各向異性擴(kuò)散程度較大時(shí),系統(tǒng)斑圖呈條紋斑圖,其空間取向由偏離分岔點(diǎn)程度大的擴(kuò)散系數(shù)所決定;各向異性擴(kuò)散程度較小時(shí),斑圖與各向同性擴(kuò)散情況相同.當(dāng)擴(kuò)散系數(shù)在空間中線性增長時(shí),由于空間異質(zhì)性的影響,不同區(qū)域可以選擇不同的斑圖類型,斑圖之間存在競爭.當(dāng)擴(kuò)散系數(shù)受到一維周期性函數(shù)調(diào)控時(shí),斑圖類型和波長由調(diào)控函數(shù)的波長以及系統(tǒng)本征波長所決定.當(dāng)調(diào)制波長大于本征波長時(shí),系統(tǒng)可呈現(xiàn)不同類型交替的雙尺度混合斑圖.此外,還發(fā)現(xiàn)非對(duì)角擴(kuò)散系數(shù)D 不僅影響了斑圖的選擇機(jī)制,還擴(kuò)展了圖靈空間的參數(shù)范圍.

1 引言

自組織斑圖是在空間和時(shí)間上具有某些規(guī)律性的非均勻宏觀結(jié)構(gòu),其普遍存在于自然界以及各種非線性系統(tǒng)[1?4].圖靈機(jī)制是理解斑圖形成的重要理論框架之一,該機(jī)制不僅對(duì)于理解各個(gè)領(lǐng)域內(nèi)的自組織斑圖現(xiàn)象具有重要作用[5],還在很多領(lǐng)域內(nèi)具有潛在的實(shí)際應(yīng)用價(jià)值,例如新型材料設(shè)計(jì)[1]、粒子運(yùn)輸輔助[6]、水質(zhì)凈化[7]以及蝗蟲集體運(yùn)動(dòng)預(yù)測[8]等.

圖靈不穩(wěn)定性本質(zhì)上是由于擴(kuò)散所引起的不穩(wěn)定性,因而擴(kuò)散過程是影響圖靈機(jī)制的一個(gè)重要因素[9].自20 世紀(jì)90 年代初期歐陽頎等[10]首次在實(shí)驗(yàn)上驗(yàn)證圖靈斑圖以來,研究者們對(duì)于正常擴(kuò)散下圖靈斑圖的形成機(jī)理已經(jīng)有了較為全面的認(rèn)識(shí)[9?12],在正常擴(kuò)散下,粒子遷移滿足完全隨機(jī)行走理論,擴(kuò)散過程與方向無關(guān),呈現(xiàn)各向同性的性質(zhì).實(shí)驗(yàn)和理論研究表明,現(xiàn)實(shí)中許多非線性物理系統(tǒng)和生物系統(tǒng)的擴(kuò)散過程大多不滿足完全隨機(jī)行走理論,其擴(kuò)散具有方向性,如分子在多孔介質(zhì)[13]、膠體[14]以及器官組織[15]中均表現(xiàn)為各向異性擴(kuò)散.各向異性擴(kuò)散作為調(diào)節(jié)斑圖形態(tài)的內(nèi)在因素之一,具有豐富的研究內(nèi)涵[16?18].Bar 等[19]發(fā)現(xiàn)各向異性雙穩(wěn)態(tài)反應(yīng)擴(kuò)散介質(zhì)中存在一種新的動(dòng)力學(xué)狀態(tài)—分層時(shí)空混沌.He 和Liu[20]研究了三維各向異性雙穩(wěn)態(tài)介質(zhì)中化學(xué)波前的動(dòng)力學(xué),發(fā)現(xiàn)足夠強(qiáng)的各向異性會(huì)引起動(dòng)力學(xué)不穩(wěn)定性并導(dǎo)致波前破裂.最近,Cui 等[21]在粗粒化模型中應(yīng)用化學(xué)波動(dòng)力學(xué)性質(zhì)成功地實(shí)現(xiàn)了納米顆粒的分離.在存在隨機(jī)性和各向異性的情況下,圖靈模式形成機(jī)制具有較強(qiáng)的魯棒性[22],并且單一擴(kuò)散方向滿足圖靈不穩(wěn)定性的發(fā)生條件時(shí),圖靈模式仍然存在[23].圖靈模式中條紋的取向由物質(zhì)各向異性擴(kuò)散的程度所決定,且各向異性提高條紋的穩(wěn)定性[17,24?26].Gao 等[18]發(fā)現(xiàn)當(dāng)各向異性擴(kuò)散在與局部反應(yīng)耦合過程中起主導(dǎo)作用時(shí)可以產(chǎn)生一種新的圖靈分岔并導(dǎo)致超晶格斑圖的形成.

此外,現(xiàn)實(shí)斑圖形成系統(tǒng)一般都具有非均勻性(異質(zhì)性)的特點(diǎn).研究表明,空間異質(zhì)性是造成時(shí)空斑圖多樣性的主要原因之一[27?29].非均勻擴(kuò)散系數(shù)是反應(yīng)擴(kuò)散系統(tǒng)空間異質(zhì)性最為常見的一種表現(xiàn)形式.2005 年,Page 等[30,31]發(fā)現(xiàn)擴(kuò)散系數(shù)的非均勻性可以擴(kuò)展圖靈空間的參數(shù)范圍,增加斑圖的復(fù)雜度以及穩(wěn)定性.胡文勇等[32,33]分析了在局域濃度擴(kuò)散系數(shù)調(diào)控下圖靈斑圖的形成機(jī)理,發(fā)現(xiàn)通過改變調(diào)控系數(shù)可實(shí)現(xiàn)六角斑點(diǎn)和條紋之間的轉(zhuǎn)化,且點(diǎn)條比率強(qiáng)烈依賴于調(diào)控系數(shù).Chen 等[34]在圖靈系統(tǒng)中施加一個(gè)圓偏振電場,也實(shí)現(xiàn)了六角斑點(diǎn)和條紋之間的轉(zhuǎn)化.最近,Van Gorder[27]詳細(xì)闡述了反應(yīng)擴(kuò)散系統(tǒng)中空間異質(zhì)性對(duì)圖靈斑圖的作用機(jī)理.在上述研究基礎(chǔ)上,本文采用雙變量反應(yīng)擴(kuò)散模型,系統(tǒng)研究了異質(zhì)環(huán)境下各向異性擴(kuò)散對(duì)圖靈斑圖的影響.

2 物理模型

將各向異性擴(kuò)散引入到最為常用的一種反應(yīng)擴(kuò)散模型—Brusselator 模型[35],在無量綱的情況下,該模型具有以下形式:

其中,u和v分別為活化子和阻滯子濃度;a和b為系統(tǒng)局部動(dòng)力學(xué)的控制參數(shù);和分別為物質(zhì)i沿x方向和y方向的擴(kuò)散系數(shù),i代表活化子u和阻滯子v,當(dāng)時(shí),系統(tǒng)為各向同性擴(kuò)散.由于反應(yīng)與擴(kuò)散的耦合,在參數(shù)空間的不同區(qū)域內(nèi)系統(tǒng)可以經(jīng)歷不同的分岔,從而產(chǎn)生不同的時(shí)空斑圖.均勻定態(tài)解(u0,v0)=(a,b/a)經(jīng)歷圖靈不穩(wěn)定的條件是阻滯子的擴(kuò)散系數(shù)大于活化子的擴(kuò)散系數(shù),相應(yīng)的圖靈斑圖具有本征波數(shù).固定參數(shù)a=3,b=9,Dv=10 ,選擇Du作為控制參數(shù),當(dāng)其小于臨界值Du≤Dc=4.44 時(shí),滿足圖靈分岔?xiàng)l件.圖1(a)為各向同性擴(kuò)散中不同擴(kuò)散系數(shù)下的色散關(guān)系曲線,虛線為特征值的實(shí)部 Re(λ),代表了圖靈模的線性增長率,Re(λ) 越大,圖靈模增長越快;實(shí)線為特征值的虛部 Im(λ),對(duì)應(yīng)模式的振蕩頻率,由于不同擴(kuò)散系數(shù)下特征值的虛部Im(λ)相差非常小,因此圖中僅給出了Du=4 時(shí)的情況作為參考.隨著擴(kuò)散系數(shù)Du的減小,系統(tǒng)逐漸遠(yuǎn)離圖靈分岔點(diǎn),圖靈模波數(shù)及其對(duì)應(yīng)的 Re(λ) 也都隨之增大.

圖1 各向同性擴(kuò)散下的色散曲線圖及其相應(yīng)的圖靈斑圖 (a)色散曲線;(b)點(diǎn)狀六邊形斑圖,Du=2;(c)條紋斑圖,Du=3 ;(d)蜂窩狀六邊形斑圖,Du=4Fig.1.Dispersion curves under isotropic diffusion and their corresponding Turing patterns: (a) Dispersion curves;(b) hexagonal spot pattern,Du=2;(c) stripe pattern,Du=3;(d) honeycomb hexagonal pattern,Du=4 .

所有模擬采用歐拉向前差分的方法進(jìn)行積分,模擬區(qū)域包含 128×128 個(gè)空間格點(diǎn),擴(kuò)散項(xiàng)在數(shù)值計(jì)算中選擇五點(diǎn)差分格式,在顯式格式下方程(1)的離散形式為

取時(shí)間步長 ?t=0.01 個(gè)時(shí)間單位,空間步長?x=?y=1.0個(gè)空間單位.采用更小的空間步長和時(shí)間步長獲得了相同的模擬結(jié)果.采用均勻定態(tài)附近的隨機(jī)擾動(dòng)為初始條件,邊界條件選用周期性邊界條件.

3 模擬結(jié)果與討論

圖1(b)—(d)給出了各向同性擴(kuò)散下獲得的3 種圖靈斑圖.偏離分岔點(diǎn)的程度是模式選擇過程中的一個(gè)關(guān)鍵因素[36].當(dāng)Du=2 時(shí),圖靈模的Re(λ)足夠高,系統(tǒng)遠(yuǎn)離主分岔點(diǎn),產(chǎn)生一種總相位為零的六邊形斑圖—點(diǎn)狀六邊形斑圖(圖1(b)),其波長λ約為7.3.隨著參數(shù)Du的增大,系統(tǒng)接近主分岔點(diǎn),產(chǎn)生波長λ=7.6 的條紋斑圖(圖1(c)),而不是點(diǎn)狀六邊形斑圖.當(dāng)Du=4 時(shí),在圖靈分岔點(diǎn)附近獲得總相位為 π 的六邊形圖案(圖1(d)),表示為蜂窩六邊形圖案,此斑圖波長λ=10.7 .

各向異性擴(kuò)散對(duì)圖靈斑圖具有影響.研究發(fā)現(xiàn)活化子u和阻滯子v的各向異性擴(kuò)散對(duì)斑圖的影響規(guī)律一致,為簡便起見,本文只給出了活化子的情況,即只考慮了活化子為各向異性擴(kuò)散的情況.保持阻滯子的擴(kuò)散系數(shù)Dv=10 不變,通過改變活化子在不同方向上的擴(kuò)散系數(shù)和的大小來研究各向異性擴(kuò)散對(duì)圖靈斑圖的影響.為了系統(tǒng)研究各向異性擴(kuò)散系數(shù)對(duì)圖靈斑圖的影響,分3 種情況來討論: 1)空間均勻分布;2)空間線性增長;3)空間周期性變化.

圖2 給出了均勻環(huán)境中兩種不同各向異性擴(kuò)散情況下獲得的條紋斑圖,其中圖2(a)為擴(kuò)散系數(shù)滿足圖靈分岔?xiàng)l件,而不滿足時(shí)的圖靈條紋斑圖.當(dāng)水平x方向上滿足圖靈分岔?xiàng)l件,而豎直y方向不滿足圖靈分岔?xiàng)l件時(shí),系統(tǒng)僅在x方向上激發(fā)圖靈模,因而產(chǎn)生豎直取向的條紋斑圖(圖2(a)),其本征波長λ=8,略大于各向同性擴(kuò)散情況下的本征波長(圖1(b)).圖2(b)為水平和豎直兩個(gè)方向均滿足圖靈分岔時(shí)的條紋斑圖,由于,即水平方向上偏離圖靈分岔點(diǎn)較遠(yuǎn),相應(yīng)的圖靈模的增長率較大,因此圖靈模的波數(shù)沿著x方向,即條紋取向?yàn)樨Q直方向.值得說明的是,此時(shí)條紋的本征波長λ=6.4,與圖2(a)相比其波長較短,這意味著y方向上的擴(kuò)散系數(shù)也會(huì)影響本征波數(shù)的大小.總而言之,條紋斑圖的取向取決于擴(kuò)散系數(shù)在不同方向上偏離分岔點(diǎn)的程度,擴(kuò)散系數(shù)越小偏離程度越高,則此方向即為條紋斑圖的取向方向.

在現(xiàn)實(shí)世界中,絕大多數(shù)非線性系統(tǒng)都處于非均勻的環(huán)境中,斑圖形成過程受到空間異質(zhì)性的影響.因此異質(zhì)性對(duì)圖靈斑圖的影響機(jī)理一直是人們研究的熱點(diǎn)之一[27].圖3 給出了擴(kuò)散系數(shù)在x方向單調(diào)增加時(shí)的混合斑圖.隨著x的不斷增加,從低于逐步高于.在圖3(a)的中間區(qū)域,當(dāng)在1.89—2.02 范圍時(shí)表現(xiàn)為各向同性情況下的點(diǎn)狀六邊形,2.02—2.13 范圍內(nèi)為斑點(diǎn)與條紋的過渡區(qū);左側(cè)區(qū)域表現(xiàn)為豎直取向條紋;右側(cè)區(qū)域,表現(xiàn)為水平取向條紋.由于圖3(a)左右兩側(cè)大小不同,所以造成了左右兩側(cè)條紋的波長不同.在圖3(b)中,在3.97—4.17 范圍時(shí)表現(xiàn)為各向同性情況下的蜂窩狀六邊形,左側(cè)區(qū)域?yàn)樨Q直條紋,右側(cè)區(qū)域?yàn)樗綏l紋.圖3 結(jié)果表明,由于空間異質(zhì)性的影響,不同區(qū)域可以選擇不同的斑圖類型.在每個(gè)區(qū)域內(nèi),斑圖由滿足圖靈條件的圖靈模式之間相互作用而成.斑圖的類型由各向異性擴(kuò)散系數(shù)的相對(duì)大小決定,在各向異性程度較小的區(qū)域內(nèi)呈現(xiàn)六邊形斑圖,而在各向異性程度較大的區(qū)域內(nèi)呈現(xiàn)條紋斑圖.

周期性變化的各向異性擴(kuò)散系數(shù)是調(diào)控圖靈斑圖的一種重要方式.引入周期性調(diào)控函數(shù)后活化子在x方向上擴(kuò)散系數(shù)變?yōu)?γcos(2πx/λF)+Du0,其中γ為周期性調(diào)制函數(shù)的振幅,調(diào)控波長λF=L/n,n為整數(shù),Du0為無周期調(diào)控時(shí)的擴(kuò)散系數(shù).當(dāng)=Du0時(shí),系統(tǒng)斑圖的本征波長為λ0(=3時(shí)λ0=8).圖4 給出了周期性調(diào)制擴(kuò)散系數(shù)在不同調(diào)制參數(shù)下的圖靈斑圖.其中圖4(a)—(c)研究了調(diào)制波長對(duì)斑圖的影響.當(dāng)n=16 時(shí),此時(shí)調(diào)控波長(λF=8)與系統(tǒng)的本征波長相同,在周期函數(shù)調(diào)控下,在不同位置的大小不同,即各位置的各向異性程度不同,但由于調(diào)控函數(shù)周期較短,各向異性程度較小的區(qū)域太窄,因此只表現(xiàn)出濃度相間的條紋斑圖,如圖4(a)所示.當(dāng)n=8 時(shí),調(diào)控函數(shù)的波長增大,各向異性程度不同的區(qū)域隨之變寬,在各向異性程度較小的區(qū)域呈現(xiàn)點(diǎn)狀斑圖,而在各向異性程度較大的區(qū)域呈現(xiàn)條紋斑圖,因此得到了雙尺度點(diǎn)線相間斑圖,如圖4(b)所示,此時(shí)λa=8.4(λa為點(diǎn)狀斑圖波長),λb=16 (λb為條紋斑圖波長),其中λb=λF=2λ0.圖4(b)與地松鼠表皮上的點(diǎn)線相間斑圖十分相似,因此空間周期異質(zhì)性可能是動(dòng)物體表點(diǎn)線斑圖的形成原因[31].當(dāng)n=4時(shí),此時(shí)λF=4λ0,調(diào)控函數(shù)的周期相較于圖4(b)進(jìn)一步增大,各向異性程度不同的區(qū)域進(jìn)一步變寬,在各向異性程度較小的區(qū)域呈現(xiàn)出橫向拉長的點(diǎn)狀斑圖,而在各向異性程度較大的區(qū)域呈現(xiàn)濃度不同的條紋斑圖,得到了圖4(c)中的點(diǎn)線相間斑圖.保持圖4(a)中的其他條件不變,圖4(d)—(f)分別給出了參數(shù)Du0,γ和對(duì)圖靈斑圖的影響.在圖4(d)中Du0的增大使得不滿足圖靈條件,系統(tǒng)僅在y方向上激發(fā)圖靈模,此時(shí)調(diào)制不起作用,系統(tǒng)產(chǎn)生了水平取向的條紋斑圖.在圖4(e)中調(diào)制幅度γ的增大擴(kuò)大了的參數(shù)范圍,但是周期大小并沒有改變,因此表現(xiàn)出另外一種濃度周期變化的條紋斑圖.在圖4(f)中,改變使其不滿足圖靈條件,系統(tǒng)僅在x方向上激發(fā)圖靈模,因而產(chǎn)生豎直取向條紋斑圖.

圖5 為不滿足圖靈條件,為周期性函數(shù)時(shí)的復(fù)雜條紋斑圖.當(dāng)n=4 時(shí),雖然調(diào)控函數(shù)周期較大,各向異性程度不同的區(qū)域較寬,但由于不滿足圖靈條件,系統(tǒng)僅在x方向上激發(fā)圖靈模,因此產(chǎn)生濃度不同的豎直取向條紋斑圖,如圖5(a)所示.當(dāng)n=8 時(shí),調(diào)控函數(shù)周期變小,各向異性程度不同的區(qū)域變窄,周期數(shù)變多,因而得到了濃度相間的周期性條紋斑圖,如圖5(b)所示.當(dāng)n=16時(shí),調(diào)控函數(shù)周期較短,因此只能得到周期內(nèi)為單一濃度的條紋斑圖,如圖5(c)所示.

在一些特殊的各向異性介質(zhì)中,例如特殊形狀的多孔介質(zhì)[13]、大腦皮層[15]、磁場中的等離子體[37]、生物種群[38]等,物質(zhì)的擴(kuò)散系數(shù)為張量,在二維情況下其形式可以表示為矩陣形式.同樣為簡便起見[這里只考慮活化]子的各向異性擴(kuò)散情況,即Du=一般來說,擴(kuò)散矩陣中的非對(duì)角元素都是相等的,令.為便于研究非對(duì)角元素D對(duì)圖靈斑圖的影響,在x方向和y方向上取各向同性擴(kuò)散系數(shù),即研究各向同性擴(kuò)散下的情況.圖6(a)和圖6(b)是相同D值、不同Du值得到的斜條紋斑圖.與無非對(duì)角元素情況相比,當(dāng)非對(duì)角擴(kuò)散系數(shù)D=1 時(shí),系統(tǒng)選擇了斜條紋斑圖(圖6(a))而非點(diǎn)狀六邊形斑圖(圖1(b)),這意味著非對(duì)角擴(kuò)散系數(shù)D影響了斑圖的選擇機(jī)理.在各向同性擴(kuò)散中,當(dāng)Du=5時(shí),此時(shí)不滿足圖靈分岔?xiàng)l件,系統(tǒng)為均勻態(tài),但令D=1 后,系統(tǒng)由均勻態(tài)轉(zhuǎn)變?yōu)樾睏l紋斑圖,如圖6(b)所示,表明各向異性擴(kuò)散擴(kuò)展了圖靈空間的參數(shù)范圍.在各向同性擴(kuò)散基礎(chǔ)上,當(dāng)非對(duì)角擴(kuò)散系數(shù)D超過一個(gè)臨界值Dc時(shí),系統(tǒng)會(huì)經(jīng)歷圖靈分岔點(diǎn).圖6(c)給出了臨界擴(kuò)散系數(shù)Dc與各向同性擴(kuò)散系數(shù)Du的變化.隨著Du的增大,Dc隨之線性增大.值得說明的是,非對(duì)角擴(kuò)散引發(fā)的圖靈斑圖總是選擇斜條紋斑圖,其本征波長隨著非對(duì)角擴(kuò)散系數(shù)D的增大而不斷減小,但是臨界擴(kuò)散系數(shù)Dc所對(duì)應(yīng)的本征波長與Du無關(guān),如圖6(d)所示.

圖6 非對(duì)角擴(kuò)散系數(shù) D 對(duì)斑圖的影響 (a)斜條紋,Du=2,D=1;(b)斜條紋,Du=5,D=1;(c)擴(kuò)散系數(shù) Du 與臨界擴(kuò)散系數(shù) Dc的關(guān)系曲線圖;(d)擴(kuò)散系數(shù) D與波長 λ 的關(guān)系曲線圖Fig.6.Effect of off-diagonal diffusion coefficient D on the pattern: (a) Oblique stripes,Du=2,D=1 ;(b) oblique stripes,Du=5,D=1;(c) the relationship between the diffusion coefficient Duand the critical diffusion coefficient Dc ;(d) the relationship between the diffusion coefficient D and the wavelength λ .

4 結(jié)論

本文通過改變活化子在不同方向上的擴(kuò)散系數(shù)和研究了各向異性擴(kuò)散對(duì)圖靈斑圖的影響.模擬結(jié)果表明,當(dāng)單方向上的擴(kuò)散系數(shù)滿足圖靈條件即可產(chǎn)生圖靈條紋斑圖,其空間取向取決于擴(kuò)散系數(shù)在不同方向上偏離分岔點(diǎn)的程度,擴(kuò)散系數(shù)越小偏離程度越高,則此方向即為條紋斑圖的取向方向.當(dāng)擴(kuò)散系數(shù)線性增長時(shí),由于空間異質(zhì)性的影響,不同區(qū)域可以選擇不同的斑圖類型.在每個(gè)區(qū)域內(nèi),斑圖由滿足圖靈條件的圖靈模式之間相互作用而成.當(dāng)擴(kuò)散系數(shù)周期性變化時(shí),可形成周期性復(fù)雜斑圖.若調(diào)控函數(shù)的周期較小,各向異性程度較小的區(qū)域太窄,則僅呈現(xiàn)出條紋斑圖而不會(huì)呈現(xiàn)各向同性擴(kuò)散時(shí)的斑圖類型.非對(duì)角擴(kuò)散系數(shù)D不僅可以影響斑圖的選擇機(jī)理而且擴(kuò)展了圖靈空間的參數(shù)范圍.隨著擴(kuò)散系數(shù)Du的增大,臨界擴(kuò)散系數(shù)Dc線性增大.斜條紋斑圖的本征波長均隨著非對(duì)角擴(kuò)散系數(shù)D的增大而不斷減小.