CFETR参数下α粒子慢化过程的数值模拟.pdf
《CFETR参数下α粒子慢化过程的数值模拟.pdf》由会员分享,可在线阅读,更多相关《CFETR参数下α粒子慢化过程的数值模拟.pdf(9页珍藏版)》请在咨信网上搜索。
1、专题:磁约束等离子体中的高能量粒子CFETR 参数下 粒子慢化过程的数值模拟*吴相凤王丰林展宏陈罗玉于召客吴凯邦王正汹(大连理工大学物理学院,大连116024)(2023年 4月 29 日收到;2023年 6月 8 日收到修改稿)氘氚聚变产生的高能量 粒子是维持未来托卡马克反应堆等离子体高温的主要加热源,良好的 粒子约束对于维持稳态燃烧等离子体至关重要.在持续发生聚变反应的系统中,粒子远离热平衡,呈现非麦克斯韦分布.如果忽略轨道效应,基于局域库仑碰撞的假设可以得到 粒子的经典慢化分布,然而由于 粒子存在较大的漂移轨道宽度,空间输运不容忽视,为得到更为准确的 粒子分布函数,需要开展相关的数值计算
2、.本文使用模拟程序 PTC(particletracercode)在中国聚变工程试验堆(CFETR)不同的放电模式下,采用粒子轨道跟踪和蒙特卡罗碰撞方法,对 粒子慢化过程进行了数值模拟,获得了更为真实的 粒子分布函数,并将其与经典慢化分布进行了对比.结果显示分布函数在中等能量附近和经典慢化分布存在较大差异.进一步的分析表明,这是由于中等能量下 粒子的较强的径向输运引起的.本文的研究结果对准确评估 粒子加热背景等离子体的能力具有重要参考价值.关键词:托卡马克,粒子,慢化分布PACS:52.50.Gj,52.55.Fa,52.55.Pi,52.65.yDOI:10.7498/aps.72.2023
3、07001引言氘氚聚变反应是目前最有希望实现商业化应用的聚变反应之一1,2.随着国际热核聚变实验堆(ITER)工程的持续推进,人们开始关注在ITER 之后的聚变装置的设计和建造问题,以便尽早地实现聚变能源商业化3,4,中国聚变工程试验堆(CFETR)将作为 ITER 和聚变电站商用示范堆(DEMO)之间的关键衔接点,以推进中国聚变能源发展,推动聚变能源商业化进程,以及实现可持续地清洁能源方面发挥重要作用5,6.200 MW在 CFETR 燃烧等离子体中,因聚变生成的携带有 聚变功率的 粒子是加热等离子体的主要方式之一,然而,这些高能量的 粒子只有被有效慢化才能起到维持等离子体自持燃烧的作用71
4、0.在氘氚等离子体中关于 3.5MeV 的 粒子在等离子体中的慢化分布或能量沉积的研究已经进行了半个多世纪1115.粒子慢化过程中,在非均匀磁场中漂移并受到电子和离子的摩擦以及与背景离子的角度散射,基于局域库仑碰撞假设求解福克-普朗克方程的稳态解可以得到经典慢化分布16,17,它被广泛应用于包括 粒子、中性束注入粒子等在内的快离子相关物理研究,如加热、电流驱动、快离子与 Alfvn 本征模的相互作用及快离子输运等1821.然而这种慢化分布具有一定局限性,在真实情况下,轨道效应和空间输运是不可忽略的,同时,粒子漂移轨道宽度也会影响 粒子在等离子体中的输运2022.这种简化模型无法准确描述 粒子的
5、行为,有必要采用更丰富的物理模型描述非局域的 粒子慢化过程23,24.为了探索 粒子的慢化过程,本文回顾了经*国家磁约束核聚变能发展研究专项(批准号:2022YFE03090000)、国家自然科学基(批准号:11975068)和大连理工大学基本科研业务费(批准号:DUT22LK18)资助的课题.通信作者.E-mail:2023中国物理学会ChinesePhysicalSocietyhttp:/物理学报ActaPhys.Sin.Vol.72,No.21(2023)215209215209-1典慢化分布理论模型并采用数值模拟方法进行对比分析.在数值模拟方面使用粒子轨道追踪耦合蒙特卡罗模拟程序 PT
6、C(particletracercode)25在CFETR 稳态运行模式和混杂运行模式下对 粒子的慢化过程进行了计算,模型考虑了粒子的有限轨道效应以及新经典输运,采用蒙特卡罗方法计算 粒子与背景等离子体之间的库仑碰撞,得到了等离子体平衡条件下 粒子的慢化分布函数.同时,统计了 粒子热化和损失的比例,计算了 粒子对背景等离子体的加热功率,并得到了 粒子的总体粒子约束时间和能量约束时间.本研究通过在不同运行模式下对 粒子整个慢化过程中的 粒子的产生、约束、对背景等离子体的加热以及整个分布演化的数值模拟分析,还有理论和程序模拟的对比分析,为 CFETR 的物理设计提供了有价值的参考依据.本文第 2
7、节描述了经典慢化分布理论模型、CFETR稳态和混杂运行模式下的等离子体参数、PTC 程序物理模型.第 3 节是两种运行模式下 粒子慢化过程的模拟结果,包括 粒子的粒子数变化、对背景等离子体的加热功率、粒子的能量变化、粒子的慢化分布函数等.第 4 节是讨论分析,包括两种运行模式下慢化过程的分析、约束时间的分析、理论与程序模拟的慢化分布函数的分析.第 5节是总结.2经典慢化分布与 PTC 模型2.1 经典慢化分布经典慢化分布是指高能粒子与热背景等离子体碰撞导致高能粒子的经典减速形成的速度或能量分布,它可以作为描述 粒子慢化过程的一种理论模型.经典慢化分布的推导基于忽略速度扩散项和轨道效应的假设,根
8、据这个假设,可以简化福克-普朗克方程,得到如下式子16:ft=Cf+S,(1)fSvi v vevivevTi=Te其中 为 粒子速度分布函数,为 粒子源项.此外,假设 ,其中 和 为背景离子和背景电子的热速度,为典型的 粒子速度,并假设背景电子温度与背景离子温度相等,即,有高能 粒子的慢化碰撞算子:Cf s1v2(/v)(v3c+v3)f,(2)s=3(2)3/220mTe3/2Z2e4m1/2eneln0mTeZemenelnvc=(3 me4mZ)1/3vt,eZvt,e其中 Spitzer 慢化时间 ;为真空介电常数;为 粒子质量;为背景电子温度;为 粒子核电荷数;为单位电荷量;为电子
9、质量;为电子密度;为碰撞参数;临界慢化速度 ,为有效电荷数,为电子热速度.接着,基于源项是各向同性的假设,同时 粒子都以初始能量 3.5MeV 诞生,形成一个以初始速度为中心的 函数,给出源项:S=S0(vb v)4 v2,(3)S0=nDnTvnDnTvvb=1.3 107m/s其中聚变反应速率 ,为氘的粒子数密度,为氚的粒子数密度,为反应截面,为相对速度大小.对于能量为 3.5MeV 的 粒子,初始速度 .将(2)式和(3)式代入(1)式并进行计算,得到稳态时三维的 粒子速度慢化分布函数:f3D(v)=S0s4H(vb v)v3+v3c.(4)H(vb v)4 v2其中 为阶跃函数,将(4
10、)式乘以 ,得到一维的 粒子速度慢化分布:f(v)=S0sv2H(vb v)v3+v3c.(5)Ev根据能量 与速度 之间的变换关系以及阶跃函数的定义,得到一维的 粒子能量慢化分布:f(E)=S0s2E1/2H(Eb E)E3/2+E3/2c.(6)EbEc其中,为 粒子初始能量,为临界能量,根据(6)式,可得到 CFETR 稳态运行模式下的局域稳态能量慢化分布如图 1 所示.接下来,利用各向同性的 粒子慢化分布计算出 粒子的密度 n 和总能量 nE,对分布函数(4)在三维空间积分,得到n=S0s3ln1+(Eb/Ec)3/2,(7)nE=S0sEc16ln(1 Eb/Ec+Eb/Ec(1+E
11、b/Ec)2)13tan1(Eb/Ec 13)13tan1(13)+S0sEb2.(8)物理学报ActaPhys.Sin.Vol.72,No.21(2023)215209215209-2根据粒子平衡方程和能量平衡方程:nt=S0np,(9)(nE)t=S0EbnEE,(10)pEnt=0(nE)t=0其中 和 为 粒子的粒子约束时间和能量约束时间,稳态时由 ,可计算得到p=s3ln(1+(Eb/Ec)3/2),(11)E=sEcEb16ln(1 Eb/Ec+Eb/Ec(1+Eb/Ec)2)13tan1(Eb/Ec 13)13tan1(13)+s2.(12)接下来的模拟中,将使用 PTC 程序对
12、 粒子的慢化过程进行模拟,并与理论模型进行比较分析.2.2 CFETR 参数与 PTC 物理模型基于 CFETR 稳态运行模式和混杂运行模式(v201903)参数下进行了 粒子慢化过程的模拟,具体等离子体参数如图 2 所示.由图 2 可看到,在稳态运行参数下,芯部区域和边界区域等离子体参数相对比较均匀,混杂运行模式下的温度略高于稳态模式,密度量级与稳态参数一致.CFETR 稳态运行模式是一种基于等离子体稳定性和安全性的非感应电流驱动的运行模式,混杂运行模式是感应和非感应电流混杂驱动的脉冲式的运行模式.这两种运行模式不是独立的,它们可以相互转换和组合,以实现更优秀的聚变反应效果和能量输出.1.6
13、1.2/(1014 m-3SeV-1)0.80.031.002.003.00/MeV4.000.40123f1,f2,f3图1电子温度分别为 27.78,14.4 和 6.7keV,对应电子密度分别为 1.141020,9.341019和 7.471019m3参数下得到的经典能量慢化分布 Fig.1.Classicalenergyslowing-downdistributionsf1,f2andf3obtainedfortheelectrontemperaturesof27.78,14.4and6.7keV,and their corresponding electron densities
14、of1.141020,9.341019and7.471019m3.061218/keV243000.30.60.91.2/(1020 m-3)00.20.40.60.81.0246810(a)eD=TeD=T061218/keV243000.30.60.91.2/(1020 m-3)00.20.40.60.81.0246810(b)eD=TeD=T图2CFETR 中的背景等离子体参数(a)稳态运行模式下的密度、温度和安全因子剖面;(b)混杂运行模式下的密度、温度和安全因子剖面Fig.2.Background plasma profiles in CFETR:(a)Density,tempera
15、ture,and safety factor profiles in steady-state scenario;(b)density,temperatureandsafetyfactorprofilesinhybridscenario.物理学报ActaPhys.Sin.Vol.72,No.21(2023)215209215209-3PTC 程序是一个拥有自主知识产权的粒子追踪程序,可以模拟聚变反应中产生的 粒子在等离子体中的慢化过程25.PTC 程序可以利用全轨道模型和漂移轨道模型来追踪粒子运动,在 粒子慢化过程模拟中,由于全轨道计算较为耗时,而且根据 粒子的能量和位置估算出 粒子的回旋半径
16、在 0.55cm,回旋半径对慢化过程的影响是可以忽略不计的,因此,可以采用漂移轨道的导心运动方程追踪 粒子.在模拟过程中,粒子源由单位体积的聚变速率来定义25,粒子运动方程采用经典的四阶龙格库塔方法求解,粒子与背景等离子的碰撞通过蒙特卡罗方法求解.PTC 程序将慢化过程看作一系列弹性散射,在每个弹性散射中,PTC 先根据 粒子的初始速度和位置,利用含时演化的轨道方程计算 粒子的运动轨迹和所处的空间位置;然后,根据该位置区域内的背景电子和离子的密度、温度,计算出散射截面;接下来,使用概率分布函数来模拟碰撞过程中的能量和动量转移.同时,在每个小时间段内会随机采样一定数量的 粒子进行统计,记录 粒子
17、的位置、能量等信息,通过大量的随机采样和弹性碰撞模拟,就可以得到 粒子慢化过程中的分布函数.相较于经典慢化分布理论模型,PTC 程序在轨道运动的基础上,考虑了粒子的新经典输运过程,包括考虑径向扩散和粒子漂移轨道宽度等因素,从而能够更全面、准确地模拟 粒子的慢化过程.3模拟结果3.1 粒子慢化过程首先,在 PTC 程序的输入模块中输入 CFETR稳态运行模式下的平衡磁场和热等离子体剖面、扰动电磁场以及 粒子初始分布;然后在极向截面划分三角形网格,将背景等离子体信息插值到网格,载入 粒子源;接着,使用漂移轨道模型追踪粒子,使用蒙特卡罗方法求解粒子与背景等离子体的弹性碰撞;在模拟过程中,程序会统计
18、粒子数量、计算每个 粒子的能量损失率、根据碰撞截面计算 粒子对背景的加热功率.同样利用 CFETR混杂运行模式下的平衡磁场、等离子体剖面等参数模拟 粒子的慢化过程,得到两种模式下慢化过程中的物理量变化对比如图 3 所示.heating power/MW(c)02040608010012000.20.40.60.81.0/s_e+_D+_T+00.20.40.60.81.010211020101910181017101610151014(a)/s_total_thermal_lost/MeV(d)1.52.02.53.03.500.20.40.60.81.0/s0.20.30.4Lost rat
19、e/%0.50.60.7(b)00.20.40.60.81.0/sSteady-stateHybrid图3CFETR 稳态运行模式(实线)和混杂运行模式(虚线)下的各个物理量随时间的变化(a)粒子数量;(b)粒子损失率;(c)粒子对背景等离子体的加热功率;(d)粒子平均能量Fig.3.TimeevolutionofvariousphysicalquantitiesinCFETRsteady-statescenario(solidlines)andhybridscenario(dashedlines):(a)Numberofparticles;(b)lossrateofparticles;(c)
20、heatingpowerofparticlestothebackgroundplasma;(d)aver-ageenergyofparticles.物理学报ActaPhys.Sin.Vol.72,No.21(2023)215209215209-41.795 10201.709 10200.41%0.38%由图 3 可以看到,在两种运行模式下,被统计的在等离子体中的 粒子数量、热化的等离子体数量、损失的 粒子数量变化趋势几乎一致,分别在 0.68s 和 0.73s 时,稳态运行模式和混杂运行模式下在等离子体中的 粒子数量趋于平稳,达到稳态时的 粒子数量在稳态运行模式和混杂运行模式下分别为 和 ,
21、稳态运行模式下的数量略高于混杂运行模式.稳态运行模式下 粒子损失率最后维持在 ,混杂运行模式下 粒子损失率最后维持在 .稳态运行模式和混杂运行模式下对背景等离子体的加热功率分别在 0.57s 和 0.59s 时开始保持不变,稳态运行模式下 粒子对背景电子的加热功率达到110MW,对背景氘离子和氚离子的加热功率分别为 45.7MW 和 31.1MW,整体略高于混杂运行模式结果.稳态运行模式和混杂运行模式下 粒子的平均能量大约经过 0.64s 和 0.72s 由初始的3.5MW 趋于稳定的 1.61MW,即 粒子经过慢化过程达到稳态.此外,也在归一化极向磁通 空间对 粒子的加热功率密度进行了模拟计
22、算,结果如图 4 所示,根据图 4 发现,在两种运行模式下 粒子对背景电子和背景离子的加热功率密度量级一致,在不同的 位置有所不同.根据粒子追踪结果也得到了在两种运行模式下达到稳态时 粒子密度分布和在能量空间、极向磁通空间的分布,如图 5 和图 6所示.00.20.40.60.81.000.51.01.52.02.53.0 heating power density/(104 MWSm-3)e-_Steady-stateion+_Steady-statee-_Hybridon+_Hybrid图4 空间的加热功率密度Fig.4.Heatingpowerdensityinthe space.6(a
23、)420/m-2-4/m9876543210T10176(b)420/m-2-456789105678910/mT10181.351.201.050.900.750.600.450.300.150图5稳态时 粒子的密度分布(a)CFETR 稳态运行模式;(b)CFETR 混杂运行模式Fig.5.Theparticledensityinsteady-state:(a)CFETRsteady-statescenario;(b)CFETRHybridscenario.00.51.0normalized1.52.02.5(a)0.031.002.003.004.00/MeVSteady-stateHy
24、brid00.3normalized0.60.91.2(b)00.20.40.60.81.0Steady-stateHybrid图6 粒子分布函数(a)能量空间;(b)归一化极向磁通空间Fig.6.Theparticledistributionfunction:(a)Energyspace;(b)normalizedpoloidalmagneticfluxspace.物理学报ActaPhys.Sin.Vol.72,No.21(2023)215209215209-53.2 约束时间的计算3.5 10203.24 10201.795 10201.709 10200.51 s0.52 s0.24 s
25、0.22 s根据模拟结果,由粒子平衡方程和能量平衡方程可以得到程序模拟得到总体的 粒子约束时间和能量约束时间,粒子源项每秒钟产生的总粒子数在稳态运行模式和混杂运行模式下分别为 和 ,稳态时在等离子体中的粒子总数在稳态运行模式和混杂运行模式下分别为 和 ,计算得到 粒子总体粒子约束时间分别为 和 .同样,根据 粒子初始总能量和稳态时总能量计算得到 粒子的总体能量约束时间在稳态运行模式和混杂运行模式下分别为 和 .4讨论4.1 慢化分布对比分析E=0.22 MeV 3.5 MeV当 粒子在非均匀磁场中运动时,粒子在轨道上的漂移会导致粒子的输运方向发生变化,对于高能量 粒子,由轨道效应引起的输运是可
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- CFETR 参数 粒子 慢化 过程 数值 模拟
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。