第3章 磁流体__力学方程

《第3章 磁流体__力学方程》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第三章磁流体力学方程（MHD）§3.1引言由上一章的讨论可以看出，等离子体动力学理论是在位形及速度空间中讨论带电粒子的分布函数随时间的演化规律。由于动力学方程是一个非线性的积分微分方程，数学处理较复杂，在一般情况下很难求解。实际上，我们可以把等离子体看成为是一种电磁流体，它的宏观状态可以用密度、流速、温度等状态变量及电磁场来描述。这些状态参量及电磁场是在三维位形空间中随时间演化的。建立电磁流体状态参置随时间的演化方程称为磁流体力学（Magnetohydrodynamics-MHD）。与动力学理论相比，磁流体力学在数学处理上简单的多，而且等离子体中的

2、许多过程，如等离子体的宏观平衡与稳定，波动过程均可以用MHD理论来描述。但对于等离子体中的另外一些现象，如Landau阻尼、速度空间中的不稳定性等则MHD理论却无能力描述。下面我们从动力学方程出发，建立MHD方程。§3.2二份量MHD方程设等离子体是由电子成份和一种离子成份组成的二份量电磁流体。首先我们引入二份量磁流体的宏观状态变量，我们知道，对于一个多粒子系统，其宏观变量是对应的微观变量的统计平均值。这样，第类成份流体的密度、流速火及温度的定义为：（3-1）（3-2）15下面我们利用上章给出的等离子体运动学方程来建立MHD方程。动力学方程可以写成

3、：(3-3)首先定义等离子体矩方程：将（3-3）两边乘以并对积分，（1）（2）（3）其中用到了分部积分和在时为零的条件。（4）其中利用了关系：这样得矩方程：其中：为统计平均。1．连续性方程设，并对积分,则（3-4）其中利用到，粒子数守恒。引入电荷密度：（3-5）和电流密度：（3-6）15将（3-4）两边乘以可以得到电荷守恒方程（3-7）将（3-4）两边乘以可以得到质量连续性方程（3-8）其中是质量密度。1．动量平衡方程设，并对积分,则可得（3-9）其中(3-10)为压强张量。而（3-11）利用连续方程（3-4），方程（3-9）可以化成为（3-12）

4、该方程中各项的物理意义是：---流体元的动量变化率；其中--为对流项；--压强梯度产生的力；--电场力；--洛仑兹力，是由电流穿越磁场而产生的力；--为第类粒子与第类粒子碰撞时，其动量的变化率。方程（3．12）是一个不封闭的方程，因为涉及到高阶矩函数及，，只有通过求解动力学方程，才能严格地计算出及。在研究等离子体的磁流体状态时，通常假定等离子体中带电粒子的速度分布基本上为各向同性分布，因此有：（3-13）15其中为静压强。的非对角部分仅与等离子体中的粘滞现象有关。另外，对于不同种类的带电粒子之间的碰撞，动量的变化率可以写成摩擦阻力的形式：（3-14

5、）其中为动量输运的平均碰撞频率.3．能量平衡方程设，并对积分,并利用连续方程（3-4），动量平衡方积（3.12），最后可以得到能量平衡方程为：(3-15)其中：（3-16）为热流矢量，而，（3-17）为不同种类带电粒子之间的碰撞产生的能量变化。在高温等离子体系统中，人们对等离子体中带电粒子的能量输运并不是太感兴趣。在研究等离子体的平衡、稳定及波动过程时，可以认为带电粒子在速度空间的分布基本上趋于各向同性的Maxwell分布。因此，通常用状态方程来确定等离子体的压强，从而取代了能量平衡方程。对于等温过程，有：（3-18）其中c1是常数。对于绝热过程，

6、压强为（3-19）其中c2是常数。这样，对于双流体等离子体，其MHD方程为：（3-20）（3-21）15（3-22）（3-23）（3-24）（3-25）（3-26）（3-27）它们与状态方程耦合，即构成一套封闭的方程组。后面几章，我们将用这套方程组研究等离于体中的波动过程及稳定性。§3.3单MHD方程在上节中，我们是把等离子体看作是由电子流体和离子流体组成的双流体。实际上，在研究等离子体中某些现象时，也可以把等离子体看成为单一的磁流体。本节我们的任务就是给出这种单一磁流体的MHD方程。首先引入单一磁流体的宏观状态参量：质量密度：（3-28）电荷密度

7、：（3-29）流速：（3-30）温度：电流密度；（3-31）总压强：（3-32）下面建立单流体的流体力学方程（1）连续性方程将电子成分的质量连续性方程（3-33）15与离子成分的质量连续性方程（3-34）相加，并利用（3-28）及（3-30），则单流体的质量连续性方程为（3-35）（1）动量平衡方程为了得到单一流体的动量平衡方程，我们假定：等离子体是准中性的，即。这样根据电子和离子的动量平衡方程，（3-36）（3-37）得到单一流体的动量平衡方程为（3-38）其中利用了如下简化假设：由于电子的质量比离子质量小的多，略去了电子的惯性项，和电中性条件：

8、，及，，，。（2）广义欧姆定律由关系：，，得：，由关系：，得：略去（3-36）中的对流项，得：(3-40)这就是广义欧姆定