Jeff Candy

Bio: Jeff Candy is an academic researcher from General Atomics. The author has contributed to research in topics: Turbulence & Tokamak. The author has an hindex of 54, co-authored 205 publications receiving 8220 citations. Previous affiliations of Jeff Candy include Joint European Torus & University of California, San Diego.

Anomalous Transport Scaling in the DIII-D Tokamak Matched by Supercomputer Simulation

Abstract: Gyrokinetic simulation of tokamak transport has evolved sufficiently to allow direct comparison of numerical results with experimental data. It is to be emphasized that only with the simultaneous inclusion of many distinct and complex effects can this comparison realistically be made. Until now, numerical studies of tokamak microturbulence have been restricted to either (a) flux tubes or (b) electrostatic fluctuations. Using a newly developed global electromagnetic solver, we have been able to recover via direct simulation the Bohm-like scaling observed in DIII-D L-mode discharges. We also match, well within experimental uncertainty, the measured energy diffusivities.

Kinetic calculation of neoclassical transport including self-consistent electron and impurity dynamics

Abstract: Numerical studies of neoclassical transport, beginning with the fundamental drift-kinetic equation (DKE), have been extended to include the self-consistent coupling of electrons and multiple ion species. The code, NEO, provides a first-principles based calculation of the neoclassical transport coefficients directly from solution of the distribution function by solving a hierarchy of equations derived by expanding the DKE in powers of ρ*i, the ratio of the ion gyroradius to system size. This includes the calculation of the first-order electrostatic potential via the Poisson equation, although this potential has exactly no effect on the steady-state transport. Systematic calculations of the second-order particle and energy fluxes and first-order plasma flows and bootstrap current and comparisons with existing theories are given for multi-species plasmas. The ambipolar relation ∑azaΓa = 0, which can only be maintained with complete cross-species collisional coupling, is confirmed, and finite mass-ratio corrections due to the collisional coupling are identified. The effects of plasma shaping are also explored, including a discussion of how analytic formulae obtained for circular plasmas (i.e. Chang–Hinton) should be applied to shaped cases. Finite-orbit-width effects are studied via solution of the higher-order DKEs and the implications of non-local transport on the validity of the δf formulation are discussed.

Integrated modeling applications for tokamak experiments with OMFIT

Abstract: One modeling framework for integrated tasks (OMFIT) is a comprehensive integrated modeling framework which has been developed to enable physics codes to interact in complicated workflows, and support scientists at all stages of the modeling cycle. The OMFIT development follows a unique bottom-up approach, where the framework design and capabilities organically evolve to support progressive integration of the components that are required to accomplish physics goals of increasing complexity. OMFIT provides a workflow for easily generating full kinetic equilibrium reconstructions that are constrained by magnetic and motional Stark effect measurements, and kinetic profile information that includes fast-ion pressure modeled by a transport code. It was found that magnetic measurements can be used to quantify the amount of anomalous fast-ion diffusion that is present in DIII-D discharges, and provide an estimate that is consistent with what would be needed for transport simulations to match the measured neutron rates. OMFIT was used to streamline edge-stability analyses, and evaluate the effect of resonant magnetic perturbation (RMP) on the pedestal stability, which have been found to be consistent with the experimental observations. The development of a five-dimensional numerical fluid model for estimating the effects of the interaction between magnetohydrodynamic (MHD) and microturbulence, and its systematic verification against analytic models was also supported by the framework. OMFIT was used for optimizing an innovative high-harmonic fast wave system proposed for DIII-D. For a parallel refractive index , the conditions for strong electron-Landau damping were found to be independent of launched and poloidal angle. OMFIT has been the platform of choice for developing a neural-network based approach to efficiently perform a non-linear multivariate regression of local transport fluxes as a function of local dimensionless parameters. Transport predictions for thousands of DIII-D discharges showed excellent agreement with the power balance calculations across the whole plasma radius and over a broad range of operating regimes. Concerning predictive transport simulations, the framework made possible the design and automation of a workflow that enables self-consistent predictions of kinetic profiles and the plasma equilibrium. It is found that the feedback between the transport fluxes and plasma equilibrium can significantly affect the kinetic profiles predictions. Such a rich set of results provide tangible evidence of how bottom-up approaches can potentially provide a fast track to integrated modeling solutions that are functional, cost-effective, and in sync with the research effort of the community.

Zonal flows in plasma—a review

Abstract: A comprehensive review of zonal flow phenomena in plasmas is presented. While the emphasis is on zonal flows in laboratory plasmas, planetary zonal flows are discussed as well. The review presents the status of theory, numerical simulation and experiments relevant to zonal flows. The emphasis is on developing an integrated understanding of the dynamics of drift wave–zonal flow turbulence by combining detailed studies of the generation of zonal flows by drift waves, the back-interaction of zonal flows on the drift waves, and the various feedback loops by which the system regulates and organizes itself. The implications of zonal flow phenomena for confinement in, and the phenomena of fusion devices are discussed. Special attention is given to the comparison of experiment with theory and to identifying directions for progress in future research.

Chapter 1: Overview and summary

Abstract: The ITER Physics Basis presents and evaluates the physics rules and methodologies for plasma performance projections, which provide the basis for the design of a tokamak burning plasma device whose goal is to demonstrate the scientific and technological feasibility of fusion energy for peaceful purposes. This Chapter summarizes the physics basis for burning plasma projections, which is developed in detail by the ITER Physics Expert Groups in subsequent chapters. To set context, the design guidelines and requirements established in the report of ITER Special Working Group 1 are presented, as are the specifics of the tokamak design developed in the Final Design Report of the ITER Engineering Design Activities, which exemplifies burning tokamak plasma experiments. The behaviour of a tokamak plasma is determined by the interaction of many diverse physics processes, all of which bear on projections for both a burning plasma experiment and an eventual tokamak reactor. Key processes summarized here are energy and particle confinement and the H-mode power threshold; MHD stability, including pressure and density limits, neoclassical islands, error fields, disruptions, sawteeth, and ELMs; power and particle exhaust, involving divertor power dispersal, helium exhaust, fuelling and density control, H-mode edge transition region, erosion of plasma facing components, tritium retention; energetic particle physics; auxiliary power physics; and the physics of plasma diagnostics. Summaries of projection methodologies, together with estimates of their attendant uncertainties, are presented in each of these areas. Since each physics element has its own scaling properties, an integrated experimental demonstration of the balance between the combined processes which obtains in a reactor plasma is inaccessible to contemporary experimental facilities: it requires a reactor scale device. It is argued, moreover, that a burning plasma experiment can be sufficiently flexible to permit operation in a steady state mode, with non-inductive plasma current drive, as well as in a pulsed mode where current is inductively driven. Overall, the ITER Physics Basis can support a range of candidate designs for a tokamak burning plasma facility. For each design, there will remain a significant uncertainty in the projected performance, but the projection methodologies outlined here do suffice to specify the major parameters of such a facility and form the basis for assuring that its phased operation will return sufficient information to design a prototype commercial fusion power reactor, thus fulfilling the goal of the ITER project.

Foundations of nonlinear gyrokinetic theory

Abstract: Nonlinear gyrokinetic equations play a fundamental role in our understanding of the long-time behavior of strongly magnetized plasmas. The foundations of modern nonlinear gyrokinetic the- ory are based on three important pillars: (1) a gyrokinetic Vlasov equation written in terms of a gyrocenter Hamiltonian with quadratic low-frequency ponderomotive-like terms; (2) a set of gyrokinetic Maxwell (Poisson-Ampere) equations written in terms of the gyrocenter Vlasov dis- tribution that contain low-frequency polarization (Poisson) and magnetization (Ampere) terms derived from the quadratic nonlinearities in the gyrocenter Hamiltonian; and (3) an exact energy conservationlaw for the gyrokineticVlasov-Maxwell equations that includes all the relevant linear and nonlinear coupling terms. The foundations of nonlinear gyrokinetic theory are reviewed with an emphasis on the rigorous applications of Lagrangian and Hamiltonian Lie-transform perturba- tion methods used in the variationalderivationof nonlineargyrokineticVlasov-Maxwell equations. The physical motivations and applications of the nonlinear gyrokinetic equations, which describe the turbulent evolution of low-frequency electromagnetic fluctuations in a nonuniform magnetized plasmas with arbitrary magnetic geometry, are also discussed.

Calculation of Gauss quadrature rules.

Abstract: Most numerical integration techniques consist of approximating the integrand by a polynomial in a region or regions and then integrating the polynomial exactly. Often a complicated integrand can be factored into a non-negative ''weight'' function and another function better approximated by a polynomial, thus $\int_{a}^{b} g(t)dt = \int_{a}^{b} \omega (t)f(t)dt \approx \sum_{i=1}^{N} w_i f(t_i)$. Hopefully, the quadrature rule ${\{w_j, t_j\}}_{j=1}^{N}$ corresponding to the weight function $\omega$(t) is available in tabulated form, but more likely it is not. We present here two algorithms for generating the Gaussian quadrature rule defined by the weight function when: a) the three term recurrence relation is known for the orthogonal polynomials generated by $\omega$(t), and b) the moments of the weight function are known or can be calculated.

Chapter 2: Plasma confinement and transport

Abstract: The understanding and predictive capability of transport physics and plasma confinement is reviewed from the perspective of achieving reactor-scale burning plasmas in the ITER tokamak, for both core and edge plasma regions. Very considerable progress has been made in understanding, controlling and predicting tokamak transport across a wide variety of plasma conditions and regimes since the publication of the ITER Physics Basis (IPB) document (1999 Nucl. Fusion 39 2137-2664). Major areas of progress considered here follow. (1) Substantial improvement in the physics content, capability and reliability of transport simulation and modelling codes, leading to much increased theory/experiment interaction as these codes are increasingly used to interpret and predict experiment. (2) Remarkable progress has been made in developing and understanding regimes of improved core confinement. Internal transport barriers and other forms of reduced core transport are now routinely obtained in all the leading tokamak devices worldwide. (3) The importance of controlling the H-mode edge pedestal is now generally recognized. Substantial progress has been made in extending high confinement H-mode operation to the Greenwald density, the demonstration of Type I ELM mitigation and control techniques and systematic explanation of Type I ELM stability. Theory-based predictive capability has also shown progress by integrating the plasma and neutral transport with MHD stability. (4) Transport projections to ITER are now made using three complementary approaches: empirical or global scaling, theory-based transport modelling and dimensionless parameter scaling (previously, empirical scaling was the dominant approach). For the ITER base case or the reference scenario of conventional ELMy H-mode operation, all three techniques predict that ITER will have sufficient confinement to meet its design target of Q = 10 operation, within similar uncertainties.