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 361 records · Page 20

Catalog of far-ultraviolet objective-prism spectrophotometry: Skylab experiment S-019, ultraviolet steller astronomy

Ultraviolet stellar spectra in the wavelength region from 1300 to 5000 A (130 to 500) were photographed during the three manned Skylab missions using a 15 cm aperture objective-prism telescope. The prismatic dispersion varied from 58 A mm/1 at 1400 A to 1600 A mm/1 at 3000 A. Approximately 1000 spectra representing 500 stars were measured and reduced to observed fluxes. About 100 stars show absorption lines of Si IV, C IV, or C II. Numerous line features are also recorded in supergiant stars, shell stars, A and F stars, and Wolf-Rayet stars. Most of the stars in the catalog are of spectral class B, with a number of O and A type stars and a sampling of WC, WN, F and C type stars. Spectrophotometric results are tabulated for these 500 stars.

Henize, K. G.↗

Two dimensional nonlinear analysis of sound transmission through a near-sonic throat flow

A two-dimensional nonlinear theory of sound transmission through a nonuniform duct carrying a near-sonic throat flow is described. The mean flow in the duct is treated according to a generalized quasi-one dimensional model, and the analysis of the unsteady perturbations is carried out using the Method of Matched Asymptotic Expansions. The linearized acoustic field in the subsonic regions of the duct is approximated by a Wave Envelope model which is matched asymptotically with a nonlinear inner solution valid in the near-sonic throat region. Numerical results are presented to illustrate the predictions of the theory. It is found that, in general, shock waves develop in the acoustic field during transmission through the near-sonic flow. Unlike previous one-dimensional theories, the current study shows that dispersion can play a major role in the propagation process.

Myers, M. K.↗

Large-scale structure after COBE: Peculiar velocities and correlations of cold dark matter halos

Large N-body simulations on parallel supercomputers allow one to simultaneously investigate large-scale structure and the formation of galactic halos with unprecedented resolution. Our study shows that the masses as well as the spatial distribution of halos on scales of tens of megaparsecs in a cold dark matter (CDM) universe with the spectrum normalized to the anisotropies detected by Cosmic Background Explorer (COBE) is compatible with the observations. We also show that the average value of the relative pairwise velocity dispersion sigma(sub v) - used as a principal argument against COBE-normalized CDM models-is significantly lower for halos than for individual particles. When the observational methods of extracting sigma(sub v) are applied to the redshift catalogs obtained from the numerical experiments, estimates differ significantly between different observation-sized samples and overlap observational estimates obtained following the same procedure.

Zurek, Wojciech H.↗

Development of a Numeric Predictor-Corrector Aerocapture Guidance for Direct Force Control

Direct force control, where the angle of attack and sideslip angle are modulated, has been proposed as an alternative to bank angle control for aero-assist maneuvers. This paper reimplements the current state-of-the-art aerocapture guidance for bank angle control, Fully Numeric Predictor-corrector Aerocapture Guidance (FNPAG), for direct force control. The optimal control theory underlying the structure of FNPAG is shown to not be applicable to the direct force control approach. Several solution structures for the longitudinal channel are compared by simulating dispersed three-degree-of-freedom trajectories for a reference mission consisting of a low lift-to-drag vehicle and a highly elliptical, 1-sol target orbit around Mars. The equations of motion for the lateral channel are derived, and a controller is designed to target a specified orbital plane. Finally, a Monte Carlo is used to demonstrate the performance of the new guidance.

Matz, Daniel A.↗

Velocities and rotational excitation of interstellar H2 toward Pi Scorpii

A spectrum of Pi Sco showing numerous atomic lines and 70 absorption features from the Lyman and Werner transitions of interstellar H2 in rotational level J from zero to five is presented. Their shapes of the composite column density profiles are very nearly Gaussian with a one-dimensional rms velocity dispersion of 3 km/s. The behavior of shifts in the inferred N(H2) as a function of velocity are consistent with the overall profiles being composed of nearly symmetrical, tightly paced assemblies of about seven unresolved components. The relative overall column densities in the higher J levels of H2 are consistent with a model where these states are populated by optical pumping through the Lyman and Werner transitions, powered by UV radiation from nearby stars. The slight narrowing of the high-J profiles may be due to small clumps of H2 at radial velocities some 5-8 km/s from the core of the profile are exposed to a pumping flux about 10 times lower than that for the material near the profile's center.

Jenkins, Edward B.↗

Dispersion-relation-preserving finite difference schemes for computational acoustics

Time-marching dispersion-relation-preserving (DRP) schemes can be constructed by optimizing the finite difference approximations of the space and time derivatives in wave number and frequency space. A set of radiation and outflow boundary conditions compatible with the DRP schemes is constructed, and a sequence of numerical simulations is conducted to test the effectiveness of the DRP schemes and the radiation and outflow boundary conditions. Close agreement with the exact solutions is obtained.

Tam, Christopher K. W.↗

Modeling the Complete Planetary Subsurface Radio Remote Sensing Problem

Ground penetrating radar (GPR) methods are in the forefront of the search for subsurface water on Mars largely because of their platform versatility. In theory these instruments can operate from orbit, on the ground, and any altitude in between. Recent review papers and white papers make the strong case for using GPR methods in subsurface water exploration. But important questions have been raised recently about effect of unknown parameters on the ability to successfully obtain planetary GPR measurements. The almost completely unknown lower ionosphere on Mars may cause serious signal absorption problems for GPR instruments on orbiting platforms, and lossy upper layers of the ground may reduce the returned signal for any platform. What is clearly needed is a minimum approximations, full-wave model of the complete GPR problem, including ionospheric dispersion and absorption, surface transmission, and subsurface scattering. We are developing such a model that is as general purpose as possible, allowing arbitrary ionospheric parameters, surface roughness, and subsurface inhomogeneities. We present some of the details of this model, and highlight some of its capabilities with numerical examples. Additional information is contained in the original extended abstract.

Cummer, S. A.↗

Scalable High Performance Computing: Direct and Large-Eddy Turbulent Flow Simulations Using Massively Parallel Computers

This final report contains reports of research related to the tasks "Scalable High Performance Computing: Direct and Lark-Eddy Turbulent FLow Simulations Using Massively Parallel Computers" and "Devleop High-Performance Time-Domain Computational Electromagnetics Capability for RCS Prediction, Wave Propagation in Dispersive Media, and Dual-Use Applications. The discussion of Scalable High Performance Computing reports on three objectives: validate, access scalability, and apply two parallel flow solvers for three-dimensional Navier-Stokes flows; develop and validate a high-order parallel solver for Direct Numerical Simulations (DNS) and Large Eddy Simulation (LES) problems; and Investigate and develop a high-order Reynolds averaged Navier-Stokes turbulence model. The discussion of High-Performance Time-Domain Computational Electromagnetics reports on five objectives: enhancement of an electromagnetics code (CHARGE) to be able to effectively model antenna problems; utilize lessons learned in high-order/spectral solution of swirling 3D jets to apply to solving electromagnetics project; transition a high-order fluids code, FDL3DI, to be able to solve Maxwell's Equations using compact-differencing; develop and demonstrate improved radiation absorbing boundary conditions for high-order CEM; and extend high-order CEM solver to address variable material properties. The report also contains a review of work done by the systems engineer.

Morgan, Philip E.↗

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Olsson & Oliger, Gerritsen & Olsson, and Yee et al. (15, 7, 37) is revisited. The objective of this paper is to prove for the first time that the entropy split scheme is an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for nonlinear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. The proof is to replace the spatial derivatives by summation-by-parts (SBP) difference operators in the entropy split form of the equations using the physical entropy of the Euler equations. The numerical boundary closure follows directly from the SBP operator. No additional numerical boundary procedure is required. In contrast, Tadmor-type entropy conserving schemes (31) using mathematical entropies and more recently in (35], do not naturally come with a numerical boundary closure and a generalized SBP operator has to be developed (18). Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grubber splitting.

Sjogreen, Bjorn↗

Finite beta effects on the drift-cyclotron instability

The effects of finite plasma beta and strong plasma inhomogeneities on the microinstability resulting from the unstable coupling of ion-cyclotron oscillations and an ion drift wave are examined. A local electromagnetic, kinetic dispersion equation, which takes into account finite orbit modifications of both ions and electrons as well as the effects of strong inhomogeneities on the ions and permits the continuous tracking of the instability from the weak gradient regime to the strong gradient regime, is derived and solved numerically. Finite beta is shown to increase the frequency and reduce the growth rate of the instability, but not completely stabilize it due to electromagnetic effects. When plasma inhomogeneity is sufficiently strong, finite beta reduces the growth rate of the large k sub y band of a given harmonic more than the small k sub y band; however, the most unstable bands remain the same as in the case of zero beta. Implications of these results for the Tandem Mirror Experiment are discussed.

Gladd, N. T.↗

A diagnostic analysis of the VVP single-doppler retrieval technique

A diagnostic analysis of the VVP (volume velocity processing) retrieval method is presented, with emphasis on understanding the technique as a linear, multivariate regression. Similarities and differences to the velocity-azimuth display and extended velocity-azimuth display retrieval techniques are discussed, using this framework. Conventional regression diagnostics are then employed to quantitatively determine situations in which the VVP technique is likely to fail. An algorithm for preparation and analysis of a robust VVP retrieval is developed and applied to synthetic and actual datasets with high temporal and spatial resolution. A fundamental (but quantifiable) limitation to some forms of VVP analysis is inadequate sampling dispersion in the n space of the multivariate regression, manifest as a collinearity between the basis functions of some fitted parameters. Such collinearity may be present either in the definition of these basis functions or in their realization in a given sampling configuration. This nonorthogonality may cause numerical instability, variance inflation (decrease in robustness), and increased sensitivity to bias from neglected wind components. It is shown that these effects prevent the application of VVP to small azimuthal sectors of data. The behavior of the VVP regression is further diagnosed over a wide range of sampling constraints, and reasonable sector limits are established.

Boccippio, Dennis J.↗

New Releases and Upgrades of Giant Planet Global Reference Atmospheric Models (GRAMs)

• GRAMs are engineering-oriented atmospheric models that estimate mean values and statistical variations of atmospheric properties for numerous planetary destinations – Outputs atmospheric density, temperature, pressure, winds, and chemical composition along a user-defined path – Provide mean values and variability for any point in an atmosphere – Includes seasonal, geographic, and altitude variations – Used by engineering community because of their ability to create realistic atmospheric dispersions; can be integrated into high fidelity flight dynamic simulations of launch, entry, descent, and landing (EDL), aerobraking, and aerocapture • GRAMs are not forecast models • GRAMs are available for: Earth, Mars, Venus, Neptune, Titan, Jupiter, and Uranus • Available through the NASA Software Catalog https://software.nasa.gov/software/MFS-33888-1

atmospheric models↗

Axisymmetric disturbances in a disk galaxy

The Lin-Shu dispersion relation, with corrections to match experimental conditions and in combination with a quantization rule, permits growth rates to be calculated for axisymmetric disturbances in an axisymmetric disk of stars. Growth rates calculated this way are in reasonable quantitative agreement with recently published results from a series of numerical experiments. The agreement confirms the physical basis of the theory and demonstrates that the theory remains valid even when strained somewhat in the application to galactic models.

Miller, R. H.↗

Lunar sample 14425 - Corrected analysis

A polished section of lunar sample 14425, an 8 mm glass bead, was studied using a scanning electron microscope and energy-dispersive X-ray analyzer. The silicon, aluminum, iron, magnesium, calcium, sodium, potassium, titanium, and chromium contents of the glass were determined. Two types of glass are visible in the polished section. One is clear and almost devoid of metallic spherules, while the other is cloudy and contains numerous metallic spherules, some less than one micron in size. Both glass types are homogeneous and identical in composition. This composition closely matches that of some Apollo 14 breccias or glass found in the breccias. The apparent similarity in composition between lunar sample 14425 and the high-magnesium microtektites found in a previous study was probably due to charging effects during analyses in which the sample was uncoated.

Glass, B. P.↗

A numerical simulation of galaxy subcluster mergers

We present preliminary results of a 3-D numerical simulation of two merging subclusters of galaxies. By self-consistently modelling the intracluster gas and dark matter dynamics, we hope to gain insight as to how the dynamics of both relate to such observables as the cluster x-ray emission, radio source morphology, and velocity dispersions.

Roettiger, Kurt↗

Catching Constrains on the Parent Body Genesis of Mesosiderites and a Possible Link to HED (Howardite-Eucrite-Diogenite) Meteorites - A New Hope?

Mesosiderites (MES) are a group of enigmatic stony-iron meteorites exhibiting fragmental matrix breccias and irregular textures; e.g. [1-3]. Mesosiderites contain roughly equal volumes metal (Fe-Ni) and silicates often intimately mixed together (Fig.1). The silicates mostly consist of basaltic, gabbroic, and pyroxenitic components, and appear similar to eucrites and howardites; [4-8]. But unlike HEDs - and other differentiated parent body meteorite groups e.g. ureilites - mesosiderites contain high metal abundances. Several studies have been published to reveal the processes leading to the formation of mesosiderites and attempt to classifiy them [1], [2], [10-15]. Because the silicate inclusions in mesosiderites are often strongly metamorphosed after formation, it is difficult to assess the origin of the silicates and implications for the differentiation process of their parent body [15-17]. Several workers have advanced a formation hypothesis for the mesosiderites where an impact between differentiated bodies occurred prior to 4.47 Ga ago (e.g. [13,18], which could explain the possible incomplete dispersal of the colliding bodies due to their low cosmic ray exposure ages and their special thermal history. However, [13] discuss and favor the model for formation of mesosiderites with the collision of two differentiated bodies, along with disruption events and gravitational re-assembly. The mesosiderites have numerous gabbroid melt clasts with anomalous rare-earth- element (REE) - especially positive Eu - values [19, 20]. HEDs do not show the same. However, the heating mechanisms of both mesosiderites and HED's are puzzling. Mesosiderites are remarkable, they consist of a mix of basalts, which are only found on or near planetary surfaces and undifferentiated metal [1,2]. The probable model is that an asteroid containing a metallic magma impacted onto a second asteroid covered with basalt [18,21]. The mix was then buried under an insulating regolith, and cooled slowly. During cooling and at low temperatures the redox reactions continued to occur and proceed (J.T. Wasson; in pers. comm. 2015).

Eucrite↗

On A Stabilization of the Ingard-Myers Impedance Boundary Condition

It has been well-known that the Ingard-Myers impedance condition, while simple to apply, is subject to the hydrodynamic Kelvin-Helmholtz instability due to its use of a vortex sheet in modeling the flow at the liner boundary. Recently, in the development of a time domain boundary element method for acoustic scattering by treated surfaces, it was found that by neglecting a certain second-order spatial derivative term in the Ingard-Myers formulation, the hydrodynamic instability can be avoided. The present paper aims to provide further analysis of this modified condition, hereby referred to as Truncated Ingard-Myers Impedance Boundary Condition (TIMIBC). It will be shown, based on the dispersion relations of linear waves, that the instability intrinsic to the Ingard-Myers condition is eliminated in the proposed new formulation. Quantitative assessments on the accuracy of TIMIBC for scattering of acoustic waves by lined surfaces will be carried out, and its effectiveness will be demonstrated by numerical examples. Specifically, the accuracy is assessed by comparing solutions obtained by the Ingard-Myers condition with that by the proposed TIMIBC where theoretical reflection coefficients at a lined surface are derived for cases of plane and spherical incident waves. It is found that the TIMIBC provides a good approximation to the original Ingard-Myers condition for flows of low to mid subsonic Mach numbers. Time domain implementations of TIMIBC are also discussed and illustrated with a numerical example using a finite difference scheme. As many studies have shown the Ingard-Myers condition to be the correct limit of boundary layer thickness going to zero, the proposed TIMIBC can offer a practical solution for overcoming the intrinsic instability associated with the Ingard-Myers condition.

Time Domain BEM↗

On A Stabilization of the Ingard-Myers Impedance Boundary Condition

It has been well-known that the Ingard-Myers impedance condition, while simple to apply, is subject to the hydrodynamic Kelvin-Helmholtz instability due to its use of a vortex sheet in modeling the flow at the liner boundary. Recently, in the development of a time domain boundary element method for acoustic scattering by treated surfaces, it was found that by neglecting a certain second-order spatial derivative term in the Ingard-Myers formulation, the hydrodynamic instability can be avoided. The present paper aims to provide further analysis of this modified condition, hereby referred to as Truncated Ingard-Myers Impedance Boundary Condition (TIMIBC). It will be shown, based on the dispersion relations of linear waves, that the instability intrinsic to the Ingard-Myers condition is eliminated in the proposed new formulation. Quantitative assessments on the accuracy of TIMIBC for scattering of acoustic waves by lined surfaces will be carried out, and its effectiveness will be demonstrated by numerical examples. Specifically, the accuracy is assessed by comparing solutions obtained by the Ingard-Myers condition with that by the proposed TIMIBC where theoretical reflection coefficients at a lined surface are derived for cases of plane and spherical incident waves. It is found that the TIMIBC provides a good approximation to the original Ingard-Myers condition for flows of low to mid subsonic Mach numbers. Time domain implementations of TIMIBC are also discussed and illustrated with a numerical example using a finite difference scheme. As many studies have shown the Ingard-Myers condition to be the correct limit of boundary layer thickness going to zero, the proposed TIMIBC can offer a practical solution for overcoming the intrinsic instability associated with the Ingard-Myers condition.

Time Domain BEM↗