EAST上由垂直不稳定性引发破裂的分析与预测

从EAST 装置2016 年的放电实验中，选取了119 次等离子体破裂放电数据，分析诱发等离子体破裂的原因，发现约60%的破裂是由垂直不稳定性直接引起的，其破裂后将会产生更大的晕电流，从而产生更大的电磁应力损坏装置。对由垂直不稳定性引起的破裂(简称为VID)(72 次放电)进行了研究，建立了分别基于单变量(垂直位移)和两维变量(垂直位移、垂直位移增长率)的预测模型用于对VID 破裂的预测。离线测试表明，基于两维变量的预测模型可以在破裂发生前20ms 给出破裂预警信号，预测成功率达93%。     

HL-2M等离子体垂直不稳定性的反馈控制研究

使用无源稳定导体和主动控制线圈来控制HL-2M等离子体的垂直不稳定性。计算了等离子体垂直不稳定性增长时间,构建了等离子体垂直不稳定性控制的线性模型,然后用MATLAB对采用PID算法的垂直不稳定性控制系统进行了模拟仿真。结果表明,无源稳定导体和主动控制线圈能够快速稳定等离子体的垂直不稳定性运动,这也表明设计是可行的。     

Synergetic Effects of Runaway Electron and Vertical Instability on the Wall Damage on HL-2A Tokamak

1 Generation of vertical instability Tokamak experiments show that the energy confinement time and performance are better, and the larger plasma current can be achieved for non-circular cross-sectional shape of plasma than circular cross-section. However, the external magnetic fields which are used to produce the non-circular cross sectional shape also cause the confined plasma to become unstable to small vertical displacement. In general, the ratio of elongation k is larger, the possibility of instability is more. In practice, this vertical displacement mode stabilized by feedback control system or other provided external radial magnetic field to balance out the plasma motion. Under some of situations the control system may be fail due to rapid growth rate of instability exceeding the ability of controlling. The plasma will then move vertically upwards or downwards depending upon the characteristics of instability and control failure.     

EASTװ�ô�ֱ���ȶ����о�

对EAST装置偏滤器位形的放电中发生的垂直位移不稳定性的全过程进行了模拟,计算了磁轴垂直位移及垂直位移不稳定性的增长率。模拟的结果与实验结果符合得较好。     

磁驱动飞片发射实验结构系数初步研究

为了确定磁驱动飞片发射实验结构系数的范围、影响因素、结构系数与影响因素的关系，对聚龙一号装置上的磁驱动飞片发射实验进行了数值模拟和分析。数值模拟表明，磁流体力学程序能正确模拟聚龙一号装置上各个磁驱动飞片发射实验；磁驱动双侧飞片发射实验的结构系数为0.7～0.8；磁驱动单侧飞片发射实验的结构系数为0.80～0.85。磁驱动飞片发射实验的结构系数与实验加载电流无关，仅由磁驱动飞片发射实验的负载结构决定。磁驱动飞片发射实验的结构系数取决于阴阳电极极板的初始宽度、阴阳电极之间的初始间隙以及阴阳电极上飞片厚度之和等三个因素。在磁驱动飞片发射实验中，电极初始宽度、阴阳电极之间的初始间隙不变的情况下，结构系数由阴阳电极上飞片厚度之和确定，阴阳电极上飞片厚度之和越大，结构系数越大。     

A method for calculating active feedback system to provide vertical position control of plasma in a tokamak

In designing tokamaks, the maintenance of vertical stability of plasma is one of the most important problems. Systems of the passive and active feedbacks are applied for this purpose. Role of the passive system consisting of a vacuum vessel and passive coils is to suppress fast MHD (magnetohydrodynamic) instabilities. The active feedback system is applied to control slow motions of plasma. The objective of the paper is to investigate two successive problems, solution of which allows to determine the possibility of controlling plasm a motions. One of these is the problem of vertical stability under the assumption of ideal conductivity of plasma and passive stabilizing elements. The problem is solved analytically and on the basis of the obtained solution a criterion of MHD-stability is formulated. The other problem is connected with the control of plasma vertical position with active feedback system. Calculation of feedback control parameters is formulated as an optimization problem and an approximate method to solve the problem is suggested. Numerical simulations are performed with parameters of the T-15M tokamak in order to justify the suggested method.      

约束条件下胶球间排空力的研究

用MonteCarlo方法对处于两平行硬板约束下三个浓度的大小胶球系统进行了模拟,通过对大胶球表面小胶球密度的统计,由密度积分公式获得了大胶球所受的排空力.研究结果显示,因为平行硬板的存在或当改变两平行硬板的距离时,同浓度下,排空力在硬板距离小的时候最明显;三个浓度中,浓度高的,排空力受硬板距离影响最大;有硬板约束比无该约束的时候,排空力效果更显著.     

HL-2M�������崹ֱ���ȶ����о�

用DⅢ-D的TokSys研究了HL-2M等离子体垂直不稳定性控制问题.使用真空室和主动线圈控制等离子体垂直不稳定性,分析了真空室的本征模和电流分布,极大简化真空室建模难度,建立了垂直不稳定性数学模型,验证了主动线圈和电源模型参数,计算了不同拉长比位形的等离子体增长率和最大可控垂直位移.利用TokSys建立了雪花、单零、双零偏滤器位形垂直不稳定仿真模型,然后通过simulink对采用PD算法的垂直不稳定系统进行仿真.结果表明,构建的模型能够较好控制不同位形的垂直不稳定性.     

HL-2M等离子体垂直不稳定性研究

用DⅢ-D的TokSys研究了HL-2M等离子体垂直不稳定性控制问题。使用真空室和主动线圈控制等离子体垂直不稳定性,分析了真空室的本征模和电流分布,极大简化真空室建模难度,建立了垂直不稳定性数学模型,验证了主动线圈和电源模型参数,计算了不同拉长比位形的等离子体增长率和最大可控垂直位移。利用TokSys建立了雪花、单零、双零偏滤器位形垂直不稳定仿真模型,然后通过simulink对采用PD算法的垂直不稳定系统进行仿真。结果表明,构建的模型能够较好控制不同位形的垂直不稳定性。     

Instabilities of cylindrical bubble clusters

Small bubbles in an experimental two-dimensional foam between glass plates regularly undergo a three-dimensional instability as the small bubbles shrink under diffusion or equivalently as the plate separation increases, and end up on one of the plates. The most recent experiments of Cox, Weaire, and Vaz are accompanied by Surface Evolver computer simulations and rough theoretical calculations. We show how a recent second variation formula may be used to perform exact theoretical calculations for infinitesimal perturbations for such a system, and verify results with Surface Evolver simulations.     

Modulation instability of an intense laser beam in an unmagnetized electron–positron–ion plasma

The modulation instability of an intense circularly polarized laser beam propagating in an unmagnetized, cold electron–positron–ion plasma is investigated. Adopting a generalized Karpman method, a three-dimensional nonlinear equation is shown to govern the laser field. Then the conditions for modulation instability and the temporal growth rate are obtained analytically. In order to compare with the usual electron–ion plasmas, the effect of positron concentration is considered. It is found that the increase in positron-to-electron density ratio shifts the instability region towards higher vertical wave numbers but does not cause displacement along the parallel wave number direction, and the growth rate increases as the positron-to-electron density ratio increases.     

Suppression of the Rayleigh-Taylor instability due to self-radiation in a multiablation target

A scheme to suppress the Rayleigh-Taylor instability has been investigated for a direct-drive inertial fusion target. In a high-Z doped-plastic target, two ablation surfaces are formed separately-one driven by thermal radiation and the other driven by electron conduction. The growth of the Rayleigh-Taylor instability is significantly suppressed on the radiation-driven ablation surface inside the target due to the large ablation velocity and long density scale length. A significant reduction of the growth rate was observed in simulations and experiments using a brominated plastic target. A new direct-drive pellet was designed using this scheme.     

Plasma position control in SST1 tokamak

For long duration steady state operation of SST1, it would be very crucial to maintain the plasma radial and vertical positions accurately. For designing the position controller in SST1 we have adopted the simple linear RZIP control model. While the vertical position instability is slowed down by a set of passive stabilizers placed closed to the plasma edge, a pair of in-vessel active feedback coils can adequately control vertical position perturbations of up to 1 cm. The shifts in radial position arising due to minor disruptions would be controlled by a separate pair of poloidal field (PF) coils also placed inside the vessel, however the controller would ignore fast but insignificant changes in radius arising due to edge localised modes. The parameters of both vertical and radial position control coils and their power supplies are determined based on the RZIP simulations.     

Spiral growth and step edge barriers

The growth of spiral mounds containing a screw dislocation is compared to the growth of wedding cakes by two-dimensional nucleation. Using phase field simulations and homoepitaxial growth experiments on the Pt(111) surface we show that both structures attain the same large scale shape when a significant step-edge barrier suppresses interlayer transport. The higher vertical growth rate of the spiral mounds on Pt(111) reflects the different incorporation mechanisms for atoms in the top region and can be formally represented by an enhanced apparent step-edge barrier.     

The Effective Mode Approximation for Analyzing Position Instability of a Tokamak Plasma

The effect of induced currents in the equilibrium field (EF) coils on ameliorating the instability of a small position perturbation of a rigid tokamak plasma is analyzed. A strong analogy between the position instability of the plasma and supercriticality of a fission reactor core is recognized. The position of the plasma corresponds to the neutron population, and the retardation of the displacement by the induced eddy currents to the suppression of population growth by delayed neutrons. The matrix equation of the dispersion relation for the position instability is diagonalized and factored into a form identical to the in-hour equation of fission reactor kinetics. An Effecitve Mode Approximation (EMA) similar to the one group of delayed neutron approximation has been introduced to greatly simplify the analysis of the position instability and feedback control. With this approximation the dispersion relation is reduced to a linear or cubic algebraic equation depending on the effectiveness of retardation by the eddy currents. The time constant of the unstable mode can be expressed in terms of the plasma parameters and the effective resistance and inductance of the current carriers, which can be conveniently computed. The vertical instability of a typical noncircular tokamak plasma is analyzed numerically as well as analytically by the EMA method. The results agree well within a negligible discrepancy.     

Statistics of defect trajectories in spatio-temporal chaos in inclined layer convection and the complex Ginzburg-Landau equation

For spatio-temporal chaos observed in numerical simulations of the complex Ginzburg-Landau equation (CGL) and in experiments on inclined-layer convection (ILC) we report numerical and experimental data on the statistics of defects and of defect loops. These loops consist of defect trajectories in space-time that are connected to each other through the pairwise annihilation or creation of the associated defects. While most such loops are small and contain only a few defects, the loop distribution functions decay only slowly with the quantities associated with the loop size, consistent with power-law behavior. For the CGL, two of the three power-law exponents are found to agree, within our computational precision, with those from previous investigations of a simple lattice model. In certain parameter regimes of the CGL and ILC, our results for the single-defect statistics show significant deviations from the previously reported findings that the defect dynamics are consistent with those of random walkers that are created with fixed probability and annihilated through random collisions.     

Theoretical analysis of polarization control for the stable output of multi-carrier source based on a re-circulating frequency shifter

Fluctuating polarization state-of-light in the optical loop is an important factor that seriously influences the output performance of a multi-carrier source based on re-circulating frequency shifter (RFS). The reason for output spectrum instability when no polarization controller (PC) is present in the loop is analyzed theoretically. Numerical simulations for the output spectra of the multi-carrier source with and without PC are conducted, and the trajectories of the several frequency components polarization states on the Poincare sphere with and without PC are compared. The results show that the performance of multi-carrier source based on a RFS can be improved effectively by adjusting the PC in the configuration properly.     

大气压力对激光辐照双层板接触传热的影响

考虑激光辐照下结构变形对双层板接触传热影响可以更真实地反映两板间的传热情况，通过实验和数值模拟分析了大气压力对接触传热的影响。结果表明:大气压力对接触热阻影响非常明显，当大气压力超过一定值后，双层板界面始终接触，后板中心的温升阶段成类抛物形状；当大气压力在某一范围内，大气压力与温度矩产生的靶板挠度相当，接触面时分时合，后板的温升阶段成振荡式类线性增加；当大气压力小于某一特定值时，后板温度有一突跃过程，通过适当的设计有可能观测到该现象。     

HT-6A托卡马克平衡系统的实验

本文叙述了HT-6A托卡马克平衡系统及第一期平衡实验结果.讨论了HT-6A上铜壳对平衡的贡献;观察到改变垂直场形态因子n而产生的水平位移及垂直位移不稳定现象;介绍了正常放电时的平衡实验,观察到等离子体位移与MHD扰动是紧密相连的;用沙弗拉诺夫薄壳理论计算了对多大的等离子体电流应施加多大的垂直场,结果与实验完全符合。     

Off-peak saturation effects of beam-plasma instability

For beam-plasma instability in the absence of a magnetic field, hybrid simulations produce almost identical results to those of particle simulations (both agree with linear and nonlinear analyses), but with much reduced computing cost and noise level. With the verification by simulations, it is shown that the saturation level is very strongly influenced by the discrete wave number spectrum, sometimes giving different, by an order of magnitude, results by choosing slightly different off-peak (in growth rate curve) parameters, which must be considered in most simulations (whether hybrid or particle or whether magnetized or unmagnetized)