Search NASA⌕ Search

SEARCH · Search NASA

Results for “electron cyclotron waves”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Cause and impact of low-frequency chirping modes in DIII-D hybrid discharges

Significant variations in MHD activity and fast-ion transport are observed in the DIIID high-beta, steady-state hybrid discharges with a mixture of Electron Cyclotron (EC) waves and Neutral Beam Injection (NBI). When Electron Cyclotron Heating (ECH) or Current Drive (ECCD) is applied, Alfven Eigenmodes (AEs) are usually suppressed and replaced by low-frequency bursting modes. The analysis of a recently compiled database of hybrid discharges suggests that the change of the fast-ion pressure especially the perpendicular pressure is the main factor responsible for the instability transition. The lower ratio of fast-ion injection speed v inj to Alfven speed v alfven and slight drop of safety factor q min during ECCD may also facilitate the transition. The database shows that AEs mainly occur when the fast-ion fraction P f /P total is less than 0.53 and v inj /v alfven is greater than 0.50, while low-frequency bursting modes appear in the opposite regime. In this research, P f and P total are the central fast-ion pressure from classical prediction and total plasma pressure, respectively. The reason why the instability transition correlates with P f /P total and v inj /v alfven is that they can significantly modify the drive of low-frequency bursting modes and AEs. The explanation is supported by the observation that low-frequency bursting modes are rarely seen in the hybrids with NBI only, with EC waves and counter-NBI, or with high plasma density. A careful check of the low-frequency bursting modes suggests that they are mainly chirping (neoclassical) tearing modes [referred to as chirping (N)TM], i.e. the mode frequency firstly jumps up from the steady (N)TM frequency, then chirps down, and finally returns to the steady (N)TM frequency. Sometimes, the (N)TMs are fully stabilized and replaced with pure fishbones. Lastly, the resonance condition calculation and ’Kick’ model simulations suggest that (N)TMs and fishbones can interact through modification of the fast ion distribution in phase space, which influences the drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energetic Electron Precipitation Driven by Electromagnetic Ion Cyclotron Waves from ELFIN’s Low Altitude Perspective

We review comprehensive observations of electromagnetic ion cyclotron (EMIC) wave-driven energetic electron precipitation using data collected by the energetic electron detector on the Electron Losses and Fields InvestigatioN (ELFIN) mission, two polar-orbiting low-altitude spinning CubeSats, measuring 50-5000 keV electrons with good pitch-angle and energy resolution. EMIC wave-driven precipitation exhibits a distinct signature in energy-spectrograms of the precipitating-to-trapped flux ratio: peaks at >0.5 MeV which are abrupt (bursty) (lasting ~17 s, or ΔL ~ 0.56) with significant substructure (occasionally down to sub-second timescale). We attribute the bursty nature of the precipitation to the spatial extent and structuredness of the wave field at the equator. Multiple ELFIN passes over the same MLT sector allow us to study the spatial and temporal evolution of the EMIC wave - electron interaction region. Case studies employing conjugate ground-based or equatorial observations of the EMIC waves reveal that the energy of moderate and strong precipitation at ELFIN approximately agrees with theoretical expectations for cyclotron resonant interactions in a cold plasma. Using multiple years of ELFIN data uniformly distributed in local time, we assemble a statistical database of ~50 events of strong EMIC wave-driven precipitation. Most reside at L ~ 5 – 7 at dusk, while a smaller subset exists at L ~ 8 – 12 at post-midnight. The energies of the peak-precipitation ratio and of the half-peak precipitation ratio (our proxy for the minimum resonance energy) exhibit an L-shell dependence in good agreement with theoretical estimates based on prior statistical observations of EMIC wave power spectra. The precipitation ratio’s spectral shape for the most intense events has an exponential falloff away from the peak (i.e., on either side of ~ 1.45 MeV). It too agrees well with quasi-linear diffusion theory based on prior statistics of wave spectra. It should be noted though that this diffusive treatment likely includes effects from nonlinear resonant interactions (especially at high energies) and nonresonant effects from sharp wave packet edges (at low energies). Sub-MeV electron precipitation observed concurrently with strong EMIC wave-driven >1 MeV precipitation has a spectral shape that is consistent with efficient pitch-angle scattering down to ~ 200-300 keV by much less intense higher frequency EMIC waves at dusk (where such waves are most frequent). At ~100 keV, whistler-mode chorus may be implicated in concurrent precipitation. These results confirm the critical role of EMIC waves in driving relativistic electron losses. Nonlinear effects may abound and require further investigation.

79 ASTRONOMY AND ASTROPHYSICS↗

Studies of EC pre-ionization in DIII-D to support development of ITER plasma initiation

New experiments have been conducted at DIII-D to improve the physics understanding of plasma initiation assisted by electron cyclotron (EC) wave injection, allowing better extrapolation to ITER. This has been achieved by applying an EC pulse prior to start of the inductive plasma initiation (i.e. the generation of a loop voltage). A pre-plasma was formed during the EC pulse that was characterized in terms of the maximum density and temperature. Parametric scans were performed to study the influence of the EC injected power, EC injection angle, and pre-fill gas pressure on the pre-plasma creation process. These experiments showed that pre-ionized plasma of good quality can have a significant effect on the subsequent V loop induced plasma initiation process, i.e. a high density pre-plasma, increases the plasma current rise and speed at which ionization is achieved when the V loop is applied. Here, a good quality pre-plasma is one that achieved a significant degree of ionization, mainly obtained by providing sufficient ECH power in DIII-D of the order of 1 MW. It was found that a minimum EC power of 0.5 MW was required in DIII-D to create ionization, and this would scale to a minimum power of roughly 6.5 MW for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Full wave modeling of radio-frequency beams in tokamaks in the electron cyclotron frequency range

Simulation of full wave, without paraxial approximation, high-resolution solution of wave equations in frequency domain in the electron cyclotron resonance (ECR) frequency range for realistic Tokamak plasma parameters became possible by using recently formulated hybrid iterative algorithm [Svidzinski et al., Phys. Plasmas 25, 082509 (2018)] for numerically solving discretized wave equations. This approach combines time evolution and iterative relaxation techniques into iteration cycles. This algorithm is implemented in 2D code FullWave, solving wave equations in Tokamaks in cold and hot plasma models, and it has been tested in 3D full wave iterative RF beams simulation tool, which is presently being developed to model 3D ECRH RF beams in fusion devices using dynamic grid adaptation. The results of 2D full wave modeling, assuming specified toroidal mode number, of ECRH RF beams in DIII-D plasma, performed in the cold and hot plasma models for outboard and top launch scenarios using FullWave are presented. Nonlocal hot plasma response model, based on accurate numerical solution of linearized Vlasov equation, is used to model beam propagation and absorption in the 2nd electron cyclotron harmonic region. Demonstration of capability of the hybrid iterative algorithm to model ECRH RF beams in 3D is made by simulating a substantial part of realistic beam in DIII-D, launched from outboard side of the machine. All relevant physics of RF beam propagation, most of which is not captured in paraxial approximation, such as beam's divergence, interference between the X and O modes in the beam, X-O mode conversion, beam splitting into the X and O mode beams, transformation of beam's cross section, and absorption at the 2nd electron cyclotron harmonic, is captured in the simulations. A numerical technique to find an optimal beam polarization at the launcher to launch a nearly pure X or O mode beam in plasma is developed and tested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Whistler waves generated by nongyrotropic and gyrotropic electron beams during asymmetric guide field reconnection

Using a two-dimensional particle-in-cell simulation of asymmetric reconnection with a guide field whose strength is 0.3 times the reconnecting magnetic field, we study electron distribution functions and wave intensities in the diffusion region, focusing on the electron diffusion region (EDR). Wave activities with frequencies below the electron cyclotron frequency are observed, and these are whistler waves propagating almost anti-parallel to the magnetic field. The waves are concentrated near the magnetospheric separatrix away from the X line, but the wave activity also spreads through the EDR near the X line. The reconnection outflows are asymmetric in the outflow direction in the magnetospheric side, and the wave intensity is stronger in the side of the faster electron outflow. We study the whistler waves using the fast Fourier transform, analyses of electron velocity distribution functions, and the dispersion solver calculation. Along the magnetospheric separatrix in the stronger outflow side, highly anisotropic electron beams exist with super-Alfvénic drift speeds. The dispersion analysis shows that there are two modes: a temperature anisotropy mode and a beam mode. Outside the EDR, the whistler wave intensity is highest near the separatrix, but the wave intensity decreases if we move away from the separatrix toward the magnetic neutral line because of the increase in the electron population near zero parallel velocity. In the EDR, in the velocity plane perpendicular to the magnetic field, ring/crescent electron distribution functions are observed. Near the X-line, the wave power is enhanced where nongyrotropic electrons contribute to increase the perpendicular temperature anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resolving ECRH deposition broadening due to edge turbulence in DIII-D

Microwave heat pulse propagation experiments have demonstrated a correlation between millimeter-scale turbulence and deposition profile broadening of electron cyclotron (EC) waves on the DIII-D tokamak. In a set of discharges in DIII-D, a variation in edge density fluctuations on the mm-scale is associated with 40%–150% broader deposition profiles, expressed in terms of normalized minor radius, as compared with equilibrium ray tracing. The 1D power profile is determined from transport analysis of the electron temperature response to EC power modulation using perturbative analysis with a square wave power modulation at 20–70 Hz, producing a series of Fourier harmonics that are fit collectively to resolve transport. Fitting an integrated heat flux expressed in the Fourier basis of the modulation to diffusive, convective, and coupled transport terms in a linear model can resolve the broadened EC deposition width from the power perturbation to resolve a broadening in each case. The best fit degree of beam broadening observed scales approximately linearly with the Doppler backscattering measured fluctuation level in the steep gradient region. Quantifying the effect of edge fluctuation broadening on EC current drive power needs of future devices will require 3D full-wave codes that can be validated on the current generation of machines. These DIII-D experiments provide a quantitative measure of fluctuation effects and a dataset to benchmark full-wave simulations that can model and eventually predict nonlinear effects neglected by 1D equilibrium beam and ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Giant Undulations Driven by Pitch‐Angle Scattering of Time Domain Structures Modulated by Plasmapause Surface Wave

Abstract Plasmapause surface waves (PSWs) near the plasmapause boundary are regarded to be the magnetospheric source of ionospheric auroral giant undulations (GUs) located at the equatorward boundary of diffuse aurora. However, the observational evidence of wave‐particle interaction connecting PSWs and GUs is absent. In this letter, we demonstrate GUs are driven by pitch‐angle scattering of time domain structures modulated by the PSWs, based on the conjugated ionospheric and magnetospheric observations. Specifically, ionospheric GUs are lighted by the pitch‐angle scattering of <1 keV thermal electron and ions and energetic ions with energy up to dozens of keV near the plasmapause. Further, the total fluxes during one PSW period and energy of scattered electron and ions determine the size and luminosity of GUs. Our research provides observational evidence that PSWs cause periodic electron precipitation via modulating the time domain structures rather than the previously predicted chorus or electron cyclotron harmonic waves.

Zhou, Yi‐Jia [Weihai Institute for Interdisciplina↗

Neural-Network-Enhanced COTSIM: Advancing Predictive Capabilities for Fast DIII-D Simulations

Sustaining fusion reactions in tokamaks requires heating plasma to thermonuclear temperatures while maintaining confinement and stability. Neutral beam injection (NBI) provides heating, current drive, torque, and fueling, while electron cyclotron (EC) waves are widely used for heating and current drive; together, these actuators shape the plasma current, temperature, and density profiles. The control-oriented tokamak simulator (COTSIM), a predictive, control-oriented code, has been enhanced with neural-network surrogates for transport and sources. Turbulent transport is predicted by MMMnet—a neural-network version of the updated multimode model (MMM 9.0.10)—with significantly reduced computation time relative to MMM; neoclassical transport follows the Chang–Hinton model. NUBEAMnet, a surrogate of the Monte Carlo NUBEAM module, predicts beam-driven heating, current, and torque. EC heating and current drive use a control-oriented, empirically scaled source model; plasma resistivity follows the Spitzer formulation; bootstrap current uses the Sauter model. Equilibrium is computed using both prescribed and fixed-boundary solvers (FBSs), and the pedestal structure is modeled with an empirical pedestal model. For a representative DIII-D discharge, COTSIM predicts electron and ion temperature and safety-factor profiles in close agreement with TRANSP predictive and interpretive simulations while extending predictions through the pedestal region to the plasma edge (versus 80% of the minor radius in TRANSP). Furthermore, the equivalent COTSIM simulation runs in under 3 min compared to about 2 h for TRANSP, enabling rapid scenario planning, optimization of tokamak operation, and between-pulse control design.

Control-oriented tokamak simulator (COTSIM)↗

Quantifying electron cyclotron power deposition broadening in DIII-D and the potential consequences for the ITER EC system

Abstract The injection of electron cyclotron (EC) waves fulfills a number of important tasks in nuclear fusion devices for which detailed knowledge of the spatial power deposition profile is critical. This deposition profile is commonly determined using forward models such as beam or ray tracing. Recent numerical and experimental studies have shown that small-angle scattering of the EC beam as it passes through the turbulent plasma edge can cause significant broadening of the effective deposition profile, leading to considerable underestimation of the deposition width by forward methods. However, traditional inverse methods to determine the deposition profile from measurements overestimate the deposition profile width due to transport broadening. In this work, we implement three novel methods to resolve the EC power deposition profile from measurements that counteract transport broadening by simultaneously resolving transport and power deposition. We validate their assumptions and compare the results from these methods to the traditional break-in-slope method as well as to the TORAY ray-tracing code in a set of DIII-D discharges spanning five different confinement modes. We show that the four different inverse methods, novel and established, paint a consistent picture of deposition broadening. Specifically, we show that the measured power deposition profile is between 1.6 and 3.6 times wider than the TORAY profiles. Moreover, we show the considerable consequences that this level of broadening can have for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Designing a validation experiment for radio frequency condensation

Abstract Theoretical studies have suggested that nonlinear effects can lead to ‘radio frequency (RF) condensation’, where an initially broad current profile can coalesce in islands when they reach sufficient width. In suitable conditions, RF condensation can ‘self-focus’ the driven current to the center of an island, improving stabilization efficiency and reducing control complexity. In unsuitable conditions, the effect can prematurely deplete the RF energy before it reaches the island center, impairing stabilization. It is predicted that the RF condensation effect can significantly impact reactor-scale tokamaks. This paper presents a set of simulations investigating the conditions under which RF condensation might be encountered in present-day tokamaks. For concreteness, the calculations use equilibrium reconstructions for two shots from DIII-D and AUG. The Current Condensation Amid Magnetic Islands (OCCAMI) simulation code has been used for this investigation. The code takes as its input a numerically specified axisymmetric EFIT equilibrium solution, and it perturbatively constructs a 3D field with an island embedded at the appropriate rational surface. In the OCCAMI code, the GENRAY code is used for ray tracing and for calculating the power deposition along a ray trajectory, and GENRAY is coupled self-consistently to a solution of the thermal diffusion equation in the island. The simulation results described in the paper illuminate the conditions required for experimental validation of the theory of RF condensation. The simulations also provide an explanation of why the effect was not noticed in experiments prior to the publication of theoretical papers on the subject.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of runaway electron driven whistler wave instability experiments

Data acquired on the DIII-D tokamak were analyzed. The data are from experiments that were conducted to study an instability that is driven unstable by intense populations of electrons with MeV energies that are known as runaway electrons. The instability is a type of plasma wave called a whistler wave that occurs at frequencies above the ion cyclotron frequency but well below the electron cyclotron frequency. The waves were measured by magnetic fluctuation coils that are embedded in the DIII-D vacuum vessel wall. After upgrades to this diagnostic were completed, new experiments were conducted on July 13, 2020 in order to measure the toroidal mode number of the whistler waves and to extend the frequency of the detected waves. Through the use of mixers, instability between 600-700 MHz was detected. (The initial experiments only measured up to 200 MHz.) Analysis of the data was led by Hari Choudhury, a PhD student at Columbia University. Mr. Choudhury has submitted two papers for publication that include contributions by UC Irvine (UCI) Professor Heidbrink and his graduated PhD student Genevieve DeGrandchamp: “Detailed Characterization of Runaway Electron Driven Whistler Waves in Low-Density DIII-D Discharges” and “First Demonstration of Resonant Pitch-Angle Scattering of Relativistic Electrons by Externally-Launched Helicon Waves.” The first paper, which has been submitted to Physics of Plasmas, has significant contributions to both the data and the interpretation by UCI scientists. In contrast, UCI contributions to the second paper, which has been submitted to Physical Review Letters, are relatively minor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Two-fluid model of rf current condensation in magnetic islands

The stabilization of tearing modes with rf waves is subject to a nonlinear effect, termed rf current condensation, that has the potential to greatly enhance and localize current driven within magnetic islands. Here we extend previous investigations of this effect with a two fluid model that captures the balance of diffusive and thermal equilibration processes within the island. We show that the effective power and resulting strength of the condensation effect can be greatly enhanced by avoiding collisional heat loss to the ions. The relative impact of collisions on the overall power balance within the island depends on the ratio of the characteristic diffusion timescale and the electron–ion equilibration time, rather than the latter alone. Although relative heat loss to ions increases with island size, the heating efficiency does as well. In particular, we show that the latter safely dominates for large deposition profiles, as is typically the case for lower hybrid current drive. This supports the possibility of passive stabilization of neoclassical tearing modes without the precise aiming of the rf waves required for electron cyclotron current drive stabilization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Case Study of Nonresonant Mode 3-s ULF Waves Observed by MMS

The nature of the 3-s ultralow frequency (ULF) wave in the Earth's foreshock region and the associated wave-particle interaction are not yet well understood. In this work, we investigate the 3-s ULF waves using Magnetospheric Multiscale (MMS) observations. By combining the plasma rest frame wave properties obtained from multiple methods with the instability analysis based on the velocity distribution in the linear wave stage, the ULF wave is determined to be due to the ion/ion nonresonant mode instability. The interaction between the wave and ions is analyzed using the phase relationship between the transverse wave fields and ion velocities and using the longitudinal momentum equation. During the stage when ULF waves have sinusoidal waveforms up to |dB|/|B 0 | ~ 3, where dB is the wave magnetic field and B0 is the background magnetic field, the wave electric fields perpendicular to B 0 do negative work to solar wind ions; along B 0 , a longitudinal electric field develops, but the V × B force is stronger and leads to solar wind ion deceleration. During the same wave stage, the backstreaming beam ions gain energy from the transverse wave fields and get deceleration along B 0 by the longitudinal electric field. The ULF wave leads to electron heating, preferentially in the direction perpendicular to the local magnetic field. Secondary waves are generated within the ULF waveforms, including whistler waves near half of the electron cyclotron frequency, high-frequency electrostatic waves, and magnetosonic whistler waves. The work improves the understanding of the nature of 3-s ULF waves and the associated wave-particle interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Differentiating Between Simultaneous Loss Drivers in Earth's Outer Radiation Belt: Multi‐Dimensional Phase Space Density Analysis

Abstract We analyzed the contribution of electromagnetic ion cyclotron (EMIC) wave driven electron loss to a flux dropout event in September 2017. The evolution of electron phase space density (PSD) through the dropout showed the formation of a radially peaked PSD profile as electrons were lost at high L* , resembling distributions created by magnetopause shadowing. By comparing 2D Fokker Planck simulations of pitch angle diffusion to the observed change in PSD, we found that the μ and K of electron loss aligned with maximum scattering rates at dropout onset. We conclude that, during this dropout event, EMIC waves produced substantial electron loss. Because pitch angle diffusion occurred on closed drift paths near the last closed drift shell, no radial PSD minimum was observed. Therefore, the radial PSD gradients resembled solely magnetopause shadowing loss, even though the local pitch angle scattering produced electron losses of several orders of magnitude of the PSD.

Staples, F. A.↗

Model predictive current profile control in tokamaks by exploiting spatially moving electron cyclotron current drives

The plasma control systems in next-generation tokamaks like ITER will balance competing control objectives to achieve the desired level of performance in advanced scenarios while preventing magnetohydrodynamic instabilities and disruptions. During normal tokamak operation, the points of incidence of the electromagnetic waves generated by the electron cyclotron heating and current drives (EC H&CDs) are usually fixed in space. However, the points of incidence can be modified in real-time by changing the angles of the mirrors that reflect the EC H&CD waves. Altering the points of incidence, in turn, varies the ability of the plasma control system to regulate a plasma property. For instance, changing the EC H&CD wave incidence location may place the power demands necessary to achieve a particular plasma target within saturation limits. Therefore, using the EC H&CD deposition location, which is related to the EC H&CD mirror angle, as a supplementary controllable variable may facilitate access to a given target scenario. However, active scenario-control algorithms have not been designed so far to fully exploit this capability in real time. In this work, a model predictive controller that can handle actuation locations as control inputs is developed. In particular, the controller is designed to regulate both the auxiliary powers and the EC H&CD deposition locations in a pre-defined optimal sense to achieve the control objective of attaining and sustaining a target current profile. Further, the proposed controller is tested for a DIII-D tokamak scenario in nonlinear simulations using the Control Oriented Transport SIMulator (COTSIM).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient electron cyclotron current drive regime for plasma current start-up in fusion reactors

We report an accessibility enhanced efficient fundamental X-mode electron cyclotron heating (ECH) current start-up regime was identified for a reactorlike toroidal magnetic field range which has more than 100 times higher current drive efficiency compared to more conventional ECH methods for the relevant start-up temperature range. Very high current drive efficiency is possible due to the strong cyclotron interaction only with unidirectional passing electrons constrained by the wave accessibility conditions. This efficient electron cyclotron current drive regime may help facilitate the design of innovative economical solenoid-free tokamak fusion reactor systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First Demonstration of Resonant Pitch-Angle Scattering of Relativistic Electrons by Externally Launched Helicon Waves

Helicon waves (a.k.a whistler waves) satisfying the normal wave-particle cyclotron resonance are observed to limit the growth and maximum energy of relativistic electrons (REs) in low-density Ohmic DIII-D tokamak plasmas. Following the application of helicon waves, pitch-angle scattering of high-energy REs causes an increase in both synchrotron and electron-cyclotron emissions. The hard x-ray emission, a proxy for the RE population, ceases to grow. Energy-resolved hard x-ray measurements also show a striking decrease in the number of high-energy REs (above the resonance at approximately 8MeV) to below the noise floor and an increase in low-energy (∼ 4 MeV) REs. This occurs despite the toroidal electric field remaining high enough to drive exponential RE growth in the absence of helicon waves. Furthermore, these results open new directions for limiting the maximum energy of RE populations in laboratory and fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗