《第3章 磁流体力学方程.doc》由会员分享，可在线阅读，更多相关《第3章 磁流体力学方程.doc（15页珍藏版）》请在金锄头文库上搜索。
1、第三章 磁流体力学方程（MHD）3.1引言 由上一章的讨论可以看出，等离子体动力学理论是在位形及速度空间中讨论带电粒子的分布函数随时间的演化规律。由于动力学方程是一个非线性的积分微分方程，数学处理较复杂，在一般情况下很难求解。实际上，我们可以把等离子体看成为是一种电磁流体，它的宏观状态可以用密度、流速、温度等状态变量及电磁场来描述。这些状态参量及电磁场是在三维位形空间中随时间演化的。建立电磁流体状态参置随时间的演化方程称为磁流体力学（Magnetohydrodynamics-MHD）。与动力学理论相比，磁流体力学在数学处理上简单的多，而且等离子体中的许多过程，如等离子体的宏观平衡与稳定，波动过
2、程均可以用MHD理论来描述。但对于等离子体中的另外一些现象，如Landau阻尼、速度空间中的不稳定性等则MHD理论却无能力描述。下面我们从动力学方程出发，建立MHD方程。3.2二份量MHD方程设等离子体是由电子成份和一种离子成份组成的二份量电磁流体。首先我们引入二份量磁流体的宏观状态变量，我们知道，对于一个多粒子系统，其宏观变量是对应的微观变量的统计平均值。这样，第类成份流体的密度、流速火及温度的定义为： （3-1） （3-2） 下面我们利用上章给出的等离子体运动学方程来建立MHD方程。动力学方程可以写成： (3-3)首先定义等离子体矩方程：将（3-3）两边乘以并对积分，（1） （2） （3）
3、 其中用到了分部积分和在时为零的条件。（4） 其中利用了关系：这样得矩方程：其中：为统计平均。1 连续性方程设，并对积分,则 （3-4）其中利用到，粒子数守恒。引入电荷密度： （3-5）和电流密度： （3-6）将（3-4）两边乘以可以得到电荷守恒方程 （3-7）将（3-4）两边乘以可以得到质量连续性方程 （3-8）其中是质量密度。2 动量平衡方程设，并对积分,则可得 （3-9）其中 (3-10)为压强张量。而 （3-11）利用连续方程（3-4），方程（3-9）可以化成为 （3-12）该方程中各项的物理意义是：-流体元的动量变化率；其中 -为对流项； -压强梯度产生的力； -电场力； -洛仑兹力
4、，是由电流穿越磁场而产生的力； -为第类粒子与第类粒子碰撞时，其动量的变化率。方程（312）是一个不封闭的方程，因为涉及到高阶矩函数及，只有通过求解动力学方程，才能严格地计算出及。在研究等离子体的磁流体状态时，通常假定等离子体中带电粒子的速度分布基本上为各向同性分布，因此有： （3-13）其中为静压强。的非对角部分仅与等离子体中的粘滞现象有关。另外，对于不同种类的带电粒子之间的碰撞，动量的变化率可以写成摩擦阻力的形式： （3-14）其中为动量输运的平均碰撞频率.3能量平衡方程设，并对积分,并利用连续方程（3-4），动量平衡方积（3.12），最后可以得到能量平衡方程为： (3-15)其中： （3
5、-16）为热流矢量，而， （3-17）为不同种类带电粒子之间的碰撞产生的能量变化。 在高温等离子体系统中，人们对等离子体中带电粒子的能量输运并不是太感兴趣。在研究等离子体的平衡、稳定及波动过程时，可以认为带电粒子在速度空间的分布基本上趋于各向同性的Maxwell分布。因此，通常用状态方程来确定等离子体的压强，从而取代了能量平衡方程。对于等温过程，有： （3-18）其中c1是常数。对于绝热过程，压强为 （3-19）其中c2是常数。这样，对于双流体等离子体，其MHD方程为： （3-20） （3-21） （3-22） （3-23） （3-24） （3-25） （3-26） （3-27）它们与状态方程
6、耦合，即构成一套封闭的方程组。后面几章，我们将用这套方程组研究等离于体中的波动过程及稳定性。3.3单MHD方程 在上节中，我们是把等离子体看作是由电子流体和离子流体组成的双流体。实际上，在研究等离子体中某些现象时，也可以把等离子体看成为单一的磁流体。本节我们的任务就是给出这种单一磁流体的MHD方程。首先引入单一磁流体的宏观状态参量：质量密度： （3-28）电荷密度： （3-29）流速： （3-30）温度： 电流密度； （3-31）总压强： （3-32）下面建立单流体的流体力学方程（1） 连续性方程将电子成分的质量连续性方程 （3-33）与离子成分的质量连续性方程 （3-34）相加，并利用（3-
7、28）及（3-30），则单流体的质量连续性方程为 （3-35）（2） 动量平衡方程为了得到单一流体的动量平衡方程，我们假定：等离子体是准中性的，即。这样根据电子和离子的动量平衡方程， （3-36） （3-37）得到单一流体的动量平衡方程为 （3-38）其中利用了如下简化假设： 由于电子的质量比离子质量小的多，略去了电子的惯性项，和电中性条件：，及，。（3） 广义欧姆定律由关系：， ，得：，由关系：，得：略去（3-36）中的对流项，得： (3-40)这就是广义欧姆定律。对于简单的欧姆定律有 （3-41）是等离子体的电导率。因此，广义欧姆定律中，多了如下几项：（1），磁流体运动引起电流；（2）：等
8、离子体受到洛兹力作用而运动产生的电流。（3）：由于压力梯度而产生的电流变化。这样采用单流体模型，等离子体的MHD方程为： （3-43）再加上Maxwell 方程组 （3-44）构成了一套封闭的方程组（设P已知，由状态方程给出）。单流体MHD方程常用于描述等离子体装置的平衡与稳定，也可以用于描述等离子体中的波动现象。用单流体的MHD方程描述等离子体中的波动现象比用双流体的MHD方程描述精确性差，这是因为在推导单流MHD方程时，做了一些简化假定。对于一些低频过程，可以略去Maxwell方程组中的 ；另外对于低温过程，还可以略去广义欧姆定律中的压力梯度项。 项与项相比，也可以略去。这样简化的单流体M
9、HD方程为： （3-45） （3-46）对于理想等离子体，还可以略去项，这样理想等离子体的MHD方程为： (3-47)下节我们将利用上述方程研究等离子体的MHD平衡与稳定。 3.4 MHD平衡与稳定1平衡与稳定的概念 对于一个磁约束的等离子体系统，人们所关心的首要问题是系统的受力是否平衡？如果受力能达到平衡，接下来的问题的就是这种平衡状态是否稳定，即系统受到小的扰动后，对平衡状态的偏离幅度是否大？下面我们用数学的语言来描述平衡与稳定问题。设x0是系统处在平衡状态下的某一物理量，满足如下方程： （3-48）如果对系统施加一个扰动，使得物理量x0变为，且满足如下运动方程： （3-49）设扰动是微扰
10、动，则可以将方程（3.49）线性化： （3-50）设扰动量可以表示成山的形式，由此可以求出扰动频率为： （3-51）如果，则表明扰动会很快地被衰减掉，则系统对应的平衡状态是稳定的。反之，则表明平衡状态是不稳定的。 我们也可以用下图4形说明平衡与稳定问题。图4a对应的是稳定平衡，因为小球偏离平衡位置后，仍能恢复到原来的状态；而图4b则对应的是不稳定的平衡，小球一旦受到扰动，将会逐渐地偏离其平衡位置，不会自动复原。下面利用上节得到的简单的MHD方程，来分析以下等离子体体系的平衡问题。2. MHD平衡 对于一般的磁约束系统，研究其平衡问题是极其复杂的。下面我们从稳态的MHD方程出发，给出等离子体平衡
11、状态的一些性质。 对稳定的状态，等离子体系统的MHD方程为(见3-7)： (3-52)Maxwell方程为： (3-53)由方程（2-52）可以看出：（a）抗磁性电流压力梯度和洛仑兹力是维持系统处在平衡状态的必要条件。为了说明这个问题，我们以一个柱状的等离子体系统为例，压强梯度指向轴心，见图5，为了消除等离子体向外膨胀，就必须在极向上施加一电流j，它与磁场作用产生的洛仑兹力正好与膨胀力相抵消。用叉乘方程（3-52）两边，可以得到： （3-54）这个电流称为抗磁性电流。从单粒子的观点看，抗磁性电流是由于粒子在磁场中做拉莫回旋运动时，受到密度梯度的作用使导向中心发生漂移而形成电流的结果。从流体力学的观点看，抗磁性电流是由于穿越磁场的压强梯度形成的，产生的电流与磁场作用形成的洛仑兹力恰好与压强力相抵消。（b）等压面由方程（3-52），说明沿j和B方向压强是不变的，电流和磁场的方向都在等压面上，j和B都与垂直。当人们在设计复杂的磁场位形时，必顺考虑到这一点。设想有一个环形的等离子体系统，见图6，存在一个径向的压强梯度，即等离子体由一层层等压面所包围。一般地讲，磁力线和电流线可以这样和那样的弯曲，但不能穿越等压面。（c）
