Search NASA⌕ Search

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.

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

At least 307 records · Page 17

Analysis of a Bank Control Guidance for Aerocapture at Uranus

We apply the Fully Numerical Predictor-corrector Aerocapture Guidance to capture a vehicle into orbit for a mission at Uranus. Using the Genesis flight mechanics simulation, we analyze both undispersed and dispersed trajectories in order to tune the guidance parameters. We then assess the performance of the guidance using Monte Carlo analyses. Properly accounting for the oblateness of Uranus within the guidance proves to be critical. We identify high-to-low density gradients as the cause of large target orbit misses. Intentionally targeting a steeper flight-path angle at entry interface mitigates the risk of exiting the atmosphere on a hyperbolic orbit. Finally, We deliver the guidance code and tuned parameters for integration into the Program to Optimize Simulated Trajectories II.

Daniel A Matz↗

Electromagnetic wave propagation in a magneto-plasma filled coaxial structure. I - Theoretical. II - Experimental

This study is concerned with the problem of electromagnetic wave propagation in a magneto-plasma filled coaxial structure. The problem is formulated using the classical boundary value problem approach. A numerical investigation shows the existence of propagating slow modes, backward modes, a quasi-TEM mode, and waveguide-type modes in a magneto-plasma filled coaxial structure. Dispersion curves for these different modes are presented. Measurements have been made of electromagnetic propagation in a coaxial electrode structure filled with longitudinally magnetized plasma. The annular plasma region had a 9.55 cm outer diameter, a 3.82 cm inner diameter and was approximately 60 cm long. A magnetic field of 300 gauss was employed. Electromagnetic wave frequencies were in the range .5 to 2.4 GHz. The plasma was generated by a continuous glow discharge. The resulting dispersion curves closely follow the predicted curves for the quasi-TEM mode.

Askins, H. W., Jr.↗

MPD thruster plasma instability studies

Results of an ongoing experimental/theoretical effort to understand the influence of plasma waves and instabilities upon coaxial plasma thruster performance and longevity are discussed. An experiment which permits measurement of the plasma dispersion relation plus a suite of plasma parameters for each one millisecond quasi-steady discharge is described. This investigation is guided by a linearized Boltzmann-Poisson kinetic numerical model which can accommodate the effects of collisions, inhomogeneities and the magnetic field. Preliminary tests reveal the presence of the low-frequency branch of an ion cyclotron wave. Although the stability of these waves still has to be established through spatial growth rate measurements, it is well known that they are prone to 'current-driven' instabilities. 'Current-driven' instabilities in the collisional tensor conducting plasma of the MPD thruster are analytically discussed and their relevance to the MPD plasma is illustrated by estimating a critical total current at which their onset occurs.

Kelly, A. J.↗

A globally well-posed finite element algorithm for aerodynamics applications

A finite element CFD algorithm is developed for Euler and Navier-Stokes aerodynamic applications. For the linear basis, the resultant approximation is at least second-order-accurate in time and space for synergistic use of three procedures: (1) a Taylor weak statement, which provides for derivation of companion conservation law systems with embedded dispersion-error control mechanisms; (2) a stiffly stable second-order-accurate implicit Rosenbrock-Runge-Kutta temporal algorithm; and (3) a matrix tensor product factorization that permits efficient numerical linear algebra handling of the terminal large-matrix statement. Thorough analyses are presented regarding well-posed boundary conditions for inviscid and viscous flow specifications. Numerical solutions are generated and compared for critical evaluation of quasi-one- and two-dimensional Euler and Navier-Stokes benchmark test problems.

Iannelli, G. S.↗

FAST Observations of Acceleration Processes in the Cusp--Evidence for Parallel Electric Fields

The existence of precipitating keV ions in the Earth's cusp originating at the magnetosheath provide unique means to test our understanding of particle acceleration and parallel electric fields in the lower altitude acceleration region. On numerous occasions, the FAST (The Fast Auroral Snapshot) spacecraft has encountered the Earth's cusp regions near its apogee of 4175 km which are characterized by their signatures of dispersed keV ion injections. The FAST instruments also reveal a complex microphysics inherent to many, but not all, of the cusp regions encountered by the spacecraft, that include upgoing ion beams and conics, inverted-V electrons, upgoing electron beams, and spikey DC-coupled electric fields and plasma waves. Detailed inspection of the FAST data often show clear modulation of the precipitating magnetosheath ions that indicate that they are affected by local electric potentials. For example, the magnetosheath ion precipitation is sometimes abruptly shut off precisely in regions where downgoing localized inverted-V electrons are observed. Such observations support the existence of a localized process, such as parallel electric fields, above the spacecraft which accelerate the electrons downward and consequently impede the precipitating ion precipitation. Other acceleration events in the cusp are sometimes organized with an apparent cellular structure that suggests Alfven waves or other large-scale phenomena are controlling the localized potentials. We examine several cusp encounters by the FAST satellite where the modulation of energetic session on acceleration particle populations reveals evidence of localized acceleration, most likely by parallel electric fields.

Pfaff, R. F.. Jr.↗

Effects of Gravity on Sheared Turbulence Laden with Bubbles or Droplets

The objective of this numerical/experimental study is to improve the understanding of the effects of gravity on the two-way interaction between dispersed particles (bubbles or liquid droplets) and the carrier turbulent flow. The first phase of the project considers isotropic turbulence. Turbulent homogeneous shear flows laden with droplets/bubbles will be studied in the next phase. The experiments reported here are concerned with the dispersion of liquid droplets by homogeneous turbulence under various gravitational conditions and the effect of these droplets on the evolution of the turbulence of the carrier fluid (air). Direct numerical simulations (DNS) of bubble - laden isotropic decaying turbulence are performed using the two-fluid approach (TF) instead of the Eulerian-Lagrangian approach (EL). The motivation for using the TF formulation is that EL requires considerable computational resources especially for the case of two-way coupling where the instantaneous trajectories of a large number of individual bubbles need to be computed. The TF formulation is developed by spatially averaging the instantaneous equations of the carrier flow and bubble phase over a scale of the order of the Kolmogorov length scale which, in our case, is much larger than the bubble diameter. On that scale, the bubbles are treated as a continuum (without molecular diffusivity) characterized by the bubble phase velocity field and concentration (volume fraction). The bubble concentration, C, is assumed small enough to neglect the bubble-bubble interactions.

Elghobashi, Said↗

WRLES: Wave Resolving Large-Eddy Simulation Code, Theory and Usage

A computational fluid dynamics code has been developed for large-eddy simulations (LES) of turbulent flow. The code uses high-order of accuracy and high-resolution numerical methods to minimize solution error and maximize the resolution of the turbulent structures. Spatial discretization is performed using explicit central differencing. The central differencing schemes in the code include 2nd- to 12th-order standard central difference methods as well as 7-, 9-, 11- and 13-point dispersion relation preserving schemes. Solution filtering and high-order shock capturing are included for stability. Time discretization is performed using multistage Runge-Kutta methods that are up to 4th order accurate. Several options are available to model turbulence including: Baldwin-Lomax and Spalart-Allmaras Reynolds-averaged Navier-Stokes turbulence models, and Smagorinsky, Dynamic Smagorinsky and Vreman sub-grid scale models for LES. This report presents the theory behind the numerical and physical models used in the code and provides a user's manual to the operation of the code.

High-Order↗

Intrinsic Carrier Concentration and Electron Effective Mass in Hg(1-x) Zn(x) Te

In this work, the intrinsic carrier concentration and electron effective mass in Hg(l-x)Zn(x)Te were numerically calculated. We adopt the procedures similar to those used by Su et. al. for calculating intrinsic carrier concentrations in Hg(1-x)Cd(x)Te which solve the exact dispersion relation in Kane model for the calculation of the conduction band electron concentrations and the corresponding electron effective masses. No approximation beyond those inherent in the k centered dot p model was used here.

Sha, Yi-Gao↗

Outer Solar System Material in Inner Solar System Regolith Breccias

There is excellent evidence that a dynamical instability in the early solar system led to gravitational interactions between the giant planets and trans-Neptunian objects (TNOs). Giant planetary migration triggered by the instability dispersed a disk of primordial TNOs and created a number of small body reservoirs (e.g. the Kuiper Belt, scattered disk, irregular satellites, and the Jupiter/Neptune Trojan populations). It also injected numerous bodies into the main asteroid belt, where modeling shows they can successfully reproduce the observed P and D-type asteroid populations. During the injection era and after implantation, some of these “main belt TNOs” would have collided with S-class asteroids. Some of this material may have survived as a component of asteroid regolith breccias. Thus, we have been searching for evidence of these impact events in the form of carbonaceous xenoliths in brecciated ordinary chondrites. These xenoliths would have experienced a wide range of impact velocities, and therefore we should expect to see everything between relatively unaltered material to completely shock-melted lithologies. This material might also be different from the carbonaceous chondrites that represent standard C-complex asteroids. A goal of this research is to define useful criteria for distinguishing between these two classes of materials, including O, Cr, N and C isotopes, petrographic characteristics, and chronology.

Zolensky, M.↗

Vlasov equation eigenvalues and eigenvectors for Fourier-Hermite dispersion matrices of order greater than 1,000

The connection between the Van Kampen and Landau representations of the Vlasov equations has been extended to Fourier-Hermite expansions containing more than 1000 terms by taking advantage of the properties of tridiagonal matrices. These numerical results are regarded as conclusive indications of the nonuniformly convergent behavior of the approximation curve in the limit of an infinite number of terms and represent an extension of work begun by Grant (1967) and by Grant and Feix (1967).

Grant, F. C.↗

Analysis of the error associated with grid representation of point sources

The steady-state, three-dimensional transport equation for inert pollutant dispersion in the atmosphere is solved analytically to obtain expressions for pollutant concentrations from both a point source and a rectangular area source. The results of this analysis have application to numerical grid models in which the smallest resolvable source is the size of a single grid. The concentrations at ground level along the plume centerline from the two types of sources are in close agreement with each other at distances greater than 10 km downwind when the horizontal dimension of the rectangular source is not less than 100 m. The dimension of the source in the vertical direction is held constant at 50 m. Similar results are obtained for concentrations along the plume centerline at the effective stack height. When the size of the source in the crosswind direction is increased to 1 km, good agreement between the two concentrations is obtained for downwind distances over 60 km for the ground level case and over 140 km for the concentration at stack height level. For the horizontal dimension equal to 10 km, good agreement between the two cases is only obtained at extremely large downwind distances at both the stack height and ground levels.

Karamchandani, P.↗

Instability of rectangular jets

The instability of rectangular jets is investigated using a vortex-sheet model. It is shown that such jets support four linearly independent families of instability waves. Within each family there are infinitely many modes. A way to classify these modes according to the characteristics of their mode shapes or eigenfunctions is proposed. It is demonstrated that the boundary element method can be used to calculate the dispersion relations and eigenfunctions of these instability wave modes. The method is robust and efficient. A parametric study of the instability wave characteristics has been carried out. A sample of the numerical results is reported here. It is found that the first and third modes of each instability wave family are corner modes. The pressure fluctuations associated with these instability waves are localized near the corners of the jet. The second mode, however, is a center mode with maximum fluctuations concentrated in the central portion of the jet flow. The center mode has the largest spatial growth rate. It is anticipated that as the instability waves propagate downstream the center mode would emerge as the dominant instability of the jet.

Tam, Christopher K. W.↗

Failure Mode Analysis of V-Shaped Pyrotechnically Actuated Valves

Current V-shaped stainless steel pyrovalve initiators have rectified many of the deficiencies of the heritage Y-shaped aluminum design. However, a credible failure mode still exists for dual simultaneous initiator (NSI) firings in which low temperatures were detected at the booster cap and less consistent ignition was observed than when a single initiator was fired. In order to asses this issue, a numerical framework has been developed for predicting the flow through pyrotechnically actuated valves. This framework includes a fully coupled solution of the gas-phase equation with a non-equilibrium dispersed phase for solid particles as well as the capability to model conjugate gradient heat transfer to the booster cap. Through a hierarchy of increasingly complex simulations, a hypothesis for the failure mode of the nearly simultaneous dual NSI firings has been proven. The simulations indicate that the failure mode for simultaneous dual NSI firings may be caused by flow interactions between the flame channels. The shock waves from each initiator interact in the booster cavity resulting in a high pressure that prevents the gas and particulate velocity from rising in the booster cap region. This impedes the bulk of the particulate phase from impacting the booster cap and reduces the heat transfer to the booster cap since the particles do not impact it. Heat transfer calculations to the solid metal indicate that gas-phase convective heat transfer may not be adequate by itself and that energy transfer from the particulate phase may be crucial for the booster cap burn through.

Sachdev, Jai S.↗

Comparison of Aerocapture Performance Using Bank Control and Direct Force Control with Two Human-Scale Vehicles at Mars

Recent studies of human-scale missions to Mars have included a wide trade space of vehicle configurations and control schemes. Some configurations fly at a low angle of attack with a low L/D,while others fly at a high angle of attack with a mid L/D. Some use bank angle control, while others use direct force control, where the angle of attack and sideslip angle are independently modulated. This paper compares the performance of three vehicle configurations: a low-L/D vehicle using direct force control, a low-L/D vehicle using bank control, and a mid-L/D vehicle using bank control. The reference mission is aerocapture at Mars into a highly elliptical, 1-sol orbit. The trajectories are integrated in three degrees of freedom. All three cases utilize numeric predictor-corrector guidances and emulate control system responses with rate and acceleration limits. The configurations are compared using a Monte Carlo analysis. The robustness of each configuration to increased dispersions is also compared.

Matz, Daniel A.↗

Entry vehicle performance analysis and atmospheric guidance algorithm for precision landing on Mars

Future missions to Mars may require pin-point landing precision, possibly on the order of tens of meters. The ability to reach a target while meeting a dynamic pressure constraint to ensure safe parachute deployment is complicated at Mars by low atmospheric density, high atmospheric uncertainty, and the desire to employ only bank angle control. The vehicle aerodynamic performance requirements and guidance necessary for 0.5 to 1.5 lift drag ratio vehicle to maximize the achievable footprint while meeting the constraints are examined. A parametric study of the various factors related to entry vehicle performance in the Mars environment is undertaken to develop general vehicle aerodynamic design requirements. The combination of low lift drag ratio and low atmospheric density at Mars result in a large phugoid motion involving the dynamic pressure which complicates trajectory control. Vehicle ballistic coefficient is demonstrated to be the predominant characteristic affecting final dynamic pressure. Additionally, a speed brake is shown to be ineffective at reducing the final dynamic pressure. An adaptive precision entry atmospheric guidance scheme is presented. The guidance uses a numeric predictor-corrector algorithm to control downrange, an azimuth controller to govern crossrange, and analytic control law to reduce the final dynamic pressure. Guidance performance is tested against a variety of dispersions, and the results from selected tests are presented. Precision entry using bank angle control only is demonstrated to be feasible at Mars.

Dieriam, Todd A.↗

Far-ultraviolet fluorescence of carbon monoxide in the red giant Arcturus. II - Analysis of high-dispersion IUE spectra

Faint, diffuse emissions near 1380 A in deeply exposed IUE spectrograms of the red giant Arcturus very likely are associated with bands of the A-X fourth-positive system of carbon monoxide, fluoresced by multiplet UV2 of neutral oxygen near 1305 A. Numerical simulations indicate that the strength of the CO bands is exceedingly sensitive, in the best available one-dimensional model of the chromosphere of Arcturus, to a delicate balance between the rapid inward attenuation of the oxygen radiation field and the rapid outward decline of the molecular absorptivity. The fortuitous character of the overlap region in the single-component model argues that one should also consider the possibility that the pumping occurs in a highly inhomogeneous chromosphere, of the type proposed in previous studies of Arcturus based on observations of the infrared absorption bands of CO.

Ayres, T. R.↗

Non-symmetric two-stream instability

A theoretical investigation is performed concerning the instability spectrum of quasi-electrostatic waves at shifted half-odd-integer values of the cyclotron frequency due to nonsymmetric counterstreaming electron beams. The beam velocities parallel and perpendicular to the static magnetic field are represented by double Dirac delta functions with no imposition of parameter limitations. Systematic consideration is given to the coupling between plasma modes and cyclotron modes as well as the coupling between cyclotron modes of the two beams that result in shifted half-odd-integer multiples of the cyclotron frequency. The general dispersion equation for quasi-electrostatic waves is analyzed, and plasma-cyclotron coupling in the nonsymmetric case is treated by deriving approximate analytical expressions for maximum growth rates and marginal stability. Exact numerical solutions in both frequency and wavenumber space are obtained and compared with the analytical expressions. Cyclotron-cyclotron coupling modes are treated in the same way, and the results for both types of coupling are compared. It is found that certain modes may be weakened or completely suppressed when the Bessel functions for specific plasma parameters vanish.

Cuperman, S.↗

Analysis of Jovian low frequency radio emissions

The density of ions in the Io plasma torus and the scattering of these ions by low frequency electromagnetic emissions detected by Voyager 1 were studied. The ion density profile was investigated using whistler dispersion measurements provided by the Voyager plasma instrument. The scale height and absolute density of H+ ions in the vicinity of the plasma torus were determined by combining the measured plasma densities with the whistler dispersion measurements. A theoretical analysis of the modes of propagation of low frequency electromagnetic emissions in the torus was undertaken. Polarization reversal effects and rough estimates of the ion diffusion coefficient were utilized. Numerical evaluation of the ion diffusion coefficients in the torus were made using the observed Voyager 1 wave intensities. Results show that the observed wave intensities produce significant ion diffusion effects in the ion torus.

Gurnett, D. A.↗