Engineering Papers⌕ Search

Engineering topics

Wang, Y. F.

Publications and source records attributed to Wang, Y. F..

At least 19 records

Turbulence spreading effects on the ELM size and SOL width

BOUT++ turbulence simulations were performed to investigate the impact of turbulence spreading on the edge localized mode (ELM) size and divertor heat flux width (λ q ) broadening in small ELM regimes. Here, this study is motivated by EAST experiments. BOUT++ linear simulations of a pedestal radial electric field (E r ) scan show that the dominant toroidal number mode (n) shifts from high-n to low-n, with a narrow mode spectrum, and the maximum linear growth rate increases as the pedestal E r well deepens. The nonlinear simulations show that as the net E × B pedestal flow increases, the pressure fluctuation level and its inward penetration beyond the top of the pedestal both increase. This leads to a transition from small ELMs to large ELMs. Both inward and outward turbulence spreading are sensitive to the scrape-off-layer (SOL) plasma profiles. The inward turbulence spreading increases for the steep SOL profiles, leading to increasing pedestal energy loss in the small ELM regime. The SOL width (λ q ) is significantly broadened progressing from the ELM-free to small ELM regime, due to the onset of strong radial turbulent transport. The extent of the SOL width (λ q ) broadening depends strongly on outward turbulence spreading. The fluctuation energy intensity flux Γ ε at the separatrix can be enhanced by increasing either pedestal E r flow shear or local SOL pressure gradient. The λ q is broadened as the fluctuation energy intensity flux Γ ε at the last close flux surface (LCFS) increases. Local SOL E × B flow shear will restrain outward turbulence spreading and the associated heat flux width broadening. Operating in H-mode with small ELMs has the potential to solve two critical problems: reducing the ELM size and broadening the SOL width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Long-pulse high-performance H-mode plasmas achieved on EAST

A record duration of a 310 s H-mode plasma (H 98y2 ~ 1.3, n e /n GW ~ 0.7, f BS > 50%) has been recently achieved on experimental advanced superconducting tokamak (EAST) with metal walls, exploiting the device's improved long-pulse capabilities. The experiment demonstrates good control of tungsten concentration, core/edge MHD stability, and particle and heat exhaust with an ITER-like tungsten divertor and zero injected torque, establishing a milestone on the path to steady-state long-pulse high-performance scenarios in support of ITER and CFETR. Important synergistic effects are leveraged toward this result, which relies purely on radio frequency (RF) powers for heating and current drive (H&CD). On-axis electron cyclotron heating enhances the H&CD efficiency from lower hybrid wave injection, increasing confinement quality and enabling fully non-inductive operation at high density (n e /n GW ~ 70%) and high poloidal beta (β P ~ 2.5). A small-amplitude grassy edge localized mode regime facilitates the RF power coupling to the H-mode edge and reduces divertor sputtering/erosion. The high energy confinement quality (H 98y2 ~ 1.3) is achieved with the experimental and simulated results pointing to the strong effect of Shafranov shift on turbulence. Transport analysis suggests that trapped electron modes dominate in the core region during the record discharge. The detailed physics processes (RF synergy, core-edge integration, confinement properties, etc.) of the steady-state operation will be illustrated in the content. In the future, EAST will aim at accessing more relevant dimensionless parameters to develop long-pulse high-performance plasma toward ITER and CFETR steady-state advanced operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of aspect ratio on plasma response to resonant magnetic perturbations in tokamak devices

A systematic numerical study is carried out, computing and comparing the plasma response to the resonant magnetic perturbation (RMP) field, applied for controlling edge localized modes (ELMs), in a series of tokamak plasmas with varying aspect ratio and utilizing the MARS-F code. The aspect ratio is scanned either by varying the plasma major radius at a fixed minor radius or by varying the latter while fixing the former. Both approaches yield similar results when compared in terms of quantities with proper normalizations. In general, a non-monotonic dependence of the resonant response field (normalized by the vacuum counterpart) near the plasma edge is found with varying aspect ratio, indicating that a given ELM control coil current configuration strongly favors plasmas with a certain aspect ratio. This optimal aspect ratio, on the other hand, depends on the toroidal as well as poloidal (i.e., coil phasing) spectra of the applied RMP field. The equilibrium (edge) safety factor, the plasma shape, and the plasma toroidal flow are all fixed to ensure that the effects identified here are predominantly due to the plasma aspect ratio.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Impact of pedestal density gradient and collisionality on ELM dynamics

BOUT++ turbulence simulations are conducted to capture the underlying physics of small ELM characteristics achieved by increasing separatrix density via controlling strike points from vertical to horizontal divertor plates for three EAST discharges. BOUT++ linear simulations show that the most unstable modes change from high-n ideal ballooning modes to intermediate-n peeling–ballooning modes and eventually to peeling–ballooning stable plasmas in the pedestal. Nonlinear simulations show that the fluctuation is saturated at a high level for the lowest separatrix density. The ELM size decreases with increasing separatrix density, until the fraction of this energy lost during the ELM crash becomes less than 1% of the pedestal stored energy, leading to small ELMs. Simulations indicate that small ELMs can be triggered either by the marginally peeling–ballooning instability near the peak pressure gradient position inside the pedestal or by a local instability in the pedestal foot with a larger separatrix density gradient. The pedestal collisionality scan for type-I ELMs with steep pedestal density gradient shows that both linear growth rate and ELM size decrease with increasing collisionality. On the contrary, the pedestal collisionality and pedestal density width scan with a weak pedestal density gradient indicate small ELMs can either be triggered by a high-n ballooning mode or by a low-n peeling mode in a low collisionality region 0.04–0.1. We report the simulations indicate the weaker the linear unstable modes near marginal stability with small linear growth rate, the lower nonlinearly saturated fluctuation intensity and the smaller turbulence spreading from the linear unstable zone to stable zone in the nonlinear saturation phase, leading to small ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characteristics of grassy ELMs and their impact on the divertor heat flux width

We report BOUT++ turbulence simulations are conducted for a 60 s steady-state long pulse high β p EAST grassy ELM discharge. BOUT++ linear simulations show that the unstable mode spectrum covers a range of toroidal mode numbers from low-n (n = 10–15) peeling–ballooning modes (P–B) to high-n (n = 40–80) drift-Alfvén instabilities. Nonlinear simulations show that the ELM crash is triggered by low-n peeling modes and fluctuation is generated at the peak pressure gradient position and radially spread outward into the scrape-off-layer, even though the drift-Alfvén instabilities dominate the linear growth phase. However, drift-Alfvén turbulence delays the onset of the grassy ELM and enhances the energy loss with the fluctuation extending to pedestal top region. Simulations further show that if the peeling drive is removed, the fluctuation amplitude drops by an order of magnitude and the ELM crashes disappear. The divertor heat flux width is ~2 times larger than the estimates based on the HD model and the Eich's ITPA multi-tokamak scaling (or empirical Eich scaling) due to the strong radial turbulence transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Diagnostic development for parallel wave-number measurement of lower hybrid waves in EAST

An eight-channel magnetic probe diagnostic system has been designed and installed adjacent to the 4.6 GHz lower hybrid (LH) grill antenna in the low-field side of the EAST tokamak in order to study the n|| evolution of lower hybrid waves in the first pass from the launcher to the core plasma. The magnetic probes are separated by 6.6 mm, which allows measurement of the dominant parallel refractive index n|| up to n|| =5 for 4.6GHz LH waves. The magnetic probes are designed to be sensitive to the magnetic field component perpendicular to the background magnetic field with a slit on the casing that encloses the probe. The intermediate frequency (IF) stage, which consists of two mixing stages, down-coverts the frequency of the measured wave signals at 4.6 GHz to 20 MHz. A bench test demonstrates the phase stability of the magnetic probe diagnostic system. By evaluating the phase variation of the measured signals along the background magnetic field, the dominant n|| of the LH wave in the scrape-off layer (SOL) has been deduced during the 2019 experimental campaign. In the low density plasma, the measured dominant n|| of the lower hybrid waves is about 2.1, corresponding to the main peak 2.04 of the launched n|| spectrum. The n|| deduced by the least square linear fit method remains near this value in the low density plasma with a high spatial correlation magnitude of 0.9. With an eight-channel probe system a wave-number spectrum has also been deduced, which has a peak near to the measured dominant n||.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of lithium wall conditioning and wave-frequency on high density lower hybrid current drive experiment on EAST

A series of dedicated lower hybrid current drive (LHCD) experiments on EAST shows that lithium wall conditioning extends LH current drive and heating up to the line-averaged density of n ̅_e≈ 4x1019 m-3 for both 2.45 and 4.6 GHz. Current drive at such a high density is crucial for the development of long-pulse non-inductive scenarios on EAST. With lithiation, the LH power injection of 1.5 MW at 2.45 GHz resulted in a drop of loop voltage of ~ 0.3 V, which is a comparable loop voltage drop observed with 1.1 MW at 4.6 GHz. The observed decrease in loop voltage is attributed mostly to the RF heating effect. Another LHCD experiment suggests that lithium wall coating has a more significant impact on the scrape-off-layer (SOL) properties than changes in the Greenwald fraction. LHCD at 2.45 GHz still suffers from a loss of efficiency. Enhanced power ionization in front of the launcher may cause the onset of density-dependent wave instabilities. The rise in the midplane SOL density may also accelerate a transition in the divertor regime, leading to additional ionization and collisional losses in the X-point divertor plasma. Ray-tracing modeling supports that a lower wave frequency is more prone to collisional power loss. The experiments confirm that lithiation is a useful tool to control the SOL plasma, and suggest that density control in front of the launcher may be critical to mitigating power loss mechanisms in the plasma boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physical mechanisms for the transition from type-III to large ELMs induced by impurity injection on EAST

Transition from type-III to large-amplitude ELMs induced by neon injection has been observed in the EAST tokamak at overlapping q 95 space between large and small ELMs. With neon injection, pedestal density gradient shows a remarkable increase accompanied by some decrease of pedestal electron temperature, and consequently the pressure gradient increases moderately and edge bootstrap current has minimal change. Further experiment demonstrates that the occurrence of large ELMs after neon injection is highly correlated with the change in edge density. Linear peeling-ballooning stability analysis indicates that the large ELM case is more unstable than the type-III ELM case during the ELM transition. A scan of pedestal density gradient in linear stability analysis shows that the direct destabilizing effect of steep pedestal density gradient on peeling-ballooning instabilities via two-fluid effects could also facilitate the transition to large ELMs. In conclusion, these results could provide more insight into the role of pedestal density gradient on pedestal stability and ELM behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sustained edge-localized-modes suppression and radiative divertor with an impurity-driven instability in tokamak plasmas

We report simultaneous control of the large edge localized modes (ELMs) and divertor heat fluxes in a metal wall environment is a critical issue for steady-state operation of a tokamak fusion reactors. Here we report a sustained ELM suppression scenario achieved in the EAST tokamak compatible with radiative divertor using different seeding impurity species over a wide range of conditions. A low-n mode appears, as manifested by the oscillations of a radiation front near the X-point. This mode appears to drive strong particle transport and tungsten exhaust, which is essential to the maintenance of the ELM-stable state. We have developed a model to explain the mode excitation, by coupling the impurity radiative condensation instability to drift waves, which could explain some characteristics of the low-n mode well. The low-n mode may offer a new ELM-stable scenario compatible with radiative divertor for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A model investigation of the impact of lower hybrid wave scattering angle on current drive profile in EAST and Alcator C-Mod

Lower hybrid current drive (LHCD) is beneficial for developing a steady-state operation scenario in a tokamak. This paper conducts a modeling investigation to identify an optimum rotation angle of the initial lower hybrid perpendicular (to the background magnetic field $\overrightarrow{B}$) wavevector for best matching the experimental RF current profile. It is hypothesized that central RF power deposition widely observed in the present-day LHCD experiments arises from wave scattering by turbulence. In a standard model without considering such interactions, the predicted power deposition profile is generally broad with off-axis peaking, not in agreement with experimental observations. A heuristic approach is adopted by introducing a spectral broadening mechanism by modifying the initial orientation of the perpendicular wavevector. The ray-tracing/Fokker–Planck solver GENRAY/CQL3D is utilized within the python-based π-scope framework. A focus is given to identify the perpendicular wavenumber orientation angle with respect to the magnetic surface normal vector at the initial ray location. Our modeling study shows that rotating the perpendicular wavevector in such a way as to increase the initial poloidal component is effective in reproducing the centrally peaked current profile observed in normal shear plasmas on both EAST and C-Mod. These waves can readily be absorbed to the central plasma, which reduces the sensitivity of the power deposition profile to a slight change of the plasma condition. The same approach is also found to help broaden the off-axis power deposition profile in a reverse-shear EAST plasma, leading to a better agreement with the experiment. The results presented here suggest that spectral modification arising from edge density fluctuations in a tokamak may need to be considered in understanding wave propagation and absorption. A further experimental and theoretical/modeling study is vital as a reverse approach is adopted in this study. Finally, our work suggests that mitigation or control measures are critical for parasitic effects occurring on the first pass in a reactor regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhancement of edge turbulence concomitant with ELM suppression during boron powder injection in EAST

A reproducible, quasi-stationary edge localized mode (ELM)-suppressed scenario was obtained over a wide range of plasma parameters by continuous injection of boron (B) powder into an upper-single null discharge in the experimental advanced superconducting tokamak [Sun et al., Nucl. Fusion 61, 014002 (2021)]. This powder-induced ELM-absent regime is associated with an edge harmonic mode (EHM) that provides continuous particle exhaust to maintain constant density without confinement degradation and impurity accumulation, the latter of which is often observed in ELM-free regimes. A flow rate threshold of B powder injection, leading to a threshold intensity of the EHM, is necessary for full ELM suppression. The fundamental harmonic of the EHM exhibits a toroidal mode number n = 1. The mode is observable in the entire poloidal cross section with a peak near the upper X-point in an upper-single null configuration. Here, the EHM spans radially across the pedestal and scrape-off layer, peaking inside the separatrix. The EHM appears to be insensitive to q 95 , heating power, plasma toroidal rotation, and pedestal collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Screening of resonant magnetic perturbation fields assuming various plasma flow models

A recently updated version of the MARS-F code [Y. Q. Liu et al., Phys. Plasmas 7, 3681 (2000); L. Li et al., Phys. Plasmas 25, 082512 (2018); and G. L. Xia et al., Nucl. Fusion 59, 126035 (2019)] is utilized to numerically investigate the plasma screening effect on the applied resonant magnetic perturbation (RMP) field, assuming various equilibrium flow models, including the toroidal flow, the parallel flow and their combinations, and poloidal and toroidal projections of the parallel flow. A parallel equilibrium flow with a uniform radial profile is found to have no effect on plasma screening of the RMP field. A sheared parallel flow, however, does change plasma screening. The poloidal projection of the parallel flow weakens plasma screening in the resistive-inertial regime. The effect on the favorable average curvature regime is found, however, to be non-monotonic. With the increasing flow speed, the poloidal projection first weakens Glasser-Green-Johnson (GGJ)-screening. Further increase in the flow speed results in enhanced GGJ-screening again. This non-monotonic behavior is related to the perturbed parallel shielding current, which appears also off the mode rational surface at fast flow due to additional resonances between the RMP perturbation and the sound wave continuum. These results indicate that flow induced plasma screening to the RMP field can have complicated characteristics, which, in turn, can have implications on the RMP field penetration into the plasma in experiments for controlling the edge localized modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design trades for mission/DSN antenna supportability using TIGRAS

In this paper, we describe the capability and usage of TIGRAS. To illustrate the usage of TIGRAS, a hypothetical mission called Jovian Moon Explorer (JME) was created. We discuss how DSN supportability affects mission and spacecraft design and how a project such as JME can use TIGRAS in an early design phase to design for better DSN coverage for its future operation.

TIGRAS↗