固定边界平衡方程的数值求解-南华大学学报.pdf

第28卷第3期 南华大学学报(自然科学版) Vol郾28 No郾3 2014年9月 Journal of University of South China(Science and Technology) Sep郾2014 文章编号:1673-0062(2014)03-0001-04 固定边界平衡方程的数值求解 黄千红 (南华大学 核科学技术学院,湖南 衡阳421001) 摘摇 要:采用基于能量原理的变分法求解等离子体平衡方程,在给定的压强和环向电 流密度分布情况下,采用变步长默森法对平衡方程进行了数值求解,得到了D 型截 面下极向和环向磁场分布.采用NSTX装置参数计算所得结果与文献结果一致.开发 的固定边界平衡方程求解程序具有通用性特点,可以模拟不同装置的平衡位型. 关键词:Grad鄄Shafranov方程;变分法;默森法 中图分类号:TL631摇 摇 摇 文献标识码:A The Solution of Equilibrium Equation with Fixed Boundary HUANG Qian鄄hong (School of Nuclear Science and Technology,University of South China,Hengyang,Hunan421001,China) Abstract:A variational method based on an energy principle is applied to solve plasma e鄄 quilibrium equation.When the plasma pressure and current density profile is given,the balance equation is solved using a variable鄄step鄄size Merson method.The poloidal and to鄄 roidal magnetic field with dee鄄shaped plasma are obtained. The results which adopted NSTX device parameters are achieved agreeable to those in the references.The program is universal,which can be applied to simulate the balance equilibrium configuration in the different devices. key words:Grad鄄Shafranov equation;variational method;Merson method 磁场把一定密度的高温等离子体稳定的约束足够 0摇 引摇 言 长时间,达到劳逊条件,进而产生净能量增益.等 磁约束受控核聚变的一个关键问题,是利用 离子体平衡是指磁压强与等离子体热压强的平 收稿日期:2014-04-25 基金项目:研究生创新课题基金资助项目(2013XCX02);衡阳市科技局基础课题基金资助项目(2013KJ24);湖南省 核聚变与等离子体物理创新团队建设基金资助项目(NHXTD03) 作者简介:黄千红(1977-),男,湖南湘潭人,南华大学核科学技术学院讲师,博士研究生.主要研究方向:核聚变与 等离子体. 2