中国聚变工程试验堆上新经典撕裂模和纵场波纹扰动叠加效应对alpha粒子损失影响的数值模拟

托卡马克聚变堆中高能量alpha粒子的良好约束是获得稳态燃烧等离子体的前提,除了磁场波纹损失,不稳定性也会引起额外的损失.本文基于中国聚变工程试验堆(CFETR)参数, alpha粒子初始分布和新经典撕裂模(NTM)扰动分布,利用粒子导心跟踪方法分别对磁场波纹和NTM两种扰动及叠加下的粒子损失进行了详细的数值模拟.结果显示粒子损失份额不随NTM扰动幅度增大而增大,两种扰动的叠加效应不明显.通过扫描装置波纹度大小以及分析相空间粒子密度和波纹损失区分布,确认原因是CFETR波纹损失区较小,没有覆盖高能量粒子(EP)密度和NTM扰动主要分布区.此外, NTM没有引起直接粒子损失和俘获粒子剖面坍塌.显著的俘获粒子密度剖面展平,并扩展到波纹损失区是两种扰动叠加效应显著的前提.无碰撞波纹随机扩散是CFETR初始分布alpha粒子的主要损失通道,通行粒子约束不受磁场波纹影响.本文研究结果对CFETR概念设计中alpha粒子物理和低频不稳定性下的EP行为具有重要意义.     

托卡马克中局域磁扰动引起的高能量离子损失研究

基于单粒子导心运动代码ORBIT，采用测试粒子模拟方法，研究了托卡马克等离子体内部不同径向位置处局域磁场扰动对高能量离子的损失的影响。研究表明，在局域磁扰动主要分布在某磁面附近、其环向具有类似纹波场形式下，可造成一些靠近等离子体中心区域的高能量离子损失，但对靠近等离子边界的离子损失影响相对不大。这些损失的高能量离子均为捕获离子，离子的投掷角越大就越容易损失。此外，造成高能量离子最大损失率的局域场径向位置与这些损失离子的初始径向位置通常存在一定的偏移，而且这个偏移与这些离子的能量密切相关。当局域场出现在某些位置时，能量较低的离子会有一定的损失，能量较高的离子反而不会损失。     

托卡马克等离子体中共振磁扰动场放大效应对离子轨道特性的作用

首先,利用MARS-F程序模拟了HL-2A装置中不同电阻率、环向旋转频率和环向模数条件下等离子体对共振磁扰动(RMP)的线性响应过程,分析了不同情况下共振场放大效应对三维扰动磁场的影响;然后,利用Boris算法追踪了三维场作用下的离子轨道,并详细探究了不同扰动磁场改变离子轨道特性的物理机制.研究发现,考虑等离子体响应后的扰动磁场可增强离子轨道径向展宽,且轨道最大径向展宽随轨道上扰动磁场平均值的增大而增大;同时,离子在通过扰动磁场被强烈放大的区域时轨道径向展宽会显著增加.该物理机制可用于解释RMP缓解边界局域模实验中离子直接损失增加和等离子体径向输运增强现象.     

聚变等离子体中高能量粒子驱动的鱼骨模不稳定性研究

在考虑高能量粒子密度的空间分布及箍缩角分布的条件下,建立了研究鱼骨模的色散关系,并作了数值研究。结果表明：当高能量粒子的密度高于某个阈值时,鱼骨模会在高能量粒子密度梯度最大处被激发,其频率与高能量粒子的环向进动频率一致。在高β_h 区间,高能量勉强通行粒子将驱动鱼骨模进入第二稳定区。高能量捕获粒子能激发非共振鱼骨模,与勉强通行粒子激发的鱼骨模类似,在高β_h 区存在第二稳定区。     

环位形下鱼骨模驱动快离子输运的数值模拟研究

考虑包含动理学效应的鱼骨模结构,使用导心轨道程序 ORBIT,在磁面坐标下研究了不同的扰动模 幅度、频率对快离子再分布的影响,并分析了粒子与扰动发生共振的条件。模拟得出,鱼骨模扰动会使快离子在 实空间、相空间中发生再分布,芯部( ψ_p/ψ _w≤0.2 ,ψ_p 为极向磁通,ψ_w 为最后一个闭合磁面的磁通)快离子密 度下降约 20%,中间磁面位置上( ψ_p/ψ _w≥0.2,ψ_p/ψ _w≤ 0.6)的快离子密度增加约 7%;通过扫描频率发现,相空间中快 离子的再分布对模式频率敏感,并分析了快离子与鱼骨模扰动共振的条件。     

CFETR阿尔法粒子波纹损失的数值模拟

使用导心轨道程序ORBIT,在平衡程序EFIT给出的中国聚变工程实验堆 (CFETR)平衡位型下,结合不同的阿尔法(α)粒子分布模型,计算了氘氚聚变产生的α粒子波纹损失情况。计算结果表明：在不考虑锯齿模不稳定性的α粒子分布下,ITER-like和super-X位型下的α粒子波纹损失份额为0.1%,snowflake位型在0.4%,反磁剪切位形在0.6%;在较平缓的α粒子分布下,损失份额增大,损失的高能量α粒子有局域性。     

电子回旋共振放电的电离特性PIC/MCC模拟(Ⅱ)——数值模拟与结果讨论

采用粒子模拟与蒙特卡罗相结合(PIC/MCC)的方法，应用电磁模型，编写了准三维的电子回旋共振(ECR)放电电离过程的模拟程序，得到了ECR放电过程中电子与离子的相空间分布、电磁场分布.通过对这些分布随时间演化的分析，得出ECR加热发生在ω≈ω_c0且垂直于轴向的区域；ECR区域，微波能量几乎全部耦合给电子，获得能量的电子通过与中性粒子的电离碰撞产生了大量的带电粒子；随着放电的进行，大量带电粒子通过频繁的碰撞，分布由各向异性逐渐趋于各向同性. 关键词： 电子回旋共振放电 粒子模拟 蒙特卡罗 电离     

基于双相位编码联合变换相关器的图像位移探测算法

提出了基于双相位编码联合变换相关器(JTC)的图像位移探测算法.随机产生一相位码对参考图像编码,与目标图像叠加作为输入图像,通过傅里叶变换得到联合功率谱,用同样的相位码对其编码并进行滤波,傅里叶逆变换后得到涉及位移矢量的互相关峰.给出了多帧图像间位移矢量的探测技术.相比传统的相关器探测技术,此算法能更有效地利用输入平面...     

共轴双涡旋声束形成离轴多涡旋的定位

涡旋声束具有螺旋的相位波前,中轴线上形成声强为零的相位奇点,其所携带的轨道角动量在粒子操控领域有着广阔的应用前景。传统声涡旋只在传播轴线上形成一个拓扑荷可控的涡旋波束,这限制了声涡旋的应用灵活性。基于环形点声源阵列和相位编码技术,利用奇偶声源分别产生共轴双涡旋声束的声场叠加,在传播截面上形成了具有中心涡旋和子涡旋的离轴多涡旋声场;研究了双涡旋拓扑参数对离轴涡旋的个数、位置及拓扑荷的影响,基于声涡旋的径向声压和相位分布,确定了离轴涡旋的离轴半径,并结合声源位置推导子涡旋中心方位角的计算公式,实现离轴涡旋的精确定位。本研究突破了沿轴分布的涡旋声场只能形成单点涡旋势阱的操控局限,为利用离轴多涡旋实现多点粒子捕获提供了理论依据,促进涡旋声场在精确粒子操控和传输方面的高效应用。      

HL-2A先进二维硬X射线成像CCD系统设计

为HL-2A装置设计了一种二维硬X射线成像CCD系统,用它能直接测量硬X射线辐射强度的二维时空分布I(z,y,t),并能得到电子速率分布fe(v,r,t)。为获得二级加热和无感电流驱动如ECH／ECCD,LHH／LHCD,ICRH／ICCD,NBIH／NBICD的效果、波可近性、共振层位置、能量沉积范围等方面的实验研究提供重要诊断手段;也可用于内破裂、大破裂、MHD,特别是q=2附近高模的MHD行为、热电子和超热电子的输运、被捕获的高能粒子图像以及与模式波(如m=1,n=1)共振引起的鱼骨模不稳定性、辅助加料引起和诱发逃逸电子雪崩现象的实验观察与研究。     

Fractional Resonances between Waves and Energetic Particles in Tokamak Plasmas

From numerical simulation and analytical modeling it is shown that fast ions can resonate with plasma waves at fractional values of the particle drift-orbit transit frequency when the plasma wave amplitude is sufficiently large. The fractional resonances, which are caused by a nonlinear interaction between the particle orbit and the wave, give rise to an increased density of resonances in phase space which reduces the threshold for stochastic transport. The effects of the fractional resonances on spatial and energy transport are illustrated for an energetic particle geodesic acoustic mode but they apply equally well to other types of MHD activity.     

Energetic Particle Transport Prediction for CFETR Steady State Scenario Based on Critical Gradient Model

The critical gradient mode(CGM) is employed to predict the energetic particle(EP) transport induced by the Alfven eigenmode(AE).To improve the model,the normalized critical density gradient is set as an inverse proportional function of energetic particle density;consequently,the threshold evolves during EP transport.Moreover,in order to consider the EP orbit loss mechanism in CGM,ORBIT code is employed to calculate the EP loss cone in phase space.With these improvements,the AE enhances EPs radial transport,pushing the particles into the loss cone.The combination of the two mechanisms raises the lost fraction to 6.6%,which is higher than the linear superposition of the two mechanisms.However,the loss is still far lower than that observed in current experiments.Avoiding significant overlap between the AE unstable region and the loss cone is a key factor in minimizing EP loss.     

托卡马克装置中高能离子的直接损失

通过严格求解导心坐标系下的哈密顿方程,提出了托卡马克装置中离子的真实轨道理论,并利用此理论对国际热核聚变实验堆(ITER)中α离子的真实轨道进行了数值模拟研究,获得了驻点轨道和近期理论预言的半俘获轨道.根据此真实轨道理论计算了ITER中α离子的直接损失率.结果发现,与所提出的真实轨道理论相比,以前的回旋平均理论大幅度低估了α离子的直接损失率,两种理论的损失率差值可达14%,对此差异进行了详细的分析并给出了物理上的解释.     

中型tokamak中快离子的约束

利用三个运动常数发展了导心轨道程序GCORBIT. 该程序在计算轴对称tokamak中带电粒子导心轨道的同时，也考虑了纵场纹波和螺距角散射的效应. GCORBIT可以用来计算轴对称tokamak中的导心轨道，计算速度空间的损失区，计算快离子的损失份额. 与已有的工作相比，GCORBIT计算出的速度空间损失区包括了多种效应：偏滤器位形中的非闭合轨道、由于导心轨道与第一壁相交引起的首次轨道损失、纵场纹波局部磁阱俘获、纹波随机扩散以及速度空间不连续约束区导致的螺距角散射效应增强. GCORBIT程序已被用于分析“先进超导托卡马克实验装置”(experimental advanced superconducting tokamak,EAST) tokamak中快离子的约束. 并讨论其中的某些数值计算结果. 关键词： 导心轨道程序 快离子约束 tokamak     

Near fields scattered by an underwater finite cylindrical baffle

In this work,acoustic vector characteristics of near fields scattered by an underwater finite cylindrical baffle are investigated theoretically and experimentally.The analytic expressions for the scattered pressure and particle velocity are derived using the elastic thin shell theory.Calculations are presented for the scattered near fields of the pressure,the particle velocity and the intensity.It is found that the pressure and the particle velocity fields near the surface of the cylindrical baffle are characterized by complex interference structure,particle velocity directions and the source bearings are not consistent.The phase difference between the pressure and the particle velocity is not zero and the intensity vector does not reflect the sound bearings.It can be noted that the distortions of the fields will make the original vector signal processing method based on the free space assumption be no longer applicable in the presence of the cylindrical baffle.These results can serve as a basis of the application for the acoustic vector sensor on board.     

Numerical study of fast ion behavior in the presence of magnetic islands and toroidal field ripple in the EAST tokamak

A peculiar first orbit loss type was found apart from the normal ones when we use ORBIT code to simulate fast ion orbits in the EAST tokamak. Fast ion orbits were studied in the presence of toroidal field （TF） ripple and magnetohydro- dynamic （MHD） perturbations. We analyzed the properties of the drifted orbits in detail and compared their differences, finding that the combined effects of ripple and magnetic islands are much greater than the effects of either one of them alone. Then we investigated the orbitdeviations as a function of pitch angle in different radial positions. The modeling results demonstrate that the loss of trapped particles is mainly caused by the ripple, while MHD＇perturbation mainly plays an important role in the passing particles. Furthermore we modeled the loss rate using different equilibriums. Results prove that a higher beta can indeed improve the confinement of fast ions, while a little change in the q profile can make the topologies of magnetic islands become quite different and results in quite different total particle losses.     

�п�����о�����Ŷ�����ĸ�����������ʧ�о�

Based on ORBIT code of a guiding center motion of single particle, the loss of energetic ions in different radial positions of tokamak plasma is studied by using test particle simulation method. The results show that the local magnetic perturbations can cause loss of many energetic ions near the central region of the plasma, but they have little effect on the ion loss near the plasma boundary, assuming that the local field is mainly located near a magnetic surface and its toroidal field is similar to the ripple field. These energetic ions are trapped ions, and the greater their pitch angle is, the easier they lose. In addition, the radial position of the local field that causes the maximum loss rate of energetic ions is usually offset from the initial radial position of these loss ions, and this shift is closely related to the energy of these ions. When the local field appears in certain locations, the ions of lower energy have some loss, but the ions higher energy does not lose.     

Three-wave interactions between fast-ion modes in the National Spherical Torus Experiment

Simultaneous bursts of energetic particle mode (EPM) and toroidicity-induced Alfvén eigenmode (TAE) activity that correlate with significant fast-ion loss are observed in beam heated plasmas. Three-wave interactions between these modes are conclusively identified, indicating fixed phase relationships. This nonlinear coupling concentrates the energy of the TAEs into a toroidally localized perturbation frozen in the frame of a rigid, toroidally rotating structure formed by the EPMs. This redistribution of energy is significant because it will modify the effect of the TAEs on fast-ion loss.     

ELECTRON ACCELERATION BY A LOCALIZED,HLGH FREQUENCY ELECTRIC FIELD OCDDRRING IN LASER PLASMA

The problem of electron acceleration by a localized,high frequency electric field has been further investigated by using single particle orbit computations and the resonant field. The computational results obtained can give a better explanation to the phenomena observed in experiments and in particle simulations,it indicates that,the Localized field occurring in resonant absorption behaves more Like a travelling wave, and the resonance of the field with electrons in the critical density region can account for the generation of energetic electrons.