1、EAST放电模拟及自举电流参数分析,吴 斌2004.9.14,主要内容,TSC程序的物理模型EAST装置放电模拟EAST自举电流参数分析,平衡程序与放电模拟程序的区别,平衡程序EQ EFIT某一时刻的等离子体位形放电模拟程序TSC等离子体位形的时间演化过程,放电模拟的不同途径,1.5维平衡演化程序二维平衡程序一维输运程序TSC程序Tokamak Simulation Code引入处理过的惯性项和粘滞项不再求解平衡方程转而求解等离子体运动方程,TSC发展概况,S-1 - 球马克位形的形成过程PBX - 位形控制、等离子体破裂、伏秒数分析TCV - 灵活等离子体形状控制系统的设计 CIT - 伏秒
2、数分析、等离子体破裂DIIID - 位形控制、 等离子体破裂 VDE 、伏秒数分析BPX - 燃烧控制、偏滤器分离磁面摆动TPX - 位形控制、垂直位移控制、等离子体破裂 ITER - 伏秒数分析、位形控制、等离子体破裂TFTR - 伏秒数分析、杂质分析NSTX - 1MA球环伏秒数分析,TSC物理模型,TSC程序是PPPL开发的一个自由边界轴对称托卡马克等离子体随时间的演变以及和它相关控制系统的数值模拟程序它通过在矩形计算网格上解磁流体方程组来模拟自由边界等离子体随时间的演变过程，MHD方程组通过边界条件和外部极向场线圈中的电流相互耦合，同时在磁面坐标中求解压强和密度面平均输运方程程序考虑了
3、反常输运、新经典电阻率、辅助加热、电流驱动、 自举电流、 粒子加热、气泡模稳定性等许多物理模型，可以用于模拟控制系统、辅助加热、电流驱动等问题。,Surface-Averaged Transport Equations,TSC主要方程,围绕TSC程序开展的工作,TSC程序的移植和开发三个版本Pub5、Pub7、Unx10.2位形演变过程的放电模拟模拟了EAST上升段和平顶段等离子体由园截面到所预定的拉长截面，由限制器位形到偏滤器位形的建立、演变、达到稳态的过程垂直位移控制系统的数值模拟探索了TSC程序进行垂直位移稳定性分析的方法模拟了EAST垂直位移系统的控制响应过程,EAST放电模拟研究,E
4、AST(HT-7U)装置特点放电模拟方案的设计放电模拟结果及分析,托卡马克磁体系统组成,环向场系统产生强的环向磁场,与等离子体电流产生的较弱极向磁场形成具有旋转变换的磁场位形极向场系统加热场 变压器原理 感应等离子体电流 焦耳热平衡场 抵消等离子体向外扩张成形场 多极场 形成拉长截面反馈场 配合被动致稳部件 控制等离子体的垂直,1, 2, 3 多极场线圈4 真空室 5 环向场线圈 6 - 欧姆线圈 7 等离子体,HL-2A 极向场系统,EAST超导托卡马克装置,EAST采用一体化设计方案的 极向场系统极向场电流是以上各种场电流等的叠加极向场系统的运行和控制变得更为复杂长脉冲及稳态运行考虑输运过
5、程对等离子体行为的影响将磁流体力学方程组和输运方程组结合起来进行放电过程的动态平衡演化模拟,放电模拟方案的设计,装置几何参数的确定 （二维模型）放电程序场的设计 PF 波形反馈控制系统的设计PID控制等离子体位形控制等离子体电流控制,EAST装置基本几何参数,EAST放电程序场的设计,放电程序场波形,加热场+平衡场,等离子体形状反馈的模型,观察对 (x1 x2)PID控制极向场线圈,EAST的等离子体形状反馈设计,随时间变更观察对合适的PID参数作用于PF2,PF3,PF4和PF6 七套PID控制系统,EAST等离子体形状反馈参数设置,Input File,EAST放电模拟结果（磁面演变）,电
6、流分布演变,等离子体位形演变,上升段时间为4秒的放电模拟,EAST自举电流参数分析,固定边界平衡程序JSOLVEREAST的六种运行模式不同运行模式的比较EAST的自举电流计算和分析,JSOLVER程序,JSOLVER,li,q,IBS,Standard Tokamak是一种标准的高值、大电流、第一稳定区的运行方式，其中心q值q0 1，等离子体压强和电流分布是一般变压器感应驱动托卡马克典型的剖面分布，这种运行模式的特点是高值、低P值、P4.30,ARIES-I是在ARIES反应堆研究计划中的一种运行模式，它仍处于第一稳定区，是一种在高值和高P值之间的折衷方案。等离子体电流减小，它的中心q值增加
7、到1.3左右，使q* /q0上升到3附近。这种模式使仍处于第一稳定区的托卡马克可以获得很高的自举电流份额和较好的自举电流分布。,Supershot 最早是在TFTR装置上发现的一种约束改善的运行模式，它的特点是具有较低的等离子体电流和拉长比，有非常峰化的等离子体电流和压强分布，仍处于在第一稳定区。在TFTR上使用限制器降低了粒子的再循环，从而进入supershot运行模式，约束改善因子可以达到3.5,具有最高的值和自举电流份额，通过改变等离子体电流分布使在中心区域q的分布出现可以反剪切dq/d 2和非常峰化的等离子体电流和压强分布而可以穿越第一稳定区进入第二稳定区。由于峰化程度很高，在等离子体
8、边界附近，等离子体的电流密度、压强和压强梯度都几乎为零。,不同运行模式的比较,EAST中的自举电流,自举电流的研究对于托卡马克的稳态运行具有非常重要的意义自举电流的计算过程-首先求解固定边界平衡方程，得到磁面以及有关物理量和它们的磁面平均值，再利用有关自举电流理论计算自举电流的磁面平均量JBbs的分布，最后计算自举电流占全部等离子体电流的份额。等离子体的密度、温度、压强等分布参数对自举电流的大小和分布有较大的影响，通过改变这些分布参数，可以改变自举电流的大小和分布形状。,等离子体 P 和中心温度对自举电流大小影响较大 P越高，自举电流越大。，中心温度越高，自举电流份额越大。,等离子体的压强分布
9、对自举电流的分布形状有较大的影响 峰化的压强分布使自举电流的分布的峰值靠近等离子体中心 平坦的分布使自举电流的分布的峰值靠近等离子体边缘,等离子体的密度分布影响自举电流的大小 密度分布越峰化，自举电流越大,反剪切位形下自举电流的优化,总结,建立和逐步完善了以TSC程序为中心的放电模拟软件。利用TSC程序模拟了ESAT等离子体由园截面到所预定的拉长截面；由限制器位形到偏滤器位形的建立、演变、达到稳态的过程，并将X点稳定在所需的位置研究了EAST装置未来可能采用的标准、反剪切模式等六种不同的运行模式，对这些运行模式约束性质进行了比较分析。对EAST的自举电流进行了研究，系统计算和分析了影响自举电流大小和分布的各种等离子体参数与自举电流份额的关系。,研究计划,考虑低杂波驱动的TSC放电模拟工作，连接LSC程序和TSC程序开展垂直位移事件（VDE）的数值模拟，模拟由VDE事件继而导致等离子体破裂的放电过程，研究Halo电流的形成，探索破裂的防护与预防措施。,谢 谢,
