Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dispersion relation”

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 343 records · Page 19

Vorticity Transport and Wave Emission In A Protoplanetary Disk

Higher order numerical algorithms (4th order in time, 3rd order in space) are applied to the Euler equations and are used to examine vorticity transport and wave motion in a non-self gravitating, initially isentropic Keplerian disk. In this talk we will examine the response of the disk to an isolated vortex with a circulation about equal to the rotation rate of Jupiter. The vortex is located on the 4 AU circle and the nebula is simulated from 1 to 24 AU. We show that the vortex emits pressure-supported density and Rossby-type wave packets before it decays within a few orbits. The acoustic density waves evolve into weak (non entropy preserving) shock waves that propagate over the entire disk. The Rossby waves remain in the vicinity of the initial vortex disturbance, but are rapidly damped. Temporal frequencies and spatial wavenumbers are derived from the nonlinear simulation data and correlated with analytical dispersion relations from the linearized Euler and energy equations.

Davis, S. S.↗

Fractal fluctuations in cardiac time series

Human heart rate, controlled by complex feedback mechanisms, is a vital index of systematic circulation. However, it has been shown that beat-to-beat values of heart rate fluctuate continually over a wide range of time scales. Herein we use the relative dispersion, the ratio of the standard deviation to the mean, to show, by systematically aggregating the data, that the correlation in the beat-to-beat cardiac time series is a modulated inverse power law. This scaling property indicates the existence of long-time memory in the underlying cardiac control process and supports the conclusion that heart rate variability is a temporal fractal. We argue that the cardiac control system has allometric properties that enable it to respond to a dynamical environment through scaling.

NASA Discipline Cardiopulmonary↗

Bose Condensation at He-4 Interfaces

Path Integral Monte Carlo was used to calculate the Bose-Einstein condensate fraction at the surface of a helium film at T = 0:77 K, as a function of density. Moving from the center of the slab to the surface, the condensate fraction was found to initially increase with decreasing density to a maximum value of 0.9, before decreasing. Long wavelength density correlations were observed in the static structure factor at the surface of the slab. A surface dispersion relation was calculated from imaginary-time density-density correlations. Similar calculations of the superfluid density throughout He-4 droplets doped with linear impurities (HCN)(sub n) are presented. After deriving a local estimator for the superfluid density distribution, we find a decreased superfluid response in the first solvation layer. This effective normal fluid exhibits temperature dependence similar to that of a two-dimensional helium system.

Draeger, E. W.↗

Direct Measurement of Wave Kernels in Time-Distance Helioseismology

Solar f-mode waves are surface-gravity waves which propagate horizontally in a thin layer near the photosphere with a dispersion relation approximately that of deep water waves. At the power maximum near 3 mHz, the wavelength of 5 Mm is large enough for various wave scattering properties to be observable. Gizon and Birch (2002,ApJ,571,966)h ave calculated kernels, in the Born approximation, for the sensitivity of wave travel times to local changes in damping rate and source strength. In this work, using isolated small magnetic features as approximate point-sourc'e scatterers, such a kernel has been measured. The observed kernel contains similar features to a theoretical damping kernel but not for a source kernel. A full understanding of the effect of small magnetic features on the waves will require more detailed modeling.

Duvall, T. L., Jr.↗

3D Relativistic Magnetohydrodynamic Simulations of Magnetized Spine-Sheath Relativistic Jets

Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor gamma = 2.5 jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are less than or equal to c/the square root of 3 in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The Alfven wave speed is less than or equal to 0.07 c in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfven wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.

Mizuno, Yosuke↗

Helicon Wave Physics Impacts on Electrodeless Thruster Design

Effective generation of helicon waves for high density plasma sources is determined by the dispersion relation and plasma power balance. Helicon wave plasma sources inherently require an applied magnetic field of .01-0.1 T, an antenna properly designed to couple to the helicon wave in the plasma, and an rf power source in the 10-100 s of MHz, depending on propellant choice. For a plasma thruster, particularly one with a high specific impulse (>2000 s), the physics of the discharge would also have to address the use of electron cyclotron resonance (ECR) heating and magnetic expansion. In all cases the system design includes an optimized magnetic field coil, plasma source chamber, and antenna. A preliminary analysis of such a system, calling on experimental data where applicable and calculations where required, has been initiated at Glenn Research Center. Analysis results showing the mass scaling of various components as well as thruster performance projections and their impact on thruster size are discussed.

Gilland, James↗

Helicon Wave Physics Impacts on Electrodeless Thruster Design

Effective generation of helicon waves for high density plasma sources is determined by the dispersion relation and plasma power balance. Helicon wave plasma sources inherently require an applied magnetic field of .01-0.1 T, an antenna properly designed to couple to the helicon wave in the plasma, and an rf power source in the 10-100 s of MHz, depending on propellant choice. For a plasma thruster, particularly one with a high specific impulse (>2000 s), the physics of the discharge would also have to address the use of electron cyclotron resonance (ECR) heating and magnetic expansion. In all cases the system design includes an optimized magnetic field coil, plasma source chamber, and antenna. A preliminary analysis of such a system, calling on experimental data where applicable and calculations where required, has been initiated at Glenn Research Center. Analysis results showing the mass scaling of various components as well as thruster performance projections and their impact on thruster size are discussed.

Gilland, James H.↗

Small-Scale Gravity Waves in ER-2 MMS/MTP Wind and Temperature Measurements during CRYSTAL-FACE

Lower stratospheric wind and temperature measurements made from NASA's high-altitude ER-2 research aircraft during the CRYSTAL-FACE campaign in July 2002 were analyzed to retrieve information on small scale gravity waves (GWs) at the aircraft's flight level (typically approximately 20 km altitude). For a given flight segment, the S-transform (a Gaussian wavelet transform) was used to search for and identify small horizontal scale GW events, and to estimate their apparent horizontal wavelengths. The horizontal propagation directions of the events were determined using the Stokes parameter method combined with the cross S-transform analysis. The vertical temperature gradient was used to determine the vertical wavelengths of the events. GW momentum fluxes were calculated from the cross S-transform. Other wave parameters such as intrinsic frequencies were calculated using the GW dispersion relation. More than 100GW events were identified. They were generally high frequency waves with vertical wavelength of approximately 5 km and horizontal wavelength generally shorter than 20 km. Their intrinsic propagation directions were predominantly toward the east, whereas their ground-based propagation directions were primarily toward the west. Among the events, approximately 20% of them had very short horizontal wavelength, very high intrinsic frequency, and relatively small momentum fluxes, and thus they were likely trapped in the lower stratosphere. Using the estimated GW parameters and the background winds and stabilities from the NCAR/NCEP reanalysis data, we were able to trace the sources of the events using a simple reverse ray-tracing. More than 70% of the events were traced back to convective sources in the troposphere, and the sources were generally located upstream of the locations of the events observed at the aircraft level. Finally, a probability density function of the reversible cooling rate due to GWs was obtained in this study, which may be useful for cirrus cloud models.

Wang, L.↗

A High-Resolution Capability for Large-Eddy Simulation of Jet Flows

A large-eddy simulation (LES) code that utilizes high-resolution numerical schemes is described and applied to a compressible jet flow. The code is written in a general manner such that the accuracy/resolution of the simulation can be selected by the user. Time discretization is performed using a family of low-dispersion Runge-Kutta schemes, selectable from first- to fourth-order. Spatial discretization is performed using central differencing schemes. Both standard schemes, second- to twelfth-order (3 to 13 point stencils) and Dispersion Relation Preserving schemes from 7 to 13 point stencils are available. The code is written in Fortran 90 and uses hybrid MPI/OpenMP parallelization. The code is applied to the simulation of a Mach 0.9 jet flow. Four-stage third-order Runge-Kutta time stepping and the 13 point DRP spatial discretization scheme of Bogey and Bailly are used. The high resolution numerics used allows for the use of relatively sparse grids. Three levels of grid resolution are examined, 3.5, 6.5, and 9.2 million points. Mean flow, first-order turbulent statistics and turbulent spectra are reported. Good agreement with experimental data for mean flow and first-order turbulent statistics is shown.

DeBonis, James R.↗

Solutions of the Taylor-Green Vortex Problem Using High-Resolution Explicit Finite Difference Methods

A computational fluid dynamics code that solves the compressible Navier-Stokes equations was applied to the Taylor-Green vortex problem to examine the code s ability to accurately simulate the vortex decay and subsequent turbulence. The code, WRLES (Wave Resolving Large-Eddy Simulation), uses explicit central-differencing to compute the spatial derivatives and explicit Low Dispersion Runge-Kutta methods for the temporal discretization. The flow was first studied and characterized using Bogey & Bailley s 13-point dispersion relation preserving (DRP) scheme. The kinetic energy dissipation rate, computed both directly and from the enstrophy field, vorticity contours, and the energy spectra are examined. Results are in excellent agreement with a reference solution obtained using a spectral method and provide insight into computations of turbulent flows. In addition the following studies were performed: a comparison of 4th-, 8th-, 12th- and DRP spatial differencing schemes, the effect of the solution filtering on the results, the effect of large-eddy simulation sub-grid scale models, and the effect of high-order discretization of the viscous terms.

DeBonis, James R.↗

The Role of Higher-Order Modes on the Electromagnetic Whistler-Cyclotron Wave Fluctuations of Thermal and Non-Thermal Plasmas

Two fundamental challenging problems of laboratory and astrophysical plasmas are the understanding of the relaxation of a collisionless plasmas with nearly isotropic velocity distribution functions and the resultant state of nearly equipartition energy density with electromagnetic plasma turbulence. Here, we present the results of a study which shows the role that higher-order-modes play in limiting the electromagnetic whistler-like fluctuations in a thermal and non-thermal plasma. Our main results show that for a thermal plasma the magnetic fluctuations are confined by regions that are bounded by the least-damped higher order modes. We further show that the zone where the whistler-cyclotron normal modes merges the electromagnetic fluctuations shifts to longer wavelengths as the beta(sub e) increases. This merging zone has been interpreted as the beginning of the region where the whistler-cyclotron waves losses their identity and become heavily damped while merging with the fluctuations. Our results further indicate that in the case of nonthermal plasmas, the higher-order modes do not confine the fluctuations due to the effective higher-temperature effects and the excess of suprathermal plasma particles. The analysis presented here considers the second-order theory of fluctuations and the dispersion relation of weakly transverse fluctuations, with wave vectors parallel to the uniform background magnetic field, in a finite temperature isotropic bi-Maxwellian and Tsallis-kappa-like magnetized electron-proton plasma. Our results indicate that the spontaneously emitted electromagnetic fluctuations are in fact enhanced over these quasi modes suggesting that such modes play an important role in the emission and absorption of electromagnetic fluctuations in thermal or quasi-thermal plasmas.

Vinas, Adolfo F.↗

Cryogenic Refractive Indices of S-LAH55, S-LAH55V, S-LAH59, S-LAM3, S-NBM51, S-NPH2, S-PHM52, and S-TIH14 Glasses

The Transiting Exoplanet Survey Satellite (TESS) is an explorer-class planet finder, whose principal goal is to detect small planets with bright host starts in the solar neighborhood. The TESS payload consists of four identical cameras and a Data Handling Unit (DHU) fitted with CCD detectors and associated electronics. Each camera consist of a lens assembly with seven optical elements that include various types of Ohara glass substrates. The successful implementation of a panchromatic and a thermal lens assembly design for these cameras requires a fairly accurate (up to 0.000001 (1e-6)) knowledge of the temperature- and wavelength-dependent of the refractive index in the wavelength and temperature range of operation. Hence, this paper is devoted to report on measurements of the refractive index over the wavelength range of 0.42-1.15 micrometers and temperature range of 110-300 K for the following Ohara glasses: S-LAH55, S-LAH55V, S-LAH59, S-LAM3, S-NBM51, S-NPH2, S-PHM52, and S-TIH14. The measurements were performed utilizing the Cryogenic High Accuracy Refraction Measuring System (CHARMS) facility at NASA's Goddard Space Flight Center. A dense coverage of the absolute refractive index for all these substrates in the aforementioned wavelength and temperature ranges was used to determine the thermo-optic coefficient (dndT) and dispersion relation (dnd) as a function of wavelength and temperature. A comparison of the measured indices with literature values, specifically the temperature-dependent refractive indices of S-PHM52 and S-TIH14 reported by Yamamuro et al. [Yamamuro et al., Opt. Eng. 45(8), 083401 (2006)], will be presented.

cryogenic temperatures↗

Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer

The objective of this work is to compare a high-order solver with a low-order solver for performing Large-Eddy Simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the highorder method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.

large eddy simulations↗

Cryogenic Refractive Indices of S-LAH55, S-LAH55V, S-LAH59, S-LAM3, S-NBM51, S-NPH2, S-PHM52, and S-TIH14 Glasses

The Transiting Exoplanet Survey Satellite (TESS) is an explorer-class planet finder, whose principal goal is to detect small planets with bright host starts in the solar neighborhood. The TESS payload consists of four identical cameras with seven optical elements each that include various types of Ohara glass substrates. The successful implementation both panchromatic and thermal lens assembly designs for these cameras requires a fairly accurate (up to 1E-6) knowledge of the temperature and wavelength dependence of the refractive index in the wavelength and temperature range of operation. Hence, this paper is devoted to report on measurements of the refractive index over the wavelength range of 0.42-1.15 um and temperature range of 110-310 K for the following Ohara glasses: S-LAH55, S-LAH55V, SLAH59, S-LAM3, S-NBM51, S-NPH2, S-PHM52, and S-TIH14. The measurements were performed utilizing the Cryogenic High Accuracy Refraction Measuring System (CHARMS) facility at NASA's Goddard Space Flight Center. A dense coverage of the absolute refractive index for the title substrates in the aforementioned wavelength and temperature ranges was used to determine the thermo-optic coefficient (dn/dT) and dispersion relation (dn/d lambda) as a function of wavelength and temperature. A comparison of the measured indices with literature values, specifically the temperature-dependent refractive indices of S-PHM52 and S-TIH14, will be presented.

geometrical optics↗

Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer

The objective of this work is to compare a high-order solver with a low-order solver for performing large-eddy simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the high-order method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.

large eddy simulations↗

Wide-Stopband Aperiodic Phononic Filters

We demonstrate that a phonon stopband can be synthesized from an aperiodic structure comprising a discrete set of phononic filter stages. Each element of the set has a dispersion relation that defines a complete bandgap when calculated under a Bloch boundary condition. Hence, the effective stopband width in an aperiodic phononic filter (PnF) may readily exceed that of a phononic crystal with a single lattice constant or a coherence scale. With simulations of multi-moded phononic waveguides, we discuss the effects of finite geometry and mode-converting junctions on the phonon transmission in PnFs. The principles described may be utilized to form a wide stopband in acoustic and surface wave media. Relative to the quantum of thermal conductance for a uniform mesoscopic beam, a PnF with a stopband covering 1.6-10.4 GHz is estimated to reduce the thermal conductance by an order of magnitude at 75 mK.

bandgap↗

Preferential Heating of Oxygen 5+ Ions by Finite-Amplitude Oblique Alfven Waves

Minor ions in the fast solar wind are known to have higher temperatures and to flow faster than protons in the interplanetary space. In this study we combine previous research on parametric instability theory and 2.5D hybrid simulations to study the onset of preferential heating of Oxygen 5+ ions by large-scale finite-amplitude Alfven waves in the collisionless fast solar wind. We consider initially non-drifting isotropic multi-species plasma, consisting of isothermal massless fluid electrons, kinetic protons and kinetic Oxygen 5+ ions. The external energy source for the plasma heating and energization are oblique monochromatic Alfven-cyclotron waves. The waves have been created by rotating the direction of initial parallel pump, which is a solution of the multi-fluid plasma dispersion relation. We consider propagation angles theta less than or equal to 30 deg. The obliquely propagating Alfven pump waves lead to strong diffusion in the ion phase space, resulting in highly anisotropic heavy ion velocity distribution functions and proton beams. We discuss the application of the model to the problems of preferential heating of minor ions in the solar corona and the fast solar wind.

Maneva, Yana G.↗

Cryogenic Refractive Index of Heraeus Homosil Glass

This paper reports measurements of the refractive index of Homosil (Heraeus) over the wavelength range of 0.34-3.16 microns and temperature range of 120-335 K. These measurements were performed by using the Cryogenic High Accuracy Refraction Measuring System (CHARMS) facility at the NASAs Goddard Space Flight Center. These measurements were in support of an integrated Structural-Thermal-Optical-Performance (STOP) model that was developed for a field-widened Michelson interferometer that is being built and tested for the High Spectral Resolution Lidar (HSRL) project at the NASA Langley Research Center (LaRC). The cryogenic refractive index measurements were required in order to account for the highly sensitivity performance of the HSRL instrument to changes in refractive index with temperature, temperature gradients, thermal expansion, and deformation due to mounting stresses. A dense coverage of the absolute refractive index over the aforementioned wavelength and temperature ranges was used to determine the thermo-optic coefficient (dn/dT) and dispersion relation (dn/d(lamda)) as a function of wavelength and temperature. Our measurements of Homosil will be compared with measurements of other glasses from the fused silica family studied in CHARMS as well as measurements reported elsewhere in literature.

cryogenic↗