胡亞軍,閆廣武

(吉林大學(xué) 數(shù)學(xué)學(xué)院,長春 130012)

大規(guī)模氣泡湮滅的元胞自動機(jī)模擬

胡亞軍,閆廣武

(吉林大學(xué) 數(shù)學(xué)學(xué)院,長春 130012)

使用元胞自動機(jī)模型,并用正六邊蜂房結(jié)構(gòu)代替二維均勻氣泡,模擬氣泡載荷超過臨界載荷時大規(guī)模二維氣泡湮滅過程. 結(jié)果表明,該過程存在Zipf律.

元胞自動機(jī); 二維氣泡; Zipf律

自Bak等[1-2]提出自組織臨界概念以來,人們已經(jīng)發(fā)現(xiàn)這種現(xiàn)象具有某種普適性,如地震、 森林火災(zāi)、 雪崩和河流的形成等都是自組織臨界的[3-5]. 研究表明,對于臨界的自組織系統(tǒng),當(dāng)能量均勻緩慢地連續(xù)加入系統(tǒng)中時,所引起的能量耗散是突然的、 雪崩式的； 而對于處于臨界狀態(tài)的系統(tǒng),存在一個表示頻率與尺度關(guān)系的對數(shù)直線關(guān)系,即遵從Zipf律,分形結(jié)構(gòu)即為這樣的系統(tǒng). 自組織臨界系統(tǒng)還有很多,如對于水的成核汽化相變,大規(guī)模氣泡湮滅中可能存在自組織臨界現(xiàn)象. 基于此,本文應(yīng)用元胞自動機(jī)模型,模擬處于臨界狀態(tài)的氣泡破裂過程. 用正六角形的蜂房結(jié)構(gòu)(實(shí)際是正六角形的內(nèi)切圓)代替均勻氣泡. 外界能量的加入采用選定氣泡的微小重量增加,研究臨界重量與氣泡破裂數(shù)的關(guān)系. 結(jié)果表明,該過程是自組織臨界的.

1 元胞自動機(jī)模型

元胞自動機(jī)(cellular automata,又稱細(xì)胞自動機(jī),點(diǎn)格自動機(jī),簡稱CA)是空間、 時間和狀態(tài)變量完全離散的動力系統(tǒng)[6]. 許多相同的元胞(即點(diǎn)格)以均勻方式排列. 在每個格點(diǎn)上有一個處理器,所有處理器中都含有同樣多的狀態(tài),每個元胞都處于一系列有限狀態(tài)中的某種狀態(tài). 時間被離散成均勻的時刻,稱為時鐘. 當(dāng)時鐘滴答一下時,所有元胞的狀態(tài)都發(fā)生改變. 元胞自動機(jī)可用于模擬許多復(fù)雜系統(tǒng),特別是處于自組織臨界狀態(tài)的系統(tǒng),例如植被的演化[7]、 生態(tài)演化[8]和交通流動[9]等.

本文考慮二維面上的一層氣泡. 將氣泡所在空間用正六角形蜂房結(jié)構(gòu)代替,如圖1所示,給出二維空間的蜂房剖分,每個正六角形的內(nèi)切圓是一個氣泡. 每個氣泡的中心將組成與蜂房的對偶網(wǎng)格Di,j,每點(diǎn)(i,j)存在不同的水重級別,稱為該點(diǎn)的狀態(tài)值. 氣泡破裂時的狀態(tài)值稱為臨界狀態(tài)值,其大于各級別的狀態(tài)值. 當(dāng)狀態(tài)值為零時表示氣泡已破裂狀態(tài). 在下一時刻,由于水重的增加使得狀態(tài)值增加,超過時氣泡發(fā)生破裂,并將水重均分給周圍可能存在的氣泡,從而構(gòu)成一個多狀態(tài)的元胞自動機(jī).

圖1 蜂房網(wǎng)格(A)和單胞及其對偶網(wǎng)格(B)Fig.1 Honeycomb grid (A) and unit cell and its dual grid (B)

1.1 非飽和狀態(tài) 非飽和狀態(tài)是指某點(diǎn)周圍當(dāng)前時刻沒有破裂氣泡,則此時元胞自動機(jī)的鄰居數(shù)為6. 本文選擇6行5列的蜂房,如圖2所示. 取氣泡的臨界狀態(tài)值為11,以3行3列位置上的氣泡為起始點(diǎn),當(dāng)起始點(diǎn)的狀態(tài)數(shù)增加2時,某個氣泡的狀態(tài)數(shù)達(dá)到11發(fā)生破裂,狀態(tài)值立刻變?yōu)?. 然后將其狀態(tài)值等分添加到周圍,如圖3所示.

圖2 6×5蜂房及級別點(diǎn)陣Fig.2 Honeycomb lattice 6 rows×5 columns and level

圖3 下一時刻狀態(tài)(A)和破裂氣泡(B)(實(shí)心圓表示破裂的氣泡)Fig.3 State (A) and burst bubbles (B) in the next moment (including solid circle said burst bubbles)

1.2 飽和狀態(tài) 當(dāng)某氣泡的周圍存在破裂氣泡時,稱該點(diǎn)處于飽和狀態(tài). 此時需要確定非破裂氣泡數(shù)量. 例如,當(dāng)周圍有一個破裂氣泡時,則需要將狀態(tài)5均分.

1.3 氣泡塌陷 氣泡塌陷是指當(dāng)超過臨界狀態(tài)時,氣泡沒有破裂而容易破裂,即它的臨界狀態(tài)值變小. 由于本文考慮的氣泡是均勻且體積不變的,因此,從體積上看不出該過程,可以理解成氣泡變薄,不能承受原來(上一步)時的水重,破裂狀態(tài)的臨界值變小. 一種極端的情況是沒有塌陷過程,當(dāng)前狀態(tài)超過臨界狀態(tài)值時,氣泡破裂,狀態(tài)值變?yōu)榱? 本文采用氣泡無塌陷破裂進(jìn)行模擬.

2 大規(guī)模氣泡湮滅模擬

2.1 42×50氣泡群 在二維平面上均勻分布42×50個氣泡. 初始時刻,隨機(jī)給出氣泡狀態(tài). 圖4為初始時刻氣泡群的狀態(tài)點(diǎn)陣,臨界狀態(tài)值為16.

圖4 初始時刻氣泡群的狀態(tài)點(diǎn)陣Fig.4 State of bubble group in the initial moment

選擇(i,j)=(20,24)為啟動點(diǎn),即在該點(diǎn)是增加水重的氣泡,在模擬過程中,該點(diǎn)始終在增加狀態(tài)值. 當(dāng)該點(diǎn)加入的狀態(tài)值為12時,有1個氣泡破裂,如圖5所示. 當(dāng)輸入狀態(tài)值為56時,氣泡破裂10個,如圖6所示. 當(dāng)輸入狀態(tài)值為62時,氣泡破裂的圖形形狀有改變,但是破裂數(shù)目仍為10個,如圖7所示. 繼續(xù)加大輸入狀態(tài)值,當(dāng)輸入狀態(tài)值增加到86時,氣泡破裂17個. 破裂情況如圖8所示.

圖5 氣泡在(i,j)=(20,24)處破裂快照Fig.5 Snapshot of the bubbles burst at location(i,j)=(20,24)

圖6 輸入狀態(tài)值為56時氣泡破裂快照Fig.6 Snapshot of bubble burst at input status value to 56

圖7 輸入狀態(tài)值為62時氣泡破裂快照Fig.7 Snapshot of bubble burst at input status value to 62

圖8 輸入狀態(tài)值為86時氣泡破裂快照Fig.8 Snapshot of bubble burst at input status value to 86

繼續(xù)增加輸入狀態(tài)值,當(dāng)輸入狀態(tài)值超過133時,氣泡破裂的范圍不再發(fā)生變化,此時整個結(jié)構(gòu)已處于自組織臨界狀態(tài),當(dāng)輸入狀態(tài)值增加1時,所有氣泡全部破裂,出現(xiàn)氣泡群的湮滅現(xiàn)象. 圖9給出了氣泡破裂個數(shù)與輸入狀態(tài)值間的關(guān)系,由于在演化過程中,每一步增加一個狀態(tài)值,因此該狀態(tài)數(shù)實(shí)際上相當(dāng)于時間步.

2.2 292×232氣泡群 模擬292×232個氣泡群的破裂情況,圖10給出了氣泡破裂個數(shù)與輸入狀態(tài)值間的關(guān)系. 模擬結(jié)果表明,該系統(tǒng)氣泡破裂的個數(shù)只有如下6種情況: 4,5,7,19,21,83,稱為系統(tǒng)的狀態(tài)級,記K=1,2,…,6. 為了考察可能出現(xiàn)的Zipf定律,將上述狀態(tài)級轉(zhuǎn)化成Richter級:L=7-K,并給出它們之間的對數(shù)關(guān)系,如圖11所示,可得一個大規(guī)模氣泡湮滅的Zipf律: log(N)=Alog(L)+C,其中:A=-1.69 4;C=1.92. 結(jié)果表明,大規(guī)模氣泡湮滅是自組織臨界的.

圖9 42×50氣泡群氣泡破裂個數(shù)與輸入狀態(tài)值間的關(guān)系Fig.9 Relationship between the 42×50 input state value and the number of bubbles burst

圖10 292×232氣泡群氣泡破裂個數(shù)與輸入狀態(tài)值間的關(guān)系Fig.10 Relationship between the 292×232 input state value and the number of bubbles burst

圖11 氣泡群Richter級L與相應(yīng)氣泡破裂數(shù)N間的對數(shù)關(guān)系Fig.11 Corresponding logarithmic relationship between bubble burst number N and magnitude on the Richter scale L

綜上可見,本文模型采用了最簡單假設(shè),僅給出了大規(guī)模氣泡的湮滅現(xiàn)象. 真實(shí)的氣泡是三維的、 大小不均勻分布的,并且具有塌陷過程. 模擬真實(shí)氣泡湮滅過程需要結(jié)合實(shí)驗(yàn),反復(fù)矯正模型的狀態(tài)級別和臨界值,同時要考慮三維非均勻分布.

[1]Bak P,TANG Chao,Wiesenfeld K,Self-organized Criticality: An Explanation of 1/fNoise [J]. Phys Rev Lett,1987,59: 381-384.

[2]Bak P,TANG Chao,Wiesenfeld K. Self-organized Criticality [J]. Phys Rev A,1988,38(3): 364-374.

[3]Bak P. How Nature Works: The Science of Self-organized Criticality [M]. New York: Springer-Verlag,1996.

[4]Seck-Tohu-Mora J C,Medina-Marino J,Martínez G J,et al. Emergence of Density Dynamics by Surface Interpolation in Elementary Cellular Automata [J]. Communications in Nonlinear Science and Numerical Simulation,2014,19(4): 941-966.

[5]Yan G W,Zhang J Y,Wang H M,et al. Simple Stochastic Lattice Gas Automaton Model for Formation of River Networks [J]. Phys Rev E,2008,78(6): 066102.

[6]閆廣武. 元胞自動機(jī)與人工生命研究進(jìn)展 [J]. 吉林大學(xué)學(xué)報: 理學(xué)版,2003,41(1): 40-44. (YAN Guangwu. Advances of the Studies on the Cellular Automata and Artificial Life [J]. Journal of Jilin University: Science Edition,2003,41(1): 40-44.)

[7]Ye F,Chen Q W,Li R N. Modelling the Riparian Vegetation Evolution Due to Flow Regulation of Lijiang River by Unstructured Cellular Automata [J]. Ecological Informatics,2010,5: 108-114.

[8]Ferreri L,Venturino E. Cellular Automata for Contact Ecoepidemic Processes in Predator-Prey Systems [J]. Ecological Complexity,2013,13: 8-20.

[9]XUN Jing,NING Bin,LI Keping,et al. The Impact of End-to-End Communication Delay on Railway Traffic Flow Using Cellular Automata Model [J]. Transportation Research Part C: Emerging Technologies,2013,35: 127-140.

(責(zé)任編輯： 趙立芹)

CellularAutomataSimulationforBubblesAnnihilationinaLargeNumbers

HU Yajun,YAN Guangwu
(CollegeofMathematics,JilinUniversity,Changchun130012,China)

With the help of a proposed cellular automata model for the two-dimensional bubble annihilation and a honeycomb structure instead of the two-dimensional bubbles,we simulated the process of the bubble collapse at a load beyond the critical load. This result shows the process meets Zipf’s law.

cellular automata; two-dimensional bubble; Zipf’s law

2014-02-25.

胡亞軍(1990—),男,漢族,碩士研究生,從事流體力學(xué)與復(fù)雜系統(tǒng)的研究,E-mail： h76101@gmail.com. 通信作者: 閆廣武(1964—),男,漢族,博士,教授,博士生導(dǎo)師,從事流體力學(xué)與復(fù)雜系統(tǒng)的研究,E-mail: yangw_jlu@126.com.

國家自然科學(xué)基金(批準(zhǔn)號: 11272133).

O231.5

A

1671-5489(2014)05-0975-04