托卡马克物理基础-3[].docVIP

III. 托卡马克平衡位形的求解 这一部分主要介绍Grad-Shafranov方程的求解方法及相关的物理内涵。重点是内部解的求解方法。在磁面截面形状接近为圆时（磁面簇为套圆环系统），如果环境比，我们可以把倒环境比作为展开小量，获得近似解。在另外一些特殊电流剖面下，还可以获得精确的位形解。 III.1 圆截面托卡马克 圆截面托卡马克不仅是一种数学近似，而且也是一些真实托卡马克装置的位形。早期的托卡马克外部有一个很厚的铜壳，在等离子体环发生向外位移时会自动产生象电流，这个象电流的垂直磁场用以维持环的平衡以及稳定一些快速增长的磁流体不稳定模（扭曲模）。Shafranov本人对这类平衡位形做了大量的研究工作，奠定了托卡马克位形的平衡理论基础。所以，平衡位形中各个磁面之间的相对位移被称为Shafranov位移。 在圆柱坐标系中写出Grad-Shafranov方程： （3-1） 先介绍大环境比近似展开法。 III.1.1大环境比展开法 设，这时， 以磁轴为原点，引入局部极坐标： （3-1-1） 图 局部极坐标 取为平均磁面半径 （3-1-2 将展开为Fourier级数： （3-1-3 对于上下对称的位形，上式中右方第二项为零。系数为Shafranov位移函数，与截面的椭圆变形有关，与三角变形有关，等等。在坐标系中，方程（3-1）化为： （3-1-4） （3-1-5） （3-1-6） 式中下标表示对该变量求偏导数。将（3-1-4）对求平均，得到 （3-1-7） 这里。在（3-1-7）式中，三个平均量中含有未知函数（n=1,2,3,..）,因此，必须将（3-1-4）式再乘以后对求平均（下面仅考虑上下对称位形），得到足够的补充方程。这些方程虽然是常微分方程，但是非现性的而且是相互偶合的。为研究解析解，我们在大环径比近似下采用展开方法。为使展开方法有效，需假定，这里。最低一级给出关于磁通函数的方程： （3-1-8） 上式其实就是直圆柱中沿圆柱对称轴流动的电流产生的极向磁通方程。不象原先的二阶非线性偏微分方程，现在对任意的电流剖面分布，求解都非常简单。一次积分给出： (3-1-9) 特别是，当电流为均匀时，得到 。 求出后，再利用环向电流密度的表示式 (3-1-10) 物理实验中压强剖面是原则上可以直接测量的量，环向电流剖面也是原则上可以测定的量（虽然测量技术非常高级）于是，极向电流通量函数F也就可以确定。 下一级近似是关于变形参量的方程，这些方程可以通过对（3-1-4）乘再求平均得到。其中实际代表各个近似的圆环面之间的位移，有时直接写成（误差为高级小量），与截面的拉长比有关，而与三角变形有关。与文献中常用的展开形式对比，如下面的展开形式： 这两种变换的关系是 ， 位移函数所满足的二阶常微分方程是： （3-1-11） 这里需要压强对平均小半径的剖面函数的信息。原则上，这也是可以由实验测定的。这个方程具有物理意义的解应该满足初始条件 ，一般解可以表示成左方非齐次项的积分形式。如果电流我准均匀分布，这个方程简化为： 这里引入安全因子 （3-1-12） 其解为 上式右方第一项与等离子体是否存在无关，代表真空位移效应；第二项是压强梯度引起的，称为Shafranov位移。主要总位移是负值，表示内部磁面相对于磁轴位置是向内测移动。内部电流剖面也随磁面位移有所位移，这要影响等离子体的内电感值；磁面边界与外导体壁之间也有相对位移，这是实验可以测定之值。在气球模稳定性的研究中，上面出现的参量经常记为 （3-1-13） 这个参量是稳定性的标志性参量。 关于变形参数的方程为 () （3-1-14） 这是齐次常微分方程，我们需求出内部区域中为正则的解。对准均匀分布电流，可以得到下面形式的解 这表明，拉长变形系数（n=2）基本上是对所有磁面都接近的；而三角变形系数（n=3）随小半径减小。靠近磁轴处三角变形消失。这些虽然是对准均匀电流剖面情况下得出的，但对其他任意电流剖面，也应具有类似的性质。 （3-1-12） 中定义的安全因子一般称为圆柱安全因子，记为。真实环形等离子体的安全因子不仅与电流剖面有关，而且与环境比有