SEARCH · Engineering Papers

Results for “eigenvalues”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,729 records · Page 96

All sources NASA NTRS DOE OSTI

NASA NTRS1995

On pseudosupersymmetric oscillators and reality of relativistic energies for vector mesons

Specific oscillators - hereafter called pseudosupersymmetric oscillators - appear as interesting nonrelativistic concepts in connection with the study of relativistic vector mesons interacting with an external constant magnetic field when the real character of the energy eigenvalues is required as expected. A new pseudosupersymmetric quantum mechanics can then be developed and the corresponding pseudosupersymmetries can be pointed out.

Beckers, Jules↗

NASA NTRS1994

Benchmark problems in computational aeroacoustics

A recent directive at NASA Langley is aimed at numerically predicting principal noise sources. During my summer stay, I worked with high-order ENO code, developed by Dr. Harold Atkins, for solving the unsteady compressible Navier-Stokes equations, as it applies to computational aeroacoustics (CAA). A CAA workshop, composed of six categories of benchmark problems, has been organized to test various numerical properties of code. My task was to determine the robustness of Atkins' code for these test problems. In one category, we tested the nonlinear wave propagation of the code for the one-dimensional Euler equations, with initial pressure, density, and velocity conditions. Using freestream boundary conditions, our results were plausible. In another category, we solved the linearized two-dimensional Euler equations to test the effectiveness of radiation boundary conditions. Here we utilized MAPLE to compute eigenvalues and eigenvectors of the Jacobian given variable and flux vectors. We experienced a minor problem with inflow and outflow boundary conditions. Next, we solved the quasi one dimensional unsteady flow equations with an incoming acoustic wave of amplitude 10(exp -6). The small amplitude sound wave was incident on a convergent-divergent nozzle. After finding a steady-state solution and then marching forward, our solution indicated that after 30 periods the acoustic wave had dissipated (a period is time required for sound wave to traverse one end of nozzle to other end).

Porter-Locklear, Freda↗

NASA NTRS1993

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗

NASA NTRS1995

Separation analysis, a tool for analyzing multigrid algorithms

The separation of vectors by multigrid (MG) algorithms is applied to the study of convergence and to the prediction of the performance of MG algorithms. The separation operator for a two level cycle algorithm is derived. It is used to analyze the efficiency of the cycle when mixing of eigenvectors occurs. In particular cases the separation analysis reduces to Fourier type analysis. The separation operator of a two level cycle for a Schridubger eigenvalue problem, is derived and analyzed in a Fourier basis. Separation analysis gives information on how to choose performance relaxations and inter-level transfers. Separation analysis is a tool for analyzing and designing algorithms, and for optimizing their performance.

Costiner, Sorin↗

NASA NTRS1994

Control theory and splines, applied to signature storage

In this report the problem we are going to study is the interpolation of a set of points in the plane with the use of control theory. We will discover how different systems generate different kinds of splines, cubic and exponential, and investigate the effect that the different systems have on the tracking problems. Actually we will see that the important parameters will be the two eigenvalues of the control matrix.

Enqvist, Per↗

NASA NTRS1995

On the nonlinear stability of the unsteady, viscous flow of an incompressible fluid in a curved pipe

The stability of the flow of an incompressible, viscous fluid through a pipe of circular cross-section curved about a central axis is investigated in a weakly nonlinear regime. A sinusoidal pressure gradient with zero mean is imposed, acting along the pipe. A WKBJ perturbation solution is constructed, taking into account the need for an inner solution in the vicinity of the outer bend, which is obtained by identifying the saddle point of the Taylor number in the complex plane of the cross-sectional angle co-ordinate. The equation governing the nonlinear evolution of the leading order vortex amplitude is thus determined. The stability analysis of this flow to periodic disturbances leads to a partial differential system dependent on three variables, and since the differential operators in this system are periodic in time, Floquet theory may be applied to reduce this system to a coupled infinite system of ordinary differential equations, together with homogeneous uncoupled boundary conditions. The eigenvalues of this system are calculated numerically to predict a critical Taylor number consistent with the analysis of Papageorgiou. A discussion of how nonlinear effects alter the linear stability analysis is also given, and the nature of the instability determined.

Shortis, Trudi A.↗

NASA NTRS1995

Analysis of gas absorption to a thin liquid film in the presence of a zero-order chemical reaction

The paper presents a detailed theoretical analysis of the process of gas absorption to a thin liquid film adjacent to a horizontal rotating disk. The film is formed by the impingement of a controlled liquid jet at the center of the disk and subsequent radial spreading of liquid along the disk. The chemical reaction between the gas and the liquid film can be expressed as a zero-order homogeneous reaction. The process was modeled by establishing equations for the conservation of mass, momentum, and species concentration and solving them analytically. A scaling analysis was used to determine dominant transport processes. Appropriate boundary conditions were used to solve these equations to develop expressions for the local concentration of gas across the thickness of the film and distributions of film height, bulk concentration, and Sherwood number along the radius of the disk. The partial differential equation for species concentration was solved using the separation of variables technique along with the Duhamel's theorem and the final analytical solution was expressed using confluent hypergeometric functions. Tables for eigenvalues and eigenfunctions are presented for a number of reaction rate constants. A parametric study was performed using Reynolds number, Ekman number, and dimensionless reaction rate as parameters. At all radial locations, Sherwood number increased with Reynolds number (flow rate) as well as Ekman number (rate of rotation). The enhancement of mass transfer due to chemical reaction was found to be small when compared to the case of no reaction (pure absorption), but the enhancement factor was very significant when compared to pure absorption in a stagnant liquid film. The zero-order reaction processes considered in the present investigation included the absorption of oxygen in aqueous alkaline solutions of sodiumdithionite and rhodium complex catalyzed carbonylation of methanol. Present analytical results were compared to previous theoretical results for limiting conditions, and were found to have very good agreement.

Rajagopalan, S.↗

NASA NTRS1995

Aeroelastic stability of wind turbine blade/aileron systems

Aeroelastic stability analyses have been performed for the MOD-5A blade/aileron system. Various configurations having different aileron torsional stiffness, mass unbalance, and control system damping have been investigated. The analysis was conducted using a code recently developed by the General Electric Company - AILSTAB. The code extracts eigenvalues for a three degree of freedom system, consisting of: (1) a blade flapwise mode; (2) a blade torsional mode; and (3) an aileron torsional mode. Mode shapes are supplied as input and the aileron can be specified over an arbitrary length of the blade span. Quasi-steady aerodynamic strip theory is used to compute aerodynamic derivatives of the wing-aileron combination as a function of spanwise position. Equations of motion are summarized herein. The program provides rotating blade stability boundaries for torsional divergence, classical flutter (bending/torsion) and wing/aileron flutter. It has been checked out against fixed-wing results published by Theodorsen and Garrick. The MOD-5A system is stable with respect to divergence and classical flutter for all practical rotor speeds. Aileron torsional stiffness must exceed a minimum critical value to prevent aileron flutter. The nominal control system stiffness greatly exceeds this minimum during normal operation. The basic system, however, is unstable for the case of a free (or floating) aileron. The instability can be removed either by the addition of torsional damping or mass-balancing the ailerons. The MOD-5A design was performed by the General Electric Company, Advanced Energy Program Department under Contract DEN3-153 with NASA Lewis Research Center and sponsored by the Department of Energy.

Strain, J. C.↗

NASA NTRS1995

Combustion-acoustic stability analysis for premixed gas turbine combustors

Lean, prevaporized, premixed combustors are susceptible to combustion-acoustic instabilities. A model was developed to predict eigenvalues of axial modes for combustion-acoustic interactions in a premixed combustor. This work extends previous work by including variable area and detailed chemical kinetics mechanisms, using the code LSENS. Thus the acoustic equations could be integrated through the flame zone. Linear perturbations were made of the continuity, momentum, energy, chemical species, and state equations. The qualitative accuracy of our approach was checked by examining its predictions for various unsteady heat release rate models. Perturbations in fuel flow rate are currently being added to the model.

Darling, Douglas↗

NASA NTRS1990

Partitioning sparse matrices with eigenvectors of graphs

The problem of computing a small vertex separator in a graph arises in the context of computing a good ordering for the parallel factorization of sparse, symmetric matrices. An algebraic approach for computing vertex separators is considered in this paper. It is shown that lower bounds on separator sizes can be obtained in terms of the eigenvalues of the Laplacian matrix associated with a graph. The Laplacian eigenvectors of grid graphs can be computed from Kronecker products involving the eigenvectors of path graphs, and these eigenvectors can be used to compute good separators in grid graphs. A heuristic algorithm is designed to compute a vertex separator in a general graph by first computing an edge separator in the graph from an eigenvector of the Laplacian matrix, and then using a maximum matching in a subgraph to compute the vertex separator. Results on the quality of the separators computed by the spectral algorithm are presented, and these are compared with separators obtained from other algorithms for computing separators. Finally, the time required to compute the Laplacian eigenvector is reported, and the accuracy with which the eigenvector must be computed to obtain good separators is considered. The spectral algorithm has the advantage that it can be implemented on a medium-size multiprocessor in a straightforward manner.

Pothen, Alex↗

NASA NTRS1994

Nonlinear equilibration of localized instabilities on a baroclinic jet

Dynamical mechanisms underlying the equilibration of absolute instability are examined in a nonlinear, quasigeostrophic, two-layer model. The key to understanding the nonlinear equilibration is in recognizing that linear absolute instabilities can be stabilized both by a reduction of the vertical shear and by enhancement of the mean barotropic velocity. In a localized domain, the equilibration process proceeds with the creation of locally convectively unstable regions downstream, which encroach onto the locally absolutely unstable region until the local instability is suppressed. That local instabilities exist only if absolutely unstable regions span a minimum size is verified by eigenvalue calculations of three-dimensional flows. Numerical examples suggest that this critical size is at least 9000 km for a wide range of parameter values chosen to investigate the midlatitude storm tracks. Fluctuations arising from local absolute instability obtain maximum amplitude in the downstream convectively unstable regions rather than in the absolutely unstable regions themselves. Together, these results suggest that if an equilibrated absolute instability were to occur in midlatitudes, a zonal band of surface easterlies exceeding 900 km would be required and the associated enhanced variances would not be found coincident with the regions of absolute instability.

Delsole, Timothy M.↗

NASA NTRS1994

Nonlinear stability of Halley comethosheath with transverse plasma motion

Weakly nonlinear Magneto Hydrodynamic (MHD) stability of the Halley cometosheath determined by the balance between the outward ion-neutral drag force and the inward Lorentz force is investigated including the transverse plasma motion as observed in the flanks with the help of the method of multiple scales. The eigenvalues and the eigenfunctions are obtained for the linear problem and the time evolution of the amplitude is obtained using the solvability condition for the solution of the second order problem. The diamagnetic cavity boundary and the adjacent layer of about 100 km thickness is found unstable for the travelling waves of certain wave numbers. Halley ionopause has been observed to have strong ripples with a wavelength of several hundred kilometers. It is found that nonlinear effects have stabilizing effect.

Srivastava, Krishna M.↗

NASA NTRS1994

Large-angle cosmic microwave background anisotropies in an open universe

If the universe is open, scales larger than the curvature scale may be probed by observation of large-angle fluctuations in the cosmic microwave background (CMB). We consider primordial adiabatic perturbations and discuss power spectra that are power laws in volume, wavelength, and eigenvalue of the Laplace operator. Such spectra may have arisen if, for example, the universe underwent a period of `frustated' inflation. The resulting large-angle anisotropies of the CMB are computed. The amplitude generally increases as Omega is decreased but decreases as h is increased. Interestingly enough, for all three Ansaetze, anisotropies on angular scales larger than the curvature scale are suppressed relative to the anisotropies on scales smaller than the curvature scale, but cosmic variance makes discrimination between various models difficult. Models with 0.2 approximately less than Omega h approximately less than 0.3 appear compatible with CMB fluctuations detected by Cosmic Background Explorer Satellite (COBE) and the Tenerife experiment and with the amplitude and spectrum of fluctuations of galaxy counts in the APM, CfA, and 1.2 Jy IRAS surveys. COBE normalization for these models yields sigma(sub 8) approximately = 0.5 - 0.7. Models with smaller values of Omega h when normalized to COBE require bias factors in excess of 2 to be compatible with the observed galaxy counts on the 8/h Mpc scale. Requiring that the age of the universe exceed 10 Gyr implies that Omega approximately greater than 0.25, while requiring that from the last-scattering term in the Sachs-Wolfe formula, large-angle anisotropies come primarily from the decay of potential fluctuations at z approximately less than 1/Omega. Thus, if the universe is open, COBE has been detecting temperature fluctuations produced at moderate redshift rather than at z approximately 1300.

Kamionkowski, Marc↗

NASA NTRS1995

The use of optimal polarizations for studying the microphysics of precipitation: Nonattenuating wavelengths

The objective of this work is to explore relationships between the microphysical properties of precipitation and optimal polarizations. The dependence of three optimal polarization parameters (asymmetry ratio A, optimal tilt tau(sub op), and optimal ellipticity epsilon(sub op) on the reflectivity-weighted mean drop shape, mean canting angle, and standard deviation of a Gaussian canting angle distribution is studied. This is accomplished by using computer simulations that provide the rms scattering matrix for an ensemble of canted drops with a prescribed two-parameter canting angle distribution. Also examined are the effects of propagation on the polarization parameters for nonattenuating wavelengths. The asymmetry ratio A is simply the ratio of the maximal to minimal total backscattered energy (ratio of the largest and smallest eigenvalue of the Graves power matrix G is identically equal to S dagger S). Similar to Z(sub DR), this ratio decreases with increasing mean axial ratio, but unlike Z(sub DR), it is not affected by canting (for a single drop). The dependence of A on the reflectivity-weighted mean drop shape is examined, and a power-law relationship similar to that which exists for Z(sub DR) is established. The asymmetry ratio A can be regarded as a generalization of Z(sub DR) because it requires only a measurement of linear depolarization ratio (in addition to Z(sub DR)), is independent of the propagation phase, and is less sensitive to canting. In a similar manner, the dependence of optimal ellipticity and tilt on the microphysical parameters is studied. In particular, it appears that the rms tilt of the optimal polarization ellipse is proportional to the variance of the canting angle distribution. Several other promising relationships between optimal polarizations and the microphysical variables of an ensemble of hydrometeors are also discussed.

Kwiatkowski, John M.↗

NASA NTRS1993

Linear static structural and vibration analysis on high-performance computers

Parallel computers offer the oppurtunity to significantly reduce the computation time necessary to analyze large-scale aerospace structures. This paper presents algorithms developed for and implemented on massively-parallel computers hereafter referred to as Scalable High-Performance Computers (SHPC), for the most computationally intensive tasks involved in structural analysis, namely, generation and assembly of system matrices, solution of systems of equations and calculation of the eigenvalues and eigenvectors. Results on SHPC are presented for large-scale structural problems (i.e. models for High-Speed Civil Transport). The goal of this research is to develop a new, efficient technique which extends structural analysis to SHPC and makes large-scale structural analyses tractable.

Baddourah, M. A.↗

NASA NTRS1994

Signatures of bulge triaxiality from kinematics in Baade's window

We study a sample of 62 Baade's Window, (l,b) = (1, -4)deg, K giants that have published proper motions, radial velocity, and metallicity. Using R(sub 0) = 8 kpc, we construct the velocity ellipsoids, namely the 3x3 velocity dispersion tensors, for the metal rich stars ((Fe/H) greater than or equal to 0) and metal poor stars ((Fe/H) less than or equal to -0.2). After diagonalizing the tensor, we find a vertex deviation characteristic of a nonaxisymmetric system. Eigenvalues for the two velocity ellipsoids (sigma(sub 1), sigma(sub 2), sigma(sub 3)) are (126, 89, 65) +/- 13 km/s for the metal rich sample and (154, 77, 83) +/- 25 km/s for the metal poor sample with their long axes pointing to two nearly perpendicular directions (l(sub v), b(sub v)) = (-65 +/- 9 deg, +14 +/- 9 deg) and (l(sub v), b(sub v)) = (25 +/- 14 deg, -11 +/- 14 deg), respectively. The vertex deviations of the velocity ellipsoids cannot be consistently explained by any oblate model. We are able to reject the hypothesis that the metal poor and metal rich populations are drawn from the same distribution at better than the 97% confidence level. We populate orbits in a realistic bar potential with a Gaussian velocity distribution, allowing us to simulate and interpret observations. We conclude that the data are consistent with a triaxial bulge pointing towards (l,b) with l less than 0 deg and b = 0 deg as suggested by earlier work on gas dynamics and the observed light distribution. We also predict that low latitude (absolute value of b less than or equal to 4 deg) bulge fields should show the vertex deviation more strongly and would therefore be the best locations for future proper motion studies. In the classification scheme of Athanassoula et al. (1983) the metal rich stars appear to occupy the B-family orbits which rotate in the prograde sense in the rest frame and have boxy shapes that are aligned with and supporting the bar. The metal poor stars in the sample lag behind the metal rich bulge and appear to occupy R-family orbits which rotate in the retrograde sense in the rest frame. They have nearly round loop shapes and are aligned perpendicularly to the bar, hence limit the triaxiality of the bar potential. The correlations between the metallicity and the orbit families can develop if the bulge forms dissipatively on a sufficiently long time scale. However, it is difficult to explain such correlations if most stars in the inner Galaxy form during the violent relaxation phase.

Zhao, Hongsheng↗

NASA NTRS1994

Buckling of imperfect, anisotropic, ring-stiffened cylinders under combined loads

The objective of this study is to develop an anlysis to predict buckling loads of ring-stiffened anisotropic cylinders subject to axial compression, torsion, and internal pressure. This structure is modeled as a branched shell. A nonlinear axisymmetric prebuckling equilibrium state is assumed which is amenable to an exact solution within each branch. Axisymmetric geometric imperfections are included. Buckling displacements are represented by a Fourier series in the circumferential coordination and the finite-element method in the axial coordinate. Application of the Trefftz criterion to the second variation of the total potential energy leads to a nonlinear eigenvalue problem for the buckling load and mode. Results are presented for both unstiffened and ring-stiffened cylinders in the form of buckling interaction diagrams. Imperfections can cause an unexpected buckling mode in the ring web which would not occur for the perfect structure, and pressurization diminishes the benefit of adding rings to the unstiffened shell to increase the buckling load. The implementation of the analysis methodology into a structural sizing algorithm is discussed.

Ley, Robert P.↗

NASA NTRS1995

Temporal variation of meandering intensity and domain-wide lateral oscillations of the Gulf Stream

The path of the Gulf Stream exhibits two modes of variability: wavelike spatial meanders associated with instability processes and large-sale lateral shifts of the path presumably due to atmospheric forcing. The objectives of this study are to examine the temporal variation of the intensity of spatial meandering in the stream, to characterize large-scale lateral oscillations in the stream's path, and to study the correlation betwen these two dynamically distinct modes of variability. The data used for this analysis are path displacemets ofthe Gulf Stream between 75 deg and 60 deg W obtained from AVHRR-derived (Advanced Very High Resolution Radiometer) infrared images for the period April 1982 through December 1989. Meandering intensity, measured by the spatial root-mean-sqaure displacement of the stream path, displays a 9-month dominant periodicity which is persistent through the study period. The 9-month fluctuation in meandering intensity may be related to the interaction of Rosseby waves with the stream. Interannual variation of meandering intensity is also found to be significant, with meandering being mich more intense during 1985 than it was in 1987. Annual variation, however,is weak and not well-defined.The spatially averaged position of the stream, which reflects nonmeandering large-scale lateral oscillations of the stream path, is dominated by an annual cycle. On average, the mean position is farthest north in November and farthest south in April. The first empirical orthogonal function mode of the space-time path displacements represents lateral oscillatins that are in-phase over the space-time domain. Interannual oscillations are also observed and are found to be weaker than the annual oscillation. The eigenvalue of the first mode indicates that about 21.5% of the total space-time variability of the stream path can be attibuted to domain-wide lateral oscillation. The correlation between meandering intensity and domain-wide lateral oscillations is very weak.

Lee, Tong↗