Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dispersion relation”

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 109 records · Page 6

Description of global EGAM in the maximum of local frequency during current ramp-up discharges in DIII-D

Energetic-particle-induced geodesic acoustic modes, EGAMs (Fu, Phys. Rev. Let., vol. 101, 2008, pp. 185002), driven by neutral beam injection (NBI), have been observed in many DIII-D tokamak experiments (Nazikian et al., Phys. Rev. Lett., vol. 101, 2008, pp. 185001). This mechanism has been theoretically investigated in (Qiu et al., Plasma Phys. Control. Fusion, vol. 52, 2010, pp. 095003), using a sharp energetic particle distribution function, and in (Qu et al., Plasma Phys. Control. Fusion, vol. 59, 2017, pp. 055018), where the dispersion relation and eigenmode behaviour were obtained for the situation of early beam scenario, that is, for times smaller than the beam slowing down time. In this work, we extend these studies determining the eigenmode for beyond the slowing down time, in a scenario with reverse safety factor q profile, where a small concentration of energetic ions can produce an off-axis maximum in the GAM dispersion relation. Here, the characteristics of EGAM are analytically studied with the drift kinetic equation together with the MHD code NOVA. The toroidal energetic ion transit frequency, coupled with the GAM frequency, produces the maximum in the dispersion relation where the eigenmode can be found. The quantitative correspondence of experimental results with the predictions of the proposed model is analysed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Finite boundary effects on the spherical Rayleigh–Taylor instability between viscous fluids

For the Rayleigh–Taylor unstable arrangement of a viscous fluid sphere embedded in a finite viscous fluid spherical shell with a rigid boundary and a radially directed acceleration, a dispersion relation is developed from a linear stability analysis using the method of normal modes. [Formula: see text] is the radially directed acceleration at the interface. ρ i denotes the density, μ i is the viscosity, and R i is the radius, where i = 1 is the inner sphere and i = 2 is the outer sphere. The dispersion relation is a function of the following dimensionless variables: viscosity ratio [Formula: see text], density ratio [Formula: see text], spherical harmonic mode n, [Formula: see text], [Formula: see text], and the dimensionless growth rate [Formula: see text], where σ is the exponential growth rate. We show that the boundedness provided by the outer spherical shell has a strong influence on the instability behavior, which is reflected not only in the modulation of the growth rate but also in the selection of the most unstable modes that are physically possible. This outer boundary effect is quantified by the relative magnitude of the radius ratio H. We find that when H is close to unity, lower order harmonics are excluded from becoming the most unstable within a vast region of the parameter space. In other words, the effect of H has precedence over the other controlling parameters d, B, and a wide range of s in establishing what the lowest most unstable mode can be. When H ~ 1, low order harmonics can become the most unstable only for s >> 1. However, in the limit when s → ∞, we show that the most unstable mode is n = 1 and derive the dispersion relation in this limit. The exclusion of most unstable low order harmonics caused by a finite outer boundary is not realized when the outer boundary extends beyond a certain threshold length-scale in which case all modes are equally possible depending on the value of B.

74 ATOMIC AND MOLECULAR PHYSICS↗

Improved analytic modeling of the linear Rayleigh–Taylor instability with plasma transport

The Rayleigh–Taylor (R–T) instability with transport effects has long been studied analytically in the neutral fluid regime. How this picture may differ in plasmas has been given less consideration, however. Leveraging previous numerical simulations with plasma viscosity and diffusion of R–T at a binary plasma interface, we show how plasma transport and kinetic effects alter the linear R–T dispersion relation. Our insight informs a much improved analytical dispersion relation, which better matches the simulation data than previously established models. Additionally, we demonstrate the importance of employing asymptotically correct viscosity coefficients in fluid codes for accurately capturing plasma R–T growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dispersion and the speed-limited particle-in-cell algorithm

This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. Here, these results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Contact solitons as spin pistons in one-dimensional ferromagnetic channels

Ferromagnetic channels subject to spin injection have been theoretically shown to sustain dissipative exchange flows (DEFs). In the strong injection regime, a soliton is stabilized at the injection site, which has been termed a contact soliton DEF or CS-DEF. Here, in this work, we investigate the modulation of CS-DEFs as a mechanism to inject magnons into a DEF. By varying the injection current, the parameters connecting the soliton with the DEF via a boundary layer are varied. These lead to a modification in the DEF which is interpreted as a spin piston. The injected magnons follow the expected dispersion relation of DEFs, akin to the Bogoliubov–de Gennes dispersion relation. As a result, this work demonstrates that changes in the injected current can pump magnons along a ferromagnetic channel via dissipative exchange flows.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cold-hot coupled waves in a flowing magnetized plasma

Abstract Nonlinear coupling of cold and hot waves in a flowing magnetized plasma is analyzed with the Vlasov equation. An analytical solution is obtained for cold waves of a small amplitude (weak flow) and a long wavelength. The distribution function is obtained by integrating the kinetic equation along a perturbed phase-space trajectory for a time-varying plasma flow. The kinetic description presents a generalized dispersion relation that involves resonances depending on cold and hot wave dispersions. Coherent fluid motion leads to radiation peaks in addition to the cyclotron harmonics, where the wavenumber of the cold wave determines the peak frequencies. The peaks appear narrow when the wave propagates perpendicular to the time-averaged flow while they become broad due to the Doppler effect when the wave propagates parallel to the flow. Fully kinetic particle-in-cell simulations corroborate the theoretical predictions. The dispersion relation and resulting wave spectra provide information about plasma parameters and flow properties.

Physics↗

Doubly diffusive magnetic buoyancy instability in the solar interior

An investigation of the buoyancy of diffuse magnetic fields has shown that in the presence of rotation, static equilibrium configurations of the toroidal magnetic field and ambient plasma can exist. In that case, the escape of toroidal magnetic flux from the solar interior may be determined by the growth of instabilities which the equilibrium configuration may be subject to. In connection with the present investigation, it is assumed that in the region of toroidal magnetic flux amplification, the magnetic field has not as yet filamented into flux ropes, and is therefore 'diffuse'. A study is conducted of the MHD stability of an electrically conducting and differentially rotating gas in the presence of a toroidal magnetic field, an external constant gravitational field, and radiance pressure. The full dispersion relation for the magnetic buoyancy problem is developed, and the solutions of the dispersion relation are discussed.

Schmitt, J. H. M. M.↗

Causality constraints on corrections to Einstein gravity

We study constraints from causality and unitarity on 2 → 2 graviton scattering in four-dimensional weakly-coupled effective field theories. Together, causality and unitarity imply dispersion relations that connect low-energy observables to high-energy data. Using such dispersion relations, we derive two-sided bounds on gravitational Wilson coefficients in terms of the mass M of new higher-spin states. Our bounds imply that gravitational interactions must shut off uniformly in the limit G → 0, and prove the scaling with M expected from dimensional analysis (up to an infrared logarithm). We speculate that causality, together with the non-observation of gravitationally-coupled higher spin states at colliders, severely restricts modifications to Einstein gravity that could be probed by experiments in the near future.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Magnetorotational instability in a swirling partially ionized gas

ABSTRACT The magnetorotational instability (MRI) has been proposed as the method of angular momentum transport that enables accretion in astrophysical discs. However, for weakly ionized discs, such as protoplanetary discs, it remains unclear whether the combined non-ideal magnetohydrodynamic (MHD) effects of Ohmic resistivity, ambipolar diffusion, and the Hall effect make these discs MRI stable. While much effort has been made to simulate non-ideal MHD MRI, these simulations make simplifying assumptions and are not always in agreement with each other. Furthermore, it is difficult to directly observe the MRI astrophysically because it occurs on small scales. Here, we propose the concept of a swirling gas experiment of weakly ionized argon gas between two concentric cylinders threaded with an axial magnetic field that can be used to study non-ideal MHD MRI. For our proposed experiment, we derive the hydrodynamic equilibrium flow and a dispersion relation for MRI that includes the three non-ideal effects. We solve this dispersion relation numerically for the parameters of our proposed experiment. We find it should be possible to produce a non-ideal MRI in such an experiment because of the Hall effect, which increases the MRI growth rate when the vertical magnetic field is anti-aligned with the rotation axis. As a proof of concept, we also present experimental results for a hydrodynamic flow in an unmagnetized prototype. We find that our prototype has a small, but non-negligible, α-parameter that could serve as a baseline for comparison to our proposed magnetized experiment, which could be subject to additional turbulence from the MRI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-harmonic Rutherford island theory

Rutherford island theory, which governs the nonlinear evolution of tearing modes in tokamak plasmas, is generalized to take into account situations in which the conventional one-harmonic approximation is not valid. The analysis incorporates non-inductive currents driven by radio frequency (RF) electromagnetic waves injected into the plasma. A multi-harmonic tearing mode dispersion relation is derived that takes the form of a nonlinear inhomogeneous matrix eigenvalue problem. The dispersion relation is solved in the so-called two-harmonic approximation, in which only the principal Fourier harmonic of the perturbed magnetic flux and its first overtone are included in the calculation. In the absence of RF current drive, the nonlinear behavior of a tearing mode predicted in the two-harmonic approximation does not differ substantially from that predicted in the one-harmonic approximation. On the other hand, RF current drive that is sufficiently localized in the vicinity of the O-points of the mode's magnetic island chain is capable of triggering bifurcations of the O-points (which is impossible in the one-harmonic approximation). However, the current drive is incapable of triggering bifurcations of the island X-points. This finding is significant because Bardóczi and Evans recently observed bifurcations of magnetic island chain O-points in the presence of RF current drive in the DIII-D tokamak but did not observe bifurcations of the X-points. Finally, the changes in the topology of the magnetic island flux-surfaces induced by RF current drive are found to facilitate the stabilization of the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Wave modes of the Io plasma torus

The corotating magnetospheric convection model is analyzed with the objective of deriving the fundamental wave modes present in the Io torus. The time-dependent equations of motion are linearized for the case of the absence of plasma mass and momentum sources. By separating the magnetosphere-ionosphere coupling equation into two equations and by adopting a specific analytic form for the equilibrium plasma density, a cubic equation is derived for the dispersion relation defining the angular frequencies of the fundamental wave modes. A detailed analysis of the dispersion relation yields three fundamental wave modes. One of the modes is identified as the previously studied inertial interchange (unstable) mode, while the other two are newly derived decaying modes. The results obtained are discussed in the context of plasma transport in the Io torus.

Summers, D.↗

The cosmic Doppler instability

The equations governing the behavior of perturbations of a mixture of nearly homogeneous and isotropic matter and radiation are derived, using a diffusion approximation where spatial perturbations in the radiation spectrum are allowed to vary with frequency. A simple model of line opacity leads to dispersion relations which display a new bulk instability. The model is used to derive an approximate dispersion relation for radiation interacting via resonance scattering opacity in atomic hydrogen at low density and low temperature. Possible applications to cosmology are briefly discussed.

Hogan, Craig J.↗

A revised parametrization of the dispersion constants for high dispersion IUE images

An economical system of parameterizing dispersion constants spectra is derived and applied to International Ultraviolet Explorer (IUE) spectra by equating specific terms to geometric coordinates of an echellogram. Because of the virtual property of the IUE camera pixels, several geometrical rectifications can be performed analytically before the dispersion relations are solved for line and sample position. This permits the relations to be decoupled and reduced to two terms apiece, the simplest representation possible. A generalization of these principles allows the dispersion relations to be determined neatly in other echelle systems having certain well understood distortions. Prismatic cross dispersion, S wave distortion, and optical aberrations are some areas in which this system can be used.

Smith, Myron A.↗

The radiation pressure-driven Rayleigh-Taylor instability - Analysis and application to QSO emission line clouds

The growth of perturbations in a photoionized gas slab accelerated by radiation pressure under conditions like those in QSO emission-line regions is analyzed. A linear dispersion relation is derived, and numerical radiation-transfer calculations are performed to evaluate the coefficients and roots of that dispersion relation for L-alpha optical depths between 0.1 and 10,000. The nonlinear growth of the waves is estimated, and it is concluded that complete dispersal of the clouds is likely. On the basis of qualitative arguments it is suggested that clouds of greater optical depth are linearly stable.

Krolik, J. H.↗

Fourier Analyses of High-Order Continuous and Discontinuous Galerkin Methods

In this paper, we present a Fourier analysis of wave propagation problems subject to a class of continuous and discontinuous discretizations using high-degree Lagrange polynomials. This allows us to obtain explicit analytical formulas for the dispersion relation and group velocity and, for the first time to our knowledge, characterize analytically the emergence of gaps in the dispersion relation at specific wavenumbers, when they exist, and compute their specific locations. Wave packets with energy at these wavenumbers will fail to propagate correctly, leading to significant numerical dispersion. We also show that the Fourier analysis generates mathematical artifacts, and we explain how to remove them through a branch selection procedure conducted by analysis of eigenvectors and associated reconstructed solutions. The higher frequency eigenmodes, named erratic in this study, are also investigated analytically and numerically.

97 MATHEMATICS AND COMPUTING↗

Beat, modulational, and decay instabilities of a circularly polarized Alfven wave

A circulary polarized low-frequency electomagnetic pump wave propoagating along an ambient magnetic field is known to be unstable to the growth of several parallel-propagating parametric instabilties. If ion-cyclotron effects are retained in a two-fluid description, the dispersion relation is a sixth-order polynomial. We present a series of new analytical approximations to this dispersion relation. We emphasize new results for the beat instability that occurs as an interaction of the forward prpagating upper sideband with the backward propagating lower sideband. The nature of the beat instabitlity depends on beta = (v(sub sound)/v(sub A)(exp 2) and on the sense of polarization of the pump wave. The beat and decay instabilities can occur together if the pump is left-handed (i.e., ion resonant) and if beta is less than or approximately 1, but they cannot occur together if the pump is right-handed. For a left-handed pump the beat mode is the only instability if beta is greater than or approximately 1. If the pump is right-handed and beta is greater than or approximately 1, then the beat instability exists only when the pump amplitude exceeds a threshold value, and the beat will be the only instability if the pump amplitude is large enough to stabilize the modulational instability. If the pump is left-handed and beta is less than or approximately 1, then the beat mode is stabilized when the pump amplitude becomes sufficiently large. The beat instability primarily produces a forward propagating transverse wave in the upper sideband. Thus if beta is greater than or approximately 1, the instabilities considered here do not produce the backward propagating waves which are thought to affect turbulence and the evolution of cross helicity in the solar wind. New analytical results are presented also for the decay and modulational instabilites when beta is approximately equal to 1.

Hollweg, Joseph V.↗

Three-dimensional coupling of electron cyclotron drift instability and ion–ion two stream instability

Electron cyclotron drift instability (ECDI) and ion–ion two stream instability (IITSI) are both kinetic instabilities that can be present in low-temperature, partially magnetized plasmas. The coupling of instabilities in a three-dimensional configuration leads to the existence of more than one unstable roots to the kinetic dispersion relation. In this paper, a generalized method has been developed for numerically evaluating solutions to the three-dimensional dispersion relation for coupled ECDI and IITSI, assuming cold singly and doubly charged ions and a Maxwellian velocity distribution function for the electrons. The present study demonstrates the coupling between ECDI and IITSI that affects the most unstable mode as a function of the wavenumbers in three dimensions and various plasma properties, including the applied electric field, magnetic field, electron temperature, ion velocities, and plasma density. One of the most notable results is that, while the most unstable mode with the largest growth rate is in the direction of the E×B drift in the two-dimensional cases, the most unstable mode for the three-dimensional configuration occurs in the oblique direction between the applied electric field and the E×B drift. Here, this agrees with experimental observations in cross field plasma sources using coherent Thomson scattering.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ULF waves in the Io torus: Ulysses observations

Throughout the Io torus, Ulysses has observed intense ultralow frequency (ULF) wave activity in both electric and magnetic components. Such ULF waves have been previously suggested as the source of ion precipitation leading to Jovian aurorae. The peaks of the wave spectra are closely related to the ion cyclotron frequencies, which is evidence of the waves being ion cyclotron waves (ICWs). Analysis of the dispersion relation using a multicomponent density model shows that at high latitudes (approximately 30 deg), peak frequencies of the waves fall into L mode branches of guided or unguided ICWs. Near the equator, in addition to the ICWs below f(sub cO(2+)), there are strong signals at approximately 10 Hz which require an unexpectedly large energetic ion temperature anistropy to be explained by the excitation of either convective or nonconvective ion cyclotron instabilities. Their generation mechanism remains open for the future study. Evaluation of the Poynting vector and the dispersion relation analysis suggest that the waves near the equator had a small wave angle relative to the magnetic field, while those observed at high latitudes were more oblique. The polarization of the waves below f(sub cH(+)) is more random than that of the whistler mode waves, but left-hand-polarized components of the waves can still be seen. The intensity of the ICWs both near the equator and at high latitudes are strong enough to meet the requirement for producing strong pitch angle scattering of energetic ions.

Lin, Naiguo↗