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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Fourier analysis of numerical algorithms for the Maxwell equations

The Fourier method is used to analyze the dispersive, dissipative, and isotropy errors of various spatial and time discretizations applied to the Maxwell equations on multi-dimensional grids. Both Cartesian grids and non-Cartesian grids based on hexagons and tetradecahedra are studied and compared. The numerical errors are quantitatively determined in terms of phase speed, wave number, propagation direction, gridspacings, and CFL number. The study shows that centered schemes are more efficient than upwind schemes. The non-Cartesian grids yield superior isotropy and higher accuracy than the Cartesian ones. For the centered schemes, the staggered grids produce less errors than the unstaggered ones. A new unstaggered scheme which has all the best properties is introduced. The study also demonstrates that a proper choice of time discretization can reduce the overall numerical errors due to the spatial discretization.

Liu, Yen↗
NASA NTRS1991

Pulmonary function in microgravity: KC-135 experience

We have commenced a KC-135 program that parallels and proceeds our Spacelab (SLS-1) pulmonary function experiment. Our first task was to elucidate the affect of normal gravitation on the shape of the maximum expiratory flow volume (MEFV) curve. Nine normal subjects performed multiple MEFV maneuvers at 0-G, 1-G, and approximately 1.7-G. The MEFV curves for each subject were filtered, aligned at RV, and ensemble-averaged to produce an average MEFV curve for each state, allowing differences to be studied. Most subjects showed a decrease in the FVC at 0-G, which we attribute to an increased intrathoracic blood volume. In most of these subjects, the mean lung volume associated with a given flow was lower at 0-G, over about the upper half of the vital capacity. This is similar to the change previously reported during heat out immersion and is consistent with the known affect of engorgement of the lung with blood, on elastic recoil. There were also consistent but highly individual changes in the position and magnitude of detailed features of the curve, the individual patterns being similar to those previously reported on transition from the erect to the supine position. This supports the idea that the location and motion of choke points which determine the detailed individual configuration of MEFV curves, can be significantly influenced by gravitational forces, presumably via the effects of change in longitudinal tension on local airway pressure-diameter behavior and wave speed. We have developed a flight mass spectrometer and have commenced a study of single breath gradients in gas exchange, inert gas washouts, and rebreathing cardiac outputs and lung volumes at 0-G, 1-G, and 1.7-G. Comparison of our results with those from SLS-1 should identify the opportunities and limitations of the KC-135 as an accessible microgravity resource.

Guy, Harold J.↗
NASA NTRS1995

Modeling of nonequilibrium space plasma flows

Godunov-type numerical solution of the 20 moment plasma transport equations. One of the centerpieces of our proposal was the development of a higher order Godunov-type numerical scheme to solve the gyration dominated 20 moment transport equations. In the first step we explored some fundamental analytic properties of the 20 moment transport equations for a low b plasma, including the eigenvectors and eigenvalues of propagating disturbances. The eigenvalues correspond to wave speeds, while the eigenvectors characterize the transported physical quantities. In this paper we also explored the physically meaningful parameter range of the normalized heat flow components. In the second step a new Godunov scheme type numerical method was developed to solve the coupled set of 20 moment transport equations for a quasineutral single-ion plasma. The numerical method and the first results were presented at several national and international meetings and a paper describing the method has been published in the Journal of Computational Physics. To our knowledge this is the first numerical method which is capable of producing stable time-dependent solutions to the full 20 (or 16) moment set of transport equations, including the full heat flow equation. Previous attempts resulted in unstable (oscillating) solutions of the heat flow equations. Our group invested over two man-years into the development and implementation of the new method. The present model solves the 20 moment transport equations for an ion species and thermal electrons in 8 domain extending from a collision dominated to a collisionless region (200 km to 12,000 km). This model has been applied to study O+ acceleration due to Joule heating in the lower ionosphere.

Gombosi, Tamas↗
NASA NTRS1994

A preliminary evaluation of ocean topography from the TOPEX/POSEIDON mission

We have analyzed 50 ten-day cycles of TOPEX/POSEIDON (T/P) altimeter data to evaluate the ocean dynamic topography and its temporal variations. We have employed data from both the U.S. and French altimeters along with the NASA precision orbits in this analysis. Errors in the diurnal and semidiurnal components of the Cartwright-Ray tide model have been significantly reduced using a correction developed from a combination of JGM-2 and OSU91A was employed, as well as a geoid model based solely on OSU91A. The long wavelengths of the comparisons to historical data, although geoid error still corrupts the dynamic topography for wavelengths shorter than 2500 km. The root mean square (RMS) variability is similar to previous results from Geosat, with bakground 'noise' approaching 3 cm RMS. The computed annual and semiannual variations are also similar to previous Geosat results, although the hemispheric distribution of the annual heating cycle is much better presented in the T/P results. They also compare reasonably well with the Levitus hydrographic compilation in the northern hemisphere, although the T/P variations generally have larger amplitudes. Ten-day average maps of variations in sea level compare well with simulations measurements at ocean tide gauges, with RMS differences of less than 4 cm and correlations greater than 0.6 for most of the island gauges. Time-longitude plots of these sea level variations at different latitudes in the Pacific clearly show the presence of equatorial Kelvin waves and Rossby waves, with the wave speeds agreeing well with theoretical and observed values. Measurement of variations in global sea level over cycles 2-51 have an RMS variability of 6.3 mm and a rate of change of -3.5 +/- 8 mm/yr, the uncertainty primarily due to insufficient averaging of the interannual and periodic sea level variations. These results show that the accuracy of the T/P measurements of sea level has dramatically improved over previous missions, with estimated time variable errors of 4 cm or less. Although geographically correlated orbit errors have also been reduced to the few centimeter level, further improvement in determinations of the mean dynamic topography will be difficult to obtain until a more accurate model of the marine geoid is available.

Nerem, R. S.↗
NASA NTRS1996

Inhomogeneous Radiation Boundary Conditions Simulating Incoming Acoustic Waves for Computational Aeroacoustics

A set of nonhomogeneous radiation and outflow conditions which automatically generate prescribed incoming acoustic or vorticity waves and, at the same time, are transparent to outgoing sound waves produced internally in a finite computation domain is proposed. This type of boundary condition is needed for the numerical solution of many exterior aeroacoustics problems. In computational aeroacoustics, the computation scheme must be as nondispersive ans nondissipative as possible. It must also support waves with wave speeds which are nearly the same as those of the original linearized Euler equations. To meet these requirements, a high-order/large-stencil scheme is necessary The proposed nonhomogeneous radiation and outflow boundary conditions are designed primarily for use in conjunction with such high-order/large-stencil finite difference schemes.

Tam, Christopher K. W.↗
NASA NTRS1997

Condensate Accretion in Shock Tube's Expansion Fan

It has been shown that turbulence and temperature influence the droplet sizes in expansion fan induced condensation by studying the Rayleigh scattering from one port in our shock tube's test section. We have modified our set-up so as to allow, using two ports, the real time measurement of the influence of turbulence and temperature on the rate at which these droplets grow. To do this, we looked at the Rayleigh scattering from two different ports for ten Reynolds numbers at five different temperatures. We modeled the time of flight of droplets, using the equations of one-dimensional gas dynamics and the measured shock wave speed in shock tube's driven section.

Mezonlin, Ephrem-Denis↗
NASA NTRS1997

Numerical Boundary Conditions for Computational Aeroacoustics Benchmark Problems

Category 1, Problems 1 and 2, Category 2, Problem 2, and Category 3, Problem 2 are solved computationally using the Dispersion-Relation-Preserving (DRP) scheme. All these problems are governed by the linearized Euler equations. The resolution requirements of the DRP scheme for maintaining low numerical dispersion and dissipation as well as accurate wave speeds in solving the linearized Euler equations are now well understood. As long as 8 or more mesh points per wavelength is employed in the numerical computation, high quality results are assured. For the first three categories of benchmark problems, therefore, the real challenge is to develop high quality numerical boundary conditions. For Category 1, Problems 1 and 2, it is the curved wall boundary conditions. For Category 2, Problem 2, it is the internal radiation boundary conditions inside the duct. For Category 3, Problem 2, they are the inflow and outflow boundary conditions upstream and downstream of the blade row. These are the foci of the present investigation. Special nonhomogeneous radiation boundary conditions that generate the incoming disturbances and at the same time allow the outgoing reflected or scattered acoustic disturbances to leave the computation domain without significant reflection are developed. Numerical results based on these boundary conditions are provided.

Tam, Chritsopher K. W.↗
NASA NTRS1996

Spectra of Baroclinic Inertia-Gravity Wave Turbulence

Baroclinic inertia-gravity (IG) waves form a persistent background of thermocline depth and sea surface height oscillations. They also contribute to the kinetic energy of horizontal motions in the subsurface layer. Measured by the ratio of water particle velocity to wave phase speed, the wave nonlinearity may be rather high. Given a continuous supply of energy from external sources, nonlinear wave-wave interactions among IG waves would result in inertial cascades of energy, momentum, and wave action. Based on a recently developed theory of wave turbulence in scale-dependent systems, these cascades are investigated and IG wave spectra are derived for an arbitrary degree of wave nonlinearity. Comparisons with satellite-altimetry-based spectra of surface height variations and with energy spectra of horizontal velocity fluctuations show good agreement. The well-known spectral peak at the inertial frequency is thus explained as a result of the inverse cascade. Finally, we discuss a possibility of inferring the internal Rossby radius of deformation and other dynamical properties of the upper thermocline from the spectra of SSH (sea surface height) variations based on altimeter measurements.

Glazman, Roman E.↗
NASA NTRS2001

Low Dissipative High Order Numerical Simulations of Supersonic Reactive Flows

The objective of this paper is to evaluate the performance of a newly developed low dissipative sixth-order spatial and fourth-order temporal scheme for viscous reactive flows interacting with shock waves that contain fine scale flow structures. The accuracy and efficiency of the scheme, and to what degree the scheme can capture the correct physical wave speeds of stiff reactive flows will be included.

Sjoegreen, B.↗
NASA NTRS2002

Advances in Time-Distance Helioseismology

Time-distance helioseismology is a way to measure travel times between surface locations for waves traversing the solar interior. Coupling the travel with an extensive modeling effort has proven to be a powerful tool for measuring flows and other wave speed inhomogeneities in the solar interior. Problems receiving current attention include studying the time variation of the meridional circulation and torsional oscillation and active region emergence and evolution, current results on these topics will be presented.

Duvall, Thomas L., Jr.↗
NASA NTRS2002

Computational Study of Near-limit Propagation of Detonation in Hydrogen-air Mixtures

A computational investigation of the near-limit propagation of detonation in lean and rich hydrogen-air mixtures is presented. The calculations were carried out over an equivalence ratio range of 0.4 to 5.0, pressures ranging from 0.2 bar to 1.0 bar and ambient initial temperature. The computations involved solution of the one-dimensional Euler equations with detailed finite-rate chemistry. The numerical method is based on a second-order spatially accurate total-variation-diminishing (TVD) scheme, and a point implicit, first-order-accurate, time marching algorithm. The hydrogen-air combustion was modeled with a 9-species, 19-step reaction mechanism. A multi-level, dynamically adaptive grid was utilized in order to resolve the structure of the detonation. The results of the computations indicate that when hydrogen concentrations are reduced below certain levels, the detonation wave switches from a high-frequency, low amplitude oscillation mode to a low frequency mode exhibiting large fluctuations in the detonation wave speed; that is, a 'galloping' propagation mode is established.

Yungster, S.↗
NASA NTRS2002

The Delta x B = 0 Constraint Versus Minimization of Numerical Errors in MHD Simulations

The MHD equations are a system of non-strictly hyperbolic conservation laws. The non-convexity of the inviscid flux vector resulted in corresponding Jacobian matrices with undesirable properties. It has previously been shown by Powell et al. (1995) that an 'almost' equivalent MHD system in non-conservative form can be derived. This non-conservative system has a better conditioned eigensystem. Aside from Powell et al., the MHD equations can be derived from basic principles in either conservative or non-conservative form. The Delta x B = 0 constraint of the MHD equations is only an initial condition constraint, it is very different from the incompressible Navier-Stokes equations in which the divergence condition is needed to close the system (i.e., to have the same number of equations and the same number of unknown). In the MHD formulations, if Delta x B = 0 initially, all one needs is to construct appropriate numerical schemes that preserve this constraint at later time evolutions. In other words, one does not need the Delta x B condition to close the MHD system. We formulate our new scheme together with the Cargo & Gallice (1997) form of the MHD approximate Riemann solver in curvilinear grids for both versions of the MHD equations. A novel feature of our new method is that the well-conditioned eigen-decomposition of the non-conservative MHD equations is used to solve the conservative equations. This new feature of the method provides well-conditioned eigenvectors for the conservative formulation, so that correct wave speeds for discontinuities are assured. The justification for using the non-conservative eigen-decomposition to solve the conservative equations is that our scheme has a better control of the numerical error associated with the divergence of the magnetic condition. Consequently, computing both forms of the equations with the same eigen-decomposition is almost equivalent. It will be shown that this approach, using the non-conservative eigensystem when solving the conservative equations, also works well in the context of standard shock-capturing schemes.

Yee, H. C.↗
NASA NTRS2002

The Delta(dot) B = O Constraint vs. Minimization of Numerical Errors in MHD Simulations

The MHD equations are a system of non-strictly hyperbolic conservation laws. The non-convexity of the inviscid flux vector resulted in corresponding Jacobian matrices with undesirable properties. On the other hand, the MHD equations can be derived from basic principles in either conservative or non-conservative form. The non-conservative system has a better conditioned eigensystem. The Delta(dot)B = 0 constraint of the A4HD equations is only an initial condition constraint. One does not need the Delta(dot)B condition to close the MHD system. We formulate our new low dissipative high order scheme together with the Cargo & Gallice (1997) form of the MHD approximate Riemann solver in curvilinear grids for both versions of the MHD equations. A novel feature of our new method is that the well-conditioned eigen-decomposition of the non-conservative MHD equations is used to solve the conservative equations. This new feature of the method provides well-conditioned eigenvectors for the conservative formulation, so that correct wave speeds for discontinuities are assured. The justification for using the non-conservative eigen-decomposition to solve the conservative equations is that our scheme has a better control of the numerical error associated with the Delta(dot)B condition. Consequently, computing both forms of the equations with the same eigen-decomposition is almost equivalent. It will be shown that this approach, using the non-conservative eigensystem when solving the conservative equations, also works well in the context of standard shock-capturing schemes.

Yee, H. C.↗
NASA NTRS2003

Comparison of Time-Distance Local Helioseismology on GONG and MDI Data Sets

We show first results derived from one rotation of GONG++ and MDI data analyzed independently by different groups with time-distance techniques. We focus on observations obtained during spring 2002 and especially on Carrington rotation 1988 (2002/3/30 - 2002/4/26) and measure flow components and wave speed inhomogeneities over a range of depths for different active regions.

Duvall, T. L., Jr.↗
NASA NTRS2004

Solar Oscillations

Oscillations were first detected in the solar photosphere in 1962 by Leighton and students. In 1970 it was calculated that these oscillations, with a period near five minutes, were the manifestations of acoustic waves trapped in the interior. The subsequent measurements of the frequencies of global oscillation modes from the spatio-temporal power spectrum of the waves made possible the refinement of solar interior models. Over the years, increased understanding of the nuclear reaction rates, the opacity, the equation of state, convection, and gravitational settling have resulted. Mass flows shift the frequencies of modes leading to very accurate measurements of the interior rotation as a function of radius and latitude. In recent years, analogues of terrestrial seismology have led to a tomography of the interior, including measurements of global north-south flows and flow and wave speed measurements below features such as sunspots. The future of helioseismology seems bright with the approval of NASA's Solar Dynamics Observatory mission, to be launched in 2008.

Duvall, Thomas↗
NASA NTRS2006

Semi-Analytic Reconstruction of Flux in Finite Volume Formulations

Semi-analytic reconstruction uses the analytic solution to a second-order, steady, ordinary differential equation (ODE) to simultaneously evaluate the convective and diffusive flux at all interfaces of a finite volume formulation. The second-order ODE is itself a linearized approximation to the governing first- and second- order partial differential equation conservation laws. Thus, semi-analytic reconstruction defines a family of formulations for finite volume interface fluxes using analytic solutions to approximating equations. Limiters are not applied in a conventional sense; rather, diffusivity is adjusted in the vicinity of changes in sign of eigenvalues in order to achieve a sufficiently small cell Reynolds number in the analytic formulation across critical points. Several approaches for application of semi-analytic reconstruction for the solution of one-dimensional scalar equations are introduced. Results are compared with exact analytic solutions to Burger s Equation as well as a conventional, upwind discretization using Roe s method. One approach, the end-point wave speed (EPWS) approximation, is further developed for more complex applications. One-dimensional vector equations are tested on a quasi one-dimensional nozzle application. The EPWS algorithm has a more compact difference stencil than Roe s algorithm but reconstruction time is approximately a factor of four larger than for Roe. Though both are second-order accurate schemes, Roe s method approaches a grid converged solution with fewer grid points. Reconstruction of flux in the context of multi-dimensional, vector conservation laws including effects of thermochemical nonequilibrium in the Navier-Stokes equations is developed.

Gnoffo, Peter A.↗
NASA NTRS2011

Time-Distance Helioseismology Data-Analysis Pipeline for Helioseismic and Magnetic Imager Onboard Solar Dynamics Observatory (SDO-HMI) and Its Initial Results

The Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) provides continuous full-disk observations of solar oscillations. We develop a data-analysis pipeline based on the time-distance helioseismology method to measure acoustic travel times using HMI Doppler-shift observations, and infer solar interior properties by inverting these measurements. The pipeline is used for routine production of near-real-time full-disk maps of subsurface wave-speed perturbations and horizontal flow velocities for depths ranging from 0 to 20 Mm, every eight hours. In addition, Carrington synoptic maps for the subsurface properties are made from these full-disk maps. The pipeline can also be used for selected target areas and time periods. We explain details of the pipeline organization and procedures, including processing of the HMI Doppler observations, measurements of the travel times, inversions, and constructions of the full-disk and synoptic maps. Some initial results from the pipeline, including full-disk flow maps, sunspot subsurface flow fields, and the interior rotation and meridional flow speeds, are presented.

Sun: helioseismology↗
NASA NTRS2014

Resonant Tidal Excitation of Internal Waves in the Earth's Fluid Core

It has long been speculated that there is a stably stratified layer below the core-mantle boundary, and two recent studies have improved the constraints on the parameters describing this stratification. Here we consider the dynamical implications of this layer using a simplified model. We first show that the stratification in this surface layer has sensitive control over the rate at which tidal energy is transferred to the core. We then show that when the stratification parameters from the recent studies are used in this model, a resonant configuration arrives whereby tidal forces perform elevated rates of work in exciting core flow. Specifically, the internal wave speed derived from the two independent studies (150 and 155 m/s) are in remarkable agreement with the speed (152 m/s) required for excitation of the primary normal mode of oscillation as calculated from full solutions of the Laplace Tidal Equations applied to a reduced-gravity idealized model representing the stratified layer. In evaluating this agreement it is noteworthy that the idealized model assumed may be regarded as the most reduced representation of the stratified dynamics of the layer, in that there are no non-essential dynamical terms in the governing equations assumed. While it is certainly possible that a more realistic treatment may require additional dynamical terms or coupling, it is also clear that this reduced representation includes no freedom for coercing the correlation described. This suggests that one must accept either (1) that tidal forces resonantly excite core flow and this is predicted by a simple model or (2) that either the independent estimates or the dynamical model does not accurately portray the core surface layer and there has simply been an unlikely coincidence between three estimates of a stratification parameter which would otherwise have a broad plausible range.

Geodynamo↗