1、A Kinetic-Fluid Model for Studying Thermal and Fast Particle Kinetic Effects on MHD Instabilities,C. Z. Cheng, N. Gorelenkov and E. BelovaPrinceton Plasma Physics Laboratory Princeton University,Outline,Energetic Particle Physics Issues Kinetic-MHD Model Advantages Limitations Linear and Nonlinear K
2、inetic-MHD codes Particle Characteristics and Kinetic Effects Nonlinear Kinetic-Fluid Model Summary,Why is Energetic Particle Physics Important?,Fast ions exist in all magnetic fusion devices and play essential roles in heating and current drive:- Fast ions in NBI, N-NBI, ICRH- Alpha particles produ
3、ced in D-T fusion reaction Significant loss of fast ions can lead to degradation of heating and current drive efficiency Lost fast ions tend to localize near outer midplane and may cause localized damage in first wall of toroidal reactors In Q 5 burning plasmas a-particles are dominant heating sourc
4、e because Pa Paux Control of fast ion pressure profile is important in controlling thermal plasma profiles, which affects global plasma stability and confinement Need to integrate energetic particle physics with global stability, confinement, and heating physics,Modeling Energetic Particle Physics,T
5、he difficulty of theoretical modeling stems from the disparatescales which traditionally are analyzed separately: global-scale phenomena are generally studied using MHD model, while microscale phenomena are described by kinetic theories. The kinetic-MHD model was developed by treating thermal partic
6、les by MHD model and fast particles by kinetic theories. Kinetic physics of both thermal and fast particles involve small spatial and fast temporal scales and can strongly affect the global structure and long time behavior of thermal plasmas and fast particles. A kinetic-fluid model has been develop
7、ed to treat kinetic physics of both thermal and fast particles, but also retains the framework of kinetic-MHD model, on which all present energetic particle codes are based.,Kinetic-MHD Model,Momentum Equation (Pc Ph):r / t + Vr V = rPc rPh + J B Continuity Equation (n nc, nh nc) :/ t + Vr r + rrV =
8、 0 Maxwells Equations:B/ t = rE, J = rB , rB = 0 Ohms Law: E + VB = 0, EB = 0 Adiabatic Pressure Law: / t + Vr (Pc/r5/3) = 0 Hot Particle Pressure Tensor:Ph = mh/2 s d3v vv fh(x,v) where fh is governed by gyrokinetic or Vlasov equations.,Advantages of Kinetic-MHD Model,Retains properly global geomet
9、rical effects such as gradients in P, B, etc. Covers most low-frequency waves and instabilities: 3 Branches of waves and instabilities:- Fast Magnetosonic Branch: compressional wvaes, mirror modes, etc.- Shear Alfven Branch: shear Alfven waves, ballooning, tearing, K-H instabilities, etc.- Slow Magn
10、etosonic Branch: sound waves, drift wave instabilities, etc. Retains hot particle kinetic physics.,Limitations of Kinetic-MHD Model,Assumes that fast particle density is negligible. Thermal particle dynamics is governed by MHD model.- Ohms law: plasma is frozen in B and moves with EB drift velocity
11、and parallel electric field vanishes. - Adiabatic pressure law: thermal plasma pressure changes adiabatically through plasma convection and compression. - Gyroviscosity, that contains ion gyroradius effects, and pressure anisotropy are ignored.- Thermal particle kinetic effects of gyroradii, trapped
12、 particle dynamics (transit, bounce and magnetic drift motions), and wave-particle resonances are ignored. Kinetic-MHD model for thermal plasmas is valid only when (a) wci w wt, wb, w*, wd(b) kL 1 and kri 1,PPPL Kinetic-MHD Codes,Linear Stability Codes- NOVA-K code: global TAE stability code with pe
13、rturbative treatment of non-MHD physics of thermal and fast particles- NOVA-2 code: global stability code with non-perturbative treatment of fast particle kinetic effects- HINST code: high-n stability code with non-perturbative treatment of fast particle kinetic effects Nonlinear Simulation Codes- M
14、3D-K code: global simulation code with fast particle kinetic physics determined by gyrokinetic equation. - HYM-1 code: global simulation code with fast particle kinetic physics determined by full equation of motion. - HYM-2 code: global hybrid simulation code with ions treated by full equation of mo
15、tion and electrons treated as massless fluid.,Kinetic Coupling Processes,Spatial scale coupling:- For k? ri O(1), ion motion is different from electron E B drift motion and large dEk can be produced.- For L? ri, particle magnetic moment is not an adiabatic invariant, ion motion is stochastic.- Banan
16、a orbit rB boundary layer width db and ri Temporal scale coupling:- If wb w, trapped particles will respond to an bounce orbit- averaged field - If wb, wt w, transit or bounce resonances are important for energy dissipation- If wd w, wave-particle drift resonance effects are important for energy dis
17、sipation- If wd w, particle magnetic drift motion dominates over E B drift,Typical Fusion Plasma Parameters,Typical Parameters of Magnetic Fusion DevicesB 5 T, ne 1014 cm-3, Ti,e 10 keV,LB, Lp a ( 1m), R/a 3, bc bh 10-2 Characteristic Scales of Core Particle Dynamicsri 3 mm, wci 3 108 sec-1,wte, wbe
18、 107 sec-1, wti, wbi 105-106 sec-1,w*i n 105 sec-1, wdi, wde n (104-105) sec-1 Fast Ions (nh nc, bh bc)a-Particles: 3.5 MeV; NBI-Particles: 100 keVN-NBI-Particles: 350 keV; ICRF Tail Ions: 1 MeVrh 1 - 5 cm, wth, wbh 106 - 107 sec-1,wdh n (105 - 106) sec-1,Temporal and Spatial Scale Orderings:- TAE M
19、odes: w VA/2qR 106 sec-1* For low-n modes: wte, wbe, w*h w wth, wbh wti, wbi, wdh k? rBh 1, k? rh k? rBi 1, k? ri wti, wbi k? rBh 1, k? rh k? rBi 1, k? ri 1- Internal Kink and Fishbone Modes: n = 1 and $w w*i, wdh 105 sec-1wte, wbe wth, wbh wti, wbi, wdh w wdi,e In the inertial layer: k? rh 1, k? ri
20、 1 Both thermal and fast particle kinetic effects are important in determining energetic particle physics.,Kinetic-Fluid Model Cheng & Johnson, J. Geophys. Res., 104, 413 (1999),Consider high-b multi-ion species plasmas in general magnetic field geometry Consider w wci, k?ri O(1) Mass Density Contin
21、uity Equation:/ t + Vr r + rrV = 0 Momentum Equation:(/t+ Vr) V = J B r j PjcmPjcm = mj s d3v (v V)(v V) fj Particle distribution functions f_j are determined from gyrokinetic (for w wci) or Vlasov (for w wci) equation. Maxwells equations in magnetostatic limit are employed.,Pressure Tensor and Gyro
22、viscosity:P = P? (I - bb) + Pk bb + Pwhere I is the unit dyadic and b = B/B.Pk = m s d3v vk2 f, P? = (m/2) s d3v v?2 fFor k? kk, gyroviscosity tensor contribution rP b (rdPc b) + b r?dPs dPc = dPc1 + dPc2 , dPc1 = sd3v (m v?2/2) g0 (J0 2 J10) dPc2 = s d3v (m v?2 /2) (q/mB) F/m (F vk Ak)(2J0J10 J02)
23、(v?d Bk /k?)(J0 J1 2 J1 J10)dPs = s d3v (i mv?2 /l2) (qF/T)(w0 - w*T)/wc (q/mB) F/m (w- kk vk - wd)/wc(F vkAk) (l J0 J1 + J02 - 1) (v?dBk/2k?)l(1 2 J12) 2 J0 J1w0 = -(Tw/m) ln F/e, l = k? v? / wc,Low-Frequency Ohms LawE + VB = (1/nee) JB r( Pecm i (qi me/e mi) Picm) + i (mi/rqi 1/nee)(B/B) (r Pi0 B/
24、B) + (me/nee2) J/ t + r(JV + VJ) + hJwhere Pi0 = mi s d3v vv fi Main Features:- The kinetic-fluid model retains most essential particle kinetic effects in low frequency phenomena (w wci) for all particle species- Gyroviscosity is included so that ion Larmor radius effects are properly retained- A ne
25、w Ohms law for multi-ion species - No assumption on nh/nc ratio- Nonlinear,Kinetic-Fluid Codes,Linear Stability Codes- Extend non-perturbative global NOVA-2 code to include both thermal and fast particle kinetic effects- Extend high-n HINST code with non-perturbative treatment of both thermal and fa
26、st particle kinetic effects Nonlinear Global Simulation Codes- Extend M3D-K and HYM codes to include both thermal and fast particle kinetic effects.,Integration of Burning Plasmas Physics,Auxiliary Heating Fueling Current Drive,P(r), n(r), q(r),Confinement, Disruption Control MHD Stability,Pa Paux,F
27、ast Ion Driven Instabilities Alpha Transport,a interaction with thermal plasmas is a strongly nonlinear process.,Must develop efficient methods to control profiles for burn control!,Fusion Output,a-Heating a-CD,Summary,A nonlinear kinetic-fluid model has been developed for high-b plasmas with multi-
28、ion species for w wci. Physics of wave-particle interaction and geometrical effects are properly included, and the kinetic-fluid model includes kinetic effects of both thermal and fast particles. Eigenmode equations for dispersive shear Alfven waves and kinetic ballooning modes derived from the kine
29、tic-fluid model agree with those derived from gyrokinetic equations for w wci. Based on the kinetic-MHD model global and high-n linear stability codes (e.g., NOVA-K, NOVA-2, HINST, etc.) and nonlinear simulation codes (e.g., M3D-K, HYM codes) have been developed to study effects of energetic particles on MHD modes such as TAEs, internal kinks, etc. Linear stability and nonlinear simulation codes based on the kinetic-fluid model can be constructed by extending these existing kinetic-MHD codes.,
