论文部分内容阅读
托卡马克中随机磁场会导致等离子体的反常输运,从而降低托卡马克的约束性能。回旋动理学数值模拟是研究湍流输运很有用的工具,模拟可以给出粒子流和能量流。一般来说这些流包括扩散项和对流项,区分每一项的贡献对理解湍流输运是重要的,但现有的模拟很难区分这两种流。利用相空间输运矩阵计算流的方法可以区分这两种流。本文基于李变换扰动方法研究了电磁扰动下的回旋动理学二阶哈密顿量,改进了导心轨道程序GYCAVA;然后基于GYCAVA程序发展了一个通过积分粒子导心的全轨道数值计算相空间输运矩阵的程序。 第一个工作中,我们考虑低频电磁扰动,用两步李变换方法详细计算了回旋动理论的二阶哈密顿量,指出长波近似下,磁扰动δA对应的回旋动理论二阶哈密顿量主要是1/2δA2‖;这是因为回旋中心变换导致的回旋中心坐标下的磁矩与导心坐标下的磁矩存在差别Δμ,这一差别对应的哈密顿量ΔμB0抵消了磁扰动垂直部分的能量1/2δA2⊥,其中B0为平衡磁场。基于该理论,我们将导心轨道程序GYCAVA加以改进,并用来模拟TFTR纵向波纹磁场中高能粒子的导心轨道及其随机输运,得到了ripple场中高能粒子随机输运的区域与GWB理论符合。 第二个工作中,我们基于改进的导心轨道程序GYCAVA发展了一个程序,该程序通过积分导心全轨道数值计算了相空间输运矩阵。我们用两个算例来验证数值计算输运矩阵的正确性:一个是计算随机磁场中无碰撞的电子扩散系数De与电子平行速度v‖和扰动场幅度|δA|的关系De∝|v‖|;De∝|δA|2;另一个是计算随机磁场中存在均匀径向电场时,电子导心径向位置的变化δr和平行速度的变化δv‖0间的交叉关联Drv‖,以及δv‖的自关联Dv‖v‖,数值计算了Drv‖,Dv‖v‖与Drr间的关系,以及Drr和Dv‖v‖与平行速度的依赖关系。数值结果与理论相符。最后我们数值计算了拉格朗日关联函数,得到托卡马克随机磁场的关联时间和平行关联长度,结果表明随机磁场的平行关联长度可以用qR0来估计,q和R0分别表示安全因子和大半径。