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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

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.↗
NASA NTRS2015

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↗
NASA NTRS2015

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↗
NASA NTRS2015

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↗
NASA NTRS2015

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↗
NASA NTRS2016

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↗
NASA NTRS2016

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.↗
NASA NTRS2017

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↗
NASA NTRS2017

Cryogenic Refractive Index of Heraeus Homosil Glass

This paper reports measurements of the refractive index of Homosil (Heraeus) over the wavelength range of 0.343.16 m and temperature range of 120335 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 (dndT) and dispersion relation (dnd) 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.

Homosil↗
NASA NTRS2019

WRLES: Wave Resolving Large-Eddy Simulation Code, Theory and Usage

A computational fluid dynamics code has been developed for large-eddy simulations (LES) of turbulent flow. The code uses high-order of accuracy and high-resolution numerical methods to minimize solution error and maximize the resolution of the turbulent structures. Spatial discretization is performed using explicit central differencing. The central differencing schemes in the code include 2nd- to 12th-order standard central difference methods as well as 7-, 9-, 11- and 13-point dispersion relation preserving schemes. Solution filtering and high-order shock capturing are included for stability. Time discretization is performed using multistage Runge-Kutta methods that are up to 4th order accurate. Several options are available to model turbulence including: Baldwin-Lomax and Spalart-Allmaras Reynolds-averaged Navier-Stokes turbulence models, and Smagorinsky, Dynamic Smagorinsky and Vreman sub-grid scale models for LES. This report presents the theory behind the numerical and physical models used in the code and provides a user's manual to the operation of the code.

High-Order↗
NASA NTRS2019

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

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Ols- son & Oliger (1994), Gerritsen & Olsson (1996), and Yee et al. (2000) is revisited. The entropy splitting is a form of skew-symmetric splitting of the nonlinear Euler flux derivatives. Central differencing applied to the entropy splitting form of the Euler flux derivatives together with SBP difference operators will, hereafter, be referred to as entropy split schemes. This study is prompted by the recent growing interest in numerical methods for which a discrete entropy conservation law holds, a discrete global entropy conservation can be proved and/or the numerical method possesses a stable entropy in the framework of SBP difference operators and L2-energy norm estimate. The objective of this paper is to recast the entropy split scheme as the re- cent definition of an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for non- linear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grub- ber splitting.

Sjogreen, Bjorn↗
NASA NTRS2019

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

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

Sjogreen, Bjorn↗
NASA NTRSDate not supplied

Methodology and Development of SLS Ascent Loads

The methodology and development of the Space Launch System (SLS) Block 1 Boost and Core phase Ascent loads will be presented in this paper. As is common practice, a series of fixed time/mass Finite Element Models of the vehicle were developed for the purpose of evaluating the structural loads and accelerations for particular Mach or time ranges associated with the Boost and Core phases of flight. Load contributors, which vary between bins, include the flight mechanic induced portion of the loads (e.g., Static-Elastic (STEL) or Vehicle Load Indicator (VLI) tool results), gust loads, buffet loads, Solid Rocket Booster thrust oscillations, Core Stage Engine thrust oscillations, maneuvering loads, thrust vector-related dispersions, and Programmed Test Input induced loads. Each of these load contributors will be examined in this paper, as well as the Loads Combination Equations used to generate the final loads. In addition to the typical coupled loads analysis, an Ascent loads evaluation is included in the SLS Artemis 1 Day of Launch process; a summary of this evaluation will also be presented in this paper.

SLS↗
NASA NTRS2020

Cross-field Streaming Instabilities in the Lower Hybrid Frequency Near the Front Pole Regions of a Magnetically Shielded Hall Thruster

The mechanism that led to the erosion of the inner front pole cover in wear tests of the Hall-effect Rocket with Magnetic Shielding (HERMeS) at 300 V and 12.5 kW has been elusive. Several classical mechanisms have been considered yet none has explained conclusively the source of the discrepancy between wear measurements and simulations. In this article we consider a new mechanism. The electrostatic dispersion relation for an unbounded, homogeneous plasma in a fixed magnetic field is solved numerically in the near-plume region of this thruster. The plasma conditions have been obtained from Hall2De (r-z) multifluid-PIC simulations. We find that the modified two-stream instability can be excited in most regions of the domain, with the fastest growth approaching the fluid-limit value of ωLH/2 near the magnetic poles where ωpe/ωce~1. In this region the relative drift between magnetized electrons and beam ions perpendicular to the magnetic field can exceed 3uTi, while that between counter-streaming beam and cathode ions can exceed 4uTi. The latter suggests that the ion-ion cross-field lower hybrid instability may also be active here. Unlike the classical ion-acoustic instability, lower hybrid instabilities can heat ions to temperatures that are comparable to the electron temperature. Such anomalous heating would enhance the average energy of these ions and, in turn, the erosion of the nearby surfaces they sputter. An idealized model that scales the anomalous ion heating with ωLH has been incorporated in Hall2De and the simulations now yield significantly improved agreement with the measured wear rates and IVDFs, as reported in our companion paper.

Ortega, Alejandro Lopez↗
NASA NTRS2022

On the Existence of Fast Modes in Compressible Magnetohydrodynamic Turbulence

We study the existence and properties of fast magnetosonic modes in 3D compressible MHD turbulence by carrying out a number of simulations with compressible and incompressible driving conditions. We use two approaches to determine the presence of fast modes: mode decomposition based on spatial variations only and spatio-temporal 4D fast Fourier transform (4D FFT) analysis of all fluctuations. The latter method enables us to quantify fluctuations that satisfy the dispersion relation of fast modes with finite frequency. Overall, we find that the fraction of fast modes identified via the spatio-temporal 4D FFT approach in total fluctuation power is either tiny with nearly incompressible driving or ∼2% with highly compressible driving. We discuss the implications of our results for understanding the compressible fluctuations in space and astrophysical plasmas.

Zhaoming Gan↗
NASA NTRS2022

Rapid Conjugate Appearance of the Giant Ionospheric Lamb Wave Signatures in the Northern Hemisphere After Hunga-Tonga Volcano Eruptions

The explosive eruption of the Hunga-Tonga volcano in the southwest Pacific at 0415UT on 15 January 2022 triggered gigantic atmospheric disturbances with surface air pressure waves propagating around the globe in Lamb mode. In space, concentric traveling ionosphere disturbances (CTIDs) are also observed as a manifestation of air pressure waves in New Zealand ∼0500UT and Australia ∼0630UT. As soon as the air pressure waves reached central Australia ∼0800UT, conjugate CTIDs appeared almost simultaneously in the northern hemispheres through interhemispheric coupling, much earlier than the arrival of the surface air pressure waves to Japan after 1100UT. Combining observations over Australia and Japan between 0800 and 1000UT, both direct and conjugate CTIDs show similar horizontal phase velocities of 320–390 m/s, matching with the dispersion relation of Lamb mode. The arrival of atmospheric Lamb wave to Japan later created in situ CTIDs showing the same Lamb mode characteristics as the earlier conjugate CTIDs.

Concentric traveling ionospheric disturbances↗
NASA NTRSDate not supplied

On A Stabilization of the Ingard-Myers Impedance Boundary Condition

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

Time Domain BEM↗