Direct multiscale coupling of a transport code to gyrokinetic turbulence codes

Direct coupling between a transport solver and local, nonlinear gyrokinetic calculations using the multiscale gyrokinetic code TRINITY [M. Barnes, Ph.D. thesis, arxiv:0901.2868] is described. The coupling of the microscopic and macroscopic physics is done within the framework of multiscale gyrokinetic theory, of which we present the assumptions and key results. An assumption of scale separation in space and time allows for the simulation of turbulence in small regions of the space-time grid, which are embedded in a coarse grid on which the transport equations are implicitly evolved. This leads to a reduction in computational expense of several orders of magnitude, making first-principles simulations of the full fusion device volume over the confinement time feasible on current computing resources. Numerical results from TRINITY simulations are presented and compared with experimental data from JET and ASDEX Upgrade plasmas.Comment: 12 pages, 13 figures, invited paper for 2009 APS-DPP meeting, submitted to Phys. Plasma

Melanoma in Buckinghamshire: Data from the Inception of the Skin Cancer Multidisciplinary Team

Background. Melanoma incidence is increasing faster than any other cancer in the UK. The introduction of specialist skin cancer multidisciplinary teams intends to improve the provision of care to patients suffering from melanoma. This study aims to investigate the management and survival of patients diagnosed with melanoma around the time of inception of the regional skin cancer multidisciplinary team both to benchmark the service against published data and to enable future analysis of the impact of the specialisation of skin cancer care. Methods. All patients diagnosed with primary cutaneous melanoma between January 1, 2003 and December 3, 2005 were identified. Data on clinical and histopathological features, surgical procedures, complications, disease recurrence and 5-year survival were collected and analysed. Results. Two hundred and fourteen patients were included, 134 female and 80 males. Median Breslow thickness was 0.74 mm (0.7 mm female and 0.8 mm male). Overall 5-year survival was 88% (90% female and 85% male). Discussion. Melanoma incidence in Buckinghamshire is in keeping with published data. Basic demographics details concur with classic melanoma distribution and more recent trends, with increased percentage of superficial spreading and thin melanomas, leading to improved survival are reflected

Effects of q(r) on the Alpha Particle Ripple Loss in TFTR

An experiment was done with TFTR DT plasmas to determine the effect of the q(r) profile on the alpha particle ripple loss to the outer midplane. The alpha particle loss measurements were made using a radially movable scintillator detector 20 degrees below the outer midplane. The experimental results were compared with TF ripple loss calculations done using a Monte Carlo guiding center orbit following code, ORBIT. Although some of the experimental results are consistent with the ORBIT code modeling, the variation of the alpha loss with the q(r) profiles is not well explained by this code. Quantitative interpretation of these measurements requires a careful analysis of the limiter shadowing effect, which strongly determines the diffusion of alphas into the detector aperture

Global Confinement, Sawtooth Mixing, and Stochastic Diffusion Ripple Loss of Fast ICRF-driven H+ Minority Ions in TFTR

This paper presents studies of ICRF-driven H+ minority ions in TFTR (Tokamak Fusion Test Reator) deuterium plasmas using primarily passive Ho flux detection in the energy range of 0.2-1.0 MeV with some corroborating active (lithium pellet charge exchange) measurements. It is shown that in the passive mode the main donors for the neutralization of H+ ions in this energy range are C5+ ions. The measured effective H+ tail temperatures range from 0.15 MeV at an ICRF power of 2 MW to 0.35 MeV at 6 MW. Analysis of the ICRF-driven H+ ion energy balance has been performed on the basis of the dependence of effective H+ temperatures on the plasma parameters. The analysis showed that H+ confinement times are comparable with their slowing-down times and tended to decrease with increasing ICRF power. Radial redistribution of ICRF-driven H+ ions was detected when giant sawtooth crashes occurred during the ICRF heating. The redistribution affected ions with energy below 0.7-0.8 MeV. The sawtooth crashes displace H+ ions outward along the plasma major radius into the stochastic ripple diffusion domain were those ions are lost in about 10 milliseconds. These observations are consistent with the model of the redistribution of energetic particles developed previously to explain the results of deuterium-tritium alpha-particle redistribution due to sawteeth observed in TFTR. The experimental data are also consistent with ORBIT code simulations of H+ stochastic ripple diffusion losses

Effect of plasma rotation on sawtooth stabilization by beam ions

The sawtooth period in JET ELM-free H-Mode plasmas is increasing with Neutral Beam Injection (NBI) power. For injected power PNBI 12MW no large sawtooth crash is observed during the ELM-free period. However, as the edge stability is improved and external kink modes and ELMs are delayed, a possible sawtooth crash at a high plasma beta becomes a concern. In JET DT experiments, delaying sawteeth was found to be crucial in the quest for high fusion power. Fast particles are known to provide stabilizing effect on sawteeth, however, sawtooth stabilization by NBI ions is not clearly understood, since NBI ions are usually not ''fast'' enough to stabilize the m/n = 1/1 internal kink mode which is believed to cause the crash. In order to understand the observed sawteeth stabilization in tokamak experiments with NBI heating, the internal kink m/n = 1/1 mode stability of JET plasmas was modeled using the NOVA-K code, which is also benchmarked with the nonperturbative version of NOVA and the M3D code. Comparison of m/n = 1/1 mode stabilization by NBI ions in JET and TFTR and application of the nonlinear stabilization criteria is given