Steven Anthony Sabbagh

Simulations of deuterium-tritium experiments in TFTR

A transport code (TRANSP) is used to simulate future deuterium-tritium (DT) experiments in TFTR. The simulations are derived from 14 TFTR DD discharges, and the modelling of one supershot is discussed in detail to indicate the degree of accuracy of the TRANSP modelling. Fusion energy yields and alpha particle parameters are calculated, including profiles of the alpha slowing down time, the alpha average energy, and the Alfven speed and frequency. Two types of simulation are discussed. The main emphasis is on the DT equivalent, where an equal mix of D and T is substituted for the D in the initial target plasma, and for the D0 in the neutral beam injection, but the other measured beam and plasma parameters are unchanged. This simulation does not assume that alpha heating will enhance the plasma parameters or that confinement will increase with the addition of tritium. The maximum relative fusion yield calculated for these simulations is QDT ~ 0.3, and the maximum alpha contribution to the central toroidal β is βα(0) ~ 0.5%. The stability of toroidicity induced Alfven eigenmodes (TAE) and kinetic ballooning modes (KBM) is discussed. The TAE mode is predicted to become unstable for some of the simulations, particularly after the termination of neutral beam injection. In the second type of simulation, empirical supershot scaling relations are used to project the performance at the maximum expected beam power. The MHD stability of the simulations is discussed

Equilibrium properties of spherical torus plasmas in NSTX

Research in NSTX has been conducted to establish spherical torus plasmas to be used for high ?, auxiliary heated experiments. This device has a major radius R0 = 0.86?m and a midplane halfwidth of 0.7?m. It has been operated with toroidal magnetic field B0 ? 0.3?T and Ip ? 1.0?MA. The evolution of the plasma equilibrium is analysed between discharges with an automated version of the EFIT code. Limiter, double null and lower single null diverted configurations have been sustained for several energy confinement times. The plasma stored energy reached 92?kJ (?t = 17.8%) with neutral beam heating. A plasma elongation in the range 1.6 ? ? ? 2.0 and a triangularity in the range 0.25 ? ? ? 0.45 have been sustained, with values of ? = 2.6 and ? = 0.6 being reached transiently. The reconstructed magnetic signals are fitted to the corresponding measured values with low errors. Aspects of the plasma boundary, pressure and safety factor profiles are supported by measurements from non-magnetic diagnostics. Plasma densities have reached 0.8 and 1.2 times the Greenwald limit in deuterium and helium plasmas, respectively, with no clear limit encountered. Instabilities including sawteeth and reconnection events, characterized by Mirnov oscillations, and a perturbation of the Ip, ? and li evolutions, have been observed. A low q limit was observed and is imposed by a low toroidal mode number kink instability.

Advances in global MHD mode stabilization research on NSTX

Stabilizing modes that limit plasma beta and reduce their deleterious effect on plasma rotation are key goals for the efficient operation of a fusion reactor. Passive stabilization and active control of global kink/ballooning modes and resistive wall modes (RWMs) have been demonstrated on NSTX and research is now advancing towards understanding the stabilization physics and reliably maintaining the high beta plasma for confident extrapolation to ITER and a fusion component test facility based on the spherical torus. Active n = 1 control experiments with an expanded sensor set, combined with low levels of n = 3 field phased to reduce error fields, reduced resonant field amplification and maintained plasma rotation, exceeded normalized beta = 6 and produced record discharge durations limited by magnet system constraints. Details of the observed RWM dynamics during active control show the mode being converted to a rotating kink that stabilizes or saturates and may lead to tearing modes. Discharges with rotation reduced by n = 3 magnetic braking suffer beta collapse at normalized beta = 4.2 approaching the no-wall limit, while normalized beta greater than 5.5 has been reached in these plasmas with n = 1 active control, in agreement with the single-mode RWM theory. Advanced state-space control algorithms proposed for RWM control in ITER theoretically yield significant stabilization improvements. Values of relative phase between the measured n = 1 mode and the applied correction field that experimentally produce stability/instability agree with RWM control modelling. Experimental mode destabilization occurs over a large range of plasma rotation, challenging the notion of a simple scalar critical rotation speed defining marginal stability. Stability calculations including kinetic modifications to the ideal MHD theory are applied to marginally stable experimental equilibria. Plasma rotation and collisionality variations are examined in the calculations. Intermediate rotation levels are less stable, consistent with experimental observations. Trapped ion resonances play a key role in this result. Recent experiments have demonstrated magnetic braking by non-resonant n = 2 fields. The observed rotation damping profile is broader than found for n = 3 fields. Increased ion temperature in the region of maximum braking torque increases the observed rate of rotation damping, consistent with the theory of neoclassical toroidal viscosity at low collisionality.

Cross-machine comparison of resonant field amplification and resistive wall mode stabilization by plasma rotation

Dedicated experiments in the DIII-D tokamak [J. L. Luxon, Nucl. Fusion, 42, 614 (2002)], the Joint European Torus (JET) [P. H. Rebut, R. J. Bickerton, and B. E. Keen, Nucl. Fusion 25, 1011 (1985)], and the National Spherical Torus Experiment (NSTX) [M. Ono, S. M. Kaye, Y.-K. M. Peng et al., Nucl. Fusion 40, 557 (2000)] reveal the commonalities of resistive wall mode (RWM) stabilization by sufficiently fast toroidal plasma rotation in devices of different size and aspect ratio. In each device the weakly damped n=1 RWM manifests itself by resonant field amplification (RFA) of externally applied n=1 magnetic fields, which increases with the plasma pressure. Probing DIII-D and JET plasmas with similar ideal magnetohydrodynamic (MHD) stability properties with externally applied magnetic n=1 fields, shows that the resulting RFA is independent of the machine size. In each device the drag resulting from RFA slows the toroidal plasma rotation and can lead to the onset of an unstable RWM. The critical plasma rotati...

Physics design of the national spherical torus experiment

The mission of the National Spherical Torus Experiment (NSTX) is to prove the principles of spherical torus physics by producing high-beta toroidal plasmas that are non-inductively sustained, and whose current profiles are in steady-state. NSTX will be one of the first ultra low a[P(input) up to 11 MW] in order to produce high-beta toroidal (25 to 40%), low collisionality, high bootstrap fraction (less than or equal to 70%) discharges. Both radio-frequency (RF) and neutral-beam (NB) heating and current drive will be employed. Built into NSTX is sufficient configurational flexibility to study a range of operating space and the resulting dependences of the confinement, micro- and MHD stability, and particle and power handling properties. NSTX research will be carried out by a nationally based science team.

An approximate analytic expression for neoclassical toroidal plasma viscosity in tokamaks

An approximate analytic expression for neoclassical toroidal plasma viscosity in tokamaks that have error fields or magnetohydrodynamic activities is presented. The expression smoothly joins transport fluxes or plasma viscosity in all the known collisionality regimes derived from the solution of the bounce averaged drift kinetic equation and should be useful in modelling results of existing and future tokamak experiments. It also incorporates some of the extensions of the known expressions to include the effects of finite ∇B drift in the non-resonant transport processes. Here, B is the magnitude of the magnetic field. The toroidal momentum balance equation is a nonlinear function of the radial electric field when the neoclassical plasma viscosity is dominant. It can have bifurcated solutions for the radial electric field and may lead to better plasma confinement as a result.

Resistive wall stabilized operation in rotating high beta NSTX plasmas

The National Spherical Torus Experiment (NSTX) has demonstrated the advantages of low aspect ratio geometry in accessing high toroidal and normalized plasma beta, and βN ≡ 10 8〈βt〉 aB0/Ip. Experiments have reached βt = 39% and βN = 7.2 through boundary and profile optimization. High βN plasmas can exceed the ideal no-wall stability limit, βNno-wall, for periods much greater than the wall eddy current decay time. Resistive wall mode (RWM) physics is studied to understand mode stabilization in these plasmas. The toroidal mode spectrum of unstable RWMs has been measured with mode number n up to 3. The critical rotation frequency of Bondeson-Chu, Ωcrit = ωA/(4q2), describes well the RWM stability of NSTX plasmas when applied over the entire rotation profile and in conjunction with the ideal stability criterion. Rotation damping and global rotation collapse observed in plasmas exceeding βNno-wall differs from the damping observed during tearing mode activity and can be described qualitatively by drag due to neoclassical toroidal viscosity in the helically perturbed field of an ideal displacement. Resonant field amplification of an applied n = 1 field perturbation has been measured and increases with increasing βN. Equilibria are reconstructed including measured ion and electron pressure, toroidal rotation and flux isotherm constraint in plasmas with core rotation ω/ωA up to 0.48. Peak pressure shifts of 18% of the minor radius from the magnetic axis have been reconstructed.

The role of kinetic effects, including plasma rotation and energetic particles, in resistive wall mode stability

particles, in resistive wall mode stability J.W. Berkery1, S.A. Sabbagh1, R. Betti2, B. Hu2, R.E. Bell3, S.P. Gerhardt3, J. Manickam3, and M. Podesta3 1Department of Applied Physics and Applied Mathematics, Columbia University, New York, New York 10027, USA 2Laboratory for Laser Energetics, University of Rochester, Rochester, New York 14623, USA 3Princeton Plasma Physics Laboratory, Princeton, New Jersey 08543, USA

High poloidal beta equilibria in the Tokamak Fusion Test Reactor limited by a natural inboard poloidal field null

Recent operation of the Tokamak Fusion Test Reactor (TFTR) [Plasma Phys. Controlled Nucl. Fusion Research 1, 51 (1986)] has produced plasma equilibria with values of Λ≡βp eq+li/2 as large as 7, eβp dia≡2μ0e〈p⊥〉/〈〈Bp〉〉2 as large as 1.6, and Troyon normalized diamagnetic beta [Plasma Phys. Controlled Fusion 26, 209 (1984); Phys. Lett. 110A, 29 (1985)], βNdia≡108〈βt⊥〉aB0/Ip as large as 4.7. When eβp dia≳1.25, a separatrix entered the vacuum chamber, producing a naturally diverted discharge that was sustained for many energy confinement times, τE. The largest values of eβp and plasma stored energy were obtained when the plasma current was ramped down prior to neutral beam injection. The measured peak ion and electron temperatures were as large as 24 and 8.5 keV, respectively. Plasma stored energy in excess of 2.5 MJ and τE greater than 130 msec were obtained. Confinement times of greater than 3 times that expected from L‐mode predictions have been achieved. The fusion power gain QDD reached a value of 1.3×10−...

Overview of TFTR transport studies

A review of TFTR plasma transport studies is presented. Parallel transport and the confinement of suprathermal ions are found to be relatively well described by theory. Cross-field transport of the thermal plasma, however, is anomalous with the momentum diffusivity being comparable to the ion thermal diffusivity and larger than the electron thermal diffusivity in neutral beam heated discharges. Perturbative experiments have studied nonlinear dependencies in the transport coefficients and examined the role of possible nonlocal phenomena. The underlying turbulence has been studied using microwave scattering, beam emission spectroscopy and microwave reflectometry over a much broader range in k perpendicular to than previously possible. Results indicate the existence of large-wavelength fluctuations correlated with enhanced transport.