The European Integrated Tokamak Modelling Effort:Achievements and First Physics Results

This article compares both new and commonly used boundary conditions for generating pressure-driven water flows through carbon nanotubes in molecular dynamics simulations. Three systems are considered: (1) a finite carbon nanotube membrane with streamwise periodicity and ‘gravity’-type Gaussian forcing, (2) a non-periodic finite carbon nanotube membrane with reservoir pressure control, and (3) an infinite carbon nanotube with periodicity and ‘gravity’-type uniform forcing. Comparison between these focuses on the flow behaviour, in particular the mass flow rate and pressure gradient along the carbon nanotube, as well as the radial distribution of water density inside the carbon nanotube. Similar flow behaviour is observed in both membrane systems, with the level of user input required for such simulations found to be largely dependent on the state controllers selected for use in the reservoirs. While System 1 is simple to implement in common molecular dynamics codes, System 2 is more complicated, and the selection of control parameters is less straightforward. A large pressure difference is required between the water reservoirs in these systems to compensate for large pressure losses sustained at the entrance and exit of the nanotube. Despite a simple set-up and a dramatic increase in computational efficiency, the infinite length carbon nanotube in System 3 does not account for these significant inlet and outlet effects, meaning that a much smaller pressure gradient is required to achieve a specified mass flow rate. The infinite tube set-up also restricts natural flow development along the carbon nanotube due to the explicit control of the fluid. Observation of radial density profiles suggests that this results in over-constraint of the water molecules in the tube

Overview of ASDEX Upgrade results

The ASDEX Upgrade programme is directed towards physics input to critical elements of the ITER design and the preparation of ITER operation, as well as addressing physics issues for a future DEMO design. After the finalization of the tungsten coating of the plasma facing components, the re-availability of all flywheel-generators allowed high-power operation with up to 20 MW heating power at I(p) up to 1.2 MA. Implementation of alternative ECRH schemes (140 GHz O2- and X3-mode) facilitated central heating above n(e) = 1.2 x 10(20) m(-3) and low q(95) operation at B(t) = 1.8 T. Central O2-mode heating was successfully used in high P/R discharges with 20 MW total heating power and divertor load control with nitrogen seeding. Improved energy confinement is obtained with nitrogen seeding both for type-I and type-III ELMy conditions. The main contributor is increased plasma temperature, no significant changes in the density profile have been observed. This behaviour may be explained by higher pedestal temperatures caused by ion dilution in combination with a pressure limited pedestal and hollow nitrogen profiles. Core particle transport simulations with gyrokinetic calculations have been benchmarked by dedicated discharges using variations of the ECRH deposition location. The reaction of normalized electron density gradients to variations of temperature gradients and the T(e)/T(i) ratio could be well reproduced. Doppler reflectometry studies at the L-H transition allowed the disentanglement of the interplay between the oscillatory geodesic acoustic modes, turbulent fluctuations and the mean equilibrium E x B flow in the edge negative E(r) well region just inside the separatrix. Improved pedestal diagnostics revealed also a refined picture of the pedestal transport in the fully developed H-mode type-I ELM cycle. Impurity ion transport turned out to be neoclassical in between ELMs. Electron and energy transport remain anomalous, but exhibit different recovery time scales after an ELM. After recovery of the pre-ELM profiles, strong fluctuations develop in the gradients of n(e) and T(e). The occurrence of the next ELM cannot be explained by the local current diffusion time scale, since this turns out to be too short. Fast ion losses induced by shear Alfven eigenmodes have been investigated by time-resolved energy and pitch angle measurements. This allowed the separation of the convective and diffusive loss mechanisms