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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
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.↗
NASA NTRSDate not supplied

A New Approach for a Wider Class of Entropy Split Methods for Compressible Gas Dynamics and MHD

The high order entropy split methods of Sjögreen & Yee [1, 2] by entropy splitting of the compressible Euler (inviscid) flux derivatives for a thermally-perfect gas are based on Harten’s entropy function [3, 4, 5]. Their derivation takes advantage of the homogeneity property of Euler flux, symmetrizable Euler flux derivatives and energy-norm stability in conjunction with high order classical spatial central, DRP (dispersion relation-preserving) [6, 7, 8] or Padé (compact) spatial discretizations [9] with summation-by-parts (SBP) operators [10]. Our entropy split methods have been proven entropy conserving and stable [1, 11, 12]. Our proofs do not rely on a two-point numerical flux, but rather only a linear difference operator is required to derive these methods. To extend the entropy split method for the MHD, we used the Godunov symmetrizable non-conservative MHD form [12, 13, 14]. These high order entropy split methods not only preserve certain physical properties of the chosen governing equations but are also known to either improve numerical stability, and/or minimize aliasing errors in long time integration of turbulent flow computations without the aid of added numerical dissipation. In our previous published work, extensive error norm comparison with grid refinement was performed to show the high accuracy performance of these methods. These studies also showed how well the entropy split methods conserve the entropy, momentum and mass, and preserve the kinetic energy for long time integration of the various flows [1, 2, 12, 13, 14]. The objective of the present work is to use a new approach to obtain a wider class of entropy split methods consisting of a two-point numerical flux portion and a non-conservative portion in such a way that the homogeneity property of the compressible Euler flux is not required. For high order classical spatial central, DRP (dispersion relation-preserving) or Padé (compact) spatial discretizations, this new approach can be proven to be entropy conservative with conservative spatial dsicretizations while at the same time allowing a wider class of symmetrizable inviscid flux derivatives. We also use this generalization to derive an entropy split scheme that is entropy conserving for the equations of MHD without the homogeneity property using the Godunov symmetrizable ideal MHD formulation [15].

High Order Physical Preserving Methods↗
NASA NTRS2004

Characteristics of Electromagnetic Pulse Propagation in Metal

It is well known that the solution of the diffusion equation for an electromagnetic field with a time harmonic term, e(sup iwt), is in the form of a traveling wave whose amplitude attenuates over distance into a conducting medium. As the attenuation is an increasing function of frequency, the high frequency components attenuate more rapidly than those of low ones upon entering a well conducting object. At the same time, the phase velocity of an individual component is also an increasing function of frequency causing a broadening of the pulse traveling inside a conductor. In the results of our previous study of numerical simulations, the problem of using a gaussian input pulse was immediately clear. First, having the dominant frequency components distributed around zero, the movement of the peak was not well defined. Second, with the amplitude of fourier components varying slowly over a wide range, the dispersion-induced blurring of the peak position was seen to be severe. For the present study, we have used a gaussian modulated single frequency sinusoidal wave, i. e., the carrier, as an input pulse in an effort to improve the issues related to the unclear movement of peak and dispersion as described above. This was based on the following two anticipated advantages: First, the packet moves in a conductor at the group velocity calculated at the carrier frequency, which means it is well controllable. Second, the amplitude of frequency components other than that of the carrier can be almost negligible, such that the effect of dispersion can be significantly reduced. A series of experiments of transmitting electromagnetic pulses through aluminum plates of various thickness was performed to test the validity of the above points. The results of numerical simulation based on wave propagation are discussed with respect to the experimental results. Finally, a simple simulation was performed based on diffusion of a continuous sine wave input and the results are compared with those of a single frequency sinusoidal wave observed over time at difference locations inside a conductor.

Namkung, M.↗
NASA NTRS1973

Principal sources and dispersal patterns of suspended particulate matter in nearshore surface waters of the northeast Pacific Ocean and the Hawaiian Islands

The author has identified the following significant results. Numerous geologic features can be discerned on ERTS-1 imagery of the northwestern Olympic Peninsula, Washington, and southern Vancouver Island, British Columbia. A thick homoclinal sequence of north-dipping Tertiary marine strata along the northwestern Olympic Peninsula is readily discernible because of the banded nature of its outcrop. The submarine basalt on which the sequence rests shows as a high, rugged ridge. Within the sedimentary sequence, alternating sandstone and siltstone members 100 m or more thick show in the southeast corner of the image. Broad folds in this banded sequence and faults that cut it can be detected. One large fault, the Pysht River fault, and a large syncline in sandstone and siltstone east of it show prominently on the image. Numerous lineations on the image of southern Vancouver Island correlate with faults shown on Muller's map. Particularly prominent are the Leech River and San Juan faults along the southernmost part of the island.

Carlson, P. R.↗
NASA NTRS1972

Experimental investigation of transient stimulated Raman scattering in a linearly dispersionless medium.

The results of numerical calculations of the transient stimulated Raman scattering reported previously have been verified experimentally under conditions where both linear dispersion and self-focusing effects were negligible. The existence of a delay between maxima of the laser and Stokes pulses is experimentally demonstrated for the first time, while the pulse shortening in time via Raman scattering is established more firmly than in previous work. The incident-laser-pulse duration, generated-Stokes-pulse duration, and delay between intensity maxima for the laser and Stokes pulses were measured with the two-photon absorption-fluorescence technique. The effective phonon-dephasing time is determined via spontaneous Raman scattering. By using these measured quantities, inferences are made as to the magnitude of the transient gain and the shape of the exciting picosecond laser pulse.

Carman, R. L.↗
NASA NTRS1993

Research on Earth's rotation and the effect of atmospheric pressure on vertical deformation and sea level variability

The work done under NASA grant NAG5-485 included modelling the deformation of the earth caused by variations in atmospheric pressure. The amount of deformation near coasts is sensitive to the nature of the oceanic response to the pressure. The PSMSL (Permanent Service for Mean Sea Level) data suggest the response is inverted barometer at periods greater than a couple months. Green's functions were constructed to describe the perturbation of the geoid caused by atmospheric and oceanic loading and by the accompanying load-induced deformation. It was found that perturbation of up to 2 cm are possible. Ice mass balance data was used for continental glaciers to look at the glacial contributions to time-dependent changes in polar motion, the lod, the earth's gravitational field, the position of the earth's center-of-mass, and global sea level. It was found that there can be lateral, non-hydrostatic structure inside the fluid core caused by gravitational forcing from the mantle, from the inner core, or from topography at the core/mantle or inner core/outer core boundaries. The nutational and tidal response of a non-hydrostatic earth with a solid inner core was modeled. Monthly, global tide gauge data from PSMSL was used to look at the 18.6-year ocean tide, the 14-month pole tide, the oceanic response to pressure, the linear trend and inter-annual variability in the earth's gravity field, the global sea level rise, and the effects of post glacial rebound. The effects of mantle anelasticity on nutations, earth tides, and tidal variation in the lod was modeled. Results of this model can be used with Crustal Dynamics observations to look at the anelastic dissipation and dispersion at tidal periods. The effects of surface topography on various components of crustal deformation was also modeled, and numerical models were developed of post glacial rebound.

Wahr, John↗
NASA NTRS1978

The velocity dispersion in Saturn's rings

The collisional dynamics of a differentially rotating disk of particles are considered. The investigation is to some extent an analytic counterpart of the numerical studies of Brahic (1977). Following Brahic, it is assumed that the particles are identical, indestructible, imperfectly elastic, smooth spheres. Gravitational interactions between the particles are neglected. The main result of the reported analysis is the derivation of the relation between the coefficient of restitution and the optical depth for a differentially rotating disk of identical, inelastic, smooth spheres. Questions concerning the applicability of the obtained results to Saturn's rings are studied. It is felt that the calculations should give accurate quantitative results if the rings' filling factor is small. If this assumption is not justified, the calculations should be at least qualitatively correct.

Goldreich, P.↗
NASA NTRS2001

Computational Aeroacoustics by the Space-time CE/SE Method

In recent years, a new numerical methodology for conservation laws-the Space-Time Conservation Element and Solution Element Method (CE/SE), was developed by Dr. Chang of NASA Glenn Research Center and collaborators. In nature, the new method may be categorized as a finite volume method, where the conservation element (CE) is equivalent to a finite control volume (or cell) and the solution element (SE) can be understood as the cell interface. However, due to its rigorous treatment of the fluxes and geometry, it is different from the existing schemes. The CE/SE scheme features: (1) space and time treated on the same footing, the integral equations of conservation laws are solve( for with second order accuracy, (2) high resolution, low dispersion and low dissipation, (3) novel, truly multi-dimensional, simple but effective non-reflecting boundary condition, (4) effortless implementation of computation, no numerical fix or parameter choice is needed, an( (5) robust enough to cover a wide spectrum of compressible flow: from weak linear acoustic waves to strong, discontinuous waves (shocks) appropriate for linear and nonlinear aeroacoustics. Currently, the CE/SE scheme has been developed to such a stage that a 3-13 unstructured CE/SE Navier-Stokes solver is already available. However, in the present paper, as a general introduction to the CE/SE method, only the 2-D unstructured Euler CE/SE solver is chosen as a prototype and is sketched in Section 2. Then applications of the CE/SE scheme to linear, nonlinear aeroacoustics and airframe noise are depicted in Sections 3, 4, and 5 respectively to demonstrate its robustness and capability.

Loh, Ching Y.↗
NASA NTRS1975

Development of a three-dimensional time-dependent flow field model

A three-dimensional, time-dependent mathematical model to represent Mobile Bay was developed. Computer programs were developed which numerically solve the appropriate conservation equations for predicting bay and estuary flow fields. The model is useful for analyzing the dispersion of sea water into fresh water and the transport of sediment, and for relating field and physical model data.

Farmer, R. C.↗
NASA NTRS1988

Approximation of the dielectric properties of Maxwellian plasmas - Dispersion functions and physical constraints

The dielectric properties of Maxwellian plasmas are approximated using high Pade approximants to the dispersion function and direct approximation of the distribution. Physical constraints on permissible approximations are discussed, and it is found that some previously published results can lead to predictions of qualitatively incorrect wave properties, including unphysical negative damping. Approximate dispersion functions for Maxwellian distributions are given explicitly, and some of the effects of these approximations on the resulting dispersion are discussed. The approximations discussed here are of use both in analytical work and in accelerating large-scale numerical computations.

Robinson, P. A.↗
NASA NTRS1991

A new approach to the determination of the striking distance from the lightning channel photos

The possible influence the stepped leader movement has on the value of the striking distance is studied. The numerical values of the striking distances are most frequently verified from the lightning channel photos. There is some explanations for dispersion of the striking distance which usually appear on tall structures. Unusually long striking distance may be influenced by the positive space charge or prior upward flashes. One more parameter which may have influence in estimation of the striking distance is a complex charge motion in the leader channel. Through introduction of two types of charge motions, possible shortening is shown of the striking distance seen on the lightning channel photos.

Vorgucic, A.↗
NASA NTRS1992

Stochastic gyroresonant electron acceleration in a low-beta plasma. I - Interaction with parallel transverse cold plasma waves

The gyroresonance of electrons with parallel transverse cold plasma waves is considered, and the Fokker-Planck equation describing the evolution of the electron distribution function in the presence of a spectrum of turbulence is derived. A new resonance which produces a divergence in the Fokker-Planck coefficients is identified; it results when the electron is in gyroresonance with a wave that has a group velocity equal to the velocity of the electron along the magnetic field. Under the assumption of a power-law spectral density, the Fokker-Planck coefficients are calculated numerically, and their complicated momentum and pitch-angle dependence, as well as the influence of various approximations to the dispersion relation, gyroresonance condition, and spectral density are discussed. It is found that there is no resonance gap at any pitch angle as long as the full gyroresonance condition is used and waves propagating on both directions are present.

Steinacker, Juergen↗
NASA NTRS1995

ICASE/LaRC Workshop on Benchmark Problems in Computational Aeroacoustics (CAA)

The proceedings of the Benchmark Problems in Computational Aeroacoustics Workshop held at NASA Langley Research Center are the subject of this report. The purpose of the Workshop was to assess the utility of a number of numerical schemes in the context of the unusual requirements of aeroacoustical calculations. The schemes were assessed from the viewpoint of dispersion and dissipation -- issues important to long time integration and long distance propagation in aeroacoustics. Also investigated were the effect of implementation of different boundary conditions. The Workshop included a forum in which practical engineering problems related to computational aeroacoustics were discussed. This discussion took the form of a dialogue between an industrial panel and the workshop participants and was an effort to suggest the direction of evolution of this field in the context of current engineering needs.

Hardin, Jay C.↗
NASA NTRS1971

Numerical solution for propagation of coupled longitudinal and transverse waves normal to the applied magnetic field in a three-fluid medium.

The solution was obtained without approximating the dispersion relation. This was made possible by the use of a computer method for carrying out the extremely large number of algebraic manipulations involved. The three-fluid model employed accounts for the interaction of the electron, ion, and neutral species. The solution for the e-folding distance describes the damping characteristics of each wave mode. Each of the e-folding solutions corresponds to one of the phase velocity solutions.

Mcclendon, R. W.↗
NASA NTRS1987

Comets

Observations of comets obtained with the IUE satellite since its launch in 1978 are reviewed. The status of UV observation of comets prior to IUE is discussed, and particular attention is given to low-resolution UV spectroscopy of cometary comae, the detection of new species in the UV emission, high-dispersion spectroscopy, spatial mapping of the emissions, abundance determinations, and short-term variability. Diagrams, graphs, sample spectra, and tables of numerical data are provided.

Festou, M. C.↗
NASA NTRS1999

Acceleration Processes in the Cusp -- Observations by the FAST Satellite

The FAST spacecraft has encountered the Earth's cusp regions near its apogee of 4175 km on numerous occasions during its first two and half years of operations. The cusp encounters are identified by their signatures of keV dispersed ion injections of solar wind origin. The FAST instruments reveal a complex microphysics inherent to many, but not all, of the cusp regions encountered by the spacecraft, that often include upgoing ion beams within regions of downgoing electrons that may appear as series of inverted-V features with energies near a few hundred eV. In many instances, upgoing electron beams have also been observed. Intense (> 100 mV/m) spikey DC-coupled electric fields and plasma waves are common features of the cusp encounters which also provide evidence for the presence of such local acceleration processes. In some cases, the FAST data show clear modulation of the precipitating magnetosheath ions indicative that they are affected by local electric potentials, as evidenced by simultaneous electron acceleration within such intervals. Furthermore, the acceleration events are sometimes organized with an apparent cellular structure that suggest Alfven waves or other large scale phenomena are controlling the localized potentials. We examine several cusp encounters in detail in order to study the complex relation of the cusp energetic particle populations with the plasma waves and DC electric fields.

Pfaff, R. F.↗
NASA NTRS1999

Acceleration Processes in the Cusp: Observations by the FAST Satellite

The Fast Auroral Snapshot (FAST) spacecraft has encountered the Earth's cusp regions near its apogee of 4175 km on numerous occasions during its first two and half years of operations. The cusp encounters are identified by their signatures of keV dispersed ion injections of solar wind origin. The FAST instruments reveal a complex microphysics inherent to many, but not all, of the cusp regions encountered by the spacecraft, that often include upgoing ion beams within regions of downgoing electrons that may appear as series of inverted-V features with energies near a few hundred eV. In many instances, upgoing electron beams have also been observed. Intense (> 100 mV/m) spikey DC-coupled electric fields and plasma waves are common features of the cusp encounters which also provide evidence for the presence of such local acceleration processes. In some cases, the FAST data show clear modulation of the precipitating magnetosheath ions indicative that they are affected by local electric potentials, as evidenced by simultaneous electron acceleration within such intervals. Furthermore, the acceleration events are sometimes organized with an apparent cellular structure that suggest Alfv6n waves or other large scale phenomena are controlling the localized potentials. We examine several cusp encounters in detail in order to study the complex relation of the cusp energetic particle populations with the plasma waves and DC electric fields.

Pfaff, R. F.↗
NASA NTRS2016

Vibroacoustic Characterization of Corrugated-Core and Honeycomb-Core Sandwich Panels

The vibroacoustic characteristics of two candidate launch vehicle fairing structures, corrugated- core and honeycomb-core sandwich designs, were studied. The study of these structures has been motivated by recent risk reduction efforts focused on mitigating high noise levels within the payload bays of large launch vehicles during launch. The corrugated-core sandwich concept is of particular interest as a dual purpose structure due to its ability to harbor resonant noise control systems without appreciably adding mass or taking up additional volume. Specifically, modal information, wavelength dispersion, and damping were determined from a series of vibrometer measurements and subsequent analysis procedures carried out on two test panels. Numerical and analytical modeling techniques were also used to assess assumed material properties and to further illuminate underlying structural dynamic aspects. Results from the tests and analyses described herein may serve as a reference for additional vibroacoustic studies involving these or similar structures.

Allen, Albert↗