托卡马克聚变等离子体参数演化的数值模拟
王时佳;王少杰
【摘 要】基于一维粒子-能量输运方程组,采用差分方法编写了模拟托卡马克聚变等离子体的密度和温度剖面演化的数值求解程序.在其中使用预估-校正方法处理了方程中的非线性项并利用能量约束时间定标律求得输运系数.模拟得到的聚变性能与之前文献采用其他方法计算的结果符合较好.该程序可用于研究托卡马克反应堆的聚变性能以及温度和密度的剖面形状.%The evolution of plasma temperature and density profiles is simulated with a code which solving the one-dimensional transport formulas with difference method.The nonlinear terms in the formulas are treated with a predictor-corrector method and the transport coefficient is evaluated by matching the scaling law of energy confinement time.The fusion performance calculated with the modeling results accord well with the results calculated with other models in the references.The code in this paper can be ud to study the fusion performance and the profiles of temperature and density in tokamak fusion reactors.
【期刊名称】《电子科技》
顾思妤
【年(卷),期】2017(030)006
【总页数】4页(P5-8)
【关键词】数值模拟;预估校正算法;非线性方程;定标律
【作 者】王时佳;王少杰
【作者单位】中国电子科技集团公司第20研究所 雷达事业部, 陕西 西安 710068;中国科学技术大学 近代物理系,安徽 合肥 230026
【正文语种】中 文
【中图分类】TN01;O241.82
国际热核聚变实验反应堆(ITER)将成为人类历史上第一个能长时间维持聚变增益因子Q(聚变释放功率与装置消耗的辅助加热功率之比)大于10的聚变实验装置[1-3]。为了达到这样的高性能状态，ITER中的反应燃料必须成为燃烧等离子体状态，即装置中的等离子体主要由聚变反应所释放的能量进行加热。研究表明聚变反应的发生率与燃料的密度呈非线性正相
关关系，而与温度呈非线性非单调关系[4]。因此在燃烧等离子体中，温度的演化将是高度非线性的。
作为一种托卡马克类型的磁约束聚变实验装置，ITER的主要设计参数是由现有的托卡马克装置实验运行参数进行统计和外推而得到的[2,5]。ITER 要实现设计目标中的高参数运行仍具有较大的不确定性和挑战性。因此，在实验装置建造之前进行可靠的数值模拟验证工作仍具有重要的价值。目前主流的关于ITER 聚变性能计算的模拟工作主要使用粒子输运和能量输运相互独立的模型来进行[3]，对于密度和温度演化的自洽性仍有欠缺。
本文从一维柱位形下的等离子体流体输运方程出发，使用能量约束定标律建立了一种约化的等离子体输运模型。其中利用预估-校正算法处理了能量输运方程中的非线性源项，模拟了等离子体密度剖面和温度剖面的自洽演化，计算了反应堆的聚变性能。
基于电中性等离子体条件以及等温等离子体假设(电子和离子具有相同的温度T)，燃烧等离子体的输运过程可由以下一维柱坐标输运方程组描述[6]
其中，n是电子密度；qr是热传导导致的径向热流密度；Γr是径向粒子流密度；ur是粒子的
径向对流速度；D和χ分别是粒子与热量的输运系数；Sp代表粒子源，Sα和Saux分别是氘氚聚变反应和辅助加热设备带来的能量源项；Srad代表因辐射带走能量而导致的能量损失项；r代表等离子体小半径方向上的位置。
与传统的能量输运方程不同[7-8]，燃烧等离子体中的物理过程是高度非线性的：一方面等离子体主要由核聚变反应生成的阿尔法粒子碰撞加热，而局域内的阿尔法粒子生成率是密度以及温度的非线性函数；另一方面在等离子体能量损失中占据相当比重的辐射损失的单位体积发射强度也是温度和密度的非线性函数。具体表现在输运方程中，有
其中，k是由氘氚燃料密度与电子密度相对比例决定的系数；<σv>是氘氚反应的碰撞截面，根据文献[4]，<σv>=3.68×10-12T-2/3exp(-19.94T-1/3)cm3c-1(T的单位为keV)；Eα=3.5 MeV是聚变产生的阿尔法粒子携带的能量；SBrm，SSyn和SL分别代表韧致辐射，回旋辐射和线辐射带来的能量损失；CB，CS，CL和A是由装置参数决定的常数[9]。
在此说明对于输运方程组的数值求解方法。记[0，T]为模拟所用的时间区间，将其分为M等分，则有Δt=T/M为时间步长，tm=mΔt，m=0,1,…,M；将径向空间分成J等分，空间步长为Δr，rj=jΔr，j=0,1,…,J。记函数f(r,t)在网格点(rj,tm)处的值为对于输运方程组中的时间偏导
sat阅读分数表采用隐式差分，对扩散项和对流项采用中心差分。假定式(1)中粒子源项分布稳定，Sp只是径向位置的函数。对于能量输运方程(2)中的源项采用预估- 校正方法，具体做法为：(1)预估：将Tm带入Sh并使用Δt/2作为推进步长，推进方程(2)，求得Tm+1/2；(2)校正：再将Tm+1/2带入Sh并使用Δt作为推进步长，推进方程(2)，求得Tm+1。在模拟中每一次Tm→Tm+1的演化都使用上述步骤。
在以往的输运模拟工作中，计算rj处物理量从tm时刻演化到tm+1时刻所需的输运系数和一般都是根据已有的输运理论由当时当地的物理量(温度，密度等)计算得出。这种“第一性原理”方法虽然很符合经典的物理学逻辑，但由于实际托卡马克等离子体的复杂性以及现有等离子体物理理论本身的限制，计算出的结果往往与实验结果偏差很大。
本文使用了一种约化方法来求出输运系数：在任意时刻根据ITER能量约束时间定标律，由系统的各宏观物理量计算出此时的能量约束时间并据此计算出下一时刻的系统中等离子体的总能量；固定输运系数的剖面形状，而对于其绝对大小数值做猜测值并带入输运方程组，计算出相应的下一时刻等离子体总能量，将其同根据能量约束时间计算出的数值进行比较并据此修正输运系数大小的猜测值，再将新的输运系数大小猜测值带入输运方程组并
迭代计算猜测值，直到输运系数对应的下一时刻等离子体总能量与根据定标律计算得到的结果相符，即得到欲求得的输运系数。具体来说，假设热输运系数有以下形式
其中，g(r)为给定的输运系数剖面形状因子；χ0为待求数值。χ0在tm时刻到tm+1时刻演化时所用的数值通过迭代求解数值方程来确定，符合要求的χ0可使得tm时刻和tm+1时刻之间等离子体总能量W的变化满足
其中，P=(Sα+Saux-Srad)d3x是等离子体的有效总加热功率；W=3nTd2x；τE为能量约束时间。目前ITER工作组给出的能量约束定标律形式为[2]
其中，HH98(y,2)为约束改善因子；Ip是等离子体电流；BT是环向磁场强度；是电子线平均密度；R为等离子体大半径；M是平均离子质量数；κα为装置的拉长比；ε为装置的逆环径比。
edc是什么
本文设定粒子输运系数D与χ有以下关系
这是一个被广泛采用的近似，并被许多实验验证过[10-11]。
根据上述数值方法，这里给出在ITER典型长脉冲放电参数下，对于燃烧等离子体演化过程的模拟过程和结果。
vere比较级
对于输运系数的剖面形状，这里假设(10)式中的输运系数形状因子g(r)有如图1中所示的剖面，注意到图中的输运系数形状具有典型的外部输运垒(ETB)特征，输运垒的宽度为0.05 个等离子体小半径。
对于方程(1)中的粒子源项，考虑到在反应堆托卡马克中粒子加料只能渗透进入很浅的区域，本文默认粒子源项Sp形式为Sp=Cexp[]。其中r0=0.95a，dr=0.1a，常数C的数值根据输出的平均密度调整，以使得稳定状态下平均密度达到ITER的预期密度( 1×1020 m-3)。
isad对于方程(2)中的辅助加热项Saux，假定其形状为Saux=Saux,0exp[-(r/raux)2]，其中raux=0.4a。Saux,0由公式Sauxd3x=Paux计算，根据文献[12]，本文取Paux=40 mW。
边界条件设置为
在设置了一系列初值条件测试后我们发现，最后的稳态剖面对于初值并不敏感，这与预期的由燃烧等离子体自主演化决定剖面的性质相符。这里设置温度初始分布T0(r)和密度分布
初值n0(r)如图2所示。
宗教笑话对于能量定标律中涉及到的装置宏观物理量，这里取值与文献[5]相同，其具体取值参见表1。对于径向粒子流密度中包含的粒子对流速度ur，假设为粒子箍缩速度，根据文献[13],其大小和分布为
其中，q=q(r)为托卡马克的安全因子，这里取q为ITER感应放电运行下的典型形状[12]，其大小和形状如图6所示。
woodstock由以上条件，计算得到了稳定的等离子体各参数剖面，其形状如图3～图7所示。由稳定剖面得到的聚变增益因子Q=10.84，与之前文献中通过其他模型计算得到的ITER性能基本相同[5,12,14-15]。
对于托卡马克聚变装置来说，其聚变性能主要依赖于密度剖面和温度剖面的维持。本文采用数值求解等离子体粒子和能量输运方程组的方法模拟了ITER参数下托卡马克聚变反应堆的聚变性能。其中考虑到聚变反应加热源项的非线性，对其采用了预估-校正方法处理。模拟中需要使用的输运系数通过能量约束定标律进行了计算。模拟得到的结果与之前文献用
其他方法得到的结果符合较好，可以用来评价聚变反应堆性能。此外，该结果还可以用于进一步研究加料和辅助加热对于聚变反应堆性能的影响。
【相关文献】
tramp是什么意思
[1] ITER Physics Expert Group.Chapter 1: Overview and summary[J].Nucl Fusion,1999,39(12):2137-2174.
customerrvice
[2] Shimada M,Campbell D J,Mukhovatov V,et al.Chapter 1:Overview and summary[J].Nucl Fusion, 2007,47(6):S1-S17.
