J-TEXT铁芯模型欧姆放电数值模拟.pdf
《J-TEXT铁芯模型欧姆放电数值模拟.pdf》由会员分享,可在线阅读,更多相关《J-TEXT铁芯模型欧姆放电数值模拟.pdf(8页珍藏版)》请在咨信网上搜索。
1、June2023NuclearFusionand Plasma Physics2023年6 月Vol.43,No.2第43 卷第2 期子体物理核聚变与等离文章编号:0 2 54-6 0 8 6(2 0 2 3)0 2-0 2 19-0 8D0I:10.16568/j.0254-6086.202302016J-TEXT铁芯模型欧姆放电数值模拟张顺1,桂腾,杨锦宏 2,陈凯杰陆野,宋强,江中和3,汪卫华2(1.安徽大学计算机科学与技术学院,合肥2 3 0 6 0 1;2.安徽大学物质科学与信息技术研究院,合肥2 3 0 6 0 1;3.华中科技大学电气与电子工程学院,武汉43 0 0 7 4)摘要
2、:本文采用有限元软件分析了J-TEXT的完整铁芯模型、等离子体电流丝模型和简化中心螺线管模型的磁场特征,确定了简化中心螺线管的几何参数与电流,之后把有限长度中心通电螺线管参数加到托卡马克模拟程序TSC模型中,模拟J-TEXT的欧姆放电过程,获得了欧姆放电电流、电子温度、等离子体密度、等离子体电流,磁轴位置等参数。模拟结果与实验结果数据相符,为具有铁芯结构的托卡马克放电模拟提供了有效手段关键词:J-TEXT;铁芯托卡马克;TSC;欧姆放电中图分类号:TL62+2文献标志码:A1引言J-TEXT装置是一个中型带铁芯常规磁体的托卡马克装置,其前身是TEXT-Ul,2 0 0 6 年在华中科技大学安装
3、调试成功 2 。在该装置内,靠近等离子体的铁芯改变了磁场拓扑结构,对磁场结构和等离子体产生显著影响,不能忽略。TSC是一个二维计算程序,它通过在矩形计算网格上解磁流体方程组来模拟自由边界等离子体随时间的演变过程 3 。用 TSC对 TFTR4,ADITYA5等装置的常规放电进行了模拟,模拟结果通常能与实验数据吻合较好,并且可以对等离子体破裂进行模拟 6 。另外,该程序还可对正在设计建造的装置进行模拟,给出如等离子体磁轴演变、电流、环压、温度、密度、压力、比压、反磁等与等离子体有关的随时间演变的海量数据信息 7 10 1,为将来的装置开展等离子体物理实验研究提供参考。研究铁芯对装置磁场的影响方法
4、有简化铁芯模型理论分析、计算机建模分析和实验测量等。袁保山等将把中心柱近似成无限长,研究了铁芯对极向场的影响。Solano等建立了卷轴模型,计算了铁芯在非饱和情况下对极向场的影响,并给出了极向场的解析表达式 12 。何勇等采用无限长和卷轴模型与J-TEXT装置在指定位置测量出的“dryshot条件下的数据进行对比(只有VF线圈放电),将卷轴模型代入空气芯托卡马克求解程序,得到铁芯托卡马克的平衡求解程序 13 。程际等将无限长和卷轴模型归类成二维模型,建立了以J-TEXT装置为原型的三维有限元模型,提出二维模型忽略了铁芯的环向不对称性,无限长模型只考虑中心柱的影响,卷轴模型假设磁导率无限大,将铁
5、芯抽象为封闭卷轴状且计算时不考虑外壁的影响。研究结果表明在净电流不为0 的情况下,两类模型计算的极向磁场相收稿日期:2 0 2 1-0 4-14;修订日期:2 0 2 2-12-2 5基金项目:国家自然科学基金(119 7 50 2 2);国家磁约束核聚变能发展研究专项(2 0 18 YFE0309101,2 0 17 Y FE0 3 0 0 6 0 3,2015GB120006);国家国际科技合作专项(2 0 15DFA61760)作者简介:张顺(19 8 2-),男,博士,副教授,从事数据差分隐私保护、移动计算优化算法等方面研究。*通讯作者:杨锦宏(19 7 7-),男,安徽桐城人,硕士研
6、究生导师,从事等离子体数值模拟研究。220第43 卷核聚变与等离子体物理差较大,具有明显的环向不对称性 14。杨锦宏等采用TSC模拟了不考虑J-TEXT装置铁芯时异常输运导致的破裂情况并与实验数据进行对比,根据模拟结果计算了面向等离子体的实验包层模块的动态响应,通过数值计算结果表明该模块的最大应力在安全范围内 15。综上所述,J-TEXT装置铁芯结构的研究比较丰富,但是J-TEXT铁芯托卡马克装置在模拟程序中进行欧姆放电研究还未开展本文采用有限元软件建立了J-TEXT装置完整铁芯模型、等离子体电流丝模型与中心螺线管模型,利用TSC代码模拟有铁芯的J-TEXT放电过程,为具有铁芯结构的托卡马克放
7、电模拟提供了有效手段。2TSC理论模型TSC程序是轴对称、自由边界、非线性的随时间演化的托卡马克模拟程序,程序在二维尺度上求解速度和磁通函数的分布演化,在一维尺度上求解磁面平均的绝热不变量的时间演化 16 。TSC 模型单流体电阻性磁流体运动方程为:(Om)/at+V.(1/p)mm=JB-Vp(1)式中,J为电流密度;B为磁场;p为压强;p为质量密度;m为动量密度,是位置和流函数的函数。TSC中时间尺度下平衡方程中惯性项相比于电磁力要小得多,因此程序中等价的增加等离子体质量,省略了对流微商项(1/p)mm,相应的增加了一个粘性阻力项F(m),以此来修正惯性项 17 。TSC中将等离子体中的粘
8、度力取为:F,(m)=-%Vm-VVm)2 V(V m)(1)式中,2 为可压缩粘滞系数和不可压缩体积粘滞系数,此时修正后的等离子体平衡方程为:(om)/at+F(m)=JxB-Vp(2)TSC模型中将等离子体之外的导体模拟为一组环向连续的电流丝 18 。每个电流丝的几何参数即是其所在计算网格的参数(x,zi,dx,d,)。对于某一个位于(x,z)处的电流丝,其电流满足回路方程:dt(4)式中,r,为电流丝电阻;I.,为电流丝电流;I,为目标电流丝外包括等离子体电流和其他导体电流在内的所有电流;V(t)为外加电压;M1.为电流丝互感,MA.jij=MoG(,),其中G(x,)是)一-,L,是电
9、流丝格林函数;Li,=ox8XlnAx的自感,且假设dx,=dz。穿过导体回路内部的磁通包括自感磁通和与其他导体的互感磁通,即:-2w.,=L,l,+Zi,M/l.j(5)导体电阻为:r(i)=(2 元x,nN(i.)/A(6)式中,n(i,j)是是等效电阻率,A=dx,=dz,是导体横截面积。根据VB=J,得到导体环向电流为:Ii.,=(AAAy)/(ox,)(7)导体极向磁通的演化:0y/at=(ny)/o+V(t)/2元(8)根据TSC导体理论模型,本文将中心螺线管简化为导体,并将其细分为数匝环向连续的电流丝,电流丝电流满足回路方程3有限元软件三维模型在有限元分析软件中建立的三维模型由P
10、F线圈、真空室、铁芯等实验装置主要组成部分构成,其中PF线圈系统是由加热场(OH)线圈、垂直场(VF)线圈、水平场(HF)线圈组成,如图1所示。OH线圈起到提供伏秒数、加热维持等离子体的作用,放电过程中可以通过OH线圈的电流来控制等离子体电流。VF 线圈主要在真空室内产生垂直方向的磁场,该磁场与等离子体电流相互作用,产生水平方向的作用力以维持等离子体在水平方向上的平衡,VF线圈由基于中平面上下对称的三组铜线圈组成,共有16 匝。J-TEXT装置采用带铁芯的结构,铁芯主要由中心柱、外柱和上下横梁组成,中心柱是直径为0.8 m、高为1.3 m的圆柱,材料是硅钢片,其与不带铁芯的托卡马克装置结构相比
11、,有2 个很221张顺等:J-TEXT铁芯模型欧姆放电数值模拟第2 期明显的优点,一是铁芯可以减小欧姆加热场的物理设计,二是可以使等离子体区域内的杂散磁场变得很小。但是铁芯本身的非线性特性会大大增加等离子体平衡计算的难度,其平衡情况在很大程度上将依赖于铁芯的饱和程度。真空室是托卡马克装置的另一个重要组成部分,它为高温等离子体提供真空环境,J-TEXT装置的真空室是由单层3 0 4不锈钢建造而成。a(5)4212019.18171615R=6.4mH-2.64m0Z=Om43)(1)中心柱(2)外柱(3)TF线圈(4)OH线圈(5)VF线圈(6)HF线圈(7)真空室图1J-TEXT装置主要组成模
12、块的右半部分截面图(a)及J-TEXT实验装置图(b)装置基于中平面Z=0上下对称;真空室中1-2 1为磁探针的标号,磁探针左右间距为6 0 mm,上下间距为12 6 mm;中心柱高度H=2.64m,半径R=0.4m。为了确定在TSC中使用的中心螺线管线圈的长度与匝数,本文通过三类三维模型来计算铁芯内以及真空室内磁感应强度以便进行对比,获得中心螺线管线圈模型的参数。图2 a为第一类模型,由三组共6 3 个观测点、完整的铁芯(其中心柱高度为1.3m),PF线圈以及真空室组成。图2 b为第二类模型,该模型通过3 3 组电流丝模拟了等离子体电流,通过对电流丝添加相应的电流密度可以模拟出某时刻的等离子
13、体电流值;图2 c为第三类模型,该模型铁芯无外柱及上下横梁仅有中心螺线管,其他组成部分与第一类模型相同。本文通过计算对比铁芯模型、等离子体模型和中心螺线管三种模型的中心点处磁感应强度、真空室内6 3 个观测点处极向磁感应强度以及实验数据,确定中心螺线管线圈的长度和匝数等参数。aCPointasmaoint.43-63poHnt中心螺线管图2 三类三维模型J-TEXT装置铁芯(蓝色模块)、PF线圈(透明模块)、真空室(绿色模块)尺寸相同的三维模型;b等离子体电流丝(紫色plasma模块)模型;a-高度1.8 m的中心螺线管(紫色模块)模型。图a中真空室内等离子1-2 1号观测点与铁芯垂直,2 2
14、-42 与43-6 3 号观测点和铁芯平行且与Z轴对称。4有限元模型分析结果4.1铁芯中心点处磁感应强度等离子体上升阶段的伏秒数主要由OH场线圈提供,铁芯能够改变真空室中的杂散磁场,但是铁芯改变磁场强度和铁芯饱和程度相关,因此将等离子体上升阶段末,平顶阶段初始时刻的OH场线圈222第43 卷核聚变与等离子体物理电流作为ANSYS模型中OH线圈输入参数,通过调节铁芯不同相对磁导率来计算出铁芯接近饱和时中心点处的磁感应强度大小。当铁芯相对磁导率为1时,铁芯相当于真空区域,此时铁芯中磁场很小,如图3 所示。在相对磁导率1 10 0 时,随着铁芯相对磁导率的增加,铁芯磁场越来越大,铁芯也逐渐饱和,在相
15、对磁导率300500时,虽然相对磁导率增加了2 0 0,但是铁芯磁场处于饱和状态,磁场增长比较缓慢,增长率约7%。因此将铁芯相对磁导率50 0 时中心点处磁场强度作为基准值。21.510.50050100150200250300350400450500u图3铁芯不同相对磁导率时中心点处磁场大小u铁芯相对磁导率;B-铁芯中心点处磁场强度。OH线圈共有40 匝,每匝线圈电流为4kA,OH线圈总安匝数约16 0 kA匝。在中心螺线管模型中,将等离子体电流设置为140 kA、16 0 k A 两组算例,仅求解只有等离子体电流、铁芯相对磁导率为500时铁芯中心点处磁感应强度。由图4可知,当仅有等离子体电
16、流为140 kA和16 0 kA时,铁芯中心点处磁感应强度均不超过1.6 T,这与图3 中完整铁芯模型计算数据相符。1.61.41.210.80.60.40.2140145150155160lasma/kAplasma图4仅有等离子体电流时,铁芯中心点处感应磁场lplasma等离子体电流;B铁芯中心点处极向磁场强度。第三种模型是将中心螺线管代替完整铁芯,中心螺线管厚度为0.0 4m,中间由空气体填充,因此中心螺线管中心点处磁感应强度与中心螺线管线圈匝数相关。本文将线圈密度用每米中心螺线管含有线圈匝数表示,这些线圈电流都与OH线圈相同。由图5可知,中心螺线管高度从1.3 1.8 m,电流单位长度
17、匝数为150 nml、3 0 0 n ml、3 50 n m时,中心螺线管中心点处在相同线圈密度磁场强度相近,由此可知,中心点处磁场强度与中心螺线管高度相关性小。21.5上10.501TVVVIVIIX图5不同高度和线圈密度中心螺线管中心点处磁场-1.3m-300nm-l;II-1.3m-350nml;III1.5m-300nm;IV-1.6m-250nm*;V-1.6m-300nml;VI1.6m-350nml;VI-1.8m-250nml;VI-1.8m-270nm-l;IX-1.8m-300nm1。横坐标表示不同高度、不同线圈密度(nm表示为每米含有OH线圈电流的匝数)中心螺线管;B一中
18、心点处极向磁场强度。随着电流单位长度匝数的增加,在相同的中心螺线管高度下,中心点处磁场均呈线性增加趋势,在2 50 3 0 0 nm-l之间时,中心点处磁场与前两种模型较为接近,由此将中心螺线管设置为250300nm-1之间的值。4.2真空室内观测点处极向磁场随着铁芯相对磁导率的增加,铁芯对极向磁场影响越来越大,则真空室内由OH线圈提供的极向磁场也会随着铁芯相对磁导率的增加而相对减少。第一种J-TEXT模型中设置了三组观测点。如图2所示,覆盖了整个真空室以便测量等离子体区域的极向磁场。当铁芯相对磁导率为1时,相当于没有铁芯存在,此时,等离子体区域极向磁场最大,且分布较为规律,如图6 所示。第一
19、组观测点位置分布垂直于铁芯,受到铁芯影响相对其他两组较弱,磁感应223张顺等:J-TEXT铁芯模型欧姆放电数值模拟第2 期强度相对较大,体现了铁芯对磁场的影响具有环向不对称性。随着铁芯相对磁导率的增加,真空室中极向磁场逐渐减小,在铁芯相对磁导率设置为3 0 0 时,为接近于0 的状态,当相对磁导率设置为50 0 时,虽然磁场也有所减少,但此时真空室中极向磁场已经很弱,铁芯已经吸收了绝大部分磁场,相较于相对磁导率3 0 0 时几乎没有变化。如图7 所示,将中心螺线管线圈密度设置为2 50 3 0 0 nml,通过观测点测量不同中心螺线管高度时真空室的磁场强度。0.030.02Hu-50point
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- TEXT 模型 欧姆 放电 数值 模拟
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。