Source author record

S. Ku

S. Ku appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2021arXiv

Constructing a new predictive scaling formula for ITER's divertor heat-load width informed by a simulation-anchored machine learning

Understanding and predicting divertor heat-load width $λ_q$ is a critically important problem for an easier and more robust operation of ITER with high fusion gain. Previous predictive simulation data for $λ_q$ using the extreme-scale edge gyrokinetic code XGC1 in the electrostatic limit under attached divertor plasma conditions in three major US tokamaks [C.S. Chang et al., Nucl. Fusion 57, 116023 (2017)] reproduced the Eich and Goldston attached-divertor formula results [formula #14 in T. Eich et al., Nucl. Fusion 53, 093031 (2013); R.J. Goldston, Nucl. Fusion 52, 013009 (2012)], and furthermore predicted over six times wider $λ_q$ than the maximal Eich and Goldston formula predictions on a full-power (Q = 10) scenario ITER plasma. After adding data from further predictive simulations on a highest current JET and highest-current Alcator C-Mod, a machine learning program is used to identify a new scaling formula for $λ_q$ as a simple modification to the Eich formula #14, which reproduces the Eich scaling formula for the present tokamaks and which embraces the wide $λ_q^X{GC}$ for the full-current Q = 10 ITER plasma. The new formula is then successfully tested on three more ITER plasmas: two corresponding to long burning scenarios with Q = 5 and one at low plasma current to be explored in the initial phases of ITER operation. The new physics that gives rise to the wider $λq_^{XGC} is identified to be the weakly-collisional, trapped-electron-mode turbulence across the magnetic separatrix, which is known to be an efficient transporter of the electron heat and mass. Electromagnetic turbulence and high-collisionality effects on the new formula are the next study topics for XGC1.

preprint2021arXiv

Effects of collisional ion orbit loss on neoclassical tokamak radial electric fields

Ion orbit loss is considered important for generating the radially inward electric field $E_r$ in a tokamak edge plasma. In particular, this effect is emphasized in diverted tokamaks with a magnetic X point. In neoclassical equilibria, Coulomb collisions can scatter ions onto loss orbits and generate a radially outward current, which in steady state is balanced by the radially inward current from viscosity. To quantitatively measure this loss-orbit current in an edge pedestal, an ion-orbit-flux diagnostic has been implemented in the axisymmetric version of the gyrokinetic particle-in-cell code XGC. As the first application of this diagnostic, a neoclassical DIII-D H-mode plasma is studied using gyrokinetic ions and adiabatic electrons. The validity of the diagnostic is demonstrated by studying the collisional relaxation of $E_r$ in the core. After this demonstration, the loss-orbit current is numerically measured in the edge pedestal in quasisteady state. In this plasma, it is found that the radial electric force on ions from $E_r$ approximately balances the ion radial pressure gradient in the edge pedestal, with the radial force from the plasma flow term being a minor component. The effect of orbit loss on $E_r$ is found to be only mild. How ion orbit loss will affect $E_r$ in the full-current ITER plasma pedestal is left as a subsequent study topic.

preprint2020arXiv

Comparison of edge turbulence characteristics between DIII-D and C-Mod simulations with XGC1

The physical processes taking place at the edge region are crucial for the operation of tokamaks as they govern the interaction of hot plasma with the vessel walls. Numerical modeling of the edge with state-of-the-art codes attempts to elucidate interactions between neoclassical drifts, turbulence, poloidal and parallel flows that control the physical set-up of the SOL region. Here, we present post-processing analysis of simulations from the gyrokinetic code XGC1, comparing edge turbulence characteristics from a simulation of DIII-D against one of C-Mod. We find that the equilibrium $E \times B$ flux across the separatrix has a similar poloidal pattern in both discharges which can be explained by magnetic drifts and trapped ion excursions. However, collisionality is noted to play a major role in that it prevents local charge accumulations from having global effects in C-Mod. In both cases, turbulent electron heat flux is higher than the ion one. This seems to be a universal characteristic of the tokamak edge. We identify turbulent frequencies and growth rates of the dominant mode in both simulations. In C-Mod, these numbers point to the presence of a drift wave. In DIII-D, linear simulations with Gene reveal a trapped electron mode. Furthermore, we present the amplitude and size distributions of the blobs from both simulations. Amplitude distributions are in qualitative agreement with experimental observations while size distributions are consistent with the fact that most blobs are not connecting to the divertor plates and suggest that they are generated by the shearing of the turbulent modes.

preprint2015arXiv

Kinetic modeling of divertor heat load fluxes in the Alcator C-Mod and DIII-D tokamaks

The guiding-center kinetic neoclassical transport code, XGC0, [C.S. Chang et. al, Phys. Plasmas 11, 2649 (2004)] is used to compute the heat fluxes and the heat-load width in the outer divertor plates of Alcator C-Mod and DIII-D tokamaks. The dependence of the width of heat-load fluxes on neoclassical effects, neutral collisions and anomalous transport is investigated using the XGC0 code. The XGC0 code includes realistic X-point geometry, a neutral source model, the effects of collisions, and a diffusion model for anomalous transport. It is observed that width of the XGC0 neoclassical heat-load is approximately inversely proportional to the total plasma current $I_{\rm p}$. The scaling of the width of the divertor heat-load with plasma current is examined for an Alcator C-Mod discharge and four DIII-D discharges. The scaling of the divertor heat-load width with plasma current is found to be weaker in the Alcator C-Mod discharge compared to scaling found in the DIII-D discharges. The effect of neutral collisions on the $1/I_{\rm p}$ scaling of heat-load width is shown not to be significant. Although inclusion of poloidally uniform anomalous transport results in a deviation from the $1/I_{\rm p}$ scaling, the inclusion of the anomalous transport that is driven by ballooning-type instabilities results in recovering the neoclassical $1/I_{\rm p}$ scaling. The Bohm or Gyro-Bohm scalings of anomalous transport does not strongly affect the dependence of the heat-load width on plasma current. The inclusion of anomalous transport, in general, results in widening the width of neoclassical divertor heat-load and enhances the neoclassical heat-load fluxes on the divertor plates. Understanding heat transport in the tokamak scrape-off layer plasmas is important for strengthening the basis for predicting divertor conditions in ITER.

preprint2010arXiv

Backward Raman compression of x-rays in metals and warm dense matters

Experimentally observed decay rate of the long wavelength Langmuir wave in metals and dense plasmas is orders of magnitude larger than the prediction of the prevalent Landau damping theory. The discrepancy is explored, and the existence of a regime where the forward Raman scattering is stable and the backward Raman scattering is unstable is examined. The amplification of an x-ray pulse in this regime, via the backward Raman compression, is computationally demonstrated, and the optimal pulse duration and intensity is estimated.

preprint2010arXiv

Photonic band gap and x-ray optics in warm dense matter

Photonic band gaps for the soft x-rays, formed in the periodic structures of solids or dense plasmas, are theoretically investigated. Optical manipulation mechanisms for the soft x-rays, which are based on these band gaps, are computationally demonstrated. The reflection and amplification of the soft x-rays, and the compression and stretching of chirped soft x-ray pulses are discussed. A scheme for lasing with atoms with two energy levels, utilizing the band gap, is also studied.

preprint2010arXiv

Theory of plasmon decay in dense plasmas and warm dense matter

The decay of the Langmuir waves in dense plasmas is not accurately predicted by the prevalent Landau damping theory. A dielectric function theory is introduced, predicting much higher damping than the Landau damping theory. This strong damping is in better agreement with the experimentally observed data in metals. It is shown that the strong plasmon decay leads to the existence of a parameter regime where the backward Raman scattering is unstable while the forward Raman scattering is stable. This regime may be used to create intense x-ray pulses, by means of the the backward Raman compression. The optimal pulse duration and intensity is estimated.