J.A. Leuer

Real time equilibrium reconstruction for tokamak discharge control

A practical method for performing a tokamak equilibrium reconstruction in real time for arbitrary time varying discharge shapes and current profiles is described. An approximate solution to the Grad-Shafranov equilibrium relation is found which best fits the diagnostic measurements. Thus, a solution for the spatial distribution of poloidal flux and toroidal current density is available in real time that is consistent with plasma force balance, allowing accurate evaluation of parameters such as discharge shape and safety factor profile. The equilibrium solutions are produced at a rate sufficient for discharge control. This equilibrium reconstruction algorithm has been implemented on the digital plasma control system for the DIII-D tokamak. The first application of real time equilibrium reconstruction to discharge shape control is described.

Experimental vertical stability studies for ITER performance and design guidance

United States Department of Energy (DE-FC02-04ER54698, DEAC52- 07NA27344, and DE-FG02-04ER54235)

Development of ITER 15 MA ELMy H-mode inductive scenario

The poloidal field (PF) coil system on ITER, which provides both feedforward and feedback control of plasma position, shape, and current, is a critical element for achieving mission performance. Analysis of PF capabilities has focused on the 15 MA Q = 10 scenario with a 300-500 s flattop burn phase. The operating space available for the 15 MA ELMy H-mode plasma discharges in ITER and upgrades to the PF coils or associated systems to establish confidence that ITER mission objectives can be reached have been identified. Time dependent self-consistent free-boundary calculations were performed to examine the impact of plasma variability, discharge programming, and plasma disturbances. Based on these calculations a new reference scenario was developed based upon a large bore initial plasma, early divertor transition, low level heating in L-mode, and a late H-mode onset. Equilibrium analyses for this scenario indicate that the original PF coil limitations do not allow low li (<0.8) operation or lower flux states, and the flattop burn durations were predicted to be less than the desired 400 s. This finding motivates the expansion of the operating space, considering several upgrade options to the PF coils. Analysis was also carried out to examine the feedback current reserve required in the CS and PF coils during a series of disturbances and a feasibility assessment of the 17 MA scenario was undertaken. Results of the studies show that the new scenario and modified PF system will allow a wide range of 15 MA 300-500 s operation and more limited but finite 17 MA operation.

Anomalies in the applied magnetic fields in DIII-D and their implications for the understanding of stability experiments

Small non-axisymmetric magnetic fields are known to cause serious loss of stability in tokamaks, leading to loss of confinement and abrupt termination of plasma current (disruptions). The best known examples are the locked mode and the resistive wall mode. Understanding of the underlying field anomalies (departures in the hardware-related fields from ideal toroidal and poloidal fields on a single axis) and the interaction of the plasma with them is crucial to tokamak development. Results of both locked mode experiments (Scoville J.T. and La Haye R.J. 2003 Nucl. Fusion 43 250) and resistive wall mode experiments (Garofalo A.M., La Haye R.J. and Scoville J.T. 2002 Nucl. Fusion 42 1335) done in DIII-D tokamak plasmas have been interpreted to indicate the presence of a significant anomalous field. New measurements of the magnetic field anomalies of the hardware systems have been made in DIII-D. The measured field anomalies due to the plasma shaping coils in DIII-D are smaller than previously reported (La Haye R.J. and Scoville J.T. 1991 Rev. Sci. Instrum. 61 2146). Additional evaluations of systematic errors have been made. New measurements of the anomalous fields of the Ohmic heating and toroidal coils have been added. Such detailed in situ measurements of the fields of a tokamak are unique. The anomalous fields from all the coils are one-third the values indicated from the stability experiments (Garofalo et al 2002, Scoville and La Haye 2003). These results indicate limitations in the understanding of the interaction of the plasma with the external field. They indicate that it may not be possible to deduce the anomalous fields in a tokamak from plasma experiments and that we may not have the basis needed to project the error field requirements of future tokamaks.

Development of ITER-relevant plasma control solutions at DIII-D

The requirements of the DIII-D physics program have led to the development of many operational control results with direct relevance to ITER. These include new algorithms for robust and sustained stabilization of neoclassical tearing modes with electron cyclotron current drive, model-based controllers for stabilization of the resistive wall mode in the presence of ELMs, coupled linear–nonlinear algorithms to provide good dynamic axisymmetric control while avoiding coil current limits, and adaptation of the DIII-D plasma control system (PCS) to operate next-generation superconducting tokamaks. Development of integrated plasma control (IPC), a systematic approach to modelbased design and controller verification, has enabled successful experimental application of high reliability control algorithms requiring a minimum of machine operations time for testing and tuning. The DIII-D PCS hardware and software and its versions adapted for other devices can be connected to IPC simulations to confirm control function prior to experimental use. This capability has been important in control system implementation for tokamaks under construction and is expected to be critical for ITER.

Fusion Nuclear Science Facility Candidates

Abstract To move to a fusion DEMO power plant after ITER, a Fusion Nuclear Science Facility (FNSF) is needed in addition to ITER and research in operating tokamaks and those under construction. The FNSF will enable research on how to utilize and deal with the products of fusion reactions, addressing such issues as how to extract the energy from neutrons and alpha particles into high-temperature process heat streams to be either used directly or converted to electricity, how to make tritium from the neutrons and lithium, how to deal with the effects of the neutrons on the blanket structures, and how to manage the first wall surface erosion caused by the alpha particle heat appearing as low-energy plasma fluxes to those surfaces. Two candidates for the FNSF are considered in this paper: normal and low aspect ratio copper magnet tokamaks. The methods of selecting optimum machine design points versus aspect ratio are fully presented. The two options are compared and contrasted; both options appear viable.

Evaluation of CFETR as a Fusion Nuclear Science Facility using multiple system codes

This paper presents the results of a multi-system codes benchmarking study of the recently published China Fusion Engineering Test Reactor (CFETR) pre-conceptual design (Wan et al 2014 IEEE Trans. Plasma Sci. 42 495). Two system codes, General Atomics System Code (GASC) and Tokamak Energy System Code (TESC), using different methodologies to arrive at CFETR performance parameters under the same CFETR constraints show that the correlation between the physics performance and the fusion performance is consistent, and the computed parameters are in good agreement. Optimization of the first wall surface for tritium breeding and the minimization of the machine size are highly compatible. Variations of the plasma currents and profiles lead to changes in the required normalized physics performance, however, they do not significantly affect the optimized size of the machine. GASC and TESC have also been used to explore a lower aspect ratio, larger volume plasma taking advantage of the engineering flexibility in the CFETR design. Assuming the ITER steady-state scenario physics, the larger plasma together with a moderately higher BT and Ip can result in a high gain Qfus ~ 12, Pfus ~ 1 GW machine approaching DEMO-like performance. It is concluded that the CFETR baseline mode can meet the minimum goal of the Fusion Nuclear Science Facility (FNSF) mission and advanced physics will enable it to address comprehensively the outstanding critical technology gaps on the path to a demonstration reactor (DEMO). Before proceeding with CFETR construction steady-state operation has to be demonstrated, further development is needed to solve the divertor heat load issue, and blankets have to be designed with tritium breeding ratio (TBR) >1 as a target.

KSTAR equilibrium operating space and projected stabilization at high normalized beta

Along with an expanded evaluation of the equilibrium operating space of the Korea Superconducting Tokamak Advanced Research, KSTAR, experimental equilibria of the most recent plasma discharges were reconstructed using the EFIT code. In near-circular plasmas created in 2009, equilibria reached a stored energy of 54kJ with a maximum plasma current of 0.34MA. Highly shaped plasmas with near double-null configuration in 2010 achieved H-mode with clear edge localized mode (ELM) activity, and transiently reached a stored energy of up to 257kJ, elongation of 1.96 and normalized beta of 1.3. The plasma current reached 0.7MA. Projecting active and passive stabilization of global MHD instabilities for operation above the ideal no-wall beta limit using the designed control hardware was also considered. Kinetic modification of the ideal MHD n = 1 stability criterion was computed by the MISK code on KSTAR theoretical equilibria with a plasma current of 2MA, internal inductance of 0.7 and normalizedbetaof4.0withsimpledensity,temperatureandrotationprofiles. Thesteepedgepressuregradientofthis equilibrium resulted in the need for significant plasma toroidal rotation to allow thermal particle kinetic resonances to stabilize the resistive wall mode (RWM). The impact of various materials and electrical connections of the passive stabilizing plates on RWM growth rates was analysed, and copper plates reduced the RWM passive growth rate by a factor of 15 compared with stainless steel plates at a normalized beta of 4.4. Computations of active RWM control using the VALEN code showed that the n = 1 mode can be stabilized at normalized beta near the ideal wall limit via control fields produced by the midplane in-vessel control coils (IVCCs) with as low as 0.83kW control power using ideal control system assumptions. The ELM mitigation potential of the IVCC, examined by evaluating the vacuum island overlap created by resonant magnetic perturbations, was analysed using the TRIP3D code. Using a combinationofallIVCCswithdominant n = 2fieldandupper/lowercoilsinanevenparityconfiguration,aChirikov parameter near unity at normalized poloidal flux 0.83, an empirically determined condition for ELM mitigation in DIII-D, was generated in theoretical high-beta equilibria. Chirikov profile optimization was addressed in terms of coil parity and safety factor profile. (Some figures in this article are in colour only in the electronic version)

Status of DIII-D plasma control

A key component of the DIII-D Advanced Tokamak and Radiative Divertor Programs is the development and implementation of methods to actively control a large number of plasma parameters. These parameters include plasma shape and position, total stored energy, density, RF loading resistance, radiated power and more detailed control of the current profile. To support this research goal, a flexible and easily expanded digital control system has been developed and implemented. We have made parallel progress in modeling of the plasma, poloidal coils, vacuum vessel, and power system dynamics and in ensuring the integrity of diagnostic and command circuits used in control. Recent activity has focused on exploiting the mature digital control platform through the implementation of simple feedback controls of more exotic plasma parameters such as enhanced divertor radiation, neutral pressure and Marfe creation and more sophisticated identification and digital feedback control algorithms for plasma shape, vertical position, and safety factor on axis (q/sub 0/). A summary of recent progress in each of these areas will be presented.

Physics Basis of a Fusion Development Facility Utilizing the Tokamak Approach

Abstract The objective of the Fusion Development Facility (FDF) under consideration is to carry forward advanced tokamak physics for optimization of fusion reactors and enable development of fusion’s energy applications. A concept of FDF based on the tokamak approach with conservative expressions of advanced physics and nonsuperconducting magnet technology is presented. It is envisioned to nominally provide 2 MW/m2 of neutron wall loading and operate continuously for up to 2 weeks as required for fusion nuclear component research and development. FDF will have tritium breeding capability with a goal of addressing the tritium self-sufficiency issue for fusion energy. A zero-dimensional system study using extrapolations of current physics and technology is used to optimize FDF for reasonable power consumption and moderate size. It projects a device that is between the DIII-D tokamak (major radius 1.8 m) [J. L. Luxon, Nucl. Fusion, Vol. 42, p. 614 (2002)] and the Joint European Torus (major radius 3 m) [P. H. Rebut, R. J. Bickerton, and B. E. Keen, Nucl. Fusion, Vol. 25, p. 1011 (1985)] in size, with an aspect ratio A of 3.5 and a fusion gain Q of 2 to 5. Theory-based stability and transport modeling is used to complement the system study and to address physics issues related to specific design points. It is demonstrated that the FDF magnetohydrodynamic stability limits can be readily met with conservative stabilizing conducting wall placement. Transport analysis using a drift-wave-based model with an edge boundary condition consistent with the pedestal stability limit indicates that the FDF confinement requirement can also be readily satisfied. A surprising finding is that the toroidal Alfvén eigenmodes are stabilized by strong ion Landau damping. Analysis of vertical stability control indicates that the basis configuration with an elongation κx [approximately] 2.35 can be controlled using a power supply technology similar to that used in DIII-D. Peak heat fluxes to the divertor are somewhat lower than those of ITER [R. Aymar, P. Barabaschi, and Y. Shimomura, Plasma Phys. Control. Fusion, Vol. 44, p. 519 (2002)], but FDF will operate with a higher duty factor.