Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dispersive wave relations”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Runge-Kutta methods combined with compact difference schemes for the unsteady Euler equations

Recent development using compact difference schemes to solve the Navier-Stokes equations show spectral-like accuracy. A study was made of the numerical characteristics of various combinations of the Runge-Kutta (RK) methods and compact difference schemes to calculate the unsteady Euler equations. The accuracy of finite difference schemes is assessed based on the evaluations of dissipative error. The objectives are reducing the numerical damping and, at the same time, preserving numerical stability. While this approach has tremendous success solving steady flows, numerical characteristics of unsteady calculations remain largely unclear. For unsteady flows, in addition to the dissipative errors, phase velocity and harmonic content of the numerical results are of concern. As a result of the discretization procedure, the simulated unsteady flow motions actually propagate in a dispersive numerical medium. Consequently, the dispersion characteristics of the numerical schemes which relate the phase velocity and wave number may greatly impact the numerical accuracy. The aim is to assess the numerical accuracy of the simulated results. To this end, the Fourier analysis is to provide the dispersive correlations of various numerical schemes. First, a detailed investigation of the existing RK methods is carried out. A generalized form of an N-step RK method is derived. With this generalized form, the criteria are derived for the three and four-step RK methods to be third and fourth-order time accurate for the non-linear equations, e.g., flow equations. These criteria are then applied to commonly used RK methods such as Jameson's 3-step and 4-step schemes and Wray's algorithm to identify the accuracy of the methods. For the spatial discretization, compact difference schemes are presented. The schemes are formulated in the operator-type to render themselves suitable for the Fourier analyses. The performance of the numerical methods is shown by numerical examples. These examples are detailed. described. The third case is a two-dimensional simulation of a Lamb vortex in an uniform flow. This calculation provides a realistic assessment of various finite difference schemes in terms of the conservation of the vortex strength and the harmonic content after travelling a substantial distance. The numerical implementation of Giles' non-refelctive equations coupled with the characteristic equations as the boundary condition is discussed in detail. Finally, the single vortex calculation is extended to simulate vortex pairing. For the distance between two vortices less than a threshold value, numerical results show crisp resolution of the vortex merging.

Yu, Sheng-Tao↗
DOE OSTI2021

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↗
NASA NTRS1992

The mirror and ion cyclotron anisotropy instabilities

The linear dispersion equation for fully electromagnetic waves and instabilities at arbitrary directions of propagation relative to a background magnetic field B(0) in a homogeneous Vlasov plasma is solved numerically for bi-Maxwellian particle distributions. For isotropic plasmas the dispersion and damping of the three modes below the proton cyclotron frequency are studied as functions of Beta(i) and T(e)/T(i). The transport ratios of helicity, cross-helicity, Alfven ratio, compressibility, and parallel compressibility are defined. Under the condition that the proton temperature perpendicular to B(0) is greater than the parallel temperature, the growth rates and transport ratios of the mirror instability and the ion cyclotron anisotropy instability are examined and compared. Both the proton parallel compressibility and the proton Alfven ratio are significantly different for the two growing modes.

Gary, S. P.↗
DOE OSTI2021

Ion-temperature-gradient stability near the magnetic axis of quasisymmetric stellarators

Abstract The stability of the ion-temperature gradient mode in quasisymmetric stellarators is assessed. This is performed using a set of analytical estimates together with linear gyrokinetic simulations. The peak growth rates, their corresponding real frequencies and wave-vectors are identified. A comparison is made between a first-order near-axis expansion model and eleven realistic designs obtained using numerical optimization methods. It is found that while the near-axis expansion is able to replicate the growth rates, real frequencies and perpendicular wave-vector at the inner core (both using simplified dispersion relations and first-principle gyrokinetic simulations), it leads to an overestimation of the growth rate at larger radii. An approximate analytic solution of the ITG dispersion relation for the non-resonant limit suggests growth rates could be systematically higher in quasi-axisymmetric (QA) configurations compared to quasi-helically (QH) symmetric ones. However except for very close to the axis, linear gyrokinetic simulations do not show systematic differences between QA and QH configurations.

Physics↗
DOE OSTI2024

Quasilinear theory: the lost ponderomotive effects and why they matter

Quasilinear theory (QLT) has been used for modeling wave–plasma interactions for decades but remains largely heuristic. Plasma inhomogeneity, ponderomotive effects, microscopic fluctuations, and collisions are not easily accommodated from first principles in QLT, and typically are ignored entirely, due to the limitations of the standard Fourier–Laplace global-mode approach. This results in inconsistencies, for example, violation of the action conservation for nonresonant waves. However, these issues can be avoided, and the theory can be substantially generalized and corrected, if QLT is formulated using more suitable analytical tools, particularly, the Weyl symbol calculus. Here, an attempt is made to deliver an accessible review of this modern formulation, provide intuitive calculations for special cases, and elaborate on the connection with the ‘oscillation-center QLT’ originally proposed by Dewar (Phys Fluids 16:1102, 1973). A Fokker–Planck equation for a ‘dressed’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, ponderomotive forces, and interactions with background fields for a generic Hamiltonian, so many known formulations of QLT for specific plasma models become corollaries of a single unifying theory. Also, waves are allowed to be off-shell (not constrained by a dispersion relation), which allows them to accommodate microscopic fluctuations. This leads to a collision integral of the Balescu–Lenard type that has all the usual properties but is not restricted to any specific plasma model. For on-shell waves, a generalized version of the classic oscillation-center QLT is obtained. Finally, combined with the wave-kinetic equation, this formulation not only conserves particles, momentum, and energy, like the classic QLT but also reinstates the action conservation for nonresonant waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1974

Electromagnetic processes in the atmosphere of pulsars

The work consists of two parts. The first deals with the fine structure of radio pulses. Based on kinetic theory, processes occurring in the plasma shell of a pulsar when external electromagnetic radiation is present are investigated. It is shown that electromagnetic waves cause electrons to drift relative to ions, and initiate longitudinal oscillations. A dispersion equation describing the longitudinal oscillations in magnetized plasma is derived. Conditions for excitation of oscillations are found. Correlation functions of electron density are calculated, along with the coefficients of electromagnetic wave scattering. It is shown that variations in the amplitude of pulsar pulses are associated with scintillations caused by fluctuations in the plasma electron density. The second part of the study presents a mechanism for the radio emission of pulsars. The model of a rotating and a pulsating star, a neutron star with dipolar or more complex magnetic field, is examined.

Yukhimuk, A. K.↗
DOE OSTI2020

Experimental parameters, combined dynamics, and nonlinearity of a magnonic-opto-electronic oscillator (MOEO)

We report the construction and characterization of a comprehensive magnonic-opto-electronic oscillator (MOEO) system based on 1550-nm photonics and yttrium iron garnet (YIG) magnonics. The system exhibits a rich and synergistic parameter space because of the ability to control individual photonic, electronic, and magnonic components. Taking advantage of the spin wave dispersion of YIG, the frequency self-generation as well as the related nonlinear processes becomes sensitive to the external magnetic field. Besides being known as a band-pass filter and a delay element, the YIG delay line possesses spin wave modes that can be controlled to mix with the optoelectronic modes to generate higher-order harmonic beating modes. With the high sensitivity and external tunability, the MOEO system may find usefulness in sensing applications in magnetism and spintronics beyond optoelectronics and photonics.

47 OTHER INSTRUMENTATION↗
NASA NTRS1988

Ion-ion waves in the auroral region - Wave excitation and ion heating

The properties of the ion-ion mode which is excited in plasmas when two or more cold ion beams are streaming relative to one another are investigated assuming a warm electron distribution at rest (the model consistent with particle distributions observed in auroral plasma cavities). Numerical solutions are derived for the generalized electrostatic dispersion equation for parallel propagation and for oblique propagation of ion waves. It is shown that the relative ion temperature and concentration have significant effects on the stability of the accelerated plasma and the expected ion heating. Finally, a relationship between ion drift and thermal speed is derived using the marginal stability of ion-ion waves for nonzero ion temperature; the relationship was found to compare favorably with DE-1 ion observations at the high-altitude boundary of the auroral cavity.

Dusenbery, P. B.↗
DOE OSTI2021

Geometric electrostatic particle-in-cell algorithm on unstructured meshes

We present a geometric particle-in-cell (PIC) algorithm on unstructured meshes for studying electrostatic perturbations with frequency lower than electron gyrofrequency in magnetized plasmas. In this method, ions are treated as fully kinetic particles and electrons are described by the adiabatic response. The PIC method is derived from a discrete variational principle on unstructured meshes. To preserve the geometric structure of the system, the discrete variational principle requires that the electric field is interpolated using Whitney 1-forms, the charge is deposited using Whitney 0-forms and the electric field is computed by discrete exterior calculus. The algorithm has been applied to study the ion Bernstein wave (IBW) in two-dimensional magnetized plasmas. The simulated dispersion relations of the IBW in a rectangular region agree well with theoretical results. In a two-dimensional circular region with fixed boundary condition, the spectrum and eigenmode structures of the IBW are obtained from simulations. We compare the energy conservation property of the geometric PIC algorithm derived from the discrete variational principle with that of previous PIC methods on unstructured meshes. The comparison shows that the new PIC algorithm significantly improves the energy conservation property.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

0n the dynamics of warped accretion disks

We present some analyses of warping modes in inviscid near Keplerian disks taking the three-dimensional structure fully into account. The results of this investigation verify the validity of a vertical averaging approximation for thin disks when the radial wavelength is significantly longer than the disk thickness. They also indicate that long wavelength disturbances may persist for long times. Shorter wavelength disturbances in non-self gravitating disks are found to propagate with little dispersion at a speed related to the sound speed. When self-gravity becomes important, fast and slow waves are found which also propagate with little dispersion. When a small viscosity is included, the evolution of the disturbances becomes more diffusive in character.

Papaloizou, J. C. B.↗
NASA NTRS2024

Hopkinson Bar Impact Force Measurement for Application to Development of an Artificial Bird

A test method has been published for measuring forces from bird cadavers and artificial birds to demonstrate equivalency to support the use of artificial birds in aircraft certification testing. The test method involves a large diameter Hopkinson bar as the technique for measuring and comparing forces. The method involves the assumption that a force on one end of the bar results in a one-dimensional plane wave that travels down the bar. To avoid reflected waves from the end of the bar overlapping the input pulse, and due to practical limitations on the length of the bar, strain transducers are located relatively close to the impacted face of the bar, raising questions about the validity of the assumption of one-dimensional wave propagation. In addition, the large diameter of the bar can introduce dispersion effects in the propagated pulse. These must be assessed to determine the impact on the accuracy of impact force measurements and, if significant, dispersion correction techniques must be used. In this study analytical, experimental, and computational methods are used to analyze the wave propagation characteristics for impact forces expected from realistic bird impacts. While wave dispersion and effects of non-concentric impacts are present, the resulting errors are relatively small. For the response we expect to see in bird impacts up to velocities of approximately 300 m/sec the large diameter Hopkinson bar test appears to be an acceptable method for measuring the impact forces and comparing the response of real and proposed artificial birds.

Birdstrike↗
DOE OSTI2026

Connection between classical and quantum descriptions of spin waves using quantum circuits

A quantum computing circuit is presented that approximates a single spin wave quantum on a linear chain of spin 1/2 particles described by a Heisenberg Hamiltonian. The circuit is a product state where each qubit represents a spin. The spin wave motion is represented by opening the cone angle using Y rotations and then adding progressive Z rotations along the chain to represent wave propagation. We show analytically that this product state yields the correct dispersion relation in the limit of an unbounded chain. This observation is confirmed using both a simulator and various quantum processors. The use of the quantum computing paradigm in this case does not lead to a computational advantage, but rather leads to a novel conceptual connection between classical and quantum descriptions of spin waves, and may also be useful for characterizing the error in quantum processors.

magnons↗
NASA NTRS1979

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.↗
NASA NTRS2003

Advanced Computational Aeroacoustics Methods for Fan Noise Prediction

Direct computation of fan noise is presently not possible. One of the major difficulties is the geometrical complexity of the problem. In the case of fan noise, the blade geometry is critical to the loading on the blade and hence the intensity of the radiated noise. The precise geometry must be incorporated into the computation. In computational fluid dynamics (CFD), there are two general ways to handle problems with complex geometry. One way is to use unstructured grids. The other is to use body fitted overset grids. In the overset grid method, accurate data transfer is of utmost importance. For acoustic computation, it is not clear that the currently used data transfer methods are sufficiently accurate as not to contaminate the very small amplitude acoustic disturbances. In CFD, low order schemes are, invariably, used in conjunction with unstructured grids. However, low order schemes are known to be numerically dispersive and dissipative. dissipative errors are extremely undesirable for acoustic wave problems. The objective of this project is to develop a high order unstructured grid Dispersion-Relation-Preserving (DRP) scheme. would minimize numerical dispersion and dissipation errors. contains the results of the funded portion of the project. scheme on an unstructured grid has been developed. constructed in the wave number space. The characteristics of the scheme can be improved by the inclusion of additional constraints. Stability of the scheme has been investigated. Stability can be improved by adopting the upwinding strategy.

Envia, Edmane↗
DOE OSTI2022

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↗
DOE OSTI2025

Turbulence simulations of drift-Alfvén wave effects on edge-localized modes and divertor target heat flux width

A series of BOUT++ turbulence simulations are performed for two Pre-Fusion Power Operation (PFPO-1,2) phases of the ITER Research plan proclaimed in 2019. Linear simulation results show that PFPO-1 is unstable to Peeling-ballooning modes (PB), while the PFPO-2 is unstable to the coupling of PB and Drift-Alfvén wave (DAW) instabilities. The linear results are qualitatively consistent with the dispersion relation of DAW. To study the influence of the DAW on the edge-localized mode (ELM) crash, the PFPO-2 is adopted in this section. Different from the grassy ELM in [Li et al 2022 Nucl. Fusion 62, 096030], nonlinear simulations show that the ELM size of PFPO-2 is almost one-third of the grassy ELM, representing a distinct small ELM. However, simulations then show that if the PB instability is removed, the fluctuation amplitude drops by an order of magnitude and the ELM crash disappears, which is in accord with the theory in [Xu et al 2010 Physical Review Letters, Vol. 105, 175005] and the results in [Li et al 2022 Nucl. Fusion 62, 096030], confirming that the PB instability is a necessary condition for an ELM crash. Furthermore, removing the DAW drive also suppresses ELM crashes, implying that PB instability is necessary but insufficient for PFPO-2 ELM and that DAW could amplify PB-driven turbulence. In addition, by integrating results of both PFPO phases by BOUT++ turbulence analysis with heat flux width (λ q ) and electron thermal diffusivity (χ e ), these simulations establish a predictive framework for heat flux width and its dominant mechanisms, thereby providing a comprehensive scheme for the prediction of the heat flux width. Moreover, simulations indicate that DAW driving can increase the transport coefficient by enhancing the turbulent transport, leading to a broadened heat flux width once the transport coefficient exceeds

Drift-Alfvén wave↗
NASA NTRS1989

The electron beam instability revisited: Growth above and below f sub p

The growth of electrostatic waves near the plasma frequency (fp) due to an unstable electron beam is investigated by solving the unmagnetized electrostatic dispersion equation numerically. The numerical solutions are compared with analytic theories for reactive (or fluid-type) and kinetic versions of the beam instability, and for the O'Neil/Malmberg connection of the beam and Langmuir modes. Conditions for growth significantly above or below fp are given. Three general results are found: (1) The unstable waves do not grow on a mode with Langmuir dispersion except in the limit of a very dilute beam with growth on O'Neil/Malmberg's connected mode; (2) The properties of the unstable mode depend strongly on beam parameters such as beam density, speed and temperature; and (3) The frequency of maximum growth frequently lies significantly above or below fp, and differs significantly from that predicted by the Langmuir dispersion relation. Results imply important consequences for theories of strong turbulence and nonlinear wave-wave processes, and observational identification of fp from observed wave frequencies.

Cairns, I. H.↗