The non-linear reduced four-field RMHD model in cylindrical geometry was extended to include plasma rotation, neoclassical poloidal viscosity and two fluid diamagnetic effects. Interaction of the static resonant magnetic perturbations (RMPs) with the rotating plasmas in tokamaks was studied. The self-consistent evolution of equilibrium electric field due to RMP penetration is taken into account in the model. It is demonstrated that in the pedestal region with steep pressure gradients, mean flows perpendicular to the magnetic field, which includes E × B and electron diamagnetic components plays an essential role in RMP screening by plasma. Generally, the screening effect increases for lower resistivity, stronger rotation and smaller RMP amplitude. Strong screening of central islands was observed limiting RMP penetration to the narrow region near the separatrix. However, at certain plasma parameters and due to the non-linear evolution of the radial electric field produced by RMPs, the E × B rotation can be compensated by electron diamagnetic rotation locally. In this case, RMPs can penetrate and form magnetic islands. Typical plasma parameters and RMPs spectra on DIII-D, JET and ITER were used in modelling examples presented in the paper.

In the perspective of operating tungsten monoblocks in WEST, the ongoing major upgrade of the Tore Supra tokamak, a dedicated modelling effort has been carried out to simulate the interaction between the edge plasma and the tungsten wall. A new transport code, SolEdge2D-EIRENE, has been developed with the ability to simulate the plasma up to the first wall. This is especially important for steady state operation, where thermal loads on all the plasma facing components, even remote from the plasma, are of interest. Moreover, main chamber tungsten sources are thought to dominate the contamination of the plasma core. We present here in particular new developments aimed at improving the description of the interface between the plasma and the wall, namely a way to treat sheath physics in a more faithful way using the output of 1D particle in cell simulations. Moreover, different models for prompt redeposition have been implemented and are compared. The latter is shown to play an important role in the balance between divertor and main chamber sources.

An integrated model capable of self-consistent Ion Cyclotron Resonant Heating (ICRH) simulations has been developed. This model includes both full shaping and pressure effects, warm contributions to the dielectric tensor, pressure anisotropy and finite orbit width. It evolves the equilibrium, wave field and full hot particle distribution function until a self-consistent solution is found. This article describes the workings of the three codes VMEC, LEMan and VENUS and how they are linked for iterated computations in a code package we have named SCENIC. The package is thoroughly tested and it is demonstrated that a number of iterations have to be performed in order to find a consistent solution. Since the formulation of the problem can treat general 3D systems, we show a quasi-axisymmetric stellarator low power test case, and then concentrate on experimentally relevant Joint European Torus (JET) 2D configurations.