我国科学院等离子体物理研究所博士毕业论文答辩.ppt

中国科学院等离子体物理研究所博士论文答辩 博士论文题目 HT-7U托卡马克偏滤器数值模拟研究 研究生姓名： 陈一平 指导导师： 虞清泉 研究员 朱思铮 研究员 报告内容 HT-7U边界等离子体和偏滤器数值模拟研究意义 边界等离子体流体模拟程序B2.5的移植 HT-7U高再循环与低再循环偏滤器的数值模拟研究 HT-7U辐射偏滤器的数值模拟研究 流体模型与动力学模型混合模拟 B2.5-EIRENE连接与耦合 HT-7U偏滤器B2.5-EIRENE数值模拟研究 论文结论和讨论 发表的论文 致谢 HT-7U托卡马克偏滤器数值模拟研究意义国际上边界等离子体和偏滤器数值模拟程序的发展 边界等离子体流体模拟程序B2.5的移植 B2.5程序所利用的物理模型 二维等离子体数值模拟程序 模拟多流体等离子体 多流体等离子体中包括中性气体 考虑原子分子过程，包括电离、复合、电荷交换等微观物理过程 考虑不同粒子种类之间的相互摩擦 流体方程中包括了电流和电势方程并利用了欧姆定律 模拟程序包包括266子程序，共20000余条语句 边界等离子体流体方程 边界等离子体流体方程 边界等离子体流体方程在准柱坐标系中的表达 边界等离子体流体方程在准柱坐标系中的表达 流体方程在准柱坐标系中的表达 边界等离子体流体模拟程序B2.5的移植 B2.5程序流体方程 边界等离子体流体模拟程序B2.5的移植 B2.5程序流体方程 边界等离子体流体模拟程序B2.5的移植 B2.5程序流体方程 边界等离子体流体模拟程序B2.5的移植 B2.5程序对流体方程的求解 流体方程 猜一套解 计算误差 离散化效正方程 求解线性代数方程 矩阵 由强隐方法求逆 边界等离子体流体模拟程序B2.5的移植 B2.5程序流体方程求解边界条件 B2.5边界等离子体数值模拟程序是世界三大边界等离子体流体模拟程序(B2.5, EDGE2D, UEDGE)之一。 将B2.5 UNIX 工作站版本移植为微机版本，利用FORTRAN POWER STATION 软件系统。充分利用国内较为流行的微机资源在边界等离子体数值模拟研究中的作用。 将B2.5应用于HT-7U托卡马克偏滤器的数值模拟研究。 移植动力学Monte-Carlo模拟程序EIRENE,将B2.5耦合到边界等离子体动力学Monte-Carlo模拟程序EIRENE (世界三大边界等离子体动力学Monte-Carlo模拟程序, EIRENE, DEGAS, NIMBUS, 之一），形成微机版本的耦合程序B2.5 -EIRENE。 将B2.5 -EIRENE应用于HT-7U偏滤器的数值模拟。 HT-7U高再循环与低再循环偏滤器数值模拟研究模拟参数和几何 HT-7U高再循环与低再循环偏滤器数值模拟研究模拟参数和几何 HT-7U高再循环与低再循环偏滤器数值模拟研究偏滤器靶板附近数值模拟网格的产生 HT-7U高再循环与低再循环偏滤器数值模拟研究模拟计算网格在笛卡尔坐标中的表示 HT-7U高再循环与低再循环偏滤器数值模拟研究模拟计算边界条件 1.芯-SOL界面：所有种类粒子的平行流速度设为零，正如径向电流。两个温度，电子温度和离子温度设置为共同值，两个密度，中性氢H0密度和氢离子H+密度分别设置为固定值。两个温度和两个密度被设置作为控制参数。 2.外壁和隐闭流区的内边界：热损失包括流损失和附加损失，这个附加损失由边界上的梯度引起，与这个边界条件相对应的边界梯度长度的量级约为0.01米。边界上有很低的温度。朝着边界的带电粒子损失流可预设，这个损失流正比于密度与热速度的乘积。零电流，平行速度的零剪切。粒子再循环系数和能量再循环系数分别等于1.0 和0.3。 HT-7U高再循环与低再循环偏滤器数值模拟研究模拟计算边界条件 3.赤道面假设为对称性边界条件。当考虑漂移流时，模拟结果将做适当修正。 4.偏滤器靶板：所有带电粒子以共同的速度外流，中性粒子的局部产生速率等于局部带电粒子外流的速率。电流由Bohm鞘条件决定但同时考虑二次电子发射。这些流平行于磁场方向，但是为了得到到达靶板上的流，这些流必须考虑磁场的局部倾斜。粒子再循环系数和能量再循环系数分别等于1.0 和0.3。 HT-7U高再循环与低再循环偏滤器数值模拟研究高再循环偏滤器数值模拟结果 热流分布 HT-7U高再循环与低再循环偏滤器数值模拟研究高再循环偏滤器数值模拟结果 偏滤器靶板上总热流分布 HT-7U高再循环与低再循环偏滤器数值模拟研究高再循环偏滤器数值模拟结果 外壁上的