SEARCH · Search NASA

Results for “numerical dispersion”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Electromagnetic electron beam instabilities - Hot, isotropic beams

This paper considers the linear theory of electromagnetic instabilities driven by an electron beam in a homogeneous, nonrelativistic, Vlasov plasma. The beam is relatively hot, isotropic in its own frame, and streams parallel or antiparallel to a magnetic field B. Numerical solutions of the full dispersion equation for propagation parallel or antiparallel to B are presented, and the linear properties of the whistler heat flux and electron beam firehose instabilities are exhibited and compared. Under a broad range of parmameters the former mode has the lower beam speed threshold, and the larger maximum growth rate. In addition, it is demonstrated that, for a sufficiently large relative beam density, relative beam temperature, and plasma beta the whistler heat flux instability has a much lower beam speed threshold than the electrostatic electron beam instability. The application of these instabilities to first-order Fermi acceleration of electrons at space plasma shocks is discussed.

Gary, S. P.↗

NASA NTRS1986

Low-frequency waves in a high-beta collisionless plasma Polarization, compressibility and helicity

This paper considers the linear theory of waves near and below the ion cyclotron frequency in an isothermal electron-ion Vlasov plasma which is isotropic, homogeneous and magnetized. Numerical solutions of the full dispersion equation for the magnetosonic/whistler and Alfven/ion cyclotron modes at beta(i) = 1.0 are presented, and the polarizations, compressibilities, helicities, ion Alfven ratios and ion cross-helicities are exhibited and compared. At sufficiently large beta(i) and theta, the angle of propagation with respect to the magnetic field, the real part of the polarization of the Alfven/ion cyclotron wave changes sign, so that, for such parameters, this mode is no longer left-hand polarized. The Alfven/ion cyclotron mode becomes more compressive as the wavenumber increases, whereas the magnetosonic/whistler becomes more compressive with increasing theta.

Gary, S. P.↗

NASA NTRS1986

Oscillations of a viscous compound drop

The effect of viscosity on normal mode oscillations of a compound drop is determined for a range of values of the viscosity, and the effect of a host medium is considered. The general form of the dispersion relation is derived, and numerical values are obtained for the case of an annular liquid region composed of silicon oil and with the core and host regions taken to be air. Results indicate that the sloshing mode of the compound drop is more damped than the bubble mode.

Lyell, M. J.↗

NASA NTRS1987

On the thermal instability of galactic and cluster halos

The paper presents a detailed study of thermal instabilities in cooling flows associated with galaxies and clusters of galaxies. In the case of purely radiation-driven accretion onto a central object such as the cD galaxy M87, it is found that the gas is largely subject to overstability, rather than to monotonic instability. If thermal conductivity is taken into account, the flow is stabilized on scales of several kiloparsecs, even if the conductivity is appreciably reduced (e.g., about 1 percent) with respect to the Spitzer value. In no case are the globular perturbations (i.e., perturbations with comparable radial and azimuthal dimensions) found to be monotonically unstable. The paper presents numerical solutions of the local dispersion relation for the cooling flow in M87 and discusses the possible consequences of the results for a correct understanding of cooling flows.

Malagoli, A.↗

NASA NTRS1988

Comet-solar wind interaction through ion-proton beam instability

The acceleration of cometary ions in the distant cometary tail by an electromagnetic beam instability, which is caused by the relative drift velocity between the cometary ions and solar protons is studied. The linear phase of the instability is analyzed by numerical evaluation of the Vlasov dispersion relation. A particle simulation is used to study the nonlinear phase of the instability. Right-hand polarized cyclotron waves are made unstable by the presence of ion beams. The cometary tail ions are accelerated mainly during the reduced-growth phase of the instability. The cyclotron waves cause momentum transfer from solar wind protons to cometary tail ions. The acceleration of cometary ions through the electromagnetic beam instability is found to be 500-5000 cm/s-squared.

Price, Channon P.↗

NASA NTRS1988

Nature and the nonlinear evolution of electrostatic waves associated with the AMPTE solar wind releases

The nonlinear evolution of the electrostatic waves observed during the AMPTE (Active Magnetosphere Particle Tracer Explorers) solar wind releases is investigated. Previous linear studies indicate two distinct sets of instabilities may be responsible for the generation. One set consists of ion-acoustic type instabilities which are insensitive to the presence of a background magnetic field, while the other group corresponds to the modified two stream instabilities and requires the solar wind flow to be across the ambient magnetic field. To establish which set of instabilities is more viable a detailed linear Vlasov theory has been conducted by numerically solving the full electromagnetic dispersion relation. In addition both the plasma wave as well as the magnetic field measurements by the IRM (Ion Release Module) spacecraft were used to correlate the frequency and the power of the observed waves with the magnitude and the direction of the solar wind magnetic field. The results of these analyses indicate that the ion-acoustic type instabilities have growth rates that are an order of magnitude or more larger than those of the modified two stream instabilities. Results show that both the solar wind protons and the released ions may be heated and accelerated in the directions oblique to the solar wind flow velocity.

Omidi, N.↗

NASA NTRS1989

The oblique behavior of low-frequency electromagnetic waves excited by newborn cometary ions

The free energy in oxygen or hydrogen ions freshly created in the solar wind stimulates low-frequency electromagnetic waves whose growth does not always maximize at parallel propagation. Exploration of the wave vector plane discloses the frequent occurrence of islets of oblique growth unconnected to the unstable parallel modes. Contour plots of the growth rate, real frequency, polarization, and magnetic compression characterize the oblique wave behavior for large values of the initial pitch angle of the cometary particles. Although wave-particle (Landau and cyclotron) resonances feed most of the surveyed oblique instabilities, some are seemingly fluidlike. The results, obtained from the numerical solution of the kinetic dispersion and wave equations, imply that newborn ions can easily excite significant oblique hydromagnetic wave activity. Cometary environments provide the adopted plasma model, but the study is helpful in the interpretation of other low-frequency wave observations in space.

Brinca, Armando L.↗

NASA NTRS1991

Flux-corrected transport techniques for multidimensional compressible magnetohydrodynamics

A prescription is given for conservatively integrating generalized hydromagnetic equations using flux-corrected transport (FCT) techniques. By placing the magnetic-field components at the interface locations of the finite-difference grid, the field is kept divergence-free to within machine roundoff error. The use of FCT techniques allows an integration scheme of high accuracy to be employed, while the numerical ripples associated with large dispersion errors are avoided. The method is particularly well suited for problems involving magnetohydrodynamic shocks and other discontinuities.

Devore, C. Richard↗

NASA NTRS1991

Feedback instability of the ionospheric resonant cavity

A model is developed that provides a theoretical basis for previous numerical results showing a feedback instability with frequencies characteristic of Alfven travel times within the region of the large increase of Alfven speed above the ionosphere. These results have been extended to arbitrary ionospheric conductivity by developing a numerical solution of the cavity dispersion relation that involves Bessel functions of complex order and argument. It is concluded that the large contrast between the magnetospheric and ionospheric Alfven speed leads to the formation of resonant cavity modes with frequencies ranging from 0.1 to 1 Hz. The presence of the cavity leads to a modification of the reflection characteristics of Alfven waves with frequencies that compare to the cavity's normal modes.

Lysak, Robert L.↗

NASA NTRS1991

ISEE 1 observations of electrostatic ion cyclotron waves in association with ion beams on auroral field lines from about 2.5 to 4.5 R(E)

Quasi-monochromatic waves at about the hydrogen cyclotron frequency were observed as the ISEE 1 satellite traversed auroral field lines at radial distances of about 2.5-4.5 R(E) near midnight on June 19, 1981. Waves and both lower and higher frequencies were observed at higher altitudes, and possible electrostatic helium cyclotron and oxygen cyclotron waves occurred at lower altitudes. Upflowing hydrogen and oxygen beams and field-aligned currents occurred simultaneously. The features of the waves are most consistent with the current-driven mode. In addition, numerical studies of the linear dispersion relation, using parameters based on the observations, show that both the parallel and oblique two-stream modes and the ion-beam-driven modes were stable while oblique current-driven modes were unstable. The O(+) and H(+) distributions provide evidence for interactions with local electrostatic ion cyclotron waves and for the H(+)-O(+) two-stream instability at altitudes below the satellite.

Catell, C. A.↗

NASA NTRS1992

Computing high-speed flows past an oscillating cylinder near a vertical wall

A computational method to simulate unsteady flows involving moving rigid boundaries and interference has been developed. The method is used to solve inviscid equations governing the fluid flow and the dynamic equations governing the motion of rigid bodies. A second-order accurate, upwind-biased, and alterating-direction-implicit method is used to solve the governing equations of the flow. A kinematic domain decomposition (KDD) procedure is extended to treat 3D problems with a high degree of accuracy and generality. The method under consideration is applied to both transonic and supertransonic flows. Both cases involve flow past a cylinder which is forced to pitch sinusoidally near a vertical wall. Benefits of the proposed approach include accurate calculation of the flow around 3D moving multiple bodies with interference; reduction of a numerical error; in particular, the dispersion error which strongly affects wave propagation; and minimization of the phase error which is accumulated according to the time advance procedure.

Yen, Guan-Wei↗

NASA NTRS1993

Strongly nonlinear evolution of low-frequency wave packets in a dispersive plasma

The evolution of strongly nonlinear, strongly modulated wave packets is investigated in a dispersive plasma using a hybrid numerical code. These wave packets have amplitudes exceeding the strength of the external magnetic field, along which they propagate. Alfven (left helicity) wave packets show strong steepening for p < 1, while fast (fight heIicity) wave packets hardly steepen for any beta. Substantial regions of opposite helicity form on the leading side of steepened Alfven wave packets. This behavior differs qualitatively from that exhibited by the solutions to the derivative nonlinear Schrodinger (DNLS) equation.

Vasquez, Bernard J.↗

NASA NTRS2023

Oblique Propagation and Temperature Effects on the Resonant Right-Hand Ion Beam Instability

The resonant right-hand instability (RHI) is often the dominant mode driven by reflected ions upstream of Earth’s quasi-parallel bow shock. In the tradition of Peter Gary, this paper further explores the right-hand instability using numerical solutions of the plasma dispersion relation and non-linear kinetic simulations, with parameters inspired by observations from NASA’s Magnetospheric Multiscale (MMS) mission. Agreement is found between the ion distributions in the particle-in-cell simulations and Magnetospheric Multiscale spacecraft data, which show the gyrophase bunching characteristic of the instability. The non-linear structures created by right-hand instability tend to be stronger when the plasma beta is lower. These structures have sizes of around 100 to 200 ion inertial lengths perpendicular to the magnetic field, presenting planet-sized disturbances to the magnetosphere. 2d and 3D hybrid particle-in-cell simulations show that modes with a range of propagation angles oblique to the magnetic field are excited, providing a ground to understand previous statistical studies of observed foreshock waves.

Ari Le↗

NASA NTRS1993

Dispersed ion structures at the poleward edge of the auroral oval: Low-altitude observations and numerical modeling

We have compared the AUREOL 3 (A3) observations of auroral ion precipitation, particularly ion beams, with the results from the global kinetic model of magnetotail plasma of Ashour-Abdalla et al. (1993). We have identified 101 energetic keV H(+) velocity dispersed precipitating ion structures (VDIS) with fluxes above 10(exp -3) ergs./sq cm./s in the A3 record between the end of 1981 and mid-1984. These beams display a systematic increase in energy with increasing latitude and were observed in a narrow region within less than 1 deg in latitude of the polar cap boundary. The VDIS are the most distinctive feature in the auroral zone of the plasma sheet boundary layer. We report first on a statistical analysis of the possible ralationships between magnetic activity or substorm phase and the VDIS properties. Our particle simulations of the precipitating ions have been extended by using a series of modified versions of the Tsyganenko (1989) magnetic field model and by varying the cross-magnetosphere electric field. In the simulations, plasma from a mantle source is subject to strong nonlinear acceleration, forming beams which flow along the PSBL. Only 3 to 4% of these beams precipitate into the ionosphere to form the VDIS while the majority return to the equatorial plane after mirroring and form the thermalized central plasma sheet. The final energy and the dispersion of the beams in the model depend on the amplitude of the cross-tail electric field. Two unsual observations of low-energy (less than 5 keV) O(+) VDIS, shifted by 4 deg 5 deg in invariant latitude equatorward of H(+) VDIS are analyzed in detail. The sparsity of such O(+) events and the absence of the changes in the flux and frequency of occurrence indicate a solar wind origin for the plasma. Finally, large-scale kinetic modeling, even with its simplifications and assumptions (e.g., static magnetic field, solar wind source), reproduces low-altitude auroral ion features fairly well; it may therefore be presented as an appropriate framework into which data on energization and transport of the hot plasma, obtained in the equatorial plane, could be inserted in the near future.

Bosqued, J. M.↗

NASA NTRS1997

Computations of Complex Three-Dimensional Turbulent Free Jets

Three-dimensional, incompressible turbulent jets with rectangular and elliptical cross-sections are simulated with a finite-difference numerical method. The full Navier- Stokes equations are solved at low Reynolds numbers, whereas at high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate the effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporary discretization and a fourth-order compact scheme is used for spatial discretization. Although such methods are widely used in the simulation of compressible flows, the lack of an evolution equation for pressure or density presents particular difficulty in incompressible flows. The pressure-velocity coupling must be established indirectly. It is achieved, in this study, through a Poisson equation which is solved by a compact scheme of the same order of accuracy. The numerical formulation is validated and the dispersion and dissipation errors are documented by the solution of a wide range of benchmark problems. Three-dimensional computations are performed for different inlet conditions which model the naturally developing and forced jets. The experimentally observed phenomenon of axis-switching is captured in the numerical simulation, and it is confirmed through flow visualization that this is based on self-induction of the vorticity field. Statistical quantities such as mean velocity, mean pressure, two-point velocity spatial correlations and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stress equations are presented to aid in the turbulence modeling of complex jets. Simulations of circular jets are used to quantify the effect of the non-uniform curvature of the non-circular jets.

Wilson, Robert V.↗

NASA NTRS1989

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗

NASA NTRSDate not supplied

The Spacecraft Mass Balance as a Diagnostic Tool for Cabin Air Quality

The spacecraft-level mass balance is a powerful tool for diagnosing and decoupling factors influencing cabin air quality, especially with respect to trace contaminant propagation. Successful implementation of this approach relies upon accurate and temporally relevant air quality measurements, which are often challenging to attain for emerging trace contaminants. In addition to the proper interpretation of air quality data trends, a thorough understanding of subsystem-level mass transfer is required to characterize process performance. For many Environmental Control and Life Support System processes, subsystem mass transfer may be predicted based on physicochemical properties and classic unit operation design approaches. Gaps in understanding of process performance should be supplemented with thorough ground testing. A high confidence in the fidelity of one of these two aspects, air quality data or process performance, may help offset uncertainties in the other. At minimum, thresholding assumptions and numerical methods can be imposed to constrain unknown parameters within the physical envelope of the situation at hand. Finally, the complexity of cabin-wide integrated mass balances may be reduced by implementing a phenomenological approach towards subsystem discretization.

Cabin Air-quality↗

NASA NTRS1996

Decoupling the Role of Inertia and Gravity on Particle Dispersion

Turbulent gas flows laden with small, dense particles are encountered in a wide number of important applications in both industrial settings and aerodynamics applications. Particle interactions with the underlying turbulent flow are exceedingly complex and, consequently, difficult to accurately model. The difficulty arises primarily due to the fact that response of a particle to the local environment is dictated by turbulence properties in the reference frame moving with the particle (particle-Lagrangian). The particle-Lagrangian reference frame is in turn dependent upon the particle relaxation time (time constant) as well as gravitational drift. The combination of inertial and gravitational effects in this frame complicates our ability to accurately predict particle-laden flows since measurements in the particle-Lagrangian reference frame are difficult to obtain. Therefore, in this work we will examine separately the effects of inertia and gravitational drift on particle dispersion through a combination of physical and numerical experiments. In this study, particle-Lagrangian measurements will be obtained in physical experiments using stereo image velocimetry. Gravitational drift will be varied in the variable-g environments of the NASA DC-9 and in the zero-g environment at the drop tower at NASA-Lewis. Direct numerical simulations will be used to corroborate the measurements from the variable-g experiments. We expect that this work will generate new insight into the underlying physics of particle dispersion and will, in turn, lead to more accurate models of particle transport in turbulent flows.

Rogers, Chris↗