5 research outputs found
Linear and nonlinear verification of gyrokinetic microstability codes
Verification of nonlinear microstability codes is a necessary step before comparisons or predictions of turbulent transport in toroidal devices can be justified. By verification we mean demonstrating that a code correctly solves the mathematical model upon which it is based. Some degree of verification can be accomplished indirectly from analytical instability threshold conditions, nonlinear saturation estimates, etc., for relatively simple plasmas. However, verification for experimentally relevant plasma conditions and physics is beyond the realm of analytical treatment and must rely on code-to-code comparisons, i.e., benchmarking. The premise is that the codes are verified for a given problem or set of parameters if they all agree within a specified tolerance. True verification requires comparisons for a number of plasma conditions, e.g., different devices, discharges, times, and radii. Running the codes and keeping track of linear and nonlinear inputs and results for all conditions could be prohibitive unless there was some degree of automation. We have written software to do just this and have formulated a metric for assessing agreement of nonlinear simulations. We present comparisons, both linear and nonlinear, between the gyrokinetic codes GYRO J. Candy and R. E. Waltz, J. Comput. Phys. 186, 545 (2003) and GS2 W. Dorland, F. Jenko, M. Kotschenreuther, and B. N. Rogers, Phys. Rev. Lett. 85, 5579 (2000). We do so at the mid-radius for the same discharge as in earlier work C. Holland, A. E. White, G. R. McKee, M. W. Shafer, J. Candy, R. E. Waltz, L. Schmitz, and G. R. Tynan, Phys. Plasmas 16, 052301 (2009). The comparisons include electromagnetic fluctuations, passing and trapped electrons, plasma shaping, one kinetic impurity, and finite Debye-length effects. Results neglecting and including electron collisions (Lorentz model) are presented. We find that the linear frequencies with or without collisions agree well between codes, as do the time averages of the nonlinear fluxes without collisions. With collisions, the differences between the time-averaged fluxes are larger than the uncertainties defined as the oscillations of the fluxes, with the GS2 fluxes consistently larger (or more positive) than those from GYRO. However, the electrostatic fluxes are much smaller than those without collisions (the electromagnetic energy flux is negligible in both cases). In fact, except for the electron energy fluxes, the absolute magnitudes of the differences in fluxes with collisions are the same or smaller than those without. None of the fluxes exhibit large absolute differences between codes. Beyond these results, the specific linear and nonlinear benchmarks proposed here, as well as the underlying methodology, provide the basis for a wide variety of future verification efforts. © 2011 American Institute of Physics
Recommended from our members
Compatibility of internal transport barrier with steady-state operation in the high bootstrap fraction regime on DIII-D
Recent EAST/DIII-D joint experiments on the high poloidal beta tokamak regime in DIII-D have demonstrated fully noninductive operation with an internal transport barrier (ITB) at large minor radius, at normalized fusion performance increased by 30% relative to earlier work (Politzer et al 2005 Nucl. Fusion 45 417). The advancement was enabled by improved understanding of the 'relaxation oscillations', previously attributed to repetitive ITB collapses, and of the fast ion behavior in this regime. It was found that the 'relaxation oscillations' are coupled core-edge modes amenable to wall-stabilization, and that fast ion losses which previously dictated a large plasma-wall separation to avoid wall over-heating, can be reduced to classical levels with sufficient plasma density. By using optimized waveforms of the plasma-wall separation and plasma density, fully noninductive plasmas have been sustained for long durations with excellent energy confinement quality, bootstrap fraction ≥80%, βN ≤ 4, βP ≥ 3, and βT ≥ 2%. These results bolster the applicability of the high poloidal beta tokamak regime toward the realization of a steady-state fusion reactor
Recommended from our members
Compatibility of internal transport barrier with steady-state operation in the high bootstrap fraction regime on DIII-D
Recent EAST/DIII-D joint experiments on the high poloidal beta tokamak regime in DIII-D have demonstrated fully noninductive operation with an internal transport barrier (ITB) at large minor radius, at normalized fusion performance increased by 30% relative to earlier work (Politzer et al 2005 Nucl. Fusion 45 417). The advancement was enabled by improved understanding of the 'relaxation oscillations', previously attributed to repetitive ITB collapses, and of the fast ion behavior in this regime. It was found that the 'relaxation oscillations' are coupled core-edge modes amenable to wall-stabilization, and that fast ion losses which previously dictated a large plasma-wall separation to avoid wall over-heating, can be reduced to classical levels with sufficient plasma density. By using optimized waveforms of the plasma-wall separation and plasma density, fully noninductive plasmas have been sustained for long durations with excellent energy confinement quality, bootstrap fraction ≥80%, βN ≤ 4, βP ≥ 3, and βT ≥ 2%. These results bolster the applicability of the high poloidal beta tokamak regime toward the realization of a steady-state fusion reactor