Nonlinear magnetohydrodynamic effects on Alfvén eigenmode evolution and zonal flow
|
|
- Arline Griffith
- 5 years ago
- Views:
Transcription
1 Home Search Collections Journals About Contact us My IOPscience Nonlinear magnetohydrodynamic effects on Alfvén eigenmode evolution and zonal flow generation This article has been downloaded from IOPscience. Please scroll down to see the full text article Nucl. Fusion ( View the table of contents for this issue, or go to the journal homepage for more Download details: IP Address: The article was downloaded on 07/10/2010 at 22:01 Please note that terms and conditions apply.
2 IOP PUBLISHING and INTERNATIONAL ATOMIC ENERGY AGENCY Nucl. Fusion 50 (2010) (9pp) NUCLEAR FUSION doi: / /50/8/ Nonlinear magnetohydrodynamic effects on Alfvén eigenmode evolution and zonal flow generation Y. Todo 1,2, H.L. Berk 3 and B.N. Breizman 3 1 National Institute for Fusion Science, Toki, Gifu , Japan 2 The Graduate University for Advanced Studies (SOKENDAI), Toki, Gifu , Japan 3 Institute for Fusion Studies, University of Texas at Austin, Austin, TX 78712, USA todo@nifs.ac.jp Received 13 November 2009, accepted for publication 16 February 2010 Published 28 July 2010 Online at stacks.iop.org/nf/50/ Abstract Nonlinear magnetohydrodynamic (MHD) effects on Alfvén eigenmode evolution were investigated via hybrid simulations of an MHD fluid interacting with energetic particles. The investigation focused on the evolution of an n = 4 toroidal Alfvén eigenmode (TAE) which is destabilized by energetic particles in a tokamak. In addition to fully nonlinear code, a linear-mhd code was used for comparison. The only nonlinearity in that linear code is from the energetic-particle dynamics. No significant difference was found in the results of the two codes for low saturation levels, δb/b In contrast, when the TAE saturation level predicted by the linear code is δb/b 10 2, the saturation amplitude in the fully nonlinear simulation was reduced by a factor of 2 due to the generation of zonal (n = 0) and higher-n (n 8) modes. This reduction is attributed to the increased dissipation arising from the nonlinearly generated modes. The fully nonlinear simulations also show that geodesic acoustic mode is excited by the MHD nonlinearity after the TAE mode saturation. PACS numbers: Bj, Mw, Pi, Ww (Some figures in this article are in colour only in the electronic version) 1. Introduction The destabilization of Alfvén eigenmodes by energetic particles is an important concern for burning plasmas since the excited modes can enhance transport and losses of energetic ions. Computer simulation is a powerful tool to investigate the interaction between Alfvén eigenmodes and energetic particles. We performed the first numerical demonstration of toroidal Alfvén eigenmode (TAE) bursts [1] with parameters similar to a TFTR experiment [2] and many of the experimental characteristics were reproduced. These include (a) the synchronization of multiple TAEs, (b) the modulation depth of the drop in the stored beam energy and (c) the stored beam energy. However, the saturation amplitude was δb/b , which is higher than the value δb/b 10 3 inferred from the experimental plasma displacement measurements [1, 3]. In the simulation of [1], the only nonlinearity retained was the nonlinearity in the energetic-particle orbits, while the nonlinear magnetohydrodynamic (MHD) effects were neglected. Thus the spatial profiles and damping rates of the TAEs were assumed to be independent of the mode amplitude. Such a simulation is justified when particle trapping is the most important saturation mechanism for the TAE instability [4, 5]. This saturation mechanism has been demonstrated in several previous simulations [6 9]. However, the linear simulation tends to give too large a saturation level. This suggests that the MHD nonlinearity neglected in the TAE burst simulation may be important and is well worth investigating carefully. In another simulation study of TAE bursts, where the MHD nonlinear effects are taken into account but the parameters are not very close to the experiment, the saturation level is roughly δb/b [10]. The effects of MHD nonlinearities on the increase in TAE damping rate were investigated theoretically [11] and numerically [12]. These studies further motivated us to investigate the nonlinear MHD effects on the TAE evolution to better understand the physics mechanism of saturation. In this work we study the nonlinear evolution of a single linearly unstable TAE mode with toroidal mode number n = 4 using two versions of the MEGA code [13 15] which is a hybrid simulation code for an MHD fluid interacting with energetic particles. With the standard version of the MEGA code, the full nonlinear dynamics of both the MHD fluid and the energetic particles is simulated. In the other version of /10/ $ IAEA, Vienna Printed in the UK & the USA
3 the MEGA code, only linear-mhd equations are used while the nonlinear particle dynamics is retained. We shall refer to this version as the linear-mhd simulation. A comparison of the results between the nonlinear MHD and linear-mhd simulations is used to clarify how the MHD nonlinearities affect the TAE evolution. We demonstrate that the nonlinearly generated modes, both the zonal (n = 0) and the higher-n (n 8) modes, reduce the TAE saturation level. To understand the underlying physics mechanism, we analysed the energy dissipation from each toroidal mode number. We found that the total energy dissipation is significantly increased by the nonlinearly generated modes. The increase in the total energy dissipation reduces the TAE saturation level. The spatial profiles and the evolution of the nonlinearly generated n = 0 and n = 8 modes are investigated. After the saturation of the TAE instability a geodesic acoustic mode (GAM) is observed, driven by the MHD nonlinearity. 2. The simulation model Several hybrid simulation models have been constructed [7, 16 18] to study the evolution of Alfvén eigenmodes destabilized by energetic particles. In the MEGA code, the bulk plasma is described by the nonlinear MHD equations and the energetic ions are simulated with the δf particle method. The MHD equations with the energetic-ion effects are given by ρ = (ρv) + ν n ( ) ρ ρ eq, (1) ρ ( v 2 ) v = ρ ω v ρ p + ( j j 2 h) B B p + 4 (νρ v) (νρ ω), (2) 3 = E, (3) = (pv) (γ 1) p v + (γ 1) [νρω νρ ( v)2 + ηj (j ) ] j eq + ν n ( p p eq ), (4) E = v B + η ( j j eq ), (5) j = 1 µ 0 B, (6) ω = v, (7) where µ 0 is the vacuum magnetic permeability, γ is the adiabatic constant, ν and ν n are artificial viscosity and diffusion coefficients chosen to maintain numerical stability and all the other quantities are conventional. The subscript eq represents the equilibrium variables. The energetic-ion contribution is included in the MHD momentum equation (equation (2)) as the energetic-ion current density. The quantity j h is the energetic-ion current density without E B drift. We see that electromagnetic field is given by the standard MHD description. This model is accurate under the condition that Figure 1. Spatial profiles of energetic-ion beta, bulk plasma beta and safety factor. the energetic-ion density is much less than the bulk plasma density. The MHD equations are solved using a fourth order (in both space and time) finite difference scheme. The driftkinetic description [19] is employed for the energetic ions. The energetic-ion current density j h in equation (2) includes the contributions from parallel velocity, magnetic curvature and gradient drifts and magnetization current. The E B drift disappears in j h due to the quasi-neutrality [13]. The computational particles are initially loaded uniformly in the phase space. For the purpose of clarifying the nonlinear MHD effects, the linear-mhd calculations were performed to solve the following equations: ρ = (ρ eq v ) + ν n ( ) ρ ρ eq, (8) ( ) ρ eq v = p + j eq j heq δb + ( δj δj ) h Beq ( νρ eq v ) ( νρ eq ω ), (9) B = E, (10) p = (p eq v ) (γ 1) p eq v + (γ 1) ηδj j eq + ν n ( p p eq ), (11) E = v B eq + ηδj. (12) Here the variables with δ such as δb represent the fluctuations, for example, δb = B B eq. A tokamak plasma with the aspect ratio of R 0 /a = 3.2 was investigated where R 0 and a represent the major radius of the geometrical centre of the simulation domain and the plasma minor radius, respectively. The cylindrical coordinates (R,ϕ,z)are employed. The shape of the outermost magnetic surface is circular. The spatial profiles of the energetic-ion beta, bulk plasma beta and safety factor are shown in figure 1. The bulk plasma density is uniform. The different initial values of the central energetic-ion beta β h0 were investigated with the profile kept constant. The central energetic-ion beta β h0 shown in figure 1 is 2.0%. The initial velocity-space distribution 2
4 Figure 3. Comparison of radial velocity evolution for the linear-mhd and the nonlinear MHD runs using the cosine part of m/n = 6/4 harmonics at r/a = 0.42 for β h0 = 1.5%. Figure 2. Spatial profiles of each poloidal harmonic of the TAE with toroidal mode number n = 4 for (a) radial velocity and (b) radial magnetic field. Solid (dashed) lines show cos(mϑ + nϕ) [sin(mϑ + nϕ)] harmonics with poloidal mode number m labelled in the figure. of the energetic ions is a slowing-down distribution with a maximum velocity 1.2v A and the critical velocity 0.5v A, where v A denotes the Alfvén velocity at the plasma centre. The ratio of the energetic-ion Larmor radius to the minor radius is 1/16 for the energetic-ion velocity equal to the Alfvén velocity. The number of grid points for the cylindrical coordinates (R,ϕ,z) is and the number of computational particles is We investigate the evolution of the n = 4 mode, which is an exact solution of the equations of a quarter of the tokamak domain with the toroidal angle taken from 0 ϕ π/2. The MHD nonlinearities generate fluctuations with toroidal mode numbers that are multiples of 4 (n = 0, 4, 8, 12, 16,...). In contrast, in the linear-mhd code only the n = 4 component is excited. The nonlinear evolution of the energetic ions generates the current density fluctuation δj h with a harmonic content in ϕ, with all harmonics being a multiple of 4. This leads to the generation of MHD fluctuations with toroidal mode numbers that are multiples of n = 4 and obscures the effects of the MHD nonlinearity. Then, we retain only the n = 4 harmonic of the hot particle current, while we artificially remove the energetic-ion current density fluctuation δj h if n 4. The viscosity, diffusion and resistivity coefficients in the MHD equations are chosen to be ν = ν n = 10 6 v A R 0 and η = 10 6 µ 0 v A R 0, respectively. The spatial profile of the unstable n = 4 TAE mode observed at the linearly growing phase of the instability for β h0 = 1.0% is shown in figure 2. The phase of the mode in the figure is chosen so that the cosine part of the dominant harmonic m/n = 6/4 is maximized at the peak location. The frequency and growth rate of the TAE are ω = 0.318ω A,γ = ω A with ω A = v A /R axis, where R axis is the major radius of the magnetic axis. The maximum values of v r /v A and δb r /B are close to each other for the TAE, Figure 4. Comparison of radial velocity evolution for the linear-mhd and the nonlinear MHD runs using the cosine part of m/n = 6/4 harmonics at r/a = 0.42 for β h0 = 2.0%. because a relation between velocity and magnetic fluctuations v r /v A = δb r /B holds for shear Alfvén waves in a uniform plasma. The major harmonics (m/n = 5/4 and 6/4) of the radial magnetic field profiles have opposite signs whereas the signs of the radial velocity harmonics are the same. The parallel wave number k = (n m/q)/r 0 is positive for m/n = 5/4 and negative for m/n = 6/4 in the region 5/4 < q < 6/4 (0.35 < r/a < 0.50). The opposite sign of the parallel wave number leads to the opposite sign of the magnetic field fluctuation. For the purpose of the data analysis, magnetic flux coordinates (r, ϕ, ϑ), where r is the radial coordinate and ϑ is the poloidal angle, were constructed for the MHD equilibrium. 3. Simulation results 3.1. Comparison of linear-mhd and nonlinear MHD simulations The evolution of the MHD radial velocity is compared in figures 3 and 4 for the linear-mhd simulation and the nonlinear MHD simulation. A comparison of the radial velocity excitation is made for two different hot particle beta values, β h0 = 1.5% and β h0 = 2.0%. Figures 3 and 4 show the evolution of the m/n = 6/4 harmonics of radial velocity v r /v A measured at the TAE peak location r = 0.42a. The frequency and growth rate are ω = 0.310ω A,γ= ω A for β h0 = 1.5%, and ω = 0.303ω A,γ= ω A for β h0 = 2.0%. We see in figure 3 that the saturation levels 3
5 3.2. Physics mechanism of saturation level reduction In this subsection, we discuss the physics mechanism that reduces the TAE saturation level. We focus on the energy dissipation that appears in equation (4) from the viscous heating and the Joule heating terms. They are, respectively, in proportion to viscosity ν and resistivity η. We analysed the time evolution of the energy and energy dissipation for each toroidal mode number n for the initial central energetic-ion beta β h0 = 1.7%. For this case, the TAE saturation level in the nonlinear MHD simulation is reduced by nearly half to v r /v A , while it is v r /v A for the linear-mhd simulation. The MHD fluctuation energy is given by ( ) 1 E = 2 ρv2 + B2 Beq 2 + δp dv. (13) 2µ 0 γ 1 We decompose the MHD fluctuations for each toroidal mode number n 0. For example, the velocity fluctuation is decomposed into v n (R,ϕ,z) = v nc (R,z) cos (nϕ) + v ns (R,z) sin (nϕ), (14) Figure 5. Comparison of radial velocity evolution for cosine part of m/n = 6/4 harmonics at r/a = 0.42 for β h0 = 2.0% between the linear-mhd and the nonlinear MHD runs (a) with only n = 0 and 4 modes retained and (b) without n = 0, but with n = 4, 8, and higher-n modes retained. of v r /v A for both the linear-mhd and nonlinear MHD runs are and in the two calculations the phases of the excitation remain synchronized in time, indicating in this case a robust accuracy of the reduced model. On the other hand, for β h0 = 2.0%, a significant reduction in the saturation level can be seen for the nonlinear MHD run in figure 4. The saturation level of v r /v A is for the nonlinear MHD simulation while it is for the linear-mhd simulation. The MHD nonlinear effects reduce the TAE saturation level by half when it reaches v r /v A δb r /B For cases where the instability growth is lower, the saturation level is v r /v A δb r /B 10 3 and the MHD nonlinearity does not play any important role. Then the saturation mechanism is dominated by the particle nonlinear dynamics, i.e. the particle trapping by the TAE causes the saturation. In the evolution of the TAE mode, the nonlinear terms in the MHD equations generate the fluctuations with toroidal mode numbers multiples of 4 (n = 0, 4, 8, 12, 16,...). It is interesting to investigate which toroidal mode number is important for the TAE saturation level reduction. We conducted two types of nonlinear MHD simulations. In the first type, only the n = 0 and 4 modes are retained while the higher-n (n 8) modes are artificially removed. In the second type, only the n = 4 and higher-n modes are retained without n = 0 fluctuations. The results are compared in figure 5 with that of the linear-mhd simulation. We can see that both the n = 0 and higher-n modes reduce the TAE saturation level. v nc (R,z) = 4 π v ns (R,z) = 4 π π/2 0 π/2 0 π/2 v (R,ϕ,z) cos (nϕ) dϕ v (R,ϕ,z) sin (nϕ) dϕ (for n 0) (15) v nc (R,z) = 2 π 0 v (R,ϕ,z) dϕ, v ns (R,z) = 0 (for n = 0). (16) The MHD kinetic energy is expressed in terms of the harmonics of density and velocity ρv2 dv = 2 ρ kv l v m dv l+m=k or l m =k (k, l, m 0). (17) We neglect the terms with k 0 because ρ 0 = ρ eq + δρ 0 δρ k 0 holds. Then, the MHD fluctuation energy is written as the sum of energy over the toroidal mode numbers E = E n, (18) n ( 1 E n 2 ρ 0vn 2 + 2δB n B eq + δbn 2 + δp ) n dv. (19) 2µ 0 γ 1 For n 0, E n is reduced to ( 1 E n 2 ρ 0vn ) δbn 2 dv, (20) 2µ 0 because the terms δb n B eq and δp n vanish after integration over the toroidal angle. We also analyse the evolution of the energy dissipation for each toroidal mode number n [νρ0 D n = ωn νρ 0 ( v n ) 2 ] + ηδj n j n dv. (21) 4
6 Figure 7. Evolution of damping rate for each toroidal mode number and total damping rate in the nonlinear MHD simulation, and damping rate in the linear-mhd simulation (γ d lin ) for β h0 = 1.7%. mode number is given by Figure 6. Evolution of (a) energy and (b) energy dissipation with toroidal mode number n = 4 denoted as E 4 and D 4 along with the total energy and total energy dissipation of toroidal mode numbers n = 0, 4, 8, 12, 16 denoted as E ALL and D ALL for β h0 = 1.7%. For n = 0, we investigate an alternative form of the fluctuation energy, S and E 0, defined by d dt S D ALL n E 0 E 0 S. D n, S t=0 = 0 (22) For the purpose of the saturation analysis, E 0 is not a useful variable because E 0 monotonically increases in time, as all the dissipated energy is transferred to the n = 0 mode. We have found that E 0 is a useful variable to understand the saturation process of the instability, as is discussed below. The evolution of energy with toroidal mode number n = 4 (E 4 ) and the sum of energy over toroidal mode numbers n = 0, 4, 8, 12, 16 (E ALL = E 0 + E 4 + E 8 + E 12 + E 16 ) is shown in figure 6(a). At saturation time, ω A t = 380, of the TAE instability, the energy in the toroidal mode numbers n = 12 (E 12 ) and 16 (E 16 ) are, respectively, 3% and 1% of the energy of the n = 4(E 4 ) component. The energy of the modes with n 20 is negligibly small and the energy and dissipation were not analysed, although such higher-n modes are included in the simulation. We also see in figure 6(a) that at saturation the total energy E ALL is greater than the energy of the n = 4 mode E 4 by 40%. The total energy dissipation of toroidal mode numbers n = 0, 4, 8, 12, 16 (D ALL = D 0 + D 4 + D 8 + D 12 + D 16 ) is compared with that of n = 4(D 4 ) in figure 6(b). The total dissipation is two times greater than the dissipation in n = 4 alone. Let us consider why the saturation level is reduced by the MHD nonlinearity. The energy evolution for each toroidal d dt E 4 = 2γ h (t) E 4 D 4 + P 4, d dt E 0 = D 0 + P 0, d dt E n = D n + P n (n 0, 4), (23) where γ h (t) represents the energetic-particle drive and P n is the energy transfer from the other toroidal mode numbers. The factor of 2 arises for γ h (t), which denotes the amplitude growth rate. The energetic-particle drive is restricted only to the n = 4 mode in this work. When the nonlinear coupling is so weak that P n can be neglected, the saturation condition for the n = 4 mode is given by γ h (t) = D 4 /2E 4. As the TAE amplitude grows, the particle trapping reduces γ h (t) to γ d4 = D 4 /2E 4 that is the damping rate of the n = 4 mode, leading to the saturation of the instability. For the case investigated here, the nonlinear coupling cannot be neglected. The summation of all equations (23) gives d dt E ALL = d E 0 dt + E n = 2γ h (t) E 4 D ALL. n 0 (24) Here, we used the relation n P n = 0 for energy conservation. As we see in figure 6(a), the saturation of E 4 and E ALL takes place roughly at the same time. Then, we investigate the saturation condition of E ALL instead of E 4. The saturation condition of E ALL is given by γ h (t) = D ALL /2E 4, which is a higher level of the energetic-particle drive than that for the linear-mhd run, corresponding to the earlier saturation of the instability. We can say that D ALL /2E 4 is an effective damping rate for the cases where nonlinear coupling between the different toroidal mode numbers is essential. In figure 7 the evolution of γ d n = D n /2E n and γ d ALL = D ALL /2E 4 are compared with γ d lin = D 4 /2E 4 in the relevant linear- MHD simulation. The total damping rate γ d ALL, in the nonlinear MHD simulation, is clearly greater than the n = 4 TAE damping rate γ d lin in the linear-mhd simulation. This explains why the saturation level is reduced by the MHD nonlinearity. The nonlinear coupling increases the total energy dissipation leading to the lower saturation level. 5
7 3.3. Effects of weak dissipation As we mentioned in section 2, the viscosity, diffusion and resistivity coefficients employed in the previous subsections and in the next subsection are chosen to be ν = ν n = 10 6 v A R 0 and η = 10 6 µ 0 v A R 0. Of course, in realistic fusion experiment plasmas the dissipation coefficients are lower by orders of magnitude. For example, normalized resistivity is η 10 9 µ 0 v A R 0 for hydrogen plasma with electron temperature 1 kev, ion density m 3, magnetic field 1 T and major radius 3 m. Thus, it is important to examine the effects of weak dissipation. On the other hand, simulation with lower dissipation coefficients needs a larger number of grid points. As the dissipation terms are second-order spatial derivatives, dissipation coefficients that are three orders of magnitude lower need times larger number of numerical grid points in one direction. This leads to times larger number of total grid points. Since such a simulation is not possible at the present time, we investigate the effects of weak dissipation by simulating a case with diffusion coefficients reduced by 1/4 to ν = ν n = v A R 0 and η = µ 0 v A R 0. The number of grid points for each of the cylindrical coordinates is doubled leading to grid points for (R,ϕ,z). The number of computational particles is increased by a factor of8to The radial velocity evolution is compared in figure 8(a) for the linear-mhd and the nonlinear MHD runs for β h0 = 1.7%. The saturation level of v r /v A is for the nonlinear MHD simulation while it is for the linear-mhd simulation. We thereby demonstrate again the reduction in the saturation amplitude by the nonlinear MHD effects but now with weaker dissipation. For the linear-mhd run, another unstable TAE with major harmonics m/n = 6/4 and 7/4 affects the evolution of the TAE focused on in this paper, leading to the second amplitude peak higher than the first peak. Shown in figure 8(b) is both the evolution of the damping rate for each toroidal mode number together with the total damping rate in the nonlinear MHD simulation and the damping rate in the linear-mhd simulation. Comparing the results for the lower with the higher dissipation case shown in figure 7, we find that in the lower dissipation case we have a lower damping rate for the linearly growing phase prior to ω A t = 250. However, the total damping rate γ d ALL increases in the nonlinear phase, just in the same way as the stronger dissipation case illustrated in figure 7. As seen from figure 8(b), the damping rate for the higher-n modes increases significantly in the nonlinear phase leading to the high level of the total damping rate. This result suggests that the saturation level would similarly reduce by the nonlinear MHD effects even if the parameters for more realistic dissipation were taken into account. We would also like to point out that a mechanism for substantial dissipation of higher-n modes might be arising from the continuum damping that does not depend on the dissipation coefficients [20, 21] yet exists within an MHD model. Further work to verify this conjecture is needed Spatial profile and evolution of nonlinearly generated modes Let us examine the spatial profiles and the evolution of the nonlinearly generated modes. We have analysed the spatial Figure 8. Comparison of the radial velocity evolution for the linear-mhd and the nonlinear MHD runs using the cosine part of m/n = 6/4 harmonics at r/a = 0.42 (a), and evolution of damping rate for each toroidal mode number and total damping rate in the nonlinear MHD simulation, and damping rate in the linear-mhd simulation (γ d lin )(b) for β h0 = 1.7% and ν = η/µ 0 = ν n = v A R 0. profiles of the n = 0 poloidal flow for different times of the nonlinear MHD simulation with β h0 = 2.0%. The radial velocity evolution for the n = 4 mode was already shown in figure 4. The spatial profiles of the n = 0 poloidal flow at ω A t = 189 and ω A t = 227 are shown in figure 9. When we compare figures 9(a) and (b), we see that the n = 0 poloidal flow profile remains constant during the linear phase of the TAE instability. The profile peaks at the maximum position of the TAE mode amplitude. The largest harmonics are the m/n = 0/0 and 1/0 cosine ( cos ϑ) harmonics. The evolution of the 0/0 harmonic of the poloidal flow, namely the zonal flow, is shown in figure 10(a) along with the 1/0 sine harmonic of the pressure fluctuation. The amplitude evolution of the zonal flow and the principal harmonic of the TAE radial velocity at the same magnetic surface are shown in figure 10(b). As can be seen in figure 10(a), the zonal flow and the pressure oscillate with the same frequency after the saturation of the instability. This oscillation is a GAM because the frequency of the oscillation ω = 0.11ω A is close to the theoretical GAM frequency ω GAM = ω A γβbulk (1+1/2q 2 ) = 0.12ω A [22, 23]. Furthermore, the coupling between the zonal flow and the 1/0 sine harmonic of the pressure fluctuation is consistent with the theory of GAM. We see in figure 10(b) that the zonal flow amplitude grows at twice the TAE growth rate. This suggests that the n = 0 modes are primarily generated by the coupling of two n = 4 harmonics of the TAE mode. Let us discuss the n = 0 mode in its exponentially growing phase and how the GAM is excited. When the n = 0 mode amplitude is sufficiently small and the source generated by the nonlinear coupling is neglected, the n = 0 mode evolution can be described by the linear MHD equations (8) (11). This can 6
8 Figure 11. Spatial profile of each poloidal harmonic of the radial velocity with n = 8 during the linearly growing phase of the TAE instability. Solid (dashed) lines show cos(mϑ + nϕ) [sin(mϑ + nϕ)] harmonics with poloidal mode number m labelled in the figure. where z is the n = 0 mode fluctuations z = t (δρ, δv,δb,δp) and M eq is the linear-mhd operator which is a function of the equilibrium variables. Next we consider the source vector s which is generated by the nonlinear coupling of the TAE eigenfunction. The evolution of z at the linear phase is described by Figure 9. Spatial profile of each poloidal harmonic of the poloidal flow with n = 0 during the linearly growing phase of the TAE instability at (a) ω A t = 189 and (b) ω A t = 227 for β h0 = 2.0%. Solid (dashed) lines show cos(mϑ) [sin(mϑ)] harmonics with poloidal mode number m labelled in the figure. z + M eq (z) = s. (26) The source vector can be expressed as s = s 0 e 2γTAEt using the TAE growth rate γ TAE because the nonlinear coupling of the TAE eigenvectors generates the n = 0 fluctuations with zero frequency and growth rate 2γ TAE. We also see from the simulation results shown in figures 9 and 10 that the n = 0 mode fluctuations can be expressed as z = z 0 e 2γTAEt. Here, s 0 and z 0 are independent of time. Equation (26) is reduced to 2γ TAE I (z 0 ) + M eq (z 0 ) = [ 2γ TAE I + M eq ] (z0 ) = s 0, (27) Figure 10. Time evolution of (a) the zonal flow and pressure fluctuation and (b) the zonal flow amplitude and the TAE amplitude for β h0 = 2.0%. be expressed in the following form z + M eq (z) = 0, (25) where I is the unit operator and both I and M eq are linear operators of z 0. Then, the spatial profile of the n = 0 mode fluctuations is given by z 0 = (2γ TAE I + M eq ) 1 (s 0 ). The zonal flow profile shown in figure 9 is a part of the solution z 0 that contains the n = 0 fluctuations of all the MHD variables. It should be emphasized that the spatial profile z 0 depends on 2γ TAE, the growth rate of the nonlinear source during the TAE linear phase. When the instability saturates, the growth rate of the nonlinear source approaches zero. Then, the n = 0 field z 0 that is matched to the TAE linear growth is no longer a solution of equation (26) during the nonlinear phase. This leads to the excitation of the n = 0 MHD waves that include the GAM. The spatial profiles of the n = 8 harmonics have also been investigated. The radial velocity profile is shown in figure 11. The profile peaks at the TAE peak location. The largest harmonics are m/n = 10/8 and 11/8. The m/n = 11/8 harmonic is consistent with the nonlinear source generated by the coupling of the TAE harmonics m/n = 5/4 and 6/4. The m/n = 10/8 harmonic may be generated from the nonlinear source with m/n = 11/8 through toroidal coupling. The evolution of the 11/8 harmonic and the frequency is shown in figure 12. We see in figure 12(b) that the real frequency of the n = 8 mode is twice the TAE frequency. 7
9 Figure 12. Time evolution of (a) the cosine part of m/n = 11/8 radial velocity harmonics and (b) the frequency for β h0 = 2.0%. 4. Discussion and summary The nonlinear MHD effects on the evolution of the Alfvén eigenmode were investigated with hybrid simulations of an MHD fluid interacting with energetic particles. To clarify the role of the MHD nonlinearity, the nonlinear MHD results were compared with results from a reduced model, where only linear-mhd equations were solved together with a nonlinear response of the energetic particles. Specifically, we studied the evolution of an n = 4 TAE mode destabilized by its resonant interaction with energetic particles in a tokamak plasma. When the TAE saturation level is δb/b 10 3 no significant difference was found between the results of the linear-mhd simulation and the nonlinear MHD simulations. On the other hand, when in the linear-mhd simulation the TAE saturation level is δb/b 10 2, the saturation level in the nonlinear MHD case is found to be reduced to half the result of the linear- MHD simulation. We found that the nonlinearly generated n = 0 and the higher-n (n 8) modes provide increased energy dissipation that appears crucial for achieving a reduced TAE saturation level. The effects of MHD nonlinearities on the increase in TAE damping rate were discussed in previous works [11, 12]. We found in this work that the total damping rate of all the toroidal mode numbers increases before the saturation of the instability. The increased dissipation leads to the reduction in the saturation amplitude even though the damping rate of the n = 4 component, which was the original TAE carrier, decreases slightly before saturation. We emphasize that this is a new picture for the mechanism by which nonlinear MHD effects cause saturation that is different from the previous works. In the simulation results of the TAE bursts presented in [1], the energetic-particle loss takes place when the TAE amplitude reaches δb/b while the saturation amplitude is δb/b It was also demonstrated in another study [10] that synchronized bursts of multiple TAEs take place when the MHD nonlinearities are taken into account. In contrast, in this work, we have focused on the evolution of a single TAE and demonstrated that the nonlinear MHD effects reduce the TAE saturation level. Computer simulations of multiple Alfvén eigenmodes have been conducted [6, 10, 12, 24]. There is still a need to investigate the nonlinear MHD effects in multiple mode evolution by comparing the linear and nonlinear MHD simulations as has been performed in this work. For future work, it is important to reproduce the TAE bursts, using realistic parameters. The challenge is to demonstrate saturation amplitudes that are closer to the value inferred from the experimental plasma displacement measurements. Another interesting discovery of this work is the excitation of the GAM after the saturation of the TAE instability. In the linearly growing phase of the instability, the zonal fluctuations are matched to the growth of the nonlinear source, with a spatial profile kept constant. When the instability saturates, however, the spatial profile of the zonal fluctuations is no longer matched to the evolution of the nonlinear source. This excites the n = 0 MHD waves that includes the GAM. It is evident that the GAM is excited through the nonlinear MHD effect because the n = 4 fluctuations of the energetic-particle current density retained in the simulations cannot directly drive the GAM. This excitation mechanism of the GAM is different from the direct destabilization by the energetic particles [25, 26]. We note that the excitation of zonal flow and GAM through the TAE nonlinearity is an interesting phenomenon for burning plasmas since its excitation may lead to improved plasma confinement and to bulk plasma heating through the damping of the GAM. Acknowledgments Numerical computations were performed at the Plasma Simulator (HITACHI SR16000) of the National Institute for Fusion Science. This work was supported by a Grant-in- Aid for Scientific Research from the Japan Society for the Promotion of Science (No ). References [1] Todo Y., Berk H.L. and Breizman B.N Phys. Plasmas [2] Wong K.L. et al 1991 Phys. Rev. Lett [3] Durst R.D., Fonck R.J., Wong K.L., Cheng C.Z., Fredrickson E.D. and Paul S.F Phys. Fluids B [4] Berk H.L. and Breizman B.B Phys. Fluids B [5] Berk H.L., Breizman B.N. and Ye H Phys. Fluids B [6] Fu G.Y. and Park W Phys. Rev. Lett [7] Todo Y., Sato T., Watanabe K., Watanabe T.H. and Horiuchi R Phys. Plasmas [8] Wu Y., White R.B., Chen Y. and Rosenbluth M.N Phys. Plasmas [9] Briguglio S., Zonca F. and Vlad G Phys. Plasmas [10] Todo Y., Watanabe T.-H., Hyoung-Bin Park and Sato T Nucl. Fusion [11] Zonca F., Romanelli F., Vlad G. and Kar C Phys. Rev. Lett [12] Spong D.A., Carreras B.A. and Hedrick C.L Phys. Plasmas
10 [13] Todo Y. and Sato T Phys. Plasmas [14] Todo Y., Shinohara K., Takechi M. and Ishikawa M Phys. Plasmas [15] Todo Y Phys. Plasmas [16] Park W. et al 1992 Phys. Fluids B [17] Spong D.A., Carreras B.A. and Hedrick C.L Phys. Fluids B [18] Briguglio S., Vlad G., Zonca F. and Kar C Phys. Plasmas [19] Littlejohn R.G Phys. Fluids [20] Zonca F. and Chen L Phys. Rev. Lett [21] Rosenbluth M.N., Berk H.L., Van Dam J.W. and Lindberg D.M Phys. Rev. Lett [22] Winsor N., Johnson J.L. and Dawson J.M Phys. Fluids [23] Hassam A.B. and Drake J.F Phys. Fluids B [24] Vlad G. et al 2009 Nucl. Fusion [25] Berk H.L. et al 2006 Nucl. Fusion 46 S888 [26] Fu G.Y Phys. Rev. Lett
Nonlinear MHD effects on TAE evolution and TAE bursts
Nonlinear MHD effects on TAE evolution and TAE bursts Y. Todo (NIFS) collaborating with H. L. Berk and B. N. Breizman (IFS, Univ. Texas) GSEP 3rd Annual Meeting (remote participation / Aug. 9-10, 2010)
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 informationNonperturbative Effects of Energetic Ions on Alfvén Eigenmodes
1 TH/3-1Ra Nonperturbative Effects of Energetic Ions on Alfvén Eigenmodes Y. Todo 1), N. Nakajima 1), K. Shinohara 2), M. Takechi 2), M. Ishikawa 2), S. Yamamoto 3) 1) National Institute for Fusion Science,
More informationMHD-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 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 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 informationNonlinear Evolution and Radial Propagation of the Energetic Particle Driven GAM
Nonlinear Evolution and Radial Propagation of the Energetic Particle Driven GAM by R. Nazikian In collaboration with G.Y. Fu, R.V. Budny, G.J. Kramer, PPPL G.R. McKee, U. Wisconsin T. Rhodes, L. Schmidt,
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 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 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 informationStabilization of sawteeth in tokamaks with toroidal flows
PHYSICS OF PLASMAS VOLUME 9, NUMBER 7 JULY 2002 Stabilization of sawteeth in tokamaks with toroidal flows Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College
More informationENERGETIC PARTICLES AND BURNING PLASMA PHYSICS
ENERGETIC PARTICLES AND BURNING PLASMA PHYSICS Reported by J. Van Dam Institute for Fusion Studies The University of Texas at Austin US-Japan JIFT Workshop on Theory-Based Modeling and Integrated Simulation
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 informationCharacteristics of Energetic-Ion-Driven Geodesic Acoustic Modes in the Large Helical Device(LHD)
O-4 12 th IAEA TM on Energetic Particles in Magnetic Confinement Systems, 7-10 Sep, Austin, USA Characteristics of Energetic-Ion-Driven Geodesic Acoustic Modes in the Large Helical Device(LHD) K. Toi,
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 informationPresentation by Herb Berk University of Texas at Austin Institute for Fusion Studies in Vienna, Austria Sept. 1-4, 2015
Review of Theory Papers at 14 th IAEA technical meeting on Engertic Particles in Magnetic Confinement systems Presentation by Herb Berk University of Texas at Austin Institute for Fusion Studies in Vienna,
More informationMHD Simulation of High Wavenumber Ballooning-like Modes in LHD
1 TH/P9-16 MHD Simulation of High Wavenumber Ballooning-like Modes in LHD H. Miura and N. Nakajima National Institute for Fusion Science, 322-6 Oroshi, Toki, Gifu 509-5292, JAPAN e-mail contact of main
More informationPlasmoid Motion in Helical Plasmas
Plasmoid Motion in Helical Plasmas Ryuichi ISHIZAKI and Noriyoshi NAKAJIMA National Institute for Fusion Science, Toki 509-5292, Japan (Received 12 December 2009 / Accepted 18 May 2010) In order to explain
More informationNonlinear Consequences of Weakly Driven Energetic Particle Instabilities
2008 International Sherwood Fusion Theory Conference March 30 - April 2, 2008, Boulder, Colorado Nonlinear Consequences of Weakly Driven Energetic Particle Instabilities Boris Breizman Institute for Fusion
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 informationGeodesic Acoustic and related modes
Geodesic Acoustic and related modes A. Smolyakov* Acknowledgements X. Garbet, C. Nguyen (CEA Cadarache) V.I. Ilgisonis, V.P. Lakhin, A.Melnikov (Kurchatov Institute) * University of Saskatchewan, Canada
More informationMHD instabilities and fast particles
ENEA F. Zonca 1 MHD instabilities and fast particles Fulvio Zonca Associazione Euratom-ENEA sulla Fusione, C.R. Frascati, C.P. 65-44 - Frascati, Italy. July 11.th, 25 : Turbulence overshoot and resonant
More informationToroidal flow stablization of disruptive high tokamaks
PHYSICS OF PLASMAS VOLUME 9, NUMBER 6 JUNE 2002 Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College Park, Maryland 20742-3511 Received 4 February 2002; accepted
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 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 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 informationMHD instability driven by supra-thermal electrons in TJ-II stellarator
MHD instability driven by supra-thermal electrons in TJ-II stellarator K. Nagaoka 1, S. Yamamoto 2, S. Ohshima 2, E. Ascasíbar 3, R. Jiménez-Gómez 3, C. Hidalgo 3, M.A. Pedrosa 3, M. Ochando 3, A.V. Melnikov
More informationConfiguration Optimization of a Planar-Axis Stellarator with a Reduced Shafranov Shift )
Configuration Optimization of a Planar-Axis Stellarator with a Reduced Shafranov Shift ) Shoichi OKAMURA 1,2) 1) National Institute for Fusion Science, Toki 509-5292, Japan 2) Department of Fusion Science,
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 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 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 informationMHD Linear Stability Analysis Using a Full Wave Code
US-Japan JIFT Workshop on Progress of Extended MHD Models NIFS, Toki,Japan 2007/03/27 MHD Linear Stability Analysis Using a Full Wave Code T. Akutsu and A. Fukuyama Department of Nuclear Engineering, Kyoto
More informationBursty Transport in Tokamaks with Internal Transport Barriers
Bursty Transport in Tokamaks with Internal Transport Barriers S. Benkadda 1), O. Agullo 1), P. Beyer 1), N. Bian 1), P. H. Diamond 3), C. Figarella 1), X. Garbet 2), P. Ghendrih 2), V. Grandgirard 1),
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 informationSimulation Study of Energetic Particle Driven Alfvén Eigenmodes and Geodesic Acoustic Modes in Toroidal Plasmas
Simulation Study of Energetic Particle Driven Alfvén Eigenmodes and Geodesic Acoustic Modes in Toroidal Plasmas WANG, Hao DOCTOR OF PHILOSOPHY Department of Fusion Science School of Physical Sciences The
More informationHigh-m Multiple Tearing Modes in Tokamaks: MHD Turbulence Generation, Interaction with the Internal Kink and Sheared Flows
TH/P3-3 High-m Multiple Tearing Modes in Tokamaks: MHD Turbulence Generation, Interaction with the Internal Kink and Sheared Flows A. Bierwage 1), S. Benkadda 2), M. Wakatani 1), S. Hamaguchi 3), Q. Yu
More informationThe Linear Theory of Tearing Modes in periodic, cyindrical plasmas. Cary Forest University of Wisconsin
The Linear Theory of Tearing Modes in periodic, cyindrical plasmas Cary Forest University of Wisconsin 1 Resistive MHD E + v B = ηj (no energy principle) Role of resistivity No frozen flux, B can tear
More informationGTC Simulation of Turbulence and Transport in Tokamak Plasmas
GTC Simulation of Turbulence and Transport in Tokamak Plasmas Z. Lin University it of California, i Irvine, CA 92697, USA and GPS-TTBP Team Supported by SciDAC GPS-TTBP, GSEP & CPES Motivation First-principles
More informationHybrid Kinetic-MHD simulations with NIMROD
simulations with NIMROD 1 Yasushi Todo 2, Dylan P. Brennan 3, Kwang-Il You 4, Jae-Chun Seol 4 and the NIMROD Team 1 University of Washington, Seattle 2 NIFS, Toki-Japan 3 University of Tulsa 4 NFRI, Daejeon-Korea
More informationGlobal gyrokinetic particle simulations with kinetic electrons
IOP PUBLISHING Plasma Phys. Control. Fusion 49 (2007) B163 B172 PLASMA PHYSICS AND CONTROLLED FUSION doi:10.1088/0741-3335/49/12b/s15 Global gyrokinetic particle simulations with kinetic electrons Z Lin,
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 informationTheory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus. K. C. Shaing
Theory for Neoclassical Toroidal Plasma Viscosity in a Toroidally Symmetric Torus K. C. Shaing Plasma and Space Science Center, and ISAPS, National Cheng Kung University, Tainan, Taiwan 70101, Republic
More informationTwo Fluid Dynamo and Edge-Resonant m=0 Tearing Instability in Reversed Field Pinch
1 Two Fluid Dynamo and Edge-Resonant m= Tearing Instability in Reversed Field Pinch V.V. Mirnov 1), C.C.Hegna 1), S.C. Prager 1), C.R.Sovinec 1), and H.Tian 1) 1) The University of Wisconsin-Madison, Madison,
More informationSimple examples of MHD equilibria
Department of Physics Seminar. grade: Nuclear engineering Simple examples of MHD equilibria Author: Ingrid Vavtar Mentor: prof. ddr. Tomaž Gyergyek Ljubljana, 017 Summary: In this seminar paper I will
More informationShear Alfvén Wave Continuous Spectrum in the Presence of a Magnetic Island
1 TH/P3-5 Shear Alfvén Wave Continuous Spectrum in the Presence of a Magnetic Island A. Biancalani 1), L. Chen 2) 3), F. Pegoraro 1), F. Zonca 4), S. V. Annibaldi 5), A. Botrugno 4), P. Buratti 4) and
More informationFinite-Orbit-Width Effect and the Radial Electric Field in Neoclassical Transport Phenomena
1 TH/P2-18 Finite-Orbit-Width Effect and the Radial Electric Field in Neoclassical Transport Phenomena S. Satake 1), M. Okamoto 1), N. Nakajima 1), H. Sugama 1), M. Yokoyama 1), and C. D. Beidler 2) 1)
More informationCollisionless nonideal ballooning modes
PHYSICS OF PLASMAS VOLUME 6, NUMBER 1 JANUARY 1999 Collisionless nonideal ballooning modes Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College Park, Maryland
More informationHybrid Kinetic-MHD simulations Status and Updates
in NIMROD simulations Status and Updates Charlson C. Kim 1,2 Yasushi Todo 2 and the NIMROD Team 1. University of Washington, Seattle 2. National Institute for Fusion Science NIMROD Team Meeting Austin,
More informationActive and Fast Particle Driven Alfvén Eigenmodes in Alcator C-Mod
Active and Fast Particle Driven Alfvén Eigenmodes in Alcator C-Mod JUST DID IT. J A Snipes, N Basse, C Boswell, E Edlund, A Fasoli #, N N Gorelenkov, R S Granetz, L Lin, Y Lin, R Parker, M Porkolab, J
More informationW.A. HOULBERG Oak Ridge National Lab., Oak Ridge, TN USA. M.C. ZARNSTORFF Princeton Plasma Plasma Physics Lab., Princeton, NJ USA
INTRINSICALLY STEADY STATE TOKAMAKS K.C. SHAING, A.Y. AYDEMIR, R.D. HAZELTINE Institute for Fusion Studies, The University of Texas at Austin, Austin TX 78712 USA W.A. HOULBERG Oak Ridge National Lab.,
More informationNumerical calculation of the Hamada basis vectors for three-dimensional toroidal magnetic configurations
PHYSICS OF PLASMAS 12, 072513 2005 Numerical calculation of the Hamada basis vectors for three-dimensional toroidal magnetic configurations J. N. Talmadge and S. P. Gerhardt a HSX Plasma Laboratory, University
More informationPredictions of fusion α-particle transport due to Alfvén eigenmodes in ITER
Predictions of fusion α-particle transport due to Alfvén eigenmodes in ITER M. Fitzgerald, S.E. Sharapov, P. Rodrigues 2, A. Polevoi 3, D. Borba 2 2 Instituto de Plasmas e Fusão Nuclear, Instituto Superior
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 information(a) (b) (c) (d) (e) (f) r (minor radius) time. time. Soft X-ray. T_e contours (ECE) r (minor radius) time time
Studies of Spherical Tori, Stellarators and Anisotropic Pressure with M3D 1 L.E. Sugiyama 1), W. Park 2), H.R. Strauss 3), S.R. Hudson 2), D. Stutman 4), X-Z. Tang 2) 1) Massachusetts Institute of Technology,
More informationFundamentals of Magnetic Island Theory in Tokamaks
Fundamentals of Magnetic Island Theory in Tokamaks Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX, USA Talk available at http://farside.ph.utexas.edu/talks/talks.html
More informationMeasuring from electron temperature fluctuations in the Tokamak Fusion Test Reactor
PHYSICS OF PLASMAS VOLUME 5, NUMBER FEBRUARY 1998 Measuring from electron temperature fluctuations in the Tokamak Fusion Test Reactor C. Ren, a) J. D. Callen, T. A. Gianakon, and C. C. Hegna University
More informationStudies on Neutral Beam Ion Confinement and MHD Induced Fast-Ion. Loss on HL-2A Tokamak
Studies on Neutral Beam Ion Confinement and MHD Induced Fast-Ion Loss on HL-A Tokamak LIU Yi, ISOBE Mitsutaka, PENG Xiao-Dong, Wang Hao, JI Xiao-Quan, CHEN Wei, ZHANG Yi-Po, Dong Yun-Bo, MORITA Shigeru
More informationSMR/ Summer College on Plasma Physics. 30 July - 24 August, Introduction to Magnetic Island Theory.
SMR/1856-1 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 Introduction to Magnetic Island Theory. R. Fitzpatrick Inst. for Fusion Studies University of Texas at Austin USA Introduction
More informationParticle Pinch Model of Passing/Trapped High-Z Impurity with Centrifugal Force Effect )
Particle Pinch Model of Passing/Trapped High-Z Impurity with Centrifugal Force Effect ) Yusuke SHIMIZU, Takaaki FUJITA, Atsushi OKAMOTO, Hideki ARIMOTO, Nobuhiko HAYASHI 1), Kazuo HOSHINO 2), Tomohide
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 informationParticle-in-cell simulations of electron transport from plasma turbulence: recent progress in gyrokinetic particle simulations of turbulent plasmas
Institute of Physics Publishing Journal of Physics: Conference Series 16 (25 16 24 doi:1.188/1742-6596/16/1/2 SciDAC 25 Particle-in-cell simulations of electron transport from plasma turbulence: recent
More informationThe Field-Reversed Configuration (FRC) is a high-beta compact toroidal in which the external field is reversed on axis by azimuthal plasma The FRC is
and Stability of Field-Reversed Equilibrium with Toroidal Field Configurations Atomics General Box 85608, San Diego, California 92186-5608 P.O. APS Annual APS Meeting of the Division of Plasma Physics
More informationAC loop voltages and MHD stability in RFP plasmas
AC loop voltages and MHD stability in RFP plasmas K. J. McCollam, D. J. Holly, V. V. Mirnov, J. S. Sar, D. R. Stone UW-Madison 54rd Annual Meeting of the APS-DPP October 29th - November 2nd, 2012 Providence,
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 informationEffect of magnetic reconnection on CT penetration into magnetized plasmas
Earth Planets Space, 53, 547 55, 200 Effect of magnetic reconnection on CT penetration into magnetized plasmas Yoshio Suzuki, Takaya Hayashi 2, and Yasuaki Kishimoto Naka Fusion esearch Establishment,
More informationOscillating-Field Current-Drive Experiment on MST
Oscillating-Field Current-Drive Experiment on MST K. J. McCollam, J. K. Anderson, D. J. Den Hartog, F. Ebrahimi, J. A. Reusch, J. S. Sarff, H. D. Stephens, D. R. Stone University of Wisconsin-Madison D.
More informationFine-Scale Zonal Flow Suppression of Electron Temperature Gradient Turbulence
Fine-Scale Zonal Flow Suppression of Electron Temperature Gradient Turbulence S.E. Parker, J.J. Kohut, Y. Chen, Z. Lin, F.L. Hinton and W.W. Lee Center for Integrated Plasma Studies, University of Colorado,
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 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 informationMomentum transport from magnetic reconnection in laboratory an. plasmas. Fatima Ebrahimi
Momentum transport from magnetic reconnection in laboratory and astrophysical plasmas Space Science Center - University of New Hampshire collaborators : V. Mirnov, S. Prager, D. Schnack, C. Sovinec Center
More informationImpact of Localized ECRH on NBI and ICRH Driven Alfven Eigenmodes in the ASDEX Upgrade Tokamak
Impact of Localized ECRH on NBI and ICRH Driven Alfven Eigenmodes in the ASDEX Upgrade Tokamak M. Garcia-Munoz M. A. Van Zeeland, S. Sharapov, Ph. Lauber, J. Ayllon, I. Classen, G. Conway, J. Ferreira,
More informationMHD WAVES AND GLOBAL ALFVÉN EIGENMODES
MHD WVES ND GLOBL LFVÉN EIGENMODES S.E. Sharapov Euratom/CCFE Fusion ssociation, Culham Science Centre, bingdon, Oxfordshire OX14 3DB, UK S.E.Sharapov, Lecture 3, ustralian National University, Canberra,
More informationNonlinear hydrid simulations of precessional Fishbone instability
Nonlinear hydrid simulations of precessional Fishbone instability M. Faganello 1, M. Idouakass 1, H. L. Berk 2, X. Garbet 3, S. Benkadda 1 1 Aix-Marseille University, CNRS, PIIM UMR 7345, Marseille, France
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 informationControl of Neo-classical tearing mode (NTM) in advanced scenarios
FIRST CHENGDU THEORY FESTIVAL Control of Neo-classical tearing mode (NTM) in advanced scenarios Zheng-Xiong Wang Dalian University of Technology (DLUT) Dalian, China Chengdu, China, 28 Aug, 2018 Outline
More informationModelling of Frequency Sweeping with the HAGIS code
Modelling of Frequency Sweeping with the HAGIS code S.D.Pinches 1 H.L.Berk 2, S.E.Sharapov 3, M.Gryaznavich 3 1 Max-Planck-Institut für Plasmaphysik, EURATOM Assoziation, Garching, Germany 2 Institute
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationarxiv: v1 [physics.plasm-ph] 11 Mar 2016
1 Effect of magnetic perturbations on the 3D MHD self-organization of shaped tokamak plasmas arxiv:1603.03572v1 [physics.plasm-ph] 11 Mar 2016 D. Bonfiglio 1, S. Cappello 1, M. Veranda 1, L. Chacón 2 and
More informationINTERACTION OF DRIFT WAVE TURBULENCE AND MAGNETIC ISLANDS
INTERACTION OF DRIFT WAVE TURBULENCE AND MAGNETIC ISLANDS A. Ishizawa and N. Nakajima National Institute for Fusion Science F. L. Waelbroeck, R. Fitzpatrick, W. Horton Institute for Fusion Studies, University
More informationEffects of drag and diffusion on nonlinear behavior of EP-driven instabilities.
IAEA-TM EP 2011 / 09 / 07 Effects of drag and diffusion on nonlinear behavior of EP-driven instabilities. Maxime Lesur Y. Idomura, X. Garbet, P. Diamond, Y. Todo, K. Shinohara, F. Zonca, S. Pinches, M.
More informationMagnetohydrodynamic stability of negative central magnetic shear, high pressure ( pol 1) toroidal equilibria
Magnetohydrodynamic stability of negative central magnetic shear, high pressure ( pol 1) toroidal equilibria Robert G. Kleva Institute for Plasma Research, University of Maryland, College Park, Maryland
More informationTurbulent Transport due to Kinetic Ballooning Modes in High-Beta Toroidal Plasmas
1 TH/P-3 Turbulent Transport due to Kinetic allooning Modes in High-eta Toroidal Plasmas A. Ishizawa 1, S. Maeyama, T.-H. Watanabe 1, H. Sugama 1 and N. Nakajima 1 1 National Institute for Fusion Science,
More informationMagnetic Reconnection Controlled by Multi-Hierarchy Physics in an Open System
Magnetic Reconnection Controlled by Multi-Hierarchy Physics in an Open System Ritoku HORIUCHI 1,2), Shunsuke USAMI 1), Hiroaki OHTANI 1,2) and Toseo MORITAKA 3) 1) National Institute for Fusion Science,
More informationEX/4-6Rb Electrostatic fluctuation and fluctuation-induced particle flux during formation of the edge transport barrier in the JFT-2M tokamak
Electrostatic fluctuation and fluctuation-induced particle flux during formation of the edge transport barrier in the JFT-2M tokamak T. Ido 1), Y. Miura 2), K. Hoshino 2), Y. Hamada 1), Y. Nagashima 1),
More informationBunno, M.; Nakamura, Y.; Suzuki, Y. Matsunaga, G.; Tani, K. Citation Plasma Science and Technology (
Title The Finite Beta Effects on the Toro Tokamak Plasma Author(s) Bunno, M.; Nakamura, Y.; Suzuki, Y. Matsunaga, G.; Tani, K. Citation Plasma Science and Technology (2013 Issue Date 2013-02 URL http://hdl.handle.net/2433/173038
More informationCoarse-graining the electron distribution in turbulence simulations of tokamak plasmas
Coarse-graining the electron distribution in turbulence simulations of tokamak plasmas Yang Chen and Scott E. Parker University of Colorado at Boulder Gregory Rewoldt Princeton Plasma Physics Laboratory
More informationControl of Sawtooth Oscillation Dynamics using Externally Applied Stellarator Transform. Jeffrey Herfindal
Control of Sawtooth Oscillation Dynamics using Externally Applied Stellarator Transform Jeffrey Herfindal D.A. Ennis, J.D. Hanson, G.J. Hartwell, S.F. Knowlton, X. Ma, D.A. Maurer, M.D. Pandya, N.A. Roberds,
More informationIssues in Neoclassical Tearing Mode Theory
Issues in Neoclassical Tearing Mode Theory Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX Tearing Mode Stability in Tokamaks According to standard (single-fluid)
More informationSTABILIZATION OF m=2/n=1 TEARING MODES BY ELECTRON CYCLOTRON CURRENT DRIVE IN THE DIII D TOKAMAK
GA A24738 STABILIZATION OF m=2/n=1 TEARING MODES BY ELECTRON CYCLOTRON CURRENT DRIVE IN THE DIII D TOKAMAK by T.C. LUCE, C.C. PETTY, D.A. HUMPHREYS, R.J. LA HAYE, and R. PRATER JULY 24 DISCLAIMER This
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 Role of Dynamo Fluctuations in Anomalous Ion Heating, Mode Locking, and Flow Generation
The Role of Dynamo Fluctuations in Anomalous Ion Heating, Mode Locking, and Flow Generation P. W. Terry 1), R. Gatto 1), R. Fitzpatrick 2), C.C. Hegna 3), and G. Fiksel 1) 1) Department of Physics, University
More informationLecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit.
Lecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit. Steve Cowley UCLA. This lecture is meant to introduce the simplest ideas about kink modes. It would take many lectures to develop the
More informationSUMMARY OF EXPERIMENTAL CORE TURBULENCE CHARACTERISTICS IN OH AND ECRH T-10 TOKAMAK PLASMAS
SUMMARY OF EXPERIMENTAL CORE TURBULENCE CHARACTERISTICS IN OH AND ECRH T-1 TOKAMAK PLASMAS V. Vershkov, L.G. Eliseev, S.A. Grashin. A.V. Melnikov, D.A. Shelukhin, S.V. Soldatov, A.O. Urazbaev and T-1 team
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 informationImpact of Energetic-Ion-Driven Global Modes on Toroidal Plasma Confinements
Impact of Energetic-Ion-Driven Global Modes on Toroidal Plasma Confinements Kazuo TOI CHS & LHD Experimental Group National Institute for Fusion Science Toki 59-5292, Japan Special contributions from:
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 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 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 informationELECTROSTATIC ION-CYCLOTRON WAVES DRIVEN BY PARALLEL VELOCITY SHEAR
1 ELECTROSTATIC ION-CYCLOTRON WAVES DRIVEN BY PARALLEL VELOCITY SHEAR R. L. Merlino Department of Physics and Astronomy University of Iowa Iowa City, IA 52242 December 21, 2001 ABSTRACT Using a fluid treatment,
More information