EFIT 反演计算研究

EAST ASIPP EFIT反演计算研究钱金平

EAST ASIPP EFIT结构 • EFIT主程序（平衡反演计算） • EFOUND程序（互感计算） • Shape程序（生成边界）

EAST ASIPP 原理原理磁轴、边界 A X ＝ B

EAST ASIPP 基本参数，中心坐标：小半径：上拉长比：下拉长比：上三角形变：下三角形变：

EAST ASIPP 平衡计算典型位形

EAST ASIPP F coils P’ FF’ fulx loops F coils magnetic probes current P’ FF’ 反演：Matrix (AX＝B)

EAST ASIPP 辅助工作 • efit程序在windows下运行 • 嵌入matfor软件，运行结果可以直接图形显示。 • 编写了mfile生成文件。

原因：同一个探针对不同线圈的灵敏度不一致。原因：同一个探针对不同线圈的灵敏度不一致。 EAST ASIPP 信号比较

EAST ASIPP 信号筛选 • 误差设定小于10％（＋1，否则－1）（对单一f-coil)

EAST ASIPP FLUX-LOOP • Flux信号选取 • 靠近内侧的单匝环

EAST ASIPP 反演 • 信号：探针，单匝环，等离子电流 • 迭代条件： 1）matrix中 errorm<=errmin(errmin=0.01) 且abs(saisq-saiold)<=0.1 2）residu中 if (errorm.le.error) idone=1

EAST ASIPP 问题一：闭合磁面(bound) • 确定等离子体边界的条件如下： • ① • ② • ③ nerr=3 return • 以上第三个条件（删除）

EAST ASIPP 问题二：开环及不收敛 • 当以设定的psi值寻找等磁面时，可能会出现此情况 (cntour) (icntour=1（删除）) • 超过设定的25次迭代次数，不满足收敛条件（大部分）

EAST ASIPP 问题三：真空室电流 • 1）带入AX=B计算？（可能现有真空室段数较多） • 2）不考虑涡流 • 3）Vloop/Rv ? • 现在以2）处理

EAST ASIPP 反演收敛条件 • (abs(saisq-saiold).le.0.01)以此为判断条件 • （errorm<=errmin(errmin=0.01) （cancel） • 同时 residu子程序中的判断条件被跳过

EAST ASIPP 2466＃1.1s /1.2s（单匝环内侧）

EAST ASIPP F-coil电流比较

EAST ASIPP Shot:2466#1.1s(Flux-loop归一化)

EAST ASIPP F-coil电流比较

EAST ASIPP M-file • 1)单匝环在原先内侧信号基础上，删除了个别与反演偏差较大的信号选取15个信号 • 2）小探针在原先以4为标注的基础上降到2选取32个信号 • 3）电流权重降为1

EAST ASIPP 2466(选取较多探针信号）1.1s/1.2s

EAST ASIPP 2481＃1s～1.6s(mfile取较多探针）

EAST ASIPP 结论 • 1）想得到精度高的反演结果（满足errorm.le.error收敛条件），可能会比较困难。 • 2）从f-coil电流上看，归一化后计算处的电流值与测量值同样偏差比较大。 • 3）真空室电流需要考虑，尤其在VDEs 。 • 4）即使按照现在的限制条件，还有一部份无法得到收敛解。

EAST ASIPP 展望 • 1）信号处理、选择、权重等可能还需要进一步细化。 • 2）真空室电流处理还需要进一步讨论。

EAST ASIPP 谢谢大家

EFIT 反演计算研究