天文學(xué)進(jìn)展太陽活動起源研究():太陽發(fā)電機理論林元章(中國科學(xué)院國家天文臺北京對試圖解釋太陽活動起源的太陽發(fā)電機理論作了綜合評述。著重介紹了平均場運動學(xué)發(fā)電機理論,包括平均場的a效應(yīng)、運動學(xué)的發(fā)電機和遷移發(fā)電機。討論了MHD發(fā)電機和其它類型發(fā)電機的研究概況。 1引言由文(I)所述可見ll,Babcock和Leighton的太陽活動周模型分別為經(jīng)驗和半經(jīng)驗?zāi)P,而非?yán)格的動力學(xué)描述。目前認(rèn)為維持周期性太陽活動過程的物理機制是太陽等離子體自身運動感應(yīng)的磁場所表現(xiàn)的周期性現(xiàn)象。運動導(dǎo)體通過感應(yīng)能夠產(chǎn)生磁場正是自激發(fā)電機原理。 其基本過程為運動導(dǎo)體切割磁力線產(chǎn)生電場,電場將產(chǎn)生電流(Ohm定律),而電流又產(chǎn)生磁場(Ampere定律),磁場將產(chǎn)生電場(Faraday定律)和Lorentz力,后者將限制導(dǎo)體運動,使磁場達(dá)到某種平衡,形成完整的發(fā)電機循環(huán)。因此人們想到太陽磁周規(guī)律可能也是太陽本身自然形成的自激發(fā)電機產(chǎn)生的。人造自激發(fā)電機中有許多導(dǎo)體、絕緣體和導(dǎo)線以及精巧的聯(lián)結(jié)方式,其運轉(zhuǎn)過程不難理解。然而太陽上是連成一片的導(dǎo)電等離子體,并且永遠(yuǎn)處在短路狀態(tài),它能形成自激發(fā)電機嗎,換句話說,也就是感應(yīng)方程是否存在發(fā)電機解,(上式Bv和r分別為磁場、速度和磁擴散系數(shù)。)許多理論研究者進(jìn)行了艱難的探索。從觀測上以及Babcock和Leighton的模型中可見,太陽活動周的主要特征是極向磁場與環(huán)向磁場之間的相互轉(zhuǎn)換W.因此探索的主要思路就在現(xiàn)有太陽速度場(主要是較差自轉(zhuǎn)和對流)條件下,尋求可以實現(xiàn)極向磁場Sp與環(huán)向磁場St之間不斷相互轉(zhuǎn)換的磁場解。 在初期的探索過程中,曾經(jīng)提出過兩個重要的制約性定理(有時稱作反發(fā)電機定理)。其一是CowlingM發(fā)現(xiàn)發(fā)電機解必須是非軸對稱的;另一是Bullard和GellmanM發(fā)現(xiàn),若速度場只有太陽自轉(zhuǎn),不管自轉(zhuǎn)是否均勻,則只能由極向磁場產(chǎn)生環(huán)向磁場Bt,不能從轉(zhuǎn)化為Sp,從而無法完成Bp與Bt之間的相互轉(zhuǎn)換,即不存在發(fā)電機解。Parker首先提出,如果采用平均場的概念,則可繞開上述制約,找到軸對稱的平均場發(fā)電機解。Parker指出,由于Coriolis力的作用,太陽對流層中的上升和下降氣團將產(chǎn)生轉(zhuǎn)動,大量小尺度轉(zhuǎn)動的平均結(jié)果等效于存在大尺度電場(稱為效應(yīng)),從而可以導(dǎo)致風(fēng)向Sp轉(zhuǎn)化。隨后Steenbeck等人問、Steenbeck和Krause、Moffatt以及Krause和Radler間等人的深入探討,逐步完善了平均場的發(fā)電機理論。 2平均場運動學(xué)發(fā)電機理論平均場和a效應(yīng)觀測表明,太陽表面磁場存在一定的起伏范圍,因此可以把磁場寫成其中〈B〉為平均場,可以把它理解為大量隨機變化的現(xiàn)實磁場集合體的平均值;6為起伏部分,其平均值〈b〉=.同時速度場也可以寫成其中平均速度〈蚴可以理解為太陽全球尺度運動(尤其是較差自轉(zhuǎn))的平均速度,u表示不規(guī)則的湍動對流速度。通常無需知道V的細(xì)節(jié),只需知道其平均值〈V〉和U的統(tǒng)計性質(zhì)。在以下的討論中,將認(rèn)為V已給定,它不受B的影響。換句話說,將作為運動學(xué)方式,而非動力學(xué)方式處理。實際上Babcock和Leighton的太陽周模型也是運動學(xué)模型。 把(2)和(3)代入(1)式,并把平均場和起伏場分離,得到其中U為擴散系數(shù),而這里的平均電場五是關(guān)鍵參量。若E已知,就可以由(4)式求解〈B〉。原則上說,必須先從式得到一個b解,再代入(6)計算以(B〉表示的E.這實際上很困難。因此通常采用所謂“一級平滑近似”,即假定二級微量項G可以忽略,于是可由(5)式求得b,再由(6)式求得五?梢宰C明,當(dāng)滿足如下兩個條件之一時,一級平滑近似成立或其中u為u的特征值,/和T分別為u和變化的空間和時間尺度。條件(8)可保證G遠(yuǎn)小于(5)式右端最后一項,這相當(dāng)于磁Reynold數(shù)很小的情況;而若(9)式成立,則G與碧相比可忽略。遺憾的是在太陽上這兩個條件均不滿足。太陽的磁Reynold數(shù)很大,而且觀測表明盡管如此,我們?nèi)钥蛇M(jìn)行以下啟發(fā)性的討論。 若僅限于了解定性結(jié)果,則除了忽略G之外,再略去(5)的〈《>項。同時,考慮到太陽等離子體的高電導(dǎo)率,再略去(5)的擴散項,立即可得到b的解為辦二上式系對時間積分,但無需寫出起始場u(-oo),因為它肯定與u⑴無關(guān),從而對計算*無貢獻(xiàn)。若u為弱各向同性湍流場,則可以證明*可以寫為tml把(11)式代入(4)式,就得到平均場的感應(yīng)方程為把(14)式與(1)式比較,可發(fā)現(xiàn)有兩個變化。其一是擴散系數(shù)增加了一項/3,其大小為在太陽對流層中,吼《= 108l9m2S*.另一更重要的變化>即方程中增加了一項a〈B〉。正是這一項保證了平均場〈S>不受上述兩個反發(fā)電機定理的制約。通常把平均場感應(yīng)方程(14)中增加的這一關(guān)鍵項稱為a效應(yīng)。從(12)式可見,a實際上代表湍動對流場u的螺度(felicity)。太陽對流層中的a值很難估計。不過從(12)式可見,a 0要求流場必須有螺度。Kraust已證明M,太陽對流層中的密度分層與Coriolis力相結(jié)合,能夠使流場存在螺度。在太陽北半球,對流層中不斷上升膨脹的氣團由于Coriolis力作用將獲得左旋螺度,下沉氣團也將獲得左旋螺度。這樣,北半球左旋螺度占優(yōu)勢,故螺度為負(fù)值。這一點已經(jīng)得到觀測證實。而由(12)式知《符號與螺度相反,故太陽北半球a為正值。Krause等人的計算表明其中p為太陽平均角速度,其取正號和負(fù)號分別對應(yīng)于北半球和南半球。在太陽對流層中,變化的空間尺度Z的范圍很大,因而a的大小約從每秒幾cm至100m.平均場感應(yīng)方程(14)是大多數(shù)太陽發(fā)電機理論研究者的工作基礎(chǔ),并且也被研究地球磁場發(fā)電機和恒星磁場發(fā)電機理論的學(xué)者采用,甚至被應(yīng)用于研究其它天體如吸積盤和星系磁場。但是它本身包含著許多假定,其中最嚴(yán)重的就是一級平滑近似。這一近似所要求的條件(8)和(9)式在太陽上至多只能達(dá)到臨界滿足。 運動學(xué)的afi發(fā)電機由上述可見,a效應(yīng)的作用是保證平均場的存在,而a效應(yīng)與太陽自轉(zhuǎn)速度場D相結(jié)合,就可構(gòu)成aD發(fā)電機。若仍然采用運動學(xué)的處理方法,就是運動學(xué)的發(fā)電機。假定在球坐標(biāo)系(rAp)中,取給定形式的a(r)和太陽的自轉(zhuǎn)角速度分布同時根據(jù)(17)式把a取為對太陽赤道反對稱,即自轉(zhuǎn)角速度則取為對赤道對稱,即若除了自轉(zhuǎn),無其它運動,于是把平均場〈B〉分解為極向場和環(huán)向場,即其中。Sp(22)的A(r,0,i)為的矢量勢,它由和軸對稱確定。于是平均場的感應(yīng)方程也分解為極向和環(huán)向分量,對于=常數(shù)的最簡單情況,它們是W其中由(24)和(25)式可以看到a效應(yīng)的關(guān)鍵作用。當(dāng)a =0時,(24)退化為擴散方程,極向場立即指數(shù)衰減,4消失。這時(25)也變成擴散方程,B也指數(shù)衰減后消失。由(26)還可清楚地看出,環(huán)向場是原始極向場通過較差自轉(zhuǎn)和a效應(yīng)產(chǎn)生的。不過對于太陽,后者貢獻(xiàn)比前者小,故常把(25)含的項略去,其條件是其中及為太陽半徑。雖然a的大小不太清楚,但(27)可能是成立的。這樣,運動學(xué)a/發(fā)電機的原理可概括為:通過a效應(yīng)產(chǎn)生5tk2必須為負(fù)值,于是(36)和(37)中的正號可以取消。另一個解即與Im(W)=0.對應(yīng)的臨界穩(wěn)定解(38)式表明乘積炻必須超過某一數(shù)值,亦即發(fā)電機數(shù)必須超過一臨界值,發(fā)電機才能運轉(zhuǎn)。 對于臨界解,振動頻率為的實部是負(fù)值。所以平均場(34)是沿a:軸正向傳播的波。若起初假定fc<0,只要a% <0,得到的結(jié)果是一樣的。反之,如a)>0,則波將沿負(fù):c軸方向傳播。通過適當(dāng)放置坐標(biāo)系,可以證明發(fā)電機波一般將沿等角速度面遷移。同時,工作在對流層的發(fā)電機還存在磁流損失問題。因而一些研究者提出了工作于對流層底部過沖層中的過沖層發(fā)電機模型(overshootlayerdynamo)。對流層底部的過沖層是指按Schwarzchild判據(jù)確定的對流層下邊界下方的一個薄層。該層已不滿足Schwarzchild的對流判據(jù)(|f|ad< |J,然而上方對流層中的氣團運動仍然對這一層沖擊。因此過沖層中仍有對流氣團運動。經(jīng)過過沖層之后,氣團才會徹底瓦解,因此過沖層的厚度約為一個混合長距離,對于太陽估計為leoookmWi.過沖層雖然比對流層穩(wěn)定,但仍有氣團運動,因此存在a效應(yīng),afi發(fā)電機能夠工作。許多研究者已對過沖層發(fā)電機作了探討。這里的界面是指太陽對流層與內(nèi)部非對流區(qū)的分界面,因此亦指對流層底部。不過與過沖層發(fā)電機工作于過沖層不同,界面發(fā)電機的工作源區(qū)實際上有兩個,分別位于界面的上方和下方。界面下方的工作源區(qū)中主要是大尺度剪切運動,其作用是完成極向場轉(zhuǎn)化為環(huán)向場;而環(huán)向場產(chǎn)生極向場的a效應(yīng)則在界面上方完成。隨后Macgregor和Charbonneau等人134,351對界面發(fā)電機模型作了具體計算。他們在平均場運動學(xué)框架下得到的結(jié)果也能解釋太陽活動周中的主要現(xiàn)象。 此外,研究表明,MHD發(fā)電機還具有混沌行為。Weiss等人曾研究一種最簡單的發(fā)電機模型。他們發(fā)現(xiàn)當(dāng)發(fā)電機數(shù)很大和擴散很小時,發(fā)電機的周期解變成不穩(wěn)定。先是被多種周期取代,而最終變成了混純解。這時,除了相當(dāng)于太陽活動周的正常周期外,還包含有低活動的周期,很像太陽活動超長期變化規(guī)律中的世紀(jì)極小期37,38.除了目前已成為主流的發(fā)電機模型外,可能還有依賴其它效應(yīng)的發(fā)電機模型也可以解釋天體的磁場。例如,其工作原理無需對流運動僅依靠等離子體湍動波效應(yīng)和天體自轉(zhuǎn)的等離子體湍動發(fā)電機(pfi發(fā)電機)391.不過在a/2發(fā)電機已能夠基本上解釋太陽活動周主要特征的前提下,目前尚難以判定發(fā)電機效應(yīng)是否會在太陽大尺度磁場的產(chǎn)生和演化中起某種作用,哪怕是次要作用。