The role of combined NNBI and ICRH heating for burning plasma physics studies in FAST H-mode scenarios. Alessandro Cardinali a
|
|
- Rose Rogers
- 5 years ago
- Views:
Transcription
1 The role of combined NNBI and ICRH heating for burning plasma physics studies in FAST H-mode scenarios Alessandro Cardinali a With the collaboration of: G. Calabrò a, M. Marinucci a, A. Bierwage c, G. Breyiannis a, S. Briguglio a, C. Di Troia a, G. Fogaccia a, P. Mantica d, G. Vlad a, X. Wang e, F. Zonca a, B. Baiocchi d, M. Baruzzo b, V. Basiuk f, R. Bilato g, M. Brambilla g, F. Imbeaux f, S. Podda a, M. Schneider f a) ENEA C.P I Frascati, Rome, Italy b) ENEA-Consorzio RFX, Padova, Italy c) Japan Atomic Energy Agency, Naka, Japan d) Istituto di Fisica del Plasma, Euratom-ENEA-CNR Association, Milan, Italy e) Institute for Fusion Theory and Simulation, Zhejiang University, Hangzhou, PRC f) CEA, IRFM, F Saint Paul Lez Durance, France g) Max-Planck-Institut für Plasmaphysik-Euratom Association, Garching, Germany 1
2 Summary Goal of the experiment High energy particle physics Description of the FAST tokamak Scenarii Description of the external power system NNBI ICRH The Physics of particle-particle & wave-particle interaction Quasi-linear Fokker-Planck equation Vlasov-Maxwell wave equation system The numerical tools Integrated scenarii Conclusions
3 Goal of the FAST experiment: High energy particle physics Preamble: The main difference between present days experiments and ITER will be the presence of alpha particles produced in DT reactions as the main heating source. Fusion alphas, with small characteristic orbit size (compared with the machine size), will mainly heat electrons and excite meso-scale fluctuations[f. Zonca]. In present day experiments relatively low energy fast ions, such as those generated by neutral beam injection and/or by ICRH minority heating, are characterized by medium or large orbit size and transfer their energy to thermal plasma ions while exciting macro-scale collective modes. The key control physical quantity for accessing burning plasma relevant dynamic regimes are plasma current, density and temperature, Plasma current is directly connected with particle orbit widths normalized to the machine size: Higher current means smaller fast particle dimensionless orbits. Relatively high density and low temperature, means lower critical energies above which fast particle predominantly heat electrons via collisional slowing down [Zonca F 008 Int. J. Mod. Phys. A ]. 3
4 Some quantities definitions & parameter requirements I Plasma radius normalized particle orbit width of fast particles: ˆρ fast Z fast 1/ M fast E fast m p a fast B T fast To be relevant for ITER, FAST experiment must have, as first requirement: ˆρ fast Z fast 1/ M fast E fast = ˆρ ITER m p a fast B T fast Z α 1/ M α E α ITER m p a ITER B T ITER In ITER the normalized orbit are ˆρ ITER = a ITER m α m proton E ITER(eV ) 0.01 ( cm) B T ITER ( Tesla)Z α ( gramms) m α 4
5 Some quantities definitions & parameter requirements II Plasma current scaling: a fast B T fast a ITER B T ITER = ξ I p fast I p ITER The requirement E fast E ITER = M α M fast Z fast Z α ξ I p fast I p ITER Recall us that in order to have a well confined orbit width and a relevant energy (above the critical energy) the plasma current must be I p fast = 1 ξ Z α Z fast M fast E fast M α E ITER I p ITER 5
6 Scaling of the plasma current for fast particles of different species 1 10 hydrogen deuterium helium3 Plasma current in MA Energy of fast particles in kev 6
7 The critical energy I To make clearer the concept of critical energy, and to quantify this concept, we write down some equations [D.V. Sivukhin, Reviews of Plasma Physics, Ed. Leontovich, Consultants Bureau New York, 4 93 (1966)]. From the theory of binary Coulomb collisions we can write the energy decrease of a test particle due to background species s with Maxwellian distribution function as: Ion dynamics dw b dt = 4πZ b e v b ( ) ( ) erf w ln Λ b,i n i Z i e i w i m i + m b m i m i m b π e w i + ( ) 4πZ b e erf w ln Λ b,e n e e e w e m e + m b v b m e m e m b π e w e b is the test particle; i (background ions); e(background electrons) ln Λ b,s is the Coulomb logarithm test particle-ions(electrons) collisions v b is the test particle-velocity-thermal-velocity-ratio w s = v ths ( ) Electron dynamics 7
8 The critical energy II The Sivukhin s equation simplifies considerably when: v thi << v b << v the For energies of test particles(deuterium) of about 1MeV, and plasma (deuterium) temperature of 10keV we have: v thd cm/sec v beam D cm/sec v the cm/sec dw b dt W b 4πZ b e 4 n e ln Λ b,e 3 π m e m b v the 1+ π Z 3 i m e n i v the m i n e W b 3/ m b 3/ ln Λ b,i ln Λ b,e Electron dynamics Ion dynamics 8
9 The critical energy III Defyining the critical energy W c ( kev ) = m b m p m e m p /3 ( π ) 1/3 ln Λ b,i ln Λ b,e /3 =14.8 M b i n i Z i i n e m p m i n i Z i n e 1 M i /3 /3 m p v the T e = ( kev ) the slowing down time τ sd = π Z b e 4 n e m b m e v the 1 ln Λ b,e 3 = 0.01 M 3/ bt e kev ( ) Z b n e10 14 ln Λ b,e = 9
10 We obtain the following equation The critical energy IV dw b dt = W b τ sd 1+ W c W b 3/ Which can be analytically integrated W b W c = 3/ W b0 3/ W +1 c e 3t τ sd 1 /3 W c If << 1 the test particle will be slowed-down by the electron species with the characteristic time t= sd / W b0 W c W b W c = W b0 W c e t τ sd If >> 1 the test particle will be slowed-down by ion species W b0 W b = 1 e W c 3t τ sd /3 10
11 The critical energy IV The fraction of energy transferred to the electrons and ions vs the ratio W is plotted showing the role c of the critical energy W b 11
12 On this basis a set of plasma and tokamak parameters have been chosen which are relevant for accessing burning plasma relevant dynamic regimes A device operating with deuterium plasmas in a dimensionless parameter range as close as possible to that of ITER and equipped with ICRH and NNBI as a main heating scheme would be able to address a number of important burning plasma physics issue. Information about the fundamental dynamic behaviour of fast ions in burning plasmas can be obtained by experimental studies of the fast ion tail produced by Ion Cyclotron Resonant Heating (ICRH) & NNBI injection. 1
13 FAST plasma parameters and scenarii Table 1: Operating scenarios H-mode & Extreme H-Mode scenarii of FAST are relevant for burning plasma physics studies The H-mode & Ext-H-Mode scenario of FAST have been obtained from the equilibrium calculated by MAXFEA code, and using iteratively the predictive 1.5-D transport code JETTO. To this end ICRH minority heating and NNBI heating are essential tools in FAST The FAST load assembling has been conceived to accommodate 30 MW of ICRH in 3 He minority heating scheme and 10MW of NBI (NNBI - Energy=0.7MeV (H) 1MeV (D)). FAST H-mode referente HMR H-mode extreme HME AT Full NICD I p (MA) q B T (T) H <n 0 > (m -3 ) β N τ E (s) τ Res (s) T 0 (kev) Q t discharge (s) t flat-top (s) I NI /I p (%) >100 P ADD (MW)
14 THE TOKAMAK In the Fusion Advanced Studies Torus (FAST) [A. Cardinali, et al.; Minority ions acceleration by ICRH: a tool for investigating burning plasma physics, Nucl. Fusion 49 (009) 09500], the reference and extreme H-mode scenarios require MW of external heating, mainly supplied by NNBI (10MW) and ICRH (30MW). These H-mode scenarios, both reference and extreme, are characterized by high magnetic field B=7.5-8,5T and high plasma current I p =6.5-8MA for a discharge time duration of about 13-0sec, with peak density.4-5x10^0m^-3 and temperature of the order of 10keV at the plasma centre. 14
15 Cryostat External Poloidal Coil Central Solenoid Toroidal Field Coil FAST load assembly view 15
16 Layout of the heating system: ICRH antenna Reference geometry (before optimization) with some relevant dimensions [m]. The antenna structure (figure) consists of a single cavity containing eight straps, two toroidal by four poloidal, protected by a Faraday Shield (FS) made of 30 conducting bars with squared crosssection. The analysis of the structure has been made with TOPICA (a computer code developed by the Plasma Facing Antenna Group of Politecnico di Torino), and 4MW per antenna has been calculated to be coupled to the plasma. 16
17 Power Spectra analysis Power spectra for the H-mode (left) and steady state(right) plasma profiles, provided a SOL length of 3 cm Power spectra for the H-mode (left) and steady state(right) plasma profiles, provided a SOL length of 3,4,5 cm 17
18 Layout of the heating system: NBI system M. Baruzzo, et al.first NBI Configuration Study for FAST Proposal, Proceedings of the 36th European Physical Society, Conference on Plasma Physics, Dublin, Ireland, June 9-July 3, 010, P4.194 (010). M. Baruzzo et al., Requirement specification for the neutral beam injector on FAST, 6th Symposium on Fusion Technology, Porto, Portugal, Sept FAST equatorial port and Neutral Beam (vertical section). FAST Neutral Beam Injector 18 Critical points between FAST equatorial port and Neutral Beam (horizontal section).
19 THE PHYSICS I The tokamak FAST, as said before, is essentially devoted to study the physics of fast particle: i) instabilities connected to this energetic ion population, ii) the radial transport of energetic ions due to collective mode excitations, iii) the coupling of meso-scalefluctuations with micro-turbulence,etc The key parameter for these studies is ß hot (thermal/magnetic energy density ratio) of the supra-thermal population, which can reach values up to 3%, well in line with the needs for exciting meso-scale fluctuations with the same characteristics of those expected in various regimes of reactor relevant conditions. ß hot, or its parallel and perpendicular component is defined as β hot = µ 0 p B β hot = µ 0 p B 1 = 4πµ 0m α B dµ dvv 4 µ f α ( v, µ )dvdµ 1 1 = πµ 0m α B dµ dvv 4 1 µ ( ) f α v,µ ( )dvdµ 19
20 THE PHYSICS II In order to reconstruct the ß hot profile it is necessary to calculate the hot-particle distribution function, which results from NNBI and/or ICRH interaction with the plasma. Therefore before doing the physics of high energy particles, me must deal with NNBI and ICRH plasma interaction physics. To this end a SIMULATION of NNBI+ICRH heating, in FAST H-mode and extreme-h-mode scenarii, is obtained by a D full wave electromagnetic code TORIC [M. BRAMBILLA, Plasma Phys. Control. Fusion 41 1 (1999)] combined with SSFPQL [M. BRAMBILLA, Nuclear Fusion (1994)], which solves the quasi-linear Fokker-Planck equation in -D velocity space in presence of ICRH and NBI. This heating code is iteratively coupled to the transport code for four time step in the profiles evolution, in order to have the self consistent evaluation of the plasma & power deposition profiles. 0
21 Interaction between Neutral beam and plasma particles: Fokker-Planck equation and distribution function I Concerning the NBI-plasma interaction, the starting point is the Fokker-Planck equation (D.V. Sivukhin; C.C.F. Karney) In spherical geometry we have: C( f b, f e ) + C( f b, f i ) = = Γb,e v v Φ ( u 1 f ew)v b v v + 1 v the + Γb,e 4v Ψ ( u ew) 3 + Γb,i v µ 1 µ ( ) f b µ + v Φ ( u 1 f iw)v b v v + 1 v thib + Γb,i 4v Ψ ( u w 3 i ) µ 1 µ µ ( ) f b f b t m α f b m e m b m i f b + + ( ) + C f b, f s s=i,e Φ u e w = Cδ ( v v birth ) where ( ) = erf ( u w) u w i i Φ u i w Ψ u i w ( ) = erf ( u i w) erf u w i Γ b,i = 4πn q i bq i ln Λ b,i m b ( ) = erf ( u w) u w e e Ψ u e w ( ) = erf ( u e w) erf u w e Γ b,e = 4πn q e bq e ln Λ b,e m b ( )er f ( u w) i u ib w ( ) ( u w)er f ( u w) i i u i w ( )erf '( u w) e u e w ( ) ( u w)erf '( u w) e e u e w u e = ( E beam m beam ) v the u ib = ( E beam m beam ) v thi w = v ( E beam m beam ) 1
22 Interaction between Neutral beam and plasma particles: Fokker-Planck equation and distribution function II To obtain the previous collisional operator we have used the following approximations: the background electrons are isotropic and Maxwellian. the background ions are isotropic and Maxwellian. The equation finally can be written at the steady state as Γ b,e w + Γb,i w w Φ( u w 1 f b e )w w w + v b w Φ u 1 f ( iw)w b w w + v b v thi m b f b v the m e m b f b m i + Γb,e 4w Ψ u w 3 e ( ) + Γb,i 4w Ψ u iw 3 ( ) µ 1 µ ( ) f b µ 1 µ µ + ( ) f b ( ) µ = δ N δt δ v v ( birth )δ µ µ birth 4πv birth ( )
23 Interaction between Neutral beam and plasma particles: Fokker-Planck equation and distribution function III To have some insight about the solution we make further assumptions: Isotropy µ = 0 v thi << v b << v the The equation assumes the relevant form: 1 τ sd v v v 3 3 ( + v c ) f b = δn birth δt ( t) roi ( ) 1 H v v birth 4πv v The solution is the classical slowing down distribution function f b = τ sd 4π 1 v 3 3 ( + v c ) δn birth ( t) δt roi H ( v v birth ) 3
24 Interaction between Neutral beam and plasma particles: Fokker-Planck equation and distribution function IV The distribution function is very different from a Maxwellian At large velocities it decays like 1/v 3 It becomes flat below the critical velocity It is not a good approximation of NBI It is a good approximation for the alpha-particle slowing down Anyhow, once known the distribution function we are able to calculate the Power Deposition Profiles of NBI and the hot parameter as defined before. These are the ingredients to perform the subsequent high energy particle physics. 4
25 Interaction between Ion cyclotron wave (ICRH) and plasma particles: Wave + Fokker-Planck equations and distribution function Maxwell-Vlasov non-linear integro-differential equation system Funzione di distribuzione E = 1 E c t 4πq α c t E = 1 B c t t f α + v f α + q α v B E + + E ext + m α c + vf α ( r, v,t)dv v B ext c v f α = 0 Campo elettrico & magnetico Campi esterni When considering collisions and mutual interaction between wave and particles we have (quasilinear-fokker-planck equation) d t f α = s C( f b, f s ) + C( f b, f b ) + S wave E ( ) + q s E dc m s f s 5
26 The previous equation system must be simplified on the basis of physics considerations If the field amplitudes are small such that W wave E << W thermal κt Fourier analysis for the time-dependence of the field (harmonic field) Linearization of the kinetic equation around a Maxwellian In the quasi-linear approximation (collisions) Drift approximation for orbits F α f α r, v,t e iωt r, v (,t) = F α r, v,t ( ) = Fα QL distorted ( ) + ( v,µ ) f α r, v,t ( ) B( r ) = B ( r 0 ) + Δr B ( r ) + 1 ( Δ r ) T B Δr +... For ICRH the cold plasma approximation (CPA) does not apply: Warm Plasma Limit. The presence of warm modes of propagation (like Bernstein modes) does not allow this limit. The integro-differential character of the equation cannot be removed. 6
27 Interaction between ICRH and plasma particles: Fokker-Planck equation and distribution function for minority species Concerning the ICRH-plasma interaction, the starting point is the Fokker-Planck equation (T.H. Stix, C.C.F. Karney, M. Brambilla) f min t s=species C( f min, f s ) + Q( f min ) + P( f min ) = 0 C( f min, f species ) = S collisions Q( f min ) = S wave P( f min ) = S Edc 7
28 The collision operator ( ) C f min, f s = species=electrons,ions = ν min,e w + ν min,i w w Φ u 1 f ( ew)w min w w + v th min v the w Φ u w 1 f ( i )w min w w + v th min v thi m min m e m min m i f min f min + ν min,e 4w Ψ u ew 3 + ν min,i 4w Ψ u w 3 i ( ) ( ) µ 1 µ µ 1 µ ( ) f min µ ( ) f min µ + Isotropic Maxwellian background of ions and electrons w = v v th min Φ( u i w) = erf ( u iw) u i w Φ( u e w) = erf u ew ( ) u e w ( )erf '( u i w) u i w ( )erf ' u e w u e w ( ) Ψ( u i w) = erf ( u i w) erf ( u iw) u i w Ψ( u e w) = erf ( u e w) erf ( u ew) u e w ( )erf ' u i w u i w ( ) ( )erf ' u e w u e w ( ) 8
29 Simplified quasi-linear diffusion operator due to ICRH In cylindrical geometry ( ) = S w = 1 ( v S w ) v Q f i = 1 v v v v D f i v where D = D ql ( v ) = D ql ( v, µ ) ( = D 0 E ) J p k v 1 µ Ω ci ( ) 1/ 9
30 Finally the equation to be solved at the steady state is: 1 w w Φ u 1 f ( ew)w min w w + v th min v the + Z w w Φ u 1 f ( iw)w min w w + v th min v thi 1 µ w ( ) m i m e m i m i f min f min w w D ql w, µ ν min,e v w fmin th min w µ f min µ + 1 4w Ψ u ew w ( ) Z 4w Ψ ( u iw) 3 µ µ 1 µ µ 1 µ ( ) f min µ 1 µ µ ( ) f min ( ) ( ) D ql w,µ µ + = ν min,e v w fmin th min w µ f min µ 30
31 Interaction between ICRH and plasma particles: Fokker-Planck equation and distribution function for minority species II To have some insight about the solution we make further assumptions: Isotropy µ = 0 v thi << v b << v the The equation assumes the relevant form: 1 w w 1 A w ( ) + w D( v, µ ) ql ν min,e v th min w F w w The solution is the classical Maxwellian distribution function for high energy particles found by Stix ln F w ( ) lnc = dw w 0 ( ) + w F( w) B( w) ( ) wb w ( ) ql w D w, µ + A w ν min,e v th min ( ) = 0 ( ) = Φ( u e w) + ZΦ( u i w) A w B w ( ) ( ) = Φ( u e w) T min T e + ZΦ( u i w) T min T i 31
32 It is possible to define a temperature of the minority tail T mintail = T e ( 1+ D ε) D = 3 A i A min D 0 1/ 4 3 π = 3 A i A min 1 P abs 1/ 4 3 π 4ν min/i n min T min ε = n e Z i n i m e m i T i T e 1/ 3
33 The numerical procedure in order to determine the distribution function of the heated (ICRH+NNBI) species and all the relevant quantities for the burning plasma physics analysis -The transport code JETTO (formation of the plasma scenario) -The wave code TORIC (which solves the D wave equation in general equilibrium geometry) -The D Fokker-Planck solver (which solves the D Fokker Planck equation for both ICRH and NBI heating system) -The HMGC and XHMGC codes which analyze the physics of energetic particles instabilities 33
34 The iteration process JETTO code t=7sec TORIC+SSQLFP codes Quasilinear Power Deposition Profiles JETTO code t=8sec TORIC+SSQLFP codes Quasilinear Power Deposition Profiles JETTO code t=9sec TORIC+SSQLFP codes Quasilinear Power Deposition Profiles JETTO code t=10sec TORIC+SSQLFP codes Quasilinear Power Deposition Profiles Steady State Reached EXIT HGCM code 34
35 Example of integrated calculation: Summary plots Summary plots show the evolution of the parallel(a) and perpendicular(b) effective temperatures, and of the tail minority fraction(c) for all the iteration procedure, corresponding to three different times of the discharge (7s,8s,9s ) and a total 30MW ICRH power coupled to the plasma. 35
36 The calculation of the ß hot for the sequence shown before 36
37 -hot is the relevant parameter in establishing the relevance of the FAST H-mode scenarii in order to study the energetic particles-induced Alfvenic instabilities. 3.5 d e n s i t y = 5 X cm -3 FAST P ICRH =10MW/m -3 ; n min =0.01; 3 He; B=7T d e n s i t y =. 3 X cm FAST Hmode extreeme n =.5X10 0 m -3 ; n e min =0.0; 3 He; B=8.5T 7. 5 power density=50m W / m beta_hot in % 1.5 beta_hot in % d e n s i t y = 5 X m power density=5m W / m -3 (a) electron temperature in kev electron temperature in kev Figures a)-b) - a) plot of -hot vs the electron temperature for several values of the plasma density at fixed ICRH power density and toroidal magnetic field for H-mode Scenario 1% of 3 He. b) as before but several values of power density and fixed density and magnetic field for the Extreme H-mode scenario % of 3 He. 37 (b)
38 hot as calculated by the ICRH full-wave code (TORIC-SSFPQL) and using the JETTO code (part of the JAMS JET suite of integrated codes with GLF3 transport model. For ICRH heating profiles we have used the TORIC code which is run outside JETTO and requires a few iterations FAST H-mode extreme scenario evolution of T_eff n=.94x10^0 m^-3; f=78mhz; B=8.5T; I_p=8MA; 3 H e = % FAST H-mode extreme scenario evolution of β_hot n=.94x10^0 m^-3; f=78mhz; B=8.5T; I_p=8MA; 3 H e = % temperature in kev T_perp t=9s T_para t=9s T_perp t=10s T_para t=10s T_perp t=11s T_para t=11s β_hot ICRH β_hot perp t=9s β_hot para t=9s β_hot perp t=10s β_hot para t=10s β_hot perp t=11s β_hot para t=11s x x Effective temperature of the energetic ion minority ß hot profile 38
39 When including NNBI (Deuterium beam of energy=1mev and P NNBI =10MW) on the same plasma target we have a PDP as calculated by ASCOT included in the JETTO code (part of the JAMS JET suite of integrated codes with GLF3 transport model, and consequently a calculated ß hot (parallel) NNBI Power Deposition Profiles [MW/m^3] electrons ions FAST H-mode extreme scenario 3rd iteration + NNBI n=.94x10^14 cm^-3; f=78mhz; T_e=10.578;T_i=9.87 kev; B=8.5T; I_p=8MA; 3He=% beta_hot-nnbi parallel FAST H-mode extreme scenario β_hot profile P_NNBI=10MW deuterium beam E=1MeV on ICRH heated plasma at P_ICRH=30MW n=.94x10^0 m^-3; f=78mhz; B=8.5T; I_p=8MA; 3 H e = % ro_pol Power deposition profiles of the beam ro_pol parallel ß hot profile 39
40 A comparison in the determination of the hot profiles is made by using two transport models for both ICRH and NNBI β hot perpendicular BgB+rotation β hot parallel BgB+rotation β hot perpendicular GLF3 β hot parallel GLF FAST H-mode extreme scenario β_hot profile P_NNBI=10MW deuterium beam E=1MeV on ICRH heated plasma at P_ICRH=30MW n=.94x10^0 m^-3; f=78mhz; B=8.5T; I_p=8MA; 3 H e = % beta_hot Scenario Hmode extreme third iteration. P ( I C R H ) = 3 0 M W + P ( N N B I = 1 0 M W on ohmic plasma target f=78mhz, n peak=10; fraction He3=1% or % beta comparison vs transport model: BgB+ rotation & GLF3 r e s o n a n c e s : 1) r=0.5(hfs) the fundamental of deuterium ) r=0.77(llfs) the fundamental of helium 3 beta_hot-nnbi parallel BgB transport model + rotation GLF3 transport model ro_pol ro_pol The BgB transport model seems to predict a broader perpendicular hot profile with respect to the GLF3 reaching at the peak the value of 1.1%, well in line with ITER expectation. Concerning the parallel hot profile (due to the NNBI power deposition) seems to be broader and much more peaked than the hot obtained with GLF3 transport model. 40
41 Changing the operation ICRH-frequency (f=8mhz) the deposition profile is on-axis. ß hot is optimized and can reach 3% temperature in kev perpendicular temperature parallel temperature FAST H-mode extreme scenario 3rd iteration n=3.x10^14 cm^-3; f=8mhz; T_e=15.4;T_i=14.6 kev; B=8.5T; I_p=8MA; He3=% β hot beta hot perpendicular beta hot parallel Scenario Hmode extreme 3rd iteration. P(ICRH)=30MW on ohmic plasma target f=8mhz, n peak=5; resonances: 1) r=0.6(hfs) the fundamental of deuterium ) r=0.1(llfs) the fundamental of helium x 0 ρ pol Effective temperature of the energetic ion minority ß hot profile 41
42 Once the FAST extreme H-mode scenario profiles are determined self-consistently with the iterative procedure described above, these are used as initial condition for numerical studies of resonant excitations of energetic particle driven Alfvénic fluctuations and related transports [3]. The FAST extreme H-mode scenario is characterized by a dense spectrum of Alfvénic fluctuations with the same wavelength and frequency spectra that are expected in ITER (peaked at 15 < n < 5). a) Application of the HMGC&XHMGC codes for the stability study b) c) -In Figs. a)-b)-c). a) Intensity contour plot of the n=8 Alfvénic mode in the space of normalized frequency and radial position. The shear Alfvén continuous spectrum for n=8 is indicated by the black solid line; b) Intensity contour plot of the n=16 Alfvénic mode in the space of normalized frequency and radial position. The shear Alfvén continuous spectrum for n=16 is indicated by the black solid line; c) Effect of n=8 and n=16 saturated EPM on the ß hot radial profile. -Both figures demonstrate the destabilization of EPM fluctuation branches, which track the toroidal precession resonance with the ICRH induced supra-thermal ion tail, as expected in this case and as predicted by theory. -The effect of n=8 and n=16 saturated EPM on the ß hot radial profile, shown in Fig. C), suggests that significant radial redistributions of energetic particle are expected with limited global losses. 4
43 Conclusions The proposed tokamak FAST in the H-mode scenarii is particularly indicated in order to investigate (experimentally) energetic particle physics: e.g. Alfvénic mode excitation, supra-thermal particle transports etc. Numerical tools like the hybrid MHD gyro-kinetic code HMGC and XHMGC can be used for predicting and interpreting the experimental data. A wide study of the wave-particle interaction is mandatory to establish if the conditions required for this goal are fulfilled by the experiment. For example the generation of a high energy population of ions with relevant value of hot. Numerical codes like TORIC+SSFPQL are devoted to study the interaction between ICRH and NBI with the plasma particles. The use of these numerical tools has been integrated in a global numerical investigation which includes a transport codes (like JETTO) with several transport models. The result of this analysis show that, as for the FAST H-mode reference scenario, the envisaged extreme H-mode scenario not only maximizes the machine performance and neutron yield, but is as well capable of addressing integrated experiments that are relevant in the view of ITER (e.g. preparation of ITER operations). 43
44 References Modeling Transport, ICRH and NBI. [1] A. Cardinali et al, Nucl. Fusion 49 (009) [] A. Pizzuto et al, Nucl. Fusion 50 (010) [3] G. Cenacchi and A. Taroni 1988 JETTO: A free boundary plasma transport code; report JET-IR (88) 03 [4] M. Brambilla, Plasma Physics and Controlled Nuclear Fusion 41 1 (1999) [5] M. Brambilla M 1994 Nucl. Fusion [6] Brambilla M 007 Nucl. Fusion [7] Brambilla M and Bilato R, submitted to Nucl. Fusion. [8] F Mirizzi et al., Fusion Engineering and Design 84 (009) [9] M Baruzzo, Proceedings of the 37th EPS Conference on Plasma Physics, Dublin, Ireland, June 010, P5-14 Modeling Instabilities. [1] Briguglio, S, Vlad G, Zonca F, and Kar C, 1995 Phys. Plasmas [] X. Wang et al., 011, Phys. of Plasmas, [3] L. Chen, F. Zonca, 007, Nucl. Fusion 47 S7. 44
MHD-particle simulations and collective alpha-particle transport: analysis of ITER scenarios and perspectives for integrated modelling
MHD-particle simulations and collective alpha-particle transport: analysis of ITER scenarios and perspectives for integrated modelling G. Vlad, S. Briguglio, G. Fogaccia, F. Zonca Associazione Euratom-ENEA
More informationEffects of Alpha Particle Transport Driven by Alfvénic Instabilities on Proposed Burning Plasma Scenarios on ITER
Effects of Alpha Particle Transport Driven by Alfvénic Instabilities on Proposed Burning Plasma Scenarios on ITER G. Vlad, S. Briguglio, G. Fogaccia, F. Zonca Associazione Euratom-ENEA sulla Fusione, C.R.
More informationFAST 1 : a Physics and Technology Experiment on the Fusion Road Map
Fusion Advanced Studies Torus FAST 1 : a Physics and Technology Experiment on the Fusion Road Map Presented by A. A. Tuccillo on behalf of ENEA-Euratom Association Univ. of Rome Tor Vergata Univ. of Catania
More informationImproved Plasma Confinement by Ion Bernstein Waves (IBWs) Interacting with Ions in JET (Joint European Torus)
Improved Plasma Confinement by Ion Bernstein Waves (IBWs) Interacting with Ions in JET (Joint European Torus) PD/P-01 C. Castaldo 1), R. Cesario 1), Y, Andrew 2), A. Cardinali 1), V. Kiptly 2), M. Mantsinen
More informationL Aquila, Maggio 2002
Nonlinear saturation of Shear Alfvén Modes and energetic ion transports in Tokamak equilibria with hollow-q profiles G. Vlad, S. Briguglio, F. Zonca, G. Fogaccia Associazione Euratom-ENEA sulla Fusione,
More informationFast ion generation with novel three-ion radiofrequency heating scenarios:
1 Fast ion generation with novel three-ion radiofrequency heating scenarios: from JET, W7-X and ITER applications to aneutronic fusion studies Yevgen Kazakov 1, D. Van Eester 1, J. Ongena 1, R. Bilato
More informationAlpha Particle Transport Induced by Alfvénic Instabilities in Proposed Burning Plasma Scenarios
Alpha Particle Transport Induced by Alfvénic Instabilities in Proposed Burning Plasma Scenarios G. Vlad, S. Briguglio, G. Fogaccia and F. Zonca Associazione Euratom-ENEA sulla Fusione, C.R. Frascati C.P.
More informationNon-ohmic ignition scenarios in Ignitor
Non-ohmic ignition scenarios in Ignitor Augusta Airoldi IFP, EURATOM-ENEA-CNR Association, Milano, Italy Francesca Bombarda, Giovanna Cenacchi Ignitor Group, ENEA, Italy Bruno Coppi MIT, USA DPP1 APS Meeting
More informationNonlinear Simulation of Energetic Particle Modes in JT-60U
TH/P6-7 Nonlinear Simulation of Energetic Particle Modes in JT-6U A.Bierwage,N.Aiba 2, K.Shinohara 2, Y.Todo 3,W.Deng 4,M.Ishikawa 2,G.Matsunaga 2 and M. Yagi Japan Atomic Energy Agency (JAEA), Rokkasho,
More informationIon energy balance during fast wave heating in TORE SUPRA
Ion energy balance during fast wave heating in TORE SUPRA Thierry Hutter, Alain Bécoulet, Jean-Pierre Coulon, Vincent Saoutic, Vincent Basiuk, G.T. Hoang To cite this version: Thierry Hutter, Alain Bécoulet,
More informationIntroduction to Plasma Physics
Introduction to Plasma Physics Hartmut Zohm Max-Planck-Institut für Plasmaphysik 85748 Garching DPG Advanced Physics School The Physics of ITER Bad Honnef, 22.09.2014 A simplistic view on a Fusion Power
More informationCurrent density modelling in JET and JT-60U identity plasma experiments. Paula Sirén
Current density modelling in JET and JT-60U identity plasma experiments Paula Sirén 1/12 1/16 Euratom-TEKES Euratom-Tekes Annual Seminar 2013 28 24 May 2013 Paula Sirén Current density modelling in JET
More informationEnergetic-Ion-Driven MHD Instab. & Transport: Simulation Methods, V&V and Predictions
Energetic-Ion-Driven MHD Instab. & Transport: Simulation Methods, V&V and Predictions 7th APTWG Intl. Conference 5-8 June 2017 Nagoya Univ., Nagoya, Japan Andreas Bierwage, Yasushi Todo 14.1MeV 10 kev
More informationDirect drive by cyclotron heating can explain spontaneous rotation in tokamaks
Direct drive by cyclotron heating can explain spontaneous rotation in tokamaks J. W. Van Dam and L.-J. Zheng Institute for Fusion Studies University of Texas at Austin 12th US-EU Transport Task Force Annual
More informationFast ion physics in the C-2U advanced, beam-driven FRC
Fast ion physics in the C-2U advanced, beam-driven FRC Richard Magee for the TAE Team 216 US-Japan Workshop on the Compact Torus August 23, 216! High β FRC embedded in magnetic mirror is a unique fast
More informationBenchmarking of electron cyclotron heating and current drive codes on ITER scenarios within the European Integrated Tokamak Modelling framework
Benchmarking of electron cyclotron heating and current drive codes on ITER scenarios within the European Integrated Tokamak Modelling framework L. Figini 1,a, J. Decker 2, D. Farina 1, N. B. Marushchenko
More informationICRF Mode Conversion Flow Drive on Alcator C-Mod and Projections to Other Tokamaks
ICRF Mode Conversion Flow Drive on Alcator C-Mod and Projections to Other Tokamaks Y. Lin, J.E. Rice, S.J. Wukitch, M.J. Greenwald, A.E. Hubbard, A. Ince- Cushman, L. Lin, E.S. Marmar, M. Porkolab, M.L.
More informationImproved Plasma Confinement by Ion Bernstein Waves (IBWs) Interacting with Ions in JET
EFDA JET CP(02)07/03 C. Castaldo, R. Cesario, Y, Andrew, A. Cardinali, V. Kiptly, M. Mantsinen, F. Meo, A. Murari, A. A. Tuccillo, M. Valisa, D. Van Eester, L. Bertalot, D. Bettella, C. Giroud, C. Ingesson,
More informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
More informationAnalysis of ICRF heating and ICRF-driven fast ions in recent JET experiments
Analysis of ICRF heating and ICRF-driven fast ions in recent JET experiments M.J. Mantsinen,,a, C. Challis 3, J. Eriksson, L.-G. Eriksson 5, D. Frigione 6, D. Gallart, J. Garcia 7, C. Giroud 3, T. Hellsten
More informationGlobal particle-in-cell simulations of Alfvénic modes
Global particle-in-cell simulations of Alfvénic modes A. Mishchenko, R. Hatzky and A. Könies Max-Planck-Institut für Plasmaphysik, EURATOM-Association, D-749 Greifswald, Germany Rechenzentrum der Max-Planck-Gesellschaft
More informationINTRODUCTION TO MAGNETIC NUCLEAR FUSION
INTRODUCTION TO MAGNETIC NUCLEAR FUSION S.E. Sharapov Euratom/CCFE Fusion Association, Culham Science Centre, Abingdon, Oxfordshire OX14 3DB, UK With acknowledgments to B.Alper for use of his transparencies
More informationTurbulent Transport Analysis of JET H-mode and Hybrid Plasmas using QuaLiKiz, TGLF and GLF23
EFDA JET CP(1)/ B. Baiocchi, J. Garcia, M. Beurkens, C. Bourdelle, F. Crisanti, C. Giroud, J. Hobirk, F. Imbeaux, I. Nunes, EU-ITM ITER Scenario Modelling group and JET EFDA contributors Turbulent Transport
More informationProgressing Performance Tokamak Core Physics. Marco Wischmeier Max-Planck-Institut für Plasmaphysik Garching marco.wischmeier at ipp.mpg.
Progressing Performance Tokamak Core Physics Marco Wischmeier Max-Planck-Institut für Plasmaphysik 85748 Garching marco.wischmeier at ipp.mpg.de Joint ICTP-IAEA College on Advanced Plasma Physics, Triest,
More informationActive Control of Alfvén Eigenmodes in the ASDEX Upgrade tokamak
Active Control of Alfvén Eigenmodes in the ASDEX Upgrade tokamak M. Garcia-Munoz, S. E. Sharapov, J. Ayllon, B. Bobkov, L. Chen, R. Coelho, M. Dunne, J. Ferreira, A. Figueiredo, M. Fitzgerald, J. Galdon-Quiroga,
More informationExcitation of Alfvén eigenmodes with sub-alfvénic neutral beam ions in JET and DIII-D plasmas
Excitation of Alfvén eigenmodes with sub-alfvénic neutral beam ions in JET and DIII-D plasmas D. Borba 1,9, R. Nazikian 2, B. Alper 3, H.L. Berk 4, A. Boboc 3, R.V. Budny 2, K.H. Burrell 5, M. De Baar
More informationEnergetic Particle Physics in Tokamak Burning Plasmas
Energetic Particle Physics in Tokamak Burning Plasmas presented by C. Z. (Frank) Cheng in collaboration with N. N. Gorelenkov, G. J. Kramer, R. Nazikian, E. Fredrickson, Princeton Plasma Physics Laboratory
More informationObservation of modes at frequencies above the Alfvén frequency in JET
Observation of modes at frequencies above the Alfvén frequency in JET F. Nabais 1, D. Borba 1, R. Coelho 1, L. Fazendeiro 1, J. Ferreira 1, A. Figueiredo 1, L. Fitzgerald 2, P. Rodrigues 1, S. Sharapov
More informationThe fast-ion distribution function
The fast-ion distribution function Source Collisions Orbits RF Losses W. Heidbrink 3 MeV & 14.7 MeV protons Charge Exchange Reactivity σv Complex neutral beam sources are described by a few parameters
More informationRELATIVISTIC EFFECTS IN ELECTRON CYCLOTRON RESONANCE HEATING AND CURRENT DRIVE
RELATIVISTIC EFFECTS IN ELECTRON CYCLOTRON RESONANCE HEATING AND CURRENT DRIVE Abhay K. Ram Plasma Science and Fusion Center Massachusetts Institute of Technology Cambridge, MA 02139. U.S.A. Joan Decker
More informationComparative Transport Analysis of JET and JT-60U Discharges
EFDA JET CP(1)/13 J. Garcia, N. Hayashi, G. Giruzzi, M. Schneider, E. Joffrin, S. Ide, Y. Sakamoto, T. Suzuki, H. Urano, the JT-U Team and JET EFDA contributors Comparative Transport Analysis of JET and
More informationToward the Realization of Fusion Energy
Toward the Realization of Fusion Energy Nuclear fusion is the energy source of the sun and stars, in which light atomic nuclei fuse together, releasing a large amount of energy. Fusion power can be generated
More informationElectron Transport and Improved Confinement on Tore Supra
Electron Transport and Improved Confinement on Tore Supra G. T. Hoang, C. Bourdelle, X. Garbet, T. Aniel, G. Giruzzi, M. Ottaviani. Association EURATOM-CEA. CEA-Cadarache, 38, St Paul-lez-Durance, France
More informationGenerating of fusion plasma neutron source with AFSI for Serpent MC neutronics computing Serpent UGM 2015 Knoxville, TN,
Generating of fusion plasma neutron source with AFSI for Serpent MC neutronics computing Serpent UGM 2015 Knoxville, TN, 14.10.2015 Paula Sirén VTT Technical Research Centre of Finland, P.O Box 1000, 02044
More informationSimulation Study of Interaction between Energetic Ions and Alfvén Eigenmodes in LHD
1 Simulation Study of Interaction between Energetic Ions and Alfvén Eigenmodes in LHD Y. Todo 1), N. Nakajima 1), M. Osakabe 1), S. Yamamoto 2), D. A. Spong 3) 1) National Institute for Fusion Science,
More informationStability Properties of Toroidal Alfvén Modes Driven. N. N. Gorelenkov, S. Bernabei, C. Z. Cheng, K. Hill, R. Nazikian, S. Kaye
Stability Properties of Toroidal Alfvén Modes Driven by Fast Particles Λ N. N. Gorelenkov, S. Bernabei, C. Z. Cheng, K. Hill, R. Nazikian, S. Kaye Princeton Plasma Physics Laboratory, P.O. Box 451, Princeton,
More informationIntegrated Modelling of ITER Scenarios with ECCD
Integrated Modelling of ITER Scenarios with ECCD J.F. Artaud, V. Basiuk, J. Garcia, G. Giruzzi*, F. Imbeaux, M. Schneider Association Euratom-CEA sur la Fusion, CEA/DSM/DRFC, CEA/Cadarache, 13108 St. Paul-lez-Durance,
More informationRelevant spatial and time scale in tokamaks. F. Bombarda ENEA-Frascati, FSN-FUSPHY-SAD
Relevant spatial and time scale in tokamaks F. Bombarda ENEA-Frascati, FSN-FUSPHY-SAD PolFusion - one day discussion Meeting, 23rd of July 2015 Ferrara Ignitor News MoU of April 2010 concerned the construction
More informationNonlinear Energetic Particle Transport in the Presence of Multiple Alfvénic Waves in ITER
Nonlinear Energetic Particle Transport in the Presence of Multiple Alfvénic Waves in ITER Mirjam Schneller, Philipp Lauber, Sergio Briguglio, Antti Snicker Acknowledgement M. Schneller 1, Ph. Lauber 1,
More informationProgress in understanding W control using ICRH in the JET-ILW tokamak
Progress in understanding W control using ICRH in the JET-ILW tokamak P.Mantica 20 th US-EU TTF Workshop, Salem, MA, US 28 th April-1 st May 2015 Acknowledgements P.Mantica1, C.Angioni2, R.Bilato2, F.Casson3,
More informationHeating and current drive: Radio Frequency
Heating and current drive: Radio Frequency Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 13 th February 2012 Dr Ben Dudson Magnetic Confinement Fusion (1 of 26)
More informationIntroduction to Fusion Physics
Introduction to Fusion Physics Hartmut Zohm Max-Planck-Institut für Plasmaphysik 85748 Garching DPG Advanced Physics School The Physics of ITER Bad Honnef, 22.09.2014 Energy from nuclear fusion Reduction
More informationTURBULENT TRANSPORT THEORY
ASDEX Upgrade Max-Planck-Institut für Plasmaphysik TURBULENT TRANSPORT THEORY C. Angioni GYRO, J. Candy and R.E. Waltz, GA The problem of Transport Transport is the physics subject which studies the physical
More informationCombined LH and ECH Experiments in the FTU Tokamak
Combined LH and ECH Experiments in the FTU Tokamak V. Pericoli Ridolfini, Y. Peysson 1, R. Dumont 1, G. Giruzzi 1, G. Granucci 2, L. Panaccione, L. Delpech 1, B. Tilia, FTU team, ECH team 2 Associazione
More informationSynergetic heating of D-NBI ions in the vicinity of the mode conversion layer in H-D plasmas in JET with the ITER like wall.
Synergetic heating of D-NBI ions in the vicinity of the mode conversion layer in H-D plasmas in JET with the ITER like wall. J. Ongena 1,*, Ye.O. Kazakov 1, Y. Baranov 2, C. Hellesen 3, J. Eriksson 3,
More informationSawtooth mixing of alphas, knock on D, T ions and its influence on NPA spectra in ITER plasma
Sawtooth mixing of alphas, knock on D, T ions and its influence on NPA spectra in ITER plasma F.S. Zaitsev 1, 4, N.N. Gorelenkov 2, M.P. Petrov 3, V.I. Afanasyev 3, M.I. Mironov 3 1 Scientific Research
More informationAnalysis and modelling of MHD instabilities in DIII-D plasmas for the ITER mission
Analysis and modelling of MHD instabilities in DIII-D plasmas for the ITER mission by F. Turco 1 with J.M. Hanson 1, A.D. Turnbull 2, G.A. Navratil 1, C. Paz-Soldan 2, F. Carpanese 3, C.C. Petty 2, T.C.
More informationEstimations of Beam-Beam Fusion Reaction Rates in the Deuterium Plasma Experiment on LHD )
Estimations of Beam-Beam Fusion Reaction Rates in the Deuterium Plasma Experiment on LHD ) Masayuki HOMMA, Sadayoshi MURAKAMI, Hideo NUGA and Hiroyuki YAMAGUCHI Department of Nuclear Engineering, Kyoto
More informationSelf-consistent modeling of ITER with BALDUR integrated predictive modeling code
Self-consistent modeling of ITER with BALDUR integrated predictive modeling code Thawatchai Onjun Sirindhorn International Institute of Technology, Thammasat University, Klong Luang, Pathumthani, 12121,
More informationEX8/3 22nd IAEA Fusion Energy Conference Geneva
P.C. de Vries JET-EFDA Culham Science Centre Abingdon OX14 3DB UK EX8/3 22nd IAEA Fusion Energy Conference Geneva P.C. de Vries1, E. Joffrin2,3, M. Brix1, C.D. Challis1, K. Crombé4, B. Esposito5, N.C.
More informationGlobal gyrokinetic modeling of geodesic acoustic modes and shear Alfvén instabilities in ASDEX Upgrade.
1 EX/P1-18 Global gyrokinetic modeling of geodesic acoustic modes and shear Alfvén instabilities in ASDEX Upgrade. A. Biancalani 1, A. Bottino 1, S. Briguglio 2, G.D. Conway 1, C. Di Troia 2, R. Kleiber
More informationSTUDY OF ADVANCED TOKAMAK PERFORMANCE USING THE INTERNATIONAL TOKAMAK PHYSICS ACTIVITY DATABASE
INTERNATIONAL ATOMIC ENERGY AGENCY 20th IAEA Fusion Energy Conference Vilamoura, Portugal, 1-6 November 2004 IAEA-CN-116/ STUDY OF ADVANCED TOKAMAK PERFORMANCE USING THE INTERNATIONAL TOKAMAK PHYSICS ACTIVITY
More informationHeating and Confinement Study of Globus-M Low Aspect Ratio Plasma
EX/P5- Heating and Confinement Study of Globus-M Low Aspect Ratio Plasma N.V. Sakharov ), V.V. Dyachenko ), B.B. Ayushin ), A.V. Bogomolov ), F.V. Chernyshev ), V.K. Gusev ), S.A. Khitrov ), N.A. Khromov
More informationDT Fusion Ignition of LHD-Type Helical Reactor by Joule Heating Associated with Magnetic Axis Shift )
DT Fusion Ignition of LHD-Type Helical Reactor by Joule Heating Associated with Magnetic Axis Shift ) Tsuguhiro WATANABE National Institute for Fusion Science, 322-6 Oroshi-cho, Toki 509-5292, Japan (Received
More informationSimulation Study of High-Frequency Magnetosonic Waves Excited by Energetic Ions in Association with Ion Cyclotron Emission )
Simulation Study of High-Frequency Magnetosonic Waves Excited by Energetic Ions in Association with Ion Cyclotron Emission ) Mieko TOIDA 1),KenjiSAITO 1), Hiroe IGAMI 1), Tsuyoshi AKIYAMA 1,2), Shuji KAMIO
More informationJoint ITER-IAEA-ICTP Advanced Workshop on Fusion and Plasma Physics October Introduction to Fusion Leading to ITER
2267-1 Joint ITER-IAEA-ICTP Advanced Workshop on Fusion and Plasma Physics 3-14 October 2011 Introduction to Fusion Leading to ITER SNIPES Joseph Allan Directorate for Plasma Operation Plasma Operations
More informationVerification of gyrokinetic particle simulation of Alfven eigenmodes excited by external antenna and by fast ions
Verification of gyrokinetic particle simulation of Alfven eigenmodes excited by external antenna and by fast ions L. Chen 1,2, W. Deng 1, Z. Lin 1, D. Spong 3, G. Y. Sun 4, X. Wang 2,1, X. Q. Xu 5, H.
More informationExperimental Evidence of Inward Momentum Pinch on JET and Comparison with Theory
Experimental Evidence of Inward Momentum Pinch on JET and Comparison with Theory Tuomas Tala, Association Euratom-Tekes, VTT, Finland JET-EFDA Culham Science Centre Abingdon, UK 22nd IAEA Fusion Energy
More informationSimulation of alpha particle current drive and heating in spherical tokamaks
Simulation of alpha particle current drive and heating in spherical tokamaks R. Farengo 1, M. Zarco 1, H. E. Ferrari 1, 1 Centro Atómico Bariloche and Instituto Balseiro, Argentina. Consejo Nacional de
More informationChapter 12. Magnetic Fusion Toroidal Machines: Principles, results, perspective
Chapter 12 Magnetic Fusion Toroidal Machines: Principles, results, perspective S. Atzeni May 10, 2010; rev.: May 16, 2012 English version: May 17, 2017 1 Magnetic confinement fusion plasmas low density
More informationSpectroscopic determination of the internal amplitude of frequency sweeping TAE
INSTITUTE OF PHYSICS PUBLISHING Plasma Phys. Control. Fusion 46 (2004) S47 S57 PLASMA PHYSICS AND CONTROLLED FUSION PII: S0741-3335(04)72680-9 Spectroscopic determination of the internal amplitude of frequency
More information1 THW/P7-09 Comparison of Quasi-linear and Exact Ion Cyclotron Resonant Heating Diffusion, With and Without Finite Width Ion Orbits Abstract.
1 Comparison of Quasi-linear and Exact Ion Cyclotron Resonant Heating Diffusion, With and Without Finite Width Ion Orbits R.W. Harvey1, Yu. V. Petrov1, D.L. Green2, E.F. Jaeger3, P.T. Bonoli4, D.B. Batchelor2,
More informationStudy of High-energy Ion Tail Formation with Second Harmonic ICRF Heating and NBI in LHD
21st IAEA Fusion Energy Conference Chengdu, China, 16-21 October, 2006 IAEA-CN-149/ Study of High-energy Ion Tail Formation with Second Harmonic ICRF Heating and NBI in LHD K. Saito et al. NIFS-851 Oct.
More informationAdditional Heating Experiments of FRC Plasma
Additional Heating Experiments of FRC Plasma S. Okada, T. Asai, F. Kodera, K. Kitano, T. Suzuki, K. Yamanaka, T. Kanki, M. Inomoto, S. Yoshimura, M. Okubo, S. Sugimoto, S. Ohi, S. Goto, Plasma Physics
More informationA Method of Knock-on Tail Observation Accounting Temperature Fluctuation Using 6 Li+T/D+T Reaction in Deuterium Plasma
A Method of Knock-on Tail Observation Accounting Temperature Fluctuation Using 6 Li+T/D+T Reaction in Deuterium Plasma Yasuko KAWAMOTO and Hideaki MATSUURA Department of Applied Quantum Physics and Nuclear
More informationLower Hybrid Current Drive Experiments on Alcator C-Mod: Comparison with Theory and Simulation
Lower Hybrid Current Drive Experiments on Alcator C-Mod: Comparison with Theory and Simulation P.T. Bonoli, A. E. Hubbard, J. Ko, R. Parker, A.E. Schmidt, G. Wallace, J. C. Wright, and the Alcator C-Mod
More informationAlcator C-Mod. Double Transport Barrier Plasmas. in Alcator C-Mod. J.E. Rice for the C-Mod Group. MIT PSFC, Cambridge, MA 02139
Alcator C-Mod Double Transport Barrier Plasmas in Alcator C-Mod J.E. Rice for the C-Mod Group MIT PSFC, Cambridge, MA 139 IAEA Lyon, Oct. 17, Outline Double Barrier Plasma Profiles and Modeling Conditions
More informationIntegration of Fokker Planck calculation in full wave FEM simulation of LH waves
Integration of Fokker Planck calculation in full wave FEM simulation of LH waves O. Meneghini S. Shiraiwa R. Parker 51 st DPP APS, Atlanta November 4, 29 L H E A F * Work supported by USDOE awards DE-FC2-99ER54512
More informationIntegrated Heat Transport Simulation of High Ion Temperature Plasma of LHD
1 TH/P6-38 Integrated Heat Transport Simulation of High Ion Temperature Plasma of LHD S. Murakami 1, H. Yamaguchi 1, A. Sakai 1, K. Nagaoka 2, H. Takahashi 2, H. Nakano 2, M. Osakabe 2, K. Ida 2, M. Yoshinuma
More informationFokker-Planck Modelling of NBI deuterons in ITER
Fokker-Planck Modelling of NBI deuterons in ITER V. Yavorskij 1,2), K. Schoepf 1), V. Goloborod ko 1,2), M. Cecconello 3,4), L.G. Eriksson 5), M. Khan 1), V. Kiptily 3), A. Korotkov 3), A. Polevoi 6),
More informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 2 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
More informationGyrokinetic Transport Driven by Energetic Particle Modes
Gyrokinetic Transport Driven by Energetic Particle Modes by Eric Bass (General Atomics) Collaborators: Ron Waltz, Ming Chu GSEP Workshop General Atomics August 10, 2009 Outline I. Background Alfvén (TAE/EPM)
More informationPhysics and Modelling of a Negative Ion Source Prototype for the ITER Neutral Beam Injection
1 ITR/P1-37 Physics and Modelling of a Negative Ion Source Prototype for the ITER Neutral Beam Injection J.P. Boeuf a, G. Fubiani a, G. Hagelaar a, N. Kohen a, L. Pitchford a, P. Sarrailh a, and A. Simonin
More informationPOWER DENSITY ABSORPTION PROFILE IN TOKAMAK PLASMA WITH ICRH
Dedicated to Professor Oliviu Gherman s 80 th Anniversary POWER DENSITY ABSORPTION PROFILE IN TOKAMAK PLASMA WITH ICRH N. POMETESCU Association EURATOM-MECTS Romania, University of Craiova, Faculty of
More informationProduction of Over-dense Plasmas by Launching. 2.45GHz Electron Cyclotron Waves in a Helical Device
Production of Over-dense Plasmas by Launching 2.45GHz Electron Cyclotron Waves in a Helical Device R. Ikeda a, M. Takeuchi a, T. Ito a, K. Toi b, C. Suzuki b, G. Matsunaga c, S. Okamura b, and CHS Group
More informationMagnetic Confinement Fusion and Tokamaks Chijin Xiao Department of Physics and Engineering Physics University of Saskatchewan
The Sun Magnetic Confinement Fusion and Tokamaks Chijin Xiao Department of Physics and Engineering Physics University of Saskatchewan 2017 CNS Conference Niagara Falls, June 4-7, 2017 Tokamak Outline Fusion
More information0 Magnetically Confined Plasma
0 Magnetically Confined Plasma 0.1 Particle Motion in Prescribed Fields The equation of motion for species s (= e, i) is written as d v ( s m s dt = q s E + vs B). The motion in a constant magnetic field
More informationOn Electron-Cyclotron Waves in Relativistic Non-Thermal Tokamak Plasmas
1 On Electron-Cyclotron Waves in Relativistic Non-Thermal Tokamak Plasmas Lj. Nikolić and M.M. Škorić Vinča Institute of Nuclear Sciences, P.O.Box 522, Belgrade 11001, Serbia and Montenegro ljnikoli@tesla.rcub.bg.ac.yu
More informationChapter IX: Nuclear fusion
Chapter IX: Nuclear fusion 1 Summary 1. General remarks 2. Basic processes 3. Characteristics of fusion 4. Solar fusion 5. Controlled fusion 2 General remarks (1) Maximum of binding energy per nucleon
More informationICRF Minority-Heated Fast-Ion Distributions on the Alcator C-Mod: Experiment and Simulation
ICRF Minority-Heated Fast-Ion Distributions on the Alcator C-Mod: Experiment and Simulation A. Bader 1, P. Bonoli 1, R. Granetz 1, R.W. Harvey 2, E.F. Jaeger 3, R. Parker 1, S. Wukitch 1. 1)MIT-PSFC, Cambridge,
More informationAbstract. see Appendix to IAEA-CN-69/OV1/2, The JET Team (presented by M.L. Watkins)
Abstract ALPHA PARTICLE STUDIES DURING JET DT EXPERIMENTS The JET Team (presented by P.R. Thomas) Joint European Torus, Abingdon, Oxfordshire, United Kingdom. The 997 DT experiment (DTE) at the Joint European
More informationLectures on basic plasma physics: Introduction
Lectures on basic plasma physics: Introduction Department of applied physics, Aalto University Compiled: January 13, 2016 Definition of a plasma Layout 1 Definition of a plasma 2 Basic plasma parameters
More informationExperimental Study of the Stability of Alfvén Eigenmodes on JET
IAEA FEC, Paper EX/P-7 Experimental Study of the Stability of Alfvén Eigenmodes on JET D.Testa,, A.Fasoli,, G.Fu 4, A.Jaun 3, D.Borba, P.de Vries 6, and JET-EFDA contributors [] Plasma Science and Fusion
More informationIon Cyclotron Emission from JET D-T Plasmas
JET P(98)7 G A Cottrell et al Ion Cyclotron Emission from JET D-T Plasmas This document is intended for publication in the open literature. It is made available on the understanding that it may not be
More informationDevelopment of LH wave fullwave simulation based on FEM
Development of LH wave fullwave simulation based on FEM S. Shiraiwa and O. Meneghini on behalf of LHCD group of Alacator C-Mod PSFC, MIT 2010/03/10 San-Diego, CA Special acknowledgements : R. Parker, P.
More informationDOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING
DOPPLER RESONANCE EFFECT ON ROTATIONAL DRIVE BY ION CYCLOTRON MINORITY HEATING V.S. Chan, S.C. Chiu, Y.A. Omelchenko General Atomics, San Diego, CA, U.S.A. 43rd Annual APS Division of Plasma Physics Meeting
More informationRole of Magnetic Configuration and Heating Power in ITB Formation in JET.
Role of Magnetic Configuration and Heating Power in ITB Formation in JET. The JET Team (presented by V. Parail 1 ) JET Joint Undertaking, Abingdon, Oxfordshire, United Kingdom 1 present address: EURATOM/UKAEA
More informationTH/P6-08. CompX, Del Mar, CA 92014, USA; 2Xcel Engineering, Oak Ridge, TN 37830, USA; 3PSFC-MIT,Cambridge, MA 02139
1 Differences Between QL and Exact Ion Cyclotron Resonant Diffusion R.W. Harvey1, Yu. V. Petrov1, E.F. Jaeger2, P.T. Bonoli3, A. Bader3, and RF-SciDAC Group 1 CompX, Del Mar, CA 92014, USA; 2Xcel Engineering,
More informationHeating and Current Drive by Electron Cyclotron Waves in JT-60U
EX/W- Heating and Current Drive by Electron Cyclotron Waves in JT-6U T. Suzuki ), S. Ide ), C. C. Petty ), Y. Ikeda ), K. Kajiwara ), A. Isayama ), K. Hamamatsu ), O. Naito ), M. Seki ), S. Moriyama )
More informationINTRODUCTION TO BURNING PLASMA PHYSICS
INTRODUCTION TO BURNING PLASMA PHYSICS Gerald A. Navratil Columbia University American Physical Society - Division of Plasma Physics 2001 Annual Meeting, Long Beach, CA 1 November 2001 THANKS TO MANY PEOPLE
More informationDer Stellarator Ein alternatives Einschlusskonzept für ein Fusionskraftwerk
Max-Planck-Institut für Plasmaphysik Der Stellarator Ein alternatives Einschlusskonzept für ein Fusionskraftwerk Robert Wolf robert.wolf@ipp.mpg.de www.ipp.mpg.de Contents Magnetic confinement The stellarator
More informationDetermination of q Profiles in JET by Consistency of Motional Stark Effect and MHD Mode Localization
Determination of q Profiles in JET by Consistency of Motional Stark Effect and MHD Mode Localization R. De Angelis 1), F. Orsitto1), M. Baruzzo 2), P. Buratti 1), B. Alper 3), L. Barrera 4), A. Botrugno
More informationAdvanced Tokamak Research in JT-60U and JT-60SA
I-07 Advanced Tokamak Research in and JT-60SA A. Isayama for the JT-60 team 18th International Toki Conference (ITC18) December 9-12, 2008 Ceratopia Toki, Toki Gifu JAPAN Contents Advanced tokamak development
More informationBenchmark of gyrokinetic, kinetic MHD and gyrofluid codes for the linear calculation of fast particle driven TAE dynamics
ITR/P-34 Benchmark of gyrokinetic, kinetic MHD and gyrofluid codes for the linear calculation of fast particle driven TAE dynamics A.Könies,S.Briguglio 2,N.Gorelenkov 3,T.Fehér 4,M.Isaev 5,P.Lauber 4,A.Mishchenko,
More informationNumKin, Strasbourg, October 17 th, 2016
F. Palermo 1 A.Biancalani 1, C.Angioni 1, F.Zonca 2, A.Bottino 1, B.Scott 1, G.D.Conway 1, E.Poli 1 1 Max Planck Institut für Plasmaphysik, Garching, Germany 2 ENEA C. R. Frascati - Via E. Fermi 45, CP
More informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More informationDiamond Neutral Particle Spectrometer at JET and proposal for ITER
Diamond Neutral Particle Spectrometer at JET and proposal for ITER Krasilnikov V.A. (vkrasilnikov@triniti.ru) SRC RF TRINITI, Troitsk, Moscow, Russia A compact fast corpuscular spectrometer with a detector
More informationProgress in Transport Modelling of Internal Transport Barrier Plasmas in JET
Progress in Transport Modelling of Internal Transport Barrier Plasmas in JET T. Tala 1, C. Bourdelle, F. Imbeaux, D. Moreau, V. Parail, G. Corrigan, F. Crisanti, X. Garbet, D. Heading, E. Joffrin, L. Laborde,
More informationMHD. Jeff Freidberg MIT
MHD Jeff Freidberg MIT 1 What is MHD MHD stands for magnetohydrodynamics MHD is a simple, self-consistent fluid description of a fusion plasma Its main application involves the macroscopic equilibrium
More informationDesign of next step tokamak: Consistent analysis of plasma flux consumption and poloidal field system
Design of next step tokamak: Consistent analysis of plasma flux consumption and poloidal field system J.M. Ané 1, V. Grandgirard, F. Albajar 1, J.Johner 1 1Euratom-CEA Association, Cadarache, France Euratom-EPFL
More information