Wave propagation in a partly ionized gas. II.
Wave propagation in partially ionized gas, deriving dispersion relations satisfying hydrodynamic equations
Meskan, D. A.↗
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.
Electron-cyclotron maser instability in relativistic plasmas
The electron-cyclotron maser instability is studied for the case of an anisotropic electron velocity distribution in the regime where the relativistic corrections to the wave dispersion are significant. Solution of the linear dispersion relation reveals that when the plasma frequency-gyrofrequency ratio is less than v(te)/c, the instability is localized just below k(perpendicular)c/Omega(e) = 1. The growth rate is then strongly peaked for emission at 90 deg to the magnetic field and is considerably larger than would be the case if the cold-plasma dispersion theory were valid. These features are confirmed by EM particle simulations.
Pritchett, P. L.↗
Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer
The objective of this work is to compare a high-order solver with a low-order solver for performing Large-Eddy Simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the highorder method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.
large eddy simulations↗
Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer
The objective of this work is to compare a high-order solver with a low-order solver for performing large-eddy simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the high-order method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.
large eddy simulations↗
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.↗
The ion cyclotron dispersion relation in a proton-alpha solar wind
Solutions to the WKB warm plasma dispersion relation are investigated for the case of a proton-alpha plasma, to determine when parallel ion cyclotron waves can propagate at the alpha particle gyrofrequency. Examples are given of the behavior of the dispersion curves in this resonant regime, for several types of plasma states that are appropriate to the fast solar wind. It is found that while ion cyclotron waves can propagate at the alpha particle gyrofrequency for a bounded set of plasma parameters, use of the correct dispersion relation will probably make resonant cyclotron acceleration of solar wind plasma particles less feasible as a mechanism for the generation of the observed preferential effects.
Isenberg, P. A.↗
Resistive wave dissipation on magnetic inhomogeneities Normal modes and phase mixing
Numerical solutions of the linearized, resistive MHD equations indicate the presence of two characteristic forms of wave decay on magnetic inhomogeneities. Solutions for relatively long-wavelength disturbances and large values of the resistivity have the behavior of decaying normal modes. In the opposite limits, phase mixing becomes important, with the eventual build-up of large spatial gradients in which the resistive terms dominate. This composite behavior is shown to be conceptually consistent with analytic results, which predict that the complete solution consists of the sum of collective and noncollective contributions. The numerical simulations go beyond analytic theory by defining where each contribution prevails and computing wave decay when phase mixing is effective. The dispersion relation in the normal-mode regime is also determined analytically, with some approximations, and is in good agreement with the numerical solution. The decay time of the normal modes varies with Lundquist number at S1/6, while the phase-mixed decay time scales as S1/3.
Steinolfson, R. S.↗
Simulation studies of plasma waves in the electron foreshock - The transition from reactive to kinetic instability
Particle simulation experiments were used to analyze the electron beam-plasma instability. It is shown that there is a transition from the reactive state of the electron beam-plasma instability to the kinetic instability of Langmuir waves. Quantitative tests, which include an evaluation of the dispersion relation for the evolving non-Maxwellian beam distribution, show that a quasi-linear theory describes the onset of this transition and applies again fully to the kinetic stage. This stage is practically identical to the late stage seen in simulations of plasma waves in the electron foreshock described by Dum (1990).
Dum, C. T.↗
Algorithmic Extensions of Low-Dispersion Scheme and Modeling Effects for Acoustic Wave Simulation
Accurate computation of acoustic wave propagation may be more efficiently performed when their dispersion relations are considered. Consequently, computational algorithms which attempt to preserve these relations have been gaining popularity in recent years. In the present paper, the extensions to one such scheme are discussed. By solving the linearized, 2-D Euler and Navier-Stokes equations with such a method for the acoustic wave propagation, several issues were investigated. Among them were higher-order accuracy, choice of boundary conditions and differencing stencils, effects of viscosity, low-storage time integration, generalized curvilinear coordinates, periodic series, their reflections and interference patterns from a flat wall and scattering from a circular cylinder. The results were found to be promising en route to the aeroacoustic simulations of realistic engineering problems.
Kaushik, Dinesh K.↗
Propagation of low-frequency electrostatic waves in a magnetic field
LF electrostatic waves oblique to magnetic field in plasma, evaluating dispersion relation near ion cyclotron frequency
Dory, R. A.↗
Shuttle wave experiments
Wave experiments on shuttle are needed to verify dispersion relations, to study nonlinear and exotic phenomena, to support other plasma experiments, and to test engineering designs. Techniques based on coherent detection and bistatic geometry are described. New instrumentation required to provide modules for a variety of missions and to incorporate advanced signal processing and control techniques is discussed. An experiment for Z to 0 coupling is included.
Calvert, W.↗
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↗
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.↗
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.↗
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.↗
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.↗
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↗
Trapped low frequency waves propagating along a warm plasma cylinder
Trapped LF ion wave propagation along warm plasma cylinder imbedded in dielectric, obtaining dispersion relation
Klevans, E. H.↗