Tokamak中的Kelvin-Helmholtz不稳定性

<正> Tokamak中的一个重要问题是加热。中性束注入加热是加热的一个有效手段,它使美国PLT上的离子温度达到7.1KeV.但PLT上的中性束注入的不对称性引起等离子体的快速环向旋转,转速可达1×10~7厘米/秒。1979年5月Suckewer等在PLT上测量了速度分布。 在具有速度剪切进行旋转的等离子体中,会不会形成新的磁流体力学不稳定性?1980年

Fusion energy-production from a deuterium-tritium plasma in the jet tokamak

Describes a series of experiments in the Joint European Torus (JET), culminating in the first tokamak discharges in deuterium-tritium fuelled mixture. The experiments were undertaken within limits imposed by restrictions on vessel activation and tritium usage. The objectives were: (i) to produce more than one megawatt of fusion power in a controlled way; (ii) to validate transport codes and provide a basis for accurately predicting the performance of deuterium-tritium plasmas from measurements made in deuterium plasmas; (iii) to determine tritium retention in the torus systems and to establish the effectiveness of discharge cleaning techniques for tritium removal; (iv) to demonstrate the technology related to tritium usage; and (v) to establish safe procedures for handling tritium in compliance with the regulatory requirements. A single-null X-point magnetic configuration, diverted onto the upper carbon target, with reversed toroidal magnetic field was chosen. Deuterium plasmas were heated by high power, long duration deuterium neutral beams from fourteen sources and fuelled also by up to two neutral beam sources injecting tritium. The results from three of these high performance hot ion H-mode discharges are described: a high performance pure deuterium discharge; a deuterium-tritium discharge with a 1% mixture of tritium fed to one neutral beam source; and a deuterium-tritium discharge with 100% tritium fed to two neutral beam sources. The TRANSP code was used to check the internal consistency of the measured data and to determine the origin of the measured neutron fluxes. In the best deuterium-tritium discharge, the tritium concentration was about 11% at the time of peak performance, when the total neutron emission rate was 6.0 × 10¹⁷ neutrons/s. The integrated total neutron yield over the high power phase, which lasted about 2 s, was 7.2 × 10¹⁷ neutrons, with an accuracy of ±7%. The actual fusion amplification factor, Q_DT was about 0.15

Optimization of positions and currents of tokamak poloidal field coils using genetic algorithms

Plasma equilibrium geometry has a great influence on the confinement and magnetohydrodynamic stability in tokamaks. The poloidal field (PF) system of a tokamak should be optimized to support the prescribed plasma equilibrium geometry. In this paper, a genetic algorithm-based method is applied to solve the optimization of the positions and currents of tokamak PF coils. To achieve this goal, we first describe the free-boundary code EQT Based on the EQT code, a genetic algorithm-based method is introduced to the optimization. We apply this new method to the PF system design of the fusion-driven subcritical system and plasma equilibrium geometry optimization of the Experimental Advanced Superconducting Tokamak (EAST). The results indicate that the optimization of the plasma equilibrium geometry can be improved by using this method.

Development of the lanzhou heavy ion linac

The heavy ion linac in Lanzhou is designed as a future injector for the Cooling Storage Ring (CSR). In order to keep the total machine within 40 meters, the IH (Interdigital H-type) structure is adopted for its higher acceleration gradient compared with the traditional DTL structure. The designed minimum charge over mass ratio is 1/6, the output energy is 16MeV/u and the beam current is 1A.mu A. The RFQ and the first DTL tank will work at 100MHz, and the other DTL tanks will work at the double frequency. The design criteria, main parameters and the detailed beam dynamic design are introduced in this paper.

Preliminary study of a niobium quarter-wave prototype cavity for a heavy-ion superconducting linac

Superconducting quarter-wave resonators, due to their compactness and their convenient shape for tuning and coupling, are very attractive for low-beta beam acceleration. In this paper, two types of cavities with different geometry have been numerically simulated: the first type with larger capacitive load in the beam line and the second type of lollipop-shape for 100 MHz, beta=0.06 beams; then the relative electromagnetic parameters and geometric sizes have been compared. It is found that the second type, whose structural design is optimized with the conical stem and shaping drift-tube, can support the better accelerating performance. At the end of the paper, some structural deformation effects on frequency shifts and appropriate solutions have been discussed.

A Particle-in-cell scheme of the RFQ in the SSC-Linac

A 52 MHz Radio Frequency Quadrupole (RFQ) linear accelerator (linac) is designed to serve as an initial structure for the SSC-Linac system (injector into Separated Sector Cyclotron). The designed injection and output energy are 3.5 keV/u and 143 keV/u, respectively. The beam dynamics in this RFQ have been studied using a three-dimensional Particle-In-Cell (PIC) code BEAMPATH. Simulation results show that this R,FQ structure is characterized by stable values of beam transmission efficiency (at least 95%) for both zero-current mode and the space charge dominated regime. The beam accelerated in the RFQ has good quality in both transverse and longitudinal directions, and could easily be accepted by Drift Tube Linac (DTL). The effect of the vane error and that of the space charge on the beam parameters have been studied as well to define the engineering tolerance for RFQ vane machining and alignment.

Design optimization of the APF DTL in the SSC-linac

A linear accelerator as a new injector for the SSC (Separated Sector Cyclotron) of the HIRFL (Heavy ton Research Facility Lanzhou) is being designed. The DTL (Drift-Tube-Linac) has been designed to accelerate U-238(34+) from 0.140 MeV/u to 0.97 MeV/u. To the first accelerating tank which accelerates U-238(34+) to 0.54 MeV/u, the approach of Alternating-Phase-Focusing (APF) is applied. The phase array is obtained by coupling optimization software Dakota and beam optics code LINREV. With the hybrid of Multi-objective Genetic Algorithm (MOGA) and a pattern search method, an optimum array of asynchronous phases is determined. The final growth, both transversely and longitudinally, can meet the design requirements. In this paper, the deign optimization of the APF DTL is presented.

Direct Numerical Simulation of Passive Scalar in Decaying Compressible Turbulence

The passive scalars in the decaying compressible turbulence with the initial Reynolds number (defined by Taylor scale and RMS velocity) Re=72, the initial turbulent Mach numbers (defined by RMS velocity and mean sound speed) Mt=0.2-0.9, and the Schmidt numbers of passive scalar Sc=2-10 are numerically simulated by using a 7th order upwind difference scheme and 8th order group velocity control scheme. The computed results are validated with different numerical methods and different mesh sizes. The Batchelor scaling with k(-1) range is found in scalar spectra. The passive scalar spectra decay faster with the increasing turbulent Mach number. The extended self-similarity (ESS) is found in the passive scalar of compressible turbulence.

圆截面轴对称等离子体环平衡的分析解和极限比压〈β_r〉

本文根据Mercier的基本方程用反纵横比ε的暴级数展开方法求出了胖Tokamak的平衡分析解。为了考查ε的影响,本文计算到了ε~3的量级,并在没有负压的限制下求得了相应的极限比压〈β_r〉。求得的平衡解亦可作为不稳定分析之用。

Jet:progress in performance and understanding

In 1990 JET operated with a number of technical improvements which led to advances in performance and permitted the carrying out of experiments specifically aimed at improving physics understanding of selected topics relevant to the "NEXT STEP". The new facilities include beryllium antenna screens, a prototype lower hybrid current drive system, and modification of the NI system to enable the injection of He-3 and He-4. Continued investigation of the hot-ion H-mode produced a value of n(D)(0)tau-E(T)(i)(0) = 9 x 10(20)m-3s keV, which is near conditions required for Q(DT) = 1, while a new peaked density profile H-mode was developed with only slightly lower performance. Progress towards steady state operation has been made by achieving ELMy H-modes under certain operating conditions, while maintaining good tau-E values. Experimental simulation of He ash transport indicates effective removal of alpha-particles from the plasma core for both L and H mode plasmas. Detailed analyses of particle and energy transport have helped establish a firmer link between particle and energy transport, and have suggested a connection between reduced energy transport and reversed shear. Numerical and analytic studies of divertor physics carried out for the pumped divertor phase of JET have helped clarify the key parameters governing impurity retention, and an intensive model validation effort has begun. Experimental simulation of alpha-particle effects with beta-fast up to 8% have shown that the slowing down processes are classical, and have given no evidence of deleterious collective effects.

槽道湍流的直接数值模拟

通过直接数值模拟(DNS)研究槽道湍流的性质和机理。包含五个部分：1）湍流直接数值模拟的差分方法研究。2）求解不可压N－S方程的高效算法和不可压槽道湍流的直接数值模拟。3）可压缩槽道湍流的直接数值模拟和压缩性机理分析。4）“二维湍流”的机理分析。5）槽道湍流的标度律分析。1．针对壁湍流计算网格变化剧烈的特点，构造了基于非等距网格的的迎风紧致格式。该方法直接针对计算网格构造格式中的系数，克服了传统方法采用 Jacobian 变换因网格变化剧烈而带来的误差。针对湍流场的多尺度特性分析了差分格式的精度、网格尺度与数值模拟能分辨的最小尺度的关系，给出不同差分格式对计算网格步长的限制。同时分析了计算中混淆误差的来源和控制方法，指出了迎风型紧致格式能很好地控制混淆误差。2．将上述格式与三阶精度的Adams半隐格式相结合，构造了不可压槽道湍流直接数值模拟的高效算法。该算法利用基于交错网格的离散形式的压力Poisson方程求解压力项，避免了压力边界条件处理的困难。利用FFT对方程中的隐式部分进行解耦，解耦后的方程采用追赶法（LU分解法）求解，大大减少了计算量。为了检验该方法，进行了三维不可压槽道湍流的直接数值模拟，得到了Re＝2800的充分发展不可压槽道湍流，并对该湍流场进行了统计分析。包括脉动速度偏斜因子在内的各阶统计量与实验结果及Kim等人的计算结果吻合十分理想，说明本方法是行之有效的。3．进行了三维充分发展的可压缩槽道湍流的直接数值模拟。得到了 Re=3300，Ma=0.8的充分发展可压槽道湍流的数据库。流场的统计特征（如等效平均速度分布，“半局部”尺度无量纲化的脉动速度均方根）和他人的数值计算结果吻合。得到了可压槽道湍流的各阶统计量，其中脉动速度的偏斜因子和平坦因子等高阶统计量尚未见其他文献报道。同时还分析了压缩性效应对壁湍流影响的机理，指出近壁处的压力－膨胀项将部分湍流脉动的动能转换成内能，使得可压湍流近壁速度条带结构更加平整。4．模拟了二维不可压槽道流动的饱和态（所谓“二维湍流”），分析了“二维槽道湍流”的非线性行为特征。分析了流场中的上抛－下扫和间歇现象，研究了“二维湍流”与三维湍流的区别。指出“二维湍流”反映了三维湍流的部分特征，同时指出了展向扰动对于湍流核心区发展的重要性。5．首次对可压缩槽道湍流及“二维槽道湍流”标度律进行了分析，得出了以下结论：a）槽道湍流中，在槽道中心线附近较宽的区域，存在标度律。b）该区域流场存在扩展自相似性（ESS）。c）在Mach数不是很高时，压缩性对标度指数影响不大。本文结果同SL标度律的理论值吻合较好，有效支持了该理论。对“二维槽道湍流”也有相似的结论，但与三维湍流不同的是，“二维槽道湍流”存在标度律的区域更宽，近壁处的标度指数比中心处有所升高。

Topology of toroidal helical fields in non-circular cross-sectional tokamaks

The ordinary differential magnetic field line equations are solved numerically; the tokamak magnetic structure is studied on Hefei Tokamak-7 Upgrade (HT-7U) when the equilibrium field with a monotonic q-profile is perturbed by a helical magnetic field. We find that a single mode (m, n) helical perturbation can cause the formation of islands on rational surfaces with q = m/n and q = (m +/- 1, +/- 2, +/- 3,...)/n due to the toroidicity and plasma shape (i.e. elongation and triangularity), while there are many undestroyed magnetic surfaces called Kolmogorov-Arnold-Moser (KAM) barriers on irrational surfaces. The islands on the same rational surface do not have the same size. When the ratio between the perturbing magnetic field B-r(r) and the toroidal magnetic field amplitude B(phi)0 is large enough, the magnetic island chains on different rational surfaces will overlap and chaotic orbits appear in the overlapping area, and the magnetic field becomes stochastic. It is remarkable that the stochastic layer appears first in the plasma edge region.

面向无线射频识别应用的复合订阅匹配

传统的发布/订阅系统不支持具有丰富时序关系的复合事件,不能很好地满足无线射频识别应用的需要。针对无线射频识别应用的特点,设计了面向无线射频识别应用的复合订阅语言,研究了订阅语言的组成和使用,基于时间事件、复合订阅的匹配图表示、事件检测的消费语义,给出了相应的复合订阅匹配算法。该算法已在OncePubSub系统中实现,并与Siena改造版J、ESS改造版进行了性能比较。