轴向磁场中电弧螺旋不稳定性的简化数值分析

在静电近似下，通过数量级的分析，简化了能量方程。用数值方法研究了外加轴向磁场，电流分布、弧电流等对电弧曙旋不稳定性的影响。数值结果表明，在具有外加磁场情况下，当弧柱内均匀电流分布时，正向磁场能够使短波长的扰动稳定，反向磁场能够使长波长的扰动稳定。当弧柱内电流为抛物线形分布时，外加正向磁场对于短波情况下的稳定性影响很小，但对于长波情况，稳定性增强，外加反向磁场，同样能稳定中等和长波情况下的电弧。     

The Implementation of Magnetic Islands in Gyrokinetic Toroidal Code

The implementation of magnetic islands in gyrokinetic simulation has been verified in the gyrokinetic toroidal code (GTC). The ion and electron density profiles become partially flattened inside the islands. The density profile at the low fleld side is less flattened than that at the high field side due to toroidally trapped particles in the low field side, which do not move along the perturbed magnetic field lines. When the fraction of trapped particles decreases, the density profile at the low field becomes more flattened.     

Study on Electron Heat Diffusion Across Stochastic Magnetic Field Affected by Magnitude of Perturbed Magnetic Field in Tokomak

Electron heat diffusion across stochastic magnetic fields is studied numerically in order to find out how the magnitude of perturbed magnetic field affect the enhanced heat conductivity and its radial profile in tokomak plasma physics. For these purposes, non-local stochastic magnetic fields are chosen as research objects in our simulation work. From our numerical results, we can find that the effects of the perturbed magnetic field magnitude are dominated parameter on the enhance electron heat transport conductivity wherever the magnetic field is single island or full stochastic field. Also, a theoretical analysis is provided and compared with numerical results.     

Directional power absorption in helicon plasma sources excited by a half-helix antenna

This paper deals with the investigation of the power absorption in helicon plasma excited through a half-helix antenna driven at 13.56 MHz. The simulations were carried out by means of a code,HELIC. They were carried out by taking into account different inhomogeneous radial density profiles and for a wide range of plasma densities, from 10~(11) cm~(-3) to 10~(13) cm~(-3). The magnetic field was 200, 400, 600 and 1000 G. A three-parameter function was used for generating various density profiles with different volume gradients, edge gradients and density widths. The density profile had a large effect on the efficient Trivelpiece–Gould(TG) and helicon mode excitation and antenna coupling to the plasma. The fraction of power deposition via the TG mode was extremely dependent on the plasma density near the plasma boundary. Interestingly, the obtained efficient parallel helicon wavelength was close to the anticipated value for Gaussian radial density profile.Power deposition was considerably asymmetric when the n/B_0 ratio was more than a specific value for a determined density width. The longitudinal power absorption was symmetric at approximately n_0 =10~(11) cm~(-3), irrespective of the magnetic field supposed. The asymmetry became more pronounced when the plasma density was 10~(12) cm~(-3). The ratio of density width to the magnetic field was an important parameter in the power coupling. At high magnetic fields, the maximum of the power absorption was reached at higher plasma density widths. There was at least one combination of the plasma density, magnetic field and density width for which the RF power deposition at both side of the tube reached its maximum value.     

Development of the W-band density profile and fluctuation reflectometer on EAST

A X-mode polarized W-band reflectometer for plasma density profile and fluctuation measurement is designed and installed on EAST. In measuring the density profile, a voltage controlled oscillator (VCO) is used as the source, allowing a high temporal resolution measurement. The density profile in a plasma with high magnetic field (3.0 T) has been measured by combination of V- and W-band reflectometers. For fluctuation measurements, a frequency synthesizer is used instead of the VCO as a microwave source. The core density fluctuations during sawtooth activity are measured and analyzed.     

On the heating mechanism of electron cyclotron resonance thruster immerged in a non-uniform magnetic field

To study the heating mechanism of electron cyclotron resonance thruster(ECRT) immersed in a non-uniform magnetic field, experiments and simulations are performed based on an electron cyclotron resonance plasma source at ASIPP. It is found that the first harmonic of electron cyclotron resonance is essential for plasma ignition at high magnetic field(0.0875 T), while the plasma can sustain without the first and second harmonics of electron cyclotron resonance at low magnetic field(till 0.0170 T). Evidence of radial hollow density profile indicates that upper hybrid resonance, which has strong edge heating effect, is the heating mechanism of low-field ECRT. The heating mode transition from electron cyclotron resonance to upper hybrid resonance is also revealed. Interestingly, the evolutions of electron temperature and electron density with input power experience a ‘delayed' jump, which may be correlated with the different power levels required for cyclotron and ionization. Moreover, when the field strength decreased, the variation of electron density behaves in an opposite trend with that of electron temperature,implying a possible competition of power deposition between them. The present work is of great interest for understanding the plasma discharge in ECRT especially immersed in a non-uniform magnetic field, and designing efficient ECRT using low magnetic field for economic space applications.     

Axial profiles of argon helicon plasma by optical emission spectroscope and Langmuir probe

We present the axial profiles of argon helicon plasma measured by a local optical emission spectroscope (OES) and Langmuir RF-compensated probe. The results show that the emission intensity of the argon atom lines (750 nm, 811 nm) is proportional to the plasma density determined by the Langmuir probe. The axial profile of helicon plasma depends on the discharge mode which changes with the RF power. Excited by helical antenna, the axial distribution of plasma density is similar to that of the external magnetic field in the capacitive coupled mode (E-mode). As the discharge mode changes into the inductively coupled mode (H-mode), the axial distribution of plasma density in the downstream can still be similar to that of the external magnetic field, but becomes more uniform in the upstream. When the discharge entered wave coupled mode (W-mode), the plasma becomes nearly uniform along the axis, showing a completely different profile from the magnetic field. The W-mode is expected to be a mixed pattern of helicon (H) and Trivelpiece-Gould (TG) waves.     

Simulation of EAST off-axis neutral beam heating and current drive

The capability of off-axis neutral beam heating and current drive has been investigated with NUBEAM for Experimental Advanced Superconducting Tokamak (EAST). Three different approaches to realize off-axis Neutral Beam Injection (NBI) have been studied. Simulation results for on- and off-axis NBI are reported. The effects of the alignment of NBI relative to the magnetic field pitch on off-axis neutral beam heating and current drive are observed and discussed qualitatively. By comparing the numerical results, a most favorable off-axis NBI configuration is recommended. The capability to control sawtooth is also investigated by comparing locations of the q = 1 rational surface and the peak of the fast ion density profile.     

Three-Dimensional Numerical Calculations on Liquid-Metal Magneto-hydrodynamic Flow through Circular Pipe in Magnetic-Field Inlet-Region

Three-dimensional numerical calculations have been performed on liquid-metal magnetohydrodynamic (MHD) flows through a circular pipe in the inlet region of the applied magnetic field, including a sufficient calculation region upstream in the magnetic field section. The continuity equation, the momentum equation including the Lorentz force term, and the induction equation derived from basic equations in the electromagnetism have been solved numerically. Along the flow axis (i.e., the channel axis), the pressure decreases slightly as a normal non-MHD flow, increases once, thereafter, decreases sharply, and finally decreases as a fully-developed MHD flow. The sharp decrease in the pressure, resulting in a large pressure drop in the inlet region is due to the increase in the induced electric current in this region compared with that in the fully-developed region. The velocity distribution changes from a parabolic profile of a laminar non-MHD flow to a profile with peaks near the walls parallel to the magnetic field, and finally to a flat profile of a fully-developed MHD flow.     

Magnetic field topology modeling under resonant magnetic perturbations on EAST

In order to understand the mechanism by which the resonant magnetic perturbation (RMP) mitigates or suppresses the edge-localized mode (ELM), the topological study of the edge magnetic field in ELM mitigation or suppression phase is a critical issue. To model the three-dimensional magnetic field topology superposed RMP on Experimental Advanced Superconducting Tokamak, a numerical model using the field line tracing method for both vacuum and ideal plasma response approximations is proposed. Using the numerical model, the topological change and the penetration depth of the stochastic field lines in the edge magnetic field are studied in an RMP experiment. Comparing profiles of minimum ρ on edge stochastic field lines and the particle flux pattern, the ideal plasma response changes the field line penetration depth while remaining similar profile relative to vacuum approximation. To mitigate and suppress ELM strongly, the deep penetration of RMP fields and topological changes of the edge magnetic field is a key from our modeling.     

Laminar pipe flow at the entrance into transverse magnetic field

High-resolution numerical simulations are conducted to analyze transformation of a liquid metal flow in a pipe at the entrance into a transverse magnetic field. The case of laminar flow, perfectly insulating pipe walls, and Hartmann number up to 200 is considered. The simulations reveal detailed structure of velocity and electric current fields and distribution of forces with particular attention given to the flow with an M-shaped velocity profile. They also establish criteria for accurate computations of laminar magnetohydrodynamic flows in strong non-uniform magnetic fields.     

强磁场作用下液态锂在矩形通道中的流动研究

对高速液态金属锂在强磁场作用下沿着非扩散矩形型通道的流动进行了分析研究 ,通过解析分析 ,推导得到了控制电场和流场的由泊松 (Poisson)方程和亥姆霍兹 (Helmholtz)方程组成的方程组 ,并且编制了求解该方程组的程序 PHsolver;根据边界的壁面函数处理 ,用PHsolver求解得到了强磁场作用下在非扩散型矩形通道中电流密度场的分布和流场的分布 ,其中流速分布呈现为马鞍型分布。     

Transmission characteristics of terahertz Bessel vortex beams through a multi-layered anisotropic magnetized plasma slab

The transmission of terahertz (THz) Bessel vortex beams through a multi-layered anisotropic magnetized plasma slab is investigated by using a hybrid method of cylindrical vector wave functions (CVWFs) and Fourier transform. On the basis of the electromagnetic boundary conditions on each interface, a cascade form of expansion coefficients of the reflected and transmitted fields is obtained. Taking a double Gaussian distribution of the plasma density as an example, the influences of the applied magnetic field, the incident angle and polarization mode of the incident beams on the magnitude, OAM mode and polarization of the transmitted beams are analyzed in detail. The results indicate that the applied magnetic field has a major effect upon the polarization state of the transmitted fields but not upon the transmitted OAM spectrum. The incident angle has a powerful influence upon both the amplitude profile and the OAM spectrum of the transmitted beam. Furthermore, for multiple coaxial vortex beams, an increase of the maximum value of the plasma density causes more remarkable distortion of both the profile and OAM spectrum of the transmitted beam. This research makes a stable foundation for the THz OAM multiplexing/demultiplexing technology in a magnetized plasma environment.     

Development of atomic beam probe for tokamaks

The concept and development of a new detection method for light alkali ions stemming from diagnostic beams installed on medium size tokamak is described. The method allows us the simultaneous measurement of plasma density fluctuations and fast variations in poloidal magnetic field, therefore one can infer the fast changes in edge plasma current. The concept has been worked out and the whole design process has been done at Wigner RCP. The test detector with appropriate mechanics and electronics is already installed on COMPASS tokamak. General ion trajectory calculation code (ABPIons) has also been developed. Detailed calculations show the possibility of reconstruction of edge plasma current density profile changes with high temporal resolution, and the possibility of density profile reconstruction with better spatial resolution compared to standard Li-BES measurement, this is important for pedestal studies.     

Self-field calculation of CICC with fast direct Biot–Savart integration

ITER magnetic device (Tokamak) requires a strong magnetic field produced by charged cable conductors and external sources to arrive at stable and reliable magnetic confinement performance. Before manufacturing and assembling conductors, preliminary analysis of self-field induction is helpful for reducing the cost of varying-parameter experiments. Spatial helix shape of numerous strand elements and multi-level twist of the finalized cable, known as CICC type, make it unpractical to direct use finite-element methods and other numerical procedures for self-field calculation. An algorithm FDBS (fast direct Biot–Savart integration) is proposed to surmount this difficulty, which improves the traditional method (DBS, direct implementing Biot–Savart law for all strand sources) in terms of computational effort. As such the complexity reduces to O(N) from the original O(N²) and speed enhancement is achieved in the parallel computation environment. FDBS calculates out a detailed self-field profile for the uncompressed ITER TF conductors carrying uniform current at each cabling level; the layered self-field distribution becomes more indistinct for higher level subcable.     

SSRF波荡器的渡越辐射特性研究

通过解析分析和数值模拟计算，对上海同步辐射装置(SSRF)波荡器U9.0在红外及远红外波段(波长范围8200μm)的渡越辐射谱通量密度角分布、亮度等特性进行了研究，并探讨了电子束发射度、能量散度及波荡器磁场误差对渡越辐射谱亮度的影响。结果表明，SSRF U9.0的红外及远红外渡越辐射通量密度最高可达351013Photons/(smrad2BW)，收集角φ=0.23mrad，且通量密度随电子上的发射度、能散和波荡器磁场误差等的变化很小。     

An RF Technique for Plasma Number Density Measurement

It is well-known that a cylindrical plasma column in a transverse rf electric field is subject to a dipole resonance. The frequency at which this occurs can be used to determine the average electron density across the plasma section. This paper describes the experimental application of the technique to plasmas with and without a confining axial static magnetic field, and gives theoretical curves to assist in rapid reduction of the data. Higher order multipole resonances can also be excited. The feasibility of applying these to determination of the number density profile is discussed.     

Feasibility of applying the lower cut-off frequency for the density radial coverage extension in EAST reflectometry measurement

The extraordinary mode (X-mode) lower cut-off frequency is proposed for use in the reflectometry diagnostic on ITER for the electron density profile measurement, which is a trade-off between extreme plasma parameters and the accessible probing frequency. In contemporary experiments, the lower cut-off frequency can be identified at the probing frequency below the electron cyclotron frequency (fce) under certain plasma conditions. We provide here, for the first time, the experimental validation of the use of the lower cut-off frequency for the density profiles via the reflectometry measurement on EAST. The corresponding group delay of the lower cut-off frequency evolves continuously with the upper one, revealing a reasonable radial coverage extension of reflectometry measurement toward the plasma core. It is concluded that the lower cut-off frequency can be used as a supplement to the upper one in the density profile inversion process, which is of particular interest in the high magnetic field and/or density discharge to extend the radial coverage of reflectometry measurement.     

Semi-Analytical Estimation of Gamma-Ray Intensity from a Volume Source

The effect of a magnetic field on laminar natural convection of liquid metal was studied experimentally using NaK as conducting fluid. The magnetic field was imposed horizontally and parallel to a uniformly heated vertical plate, to act perpendicularly across the convective flow. In a low magnetic field, the temperature profile across the layer of flowing fluid acquired an η-shaped profile characterized by a valley close to the wall and a peak further away, which had the effect of raising heat transfer rate above that obtained in the absence of magnetic field. When the magnetic field was intensified, its braking effect on the flow approached the temperature profile to the case of pure heat conduction through solid, to result in dwindling heat transfer rate.     

Effects of the consistency of the fringing magnetic field on direct numerical simulations of liquid-metal flow

We investigate the effects of the consistency of strong fringing decreasing magnetic fields on numerical simulations of classical experimental data. Studies about fringing magnetic fields have attracted the attention of the fusion community in relation to the design of the liquid-metal flow blankets for fusion nuclear reactors. One-dimensional fitting functions neglecting magnetic field consistency have been adopted in previous numerical studies. Thanks to complete three-dimensional numerical simulations, the effect of the physical consistency of the magnetic field on fluid flow can now be assessed. We present a technique for generating discretely consistent magnetic fields based on classical one-dimensional fittings. With this method, key magnetic field features, such as the bending of the magnetic lines, are accurately reproduced and, therefore, the validity of the technique is established. Consistent and inconsistent magnetic fields have been tested under very strong decreasing magnetic fields with insulating and conducting walls using direct numerical simulations. The results show a moderate, but systematic, improvement of the predictions with respect to the experiments. As an example, the repeated under-prediction of the peak transverse pressure gradient, observed in the results of asymptotic methods and of direct numerical simulations, is explained by the historically neglected consistency of the fringing magnetic field.