EFIT反演计算研究教学课件.pptVIP

EFIT反演计算研究 钱金平 EFIT结构 EFIT主程序（平衡反演计算） EFOUND程序 （互感计算） Shape程序（生成边界） 原理 磁轴、边界 基本参数 平衡计算典型位形 反演：Matrix (AX＝B) 辅助工作 efit程序在windows下运行 嵌入matfor软件，运行结果可以直接图形显示。 编写了mfile生成文件。 信号比较 原因：同一个探针对不同线圈的灵敏度不一致。 信号筛选 误差设定小于10％（＋1，否则－1）（对单一f-coil) FLUX-LOOP Flux信号选取 靠近内侧的单匝环 反演 信号：探针，单匝环，等离子电流 迭代条件： 1）matrix中 errorm=errmin(errmin=0.01) 且abs(saisq-saiold)=0.1 2）residu中 if (errorm.le.error) idone=1 问题一：闭合磁面(bound) 确定等离子体边界的条件如下： ① ② ③ nerr=3 return 以上第三个条件（删除） 问题二 ：开环及不收敛 当以设定的psi值寻找等磁面时，可能会出现此情况 (cntour) (icntour=1（删除）) 超过设定的25次迭代次数，不满足收敛条件（大部分） 问题三：真空室电流 1）带入AX=B计算？（可能现有真空室段数较多） 2）不考虑涡流 3）Vloop/Rv ? 现在以2）处理 反演收敛条件 (abs(saisq-saiold).le.0.01)以此为判断条件 （errorm=errmin(errmin=0.01) （cancel） 同时 residu子程序中的判断条件被跳过 2466＃1.1s /1.2s（单匝环内侧） F-coil电流比较 Shot:2466#1.1s(Flux-loop归一化) F-coil电流比较 M-file 1)单匝环在原先内侧信号基础上，删除了个别与反演偏差较大的信号 选取15个信号 2）小探针在原先以4为标注的基础上降到2选取32个信号 3）电流权重降为1 2466(选取较多探针信号）1.1s/1.2s 2481＃1s～1.6s(mfile取较多探针） 结论 1）想得到精度高的反演结果（满足errorm.le.error收敛条件），可能会比较困难。 2）从f-coil电流上看，归一化后计算处的电流值与测量值同样偏差比较大。 3）真空室电流需要考虑，尤其在VDEs 。 4）即使按照现在的限制条件，还有一部份无法得到收敛解。 展望 1）信号处理、选择、权重等可能还需要进一步细化。 2）真空室电流处理还需要进一步讨论。 谢谢大家 * EAST ASIPP EAST ASIPP 原理 A X ＝ B EAST ASIPP 中心坐标： ， 小半径： 上拉长比： 下拉长比： 上三角形变： 下三角形变： EAST ASIPP EAST ASIPP fulx loops magnetic probes current F coils P’ FF’ F coils P’ FF’ EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP EAST ASIPP *