EAST纯欧姆放电的数值模拟-EAST纯欧姆放电的数值模拟.pdfVIP

维普资讯 第27卷 第4期 核 聚 变 与 等 离 子 体 物 理 Vo1．27。NO．4 2007年 l2月 NuclearFusionandPlasmaPhysics Dec． 20o7 文章编号：0254-6086(2007)04--0280-06 EAST纯欧姆放 电的数值模拟 刘成岳 1,2，肖炳甲 ，吴 斌 ，刘连忠 ，罗正平． (1．中国科学院等离子体物理研究所，合肥 230031；2．合肥工业大学理学院，合肥 230009) 摘 要：利用 l，2一D平衡演化程序对EAST的首轮纯欧姆放电进行了数值模拟，获得的模拟结果与实验结果 吻合得很好。同时，对放电全过程中的伏秒消耗进行了分析。 关键词：EAST；l’，2D平衡演化程序；伏秒；数值模拟 中图分类号：O51l．9 文献标识码 ：A 1 引言 区、真空区和导体区三个区。等离子体区与真空区 EAST是一台具有非圆截面的全超导托卡马克 的分界面可以随时问变化而变化，它由第—个接触 装置。实验 目标是要实现长脉冲稳态运行，以及在 到限制器或含 X．point的等磁通面确定。在每个区 稳态运行条WF改善约束，从而为未来先进聚变反 内分别求解经过一系列处理的MHD—Maxwell方程 应堆奠定工程技术和物理基础。同时其上开展的工 组，这些处理包括引入了等离子体质量放大因子和 作对未来在法国建成的 ITER也具有重要的参考价 粘滞放大因子、真空区电导等参数。由于这些特殊 值，这使得它的建成在聚变领域具有重要的意义。 的处理，方程组在穿越区域边界时形式并不发生改 本文使用 Il／2一D平衡演化程序对EAST首轮纯 变，其解仍保持连续，便于进行数值求解。这里主 欧姆限制器位形的放电进行了模拟，该程序已用于 要介绍 Il／2_D平衡演化程序的特点，详细物理模型 PBX、TFTR、DIII．D等托卡马克的放电模拟及BPX、 的数学表述见文献 【l，2]。本程序中包括：(1)单独的 KSTAR、ITER等装置的设计工作。 电子能量方程和离子能量方程；(2)新经典电阻率； (3)自举 电流 ；(4)时间平均 的锯齿模型 ；(5) 2 物理模型Il，z】 Coppi—Tang输运模型；(6)只考虑完全电离离子的辐 射。初始的自由边界等离子平衡是在—个矩形网格 I1／2一D平衡演化程序是一个 自由边界轴对称托 内形成，由极向场电流提供边界条件。等离子体电 卡马克等离子体随时间演变及其相关控制系统的 子温度和离子温度是通过求解 1，2一D面平均输运方 数值模拟程序，它通过在矩形计算网格上解磁流体 程获得。在粒子输运的模拟中，等离子体密度具有 方程组来模拟 自由边界等离子体随时间的演变过 程。MHD方程细j萤过边界条件和外部极向场线圈中 如下的剖面分布形式： (，f)=，l)Il一 l 的电流相互耦合，同时在磁面坐标中求解压力和温 +，l(f)。其中， 为归一化的磁通坐标，在磁轴处 度的面平均输运方程。程序考虑了反常输运、新经 为0，等离子体边界处为 1， (t)和 0(t)分别为 典电阻率、自举电流、辅助加热、电流驱动、G【粒 等离子体在边界及中心处的密度。在EAST第2344 =加热、气泡模稳定性等多种物理模型，它可以用