اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدStructure and overstability of resistive modes with runaway electrons
154
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We investigate the effects of runaway electron current on the dispersion relation of resistive magnetohydrodynamic modes in tokamaks. We present a new theoretical model to derive the dispersion relation, which is based on the asymptotic analysis of the resistive layer structure of the modes. It is found that in addition to the conventional resistive layer, a new runaway current layer can emerge whose properties depend on the ratio of the Alfven velocity to the runaway electron convection speed. Due to the contribution from this layer, both the tearing mode and kink mode will have a real frequency in addition to a growth rate. The derived dispersion relation has been compared with numerical results using both a simplified eigenvalue calculation and a M3D-C1 linear simulation, and good agreement is found in both cases.
قيم البحث
اقرأ أيضاً
A new fluid model for runaway electron simulation based on fluid description is introduced and implemented in the magnetohydrodynamics code M3D-C1, which includes self-consistent interactions between plasma and runaway electrons. The model utilizes t
he method of characteristics to solve the continuity equation for the runaway electron density with large convection speed, and uses a modified Boris algorithm for pseudo particle pushing. The model was employed to simulate magnetohydrodynamics instabilities happening in a runaway electron final loss event in the DIII-D tokamak. Nonlinear simulation reveals that a large fraction of runaway electrons get lost to the wall when kink instabilities are excited and form stochastic field lines in the outer region of the plasma. Plasma current converts from runaway electron current to Ohmic current, and get pinched at the magnetic axis. Given the good agreement with experiment, the simulation model provides a reliable tool to study macroscopic plasma instabilities in existence of runaway electron current, and can be used to support future studies of runaway electron mitigation strategies in ITER.
Improved understanding of runaway-electron formation and decay processes are of prime interest for the safe operation of large tokamaks, and the dynamics of the runaway electrons during dynamical scenarios such as disruptions are of particular concer
n. In this paper, we present kinetic modelling of scenarios with time-dependent plasma parameters; in particular, we investigate hot-tail runaway generation during a rapid drop in plasma temperature. With the goal of studying runaway-electron generation with a self-consistent electric-field evolution, we also discuss the implementation of a collision operator that conserves momentum and energy and demonstrate its properties. An operator for avalanche runaway-electron generation, which takes the energy dependence of the scattering cross section and the runaway distribution into account, is investigated. We show that the simplified avalanche model of Rosenbluth & Putvinskii [Nucl. Fusion 1997 37 1355] can give inaccurate results for the avalanche growth rate (either lower or higher) for many parameters, especially when the average runaway energy is modest, such as during the initial phase of the avalanche multiplication. The developments presented pave the way for improved modelling of runaway-electron dynamics during disruptions or other dynamic events.
A kinetic formalism of parametric decay of a large amplitude lower hybrid pump wave into runaway electron mode and a uppersideband mode is investigated. The pump and the sideband exert a ponderomotive force on runaway electrons, driving the runaway m
ode. The density perturbation associated with the latter beats with the oscillatory velocity due to the pump to produce the sideband. The finite parallel velocity spread of the runaway electrons turns the parametric instability into a stimulated compton scattering process where growth rate scales as the square of the pump amplitude. The large phase velocity waves thus generated can potentially generate relativistic electrons.
In this research, formation of finite two-dimensional (2D) plasma crystals was numerically simulated. The structure was investigated for systems with various Debye length and it was found there is a transition from a complete hexagonal structure to a
structure with concentric rings on the outer edge and hexagonal lattice in the interior as the Debye length increases. The vertical as well as horizontal oscillation modes thermally excited in the 2D dust coulomb clusters were investigated. The horizontal mode spectra is shown to agree with published results while the vertical mode spectra obtained from numerical simulation and analytical method agree with one another. The frequency of the vertical modes is shown to have a maximum corresponding to the whole system acting as a solid plane. For low frequency modes, the largest amplitude particle motion is concentrated in a few inner rings with the outer rings remaining almost motionless in contrast to the horizontal modes for which the strongest motion of the particles is concentrated in the inner rings at high frequencies.
The stability of the $n=1$ resistive wall modes (RWMs) is investigated using the AEGIS code for the newly designed China Fusion Engineering Test Reactor (CFETR) 1GW steady-state operating (SSO) scenario. Here, $n$ is the toroidal mode number. Due to
the large fraction of bootstrap current contribution, the profile of safety factor q is deeply reversed in magnetic shear in the central core region and locally flattened within the edge pedestal. Consequently the pressure-driven infernal components develop in the corresponding q-flattened regions of both core and edge. However, the edge infernal components dominate the $n=1$ RWM structure and lead to lower $beta_N$ limits than the designed target $beta_N$ for the CFETR 1GW SSO scenario. The edge rotation is found the most critical to the stabilization due to the dominant influence of the edge infernal components, which should be maintained above $1.5%Omega_{A0}$ in magnitude in order for the rotation alone to fully suppress the $n=1$ RWM in the CFETR 1GW SSO scenario.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
