托卡马克中的静电漂移不稳定性

合集下载

托卡马克中的静电漂移不稳定性

托卡马克中的静电漂移不稳定性

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)核聚变

受控热核聚变能的研究主要有两种--惯性约束核聚变和磁约束核聚变。前者利用超高强度的激光在极短的时间内辐照氘氚靶来实现聚变,后者则利用强磁场可很好地约束带电粒子的特性,将氘氚气体约束在一个特殊的磁容器中并加热至数亿摄氏度高温,实现聚变反应。
托卡马克(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]

托卡马克物理基础-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装置优缺点简析

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在极向场线圈上感应的高电压将一直传递到电流引线箱,如果在超导电流传输线和电流引线箱内绝缘薄弱、损坏和真空度下降均会引发起弧、放电和因此烧毁部件。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。


×
∇ϕ


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

×∇φ
+
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

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||

⎛ ⋅⎜ ⎝
∇φ
+
P0 en0
∇Pˆ
⎞ ⎟ ⎠
(17)
方程(3)化为:
∂Pˆ ∂t
+V0||bˆ

∇Pˆ
+
c B

×∇ϕ

∇ln
P0
+
Γ∇


=
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||



⎢ ⎢∇
ϕ

⎢⎣
+
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
相关文档
最新文档