HT7U快速反馈控制实时平衡重建算法简介-EAST.ppt

HT7U快速反馈控制实时平衡重建算法简介;主要内容;EFUND程序简介;开始;EFITPLOT的绘图原理;新的数据格式;软件功能;;SHAPE程序 主要是为固定边界计算提供输入的 GetSets 子程序主要是用于获得控制参量并读入 Green Table (ec6565.ddd: e_coil的格林函数 rfcoil.ddd:f_coil的格林函数 ep6565.ddd：plasma的格林函数) Data_input的主要功能 设置磁数据和权重矩阵;EFIT程序简介;描叙等离子体压力平衡的Grad-Shafranov平衡方程可以通过例如环向电流密度上的限制条件等有用的测量信号得到它的解。 这样就使可以用两个磁通值的函数和来表示,这里P和F分别是等离子体压强和角向电流;Jt表示等离子体环向电流密度分量。 在以上方程的基础上寻找一个使方程 最小的解从而得到 和 x是归一化的磁面函数;; 1假设初始等离子体电流分布 2 计算由在计算边界上产生的磁通值做为边界条件 3 求解Grad-Shafranov方程，求出 ，确定等离子体边界 4 求解最小二乘问题，确定参数 5计算等离子体电流密度分布 6 重复步骤2,3,4,5，直到满足一定的收敛条件;; EFIT程序的相应子程序 FIT: Equilibrium and fitting iterations. GREEN: 重新计算格林函数值. PRESUR: 重新计算压力剖面值. MATRIX:重新计算电流剖面值. CURRNT:计算等离子体电流值. FCURRT: 得到外部成形场线圈电流值. PFLUX: 得到等离子体磁通值. BOUND:得到等离子体边界. FINDAX:得到等离子体磁轴. RESIDU: 检察是否满足平衡. ;进展情况FIT;currnt computes the current density on the r-z mesh. currnt程序可以看出如果是平衡模式 + 多项式 的话就跳到1100; 如果是繁衍模式第一次跳到3100 执行GAQ电流模型 do 3300 i=1,nw do 3300 j=1,nh kk=(i-1)*nh+j pcurrt(kk)=0.0 if ((xpsi(kk).lt.0.0).or.(xpsi(kk).gt.1.0)) go to 3300 rdiml=rgrid(i)/rzero pp0 = (1.-xpsi(kk)**enp)**emp*(1.-gammap)+gammap pcurrt(kk)= rbetap/rdiml*((1.-xpsi(kk)**enf)**emf *(1.-gammaf)+gammaf) pcurrt(kk)=pcurrt(kk)+pp0*rdiml pcurrt(kk)=pcurrt(kk)*www(kk) tcurrt=tcurrt+pcurrt(kk) 3300 continue do 4050 i=1,kppcur brsp(nfcoil+i)=0.0 4050 continue do 4060 i=1,kffcur brsp(nbase+i)=0.0 4060 continue brsp(nfcoil+1)=cratio/rzero brsp(nbase+1)=cratio*rbetap*rzero brsp(nfcoil+2)=-brsp(nfcoil+1) brsp(nbase+2)=-brsp(nbase+1) 4100 continue return;fcurrt computes the currents in the f coils. 如果是固定极向场的时候if (ifcurr.gt.0) return 该程序不执行； 该程序首先从\green_table 目录下读入 fc6565.ddd 等文件中的值，然后改变 brsp(i) 变量的数值； 然后调用call solve(nfcoil,nfcoil,afma,brsp,ifmatr)。 solve一方面是采用高斯消元法分解线形代数方程组的矩