Tokamak profile prediction using direct gyrokinetic and neoclassical simulation

Candy, J.; Holland, C.; Waltz, R. E.; Fahey, M. R.; Belli, E.
June 2009
Physics of Plasmas;Jun2009, Vol. 16 Issue 6, p060704
Academic Journal
Tokamak transport modeling scenarios, including ITER [ITER Physics Basis Editors, Nucl. Fusion 39, 2137 (1999)] performance predictions, are based exclusively on reduced models for core thermal and particle transport. The reason for this is simple: computational cost. A typical modeling scenario may require the evaluation of thousands of individual transport fluxes (local transport models calculate the energy and particle fluxes across a specified flux surface given fixed profiles). Despite continuous advances in direct gyrokinetic simulation, the cost of an individual simulation remains so high that direct gyrokinetic transport calculations have been avoided. By developing a steady-state iteration scheme suitable for direct gyrokinetic and neoclassical simulations, we can now compute steady-state temperature profiles for DIII-D [J. L. Luxon, Nucl. Fusion 42, 614 (2002)] plasmas given known plasma sources. The new code, TGYRO, encapsulates the GYRO [J. Candy and R. E. Waltz, J. Comput. Phys. 186, 545 (2003)] code, for turbulent transport, and the NEO [E. A. Belli and J. Candy, Plasma Phys. Controlled Fusion 50, 095010 (2008)] code, for kinetic neoclassical transport. Results for DIII-D L-mode discharge 128913 are given, with computational and experimental results consistent in the region 0≤r/a≤0.8.


Related Articles

  • Rotation velocities in the plasma edge driven viscously by scrape-off layer flows. Stacey, W. M. // Physics of Plasmas;Jun2009, Vol. 16 Issue 6, p062505 

    Scrape-off layer parallel flows and the viscous fluxes in the plasma edge driven thereby are calculated from neoclassical theory for a model problem representative of a present experiment, using an analytical model for elongated flux surface geometry with a Shafranov shift to provide a realistic...

  • Shock Hugoniot and temperature data for polystyrene obtained with quartz standard. Ozaki, N.; Sano, T.; Ikoma, M.; Shigemori, K.; Kimura, T.; Miyanishi, K.; Vinci, T.; Ree, F. H.; Azechi, H.; Endo, T.; Hironaka, Y.; Hori, Y.; Iwamoto, A.; Kadono, T.; Nagatomo, H.; Nakai, M.; Norimatsu, T.; Okuchi, T.; Otani, K.; Sakaiya, T. // Physics of Plasmas;Jun2009, Vol. 16 Issue 6, p062702 

    Equation-of-state data, not only pressure and density but also temperature, for polystyrene (CH) are obtained up to 510 GPa. The region investigated in this work corresponds to an intermediate region, bridging a large gap between available gas-gun data below 60 GPa and laser shock data above 500...

  • Fokker–Planck simulations for core heating in subignition cone-guiding fast ignition targets. Johzaki, Tomoyuki; Nakao, Yasuyuki; Mima, Kunioki // Physics of Plasmas;Jun2009, Vol. 16 Issue 6, p062706 

    On the basis of two-dimensional coupled relativistic Fokker–Planck and fluid simulations, the core heating by fast electron beam in subignition cone-guiding fast ignition targets is investigated. It was found that the magnetic field due to the thermoelectric current reduces the beam...

  • Wave coupling in sheet- and multiple-beam traveling-wave tubes. Nusinovich, Gregory S.; Cooke, Simon J.; Botton, Moti; Levush, Baruch // Physics of Plasmas;Jun2009, Vol. 16 Issue 6, p063102 

    To increase the power level of the sources of coherent electromagnetic radiation at frequencies from 100 GHz up to the terahertz range it makes sense to develop devices with a spatially extended interaction space. Sheet-beam and multiple-beam devices belong to the category. In the present paper...

  • On the reliability of the Bhatnagar–Gross–Krook collision model in weakly ionized plasmas. Else, D.; Kompaneets, R.; Vladimirov, S. V. // Physics of Plasmas;Jun2009, Vol. 16 Issue 6, p062106 

    The ion velocity distribution in weakly ionized plasma with mobility-limited ion flow driven by an electric field is studied numerically, with the assumption that the charge exchange is the dominant mechanism of ion scattering. The calculations are performed using the realistic constant mean...

  • Charging and E×B rotation of ablation clouds surrounding refueling pellets in hot fusion plasmas. Parks, P. B.; Lu, T.; Samulyak, R. // Physics of Plasmas;Jun2009, Vol. 16 Issue 6, p060705 

    The finite resistivity magnetohydrodynamic code FRONTIER-MHD [R. Samulyak et al., Nucl. Fusion 47, 103 (2007)] is used to simulate the ablation rate of refueling pellets, including the novel effect of electrostatically induced E×B rotation of the ablation cloud about its symmetry axis...

  • Anomalous resistivity due to kink modes in a thin current sheet. Moritaka, Toseo; Horiuchi, Ritoku; Ohtani, Hiroaki // Physics of Plasmas;Oct2007, Vol. 14 Issue 10, p102109 

    The roles of microscopic plasma instabilities on the violation of the frozen-in constraint are investigated by examining the force balance equation based on explicit electromagnetic particle simulation for a thin current sheet. Wave-particle interactions associated with lower hybrid drift...

  • The transition from plasma gratings to cavitons in laser-plasma interactions. Liu, Z. J.; He, X. T.; Zheng, C. Y.; Wang, Y. G. // Physics of Plasmas;Sep2009, Vol. 16 Issue 9, p093108 

    One-dimensional Vlasov–Maxwell simulations of laser-plasma interactions are presented. It is shown that plasma gratings and density cavitons are formed sequentially. There are strong electromagnetic fields in the cavitons and the electromagnetic structures are nearly standing and...

  • Landau and non-Landau linear damping: Physics of the dissipation. Chust, T.; Belmont, G.; Mottez, F.; Hess, S. // Physics of Plasmas;Sep2009, Vol. 16 Issue 9, p092104 

    For linear Langmuir waves, it is well known that the energy exchanges generally lead to a continuous dissipation, on average, from the electric form to the kinetic one. Many papers have estimated these exchanges and indeed shown that the classical Landau value γL, characterizing the electric...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics