理工論文三維可激發(fā)系統(tǒng)中Ｗｉｎｆｒｅｅ湍流態(tài)的控制

1、三維可激發(fā)系統(tǒng)中湍流態(tài)的控制 三維可激發(fā)系統(tǒng)中湍流態(tài)的控制是小柯論文網(wǎng)通過網(wǎng)絡(luò)搜集，并由本站工作人員整理后發(fā)布的，三維可激發(fā)系統(tǒng)中湍流態(tài)的控制是篇質(zhì)量較高的學(xué)術(shù)論文，供本站訪問者學(xué)習(xí)和學(xué)術(shù)交流參考之用，不可用于其他商業(yè)目的，三維可激發(fā)系統(tǒng)中湍流態(tài)的控制的論文版權(quán)歸原作者所有，因網(wǎng)絡(luò)整理，有些文章作者不詳，敬請諒解，如需轉(zhuǎn)摘，請注明出處小柯論文網(wǎng)，如果此論文無法滿足您的論文要求，您可以申請本站幫您代寫論文，以下是正文。 摘 要 三維可激發(fā)系統(tǒng)中的winfree湍流態(tài)被認(rèn)為與心室纖維性顫動可能有著密切的關(guān)系.如何快速而有效地消除winfree湍流是理論工作者非常關(guān)心的問題.文章作者發(fā)展了兩種可有效

2、地控制三維可激發(fā)系統(tǒng)中winfree湍流態(tài)的方法:(1)局域的周期刺激,即通過小區(qū)域內(nèi)的周期刺激所產(chǎn)生的激發(fā)波來消除winfree湍流;(2)全局的時空微擾,通過對參數(shù)的時空調(diào)制來消除wnifree湍流.關(guān)鍵詞 可激發(fā)系統(tǒng),winfree湍流態(tài),心顫,控制controlling winfree turbulence in excitable mediazhang hong1 cao zhou-jian2 li bing-wei1 wu ning-jie1ying he-ping1 hu gang2,(1 zhejiang institute of modern physics and depa

3、rtment of physics, zhejiang university, hangzhou 310027, china)(2 department of physics, beijing normal university, beijing 100875, china)abstract it is believed that ventricular fibrillation is related to winfree turbulence in three-dimensional excitable media and much attention has been directed t

4、owards how to eliminate effectively this kind of spatiotemporal chaos. we present two schemes to suppress winfree turbulence: (1) localized stimulations; these generate control waves which can successfully suppress existing turbulence and dominate the whole system; (2) spatiotemporal modulation (tra

5、veling modulation) of the excitability; in this case, winfree turbulence can also be eliminated by choosing suitable control parameters.keywords excitable media, winfree turbulence, ventricular fibrillation, control1 引言可激發(fā)介質(zhì)中非線性波的傳播是當(dāng)前非線性動力學(xué)的一個重要的研究領(lǐng)域,螺旋波則是該領(lǐng)域中一個研究熱點(diǎn)1.螺旋波廣泛地存在于各類激發(fā)介質(zhì),包括心臟組織2一氧化碳在鉑金表

6、面的氧化反應(yīng)3以及著名的belousov-zhabotinsky(簡稱bz)化學(xué)反應(yīng)之中46.第一個在bz化學(xué)反應(yīng)中觀測到會旋轉(zhuǎn)的螺旋波的是美國著名生物化學(xué)家winfree4.在過去的30多年中,對螺旋波的基本動力學(xué)以及如何控制螺旋波取得了相當(dāng)大的進(jìn)展710.三維螺旋波又稱回卷波,可以認(rèn)為是二維螺旋波在三維空間的延伸.對應(yīng)于螺旋波的相奇點(diǎn)(又稱拓?fù)淙毕?,回卷波的中心是一條相奇異線,波繞著這條奇異線旋轉(zhuǎn).如果沿著奇異線的垂直方向做一個截面,這個截面便是螺旋波.值得指出的是,回卷波的奇異線并不是總是成一條直線的,它還可以彎曲形成一個環(huán),甚至本身還有非常復(fù)雜的拓?fù)浣Y(jié)構(gòu)11.因此從這個意義上講,無論

7、是回卷波的幾何結(jié)構(gòu)和動力學(xué)性質(zhì)都要比二維螺旋波復(fù)雜得多.winfree在bz反應(yīng)中發(fā)現(xiàn)螺旋波一年以后,又在實(shí)驗(yàn)中發(fā)現(xiàn)了回卷波12.更為重要的是,他還把心臟節(jié)律不齊,特別是心臟纖維性顫動,與此相聯(lián)系起來13.據(jù)美國心臟協(xié)會估計(jì),在美國,每年死于心臟病的人數(shù)大概為480,400,而由于心室纖維性顫動導(dǎo)致死亡的約為335,00014.到目前為止,科學(xué)家為了解釋這種心室纖維性顫動,提出了各種各樣的機(jī)制15.其中有一種觀點(diǎn)認(rèn)為,發(fā)生心臟纖維性顫動很可能與三維回卷波失穩(wěn)有著密切的關(guān)系13, 1618.在三維可激發(fā)系統(tǒng)中,由于奇異線負(fù)張力引起的回卷波失穩(wěn)而導(dǎo)致的時空混沌稱之為winfree湍流.目前,在臨床

8、上使用比較廣泛的除顫方法是利用高壓電(5kv ).當(dāng)發(fā)生心顫時,給心臟施加一個強(qiáng)電擊,瞬間全面消除心臟中的電活動,從而使竇房結(jié)信號可以重新控制心臟系統(tǒng),這樣心臟就可以恢復(fù)到正常起搏的狀態(tài)19.這種方法雖然說比較有效,但也存在著一些不足之處.比如,強(qiáng)電擊會損害一些正常的心肌細(xì)胞,甚至造成死亡20.此外,這種方法會給患者帶來極大的痛苦.因此,如何利用弱電擊(比如幾毫伏的電壓)來控制甚至是消除心顫是臨床醫(yī)學(xué)工作者所期望的,這也促使了理論工作者提出各種有效的方法來控制三維時空湍流.本文介紹了三維可激發(fā)系統(tǒng)中winfree湍流態(tài)控制的最新進(jìn)展.我們在alonso等人的工作基礎(chǔ)上18,發(fā)展了兩種有效控制w

9、infree湍流態(tài)的方法:(1)局域的周期刺激:通過在一個非常小的區(qū)域內(nèi)施加較弱的周期信號,產(chǎn)生新的激發(fā)波來消除已存在的三維時空湍流21;(2)全局的時空微擾:通過對系統(tǒng)參數(shù)的時空調(diào)制(行波調(diào)制)來消除winfree湍流態(tài)22.這些結(jié)果對于臨床上用低幅控制信號來消除心室纖維性顫動將有一定的啟發(fā)意義.2 全空間周期調(diào)制alonso等人在science周刊上報(bào)道了通過對參數(shù)的全空間的周期調(diào)制,消除三維可激發(fā)系統(tǒng)中的winfree湍流態(tài)18.他們采用的數(shù)值模型為典型的可激發(fā)系統(tǒng)23:ut=1u(1-u)(u-+ba)+d2u,(1)t=u-,(2)圖1 (a)三維可激發(fā)模型(1-2)中的winfre

10、e湍流態(tài),曲線代表回卷波的奇異線;(b)至(d)為施加局域周期刺激后的各種漸近狀態(tài),其中bf=0.75,f=1.6;(b)小立方體(666格點(diǎn))內(nèi)施加周期刺激;(c)條狀區(qū)域內(nèi)(16150格點(diǎn))內(nèi)施加周期刺激;(d)平面上(1150150格點(diǎn))施加周期刺激;(e)奇異線長度隨時間的變化:曲線a為無任何周期刺激,曲線b為小立方體內(nèi)的周期刺激(對應(yīng)于(b)的情況),曲線c為條狀區(qū)域內(nèi)周期刺激(對應(yīng)(c)的情況),曲線d為平面上的局域刺激(對應(yīng)(d)的情況).從圖(a)至(e)中,系統(tǒng)尺寸大小為606060 (150150150格點(diǎn)),圖(b),(c),(d)的初始條件為(a) 這里的參數(shù)代表快變量

11、u(相應(yīng)于心臟細(xì)胞的膜電壓)和慢變量(相應(yīng)于跨膜電流)的時間尺度之比,是一個小量.對于該模型,盡管在二維系統(tǒng)中不存在螺旋波湍流,但在三維系統(tǒng)中,存在由于奇異線負(fù)張力失穩(wěn)而引起回卷波向winfree湍流態(tài)的轉(zhuǎn)變.為了抑制winfree湍流態(tài),他們采取的控制方法是對參數(shù)b進(jìn)行全局周期調(diào)制,即b=b0+bfcos(ft).(3)研究結(jié)果表明,當(dāng)f0時(0是在同樣參數(shù)下,二維系統(tǒng)中螺旋波的旋轉(zhuǎn)頻率),winfree湍流態(tài)最終會被消除.在此基礎(chǔ)上,我們發(fā)展了兩種較為有效的方法來控制winfree湍流態(tài),這包括局部的周期刺激和全局的時空微擾.下面予以詳細(xì)的介紹.3 局部周期刺激圖2 在周期刺激區(qū)域?yàn)樾×⒎?/p>

12、體的情況下,消除湍流態(tài)的速度r對各種外加參數(shù)的依賴關(guān)系 (a)n=6,bf=0.75 ,r對f 依賴關(guān)系;(b)bf=0.75 ,f =1.6,r對n依賴關(guān)系;(c)f=1.6,n=6,r對bf依賴關(guān)系;(d)對于f=1.6的情況,在相平面bf-n下的可控制區(qū)域.在數(shù)值模擬中,可控是指奇異線長度在t500時縮短到0的事實(shí)我們采用的仍舊是由(1)式和(2)式描述的可激發(fā)系統(tǒng),也同樣采取b的周期調(diào)制:b=b0+bfcos(ft).但有一點(diǎn)不同的是控制區(qū)域是集中在一個小的范圍內(nèi)而并不是在全空間(即整個系統(tǒng)).圖1展示了通過各種不同形式的局部刺激來控制(消除)winfree湍流態(tài).在圖1(b)中,在一

13、個小立方體上施加周期信號;在圖1(c)中,在某一個面上的條帶上施加信號;而在圖1(d)中,則是在整個平面上加信號.從圖1中不難發(fā)現(xiàn),局域的刺激將誘導(dǎo)各種新波的產(chǎn)生,其中包括球面波(圖1(b)柱面波(圖1(c)以及平面波(圖1(d).這種新產(chǎn)生的波與原來存在的湍流波相互競爭,最終將其成功地消除掉.值得指出的是,我們所采取的局域周期信號的強(qiáng)度與變量u和屬于同一數(shù)量級,即對winfree湍流態(tài)實(shí)現(xiàn)了低強(qiáng)度控制,對應(yīng)于對心臟的弱信號除顫.與全局控制的結(jié)果相比18,局部控制方法具有下面的特點(diǎn):(1)全局的方法需要將信號(或刺激)施加到整個系統(tǒng),這種方法在實(shí)際情況當(dāng)中會遇到很多不便(比如,將電脈沖信號施加

14、在心臟內(nèi)部的整個心肌組織上是非 三維可激發(fā)系統(tǒng)中湍流態(tài)的控制是小柯論文網(wǎng)通過網(wǎng)絡(luò)搜集，并由本站工作人員整理后發(fā)布的，三維可激發(fā)系統(tǒng)中湍流態(tài)的控制是篇質(zhì)量較高的學(xué)術(shù)論文，供本站訪問者學(xué)習(xí)和學(xué)術(shù)交流參考之用，不可用于其他商業(yè)目的，三維可激發(fā)系統(tǒng)中湍流態(tài)的控制的論文版權(quán)歸原作者所有，因網(wǎng)絡(luò)整理，有些文章作者不詳，敬請諒解，如需轉(zhuǎn)摘，請注明出處小柯論文網(wǎng)，如果此論文無法滿足您的論文要求，您可以申請本站幫您代寫論文，以下是正文。常困難的);相反,將信號施加在系統(tǒng)的表面部分區(qū)域是比較方便的;(2)用局域周期刺激的方法來消除winfree湍流態(tài)的速度很快(在可比的輸入信號總功率下,遠(yuǎn)遠(yuǎn)快于全局控制方法).圖

15、1(e)給出了在無任何刺激下,小立方體條狀以及整個平面刺激下的奇異線長度隨時間的變化曲線.我們發(fā)現(xiàn),局域刺激可以使奇異線的長度在2030個螺旋波周期內(nèi)減小至0(即消除),而在文獻(xiàn)18中,全局刺激則需要300個螺旋波周期左右.對于局域刺激來控制winfree湍流態(tài),效率是一個重要因素,即用盡可能小的功率來達(dá)到消除winree湍流態(tài)的目的.假設(shè)注入信號的總功率正比于mb2f,其中m是注入?yún)^(qū)域的格點(diǎn)數(shù).數(shù)值模擬發(fā)現(xiàn),注入的區(qū)域如果是小立方體形狀的,那么控制winfree湍流態(tài)的效率是最高的.盡管由于曲率的影響,平面波峰的運(yùn)動速度要比凸的波峰(比如說球面波)來的大,但是要通過局部刺激首先必須要產(chǎn)生一個

16、完整的波峰.產(chǎn)生一個大的完整的平面波顯然要比產(chǎn)生一個小的球面波要困難,因此小立方體刺激要優(yōu)于其他兩種局部控制的方法.在接下來的討論當(dāng)中,我們著重關(guān)注局部刺激是小立方體的情況.圖2給出了在立方體局部刺激下,控制速度對于外加頻率外加振幅以及刺激區(qū)域大小的依賴關(guān)系.從這些圖中不難發(fā)現(xiàn):(1)要消除winfree湍流態(tài),對外加頻率要有一定的選擇:頻率太小或太大都不能消除winfree湍流態(tài),并且外加頻率要稍大于同樣參數(shù)下二維螺旋波的特征頻率,即f0.這是因?yàn)?只有當(dāng)局部刺激產(chǎn)生的球面波的頻率大于回卷波的頻率時,球面波才能占優(yōu)勢,從而消除湍流態(tài)(頻率大的占優(yōu));另一方面,外加頻率應(yīng)有一個上限,在原有wi

17、nfree 湍流態(tài)存在的情況下,要有效地產(chǎn)生球面波,外加頻率不應(yīng)該遠(yuǎn)大于0.(2)注入信號的小立方體區(qū)域必須包括足夠多的格點(diǎn)數(shù)(但與整個系統(tǒng)相比,還是非常小的).這樣產(chǎn)生的波峰曲率才能小于波傳播的臨界曲率,球面靶波才能夠傳播,從而控制winfree湍流態(tài).此外,信號的強(qiáng)度必須大于某一個臨界值(閾值).比如在本系統(tǒng)的數(shù)值模擬中,對于n4(nnn格點(diǎn))或bf0.53時小立方體方法不再奏效.這是因?yàn)橐a(chǎn)生球面靶波,局部刺激的振幅必須要大于一個閾值.從圖2(b)和2(c)還可以發(fā)現(xiàn),一旦n和bf大于臨界值,消除湍流態(tài)的速度幾乎不依賴于n和bf. 換句話說,如果繼續(xù)增加n和bf,控制效率也不會明顯上升.

18、這就意味著在一個優(yōu)化的參數(shù)組合情況下,可以達(dá)到高效率的控制.以上結(jié)果為臨床中用低幅電信號快速除顫提供了有益的參考.為了更加清晰地表明如何通過局部刺激(小立方體)來消除winfree 湍流態(tài)的,在圖3中,我們給出了其時間演化圖.從圖3可以看到,一個在小立方體區(qū)域內(nèi)的刺激,激發(fā)出一個球面靶波.該球面靶波與湍流波相互作用,最終將回卷波的奇異線排擠到系統(tǒng)的邊界外,整個系統(tǒng)最終被球面波所占據(jù).這時如果撤離局部周期信號,整個系統(tǒng)將恢復(fù)到靜息狀態(tài).圖3 周期刺激區(qū)域?yàn)樾×⒎襟w的情況下,消除winfree湍流態(tài)的具體演化圖(相當(dāng)于圖1中的(b).從圖中可以發(fā)現(xiàn),一個局域的周期刺激將激發(fā)出一個球面靶波.借助于靶

19、波的傳播,最終將winfree湍流態(tài)趕出整個系統(tǒng)4 全局時空微擾前面講的控制方法18, 21以及以前許多控制二維系統(tǒng)的螺旋波湍流態(tài)方法都集中在對參數(shù)的純時間的周期調(diào)制.最近的研究表明,在二維激發(fā)系統(tǒng)中外加了一個較強(qiáng)的平面波控制信號,系統(tǒng)最終被激發(fā)出平面波,同時消除了系統(tǒng)中的螺旋波及湍流態(tài)24.然而在三維的情況下,通過對參數(shù)的時空調(diào)制來抑制winfree湍流態(tài),還沒人研究過.以下我們研究如何通過對參數(shù)弱的時空調(diào)制-行波調(diào)制來消除winfree湍流態(tài).對于時空調(diào)制,可以把參數(shù)b寫成下面的形式:b=b0+bfcos(kfr-ft).(4)純時間的調(diào)制可以認(rèn)為是行波調(diào)制在kf=0 的一種特殊情況.在接

20、下來的討論中,都是針對kf0,f0 的情況.下面主要討論空間共振的情況,即kf=k0 (k0=2/0,0是在同樣參數(shù)下,二維系統(tǒng)中螺旋波的特征波長).圖4展示了在參數(shù)bf=0.03,kf=k0的情況下,通過對b進(jìn)行行波調(diào)制來消除winfree湍流態(tài)的一個例子.與空間均勻周期調(diào)制(即kf=0)相比18,行波調(diào)制對于控制湍流態(tài),效率更高:圖5比較了參數(shù)做純時間和做時空調(diào)制下,奇異線長度隨時間的改變.我們發(fā)現(xiàn),在同樣的組合參數(shù)(bf,f)=(0.03,1.2)情況下,行波調(diào)制使奇異線長度收縮到0的速度要比在全局周期調(diào)制快12倍.同時也可以看到,對于bf=0.03和kf=k0 ,存在一個最優(yōu)頻率,在這

21、種情況下,通過行波調(diào)制來抑制winfree湍流態(tài)的速度比均勻周期調(diào)制快將近36倍. 圖4 通過時空調(diào)制消除winfree湍流態(tài) (a)和(b)表示t=0時,初始時回卷環(huán)以及垂直奇異線截面上的二維螺旋波斑圖;(c)t=224,奇異線由于負(fù)張力失穩(wěn)而產(chǎn)生的winfree湍流態(tài)(類似圖1(a).此時,施加的時空調(diào)制信號bf=0.03,kf=k0, f=1.2;(d)t=400,winfree湍流態(tài)最終被消除圖5 在兩種不同的全局周期刺激下,奇異線長度l隨時間的變化 (a)全局時空調(diào)制:b=b0+bfcos(kfr-ft) ,bf=0.03 ,kf=k0 ,f=1.2 ;(b)全局均勻調(diào)制:b=b0+

22、bfcos(ft),bf=0.03,f=1.2 ;(c)全局時空調(diào)制:f=1.0,其他參數(shù)同(a)為了更好的理解行波調(diào)制對winfree湍流態(tài)的高效控制,我們研究了行波調(diào)制對回卷環(huán)(圖4(a)的影響,特別考慮了角(是指kf 和z方向的夾角)對回卷環(huán)的抑制作用.圖6展示了抑制速度r對的依賴關(guān)系:當(dāng)=0,即kf平行于z軸的時候,控制速度達(dá)到最大;當(dāng)=/2,即kf垂直于z軸的時候,控制速度最小,但是不管取何值,行波調(diào)制的控制速度總是要大于空間均勻調(diào)制的速度.winfree湍流態(tài)可以近似地看成由許許多多不同取向的部分回卷環(huán)構(gòu)成的,因此,利用行波調(diào)制的方法來控制winfree湍流態(tài)要比空間均勻調(diào)制來得更

23、加有效.圖6 回卷環(huán)的收縮速度r與角(是指kf 和z方向的夾角)的依賴關(guān)系(b=b0+bfcos(kfr-ft),bf=0.03,kf=k0,f=1.2.直線表示在全空間均勻調(diào)制bf=0.03,f=1.2 的情況下的收縮速度.從圖中對比可以發(fā)現(xiàn),對于任何取向(角)在行波調(diào)制下,回卷環(huán)的收縮速度要快于全空間均勻調(diào)制的情況)5 總結(jié)和展望我們介紹了兩種有效的控制三維可激發(fā)系統(tǒng)中winfree湍流態(tài)的方法.需要強(qiáng)調(diào)的是,無論是通過局部的周期刺激,還是全局的時空微擾,所施加的信號強(qiáng)度均是比較微弱的.這對臨床上低幅除顫有著潛在的應(yīng)用價(jià)值.但要指出的是,本文采用的模型過于簡單,與真實(shí)的心臟模型有很大距離.

24、比如,我們沒有考慮心臟復(fù)雜的幾何結(jié)構(gòu),力電耦合等一些基本屬性,而這些屬性對心電信號的傳播有著重要的影響.因此,接下來要做的是如何將本文的方法直接應(yīng)用于更加實(shí)際的心臟模型.參考文獻(xiàn)1 cross m c, hohenberg p c. rev. mod. phys., 1993, 65:8512 davidenko j m, pertsov a v, salomonsz r et al. nature, 1992,355:3493 jakubith s, rotermund h h, engel w et al. phys. rev. lett., 1990, 65:30134 winfree

25、a t. science, 1972, 175: 6345 vanag v k, epstein i r. science, 2001, 294: 8356 ouyang q, swinney h l, li g. phys. rev. lett., 2000, 84: 10477 tyson j j, keener j p. physica d, 1988, 32: 3278 mikhailov a s, davydov v a, zykov v s. physica d, 1994,70: 19 meron e. phys. rep., 1992, 218: 110 mikhailov a

26、 s, showalter k. phys. rep., 2006, 425: 7911 winfree a t, strogatz s h. physica d, 1983,9: 33 三維可激發(fā)系統(tǒng)中湍流態(tài)的控制是小柯論文網(wǎng)通過網(wǎng)絡(luò)搜集，并由本站工作人員整理后發(fā)布的，三維可激發(fā)系統(tǒng)中湍流態(tài)的控制是篇質(zhì)量較高的學(xué)術(shù)論文，供本站訪問者學(xué)習(xí)和學(xué)術(shù)交流參考之用，不可用于其他商業(yè)目的，三維可激發(fā)系統(tǒng)中湍流態(tài)的控制的論文版權(quán)歸原作者所有，因網(wǎng)絡(luò)整理，有些文章作者不詳，敬請諒解，如需轉(zhuǎn)摘，請注明出處小柯論文網(wǎng)，如果此論文無法滿足您的論文要求，您可以申請本站幫您代寫論文，以下是正文。312 winfree

27、 a t. science, 1973, 181: 93713 winfree a t. science, 1994, 266: 100314 american hear association. 2005, heart disease and stroke statistics-2005 update. dallas (tx) american heart association15 myerburg r j et al. in: cardiac electrophysiology, from cell to beside, zipes d p, jalife j eds. (saunder

28、, philadelphia pa, 1990), 66616 gray r a, jalife j, panfilov a v et al. science, 1995, 270 :122217 panfilov a v, hogeweg p. science, 1995, 270: 1224718 alonso a, sagues f, mikhailov a s. science, 2003, 299: 1722; chaos, 2006, 16: 02312419 eisenberg m s et al. sci. am. 1986, 254: 2520 pool r. science

29、, 1990, 247: 129421 zhang h, cao z, wu n j, ying h p, hu g. phys. rev. lett., 2005, 94: 18830122 wu n j, zhang h, ying h p, cao z, hu g. phys. rev. e, 2006, 73: 060901 (r)23 barkley d, kness m, tuckerman l s. phys. rev. a, 1990,42: 248924 wang p y, xie p, yin h w. chinese phys., 2003, 12: 674“本文中所涉及到的圖表、注解、公式等內(nèi)容請以pdf格式閱讀原文”其他參考文獻(xiàn)baker, sheridan. the practical stylist. 6th ed. new york: harper & row, 1985.flesch, rudolf. the art of plain talk. new york: harper & brothers, 1946.gowers, ernest. the complete plain