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 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Lower-hybrid instabilities driven by an ion velocity ring

The lower-hybrid instabilities in high-beta (ratio of plasma to magnetic pressure) plasmas driven by ring-ion distributions in velocity space are investigated. A dispersion equation including electromagnetic effects is derived. In the low-beta limit, analytic expressions are obtained which illuminate the physical nature of the instabilities. The complete dispersion equation is solved numerically as a function of ring speed and plasma beta for several types of ring distribution. Electromagnetic effects are important for relatively energetic rings even in the low-beta regime, suppressing growth rates and shifting the angle of propagation to more oblique angles. Stabilization by thermal effects is also discussed. Application of these results to the earth's bow shock, AMPTE, comets and solar flares is suggested.

Akimoto, K.↗

NASA NTRS1987

A linear MHD instability analysis of solar mass ejections with gravitation

The linear MHD instability of a cylindrical plasma is used to investigate the origin of solar mass ejections, and the dispersion relation is solved numerically. The initial plasma-flow velocity is found to have a significant effect on the instability criteria and growth rate, and the instability growth-rate is shown to be larger in cases where plasma flow exists, relative to the static case. Results suggest that the plasma column may break into small pieces. Assuming a thin-tube approximation, gravity is found to have little effect on the instability of quasi-horizontal ejection, but to have considerable effect on the vertical ejection. In considering the gravitational force, an exact analytical solution is found for the vertical case, while asymptotic solutions are given for the horizontal and oblique cases.

Song, M. T.↗

NASA NTRS1988

Growth of electron plasma waves above and below f(p) in the electron foreshock

This paper investigates the conditions required for electron beams to drive wave growth significantly above and below the electron plasma frequency, f(p), by numerically solving the linear dispersion equation. It is shown that kinetic growth well below f(p) may occur over a broad range of frequencies due to the beam instability, when the electron beam is slow, dilute, and relatively cold. Alternatively, a cold or sharp feature at low parallel velocities in the distribution function may drive kinetic growth significantly below f(p). Kinetic broadband growth significantly above f(p) is explained in terms of faster warmer beams. A unified qualitative theory for the narrow-band and broad-band waves is proposed.

Cairns, Iver H.↗

NASA NTRS1989

Effects of unequal particle number densities on Alfven waves

Analytic plasma theory and numerical solutions of the dispersion equation are used to show that the assumption that the linear properties of the waves are determined by a charge-neutral plasma in the absence of the nonthermal particles, while the nonthermal particles cause growth or additional damping superposed onto the background, is seriously flawed even for stable plasmas. Even when the nonthermal particles do not contribute significantly to the dispersion equation, unequal thermal electron and ion number densities (due to the presence of the nonthermal particles) may cause fundamental low wave number modifications to the Alfven modes, including the creation of a new resonance and severely modified dispersion. These results are found for both cold and warm plasmas. Previous work on Alfven waves should be reevaluated in view of these results.

Cairns, I. H.↗

NASA NTRS1990

A current disruption mechanism in the neutral sheet - A possible trigger for substorm expansions

A linear analysis is performed to investigate the kinetic cross-field streaming instability in the earth's magnetotail neutral sheet region. Numerical solution of the dispersion equation shows that the instability can occur under conditions expected for the neutral sheet just prior to the onset of substorm expansion. The excited waves are obliquely propagating whistlers with a mixed polarization in the lower hybrid frequency range. The ensuing turbulence of this instability can lead to a local reduction of the cross-tail current causing it to continue through the ionosphere to form a substorm current wedge. A substorm expansion onset scenario is proposed based on this instability in which the relative drift between ions and electrons is primarily due to unmagnetized ions undergoing current sheet acceleration in the presence of a cross-tail electric field. The required electric field strength is within the range of electric field values detected in the neutral sheet region during substorm intervals. The skew in local time of substorm onset location and the three conditions under which substorm onset is observed can be understood on the basis of the proposed scenario.

Lui, A. T. Y.↗

NASA NTRS1993

Fuel injector: Air swirl characterization aerothermal modeling, phase 2, volume 1

A well integrated experimental/analytical investigation was conducted to provide benchmark quality relevant to a prefilming type airblast fuel nozzle and its interaction with the combustor dome air swirler. The experimental investigation included a systematic study of both single-phase flows that involved single and twin co-axial jets with and without swirl. A two-component Phase Doppler Particle Analyzer (PDPA) was used to document the interaction of single and co-axial air jets with glass beads that simulate nonevaporating spray and simultaneously avoid the complexities associated with fuel atomization processes and attendant issues about the specification of relevant boundary conditions. The interaction of jets with methanol spray produced by practical airblast nozzle was also documented in the spatial domain of practical interest. Model assessment activities included the use of three turbulence models (k-epsilon, algebraic second moment (ASM), and differential second moment (DSM)) for the carrier phase, deterministic or stochastic Lagrangian treatment of the dispersed phase, and advanced numerical schemes. Although qualitatively good comparison with data was obtained for most of the cases investigated, the model deficiencies in regard to modeled dissipation rate transport equation, single length scale, pressure-strain correlation, and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗

NASA NTRS1993

Fuel Injector: Air swirl characterization aerothermal modeling, phase 2, volume 2

A well integrated experimental/analytical investigation was conducted to provide benchmark quality data relevant to prefilming type airblast fuel nozzle and its interaction with combustor dome air swirler. The experimental investigation included a systematic study of both single-phase flows that involved single and twin co-axial jets with and without swirl. A two-component Phase Doppler Particle Analyzer (PDPA) equipment was used to document the interaction of single and co-axial air jets with glass beads that simulate nonevaporating spray and simultaneously avoid the complexities associated with fuel atomization processes and attendant issues about the specification of relevant boundary conditions. The interaction of jets with methanol spray produced by practical airblast nozzle was also documented in the spatial domain of practical interest. Model assessment activities included the use of three turbulence models (k-epsilon, algebraic second moment (ASM) and differential second moment (DSM)) for the carrier phase, deterministic or stochastic Lagrangian treatment of the dispersed phase, and advanced numerical schemes. Although qualitatively good comparison with data was obtained for most of the cases investigated, the model deficiencies in regard to modeled dissipation rate transport equation, single length scale, pressure-strain correlation, and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗

NASA NTRS1996

The Effects of Dissipation and Coarse Grid Resolution for Multigrid in Flow Problems

The objective of this paper is to investigate the effects of the numerical dissipation and the resolution of the solution on coarser grids for multigrid with the Euler equation approximations. The convergence is accomplished by multi-stage explicit time-stepping to steady state accelerated by FAS multigrid. A theoretical investigation is carried out for linear hyperbolic equations in one and two dimensions. The spectra reveals that for stability and hence robustness of spatial discretizations with a small amount of numerical dissipation the grid transfer operators have to be accurate enough and the smoother of low temporal accuracy. Numerical results give grid independent convergence in one dimension. For two-dimensional problems with a small amount of numerical dissipation, however, only a few grid levels contribute to an increased speed of convergence. This is explained by the small numerical dissipation leading to dispersion. Increasing the mesh density and hence making the problem over resolved increases the number of mesh levels contributing to an increased speed of convergence. If the steady state equations are elliptic, all grid levels contribute to the convergence regardless of the mesh density.

Eliasson, Peter↗

NASA NTRS2015

A Comparison of Three Algorithms for Orion Drogue Parachute Release

The Orion Multi-Purpose Crew Vehicle is susceptible to ipping apex forward between drogue parachute release and main parachute in ation. A smart drogue release algorithm is required to select a drogue release condition that will not result in an apex forward main parachute deployment. The baseline algorithm is simple and elegant, but does not perform as well as desired in drogue failure cases. A simple modi cation to the baseline algorithm can improve performance, but can also sometimes fail to identify a good release condition. A new algorithm employing simpli ed rotational dynamics and a numeric predictor to minimize a rotational energy metric is proposed. A Monte Carlo analysis of a drogue failure scenario is used to compare the performance of the algorithms. The numeric predictor prevents more of the cases from ipping apex forward, and also results in an improvement in the capsule attitude at main bag extraction. The sensitivity of the numeric predictor to aerodynamic dispersions, errors in the navigated state, and execution rate is investigated, showing little degradation in performance.

Matz, Daniel A.↗

NASA NTRS1986

Numerical analysis of confined recirculating gas-solid turbulent flows

Particle dispersion in confined recirculating turbulent flows has been investigated numerically. The present computational scheme utilizes Eulerian forms of the governing equations and allows two-way coupling between gas and solid phases. A recently developed two-phase closure model based on the multiple-scale turbulence model was used for the estimation of turbulent viscosities and diffusivities. For the particle size and loading considered in this study, the effect of particulate phase on the conveying gas is small, the nonequilibrium between the two phases is significant. Effects of recirculation, expanded chamber size and secondary annular jet momentum on the particle mixing rate are also investigated. In general, the present numerical results are in reasonably good agreement with the available experimental data.

Chen, C. P.↗

DOE OSTI2024

Rheology of bi-disperse dense fiber suspensions

While significant progress has been made in the modeling and simulation of uniform fiber suspensions, no existing model has been validated for industrially-relevant concentrated suspensions containing fibers of multiple aspect ratios. In the present work, we investigate bi-disperse suspensions with two fiber populations in varying aspect ratios in a steady shear flow using direct numerical simulations. Moreover, we measure the suspension viscosity by creating a controlled length bidispersity for nylon fibers suspended in a Newtonian fluid. Here the results showed good agreement between the experimentally measured and numerically predicted viscosity for bi-disperse suspensions. The ratio between the aspect ratio of large to small fibers (size ratio) and the volume fraction of large fibers (composition) in bi-disperse systems strongly affected the rheological behavior of the suspension. The increment of relative viscosity associated with size ratio and composition can be explained by the decrease in the maximum flowable limit or jamming volume fraction. Moreover, the relative viscosity of bi-disperse suspensions collapses, when plotted against the reduced volume fraction, demonstrating the controlling influence of the jamming fraction in bi-disperse fiber suspensions.

42 ENGINEERING↗

NASA NTRS1994

Electron time dispersion

Observations of electrons during a sounding rocket flight from Sondrestromfjord, Greenland, in 1985 showed numerous instances of electron time-of-flight dispersion, that is, a coherent structure in which fast electrons arrived before slow electrons. The instrumentation measured many pitch angles simultaneously with high time resolution, providing a detailed view of the temporal evolution of the electron signature at rocket altitudes which was indicative of an impulsive injection of electrons over a range of energies. A model utilizing an extended, quasi-static (static on a timescale of 1 s), parallel electric field was developed to investigate a possible mechanism for the injection of these electrons. The model is found to fit the data well only when a background is added which cannot be self-consistently generated by the model, although other parameters used in the model are in general agreement with previous reports. This leads to the conclusion that electron dispersion is more probably generated by a wave acceleration mechanism.

Kletzing, C. A.↗

NASA NTRS1995

Absorbing Boundary Conditions For Optical Pulses In Dispersive, Nonlinear Materials

This paper will present results in computational nonlinear optics. An algorithm will be described that provides absorbing boundary conditions for optical pulses in dispersive, nonlinear materials. A new numerical absorber at the boundaries has been developed that is responsive to the spectral content of the pulse. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of "light bullet" like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. Comparisons will be shown of calculations that use the standard boundary conditions and the new ones.

Goorjian, Peter M.↗

NASA NTRS1977

A simulation of air pollution model parameter estimation using data from a ground-based LIDAR remote sensor

One way to obtain estimates of the unknown parameters in a pollution dispersion model is to compare the model predictions with remotely sensed air quality data. A ground-based LIDAR sensor provides relative pollution concentration measurements as a function of space and time. The measured sensor data are compared with the dispersion model output through a numerical estimation procedure to yield parameter estimates which best fit the data. This overall process is tested in a computer simulation to study the effects of various measurement strategies. Such a simulation is useful prior to a field measurement exercise to maximize the information content in the collected data. Parametric studies of simulated data matched to a Gaussian plume dispersion model indicate the trade offs available between estimation accuracy and data acquisition strategy.

Kibler, J. F.↗

NASA NTRS1975

Development of a three dimensional numerical water quality model for continental shelf applications

A model to predict the distribution of water quality parameters in three dimensions was developed. The mass transport equation was solved using a non-dimensional vertical axis and an alternating-direction-implicit finite difference technique. The reaction kinetics of the constituents were incorporated into a matrix method which permits computation of the interactions of multiple constituents. Methods for the computation of dispersion coefficients and coliform bacteria decay rates were determined. Numerical investigations of dispersive and dissipative effects showed that the three-dimensional model performs as predicted by analysis of simpler cases. The model was then applied to a two dimensional vertically averaged tidal dynamics model for the Providence River. It was also extended to a steady state application by replacing the time step with an iteration sequence. This modification was verified by comparison to analytical solutions and applied to a river confluence situation.

Spaulding, M.↗

DOE OSTI2025

Floquet engineering of polaritonic amplification in dispersive photonic time crystals

In this study, we investigate the dynamics of dispersive photonic time crystals (PTCs) and their potential applications for controlling light-matter interaction. Using the Lorentz-Drude as dispersive material model, we analyze the effects of periodic modulation of dispersion parameters through both theoretical and numerical analysis, revealing the emergence of hybrid bandgaps from the interaction of polaritonic branches with distinctive characteristics. Our study demonstrates that dispersive PTCs offer novel excitation channels and amplification possibilities, that require lower modulation frequencies compared to nondispersive systems thus alleviating experimental challenges for the realization of PTCs in the optical regime. These findings pave the way for advancements in polaritonic lasing and resonant Raman scattering.

Classical optics↗

NASA NTRS1983

Linear density drift instabilities in very low beta plasmas A different approach

This paper reports a study of the linear Vlasov electromagnetic dispersion equation for density drift instabilities in very low beta plasmas. A uniform magnetic field and a weak uniform density gradient are assumed. This paper differs from most other studies of this topic in three important ways: (1) no low frequency or long wavelength approximations are made, (2) no gauge condition is imposed and (3) a modification of the local approximation is used which is argued to be less arbitrary than the usual local approximation. Numerical solution of the resulting dispersion equation yields stabilization of the universal density drift instability at values of beta somewhat higher than those obtained from the local approximation, and near maximum growth shows no evidence of coupling between the universal and Alfven modes.

Gary, S. P.↗

DOE OSTI2023

Computation of high-frequency magnetoelastic waves in layered materials

Here, the direct calculation of magnetoelastic wave dispersion in layered media is presented using an efficient, accurate computational technique. The governing, coupled equations for elasticity and magnetism, the Navier and Landau-Lifshitz equations, respectively, are linearized to form a quadratic eigenvalue problem that determines a complex web of wave-number–frequency dispersion branches and their corresponding mode profiles. Numerical discretization of the eigenvalue problem via a spectral collocation method (SCM) is employed to determine the complete dispersion maps for both a single, finite-thickness magnetic layer and a finite magnetic-nonmagnetic double-layer. The SCM, previously used to study elastic waves in nonmagnetic media, is fast, accurate, and adaptable to a variety of sample configurations and geometries. Emphasis is placed on the extremely high-frequency regimes being accessed in ultrafast magnetism experiments. The dispersion maps and modes provide insight into how energy propagates through the coupled system, including how energy can be transferred between elastic- and magnetic-dominated waves as well as between different layers. The numerical computations for a single layer are further understood by a simplified analytical calculation in the high-frequency, exchange-dominated regime where the resonance condition required for energy exchange (an anticrossing) between quasi-elastic and quasi-magnetic dispersion branches is determined. Nonresonant interactions are shown to be well approximated by the dispersion of uncoupled elastic and magnetic waves. The methods and results provide fundamental theoretical tools to model and understand current and future magnetic devices powering spintronic innovation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗