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

Modulation of energetic particle fluxes by a mixed mode of transverse and compressional waves

Modulation characteristics of particle fluxes in the presence of a mixed mode of compressional and transverse magnetic waves at hydromagnetic frequencies are investigated through kinetic perturbation of the distribution function. The magnetospheric medium where the particles are modulated contains both the magnetic and pressure gradients. The modulation features are found to be strongly dependent on the energy and pitch angle of the particles. Drifting particles can resonate with waves whose phase velocities are close to their drift velocities. When this occurs, the modulation amplitudes become significantly large and large phase shifts will occur. It is pointed out that resonance is important for particles with mid pitch angles (40-70 deg). The phase shift between the particle modulations and the magnetic field oscillations are strongly controlled by the combined effects of transverse and compressional wave components and/or the occurrence of drift resonance. Numerical calculations are performed using the dispersion relation of drift mirror Alfven waves as an example of waves with both compressional and transverse components.

Lin, C. S.↗

A search for spectroscopic binaries among the runaway O type stars

Numerous radial velocity measurements of medium dispersion were made for the 10 brighter stars given in Stone's list of very probable O type runaways. All plates were measured with the KPNO PDS microdensitometer, and a new iterative reductional analysis was used to derive plate velocities, which are estimated to be 1.6 times more accurate internally than those found by using the traditional method. Of thse stars, psi Per, alpha Cam, HD 188209, and 26 Cep are identified as probable velocity variables, while 9 Sge, lambda Cep, and HD 218915 are classed as possible variables. If the source of this variability is Keplerian rather than atmospheric, which cannot be established unequivocally from the observations of this paper, psi Per could be a spectroscopic binary with a black hole companion, and at least 1.2 solar mass. The detection of runaway binary systems from radial velocity measurements is discussed.

Stone, R. C.↗

Electromagnetic ion beam instabilities - Hot beams at interplanetary shocks

This paper considers Vlasov instabilities driven by a very hot ion beam streaming parallel to a magnetic field. The linear theory of electromagnetic instabilities driven by such a beam in a homogeneous, collisionless, nonrelativistic plasma is used. Numerical solutions of the full linear dispersion equation for bi-Maxwellian distribution functions are presented. If the thermal speed of the beam is greater than its drift speed, the two dominant modes are the right-hand and left-hand resonant ion beam instabilities. The parametric dependencies of the threshold drift speeds and growth rates for both modes are presented. In particular, the thresholds are shown to be sensitive functions of the beam temperature anisotropies. Arguments are presented that these two modes are driven by the suprathermal ion component at interplanetary shocks, and that the growth of these modes is sufficient to account for the amplitudes of MHD-like waves observed at such shocks.

Gary, S. P.↗

Electromagnetic electron temperature anisotropy instabilities

This paper considers electromagnetic Vlasov instabilities driven by electron temperature anisotropies in a homogeneous, nonrelativistic magnetized plasma. Numerical solutions of the full linear dispersion equation for bi-Maxwellian distribution functions and instabilities propagating parallel to the magnetic field are presented. Parametric dependences of the maximum growth rates of the electron fire hose and whistler anisotropy instabilities are given.

Gary, S. P.↗

The ion-ion acoustic instability

This paper considers the linear theory of ion-acoustic-like instabilities in a homogeneous Vlasov plasma with two ion components, a less dense beam and a more dense core, with a relative drift velocity. Numerical solutions of the full electrostatic dispersion equation are presented, and the properties of the ion-ion acoustic instability are studied in detail. At relatively cold beam temperatures, the instability is fluid-like, but it becomes a beam resonant kinetic instability as the beam temperature becomes of the order of the core temperature; if the mode is unstable, its threshold lies well below the threshold of the electron-ion acoustic instability. An electron temperature anisotropy enhances the instability, and at sufficiently large beam-core relative drift speeds, electron magnetization can either detract from or enhance the instability.

Gary, S. Peter↗

Magnetohydrodynamic waves and particle acceleration in a coronal loop

A model is presented for the acceleration of electrons in a flaring solar coronal loop. It is suggested that the primary energy release event in the flare establishes a spectrum of MHD waves within the loop which accelerate electrons by wave-particle resonant interactions. The nature of the MHD oscillations are described and numerical and analytic solution of the dispersion relation are examined. It is found that the behavior of a particle depends on whether or not it is in resonance with the wave spectrum. Resonant particles follow chaotic trajectories and can diffuse in velocity space between wave effective potentials, thereby gaining energy.

De La Beaujardiere, Jean-Francois↗

Theory of circular dielectric waveguide with anisotropic sheet cover

A circular dielectric waveguide consisting of an isotropic core covered with a thin anisotropic sheet is considered. The sheet is represented as a jump immittance and Maxwell's equations are applied. Solution of the boundary value problem yields the characteristic equation, or dispersion relation, which is then solved numerically. The results are verified for the step-index fiber and circular, metallic waveguides. Finally, examples are included to investigate the effects of the anisotropic sheet.

Kish, Jerome Thomas↗

Lagrangian and Eulerian statistics obtained from direct numerical simulations of homogeneous turbulence

Direct numerical simulation is used to study dispersion in decaying isotropic turbulence and homogeneous shear flow. Both Lagrangian and Eulerian data are presented allowing direct comparison, but at fairly low Reynolds number. The quantities presented include properties of the dispersion tensor, isoprobability contours of particle displacement, Lagrangian and Eulerian velocity autocorrelations and time scale ratios, and the eddy diffusivity tensor. The Lagrangian time microscale is found to be consistently larger than the Eulerian microscale, presumably due to the advection of the small scales by the large scales in the Eulerian reference frame.

Squires, Kyle D.↗

Particle dispersion in confined turbulent swirling flows

This paper reports a numerical investigation of confined swirling flows of gas-particle mixtures. A recently developed two-fluid, multiple-scale mixing model is applied to study the influence of particles on the intensity of the reverse flows of the gas phase and the effects of swirl on the particle dispersion in an annular expanding chamber under isothermal condition. The calculations were made for different swirl strength of injection of the annular jet into the mixing chamber. Results agree well qualitatively with experimental information available. It is also found that the calculated flow fields depend heavily on the prescription of the inlet flow conditions.

Chen, C. P.↗

Assessment of Total Variation Diminishing schemes in combustion instability problems

Numerical simulations of acoustic waves in a shear layer and in an idealized combustion chamber using high resolution Total Variation Diminishing (TVD) schemes has been carried out to study the effects of inherent scheme dissipation and dispersion errors of this class of problems. The numerical results are compared against available exact solutions to quantify these errors. Several popular TVD limiters widely used in the Computational Fluid Dynamics (CFD) community have been assessed. Osher-Chakravarthy limiters are modified so that they can be used in explicit schemes. Among all the limiters investigated, Osher-Chakravarthy third-order limiter is identified as having performed the best. It is also found that all TVD schemes have exceptionally small dispersive errors.

Wang, Z. J.↗

Higher order dispersion in the propagation of a gravity wave packet

To the first order of approximation, the complex amplitude of a wave packet in an anisotropic and dispersive medium is convected with the group of velocity. However, a gravity wave is a vector wave. Its wave packet must be formed by superposition of various wave numbers with corresponding frequencies, as is the case for scalar waves, and additionally by superposing many eigenmodes which also depend on the wave number. To represent the vector wave packet self-consistently, it is found that a gradient term must be included in the expansion. For a Guassian wave packet, this gradient term is shown to have important implications on the velocity vector as represented by its hodograph. Numerical results show that the hodograph is influenced by the location of the relative position of interest from the center of a Gaussian pulse. Higher order expansion shows that an initial Gaussian wave packet will retain its Gaussian shape as it propagates, but the pulse will spread in all directions with its major axis undergoing a rotation. Numerical results indicate that these higher order dispersive effects may be marginally observable in the atmosphere.

Yeh, K. C.↗

Computation of Transonic Nozzle Sound Transmission and Rotor Problems by the Dispersion-Relation-Preserving Scheme

The transonic nozzle transmission problem and the open rotor noise radiation problem are solved computationally. Both are multiple length scales problems. For efficient and accurate numerical simulation, the multiple-size-mesh multiple-time-step Dispersion-Relation-Preserving scheme is used to calculate the time periodic solution. To ensure an accurate solution, high quality numerical boundary conditions are also needed. For the nozzle problem, a set of nonhomogeneous, outflow boundary conditions are required. The nonhomogeneous boundary conditions not only generate the incoming sound waves but also, at the same time, allow the reflected acoustic waves and entropy waves, if present, to exit the computation domain without reflection. For the open rotor problem, there is an apparent singularity at the axis of rotation. An analytic extension approach is developed to provide a high quality axis boundary treatment.

Tam, Christopher K. W.↗

Finite difference methods for transient signal propagation in stratified dispersive media

Explicit difference equations are presented for the solution of a signal of arbitrary waveform propagating in an ohmic dielectric, a cold plasma, a Debye model dielectric, and a Lorentz model dielectric. These difference equations are derived from the governing time-dependent integro-differential equations for the electric fields by a finite difference method. A special difference equation is derived for the grid point at the boundary of two different media. Employing this difference equation, transient signal propagation in an inhomogeneous media can be solved provided that the medium is approximated in a step-wise fashion. The solutions are generated simply by marching on in time. It is concluded that while the classical transform methods will remain useful in certain cases, with the development of the finite difference methods described, an extensive class of problems of transient signal propagating in stratified dispersive media can be effectively solved by numerical methods.

Lam, D. H.↗

The whistler mode in a Vlasov plasma

In this study, properties of small-amplitude parallel and oblique whistler-mode waves are investigated for a wide range of plasma parameters by numerically solving the full electromagnetic Vlasov-dispersion equation. To investigate the cold-plasma and electrostatic approximations for the whistler mode, the results are compared with results obtained using these descriptions. For large wavelengths, the cold-plasma description is often accurate, while for short wavelengths and sufficiently oblique propagation, the electrostatic description is often accurate. The study demonstrates that in a Vlasov plasma the whistler mode near resonance has a group velocity more nearly parallel to the magnetic field than that predicted by cold-plasma theory.

Tokar, R. L.↗

A predictive guidance algorithm for Mars entry

This report presents the derivation, logic, and performance of a Mars entry guidance algorithm for the Mars Rover Sample Return Mission. This algorithm guides the Mars entry vehicle along an equilibrium glide trajectory designed to satisfy a set of equality and inequality constraints at chute deploy. An equilibrium glide assumption permits analytic prediction of terminal state conditions and eliminates the need for determining numeric state sensitivities. Trajectory deviations due to dispersions or uncertainties in Martian atmosphere, vehicle aerodynamics, or entry interface conditions are accommodated simply by predicting the compensating equilibrium glide profile and achieving this profile with the available lift control. This Mars Predictive Equilibrium Glide guidance scheme was tested on a range of vehicle configurations with a variety of aerodynamic, atmosphere, and trajectory dispersions with a high degree of success.

Tigges, M. A.↗

Coupling of Multiple Coulomb Scattering with Energy Loss and Straggling in HZETRN

The new version of the HZETRN deterministic transport code based on Green's function methods, and the incorporation of ground-based laboratory boundary conditions, has lead to the development of analytical and numerical procedures to include off-axis dispersion of primary ion beams due to small-angle multiple Coulomb scattering. In this paper we present the theoretical formulation and computational procedures to compute ion beam broadening and a methodology towards achieving a self-consistent approach to coupling multiple scattering interactions with ionization energy loss and straggling. Our initial benchmark case is a 60 MeV proton beam on muscle tissue, for which we can compare various attributes of beam broadening with Monte Carlo simulations reported in the open literature.

Mertens, Christopher J.↗

Flight Mechanics Modeling and Simulation of the Earth Entry System

Introduction: The Mars Sample Return (MSR) Campaign being planned by NASA and ESA has the ambitious goal to return Mars samples back to Earth. This international collaboration had developed a concept of operations that included a ESA-designed Earth Return Orbiter (ERO) and NASA-designed Capture, Containment, and Return System (CCRS). The Earth Entry System (EES), consisting of a protective aeroshell that houses the samples as well as sample containment vessels, would conduct entry, descent, and landing (EDL) on a direct Earth trajectory. The EES would enter on a spin-stabilized ballistic trajectory with the goal to passively achieve aerodynamic stability throughout all regions of flight. The EDL sequence would end with the EES impacting the soft playa soil of the Utah Test and Training Range (UTTR). As of the submission of this abstract, the MSR campaign is undergoing a re-architecture leading to a pause in EES development. However, the novel approaches developed in flight mechanics modeling and simulation can significantly benefit the greater IPPW community in the development of Earth return vehicles. This paper will present the latest state of EES flight mechanics modeling and simulation. The paper will highlight the simulation architecture developed and key lessons learned from understanding of EDL trajectory sensitivities. Modeling and Simulation: Figure 1 provides a high-level concept of operations for the approach, entry, descent, and landing (AEDL) phase of the CCRS-portion of MSR. The objective of EES flight mechanics is to model and simulate the EES trajectory from ERO separation to ground impact at UTTR. A variety of flight mechanics simulation models were utilized to model both exo-atmopsheric and atmospheric portions of flight. 42, a 6-DOF simulation developed at Goddard Space Flight Center, is utilized for propagating the attitude of EES during exo-atmospheric flight. 42 allows for a variety of spin eject mechanism scenarios to be simulated for analysis. 10 minutes prior to entry, the 42 states are handed off to the EDL sims. The prime EDL sim utilized by EES is the Program to Optimize Simulated Trajectories II (POST2), a 6-DOF sim developed at Langley Research Center, and the independent verification and validation EDL sim utilized is DSENDS, a 6-DOF sim developed at Jet Propulsion Laboratory. Figure 2 provides a visualization of the flight mechanics simulation model flow through various points in the AEDL phase. Due to the existence of a variety of sim models, the EES flight mechanics team developed processes for data hand-off. These processes included the development of a centralized coordinate frame document, utilization of a single, centralized simulation input document for all sims to reference, and hand-off files containing both the technical data to be ingested by other flight mechanics sims as well as annotations of modeling assumptions utilized to generate the data. Figure~\ref{fig:post2simarchitecture} provides an overview of the POST2 sim architecture wherein POST2 ingests numerous subsystem models and input files. The dispersed state file generated by MONTE provides the position/velocity state of the trajectory while the 42 Handoff file provides the attitude. The aerodynamics database, delivered by the EES aeroscience team, is utilized to simulate the aerodynamic forces and moments experienced during EDL. A custom atmosphere model, developed by EES atmosphere team, is utilized to simulate the anticipated atmosphere environment around the region of Earth through which the EES trajectory flys. These inputs and subsystem models can be varied depending on the AEDL flight mechanics scenario being simulated. Monte Carlo simulations are utilized to generate statistical AEDL performance metrics in the form of scorecards and violin plots. Furthermore, outputs from the POST2 simulation are utilized for follow-on analyses including aerothermal and landing performance. \section{Flight Mechanics Lessons Learned} Though the EES flight mechanics team uncovered a variety of lessons learned through the analysis conducted to support CCRS through preliminary design review, this paper will highlight the most important lessons. A key AEDL performance goal is to ensure the landing footprint of EES remains on the UTTR south range. A common modeling strategy used in EDL analysis is One-Variable-At-a-Time (OVAT). OVAT analysis provides insight into the key drivers that affect AEDL performance metrics. Figure 3 shows the landing ellipses for single dispersion sources as compared to the baseline aggregate of all dispersions. The figure shows that atmosphere winds alone dominate the size of the footprint ellipse (note: EES does not use a parachute unlike previous Earth-return missions and is in wind-driven free fall for ~5min). The significance of the wind led the EES flight mechanics team to pursue the development of a Custom Atmosphere Model [4], in lieu of EarthGRAM [1], built on actual radiosonde wind measurements around the UTTR-region. This decision was driven by the realism in the generated footprint ellipses and lessons-learned from Stardust [5]. These findings will be invaluable for future Earth-return missions in providing an early understanding of the key drivers affecting footprint size and modeling considerations for which to account. Another lesson learned is tied to the AEDL performance goal of achieving passive stability throughout all regions of flight. It is well understood that blunt-body aeroshells are less stable as they transition from supersonic to subsonic. Eliminating a backshell does help improvestability; however, other phenomena such as roll-induced instability during terminal descent can still arise. The EES flight mechanics team developed stability metrics as tools to better understand the causes of and better predict the onset of dynamic instability. These tools were built upon analytical models developed by Jaffe [3] and Murphy [2]. The tools were shown to both be very accurate in correlation with actual unstable cases and useful in developing stability margin policies based on the vehicle design and simulation considerations (e.g. sphere-cone angle change, mass change, wind turbulence). These tools allowed for the current EES design to demonstrate the ability to achieve passive stability and can be an invaluable tool for consideration in the design of parachute-less Earth-return vehicles.

Rohan Deshmukh↗

Velocity dispersions in galaxies: 1: The SO galaxy NGC 7332

A Coude spectrum of the SO galaxy NGC 7332 with 0.9 A resolution from 4186 to 4364 A was obtained with the SEC vidicon television camera and the Hale telescope. Comparisons with spectra of G and K giant stars, numerically broadened for various Maxwellian velocity distributions, give a dispersion velocity in the line of sight of 160 + or - 20 km/sec with the best fit at G8III. The dispersion appears to be constant within + or - 35 km/sec out to 1.4 kpc (H = 100 km/sec/mpc). After correction for projection, the rotation curve has a slope of 0.16 km/sec/pc at the center and a velocity of 130 km/sec at 1.4 kpc where it is still increasing. For an estimated effective radius of 3.5 kpc enclosing half the light, the virial theorem gives a mass of 1.4 x 10 to the 11th power solar masses if the mass-to-light ratio is constant throughout the galaxy. The photographic luminosity is 8.3 x 10 to the 9th power solar luminosities so that the M/L ratio is 17.

Morton, D. C.↗