We report the status of hybrid scenario experiments in Korea Superconducting Tokamak Advanced Research (KSTAR). The hybrid scenario is defined as stationary discharges with β N ⩾ 2.4 and H 89 ⩾ 2.0 at q 95 < 6.5 without or with very mild sawtooth activities in KSTAR. It is being developed towards reactor-relevant conditions. High performance of β N ≲ 3.0, H 89 ≲ 2.4 and G-factor (≡ β N H 89 /q 2 95 ) ≲ 0.46 has been achieved and sustained for ≳ 40τ E at n e /n GW ~0.7 with heating power of ≲5 MW. Some KSTAR hybrid discharges exhibit a unique feature of a slow transition from conventional H-mode to hybrid mode after the third neutral beam injection. The reason for the confinement enhancement is extensively studied in this transition period of a representative discharge exhibiting a common feature of KSTAR hybrid scenarios. 0D performance analysis with magnetohydrodynamic activities, 1D kinetic profile dynamics, power balance analysis, linear gyro-kinetic analysis and edge pedestal stability analysis were conducted. The enhancement is thought to be from both the core and the pedestal. The improvement in the core region of the ion energy channel is observed from the linear gyro-kinetic analysis considering the electromagnetic, the fast ion, the Shafranov shift, ω E×B , and the magnetic shear effect. The electromagnetic finite β stabilisation plays a role in the inner core region at ρ tor ∼ 0.35 together with the fast ion effect. The alpha stabilisation effect is also found at ρ tor ∼ 0.5. ω E×B , which could reduce the linear growth of the ion temperature gradient mode in the outer core region at ρ tor ∼ 0.5 − 0.7 with the highest contribution from the toroidal rotation. Regarding the improvement in the pedestal, Shafranov shift broadens the stability boundary of the pedestal in support of the diamagnetic effect. The pedestal height and width could be reproduced by the EPED model, while a realistic current profile is used to calculate the internal inductance for Shafranov shift. Based on these findings, a comprehensive confinement enhancement mechanism has been proposed by considering the core-edge interplay.

A newly developed integrated suite of codes coined as tokamak reactor integrated automated suite for simulation and computation (TRIASSIC) is reported. The suite comprises existing plasma simulation codes, including 1.5D/2D plasma transport solvers and neoclassical/anomalous transport, plasma heating/cooling, and cold neutral models. The components in TRIASSIC are fully modularized by adopting a generic data structure as its internal storage. Primary components such as the transport solver and the neutral beam or electron cyclotron wave actuator were verified to its standalone implementation. The interpretive and predictive simulations of 50 stationary plasma phases from 30 KSTAR discharges were conducted for validation, and a good agreement with experimental measurements on all simulation cases was found.

We report results of benchmarking of core particle transport simulations by a collection of codes widely used in transport modelling of tokamak plasmas. Our analysis includes formulation of transport equations, difference between electron and ion solvers, comparison of modules of the 2 pellet and edge gas fuelling on the ITER baseline scenario. During the first phase of benchmarking we address the particle transport effects in the stationary phase. Firstly, simulations are performed with identical sources, sinks, transport coefficients, and boundary conditions prescribed in the flattop H-mode phase. The transformation of ion particle transport equations is introduced so to directly compare their results to electron transport solvers. Secondly, the pellet fuelling models are benchmarked in various conditions to evaluate the dependency of the pellet deposition on the pellet volume, injection side, pedestal and separatrix parameters. Thirdly, edge gas fuelling is benchmarked to assess sensitivities of source profile predictions to uncertainties in plasma conditions and detailed model assumptions. At the second phase, we address particle transport effects in the time-evolving plasma including the current ramp-up to the ramp-down phase. The ion and the electron solvers are benchmarked together. Differences between the simulation results of the solvers are investigated in terms of equilibrium, grid resolution, radial coordinate, radial grid distribution, and plasma volume evolution term. We found that the selection of the radial coordinate can yield prominent differences between the solvers mainly due to differences in the edge grid distribution. The simulations reveal that electron and ion solvers predict noticeably different density peaking for the same diffusion and pinch velocity while with the peaked profile of helium, expected in fusion reactors. The fuelling benchmarking shows that gas puffing is not efficient for core fuelling in H-modes and density control should be done by the high field side pellet injection in contrast to present machines. mode in the current flattop phase depends sensitively on the particle balance of the mixed D-T fuels, He and impurities. In ITER, the neutral beam injection (NBI) does not play a noticeable role neither in the global particle balance [7], nor for the central fuelling. Moreover, the SOLPS modelling [7] predicts dramatic reduction of the gas penetrated from the edge, making the pellet injection the main tool for the density control in the H-mode plasmas, though the gas penetrated from the edge still can play the dominant role for the L-mode operation. Features like the recycling and penetration of He and the fuel into the core plasma are central to understanding the dilution and tritium burnup. The SOL/divertor plasma and its interactions with plasma facing components will set the boundary conditions for the core transport. Eventually, the particle transport alters the heat and the momentum transport so all these non-linear connections need to be understood simultaneously t...

No abstract

A tokamak, a torus-shaped nuclear fusion device, needs an electric current in the plasma to produce magnetic field in the poloidal direction for confining fusion plasmas. Plasma current is conventionally generated by electromagnetic induction. However, for a steady-state fusion reactor, minimizing the inductive current is essential to extend the tokamak operating duration. Several non-inductive current drive schemes have been developed for steady-state operations such as radio-frequency waves and neutral beams. However, commercial reactors require minimal use of these external sources to maximize the fusion gain, Q, the ratio of the fusion power to the external power. Apart from these external current drives, a self-generated current, so-called bootstrap current, was predicted theoretically and demonstrated experimentally. Here, we reveal another self-generated current that can exist in a tokamak and this has not yet been discussed by present theories. We report conclusive experimental evidence of this self-generated current observed in the KSTAR tokamak.

No abstract