磁流体力学平衡.docVIP

4.7.2 正交曲线坐标变换 设为两组正交曲线坐标， 对应的单位矢量为：，满足， 长度元： ， 为拉梅系数， ， 线元为： ， 微分面元： ， 体积元： ， 其它微分算符关系为： ， 例题1:柱坐标变换： 例题2：对于环形坐标：从大柱坐标（）到小柱坐标（） ， 或 ，， ， 4.7.3 二维问题的平衡方程 1．磁场、电流及安全因子的表达式 对tokamak这类的轴对称系统有 其中环向磁场为，而是不依赖于的常量．这样 ， 定义：，是引进标量函数 在大柱坐标下（）：，，代入定义中，得到关系 可见 正好满足 在小柱坐标下（）中：， ， 得到： 正好满足： 最后一项为零：，以及关系，而有角向磁通的物理意义 其中S是平行于赤道面、高为Z的而半径为R的圆面，所以最后总磁场可写成 把上式代入安培定律，可得电流的表达式分别为 （4.43） 其中是二维拉普拉斯算符 （4.44） 而角向电流为，其中 积分后一方程得 其中S是高为Z、半径为R的圆面，所以穿过此面的电流正是相应磁面之外的等离子体电流，它也是磁通的函数．由此可知环向磁场和I可以互相表示成 （4.45） 而角向电流及它的两个分量也可用它来表示 （4.46） 结果总电流可写成 （4.47） 最后来给出安全因子的表达式．用坐标，假设当从时，磁力线在大环方向上转了，按定义安全因子。而为了求，用力线方程，可以有 上面用到关系。将上面代入安全因子的定义后可得 （4.48） 2．二维 Grad－Shafranov方程 GS方程是平衡时角向磁通所满足的微分方 程，由（4.43）式有 现在只要把方程右面的环向电流设法用表示出来，就可以得到一个关于的二阶椭圆方程．而这可以通过平衡方程来做到，由得 将上面已知的表达式 ，， 代人R向平衡方程后，可以求出，即 将上式代人（4．43）式的右面，就得到二维的Grad-Shafranov方程 （4.49） 其中的定义见（4.44）式．这时如果给定了和的具体函数形式后，加上合适的边界条件，就可以从GS方程解出．而=常数的空间曲面，就是相互叠套的环形磁面．例如，往往取 然后用数值方法解出。只有在个别情况下，如=常数，或时，G-S方程才有严格解（见文献[3]）． 3 G-S方程的小参量展开，解析解 当把作为小参量，对 G-S方程作微扰展开时，其最低阶（即阶）的方程相当于无环效应时、直柱螺旋箍缩的平衡方程；而阶方程及解则是相应的环修正．对普通tokamak及反场箍缩来说，它们展开的下一阶方程是完全一样的，但反场箍缩位形的推导过程比较简单，所以下面从这种位形来导出G-S方程（以及磁通）的小参量展开表达式．对反场箍缩，有下面关系： （对tokamak相应的量级关系则是）． 1．的展开 先看小参量展开的零阶，即直柱形螺旋箍缩的情况．假定等离子体直柱外面被一导体壳包围，则的等离子体表面是最外面的一个磁面，它是半径为a，中心在真空室小截面几何中心的圆柱面．而其内部的磁面是一个个同心的圆柱面（图2.20），其中半径为的磁面，其磁通可写作 在展开的下一阶，直柱已弯成了环，这时由于在R方向的扩张力使磁面都朝外挤压．最外的磁面由于和金属壳相重而无法移动，因此只有内磁面移向外侧，而且半径愈小的磁面移动得愈多．因为。故移动量也是一阶小量，所以这时环形磁面仍可以看成是原来半径为的那个磁面，只是其圆心移至 处，变成了一个偏心圆（图2.21），在原来的极坐标中其轨迹可以表示成（图2.22） 环形磁面中心的这个位移量被称为Shafranov位移． 这时相应的磁通就变成 其中是直柱近似下的零阶磁通，仅其自变量改成。而则是磁通的一阶修正值．在处对作展开，准确到的一阶，上式变成 因为和是同一个磁面的磁通值，所以大小一定相同，（注意 ）。这样从上式可得一阶磁通的修正值 由关系 得到 ：，这样 或得Shafranov位移 把上面对特定的讨论变成对任意的磁面的讨论，就可得到一般的表达式 2．Grad-Sha