Search NASA⌕ Search

SEARCH · Search NASA

Results for “LANDAU DAMPING”

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 217 records · Page 12

Theory and simulation of the beam cyclotron instability.

A theory of plasma beam cyclotron instability is developed on the basis of computer simulation experiments. The theory holds that at certain turbulence levels, electron cross-field diffusion which supresses the electron gyroresonances is created by turbulent wave-particle interactions in a plasma beam after a period of quasi-linear exponential development of turbulence. The stabilizing effect of Landau ion damping is noted. The behavior of cold and hot ions is discussed.

Lampe, M.↗

The theory of magnetohydrodynamic wave generation by localized sources. II - Collisionless dissipation of wave packets

The dispersion equation of Barnes (1966) is used to study the dissipation of asymptotic wave packets generated by localized periodic sources. The solutions of the equation are linear waves, damped by Landau and transit-time processes, in a collisionless warm plasma. For the case of an ideal MHD system, most of the waves emitted from a source are shown to cancel asympotically through destructive interference. The modes transporting significant flux to asymptotic distances are found to be Alfven waves and fast waves with theta (the angle between the magnetic field and the characteristics of the far-field waves) of about 0 and about pi/2.

Collins, William↗

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On Stochastic Heating and Its Phase-space Signatures in Low-beta Kinetic Turbulence

We revisit the theory of stochastic heating of ions and investigate its phase-space signatures in kinetic turbulence of relevance to low- β portions of the solar wind. In particular, we retain a full scale-dependent approach in our treatment, and we explicitly consider the case in which electric-field fluctuations can be described by a generalized Ohm's law that includes Hall and thermoelectric effects. These two electric-field terms provide the dominant contributions to stochastic ion heating when the ion-Larmor scale is much smaller than the ion skin depth, ρ i << d i , which is the case at β << 1. Employing well-known spectral scaling laws for Alfvén-wave and kinetic-Alfvén-wave turbulent fluctuations, we obtain scaling relations characterizing the field-perpendicular particle-energization rate and energy diffusion coefficient associated with stochastic heating in these two regimes. Phase-space signatures of ion heating are then investigated using three-dimensional hybrid-kinetic simulations of continuously driven Alfvénic turbulence at low β (namely, β i = β e = 0.3 and β i = β e = 1/9). In these simulations, energization of ions parallel to the magnetic field is subdominant compared to its perpendicular counterpart ( Q ∥,i << Q ⊥,i ), and the fraction of turbulent energy that goes into ion heating is ≈75% at β i = 0.3 and ≈40% at β i ≃ 0.1. The phase-space signatures of ion energization are consistent with Landau-resonant collisionless damping and a ( β -dependent) combination of ion-cyclotron and stochastic heating. We demonstrate good agreement between our scale-dependent theory and various signatures associated with the stochastic portion of the heating. We discuss briefly the effect of intermittency on stochastic heating and the implications of our work for the interpretation of stochastic heating in solar-wind spacecraft data.

79 ASTRONOMY AND ASTROPHYSICS↗

Are Ring Current Ions Lost in Electromagnetic Ion Cyclotron Wave Dispersion Relation?

Electromagnetic ion cyclotron (EMIC) waves are widely observed in the inner and outer magnetosphere, at geostationary orbit, at high latitudes along the plasmapause, and at the ionospheric altitudes. Interaction of the Ring Current (RC) ions and EMIC waves causes ion scattering into the loss cone and leads to decay of the RC, especially during the main phase of storms when the RC decay times of about one hour or less are observed. The oblique EMIC waves damp due to Landau resonance with the thermal plasmaspheric electrons, and subsequent transport of the dissipating wave energy into the ionosphere below causes an ionosphere temperature enhancement. Induced scattering of these waves by the plasmaspheric thermal ions leads to ion temperature enhancement, and forms a so-called hot zone near the plasmapause where the temperature of core plasma ions can reach tens of thousands of degrees. Relativistic electrons in the outer radiation belt also interact well with the EMIC waves, and during the main and/or recovery phases of the storms these electrons can easily be scattered into the loss cone over a time scale from several hours to a day. The plasma density distribution in the magnetosphere and the ion content play a critical role in EMIC wave generation and propagation, but the wave dispersion relation in the known RC-EMIC wave interaction models is assumed to be determined by the thermal plasma distribution only. In these models, the modification of the EMIC wave dispersion relation caused by the RC ions is not taken into account, and the RC ions are only treated as a source of free energy in order to generate EMIC waves. At the same time, the RC ions can dominate the thermal magnetospheric content in the night MLT sector at great L shells during the main and/or recovery storm phase. In this study, using our self-consistent RC-EMIC wave model [Khazanov et al., 2006], we simulate the May 1998 storm in order to quantify the global EMIC wave redistribution caused by taking into account the RC ions in the EMIC wave dispersion relation. The dramatic wave pattern redistribution is observed in the postdusk-predawn MLT sector (night sector) for L greater than 5. We found the intense EMIC waves (about a few nT) there during the main and early recovery phases of the storm. The observed wave generation in this sector is caused by taking into account the EMIC wave dispersion change due to the RC ions. There are no waves at these locations in our model if the RC ions are taken into account in the wave growth rate only, and the wave dispersion relation is only governed by the thermal plasmaspheric model.

Khazanov, G. V.↗

Effect of Ring Current Ions on Electromagnetic Ion Cyclotron Wave Dispersion Relation

Electromagnetic ion cyclotron (EMIC) waves are widely observed in the inner and outer magnetosphere, at geostationary orbit, at high latitudes along the plasmapause, and at the ionospheric altitudes. Interaction of the Ring Current (RC) ions and EMIC waves causes ion scattering into the loss cone and leads to decay of the RC, especially during the main phase of storms when the RC decay times of about one hour or less are observed. The oblique EMIC waves damp due to Landau resonance with the thermal plasmaspheric electrons, and subsequent transport of the dissipating wave energy into the ionosphere below causes an ionosphere temperature enhancement. Induced scattering of these waves by the plasmaspheric thermal ions leads to ion temperature enhancement, and forms a so-called hot zone near the plasmapause where the temperature of core plasma ions can reach tens of thousands of degrees. Relativistic electrons in the outer radiation belt also interact well with the EMIC waves, and during the main and/or recovery phases of the storms these electrons can easily be scattered into the loss cone over a time scale from several hours to a day. The plasma density distribution in the magnetosphere and the ion content play a critical role in EMIC wave generation and propagation, but the wave dispersion relation in the known RC-EMIC wave interaction models is assumed to be determined by the thermal plasma distribution only. In these models, the modification of the EMIC wave dispersion relation caused by the RC ions is not taken into account, and the RC ions are only treated as a source of free energy in order to generate EMIC waves. At the same time, the RC ions can dominate the thermal magnetospheric content in the night MLT sector at great L shells during the main and/or recovery storm phase. In this study, using our self-consistent RC-EMIC wave model [Khazanov et al., 2006], we simulate the May 1998 storm in order to quantify the global EMIC wave redistribution caused by taking into account the RC ions in the EMIC wave dispersion relation. The dramatic wave pattern redistribution is observed in the postdusk-predawn MLT sector (night sector) for L greater than 5. We found the intense EMIC waves (about a few nT) there during the main and early recovery phases of the storm. The observed wave generation in this sector is caused by taking into account the EMIC wave dispersion change due to the RC ions. There are no waves at these locations in our model if the RC ions are taken into account in the wave growth rate only, and the wave dispersion relation is only governed by the thermal plasmaspheric model.

Gamayunov, K. V.↗

Self-Consistent Ring Current/Electromagnetic Ion Cyclotron Waves Modeling

The self-consistent treatment of the RC ion dynamics and EMIC waves, which are thought to exert important influences on the ion dynamical evolution, is an important missing element in our understanding of the storm-and recovery-time ring current evolution. For example, the EMlC waves cause the RC decay on a time scale of about one hour or less during the main phase of storms. The oblique EMIC waves damp due to Landau resonance with the thermal plasmaspheric electrons, and subsequent transport of the dissipating wave energy into the ionosphere below causes an ionosphere temperature enhancement. Under certain conditions, relativistic electrons, with energies 21 MeV, can be removed from the outer radiation belt by EMIC wave scattering during a magnetic storm. That is why the modeling of EMIC waves is critical and timely issue in magnetospheric physics. This study will generalize the self-consistent theoretical description of RC ions and EMIC waves that has been developed by Khazanov et al. [2002, 2003] and include the heavy ions and propagation effects of EMIC waves in the global dynamic of self-consistent RC - EMIC waves coupling. The results of our newly developed model that will be presented at the meeting, focusing mainly on the dynamic of EMIC waves and comparison of these results with the previous global RC modeling studies devoted to EMIC waves formation. We also discuss RC ion precipitations and wave induced thermal electron fluxes into the ionosphere.

Khazanov, G. V.↗

Parametric instabilities of parallel-propagating Alfven waves: Some analytical results

We consider the stability of a circularly polarized Alfven wave (the pump wave) which propagates parallel to the ambient magnetic field. Only parallel-propagating perturbations are considered, and we ignore dispersive effects due to the ion cyclotron frequency. The dissipationless MHD equations are used throughout; thus possibibly important effects arising from Landau and transit time damping are omitted. We derive a series of analytical approximations to the dispersion relation using A = (Delta B/B(sub O))(exp 2) as a small expansion parameter; Delta B is the pump amplitude, and B(sub O) is the ambient magnetic field strength. We find that the plasma beta (the square of the ratio of the sound speed to the Alfven speed) plays a crucial role in determining the behavior of the parametric instabilities of the pump. If 0 less than beta less than 1 we find the familiar result that the pump decays into a forward propagating sound wave and a backward propagating Alfven wave with maximum growth rate gamma(sub max) varies A(sup 1/2), but beta cannot be too close to 0 or to 1. If beta approx. 1, we find gamma(sub max) varies A(sup 3/4), if beta greater than 1, we find gamma(sub max) varies A(sup 3/2), while if beta approx. 0, we obtain gamma(sub max) varies A(sup 1/3); moreover, if beta approx. 0 there is a nearly purely growing instability. In constrast to the familiar decay instability, for which the backward propagating Alfven wave has lower frequency and wavenumber than the pump, we find that if beta greater than or approx. equal to 1 the instability is really a beat instability which is dominated by a transverse wave which is forward propagating and has frequency and wavenumber which are nearly twice the pump values. Only the decay instability for 0 less than beta less than 1 can be regarded as producing two recognizable normal modes, namely, a sound wave and an Alfven wave. We discuss how the different characteristics of the instabilities may affect the evolution of Alfven waves in the solar wind. However, for a solar wind in which beta approx. 1 the growth times of the instabilities are probably too long for these instabilities to have an appreciable effect inside 1 AU.

Jayanti, V.↗

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↗

Self-sustaining sound in collisionless, high- β plasma

Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $δn/n_0≳2/β$ (where $β\gg1$ is the ratio of thermal to magnetic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly to reduce the pressure anisotropy by pitch-angle scattering and trapping particles, respectively, thereby impeding the maintenance of Landau resonances that enable such waves’ otherwise potent collisionless damping. The result is wave dynamics that evince a weakly collisional plasma: the ion distribution function is near-Maxwellian, the field-parallel flow of heat resembles its Braginskii form (except in regions where large-amplitude magnetic mirrors strongly suppress particle transport), and the relations between various thermodynamic quantities are more ‘fluid-like’ than kinetic. A nonlinear fluctuation–dissipation relation for self-sustaining IAWs is obtained by solving a plasma-kinetic Langevin problem, which demonstrates suppressed damping, enhanced fluctuation levels and weakly collisional thermodynamics when IAWs with $δn/n_0≳2/β$ are stochastically driven. We investigate how our results depend upon the scale separation between the wavelength of the IAW and the Larmor radius of the ions, and discuss briefly their implications for our understanding of turbulence and transport in the solar wind and the intracluster medium of galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Langevin dynamics of generalized spins as SU( N ) coherent states

Classical models of spin systems traditionally retain only the dipole moments, but a quantum spin state will frequently have additional structure. Spins of magnitude S have N = 2⁢S + 1 levels. Alternatively, the spin state is fully characterized by a set of N 2 - 1 local physical observables, which we interpret as generalized spin components. For example, a spin with S = 1 has three dipole components and five quadrupole components. Further, these components evolve under a generalization of the classical Landau-Lifshitz dynamics, which can be extended with noise and damping terms. In this paper, we reformulate the dynamical equations of motion as a Langevin dynamics of SU(N) coherent states in the Schrödinger picture. This viewpoint is especially useful as the basis for an efficient numerical method to sample spin configurations in thermal equilibrium and to simulate the relaxation and driven motion of topological solitons. To illustrate the approach, we simulate a nonequilibrium relaxation process that creates CP 2 skyrmions, which are topological defects with both dipole and quadrupole character.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Control of magnon–photon coupling by spin torque

We demonstrate the influence of damping and field-like torques in the magnon–photon coupling process by classically integrating the generalized Landau–Lifshitz–Gilbert equation with RLC equation in which a phase correlation between dynamic magnetization and microwave current through combined Ampère and Faraday effects are considered. We show that the gap between two hybridized modes can be controlled in samples with damping parameter in the order of 10 –3 by changing the direction of the dc current density J if a certain threshold is reached. Furthermore, our results suggest that an experimental realization of the proposed magnon–photon coupling control mechanism is feasible in yttrium iron garnet/Pt hybrid structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Photogalvanic Effects in Surface States of Topological Insulators under Perpendicular Magnetic Fields

We present a theoretical study of the nonlinear magneto-optical shift conductivity in the surface states of the prototypical topological insulator Bi$_2$Se$_3$ under a perpendicular quantizing magnetic field. By describing the electronic states as Landau levels and using a perturbative approach, we derive the microscopic expression for the shift conductivity $\sigma^{(2);\alpha\beta\gamma}(-\omega,\omega)$, where $\alpha,\beta,\gamma=\pm$ stand for the circular polarization of light and $\omega$ is the light frequency; the spectra are further decomposed into contributions from the interband and intraband optical transitions, for which the selection rules are identified. Considering that the system possesses $C_3$ point group of symmetry, the nonzero components of the conductivity tensor are $\sigma^{(2);-++}=[\sigma^{(2);+--}]^\ast$. Therefore, a pure circularly polarized light generates zero shift current. In the clean limit, the conductivities are nonzero only for discrete photon energies because of the discrete Landau levels and energy conservation, and they become Lorentzian lineshapes with the inclusion of damping, which relaxes the condition of energy conservation. The dependence of the spectra on the damping parameters, the magnetic fields, and the chemical potentials is investigated in detail. Our results reveal that the shift current is highly tunable by the chemical potential and the magnetic field. Furthermore, these results underscore the potential of topological insulators for tunable, strong nonlinear magneto-optical applications.

Landau level↗

Landau Modes are Eigenmodes of Stellar Systems in the Limit of Zero Collisions

Abstract We consider the spectrum of eigenmodes in a stellar system dominated by gravitational forces in the limit of zero collisions. We show analytically and numerically using the Lenard–Bernstein collision operator that the Landau modes, which are not true eigenmodes in a strictly collisionless system (except for the Jeans unstable mode), become part of the true eigenmode spectrum in the limit of zero collisions. Under these conditions, the continuous spectrum of true eigenmodes in a collisionless system, also known as the Case–van Kampen modes, is eliminated. Furthermore, because the background distribution function in a weakly collisional system can exhibit significant deviations from a Maxwellian distribution function over long times, we show that the spectrum of Landau modes can change drastically even in the presence of slight deviations from a Maxwellian, primarily through the appearance of weakly damped modes that may be otherwise heavily damped for a Maxwellian distribution. Our results provide important insights for developing statistical theories to describe thermal fluctuations in a stellar system, which are currently a subject of great interest for N -body simulations as well as observations of gravitational systems.

79 ASTRONOMY AND ASTROPHYSICS↗

Breakdown of sound in superfluid helium

As elementary particles carry energy and momentum in the Universe, quasiparticles are the elementary carriers of energy and momentum quanta in condensed matter. And, as elementary particles, under certain conditions quasiparticles can be unstable and decay, emitting pairs of less energetic ones. Pitaevskii [Sov. Phys. JETP 9, 830 (1959)] proposed that such processes exist in superfluid helium, a quantum fluid where the very concept of quasiparticles was borne by Landau and which presented the first notable success of that concept. Pitaevskii’s decays have important consequences, including the possible breakdown of a quasiparticle [M. B. Stone et al. , Nature (London) 440, 187 (2006)]. Here, we present neutron scattering experiments, which provide evidence that such decays explain the collapsing lifetime (strong damping) of higher-energy phonon-roton sound-wave quasiparticles in superfluid helium. Importantly, this damping develops when helium is pressurized towards crystallization, or warmed towards approaching the superfluid transition. Our results resolve a number of puzzles raised by previous experiments and reveal the ubiquity of quasiparticle decays and their importance for understanding quantum matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Current-driven transverse domain wall oscillations in perpendicular spin-valve structures

Spin-transfer-driven oscillations of a transverse domain wall confined to a perpendicular spin-valve structure are investigated using a one-dimensional model. The stack consists of a polarizer, nonmagnetic spacer, soft free layer, and pinned magnetic layer. It is found that the domain-wall oscillation frequency is a nonmonotonic, highly asymmetrical function of applied electric current, showing a strong dependence on the current direction and the relative strengths of the interfacial and bulk spin-transfer torques. Micromagnetic analysis reveals that the surprising and atypical oscillator response is due to an interplay between the interfacial spin-transfer torque, the bulk spin-transfer torque, the exchange torque, and the damping torque. The underlying physical and material responses are examined, including the important role of the domain-wall twist. The competitions between the involved torques under different operating conditions suggest that the oscillator could serve as a model system to investigate magnetic and spintronic phenomena at the nanoscale. Finally, the observed current-dependent twisting of the free-layer magnetization about the axis of precession may further be found interesting for investigations of the interaction between spin-polarized current and chiral spin structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗