Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Eigen values”

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.

18 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

The Benard problem: A comparison of finite difference and spectral collocation eigen value solutions

The application of spectral methods, using a Chebyshev collocation scheme, to solve hydrodynamic stability problems is demonstrated on the Benard problem. Implementation of the Chebyshev collocation formulation is described. The performance of the spectral scheme is compared with that of a 2nd order finite difference scheme. An exact solution to the Marangoni-Benard problem is used to evaluate the performance of both schemes. The error of the spectral scheme is at least seven orders of magnitude smaller than finite difference error for a grid resolution of N = 15 (number of points used). The performance of the spectral formulation far exceeded the performance of the finite difference formulation for this problem. The spectral scheme required only slightly more effort to set up than the 2nd order finite difference scheme. This suggests that the spectral scheme may actually be faster to implement than higher order finite difference schemes.

Skarda, J. Raymond Lee↗
NASA NTRS1993

Tsien's method for generating non-Keplerian trajectories. Part 2: The question of thrust to orbit a sphere and the restricted three-body problem

Tsien's method is extended to treat the orbital motion of a body undergoing accelerations and decelerations. A generalized solution is discussed for the generalized case where a body undergoes azimuthal and radial thrust and the problem is further simplified for azimuthal thrust alone. Judicious selection of thrust could generate either an elliptic or hyperbolic trajectory. This is unexpected especially when the body has only enough energy for a lower state trajectory. The methodology is extended treating the problem of vehicle thrust for orbiting a sphere and vehicle thrust within the classical restricted three-body problem. Results for the latter situation can produce hyperbolic trajectories through eigen value decomposition. Since eigen values for no-thrust can be imaginary, thrust can generate real eigen values to describe hyperbolic trajectories. Keplerian dynamics appears to represent but a small subset of a much larger non-Keplerian domain especially when thrust effects are considered. The need for high thrust long duration space-based propulsion systems for changing a trajectory's canonical form is clearly demonstrated.

Murad, P. A.↗
NASA NTRS1987

An analysis of the vertical structure equation for arbitrary thermal profiles

The vertical structure equation is a singular Sturm-Liouville problem whose eigenfunctions describe the vertical dependence of the normal modes of the primitive equations linearized about a given thermal profile. The eigenvalues give the equivalent depths of the modes. The spectrum of the vertical structure equation and the appropriateness of various upper boundary conditions, both for arbitrary thermal profiles were studied. The results depend critically upon whether or not the thermal profile is such that the basic state atmosphere is bounded. In the case of a bounded atmosphere it is shown that the spectrum is always totally discrete, regardless of details of the thermal profile. For the barotropic equivalent depth, which corresponds to the lowest eigen value, upper and lower bounds which depend only on the surface temperature and the atmosphere height were obtained. All eigenfunctions are bounded, but always have unbounded first derivatives. It was proved that the commonly invoked upper boundary condition that vertical velocity must vanish as pressure tends to zero, as well as a number of alternative conditions, is well posed. It was concluded that the vertical structure equation always has a totally discrete spectrum under the assumptions implicit in the primitive equations.

Cohn, Stephen E.↗
NASA NTRS1989

An analysis of the vertical structure equation for arbitrary thermal profiles

The vertical structure equation is a singular Sturm-Liouville problem whose eigenfunctions describe the vertical dependence of the normal modes of the primitive equations linearized about a given thermal profile. The eigenvalues give the equivalent depths of the modes. The spectrum of the vertical structure equation and the appropriateness of various upper boundary conditions, both for arbitrary thermal profiles were studied. The results depend critically upon whether or not the thermal profile is such that the basic state atmosphere is bounded. In the case of a bounded atmosphere it is shown that the spectrum is always totally discrete, regardless of details of the thermal profile. For the barotropic equivalent depth, which corresponds to the lowest eigen value, upper and lower bounds which depend only on the surface temperature and the atmosphere height were obtained. All eigenfunctions are bounded, but always have unbounded first derivatives. It was proved that the commonly invoked upper boundary condition that vertical velocity must vanish as pressure tends to zero, as well as a number of alternative conditions, is well posed. It was concluded that the vertical structure equation always has a totally discrete spectrum under the assumptions implicit in the primitive equations.

Cohn, Stephen E.↗
NASA NTRS1991

Field Telemetry of Blade-rotor Coupled Torsional Vibration at Matuura Power Station Number 1 Unit

The quasi-modal reduction technique and finite element model (FEM) were used to construct an analytical model for the blade-rotor coupled torsional vibration of a steam turbine generator of the Matuura Power Station. A single rotor test was executed in order to evaluate umbrella vibration characteristics. Based on the single rotor test results and the quasi-modal procedure, the total rotor system was analyzed to predict coupled torsional frequencies. Finally, field measurement of the vibration of the last stage buckets was made, which confirmed that the double synchronous resonance was 124.2 Hz, meaning that the machine can be safely operated. The measured eigen values are very close to the predicted value. The single rotor test and this analytical procedure thus proved to be a valid technique to estimate coupled torsional vibration.

Isii, Kuniyoshi↗
NASA NTRS1992

Eno-Osher schemes for Euler equations

The combination of the Osher approximate Riemann solver for the Euler equations and various ENO schemes is discussed for one-dimensional flow. The three basic approaches, viz. the ENO scheme using primitive variable reconstruction, either with Cauchy-Kowalewski procedure for time integration or the TVD Runge-Kutta scheme, and the flux-ENO method are tested on different shock tube cases. The shock tube cases were chosen to present a serious challenge to the ENO schemes in order to test their ability to capture flow discontinuities, such as shocks. Also the effect of the ordering of the eigen values, viz. natural or reversed ordering, in the Osher scheme is investigated. The ENO schemes are tested up to fifth order accuracy in space and time. The ENO-Osher scheme using the Cauchy-Kowalewski procedure for time integration is found to be the most accurate and robust compared with the other methods and is also computationally efficient. The tests showed that the ENO schemes perform reasonably well, but have problems in cases where two discontinuities are close together. In that case there are not enough points in the smooth part of the flow to create a non-oscillatory interpolation.

Vandervegt, Jacobus J.↗
NASA NTRS1993

A model of concurrent flow flame spread over a thin solid fuel

A numerical model is developed to examine laminar flame spread and extinction over a thin solid fuel in lowspeed concurrent flows. The model provides a more precise fluid-mechanical description of the flame by incorporating an elliptic treatment of the upstream flame stabilization zone near the fuel burnout point. Parabolic equations are used to treat the downstream flame, which has a higher flow Reynolds number. The parabolic and elliptic regions are coupled smoothly by an appropriate matching of boundary conditions. The solid phase consists of an energy equation with surface radiative loss and a surface pyrolysis relation. Steady spread with constant flame and pyrolysis lengths is found possible for thin fuels and this facilitates the adoption of a moving coordinate system attached to the flame with the flame spread rate being an eigen value. Calculations are performed in purely forced flow in a range of velocities which are lower than those induced in a normal gravity buoyant environment. Both quenching and blowoff extinction are observed. The results show that as flow velocity or oxygen percentage is reduced, the flame spread rate, the pyrolysis length, and the flame length all decrease, as expected. The flame standoff distance from the solid and the reaction zone thickness, however, first increase with decreasing flow velocity, but eventually decrease very near the quenching extinction limit. The short, diffuse flames observed at low flow velocities and oxygen levels are consistent with available experimental data. The maximum flame temperature decreases slowly at first as flow velocity is reduced, then falls more steeply close to the quenching extinction limit. Low velocity quenching occurs as a result of heat loss. At low velocities, surface radiative loss becomes a significant fraction of the total combustion heat release. In addition, the shorter flame length causes an increase in the fraction of conduction downstream compared to conduction to the fuel. These heat losses lead to lower flame temperatures, and ultimately, extinction. This extinction mechanism differs from that of blowoff, where the flame is unable to be stabilized due to the high flow velocity.

Ferkul, Paul V.↗
NASA NTRS1993

On the cross-stream spectral method for the Orr-Sommerfeld equation

Cross-stream models are defined as solutions to the Orr-Sommerfeld equation which are propagating normal to the flow direction. These models are utilized as a basis for a Hilbert space to approximate the spectrum of the Orr-Sommerfeld equation with plane Poiseuille flow. The cross-stream basis leads to a standard eigenvalue problem for the frequencies of Poiseuille flow instability waves. The coefficient matrix in the eigenvalue problem is shown to be the sum of a real matrix and a negative-imaginary diagonal matrix which represents the frequencies of the cross-stream modes. The real coefficient matrix is shown to approach a Toeplitz matrix when the row and column indices are large. The Toeplitz matrix is diagonally dominant, and the diagonal elements vary inversely in magnitude with diagonal position. The Poiseuille flow eigenvalues are shown to lie within Gersgorin disks with radii bounded by the product of the average flow speed and the axial wavenumber. It is shown that the eigenvalues approach the Gersgorin disk centers when the mode index is large, so that the method may be used to compute spectra with an essentially unlimited number of elements. When the mode index is large, the real part of the eigenvalue is the product of the axial wavenumber and the average flow speed, and the imaginary part of the eigen value is identical to the corresponding cross-stream mode frequency. The cross-stream method is numerically well-conditioned in comparison to Chebyshev based methods, providing equivalent accuracy for small mode indices and superior accuracy for large indices.

Zorumski, William E.↗
NASA NTRS1994

Sensitivity Analysis of Flutter Response of a Wing Incorporating Finite-Span Corrections

Flutter analysis of a wing is performed in compressible flow using state-space representation of the unsteady aerodynamic behavior. Three different expressions are used to incorporate corrections due to the finite-span effects of the wing in estimating the lift-curve slope. The structural formulation is based on a Rayleigh-Pitz technique with Chebyshev polynomials used for the wing deflections. The aeroelastic equations are solved as an eigen-value problem to determine the flutter speed of the wing. The flutter speeds are found to be higher in these cases, when compared to that obtained without accounting for the finite-span effects. The derivatives of the flutter speed with respect to the shape parameters, namely: aspect ratio, area, taper ratio and sweep angle, are calculated analytically. The shape sensitivity derivatives give a linear approximation to the flutter speed curves over a range of values of the shape parameter which is perturbed. Flutter and sensitivity calculations are performed on a wing using a lifting-surface unsteady aerodynamic theory using modules from a system of programs called FAST.

Issac, Jason Cherian↗
NASA NTRS1997

Convection in an Infinite Vertical Cylinder in the Presence of a Horizontal Magnetic Field

The use of a magnetic field to suppress convection in low gravity materials processing is an area of active research. This study aims to investigate some of the basic features of an verse magnetic field on convective motion in electrically conducting melts. The geometry considered is in infinite vertical column of liquid metal with an axial temperature gradient and an imposed transverse magnetic field. The incipience and azimuthal orientation of convective roM is investigated. The governing transport equations are nondimensionalized, linearized and the eigen value problem is solved to determine the most physically significant states. The analysis shows that the initial degeneracy (with no magnetic field) is removed by the presence of a magnetic field and the primary convection roll aligns itself so that its plane is perpendicular to the applied magnetic field. A cross over to another mode is also found when the strength of the magnetic field is increased. Velocity distributions are presented corresponding to the lowest modes and the system behavior within limit s of non-linearity is discussed.

Mazuruk, Konstantin↗
NASA NTRS2001

2-D/Axisymmetric Formulation of Multi-dimensional Upwind Scheme

A multi-dimensional upwind discretization of the two-dimensional/axisymmetric Navier-Stokes equations is detailed for unstructured meshes. The algorithm is an extension of the fluctuation splitting scheme of Sidilkover. Boundary conditions are implemented weakly so that all nodes are updated using the base scheme, and eigen-value limiting is incorporated to suppress expansion shocks. Test cases for Mach numbers ranging from 0.1-17 are considered, with results compared against an unstructured upwind finite volume scheme. The fluctuation splitting inviscid distribution requires fewer operations than the finite volume routine, and is seen to produce less artificial dissipation, leading to generally improved solution accuracy.

Wood, William A.↗
NASA NTRS1983

Study of LANDSAT-D thematic mapper performance as applied to hydrocarbon exploration

Many scenes of particular interest have a light dusting of snow cover. The possible use of hue-saturation-intensity transformations to reduce the effect of snow cover is being investigated. A tape of the Greeley, Colorado scene was reviewed on the interactive system and image types to be produced (decorrelated 2,3,4; natural color 1,2,3,; hue separation value 5/2,5/7 eigen 1; 4,5,7 in two color combinations) were selected. In several instances, a 1,3,4 combination produces a more useful false color infrared version of TM data than the more common 2,3,4 arrangement, probably because band 1 is less highly correlated with the band 3 and 4 than is band 2. A review of spacecraft performance suggests that the standard corrections applied at GSFC are more complicated than necessary in some areas and insufficient in other cases. The image motion compensation device on the TM works so well that bow-typing effects are very small; there are differences in the radiometry of forward and backward scans that make additional calibration necessary.

Source record↗
NASA NTRS1974

Propagation of waves in elliptic ducts. A theoretical study

The cut-off frequencies for high order circumferential modes were calculated for various eccentricities of an elliptic duct section. The problem was studied with a view to the reduction of jet engine compressor noise by elliptic ducts, instead of circular ducts. The cut-off frequencies for even functions decrease with increasing eccentricity. The third order eigen frequencies are oscillatory as the eccentricity increases for odd functions. The eigen frequencies decrease for higher order odd functions inasmuch as, for higher orders, they assume the same values as those for even functions. Deformation of a circular pipe into an elliptic one of sufficiently large eccentricity produces only a small reduction in the cut-off frequency, provided the area of the pipe section is kept invariable.

Baskaran, S.↗
NASA NTRS2002

Effects of Resonance in Quasi-Periodic Oscillators of Neutron Star Binaries

Using a large quantity of Rossi X-Ray Timing Explorer data presented in the literature, I offer a detailed investigation into the accuracy of the quasi-periodic oscillation (QPO) frequency determination. The QPO phenomenon seen in X-ray binaries is possibly a result of the resonance of the intrinsic (eigen) oscillations and harmonic driving forces of the system. I show that the resonances, in the presence of the damping of oscillations, occur at frequencies that are systematically and randomly shifted with respect to the eigenfrequencies of the system. The shift value strongly depends on the damping rate that is measured by the half-width of the QPO feature. Taking into account this effect, I analyze the QPO data for four Z sources, Scorpius X-I, GX 340+0, GX 5-1, and GX 17+2, and two atoll sources, 4U 1728-34 and 4U 0614+09. The transition-layer model (TLM) predicts the existence of the invariant quantity delta, an inclination angle of the magnetospheric axis with respect to the normal to the disk. I calculate delta and the error bars of delta using the resonance shift, and I find that the inferred delta-values are consistent with constants for these four Z sources, in which horizontal-branch oscillation and kilohertz frequencies have been detected and correctly identified. It is shown that the inferred delta are in the range between 5.5 deg and 6.5 deg. I conclude that the TLM seems to be compatible with the data.

Titarchuk, Lev↗
NASA NTRS2000

An Interpretation of Banded Magnetospheric Radio Emissions

Recently-published Active Magnetospheric Particle Tracer Explorer/Isothermal Remanent Magnetization (AMPTE/IRM) banded magnetospheric emissions, commonly referred to as '(n + 1/2)f(sub ce)' emissions where f(sub ce) is the electron gyrofrequency, are analyzed by treating them as analogous to sounder-stimulated ionospheric emissions. We show that both individual AMPTE/IRM spectra of magnetospheric banded emissions, and a statistically-derived spectra observed over the two-year lifetime of the mission, can be interpreted in a self-consistent manner. The analysis, which predicts all spectral peaks within 4% of the observed peaks, interprets the higher-frequency emissions as due to low group-velocity Bernstein-mode waves and the lower-frequency emissions as eigen modes of cylindrical-electromagnetic-plasma-oscillations. The demarcation between these two classes of emissions is the electron plasma frequency f(sub pe), where an emission is often observed. This f(sub pe), emission is not necessarily the strongest. None of the observed banded emissions were attributed to the upper-hybrid frequency. We present Alouette-2 and ISIS-1 plasma-resonance data, and model electron temperature (T(sub e)) values, to support the argument that the frequency-spectrum of ionospheric sounder-stimulated emissions is not strongly temperature dependent and thus that the interpretation of these emissions in the ionosphere is relevant to other plasmas (such as the magnetosphere) where N(sub e) and T(sub e) can be quite different but where the ratio f(sub pe)/f(sub ce) is identical.

Benson, Robert F.↗
NASA NTRS2010

Broken Ergodicity in Two-Dimensional Homogeneous Magnetohydrodynamic Turbulence

Two-dimensional (2-D) homogeneous magnetohydrodynamic (MHD) turbulence has many of the same qualitative features as three-dimensional (3-D) homogeneous MHD turbulence.The se features include several ideal invariants, along with the phenomenon of broken ergodicity. Broken ergodicity appears when certain modes act like random variables with mean values that are large compared to their standard deviations, indicating a coherent structure or dynamo.Recently, the origin of broken ergodicity in 3-D MHD turbulence that is manifest in the lowest wavenumbers was explained. Here, a detailed description of the origins of broken ergodicity in 2-D MHD turbulence is presented. It will be seen that broken ergodicity in ideal 2-D MHD turbulence can be manifest in the lowest wavenumbers of a finite numerical model for certain initial conditions or in the highest wavenumbers for another set of initial conditions.T he origins of broken ergodicity in ideal 2-D homogeneous MHD turbulence are found through an eigen analysis of the covariance matrices of the modal probability density functions.It will also be shown that when the lowest wavenumber magnetic field becomes quasi-stationary, the higher wavenumber modes can propagate as Alfven waves on these almost static large-scale magnetic structures

Shebalin, John V.↗
NASA NTRS2018

Decadal Seasonal Shifts of Precipitation and Temperature in TRMM and AIRS Data

We present results from an analysis of seasonal phase shifts in the global precipitation and surface temperatures. We use data from the TRMM (Tropical Rainfall Measuring Mission) Multi-satellite Precipitation Algorithm (TMPA), and the Atmospheric Infrared Sounder (AIRS) on Aqua satellite, all hosted at NASA Goddard Earth Science Data and Information Services Center (GES DISC). We explore the information content and data usability by first aggregating daily grids from the entire records of both missions to pentad (5-day) series which are then processed using Singular Value Decomposition approach. A strength of this approach is the normalized principal components that can then be easily converted from real to complex time series. Thus, we can separate the most informative, the seasonal, components and analyze unambiguously for potential seasonal phase drifts. TMPA and AIRS records represent correspondingly 20 and 15 years of data, which allows us to run simple “phase learning†from the first 5 years of records and use it as reference. The most recent 5 years are then phase-compared with the reference. We demonstrate that the seasonal phase of global precipitation and surface temperatures has been stable in the past two decades. However, a small global trend of delayed precipitation, and earlier arrival of surface temperatures seasons, are detectable at 95% confidence level. Larger phase shifts are detectable at regional level, in regions recognizable from the Eigen vectors to having strong seasonal patterns. For instance, in Central North America, including the North American Monsoon region, confident phase shifts of 1-2 days per decade are detected at 95% confidence level. While seemingly symbolic, these shifts are indicative of larger changes in the Earth Climate System. We thus also demonstrate a potential usability scenario of Earth Science Data Records curated at the NASA GES DISC in partnership with Earth Science Missions.

surface temperatures↗