1、第6章 磁流體力學不穩定性§6.1概論等離子體能夠被磁場約束并處于力學平衡狀態。一個處于力學平衡狀態的等離子體位形，當它受到某種擾動，偏離平衡態時，等離子體將如何反應？是越來越偏離平衡態，最后導致平衡態被破壞呢，還是很快將擾動抑制住回到平衡態前者是不穩定平衡，后者是穩定平衡但當磁流體處在非熱力學平衡態，其內部存在著可以轉換成擾動能量的自由能時，在合適的條件下有些擾動就可能發展成為在大范圍、長時間、能量超過熱噪聲水平的大幅度集體運動這種集體運動就稱為不穩定的模式，相應現象就稱為磁流體的不穩定性研究等離子體的各種不穩定性，闡明其物理機制，探索抑制不穩定性的方法，一直是受控核聚變研究的重要

2、課題磁約束等離子體可以處于力學平衡狀態，但它不是完全的熱力學平衡態等離子體處于非熱力學平衡狀態意味著等離子體具有較高的自由能，因而必然會產生從較高能量狀態過渡到較低能量狀態的宏觀或微觀運動等離子體偏離熱力學平衡態大體有兩類方式一類是等離子體宏觀參數如密度、溫度、壓強或其它熱力學量的空間局域性和不均勻性；另一類是等離子體的速度空間分布函數偏離麥克斯韋分布由于前一種原因產生不穩定性時，等離子體通常以整體形式在空間改變其形狀，因而稱為宏觀不穩定性。由后一種原因產生的不穩定性稱為微觀不穩定性宏觀不穩定性通常用磁流體力學方程進行分析，因而也稱為磁流體力學不穩定性，而微觀不穩定性則用動力論方程進行分析，因

3、而也叫動力學不穩定性由于磁流體力學不穩定性在磁約束核聚變等離子體中具有更重要的地位，處理方法也相對地比較容易，因此本節僅討論磁流體力學不穩定性下面我們將首先從分析流體的瑞利一泰勒不穩定性（RayleighTaylor instability）入手，這樣做物理圖像清晰，易于理解然后討論在分析磁流體力學不穩定性中得到廣泛應用的能量原理在這基礎上分析幾種主要的宏觀不穩定性，最后討論等離子體電阻對不穩定性的影響 下面是幾種典型的磁流體不穩定模式 例 1瑞利一泰勒（Rayleigh-Taylor）不穩定性（圖41）； 例2開爾文一亥姆霍茲（KelvinHelmholtz）不穩定性（圖42）； 例3臘腸型

4、不穩定性（圖43）； 例4彎曲型不穩定性（圖4.4）； 例5. 磁島（圖4.5）；例6. 磁重聯（圖46）每種不穩定的擾動在其演化過程中都會依次經歷下面三個階段：線性階段、非線性階段及飽和階段在線性階段，擾動的幅度較小，不同類型的擾動彼此之間并不相互作用，擾動對它所處的平衡態也無影響，這時擾動的幅度是隨時間指數增長的在非線性階段，擾動幅度增大到會反過來使原有的平衡量作一定調整（因此改變了自己得以不穩定增長的初始條件，使饋入的自由能量減少），并達到開始和其他擾動模式相互作用（從而彼此間交換能量）的程度，從而使增長率木斷下降這時擾動幅度是依次隨時間的不同冪次（一般是從高冪到低冪次）而增長的當時間的

5、冪次最后降低到零時，就達到了演化的終點擾動的幅度不再隨時間增加，而一直保持極大值，這就是飽和本章只討論磁流體的線性不穩定性線性不穩定性的基本描述方法（1）簡正模法 先將描述所研究對象的狀態量寫成平衡量（零級量）和擾動量（一級小量）之和，然后把它們代入所用的磁流體方程組，從中減去平衡方程并略去二級小量就得到了線性化的方程組對這些方程作（時間）拉氏變換和（空間）傅氏變換后可能出現下列幾種情況：（i）全部空間坐標都能進行傅氏變換這樣線性微分方程組就變成了線性的齊次代數方程組，它的有非平凡解的條件（系數行列式為零）就給出了關于的色散關系例如上一章中平板幾何位形下的阿爾文波的色散關系正是由這種方式得到的

6、（ii）只有部分空間坐標能進行傅氏變換，剩余的坐標構成了約化的微分方程組這時要設法先得到它的通解，然后利用邊條件或連接條件也可以得到的色散關系例如上一章中，柱坐標下阿爾文波的色散關系就是這樣求得的（iii）所得出的約化微分方程如果是奇異的，如上一章中連續譜阿爾文波所滿足的方程(2)能量原理（僅對理想磁流體適用）§6.2瑞利一泰勒不穩定性 這是一種經典的流體不穩定性因為這種不穩定性是由重力驅動的，故又稱重力不穩定性讓我們來研究圖3.25所示的一個容器該容器內盛有兩種不同質量密度的液體，上面的液體質量密度大，下面的質量密度小兩種流體之間有明顯的分界線顯然，質量密度梯度由下向上，受到的重力

7、由上向下，用來表示液體的平衡方程是 (1) (2)式中是流體元的速度流體達到平衡現在假定在交界面上出現了一個微擾動，其形式為 (3)這樣，密度和流體速度便可寫成： (4)從這里開始，參數下標為0表示平衡量，參數下標為1表示擾動量將（4）式代入平衡方程（3），我們得到質量守恒方程 （5）在整理上式時，已考慮到流體是不可壓縮的，將（3）式代人（5）式便得到表達式： （6）同樣可以得到擾動后的動量方程和的表達式： （7） （8）將（6）式和（8）相結合使得到如下的方程： （9）（9）式說明，當流體的密度梯度方向跟受到的重力方向相反時就會產生不穩定性，此時，這就是說重流體在上面輕流體在下面的這種平衡是

8、不穩定的只要有微擾（輕輕晃動），就會破壞原來的平衡狀態，直到達到另一種新的平衡態為止這時重流體在下，輕流體在上，正好跟原來交換了位置，所以這種不穩定性也叫做交換不穩定性現在我們采用類比的方法來研究約束在磁場中的等離子體假定磁場與等離子體之間達到了平衡，中間有明顯的分界面就是說在等離子體中沒有磁場，在磁場中沒有等離子體這時，等離子體除了受到重力之外，還受到磁場的作用力，包括磁場梯度引起的力和磁場的彎曲引起的力當然這是指單個粒子受到的力，我們把它們當作等效重力（跟流體情況作類比），記作， (10) 將 以及粒子能量代入上式并對整個麥克斯韋速度分布函數積分，我們可以得到作為流體元的等效重力： （11

9、）對干各向同性等離子體，因此 因為在低情況下 所以 （12） 將（12）式代入（9）式便得到描述瑞利一泰勒不穩定性的方程 （13）上式說明，當磁場曲率與等離子體密度梯度方向相反，即，就會產生不穩定性這種不穩定性條件也可以表示為磁場梯度與等離子體密度梯度同向，即如圖3.26（a）所示從圖中可以看出，這時的磁力線是凹向等離子體的這種曲率被稱為“壞曲率”圖326（b）畫出了穩定的磁場位形此時，磁場曲率與等離子體壓強梯度（或密度梯度）同向磁力線凸向等離子體，這種磁場位形的曲率被稱為“好曲率”在實際的磁場位形中，曲率矢量往往不斷改變方向也就是說，在某個地方是“好曲率”，在另一個地方則變成“壞曲率”如在簡

10、單磁鏡場中，在中心部位是“壞曲率”，而在“咽喉”部位則是“好曲率”因此，有必要引入“平均曲率”的概念定義： 磁力線管的比容，它是磁力線管的幾何體積與管內的磁通量的比值：，平均曲率的定義為因此，平均曲率半徑為前面得到的穩定條件（好曲率）是曲率與同向，即，在聚變等離子體中，一般都是中心密度大，即；因此穩定條件要求這就相當于要求其中為磁面包圍的體積因此，即有極大值，其中必有磁場極小值，這相當于平均磁阱這說明位于磁阱的等離子體是穩定的與之相反，位于磁山“磁山”的等離子體是不穩定的，§6.2 等離子體的能量原理 不考慮離子和電子的效應，可將等離子體作為單流體來處理。采用理想磁流體力學方程組作為

11、出發點 （1） （2） （3） （4） （5） （6）其中 表示比熱比。設每一個變量均為平衡量和擾動量的疊加，即。為簡化起見，不考慮平衡流，即。（如果可以討論）則將方程（1）（6）線性化之后可得關于一階擾動量的微分方程組 （7） （8） （9） （10） （11）令相對于流體元平衡位置的擾動位移為一階小量，則有 （12）將上式分別代入方程（7）、（9）和（10），對時間積分，可將擾動密度、擾動壓強和擾動磁場均用擾動位移來表示 （13） （14） （15）將這些表達式代入方程（8），并利用方程（11），則可得到關于擾動位移的二階微分方程 （16） （17）顯而易見，相當于由擾動位移所引起的作用在

12、單位流體體積上的力。在適當的邊界條件下解此方程，可以確定位形的穩定性。根據能量守恒原理，擾動位移引起的系統總能量的變化為零，即動能和勢能的變化之和為零 （18）將上式對時間微商可得 （19）利用擾動方程（16）和函數的自伴性，即 （20）可將方程（19）中的第一項寫成 （21）則由方程（19）可得擾動勢能的變化 （22）從直觀上來說，線性系統在力的方向上做位移時，所做的功為。根據能量守恒原理，這個功之可能是以消耗勢能為代價，因此可得方程（22）。 假設等離子體邊界為理想導電壁，一次在邊界上垂直位移。將式（17）代入方程（22）可得等離子體內部擾動勢能變化的表達式 （23）右邊被積函數中第一項總

13、是正的，代表流體可壓縮性的穩定作用；第二項在很多情況下是負的，可導致壓強梯度驅動的不穩定性；第三項總是正的，代表磁張力的穩定作用，因為彎曲磁感應線會導致磁能增加；第四項有時是負的，可導致電流驅動的不穩定性。如果等離子體邊界為真空區，在等離子體內部，根據歐姆定律（6）的線性化形式，在固定于流體上的坐標系中必有。而電場的切向分量在等離子體和真空的交界面上連續，故在等離子體之外的真空區也有 （24）即 （25）這里是垂直于等離子體邊界、法向的單位矢量，是垂直速度。在真空中可將擾動磁場用矢勢表示成，利用法拉第電磁感應定律（5）式，則有，代入方程（25）并利用可得邊界條件為 （26）在真空之外的理想導電壁上，則有 （27）根據平衡方程，在跨越等離子體和真空的交界面上的總壓強連續， （28）相應地，在位移的邊界上總壓強也必須連續 （29）將方程（29）的各項用其在平衡位置上的值展開，且只保留一級近似，則有 （30）此時，這擾動勢能除了方程（23）所示的外，還有等離子體表面的部分 （31）利用邊界條件（25）和（31）可將分成邊界上的擾動勢能 （32）表示穩定性取決于界面內外的勢能差，和真空區中的擾動勢能 （33）代表環繞等