Tongnyeol Rhee

Gyro-fluid and two-fluid theory and simulations of edge-localized-modesa)

This paper reports on the theoretical and simulation results of a gyro-Landau-fluid extension of the BOUT++ code, which contributes to increasing the physics understanding of edge-localized-modes (ELMs). Large ELMs with low-to-intermediate-n peeling-ballooning (P-B) modes are significantly suppressed due to finite Larmor radius (FLR) effects when the ion temperature increases. For type-I ELMs, it is found from linear simulations that retaining complete first order FLR corrections as resulting from the incomplete “gyroviscous cancellation” in Braginskiis two-fluid model is necessary to obtain good agreement with gyro-fluid results for high ion temperature cases (Ti≽3  keV) when the ion density has a strong radial variation, which goes beyond the simple local model of ion diamagnetic stabilization of ideal ballooning modes. The maximum growth rate is inversely proportional to Ti because the FLR effect is proportional to Ti. The FLR effect is also proportional to toroidal mode number n, so for high n cases,...

Analysis of symmetry breaking mechanisms and the role of turbulence self-regulation in intrinsic rotation

We present analyses of mechanisms which convert radial inhomogeneity to broken k||-symmetry and thus produce turbulence driven intrinsic rotation in tokamak plasmas. By performing gyrokinetic simulations of ITG turbulence, we explore the many origins of broken k||-symmetry in the fluctuation spectrum and identify both E ? B shear and the radial gradient of turbulence intensity?a ubiquitous radial inhomogeneity in tokamak plasmas?as important k||-symmetry breaking mechanisms. By studying and comparing the correlations between residual stress, E ? B shearing, fluctuation intensity and its radial gradient, we investigate the dynamics of residual stress generation by various symmetry breaking mechanisms and explore the implication of the self-regulating dynamics of fluctuation intensity and E ? B shearing for intrinsic rotation generation. Several scalings for intrinsic rotation are reported and are linked to investigations of underlying local dynamics. It is found that stronger intrinsic rotation is generated for higher values of ion temperature gradient, safety factor and weaker magnetic shear. These trends are broadly consistent with the intrinsic rotation scaling found from experiment?the so-called Rice scaling.

Turbulent acceleration of superthermal electrons

In a recent paper it was suggested on the basis of weak turbulence theory that the collisionality of a plasma, coupled with nonlinear wave-particle interaction, is crucial for the acceleration of electrons by Langmuir turbulence to a superthermal energy level. In this Letter, fully nonlinear Vlasov and particle-in-cell (PIC) simulation techniques are employed to further verify this potentially important finding. The previous conclusion is fully confirmed by observing the expected difference between the Vlasov and PIC simulation results in the weak beam regime. However, in the strong beam regime, both the Vlasov and PIC simulations are found to produce a high-energy tail population, which indicates that there may be other mechanisms in the high beam speed situation, that are responsible for the generation the superthermal electrons.

On the mechanism for edge localized mode mitigation by supersonic molecular beam injection

We construct a diffusive, bi-stable cellular automata model to elucidate the physical mechanisms underlying observed edge localized mode (ELM) mitigation by supersonic molecular beam injection (SMBI). The extended cellular automata model reproduces key qualitative features of ELM mitigation experiments, most significantly the increase in frequency of grain ejection events (ELMs), and the decrease in the number of grains ejected by these transport events. The basic mechanism of mitigation is the triggering of small scale pedestal avalanches by additional grain injection directly into the H-mode pedestal. The small scale avalanches prevent the gradient from building-up to marginality throughout the pedestal, thus avoiding large scale transport events which span the full extent of that region. We explore different grain injection parameters to find an optimal SMBI scenario. We show that shallow SMBI deposition is sufficient for ELM mitigation.

Multiple Harmonic Plasma Emission

Electromagnetic radiation at the plasma frequency and/or its second harmonic, the so-called plasma emission, is widely accepted as the fundamental process responsible for solar type II and III radio bursts. There have also been occasional observations of higher-harmonic plasma emissions in the solar-terrestrial environment. This paper presents the first demonstration of multiple harmonic emission by means of two-dimensional electromagnetic particle-in-cell simulation. This finding indicates that under certain circumstances the traditional mechanism of fundamental-harmonic pair emission might also be accompanied by higher-harmonic components. Consequently, the present findings are highly relevant to in situ observations of third- and/or higher-harmonic plasma emission in astrophysical and solar-terrestrial environments.

Flux-driven simulations of turbulence collapse

Using three-dimensional nonlinear simulations of tokamak turbulence, we show that an edge transport barrier (ETB) forms naturally once input power exceeds a threshold value. Profiles, turbulence-driven flows, and neoclassical coefficients are evolved self-consistently. A slow power ramp-up simulation shows that ETB transition is triggered by the turbulence-driven flows via an intermediate phase which involves coherent oscillation of turbulence intensity and E×B flow shear. A novel observation of the evolution is that the turbulence collapses and the ETB transition begins when RT > 1 at t = tR (RT: normalized Reynolds power), while the conventional transition criterion ( ωE×B>γlin where ωE×B denotes mean flow shear) is satisfied only after t = tC ( >tR), when the mean flow shear grows due to positive feedback.

Simulation and theory for two-dimensional beam-plasma instability

A comparative study of the dynamics of the electron beam-plasma system in two spatial dimensions is carried out by means of particle-in-cell (PIC) simulation and quasilinear theory. In the literature, the beam-plasma instability is usually studied with one-dimensional assumption. Among the few works that include higher-dimensional effects are two- and three-dimensional quasilinear theory and two-dimensional PIC simulations. However, no efforts were made to compare the theory and simulation side by side. The present paper carries out a detailed comparative study of two-dimensional simulation and quasilinear theory. It is found that the quasilinear theory quite adequately accounts for most important features associated with the simulation result. For instance, the particle diffusion time scale, the maximum wave intensity, dynamical development of the electron distribution function, and the change in the wave spectrum all agree quantitatively. However, certain nonlinear effects such as the Langmuir condensation phenomenon are not reproduced by the quasilinear theory. Nevertheless, the present paper verifies that the simple quasilinear theory is quite effective for the study of beam-plasma instability for the present choice of parameters.

Simulation of asymmetric solar wind electron distributions

The electron distributions detected in the solar wind feature varying degrees of anisotropic high-energy tail. In a recent work the present authors numerically solved the one-dimensional electrostatic weak turbulence equations by assuming that the solar wind electrons are initially composed of thermal core plus field-aligned counterstreaming beams, and demonstrated that a wide variety of asymmetric energetic tail distribution may result. In the present paper, the essential findings in this work are tested by means of full particle-in-cell simulation technique. It is found that the previous results are largely confirmed, thus providing evidence that the paradigm of local electron acceleration to high-energy tail by self-consistently excited Langmuir turbulence may be relevant to the solar wind environment under certain circumstances. However, some discrepancies are found such that the nearly one-sided energetic tail reported in the numerical solution of the weak turbulence kinetic equation is not shown.

Effects of spontaneous thermal fluctuations on nonlinear beam-plasma interaction

In this paper, the effects of nonvanishing plasma parameter 1∕nλDe3 on the nonlinear beam-plasma interaction process are discussed on the basis of numerical solutions of weak turbulence equation. The finiteness of the plasma parameter is directly related to the single-particle spontaneous fluctuation phenomena. It is shown that spontaneous fluctuations promote the Langmuir condensation effect, lead to a finite level of ambient turbulence, and enhance nonlinear mode coupling process.

Impact of zonal flows on edge pedestal collapse

We perform a computational study of the role of zonal flows in edge pedestal collapse on the basis of a nonlinear three-field reduced magnetohydrodynamic (MHD) model. A dramatic change of dynamics takes place when ideal ballooning modes are completely stabilized. Analyses show that a new instability is developed due to a strong excitation of zonal vorticity, resulting in a series of secondary crashes. The presence of subsidiary bursts after a main crash increases the effective crash time and energy loss. These simulation results resemble the behavior of compound edge localized modes (ELMs). Analyses in this paper indicate that a complete understanding of ELM crash dynamics requires the self-consistent inclusion of nonlinear zonal flows-MHD interaction and transport physics.