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 145 records · Page 8

The possible role of high-frequency waves in heating solar coronal loops

We investigate the role of high-frequency waves in the heating of solar active region coronal loops. We assume a uniform background magnetic field, and we introduce a density stratification in a direction perpendicular to this field. We focus on ion compressive viscosity as the damping mechanism of the waves. We incorporate viscosity self-consistently into the equations, and we derive a dispersion relation by adopting a slab model, where the density inside the slab is greater than that outside. Such a configuration supports two types of modes: surface waves and trapped body waves. In order to determine under what conditions these waves may contribute to the heating of active regions, we solve our dispersion relation for a range of densities, temperatures, magnetic field strengths, density ratios, wavevector magnitudes, wavevector ratios, and slab widths. We find that surface waves exhibit very small damping, but body waves can potentially damp at rates needed to balance radiative losses. However, the required frequencies of these body waves are very high. For example, the wave frequency must be at least 5.0/s for a slab density of 10(exp 9,5)/cc, a slab temperature of 10(exp 6,5) K, a field strength of 100 G, and a density ratio of 5. For a slab density of 10(exp 10)/cc, this frequency increases to 8.8/s. Although these frequencies are very high, there in no observational evidence to rule out their existence, and they may be generated both below the corona and at magnetic reconnection sites in the corona. However, we do find that, for slab densities of 10(exp 10)/cc or less, the dissipation of high-frequency waves will be insufficient to balance the radiative losses if the magnetic field strength exceeds roughly 200 G. Because the magnetic field is known to exceed 200 G in many active region loops, particularly low-lying loops and loops emanating from sunspots, it is unlikely that high-frequency waves can provide sufficient heating in these regions.

Porter, Lisa J.↗

Rigid-Mode Limit of the Yokoya Matrix Formalism and the Burov-Lebedev Dispersion Equation

Transverse single-bunch instabilities of space-charge-dominated coasting beams with round and flat transverse geometries are studied using a unified dispersion-relation framework. The analysis combines the Burov-Lebedev formalism, which captures space-charge tune spread, Landau damping, and instability threshold behavior, with Yokoya’s projection method for representing coherent transverse mode structure and its dependence on beam aspect ratio. In the rigid-beam limit, the formulation reduces to a scalar dispersion relation of Burov-Lebedev paper. For non-rigid transverse oscillations, truncation of Yokoya’s Hermite-based expansion yields a finite-dimensional matrix eigenvalue problem in which space-charge and coupling impedance effects enter through Burov-Lebedev–type denominators. This approach provides a consistent basis for comparing rigid and non-rigid instability behavior in round and flat beams and for assessing the role of beam ellipticity in modifying coherent mode structure and stability thresholds.

43 PARTICLE ACCELERATORS↗

Theoretical and global simulation analysis of collisional microtearing modes

Microtearing modes (MTMs) are suggested as a candidate for anomalous thermal transport in tokamak H-mode discharges. This study investigates MTMs in tokamak plasmas, employing simulations in the BOUT++ framework. It simplifies and linearizes the governing equations in detailed linear simulations. The study meticulously evaluates various conductivity models under diverse plasma conditions and collision regimes. The research thoroughly assesses different conductivity models across a range of plasma conditions and collision regimes. A unified dispersion relation that includes both MTM and Drift-Alfvén Wave (DAW) instabilities is derived, showing that DAW and MTM instabilities occur at varying distances from the rational surface. Specifically, MTMs become unstable near the rational surface but stabilize farther away, while drift-Alfvén instability appears farther from the rational surface. The study also re-derives MTM dispersion relations using Ohm's law and the vorticity equation, providing a thorough analysis of electromagnetic and electrostatic interactions in tokamaks. Global simulations demonstrate an inverse correlation between MTM growth rates and collisionality, and a direct correlation with temperature gradients. The nonalignment of the rational surface with the peak ω*e stabilizes the MTMs. Nonlinear simulations highlight electron temperature relaxation as the primary saturation mechanism for MTMs, with magnetic flutter identified as the dominant mode of electron thermal transport.

Fan, K. (ORCID:0000000227518809)↗

Wave evolution in the marginal ice zone - Model predictions and comparisons with on-site and remote data

The ocean-wave dispersion relation and viscous attenuation by a sea ice cover were studied for waves in the marginal ice zone (MIZ). The Labrador ice margin experiment (Limex), conducted off the east coast of Newfoundland, Canada in March 1987, provided aircraft SAR, wave buoy, and ice property data. Based on the wave number spectrum from SAR data, the concurrent wave frequency spectrum from ocean buoy data, and accelerometer data on the ice during Limex '87, the dispersion relation has been derived and compared with the model. Accelerometers were deployed at the ice edge and into the ice pack. Data from the accelerometers were used to estimate wave energy attenuation rates and compared with the model. The model-data comparisons are reasonably good for the ice conditions observed during Limex' 87.

Liu, A. K.↗

Current driven instabilities of an electromagnetically accelerated plasma

A plasma instability that strongly influences the efficiency and lifetime of electromagnetic plasma accelerators was quantitatively measured. Experimental measurements of dispersion relations (wave phase velocities), spatial growth rates, and stability boundaries are reported. The measured critical wave parameters are in excellent agreement with theoretical instability boundary predictions. The instability is current driven and affects a wide spectrum of longitudinal (electrostatic) oscillations. Current driven instabilities, which are intrinsic to the high-current-carrying magnetized plasma of the magnetoplasmadynmic (MPD) accelerator, were investigated with a kinetic theoretical model based on first principles. Analytical limits of the appropriate dispersion relation yield unstable ion acoustic waves for T(i)/T(e) much less than 1 and electron acoustic waves for T(i)/T(e) much greater than 1. The resulting set of nonlinear equations for the case of T(i)/T(e) = 1, of most interest to the MPD thruster Plasma Wave Experiment, was numerically solved to yield a multiparameter set of stability boundaries. Under certain conditions, marginally stable waves traveling almost perpendicular to the magnetic field would travel at a velocity equal to that of the electron current. Such waves were termed current waves. Unstable current waves near the upper stability boundary were observed experimentally and are in accordance with theoretical predictions. This provides unambiguous proof of the existence of such instabilites in electromagnetic plasma accelerators.

Chouetri, E. Y.↗

The Lifshitz Regime and its Experimental Signals

We discuss the possibility of a Lifshitz regime, where the dispersion relation for Goldstone bosons and related fields has a minimum at nonzero momenta. Studies with the Functional Renormalization Group suggest that this occurs over a wide region in the plane of temperature and baryon chemical potential. Conversely, the FRG finds that the region in which fluctuations from a critical endpoint are significant is rather small. Here, we suggest that this is due generically to the narrowness of the tricritical region in the chiral limit. Even if particles are produced in thermal equilibrium, a dispersion relation which is non-monotonic in momenta produces what appears to be non-thermal behavior.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Polarization and properties of low-frequency waves in warm magnetized two-fluid plasma

This paper presents the derivation of a general wave dispersion relation for warm magnetized plasma under the two-fluid formalism. The discussion is quite general except for the assumption of low frequency and slow phase speed, for which the displacement current is negligible, under the implicit assumption that the plasma is sufficiently dense to satisfy the condition ωpe>ωce, where ωpe and ωce denote the plasma oscillation frequency and electron gyro frequency, respectively. The present discussion does not invoke charge neutrality associated with the fluctuations although it is implicitly satisfied. The resulting dispersion relation that includes the fluid thermal effects shows that there are three eigen modes, which include those corresponding to ideal MHD, namely, fast, slow, and kinetic Alfvén waves, as well as higher-frequency modes including the ion and electron cyclotron and lower-hybrid resonances. The fluid effects in the ideal MHD wave branches are influenced by the finite Larmor radius scales, and when the wave number in the cross field direction is comparable to these values, the fluid effects become significant. It is found that the Larmor radius should be interpreted in the sense as ion-acoustic gyro-radius instead of ion thermal gyro radius only. That is, it is found that the electrons also contribute to the non-ideal effect associated with the kinetic Alfvén wave. A comprehensive explanation of the polarization of each mode is also presented. The present findings indicate that the polarity may change its sign only for the kinetic Alfvén mode branch and that such a transition is based on the propagation angle. When such a change does take place, it is found that the kinetic Alfvén wave transits to an ion-acoustic mode. For each branch, it is also found that the electric field along the ambient magnetic field is purely transverse.

Physics↗

Magnetized laser–plasma interactions in high-energy-density systems: Parallel propagation

In this study, we investigate parametric processes in magnetized plasmas, driven by a large-amplitude pump light wave. Our focus is on laser–plasma interactions relevant to high-energy-density (HED) systems, such as the National Ignition Facility and the Sandia MagLIF concept. We present a self-contained derivation of a “parametric” dispersion relation for magnetized three-wave interactions, meaning the pump wave is included in the equilibrium, similar to the unmagnetized work of Drake et al., Phys. Fluids 17, 778 (1974). For this, we use a multi-species plasma fluid model and Maxwell's equations. The application of an external B field causes right- and left-polarized light waves to propagate with differing phase velocities. This leads to Faraday rotation of the polarization, which can be significant in HED conditions. Phase-matching and linear wave dispersion relations show that Raman and Brillouin scattering have modified spectra due to the background B field, though this effect is usually small in systems of current practical interest. We study a scattering process we call stimulated whistler scattering, where a light wave decays to an electromagnetic whistler wave (ω≲ω ce ) and a Langmuir wave. This only occurs in the presence of an external B field, which is required for the whistler wave to exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear waves in a thin dielectric film on the surface of a topological insulator

Based on the dispersion relation for a wave localised in a thin film of a nonlinear dielectric, which is located on the surface of a topological insulator, we have derived a system of equations that describes the propagation of a surface wave. It is shown that the longitudinal and transverse tangential components of the electric field vector are related due to the nonlinearity of the film and change periodically during propagation. It is found that the rotation period of this vector is determined by the axion charge of the topological dielectric and the nonlinear susceptibility of the thin film. (paper)

36 MATERIALS SCIENCE↗

Comparison of Non-Parabolic Hydrodynamic Simulations for Semiconductor Devices

Parabolic drift-diffusion simulators are common engineering level design tools for semiconductor devices. Hydrodynamic simulators, based on the parabolic band approximation, are becoming more prevalent as device dimensions shrink and energy transport effects begin to dominate device characteristic. However, band structure effects present in state-of-the-art devices necessitate relaxing the parabolic band approximation. This paper presents simulations of ballistic diodes, a benchmark device, of Si and GaAs using two different non-parabolic hydrodynamic formulations. The first formulation uses the Kane dispersion relationship in the derivation of the conservation equations. The second model uses a power law dispersion relation {(hk)(exp 2)/2m = xW(exp Y)}. Current-voltage relations show that for the ballistic diodes considered. the non-parabolic formulations predict less current than the parabolic case. Explanations of this will be provided by examination of velocity and energy profiles. At low bias, the simulations based on the Kane formulation predict greater current flow than the power law formulation. As the bias is increased this trend changes and the power law predicts greater current than the Kane formulation. It will be shown that the non-parabolicity and energy range of the hydrodynamic model based on the Kane dispersion relation are limited due to the binomial approximation which was utilized in the derivation.

Smith, A. W.↗

Azimuthal spoke propagation in Hall Effect Thrusters

Spokes are azimuthally propagating perturbations in the plasma discharge of Hall Effect Thrusters (HETs) that travel in the E × B direction and have been observed in many different systems. The propagation of azimuthal spokes are investigated in a 6 kW HET known as the H6 using ultra-fast imaging and azimuthally spaced probes. A spoke surface is a 2-D plot of azimuthal light intensity evolution over time calculated from 87,500 frames/s videos. The spoke velocity has been determined using three methods with similar results: manual fitting of diagonal lines on the spoke surface, linear cross-correlation between azimuthal locations and an approximated dispersion relation. The spoke velocity for three discharge voltages (300, 400 and 450 V) and three anode mass flow rates (14.7, 19.5 and 25.2 mg/s) yielded spoke velocities between 1500 and 2200 m/s across a range of normalized magnetic field settings. The spoke velocity was inversely dependent on magnetic field strength for low B-field settings and asymptoted at B-field higher values. The velocities and frequencies are compared to standard drifts and plasma waves such as E × B drift, electrostatic ion cyclotron, magnetosonic and various drift waves. The empirically approximated dispersion relation yielded a characteristic velocity that matched the ion acoustic speed for 5 eV electrons that exist in the near-anode and near-field plume regions of the discharge channel based on internal measurements. Thruster performance has been linked to operating mode where thrust-to-power is maximized when azimuthal spokes are present so investigating the underlying mechanism of spokes will benefit thruster operation.

Polk, James E.↗

Azimuthal Spoke Propagation in Hall Effect Thrusters

Spokes are azimuthally propagating perturbations in the plasma discharge of Hall Effect Thrusters (HETs) that travel in the E x B direction and have been observed in many different systems. The propagation of azimuthal spokes are investigated in a 6 kW HET known as the H6 using ultra-fast imaging and azimuthally spaced probes. A spoke surface is a 2-D plot of azimuthal light intensity evolution over time calculated from 87,500 frames/s videos. The spoke velocity has been determined using three methods with similar results: manual fitting of diagonal lines on the spoke surface, linear cross-correlation between azimuthal locations and an approximated dispersion relation. The spoke velocity for three discharge voltages (300, 400 and 450 V) and three anode mass flow rates (14.7, 19.5 and 25.2 mg/s) yielded spoke velocities between 1500 and 2200 m/s across a range of normalized magnetic field settings. The spoke velocity was inversely dependent on magnetic field strength for low B-field settings and asymptoted at B-field higher values. The velocities and frequencies are compared to standard drifts and plasma waves such as E x B drift, electrostatic ion cyclotron, magnetosonic and various drift waves. The empirically approximated dispersion relation yielded a characteristic velocity that matched the ion acoustic speed for 5 eV electrons that exist in the near-anode and near-field plume regions of the discharge channel based on internal measurements. Thruster performance has been linked to operating mode where thrust-to-power is maximized when azimuthal spokes are present so investigating the underlying mechanism of spokes will benefit thruster operation.

Cylindrical Hall Thruster (CHT)↗

Modification of ion-temperature-gradient turbulence by impurities in stellarator plasmas

Recent nonlinear gyrokinetic simulations have shown that impurities can strongly modify the turbulent heat flux in stellarator plasmas. Here, the ion-temperature-gradient (ITG) dispersion relation in a plasma containing impurities is analytically solved in certain limits and an expression for the modification of the ITG growth rate by impurities is derived. The analytical expression is the sum of three terms corresponding to three different physical causes (impurity density gradient, impurity temperature gradient and dilution) of the change in the growth rate. The scalings predicted analytically for the modification of the growth rate are shown to be reproduced by linear gyrokinetic simulations. The conditions for reduction or increase of the ITG growth by impurities are also correctly predicted by the analytical solution to the dispersion relation. Finally, a remarkable correlation is found between the analytical expression for the modification of the growth rate and the modification of the turbulent heat flux obtained from nonlinear gyrokinetic simulations.

Calvo, Iván [Research Centre for Energy, Environme↗

Effects of finite plasma beta on the lower-hybrid-drift instability

The local dispersion relation for the lower-hybrid-drift (LHD) instability is derived and analyzed, taking into account the finite-beta effects associated with transverse electromagnetic perturbations as well as with resonant and nonresonant electron-orbit modifications due to magnetic-field gradients. The influence of finite-beta effects on the LHD instability is calculated in a fully self-consistent manner for arbitrary values of electron-ion temperature ratio, local beta, cross-field ExB velocity/ion thermal speed ratio, and other plasma parameters. Stability properties are investigated analytically for the case of cold electrons, and the local dispersion relation is solved numerically in the parameter regime of most interest for high-density plasma pinches. The results show that for all parameter regimes studied, the net effect of finite plasma beta is to reduce the maximum growth rate of the LHD instability, although the details can vary, depending on the plasma parameters. Except in the limit where the electron/ion temperature ratio tends to zero, it is found that there is a critical value of plasma beta above which the LHD instability is completely stabilized.

Davidson, R. C.↗

Waves in frazil and pancake ice and their detection in Seasat synthetic aperture radar imagery

Measurements of wavelength and direction derived from SAR imagery obtained from Seasat in October 1987 are compared with a theoretical model of waves propagating into an ice cover composed of frazil and pancake ice. The dispersion relation for phase and group velocities is derived, and the existence of an upper frequency limit for propagation into the ice is found. The wave radiation pressure exerted on the ice is derived from the reflection coefficient at the ice edge; it is found to cause a slick of frazil ice backed by thicker floes to become more dense or thick with increasing penetration. When the frazil ice becomes transformed into pancake ice, the raised edges of the pancakes cause the ice to appear bright despite the fact that there are no waves present at the Bragg wavelength. Directional wave number spectra are derived outside and inside the ice cover. The change of wavelength and angle of refraction of the dominant wave entering the ice field are both found to be characteristic of the dispersion relation derived theoretically.

Wadhams, Peter↗

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗

Observation of intravalley phonon scattering of 2s excitons in MoSe 2 and WSe 2 monolayers

Abstract We present a high-resolution resonance Raman study of hBN encapsulated MoSe 2 and WSe 2 monolayers at 4 K using excitation energies from 1.6 eV to 2.25 eV. We report resonances with the WSe 2 A2s and MoSe 2 A2s and B2s excited Rydberg states despite their low oscillator strength. When resonant with the 2s states we identify new Raman peaks which are associated with intravalley scattering between different Rydberg states via optical phonons. By calibrating the Raman scattering efficiency and separately constraining the electric dipole matrix elements, we reveal that the scattering rates for k = 0 optical phonons are comparable for both 1s and 2s states despite differences in the envelope functions. We also observe multiple new dispersive Raman peaks including a peak at the WSe 2 A2s resonance that demonstrates non-linear dispersion and peak-splitting behavior that suggests the dispersion relations for dark excitonic states at energies near the 2s state are extremely complex.

Materials Science↗

Self-consistent theory of the Darrieus–Landau and Rayleigh–Taylor instabilities with self-generated magnetic fields

The Rayleigh–Taylor (RT) and Darrieus–Landau (DL) instabilities are studied in an inertial confinement fusion context within the framework of small critical-to-shell density ratio D R and weak acceleration regime, i.e., large Froude number Fr. The quasi-isobaric analysis in Sanz et al. is completed with the inclusion of non-isobaric and self-generated magnetic-field effects. The analysis is restricted to perturbation wavelengths k~1 larger than the conduction length scale at the ablation front, yet its validity ranges from wavelengths shorter and larger than the conduction layer width (distance between the ablation front and the critical surface). The use of a sharp boundary model leads to a single analytical expression of the dispersion relation encompassing both instabilities. The two new effects come into play by modifying the perturbed mass and momentum fluxes at the ablation front. The momentum flux (perturbed pressure at the spike) is the predominant stabilizing mechanism in the RT instability (overpressure) and the driving mechanism in the DL instability (underpressure). The non-isobaric effects notably modify the scaling laws in the DL limit, leading to an underpressure scaling as ~k –11/15 rather than ~k –2/5 obtained in the quasi-isobaric model. The magnetic fields are generated due to misalignment between pressure and density gradients (Biermann battery effect). They affect the hydrodynamics by bending the heat flux lines. Within the framework of this paper, they enhance ablation, resulting in a stabilizing effect that peaks for perturbation wavelengths comparable to the conduction layer width. Here, the combination of parameters D R Fr 2/3 defines the region of predominance of each instability in the dispersion relation. It is proven that the DL region falls outside of the parameter range in inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗