Search NASA⌕ Search

SEARCH · Search NASA

Results for “wave speed”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Computation of viscous blast wave solutions with an upwind finite volume method

A fully conservative, viscous, implicit, upwind, finite-volume scheme for the thin-layer Navier-Stokes equations is described with application to blast wave flow fields. In this scheme, shocks are captured without the oscillations typical of central differencing techniques and wave speeds are accurately predicted. The finite volume philosophy ensures conservation and since boundary conditions are also treated conservatively, accurate reflections of waves from surfaces are assured. Viscous terms in the governing equations are treated in a manner consistent with the finite volume philosophy, resulting in very accurate prediction of boundary layer quantities. Numerical results are presented for four viscous problems: a steady boundary layer, a shock-induced boundary layer, a blast wave/cylinder interaction and a blast wave/supersonic missile interaction. Comparisons of the results with an established boundary layer code, similarity solution, and experimental data show excellent agreement.

Molvik, Gregory A.↗
NASA NTRS1991

Computation of viscous blast wave flowfields

A method to determine unsteady solutions of the Navier-Stokes equations was developed and applied. The structural finite-volume, approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the interaction of blast-waves with stationary targets. The inviscid flux is evaluated using MacCormack's modified Steger-Warming flux or Roe flux difference splittings with total variation diminishing limiters, while the viscous flux is computed using central differences. The use of implicit boundary conditions in conjunction with a telescoping in time and space method permitted solutions to this strongly unsteady class of problems. Comparisons of numerical, analytical, and experimental results were made in two and three dimensions. These comparisons revealed accurate wave speed resolution with nonoscillatory discontinuity capturing. The purpose of this effort was to address the three-dimensional, viscous blast-wave problem. Test cases were undertaken to reveal these methods' weaknesses in three regimes: (1) viscous-dominated flow; (2) complex unsteady flow; and (3) three-dimensional flow. Comparisons of these computations to analytic and experimental results provided initial validation of the resultant code. Addition details on the numerical method and on the validation can be found in the appendix. Presently, the code is capable of single zone computations with selection of any permutation of solid wall or flow-through boundaries.

Atwood, Christopher A.↗
NASA NTRS1986

Analysis of the Giacobini-Zinner bow wave

The cometary bow wave of P/Giacobini-Zinner has been analyzed using the complete set of ICE field and particle observations to determine if it is a shock. Changes in the magnetic field and plasma flow velocities from upstream to downstream have been analyzed to determine the direction of the normal and the propagation velocity of the bow wave. The velocity has then been compared with the fast magnetosonic wave speed upstream to derive the Mach number and establish whether it is supersonic, i.e., a shock, or subsonic, i.e., a large amplitude wave. The various measurements have also been compared with values derived from a Rankine-Hugoniot analysis. The results indicate that, inbound, the bow wave is a shock with M = 1.5. Outbound, a subsonic Mach number is obtained, however, arguments are presented that the bow wave is also likely to be a shock at this location.

Smith, E. J.↗
NASA NTRS1986

Analysis of the Giacobini-Zinner bow wave

The cometary bow wave of P/Giacobini-Zinner was analyzed using the ICE field and particle observations to determine if it is a shock. Changes in the magnetic field and plasma flow velocities from upstream to downstream were analyzed to determine the direction of the normal and the propagation velocity of the bow wave. The velocity was compared with the fast magnetosonic wave speed upstream to derive the Mach number and establish whether it is supersonic, i.e., a shock, or subsonic, i.e., a large amplitude wave. The measurements were compared with values derived from a Rankine-Hugoniot analysis. The results indicate that inbound the bow wave is a shock with M = 1.5. Outbound, a subsonic Mach number is obtained. However, arguments that the bow wave is also likely to be a shock at this location are presented.

Smith, E. J.↗
NASA NTRS1989

A procedure for computing surface wave trajectories on an inhomogeneous surface

Equations are derived for computing surface waves on smooth surfaces, including surfaces with a nonuniform wave speed. The prior literature dealt primarily with the theoretical development with little consideration given to computational methods, and examples were limited to waves on surfaces of simple analytic description, such as cones, spheres, and cylinders. The computational procedure presented is a relatively general method. Sample calculations illustrate the procedure for a class of practical shapes of the type that include aerodynamic and hydrodynamic surfaces. Equations are also included for computing the spreading of rays into a surrounding medium that will support waves.

Barger, Raymond L.↗
NASA NTRS1978

Analytical study of laser-supported combustion waves in hydrogen

Laser supported combustion (LSC) waves are an important ingredient in the fluid mechanics of CW laser propulsion using a hydrogen propellant and 10.6 micron lasers. Therefore, a computer model has been constructed to solve the one-dimensional energy equation with constant pressure and area. Physical processes considered include convection, conduction, absorption of laser energy, radiation energy loss, and accurate properties of equilibrium hydrogen. Calculations for 1, 3, 10 and 30 atm were made for intensities of 10 to the 4th to 10 to the 6th W/sq cm, which gave temperature profiles, wave speed, etc. To pursue the propulsion application, a second computer model was developed to describe the acceleration of the gas emerging from the LSC wave into a variable-pressure, converging streamtube, still including all the above-mentioned physical processes. The results show very high temperatures in LSC waves which absorb all the laser energy, and high radiative losses.

Kemp, N. H.↗
NASA NTRS1992

Navier-Stokes simulations of unsteady transonic flow phenomena

Numerical simulations of two classes of unsteady flows are obtained via the Navier-Stokes equations: a blast-wave/target interaction problem class and a transonic cavity flow problem class. The method developed for the viscous blast-wave/target interaction problem assumes a laminar, perfect gas implemented in a structured finite-volume framework. The approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the blast-waves with stationary targets. The inviscid flux is evaluated using either of two upwind techniques, while the full viscous terms are computed by central differencing. Comparisons of unsteady numerical, analytical, and experimental results are made in two- and three-dimensions for Couette flows, a starting shock-tunnel, and a shock-tube blockage study. The results show accurate wave speed resolution and nonoscillatory discontinuity capturing of the predominantly inviscid flows. Viscous effects were increasingly significant at large post-interaction times. While the blast-wave/target interaction problem benefits from high-resolution methods applied to the Euler terms, the transonic cavity flow problem requires the use of an efficient scheme implemented in a geometrically flexible overset mesh environment. Hence, the Reynolds averaged Navier-Stokes equations implemented in a diagonal form are applied to the cavity flow class of problems. Comparisons between numerical and experimental results are made in two-dimensions for free shear layers and both rectangular and quieted cavities, and in three-dimensions for Stratospheric Observatory For Infrared Astronomy (SOFIA) geometries. The acoustic behavior of the rectangular and three-dimensional cavity flows compare well with experiment in terms of frequency, magnitude, and quieting trends. However, there is a more rapid decrease in computed acoustic energy with frequency than observed experimentally owing to numerical dissipation. In addition, optical phase distortion due to the time-varying density field is modelled using geometrical constructs. The computed optical distortion trends compare with the experimentally inferred result, but underpredicts the fluctuating phase difference magnitude.

Atwood, C. A.↗
NASA NTRS1977

Nonlinear stress waves in a perfectly flexible string

This paper discusses nonlinear stress-wave propagation in a perfectly flexible string obeying a quasilinear (rate-dependent) constitutive equation. Wave speeds and compatibility relations valid along various families of characteristics were determined. It was shown that the compatibility relations associated with the transverse as well as the longitudinal waves readily yield a physical interpretation when they are expressed in suitable variables and in vector form. Coding based on the present information was completed for the machine solution of a class of mixed initial- and boundary-value problems of practical interest. Computer simulation of the stress-wave interaction in the 40-foot lanyard in the Arcas 'Rocoz' system during deployment was carried out using a stress-strain relation for nylon at the strain rate of 30/second. A method for estimating the maximum tension and strain in a string during the initial loading phase is proposed.

Fan, D.-N.↗
NASA NTRS1996

Nonstationary Gravity Wave Forcing of the Stratospheric Zonal Mean Wind

The role of gravity wave forcing in the zonal mean circulation of the stratosphere is discussed. Starting from some very simple assumptions about the momentum flux spectrum of nonstationary (non-zero phase speed) waves at forcing levels in the troposphere, a linear model is used to calculate wave propagation through climatological zonal mean winds at solstice seasons. As the wave amplitudes exceed their stable limits, a saturation criterion is imposed to account for nonlinear wave breakdown effects, and the resulting vertical gradient in the wave momentum flux is then used to estimate the mean flow forcing per unit mass. Evidence from global, assimilated data sets are used to constrain these forcing estimates. The results suggest the gravity-wave-driven force is accelerative (has the same sign as the mean wind) throughout most of the stratosphere above 20 km. The sense of the gravity wave forcing in the stratosphere is thus opposite to that in the mesosphere, where gravity wave drag is widely believed to play a principal role in decelerating the mesospheric jets. The forcing estimates are further compared to existing gravity wave parameterizations for the same climatological zonal mean conditions. Substantial disagreement is evident in the stratosphere, and we discuss the reasons for the disagreement. The results suggest limits on typical gravity wave amplitudes near source levels in the troposphere at solstice seasons. The gravity wave forcing in the stratosphere appears to have a substantial effect on lower stratospheric temperatures during southern hemisphere summer and thus may be relevant to climate.

Alexander, M. J.↗
NASA NTRS1992

Film condensation in a horizontal rectangular duct

Condensation heat transfer in an annular flow regime with and without interfacial waves was experimentally investigated. The study included measurements of heat transfer rate with condensation of vapor flowing inside a horizontal rectangular duct and experiments on the initiation of interfacial waves in condensation, and adiabatic air-liquid flow. An analytical model for the condensation was developed to predict condensate film thickness and heat transfer coefficients. Some conclusions drawn from the study are that the condensate film thickness was very thin (less than 0.6 mm). The average heat transfer coefficient increased with increasing the inlet vapor velocity. The local heat transfer coefficient decreased with the axial distance of the condensing surface, with the largest change at the leading edge of the test section. The interfacial shear stress, which consisted of the momentum shear stress and the adiabatic shear stress, appeared to have a significant effect on the heat transfer coefficients. In the experiment, the condensate flow along the condensing surface experienced a smooth flow, a two-dimensional wavy flow, and a three-dimensional wavy flow. In the condensation experiment, the local wave length decreased with the axial distance of the condensing surface and the average wave length decreased with increasing inlet vapor velocity, while the wave speed increased with increasing vapor velocity. The heat transfer measurements are reliable. And, the ultrasonic technique was effective for measuring the condensate film thickness when the surface was smooth or had waves of small amplitude.

Lu, Qing↗
NASA NTRS1991

Transport of gyration-dominated space plasmas of thermal origin. I - Generalized transport equations

Grad's 20-moment set of transport equations has been examined in the limit of strong external magnetic fields. This simplified set of equations describes the transport of mass and parallel momentum as well as the transport of parallel and perpendicular energy and heat flow in the magnetic field direction. The effect of collisions was calculated assuming a modified relaxation model. Wave speeds and normal modes of the simplified set of equations were examined for an ion and electron gas. It was found that four of the ten normal modes are electron thermal-heat waves which approximately decouple from the six ion waves in the system. When low-frequency waves are considered (slow-wave approximation), this allows the electron energy and heat flow equations to be solved separately from the ion equations and in a time-independent fashion. When this was done, it was found that, under certain conditions, these equations predict an infinite electron perpendicular temperature, Te-perpendicular, in the collisionless regime. This occurs whenever Te-perpendicular is greater than the parallel temperature, at any point along collisionless and diverging magnetic field lines. When applied appropriately, the significantly simplified set of equations derived here are well suited for application to a large variety of problems in planetary ionospheres and magnetospheres.

Gombosi, Tamas I.↗
NASA NTRS1995

Molecular dynamics simulation of a piston driven shock wave in a hard sphere gas

Molecular dynamics simulation is used to study the piston driven shock wave at Mach 1.5, 3, and 10. A shock tube, whose shape is a circular cylinder, is filled with hard sphere molecules having a Maxwellian thermal velocity distribution and zero mean velocity. The piston moves and a shock wave is generated. All collisions are specular, including those between the molecules and the computational boundaries, so that the shock development is entirely causal, with no imposed statistics. The structure of the generated shock is examined in detail, and the wave speed; profiles of density, velocity, and temperature; and shock thickness are determined. The results are compared with published results of other methods, especially the direct simulation Monte-Carlo method. Property profiles are similar to those generated by direct simulation Monte-Carlo method. The shock wave thicknesses are smaller than the direct simulation Monte-Carlo results, but larger than those of the other methods. Simulation of a shock wave, which is one-dimensional, is a severe test of the molecular dynamics method, which is always three-dimensional. A major challenge of the thesis is to examine the capability of the molecular dynamics methods by choosing a difficult task.

Woo, Myeung-Jouh↗
NASA NTRS2022

How a Realistic Magnetosphere Alters the Polarizations of Surface, Fast Magnetosonic, and Alfvén Waves

System-scale magnetohydrodynamic (MHD) waves within Earth's magnetosphere are often understood theoretically using box models. While these have been highly instructive in understanding many fundamental features of the various wave modes present, they neglect the complexities of geospace such as the inhomogeneities and curvilinear geometries present. Here, we show global MHD simulations of resonant waves impulsively excited by a solar wind pressure pulse. Although many aspects of the surface, fast magnetosonic (cavity/waveguide), and Alfvén modes present agree with the box and axially symmetric dipole models, we find some predictions for large-scale waves are significantly altered in a realistic magnetosphere. The radial ordering of fast mode turning points and Alfvén resonant locations may be reversed even with monotonic wave speeds. Additional nodes along field lines that are not present in the displacement/velocity occur in both the perpendicular and compressional components of the magnetic field. Close to the magnetopause, the perpendicular oscillations of the magnetic field have the opposite handedness to the velocity. Finally, widely used detection techniques for standing waves, both across and along the field, can fail to identify their presence. We explain how all these features arise from the MHD equations when accounting for a non-uniform background field and propose modified methods that might be applied to spacecraft observations.

ULF Waves↗
NASA NTRS2019

The Effect of Redox Conditions on Seismic Waves in Iron-Bearing Olivine: Implications for Understanding Planetary Interiros Through Seismilogy

Seismic data, inclusive of velocities and attenuation, can be utilized to elucidate the physical state of planetary interiors]. However, numerous micromechanical factors have been either experimentally demonstrated, or theoretically considered, to affect the propagation and dissipation of seismic energy within crystalline solids - including, but not limited to, changes in grain size, temperature, melt fraction, pressure and dislocation density. Thus, observed variations in seismic wave speeds and attenuation may be used to ultimately ‘map’ variations in physical properties, such as those listed above, within planetary bodies. But, in order to complete a successful inversion of seismic data into representations of physical properties, a first requirement is to obtain a fundamental laboratory based understanding of how each of these possible factors individually influences seismic waves. Here we conduct an experimental study with the initial objective to further understand one of the most commonly invoked, yet least studied, mechanisms that could alter intrinsic seismic wave attenuation: water content (occurring as chemically-bound hydroxyl). The historical basis for determining the effect of water on seismic properties was established predominantly through analogy with large-strain creep experiments conducted on olivine under water-saturated conditions. While these deformation experiments routinely demonstrate a weakening of olivine in the presence of water, they represent a fundamentally different deformation regime in comparison to the microstrains experienced due to a passing seismic wave. Thus, in order to directly assess the effects of water on seismic properties, small-strain experiments are required. Substantially modified seismic properties in the presence of water have been observed previously at low strains and low frequencies, but only in a single exploratory study conducted under water-saturated conditions. Thus, to properly test the theoretical predictions we conducted a systematic study of the seismic properties of olivine using low-frequency torsional oscillation on aggregates containing varying concentrations of bound hydroxyl, for the first time at under saturated conditions.

Cline, C. J., II↗
NASA NTRS1993

Water-Group Ion Distributions in the Mid-Cometosheath of Comet Halley

In the mid-cometosheath of comet Halley (1-2x10^5 km from the nucleus) the center- of-mass plasma frame is approximately the bulk flow velocity of the cometary ions, and the Alfven wave speed is an appreciable fraction of the flow speed. Here, the peaks of the water-group ion distributions observed by the Giotto Ion Mass Spectrometer are at velocities consistently below the expected pickup speed. It is shown that this effect is consistent with the scattering of the new pickup ions onto a bispherical shell distribution. The model does not fit the data inside similar 1.2x10^5 km however, possibly as a result of the growing importance of collisions or the presence of other processes such as scattering on obliquely-propagating magnetosonic waves.

Huddleston, D. E.↗
NASA NTRS1993

Water group ion distributions in the midcometosheath of comet Halley

In the midcometosheath of comet Halley (1 x 10(exp 5) to 2 x 10(exp 5) km from the nucleus) the center-of-mass plasma frame is approximately the bulk flow velocity of the cometary ions, and the Alfven wave speed is an appreciable fraction of the flow speed. Here, the peaks of the water group ion distributions observed by the Giotto ion mass spectrometer are at velocities consistently below the expected pickup speed. It is shown that this effect is consistent with the scattering of the new pickup ions onto a bispherical shell distribution. The model does not fit the data inside approximately 1.2 x 10(exp 5) km, however, possibly as a result of the growing importance of collisions or the presence of other processes such as scattering on obliquely propagating magnetosonic waves.

Huddleston, D. E.↗
NASA NTRS1997

Correspondence between discrete and continuous models of excitable media: trigger waves

We present a theoretical framework for relating continuous partial differential equation (PDE) models of excitable media to discrete cellular automata (CA) models on a randomized lattice. These relations establish a quantitative link between the CA model and the specific physical system under study. We derive expressions for the CA model's plane wave speed, critical curvature, and effective diffusion constant in terms of the model's internal parameters (the interaction radius, excitation threshold, and time step). We then equate these expressions to the corresponding quantities obtained from solution of the PDEs (for a fixed excitability). This yields a set of coupled equations with a unique solution for the required CA parameter values. Here we restrict our analysis to "trigger" wave solutions obtained in the limiting case of a two-dimensional excitable medium with no recovery processes. We tested the correspondence between our CA model and two PDE models (the FitzHugh-Nagumo medium and a medium with a "sawtooth" nonlinear reaction source) and found good agreement with the numerical solutions of the PDEs. Our results suggest that the behavior of trigger waves is actually controlled by a small number of parameters.

Non-NASA Center↗
NASA NTRS1995

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗