托卡马克中的静电漂移不稳定性
托卡马克中的静电漂移不稳定性
E ⊥ B
练习题:
1,什么是反磁漂移?试证明反磁漂移对 密度扰动无贡献并讨论其物理意义。 2,对典型的托卡马克参数, kθ ∼ (0.1 − 1) ρ , 假设 计算电子漂移波频率和离子漂移波频率。
Ti ∼ Te ∼ 10kev, B ∼ 3T , Ln ∼ 50cm
Ⅲ、离子温度梯度不稳定性的流体研究 1,基本方程: 离子连续性方程:
4)线性理论的作用:饱和湍流幅的计算需要包括非 线性效应, 线性理论可以: i)确定可能的驱动机制,抑制机制; ii)确定不稳定性条件; iii)当湍流引起的输运占主导时,等离子体的密度和 温度梯度可能被调整到接近由线性不稳定性理论所预 言的阈值; iv)线性模的时间和空间特征可能与湍流态有一定的 联系,从而可以提供对湍流输运的粗略估算: 准线性理论 γ 混合长度计算: x ~ 2
2, Charged particle motion in tokamaks
• Parallel (lognitudinal) motion: v|| ≈ vth • Rotation: v⊥ = Ωρ ≈ vth • Drifts of guiding center i) Electric field drift:
k
Ⅱ、漂移波和漂移型不稳定性
1,磁化等离子体中的(反磁)漂移运动: 1) 漂移速度(均匀磁场):
Vdj =
− c ∇ Pj × B n jq j B
2
=
ˆ cT j b × ∇ n j q jn j B
2)粒子流强:
cT ∂n ∂n V⊥ Γ = ⋅ ⋅V⊥ = ∂x Ω qB ∂x cT ∂ n Vd = qBn ∂ x
练习题: 1,证明反磁漂移对密度扰动无贡献,讨论其物理 意义。 2,对典型的托卡马克参数,
托卡马克(Tokamak)核聚变
受控热核聚变能的研究主要有两种--惯性约束核聚变和磁约束核聚变。前者利用超高强度的激光在极短的时间内辐照氘氚靶来实现聚变,后者则利用强磁场可很好地约束带电粒子的特性,将氘氚气体约束在一个特殊的磁容器中并加热至数亿摄氏度高温,实现聚变反应。
托卡马克(Tokamak)是前苏联科学家于20世纪50年代发明的环形磁约束受控核聚变实验装置。经过近半个世纪的努力,在托卡马克上产生聚变能的科学可行性已被证实,但相关结果都是以短脉冲形式产生的,与实际反应堆的连续运行有较大距离。超导技术成功地应用于产生托卡马克强磁场的线圈上,是受控热核聚变能研究的一个重大突破。
法国研究的托卡马克核聚变装置托卡马克核聚变研究举步维艰,根本原因是轻元素原子核的聚合远比重元素原子核的分裂困难。原子核之间的吸引力是很大的,但原子核都带正电,又互相排斥,只有当两个原子核之间的距离非常接近,大约相距只有万亿分之三毫米时,它们的吸引力才大于静电斥力,两个原子核才可能聚合到一起同时放出巨大的能量。因此,首先必须使聚变物质处于等离子状态,让它们的原子核完全裸露出来。然而,两个带正电的原子核越互相接近,它们之间的静电斥力也越大。只有当带正电的原子核达到足够高的动能时,这需要几千万甚至几亿摄氏度的高温,它们的碰撞才有机会使它们非常接近,以致产生聚合。
实现梦想需要科学。经过多年大量科学实验证明在一种称为托卡马克核聚变能开发出无限而清洁的聚变能,它能帮助人类实现寻求能源的梦想。
托卡马克核聚变 - 基本原理
托卡马克核聚变原理演示核能是能源家族的新成员,包括裂变能和聚变能两种主要形式。裂变能 是重金属元素的核子通过裂变而释放的巨大能量。受控核裂变技术的发展已使裂变能的应用实现了商用化,如核(裂变)电站。裂变需要的铀等重金属元素在地球上含量稀少,而且常规裂变反应堆会产生放射性较强的核废料,这些因素限制了裂变能的发展。聚变能是两个较轻的原子核聚合为一个较重的原子核并释放出的能量。目前开展的受控核聚变研究正是致力于实现聚变能的和平利用。其实,人类已经实现了氘氚核聚变--氢弹爆炸,但那是不可控制的瞬间能量释放,人类更需要受控核聚变。维系聚变的燃料是氢的同位素氘和氚,氘在地球的海水中有极其丰富的蕴藏量。经测算,l升海水所含氘产生的聚变能等同于300升汽油所释放的能量。海水中氘的储量可使人类使用几十亿年。特别的,聚变产生的废料为氦气,是清洁和安全的。因此,聚变能是一种无限的、清洁的、安全的新能源。这就是世界各国尤其是发达国家不遗余力竞相研究、开发聚变能的根本原因。
托卡马克芯部等离子体微观不稳定性研究
托卡马克芯部等离子体微观不稳定性研究王亮;龚学余;李新霞;邹任忠【摘要】在Weiland模型的基础上,考虑碰撞效应的双流体方程,获得了托卡马克等离子体芯部区域ITG(离子温度梯度模)和TEM(捕获电子模)的色散关系,分析了归一化温度梯度和密度梯度以及归一化径向坐标对这两种模式增长率的影响。
计算结果表明等离子体密度梯度对离子温度梯度模有致稳作用,而对捕获电子模模有促进作用。
两种模的增长率都随各自温度梯度和归一化径向坐标的增大而不同程度增大。
%With the dispersion relation of the ion temperature gradient mode( ITG) and the trapped electron mode( TEM) derived from the Weiland model,where the two fluid equa-tions with the collision term are used,the effects of the normalized temperature gradient, the normalized density gradient and the normalized radial coordiante on the growth rates of these two modes are analyzed. The computation results indicate that the ion temperature gradient modes are stablized by the plasma density gradient and the traped electron mode is destabilized by the density gradient. The growth rate of both modes increase with plasma temperature gradient and the radial coordinate.【期刊名称】《南华大学学报(自然科学版)》【年(卷),期】2014(000)002【总页数】6页(P1-5,10)【关键词】芯部等离子体;增长率;离子温度梯度模;俘获电子模【作者】王亮;龚学余;李新霞;邹任忠【作者单位】南华大学核科学技术学院,湖南衡阳421001;南华大学核科学技术学院,湖南衡阳421001;南华大学核科学技术学院,湖南衡阳421001;南华大学核科学技术学院,湖南衡阳421001【正文语种】中文【中图分类】TL610 引言在托卡马克聚变装置中,等离子体芯部输运研究是目前的热点之一.等离子体芯部输运主要由各种各样的微观不稳定性[1-2]主导,而这些微观不稳定性所驱动的湍流导致等离子体出现了异常大的输运系数,严重制约着聚变装置的约束性能.等离子芯部湍流理论研究表明,离子温度梯度模(ITG),俘获电子模(TEM)和电子温度梯度模(ETG)这三种静电不稳定性模[3]所驱动.当存在径向川流(streamers)的情况时,相比于 ITG和TEM所驱动起来的湍流尺度,ETG的贡献可忽略不计[4],本文主要研究了ITG和TEM对芯部湍流的影响.关于芯部ITG和TEM的理论研究,近年来取得了积极进展.Dimits[5]和 Candy [6]分别于 2000年和2003年对ITG进行了研究,结果表明,在电子绝热近似下,当离子温度梯度满足|R▽Ti|/Ti>(R/LTi)(其中R表示托卡马克大半径,Ti表示离子温度,LTi=▽Ti/Ti表示归一化临界温度梯度)时,ITG 才能被激发.2005 年,Jenko[7]基于回旋动理学方法,在忽略碰撞影响的情况下,认为TEM由电子温度梯度和电子密度梯度共同驱动.此外,Weialnd从双流体模型出发,用 MMM95 CODE计算了kyρsH时(ky表示垂直磁场方程波数,ρsH表示氢离子拉莫尔回旋半径)ITG和TEM对输运的影响,并结合BALDUR CODE预测了TFTR,DⅢ-D,JET等聚变装置中的等离子体剖面,结果与实验吻合较好[8].本文在Weiland模型的基础上,进一步考虑了电子—离子间的库伦碰撞对ITG和TEM增长率的影响,并针对HL-2A装置芯部等离子体计算了0<kyρsH<2.0范围内这两种不稳定性模的增长率,分别分析了归一化温度梯度gT,密度梯度gn 和归一化径向坐标ρ对这两种模增长率的影响.1 物理模型1.1 基本方程采用双流体模型.对离子,连续性方程(粒子数守恒)为平行磁场方向动量方程为能量方程为其中,mH,nH,TH分别表示氢离子的质量数,密度、温度和抗磁漂移热通量,=2.5nHTH(× ∇TH)/eB2 ,速度+垂直磁场方向的速度V⊥H主要由电场、压力梯度、惯性力和有限拉莫效应引起的粘滞力这四种效应有关,可表示为[9]这里,电场∇φ,φ 是静电势表示压力张量.在托卡马克中,电子可以分成俘获电子和通行电子两部分.对于前者,反弹平均后的平行速度为零,俘获电子的密度net和温度Tet可以由连续性方程和能量守恒方程给出[10]:其中=2.5netTe(×∇Te)/eB2表示电子抗磁漂移热通量,St,f= -vth(δnet -Γneeφ/Te),Qt,f=vthTene(βeφ/Te -0.5δnet/ne),St,f,Qt,f 分别表示库伦碰撞引起的俘获粒子和通行粒子之间的粒子交换和能量交换.Γ = 1 + (dlnTe/dlnne)/(ω/ωDe -1 +iνth/ωDe),β =1.5 -2.5Γ,νth=Rνei/r,δnet,ωDe=2kyTe/(eBR)分别表示电子—离子有效碰撞频率,俘获电子扰动密度,电子抗磁漂移频率,其中r表示小半径,vei为电子离子碰撞频率.对于通行电子,在平行磁场方向的动量方程中,当忽略惯性力和热力学力时,电子的压强梯度力和电场力平衡:1.2 线性化方程假设稳态时粒子的温度、密度都只是小半径r的函数,平行方向的速度和静电势均为零.考虑上述参量在时间上随频率w变化,空间上,沿磁场方向l和垂直磁场方向y上(波数分别为k//和ky,其中ky包括垂直磁面的波数和在磁面上垂直于磁力线的波数)的一个微小扰动δf,扰动波形式为δf~ exp(-iωt+ik//l+ikyy).方程(1)、方程(2)、方程(3)线性化,有:其中=w/wDe分别为归一化扰动电势、密度、温度和频率;gx=-R*dlnx/dr.为归一化等离子体参量(下标x表示电子或离子的温度和密度);τH=TH/Te,同样,方程(5)、方程(6)线性化处理后可得到:其中,=vth/ωDe为归一化有效碰撞频率,为俘获电子份额,f1=ft(gTe/2-gne/3)+i·(5ft/3-2/3).1.3 ITG和TE模色散方程对离子,由方程(8)、方程(9)、方程(10)有对俘获电子,由方程(11)、方程(12),有:其中,=δnH/nH=eφ/Te分别为氢离子归一化扰动密度,归一化扰动电势.PH,Pet表达式见附录所示.对通行电子,由于其平行磁场方向高的热导率,温度扰动可忽略不计,由玻尔兹曼关系和(7),有:其中,=δnef/nef为俘获电子归一化扰动密度.此外,根据准中性条件:把式(13)、式(14)、式(15)代入式(16)中,则 ITG模和TE模色散关系可表示为:这里,方程(17)是关于ky和w的多项式.其中,w=wr+iγ,实部wr表示漂移波的频率,虚部表示阻尼,当γ<0时,漂移波随时间衰减;当γ>0时,漂移波随时间增长.计算中,取增长率最大的模,当γ>0时,若wr<0时,γ表示ITG的增长率;若wr>0,γ表示 TEM 的增长率.本文中,取0<kyρsH <2.0,并取该波长范围内最大波长 kmax-yρsH所对应的γmax作为该不稳定性模的增长率.2 计算结果与分析等离子体芯部区域ITG和TE不稳定性模的一个重要特征是存在阈值.其增长率与等离子体温度梯度、密度梯度、归一化径向坐标等参量密切相关,本文主要研究了上述等离子体参量对HL-2A装置芯部ITG和TEM增长率的影响.2.1 gT对ITG和TEM的影响图1和图2分别给出了gTe+gTH=5时,表1所给参数下HL-2A装置中ITG和TEM增长率γmax以及γmax 所对应的波长 kmax-yρsH 随参量θ(θ=(gTe-gTH)/(gTe+gTH))的变化情况.图1 增长率γmax与θ关系Fig.1 θ-dependences of γmax图2 kmax-yρsH与θ关系Fig.2 θ-dependences of kmax-yρsH从图1和图2可以看出,对于ITG模,增长率γmax随θ的增大而逐渐减少,当θ≥0时,ITG完全被抑制住了,这是由于此时离子温γmax度梯度gTH低于ITG 的临界温度梯度阈值,容易算出ηH=1.25(ηH=gTH/gnH) ,该结论与 Chen.L [11-12]等人早期的理论预测吻合.此外,ITG的γmax所对应波长 kmax-yρsH随θ的增大而缓慢降低.对于 TEM,在整个区间都是不稳定的,当-1<θ<-0.4时,增长率γmax随θ的增大而减少,当-0.4<θ<1时,增长率γmax 随θ的增大而增大,这是由于TEM是由温度梯度和密度梯度共同驱动,此时密度梯度gn=2,TEM已经被密度梯度激发起来了.此外,TEM的γmax所对应波长kmax-yρsH也是随θ的增加先变大后减小.表1 HL-2A装置计算参数Table 1 The calculation parameters of HL-2A注:表1中,温度Te(TH),密度ne(nH),小半径a(大半径R,径向位置r),磁场B单位分别为KeV、1019m-3、m和T.Te TH ne(nH) gne(gnH) qsf S a R r B 2 1.5 2.25 2 1 2 0.4 1.65 0.16 2.82.2 gn对ITG和TEM的影响图3和图4分别给出了表1参数下不同归一化温度梯度 gn(0.5,1.5,2.5,3.5)时 ITG 和 TEM增长率随gT的变化情况.图3 不同gn时ITG增长率与gTH关系Fig.3 gTH-dependences of the ITG growth rate in different gn从图3可以看出,ITG的温度梯度阈值随着密度梯度gn的增大而增大,这说明密度梯度对ITG有致稳作用,这是由于ITG由ηH(gTH/gnH)驱动,其值保持在一定范围内.当gn增大时,显然ITG的温度梯度阈值也会相应增大.此外,当gTH>9时,γmax几乎不随gn变化,这说明此时ITG已基本由温度梯度所驱动.从图4可以看出,TEM增长率γmax随gn的增大而增大,当gn=0.5时,TEM还存在温度梯度阈值,但之后随着gn的增大,TEM在整个温度梯度区间都被激发起来了,这是由于TEM是由温度梯度和密度梯度所共同驱动的,此时,密度梯度已足够激发TEM.图4 不同gn时TEM增长率与gTe关系Fig.4 gTe-dependences of the TE growth rate in different gn2.3 ρ对ITG和TEM的影响取表1计算中参数计算,图5和6分别给出了不同归一化径向坐标ρ(r/a)(0.2,0.4,0.6)时ITG和TEM增长率随gT的变化情况.图5 不同ρ下ITG增长率与gTH关系Fig.5 gTH-dependences of the ITG growth rate in different ρ图6 不同ρ下TEM增长率与gTe关系Fig.6 gTe-dependences of the TEM growth rate in different ρ从图5可以看出,ρ对于ITG阈值几乎没有影响,但一旦超过温度梯度阈值时,ITG模增长率γmax会随着ρ的增大而稍微增大;从图6可以看出,对于TEM,其增长率会随着ρ的增大而增大,且增加幅度随着温度梯度的增大而增大.这是由于俘获电子所占的份额ft随着ρ的增大而增大,从而促进了TEM的增长.3 结论从双流体模型出发,给出了芯部ITG和TEM的色散关系.在Weiland模型的基本上,考虑了离子平行方向速度、俘获电子和碰撞对ITG和TEM的增长率的影响,并分析了归一化温度梯度、密度梯度以及归一化径向坐标对ITG和TEM增长率的影响.对HL-2A托卡马克等离子体的计算结果表明:1)对于ITG,存在温度梯度阈值,一旦超过温度梯度阈值后,增长率γmax会随温度梯度增大而增大.离子密度梯度对ITG有致稳作用,可以使ITG温度梯度阈值增大,且使增长率γmax减小.2)对于TEM,当电子密度梯度比较小时,TEM也存在温度梯度阈值,但当电子密度梯度超过某一值后,TEM在整个温度梯度区间都是不稳定的,其增长率γmax会先随着温度梯度的增大而减小,而后又随着温度梯度的增大而增大.3)ITG和TEM增长率都会随归一化径向坐标ρ的增大而增大,TEM的增加幅度更大.参考文献:[1]Liewer P C.Measurements of microturbulence in tokamaks and comparisons with theories of turbulence and anomalous transport [J].Nucl.Fusion.,1985,25(5):543-622.[2] Horton W.Drift waves and transport[J].Rev.Mod.Phys.,1999,71(3):735-778.[3]Garbet X,Idomura Y,Villard L,et.al.Gyrokinetic simulations of turbulent transport[J].Nucl.Fusion.,2010,50(4):043002-1-30.[4]Waltz R E,Candy J,Fahey M.Coupled ion temperature gradient and trapped electron mode to electron temperature gradient mode gyrokinetic simulations[J].Phys.Plasmas.,2007,14(5):056116-1-8.[5]Dimits A M,Bateman G,Beer M A,et parisons and physics basis of tokamak transport models and turbulence simulations[J].Phys Plasmas.,2000,7(3):969-983.[6]Candy J,Waltz R E.An Eulerian gyrokinetic-Maxwell solver[J].J Comput Phys.,2003,186(2):545-581.[7]Jenko F,Dannert T,Angioni C.Heat and particle transport in a tokamak:advances in nonlinear gyrokinetics[J].PlasmaPhys.Control.Fusion.,2005,47:B195-B206.[8]Bateman G,Kritz A H,Kinsey J E.Predicting temperature and density profiles in tokamaks[J].Physics of Plasma.,1998,5(5):1793-1799. [9]Weiland J.Collective modes in inhomogeneous plasma:kinetic and advanced fluid theory[M]Bristol:Institute of Physics Publishing,2000. [10]Nilsson J,Weiland J.Fluid model for general collisionality and magnetic curvature[J].Nucl.Fusion.,1994,34(6):803-808.[11]Terry P,Anderson W,Horton W.Kinetic effects on the toroidal ion pressure gradient drift mode[J].Nucl.Fusion.,1982,22(4):487-498. [12]Guzdar P N,Chen L,Tang W M,et al.Ion-temperaturegradient instability in toroidal plasmas[J].Phys Fluids.,1983,26(3):673-677.。
托卡马克物理基础-3[1]
III. 托卡马克平衡位形的求解这一部分主要介绍Grad-Shafranov 方程的求解方法及相关的物理内涵。
重点是内部解的求解方法。
在磁面截面形状接近为圆时(磁面簇为套圆环系统),如果环境比1/>>=a R A ,我们可以把倒环境比A /1=ε作为展开小量,获得近似解。
在另外一些特殊电流剖面下,还可以获得精确的位形解。
III.1 圆截面托卡马克圆截面托卡马克不仅是一种数学近似,而且也是一些真实托卡马克装置的位形。
早期的托卡马克外部有一个很厚的铜壳,在等离子体环发生向外位移时会自动产生象电流,这个象电流的垂直磁场用以维持环的平衡以及稳定一些快速增长的磁流体不稳定模(扭曲模)。
Shafranov 本人对这类平衡位形做了大量的研究工作,奠定了托卡马克位形的平衡理论基础。
所以,平衡位形中各个磁面之间的相对位移被称为Shafranov 位移。
在圆柱坐标系中写出Grad-Shafranov 方程:),(21102202222ψμψψμψψψζR Rj d dF d dP R Z R R R -=--=∂∂+∂∂-∂∂ (3-1)先介绍大环境比近似展开法。
III.1.1大环境比展开法设1/<<=R r ε,这时, 以磁轴为原点,引入局部极坐标),,(ζθρ:θθρθθρsin ),(,cos ),(0r Z r R R =+= (3-1-1)图 局部极坐标取ρ为平均磁面半径⎰=θπρrd 21(3-1-2)将),(θρr 展开为Fourier 级数:]sin )(cos )([),(1θρθρρθρδρn b n a r n n n ++=+=∑=(3-1-3)对于上下对称的位形,上式中右方第二项为零。
系数1a 为Shafranov 位移函数,2a 与截面的椭圆变形有关,3a 与三角变形有关,等等。
在),(θρ坐标系中,方程(3-1)化为:),(),(),(22θρρψθρρψθρY d d C d d A =+(3-1-4)2222ρθr r r r A +=(3-1-5)]2)(sin cos 1[1222r r rr r r rr r r r R r r rr C θθρρθθρθρρθρθθ-++-+-= (3-1-6)式中下标表示对该变量求偏导数。
带电粒子在电磁场中几种漂移运动以及托克马克约束下的运动
带电粒子在电磁场中几种漂移运动以及托克马克约束下的运动作者:吴小锋来源:《中国科技纵横》2019年第16期摘 ;要:本文将以回旋单粒子模拟方法为基础,研究带电粒子在不同种类的电磁场中的运动状况,并带入简化的托克马克磁场中模拟观察其运动轨迹,以验证利用托克马克装置约束等离子体的可行性。
关键词:等离子体;单粒子模拟;磁约束中图分类号:O441;;文献标识码:A ;;文章编号:1671-2064(2019)16-0000-000引言托卡马克,是一种利用磁约束来实现受控核聚变的环形容器。
托卡马克的中央是一个环形的真空室,外面缠绕着线圈。
在通电的时候托卡马克的内部会产生巨大的螺旋型磁场,将其中的等离子体加热到很高的温度,以达到核聚变的目的。
描述等离子体有很多种方法,比如理想磁流体模型MHD。
而回转运动粒子模拟作为最早出现的一种模拟等离子体和描述等离子体湍流和输运的方法,其理论和实用性得到了广泛的实验证明,现如今它已成为一种强大而可靠的工具。
[1]1回旋单粒子模拟方法单粒子模拟通过将等离子体内带电粒子如电子等看做带一定量电荷和质量的质点,忽略粒子的其它属性和与其它粒子的相互作用,通过运动方程计算得到带电粒子在一定电场和磁场作用下的运动轨迹,从而由单个带电粒子的微观状态来研究等离子体的宏观状态。
1.1电子在恒定均匀的磁场下的漂移运动如果电量为q的粒子在磁场中除了受到恒定均匀磁场B作用外,还受到其他外力F的作用,则粒子除了以磁力线为轴的螺旋运动外,还要在垂直于磁场B和外力F外的方向运动,这个运动的速度矢量v D=(F×B)/(qB2)。
[1]1.2电子在均匀电场和均匀磁场的漂移运动假设电场强度为1,方向为x轴正方向;磁场强度为1,方向为z轴正方向电场项:Ex[x_, y_, z_] = 1; Ey[x_, y_, z_] = 0; Ez[x_, y_, z_] = 0;磁场项:Bx[x_, y_, z_] = 0; By[x_, y_, z_] = 0;Bz[x_, y_, z_] = 1;此时电场作为电子所受外力F的来源。
TOKAMAK装置优缺点简析
简析托卡马克装置的优缺点摘要:托卡马克装置是目前最有希望实现可控核聚变反应的途径。
虽然托卡马克存在着很多尚未解决的问题,例如磁场的起伏、热压强的影响,内部电流的不稳定等,但只要这些影响不大且发展得足够缓慢,我们仍然可以实现控核聚变。
本文简要分析了托卡马克装置的优缺点。
一、托卡马克装置简介1.1托卡马克装置的产生背景能源是社会发展的基石。
以煤炭、石油、天然气等化石能源替代柴薪的第一次能源革命带来了社会经济的飞速发展。
然而这些宝贵的资源就这样被燃烧掉,同时造成了严重的污染。
据估计,一百年后地球上的化石能源将会面临枯竭。
面对着即将来临的能源危机,人类有了一个共同的梦想——寻求一种无限而清洁的能源来实现人类的持续发展[1]。
托卡马克(Tokamak)核聚变是一种利用磁约束来实现受控的核聚变。
它的名字Tokamak来源于环形(toroidal)、真空托卡马克核聚变空室(kamera)、磁(magnit)、线圈(kotushka)。
托卡马克核聚变的中央是一个环形的真空室,外面缠绕着线圈。
在通电的时候托卡马克的内部会产生巨大的螺旋型磁场,将其中的等离子体加热到很高的温度,以达到核聚变的目的。
托卡马克(Tokamak)是前苏联科学家于20世纪50年代发明的环形磁约束受控核聚变实验装置。
经过近半个世纪的努力,在托卡马克上产生聚变能的科学可行性已被证实,但相关结果都是以短脉冲形式产生的,与实际反应堆的连续运行有较大距离。
超导技术成功地应用于产生托卡马克强磁场的线圈上,是受控热核聚变能研究的一个重大突破。
1.2托卡马克核聚变实验装置“超导托卡马克核聚变”实验包括一个具有非圆小截面的大型超导托卡马克实验装置和低温、真空、水冷、电源及控制、数据采集和处理、波加热、波驱动电流、诊断等子系统。
其中超导托卡马克装置是本项目的核心。
而超导托卡马克装置又包括超导纵场与极向场磁体系统、真空室、冷屏、外真空杜瓦及面对等离子体部件等部件[2]。
托卡马克工程
2,包层 包层(blanket)是聚变堆中置于真空室内面对等离子
体的屏蔽层,功能是吸收中子、传输能量、氚增殖。 ITER 的包层由冷却剂(氦或水)、氚增殖剂、中子
增殖剂和结构材料组成。氚增殖剂为含锂材料,由下列 反应产生氚:
6Li + n → 4He + T + 4.78MeV 7Li + n → 4He + T + n -2.47MeV 第二个反应虽是吸热的,但又产生一个中子,所以两 种增殖剂混合使用,可以使氚的增殖系数大于1。 为进一步增殖中子,使用铍或铅作中子增殖剂,发生 如下的反应
伏秒数ΔΦ决定于变压器线圈(螺线管和外线圈)的 几何尺寸或者说其电感,以及放电电流。为了增加这 个伏秒数,一般使变压器先向反向磁化到最大电流, 再使电流反向。此时使气体击穿,产生等离子体和感 应等离体电流。这样,在最大电流为一定值条件下, 使伏秒数加倍。
在一次放电中,首先启动环向场,在它达到或接近 平顶时,启动变压器反向磁化,但不使气体击穿。反 向磁化达最大电流后,在正向磁化,使气体击穿,产 生等离子体。
在B∝1/ R磁场形态中,存在一种纯张力线圈。电 动力在其中引起的应力是纯张力,无垂直方向的切 应力。从受力角度这样的线圈比较合理。它的形状 可用微分方程表示:
ky 1 y3 2
y
其中k为一常数。这一方程无解析解,但可得到一积 分式
y ln xdx C k 2 ln2 x
常数C说明它可上下任 意移动。如果要上下对称 则 C=0 。 常 数 k 决 定 曲 线 的形状。这一线圈接近于 D形。所以在一些大的装 置中,环向场线圈往往做 成这一形状或称D形。当 前所有大型装置的环向磁 体均采用D形的纯张力线 圈。
的气体,表面状况对于等离子体的纯度至关重要。 ➢ 除去真空性能外,还要求真空室材料的电阻高、无
全超导托卡马克的运行风险和需要采取的对策
全超导托卡马克的运行风险和需要采取的对策全超导托卡马克的定义:不仅纵场磁体而且所有主要的极向场磁体均为超导磁体的托卡马克。
一风险1常规托卡马克所有磁体系统在设计运行的参数范围内均是自安全的;2所有超导托卡马克的超导磁体即使在其设计的运行参数范围内均会在某些运行条件下发生失超,如不加以保护,不能及时将磁体内储存的磁能泄放,超导磁体或与其相关联的某些重要部件将被烧毁;3对于只有纵场是超导磁体的超导托卡马克:运行风险主要来源于超导纵场磁体系统失超。
引起失超的原因可以是:励磁电流过快、过大;制冷系统发生故障;磁体上局部地点温度升高;等离子体电流快速破裂;外杜瓦真空被破坏等等。
失超后的具体风险为:3.1由于纵场是托卡马克上储能最大的磁体系统,因此,发生失超时如不能及时进行失超保护,将会造成磁体系统或相关部件烧毁的严重事故;3.2如果及时进行了失超保护,但泄放回路的时间常数太小,则会使泄放时感生的电压过高而造成部件绝缘击穿,烧毁磁体或相关部件;3.3纵场系统失超时在与之有最强耦合的真空室和内冷屏的极向方向上将感应出电流,这一电流将与纵场作用产生小截面上的扩张力和大环的收缩力,它是一个将造成真空室和内冷屏剧烈振动的冲击力;3.4由于在托卡马克装置上极向场系统与纵场系统的磁耦合很弱,因此在只有纵场是超导磁体的超导托卡马克上极向场的快速变化(包括等离子体电流的破裂)只会对外杜瓦内的纵场磁体系统造成小的影响;4对于全超导托卡马克:在外杜瓦内不仅有超导纵场磁体系统而且全部极向场超导磁体系统也在同一外杜瓦内,由此,除了具有上述纵场系统失超引发的同样风险外还具有更大的运行风险,它是来源于所有极向磁场系统(包括等离子体电流)的特殊运行要求和自身相互之间的强耦合,具体是:4.1等离子体电流的建立必需极向场线圈系统提供极快速(击穿)和快速(电流爬升)的磁通变化:4.1.1在这一阶段,极向场线圈,特别是中心螺管线圈极容易发生失超;4.1.2同时所有与之强耦合的,构成环向回路的金属部件上将会引发感应电压或涡流,前者可以引发电弧,后者将引起部件发热和电动力;4.1.3电弧将会破坏外杜瓦真空,并有可能像预电离一样在外杜瓦内引发更大面积的放电,造成不可收拾的严重后果;4.1.4涡流将引发冷质部件的发热,从而也有可能引发线圈失超;4.1.5涡流引发的电动力将有可能破坏或逐渐破坏绝缘从而引发电弧4.1.6所有与之强耦合的极向场线圈上(包括中心螺管线圈自身)将引发感应高电压,如果任何地方绝缘薄弱或损坏,则一定会引发线圈及其部件烧毁的严重事故;4.1.7在极向场线圈上感应的高电压将一直传递到电流引线箱,如果在超导电流传输线和电流引线箱内绝缘薄弱、损坏和真空度下降均会引发起弧、放电和因此烧毁部件。
托卡马克装置原理2
单位长 µ 0 l 度内感 4π i
li 1 l a2 a2 ( β p + − ) < ∆ (o) < (β p + i ) 2 R0 2 2 2 R0 2
取决于电流压强轮廓
真空磁场
在真空区域解齐次平衡方程,即M方程,引进
Λ = βp +
µ0 I p µ0 I p 8 R0 r 1 a2 − 2) − ψ = R0 (ln r[ln + (Λ + )(1 − 2 )] cos θ 2π 4π 2 r a r µ0 I p µ0 I p 8 R0 8 R0 8R 1 a2 1 = − 2) + − 1) + (Λ + ) 2 ] − (ln 0 + Λ − )} cos θ R0 (ln r{[ln( 2π r 4π r 2 r a 2
jR = −
1 ∂F R ∂z
两组量的 对应关系
v A
v B
µ0
v B v j
Ψ
ψ
F
磁场量
I
电流量
Grad-Shafranov方程 方程
从磁面量的定义和M方程
1 ∂ψ BR = − R ∂z
∂BR 1 ∂ 2ψ =− R ∂z 2 ∂z
∂B z 1 ∂ 2ψ 1 ∂ψ = − 2 ∂R R ∂R 2 R ∂ψ
∂p = jR Bϕ − jϕ BR ∂z
j z / j R = B z / BR
在子午面内电流与磁场的投影同向, 即电流磁场在一个环面上
压强 p的梯度垂直于这一面,所以在面上 p=const,称为磁面
托卡马克等离子体的磁面是嵌套的环面. 最里面的磁面退化为磁轴
等离子体温度和密度在磁面上近似为常量
托卡马克装置资料
ASIPP
HT-7
等离子体在环形螺旋磁场中的约束
应该注意,如果有了磁场误差或其它的非轴对称磁场, 则磁力线在多次绕环以后,它们常常和器壁相交,粒 子也就约束不住。在托卡马克装置中一个无法避免的 缺陷是装置只能使用有限个纵场线圈,引起纵向磁场 的起伏,这就会沿磁力线产生非常浅的局部磁镜,它 们能够捕获一小部分等离子体粒子:这种粒子既不围 绕环的小截面画圆圈,也不相对于赤道平面对称地被 捕获,它们将漂移出系统:除了这些磁捕获效应以外, 一个电场本身或者和磁场一起也能引起粒子捕获,使 粒子回旋中心轨道发生很大的变化:为了减小磁场误 差,一般需要细致的进行纵场线圈的设计。
II、世界上现有g束线站的简 介 高能g光束线
光子束性能 装置 产生方法 ē Brem. coh.Brem. ē Brem. coh.Brem. ē Brem. coh.Brem. Compton Compton Compton Compton Compton 能量(MeV) 140-800 -400 500-2400 -1.2 300-2300 -1800 -1600 180-320 280-470 550-1470 1500-2400 200-870 强度 (10 s MeV-1)
50 10 5 50 2 1 33 1 1 10
Spring-8 Result
Q+
Q+
g+ C(n/p) K- K+ (n/ p) Cuts: no fK+Kno recoil p (gn only) Missing mass for n
M = 1.540.01 GeV G < 25 MeV Gaussian significance 4.6s
ASIPP
托卡马克等离子体中逃逸电子动理学的实验研究
报告内容
• • • • • 1. 2. 3. 4. 5. 基本理论; 研究意义; 选题内容; 现有诊断系统; 已做工作。
选题内容
• 1. 逃逸电子的产生机制、能量分布以及诊 断逃逸电子的可行探测手段。 • 2. 逃逸电子的两种产生机制在不同放电条 件和放电阶段的不同作用效果。 • 3.低杂波电流驱动下剩余环电压对逃逸电 子的作用效果。 • 4.利用低杂波、磁场波纹共振或通过注入 惰性气体等方法来寻求抑制逃逸电子可行 手段。( )
(correspond to 50 times the electron content)
Injection of
20.9 ×1021atoms of argon
(correspond to 200 times the plasma electron content)
1. Argon injection can cause the generation of runaways carrying up. 2. Disruptions caused by argon injection finally become runaway free for very large numbers of injected atoms (about 140 times the plasma electron content) 3. Disruptions triggered by injections of deuterium, helium or of mixture of argon with deuterium are runaway free.
Lower hybrid current drive (LHCD)
1.低杂波在0.44s注入; 2.环电压由1.52V降到0.14V; FEB 3.FEB保持在一个较高的水平;
托卡马克的基本慨念
其他情况: 如果有了磁场误差或其它的非轴对称磁场, 则磁力线在多次绕环以后, 它们常常和器壁相交,粒子也就随磁力线碰壁。 托卡马克装置只能使用有限个纵场线圈,引起纵向磁场的起伏,这就会 沿磁力线产生非常浅的局部磁镜,它们能够捕获一小部分带电粒子,它们容 易漂移出系统(见下图) 。
(利用永久磁铁降低纵场线圈纹波,提高粒子约束)
4. 托卡马克等离子体平衡、平衡(Grad Shafranvo)方程
导言: (a)托卡马克中等离子体的质量非常小,一般仅 10-4 克/m3,不考虑重力因素: (b)载流的环形等离子体柱有向外扩张的趋势,且受力大(10 吨/m3) ,如不设
法加以平衡,等离子体就会向外运动碰撞容器壁。
(c)沿磁力线方向的压力平衡:容易达到,速度很快(压力的传递以类似声波的
Triangularity 三角变形度 等离子体截面的三角形程度的一种量 度。
Vertical Displacement Event 垂直位移事件 整个等离子体在离其平衡位置向上乮或 向下乯运动期间的一种不稳定性。高拉 长度等离子体更有这种运动的倾向丆因 为产生这些需要更强的成形磁场。要是 不用反馈系统控制丆当等离子体与容器 碰撞时丆它将迅速地损失掉丆并导致等 离子体电流流过壁和其它部件。
P
R
子体在R方向受到一个向外的力.而等离子 体电流在极向磁场中,收到向内的力,正好 与之平衡所以等离子体在R方向平衡.
JB
R 描述等离子体平衡的基本方程组 前面的力学平衡方程: 再加上麦克斯韦方程:
J B P
B 0 J
B 0
这三个方程组成描述托卡马克和其他轴对称环形平衡位形的基本方程组。 它们完整地
螺旋磁场的不均匀性 大环外侧磁场最小;大环内侧最大。 造成的结果:因为有磁场梯度存在,所以带电粒子通旋中心的运动是 由沿磁力线的导引运动和磁漂移两部分合成。 两类粒子 “捕获粒子”或“约束粒子”: 托卡马克等离子体外侧磁场 比内侧磁场低。 平行于磁场的低速 度粒子没有足够的能量进入强场 (内侧)区,并且在外侧被捕获,
托卡马克装置的核心组件分析
托卡马克装置的核心组件分析揭开科学界的神秘面纱,人类对于能源领域的追求一直不曾停止。
而托卡马克装置作为核聚变研究的关键设备,正在被广泛关注。
在这篇文章中,我们将对托卡马克装置的核心组件进行分析,探寻其背后的奥秘。
首先,让我们来了解一下托卡马克装置的基本原理。
托卡马克装置是一种通过控制高温、高密度等条件,实现核聚变反应的装置。
其核心组件主要包括等离子体环状炉体、磁场线圈、导向磁场等。
等离子体环状炉体是托卡马克装置的核心,它是容纳等离子体的空心环形结构。
等离子体是物质的第四态,是由高能粒子通过碰撞、电离等过程产生的带有正电荷的气体。
等离子体环状炉体的材料通常选用耐高温、耐辐射的特殊合金,以承受高能粒子冲击和极高温度。
在等离子体环状炉体的外部,布置有磁场线圈。
磁场线圈通过产生强大的磁场,将等离子体束缚在环状炉体的中心区域,形成稳定的磁约束。
这一磁场约束是核聚变反应能够进行的关键,它能够有效减少等离子体与环状炉体之间的相互作用,延长等离子体的滞留时间,提高反应效率。
同时,托卡马克装置中还需要导向磁场来引导等离子体的运动轨迹。
导向磁场通过控制等离子体中的带电粒子的运动方向,维持等离子体在环状炉体中的循环运动。
这种稳定的循环运动是保持等离子体稳定和均匀分布的关键。
除了上述核心组件外,托卡马克装置还需要带电粒子注入系统、真空系统、热控制系统等配套设备的支持。
带电粒子注入系统用于注入带有高能量的粒子,提高等离子体的能量密度。
真空系统则用于维护装置内部的真空环境,避免气体的干扰。
热控制系统则用于控制等离子体的温度,避免装置因高温而受损。
总结起来,托卡马克装置是一种将高能粒子束缚在稳定磁场中,通过控制温度、密度等参数实现核聚变反应的设备。
其中,等离子体环状炉体、磁场线圈和导向磁场是其核心组件。
托卡马克装置的研究对于未来的能源开发具有重要意义,它有望成为可持续、清洁的能源选择。
然而,托卡马克装置的核心组件研发面临着许多挑战。
托卡马克原理简介:核聚变的未来
托卡马克原理简介:核聚变的未来核聚变是一种将轻元素聚集在一起,通过释放巨大能量的过程。
科学家们一直追求着核聚变技术的发展,因为它被认为是未来清洁和可持续能源的解决方案。
而托卡马克原理则是实现核聚变的一种重要方法。
1. 什么是托卡马克原理?托卡马克原理是一种通过在磁场中限制等离子体的运动来实现核聚变的理论和实验技术。
其核心设备被称为托卡马克装置。
所谓等离子体,是一种高温高能状态下的物质,具有高度离子化的性质,常见于太阳和恒星中。
2. 托卡马克原理的工作原理托卡马克原理基于等离子体受磁场约束的特性。
当等离子体的温度和气压达到足够高的程度时,原子核开始碰撞并融合在一起,释放出巨大的能量。
为了实现这一过程,科学家们将等离子体置于一个扁平环形状的装置中,通过外部的磁场使等离子体保持稳定。
这种装置形状被称为“托卡马克”,得名于俄罗斯单词“тор”和“焊接”。
3. 托卡马克装置的主要部件托卡马克装置由几个关键组件组成。
首先是真空室,用于创造一个低压环境,使等离子体不与外部气体发生碰撞。
其次是加热系统,用来提供外部能量来加热等离子体以达到聚变所需的温度。
同时,还需要一个磁场系统,能够提供足够强度和稳定性的磁场来限制等离子体的运动。
最后是等离子体诊断装置,用于监测等离子体的性质和行为。
4. 托卡马克原理面临的挑战虽然托卡马克原理潜力巨大,但也存在一些挑战和困难。
首先,高温高能的等离子体对材料的破坏性很大,因此装置材料的选择和应用也是一个挑战。
其次,维持长时间的高温等离子体状态需要大量的能量输入和精确的控制,这也是一个技术难题。
此外,聚变反应不稳定性、等离子体的热损失等问题也需要解决。
5. 托卡马克原理的前景尽管面临着种种挑战,托卡马克原理仍然被认为是实现核聚变的一项重要技术。
当前,许多国家和国际组织正在进行相关的研究和实验,试图克服这些难题。
如果能够实现可控的核聚变反应,将为人类提供可持续、高效、清洁的能源解决方案。
五、托卡马克中的重要问题(磁约束、平衡、加热、第一壁之外)
五、托卡马克中的重要问题(磁约束、平衡、加热、第一壁之外)五、托卡马克中的其他重要问题(磁约束、平衡、加热、第一壁之外)1. 托卡马克物理发展的重要点19世纪30年代英国的M.法拉第以及其后的J.J.汤姆孙、J.S.E.汤森德等人相继研究气体放电现象,这实际上是等离子体实验研究的起步时期。
1879年英国的W.克鲁克斯采用“物质第四态”这个名词来描述气体放电管中的电离气体。
美国的I.朗缪尔在1928年首先引入等离子体这个名词,等离子体物理学才正式问世。
1929年美国的L.汤克斯和朗缪尔指出了等离子体中电子密度的疏密波(即朗缪尔波)。
对空间等离子体的探索,也在20世纪初开始。
1902年英国的亥维赛(发现地球上电离离层的中层,E层,被称为亥维赛层)等为了解释无线电波可以远距离传播的现象,推测地球上空存在着能反射电磁波的电离层。
这个假说为英国的E.V.阿普顿用实验证实。
英国的D.R.哈特里(1931)和阿普顿(1932)提出了电离层的折射率公式,并得到磁化等离子体的色散方程。
――――以下与托卡马克密切相关(在高温等离子体书中有对应内容)―――― 从20世纪30年代起,磁流体力学及等离子体动力论逐步形成。
等离子体的速度分布函数服从福克,普朗克方程。
苏联的Л.Д.朗道在1936年给出方程中由于等离子体中的粒子碰撞而造成的碰撞项的碰撞积分形式。
1938年苏联的A.A.符拉索夫提出了符拉索夫方程,即弃去碰撞项的无碰撞方程。
朗道碰撞积分和符拉索夫方程的提出,标志着动力论的发端。
1942年瑞典的H.阿尔文指出,当理想导电流体处在磁场中,会产生沿磁力线传播的横波(即阿尔文波)。
印度的S.钱德拉塞卡在1942年提出用试探粒子模型来研究弛豫过程。
1946年朗道证明当朗缪尔波传播时,共振电子会吸收波的能量造成波衰减,这称为朗道阻尼。
朗道的这个理论,开创了等离子体中波和粒子相互作用和微观不稳定性这些新的研究领域。
从1935年延续至1952年,苏联的H.H.博戈留博夫、英国的M.玻恩等从刘维定理出发,得到了不封闭的方程组系列,名为BBGKY链。
托卡马克物理基础-4
IV 托卡马克等离子体的约束模式从上世纪六十年代以来,全世界已建造了大量的托卡马克实验装置,其中特别重要的如JET(欧洲联合环),JT-60(日本),TFTR (PPPL ,US), DIII-D(GA,US ) 等四个超大装置, 以及以ASDEX, ADDEX-U(德),TORE-S (超导,法), ALCATOR-CMOD (强磁场,US), FT (强磁场,意大利), TEXTOR(德),TCV (瑞士),T-10 (俄),TEXT( US), HL-1(中国乐山),HT-7(超导,中国合肥)等大装置,以及近期投入运行的HL-2A (中国成都),EAST(超导,中国合肥)等新装置,都对托卡马克的等离子体的约束特性做了大量的研究。
其中四个特大装置的等离子体参数都达到聚变堆所要求的水准,即形成了堆芯等离子体(等离子体温度10keV 量级)。
所有的装置都以欧姆加热所形成的初始等离子体为基准,进行了不同程度的高功率加热,包括中性束注入(Neutral Beam Injection, NBI ),离子回旋频段波加热 ( Ion Cyclotron Frequency Range Heating, ICFRH),低混杂波电流驱动和加热 ( Lower Hybrid Wave Current Drive & Heating ,LHCD&LHH) 电子回旋共振加热 ( Electron Cyclotron Resonance Heating, ECRH) 等。
所积累的数据库包含极其丰富的各类信息。
对托卡马克等离子体约束规律的研究的最有效的推动是ITER 设计。
历时两年多(1988-1990)的ITER 的物理设计及跟着进行的近十年(1991-1999)的工程设计,以及几乎又是快十年的后续更新设计(2000-2007),使得对这些规律的认识极大地得到深化。
形成了托卡马克型反应堆物理和工程设计的基础。
微观不稳定性
7.3 微观不稳定性鉴于环形等离子体的输运水平远超过新经典输运理论的预言,很早就将这一反常输运现象的来源归于湍流输运。
L-H 模转换发现以后,特别是观测到随输运的降低,等离子体的涨落水平也相应降低,进一步将输运和微观不稳定性联系起来。
微观不稳定性指系统偏离Maxwell 分布而产生的不稳定现象,例如发生在速度空间的束流不稳定性(逆朗道阻尼)和损失锥不稳定性。
而和环形装置输运问题联系的,主要是由等离子体自由能推动的微观涨落现象,尺度为离子回旋半径量。
广义微观涨落的还包括对流元(convective cell )、随机磁场等机构。
但实验研究表明,所观测的不稳定性,一般具有漂移波的特征。
1,湍流测量结果边界区 在边界区,主要是静电探针测量,得到的密度涨落在30%以上,可达100%。
电位涨落和密度涨落关系有些偏离Maxwell 分布的关系,有nnT e e ~~>φ。
用静电探针得到的温度涨落为%20~>eeT T ,在不同装置上有所不同。
图7-35是TEXT 上的温度、密度和电位的相对涨落振幅,也表示了磁涨落的相对振幅(径向涨落和环向磁场比较)。
密度涨落频谱在高频有一临界频率,处于10-100kHz 范围,与装置有关。
在这个频率以上,功率谱以负指数f -1~f -4衰减。
在频谱图的横轴“归一化”后,不同装置谱的形状非常相似(图7-36)。
扰动的极向相速<106cm/s ,在电子逆磁漂移方向,但在分支磁面以外,改变为离子逆磁方向。
而径向速度小与极向速度。
扰动的结构大致是沿磁场延伸的丝状结构,宽度为cm 量级。
极向的k 谱类似频谱。
图7-35 TEXT 上边缘区相对涨落水平 图7-36 几个装置上边缘区密度涨落功率谱核心区 可用集体散射、重离子探针、束光谱等诊断,探测的涨落水平只有1%左右。
其频率宽度ωΔ为100kHz 左右,但可能来源于极向剪切流的Doppler 效应。
波数谱主要由长波组成,径向相关长度约2-3cm ,极向波数谱的峰值为3.0≤⊥L r k ,L r 为离子回旋半径。
托卡马克物理基础-2[1]
II 磁约束聚变装置和反应堆的约束原理聚变反应必须在高温下才能进行。
按粒子动能与温度的转换关系1eV=11600 K. 通常我们简单地将 1keV 与1千万度互相对应。
氘氚聚变要在10keV 以上才能充分进行,对应于1亿度的高温。
物质在几千度温度下已经气化并在进一步的高温下形成等离子体-物质的第四态。
聚变堆堆心是几亿度的等离子体。
就算是实验室中的等离子体,温度也在几百万度以上。
这些等离子体必须放在高真空容器中。
通常的材料最高能经受3000 K 就很不错了。
显然,不可能用普通的材料来包围聚变堆中的等离子体。
最简单的想法是利用电磁场。
利用电场是最简单的,在历史上也有人试验过。
现在也有人继续研究用电场约束等离子体。
这对低温等离子体是可以的。
但对高温等离子体不行。
主要是因为,等离子体中的离子和电子在电场作用下沿相反的方向运动,于是外加电场很快就被极化的等离子体屏蔽了。
另外,也无法设计一种三维的封闭电场位形。
用磁场约束等离子体则有完全不同的图像。
带电粒子在磁场中的运动分成两部分,在垂直于磁力线方向做Larmor 运动;沿磁力线方向则可以自由运动(如果磁场是均匀的)。
因此,除非受到其他作用,带电粒子不会离开磁力线。
所以,磁场可以将高温等离子体与周围物质(真空室)隔开。
磁场的这种热绝缘本领与磁场的强度有关,也与等离子体的参数有关。
进一步的分析表明,更与磁场位形的特性有关。
II.1 磁约束位形 (magnetic configuration)磁场是一类无源矢量场,磁力线既不能产生,也不能消除。
磁场一定与电流有关。
对稳态磁场,描述磁场的两个电动力学方程是Ampere 定律和无散度矢量关系:j B 0μ=⨯∇,0=⋅∇B(在等离子体物理中,不区分磁感应强度和磁场的差别,即将等离子体当作真空介质处理。
所有的电流都看成真空电流。
)磁约束聚变研究近60年的历史表明,现有技术进展只能确保环型磁约束位形有可能建造聚变反应堆。
HL-2A托卡马克大破裂期间逃逸电子产生与扩散的模拟研究的开题报告
HL-2A托卡马克大破裂期间逃逸电子产生与扩散的模拟研究的开题报告1. 研究背景与意义HL-2A托卡马克是中国自主研发的核聚变装置,是国际上少数完整实现了托卡马克核聚变的研究装置之一。
在托卡马克大破裂等异常事件中,会有大量的高能粒子与辐射释放出来,其中就包括逃逸电子。
逃逸电子对人员和设备的安全构成一定威胁,因此需要进一步研究其产生和扩散规律,以制定相应的安全保护措施。
2. 研究目的本研究旨在通过数值模拟的方法,探究托卡马克大破裂事件中逃逸电子的产生和扩散规律,为安全保护提供理论和技术支持。
3. 研究方法本研究将采用数值计算方法,借助COMSOL Multiphysics软件建立HL-2A托卡马克逃逸电子扩散的三维模型,并结合实验数据,建立逃逸电子产生和扩散的数学模型。
通过模拟计算,得到逃逸电子在托卡马克中的分布和扩散规律,并分析不同因素对逃逸电子扩散的影响。
4. 预期结果本研究预期得到HL-2A托卡马克大破裂期间逃逸电子的产生和扩散规律,对逃逸电子对人员和设备的威胁进行定性和定量分析,并提出相应的安全保护措施,对保障核聚变装置的安全运行具有重要意义。
5. 参考文献[1] 张骞, 邵芝高, 王志新, 等. HL-2A 托卡马克核聚变装置物理实验进展[J]. 科学技术与工程, 2015,15(6):6-13.[2] 程帆, 袁小伟, 卢鑫, 等. 托卡马克设备中逃逸电子的研究进展[J]. 核动力工程, 2015,36(2):195-199.[3] H. D. Yang, L. Q. Hu, J. Li, et al. Electron cyclotron emission imaging system on the HL-2A tokamak [J]. Review of Scientific Instruments, 2010, 81(10): 10F119.。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
bˆ
×
∇ϕ
⋅
∇
ln
P0
i(ω − V0||k|| )
(19)
在平板位形中,κ = 0,把(19)代入(17)给出:
∂Vˆ|| ∂t
+ V0||bˆ ⋅ ∇Vˆ||
+
c B
bˆ × ∇ lnV0||
=
( ) − e
mV0||
bˆ ⋅
⎡ ⎢
⎢∇ϕ
⎢ ⎢⎣
+
P0 en0
⎜⎛ ∇⎜
⎜⎜⎝
c B
bˆ × ∇ϕ ⋅ ∇ ln i ω − k||V0||
τ =4
2
Rq
v⊥ε 0.5
K (κ ), κ 2
≈
R 2r
v||2 v⊥2
<1
c) Toroidal drift of trapped particles
vdϕ = 4
2
mv ⊥ 2
E (κ )
(
−
1 ),
qBθ R K (κ ) 2
3, Diamagnetic drift of plasma fluids
i) 宏观不稳定性和输运; ii) 波和等离子体相互作用;
iii) 微湍流和反常输运; iV) 边界层物理; V) 高能量粒子物理
4)线性理论的作用:饱和湍流幅的计算需要包括非
线性效应, 线性理论可以:
i)确定可能的驱动机制,抑制机制;
ii)确定不稳定性条件;
iii)当湍流引起的输运占主导时,等离子体的密度和
P0
⎟⎟⎞⎥⎥⎤ ⎟⎟⎠⎥⎥⎦
(20)
由方程(20)可得:
( ) ( ) V||
= i
1
ω−k||V0||
⎧
⎪⎨⎪⎪Bc bˆ×∇φ⋅∇lnV0||
+
e mV0||
bˆ⋅
⎡ ⎢
⎢∇φ
⎢
+
P0 en0
⎛ ⎜ ∇⎜ ⎜
Bc bˆ×∇φ⋅∇lnP0 i ω−k||V0||
⎞⎤⎫ ⎟⎥⎪⎪ ⎟⎥⎬ ⎟⎥⎪
计算电子漂移波频率和离子漂移波频率。
Ⅲ、离子温度梯度不稳定性的流体研究
1,基本方程: 离子连续性方程:
∂ ∂t
ni
+
∇
⋅ (n iν
i )=
0
(1)
离子运动方程 :
mi
ni
⎜⎛ ⎝
∂ ∂t
+Vi
⋅ ∇ ⎟⎠⎞Vi
=
eni ⎜⎜⎝⎛∇φ
+
Vi
×B c
⎟⎟⎠⎞
− ∇Pi
− ∇ ⋅ Πi
(2)
离子压强的绝热演化方程:
不稳定性;其空间尺度:λ → ρe − ρi
1)捕获粒子不稳定性:环状位型(托卡马克)
中: ⎜⎜⎝⎛
ν ν
||
⎟⎟⎠⎞ 2
<
r
R
的粒子被捕获在环的外侧而引起的不稳定性
(TEM)。
2)漂移型不稳定性:由粒子的密度,温度以及
磁场的空间不均匀性引起的不稳定性,是托卡马
克中的主要微观不稳定性(ITG,ETG)。
4,漂移型不稳定性:一切与漂移运动相关的不 稳定性。 i)密度梯度( ∇n )漂移: ii)温度梯度( ∇T)漂移:
iii)磁场梯度(∇B )漂移:
Vg
=
bˆ Ω
×
μm
m
∇B
=
V⊥2 2RΩ
nˆ × bˆ
iV)磁场曲率漂移:
∇||bˆ
=
nˆ R
,
Vd
=
bˆ Ω
×V||2
∂bˆ ∂l
=
V||2nˆ × bˆ RΩ
eni ⎜⎜⎝⎛ −
∇φ
+
Vi c
×
B ⎟⎟⎠⎞
−∇pi=0 Nhomakorabea(2a)
bˆ×(2a) 给出:
Vi(0)
=V||bˆ
+
c B
bˆ
×∇φ
+
c bˆ eBn
×∇pi
(6)
其中: bˆ ≡ B B
将(6)代入(2)的右端,取到漂移近似的一级,我们
得到:
Vi (1)
=
−
mi c 2 eB 2
⎜⎛ ⎝
∂ ∂t
+ Vi(0)
Q(r,t) = Q(x) exp i(ky y + kz z − ωt)
电子漂移波频率: ω
∼ ω∗e
≡
− ckθTe eBn
dn dr
=
ckθ Te eBLn
离子漂移波频率:ω
∼ ω∗i
≡
ckθ Te qi Bn
dn dr
=
−
ckθ Te qi BLn
练习题: 1,证明反磁漂移对密度扰动无贡献,讨论其物理 意义。
d dt
=
∂ ∂t
+ V0
⋅∇
=
∂ ∂t
+ V0||bˆ ⋅ ∇ +
c enB
bˆ × ∇P0
⋅∇
=
−iω
+
ik||V0||
−
P0c(1 + ηi
enBLn
)
∂ ∂y
其中:
Ln
=
−
1 n0
dn0 dx
ηi
=
Ln LTi
LTi
=− 1 Ti
dTi dx
Pi = n iT i
方程(22)可以整理为:
( ) ( ) ( ) ( ) {1 +
Ⅰ.引言
磁场位形 θ
φ
r
R
r=r2 r
r=r1 0
πr
2πr
rθ
Poloidal direction
2πR
πR
Toroidal direRcφtion
1,Tokamak magnetic
configuration
• Equilibrium magnetic field:
Toroidal field
V)极化漂移:
Vdi
=
−
mc qB2
⎜⎜⎝⎛
dE dt
×b ⎟⎟⎠⎞
×b
=
1 ΩB
dE dt
,
E⊥B
练习题:
1,什么是反磁漂移?试证明反磁漂移对 密度扰动无贡献并讨论其物理意义。 2,Ti 对∼ T典e ∼型10的ke托v,卡B ∼马3克T ,参Ln数∼ ,50cm
kθ ∼ (0.1−1)ρ,
假设
• Rotation: v⊥ = Ωρ ≈ vth
• Drifts of guiding center
i) Electric field drift:
VE
=
E
×
bˆ B
ii) magnetic gradient ( ∇B ) drift:
Vg
=
bˆ Ω
×
μm
m
∇B
=
V⊥2 2RΩ
nˆ × bˆ
iii) magnetic curvature drift:
−
e miV0||
bˆ
⎛ ⋅⎜ ⎝
∇φ
+
P0 en0
∇Pˆ
⎞ ⎟ ⎠
(17)
方程(3)化为:
∂Pˆ ∂t
+V0||bˆ
⋅
∇Pˆ
+
c B
bˆ
×∇ϕ
⋅
∇ln
P0
+
Γ∇
⋅
Vˆ
=
0(18)
方程(16-18)是基本 方程。取 ~f = fe−iωt+ik⋅r ,略去Γ∇ ⋅Vˆ 则(18)给出:
Pˆ =
c B
2,对典型的托卡马克参数,
Ti ∼ Te ∼ 10kev, B ∼ 3T , Ln ∼ 50cm
假设 kθ ∼ (0.1−1)ρ,
计算电子漂移波频率和离子漂移波频率。
3,漂移不稳定性: 在一定条件下,漂移波会增长,称为漂移不 稳定性。由密度梯度漂移引起的不稳定性称 为漂移不稳定性,又称普适不稳定性。在电 子漂移方向传播的称为电子模,在离子漂移 方向传播的称为离子模。
⎪⎩
⎢⎣
⎝
⎠⎥⎦⎪⎭
(21)
将(19)和(21)代入(16),我们得到:
( ) ∂nˆ
∂t
+ V0||bˆ
⋅ ∇ nˆ
+
i
ω
V 0 || − V0|| k||
bˆ ⋅ ∇
( ) ⎧
⎪⎪ c
⎨ ⎪
B
⎪⎩
bˆ × ∇ ϕ
⋅∇
ln V0||
+
e mV 0||
⎡
bˆ
⋅
⎢ ⎢∇
ϕ
⎢
⎢⎣
+
P0 en 0
⎜⎛ ∇⎜
传输的三种方式:传导,对流,辐射, 这里的反常输运指的
是反常热传导(在垂直于 i)经典输运:
B0
的方向)
χe = 4.66ρe2υee
χi
=
2
ρ
υ2
i ii
ii)新经典输运(香蕉区)
−3
χ ne e
~
ε
2
q 2 ρ e2υ ee
−3
χ ne i
~
ε
2
q 2 ρ i2υii
iii)由扰动引起的输运(反常输运)
=
cT jbˆ × ∇n j q jn jB
2)粒子流强:
Γ
=
∂n ∂x