求解有相对论效应的Fokker Plank方程的程序包的开发与应用 Numerical.ppt

 龙永兴 董家齐 焦一鸣 石秉仁 王爱科 张锦华 郦文忠 （核工业西南物理研究院） # 国家自然科学重点基金支持 程序包的功能 求解没有相对论效应的Fokker-Plank方程的程序包已于2002年完成[1,2]。 该程序包具有相对论效应，是反跳(bounce averaged)平均的，考虑了捕获电子效应和各种不同加热模式（LHCD ECRH）下的准线性扩散系数[3-5] 。 程序包由意大利引进，可用于模拟等离子体中粒子的动力学过程，适用于磁约束环形装置等离子体辅助加热与电流驱动的动力学计算,对HL-1M和HL-2A的实验数据分析非常有用。 程序包的功能 我们将程序包成功移植，增加了计算功能，验证了程序的可靠性，并利用开发的编码程序对低杂波在托卡马克等离子体中的吸收和驱动效率进行了数值模拟，数值结果与理论分析完全一致，详细分析见文[8]。 程序包的功能 计算结果包括： 1. (不同坐标系)分布函数 2.密度,能量密度 3. 各向异性参数 4.加热功率 5.电流驱动效率 6.高能成份,能量传输等值. 参 考 文 献 1.龙永兴 石秉仁 董家齐等,“Fokker-Planck程序包的开发与利用”《科学技术与工程》Vol.3No5,p257(2003). 2.石秉仁 龙永兴 董家齐等,“Fokker-Planck study of tokamak electron cyclotron resonance heating”，《中国物理Chinese Physics》,vol.12,no.11,(2003) 3．D. Farina, IFP Internal Report, FP 99/14, 2000. 4．A. Simonetto et. IFP Internal Report, FP 00/08, 2000. 5． A. Cardinali, Recent research Developments in Plasmas, Trasnworld Research Network Magazine, Vol. 1 185 (2000). 6．M. Shoucri, I. Shkarofsky, Comp. Phys. Comm. 82 287 (1994) 7．Y. Peysson, M. Shoucri, Comp. Phys. Comm. 109 55 (1998) 8. 焦一鸣 龙永兴 董家齐 石秉仁 高庆弟,”捕获电子效应对低杂波电流驱动的影响”,2003年GF报告，《物理学报》2005,n0.1待发表。 Fokker-Planck方程  Fokker-Planck方程 速度空间总流量是: 表示为两项: Fokker-Planck方程 波加热项只是扩散. DC电场只是磨擦力. 碰撞项包含了扩散与磨擦力。 Fokker-Planck方程 Fokker-Planck方程 准线性扩散系数 LHCD： 准线性扩散系数 ECRH: 数值方法 利用有限差分方法求解Fokker-Planck方程。空间差分近似为中心差分，时间轴采用Crank-Nicholson差分格式。 对差分方程采用5点（ 9点）强隐式方法,这是一个高精度算法[6,7]。 Fokker-Planck系数及其导数通过分布函数和Rosenbluth势的Legendre系列展开求得。 F-P方程(code) F-P方程的离散化 差分格式 差分方程 九点差分方程 五点差分方程 输入参数 N1: 在 P方向的离散点数. N2: cos方向的离散点数，u(2)=-1,u(N2-1)=1. DT: 时间步长. PRINT: 结果输出频率. ITS: 计算的时间的总步数。 TRAP: =1 包含了捕获电子效应. TRAP: =0 没有捕获电子效应. 输入参数 使用数值方法的选取： MODE: =5 使用五点格式. =9 使用九点格式. APARAM: 方法收敛参数 (通常取 60). 输入参数 PMAX: P 的最大值. PMIN: P 的最小值，通常=0 TRUNC: =1 包含E-E碰撞项（动量守恒） =0 没有E-E碰撞项 EEPS： Ohmic 电场参数 ? 输入参数 Ti: