The linear and nonlinear simulations are carried out using the gyrokinetic code NLT for the electrostatic instabilities in the core region of a deuterium plasma based on the International Thermonuclear Experimental Re...The linear and nonlinear simulations are carried out using the gyrokinetic code NLT for the electrostatic instabilities in the core region of a deuterium plasma based on the International Thermonuclear Experimental Reactor(ITER)baseline scenario.The kinetic electron effects on the linear frequency and nonlinear transport are studied by adopting the adiabatic electron model and the fully drift-kinetic electron model in the NLT code,respectively.The linear simulations focus on the dependence of linear frequency on the plasma parameters,such as the ion and electron temperature gradientsκ_(Ti,e)≡R=L_(Ti,e),the density gradientκ_(n)≡R/L_(n)and the ion-electron temperature ratioτ=T_(e)=T_(i).Here,is the major radius,and T_(e)and T_(i)denote the electron and ion temperatures,respectively.L_(A)=-(δ_(r)lnA)^(-1)is the gradient scale length,with denoting the density,the ion and electron temperatures,respectively.In the kinetic electron model,the ion temperature gradient(ITG)instability and the trapped electron mode(TEM)dominate in the small and large k_(θ)region,respectively,wherek_(θ)is the poloidal wavenumber.The TEMdominant region becomes wider by increasing(decreasing)κ_(T_(e))(κ_(T_(i)))or by decreasingκ_(n).For the nominal parameters of the ITER baseline scenario,the maximum growth rate of dominant ITG instability in the kinetic electron model is about three times larger than that in the adiabatic electron model.The normalized linear frequency depends on the value ofτ,rather than the value of T_(e)or T_(i),in both the adiabatic and kinetic electron models.The nonlinear simulation results show that the ion heat diffusivity in the kinetic electron model is quite a lot larger than that in the adiabatic electron model,the radial structure is finer and the time oscillation is more rapid.In addition,the magnitude of the fluctuated potential at the saturated stage peaks in the ITGdominated region,and contributions from the TEM(dominating in the higher k_(θ)region)to the nonlinear transport can be neglected.In the adiabatic electron model,the zonal radial electric field is found to be mainly driven by the turbulent energy flux,and the contribution of turbulent poloidal Reynolds stress is quite small due to the toroidal shielding effect.However,in the kinetic electron model,the turbulent energy flux is not strong enough to drive the zonal radial electric field in the nonlinear saturated stage.The kinetic electron effects on the mechanism of the turbulence-driven zonal radial electric field should be further investigated.展开更多
A conservative scheme of kinetic electrons for gyrokinetic simulations in the presence of magnetic islands has been implemented and verified in the gyrokinetic toroidal code, where zonal and nonzonal components of all...A conservative scheme of kinetic electrons for gyrokinetic simulations in the presence of magnetic islands has been implemented and verified in the gyrokinetic toroidal code, where zonal and nonzonal components of all perturbed quantities are solved together. Using this new conservative scheme, linear simulation of kinetic ballooning mode has been successfully benchmarked with the electromagnetic hybrid model. Simulations of nonlinear interactions between magnetic islands and the ion temperature gradient(ITG) mode in a tokamak show that the islands rotate at the electron diamagnetic drift velocity. The linear ITG structure shifts from the island O-point toward the X-point due to the pressure flattening effect inside the islands, and the nonlinear ITG structure peaks along the magnetic island separatrix because of the increased pressure gradient there.展开更多
The turbulence characteristics of plasmas with internal transport barriers in the HL-2A tokamak are analyzed by means of linear gyrokinetic simulations. It is found that turbulence is dominated by the ion temperature ...The turbulence characteristics of plasmas with internal transport barriers in the HL-2A tokamak are analyzed by means of linear gyrokinetic simulations. It is found that turbulence is dominated by the ion temperature gradient(ITG)mode together with large-scale modes characterized by high-frequency electromagnetic fluctuation, which are destabilized by the steep ion temperature gradient in the weak magnetic shear regime. Comparison with solutions of analytical dispersion relations shows that their linear features match well with the beta-induced Alfvén eigenmode branch of the shear Alfvénic spectrum. It is further clarified that the large population of fast ions in these plasmas plays a stabilization role through the dilution mechanism in high-n ITG mode regimes.展开更多
Field-aligned coordinates have been implemented in the gyrokinetic semi-Lagrangian code NLT, Ye et al （2016 J. Comput. Phys. 316 180）, to improve the computational efficiency for the numerical simulations of tokamak...Field-aligned coordinates have been implemented in the gyrokinetic semi-Lagrangian code NLT, Ye et al （2016 J. Comput. Phys. 316 180）, to improve the computational efficiency for the numerical simulations of tokamak turbulence and transport. 4D B-spline interpolation in field- aligned coordinates is applied to solve the gyrokinetic Vlasov equation. A fast iterative algorithm is proposed for efficiently solving the quasi-neutrality equation. A pseudo transform method is used for the numerical integration of the gyro-average operator for perturbations with a high toroidal mode number. The new method is shown to result in an improved code performance for reaching a given accuracy. Some numerical tests are presented to illustrate the new methods.展开更多
Gyrokinetic theory is arguably the most important tool for numerical studies of transport physics in magnetized plasmas.However,exact local energy–momentum conservation laws for the electromagnetic gyrokinetic system...Gyrokinetic theory is arguably the most important tool for numerical studies of transport physics in magnetized plasmas.However,exact local energy–momentum conservation laws for the electromagnetic gyrokinetic system have not been found despite continuous effort.Without such local conservation laws,energy and momentum can be instantaneously transported across spacetime,which is unphysical and casts doubt on the validity of numerical simulations based on the gyrokinetic theory.The standard Noether procedure for deriving conservation laws from corresponding symmetries does not apply to gyrokinetic systems because the gyrocenters and electromagnetic field reside on different manifolds.To overcome this difficulty,we develop a high-order field theory on heterogeneous manifolds for classical particle-field systems and apply it to derive exact,local conservation laws,in particular the energy–momentum conservation laws,for the electromagnetic gyrokinetic system.A weak Euler–Lagrange(EL)equation is established to replace the standard EL equation for the particles.It is discovered that an induced weak EL current enters the local conservation laws,and it is the new physics captured by the high-order field theory on heterogeneous manifolds.A recently developed gauge-symmetrization method for high-order electromagnetic field theories using the electromagnetic displacement-potential tensor is applied to render the derived energy–momentum conservation laws electromagnetic gauge-invariant.展开更多
Global linear gyrokinetic simulations using realistic DIII-D tokamak geometry and plasma profiles find co-existence of unstable reversed shear Alfvén eigenmodes(RSAE)with low toroidal mode number n and electromag...Global linear gyrokinetic simulations using realistic DIII-D tokamak geometry and plasma profiles find co-existence of unstable reversed shear Alfvén eigenmodes(RSAE)with low toroidal mode number n and electromagnetic ion temperature gradient(ITG)instabilities with higher toroidal mode number n.For intermediate n?=?[10,12],RSAE and ITG co-exist and overlap weakly in the radial domain with similar growth rates but different real frequencies.Both RSAE and ITG growth rates decrease less than 5%when compressible magnetic perturbations are neglected in the simulations.The ITG growth rates increase less than 7%when fast ions are not included in the simulations.Finally,the effects of trapped electrons on the RSAE are negligible.展开更多
Nonlinear electromagnetic gyrokinetic equations have been constructed without expanding the field variables into background and finite but small-amplitude fluctuating components. At the long-wavelength limit, these un...Nonlinear electromagnetic gyrokinetic equations have been constructed without expanding the field variables into background and finite but small-amplitude fluctuating components. At the long-wavelength limit, these unexpanded nonlinear gyrokinetic equations recover the wellknown drift-kinetic equations. At the expanded limit, they recover the usual nonlinear gyrokinetic equations. These equations can therefore be applied to long-term simulations covering from microscopic to macroscopic spatial scales.展开更多
Fusion-born alpha particles in burning plasmas are usually regarded as have a slowing-down distribution,which differs significantly from the Maxwellian distribution of thermal particles in velocity space.A generalized...Fusion-born alpha particles in burning plasmas are usually regarded as have a slowing-down distribution,which differs significantly from the Maxwellian distribution of thermal particles in velocity space.A generalized multi-point average method has been developed for gyrokinetic Poisson equation with slowing-down equilibrium distribution using optimization in Fourier space.Its accuracy is verified in both long and short wavelength limits.The influence of changing equilibrium distribution from Maxwellian to slowing-down on gyrokinetic Poisson equation is analyzed to illustrate the significance of the new method.The effect of critical speed in the slowingdown distribution on the field solver is also presented.This method forms an important basis for global gyrokinetic simulation of low-frequency drift Alfvénic turbulence in burning plasmas.展开更多
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 densit...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 field 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.展开更多
In a tokamak fusion reactor operated at steady state,the equilibrium magnetic field is likely to have reversed shear in the core region,as the noninductive bootstrap current profile generally peaks off-axis.The revers...In a tokamak fusion reactor operated at steady state,the equilibrium magnetic field is likely to have reversed shear in the core region,as the noninductive bootstrap current profile generally peaks off-axis.The reversed shear Alfvén eigenmode(RSAE)as a unique branch of the shear Alfvén wave in this equilibrium,can exist with a broad spectrum in wavenumber and frequency,and be resonantly driven unstable by energetic particles(EP).After briefly discussing the RSAE linear properties in burning plasma condition,we review several key topics of the nonlinear dynamics for the RSAE through both wave-EP resonance and wave-wave coupling channels,and illustrate their potentially important role in reactor-scale fusion plasmas.By means of simplified hybrid MHD-kinetic simulations,the RSAEs are shown to have typically broad phase space resonance structure with both circulating and trapped EP,as results of weak/vanishing magnetic shear and relatively low frequency.Through the route of wave-EP nonlinearity,the dominant saturation mechanism is mainly due to the transported resonant EP radially decoupling with the localized RSAE mode structure,and the resultant EP transport generally has a convective feature.The saturated RSAEs also undergo various nonlinear couplings with other collective oscillations.Two typical routes as parametric decay and modulational instability are studied using nonlinear gyrokinetic theory,and applied to the scenario of spontaneous excitation by a finite amplitude pump RSAE.Multiple RSAEs could naturally couple and induce the spectral energy cascade into a low frequency Alfvénic mode,which may effectively transfer the EP energy to fuel ions via collisionless Landau damping.Moreover,zero frequency zonal field structure could be spontaneously excited by modulation of the pump RSAE envelope,and may also lead to saturation of the pump RSAE by both scattering into stable domain and local distortion of the continuum structure.展开更多
We have developed, based on the oscillating-center transformation, a general theoretical approach for self-consistent plasma dynamics including, explicitly, effects of nonlinear(higherorder) wave-particle resonances. ...We have developed, based on the oscillating-center transformation, a general theoretical approach for self-consistent plasma dynamics including, explicitly, effects of nonlinear(higherorder) wave-particle resonances. A specific example is then given for low-frequency responses of trapped particles in axisymmetric tokamaks. Possible applications to transport as well as nonlinear wave growth/damping are also briefly discussed.展开更多
The parameter dependence of transition between electrostatic instabilities is studied using gyrokinetic simulation based on a real discharge of steady-state scenario in the Experimental Advanced Superconducting Tokama...The parameter dependence of transition between electrostatic instabilities is studied using gyrokinetic simulation based on a real discharge of steady-state scenario in the Experimental Advanced Superconducting Tokamak.The scan of radial locations shows that trapped electron mode(TEM)dominates around the core while the ion temperature gradient mode(ITG)simultaneously dominates outside.The maximum growth rate of TEM appears aroundρ=0.24,where the maximum electron temperature gradient R/LTelocates,ρis the normalized poloidal flux.Effects of the parameters on the transition between TEM and ITG instability are studied atρ=0.24.It is found that TEM dominates in the scanning with individually changing R/LTe from 2.50 to 25.02 or the density gradient R/L_(n)from 1.38 to 13.76.Meanwhile,the electron-ion temperature ratio T_(e)/T_(i)is found to destabilize TEM,the effect of Teis more sensitive than that of T_(i).The dominant instability diagrams in the(R/L_(Te),R/L_(Ti))plane at different T_(e)/T_(i)and R/Lnare numerically obtained,which clearly show the parameter range of the dominant TEM or dominant ITG instability region.It is found that the dominant TEM region becomes narrower in the plane by decreasing R/L_(n)when T_(e)/T_(i)>0.5.展开更多
The parametric decay process of a reversed shear Alfvén eigenmeode(RSAE)into a geodesic acoustic mode and a kinetic RSAE is investigated using nonlinear gyrokinetic theory.The excitation conditions mainly require...The parametric decay process of a reversed shear Alfvén eigenmeode(RSAE)into a geodesic acoustic mode and a kinetic RSAE is investigated using nonlinear gyrokinetic theory.The excitation conditions mainly require the pump RSAE amplitude to exceed a certain threshold,which could be readily satisfied in burning plasmas operated in steady-state advanced scenario.This decay process can contribute to thermal plasma heating and confinement improvement.展开更多
In order to implement large-scale and high-beta tokamak simulation, a new algorithm of the electromagnetic gyrokinetic PIC （particle-in-cell） code was proposed and installed on the Gpic-MHD code [Gyrokinetic PIC cod...In order to implement large-scale and high-beta tokamak simulation, a new algorithm of the electromagnetic gyrokinetic PIC （particle-in-cell） code was proposed and installed on the Gpic-MHD code [Gyrokinetic PIC code for magnetohydrodynamic （MHD） simulation]. In the new algorithm, the vorticity equation and the generalized Ohm＇s law along the magnetic field are derived from the basic equations of the gyrokinetic Vlasov, Poisson, and Ampere system and are used to describe the spatio-temporal evolution of the field quantities of the electrostatic potential φ and the longitudinal component of the vector potential Az. The basic algorithm is equivalent to solving the reduced-MHD-type equations with kinetic corrections, in which MHD physics related to Alfven modes are well described. The estimation of perturbed electron pressure from particle dynamics is dominant, while the effects of other moments are negligible. Another advantage of the algorithm is that the longitudinal induced electric field, ETz = -δAz/δt, is explicitly estimated by the generalized Ohm＇s law and used in the equations of motion. Furthermore, the particle velocities along the magnetic field are used （vz-formulation） instead of generalized momentums （pz-formulation）, hence there is no problem of ＇cancellation＇, which would otherwise appear when Az is estimated from the Ampere＇s law in the pz-formulation. The successful simulation of the collisionless internal kink mode by the new Gpic-MHD with realistic values of the large-scale and high-beta tokamaks revealed the usefulness of the new algorithm.展开更多
In the present paper,we first derive the eigenmode equation of the ideal ballooning mode in tokamak plasmas using a gyrokinetic equation.It is shown that the gyrokinetic eigenmode equation can be reduced to the magnet...In the present paper,we first derive the eigenmode equation of the ideal ballooning mode in tokamak plasmas using a gyrokinetic equation.It is shown that the gyrokinetic eigenmode equation can be reduced to the magnetohydrodynamic（MHD） form in the long wavelength limit when kinetic effects are ignored.Then,the global gyrokinetic toroidal code（GTC） is applied for simulations of the edge-localized ideal ballooning modes.The obtained mode structures are compared with the results of ideal MHD simulations.The observed scaling of the linear growth rate with the toroidal mode number is consistent with the ideal MHD theory.The simulation results verify the GTC capability of simulating MHD processes in toroidal plasmas.展开更多
The theoretical and numerical studies on kinetic micro-instabilities,including ion temperature gradient（ITG） driven modes,trapped electron modes（TEMs） in the presence of impurity ions as well as impurity modes（IM...The theoretical and numerical studies on kinetic micro-instabilities,including ion temperature gradient（ITG） driven modes,trapped electron modes（TEMs） in the presence of impurity ions as well as impurity modes（IMs）,induced by impurity density gradient alone,in toroidal magnetized plasmas,such as tokamak and reversed-field pinch（RFP） are reviewed briefly.The basic theory for IMs,the electrostatic instabilities in tokamak and RFP plasmas are discussed.The observations of hybrid and coexistence of the instabilities are categorized systematically.The effects of impurity ions on electromagnetic instabilities such as ITG modes,the kinetic ballooning modes（KBMs） and kinetic shear Alfvén modes induced by impurity ions in tokamak plasmas of finite β（=plasma pressure/magnetic pressure） are analyzed.The interesting topics for future investigation are suggested.展开更多
A nonlinear saturation mechanism for reversed shear Alfvén eigenmode(RSAE)is proposed and analyzed,and is shown to be of relevance to typical reactor parameter region.The saturation is achieved through the genera...A nonlinear saturation mechanism for reversed shear Alfvén eigenmode(RSAE)is proposed and analyzed,and is shown to be of relevance to typical reactor parameter region.The saturation is achieved through the generation of high-frequency quasi-mode due to nonlinear coupling of two RSAEs,which is then damped due to coupling with the shear Alfvén continuum,and leads to the nonlinear saturation of the primary RSAEs.An estimation of the nonlinear damping rate is also provided.展开更多
基金supported by the National MCF Energy R&D Program of China(No.2019YFE03060000)National Natural Science Foundation of China(Nos.12005063,12375215 and 12175034)the Collaborative Innovation Program of Hefei Science Center,CAS(No.2022HSC-CIP008).
文摘The linear and nonlinear simulations are carried out using the gyrokinetic code NLT for the electrostatic instabilities in the core region of a deuterium plasma based on the International Thermonuclear Experimental Reactor(ITER)baseline scenario.The kinetic electron effects on the linear frequency and nonlinear transport are studied by adopting the adiabatic electron model and the fully drift-kinetic electron model in the NLT code,respectively.The linear simulations focus on the dependence of linear frequency on the plasma parameters,such as the ion and electron temperature gradientsκ_(Ti,e)≡R=L_(Ti,e),the density gradientκ_(n)≡R/L_(n)and the ion-electron temperature ratioτ=T_(e)=T_(i).Here,is the major radius,and T_(e)and T_(i)denote the electron and ion temperatures,respectively.L_(A)=-(δ_(r)lnA)^(-1)is the gradient scale length,with denoting the density,the ion and electron temperatures,respectively.In the kinetic electron model,the ion temperature gradient(ITG)instability and the trapped electron mode(TEM)dominate in the small and large k_(θ)region,respectively,wherek_(θ)is the poloidal wavenumber.The TEMdominant region becomes wider by increasing(decreasing)κ_(T_(e))(κ_(T_(i)))or by decreasingκ_(n).For the nominal parameters of the ITER baseline scenario,the maximum growth rate of dominant ITG instability in the kinetic electron model is about three times larger than that in the adiabatic electron model.The normalized linear frequency depends on the value ofτ,rather than the value of T_(e)or T_(i),in both the adiabatic and kinetic electron models.The nonlinear simulation results show that the ion heat diffusivity in the kinetic electron model is quite a lot larger than that in the adiabatic electron model,the radial structure is finer and the time oscillation is more rapid.In addition,the magnitude of the fluctuated potential at the saturated stage peaks in the ITGdominated region,and contributions from the TEM(dominating in the higher k_(θ)region)to the nonlinear transport can be neglected.In the adiabatic electron model,the zonal radial electric field is found to be mainly driven by the turbulent energy flux,and the contribution of turbulent poloidal Reynolds stress is quite small due to the toroidal shielding effect.However,in the kinetic electron model,the turbulent energy flux is not strong enough to drive the zonal radial electric field in the nonlinear saturated stage.The kinetic electron effects on the mechanism of the turbulence-driven zonal radial electric field should be further investigated.
基金supported by the China National Magnetic Confinement Fusion Science Program (Grant No. 2018YFE0304100)the US Department of Energy, Office of Science, Office of Advanced Scientific Computing Research and Office of Fusion Energy Sciences, Scientific Discovery through Advanced Computing (Sci DAC) program under Award Number DE-SC0018270 (Sci DAC ISEP Center)+2 种基金the China Scholarship Council (Grant No. 201306010032)resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory (DOE Contract No. DE-AC05-00OR22725)the National Energy Research Scientific Computing Center (DOE Contract No. DE-AC02-05CH11231)
文摘A conservative scheme of kinetic electrons for gyrokinetic simulations in the presence of magnetic islands has been implemented and verified in the gyrokinetic toroidal code, where zonal and nonzonal components of all perturbed quantities are solved together. Using this new conservative scheme, linear simulation of kinetic ballooning mode has been successfully benchmarked with the electromagnetic hybrid model. Simulations of nonlinear interactions between magnetic islands and the ion temperature gradient(ITG) mode in a tokamak show that the islands rotate at the electron diamagnetic drift velocity. The linear ITG structure shifts from the island O-point toward the X-point due to the pressure flattening effect inside the islands, and the nonlinear ITG structure peaks along the magnetic island separatrix because of the increased pressure gradient there.
基金supported by the National Key Research and Development Program of China (Grant No. 2017YFE0301201)partially by the National Natural Science Foundation of China (Grant Nos. U1967206 and 11775069)。
文摘The turbulence characteristics of plasmas with internal transport barriers in the HL-2A tokamak are analyzed by means of linear gyrokinetic simulations. It is found that turbulence is dominated by the ion temperature gradient(ITG)mode together with large-scale modes characterized by high-frequency electromagnetic fluctuation, which are destabilized by the steep ion temperature gradient in the weak magnetic shear regime. Comparison with solutions of analytical dispersion relations shows that their linear features match well with the beta-induced Alfvén eigenmode branch of the shear Alfvénic spectrum. It is further clarified that the large population of fast ions in these plasmas plays a stabilization role through the dilution mechanism in high-n ITG mode regimes.
基金supported by National Natural Science Foundation of China under Grant Nos.11505240,11375196 and 11405174the National ITER program of China under Contract No.2014GB113000
文摘Field-aligned coordinates have been implemented in the gyrokinetic semi-Lagrangian code NLT, Ye et al （2016 J. Comput. Phys. 316 180）, to improve the computational efficiency for the numerical simulations of tokamak turbulence and transport. 4D B-spline interpolation in field- aligned coordinates is applied to solve the gyrokinetic Vlasov equation. A fast iterative algorithm is proposed for efficiently solving the quasi-neutrality equation. A pseudo transform method is used for the numerical integration of the gyro-average operator for perturbations with a high toroidal mode number. The new method is shown to result in an improved code performance for reaching a given accuracy. Some numerical tests are presented to illustrate the new methods.
基金supported by the Chinese Scholarship Council(CSC)(No.201806340074)Shenzhen Clean Energy Research Institute and National Natural Science Foundation of China(No.12005141)+3 种基金supported by the US Department of Energy(No.DE-AC02-09CH11466)supported by the National MC Energy R&D Program(No.2018YFE0304100)National Key Research and Development Program(Nos.2016YFA0400600,2016YFA0400601 and 2016YFA0400602)the National Natural Science Foundation of China(Nos.11905220 and 11805273)。
文摘Gyrokinetic theory is arguably the most important tool for numerical studies of transport physics in magnetized plasmas.However,exact local energy–momentum conservation laws for the electromagnetic gyrokinetic system have not been found despite continuous effort.Without such local conservation laws,energy and momentum can be instantaneously transported across spacetime,which is unphysical and casts doubt on the validity of numerical simulations based on the gyrokinetic theory.The standard Noether procedure for deriving conservation laws from corresponding symmetries does not apply to gyrokinetic systems because the gyrocenters and electromagnetic field reside on different manifolds.To overcome this difficulty,we develop a high-order field theory on heterogeneous manifolds for classical particle-field systems and apply it to derive exact,local conservation laws,in particular the energy–momentum conservation laws,for the electromagnetic gyrokinetic system.A weak Euler–Lagrange(EL)equation is established to replace the standard EL equation for the particles.It is discovered that an induced weak EL current enters the local conservation laws,and it is the new physics captured by the high-order field theory on heterogeneous manifolds.A recently developed gauge-symmetrization method for high-order electromagnetic field theories using the electromagnetic displacement-potential tensor is applied to render the derived energy–momentum conservation laws electromagnetic gauge-invariant.
基金supported by the China National Magnetic Confinement Fusion Science Program(Grant No.2018YFE0304100)the US Department of Energy,Office of Science,Office of Advanced Scientific Computing Research and Office of Fusion Energy Sciences,and the Scientific Discovery through Advanced Computing(Sci DAC)program under Award No.DE-SC0018270(Sci DAC ISEP Center)+2 种基金the China Scholarship Council(Grant No.201806010067)used resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory(DOE Contract No.DEAC05-00OR22725)the National Energy Research Scientific Computing Center(DOE Contract No.DE-AC02-05CH11231)
文摘Global linear gyrokinetic simulations using realistic DIII-D tokamak geometry and plasma profiles find co-existence of unstable reversed shear Alfvén eigenmodes(RSAE)with low toroidal mode number n and electromagnetic ion temperature gradient(ITG)instabilities with higher toroidal mode number n.For intermediate n?=?[10,12],RSAE and ITG co-exist and overlap weakly in the radial domain with similar growth rates but different real frequencies.Both RSAE and ITG growth rates decrease less than 5%when compressible magnetic perturbations are neglected in the simulations.The ITG growth rates increase less than 7%when fast ions are not included in the simulations.Finally,the effects of trapped electrons on the RSAE are negligible.
基金supported by National Natural Science Fundation of China (Nos. 11235009 and 11905097)Fundamental Research Fund for Chinese Central Universities (No. 2019FZA3003)carried out within the framework of the EUROfusion Consortium and received funding from Euratom research and training programme 2014–2018 and 2019–2020 under Grant Agreement No. 633053 (Project No. WP19-ER/ENEA-05)。
文摘Nonlinear electromagnetic gyrokinetic equations have been constructed without expanding the field variables into background and finite but small-amplitude fluctuating components. At the long-wavelength limit, these unexpanded nonlinear gyrokinetic equations recover the wellknown drift-kinetic equations. At the expanded limit, they recover the usual nonlinear gyrokinetic equations. These equations can therefore be applied to long-term simulations covering from microscopic to macroscopic spatial scales.
基金the National Magnetic Confinement Fusion Program of China(No.2015GB110000)National Natural Science Foundation of China(No.11975201).
文摘Fusion-born alpha particles in burning plasmas are usually regarded as have a slowing-down distribution,which differs significantly from the Maxwellian distribution of thermal particles in velocity space.A generalized multi-point average method has been developed for gyrokinetic Poisson equation with slowing-down equilibrium distribution using optimization in Fourier space.Its accuracy is verified in both long and short wavelength limits.The influence of changing equilibrium distribution from Maxwellian to slowing-down on gyrokinetic Poisson equation is analyzed to illustrate the significance of the new method.The effect of critical speed in the slowingdown distribution on the field solver is also presented.This method forms an important basis for global gyrokinetic simulation of low-frequency drift Alfvénic turbulence in burning plasmas.
基金supported by National Special Research Program of China for ITER(Nos.2013GB111000 and 2014GB107004)China Scholarship Council(No.2011601098)U.S.DOE Grants DE-SC0010416 and DE-FG02-07ER54916
文摘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 field 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.
基金supported by National Natural Science Foundation of China (Nos. 12205251, 12275236 and 12261131622)Italian Ministry for Foreign Affairs and International Cooperation Project (No. CN23GR02)+2 种基金the National Key Research and Development Program of China (Nos. 2019YFE03020003 and 2017YFE0301900)Users of Excellence program of Hefei Science Center CAS (No. 2021HSC-UE016)funded by the European Union via the Euratom Research and Training Programme (No. 101052200–EUROfusion)
文摘In a tokamak fusion reactor operated at steady state,the equilibrium magnetic field is likely to have reversed shear in the core region,as the noninductive bootstrap current profile generally peaks off-axis.The reversed shear Alfvén eigenmode(RSAE)as a unique branch of the shear Alfvén wave in this equilibrium,can exist with a broad spectrum in wavenumber and frequency,and be resonantly driven unstable by energetic particles(EP).After briefly discussing the RSAE linear properties in burning plasma condition,we review several key topics of the nonlinear dynamics for the RSAE through both wave-EP resonance and wave-wave coupling channels,and illustrate their potentially important role in reactor-scale fusion plasmas.By means of simplified hybrid MHD-kinetic simulations,the RSAEs are shown to have typically broad phase space resonance structure with both circulating and trapped EP,as results of weak/vanishing magnetic shear and relatively low frequency.Through the route of wave-EP nonlinearity,the dominant saturation mechanism is mainly due to the transported resonant EP radially decoupling with the localized RSAE mode structure,and the resultant EP transport generally has a convective feature.The saturated RSAEs also undergo various nonlinear couplings with other collective oscillations.Two typical routes as parametric decay and modulational instability are studied using nonlinear gyrokinetic theory,and applied to the scenario of spontaneous excitation by a finite amplitude pump RSAE.Multiple RSAEs could naturally couple and induce the spectral energy cascade into a low frequency Alfvénic mode,which may effectively transfer the EP energy to fuel ions via collisionless Landau damping.Moreover,zero frequency zonal field structure could be spontaneously excited by modulation of the pump RSAE envelope,and may also lead to saturation of the pump RSAE by both scattering into stable domain and local distortion of the continuum structure.
基金supported by the ITER-CN under Grant No. 2017YFE0301900carried out within the framework of the EUROfusion Consortium and received funding from the Euratom research and training programme 2014-2018 and 2019-2020 under Grant Agreement No. 633053 (Project No. WP19-ER/ENEA-05)
文摘We have developed, based on the oscillating-center transformation, a general theoretical approach for self-consistent plasma dynamics including, explicitly, effects of nonlinear(higherorder) wave-particle resonances. A specific example is then given for low-frequency responses of trapped particles in axisymmetric tokamaks. Possible applications to transport as well as nonlinear wave growth/damping are also briefly discussed.
基金supported by the National MCF Energy R&D Program of China(Nos.2019YFE03060000,2019YFE03050000 and 2019YFE03020004)National Natural Science Foundation of China(Nos.12005063 and 11875131)+1 种基金Users with Excellence Program of Hefei Science Center CAS(Nos.2020HSC-UE011 and 2021HSC-UE015)Anhui Provincial Natural Science Foundation(No.2008085Jo4)。
文摘The parameter dependence of transition between electrostatic instabilities is studied using gyrokinetic simulation based on a real discharge of steady-state scenario in the Experimental Advanced Superconducting Tokamak.The scan of radial locations shows that trapped electron mode(TEM)dominates around the core while the ion temperature gradient mode(ITG)simultaneously dominates outside.The maximum growth rate of TEM appears aroundρ=0.24,where the maximum electron temperature gradient R/LTelocates,ρis the normalized poloidal flux.Effects of the parameters on the transition between TEM and ITG instability are studied atρ=0.24.It is found that TEM dominates in the scanning with individually changing R/LTe from 2.50 to 25.02 or the density gradient R/L_(n)from 1.38 to 13.76.Meanwhile,the electron-ion temperature ratio T_(e)/T_(i)is found to destabilize TEM,the effect of Teis more sensitive than that of T_(i).The dominant instability diagrams in the(R/L_(Te),R/L_(Ti))plane at different T_(e)/T_(i)and R/Lnare numerically obtained,which clearly show the parameter range of the dominant TEM or dominant ITG instability region.It is found that the dominant TEM region becomes narrower in the plane by decreasing R/L_(n)when T_(e)/T_(i)>0.5.
基金supported by the National Key R&D Program of China(No.2017YFE0301900)National Natural Science Foundation of China(No.11875233)Users of Excellence Program of Hefei Science Center CAS(No.2021HSC-UE016)。
文摘The parametric decay process of a reversed shear Alfvén eigenmeode(RSAE)into a geodesic acoustic mode and a kinetic RSAE is investigated using nonlinear gyrokinetic theory.The excitation conditions mainly require the pump RSAE amplitude to exceed a certain threshold,which could be readily satisfied in burning plasmas operated in steady-state advanced scenario.This decay process can contribute to thermal plasma heating and confinement improvement.
文摘In order to implement large-scale and high-beta tokamak simulation, a new algorithm of the electromagnetic gyrokinetic PIC （particle-in-cell） code was proposed and installed on the Gpic-MHD code [Gyrokinetic PIC code for magnetohydrodynamic （MHD） simulation]. In the new algorithm, the vorticity equation and the generalized Ohm＇s law along the magnetic field are derived from the basic equations of the gyrokinetic Vlasov, Poisson, and Ampere system and are used to describe the spatio-temporal evolution of the field quantities of the electrostatic potential φ and the longitudinal component of the vector potential Az. The basic algorithm is equivalent to solving the reduced-MHD-type equations with kinetic corrections, in which MHD physics related to Alfven modes are well described. The estimation of perturbed electron pressure from particle dynamics is dominant, while the effects of other moments are negligible. Another advantage of the algorithm is that the longitudinal induced electric field, ETz = -δAz/δt, is explicitly estimated by the generalized Ohm＇s law and used in the equations of motion. Furthermore, the particle velocities along the magnetic field are used （vz-formulation） instead of generalized momentums （pz-formulation）, hence there is no problem of ＇cancellation＇, which would otherwise appear when Az is estimated from the Ampere＇s law in the pz-formulation. The successful simulation of the collisionless internal kink mode by the new Gpic-MHD with realistic values of the large-scale and high-beta tokamaks revealed the usefulness of the new algorithm.
基金supported by U.S.Department of Energy(DOE) SciDAC GSEP Center and National Special Research Program of China for ITER
文摘In the present paper,we first derive the eigenmode equation of the ideal ballooning mode in tokamak plasmas using a gyrokinetic equation.It is shown that the gyrokinetic eigenmode equation can be reduced to the magnetohydrodynamic（MHD） form in the long wavelength limit when kinetic effects are ignored.Then,the global gyrokinetic toroidal code（GTC） is applied for simulations of the edge-localized ideal ballooning modes.The obtained mode structures are compared with the results of ideal MHD simulations.The observed scaling of the linear growth rate with the toroidal mode number is consistent with the ideal MHD theory.The simulation results verify the GTC capability of simulating MHD processes in toroidal plasmas.
基金supported by National Natural Science Foundation of China(Nos.11475057 and 11575158)the National Key R&D Program of China under Grant No.2017YFE0300405
文摘The theoretical and numerical studies on kinetic micro-instabilities,including ion temperature gradient（ITG） driven modes,trapped electron modes（TEMs） in the presence of impurity ions as well as impurity modes（IMs）,induced by impurity density gradient alone,in toroidal magnetized plasmas,such as tokamak and reversed-field pinch（RFP） are reviewed briefly.The basic theory for IMs,the electrostatic instabilities in tokamak and RFP plasmas are discussed.The observations of hybrid and coexistence of the instabilities are categorized systematically.The effects of impurity ions on electromagnetic instabilities such as ITG modes,the kinetic ballooning modes（KBMs） and kinetic shear Alfvén modes induced by impurity ions in tokamak plasmas of finite β（=plasma pressure/magnetic pressure） are analyzed.The interesting topics for future investigation are suggested.
基金supported by the Strategic Priority Research Program of Chinese Academy of Sciences(No.XDB0790000)the Collaborative Innovation Program of Hefei Science Center,CAS(No.2022HSC-CIP008)National Natural Science Foundation of China(Nos.12275236 and 12261131622)。
文摘A nonlinear saturation mechanism for reversed shear Alfvén eigenmode(RSAE)is proposed and analyzed,and is shown to be of relevance to typical reactor parameter region.The saturation is achieved through the generation of high-frequency quasi-mode due to nonlinear coupling of two RSAEs,which is then damped due to coupling with the shear Alfvén continuum,and leads to the nonlinear saturation of the primary RSAEs.An estimation of the nonlinear damping rate is also provided.
