Search NASA⌕ Search

SEARCH · Search NASA

Results for “gyro”

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 19 records

Benchmark studies of the gyro-Landau-fluid code and gyro-kinetic codes on kinetic ballooning modes

A Gyro-Landau-Fluid (GLF) 3 + 1 model has been recently implemented in BOUT++ framework, which contains full Finite-Larmor-Radius effects, Landau damping, and toroidal resonance [Ma et al., Phys. Plasmas 22, 055903 (2015)]. A linear global beta scan has been conducted using the JET-like circular equilibria (cbm18 series), showing that the unstable modes are kinetic ballooning modes (KBMs). In this work, we use the GYRO code, which is a gyrokinetic continuum code widely used for simulation of the plasma microturbulence, to benchmark with GLF 3 + 1 code on KBMs. To verify our code on the KBM case, we first perform the beta scan based on “Cyclone base case parameter set.” We find that the growth rate is almost the same for two codes, and the KBM mode is further destabilized as beta increases. For JET-like global circular equilibria, as the modes localize in peak pressure gradient region, a linear local beta scan using the same set of equilibria has been performed at this position for comparison. With the drift kinetic electron module in the GYRO code by including small electron-electron collision to damp electron modes, GYRO generated mode structures and parity suggest that they are kinetic ballooning modes, and the growth rate is comparable to the GLF results. However, a radial scan of the pedestal for a particular set of cbm18 equilibria, using GYRO code, shows different trends for the low-n and high-n modes. The low-n modes show that the linear growth rate peaks at peak pressure gradient position as GLF results. However, for high-n modes, the growth rate of the most unstable mode shifts outward to the bottom of pedestal and the real frequency of what was originally the KBMs in ion diamagnetic drift direction steadily approaches and crosses over to the electron diamagnetic drift direction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Predictions for Gyro-phase Drift in MDPX

This work assesses the feasibility of observing the gyro-phase drift in the Auburn Magnetized Dusty Plasma Experiment [MDPX, described by Thomas et al., Plasma Phys. Controlled Fusion 54, 124034 (2012)]. The gyro-phase drift arises when a dust grain does not instantaneously reach the in-situ-equilibrium grain charge during gyro-synchronous grain-charge modulation. Koepke et al. [J. Plasma Phys. 79, 1099 (2013)] first suggested using MDPX to observe the gyro-phase drift, and here we use a single-particle trajectory tracker with an iterative velocity solver, using a fixed timestep for grain motion and an adaptive time step for grain charging, to consider all relevant dust grain forces to assess gyro-phase drift arising from gradual inhomogeneity. Additionally, the semi-analytic theory developed by Walker et al. [J. Plasma Phys. 80, 395 (2014)] predicts dust grain motion in abrupt inhomogeneity for MDPX-relevant conditions. We compare three grain-charging models with each other and with the single-particle trajectory tracker and found to predict distinctly different trajectories depending on the treatment of neutral drag and flowing ions. The measurement thresholds for Particle Tracking Velocimetry permit gyro-phase drift detection in MDPX for the abrupt inhomogeneity, given sufficiently large enough UV photoelectron flux (fuv/[n e v the ] > 0.01) and low enough neutral gas pressure (less than one mTorr). The Orbit-Motion-Limited charge model and the charge models developed by Patacchini et al. [Phys. Plasmas 14, 062111 (2007)] and Gatti and Kortshagen [Phys. Rev. E 78, 046402 (2008)] can, in principle, be distinguished by gyro-phase drift in the abrupt inhomogeneity, but large magnetic fields, large UV photoelectron flux, and low neutral gas pressure are required. Gyro-phase drift for a gradual inhomogeneity in the ratio ne/ni, arising from the presence of a radial electric field, is predicted to be undetectable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

3D ion gyro-orbit heat load predictions for NSTX-U

Abstract High power tokamaks operate with divertor heat loads capable of destroying the plasma facing components (PFCs). High fidelity heat load predictions are necessary to ascertain the PFC state for design and during operation. Typical heat flux calculations are 2D, time invariant, and assume that power flows directly along the magnetic field lines (the optical approximation). These assumptions neglect the complex 3D geometries employed to protect the PFCs, the time varying nature of the plasma and PFC thermal state, and the helical trajectories of ions with finite Larmor radii (the gyro-orbit approximation). An integrated software framework, the heat flux engineering analysis toolkit (HEAT), was developed to generate time varying optical heat loads applied to real engineering computer aided design (CAD) (Looby et al 2022 Fusion Sci. Technol. 78 10–27). Recently, an ion-gyro orbit module has been added to HEAT. This module calculates the helical trajectories of ions as they gyrate about the magnetic field lines using kinetic theory macro-particles to accelerate the calculation. First, the new gyro-orbit module will be presented. Next, a comparison to existing research is performed. Finally, an analysis of the gyro-orbit heat loads for NSTX-U is presented for diverted discharges using the engineering CAD models utilized for PFC fabrication. Including these gyro-orbit effects can enhance the PFC performance by ‘smearing’ out the magnetic shadows associated with the castellated fish-scaled geometry. Simultaneously, the helical trajectories can degrade performance when they load narrow regions on edges and corners with high heat fluxes. Analysis of the trade-offs between these competing effects is included, and regions for further investigation are identified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gyro-viscosity and linear dispersion relations in pair-ion magnetized plasmas

A fluid theory has been developed by taking account of gyro-viscosity to study wave propagation characteristics in a homogeneous pair-ion magnetized plasma with a cylindrical symmetry. The exact dispersion relations derived by the Hankel-Fourier transformation are shown comparable with those observed in the experiment by Oohara and co-workers. The gyro-viscosity is responsible for the change in propagation characteristics of the ion cyclotron wave from forward to backward by suppressing the effect of the thermal pressure which normally causes the forward nature of dispersion. Although the experiment has been already explained by a kinetic theory by the present authors, the kinetic derivations are so involved because of exact particle orbits in phase space, finite Lamor radius effects, and higher order ion cyclotron resonances. The present fluid theory provides a simple and transparent structure to the dispersion relations since the gyro-viscosity is renormalized into the ion cyclotron frequency which itself indicates the backward nature of dispersion. The usual disadvantage of a fluid theory, which treats only fundamental modes of eigen-waves excited in a system and is not able to describe higher harmonics that a kinetic theory does, is compensated by simple derivations and clear picture based on the renormalization of the gyro-viscosity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Alfvén eigenmode stability and critical gradient energetic particle transport using the Trapped-Gyro-Landau-Fluid model

The Trapped-Gyro-Landau-Fluid (TGLF) transport model is a physically realistic and comprehensive theory based on a local quasilinear transport model fitted to linear and nonlinear GYRO gyrokinetic simulations [Staebler et al., Phys. Plasmas 14, 55909 (2007)]. This work presents the first use of the TGLF model to treat low-n Alfvén eigenmode (AE) stability and energetic particle (EP) transport. TGLF accurately recovers the local GYRO toroidicity-induced AE (TAE) and energetic particle mode (EPM) linear growth and frequency rates for a fusion alpha case. With a very high grid resolution, TGLF can quickly find the critical EP pressure gradient profile for stiff EP transport based on an AE linear threshold given the background thermal plasma profiles in DIII-D. The TGLF critical gradient profile using the recipe γ AE =0, that is the linear AE growth rate without additional driving rates from the background plasma gradients, matches the more expensive linear GYRO results with a single worst toroidal mode number n. TGLF can easily find the minimum critical gradient profile with testing multiple ns. From a database of runs using a newly developed TGLFEP code, a rough but insightful parametric “power law” scaling for critical EP beta is demonstrated. An important toroidal stabilization condition on the EP pressure gradient p EP /L$EP\atop{p}$ drive is isolated: R/L$EP\atop{p}$ > C R ~ 3, where L$EP\atop{p}$ is the EP pressure gradient length and R is the tokamak major radius. This paper also demonstrates that relaxation of the fixed slowing down EP profile shape approximation often used to find the critical EP density profile has little effect on the resulting EP transport. The single EP species critical gradient model is generalized to handle two EP species.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Phase Measurements of a 140-GHz Confocal Gyro-Amplifier

The phase stability of a 140-GHz, 1-kW pulsed gyro-amplifier system and the phase dependence on the cathode voltage were experimentally measured. To optimize the measurement precision, the amplifier was operated at 47 kV and 1 A, where the output power was ~ 30 W. The phase was determined to be stable both pulse-to-pulse and during each pulse, so far as the cathode voltage and electron beam current are constant. The phase variation with voltage was measured and found to be 130 ± 30° /kV, in excellent agreement with simulations. The electron gun used in this device is non-adiabatic, resulting in a steep slope of the beam pitch factor with respect to cathode voltage. This was discovered to be the dominant factor in the phase dependence on voltage. The use of an adiabatic electron gun is predicted to yield a significantly smaller phase sensitivity to voltage, and thus a more phase-stable performance. To our knowledge, these are the first phase measurements reported for a gyro-amplifier operating at a frequency above W-band.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Large gyro-orbit model of ion velocity distribution in plasma near a wall in a grazing-angle magnetic field

A model is presented for the ion distribution function in a plasma at a solid target with a magnetic field $\boldsymbol {B}$ inclined at a small angle, $\alpha \ll 1$ (in radians), to the target. Adiabatic electrons are assumed, requiring $\alpha \gg \sqrt {Zm_{e}/m_{i}}$ , where $m_{e}$ and $m_{i}$ are the electron and ion mass, respectively, and $Z$ is the charge state of the ion. An electric field $\boldsymbol {E}$ is present to repel electrons, and so the characteristic size of the electrostatic potential $\phi$ is set by the electron temperature $T_{e}$ , $e\phi \sim T_{e}$ , where $e$ is the proton charge. An asymptotic scale separation between the Debye length $\lambda _{D} = \sqrt {\epsilon _0 T_{{e}} / e^{2} n_{{e}} }$ , the ion sound gyro-radius $\rho _{s} = \sqrt { m_{i} ( ZT_{e} + T_{i} ) } / (ZeB)$ and the size of the collisional region $d_{c} = \alpha \lambda _{\textrm {mfp}}$ is assumed, $\lambda _{D} \ll \rho _{s} \ll d_{c}$ . Here $\epsilon _0$ is the permittivity of free space, $n_{e}$ is the electron density, $T_{i}$ is the ion temperature, $B= |\boldsymbol {B}|$ and $\lambda _{\textrm {mfp}}$ is the collisional mean free path of an ion. The form of the ion distribution function is assumed at distances $x$ from the wall such that $\rho _{s} \ll x \ll d_{c}$ , that is, collisions are not treated. A self-consistent solution of the electrostatic potential for $x \sim \rho _{s}$ is required to solve for the quasi-periodic ion trajectories and for the ion distribution function at the target. The large gyro-orbit model presented here allows to bypass the numerical solution of $\phi (x)$ and results in an analytical expression for the ion distribution function at the target. It assumes that $\tau =T_{i}/(ZT_{e})\gg 1$ , and ignores the electric force on the quasi-periodic ion trajectory until close to the target. For $\tau \gtrsim 1$ , the model provides an extremely fast approximation to energy–angle distributions of ions at the target. These can be used to make sputtering predictions.

Physics↗

A statistical study of gyro-averaging effects in a reduced model of drift-wave transport

Here, a statistical study of finite Larmor radius (FLR) effects on transport driven by electrostatic driftwaves is presented. The study is based on a reduced discrete Hamiltonian dynamical system known as the gyro-averaged standard map (GSM). In this system, FLR effects are incorporated through the gyro-averaging of a simplified weak-turbulence model of electrostatic fluctuations. Formally, the GSM is a modified version of the standard map in which the perturbation amplitude, K 0 , becomes K 0 J 0 ($\hat{p}$), where J 0 is the zeroth-order Bessel function and $\hat{p}$ s the Larmor radius. Assuming a Maxwellian probability density function (pdf) for $\hat{p}$ , we compute analytically and numerically the pdf and the cumulative distribution function of the effective drift-wave perturba- tion amplitude K 0 J 0 ($\hat{p}$). Using these results, we compute the probability of loss of confinement (i.e., global chaos), P c provides an upper bound for the escape rate, and that P t rovides a good estimate of the particle trapping rate. Lastly. the analytical results are compared with direct numerical Monte-Carlo simulations of particle transport.

42 ENGINEERING↗

Verification of GENE and GYRO with L-mode and I-mode plasmas in Alcator C-Mod

Here, verification comparisons are carried out for L-mode and I-mode plasma conditions in Alcator C-Mod. We compare linear and nonlinear ion-scale calculations by the gyrokinetic codes GENE and GYRO to each other and to the experimental power balance analysis. The two gyrokinetic codes' linear growth rates and real frequencies are in good agreement throughout all the ion temperature gradient mode branches and most of the trapped electron mode branches of the kyρs spectra at r/a = 0.65, 0.7, and 0.8. The shapes of the toroidal mode spectra of heat fluxes in nonlinear simulations are very similar for k y ρ s ≤ 0.5, but in most cases GENE has a relatively higher heat flux than GYRO at higher mode numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Asymptotic method for constructing a dynamic frequency response of a laser gyro

We consider an equation describing the dynamics of the phase of a beat signal at the output of a ring laser gyro in the presence of an alternating-sign frequency biasing in the form of a superposition of two meanders with different amplitudes and periods. It is assumed that the amplitude of one of the meanders is much larger than that of the other as well as than the half-widths of the static lock-in and the measured difference of eigenfrequencies of a ring resonator. A method is proposed for approximately calculating the beat frequency as a function of the measured difference between the eigenrequencies of the resonator. The result of applying this method represents analytically a recursive algorithm. Its implementation on a computer has made it possible to construct a dynamic frequency (output) response of a laser gyro with a biasing of given form. (laser applications and other topics in quantum electronics)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Phase Measurements of a 140-GHz Confocal Gyro-Amplifier

The phase stability of a 140-GHz, 1-kW pulsed gyro-amplifier system and the phase dependence on the cathode voltage were experimentally measured. To optimize the measurement precision, the amplifier was operated at 47 kV and 1 A, where the output power was ~30 W. The phase was determined to be stable both pulse-to-pulse and during each pulse, so far as the cathode voltage and electron beam current are constant. The phase variation with voltage was measured and found to be 130 ± 30 deg/kV, in excellent agreement with simulations. The electron gun used in this device is non-adiabatic, resulting in a steep slope of the beam pitch factor with respect to cathode voltage. This was discovered to be the dominant factor in the phase dependence on voltage. The use of an adiabatic electron gun is predicted to yield a significantly smaller phase sensitivity to voltage, and thus a more phase-stable performance. To our knowledge, these are the first phase measurements reported for a gyro-amplifier operating at a frequency above W-band.

Rosenzweig, Guy↗

Stability optimization of energetic particle driven modes in nuclear fusion devices: the FAR3d gyro-fluid code

The development of reduced models provide efficient methods that can be used to perform short term experimental data analysis or narrow down the parametric range of more sophisticated numerical approaches. Reduced models are derived by simplifying the physics description with the goal of retaining only the essential ingredients required to reproduce the phenomena under study. This is the role of the gyro-fluid code FAR3d, dedicated to analyze the linear and nonlinear stability of Alfvén Eigenmodes (AE), Energetic Particle Modes (EPM) and magnetic-hydrodynamic modes as pressure gradient driven mode (PGDM) and current driven modes (CDM) in nuclear fusion devices. Such analysis is valuable for improving the plasma heating efficiency and confinement; this can enhance the overall device performance. The present review is dedicated to a description of the most important contributions of the FAR3d code in the field of energetic particles (EP) and AE/EPM stability. FAR3d is used to model and characterize the AE/EPM activity measured in fusion devices as LHD, JET, DIII-D, EAST, TJ-II and Heliotron J. In addition, the computational efficiency of FAR3d facilitates performing massive parametric studies leading to the identification of optimization trends with respect to the AE/EPM stability. This can aid in identifying operational regimes where AE/EPM activity is avoided or minimized. This technique is applied to the analysis of optimized configurations with respect to the thermal plasma parameters, magnetic field configuration, external actuators and the effect of multiple EP populations. In addition, the AE/EPM saturation phase is analyzed, taking into account both steady-state phases and bursting activity observed in LHD and DIII-D devices. The nonlinear calculations provide: the induced EP transport, the generation of zonal structures as well as the energy transfer towards the thermal plasma and between different toroidal/helical families. Finally, FAR3d is used to forecast the AE/EPM stability in operational scenarios of future devices as ITER, CFETR, JT60SA and CFQS as well as possible approaches to optimization with respect to variations in the most important plasma parameters.

Alfv én Eigenmodes↗

Solitary Magnetic Structures Developed From Gyro-Resonance With Solar Wind Ions at Mars and Earth

We investigate solitary magnetic structures in the foreshock of Mars and Earth. The structures exhibit pulse-like magnetic field and density enhancements along with plasma heating and local solar wind (SW) ion reflection. The structures at Mars resemble the foreshock structures developed from ultra-low-frequency electromagnetic waves gyro-resonant with SW ions at Earth, and propagate toward the magnetosphere. We perform fully kinetic simulations to reproduce the solitary structures and the nonlinear evolution of ion distribution functions, illustrating their resonance with SW ions. The structures present self-induced foreshock turbulence that can have space weather effects. Our results advance the fundamental understanding of how SW interacts with planetary magnetospheres, and have potential impact on the current picture of planet-origin ion escape.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating gyro-viscosity in the Kelvin-Helmholtz instability by kinetic simulations

In the present paper, the finite-Larmor-radius (gyro-viscous) term [K. V. Roberts and J. B. Taylor, Phys. Rev. Lett. 8, 197–198 (1962)] is evaluated by using a full kinetic Vlasov simulation result of the Kelvin-Helmholtz instability (KHI). The velocity field and the pressure tensor are calculated from the high-resolution data of the velocity distribution functions obtained by the Vlasov simulation, which are used to approximate the Finite-Larmor-Radius (FLR) term according to Roberts and Taylor [Phys. Rev. Lett. 8, 197–198 (1962)]. The direct comparison between the pressure tensor and the FLR term shows an agreement. It is also shown that the anisotropic pressure gradient enhanced the linear growth of the KHI when the inner product between the vorticity of the primary velocity shear layer and the magnetic field is negative, which is consistent with the previous FLR-magnetohydrodynamic simulation result. This result suggests that it is not sufficient for reproducing the kinetic simulation result by fluid simulations to include the FLR term (or the pressure tensor) only in the equation of motion for fluid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The effects of main-ion dilution on turbulence in low q 95 C-Mod ohmic plasmas, and comparisons with nonlinear GYRO

Recent experiments on C-mod seeding nitrogen into ohmic plasmas with q 95 = 3.4 found that the seeding greatly reduced long-wavelength (ITG-scale) turbulence. The long-wavelength turbulence that was reduced by the nitrogen seeding was localized to the region of r/a ≈ 0.85, where the turbulence is well above marginal stability (as evidenced by Q i /Q GB $\gg$ 1). The nonlinear gyrokinetic code GYRO was used to simulate the expected turbulence in these plasmas, and the simulated turbulent density fluctuations and turbulent energy fluxes quantitatively agreed with the experimental measurements both before and after the nitrogen seeding. Unexpectedly, the intrinsic rotation of the plasma was also found to be affected by the nitrogen seeding, in a manner apparently unrelated to a change in the electron-ion collisionality that was proposed by other experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of the impact of parallel magnetic fluctuations on linear gyrokinetic stability in NSTX-U and verification of gyro-fluid models

In this work, we use the CGYRO gyrokinetic code to analyze two L- and one H-mode discharges from the National Spherical Torus Experiment (NSTX) and NSTX-Upgrade (NSTX-U) selected due to their different mix of ion-scale driftwaves, ion temperature gradient (ITG) mode and trapped electron mode (TEM), and electromagnetic instabilities, kinetic ballooning mode (KBM), and micro-tearing mode (MTM) in the plasma core. It is found that the effect of parallel magnetic fluctuations is strongly destabilizing to the unstable KBMs compared to calculations with only perpendicular magnetic fluctuations. Two discharges have a mix of ITG/TEM and MTMs that are predicted to be dominant instability across the plasma radius. The parallel magnetic fluctuations are found to have little effect on the MTM stability but are destabilizing to ITG/TEM modes. To test the validity of the gyro-fluid linear stability codes TGLF and GFS at low aspect ratio, a database of linear growth rates has been created using the CGYRO gyrokinetic code. The database is comprised of various parameter scans around a standardized set of NSTX-U core parameters. It contains a group of electrostatic cases and an electromagnetic group that includes the effects of perpendicular and parallel magnetic fluctuations. Comparing the results from the GFS and TGLF models, we find that GFS exhibits the best agreement with the database of CGYRO linear growth rates. Comparing the model results for the electromagnetic scans shows that GFS captures the effects of parallel magnetic fluctuations accurately, while the TGLF model does not, as it lacks sufficient perpendicular energy resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗