Eigenmode analysis of the sheared-flow Z-pinch
Autor: | J. R. Angus, J. J. Van De Wetering, Mikhail Dorf, V. I. Geyko |
---|---|
Rok vydání: | 2020 |
Předmět: | |
Zdroj: | Physics of Plasmas. 27:122108 |
ISSN: | 1089-7674 1070-664X |
DOI: | 10.1063/5.0029716 |
Popis: | Experiments have demonstrated that a Z-pinch can persist for thousands of times longer than the growth time of global magnetohydrodynamic (MHD) instabilities such as the m = 0 sausage and m = 1 kink modes. These modes have growth times on the order of t a = a / v i, where vi is the ion thermal speed and a is the pinch radius. Axial flows with d u z / d r ≲ v i / a have been measured during the stable period, and the commonly accepted theory is that this amount of shear is sufficient to stabilize these modes as predicted by numerical studies using the ideal MHD equations. However, these studies only consider specific equilibrium profiles that typically have a modest magnitude for the logarithmic pressure gradient, q P ≡ d ln P / d ln r, and may not represent experimental conditions. Linear stability of the sheared-flow Z-pinch is studied here via a direct eigen-decomposition of the matrix operator obtained from the linear ideal MHD equations. Several equilibrium profiles with a large variation of qP are examined. Considering a practical range of k, 1 / 3 ≲ ka ≲ 10, it is shown that the shear required to stabilize m = 0 modes can be expressed as d u z / d r ≥ C γ 0 / ( k a ) α. Here, γ 0 = γ 0 ( k a ) is the profile-specific growth rate in the absence of shear, which scales approximately with | q P |. Both C and α are profile-specific constants, but C is order unity and α ≈ 1. It is further demonstrated that even a large value of shear, d u z / d r = 3 v i / a, is not sufficient to provide linear stabilization of the m = 1 kink mode for all profiles considered. This result is in contrast to the currently accepted theory predicting stabilization at much lower shear, d u z / d r = 0.1 v i / a, and suggests that the experimentally observed stability cannot be explained within the linear ideal-MHD model. |
Databáze: | OpenAIRE |
Externí odkaz: |