-
To understand the effect of given radial electric field on ion-temperature gradient driven mode (ITG) stability in tokamak plasmas, we derive the eigenmode equation for ITG including the poloidal rotation and density modulation associated with radial electric field by using nonlinear gyrokinetic theory. The equation is solved for the eigenfrequency, growth rate and parallel mode structure of ITG both in short- and long-wavelength limit with energetic-particle-induced geodesic acoustic mode (EGAM) as a specific form. The eigenmode equation is not only solved analytically, but also solved numerically to validate the analytic solutions. It is found that, the radial electric field induced poloidal rotation can significantly stabilize ITG, while the density perturbation of the radial electric field may slightly distort the ITG parallel mode structure, but has little effect on ITG stability. The result is consistent with the common picture of turbulence suppression by poloidal shear flow. The general model is also applicable to the investigation of the indirect interaction of ITG and energetic particle driven Alfvén instabilities via zonal structures generation, by introducing poloidal rotation and density modulation associated with zonal structures spontaneously excited by Alfvén instabilities. The indirect channel is supplement to the direct interaction of microturbulences and energetic particle driven Alfvén instabilities.
-
Keywords:
- ion-temperature gradient driven mode /
- radial electric field /
- energetic-particle induced geodesic acoustic mode /
- nonlinear gyrokinetic theory
[1] Horton W 1999 Rev. Mod. Phys. 71 735Google Scholar
[2] Lin Z H, Hahm T S, Lee W W, Tang W, White R 1998 Science 281 1835Google Scholar
[3] Hahm T S, Beer M, Lin Z H, Hammett G, Lee W W, Tang W 1999 Phys. Plasmas 6 922Google Scholar
[4] Chen L, Lin Z H, White R 2000 Phys. Plasmas 7 3129Google Scholar
[5] Zonca F, Chen L 2008 Europhys. Lett. 83 35001Google Scholar
[6] Nazikian R, Fu G Y, Austin M, Berk H, Budny R, Gorelenkov N, Heidbrink W, Holcomb C, Kramer G, McKee G, Makowski M, Solomon W, Shafer M, Strait E, Van Zeeland M 2008 Phys. Rev. Lett. 101 185001Google Scholar
[7] Fu G Y 2008 Phys. Rev. Lett. 101 185002Google Scholar
[8] Qiu Z Y, Zonca F, Chen L 2010 Plasma Phys. Control. Fusion 52 095003Google Scholar
[9] Berk H, Zhou T 2010 Nuclear Fusion 50 035007Google Scholar
[10] Qiu Z Y, Chen L, Zonca F, Chen W 2018 Phys. Rev. Lett. 120 135001Google Scholar
[11] Zarzoso D, Sarazin Y, Garbet X, Dumont R, Strugarek A, Abiteboul J, Cartier-Michaud T, Dif-Pradalier G, Ghendrih P, Grandgirard V, Latu G, Passeron C, Thomine O 2013 Phys. Rev. Lett. 110 125002Google Scholar
[12] Dimits A, Bateman G, Beer M, Cohen B, Dorland W, Hammett G, Kim C, Kinsey J, Kotschenreuther M, Kritz A, Lao L L, Mandrekas J, Nevins W, Parker S, Redd A, Shumaker D, Sydora R, Weiland J 2000 Phys. Plasmas 7 969Google Scholar
[13] Qiu Z Y, Chen L, Zonca F 2014 Phys. Plasmas 21 022304Google Scholar
[14] Qiu Z Y, Chen L, Zonca F 2015 Phys. Plasmas 22 042512Google Scholar
[15] Tardini G, Hobirk J, Igochine V, Maggi C, Martin P, McCune D, Peeters A, Sips A, Stäbler A, Stober J, the ASDEX Upgrade Team 2007 Nuclear Fusion 47 280Google Scholar
[16] Mantica P, Strintzi D, Tala T, Giroud C, Johnson T, Leggate H, Lerche E, Loarer T, Peeters A, Salmi A, Sharapov S, Van Eester D, de Vries P, Zabeo L, Zastrow K 2009 Phys. Rev. Lett. 102 175002Google Scholar
[17] Romanelli M, Zocco A, Crisanti F, JET-EFDA Contributors 2010 Plasma Phys. Control. Fusion 52 045007Google Scholar
[18] Garcia J, Challis C, Citrin J, Doerk H, Giruzzi G, Görler T, Jenko F, Maget P. JET Contributors 2015 Nuclear Fusion 55 053007Google Scholar
[19] Citrin J, Jenko F, Mantica P, Told D, Bourdelle C, Garcia J, Haverkort J, Hogeweij G, Johnson T, Pueschel M 2013 Phys. Rev. Lett. 11 1Google Scholar
[20] Chankin A, McCracken G 1993 Nuclear Fusion 33 1459Google Scholar
[21] Cheng C Z, Chen L, Chance M 1985 Ann. Phys. 161 21Google Scholar
[22] Chen L 1994 Phys. Plasmas 1 1519Google Scholar
[23] Qiu Z Y, Chen L, Zonca F 2016 Nuclear Fusion 56 106013Google Scholar
[24] Qiu Z Y, Chen L, Zonca F 2017 Nuclear Fusion 57 056017Google Scholar
[25] Chen L, Qiu Z Y, Zonca F 2022 Nuclear Fusion 62 094001Google Scholar
[26] Chen L, Qiu Z Y, Zonca F 2014 Europhys. Lett. 107 15003Google Scholar
[27] Romanelli F, Zonca F 1993 Phys. Fluids B 5 4081Google Scholar
[28] Chen L, Briguglio S, Romanelli F 1991 Phys. Fluids B 3 611Google Scholar
[29] Chen N F, Hu H Y, Zhang X Y, Wei S Z, Qiu Z Y 2021 Phys. Plasmas 28 042505Google Scholar
-
图 1 短波长极限下,
$ e{{\text{δ}}}\phi_{{\rm E}}/T_{{\rm i}}=0 $ (a)和$ 0.1 $ (b)时, ITG本征值的实部($ \varOmega_{\rm r} $ )和虚部($ \varOmega_{\rm i} $ )的分布. 在两种情况下, 最不稳定的都是l = 0的基态Figure 1. Distribution of the real (
$ \varOmega_{\rm r} $ ) and imaginary ($ \varOmega_{\rm i} $ ) parts of eigenvalues of ITG when$ e{{\text{δ}}}\phi_ {{\rm E}}/T_{{\rm i}}=0 $ (a) and$ 0.1 $ (b) in the short-wavelength limit. In both cases, the ground state with l = 0 is the most unstable eigenstate.图 2 短波长极限下,
$ e{{\text{δ}}}\phi_{{\rm E}}/T_{{\rm i}}=0.1 $ 和0时, ITG最不稳定的基态的平行模结构. 其中, 蓝色实线表示径向电场为零时ITG的模结构, 红色虚线表示有限径向电场时ITG的模结构Figure 2. Mode structure of the most unstable mode of ITG. The blue solid and red dashed lines represent the cases with
$ e{{\text{δ}}}\phi_{{\rm E}}/T_{{\rm i}}=0 $ and$ e{{\text{δ}}}\phi_{{\rm E}}/T_{{\rm i}}=0.1 $ , respectively.图 3 短波长极限下, ITG的增长率 (a)和实频(b)与EGAM幅度
$ e{{\text{δ}}}\phi_{{\rm E}}/T_{{\rm i}} $ 的关系. 蓝色圆点表示直接数值求解本征值方程(9)得到的数值解, 红色叉号表示求解基态色散关系(14)得到的理论解. 实频和增长率都是用$ C_{\rm s}/L_{T_{\rm i}} $ 进行归一化, 其中,$ C_{\rm s}^2=2 T_{\rm e}/m_{\rm i} $ 表示声速; 这两幅图所用的参数为$ \varepsilon_ {T_{\rm i}}= $ $ 0.2 $ ,$ b=1 $ Figure 3. Dependence of the growth rate (a) and real frequency (b) of ITG on the amplitude of EGAM
$ e{{\text{δ}}}\phi_{{\rm E}}/T_{{\rm i}} $ in short-wavelength limit. The blue circles represent the numerical value obtained by directly solving the eigenmode equation (9); while the red crosses represent the theoretical value obtained by solving the dispersion relation (14). The real frequency and growth rate of ITG are normalized to$ C_{\rm s}/L_{T_{\rm i}} $ , with$ C_{\rm s}^2=2 T_{\rm e}/m_{\rm i} $ representing sound velocity. The parameters used here are$ \varepsilon_ {T_{\rm i}}=0.2 $ and$ b=1 $ .图 4 短波长极限下, 极向旋转和密度调制共同作用和分别作用时 ITG 增长率对 EGAM 幅度的依赖关系. 图中绿色实线表示只有密度扰动的情形, 红色虚线表示只有电势扰动的情形, 蓝色实线表示两种扰动都存在的情形
Figure 4. Dependence of the growth rate of ITG in presence of poloidal rotation and/or density modulation in short-wavelength limit. The green solid and red dashed lines represent the cases with only density modulation and poloidal rotation, respectively; while the blue line represents the case with both effects.
图 5 长波长极限下,
$ e{{\text{δ}}}\phi_{{\rm E}}/T_{{\rm i}} = 0 $ (a)和$ 0.1 $ (b)时, ITG的本征值分布. 在两种情况下, 最不稳定的都是l = 0的基态Figure 5. Distribution of eigenvalues of ITG when
$ e{{\text{δ}}}\phi_ {{\rm E}}/T_{{\rm i}} = 0 $ (a) and$ 0.1 $ (b) in the long-wavelength limit. In both cases, the ground state with l = 0 is the most unstable eigenstate.图 6 长波长极限下,
$ e{{\text{δ}}}\phi_{{\rm E}}/T_{{\rm i}}=0.1 $ 和0时, ITG最不稳定的基态的平行模结构. 其中, 蓝色实线表示径向电场为零时ITG的模结构, 红色虚线表示有限径向电场时ITG的模结构Figure 6. Mode structure of the most unstable mode of ITG. The blue solid and red dashed lines represent the cases with
$ e{{\text{δ}}}\phi_{{\rm E}}/T_{{\rm i}}=0 $ and$ e{{\text{δ}}}\phi_{{\rm E}}/T_{{\rm i}}=0.1 $ , respectively.图 7 长波长极限下, ITG的增长率 (a)和实频(b)与EGAM幅度
$ e{{\text{δ}}}\phi_{{\rm E}}/T_{{\rm i}} $ 的关系. 蓝色圆点表示直接数值求解本征值方程(9) 得到的数值解, 红色叉号表示求解基态色散关系(19)得到的理论解. 这两幅图所用的参数为$ \varepsilon_ {T_{\rm i}}=0.2 $ , b = 0.01Figure 7. Dependence of the growth rate (a) and real frequency (b) of ITG on the amplitude of EGAM
$ e{{\text{δ}}}\phi_{{\rm E}}/T_{{\rm i}} $ in long-wavelength limit. The blue circles represent the numerical value obtained by directly solving the eigenmode equation (9); while the red crosses represent the theoretical value obtained by solving the dispersion relation (19). The parameters used here are$ \varepsilon_ {T_{\rm i}}=0.2 $ and b = 0.01.图 8 长波长极限下, 极向旋转和密度调制共同作用和分别作用时 ITG 增长率对 EGAM 幅度的依赖关系. 图中绿色实线表示只有密度扰动的情形, 红色虚线表示只有电势扰动的情形, 蓝色实线表示两种扰动都存在的情形
Figure 8. Dependence of the growth rate of ITG in presence of poloidal rotation and/or density modulation in long-wavelength limit. The green solid and red dashed lines represent the cases with only density modulation and poloidal rotation, respectively; while the blue line represents the case with both effects.
-
[1] Horton W 1999 Rev. Mod. Phys. 71 735Google Scholar
[2] Lin Z H, Hahm T S, Lee W W, Tang W, White R 1998 Science 281 1835Google Scholar
[3] Hahm T S, Beer M, Lin Z H, Hammett G, Lee W W, Tang W 1999 Phys. Plasmas 6 922Google Scholar
[4] Chen L, Lin Z H, White R 2000 Phys. Plasmas 7 3129Google Scholar
[5] Zonca F, Chen L 2008 Europhys. Lett. 83 35001Google Scholar
[6] Nazikian R, Fu G Y, Austin M, Berk H, Budny R, Gorelenkov N, Heidbrink W, Holcomb C, Kramer G, McKee G, Makowski M, Solomon W, Shafer M, Strait E, Van Zeeland M 2008 Phys. Rev. Lett. 101 185001Google Scholar
[7] Fu G Y 2008 Phys. Rev. Lett. 101 185002Google Scholar
[8] Qiu Z Y, Zonca F, Chen L 2010 Plasma Phys. Control. Fusion 52 095003Google Scholar
[9] Berk H, Zhou T 2010 Nuclear Fusion 50 035007Google Scholar
[10] Qiu Z Y, Chen L, Zonca F, Chen W 2018 Phys. Rev. Lett. 120 135001Google Scholar
[11] Zarzoso D, Sarazin Y, Garbet X, Dumont R, Strugarek A, Abiteboul J, Cartier-Michaud T, Dif-Pradalier G, Ghendrih P, Grandgirard V, Latu G, Passeron C, Thomine O 2013 Phys. Rev. Lett. 110 125002Google Scholar
[12] Dimits A, Bateman G, Beer M, Cohen B, Dorland W, Hammett G, Kim C, Kinsey J, Kotschenreuther M, Kritz A, Lao L L, Mandrekas J, Nevins W, Parker S, Redd A, Shumaker D, Sydora R, Weiland J 2000 Phys. Plasmas 7 969Google Scholar
[13] Qiu Z Y, Chen L, Zonca F 2014 Phys. Plasmas 21 022304Google Scholar
[14] Qiu Z Y, Chen L, Zonca F 2015 Phys. Plasmas 22 042512Google Scholar
[15] Tardini G, Hobirk J, Igochine V, Maggi C, Martin P, McCune D, Peeters A, Sips A, Stäbler A, Stober J, the ASDEX Upgrade Team 2007 Nuclear Fusion 47 280Google Scholar
[16] Mantica P, Strintzi D, Tala T, Giroud C, Johnson T, Leggate H, Lerche E, Loarer T, Peeters A, Salmi A, Sharapov S, Van Eester D, de Vries P, Zabeo L, Zastrow K 2009 Phys. Rev. Lett. 102 175002Google Scholar
[17] Romanelli M, Zocco A, Crisanti F, JET-EFDA Contributors 2010 Plasma Phys. Control. Fusion 52 045007Google Scholar
[18] Garcia J, Challis C, Citrin J, Doerk H, Giruzzi G, Görler T, Jenko F, Maget P. JET Contributors 2015 Nuclear Fusion 55 053007Google Scholar
[19] Citrin J, Jenko F, Mantica P, Told D, Bourdelle C, Garcia J, Haverkort J, Hogeweij G, Johnson T, Pueschel M 2013 Phys. Rev. Lett. 11 1Google Scholar
[20] Chankin A, McCracken G 1993 Nuclear Fusion 33 1459Google Scholar
[21] Cheng C Z, Chen L, Chance M 1985 Ann. Phys. 161 21Google Scholar
[22] Chen L 1994 Phys. Plasmas 1 1519Google Scholar
[23] Qiu Z Y, Chen L, Zonca F 2016 Nuclear Fusion 56 106013Google Scholar
[24] Qiu Z Y, Chen L, Zonca F 2017 Nuclear Fusion 57 056017Google Scholar
[25] Chen L, Qiu Z Y, Zonca F 2022 Nuclear Fusion 62 094001Google Scholar
[26] Chen L, Qiu Z Y, Zonca F 2014 Europhys. Lett. 107 15003Google Scholar
[27] Romanelli F, Zonca F 1993 Phys. Fluids B 5 4081Google Scholar
[28] Chen L, Briguglio S, Romanelli F 1991 Phys. Fluids B 3 611Google Scholar
[29] Chen N F, Hu H Y, Zhang X Y, Wei S Z, Qiu Z Y 2021 Phys. Plasmas 28 042505Google Scholar
Catalog
Metrics
- Abstract views: 2597
- PDF Downloads: 75
- Cited By: 0