Abstract.A selection of achievements and first physics results are presented of the European Integrated Tokamak Modelling Task Force (EFDA ITM-TF) simulation framework, which aims to provide a standardized platform and an integrated modelling suite of validated numerical codes for the simulation and prediction of a complete plasma discharge of an arbitrary tokamak. The framework developed by the ITM-TF, based on a generic data structure including both simulated and experimental data, allows for the development of sophisticated integrated simulations (workflows) for physics application. The equilibrium reconstruction and linear MHD stability simulation chain was applied, in particular, to the analysis of the edge MHD stability of ASDEX Upgrade type-I ELMy Hmode discharges and ITER hybrid scenario, demonstrating the stabilizing effect of an increased Shafranov shift on edge modes. Interpretive simulations of a JET hybrid discharge were performed with two electromagnetic turbulence codes within ITM infrastructure showing the signature of trapped-electron assisted ITG turbulence. A successful benchmark among five EC beam/ray-tracing codes was performed in the ITM framework for an ITER inductive scenario for different launching conditions from the Equatorial and Upper Launcher, showing good agreement of the computed * See the Appendix.

The first results of the Dynamic Ergodic Divertor in TEXTOR, when operating in the m=n 3=1 mode configuration, are presented. The deeply penetrating external magnetic field perturbation of this configuration increases the toroidal plasma rotation. Staying below the excitation threshold for the m=n 2=1 tearing mode, this toroidal rotation is always in the direction of the plasma current, even if the toroidal projection of the rotating magnetic field perturbation is in the opposite direction. The observed toroidal rotation direction is consistent with a radial electric field, generated by an enhanced electron transport in the ergodic layers near the resonances of the perturbation. This is an effect different from theoretical predictions, which assume a direct coupling between rotating perturbation and plasma to be the dominant effect of momentum transfer. Helical magnetic field perturbations are introduced in tokamak plasmas to study, on the one hand, the ergodic divertor concept [1,2] and, on the other hand, the interaction of such perturbations with the magnetohydrodynamics (MHD) stability of the plasma [3,4]. Recent experiments, for instance, suggest a control method to mitigate edge localized modes while maintaining the pedestal pressure and thus plasma confinement [5][6][7]. However, open questions remain, in particular, with regard to the influence on the momentum transport of the plasma. Indeed, one motivation to equip the tokamak TEXTOR with the Dynamic Ergodic Divertor (DED) [8] was to be able to study the interaction between helical magnetic field perturbations and plasma transport and stability.The DED consists of 16 magnetic perturbation coils (four quadruples), plus two additional coils for the compensation of the magnetic field imperfections at the feeder regions of the coils. The coils wind helically around the inner side of the torus (major radius: R 1:75 m; minor radius of the circular plasma cross section typically a 0:47 m) with a pitch corresponding to the magnetic field lines of the magnetic flux surface with a safety factor of q 3. Depending on the choice of coil connections to the power supplies, base modes with different poloidal and toroidal mode numbers can be produced. For the DED these are m=n 12=4, 6=2, and 3=1. The penetration depth into the plasma strongly depends on the mode numbers: While the m=n 12=4 affects the edge plasma only, the m=n 3=1 mode reaches into the plasma center (the maximum radial magnetic field component achievable by the DED at the q 2 surface is 10 ÿ3 of the total magnetic field).In this Letter we present results obtained by the m=n 3=1 mode operation. Covering about one-third of the poloidal cross section of the torus, the mode spectrum of the DED does not contain many sidebands. For the m=n 3=1 configuration the three dominant resonant components inside the plasma are m 1, 2, and 3. In Fig. 1 their strengths at the respective resonances are PRL 94, 015003 (2005) P H Y S I C A L

The "European Transport Simulator" (ETS) [1,2] is the new modular package for 1-D discharge evolution developed within the EFDA Integrated Tokamak Modelling (ITM) Task Force. It consists of precompiled physics modules combined into a workflow through standardized input/output data-structures. Ultimately, the ETS will allow for an entire discharge simulation from the start up until the current termination phase, including controllers and sub-systems. The paper presents the current status of the ETS towards this ultimate goal. It discusses the design of the workflow, the validation and verification of its components on the example of impurity solver and demonstrates a proof-of-principles coupling of a local gyrofluid model for turbulent transport to the ETS. It also presents first results on the application of the ETS to JET tokamak discharges with ITER like wall. It studies the correlations of the radiation from impurity to the choice of the sources and transport coefficients. ETS transport solver and workflow designThe ETS adopts a modular approach, where standalone precompiled physics modules ('actors' in the context of the workflow) are coupled into the workflow through standardised interfaces linked with the ITM data-structure. In view of allowing collective development of

Since the installation of an ITER-like wall, the JET programme has focused on the consolidation of ITER design choices and the preparation for ITER operation, with a specific emphasis given to the bulk tungsten melt experiment, which has been crucial for the final decision on the material choice for the day-one tungsten divertor in ITER. Integrated scenarios have been progressed with the re-establishment of long-pulse, high-confinement H-modes by optimizing the magnetic configuration and the use of ICRH to avoid tungsten impurity accumulation. Stationary discharges with detached divertor conditions and small edge localized modes have been demonstrated by nitrogen seeding. The differences in confinement and pedestal behaviour before and after the ITER-like wall installation have been better characterized towards the development of high fusion yield scenarios in DT. Post-mortem analyses of the plasma-facing components have confirmed the previously reported low fuel retention obtained by gas balance and shown that the pattern of deposition within the divertor has changed significantly with respect to the JET carbon wall campaigns due to the absence of thermally activated chemical erosion of beryllium in contrast to carbon. Transport to remote areas is almost absent and two orders of magnitude less material is found in the divertor.

The European Integrated Tokamak Modelling Task Force (ITM-TF) is developing a new type of fully modular and flexible integrated tokamak simulator, which will allow a large variety of simulation types. This ambitious goal requires new concepts of data structure and workflow organisation, which are described for the first time in this paper. The backbone of the system is a physics-and workflow-oriented data structure which allows for the deployment of a fully modular and flexible workflow organisation. The data structure is designed to be generic for any tokamak device and can be used to address physics simulation results, experimental data (including description of subsystem hardware) and engineering issues.

The status of the European Transport Solver, a new 1-D core transport code that is being developed by the members of Integrated Modelling Project 3 ("Transport Code and Discharge Evolution") of the EFDA Task Force on Integrated Tokamak Modelling (ITM), is described. The approach taken by the ITM is to couple codes so that the only exchange is via well-specified data structures (Consistent Physical Objects), with the aim of having the workflow managed by Kepler, a scientific workflow engine.

The results of predictive self-consistent modelling of plasma parameters in low (L) and high (H) confinement modes by the one-dimensional transport code RITM, with particular emphasis on the properties of the edge transport barrier, are presented and discussed. The same transport model is used under both Land H-mode conditions and includes contributions from ion temperature gradient (ITG), trapped electron, drift Alfvén (DA) and drift resistive ballooning instabilities described in the fluid approximation. The computations predict the formation of the edge transport barrier at a high enough heating power due to stabilization of ITG and DA modes, dominating the edge transport in the L-mode, through the effects of the density gradient and the pressure gradient at low collisionality, respectively. The calculated radial profiles and scalings for pedestal and confinement characteristics are compared with measurements on JET, DIII-D and JT-60U tokamaks.

Confinement deterioration and rollover from the radiative improved (RI) mode to the low (L) confinement mode in plasmas with a strong gas puff in the tokamak TEXTOR-94 are modelled by the one-dimensional transport code RITM (radiation from impurities in the transport model). The anomalous transport coefficients in the plasma core include contributions from the ion temperature gradient and dissipative trapped electron instabilities. This model allows us to reproduce the L-RI bifurcation with impurity seeding in good agreement with experimental observations. Whereas the transport at the plasma edge under the RI-mode conditions might be described by the electrostatic turbulence caused by electric currents in the scrape-off layer of the limiter, the present computations show that the level of the edge transport must be increased by at least a factor of three in order to reproduce the evolution of the plasma density profile and the effective ion charge during RI-to-L back transitions. An increase of this order is in agreement with reflectometer measurements and could tentatively be explained by the effect of neutrals on resistive drift modes.