Simulations of H-Mode Plasmas in Tokamak Using a Complete Core-Edge Modeling in the BALDUR Code

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "Simulations of H-Mode Plasmas in Tokamak Using a Complete Core-Edge Modeling in the BALDUR Code"

Transcription

1 Plasma Science and Technology, Vol.14, No.9, Sep Simulations of H-Mode Plasmas in Tokamak Using a Complete Core-Edge Modeling in the BALDUR Code Y. PIANROJ, T. ONJUN School of Manufacturing Systems and Mechanical Engineering, Sirindhorn International Institute of Technology, Thammasat University, Pathum Thani 12121, Thailand Abstract A theory-based model for predicting the pedestal formation in both ion and electron temperatures, and hydrogenic and impurity density is developed and implemented in the 1.5D BALDUR codes for self-consistently simulating H-mode plasma in tokamak. In the simulation, the transports around pedestal, including the electron and ion thermal, hydrogenic and impurity particle transports are calculated using an anomalous semi-empirical mixed Bohm/gyro-Bohm (Mixed B/gB) model, which is modified to include the effects of ω E B flow shear and magnetic shear. Because of the reduction of transport, the pedestal can be formed. For a preliminary test, this core-edge model is used to simulate the temporal evolution of plasma current, temperature, and density profiles for DIII-D discharges. It is found that the simulations successfully reproduce the experimental results. A statistical analysis, including RMSE and offset, is used to quantify the agreement between the prediction and the corresponding experimental results. The simulation results show an agreement with average RMSE of 11.87%, 14.53%, 7.59% and 12.21% for electron temperature, ion temperature, electron density, and deuterium density profiles, respectively. In addition, it is found that the suppression function developed is effective only in the edge region. Keywords: ITER, pedestal, H-mode, BALDUR, fusion performance, transport PACS: y, Fa, Fi DOI: / /14/9/02 1 Introduction A major advance in magnetic confinement fusion occurred with the discovery of a new operational regime, called the High confinement mode (H-mode) [1]. The H-mode operation results in a significant increase in the plasma temperature and confinement time. The radial pressure profiles for typical Low confinement mode (L-mode) and H-mode discharges are depicted in Fig. 1. It can be seen that there is a significant increase in the core pressure, from the L-mode discharge to the H-mode discharge. The significant enhancement in the plasma performance is the result of a transport barrier that forms at the edge of the plasma. This edge transport barrier (ETB) is usually referred as the pedestal. Typically, the energy content in an H-mode discharge is approximately twice the energy contained in an L- mode discharge, for the plasma heated with the same input power [2]. To understand better the physical processes that take place in tokamak plasma, many advanced computer codes have been developed. The integrated predictive modelling codes, such as BALDUR [3], TASK/TR [4], JETTO [5], ASTRA [6] and CRONOS [7], have played an important role in the simulations to predict the temporal evolution of plasma current, temperature, and density profiles. They lead to significant improvement in understanding of plasmas behaviours in tokamaks. Normally, the simulations carried out with these integrated predictive modelling codes making use of boundary conditions taken from experimental data. In simulating H-mode discharges, the evolution of the core plasma was carried out using boundary conditions taken from experimental data at the top of the pedestal [8 13], because the distance from the top of the pedestal to the center of plasma is approximately 95% of the minor radius, as shown in Fig. 1. Several years after discovery of the H-mode, however, different physical phenomena were found at the plasma edge such as the formation of a pedestal that affects the plasma core. These phenomena suggest that the integrated predictive modelling code must be extended to include the entire plasma cross section, because a reliable methodology to predict the plasma properties is needed in order to advance the predictive capability. This capability is essential in designing future experiments in existing tokamaks. Fig.1 Pressure profiles for both L-mode and H-mode plasma supported by the Commission on Higher Education and the Thailand Research Fund (No. RSA ) and the Government Annual Research Budget through Thammasat University

2 Y. PIANROJ et al.: Simulations of H-Mode Plasmas in Tokamak Using a Complete Core-Edge Modeling The requirement for the formation of pedestal is the reduction of fluctuation-driven transport. This can be achieved by stabilization or decorrelation of microturbulence in the plasma. The stabilization mechanisms, which can suppress turbulent modes, have to account for the different dynamical behaviors of the various species in the plasma. The first candidate for edge turbulence stabilization is the ω E B flow shear. The ω E B flow shear can suppress turbulence through linear stabilization of turbulent modes, as well as in particular by non-linear decorrelation of turbulence vortices [14 16], thereby reducing transport by acting on both the amplitude of the fluctuations and the phase between them [17]. The second candidate is the magnetic shear s, which reduces transport only in the region where the shear s exceeds a threshold and is unaffected elsewhere. In Ref. [18], the ω E B flow shear was used as a mechanism for reduction of ion and electron thermal diffusivity in order to initiate the formation of a pedestal. This work was carried out using ASTRA by A. Y. PANKIN et al. However, the hydrogenic and impurity transports are not reduced due to some limitations. Moreover, G. W. Pacher et al. [18] implemented the suppression function which consists of two stabilization candidate terms, namely ω E B flow shear and magnetic shear, into ASTRA to carry out a prediction for ITER. Both temperature and density profile evolutions were considered with a simple anomalous core transport model. In this work, a core transport model, namely in mixed Bohm/gyro-Bohm (mixed B/gB) model, was used, by suppressing every channel of transport coefficient, i.e., electron thermal diffusivity, ion thermal diffusivity, hydrogenic diffusivity, and impurity diffusivity to describe the pedestal formation. As a result, full core-edge simulations of tokamak plasmas can be conducted. This paper is organized as follows. The edge model with suppression function will be described in section 2. A brief introduction of the mixed B/gB anomalous core transport is given in section 3. In section 4, the calibration and the sensitivity of coefficient C x and the simulation results for typical H-mode will be validated by the statistical comparisons with the experimental data from DIII-D. The final section is for the conclusion. 2 Model for edge plasma The modelling of pedestal structure inevitably requires the full integration of plasma in both core and scrape-off layer (SOL)/divertor regions. One of the difficulties in modelling the pedestal structure is that there exist different physical mechanisms with different time scales. The pedestal structure evolves on a transport time scale. However during the profile evolution, magnetohydrodynamic (MHD) phenomena, which occur on short time scales, can develop, such as edge localized modes (ELMs) [19]. The occurrence of an ELM event burst gives rise a significant pulsed flow of both particles and energy onto the divertor target, diminishing the edge pressure gradients in the process. A paint for the pedestal transport modelling is to extend the anomalous core transport to include the effect of ω E B flow shearing rate and magnetic shear. In this work, the anomalous transport in both the core and pedestal is taken in forms of: χ is = χ i f sion, (1) χ es = χ e f selectron, (2) D Hs = D H f shydrogenic, (3) D zs = D z f simpurity, (4) where χ is and χ es are the modified anomalous ion and electron thermal diffusivity, D Hs and D zs are the modified anomalous particle and impurity diffusivity, respectively; f sion and f selectron are the suppression function for ion and electron thermal diffusivity, f shydrogenic and f simpurity are the suppression function for particle and impurity diffusivity, respectively. An appropriate transport suppression function due to ω E B flow shearing rate together with the reduction of turbulence growth rate [20 22] is represented in the first term in which the E B flow shear alone leads to a formation of pedestals which are appreciably lower than those obtained experimentally. An additional magnetic shear stabilization is realized in the second term [18,23]. The transport is reduced only in the region where the magnetic shear exceeds a threshold (in this work the threshold is set to 0.5 [18] ) and the suppression function (f sx ) can be written in 1 f sx = 1 + C x ( ω E B ) 1 2 max(1, (s 0.5) 2 ), (5) γ ITG where, C x is the coefficient for each species calibrated later in section 4.2.1, γ ITG is the growth rate of ion temperature gradient mode (ITG), estimated as V ti /qr [24] with V ti the ion thermal velocity and s the magnetic shear. 3 Model for core plasma The mixed B/gB anomalous core transport model [25] is an empirical model. It consists of Bohm and gyro- Bohm contributions. The Bohm term is linear in the gyro-radius and is a non-local transport model, in which the transport depends on a finite difference approximation to the electron temperature gradient at the edge of plasma. The gyro-bohm term is a local transport model, which is added for simulation to be able to match data from both smaller and larger tokamaks. It is proportional to the square of the gyro-radius times thermal velocity over the square of the plasma dimension. Thus, the mixed B/gB transport model can be expressed as [26] : χ e = 1.0χ gb + 2.0χ B, (6) χ i = 0.5χ gb + 4.0χ B + χ neo, (7) χ e χ i D H = [ ρ], χ e + χ i (8) 779

3 Plasma Science and Technology, Vol.14, No.9, Sep where D z = D H, (9) χ gb = T e T e, (10) B 2 φ χ B = R (n et e ) ( q 2 T e,0.8 T ) e,1.0, (11) n e B φ T e,1.0 where, χ e is the electron diffusivity, χ i is the ion diffusivity, D H is the particle diffusivity, D z is the impurity diffusivity, χ gb is the gyro-bohm contribution, χ B is the Bohm contribution, ρ is the normalized minor radius, T e is the electron temperature, B φ is the toroidal magnetic field, R is the major radius, n e is the local electron density, q is the safety factor, s is the magnetic shear, ω E B is the flow shearing rate, and γ ITG is the growth rate of ion temperature gradient (ITG). 4 Simulation results and discussions 4.1 Simulated plasma profiles The experimental data from ten DIII-D H- mode shots are taken from the international profile database [27]. These experimental data are classified into four pairs of H-mode plasmas which comprise the DIII-D discharges examined in this paper. The major plasma parameters for all of ten shots are listed in Table 1. They are used as the initial and boundary conditions for both core and core-edge simulations while the discharges include scans for plasma power (shots and 77559), density (shots and 81329), elongation (κ) (shots and 81507), and normalized gyroradius (ρ ) (shots and 82788). However, the simulation results carried out by the core model are calculated at the top of pedestal which depends on the boundary conditions of each shot, but the simulation results carried out by the core-edge model are calculated at the plasma edge with the same temperature 10 ev, and compared with the DIII-D data. Note that the change of edge temperature does not affect the simulation results. First, consider the elongation scan (shots and 81507). The plasma elongation (at the 95% flux surface) varies from 1.68 to 1.95 while the minor radius of plasma varies from 0.63 m to 0.54 m correspondingly. The simulation results using the core-edge model match the experimental data quite well for both shots in this scan, as shown in Fig. 2. The shapes of all profiles remained nearly unchanged in the two shots of this scan, in both the experimental and the simulation results except for some structures in the core region where the magnitudes from simulation are lower than those in experiment. Unfortunately, there is no result for impurity in the database for these DIII-D shots. In Fig. 3, the simulated plasma profiles are compared to the corresponding experimental ones for the DIII-D ρ scan represented (shots and 82788). It is found that the simulated temperature profiles of shot (low ρ ) can match the experimental data. However, the simulated density profiles of this shot do not match the experimental data. On the other hand, the temperature profiles in core region of shot (high ρ ) tends to over-predict both electron s and ion s, while the density profiles of both electron and deuterium match very well. In the power scan (shots and 77559), the neutral beam injection power varies while the average plasma density is kept constant. In the density scan (shots and 81329), the heating power was adjusted to hold the temperatures fixed as the plasma density varies. It is found in the simulations that the temperature profiles of these four discharges do not match the experimental well in the core region. They tend to over-predict the experimental results, except in the case of low density for the density scan. On the other hand, the density profiles of these four shots match the experimental data quite well and most simulations tend to produce higher prediction than the experimental data. The last two shots (shots and 82183) are not a pair of experimental discharges. It is found that the temperature profiles and density profiles Table 1. Plasma parameters for ten shots in DIII-D Discharges Type Low power High power Low n e High n e Low κ High κ Low ρ High ρ R (m) a (m) κ δ B T(T) I p(ma) n e( m 3 ) Z eff P NB(MW) Diagnostic time (s)

4 Y. PIANROJ et al.: Simulations of H-Mode Plasmas in Tokamak Using a Complete Core-Edge Modeling Fig.2 Profiles of electron temperature, ion temperature, electron density and deuterium density. The simulated results by using BALDUR with core-edge transport model are compared to the DIII-D data of shots (low κ; left panel) and (high κ; right panel) Fig.3 Profiles of electron temperature, ion temperature, electron density and deuterium density. The simulated results by using BALDUR with core-edge transport model are compared to the DIII-D data of shots (low ρ ; left panel) and (high ρ ; right panel) of shot do not match the experimental ones quite well and they mostly under-predict the experimental data. The simulations results carried out for shot are different when compared to those for shot In shot 82183, the temperature and density profiles match very well to the experimental data except in deuterium profiles, which significantly under-predict the experimental data. A statistical analysis (presented in section 4.2) will be carried out to quantify the agreement between the core-edge model predictions and the corre- sponding experimental data. The effect that helps to reduce the anomalous transport and to form the transport barrier at the plasma edge is described by the suppression function that is shown in Eq. (1) and Fig. 4 (shots for low κ, and for high κ) and Fig. 5 (shots for low ρ and for high ρ ). This function is composed of two terms which reduce the turbulent transport. The first term is for shear in the E B flow. In this work, it is effective in the edge region only (at about a normalized 781

5 Plasma Science and Technology, Vol.14, No.9, Sep Fig.4 Profiles of radial electric field, ω E B flow shear, magnetic shear and suppression function obtained by using BALDUR with core-edge transport model are compared to the DIII-D data of shots (low κ; left panel) and (high κ; right panel) Fig.5 Profiles of radial electric field, ω E B flow shear, magnetic shear, and suppression function obtained by using BALDUR with core-edge transport model are compared to the DIII-D data of shots (low ρ ; left panel) and (high ρ ; right panel) minor radius (r/a) = ). To demonstrate the E B flow shear stabilization effect, the radial electric field (E r ) calculated by BALDUR is shown. It seems to be the same pattern and stronger near the edge area because the pressure gradient is high there. Thus, the radial electric field will affect the ω E B flow shear [28] that can be depicted by: where, R is the major radius, B θ is the poloidal magnetic field and B φ is the toroidal magnetic field. The ω E B flow shear calculated by the predictive modelling codes is shown in Figs. 4 and 5, which show the same trend as the radial electric field with a high magnitude at the plasma edge. It is worth to note that the radial electric field, E r, is calculated by 782 ω E B = (RB θ) 2 ( ) Er, (12) B φ Ψ RB θ E r = 1 Zen i p i r v θb φ + v tor B θ, (13)

6 Y. PIANROJ et al.: Simulations of H-Mode Plasmas in Tokamak Using a Complete Core-Edge Modeling p i where, r is the pressure gradient, υ θ and υ tor are the poloidal and toroidal velocities, respectively, n i is the ion density, Z is the ion charge number, and e is the elementary charge. The calculated toroidal velocity, v tor, is assumed to be a function of ion temperature [29] ; whereas the poloidal velocity, v θ is estimated using the NCLASS code. The last term, the magnetic shear, also plays a key role in facilitating entry into enhanced confinement or low magnetic shear acts to reduce the growth rate of turbulence [19,29]. Thus, the magnetic shear profile simulated by BALDUR is shown in these figures. The magnetic shear profile increases rapidly at r/a = 0.6 to 1.0. Therefore, the inverse of maximum function between 1 to (s 0.5) 2 of this term is effective in term of suppression at the plasma edge area as well. This behaviour of these parameters for suppression is similar in all simulations in this work. 4.2 Statistical analysis To quantify the comparison between simulations and experimental results, a percentage of the root-meansquare error (%RMSE) deviation and the %offset are calculated based on the difference between simulated and experimental results. In this paper, the %RMSE and the %offset are defined respectively as RMSE(%) = 1 N offset(%) = 1 N N i=1 N i=1 ( Xsimi X expi X exp0 ) 2 100, (14) ( Xsimi X ) expi 100, (15) X exp0 where, X expi is the i th data point of the experimental profile, X simi is the corresponding data point of the simulation profile, and X exp0 is the maximum data point of the experimental profile of X as a function of radius, with N points totally. The %RMSE and the %offset are evaluated for each of the four profiles, namely: electron temperature, ion temperature, electron density, and deuterium density for the shots considered. Note that the %offset is positive if the simulated data is larger than the experimental data and negative if the simulated is smaller than the experimental data. If the %offset is zero, then the %RSME is a measure of how much the shape of the profiles differ between simulated and experimental one, thus the statistical analysis will be used to analyze the coefficients C x and the simulated and experimental data from DIII-D H-mode discharges that will be described in the next section Calibration of coefficients C x In the pedestal transport model, described in Eq. (5), the coefficient, C x for each species is chosen to optimize the agreement in the data from the simulation and the DIII-D H-mode shot for both core and edge regions. It is found that the ion coefficient (C i ) and the electron coefficient (C e ) are equal to , which yield the best agreement. For the hydrogenic coefficient and the impurity coefficient, the best choices are and , respectively. The effect of different coefficients, by a comparison to the experimental data of DIII-D is shown in Fig. 6, with the values of coefficients listed in Table 2. It depicts the set of coefficients that are used in this work C x1, the set of coefficients that are less than ten times that of C x1 are C x2, and the set of coefficients that are more than ten times that of C x1 are C x3. These coefficients affect the electron temperature, ion temperature, electron density, and deuterium density profiles. Fig.6 Profiles of electron temperature, ion temperature electron density and deuterium density. The simulated results by using BALDUR in order to show the calibration and sensitivity of coefficients C x, are compared to the DIII-D data of shot

7 Plasma Science and Technology, Vol.14, No.9, Sep Simulation results compared with 10 DIII-D H-mode discharges In this section, the simulation results carried out by BALDUR are validated with the data from ten shots in DIII-D by using statistical analysis. The results are presented in Figs and listed in Table 3. It should be noted that the statistical analysis in this work consists of two groups. The first group is named core model which is used to quantify all the simulated and experimental data when the anomalous transport mixed B/gB model is used from the center of plasma to the top of pedestal in which the boundary is taken from experiments at the top of pedestal at the diagnostic time. Another group is named core-edge model which is used to quantify all the simulated and experimental data when the anomalous transport mixed B/gB model is used with the suppression function from the center to the edge of plasma (separatrix). Table 2. The %RMSEs and %offsets for the electron temperature profiles using either the core model or using the core-edge model for ten DIII-D shots are shown in Fig. 7. The %RMSEs of the electron temperature using the core model range from 4.50% to 8.38%. However, with the core-edge model, the %RMSEs range from 6.88% to 18.71%. The %offsets using the core model are mostly positive, indicating that simulations tend to over-predict the experimental results. However, the %offsets using the core-edge model are mostly negative, indicating that simulations tend to under-predict the experimental results. In Fig. 8, the %RMSEs of the ion temperature using the core model range from 3.83% to 12.21% while those using the core-edge model range from 6.88% to 18.71%. The %offsets of simulations using both core and core-edge model are mostly positive, indicating that simulations tend to over-predict the experimental results. %RMSE and %offset for the calibration and sensitivity of coefficients Cx Profiles Cx1 Cx2 Cx3 %RMSE %offset %RMSE %offset %RMSE %offset Te Ti ne nd Fig.7 Percentage of root mean square error (%RMSE) for the electron temperature profiles obtained by simulations using both core model (mixed B/gB) and core-edge model in BALDUR, compared to the experimental data from ten H-mode shots (pedestal occurred) in DIII-D 784

8 Y. PIANROJ et al.: Simulations of H-Mode Plasmas in Tokamak Using a Complete Core-Edge Modeling Fig.8 Percentage of root mean square error (%RMSE) for the ion temperature profiles obtained by the simulation using both core model (mixed B/gB) and core-edge model in BALDUR, compared to the experimental data from ten H-mode shots (pedestal occurred) in DIII-D Fig.9 Percentage of root mean square error (%RMSE) for the electron density profiles obtained by the simulation using both core model (mixed B/gB) and core-edge model in BALDUR, compared to the experimental data from ten H-mode shots (pedestal occurred) in DIII-D In Fig. 9, the %RMSE for the simulations for electron density using the core model ranges from 5.37% to 25.51% while that of simulations using the coreedge model ranges from 4.03% to 13.21%. The %offsets in the case of the simulations using the core model are mostly negative, indicating that simulations under- predict the experimental results. However, the %offsets in the case of the simulations using the core-edge model are mostly positive. Finally in Fig. 10, the %RMSE of the simulations for deuterium density using the core model ranges from 5.66% to 21.51% and the %RMSE of the simulations using the core-edge model ranges from 785

9 Plasma Science and Technology, Vol.14, No.9, Sep % to 17.66%. The %offsets are mostly negative for the simulations using both models, indicating that simulations tend to under-predict the experimental data. The average of %RMSE and average of %offset, averaged over ten shots for electron temperature, ion temperature, electron density, and deuterium density, are shown in Fig. 11. The averaged %RMSE differs in less than 6% between the two models for four profiles, and the averaged %offsets differs in less than 8%. A com- parison between the results using the core model to the core-edge mode for the simulation of temperature profiles indicates that the average of %RMSE increase and the simulations over-predict the experimental results because of the positive averages of %offset. For density profiles, the average of %RMSE for electron density decreases in some shots (77559, 81499, and 82788) and the simulations over-predict the experimental results. Fig.10 Percentage of root mean square error (%RMSE) for the deuterium density profiles obtained by the simulation using core model (mixed B/gB) and core-edge model in BALDUR, compared to the experimental data from ten H-mode shots (pedestal occurred) in DIII-D Table 3. Discharges %RMSE and %offset for both core and core-edge models Core model Core-Edge model %RMSE Te Ti %offset %RMSE Ti ne %offset ne nd Te Ti ne nd Te nd Te Ti ne nd

10 Y. PIANROJ et al.: Simulations of H-Mode Plasmas in Tokamak Using a Complete Core-Edge Modeling Fig.11 Averaged percentage of root mean square error (%RMSE) and the averaged percentage of offset for four profiles of electron temperature (Te ), ion temperature (Ti ), electron density (ne ) and deuterium density (nd ) obtained by simulations using both the core and core-edge models, compared to the experimental data from ten H-mode shots in DIII-D 5 Conclusion A theory-based model for describing the pedestal transport is developed. An anomalous transport model for the core region is extended to the pedestal region by the inclusion of the effect of both ωe B flow shear and magnetic shear. This model is developed and tested in BALDUR for self-consistently simulating H-mode plasmas. It is found that the simulations successfully reproduce the experimental results. The core simulations carried out using BALDUR with the mixed B/gB core model and prescribed boundary conditions yield the simulations with average RMSEs between 6.38% and 12.62% for both temperature and density profiles. For the core-edge simulations, the agreement changes slightly with a range between 7.59% and 14.53%. As a result, the simulations using the core-edge model with the pedestal model developed yield an agreement comparable to those obtained from the simulations using the core model with the boundary conditions taken from experiment. It is found that the suppression function developed is effective only in the edge region. References Acknowledgments Mr. YUTTHAPONG thanks for the program Strategic Scholarships for Frontier Research Network for Thai Ph.D. Program Wagner F, Becker G, Behringer K, et al. 1982, Physical Review Letters, 49: 1408 Connor J W, Wilson H R. 2000, Plasma Physics and Controlled Fusion, 42: R1 Singer C E, Post D E, Mikkelsen D R. 1988, Comput. Phys. Commun., 49: 275 Honda M, Fukuyama A. 2006, Nuclear Fusion, 46: 580 Cenacchi G, Taroni A. 1988, JET-IR, 88, Technical Report, Joint European Torus Underlaking, UK Pereverzev G, Yushmanov P N. 2002, IPP 5/98, Technical Report, Max-Planck Institut fur Plasmaphysik, Garching, Germany Basiuk V, Artaud J F, Imbeaux F, et al. 2003, Nuclear Fusion, 43: 822 Kinsey J E, Staebler G M, Waltz R E. 2002, Physics of Plasmas, 9: 1676 Bateman G, Kritz A H, Kinsey J E. 1998, Physics of Plasmas, 5: 1793 Bateman G, Onjun T, Kritz A H. 2003, Plasma Physics and Controlled Fusion, 45: 1939 Kinsey J E, Bateman G, Onjun T, et al. 2003, Nuclear Fusion, 43: 1845 Onjun T, Bateman G, Kritz A H. 2001, Physics of Plasmas, 8: 975 Hannum D, Bateman G, Kinsey J. 2001, Physics of Plasmas, 8: 964 Biglari H, Diamond P H, Terry P W. 1990, Physics of Fluids B: Plasma Physics, 2: 1 Gohil P. 2006, Comptes Rendus Physique, 7: 606 Boedo J, Gray D, Jachmich S, et al. 2000, Nuclear 787

11 Plasma Science and Technology, Vol.14, No.9, Sep Fusion, 40: Oost G V, Adamek J, Antoni V, et al. 2003, Plasma Physics and Controlled Fusion, 45: 1 18 Pacher G W, et al. 2004, Plasma Physics and Controlled Fusion, 46: A Doyle E J, Houlberg W A, Kamada Y, et al. 2007, Nuclear Fusion, 47: S18 20 Zolotukhin O V, Igitkhanov Y, Janeschitz G, et al. 2001, Transport simulations of type I ELMs. Presented at the 28th EPS Conf. on Contr. Fusion and Plasma Phys., (Funchal, Portugal, June 2001), European Physical Society, Geneva 21 Sugihara M, Igitkhanov Y, Janeschitz G, et al. 2001, Simulation studies on H-mode pedestal behavioru during type-i ELMs under various plasma conditions. Presented at the 28th EPS Conf. on Contr. Fusion and Plasma Phys., (Funchal, Portugal, June 2001), European Physical Society, Geneva 22 Pankin A Y, Voitsekhovitch I, Bateman G, et al. 2005, Plasma Physics and Controlled Fusion, 47: Janeschitz G, Pacher G W, Zolotukhin O, et al. 2002, Plasma Physics and Controlled Fusion, 44: A Tala T. 2002, Transport Barrier and Current Profile Studies on the JET Tokamak. [Ph.D] Helsinki University of Technology, Espoo, Finland 25 Tala T J J, Heikkinen J A, Parail V V. 2001, Plasma Physics and Controlled Fusion, 43: Onjun T, Pianroj Y. 2009, Nuclear Fusion, 49: Boucher D, Connor J W, Houlberg W A, et al. 2000, Nuclear Fusion, 40: Burrell K H. 1997, Physics of Plasmas, 4: Waltz R E, Staebler G M, Dorland W, et al. 1997, Physics of Plasmas, 4: 2482 (Manuscript received 21 February 2011) (Manuscript accepted 3 June 2011) address of corresponding author T. ONJUN: 788

A THEORETICAL AND EXPERIMENTAL INVESTIGATION INTO ENERGY TRANSPORT IN HIGH TEMPERATURE TOKAMAK PLASMAS

A THEORETICAL AND EXPERIMENTAL INVESTIGATION INTO ENERGY TRANSPORT IN HIGH TEMPERATURE TOKAMAK PLASMAS A THEORETICAL AND EXPERIMENTAL INVESTIGATION INTO ENERGY TRANSPORT IN HIGH TEMPERATURE TOKAMAK PLASMAS Presented by D.P. SCHISSEL Presented to APS Centennial Meeting March 20 26, 1999 Atlanta, Georgia

More information

GA A23736 EFFECTS OF CROSS-SECTION SHAPE ON L MODE AND H MODE ENERGY TRANSPORT

GA A23736 EFFECTS OF CROSS-SECTION SHAPE ON L MODE AND H MODE ENERGY TRANSPORT GA A3736 EFFECTS OF CROSS-SECTION SHAPE ON L MODE AND H MODE ENERGY TRANSPORT by T.C. LUCE, C.C. PETTY, and J.E. KINSEY AUGUST DISCLAIMER This report was prepared as an account of work sponsored by an

More information

Plasma instability during ITBs formation with pellet injection in tokamak

Plasma instability during ITBs formation with pellet injection in tokamak Plasma instability during ITBs formation with pellet injection in tokamak P. Klaywittaphat 1, B. Chatthong 2, T. Onjun. R. Picha 3, J. Promping 3 1 Faculty of Engineering, Thaksin University, Phatthalung,

More information

EFFECT OF EDGE NEUTRAL SOUCE PROFILE ON H-MODE PEDESTAL HEIGHT AND ELM SIZE

EFFECT OF EDGE NEUTRAL SOUCE PROFILE ON H-MODE PEDESTAL HEIGHT AND ELM SIZE EFFECT OF EDGE NEUTRAL SOUCE PROFILE ON H-MODE PEDESTAL HEIGHT AND ELM SIZE T.H. Osborne 1, P.B. Snyder 1, R.J. Groebner 1, A.W. Leonard 1, M.E. Fenstermacher 2, and the DIII-D Group 47 th Annual Meeting

More information

Comparison of ITER Performance Predicted by Semi-Empirical and Theory-Based Transport Models

Comparison of ITER Performance Predicted by Semi-Empirical and Theory-Based Transport Models 1 CT/P-3 Comparison of ITER Performance Predicted by Semi-Empirical and Theory-Based Transport Models V. Mukhovatov 1), Y. Shimomura 1), A. Polevoi 1), M. Shimada 1), M. Sugihara 1), G. Bateman 2), J.G.

More information

Progress in characterization of the H-mode pedestal

Progress in characterization of the H-mode pedestal Journal of Physics: Conference Series Progress in characterization of the H-mode pedestal To cite this article: A W Leonard 2008 J. Phys.: Conf. Ser. 123 012001 View the article online for updates and

More information

Improved 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) 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 information

Core Transport Properties in JT-60U and JET Identity Plasmas

Core Transport Properties in JT-60U and JET Identity Plasmas 1 EXC/P4-12 Core Transport Properties in JT-60U and JET Identity Plasmas X. Litaudon 1, Y. Sakamoto 2, P.C. de Vries 3, A. Salmi 4, T. Tala 5, C. Angioni 6, S. Benkadda 7, M.N.A. Beurskens 8, C. Bourdelle

More information

Scaling of the Density Peak with Pellet Injection in ITER

Scaling of the Density Peak with Pellet Injection in ITER Plasma Science and Technology, Vol.14, No.12, Dec. 2012 Scaling of the Density Peak with Pellet Injection in ITER P. KLAYWITTAPHAT and T. ONJUN School of Manufacturing Systems and Mechanical Engineering,

More information

Integrated Simulation of ELM Energy Loss Determined by Pedestal MHD and SOL Transport

Integrated Simulation of ELM Energy Loss Determined by Pedestal MHD and SOL Transport 1 Integrated Simulation of ELM Energy Loss Determined by Pedestal MHD and SOL Transport N. Hayashi, T. Takizuka, T. Ozeki, N. Aiba, N. Oyama Japan Atomic Energy Agency, Naka, Ibaraki-ken, 311-0193 Japan

More information

QTYUIOP ENERGY TRANSPORT IN NEUTRAL BEAM HEATED DIII D DISCHARGES WITH NEGATIVE MAGNETIC SHEAR D.P. SCHISSEL. Presented by. for the DIII D Team*

QTYUIOP ENERGY TRANSPORT IN NEUTRAL BEAM HEATED DIII D DISCHARGES WITH NEGATIVE MAGNETIC SHEAR D.P. SCHISSEL. Presented by. for the DIII D Team* ENERGY TRANSPORT IN NEUTRAL BEAM HEATED DIII D DISCHARGES WITH NEGATIVE MAGNETIC SHEAR Presented by D.P. SCHISSEL for the DIII D Team* Presented to 38th APS/DPP Meeting NOVEMBER 11 15, 1996 Denver, Colorado

More information

Current 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 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 information

Towards Multiscale Gyrokinetic Simulations of ITER-like Plasmas

Towards Multiscale Gyrokinetic Simulations of ITER-like Plasmas Frank Jenko Max-Planck-Institut für Plasmaphysik, Garching Universität Ulm Towards Multiscale Gyrokinetic Simulations of ITER-like Plasmas 23 rd IAEA Fusion Energy Conference 11-16 October 2010, Daejeon,

More information

ITER Predictions Using the GYRO Verified and Experimentally Validated TGLF Transport Model

ITER Predictions Using the GYRO Verified and Experimentally Validated TGLF Transport Model 1 THC/3-3 ITER Predictions Using the GYRO Verified and Experimentally Validated TGLF Transport Model J.E. Kinsey, G.M. Staebler, J. Candy, and R.E. Waltz General Atomics, P.O. Box 8608, San Diego, California

More information

EX/C3-5Rb Relationship between particle and heat transport in JT-60U plasmas with internal transport barrier

EX/C3-5Rb Relationship between particle and heat transport in JT-60U plasmas with internal transport barrier EX/C-Rb Relationship between particle and heat transport in JT-U plasmas with internal transport barrier H. Takenaga ), S. Higashijima ), N. Oyama ), L. G. Bruskin ), Y. Koide ), S. Ide ), H. Shirai ),

More information

Particle Transport and Density Gradient Scale Lengths in the Edge Pedestal

Particle Transport and Density Gradient Scale Lengths in the Edge Pedestal Particle Transport and Density Gradient Scale Lengths in the Edge Pedestal W. M. Stacey Fusion Research Center, Georgia Institute of Technology, Atlanta, GA, USA Email: weston.stacey@nre.gatech.edu Abstract

More information

GA A26887 ADVANCES TOWARD QH-MODE VIABILITY FOR ELM-FREE OPERATION IN ITER

GA A26887 ADVANCES TOWARD QH-MODE VIABILITY FOR ELM-FREE OPERATION IN ITER GA A26887 ADVANCES TOWARD QH-MODE VIABILITY FOR ELM-FREE OPERATION IN ITER by A.M. GAROFALO, K.H. BURRELL, M.J. LANCTOT, H. REIMERDES, W.M. SOLOMON and L. SCHMITZ OCTOBER 2010 DISCLAIMER This report was

More information

Integrated Simulation of ELM Energy Loss and Cycle in Improved H-mode Plasmas

Integrated Simulation of ELM Energy Loss and Cycle in Improved H-mode Plasmas 1 Integrated Simulation of ELM Energy Loss and Cycle in Improved H-mode Plasmas N. Hayashi 1), T. Takizuka 1), N. Aiba 1), N. Oyama 1), T. Ozeki 1), S. Wiesen 2), V. Parail 3) 1) Japan Atomic Energy Agency,

More information

Direct drive by cyclotron heating can explain spontaneous rotation in tokamaks

Direct 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 information

C-Mod Core Transport Program. Presented by Martin Greenwald C-Mod PAC Feb. 6-8, 2008 MIT Plasma Science & Fusion Center

C-Mod Core Transport Program. Presented by Martin Greenwald C-Mod PAC Feb. 6-8, 2008 MIT Plasma Science & Fusion Center C-Mod Core Transport Program Presented by Martin Greenwald C-Mod PAC Feb. 6-8, 2008 MIT Plasma Science & Fusion Center Practical Motivations for Transport Research Overall plasma behavior must be robustly

More information

DT Fusion Power Production in ELM-free H-modes in JET

DT Fusion Power Production in ELM-free H-modes in JET JET C(98)69 FG Rimini and e JET Team DT Fusion ower roduction in ELM-free H-modes in JET This document is intended for publication in e open literature. It is made available on e understanding at it may

More information

ELM Suppression in DIII-D Hybrid Plasmas Using n=3 Resonant Magnetic Perturbations

ELM Suppression in DIII-D Hybrid Plasmas Using n=3 Resonant Magnetic Perturbations 1 EXC/P5-02 ELM Suppression in DIII-D Hybrid Plasmas Using n=3 Resonant Magnetic Perturbations B. Hudson 1, T.E. Evans 2, T.H. Osborne 2, C.C. Petty 2, and P.B. Snyder 2 1 Oak Ridge Institute for Science

More information

Recent Development of LHD Experiment. O.Motojima for the LHD team National Institute for Fusion Science

Recent Development of LHD Experiment. O.Motojima for the LHD team National Institute for Fusion Science Recent Development of LHD Experiment O.Motojima for the LHD team National Institute for Fusion Science 4521 1 Primary goal of LHD project 1. Transport studies in sufficiently high n E T regime relevant

More information

Dimensionless Identity Experiments in JT-60U and JET

Dimensionless Identity Experiments in JT-60U and JET 1 IAEA-CN-116/IT/1-2 Dimensionless Identity Experiments in JT-60U and JET G Saibene 1, N Oyama 2, Y Andrew 3, JG Cordey 3, E de la Luna 4, C Giroud 3, K Guenther 3, T Hatae 2, GTA Huysmans 5, Y Kamada

More information

Introduction to Fusion Physics

Introduction 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 information

The EPED Pedestal Model: Extensions, Application to ELM-Suppressed Regimes, and ITER Predictions

The EPED Pedestal Model: Extensions, Application to ELM-Suppressed Regimes, and ITER Predictions The EPED Pedestal Model: Extensions, Application to ELM-Suppressed Regimes, and ITER Predictions P.B. Snyder 1, T.H. Osborne 1, M.N.A. Beurskens 2, K.H. Burrell 1, R.J. Groebner 1, J.W. Hughes 3, R. Maingi

More information

GA A25853 FAST ION REDISTRIBUTION AND IMPLICATIONS FOR THE HYBRID REGIME

GA A25853 FAST ION REDISTRIBUTION AND IMPLICATIONS FOR THE HYBRID REGIME GA A25853 FAST ION REDISTRIBUTION AND IMPLICATIONS FOR THE HYBRID REGIME by R. NAZIKIAN, M.E. AUSTIN, R.V. BUDNY, M.S. CHU, W.W. HEIDBRINK, M.A. MAKOWSKI, C.C. PETTY, P.A. POLITZER, W.M. SOLOMON, M.A.

More information

DEPENDENCE OF THE H-MODE PEDESTAL STRUCTURE ON ASPECT RATIO

DEPENDENCE OF THE H-MODE PEDESTAL STRUCTURE ON ASPECT RATIO 21 st IAEA Fusion Energy Conference Chengdu, China Oct. 16-21, 2006 DEPENDENCE OF THE H-MODE PEDESTAL STRUCTURE ON ASPECT RATIO R. Maingi 1, A. Kirk 2, T. Osborne 3, P. Snyder 3, S. Saarelma 2, R. Scannell

More information

Role of Flow Shear in Enhanced Core Confinement

Role of Flow Shear in Enhanced Core Confinement PPPL-3156, Preprint: March 1996, UC-420, 427 Role of Flow Shear in Enhanced Core Confinement Regimes TS Hahm and KH Burrell The importance of the ExB fla shear in arioi s enhanced confinement regimes is

More information

Gyrokinetic Transport Driven by Energetic Particle Modes

Gyrokinetic 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 information

Divertor power deposition and target current asymmetries during type-i ELMs in ASDEX Upgrade and JET

Divertor power deposition and target current asymmetries during type-i ELMs in ASDEX Upgrade and JET Journal of Nuclear Materials 363 365 (2007) 989 993 www.elsevier.com/locate/jnucmat Divertor power deposition and target current asymmetries during type-i ELMs in ASDEX Upgrade and JET T. Eich a, *, A.

More information

Effect of the Radial Electric Field on Lower Hybrid Plasma Heating in the FT-2 Tokamak

Effect of the Radial Electric Field on Lower Hybrid Plasma Heating in the FT-2 Tokamak Plasma Physics Reports, Vol. 7, No.,, pp.. Translated from Fizika Plazmy, Vol. 7, No.,, pp. 9 9. Original Russian Text Copyright by Lashkul, Budnikov, Vekshina, D yachenko, Ermolaev, Esipov, Its, Kantor,

More information

Phase ramping and modulation of reflectometer signals

Phase ramping and modulation of reflectometer signals 4th Intl. Reflectometry Workshop - IRW4, Cadarache, March 22nd - 24th 1999 1 Phase ramping and modulation of reflectometer signals G.D.Conway, D.V.Bartlett, P.E.Stott JET Joint Undertaking, Abingdon, Oxon,

More information

STABILIZATION OF m=2/n=1 TEARING MODES BY ELECTRON CYCLOTRON CURRENT DRIVE IN THE DIII D TOKAMAK

STABILIZATION 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 information

GA A26874 ITER PREDICTIONS USING THE GYRO VERIFIED AND EXPERIMENTALLY VALIDATED TGLF TRANSPORT MODEL

GA A26874 ITER PREDICTIONS USING THE GYRO VERIFIED AND EXPERIMENTALLY VALIDATED TGLF TRANSPORT MODEL GA A26874 ITER PREDICTIONS USING THE GYRO VERIFIED AND EXPERIMENTALLY VALIDATED TGLF TRANSPORT MODEL by J.E. KINSEY, G.M. STAEBLER, J. CANDY and R.E. WALTZ NOVEMBER 20 DISCLAIMER This report was prepared

More information

Impact of neutral atoms on plasma turbulence in the tokamak edge region

Impact of neutral atoms on plasma turbulence in the tokamak edge region Impact of neutral atoms on plasma turbulence in the tokamak edge region C. Wersal P. Ricci, F.D. Halpern, R. Jorge, J. Morales, P. Paruta, F. Riva Theory of Fusion Plasmas Joint Varenna-Lausanne International

More information

J. Kesner. April Plasma Fusion Center Massachusetts Institute of Technology Cambridge, Massachusetts USA

J. Kesner. April Plasma Fusion Center Massachusetts Institute of Technology Cambridge, Massachusetts USA PFC/JA-88-38 Effect of Local Shear on Drift Fluctuation Driven T'ransport in Tokamaks J. Kesner April 1989 Plasma Fusion Center Massachusetts Institute of Technology Cambridge, Massachusetts 2139 USA Submitted

More information

Approach to Steady-state High Performance in DD and DT with Optimised Shear on JET

Approach to Steady-state High Performance in DD and DT with Optimised Shear on JET JET P(98) F X Söldner et al Approach to Steady-state High Performance in DD and DT with Optimised Shear on JET This document is intended for publication in the open literature. It is made available on

More information

proposed [6]. This scaling is found for single null divertor configurations with the VB

proposed [6]. This scaling is found for single null divertor configurations with the VB ORNL/CP-99185 Assessment of Effects of Neutrals on the Power Threshold for L to H Transitions in DIII-D ^ ^ L W Owent, B A Carrerast, R Maingit, P K Mioduszewskit, T N CarlstrfS. R J Groebner* A(jQ j toak

More information

GA A25351 PHYSICS ADVANCES IN THE ITER HYBRID SCENARIO IN DIII-D

GA A25351 PHYSICS ADVANCES IN THE ITER HYBRID SCENARIO IN DIII-D GA A25351 PHYSICS ADVANCES IN THE ITER HYBRID SCENARIO IN DIII-D by C.C. PETTY, P.A. POLITZER, R.J. JAYAKUMAR, T.C. LUCE, M.R. WADE, M.E. AUSTIN, D.P. BRENNAN, T.A. CASPER, M.S. CHU, J.C. DeBOO, E.J. DOYLE,

More information

Influence of Impurity Seeding on ELM Behaviour and Edge Pedestal in ELMy H-Mode Discharges

Influence of Impurity Seeding on ELM Behaviour and Edge Pedestal in ELMy H-Mode Discharges EFDA JET CP()-5 S.Jachmich, G.Maddison, M.N.A.Beurskens, P.Dumortier, T.Eich, A.Messiaen, M.F.F.Nave, J.Ongena, J.Rapp, J.Strachan, M. Stamp, W.Suttrop, G.Telesca, B.Unterberg and JET EFDA Contributors

More information

The Study of Correlation Properties of Geodesic Acoustic Modes in the T-10 Tokamak

The Study of Correlation Properties of Geodesic Acoustic Modes in the T-10 Tokamak 1 The Study of Correlation Properties of Geodesic Acoustic Modes in the T-1 Tokamak A.V. Melnikov 1), L.G. Eliseev 1), S.V. Perfilov 1), S.E. Lysenko 1), V.A. Mavrin 1), R.V. Shurygin 1), D.A. Shelukhin

More information

DPG School The Physics of ITER Physikzentrum Bad Honnef, Energy Transport, Theory (and Experiment) Clemente Angioni

DPG School The Physics of ITER Physikzentrum Bad Honnef, Energy Transport, Theory (and Experiment) Clemente Angioni Max-Planck-Institut für Plasmaphysik DPG School The Physics of ITER Physikzentrum Bad Honnef, 23.09.2014 Energy Transport, Theory (and Experiment) Clemente Angioni Special acknowledgments for material

More information

Rotation and Neoclassical Ripple Transport in ITER

Rotation and Neoclassical Ripple Transport in ITER Rotation and Neoclassical Ripple Transport in ITER Elizabeth J. Paul 1 Matt Landreman 1 Francesca Poli 2 Don Spong 3 Håkan Smith 4 William Dorland 1 1 University of Maryland 2 Princeton Plasma Physics

More information

International Workshop on the Frontiers of Modern Plasma Physics July On the Nature of Plasma Core Turbulence.

International Workshop on the Frontiers of Modern Plasma Physics July On the Nature of Plasma Core Turbulence. 1953-43 International Workshop on the Frontiers of Modern Plasma Physics 14-25 July 2008 On the Nature of Plasma Core Turbulence. F. Jenko Max-Planck Institute fuer Plasmaphysik Garching bei Munchen Germany

More information

Towards the construction of a model to describe the inter-elm evolution of the pedestal on MAST

Towards the construction of a model to describe the inter-elm evolution of the pedestal on MAST Towards the construction of a model to describe the inter-elm evolution of the pedestal on MAST D. Dickinson 1,2, S. Saarelma 2, R. Scannell 2, A. Kirk 2, C.M. Roach 2 and H.R. Wilson 1 June 17, 211 1

More information

Characterization of neo-classical tearing modes in high-performance I- mode plasmas with ICRF mode conversion flow drive on Alcator C-Mod

Characterization of neo-classical tearing modes in high-performance I- mode plasmas with ICRF mode conversion flow drive on Alcator C-Mod 1 EX/P4-22 Characterization of neo-classical tearing modes in high-performance I- mode plasmas with ICRF mode conversion flow drive on Alcator C-Mod Y. Lin, R.S. Granetz, A.E. Hubbard, M.L. Reinke, J.E.

More information

Effect of Biasing on Electron Temperature in IR-T1 Tokamak

Effect of Biasing on Electron Temperature in IR-T1 Tokamak Effect of Biasing on Electron Temperature in IR-T1 Tokamak Sakineh Meshkani 1, Mahmood Ghoranneviss 1 and Mansoureh Lafouti 2 1 Plasma Physics Research Center, Science and Research Branch, Islamic Azad

More information

The role of zonal flows and predator-prey oscillations in triggering

The role of zonal flows and predator-prey oscillations in triggering The role of zonal flows and predator-prey oscillations in triggering the formation of edge and core transport barriers L. Schmitz, 1 L. Zeng, 1 T.L. Rhodes, 1 J.C. Hillesheim 2 W.A. Peebles, 1 R.J. Groebner,

More information

Tungsten impurity transport experiments in Alcator C-Mod to address high priority R&D for ITER

Tungsten impurity transport experiments in Alcator C-Mod to address high priority R&D for ITER Tungsten impurity transport experiments in Alcator C-Mod to address high priority R&D for ITER A. Loarte ITER Organization, Route de Vinon-sur-Verdon, CS 90 046, 13067 St Paul Lez Durance Cedex, France

More information

Multi scale drift turbulence dynamics in an Ohmic discharge as measured at the FT 2 tokamak and modelled by full f gyrokinetic ELMFIRE code

Multi scale drift turbulence dynamics in an Ohmic discharge as measured at the FT 2 tokamak and modelled by full f gyrokinetic ELMFIRE code Multi scale drift turbulence dynamics in an Ohmic discharge as measured at the FT tokamak and modelled by full f gyrokinetic ELMFIRE code E. Gusakov1, V. Bulanin, A. Gurchenko1, J. Heikkinen3, S. Janhunen4,

More information

Internal Transport Barrier Triggering by Rational Magnetic Flux Surfaces in Tokamaks

Internal Transport Barrier Triggering by Rational Magnetic Flux Surfaces in Tokamaks EFDA JET CP(0)07/09 E. Joffrin, C.D. Challis, G.D. Conway, X. Garbet, A. Gude, S. Guenther, N. C. Hawkes, T.C. Hender, D. Howell, G.T.A. Huysmans, E. Lazarro, P. Maget, M. Marachek, A.G. Peeters, S.D.

More information

GA A27433 THE EPED PEDESTAL MODEL: EXTENSIONS, APPLICATION TO ELM-SUPPRESSED REGIMES, AND ITER PREDICTIONS

GA A27433 THE EPED PEDESTAL MODEL: EXTENSIONS, APPLICATION TO ELM-SUPPRESSED REGIMES, AND ITER PREDICTIONS GA A27433 THE EPED PEDESTAL MODEL: EXTENSIONS, APPLICATION TO ELM-SUPPRESSED REGIMES, AND ITER PREDICTIONS by P.B. SNYDER, T.H. OSBORNE, M.N.A. BEURSKENS, K.H. BURRELL, R.J. GROEBNER, J.W. HUGHES, R. MAINGI,

More information

The Physics Basis of ITER Confinement

The Physics Basis of ITER Confinement The Physics Basis of ITER Confinement F. Wagner Max-Planck-Institut für Plasmaphysik EURATOM Association Wendelsteinstr. 1, 17491 Greifswald, Germany Abstract. ITER will be the first fusion reactor and

More information

Tokamak Edge Turbulence background theory and computation

Tokamak Edge Turbulence background theory and computation ASDEX Upgrade Tokamak Edge Turbulence background theory and computation B. Scott Max Planck Institut für Plasmaphysik Euratom Association D-85748 Garching, Germany Krakow, Sep 2006 Outline Basic Concepts

More information

Measuring from electron temperature fluctuations in the Tokamak Fusion Test Reactor

Measuring 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 information

Potentials, E B drifts, and uctuations in the DIII-D boundary

Potentials, E B drifts, and uctuations in the DIII-D boundary Journal of Nuclear Materials 266±269 (1999) 1145±1150 Potentials, E B drifts, and uctuations in the DIII-D boundary R.A. Moyer a, *, R. Lehmer a, J.A. Boedo a, J.G. Watkins b,x.xu c, J.R. Myra d, R. Cohen

More information

On the Nature of ETG Turbulence and Cross-Scale Coupling

On the Nature of ETG Turbulence and Cross-Scale Coupling J. Plasma Fusion Res. SERIES, Vol. Vol. 6 6 (2004) (2004) 11 16 000 000 On the Nature of ETG Turbulence and Cross-Scale Coupling JENKO Frank Max-Planck-Institut für Plasmaphysik, EURATOM-Association, D-85748

More information

PROGRESS IN STEADY-STATE SCENARIO DEVELOPMENT IN THE DIII-D TOKAMAK

PROGRESS IN STEADY-STATE SCENARIO DEVELOPMENT IN THE DIII-D TOKAMAK PROGRESS IN STEADY-STATE SCENARIO DEVELOPMENT IN THE DIII-D TOKAMAK by T.C. LUCE, J.R. FERRON, C.T. HOLCOMB, F. TURCO, P.A. POLITZER, and T.W. PETRIE GA A26981 JANUARY 2011 DISCLAIMER This report was prepared

More information

q(0) pressure after crash 1.0 Single tearing on q=2 Double tearing on q=2 0.5

q(0) pressure after crash 1.0 Single tearing on q=2 Double tearing on q=2 0.5 EX/P-1 MHD issues in Tore Supra steady-state fully non-inductive scenario P Maget 1), F Imbeaux 1), G Giruzzi 1), V S Udintsev ), G T A Huysmans 1), H Lütjens 3), J-L Ségui 1), M Goniche 1), Ph Moreau

More information

TH/P6-14 Integrated particle simulation of neoclassical and turbulence physics in the tokamak pedestal/edge region using XGC a)

TH/P6-14 Integrated particle simulation of neoclassical and turbulence physics in the tokamak pedestal/edge region using XGC a) 1 TH/P6-14 Integrated particle simulation of neoclassical and turbulence physics in the tokamak pedestal/edge region using XGC a) 1 Chang, C.S., 1 Ku, S., 2 Adams M., 3 D Azevedo, G., 4 Chen, Y., 5 Cummings,

More information

Mechanisms for ITB Formation and Control in Alcator C-Mod Identified through Gyrokinetic Simulations of TEM Turbulence

Mechanisms for ITB Formation and Control in Alcator C-Mod Identified through Gyrokinetic Simulations of TEM Turbulence th IAEA Fusion Energy Conference Vilamoura, Portugal, 1-6 November IAEA-CN-116/TH/-1 Mechanisms for ITB Formation and Control in Alcator C-Mod Identified through Gyrokinetic Simulations of TEM Turbulence

More information

Minimal Model Study for ELM Control by Supersonic Molecular Beam Injection and Pellet Injection

Minimal Model Study for ELM Control by Supersonic Molecular Beam Injection and Pellet Injection 25 th Fusion Energy Conference, Saint Petersburg, Russia, 2014 TH/P2-9 Minimal Model Study for ELM Control by Supersonic Molecular Beam Injection and Pellet Injection Tongnyeol Rhee 1,2, J.M. Kwon 1, P.H.

More information

Tokamak Divertor System Concept and the Design for ITER. Chris Stoafer April 14, 2011

Tokamak Divertor System Concept and the Design for ITER. Chris Stoafer April 14, 2011 Tokamak Divertor System Concept and the Design for ITER Chris Stoafer April 14, 2011 Presentation Overview Divertor concept and purpose Divertor physics General design considerations Overview of ITER divertor

More information

W.M. Solomon 1. Presented at the 54th Annual Meeting of the APS Division of Plasma Physics Providence, RI October 29-November 2, 2012

W.M. Solomon 1. Presented at the 54th Annual Meeting of the APS Division of Plasma Physics Providence, RI October 29-November 2, 2012 Impact of Torque and Rotation in High Fusion Performance Plasmas by W.M. Solomon 1 K.H. Burrell 2, R.J. Buttery 2, J.S.deGrassie 2, E.J. Doyle 3, A.M. Garofalo 2, G.L. Jackson 2, T.C. Luce 2, C.C. Petty

More information

L-to-H power threshold comparisons between NBI and RF heated plasmas in NSTX

L-to-H power threshold comparisons between NBI and RF heated plasmas in NSTX Research Supported by L-to-H power threshold comparisons between NBI and RF heated plasmas in NSTX T.M. Biewer 1, R. Maingi 1, H. Meyer 2, R.E. Bell 3, C. Bush 1, S. Kaye 3, S. Kubota 3, B. LeBlanc 3,

More information

Physics of the detached radiative divertor regime in DIII-D

Physics of the detached radiative divertor regime in DIII-D Plasma Phys. Control. Fusion 41 (1999) A345 A355. Printed in the UK PII: S741-3335(99)97299-8 Physics of the detached radiative divertor regime in DIII-D M E Fenstermacher, J Boedo, R C Isler, A W Leonard,

More information

The performance of improved H-modes at ASDEX Upgrade and projection to ITER

The performance of improved H-modes at ASDEX Upgrade and projection to ITER EX/1-1 The performance of improved H-modes at ASDEX Upgrade and projection to George Sips MPI für Plasmaphysik, EURATOM-Association, D-85748, Germany G. Tardini 1, C. Forest 2, O. Gruber 1, P. Mc Carthy

More information

GA A27849 APPLICATION OF ELECTRON CYCLOTRON HEATING TO THE STUDY OF TRANSPORT IN ITER BASELINE SCENARIO-LIKE DISCHARGES IN DIII-D

GA A27849 APPLICATION OF ELECTRON CYCLOTRON HEATING TO THE STUDY OF TRANSPORT IN ITER BASELINE SCENARIO-LIKE DISCHARGES IN DIII-D GA A27849 APPLICATION OF ELECTRON CYCLOTRON HEATING TO THE STUDY OF TRANSPORT IN ITER BASELINE by R.I. PINSKER, M.E. AUSTIN, D.R. ERNST, A.M. GAROFALO, B.A. GRIERSON, J.C. HOSEA, T.C. LUCE, A. MARINONI,

More information

Effects of stellarator transform on sawtooth oscillations in CTH. Jeffrey Herfindal

Effects of stellarator transform on sawtooth oscillations in CTH. Jeffrey Herfindal Effects of stellarator transform on sawtooth oscillations in CTH Jeffrey Herfindal D.A. Ennis, J.D. Hanson, G.J. Hartwell, E.C. Howell, C.A. Johnson, S.F. Knowlton, X. Ma, D.A. Maurer, M.D. Pandya, N.A.

More information

Chapter 1: Overview and summary

Chapter 1: Overview and summary Chapter 1: Overview and summary ITER Physics Basis Editors ITER Physics Expert Group Chairs and Co-Chairs ITER Joint Central Team and Physics Integration Unit ITER EDA, Naka Joint Work Site, Mukouyama,

More information

Confinement and Transport Research in Alcator C-Mod

Confinement and Transport Research in Alcator C-Mod PSFC/JA-05-32. Confinement and Transport Research in Alcator C-Mod M. Greenwald, N. Basse, P. Bonoli, R. Bravenec 1, E. Edlund, D. Ernst, C. Fiore, R. Granetz, A. Hubbard, J. Hughes, I. Hutchinson, J.

More information

Max-Planck-Institut für Plasmaphysik, EURATOM Association POB 1533, D Garching, Germany

Max-Planck-Institut für Plasmaphysik, EURATOM Association POB 1533, D Garching, Germany DEPTH PROFILE REONSTRUTION FROM RUTHERFORD BAKSATTERING DATA U. V. TOUSSAINT, K. KRIEGER, R. FISHER, V. DOSE Max-Planck-Institut für Plasmaphysik, EURATOM Association POB 1533, D-8574 Garching, Germany

More information

Production 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 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 information

A kinetic neutral atom model for tokamak scrape-off layer tubulence simulations. Christoph Wersal, Paolo Ricci, Federico Halpern, Fabio Riva

A kinetic neutral atom model for tokamak scrape-off layer tubulence simulations. Christoph Wersal, Paolo Ricci, Federico Halpern, Fabio Riva A kinetic neutral atom model for tokamak scrape-off layer tubulence simulations Christoph Wersal, Paolo Ricci, Federico Halpern, Fabio Riva CRPP - EPFL SPS Annual Meeting 2014 02.07.2014 CRPP The tokamak

More information

Divertor Detachment on TCV

Divertor Detachment on TCV Divertor Detachment on TCV R. A. Pitts, Association EURATOM-Confédération Suisse,, CH- LAUSANNE, Switzerland thanks to A. Loarte a, B. P. Duval, J.-M. Moret, J. A. Boedo b, L. Chousal b, D. Coster c, G.

More information

Reinstallation of the COMPASS D Tokamak in IPP ASCR

Reinstallation of the COMPASS D Tokamak in IPP ASCR R. Pánek, O. Bilyková, V. Fuchs, M. Hron, P. Chráska, P. Pavlo, J. Stöckel, J. Urban, V. Weinzettl, J. Zajac and F. Žáček Institute of Plasma Physics, Academy of Sciences of the Czech Republic, Za Slovankou

More information

Tungsten impurity transport experiments in Alcator C-Mod to address high priority R&D for ITER

Tungsten impurity transport experiments in Alcator C-Mod to address high priority R&D for ITER Tungsten impurity transport experiments in Alcator C-Mod to address high priority R&D for ITER M.L. Reinke 1, A. Loarte 2, M. Chilenski 3, N. Howard 3, F. Köchl 4, A. Polevoi 2, A. Hubbard 3, J.W. Hughes

More information

Estimating the plasma flow in a recombining plasma from

Estimating the plasma flow in a recombining plasma from Paper P3-38 Estimating the plasma flow in a recombining plasma from the H α emission U. Wenzel a, M. Goto b a Max-Planck-Institut für Plasmaphysik (IPP) b National Institute for Fusion Science, Toki 509-5292,

More information

T (ρ)/t or n(ρ)/n = µ(ρ)i f I (ρ)

T (ρ)/t or n(ρ)/n = µ(ρ)i f I (ρ) PROFILE SHAPE PARAMETERIZATION OF JET ELECTRON TEMPERATURE AND DENSITY PROFILES Beatrix Schunke JET Joint Undertaking, Abingdon, Oxon, OX14 3EA, UK Kaya Imre and Kurt S. Riedel New York University, 251

More information

Evolution of the pedestal on MAST and the implications for ELM power loadings

Evolution of the pedestal on MAST and the implications for ELM power loadings Evolution of the pedestal on MAST and the implications for ELM power loadings Presented by Andrew Kirk EURATOM / UKAEA Fusion Association UKAEA authors were funded jointly by the United Kingdom Engineering

More information

INTRODUCTION TO GYROKINETIC AND FLUID SIMULATIONS OF PLASMA TURBULENCE AND OPPORTUNITES FOR ADVANCED FUSION SIMULATIONS

INTRODUCTION TO GYROKINETIC AND FLUID SIMULATIONS OF PLASMA TURBULENCE AND OPPORTUNITES FOR ADVANCED FUSION SIMULATIONS INTRODUCTION TO GYROKINETIC AND FLUID SIMULATIONS OF PLASMA TURBULENCE AND OPPORTUNITES FOR ADVANCED FUSION SIMULATIONS G.W. Hammett, Princeton Plasma Physics Lab w3.pppl.gov/ hammett Fusion Simulation

More information

Identification of ELM evolution patterns with unsupervised clustering of time-series similarity metrics

Identification of ELM evolution patterns with unsupervised clustering of time-series similarity metrics Identification of ELM evolution patterns with unsupervised clustering of time-series similarity metrics David R. Smith 1, G. McKee 1, R. Fonck 1, A. Diallo 2, S. Kaye 2, B. LeBlanc 2, S. Sabbagh 3, and

More information

A Study of Directly Launched Ion Bernstein Waves in a Tokamak

A Study of Directly Launched Ion Bernstein Waves in a Tokamak PFC-/JA-86-6 A Study of Directly Launched Ion Bernstein Waves in a Tokamak Y. Takase, J. D. Moody, C. L. Fiore, F. S. McDermott, M. Porkolab, and J. Squire Plasma Fusion Center Massachusetts Institute

More information

RESISTIVE BALLOONING MODES AND THE SECOND REGION OF STABILITY

RESISTIVE BALLOONING MODES AND THE SECOND REGION OF STABILITY Plasma Physics and Controlled Fusion, Vol. 29, No. 6, pp. 719 to 121, 1987 Printed in Great Britain 0741-3335/87$3.00+.OO 1OP Publishing Ltd. and Pergamon Journals Ltd. RESISTIVE BALLOONING MODES AND THE

More information

Review of experimental observations of plasma detachment and of the effects of divertor geometry on divertor performance

Review of experimental observations of plasma detachment and of the effects of divertor geometry on divertor performance Review of experimental observations of plasma detachment and of the effects of divertor geometry on divertor performance Alberto Loarte European Fusion Development Agreement Close Support Unit - Garching

More information

Tungsten impurity transport experiments in Alcator C-Mod to address high priority research and development for ITER

Tungsten impurity transport experiments in Alcator C-Mod to address high priority research and development for ITER Tungsten impurity transport experiments in Alcator C-Mod to address high priority research and development for ITER The MIT Faculty has made this article openly available. Please share how this access

More information

Anomalous transport in the H-mode pedestal of Alcator C-Mod discharges

Anomalous transport in the H-mode pedestal of Alcator C-Mod discharges PAPER Anomalous transport in the H-mode pedestal of Alcator C-Mod discharges To cite this article: A.Y. Pankin et al Nucl. Fusion 00 Manuscript version: Accepted Manuscript Accepted Manuscript is the version

More information

Influence of gas impuritites on plasma performance on JET tokamak

Influence of gas impuritites on plasma performance on JET tokamak 19 Influence of gas impuritites on plasma performance on JET tokamak Martin Kubič České vysoké učení technické v Praze Fakulta jaderná a fyzikálně inženýrská Katedra fyziky THERMONUCLEAR FUSION 1.1 Introduction

More information

The International Multi-Tokamak Profile Database

The International Multi-Tokamak Profile Database The International Multi-Tokamak Profile Database The ITER 1D Modelling Working Group: D. Boucher a, J.W. Connor b, W.A. Houlberg c,m.f.turner b, G. Bracco d, A. Chudnovskiy a,e,j.g.cordey b, M.J. Greenwald

More information

Effect of 3D magnetic perturbations on the plasma rotation in ASDEX Upgrade

Effect of 3D magnetic perturbations on the plasma rotation in ASDEX Upgrade Effect of 3D magnetic perturbations on the plasma rotation in ASDEX Upgrade A F Martitsch 1, S V Kasilov 1,2, W Kernbichler 1, G Kapper 1, C G Albert 1, M F Heyn 1, H M Smith 3, E Strumberger 4, S Fietz

More information

Non-linear modeling of the Edge Localized Mode control by Resonant Magnetic Perturbations in ASDEX Upgrade

Non-linear modeling of the Edge Localized Mode control by Resonant Magnetic Perturbations in ASDEX Upgrade 1 TH/P1-26 Non-linear modeling of the Edge Localized Mode control by Resonant Magnetic Perturbations in ASDEX Upgrade F.Orain 1, M.Hölzl 1, E.Viezzer 1, M.Dunne 1, M.Bécoulet 2, P.Cahyna 3, G.T.A.Huijsmans

More information

arxiv: v1 [physics.plasm-ph] 14 Jan 2009

arxiv: v1 [physics.plasm-ph] 14 Jan 2009 arxiv:0901.2043v1 [physics.plasm-ph] 14 Jan 2009 Effectsofparallel ion motion onzonal flowgeneration in ion-temperature-gradient mode turbulence J. Anderson 1, J. Li, Y. Kishimoto Department of Fundamental

More information

Confinement Studies during LHCD and LHW Ion Heating on HL-1M

Confinement Studies during LHCD and LHW Ion Heating on HL-1M Confinement Studies during LHCD and LHW Ion Heating on HL-1M Y. Liu, X.D.Li, E.Y. Wang, J. Rao, Y. Yuan, H. Xia, W.M. Xuan, S.W. Xue, X.T. Ding, G.C Guo, S.K. Yang, J.L. Luo, G.Y Liu, J.E. Zeng, L.F. Xie,

More information

On existence of resistive magnetohydrodynamic equilibria

On existence of resistive magnetohydrodynamic equilibria arxiv:physics/0503077v1 [physics.plasm-ph] 9 Mar 2005 On existence of resistive magnetohydrodynamic equilibria H. Tasso, G. N. Throumoulopoulos Max-Planck-Institut für Plasmaphysik Euratom Association

More information

Comparison of Pellet Injection Measurements with a Pellet Cloud Drift Model on the DIII-D Tokamak

Comparison of Pellet Injection Measurements with a Pellet Cloud Drift Model on the DIII-D Tokamak Comparison of Pellet Injection Measurements with a Pellet Cloud Drift Model on the DIII-D Tokamak T.C. Jernigan, L.R. Baylor, S.K. Combs, W.A. Houlberg (Oak Ridge National Laboratory) P.B. Parks (General

More information

Non-linear MHD Modelling of Rotating Plasma Response to Resonant Magnetic Perturbations.

Non-linear MHD Modelling of Rotating Plasma Response to Resonant Magnetic Perturbations. Non-linear MHD Modelling of Rotating Plasma Response to Resonant Magnetic Perturbations. M. Becoulet 1, F. Orain 1, G.T.A. Huijsmans 2, P. Maget 1, N. Mellet 1, G. Dif-Pradalier 1, G. Latu 1, C. Passeron

More information

ASCOT simulations of electron energy distribution at the divertor targets in an ASDEX Upgrade H-mode discharge

ASCOT simulations of electron energy distribution at the divertor targets in an ASDEX Upgrade H-mode discharge ASCOT simulations of electron energy distribution at the divertor targets in an ASDEX Upgrade H-mode discharge L K Aho-Mantila 1, T Kurki-Suonio 1, A V Chankin 2, D P Coster 2 and S K Sipilä 1 1 Helsinki

More information

Experimental Results on Pellet Injection and MHD from the RTP Tokamak

Experimental Results on Pellet Injection and MHD from the RTP Tokamak Experimental Results on Pellet Injection and MHD from the RTP Tokamak A.A.M. Oomens, J. de Kloe, F.J.B. Salzedas, M.R. de Baar, C.J. Barth, M.N.A. Beurskens, A.J.H. Donné, B. de Groot, G.M.D. Hogeweij,

More information

The RFP: Plasma Confinement with a Reversed Twist

The RFP: Plasma Confinement with a Reversed Twist The RFP: Plasma Confinement with a Reversed Twist JOHN SARFF Department of Physics University of Wisconsin-Madison Invited Tutorial 1997 Meeting APS DPP Pittsburgh Nov. 19, 1997 A tutorial on the Reversed

More information