普朗克碰撞项课件

1、类似地，沿-x方向的粒子流因为粒子密度不均匀，所以有净粒子流 设 随x变化缓慢，对x展开，只保留一级小量项， 代入，经简单计算得 式中由于粒子运动速度分布近似服从麦克斯韦分布，对麦克斯韦速度分布求平均后得应用于等离子体，则扩散系数以上讨论的是无磁场时的扩散系数。当有磁场存在时，平行磁场方向的输运系数与无磁场时相同。但在垂直磁场方向情况就不相同。 垂直磁场方向碰撞平均自由程与无磁场时完全不同。因为无磁场时两次碰撞间粒子运动是直线运动，有磁场时，则是绕磁力线的回旋运动。无碰撞时回旋中心无横向位移，只有发生碰撞时，回旋中心才会从一根磁力线跳到另一跟磁力线上，如图所示。（a）经碰撞，偏转180o， 这

2、时回旋中心位移最大一般情况如（b），回旋中心位移， ，平均讲，一次碰撞回旋中心横向移动近似为 现在用 代替式中的平均自由程l，则横向扩散系数 上式就是垂直磁场（横向）扩散系数， 就是上面已给出的无磁场时的扩散系数。一般在磁约束等离子体中， ，所以有磁场时横向扩散系数比无磁场时小很多。必须指出，当取磁场 时，不能回到无磁场时的结果，显然，这是不合理的。 如果将式子改写为因为 ，这样改写对 无多大影响，但当 时，这样结果就合理了。下面从输运方程出发计算横向扩散系数，就可得到这种合理的结果。 2. 横越磁场扩散（1）扩散方程 当流体元无整体平均运动，而且又无外电场，只是存在粒子数密度空间不均匀，则会

3、引起粒子流。根据动量输运方程 对于稳态情况 ，并假定温度 为常量， 则 ，式中摩擦阻力 ， 于是由输运方程式可以得到粒子流 代入连续性方程得当 近似地与r无关时，上式化为 这就是扩散方程， 为扩散系数，它可以从输运方程求解得到。 （2）平行磁场方向扩散 假定 ， 为常量，则平行磁场方向动量输运方程 利用 得 平行磁场方向扩散粒子流 平行磁场方向的扩散系数，结果与前面的式子相同，因0.51的修正因子，式中多了“2”因子。 （3）垂直磁场方向的扩散 由稳态的动量输运方程（垂直分量）设 ， ， =常量， ，则输运方程沿两个方向分量为方程组中消去 ，得 ，则 即式中 为垂直磁场方向的扩散系数 一般 ，

4、 表明有了磁场，横向扩散系数大大减小了。横越磁场扩散系数与磁场的平方成反比。式中最后一个等式有1/2因子，是由修正因子“0.51”引起的。上式表明，当磁场 时，可以还原到无磁场的结果。关于横越磁场扩散，在第4章4.5节做过计算和讨论，其结果是类似的。在4.5节是用（单）磁流体模型，而这里是用双流体模型。3. 同类粒子碰撞不会引起横越磁场扩散 因为扩散粒子流 ，实际上也就是动量流 ，由于同类粒子弹性碰撞动量是守恒的，不会改变其粒子的动量流，因而不引起横向磁场扩散。 下面考察两个同类粒子的碰撞，如图： 设在A点相碰， 和 为碰撞前后两粒子的 回旋半径， 和 为 碰撞前后两粒子的动量， 因为回旋半径

5、 碰撞前碰撞后动量守恒则得结果表明，碰撞后二粒子的回旋中心位置 、 ,回旋半径 、 无论怎么发生变化，但是 是不变的！ 这意味着碰撞前后回旋中心连线的中点M是固定不动的，因此不引起横越磁场扩散。由此可见，横越磁场扩散是由不同类粒子碰撞引起的。 4. 双极扩散 密度不均匀引起的扩散，其扩散系数与碰撞频率有密切关系。由于离子、电子质量差别很大，相应的碰撞频率也相差很大。因而电子和离子在无磁场时或有外磁场时的平行磁场方向或垂直磁场方向的扩散系数相差很远。因此当等离子体中出现密度梯度时，就会有一种粒子扩散速度快（粒子流大）、另一种粒子扩散速度慢（粒子流小），使原来准电中性的等离子体引起电荷分离，出现电

6、场，这种电场也会引起粒子的输运（称迁移），其效果是使原来扩散快的粒子速度减慢，扩散慢的粒子速度加快，达到准稳态时这两种粒子流速度相等。这种扩散称双极扩散，它是密度梯度扩散与电场迁移产生的总的效果。相应的电场称双极电场，而等效的扩散系数称双极扩散系数。离子和电子的扩散系数之比（设 ）为 结果说明：在平行磁场方向电子扩散系数比离子的快得多，而在垂直磁场方向电子的扩散系数比离子的慢得多。这样等离子体会引起电荷分离，出现电场，因此就会出现前面分析的双极扩散。可以证明，双极扩散方程为 式中 为双极扩散系数。 用输运方程研究双极扩散问题。稳态的输运方程 式中 是由等离子体空间电荷产生的，摩擦力 设 ， ，

7、 ， 则稳态的输运方程可写为 由此解得 第2式，是忽略了高级小量项 后得到的，因为此项多个因子 。由此得 现在分别讨论几种情况：（1）同类粒子碰撞，即 ， 说明同类粒子碰撞不引起横越磁场扩散，现在从输运方程证明了前面论述过的结论。（2）不同类粒子碰撞，令=p（质子），=e（电子），假定 ， ， 因为 则得 结果与第4章用磁单流体力学得到的（4.4.16）式扩散系数 相同。现在 、 给出的横越磁场扩散系数，那里没考虑等离子体空间电荷产生的电场作用。 现在结果表明，垂直磁场方向电子、离子是以相同的速度向 方向扩散，称双极扩散；双极扩散系数不依赖于电场强度，而且比 大一倍，但比 小很多，即原先扩散很

8、慢的电子，现在稍微加快了，为原来的两倍，而原先扩散很快的离子，现在减慢了很多，最终电子、离子扩散速度一致，粒子流相同，双极扩散系数主要取决于扩散较慢的粒子成分。6.8 环形磁场的新经典扩散 前面讨论的输运过程理论，通常称为经典理论，它的基础是碰撞十分频繁，碰撞频率远大于其它特征频率，平均自由程远小于系统内的其它特征长度，因而可以应用磁流体力学方法处理等离子体中的输运过程。经典输运理论可以解释一些等离子体中的输运现象，但在热核等离子体中，如环形的磁场装置，经典理论与实验结果相差甚远，经典理论的输运系数比实验测量结果要小几个量级。其主要原因是在环形装置中粒子运动轨道对输运系数有很大影响，因此输运理

9、论必需考虑粒子运动轨道的特征。这种考虑了等离子体环形效应的经典输运理论称新经典输运理论。 第4章用磁流体力学得到的扩散系数和上节应用双流体力学方程得到的双极扩散系数完全相同: 现在称 为经典扩散系数，其中 为电子特征热速度， 就是特征热速度相应的电子回旋半径。经典扩散系数可以解释为：电子每经历一次碰撞在垂直磁场方向无规地移动一个电子回旋半径，而扩散系数就是单位时间内碰撞（碰撞频率）产生的无规移动平方平均值的总和。经典扩散系数公式都只是适用于平直磁场情况。对于环形磁场系统，应当考虑粒子运动轨道的影响。托卡马克装置中粒子的轨道运动。由于环向磁场 的不均匀性和磁力线旋转变换，在环的外侧磁场较弱，环的

10、内侧磁场较强，因此在磁面上沿一根磁力线运行的粒子感受到磁场强弱变化，相当于磁镜场结构。当粒子的初始速度平行分量和垂直分量比值较大时，这种磁场强弱变化不影响粒子的“通行”，可以在整个环形等离子体中运动，称通行粒子或环形粒子。但当较小时，就会出现粒子由外侧到内侧运动时被强磁场区域反射，被捕获在磁场较弱的磁阱中，称这种粒子为捕获粒子。根据磁矩守恒，可得到划分这两类粒子的界线。设在外侧 粒子的初始速度为v，平行分量v 、垂直分量v 到内侧 粒子速度为v，平行分量v=0，垂直分量v= v 即全部速度都是垂直分量。磁矩守恒因 ，由此得通行粒子：捕获粒子： 如果在速度空间粒子分布是各向同性的，则捕获粒子所占

11、的比例 因此捕获粒子所占的比例是很小的。若n0 为粒子数密度，则捕获粒子密度 。两类粒子运动轨道特性 设由磁面上A点出发的通行粒子，由于磁力线曲率和梯度引起漂移，当v0 时，粒子的导向中心沿磁面内的一个小圆形漂移面运动，当v0时，则沿磁面外的一个大圆漂移面运动，小圆漂移面中心向外侧移动，大圆漂移面中心向内侧移动从磁面上A点出发的捕获粒子，其导向中心轨道是一个香蕉形， v0时，粒子的导向中心轨道是沿磁面内的小香蕉， v0时，是沿磁面外的大香蕉，大小香蕉的半宽度大致相等，可近似地取为 根据经典扩散系数公式， 粒子经历1次碰撞平均无规移动 现在考虑粒子运动轨道影响，即新经典扩散， 通行粒子平均无规移

12、动是 捕获粒子平均无规移动是 由于扩散系数是与平均无规移动的平方成正比，应该用新的平均无规移动来代替，因而粒子运动轨道对扩散有很大影响，而且扩散系数与碰撞频率也有密切关系。现在考察这两类粒子的特征频率 （i）通行粒子特征频率 （i）通行粒子特征频率 通行粒子绕小圆环一周时在大环上走过的距离 是旋转变换角。 粒子沿磁力线的速度为 ，所以通行粒子绕小环一周所需的时间 因为通行粒子的 较大，粒子沿磁力线运动时变化不大，上式 可用特征热速度 代替。因此通行粒子绕小环的周期和频率分别为 式中 、 为通行粒子的特征时间与特征频率。（ii）捕获粒子特征频率 捕获粒子绕香蕉轨道一周所需的时间 但是捕获粒子的

13、较小，而且在香蕉轨道一周内 变化很大，它应作为变量，而且在反射点 ，所以上式中 应该用香蕉轨道一周上的平均值 代替。即绕香蕉轨道一周的平均时间 捕获粒子满足条件 由此估计捕获粒子特征时间 与 相比 表明，捕获粒子绕一个完整的香蕉轨道所需的时间比环行粒子绕小环一周的时间要长。现在定义捕获粒子有效碰撞时间和有效碰撞频率。有效碰撞时间为捕获粒子在香蕉轨道中经历碰撞而破坏了原有捕获条件的时间，即图中捕获粒子 偏转角改变 、使 所需的时间，而不是电子碰撞使 在速度空间偏转 、 需要的时间 因此 有效碰撞频率有效碰撞频率实际就是捕获粒子逃脱一个香蕉轨道而跳到另一香蕉轨道的频率。为使一个捕获粒子能够完成一个

14、香蕉轨道的运动，有效碰撞时间 有效碰撞频率 利用 ，上式条件变为 为捕获粒子的特征频率，其意义是当电子碰撞频率 时，捕获粒子至少可以完成一个完整的香蕉轨道运动。 根据以上两类粒子的特征频率 和， 而且 ，可以把电子碰撞频率分为三个区域：（i） 称流体区 在流体区，碰撞频繁，磁流体力学方法可以适用。在这个区域，只有饶小环形的通行粒子可以存在，而绕大一点环形的通行粒子和捕获粒子都完不成一个闭合的轨道就被碰撞离开原有的环形或香蕉轨道。这部分粒子的扩散仍可用经典扩散公式，但仍有大量绕小环形粒子存在，其扩散系数应考虑轨道效应。因此经典扩散系数公式中的电子回旋半径 （电子特征热速度相应的值）应改用 即新经

15、典扩散系数 由于在这个区域还有大量粒子，在碰撞时间内没有完成任何闭合轨道，应考虑这部分粒子的经典扩散的贡献，所以在流体区的总扩散系数应为部分新经典扩散公式和部分经典扩散公式相加，即 上式称斐斯-舒鲁特（Pfirsch-Schluter）扩散系数。这个公式可以直接由磁流体力学方程导出。 因此上式结果比经典扩散的公式增加了 因子。如果设 ，则扩散系数增加了1个量级。 （ii） ,称香蕉区 在这个区域电子碰撞频率很低，捕获粒子有足够时间完成整个香蕉轨道运动。因为每一次有效碰撞无规移动 ，有效碰撞频率 还应考虑捕获粒子在总粒子数密度中占的比例 根据经典扩散系数公式，捕获粒子扩散系数 上式与碰撞频率很大

16、时的流体区新经典扩散系数相比，增加了因子 ，对大部分约束区域 。因此，尽管捕获粒子占的比例很小，但它比经典扩散系数大了2个量级。通行粒子的新经典扩散虽然比经典情况大1个量级，但还比捕获粒子的新经典扩散小1个量级，所以在托卡马克装置中，研究微观不稳定性的作用机制时，捕获粒子效应往往是起决定性作用。（iii） ，称过渡区，也称平台区。 在这个区域电子不能完成香蕉轨道，上面简单的处理方法不适用，而且磁流体力学方法也不适用。严格讲，香蕉区和平台区的输运过程都应该从动理学方程出发来研究，得到过渡区的扩散系数 在这个区域扩散系数与电子碰撞频率无关，平台区就是由此而得名。流体区的公式，当取，时和香蕉区公式，

17、当取时的扩散系数与平台区公式相等。 三个区域扩散系数随电子碰撞频率的关系如图所示。这就是托卡马克系统中的新经典扩散系数的结果。 在过渡区，扩散系数曲线为平台型，在物理上可做些定性说明: 首先在平台区两端分界点 和 扩散系数相等，这为平台区形成的前提。香蕉区流体区平台区从下端分界点 开始，随电子碰撞频率 逐渐增加，相应的捕获粒子的有效碰撞频率 也相应增加，这个因素使扩散系数增大，但 的增加，最初使许多大香蕉粒子完不成香蕉轨道运动，频率再增高，中香蕉粒子、小香蕉粒子也相继消失，这样香蕉粒子效应逐渐减弱以致完全消失，使扩散系数逐渐减小。这两种因素的消长相抵，大体保持扩散系数不随电子碰撞频率变化。从上

18、端分界点 开始，随 减小，扩散系数也相应减小，但绕环形的通行粒子数量不断增加，开始是绕小环形的增加，后来绕中环形的、大环形的通行粒子相继出现，数量也不断增长，于是绕环形的通行粒子效应增强使扩散系数增大。同样，这两种因素消长相抵，基本上保持扩散系数为常量。因此在过渡区，扩散系数出现平台型曲线。6.9动理学方程 前面讨论的等离子体输运过程和输运系数的计算比较简单、但都比较粗糙，因为输运过程是一个非平衡态过程，严格地讲，它必须从粒子分布函数的变化规律出发，用统计力学方法确定系统的全部物理状态和性质及其随时间的变化过程。设系统的单粒子分布函数为 ，则描述分布函数在相空间中演化的动理学方程 表示类粒子的

19、分布函数， 表示一个类粒子所受的力场，它包括外场和等离子体内部的平均场（自洽场） 碰撞项表示因碰撞引起的分布函数的变化，对于多种类型粒子： 表示各种类型粒子的碰撞效应总和，包括同类粒子。碰撞项可以用各种模型导出，由此得到相应的碰撞积分，因而动理学方程就有各种不同的形式。1. BGK方程（Krook碰撞项） 碰撞项取为如下形式 称BGK（P.L.Bhatnagar, E.P. Groos, M.Krook）碰撞项或称Krook碰撞项，式中 为特征时间，通常取为平均碰撞时间， 为平均碰撞频率。BGK碰撞项的意义为系统从初始的非平衡态速度分布，经过时间 趋向平衡态分布，所以 就是这一过程的弛豫时间。

20、取BGK碰撞项后的动理学方程为(省略下标 ) 一般地 或 是速度的函数，而且与分布函数有关，现在碰撞项中的 或 都取为常量，这样碰撞项计算就大为简化，因此BGK碰撞项也称弛豫时间近似。例如：如果系统分布偏离平衡态不远，则分布函数可以写成线性化形式 f0 为平衡态速度分布，f1 为偏离平衡态的小量。 f0 满足的动理学方程为 对于局域性热力学平衡情况，方程的解为 由于系统偏离平衡态不远， BGK动理学方程线性化后为 上式左方第1项 ，第2、3项忽略小量f1的贡 献。对于定态情况， ， 则方程的定态解为式中 f0 为已知，于是各种输运流的表达式： 将解 f1 代入上式，计算各种输运流，就可得到各种

21、输运过程系数。 2. 玻尔兹曼方程 这是最早的一种碰撞模型，它是建立在二体碰撞基础上的一种碰撞积分表示式。在质心系中，二粒子的碰撞相当于质量为的一个粒子被固定的散射中心的作用而改变速度方向。每秒入射到面积元 上粒子数场粒子数为 ，因此每秒单位体积中碰撞数为 粒子的减少率 逆过程每秒单位体积内的碰撞数为 上式为二粒子发生碰撞，使粒子进入速度区间 的粒子数，积分后为 粒子的增加率。因为正逆过程的作用完全对称，所以 ，由弹性碰撞动量能量守恒， ，根据刘维定理 因此，积分后，得 粒子的增加率碰撞项为（净增加率）利用微分截面定义则碰撞积分改写为 这就是玻尔兹曼碰撞积分。如果有多种粒子，碰撞项还应该对求和

22、（包括同类粒子）。用玻尔兹曼碰撞积分的动理学方程称玻尔兹曼方程，即 式中 为 ， 为 ，其余类似。微分截面如用库伦作用 也可用其他形式结果。 注意：导出玻尔兹曼碰撞积分实际含有如下假设：（i）所有碰撞都是二体碰撞；（ii）相互作用长度远小于发生显著变化的长度；（iii）碰撞持续时间远小于发生显著变化的时间。 3. 福克-普朗克方程 玻尔兹曼碰撞积分应用于等离子体时的主要缺点是假定是短程二体碰撞。等离子体中的带电粒子是库仑长程相互作用，粒子间的碰撞大部分是小角度散射，而且每个粒子同时与周围许多粒子相互作用，它的大角度偏转大部分是由小角度偏转积累的效应。因此引用20世纪初在处理布朗粒子运动时导出的

23、碰撞积分比较合适。因为布朗粒子质量大，受到周围流体分子碰撞时，每个时刻速度改变量都很小， ，因此可以把 作为小量对布朗粒子分布函数作泰勒展开，这样得到碰撞项为微分形式。用这种方法得到的动理学方程称福克-普朗克（Fokker-Planch）方程。设速度为v的粒子、由于碰撞在时间t 内速度增量为v 的几率为 ，这里假定w与时间无关，即与粒子过去的历史无关，这种过程一般称马尔科夫过程。由w的定义，可得粒子的分布函数上式对小量v做展开所有可能的增量的总几率为1, 碰撞产生的分布函数的变化率为将前面式子代入上式，保留到二阶小量项，则得式中 这就是福克-普朗克碰撞项，代入动理学方程后就是福克-普朗克方程。

24、在福克-普朗克碰撞项中包含动摩擦项（第1项）和扩散系数项（第2项）。 为动摩擦系数，它表示碰撞引起的速度慢化； 为扩散系数，它是碰撞把初始单一方向速度的粒子在速度空间扩展开来。福克-普朗克碰撞项是微分算子，因此它相应的动理学方程是微分方程，所以比玻尔兹曼微分积分方程更容易求解。求解福克-普朗克方程，首先要计算动摩擦系数和扩散系数，这就需要知道适合等离子体的转移概率 的表示式。到目前为止，已有多种表达式，但也都是假定多体相互作用当作一系列二体碰撞，这样才能写出明显表达式。如果假定二体碰撞，则和前面推导玻尔兹曼碰撞积分的方法一样，可以导得转移概率函数，并求出相应的摩擦系数和扩散系数。 4. Ros

25、enbluth势碰撞项可以用Rosenbluth势表达福克普朗克碰撞项。 设速度为v1、分布函数为 的粒子与速度为v2、 分布函数为 的粒子的相碰撞，第1个粒子被散射到立体角元 内的几率为 现在 中的一个粒子在单位时间和一群速度在 间 的粒子相碰，被散射到立体角 内的转移几率为 式中 是与散射角对应的速度改变量。福克-普朗克的2个系数可以表示为 由二体碰撞运动学给出：相对速度 ，质心速度 ，得 如图，直角坐标系三个单位矢量将上面有关各式代入，对 积分得利用 并引入两个速度函数则2个系数公式简化为函数 、 称为Rosenbluth势。因为它们满足速度空间泊松方程，和坐标空间静电势相类似。 将这2

26、个系数公式代入福克-普朗克碰撞项，就得Rosenbluth势表达的福克普朗克碰撞项。5. 朗道方程 1936年，朗道取库仑势场的散射微分截面，对玻尔兹曼方程碰撞积分作小的 展开，保留到二阶小量项，得到了朗道碰撞项： 因为朗道碰撞项还可以写成 应用朗道碰撞项的动理学方程称朗道方程。朗道碰撞项是最早导出的适合库仑长程作用条件下二体碰撞项。朗道方程与具有Rosenbluth势的福克-普朗克方程完全等价的，只是形式上更对称。 以上介绍了几种碰撞项，实际上用动理学方程求解输运问题时，无论用哪个碰撞项都可以，只不过对不同问题选择不同的碰撞项形式可能会更方便些。6. 符拉索夫方程如果在动理学方程中，忽略碰撞

27、项，即 称为符拉索夫方程或无碰撞玻尔兹曼方程。这个方程适合于讨论特征时间远小于碰撞时间或特征长度远小于平均自由程情况，因为在这种情况下，碰撞效应完全可以忽略。 但需指出，式中的力场包含外场和等离子体内部的平均场（自洽场），因此方程要与电磁场的麦克斯韦方程联合求解讨论等离子体微观不稳定性等问题时，常从这个方程出发求解。 6.10 波与粒子相互作用， 朗道阻尼与朗道增长没有碰撞情况下，也可能存在波和粒子间的能量交换机制，波的能量可能被粒子吸收而衰减，也可能出现逆过程。波与粒子作用机制对波加热等离子体或波的不稳定性都是非常重要。现在用简单的物理模型来说明波与粒子相互作用机制的物理概念。 设一群初始速度v0不同的 带电粒子，沿磁力线漂移， 并设 的静电波（纵波） 沿平行磁力线方向传播，如图。 如果粒子群的漂移速度与静电波传播速度相近，则粒子群与静电波之